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

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

基于转录组分析铜绿假单胞菌DN1降解荧蒽特性

[背景]铜绿假单胞菌DN1是一株从石油污染土壤中分离筛选到的具有广谱降解功能的菌株。[目的]深入了解荧蒽胁迫条件下铜绿假单胞菌DN1降解污染物过程中重要的降解相关基因信息。[方法]通过高通量测序技术对铜绿假单胞菌DN1进行转录组测序,对其所有的转录本进行KEGG (kyoto encyclopedia of genes and genomes)分类和Pathway注释、GO (gene ontology)分类和富集分析。[结果]转录组测序显示:与对照组相比,荧蒽诱导组检测到6 189个基因,其中1 919个基因上调表达,1 603个基因下调表达。KEGG注释分析显示差异上调表达基因匹配到了112个KEGG代谢途径,注释到"代谢途径"的1 408个基因(约占总差异基因的73.4%)中有317个基因参与了碳氢化合物代谢及含有苯环结构的异源生物质的生物降解,占"代谢途径"的16.53%,暗示了菌株DN1降解荧蒽可能与这些途径有密切关系。另外,主要代谢途径中的差异表达基因主要集中在ABC转运系统、氨基酸生物合成、双组分系统及碳代谢,这些途径大多数参与了底物的识别转运、信号转导及基因表达调控。[结论]进一步拓展了铜绿假单胞菌DN1在荧蒽胁迫条件下的代谢途径和逆境反应,也为微生物修复环境污染物研究夯实了理论基础。     

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

枫香因其树形优美,入秋后叶色红艳或橙黄,极具观赏价值,是优良的景观生态树种。为了解枫香叶片变色及其次级代谢过程的遗传基础,该文以枫香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个转录本参与。上述结果初步揭示了枫香叶片变色期转录组信息以及功能特性,为后续研究枫香叶片变色分子机制、色素代谢合成途径和调控、相关功能基因克隆以及叶色改良提供基础数据。     

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

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

朝仓花椒幼芽转录组测序数据组装及基因功能注释

采用2代Illumina Hi-Seq测序技术对朝仓花椒雌株结果期幼芽的转录组进行测序,构建了转录组数据库,获得49 672 080条Clean Reads数据,包含总长度为4 967 208 000 nt序列数据信息;经拼接组装,获得转录组基因信息长达55 902 243 nt的176 407个Contig片段,再通过进一步拼接,共获得平均长度为878 nt的96 475个Unigene片段。与Nt、Nr、Swiss-Prot、COG、GO、KEGG等数据库进行BLAST信息比对(E-value艽10-5),共获得68 315个注释基因,以及62 150个表达序列标签(EST)。与NR数据库比对,发现花椒转录组基因序列与同科柑橘属的甜橙(Citrus sinensis)和克莱门柚(Citrus clementina)具有较高的同源性,分别为43.96%及41.86%,与其他物种的同源性较低,均不足10%;可将花椒转录组的Unigene的功能通过与COG数据库进行注释比对划分为25类;根据GO数据库的注释,可将其分为生物过程、细胞组分和分子功能3大类共55分支,根据与KEGG数据库的比较分析发现朝仓花椒的转录组数据中含有代谢通路相关基因128类,其中包括多种化合物的代谢途径和次生代谢产物的生物合成途径,其中有较多的次生物质代谢途径,如萜类化合物生物合成途径、黄酮类生物合成代谢途径以及花青素的生物合成途径等。     

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

为探究不同品种宁夏枸杞果实活性成分生物合成相关基因的表达水平,筛选关键差异表达基因(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验证,结果与转录组测序数据表达趋势一致。本研究从转录水平为不同品种宁夏枸杞活性成分含量差异提供了初步证据,为进一步挖掘枸杞活性成分生物合成的关键基因及解析其表达调控机制提供了研究基础。     

建兰花色形成的分子调控机理研究

为探究建兰花色形成的分子调控机理,该研究以同株建兰黄绿色、红色花瓣为实验材料,采用高通量Highseq测序技术进行转录组文库构建,从基因水平探索花色物质的合成代谢通路及关键调控基因的转录活性。结果表明:(1)转录组测序共获得131 110 030条过滤序列(Clean Reads)、106 479条单基因(Unigenes),通过NR、GO、COG、KEGG等公共数据库比对,获得了29 748个有注释信息的Unigenes。(2)建兰黄绿色花瓣与红色花瓣中包括20个类黄酮合成代谢相关差异表达基因,与黄绿色建兰花瓣相比,红色花瓣中776个基因转录表达量上升,589个基因转录表达量下降,在KEGG数据库被注释到93条代谢通路,涉及6条类黄酮合成代谢相关的途径,共20个差异表达基因(83Unigenes)。(3)QRT-PCR分析显示,所选20个基因在建兰2种颜色花瓣中的表达量比值趋势与转录组FPKM比值趋势一致,表明该研究获得的转录组数据具有较高的参考价值;其中分支酸、苯丙氨酸解氨酶、肉桂酸4-羟化酶、4-香豆酸:辅酶A连接酶基因表达上调,有利于类黄酮前体积累;查尔酮合成酶、二氢黄酮醇还原酶、花青素合成酶基因在黄绿色花瓣中几乎不表达,在红色花瓣中表达明显上调,可能与建兰花色形成相关。     

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

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

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

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

Pigment variation and transcriptional response of the pigment synthesis pathway in the S2309 triple-color ornamental kale (Brassica oleracea L. var. acephala) line

Ornamental kale is popular because of its colorful leaves and few studies have investigated the mechanism of color changes. In this study, an ornamental kale line (S2309) with three leaf colors was developed. Analysis of the anthocyanin, chlorophyll, and carotenoid contents and RNA-seq were performed on the three leaf color types. There was less chlorophyll in the white leaves and purple leaves than in the green leaves, and the anthocyanin content was greatest in the purple leaves. All the downregulated DEGs related to chlorophyll metabolism were detected only in the S2309_G vs. S2309_W comparison, which indicated that the decrease in chlorophyll content was caused mainly by the inhibition of chlorophyll biosynthesis during the leaf color change from green to white. Moreover, the expression of 19 DEGs involved in the anthocyanin biosynthesis pathway was upregulated. These results provide new insight into the mechanisms underlying the three-color formation.     

榆树叶片形成虫瘿过程中转录组学分析

榆瘿蚜取食侵染榆树叶片形成了榆树虫瘿,本研究采用新一代的高通量Illumina Hi Seq^TM 2000技术测序平台对榆瘿蚜取食刺激的榆树叶片进行转录组测序和功能注释,利用生物学方法对基因表达和功能进行研究。测序获得23.19 Gb碱基序列信息,通过对测序数据进行序列过滤、拼接和去冗余,共获得102 017个Unigenes,通过NR与BLAST等数据库比对,其中有37 899个(37.15%) Unigense被注释。利用KOG、GO、KEGG等数据库对榆树虫瘿叶片的Unigense进行比对,按功其能将匹配的Unigenes基因划分25大类;GO注释将信息归纳为基因的3大主类,57个亚类;以KEGG数据库为参考,将Unigene定位到110个不同的代谢通路,包括氧化应激防御、植物激素信号转导、碳水化合物以及次生物代谢等代谢相关的Unigenes通路。本研究通过二代高通量转录组测序技术研究榆瘿蚜侵染下榆树虫瘿的相关基因,为今后研究榆瘿蚜侵染榆树叶片形成虫瘿的分子机理提供了基础资料。     

三叶木通转录组分析及三萜皂苷生物合成途径关键酶基因的发掘

为了解三叶木通(Akebia trifoliata(Thunb.)Koidz.)的三萜皂苷合成途径及其关键酶,本研究对其花、叶、根、茎进行转录组测序,组装获得了57.25 Gb数据,含140 859个unigenes,序列平均长度为1350 bp。KEGG代谢通路富集结果显示,517个unigenes参与三萜皂苷合成相关的3条代谢途径,其中415个unigenes编码三萜皂苷生物合成途径的19个关键酶。对三萜皂苷生物合成过程中的关键酶角鲨烯环氧酶(SE)进行序列分析和同源建模,发现其具有保守的底物结合结构域。将三叶木通茎与花、叶、根的基因表达水平进行比较,发现茎与根相比较其上调基因数目最多,其中295个差异表达基因(DEGs)与三萜皂苷生物合成途径相关。     

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

叶子花(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)以及发现与花瓣发育相关的基因,在叶片中高表达的差异化基因主要是参与光合作用。