【目的】通过NaCl胁迫下哈茨木霉(Trichoderma harzianum)ACCC32527差异转录组和代谢组数据分析,挖掘关键功能型差异表达基因(DEG)及次级代谢产物。【方法】采用RNA-seq和GC-TOF-MS技术分别完成0(T1)、0.4(T2)、0.6(T3)mol·L~(-1) NaCl胁迫下哈茨木霉ACCC32527的差异表达基因和次级代谢产物检测,利用相关软件及数据库(GO、COG、KEGG pathway等)完成对差异表达基因(|log_2fold change)|>1&FDR<0.01)和次级代谢产物(P-value≤0.05&VIP>1)的筛选、注释和分类,并进行RT-qPCR验证。【结果】NaCl胁迫处理后,ACCC32527的生长受到抑制,繁殖速率明显降低,并在该条件下分别获得T1 vs T2和T1 vs T3比较组637、1 570个DEG,获得DEG的16种表达模式,包括9种上调表达模式(950 gene)、3种下调表达模式(662 gene)和4种不规则表达模式(309 gene),共有的GO功能分类为催化活性、结合、细胞器、细胞膜部分、细胞膜、细胞部分、细胞、单组织过程和代谢过程,且占主要比例,约为61%—94%,其中处理前后基因比例差异较大的分类为催化活性。KEGG代谢通路富集分析发现,不同NaCl浓度处理下的DEG富集到的代谢通路种类及数量存在较大差异,分别有225和535个差异基因富集于20条KEGG通路,包括代谢通路、抗生素的合成和次级代谢产物合成,其中T2和T3处理下核糖体的生物合成途径分别有12个和18个相关基因受到抑制。从NaCl胁迫下差异转录组中筛选出活性氧清除、离子转运和细胞壁及胞外结构等共73个相关基因,并且多数为上调表达基因。另外,差异代谢组学分析表明,T3处理下,胞内小分子代谢物出现明显变化,累积量增加的代谢产物有30种,包括氨基酸及其衍生物、糖类及其衍生物和醇类,而减少的代谢物有53种,涉及脂肪酸和有机酸等。其中,甘油是已知的真菌重要渗透保护物,T3处理下ACCC32527胞内甘油水平提高,可参与细胞的渗透调节,维持渗透压平衡。RT-qPCR验证了7个差异表达显著基因,虽然在表达量上存在一定的差异,但表达趋势与RNA-seq分析结果基本一致,表明转录组测序结果的可靠性。【结论】NaCl胁迫引起哈茨木霉ACCC32527大量基因及次级代谢产物变化,获得了与NaCl胁迫调节相关基因和代谢产物,这些机制共同作用减少盐离子对细胞的毒害作用,研究结果可为研究木霉菌的耐盐机理提供重要信息。

采用菌物学、生物化学等方法,对从以色列引进的哈茨木霉T28进行盐碱胁迫研究,从菌株生长和营养代谢的角度研究了该菌株对盐碱胁迫的响应。结果表明:在NaCl、Na2SO4和NaHCO3胁迫条件下,菌株生长和营养代谢与盐碱浓度呈明显的相关性,在cNaCl≤0.4 mol/L和cNa2SO4≤0.2 mol/L胁迫的56 h内,T28均可长满平板(Φ=90mm),生物量为0.17~0.38 g/150 mL和0.51~0.52 g/150 mL,碳利用率为16%~43%和8.4%~8.7%,磷利用率为3.3%~4.7%和4.3%~5.7%。菌株生长和对碳、磷相对利用率受抑制程度均较低;在NaHCO3胁...

[目的]探究不同氮效率马铃薯品种氮代谢对低氮胁迫的响应及分子机制,为马铃薯氮高效育种及高产栽培提供理论依据。[方法]以氮高效马铃薯品种冀张薯12号和氮低效马铃薯品种尤佳70为材料,采用盆栽试验方式,设置低氮(3 mmol/L纯氮)和正常施氮(15 mmol/L纯氮)处理,在出苗后30 d测定植株的干质量、氮积累量及氮代谢相关酶(硝酸还原酶(NR)、亚硝酸还原酶(NiR)、谷氨酰胺合成酶(GS)、谷氨酸合成酶(GOGAT))活性,并依托Illumina HiSeq 2000平台进行转录组测序分析。利用DESeq进行差异基因分析,通过GO富集和KEGG Pathway数据库注释参与马铃薯氮胁迫响应的差异表达基因。[结果]2种施氮量下冀张薯12号茎叶和根系的干质量、氮积累量以及NR、NiR、GS、GOGAT活性均高于尤佳70。与正常施氮相比,低氮胁迫下2个马铃薯品种茎叶和根系的干质量、氮积累量及叶片的NR和NiR活性均显著降低,根系的GS活性显著提高;尤佳70叶片的GS活性、GOGAT活性和根系的GOGAT活性均显著降低,冀张薯12号根系的NiR活性显著降低。转录组分析结果表明,在2个施氮量对比组中,尤佳70叶片和根系中的差异表达基因分别有6 173个和249个,冀张薯12号叶片和根系中的差异表达基因分别有3 455个和1 990个。GO和KEGG结果显示,冀张薯12号的差异基因主要涉及氮代谢、氨基酸代谢等生物过程,尤佳70的差异基因主要涉及光合作用等生物过程。尤佳70和冀张薯12号根系中富集到氮代谢通路的差异基因分别有1个和8个,叶片中分别有11个和8个。尤佳70和冀张薯12号叶片的氨基酸相关代谢途径中涉及上调差异表达基因最多的均是苯丙氨酸代谢,其中编码苯丙氨酸解氨酶、肉桂酸-4-羟化酶、4-香豆酸-CoA连接酶的基因表达显著上调。[结论]在低氮条件下,氮高效品种冀张薯12号氮代谢相关酶活性较高,能减缓低氮胁迫对氮代谢的影响,进而保证较高的氮效率。

为了探明我国玉米主要推广品种郑单958、先玉335花粒期高温胁迫下基因表达及代谢物差异,挖掘玉米响应高温胁迫的关键基因及代谢物,利用Illumina HiSeq~(TM) 2500高通量测序技术分别对郑单958、先玉335高温胁迫7,14 d的穗位叶及对照叶进行高通量转录组测序,利用广泛靶向代谢组测序技术对样品进行高通量代谢物测序。转录组测序共获得214.81 Gb干净序列。郑单958高温胁迫7 d后检测到49个差异表达基因,其中24个为上调基因。继续高温胁迫14 d共检测到306个差异表达基因,其中130个为上调基因。高温胁迫7 d与高温胁迫14 d对比组中检测到1 462个差异表达基因,其中647个为上调基因。先玉335高温胁迫7 d后共检测到381个差异表达基因,164个上调基因。继续高温胁迫14 d后共检测到299个差异表达基因,226个为上调基因。高温胁迫7 d与高温胁迫14 d对比组中共检测到2 481个差异表达基因,其中1 275个上调基因。在郑单958、先玉335高温胁迫7 d后对比组中检测到6 646个差异表达基因,其中3 253个上调基因。在郑单958、先玉335高温胁迫14 d对比组中检测到5 958个差异表达基因,其中3 110个上调基因。代谢物测序共检测到654个代谢物,郑单958高温胁迫7 d后共检测出28个差异代谢物,共注释到5个代谢通路中。高温胁迫14 d后共检测出54个差异代谢物,共注释到13个代谢通路中,其中9个差异代谢物呈上调表达。先玉335高温胁迫7 d后检测出98个差异代谢物,共注释到24个代谢通路中,其中43个差异代谢物呈上调表达。先玉335高温胁迫14 d后共检测出38个差异代谢物,共注释到13个代谢通路中,其中14个上调表达。郑单958高温胁迫7 d与先玉335高温胁迫7 d对比组共检测出144个差异代谢物,共注释到36个代谢通路中,其中81个差异代谢物呈上调表达。郑单958高温胁迫14 d与先玉335高温胁迫14 d对比组共检测出158个差异代谢物,共注释到40个代谢通路中,其中81个差异代谢物呈上调表达。

通过“广靶代谢组+转录组”联合分析方法，对Na Cl胁迫下胡麻(Linure usitatissmum)耐盐种质资源R40和敏感资源R24根部差异表达基因与差异积累代谢物进行共表达分析，以探明代谢物与关键基因的协同变化情况，阐明胡麻的耐盐代谢过程和变化规律，发掘胡麻中可能调控抗盐代谢途径的重要基因。共鉴定出具有显著差异的代谢物741个，其中下调409个，上调332个，鉴定出差异表达基因30 233条，其中上调15 827条，下调14 406条，有注释结果的转录因子13 015个。富集前十的代谢物几乎都有相对应的功能转录因子富集，如黄酮和黄酮醇的生物合成、苯丙素的生物合成、氨基糖和核苷酸糖的代谢、花青素生物合成、类黄酮生物合成等通路均有对应的转录因子参与其中。胡麻根部在盐胁迫下调控抗盐代谢途径的重要基因主要以参与调控糖类代谢为主，增强细胞的渗透压，从而有效降低盐碱胁迫的伤害；筛选出的8个候选基因与5个代谢物具有明显的相关性，每一个代谢物与一个或多个基因相关联，色氨酸–生长素通路发现，GH3基因是一种典型的与植物生长素原初反应相关的基因，在植物的抗逆防卫反应中起重要作用。

【目的】克隆哈茨木霉几丁质合酶基因全序列,并对序列进行生物信息学分析。【方法】利用多聚酶链式反应(PCR)和染色体步移技术(genome walking)克隆几丁质合酶基因全序列;应用生物信息学软件对获得的基因序列及编码的蛋白序列进行分析;通过同源建模,对蛋白质的三维结构进行预测。【结果】从哈茨木霉(Trichoderma harzianum Th-33)中扩增得到几丁质合酶基因ThChsC(GenBank登录号HQ419000),编码区(CDS)长为2 835 bp,由4个外显子和3个内含子组成,开放阅读框(ORF)长2 688 bp,编码895个氨基酸。序列比对和同源性分析显示,该基因与串...

以流苏2年生苗木为材料,测定其叶片与茎在NaCl胁迫(0,0.4%,0.8%,1.2%)下的电阻抗图谱参数和电导率的变化,分析电阻抗图谱参数与电导率间的相关性,以讨论在NaCl处理下用电阻抗图谱法研究流苏受胁迫程度的有效性。结果表明:1)随着NaCl浓度的升高,叶片电导率显现逐渐升高或先升高后下降的趋势,茎的表现不明显;随着胁迫时间的延长叶片和茎电导率,整体均显现逐渐升高或先升高后降低的趋势。流苏茎的耐盐性强于叶片。2)随着NaCl浓度的升高,流苏叶片和茎的高频电阻率r、胞内电阻率ri和弛豫时间τ呈现先增大后减小或逐渐减小的特征;随着胁迫时间的延长叶片和茎的各处理的电阻抗图谱参数均在胁迫15天...

沙枣（ＥｌａＥａｇｎｕｓ ａｎｇｕｓｔｉｆｏｌｉａ ｌ．）是胡颓子科（ＥｌａＥａｇｎａｃＥａＥ）胡颓子属（ＥｌａＥａｇｎｕｓ ｌｎｎ．）的落叶小乔木，主要分布于我国西北五省，在新疆塔里木河、玛纳斯河、甘肃疏勒河、内蒙古的额济纳河两岸有天然的沙枣林分布，在北京、天津和山东等地也开展了引种栽培［１－２］。沙枣的果实、叶片、花和种子等在食品加工、药物提炼和饲用开发等领域具有较高的价值，沙枣木材也是高级家具和工艺品制作的良材［１］；同时，沙枣具有耐旱、耐盐碱和抗风沙的特性，并且沙枣根系上有固氮根瘤菌，能提高土壤肥力，改善土壤环境，因此，沙枣已经被列为半干旱地区和盐土地

bZIP转录因子广泛存在于植物中，在调节植物生长发育和非生物胁迫应答中起着重要作用。为探究bZIP转录因子在玉米干旱胁迫响应中的功能，在玉米苗期干旱胁迫处理5 d和复水3 d时，利用转录组测序技术对转录因子基因的表达变化进行分析，筛选到一个响应干旱和复水处理的bZIP转录因子(ZmbZIP26),共表达网络分析发现，ZmbZIP26处于网络调节的核心节点位置。ZmbZIP26基因含有558 bp的开放阅读框，编码185个氨基酸，属于亲水性蛋白。蛋白质进化树及保守序列分析发现，ZmbZIP26蛋白与高粱和芒草同源蛋白的同源性较高，而且在同一氨基酸位置上的保守基序相同。启动子顺式作用元件分析表明，ATG上游2 000 bp区域内含有干旱响应元件、激素响应元件和光响应元件等。qRT-PCR分析发现，ZmbZIP26是组成型表达基因，在幼茎、雌穗和根中高表达，且ZmbZIP26基因积极响应干旱、高温、高盐、氮胁迫及恢复正常的过程，可能在植物的抗逆过程中发挥着重要作用。亚细胞定位分析发现，ZmbZIP26属于核蛋白，定位在细胞核上。蛋白质互作预测发现，ZmbZIP26可能与锌指蛋白、丝氨酸蛋白、钙依赖性蛋白、谷胱甘肽转移蛋白等互作构建了一张调控网络，协同调控玉米生长发育和胁迫应答过程。

以西瓜品种和欣(HX)和黑媚娘(HMN)的幼苗叶片为材料,经60 mg/kg氯化镉溶液处理(以未经氯化镉溶液处理的为对照)后,采用Illumina测序技术对转录组进行高通量测序,探讨西瓜镉胁迫的分子机理。结果表明:在黑媚娘处理组(HMN)和对照组(HMNCK)中共检测到8 579条差异基因,其中有4 297条基因表达显著上调,4 282条基因表达明显下调;在和欣处理组(HX)和对照组(HXCK)中共检测到7 108条差异基因,其中有3 662条基因表达显著上调,3 446条基因表达显著下调;不同组间差异基因的维恩图表明,HMN和HX组间有5 656条差异基因;根据不同数据库的注释信息,发现西瓜体内与镉胁迫相关的差异表达基因主要包括响应光合作用、信号转导、次生物质代谢、转录因子等相关基因;转录因子MYB、AP2/EREBP、bHLH、NAC和WRKY家族可能在西瓜响应镉胁迫中具有重要作用。用q RT–PCR技术对随机挑选的7条基因进行荧光定量验证,结果与Illumina测序数据一致,证实了差异表达基因数据的有效性。

[目的]本文旨在探明2个大豆品种根系对早期缺铁胁迫的响应差异。[方法]以铁高效品种‘吉育99’和铁低效品种‘吉育93’大豆为材料进行水培试验，将对照(25μmol·L~(-1) Fe)和早期(1和6 h)缺铁胁迫(1μmol·L~(-1) Fe)处理的大豆根系进行转录组测序。[结果]与对照相比，铁高效大豆品种差异基因数随着缺铁胁迫时间的延长而增加，而铁低效品种差异基因数随着胁迫时间的延长而降低。差异表达基因(differential expression gene, DEG)的KEGG通路富集结果表明，缺铁胁迫处理1 h, 2个大豆品种在苯丙烷生物合成、倍半萜和三萜生物合成、植物病原体相互作用、植物MAPK信号通路、淀粉和蔗糖代谢等代谢通路中的相关基因表达有显著差异。缺铁胁迫处理6 h, 2个大豆品种在苯丙烷生物合成、倍半萜和三萜生物合成、类黄酮生物合成、谷胱甘肽代谢等代谢通路中的相关基因表达有显著差异。差异表达基因转录因子分析结果表明，铁高效大豆品种受早期缺铁胁迫调控的转录因子家族主要有AP2/ERF-ERF、C2H2、MYB、WRKY、bHLH、NAC,铁低效大豆品种中转录因子家族主要有AP2/ERF-ERF、MYB、WRKY、bHLH、NAC。[结论]缺铁胁迫1 h大豆就能对缺铁胁迫作出响应，且不同铁效率大豆品种在苯丙烷生物合成、转录因子调控等相关基因的表达存在显著差异。

低温冷害是限制我国东北地区花生产量和品质提升的主要环境因素,鉴定花生耐寒的关键功能基因,并揭示其调控机制是创制耐寒花生种质,提高花生耐冷性的重要途径之一。以耐寒型花生品种农花5号为试验材料,对6℃低温胁迫下的花生叶片进行RNA-Seq转录组测序分析。结果表明:与对照相比,共有3910个基因持续差异表达;GO功能注释将差异表达基因富集到43个功能组,其中生物学过程中的代谢过程富集到的基因数目最多;KEGG代谢通路分析将持续差异表达基因富集到116个代谢通路中,其中植物激素信号转导、淀粉和蔗糖代谢、植物-病原菌互作、植物昼夜节律和α-亚麻酸代谢途径富集到的基因数目最多,且大部分为上调表达,在低温胁迫中具有重要作用。从转录组数据库中随机挑选6个差异表达基因进行qRT-PCR验证,其表达趋势与高通量测序结果相一致,证明了转录组数据的可靠性。研究结果将为揭示花生响应低温胁迫的分子机制提供理论依据。

分析干旱胁迫下枇杷叶片的转录组,挖掘功能基因并对其差异表达基因进行筛选和分析,为枇杷抗旱提供理论依据。利用新一代高通量测序技术测序,对测序结果进行de novo拼接、功能注释和ORF预测,将差异表达基因在COG、GO和KEGG数据库中进行比对注释。测序结果表明,获得转录本共88 530个,平均长度为740.64 bp,ORF41 748条。COG、GO和KEGG数据库将转录本分别划分为24,54个功能类别及291条代谢通路中。25 197个差异表达的基因在30条代谢通路中显著富集,与酶活性、激素合成代谢和

玉米籽粒发育时期对高温胁迫十分敏感，严重影响玉米的产量和品质。为研究不同耐高温玉米自交系籽粒响应高温胁迫的基因表达差异，在基因水平解析不同玉米种质响应高温胁迫的分子机制，以前期筛选的耐高温自交系XN202和敏感自交系CT110为材料，利用RNA-Seq技术挖掘正常和高温胁迫下授粉后15 d玉米籽粒的差异表达基因(DEG)。与对照相比，XN202和CT110分别检测到1 517,1 012个DEGs,共同的DEGs有142个，包含7个转录因子。不同耐高温玉米自交系的籽粒对高温胁迫的响应存在明显差异，通过基因本体(GO)功能注释和全基因组及代谢途径数据库(KEGG)信号通路富集分析筛选出DNA复制、核小体、微小染色体维持复合物、DNA引物酶-多聚酶α复合体、营养库活性、丙氨酸、天冬氨酸和谷氨酸代谢等共同参与响应高温胁迫的应答过程。通过生物信息学分析，共有374个DEGs位于已报道的玉米耐高温相关性状QTL区间，其中42个DEGs为已报道的耐高温相关基因。综上所述，高温胁迫下玉米籽粒会形成复杂的细胞保护和防御系统，AP2/ERF、MYB、bHLH、NAC、HSF、HSP等耐高温相关的DEGs可能在分子调控网络中发挥重要作用。

【目的】分析镉胁迫下中辽1号杨的转录水平,筛选差异表达基因,挖掘与镉胁迫相关的功能基因,为深入探索中辽1号杨对镉胁迫响应的分子机制提供理论依据。【方法】以中辽1号杨为试验材料,采用盆栽试验的方法将其扦插于土壤Cd含量为20 mg/kg(M20)的花盆中,以不加Cd为对照(CK),80 d后采集不同处理中辽1号杨叶片进行转录组测序,使用DESeq2软件筛选CK与M20处理中辽1号杨的差异表达基因,将得到的差异基因在基因本体数据库(gene ontology,GO)、京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)、蛋白相邻类聚簇数据库(clusters of orthologous groups of proteins,COG)中进行注释,分析差异表达基因在不同数据库中的注释信息。随机挑选6个差异表达基因(超氧化物歧化酶(SOD)、丝氨酸/苏氨酸蛋白激酶(STK)、脱落酸(ABA)、T复合蛋白1(TCP1)、MYB基因(MYB)、丝裂原活蛋白激酶12(MAPK12)),利用Primer 5软件设计特异性引物进行实时荧光定量PCR(RT-q PCR)试验,检测这6个基因表达水平的变化,验证转录组结果的准确性。【结果】在CK和M20处理的中辽1号杨中,共发现3 812个差异表达基因,其中上调表达的有2 209个,下调表达的有1 603个。GO分析发现,差异表达基因注释到3大类49个功能组中,其中与代谢过程有关基因最多。转录组KEGG代谢通路分析发现,差异表达基因主要在信号转导途径、代谢途径和生物合成途径富集。COG分析发现,1 455个差异基因被注释到22种分类中,其中一般功能预测基因注释最多,其次是信号转导机制基因。结合各数据库分析结果发现,Cd胁迫下中辽1号杨转录组的注释结果中代谢过程和信号转导过程相关基因较多。镉胁迫下差异表达基因最多的家族是ABC、MYB、WRKY、b HLH和NAC基因家族,挖掘出与镉胁迫相关基因WRKY家族基因47(WRKY47)、硝酸盐转运蛋白基因(NRT)、ABC家族转运蛋白基因2(ABC2)、苹果酸脱氢酶(MDH)、谷胱甘肽S-转移酶基因(GSTs)、NAC家族基因2(NAC2)的表达量显著上调,MYB家族基因44(MYB44)、重金属相关异戊二烯化植物蛋白基因39(HIPP39)的表达量显著下调。RT-q PCR结果显示,随机挑选6个差异表达基因表达量的变化趋势与转录组测序结果一致。【结论】镉胁迫后中辽1号杨转录组的差异表达基因主要富集在代谢过程和信号转导过程,挖掘出与镉胁迫相关基因,分别是WRKY47、NRT、ABC2、MDH、GSTs、NAC2、MYB44、HIPP39,为深入探索杨树对镉胁迫响应的分子机制提供了理论依据。