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

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

基于转录组测序的羽衣甘蓝叶色相关基因分析

叶色是羽衣甘蓝重要的观赏性状之一。本研究采用Illumina Hi-Seq2500高通量测序技术,对基因型纯合的紫叶和白叶羽衣甘蓝叶片进行转录组测序,筛选差异基因并与GO和KEGG数据库比对进行注释分析,分析羽衣甘蓝叶色形成相关基因。结果显示,获得高质量短读序共104 608 770条,筛选出紫叶相对白叶的差异表达基因1 993个,其中上调表达基因1 094个,下调表达基因899个。根据GO功能分类可分为生物过程、细胞组分和分子功能3大类64功能组。根据KEGG代谢通路分析可以分为171类,在叶色相关的类黄酮生物合成途径中黄酮醇合成酶(FLS)、二氢黄酮醇4-还原酶(DFR)上调以及类胡萝卜素生物合成途径中的类胡萝卜素β-环化酶上调与紫叶形成关系密切。本研究丰富了羽衣甘蓝的转录组信息,获得了一些差异表达基因,为进一步研究羽衣甘蓝叶色形成的遗传机制提供了有价值的信息。     

多年宿根甘蔗的转录组分析及差异基因挖掘

甘蔗(Saccharum hybirds)宿根性直接关系到甘蔗生产成本和种植效益,品种、环境和栽培措施均会影响甘蔗宿根能力,但品种的种性是影响宿根性最关键因素。国内外从分子生物学层面解释甘蔗宿根性差异的文献鲜见报道,从转录水平分析不同宿根年限GR2号和ROC22号基因表达的差异。结果表明：转录组测序共得到100 558条转录本和25 582条Unigene。获得53 790条Unigene的注释结果,分别在NR、Swiss-Prot、KEGG、COG、KOG、GO和Pfam数据库进行比对。GO功能注释共分成三大类,51小类,6年宿根蔗注释到1 029个差异基因,3年宿根蔗注释到3 391个差异基因。主要KEGG代谢通路有8条,分别涉及植物激素信号转导,淀粉和蔗糖代谢,甘氨酸、丝氨酸和苏氨酸代谢,苯丙素生物合成,同源重组,DNA复制,错配修复及植物病原体相互作用。筛选出脱落酸(ABA)相关差异基因6个,分别为脱落酸不敏感蛋白2基因(ABI2)、bZIP转录因子超家族蛋白、碱性亮氨酸拉链型转录因子基因(ABI5)、aba响应元件结合因子1基因(ABF1)、G-box结合因子基因(GBFs)...     

基于转录组测序的宁夏枸杞不同品种果实活性成分合成差异表达基因分析

为探究不同品种宁夏枸杞果实活性成分生物合成相关基因的表达水平,筛选关键差异表达基因(differentially expressed genes,DEGs),揭示宁夏枸杞品种间活性成分含量差异的分子机制,本研究采用Illumina NovaSeq 6000高通量测序技术,对宁夏枸杞‘宁杞1号’和‘宁杞7号’青果期、转色期及成熟期果实进行转录组测序,比较2个品种果实不同发育期相关基因表达谱的变化。结果显示:转录组测序共获得811818178条clean reads,有121.76 Gb有效数据。‘宁杞1号’和‘宁杞7号’在青果期、转色期和成熟期差异表达基因分别有2827、2552和2311个;分别有2153、2050和1825个差异基因在基因本体论(gene ontology,GO)、京都基因与基因组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)富集分析和同源蛋白簇(clusters of orthologous groups of proteins,KOG)分析等6个数据库中被成功注释。青果期、转色期和成熟期果实的差异表达基因,在GO数据库分别有1307、865和624个被富集到生物学过程、细胞组分及分子功能3个部分中;KEGG通路富集结果均集中在代谢途径、次生代谢物生物合成和植物-病原互作过程;在KOG数据库,3个发育期分别注释了1775、1751和1541个差异表达基因。对注释的基因进行PubMed数据库检索,在青果期、转色期和成熟期分别筛选到与枸杞活性成分合成相关的差异表达基因18、26和24个,这些基因主要参与类胡萝卜素、类黄酮、萜类、生物碱和维生素等代谢途径。选取7个差异表达基因进行RT-qPCR验证,结果与转录组测序数据表达趋势一致。本研究从转录水平为不同品种宁夏枸杞活性成分含量差异提供了初步证据,为进一步挖掘枸杞活性成分生物合成的关键基因及解析其表达调控机制提供了研究基础。     

基于转录组分析酸枣仁斯皮诺素生物合成途径的相关基因

采用转录组测序技术分析酸枣仁(Ziziphi spinosae semen)不同发育时期基因表达差异,识别酸枣仁黄酮生物合成途径的相关基因,为进一步分析酸枣仁中黄酮类成分生物合成相关酶基因提供研究基础.以3个不同生长发育时期的酸枣仁样本为供试材料,采用Illumina Hiseq 4000进行转录组测序,并从头拼接,对转录本进行功能注释和差异性分析.经StringTie软件组装拼接得到74472个转录本,其中有26.83％在NR、Swiss-Prot、Pfam、GO、KOG、KEGG数据库中得到注释;代谢通路分析结果发现,在酸枣仁转录组中与黄酮类合成途径相关的11条差异表达基因可能编码斯皮诺素合成途径的一些酶类.酸枣仁中黄酮类成分斯皮诺素相关转录本信息为阐明斯皮诺素在生长期的积累提供重要依据,为提高酸枣仁品质研究提供参考.     

草菇耐低温突变菌株Vtlt-1和原始菌株V23转录组差异的研究

为了培育草菇耐低温菌株与解析其耐低温的分子机制,采用紫外诱变的方法,选育出耐低温草菇菌株Vtlt-1,并利用表达谱芯片技术,比较突变菌株Vtlt-1与原始菌株V23的表达差异基因,筛选差异显著基因后利用GO（gene ontology）功能注释和KEGG（kyoto encyclopedia of genes and genomes）通路富集等方法分析,结果发现与V23相比Vtlt-1表达显著差异基因共有1 600个,其中704个基因上调表达,896个基因下调表达。针对差异基因的GO功能分类结果发现：生物学过程方面,差异基因主要分布在金属离子结合、氧化还原过程、碳水化合物的代谢与核酸的结合。细胞组分方面主要与细胞核相关;分子功能中,氧化还原活性以及解旋酶活性,依赖于ATP 的解旋酶活性、DNA指导的RNA聚合酶活性。KEGG注释结果发现差异表达基因主要富集在氨基酸与氮类物质的代谢、脂肪酸与生物碱的合成这两方面,此外还富集到核糖体的生物合成、细胞色素P450、RNA聚合酶、硫和氨基酸的代谢、核酸的修复等通路。这些结果为解析草菇耐低温机制提供分子依据和理论基础。     

宁乡猪皮下脂肪与肌内脂肪组织转录组差异分析

为了比较分析宁乡猪皮下脂肪(subcutaneous fat,SAF)和肌内脂肪(intramuscular fat,IMF)组织脂肪沉积的分子机制,本文利用RNA-seq技术鉴定和分析宁乡猪SAF和IMF组织中差异基因表达谱。选取6头250日龄健康状况良好、种内个体体重相近(约85 kg)的雄性宁乡猪为实验材料,采集SAF与IMF组织样品。通过对两个脂肪组织转录组测序并进行GO(Gene Ontology)功能注释及KEGG(Kyoto Encyclopedia of Genes and Genomes)信号通路富集分析,得到与脂肪沉积和脂质代谢相关的差异基因。为验证测序数据结果的可靠性,本研究随机选取6个差异基因进行qRT-PCR验证。结果表明,宁乡猪SAF和IMF组织中有2406个基因差异表达,其中1422个基因表达上调,984个基因表达下调。通过GO功能注释分析发现,差异表达基因主要通过参与类固醇生物合成、不饱和脂肪酸的生物合成、甘油磷脂代谢及自噬途径等与脂质代谢相关途径,调控宁乡猪SAF和IMF的沉积。KEGG通路富集结果显示,差异基因主要富集在脂质结合、脂肪酸代谢过程、甘油...     

京海黄鸡柔嫩艾美耳球虫感染后盲肠转录组分析

为了探究京海黄鸡感染柔嫩艾美尔球虫(E. tenella)后盲肠组织差异表达基因,以及球虫感染分子应答过程和免疫应答机制,试验采用RNA-seq技术对E. tenella感染和非感染组第7天的盲肠组织进行转录组测序,筛选差异表达基因,并进行差异基因的功能、通路富集分析。结果表明,在感染和非感染组中有显著差异的表达基因2 830个(P0.05),其中1 419个基因上调,1 411个基因下调。随机选取10个差异基因进行qRT-PCR验证,结果显示差异基因的表达倍数与RNA-seq检测结果显著相关(r=0.988,P0.000),决定系数达0.975。GO分析表明,有2 356个差异基因获得GO功能注释,显著富集的前30个GO terms主要涉及细胞交流、信号转导、血管生成、氧化还原酶活性等。KEGG分析发现差异基因显著富集的信号通路有黏着斑、细胞外基质-受体相互作用、过氧化物酶体增殖物激活受体等。这些通路中的差异基因有ANGPTL4、ACSL5、VEGFC、CD44和MAKP10等,提示这些基因在宿主柔嫩艾美耳球虫感染过程中发挥重要作用。     

枫香叶片变色期全长转录组测序及分析

枫香因其树形优美,入秋后叶色红艳或橙黄,极具观赏价值,是优良的景观生态树种。为了解枫香叶片变色及其次级代谢过程的遗传基础,该文以枫香5个叶片变色期叶片混合样品为材料,利用单分子实时测序技术(PacBio平台)对其进行全长转录组测序。结果表明:(1)全长转录组测序共获得41.04 Gb的高质量数据,从中鉴定出全长非嵌合序列563 180条,通过聚类和去冗余,获得27 269条高质量全长转录本。在27 269条全长转录本中预测到2 035条长链非编码RNA(lncRNA),并检测出14 892个简单重复序列(SSR)位点和1 856个转录因子。(2)基因注释结果表明,NR、GO、COG、KEGG 等8个数据库共注释了24 857条转录本,KEGG数据库共获得了124个条代谢途径,主要有核糖体、碳代谢、氨基酸生物合成等,在类黄酮和叶绿素代谢途径中分别有49和71个转录本参与。上述结果初步揭示了枫香叶片变色期转录组信息以及功能特性,为后续研究枫香叶片变色分子机制、色素代谢合成途径和调控、相关功能基因克隆以及叶色改良提供基础数据。     

侧孢短芽孢杆菌S2-31拮抗下细辛叶枯病菌转录组差异表达分析

通过分析侧孢短芽孢杆菌拮抗作用下叶枯病菌的转录组学特征,研究差异表达基因(DEGs)和代谢通路的富集情况,初步探索侧孢短芽孢杆菌拮抗叶枯病菌的分子机制。首先利用S2-31与叶枯病菌的对峙培养观察其拮抗作用,然后利用转录组测序探究侧孢短芽孢杆菌拮抗下和正常生长下的叶枯病菌的基因表达水平差异,并进行RT-qPCR验证,最后,利用相关数据库对DEGs进行注释和富集分析。对峙培养后,叶枯病菌菌丝出现皱缩、扭曲等现象。测序后共得到3681个DEGs,其中2224个相对上调,1457个相对下调。GO功能注释分析表明,上调和下调表达的DEGs均注释到生物过程8个亚群、细胞组分8个亚群和分子功能4个亚群。KEGG富集分析中,共有732个unigenes定位到115条生物学通路中,其中富集程度相对较高的通路有氨基糖和核苷糖代谢、糖基磷脂酰肌醇(GPI)锚定生物合成和过氧化物酶体;涉及差异基因较多的通路有淀粉和蔗糖代谢、剪接体和胞吞作用;上调表达的DEGs主要富集在代谢途径,下调表达的DEGs主要富集在遗传信息处理和细胞过程途径。从显著富集的通路中筛选出与菌丝生长发育相关的差异表达基因,结果表明侧孢短芽孢杆菌主要抑制病原菌菌丝过氧化物酶体的形成、胞吞过程和GPI锚定的合成等过程。差异表达基因的RT-qPCR验证结果与转录组测序的结果一致。在侧孢短芽孢杆菌S2-31的拮抗下抑制了叶枯病菌的生长,其转录组特征也发生显著变化,差异基因主要涉及代谢、细胞过程以及遗传信息处理等途径的相关通路。     

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

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

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

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

基于转录组测序的铁皮石斛植物甾醇生物合成相关基因分析

为了解铁皮石斛(Dendrobium officinale)植物甾醇的生物合成途径,利用Illumina HiSeq 4000高通量测序技术对茎和叶进行转录组测序,比较植物甾醇生物合成关键酶基因的表达。结果表明,共获得43 085条Unigenes,其中24 459条在Nr、Swiss-prot、KOG和KEGG数据库获得注释,7 333条获得共同注释。KEGG代谢途径分析表明,铁皮石斛植物甾醇生物合成分为3个阶段,共有50个Unigenes (30种酶)参与。表达分析表明,DXR和HMED的表达量明显高于MK和MVD;成熟期茎、叶SMT1的表达量比生长期高,SMT2则生长期高于成熟期;同一时期,SMT1和SMT2在叶的表达量都比茎高。这为铁皮石斛植物甾醇的开发利用和调控植物甾醇生物合成研究提供了科学依据。     

樟树5种化学类型叶片转录组分析

樟树(Cinnamomum camphora )是樟科植物的一个代表种, 具有材用、药用、香料、油用和生态环境建设等多种用途。叶精油中富含利用价值极高的樟脑、芳樟醇、1,8-桉叶油素、异-橙花叔醇和右旋龙脑等萜类化合物。依据叶精油中主要成分的种类和含量, 可将樟树划分为脑樟、芳樟、油樟、异樟、龙脑樟5种化学类型。文章采用Illumina HiSeq™ 2000高通量测序技术, 对5种化学类型叶片转录组进行测序, 对测序得到的所有Unigene进行GO(Gene Ontology)、COG(Clusters of Orthologous Groups)和KEGG(Kyoto Encyclopedia of Genes and Genomes)分类, 给出功能注释和Pathway注释, 并预测Unigene蛋白编码区(Coding sequence, CDS)。De novo组装共获得156 278个Unigene, 序列平均长度584 bp, N50(覆盖50%所有核苷酸的最大Unigene长度)为1 023 bp。通过与其他核酸、蛋白数据库的Blast搜索比对, 共有55 955条Unigene获得了基因注释, 占所有Unigene的35.80%。其中, 有24 717条Unigene得到GO注释, 有21 806条Unigene得到COG注释。KEGG pathways分析结果表明, 共有3 350条基因(10.19%)注释到次生代谢生物合成途径, 其中参与单萜、二萜、倍半萜和萜类骨架合成的Unigene有424个。在单萜合成的代谢通路中, 有9条Unigene可能编码芳樟醇合成酶基因, 且表达分析结果显示, 芳樟醇合成酶基因在芳樟化学类型中优势表达, 在油樟化学类型中表达水平较低。这些注释信息的完成为樟树功能基因及相关候选基因的发掘提供了基础数据和重要依据。     

Transcriptome analysis of Verticillium lecanii based on RNA-Seq technology

[目的]蜡蚧轮枝菌是一种防治植物病原物和害虫的生防菌,弄清其致病的分子基础具有重要的意义。[方法]利用Illumina HiSeq技术测序蜡蚧轮枝菌的cDNA文库,并进行RNA-Seq序列拼接和功能注释。[结果]共获得1634670条高品质reads,碱基GC含量48.50%。14856条Unigene经组装拼接后,在Nr、Swiss-Prot、COG和KEGG数据库中分别比对得到1467、9379、6518和4188条Unigenes。采用GO分类工具将21403条Unigenes分成24个功能组,主要包含在细胞组成（1369）、分子功能（1679）和生物学过程（1928）3个大类里,在COG数据库注释的基础上,把6518条Unigenes分成38个功能组,通过GO和COG丰度识别分类,聚类的主要是一般功能预测、次生代谢产物的合成、运输和分解代谢以及信号转导机制起作用的独立基因。能与KEGG路径匹配的108条Unigenes中大部分和新陈代谢途径及次生代谢产物的合成有关。[结论]这些结果将有助于找寻蜡蚧轮枝菌的致病因子,从而丰富其致病机理。     

基于高通量测序的阿尔泰蝠蛾(鳞翅目:蝙蝠蛾科)幼虫转录组分析

本研究采用Illumina HiSeq TM 2500测序平台对阿尔泰蝠蛾Hepialus altaicola Wang幼虫进行转录组测序及生物信息学分析.经序列拼接后共获得100133个Unigenes,总长度86319112 bp,平均长度862 bp,N50长度1628 bp.将Unigenes与NR、COG/KOG、Pfam、Swiss-Prot、GO、KEGG数据库比对,共获得38198条Unigenes,其中Nr数据库注释的Unigenes最多,为32381条,占32.34％.通过GO功能分类,共有13216个Unigenes在GO数据库中细胞组分、分子功能和生物学过程等3大类57个分支中找到注释;KEGG通路分析,共有15058条Unigenes被注释,归属于305条代谢通路.CDS预测发现54002条序列可被编码,占全部基因的53.93％.基因注释进一步获得311个与冷适应相关的代谢调节基因,并用FPKM值对基因表达量进行评估.本研究获得的转录组信息及分析结果,为进一步研究阿尔泰蝠蛾的基因功能及低温生态适应性奠定分子基础.