黔北麻羊毛色差异主要调控基因及其通路分析

【目的】明确调控黔北麻羊毛色差异的候选基因及其通路,为揭示黔北麻羊特征毛色的形成机制提供理论依据。【方法】以1岁黔北麻羊公羊为研究对象,分别采集其肩部黑色羊毛下的皮肤组织块和腹部白色羊毛下的皮肤组织块,以RNA为模板合成双链cDNA,经末端修复等处理后进行PCR富集构建cDNA文库,采用Illumina HiSeqTM...     

藏绵羊不同毛色皮肤组织miRNA表达谱及靶基因分析

【目的】藏绵羊又称藏系绵羊,是中国三大粗毛羊品种之一,主要分布在西藏及其毗邻的高寒牧区,如青海、甘肃、四川、云南等地。藏绵羊由于其被毛纤维粗长,毛质细腻且保暖性好,是制造藏式地毯的上好原料。而毛色作为一种重要的经济性状,同时制约羊毛的经济价值。然而,目前对于绵羊毛色的分子调控机制尚不明确,以藏绵羊为研究对象,对其不同颜色皮肤组织（黑色和白色）进行转录组测序,旨在探讨miRNA在藏绵羊不同毛色皮肤组织转录后水平中发挥的作用及其可能参与的调控通路。【方法】选取2只具有黑白花毛色的健康藏绵羊,减去体侧被毛使皮肤完全裸露并用75%酒精消毒处理,屠宰后快速切取皮肤组织并保存在组织样品保存液中防止RNA降解,提取RNA后通过miRNA测序和生物信息学分析,获得藏绵羊不同毛色（黑色和白色）皮肤组织miRNA表达谱,筛选组织差异表达miRNA并预测相关靶基因。通过靶基因分析,获得与调控藏绵羊毛色性状可能相关的信号通路。【结果】黑色和白色皮肤组织共获得85.76兆条原始读长,经过滤后得到85.08兆条clean reads;共鉴定出334个已知的miRNA和59个新的miRNA;从中获得23个差异表达miRNA,其中14个差异表达miRNA在白色皮肤组织中表达显著上调,9个表达下调。miR-2284b和miR-744仅在藏绵羊白色皮肤中表达,而miR-23b、miR-411a-5p、miR-30c、miR-423-3p和miR-324-5p则仅在藏绵羊黑色皮肤中表达,但表达量相对较低。表达量最高的miRNA为miR-10a,而倍数差异最大的为miR-23b。其中,miR-10a可参与调节多种信号通路,其作用涉及增殖、分化、凋亡、细胞粘附等各个方面,目前并无调控绵羊毛色方面的相关报道,但其功能及作用机制的研究仍在不断扩展。miR-23b则与Wnt、Notch等信号通路有关,这些信号通路中的部分基因也和黑色素的生成有密切关系。结合倍数差异和miRNA表达量分析,miR-411a-5p、miR103和miR-200b可能也和毛色调控密切相关。本研究共预测出981个靶基因可能参与毛色调控,多数基因和WNT信号通路、MAPK通路、EDNRB信号通路及黑素原生成通路（melanogenesis pathway）相关,其中黑素原生成通路显示和WNT信号通路、KIT信号通路及EDNRB信号通路相互关联。黑素原生成通路显示藏绵羊不同毛色可能最终通过MITF基因调控下游的TYR基因（酪氨酸酶基因）影响黑色素的合成。同时,筛选7个差异表达miRNA进行定量验证,结果表明,5个miRNA在黑色皮肤组织中上调表达,2个miRNA在白色组织中上调表达,定量结果与测序结果一致。【结论】获得了藏绵羊皮肤组织miRNA表达谱,并通过生物信息学分析得到可能和调控毛色性状相关的miRNA和信号通路。这些结果表明,藏绵羊毛色性状可能受到多种miRNA的调控并涉及多个信号通路,这有助于更进一步了解毛色转录后水平的调控过程,并对进一步的功能验证奠定基础。     

持续感染口蹄疫病毒VP1蛋白对细胞转录组的影响

【目的】通过转录组测序（RNA-Seq）方法分析持续感染口蹄疫病毒结构蛋白VP1诱导的PK-15细胞基因表达谱的变化,探索VP1蛋白对宿主细胞信号通路的影响。为后续开展VP1影响FMDV持续感染相关机制研究提供基础数据和理论支持。【方法】PK-15细胞分别转染pCAGGS-HA（对照载体）、pCAGGS-HA-Pi-VP1质粒,转染36 h后,Western blot检测VP1蛋白表达,同时提取两组细胞总RNA进行转录组测序（每组3个重复,2组6个样本）。将测序所得的序列进行过滤比对,获得差异表达基因后,对差异表达基因进行了生物信息学分析。采用qRT-PCR验证关键差异表达基因的表达。【结果】RNA-Seq共鉴定出3 376个差异表达基因（DEGs）|log2(Fold Change)|>0和padj<0.05,其中上调基因1 744个,下调基因1 632个。许多参与免疫反应和炎症的效应基因差异表达。GO和KEGG富集分析显示,差异基因GO条目绝大部分与核糖体有关,而KEGG富集中涉及了与天然免疫相关的信号通路,如TNF信号通路、Toll受体信号通路。【结论】pCAGGS-...     

鸡白痢沙门氏菌感染对雏鸡脾脏miRNA表达谱的影响

【目的】鉴定感染鸡白痢沙门氏菌后雏鸡脾脏差异表达的miRNA，为阐明miRNA在鸡白痢沙门氏菌感染中的调控机制提供理论基础。【方法】采集感染和未感染鸡白痢沙门氏菌的SPF雏鸡脾脏，提取脾脏RNA，构建6个miRNA文库；利用Illumina Hiseq 2500测序技术和生物信息学分析筛选差异表达miRNA，通过TargetScan和miRanda算法预测差异表达miRNA的靶基因，并进行GO功能富集和KEGG通路富集分析，最后利用RT-qPCR方法验证测序结果。【结果】鉴定到29个差异表达miRNAs，其中15个显著上调、14个显著下调。GO分析表明：脂多糖介导的信号通路正调控以及NIK/NF-κB介导的信号通路负调控等免疫相关生物学过程显著富集；KEGG分析显示：差异表达的miRNA主要参与溶酶体和内吞作用等免疫相关通路。高通量测序结果与RT-qPCR结果一致。【结论】成功鉴定到鸡白痢沙门氏菌感染雏鸡脾脏miRNA表达谱，获得29个鸡白痢沙门氏菌感染潜在相关miRNAs。     

枯草芽孢杆菌NCD-2对棉花根系分泌物L-脯氨酸响应的转录-蛋白质组学联合分析

【目的】棉花根系分泌物L-脯氨酸是影响生防菌株定殖的关键因素之一,前期研究发现L-脯氨酸能够提高枯草芽孢杆菌（Bacillus subtilis）NCD-2生物膜形成能力。通过高通量测序技术分析L-脯氨酸调控枯草芽孢杆菌NCD-2生物膜形成和生防潜力相关的基因,为深入认识棉花根系分泌物与生防菌株的分子互作关系打下基础。【方法】以枯草芽孢杆菌NCD-2菌株为供试材料,外源添加浓度为10 mg·mL^-1的L-脯氨酸共培养24 h后,分别进行转录组（RNA-seq）和同位素标记相对定量蛋白质组技术（iTRAQ）分析,对所获得的转录组和蛋白质组数据进行生物信息学分析,利用实时荧光定量PCR（RT-qPCR）验证不同代谢通路中部分差异基因的表达情况。 【结果】转录组分析发现,L-脯氨酸和NCD-2菌株共培养后,获得1 071个差异表达基因（DEG）,其中602个基因上调,469个基因下调。GO分析发现在生物过程、细胞组分和分子功能方面分别有49、14和30个功能条目显著富集。KEGG代谢通路主要富集在化合物代谢、鞭毛组装、细菌运动和趋化作用。蛋白质组学分析发现,筛选到211个差异表达蛋白（DEP）,其中118个蛋白上调,93个蛋白下调。GO分析发现在生物过程和分子功能方面分别有13和8个功能条目显著富集。KEGG代谢通路主要富集在氨基酸代谢、碳水化合物类代谢、鞭毛组装和ABC转运蛋白。进一步转录-蛋白质组学联合分析发现,在L-脯氨酸作用下,检测到共有的显著差异表达基因（或蛋白）112个,其中38个基因（或蛋白）下调,74个基因（或蛋白）上调。GO功能显著富集主要集中在营养库活性、催化活性、细胞膜、定位、细胞脂质代谢过程、氧化还原过程、sigma因子活性、转运活性、芽孢形成9个方面。KEGG代谢通路主要富集在能量代谢、ABC转运蛋白、抗生素生物合成、鞭毛组装、运动或趋化作用和双组分系统方面。RT­qPCR验证了26个差异表达显著基因,结果发现在表达量上存在一定的差异,但表达趋势与RNA­seq和iTRAQ组学分析结果基本一致。【结论】棉花根系分泌物L-脯氨酸与枯草芽孢杆菌NCD-2之间的互作存在一个复杂的生物过程,依赖于不同代谢通路网络中的多个基因。双组分系统、抗生素生物合成、能量代谢、运动或趋化作用、鞭毛组装和ABC转运蛋白通路中的差异基因（或蛋白）可能在棉花根系分泌物与枯草芽孢杆菌互作过程方面发挥着重要作用。     

西方蜜蜂工蜂不同虫态发育的转录组学分析

【目的】基于转录组学对西方蜜蜂工蜂不同虫态间的差异表达基因（DEGs）进行筛选和功能注释分析,揭示与工蜂生长发育相关的信号通路,为深入解析工蜂生长发育的分子调控机理提供基础数据。【方法】以西方蜜蜂工蜂的3日龄幼虫、1日龄蛹和1日龄羽化工蜂3个虫态为研究对象,利用llumina NovaSeq 6000平台进行转录组测序,采用DESeq2筛选不同虫态样品间的表达差异基因,然后分别进行GO功能注释分析及KEGG信号通路富集分析,并通过实时荧光定量PCR进行验证。【结果】经转录组测序,在西方蜜蜂工蜂3日龄幼虫与1日龄蛹间筛选出4823个差异表达基因（51.86%上调,48.14%下调）,在1日龄蛹与1日龄羽化工蜂间筛选出3295个差异表达基因（57.51%上调,42.49%下调）,在3日龄幼虫与1日龄羽化工蜂间筛选出5267个差异表达基因（52.95%上调,47.05%下调）。GO功能注释分析结果显示,3日龄幼虫与1日龄蛹间的差异表达基因注释到43个GO功能条目,1日龄蛹与1日龄羽化工蜂间的差异表达基因注释到45个GO功能条目,3日龄幼虫与1日龄羽化工蜂间的差异表达基因注释到44个GO功能条目,主要涉及细胞过程、细胞部分及结合等。KEGG信号通路富集分析发现,3日龄幼虫与1日龄蛹间有2905个差异表达基因富集到332条KEGG信号通路上,其中17条KEGG信号通路呈显著富集,涉及核糖体、氧化磷酸化和昆虫激素生物合成等;1日龄蛹与1日龄羽化工蜂间有1644个差异表达基因富集到331条KEGG信号通路上,其中45条KEGG信号通路呈显著富集,涉及氧化磷酸化、生热作用和胰岛素分泌等;3日龄幼虫与1日龄羽化工蜂间有2958个差异表达基因富集到337条KEGG信号通路上,其中14条KEGG信号通路呈显著富集,涉及核糖体、蛋白酶体和胰岛素分泌等。6个随机挑选差异表达基因的实时荧光定量PCR检测结果与转录组测序结果相符,即转录组测序结果可靠。【结论】昆虫激素生物合成通路相关差异表达基因调控与西方蜜蜂工蜂各虫态JH滴度变化规律一致,氧化磷酸化信号通路则与各虫态的营养摄入和活动行为相关,而胰岛素分泌通路涉及各虫态的营养调控、脂肪体合成及细胞凋亡。可见,昆虫激素生物合成、胰岛素分泌和氧化磷酸化3种信号通路在西方蜜蜂工蜂幼虫、蛹和成虫的发育调控中发挥着重要作用。     

转录组测序揭示不同生长速度番鸭胸肌组织的基因表达差异

【目的】利用转录组测序技术分析不同生长速度番鸭的基因表达差异，挖掘影响番鸭生长性能的基因与信号通路，为阐明生长性能差异的遗传机制提供理论基础。【方法】选取不同生长速度两尾样番鸭的胸肌组织进行转录组测序，采用无参考基因组的分析手段对测序结果进行分析，获得差异表达基因（Differentially expressed genes，DEGs）。利用 Blast2 GO 软件对 DEGs 进行功能富集分析，利用 KAAS 网站进行 DEGs 的通路富集分析。【结果】转录组测序分析结果显示，快速生长型番鸭与缓慢生长型番鸭相比共有 186 个显著 DEGs，其中 49 个表达量上调，137 个表达量下调。这些 DEGs 主要富集在磷脂酰肌醇三激酶 - 丝氨酸 / 苏氨酸激酶（PI3K-AKT）信号通路、丝裂原活化蛋白激酶（MAPK）信号通路和腺苷酸激活蛋白激酶（AMPK）信号通路。【结论】PI3K-AKT、MAPK、AMPK 这 3 条信号通路是控制番鸭机体生长发育的重要调节信号通路。     

基于TMT蛋白组学及生物信息学分析牦牛抗冻差异蛋白

【目的】从蛋白水平阐明牦牛抗寒性能机理,并进一步从营养学角度提高其代谢性能,为牦牛有效抵御外界恶劣气候条件提供科学依据。【方法】应用TMT蛋白组学技术对寒冷季节（1月）和温暖季节（8月）的牦牛抗冻蛋白进行挖掘,并对鉴定到的牦牛抗冻蛋白进行亚细胞定位、结构域、GO功能富集、KEGG信号通路注释、蛋白相互作用等生物信息学分析。【结果】从牦牛耳组织中共鉴定获得21856个肽段（Peptide）,其中特有肽段（Unique peptide）序列为18452个,定量获得4519个蛋白,最终筛选出144个差异蛋白,其中上调蛋白89个、下调蛋白55个。144个牦牛抗冻差异蛋白亚细胞定位到7个条目上,分别是细胞核蛋白56个、细胞质蛋白51个、质膜蛋白24个、细胞外蛋白23个、线粒体蛋白18个、细胞骨架蛋白1个和溶酶体蛋白1个;共鉴定到194个结构域。GO功能富集分析结果显示,生物过程主要富集到细胞过程蛋白79个、代谢过程蛋白70个和生物调控蛋白42个等,分子功能主要富集到结合功能蛋白75个和催化活性蛋白64个等,细胞组分主要富集到细胞部分蛋白89个和细胞蛋白89个等。144个牦牛抗冻差异蛋白在KEGG数据库中注释到205条KEGG信号通路,主要涉及核糖体、氮代谢、胞质DNA感受、动物体内生热作用、氧化磷酸化、白细胞介素-17及钙离子信号等通路。牦牛抗冻差异蛋白相互作用网络分析发现L8IHE5的关联度最高,且冷诱导RNA结合蛋白（CIRP）和HSP70结合蛋白在蛋白相互作用网络中具有更多的相互作用关系。【结论】基于TMT蛋白组学对牦牛抗冻差异蛋白进行挖掘,结果鉴定获得144个抗冻差异蛋白（上调蛋白89个,下调蛋白55个）,其中CIRP和HSP70在冷应激条件下呈上调趋势,能促使牦牛肌体适应低温环境,可作为牦牛抗冻性育种的候选分子标记。     

POMC在不同毛色绵羊皮肤中的表达与定位

【目的】对不同毛色绵羊皮肤中的阿黑皮素原（proopiomelanocortin,POMC）进行定位和表达差异分析,确定POMC与毛色形成的关系,进一步丰富有关哺乳动物毛色形成机制的理论知识。【方法】选取年龄相近的雄性黑色和白色成年绵羊各3只,用取皮器分别在臀部采集皮肤组织3块,对其中一块制作石蜡切片后采用免疫组织化学法确定POMC在不同绵羊皮肤中的定位,另两块置于液氮中冻存后分别提取皮肤总RNA和总蛋白,总RNA经反转录成cDNA后,经NCBI分析比对后设计POMC特异性PCR引物,采用QRT-PCR方法检测POMC基因在mRNA水平的表达差异;总蛋白经SDS-PAGE凝胶电泳后,采用Western blotting方法加入POMC多克隆抗体检测POMC在蛋白水平的表达差异,得出的数据用SPSS19.0软件进行差异显著性分析,探讨POMC对毛色形成机制的调控作用。【结果】免疫组织化学结果显示POMC在黑色和白色绵羊皮肤中均有表达,且表达部位都位于毛囊根鞘周围和表皮的皮肤角质形成层细胞内;QRT-PCR结果显示POMC在黑色绵羊皮肤中mRNA的相对表达量为2.5004±0.2577^**,而在白色绵羊皮肤中POMC的mRNA的相对表达量为1.0032±0.0350,黑色绵羊皮肤中POMC的mRNA的相对表达量为白色绵羊皮肤的2.5倍,二者表达差异极显著（P＜0.01）;Western boltting结果显示黑色绵羊皮肤POMC蛋白水平相对表达量为1.5253±0.0426^*,而白色皮肤中POMC蛋白水平相对表达量为1.0005±0.0180,黑色皮肤POMC蛋白表达量是白色皮肤的1.55倍（P＜0.05）,两者差异显著。【结论】试验表明POMC可在黑色和白色绵羊皮肤中正常表达,且其表达位置一致,但POMC在不同毛色绵羊皮肤组织中mRNA和蛋白的表达水平存在一定的差异,表明POMC的表达量变化会在一定程度上影响毛色的生成。     

中国美利奴羊羊毛弯曲相关microRNA的筛选

【目的】 研究miRNA在不同羊毛弯曲中国美利奴羊皮肤组织中的表达差异及其可能参与的调控通路,筛选出与中国美利奴羊羊毛弯曲相关的miRNA,为进一步探究细毛羊羊毛弯曲性状的调控机理提供理论依据。【方法】 运用高通量测序技术对羊毛弯曲频率有极端差异表型值(大弯曲组与小弯曲组)细毛羊个体皮肤的miRNA组进行测序与比较分析。【结果】 共鉴定出425个miRNAs,其中包括143个已知miRNAs和282个新发现的miRNAs。在大弯曲组与小弯曲组共筛选到8个上调和17个下调的miRNAs。对差异miRNA的预测靶基因进行简单的Gene Ontology与KEGG Pathway富集分析。差异表达的miRNA的靶基因主要参与跨膜信号受体活动、信号转导活动、分子转导活动、膜的组成等过程。15 761个靶基因注释到252个信号通路。【结论】 筛选出了25个与细毛羊羊毛弯曲性状相关的候选miRNAs。     

福清山羊与努比亚黑山羊发情期卵巢组织RNA-Seq分析

【目的】 研究福清山羊与努比亚黑山羊发情期卵巢组织转录组差异表达水平,一方面为山羊繁殖性状形成的分子机制提供理论依据,另一方面为利用努比亚黑山羊杂交改良福清山羊提供可加快遗传进展的分子标记。【方法】 利用转录组测序方法对福清山羊和努比亚黑山羊发情期卵巢组织进行研究,筛选品种间的差异表达基因（differentially expressed genes,DEGs）,并对DEGs功能进行注释和若干基因荧光定量PCR（quantitative real-time PCR,qRT-PCR）验证;同时通过与参考基因组比对,分析和筛选测序样品中存在的SNP/InDel。【结果】 6个样品共得到46.68Gb Clean Data,DESeq分析发现了福清山羊和努比亚黑山羊发情期卵巢组织的DEGs149个（包含25个新转录本）,其中表达上调53个,表达下调96个;初步认为输卵管素基因（oviductin,OVN）、类固醇合成快速调节蛋白基因(steroidogenic acute regulatory protein,STAR)、早期生长应答1基因(early growth response 1,EGR1)可作为福清山羊繁殖性能的候选基因。149个DEGs中的108个基因能被GO（gene ontology）数据库注释,30个DEGs能够被COG（Cluster of Orthologous Groups of proteins）数据库注释,91个DEGs能够被KEGG（kyoto encyclopediaof genes and genomes）数据库注释。KEGG通路分析表明DEGs共富集到21条信号通路中,6条通路显著富集。经过进一步的与参考基因组序列比对分析,共发掘1 506个新转录本。经qRT-PCR验证,所选基因（转录本）表达变化模式与转录组测序结果一致,表明测序结果可靠。【结论】 在转录组水平上筛选出了福清山羊和努比亚黑山羊发情期卵巢组织的149个DEGs,发掘了1 506个新转录本,初步揭示了OVN、STAR和EGR1在山羊繁殖过程中可能发挥重要作用,为进一步探索山羊繁殖性状相关机理提供参考依据。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

抗氟苯虫酰胺小菜蛾差异表达基因及其通路

【目的】随着氟苯虫酰胺等双酰胺类杀虫剂的广泛应用,小菜蛾（Plutella xylostella）对该类药剂的抗性愈发突出。研究旨在从转录组水平研究小菜蛾对氟苯虫酰胺的抗性机制,明确小菜蛾对氟苯虫酰胺产生抗性的关键基因及通路,为揭示小菜蛾对氟苯虫酰胺的抗性机制打下基础。【方法】以室内筛选的小菜蛾抗氟苯虫酰胺品系（Rh28）、田间抗性种群（Rz36）和敏感品系（S）为材料,通过转录组测序技术（RNA-Seq）,获得抗性小菜蛾品系/种群在转录水平上的差异表达基因及抗性显著上调表达基因,采用基因本体（GO）显著性富集分析及京都基因与基因组百科全书数据库（KEGG）富集分析,确定差异表达基因具有的主要生物学功能及其参与的主要生化代谢和信号转导途径。【结果】通过将不同抗性品系/种群与敏感品系的RNA-Seq数据对比,获得数量不等的差异表达基因。将差异基因富集的生物学过程（biological process）和信号通路（pathway）进行GO分析,发现其主要集中于对刺激物的反应（response to stimulus）、催化活性（catalytic activity）和代谢过程（metabolic process）等条目中;而KEGG分析发现,差异基因主要集中于代谢通路（metabolic pathway）中,基因数量占比最高,为17.84%。通过统计分析进一步获得了抗性显著上调表达基因,并对其中218个基因开展了GO分析,其仍然主要集中在代谢过程、对刺激物的反应和生物调节（biological regulation）等条目中;通过对抗性显著上调基因进行KEGG分析,发现其仍主要集中在代谢通路中,其中谷胱甘肽-S-转移酶、细胞色素P450解毒酶及参与昆虫多种生理反应的热激蛋白（Hsp）超基因家族等基因的表达量均显著上调。通过聚类热图分析发现,显著表达上调的基因较多集中在内质网上蛋白质的加工（protein processing in endoplasmic reticulum）、酪氨酸代谢（tyrosine metabolism）、咖啡因代谢（caffeine metabolism）和Wnt信号通路（Wnt signaling pathway）中。【结论】获得了小菜蛾抗氟苯虫酰胺差异表达基因及抗性显著上调表达基因,其主要富集于代谢过程、应激反应及对刺激物的反应等条目中,这些基因的相互协同调控是小菜蛾对氟苯虫酰胺产生抗性的重要机制。