肺形侧耳低温胁迫时期的转录组分析

本文以肺形侧耳栽培菌株X57为研究对象,利用高通量RNA测序技术对4℃低温处理0h、6h和12h后的肺形侧耳进行基因表达分析来探讨肺形侧耳低温应答机制。分析结果筛选出低温处理6h差异表达基因742个,其中上调表达基因374个,下调表达基因368个;低温处理12h差异表达基因1 489个,其中上调表达基因占53%,下调表达基因占47%。Gene Ontology（GO）功能聚类分析表明,差异表达基因主要富集在结合、催化分子功能组和代谢过程、细胞过程生物学过程中。KEGG功能富集分析结果显示,差异基因主要富集在氨基酸、核糖体和类固醇生物合成,及氮类物质代谢的通路上,富集到与低温胁迫相关通路MAPK signaling pathway-yeast上的基因表达随低温处理时间的增加呈上调趋势。已报道HOG-MAPK通路是MAPK途径研究较为明确、信号传递单一的真菌低温胁迫中重要的通路,本文运用生物信息学软件构建肺形侧耳的HOG-MAPK通路,并利用荧光定量PCR对相关基因表达进行验证,结果显示以低温处理0h为参照,随着低温胁迫时间增加通路上大部分基因的表达量持续上升,且差异显著,与RNA‐Seq分析结果一致。本文通过全面分析肺形侧耳低温胁迫时期基因表达情况,为进一步研究肺形侧耳低温胁迫相关重要基因的功能鉴定与信号通路调控机理提供基础。     

基于转录组学的蓖麻耐盐基因的挖掘

为获取蓖麻耐盐基因的序列信息,挖掘盐胁迫下差异表达基因及相关代谢途径,本研究以盐胁迫（300 mmol/L NaCl）处理0、12和24 h后通篦5号的幼苗真叶为试验材料,借助高通量测序技术进行转录组测序分析。结果表明,在盐胁迫12 h和24 h分别有4822和3103个差异表达基因。对共有的1872个差异表达基因进行共表达模式聚类分析,发现这些基因共有3种表达模式。KEGG代谢通路分析结果显示,缬氨酸、亮氨酸和异亮氨酸降解（ko00280）、植物昼夜节律调控（ko04712）以及淀粉和蔗糖代谢（ko00500）3个通路在盐胁迫适应过程中显著富集;GO功能富集分析结果显示,多数差异表达基因被富集到生物过程中,其中细胞进程（GO:0009987）和响应非生物胁迫（GO:0009628）过程富集到差异表达基因的数目最多。另外,共有19个转录因子参与蓖麻的盐胁迫响应。植物激素信号转导通路中有42个差异表达基因,其中有97.6%的基因分别在12 h和24 h是上调表达的。此外,还筛选出包括光合作用途径、抗氧化调节、Na+、K+和Ca2+转运相关参与蓖麻幼苗盐胁迫的差异表达基因。qRT-PCR结...     

甘薯响应冷胁迫的转录组分析

甘薯是世界第七大粮食作物,虽然甘薯对盐、旱等胁迫的抗逆性较强,但因其起源热带而对低温耐受性差,目前低温已经成为影响甘薯产量及限制甘薯种植面积的重要非生物胁迫。本研究以耐冷品种辽薯36为试验材料,在低温处理0 h、6 h、24 h和48 h后,利用转录组测序技术筛选与冷胁迫紧密相关的代谢通路和差异表达基因。结果共获得甘薯转录组数据75.01 Gb,与参考基因组比对效率达76.33%。差异表达基因分析结果表明,甘薯响应冷胁迫的过程中有6282个基因上调表达,3881个基因下调表达,其中283个共有基因在低温处理6 h、24 h和48 h上调表达,26个基因在3个处理中均下调表达。GO功能富集分析发现,差异表达基因被大量富集到蛋白激酶活性、过氧化物酶活性、氧化还原酶活性、蛋白质磷酸化、氨基酸跨膜转运、蛋白质生色团连锁及果胶分解代谢分子功能和生物学过程等分类条目中。KEGG代谢通路富集分析表明,苯丙烷生物合成、淀粉和蔗糖代谢、植物激素信号转导、植物-病原相互作用、谷胱甘肽代谢这5个代谢通路富集了大量的差异表达基因。本研究为甘薯耐冷基因克隆及耐冷机制解析提供了理论依据。     

金钱鱼下丘脑转录组分析及其对雌激素在体注射的响应

为探究XX雌性和XY雄性金钱鱼(Scatophagus argus)下丘脑中的差异表达基因及雌激素(17β-Estradiol,E2)在体注射对XY雄鱼下丘脑中基因表达的影响,开展了转录组测序分析,包括测序数据质控、基因功能注释,差异基因筛选、鉴定和差异基因功能富集分析等,并利用实时定量PCR(Real-time quantitative PCR, qPCR)检测了金钱鱼下丘脑中相关基因的表达。金钱鱼下丘脑转录组共测得Clean reads 275833710,测序质量30(Q30)大于95%, GC含量值大于48%。其中, Ctrl-XX-H vs. Ctrl-XY-H组,共筛选和鉴定了91个差异表达基因,包括36个上调基因和55个下调基因。在Ctrl-XY-H vs. E2-XY-H组,筛选和鉴定了28个差异表达基因,包括11个上调基因和17个下调基因。GO和KEGG富集分析发现,差异表达基因主要显著富集在细胞、单生物过程,膜、膜组分,结合等生物学功能及泛醌和其他萜类-醌生物合成及类固醇激素生物合成通路等。转录组和qPCR结果表明, XX和XY金钱鱼下丘脑中prl、ccna1和hs...     

金针菇单、双核菌丝差异表达基因分析

对金针菇Flammulina velutipes单核菌丝W23的菌丝体以及与L11质配后的双核菌丝H1123菌丝体进行了转录组测序,以本实验室已获得的W23基因组为参考基因组研究两样本间差异基因,并对这些差异基因进行了GO功能和Pathway显著性富集分析。差异基因分析显示,两个样本中共有显著性差异表达的基因3 504个,其中在双核菌丝中上调、下调的基因数分别为2 151和1 353个。研究发现差异表达基因含有很多的转录因子基因、蛋白激酶以及WD40 repeat-like蛋白。Gene Ontology（GO）功能分析结果表明,extracellular region和membrane-enclosed lumen条目下的差异基因全部为上调表达,而envelope下的差异基因全部为下调表达,以利于双核菌丝分裂时锁状联合的形成而便于核的迁移。Pathway 功能富集分析结果表明,脂肪酸、氨基酸以及大部分糖类合成相关基因具有比较活跃的上调表达。说明双核菌丝主要进行营养物质的富集,为下一步在合适条件下分化成原基,进入生殖生长阶段储备物质基础。     

苗期水稻响应褐飞虱取食的基因差异表达分析

【目的】褐飞虱Nilaparvata lugens (St?l)是水稻的重要害虫之一,主要刺吸水稻韧皮部汁液为害,对水稻生产造成严重的产量和经济损失。为研究水稻响应褐飞虱取食的分子机制,对褐飞虱取食6 h后的苗期水稻进行转录组测序及基因差异表达分析。【方法】采用illumina二代测序技术获得褐飞虱取食前后水稻组织的转录组数据,利用RSEM软件进行基因表达定量和DEseq2进行差异表达分析;从差异表达基因中随机选取20个基因采用荧光定量PCR技术进行验证;采用GeneMerge软件对差异表达基因进行KEGG和GO富集分析。【结果】褐飞虱取食后,水稻转录组中的1 104个基因出现了差异表达,其中435个基因表达上调,669个基因表达下调。荧光定量PCR结果显示,20个差异表达基因中18个基因的表达变化趋势和测序结果一致,证明了转录组分析结果可靠。GO和KEGG富集分析表明,表达上调基因主要与水稻氧化应激、海藻糖合成及次生化合物代谢有关,显著富集在14个KEGG通路和30个GO功能分类中;而表达下调基因主要参与水稻纤维素、蛋白质及脂肪酸合成过程,显著富集在29个KEGG通路和26个GO功能分类。在差异表达基因中,分别有61个转录因子和13个水杨酸和茉莉酸信号通路相关基因。【结论】褐飞虱取食激发了水稻的应激反应和保护机制,同时还降低了营养合成的过程,是飞虱为害造成水稻减产的原因之一。本研究初步揭示了苗期水稻响应褐飞虱取食的差异表达基因,为研究水稻-褐飞虱互作机制以及褐飞虱抗性水稻品种培育提供了参考和依据。     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

盐度胁迫对幼体大海马基因转录表达的影响

为探究大海马(Hippocampus kuda Bleeker)幼体在高盐、低盐胁迫条件下的基因表达水平的变化规律,对实验条件下的大海马幼体的肝脏样品进行了转录组测序。对照组(CK, 25‰)、高盐(HS-test, 31‰)和低盐(LS-test, 17‰)胁迫组共获得71794个单基因簇(Unigenes), N50为1780 bp,平均长度为820.71 bp。高盐胁迫组与对照组比较,共获得2740个差异表达基因(DEGs),其中495个DEGs上调, 2245个DEGs下调;与对照组相比,低盐胁迫组共获得3715个DEGs,其中1854个DEGs上调, 1861个DEGs下调。高/低盐胁迫组DEGs经KEGG数据库富集发现,高/低盐度胁迫均能导致大海马幼体体内氨基酸代谢、免疫代谢、能量和脂肪酸代谢相关基因受到影响。其中,低盐胁迫时能量代谢和氨基酸代谢的相关基因显著上调,高盐胁迫时脂肪酸代谢的相关基因显著下调,而高/低盐胁迫时免疫代谢的相关基因都显著上调。从经过盐度胁迫的大海马幼体的肝脏转录组中分别筛选到免疫相关基因Gst、Hsp70、Hsp90、Sod、Bcl-2、Gadd45...     

基于转录组测序的花烟草低温胁迫响应转录因子挖掘

为了鉴定参与花烟草低温胁迫的转录因子,对花烟草幼苗进行了4℃低温处理,并在处理后12 h采集其幼苗样本,提取总RNA后,采用高通量测序技术,进行了转录组测序。在对差异表达基因进行GO及KEGG分析的基础上,对参与其中的转录因子进行了挖掘,并采用qRT-PCR的方法,对转录组测序的结果进行了验证。结果表明,低温处理后,花烟草基因表达量变化在2倍以上的基因有8388个(P<0.01),其中,上调表达4229个,下调表达4159个。这些差异表达基因的功能归类于生物过程、细胞组分及分子功能3大类69个GO条目,并显著富集在40条KEGG代谢通路中。同时,在低温胁迫下,花烟草有118个转录因子的表达发生了显著改变,其中,上调表达82个,下调表达36个。这些转录因子属于28个家族,其中,数量最多的为NAC家族19个,其次为ERF家族16个、MYB家族15个、WRKY家族15个。本研究的结果为花烟草低温响应分子机制的研究提供了借鉴。     

阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫低温胁迫转录组分析

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾Hepialus altaicola Wang幼虫基因转录差异及上调表达基因的主要功能类群和参与的主要代谢通路.[方法]采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾幼虫进行转录组测序、组装,利用Blast软件进行数据库比对和基因功能注释,用DEGSeq R软件包分析4℃低温处理与室内适温饲养试虫的差异表达基因,并对上调表达基因进行GO和KEGG代谢途径富集分析.[结果]经序列拼接后共获得100300个unigenes,总长度81600309 bp,平均长度813 bp,N50长度1719 bp.与7大数据库同源比对,共获得34691(34.59％)条unigenes.低温胁迫转录组分析得到11569个差异表达基因(DEGs),7158条基因上调,4411条基因下调.富集到47个GO类群,217个KEGG途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

圆口铜鱼亲本与子代肝脏转录组特征及对比分析

通过Illumina HiseqTM 2000测序平台首次开展了圆口铜鱼(Coreius guichenoti)亲本与子代肝脏转录组测序并对比分析其测序结果。对转录组进行拼接和组装, 共获得80688个Unigene, Unigene的长度主要分布在401—600 bp, 占总Unigene的43.56%。与Nr、GO、COG、KEGG等公共数据库比对并进行功能注释, 经分析圆口铜鱼亲本与子代差异表达基因共2701个, 其中上调1240个, 下调1461个。差异基因表达模式聚类热图显示, 同一样本不同重复基因表达相似。将差异基因映射到代谢通路KEGG数据库进行富集分析, 并对37个KEGG显著富集的代谢通路(P<0.05)构建其基因共表达网络, 结果显示丙酸盐代谢、丙酮酸代谢、原核生物固碳途径、乙醛和二羧酸代谢、脂肪酸代谢、萜类骨架生物合成是KEGG差异表达基因的核心代谢途径, 且这些基因在圆口铜鱼子代中表达量均显著上调。该结果丰富了圆口铜鱼基因组数据, 并首次从分子水平比较了圆口铜鱼亲本和子代的代谢差异, 同时解释了在同一循环水系统中子代不易患病的现象。     

烟草脉斑驳病毒侵染本氏烟的转录组研究

为了研究烟草脉斑驳病毒（TVMV）侵染本氏烟的抗病分子机制,本文建立了TVMV转录组数据库,挖掘抗病相关基因以及代谢和信号通路,采用Illumina Novaseq 6000高通量测序平台对侵染TVMV的本氏烟进行转录组测序,并进行生物信息学分析;利用R package:edge R等软件分析了有关抗病的差异表达基因,对差异表达基因进行基因本体论（GO）及京都基因和基因组百科全书（KEGG）基于途径的通路富集分析。基于TVMV侵染本氏烟转录组的测序结果中共获得4 593个差异表达基因,其中3 564个基因上调表达,1 029个基因下调表达。共筛选出10个抗病相关的差异表达显著基因,6个上调,4个下调。GO数据库中注释到生物学过程、细胞组分及分子功能等三大类共50个功能组。其中,在第三大类分子功能中,蛋白质结合功能以及ATP结合功能所涉及的差异表达基因最显著,分别为501个和453个差异基因。KEGG共富集20条通路,注释到核糖体途径的差异基因富集程度最高,基因最显著,差异表达基因数为307个。蛋白质结合、ATP结合、核糖体、真核生物核糖体的生物反应、DNA复制（DNA replicat...     

NaCl胁迫下黑果枸杞转录组测序分析

研究黑果枸杞在不同浓度盐胁迫下基因表达谱变化情况,为进一步研究黑果枸杞抗盐分子机制奠定研究基础。对0(CK)、50、250 mmol/L NaCl溶液胁迫的黑果枸杞组培苗的根和叶在胁迫时间为0、1、12 h时分别取样,采用转录组测序(RNASeq)技术进行测序分析。结果表明,转录组测序共产生222.49 Gb原始数据,拼接出Unigenes 86 037条,注释到7大功能数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG)上的Unigenes总数为46 594个,占总Unigenes的54.76%,还有38 929个Unigenes在这些数据库中没有得到注释。通过GO分类和KEGG Pathway富集性分析,分别归于51个GO类别和211条代谢途径。差异表达基因分析显示,黑果枸杞叶片和根的上调基因和下调基因数随着NaCl浓度的增大和处理时间的延长均呈增加趋势,叶片中的上调基因数(7 514)小于下调基因数(9 032),根中的上调基因数(12 347)大于下调基因数(11 559)。在黑果枸杞盐胁迫下转录组中发现28 325个SSR位点,最多的为单核苷酸SSR,占70.47%。综合分析表明,黑果枸杞对盐胁迫的反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是基因调控的主要方式。     

基于cDNA-AFLP技术的玉米种子劣变相关基因分析

以人工控制劣变处理（温度45℃,相对湿度100%）3 d的玉米杂交种‘郑单958’种子为材料,通过cDNA-AFLP技术,共获得约5 400个转录本。对差异片段进行测序分析,最终获得122个种子劣变相关基因序列（DRS）,其中,上调表达74个,下调表达48个,功能涉及初生代谢、次生代谢、蛋白代谢、转录与调控、胁迫响应、信号转导和运输。采用qRT-PCR技术分析了11个上调表达差异基因的表达模式,结果显示不同基因的表达模式存在较大差异,说明不同基因在种子劣变中可能扮演不同角色。     

眼优角蚱两性成体转录组差异分析

为了解蚱类昆虫两性间基因转录水平上的整体差异,本研究以眼优角蚱Eucriotettix oculatus为研究对象,通过高通量转录组(RNA-Seq)测序和分析技术对其雌、雄成体的转录组进行了测定和分析。结果显示,在雌、雄转录组中共获得13 432条转录本,两性间差异表达基因共3 597个,其中雌性相对于雄性显著上调基因1 295个,显著下调基因2 302个。对其中9个差异基因进行实时荧光定量PCR验证,其结果与RNA-Seq分析一致。差异基因显著富集的GO条目大多是与细胞的形成有关,结合、细胞组分和细胞过程等条目相关基因显著高表达。KEGG富集显示差异基因主要是与性腺发育、遗传信息加工、能量代谢和消化等通路有关。研究结果对从分子水平进一步研究眼优角蚱两性差异提供一定的理论基础。     

基于高通量测序的都匀地区福鼎大白种茶树根茎叶分析

为探究茶树中茶多酚等产物代谢途径的相关基因,该研究以贵州都匀地区福鼎大白种茶树的根茎叶为对象,利用高通量测序技术构建茶的转录组数据库并筛选其根茎叶差异表达基因。结果表明:共获得70.88 Gb Clean Data,各样品Clean Data均达到6.33 Gb,Q30碱基百分比在93.22%以上。将Clean Reads与中国种茶树参考基因组进行序列比对,比对效率从87.83%到91.14%。基于比对结果,进行可变剪接预测分析和基因结构优化分析,发掘新基因13 531个,其中10 244个得到功能注释。利用FPKM进行基因表达量分析,根据基因在不同样品中表达量识别差异表达基因。叶与茎的差异基因有5 595个,其中2 769个在茎中上调,2 826个下调,叶与根有9 650个差异基因,5 056个上调,4 594个下调,茎与根中有5 644个差异基因,2 938个上调,2 706个下调,并通过GO和KEGG分析,将差异基因进行功能注释和富集分析。上述结果为揭示都匀地区福鼎大白种茶参与类黄酮、茶氨酸和咖啡碱等代谢途径相关的基因提供了参考,为选育优良品种等提供了理论依据。     

茶树种子发育过程的转录组分析

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的3个时期(幼果期、膨大期、成熟期)的表达差异,探究茶树种子油脂代谢的分子机制。结果表明:(1)经转录组测序、组装后共获得30 940 581个clean reads,经数据合并拼接最终得到36 951条非冗余Unigene序列,其中28 476个Unigene可得到功能注释;在转录本中能够被注释到GO分类的Unigene有11 201条(30.3%),KEGG分析发现共有17 172个基因参与了127个代谢通路。(2)经KEGG通路筛选出14条与脂肪酸代谢相关的通路,且随着茶籽的发育,大部分脂肪酸调控途径相关基因呈下调趋势,其中上调基因数最多的有α-亚麻酸代谢途径和脂肪酸降解途径(有17个基因表达量上调),下调基因数最多的是甘油磷脂代谢途径(有58个基因表达量下调);在茶籽发育幼果期α-亚麻酸代谢途径中表达量上调的基因数超过表达量下调的基因数。(3)研究发现茶籽脂肪酸合成相关的基因涉及14个脂类调控途径,共409条差异基因;随着茶树种子发育到成熟期,上调的差异表达基因数量在减少,下调的差异表达基因数量增加,其中α-亚麻酸途径中的基因PLA2G16、DAD1、pldA、FabF、FabI表达量上调显著,随后表达量下调。(4)qRT-PCR检测结果表明,7个茶树FAD和1个ACP差异表达基因的水平与转录组测序结果基本一致;随着茶籽的发育,基因CsFAD7和Δ6-CsFAD从幼果期、果实膨大期至果实成熟期都为差异下调表达,CsFAD2、CsFAD6和Δ7-CsFAD为差异上调表达,CsFAD8、Δ8-CsFAD和CsACP在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

鲤在低温胁迫下肝胰腺转录组测序分析

检测低温耐受与低温敏感鲤品种低温胁迫下的差异表达基因,探究鲤低温适应的分子应答过程。低温胁迫后对松浦镜鲤及荷包红鲤肝胰腺进行转录组测序,并对表达差异基因进行功能、通路富集分析。通过差异表达分析,共获得10 521个差异表达基因,其中5 246个基因仅在低温耐受品种发生表达变化。功能、通路富集分析将差异基因显著富集在包括糖酵解/糖质新生代谢途径在内的144个代谢途径。选取5个基因进行qRT-PCR验证,定量结果与RNA-seq结果相关性达0.89。研究结果表明低温胁迫会显著改变鲤肝胰腺的转录活动,编码糖酵解/糖质新生途径的多个关键酶基因在低温耐受品种的表达量显著高于低温敏感品种,推测这些基因功能有助于其长期低温适应。     

盐胁迫下两个甜瓜品种转录因子的转录组分析

利用新一代高通量测序手段——转录组测序（RNA-Seq）技术研究在300 mmol.L-1NaCl胁迫下两个甜瓜（Cucumis melo L.）品种‘玉露’和‘冰雪脆’的转录因子基因表达变化。这两个甜瓜品种在叶绿素荧光参数上的差异表明其在盐胁迫下有不同的生理反应。转录组测序结果表明,盐胁迫与对照相比,‘玉露’共有属于19个转录因子家族的56个转录因子基因表达发生变化（在转录水平,属于7个转录因子家族的22个转录因子上调表达,属于14个转录因子家族的34个转录因子下调表达）。‘冰雪脆’有属于20个转录因子家族的47个转录因子基因表达发生变化（在转录水平上,属于5个转录因子家族的17个转录因子上调表达,属于17个转录因子家族的30个转录因子下调表达）。盐胁迫下,两个甜瓜品种差异表达的转录因子既表现特异性,也存在部分重叠。‘玉露’有29个转录因子特异响应,‘冰雪脆’有20个特异响应。盐胁迫响应重叠的转录因子有27个,其中9个上调表达,18个下调表达。采用实时荧光定量PCR对几个转录因子进行了盐胁迫下的表达检测,其趋势与转录组分析结果基本一致。