烟草幼苗响应弱光胁迫的转录组分析

为了研究烟草幼苗对弱光胁迫反应的分子机制,以'云烟87'为材料,构建全光照和弱光处理条件下大十字期烟苗cDNA文库,利用Illumina测序技术进行转录组测序,筛选差异表达基因.结果表明:借助RNA-seq技术共筛选到2 956个DEGs,其中弱光相对于全光照表达上调的DEGs有691个,表达下调的DEGs为2 265个;通过GO功能分类和KEGG富集分析,将差异基因划分为30个功能类别和113个代谢途径,包括代谢进程、碳水化合物结合、核酸结合转录因子活性、水解酶活性、转录因子活性、催化活性、内质网蛋白质加工、植物-病原体互作、苯丙烷生物合成、苯丙氨酸代谢、胁迫反应、激素信号转导以及光合作用等.本研究将为揭示烟草对弱光胁迫反应的分子机制提供理论依据.     

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

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾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途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

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

[目的]本研究旨在通过转录组测序技术分析低温胁迫引起的阿尔泰蝠蛾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途径.其中代谢过程、催化、结合活性类群占有重要比例,核糖体、碳代谢、剪接体等途径显著富集.另外,热激蛋白、昆虫表皮蛋白、海藻糖酶、超氧化物歧化酶等非生物胁迫相关基因显著上调表达.[结论]阿尔泰蝠蛾幼虫低温胁迫转录组分析揭示,代谢过程、细胞过程、生物调节、对刺激的反应等生物学过程相关基因和部分非生物胁迫响应基因显著上调表达,提示蝠蛾幼虫可能从抗氧化防御、分子伴侣、体温调节和维持细胞的渗透平衡等多方面应对低温胁迫.     

东乡野生稻苗期响应低温胁迫的转录组分析

为了解东乡野生稻（Oryza rufipogon）对低温胁迫的响应机制,对苗期的RNA-seq转录表达谱进行了研究。结果表明,与对照相比,共检测到10 200个差异表达基因（DEGs）,其中5 201个上调表达,4 999个下调表达,其中有426个DEGs位于已报道的水稻耐冷QTL区间,且37个为耐冷调控相关的家族基因。GO功能分类和KEGG代谢路径分析表明,核酸结合转录因子活性、氨基酸生物合成以及光合作用代谢等均参与响应低温胁迫过程。实时荧光定量分析表明,ABA响应蛋白基因、MYB转录因子和40S核糖体蛋白SA基因等12个可能与低温胁迫响应相关的DEGs表达模式与RNA-seq的一致。可见,植物激素传导途径和转录因子相关调控基因在东乡野生稻苗期响应低温胁迫过程中起重要作用。     

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

利用新一代高通量测序手段——转录组测序（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对几个转录因子进行了盐胁迫下的表达检测,其趋势与转录组分析结果基本一致。     

准噶尔小胸鳖甲短时低温胁迫响应的转录组分析

【目的】拟步甲科昆虫小胸鳖甲Microdera punctipennis是分布于中国新疆古尔班通古特沙漠的特有物种,具有很强的耐寒性,但其低温响应的分子机制尚不明确。本研究旨在利用转录组测序技术丰富小胸鳖甲的已知基因信息,分析在短时低温胁迫下,上调表达基因的主要类群及参与的主要代谢通路。【方法】利用Trinity软件对分别在4℃低温和25℃(对照)下处理3 h的小胸鳖甲成虫RNA-Seq数据进行从头组装,利用Blast2go软件进行基因注释。通过DESeq和GFOLD软件分析差异表达基因,并对上调表达基因进行GO(Gene Ontology)和KEGG(Kyoto Encyclopedia of Genes and Genomes)代谢途径富集分析。【结果】测序过滤后得到68 165 001个有效读序,51 712个平均长度为984 bp的非冗余基因,其中29 925个基因(约57.87%)有同源序列(E-value10-5),与拟步甲科的赤拟谷盗Tribolium castaneum相应基因相似性最高(核苷酸序列一致性为72.75%)。低温响应上调表达基因有514个,富集到35个GO类群,18个KEGG途径。其中胁迫响应类群、生物调节类群和免疫系统过程类群都占有重要比例,嘌呤代谢、硫胺素代谢、糖酵解/糖异生等代谢通路显著富集。在胁迫响应类群中,热激蛋白Hsp90基因、超氧化物歧化酶(SOD)基因和钙联接蛋白Calnexin基因等8个非生物胁迫相关基因显著上调表达。【结论】荒漠昆虫小胸鳖甲低温转录组揭示,在4℃冷驯化过程中,代谢过程、胁迫响应、生物调节和免疫系统等生物学过程相关基因表达显著上调,其中几个非生物胁迫响应基因提示小胸鳖甲可能从分子伴侣、细胞周期阻滞、线粒体稳定性、抗氧化胁迫等多方面应对低温胁迫。KEGG富集分析所揭示的小胸鳖甲在低温下的转录组代谢通路与其他物种有较大差异。研究结果为后续生物信息学分析及小胸鳖甲耐寒关键基因的发掘及耐寒机制的揭示提供了基础数据。     

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%。综合分析表明,黑果枸杞对盐胁迫的反应是一个多基因参与、多个生物过程协同调控的过程,基因表达量的变化可能是基因调控的主要方式。     

萱草叶片响应低温胁迫的转录组分析

低温胁迫是萱草（Hemerocallis fulva）生长过程中经常会遭遇的一种非生物胁迫。比较了萱草叶片在低温处理（10、5、0 ℃）下转录组与对照（15 ℃）数据的差异,共筛选出差异表达基因2 457个,其中上调基因1 253个,下调基因1 204个。差异表达基因主要富集在细胞过程、代谢过程和催化活性等49个GO过程,代谢途径、次生代谢产物的生物合成、植物激素信号转导等42条KEGG代谢途径中。其中参与植物激素信号转导通路的差异表达基因发生了不同程度的变化,GH3.10基因上调至对照组的13.624倍,IAA1基因下调0.120倍;参与可溶性糖合成通路的差异基因发生了0.076~28.114倍不同程度的变化。随后对3个低温处理组共有的29个差异表达基因进行热图和网络调控分析,基于基因在网络调控中的位置,对ABCF5、OFPs和SWEETs等基因在冷应答的作用进行了分析。本研究结果为进一步挖掘萱草低温响应的关键基因及耐寒萱草种质开发、分子育种提供了一定的理论支撑。     

基于转录组测序的扎龙湿地野大麦耐盐性分析

通过转录组测序技术对扎龙湿地野大麦幼苗在200 mmol/L NaCl胁迫2 d后的叶片进行分析。对测序结果进行De novo拼接后,将差异表达基因在GO、KEGG数据库中进行比对注释。结果表明:NaCl胁迫2 d后,对照组与处理组分别获得Unigene序列61038个和46754个,检测到差异表达基因25465个,差异基因GO功能注释到3个大类的55个功能组。差异表达基因被注释到135个Pathway上,直观地显示出NaCl胁迫下野大麦幼苗体内发生调节及改变的代谢过程和信号通路,其中主要涉及光合作用、脯氨酸代谢、叶绿素代谢等途径。通过对野大麦幼苗叶片转录组分析,发掘相关耐盐基因,有助于培育野大麦及其近缘作物新品种,并对揭示植物耐盐性分子机制及相关代谢途径具有重要意义。     

独行菜种子bHLH类转录因子基因家族及幼苗laICE1表达对冷胁迫的响应

该研究在转录组数据基础上,以独行菜(Lepidium apetalum Willd.)种子耐受和不能耐受低温萌发的2组样本为材料,对其bHLH转录因子家族成员进行分析,并对该家族成员表达差异极显著的laICE1基因进行克隆、序列分析、结构预测,以探讨独行菜幼苗中laICE1基因表达对低温胁迫的响应特征。结果表明:(1)萌发中的独行菜种子,至少有83个bHLH转录因子序列表达,相对于低温萌发停滞组,在低温萌发耐受组的独行菜种子中,有10个下调、13个上调、60个表达差异不显著。其中表达量极显著下调的序列c20009_g1具有1 503bp开放阅读框,GO注释到ICE1基因,该基因命名为laICE1。(2)laICE1基因编码500aa,蛋白分子量为54 635.19kD,理论等电点为5.45,分子式为C2364H3742N688O758S22;该蛋白具有保守结构域bHLH。(3)定量分析表明,laICE1基因在非低温胁迫的独行菜种子中的表达量显著低于一直处于低温胁迫中不能进行萌发的独行菜种子,这与转录组数据库中该序列表达情况一致;而laICE1基因在独行菜幼苗期经低温处理后,其表达量显著上调,表明laICE1基因可能在独行菜幼苗耐受低温生长中具有一定的作用。     

杨树MYB转录因子家族基因应答盐胁迫表达特性分析

MYB转录因子家族是植物中数量最多的转录因子家族之一,在植物次生代谢调节、信号转导和抗逆等生物过程起重要作用。根据MYB转录因子结构域组成差异可分4个亚家族：即1R-MYB（MYB-relaed）、R2R3-MYB、3R-MYB和4R-MYB。其中,R2R3-MYB亚家族数量最多,可进一步分为22个亚组;利用生物信息学分析杨树MYB转录因子蛋白序列的保守结构域、系统发生、基因组定位、氨基酸组成和理化性质等;参照拟南芥MYB转录因子功能,预测杨树MYB转录因子功能;基于84K杨转录组测序和RT-qPCR分析,从301个杨树MYB转录因子基因中筛选出69个应答盐胁迫基因（P≤0.05）。其中,上调表达基因32个,下调表达基因37个。该研究可为进一步研究杨树MYB家族基因功能提供参考依据。     

基于转录组的八角挥发油合成相关基因挖掘与分析

为更好地挖掘八角(Illicium verum)挥发油合成相关基因，该文对挥发油性状差异显著的优良无性系桂角69号及普通品种砧01号叶片进行了转录组测序及组装注释，并对差异表达基因进行了GO分类和KEGG通路分析。结果表明：(1)转录本经组装后获得84 182条序列，使用NR、NT、Swiss-Prot、KEGG、KOG、GO和Pfam数据库进行序列比对，共注释了59 161条序列，筛选出30 572个差异表达基因。与砧01号相比，桂角69号叶片中上调基因有15 025个，下调基因有15 547个。(2)GO分类结果显示共有20 287个差异基因被注释。KEGG分析结果表明，有21 600个差异基因被注释到133条KEGG通路上，其中挥发油合成相关的单萜生物合成通路、萜类骨架生物合成通路、苯丙素合成通路中的芳樟醇合酶、月桂烯合酶、香叶基香叶基焦磷酸合酶、肉桂酰辅酶A还原酶、咖啡酸3-O-甲基转移酶、肉桂醇脱氢酶等关键酶基因呈差异表达。(3)转录因子分析发现差异表达基因分布于31个转录因子家族，其中MYB家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

脂多糖调控对大鼠心脏微血管内皮细胞的转录组分析

目的探索脂多糖(LPS)对大鼠心脏微血管内皮细胞(rCMECs)转录组的调节作用。方法对照组(正常培养rCMECs),LPS组(100 ng/mL LPS处理6 h的rCMECs),每组进行3个生物学重复转录组测序。得到差异基因后,使用实时定量PCR对部分差异基因mRNA的表达进行验证。分别对上调和下调基因进行GO和KEGG富集,并对差异基因进行共表达网络分析。采用独立t检验进行统计学分析。结果LPS处理后,265个基因表达上调,118个基因的表达下调。前10个最显著上调基因为:Mt2a、Cyp7b1、Sod2、Icam1、Ccl2、AC128848.1、Mt1、Cebpd、Serpinb2和Tnfrsf11b。前10个最显著下调基因为:Cavin2、Ankrd1、Edn1、Prss35、Lmod1、Dhrs3、Ttc22、Sema6a、Map2k3和Sema7a。定量PCR的结果表明Mt2a、Sod、Ccl2、Cxcl1、Icam1和Vcaml基因的表达得到了上调(P均<0.01);而Cavin2、Ankrd1、Edn1和Prss35基因表达下调(P均<0.05)。GO和KEGG富集的结果表明,上调基因与内皮细胞对炎性免疫细胞的趋化作用和黏附作用密切相关;而下调基因则是与钙离子信号和G蛋白相关通路以及内皮通透性增加有关。此外,差异基因进行共表达网络分析发现Sod2处于核心位置,提示其可能与LPS诱导的rCMECs的各种变化密切相关。结论LPS调控了rCMECs中大量与炎性免疫细胞进入心肌组织相关基因的表达。     

葡萄糖诱导绿色杜氏藻转录组及相关通路差异分析

绿色杜氏藻是研究耐盐机理的模式绿藻．葡萄糖不仅是营养物质,而且还是信号物质．目前,对绿色杜氏藻转录组、糖处理后差异表达基因和β-胡萝卜素生物合成途径关键基因表达还不清楚．本研究通过Illumina HiSeqTM 2000高通量测序,获得葡萄糖处理和未处理绿色杜氏藻转录组信息．利用P value值和差异倍数对样本进行差异表达分析,共111条转录本存在差异表达,3条为上调转录本,108条为下调转录本．利用RT-qPCR检验差异表达分析的准确性. 结果表明,转录本表达结果与转录组分析结果一致．GO功能富集结果表明,71条下调转录本与代谢相关,占所有下调转录本的65.74%．KEGG富集分析结果表明,21条KEGG通路含89条下调转录本,14条通路与代谢相关．代谢中通路最多的为能量代谢（6条）,含63条下调转录本．能量代谢中与光合作用相关的下调转录本最多,为29条．通过分析找到2条与β-胡萝卜素生物合成相关通路（MVA/MEP途径及β-胡萝卜素合成途径）,并发现通路的关键基因hmgs、dxs、dxr、psy、pds、chyb,对其进行差异表达分析,均不存在差异表达．研究表明,葡萄糖抑制了绿色杜氏藻光合作用,代谢受阻,但未影响β-胡萝卜素生物合成相关通路及关键基因．     

白化菠萝蜜茎2个蛋白激酶和4个转录因子家族分析

为探究蛋白激酶(PKs)和转录因子(TFs)在白化菠萝蜜(Artocarpus heterophyllus)幼苗茎次生生长中的表达变化,基于转录组数据对其差异表达基因(DEGs)进行预测及分类,并对挑选出的2个PKs和4个TFs家族构建系统进化树。结果表明,胞质类受体激酶(RLCK)-VIII家族的DEGs上下调表达各4个,亮氨酸富集重复类受体激酶(LRR-RLK)-X家族Xa和Xb-2分支中的DEGs均下调表达,Xb-1中的均上调,TCP家族的20个DEGs中有15个上调表达,zf-HD和GRF家族中的大多数DEGs上调表达,Alfin-like家族中的DEGs均下调表达。因此,这表明6个家族可能在菠萝蜜茎的次生生长过程和应对非生物胁迫中发挥重要作用。     

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

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

不同温度条件下鲁氏耶尔森氏菌的链特异性转录组分析

为鉴定鱼源鲁氏耶尔森氏菌(Yersinia ruckeri) SC09菌株水生环境中不同温度的转录组水平上的差异, 研究采用链特异性转录组测序(Strand-specific RNA-seq)技术对菌体生理温度(28℃)和实验培养温度(37℃)下进行链特异性测序, 原始数据质控后, 筛选得到差异表达基因, 通过KEGG (Kyoto Encyclopedia of Genes and Genomes)数据库对差异表达基因进行富集分析, 并利用Rockhopper软件筛选出的重要原核生物基因簇进行验证。结果显示, 共筛选获得173个显著差异表达基因(P<0.05), 其中包括58个上调基因, 主要富集到一些特殊的碳水化合物代谢相关的通路中; 以及115个下调基因, 主要富集到双组份信号系统中与三羧酸循环相关的代谢通路上, 同时部分基因富集到编码鞭毛素相关的基因簇中。结果表明, 相对于37℃的实验室培养温度, 在水生环境的生理温度条件下(28℃) SC09菌株拥有较高的运动性和较强的葡萄糖代谢, 但相对的SC09菌株代谢一些特殊糖类的能力减弱。     

白桦BpNAC012基因调控白桦木质部发育表达谱分析

NAC转录因子成员被认为是调控植物次生壁合成的开关。前期研究结果显示BpNAC012基因的表达能够调控白桦次生细胞壁的合成。为研究BpNAC012调控的下游靶基因,本研究分别以该基因的过表达,抑制表达株系茎为材料构建转录组,以野生型为对照分析差异表达基因。结果显示:与对照相比,过表达株系OE中上调的基因有627条,下调的基因有229条,抑制表达株系中上调的基因有299条,下调表达的基因有207条。过表达BpNAC012相对于抑制表达能够调控更多的基因表达变化。而抑制表达BpNAC012更多的是影响蛋白修饰和转运类基因的表达变化。在差异表达基因中,涉及受体信号通路,营养代谢,氨基酸合成,及苯丙烷生物合成相关代谢通路基因比较富集。BpNAC012能够调控纤维素、木质素合成及木质部发育相关基因的表达变化,同时能够调控多种转录因子的表达变化。该研究为深入分析BpNAC012在白桦次生细胞壁合成的分子调控机制奠定基础。