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

该研究以‘铁观音’茶树品种的种子为试验材料,采用转录组测序技术分析种子发育的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在幼果期至果实膨大期差异上调表达,在果实膨大期至果实成熟期差异下调表达。     

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

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

沙葱萤叶甲卵滞育的转录组学分析

【目的】建立沙葱萤叶甲Galeruca daurica滞育卵转录组数据库,挖掘卵滞育相关的基因以及代谢和信号通路,在转录组水平探讨卵滞育的分子机制。【方法】采用Illumina NovaSeq6000高通量测序平台对沙葱萤叶甲滞育卵与解除滞育卵进行转录组测序,并进行生物信息学分析;利用DESeq软件分析沙葱萤叶甲滞育卵与解除滞育卵中的差异表达基因,对差异表达基因进行KEGG通路富集分析;利用qRT PCR技术对10个差异表达基因的表达模式进行验证。【结果】基于沙葱萤叶甲滞育卵与解除滞育卵转录组测序结果,共获得53 389个unigene,其中差异表达基因2 145个,24个差异表达基因与保幼激素信号及脂肪酸生物合成和降解相关。与解除滞育卵相比,滞育卵转录组中1 297个基因上调表达,富集于124条KEGG通路,其中核糖体通路显著富集;848个基因下调表达,富集于73条KEGG通路,其中MAPK信号通路和糖胺聚糖生物合成通路显著富集。qRT-PCR结果表明,随机选取的10个差异表达基因的表达趋势与RNA-Seq转录组测序结果完全一致。【结论】保幼激素,脂肪酸生物合成和降解,核糖体,MAPK信号及糖胺聚糖生物合成等通路可能在沙葱萤叶甲卵滞育调控中起着重要的作用。     

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

为了研究烟草脉斑驳病毒（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...     

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

为探究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...     

黑果枸杞果实成熟发育过程表达谱差异分析

对黑果枸杞(Lycium ruthenicum Murr.)3个不同发育时期的果实进行转录组测序,分析对比黑果枸杞果实不同发育时期相关基因表达谱的变化。提取果实RNA进行Illumina高通量测序,利用GO、KEGG等公共数据库进行功能注释、分类,利用数字基因表达谱技术对比分析。结果显示,KEGG pathway富集分析表明亮氨酸生物合成途径在变色期对比青果期有7个差异表达基因,黑熟期对比变色期有35个差异表达基因,且全部为上调表达。在糖酵解途径中3个不可逆的反应步骤限速酶基因在黑熟期均为上调表达。对黑果枸杞亮氨酸合成和糖酵解涉及到的基因进行了功能和代谢通路的分析,为黑果枸杞种质资源探索提供理论基础。     

伞裙追寄蝇滞育关联基因的转录组学分析

【目的】本研究旨在探讨调控伞裙追寄蝇Exorista civilis滞育的重要关联基因以及代谢途径,为明确该虫在转录组水平下的滞育分子机制提供理论基础。【方法】采用新一代高通量测序平台Illumina HiSeq^TM2000对伞裙追寄蝇非滞育以及滞育的蛹进行转录组测序分析以及生物信息学分析;应用KAAS在线pathway比对分析工具对筛选出的满足padj＜0.05且fold change≥32或padj＜0.05且fold change≤1/32条件的差异表达基因进行KEGG通路富集分析。【结果】根据测序结果,共获取58 050个基因。差异倍数在32倍以上的滞育关联基因(diapause-associated genes, DAGs)有454个,其中406个表达上调,共涉及134条通路,包括氧化磷酸化、柠檬酸循环及其他重要通路;48个表达下调,涉及32条通路。KEGG通路富集分析结果表明,滞育关联基因主要集中在信号转导、内分泌系统、碳水化合物代谢等途径中。【结论】本研究获得的伞裙追寄蝇转录组数据揭示了其滞育调控的重要代谢途径与关联基因。     

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

绿色杜氏藻是研究耐盐机理的模式绿藻．葡萄糖不仅是营养物质,而且还是信号物质．目前,对绿色杜氏藻转录组、糖处理后差异表达基因和β-胡萝卜素生物合成途径关键基因表达还不清楚．本研究通过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,对其进行差异表达分析,均不存在差异表达．研究表明,葡萄糖抑制了绿色杜氏藻光合作用,代谢受阻,但未影响β-胡萝卜素生物合成相关通路及关键基因．     

镧对铜胁迫下水稻转录组模式的调控分析

探究稀土镧对铜胁迫下水稻转录组的影响,鉴定镧在调控水稻铜胁迫应答中的关键基因和功能。利用抗氧化酶活性筛选最适铜胁迫及镧处理浓度,进行转录组测序和差异表达分析,并用qRT-PCR验证差异表达水平。对测序数据进行差异表达分析发现3 222个基因显著上调,3 798个显著下调,qRT-PCR验证了4个细胞壁防御相关基因CHIT13、Laccase、Expansin和GH3.4。功能和通路分析获得95组GO条目和112条KEGG通路,富集于激酶和氧分子结合等分子功能、细胞壁及质膜等细胞组分、次级代谢和脂质代谢等生物学过程以及苯丙烷生物合成和植物激素信号转导通路。镧通过调节细胞壁形成和组分提高水稻铜胁迫耐受性。     

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

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

基于油桐种子3个不同发育时期转录组的油脂合成代谢途径分析

油桐是我国重要的木本油料植物, 过去对油桐的研究主要集中于栽培和常规育种, 与油桐种仁油脂合成相关的分子机理研究还未见报道。文章采用RNA-seq技术对油桐种子油脂合成的3个不同时期的转录组进行比较, 获得了大量差异表达的Unigene序列。在此基础上, 通过GO分类和Pathway富集性分析将这些差异表达Unigene归类于128个代谢途径, 其中包含与油脂合成相关的脂肪酸生物合成和甘油磷脂代谢途径。桐酸经脂肪酸生物合成途径合成后通过甘油磷脂代谢途径以桐油的形式贮存。将这两个代谢途径的Unigene序列在KEGG数据库中进行比对, 获得了一些关键酶的同源蛋白质。文章通过对编码这些同源蛋白质的基因在油桐种子油脂合成期的表达模式进行分析, 以期为油桐油脂合成, 尤其是桐酸合成机理的解析提供理论参考, 并为油桐的遗传改良提供潜在的基因资源, 从而提高油桐的单位面积产量。     

芒果胶孢炭疽病菌应答菌丝机械损伤产生无性孢子的分子机制

胶孢炭疽菌(Colletotrichum gloeosporioides)是引发芒果(Mangifera indica)炭疽病的主要病原体。室内平板培养胶孢炭疽菌不产生或产生很少分生孢子的情况时有发生, 但菌丝在机械损伤后24-48小时会产生大量分生孢子。胶孢炭疽菌应答机械损伤诱导产孢的核心基因及关键代谢通路尚未见报道。基于转录组测序(RNA-seq)技术检测了芒果胶孢炭疽菌菌丝在机械损伤处理后2小时内5个时间点的基因表达变化, 对差异表达基因进行GO富集和KEGG代谢通路富集分析, 并对菌丝响应胁迫的基因动态表达数据进行分析。基于常微分方程ODE模型结合变量选择技术, 构建了动态基因调控网络。结果表明, 有417个差异表达基因参与应答胶孢炭疽菌菌丝机械损伤, 分属12个聚类模块, 有4条通路存在显著富集, 分别是丙酮酸代谢、硫代谢、黄曲霉素合成途径和二萜合成途径。结合功能注释筛选出12个应答菌丝损伤胁迫的核心基因。研究结果为后续深入开展芒果胶孢炭疽菌产孢和致病机理研究奠定了重要基础。     

低温胁迫下两种不同形态群体拟穴青蟹转录组分析

为研究拟穴青蟹(Scylla paramamosain)不同形态群体对温度变化的适应差异,以Sp1群体(腕节外缘中部都不具有发达的刺)和Sp2群体(至少一只螯足腕节外缘中部具有发达的刺)为研究对象,选用低温胁迫组(8℃)和对照组(20℃)拟穴青蟹肝胰腺为实验材料,进行转录组测序。结果显示,测序共获得81853个Unigene,平均长度为420 bp, N₅₀为1460 bp,其中22.33%为已知基因。差异表达基因筛选发现Sp1群体在低温胁迫下共有15773个Unigene显著差异表达,其中15628个Unigene表达下调, 145个Unigene表达上调,这些Unigene显著富集于细胞质基质等生物学过程和RNA结合、核糖体结构组分等分子功能,并进一步富集在核糖体、剪接体和内质网蛋白质加工等信号通路。Sp2群体在低温胁迫下共有323个Unigene显著差异表达,其中上调114个,下调209个。这些Unigene显著富集于DNA整合等生物学过程和甾体羟化酶活性、碳酸氢盐跨膜转运体活性等分子功能,并进一步富集在赖氨酸降解、甘氨酸、丝氨酸和苏氨酸的代谢、精氨酸和脯氨...     

RNA-seq用于获得共轭亚油酸对贵妃鸡肌内脂肪代谢调控的关键基因

通过转录物组测序获得在贵妃鸡基础日粮中添加共轭亚油酸(CLA)对肌内脂肪代谢的差异表达基因,经生物信息学分析获得相关的信号通路及可能发挥重要作用的候选基因,为CLA对肌内脂肪沉积的分子机制奠定基础。本研究选用55日龄健康的贵妃鸡为试验动物,在基础日粮中添加CLA 0%、1%和2%,预饲期1周,正饲期6周。屠宰采集胸肌组织进行转录物组测序,对测序数据进行差异表达分析,差异表达基因GO功能和差异表达基因KEGG通路富集分析,筛选出与胸肌脂类代谢相关的差异表达基因,利用qRT-PCR对差异表达基因进行验证。结果显示,共获得1 065个差异表达基因,其中上调基因703个,下调基因362个。GO富集结果显示,差异表达基因主要富集在生物过程的细胞过程、单一生物过程、生物调节和代谢过程。KEGG信号通路富集显示,差异表达基因显著富集在黏着斑、不饱和脂肪酸生物合成、脂肪酸生物合成和类固醇生物合成等信号通路中,发现11个主要与肌内脂肪代谢相关的候选基因,分别是FADS1、FADS2、ELOVL5、ACOX2、SLC27A1、FABP5、LPL、LOC107050163、ENSGALG00000030996、ENSGALG00000005043和ENSGALG00000048882。并随机选取6个基因进行qRT-PCR验证,其相对表达量变化趋势与测序结果一致。本研究筛选到CLA影响贵妃鸡胸肌脂类代谢相关的差异表达基因,并对11个主要参与脂肪代谢相关的基因进行分析,为揭示CLA调控肌内脂肪沉积的分子机制奠定基础。     

转录组测序揭示油桐脂肪酸的合成途径

油桐是我国重要的木本油料植物,过去对油桐的研究主要集中于栽培和常规育种,与油桐种仁油脂合成相关的分子机理研究还未见报道。本研究采用转录组测序（RNA-Seq）技术,比较了油桐种子油脂合成3个不同时期的转录组,获得了大量差异表达的Unigene序列。通过GO（GeneOntology）分类和代谢途径富集性分析将这些差异表达的Unigene归类于包含脂肪酸生物合成途径在内的128个代谢途径。随后将脂肪酸生物合成途径中的54个Unigene序列在KEGG（KyotoEncyclopediaofGenesandGenomes）数据库中进行比对,获得了14个关键酶的同源蛋白质。本研究通过对编码这些同源蛋白质的基因在油桐种子油脂合成期的表达模式进行分析,以期为油桐油脂合成,尤其是桐酸合成机理的解析提供理论参考,并为进一步的理论研究和油桐的遗传改良提供潜在的基因资源,从而提高油桐的单位面积产量。     

利用RNA-seq数据分析肝细胞癌差异表达基因及蛋白相互作用

RNA-Seq已成为当前转录组学研究的强有力工具,尤其在肿瘤差异表达基因的筛选方面有重要的应用价值。为进一步阐明肝细胞癌(HCC)的分子机制,本研究对GEO中1个包括12对HCC组织标本的RNA-Seq数据集(GSE63863)进行了生物信息学分析。采用edgeR、DESeq2、voom等3种不同算法的软件进行统计分析,共获得976个差异表达基因(adj. p-value<0.01或FDR<0.01,|logFC|≥2),其中上调表达422个(43.2%),下调554个(56.8%)。GO富集分析显示这些差异表达基因主要涉及离子结合、氧化还原酶活性等分子功能以及氧化还原、细胞分裂等生物学过程;KEGG通路分析显示,这些差异表达基因主要涉及细胞周期、视黄醇等代谢通路。STRING分析显示,共有654个基因编码的蛋白质存在相互作用,进一步利用MCODE分析显示,169个基因编码蛋白构成4个子网络,相应的中心节点基因分别为UBE2C、GNG4、TTR、FOS,这些基因的异常表达可能在HCC的发生发展过程中具有重要作用。上述研究结果将为进一步阐明HCC分子发病机制、寻找新型生物标志物提供初步的依据。     

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

为更好地挖掘八角(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家族序列数量最多。该文利用转录组测序技术分析八角优良无性系与普通品种叶片的差异基因及...     

叶子花花瓣状苞片和叶片的比较转录组学研究

叶子花(Bougainvillea spectabilis Willd.)的苞片色彩鲜艳,呈花瓣状,在景观园艺和环境美化有着广泛应用。为了进一步了解花瓣状苞片的演化机制,该研究采用RNA-seq技术分别对苞片和叶片进行测序,比较苞片和叶片在基因表达水平的差异。最终获得54.48 M对长度为150bp的双向clean reads,从头拼接出55 782条Unigenes,平均长度为1 003.79bp,N50为1 402bp。共有30 874条Unigenes能被注释到公共数据库,其中有12 631条Unigenes能被注释到基因本体。差异化表达分析确定苞片和叶片差异表达的基因共有456个。该研究分别对不同组织差异表达的基因进行GO富集分析和KEGG富集分析,结果表明,在苞片中高表达的差异化基因有黄酮类合成相关的基因,花色合成关键酶多巴双加氧酶(DODA)以及发现与花瓣发育相关的基因,在叶片中高表达的差异化基因主要是参与光合作用。