以西瓜品种和欣(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测序数据一致,证实了差异表达基因数据的有效性。

本文采用高通量Illumina Hiseq测序平台对紫花苜蓿(Medicago sativa)的铅胁迫处理组(Pb96)和对照组(CK)进行转录组测序,并对测序数据进行分析,探究紫花苜蓿抗铅应答的分子机制。结果表明,在铅胁迫处理下共检测到2 242个差异表达基因(Differentially expressed genes,DEGs)。其中,1 321个DEGs上调表达,921个DEGs下调表达。GO (Gene ontology,http://www.geneontology.org/)富集分析,差异表达基因主要涉及物质代谢、蛋白结合、催化活性等。KEGG (Kyoto encyclopedia of genes and genomes,http://www.genome.jp/kegg/)分析,差异基因主要涉及代谢途径、光合作用、淀粉和蔗糖代谢等。COG (Clusters of orthologous groups of proteins/orthologous groups of genes)分析,差异基因主要涉及"碳水化合物的运输和新陈代谢"、"翻译、核糖体结构和生物发生"、"次生代谢产物的生物合成、转运和分解代谢"等。本研究通过对紫花苜蓿根系转录组分析,为研究紫花苜蓿耐铅分子机制提供参考。

【目的】分析镉胁迫下中辽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,为深入探索杨树对镉胁迫响应的分子机制提供了理论依据。

【目的】探究盐碱胁迫对玉米中木质素、叶绿素含量的变化，以及其生物合成途径上关键酶基因表达水平的变化，为揭示玉米响应盐碱胁迫的分子机制、培育抗盐碱玉米新品种提供理论基础。【方法】对玉米苗期进行盐碱胁迫处理，并进行高通量转录组测序、叶绿素和木质素含量测定。【结果】转录组测序共获得39 722个基因，其中差异表达基因12 432个。KEGG注释分析表明，差异基因主要富集在碳代谢、苯丙烷生物合成通路、淀粉和蔗糖代谢通路。生化分析表明，盐碱胁迫下玉米叶绿素含量和总木质素含量分别降低30.45%和31.26%，其中H和S型木质素单体含量降低，而G型木质素单体含量升高。对基因表达水平和成分含量之间的关系分析表明，叶绿素合成酶基因的下调和叶绿素b还原酶基因的上调导致了叶绿素含量的降低，苯丙素生物合成途径中的161个差异表达基因导致木质素和木质素单体含量的变化。【结论】盐碱胁迫抑制玉米木质素的生成，且主要抑制H和S型木质素单体的生成；盐碱胁迫通过调控叶绿素合成酶基因的下调和叶绿素b还原酶基因的上调，致使叶绿素的生成量减少。本研究为采用基因工程进行分子育种提高玉米耐盐碱性和调控玉米木质素含量提供了理论依据。

【目的】皂苷是苦瓜果实苦味的主要来源，具有降血糖、抗癌等多种药用价值。鉴定和挖掘参与调控苦瓜皂苷生物合成的代谢通路及基因，为进一步深入解析苦瓜皂苷形成的分子机制提供参考。【方法】以苦瓜高代自交系GK24为试验材料，分别在子房期（T1）、幼果期（T2）、商品果期（T3）和完全成熟期（T4）取果实组织，通过香草醛-冰醋酸法测定不同时期的苦瓜皂苷含量，利用转录组测序的方法鉴定差异表达基因。【结果】从T1到T4，苦瓜内源皂苷含量随着果实发育而呈现显著下降的趋势。转录组测序共鉴定到17 504个基因，在T1-vs-T2、T2-vs-T3和T3-vs-T4三组比较中，下调表达基因数量均高于上调表达基因数量。GO富集分析表明含磷复合物代谢过程、驱动蛋白复合体和2-琥珀酰-6-羟基-环己二烯-1-羧酸合成酶活性分别是生物过程、细胞组分及分子功能模块最显著富集的条目。KEGG分析显示全局与概述图谱、碳水化合物代谢和氨基酸代谢是3组比较中共同的最显著富集的代谢通路。根据表达模式可将基因划分为20个模块，其中profile 0中基因的表达量从T1到T4逐渐降低，与苦瓜果实发育过程中皂苷含量变化规律一致。从profile 0中鉴定出与苦瓜皂苷生物合成相关的基因14个，包括萜类骨架生物合成通路上的基因AAT1、HMG1、MVK、PMK、MVD2和FPS1，倍半萜与三萜生物合成途径中的基因SS12、SQE1和CPQ及后修饰酶基因CYP97A3、CYP71AN24、UGT94E5及UGT73C6。实时荧光定量PCR(qRT-PCR)结果与RNA-seq结果一致。【结论】苦瓜果实中皂苷含量随着果实不断成熟而逐渐降低。苦瓜皂苷生物合成首先通过上游的萜类骨架生物合成（ko00900）代谢途径完成萜类皂苷骨架的积累，随后通过倍半萜与三萜生物合成（ko00909）代谢途径以及下游的氧化还原酶、糖基转移酶的修饰最终生成皂苷，本研究通过转录组测序共鉴定到14个与苦瓜皂苷生物合成相关的基因。

镉(Cd)作为生物非必需、毒性极强的蓄积性重金属,易通过食物链进入人体,严重威胁人类健康。扇贝相对于其他贝类具有特异性蓄积镉的能力,成为水产品安全领域关注的焦点。为了阐释扇贝高蓄积镉的分子机制,本研究以镉胁迫的栉孔扇贝(Chlamys farreri)消化盲囊组织为研究对象,通过转录组测序技术进行基因转录水平分析。通过比较转录组拼接获得105071个unigene和3800个差异基因,对所得unigene进行功能注释,发现这些基因集中在蛋白绑定、细胞黏附、免疫应答、细胞凋亡和能量代谢等方面。对这些蛋白的分

【目的】通过转录组测序和分析,发掘美洲黑杨无性系‘中菏1号’耐镉关键功能基因,揭示美洲黑杨在镉胁迫下的响应机制。【方法】以‘中菏1号’为试验材料,通过沙培试验法对处在不同质量浓度(0 mg/L(ZH1C0),10 mg/L(ZH1C1),20 mg/L(ZH1C2))Cd~(2+)环境下的无性系进行胁迫处理,而后采集叶片,提取RNA,并使用Illumina HiSeq TM2500系统对转录组进行高通量测序。【结果】主成分分析确定PC1作为区分不同Cd~(2+)浓度样本的标准;经过差异表达分析,在ZH1C0和ZH1C2之间共检测到4 109条差异基因,其中2 741条基因显著上调表达,1 368条基因显著下调表达;在ZH1C1和ZH1C2之间检测到3 710条差异基因,其中1 969条基因显著上调表达,1 741条基因显著下调表达;不同数据库之间的注释表明Cd~(2+)不仅会影响叶绿体、线粒体等细胞器的微结构,而且会干扰光合作用和呼吸作用中的电子传递。转录因子分析发现,WRKY、MYB、ZIP、ERF、bHLH在镉胁迫响应和耐受中具有重要作用。【结论】在‘中菏1号’中,异戊二烯化植物蛋白在阻止Cd~(2+)进入细胞内部发挥了重要作用;谷胱甘肽不仅可以促进细胞内Cd~(2+)的螯合,还可以缓解细胞内的氧化胁迫;ABC家族蛋白和MATE家族蛋白是‘中菏1号’体内转移Cd~(2+)螯合蛋白的关键载体。

为探讨盐度变化对大黄鱼(Larimichthys crocea)肝脏基因转录表达的影响，将体质量为(59.68±1.32) g的大黄鱼分别在盐度25（对照组）、盐度36（高盐组）、盐度10（低盐组）的水体中暴露24 h，采用Illumina Hiseq 2500对肝脏样本进行转录组测序。结果显示，从高盐vs.对照组(H vs. C)、低盐vs.对照组(L vs. C)和高盐vs.低盐(H vs. L)中分别筛选出71个、180个和438个差异基因(DEGs)。GO (Gene Ontology)和KEGG (Kyoto Encyclopedia of Genes and Genomes)分析发现，H vs. C中的DEGs显著富集在三羧酸循环、MAPK信号通路、内质网蛋白加工和Toll样受体信号通路等；L vs. C中的DEGs显著富集在FoxO信号通路、戊糖磷酸途径、mTOR信号通路、RNA降解、P53信号通路和内质网蛋白加工等，表明大黄鱼对高盐和低盐胁迫的响应均涉及到糖类代谢、内质网应激、细胞凋亡和自噬等生理功能。H vs. L中的DEGs显著富集在内质网蛋白加工、FoxO信号通路、甘油脂代谢、糖酵解/糖异生、吞噬体和P53信号通路等，表明大黄鱼在高盐和低盐适应过程中，内质网应激、脂类代谢、糖类代谢、细胞凋亡和自噬等生理功能存在差异。研究结果可为深入研究鱼类对盐度胁迫的响应机制奠定基础。

为了探究低氧胁迫下中华圆田螺(Cipangopaludina cathayensis)肝脏组织基因的差异表达，本研究通过高通量测序技术，分析中华园田螺低氧胁迫组(2.5 mg/L)和常氧组(6.9 mg/L)某些基因的差异表达，并对差异基因进行生物信息学分析，进一步采用实时荧光定量PCR (RT-qPCR)对关键差异表达基因进行验证。结果显示，测序共获得232 379条基因(unigene)，与对照组比较，低氧胁迫组筛到176个差异基因，包含64个上调基因和112个下调基因。GO功能注释分析显示，差异基因主要富集在生物学过程中的几丁质代谢过程和含氨基葡萄糖的复合代谢过程，细胞组分中的胶原三聚体成分，分子功能中的几丁质结合功能和糖衍生物结合功能。KEGG通路富集分析显示，差异基因主要集中于环境信息处理、遗传信息处理、代谢和生物系统这4大类通路。6个关键差异基因的RT-qPCR结果显示，热休克蛋白70B2和热休克蛋白β-6基因表达量上调，几丁质酶蛋白4、α-1胶原蛋白(XIV)、α-4胶原蛋白(XIV)、5-磷酸酶蛋白基因表达量下调，证实了转录组测序结果的可靠性。本研究发现，低氧胁迫激活了中华圆田螺适应缺氧的生理活动，并获得了低氧胁迫下中华圆田螺肝脏组织中相关功能基因的表达信息，为深入研究中华圆田螺响应低氧胁迫的调控机制提供了基础数据和理论依据。

为探讨盐度变化对大黄鱼(Larimichthys crocea)肝脏基因转录表达的影响，将体质量为(59.68±1.32) g的大黄鱼从盐度25转入盐度36或10的水体中暴露24 h，采用Illumina Hiseq 2500对肝脏样本进行转录组测序。结果显示，从高盐vs.对照组(H vs. C)、低盐vs.对照组(L vs. C)和高盐vs.低盐(H vs. L)中分别筛选出71个、180个和438个差异基因(DEGs)。GO (Gene Ontology)和KEGG (Kyoto Encyclopedia of Genes and Genomes)分析发现，H vs. C中的DEGs显著富集在三羧酸循环、MAPK信号通路、内质网蛋白加工和Toll样受体信号通路等；L vs. C中的DEGs显著富集在FoxO信号通路、戊糖磷酸途径、mTOR信号通路、RNA 降解、P53信号通路和内质网蛋白加工等，表明大黄鱼对高盐和低盐胁迫的响应均涉及到糖类代谢、内质网应激、细胞凋亡和自噬等生理功能。H vs. L中的DEGs显著富集在内质网蛋白加工、FoxO信号通路、甘油脂代谢、糖酵解/糖异生、吞噬体和P53信号通路等，表明大黄鱼在高盐和低盐适应过程中，内质网应激、脂类代谢、糖类代谢、细胞凋亡和自噬等生理功能存在差异。本结果可为深入研究鱼类对盐度胁迫的响应机制奠定基础。

利用农杆菌介导的花序浸润法,将金柑(Fortunella margarita(Lour.)Swingle)MLP2–1基因转入野生型拟南芥(Col–1),对转MLP2–1基因拟南芥在冷冻胁迫下目的基因的转录组数据进行分析。结果显示:转MLP2–1基因拟南芥的抗冻能力强于野生型拟南芥;转录组测序分析表明,转MLP2–1拟南芥与野生型拟南芥植株差异表达基因为1 422个,其中551个表达上调,871个表达下调;差异表达基因数目在GO数据库注释到1 208个,在KEGG数据库中注释到474个,在COG数据库中注释到603个;KEGG代谢通路富集结果表明,差异表达基因主要集中在植物激素与信号转导代谢通路。

为探究花丝对干旱胁迫的分子响应，以安农591和先玉335为材料，设置土壤含水量为80%田间持水量(CK)和60%田间持水量(DS)的2个处理，在玉米花丝伸长速增期取样，进行转录组测序。通过对比2个品种干旱组和对照组的基因表达差异，明确玉米花丝响应干旱的关键通路和相关候选基因。结果表明：苯丙烷类生物合成途径(ko00940)和淀粉和蔗糖代谢途径(ko00500)为玉米花丝响应干旱的关键通路。在苯丙烷类生物合成途径通路中编码4-香豆酸-辅酶A连接酶(4CL)、肉桂酰辅酶A还原酶(CCR)、肉桂醇脱氢酶(CAD)和过氧化物酶(POD)的基因在干旱的玉米花丝中低表达。相关分析表明，编码4CL、CCR、CAD和POD的基因下调与木质素含量降低有关。在淀粉和蔗糖代谢途径中转化酶和蔗糖合酶基因在干旱的玉米花丝中低表达；相关分析表明，转化酶和蔗糖合酶基因下调与蔗糖代谢有关。综上所述，本研究确定了参与蔗糖和木质素代谢的关键基因和代谢物，木质素合成可能与玉米花丝耐旱相关，4CL、CCR、CAD和POD是木质素合成通路中玉米花丝耐旱的候选基因；此外，淀粉和蔗糖代谢途径可能与花丝伸长和花丝耐旱有关，转化酶和蔗糖合酶可能是调控花丝耐旱的关键酶。

[目的]通过对红豆杉科穗花杉属的穗花杉进行转录组测序,为穗花杉的萜类合成途径和分类学研究提供支持。[方法]利用Hi Seq2500技术对穗花杉的茎叶进行转录组测序。[结果]穗花杉转录组测序共得到8.14 Gb的clean data。从头组装共获得82 884条unigene,其中有27 495条unigene被注释。此外,在82 884条unigene中共搜索到2 827个SSR位点,单核苷酸重复SSR出现频率最高(60.1%),其次为三核苷酸重复SSR(25.4%)。在穗花杉unigene中挖掘出了1个

【目的】WRKY转录因子是一类植物响应生物、非生物胁迫,对生长发育都起重要调控作用的转录因子。在甜菜全基因组信息分析的基础上,鉴定WRKY家族基因(Bv WRKYs),解析其组织特异性及盐、热胁迫下的表达情况,为该类基因的功能研究提供参考,为观赏甜菜和石竹目其他观赏植物的基因工程打下基础。【方法】以75条拟南芥WRKY蛋白为参考,根据WRKY保守蛋白序列(PF03106)利用hmm和BLAST同源性搜索对甜菜WRKY家族基因进行鉴定。利用Map Inspect、GSDS2.0、MEGA5.0、DNAMAN

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