野生矮牡丹高通量转录组测序分析

为探究野生稷山矮牡丹(Paeonia jishanensis)花青素等抗氧化活性物质生物合成的遗传基础,本文采用高通量测序技术Illumina NextSeqTM500对其幼嫩叶片、嫩茎和花芽混合池样本进行转录组测序分析。共获得39450个unigene,其中3563个unigene中有4106个SSR位点。在Swiss Prot数据库中注释到29420个unigene,其中11个与花青素生物合成过程相关,22个与类黄酮生物合成过程相关,15个与异黄酮生物合成过程相关;NR数据库中注释高达100%,有12个与类黄酮生物合成相关,7个与花青素生物合成相关;GO数据库中注释到24291个unigene,分为生物学过程、细胞组分及分子功能等三大类53个功能组,有21个unigene与黄酮、黄酮醇及异黄酮生物合成有关,4个unigene与花青素合成有关;eggNOG数据库中注释到38238个unigene,涉及植物生长发育过程绝大多数生命活动过程;KEGG数据库注释到5709个unigene,分别定位到6大类别42个代谢途径,其中25个unigene与花青素生物合成相关。研究表明,在稷山矮牡丹转录组序列中,共鉴定和发掘出多个涉及类黄酮化合物生物合成、花青素生物合成过程的相关基因序列及SSR位点,为下一步开展相关代谢合成途径研究奠定了基础。     

基于高通量测序的“紫娟”茶树叶片miRNA分析

"紫娟"是特异的茶树资源,含有丰富的次生代谢物儿茶素、花青素、黄酮等,其生物合成是由多条代谢途径通过相应的节点连在一起的网络途径,受到多种结构基因和调控基因的控制。Micro RNA又称miRNA,是一种非编码RNAs,能通过对基因表达的调控来调节植物生长、发育以及次生代谢等许多方面。本研究利用高通量测序技术分别构建了"紫娟"茶树芽、第二叶、开面叶、成熟叶的miRNA文库并进行测序,鉴定出已知的miRNA 126个,分为26个家族,预测到的新miRNA 119个。基于"紫娟"茶树转录组数据分析已知miRNA和新miRNA分的靶基因,分别预测到靶基因724条和2 285条。预测的靶基因大多为转录因子,包括调控次生代谢物花青素、黄酮类生物合成的MYB、b HLH转录因子等。这些注释信息的完成为后期研究miRNA在调控茶树叶片发育和次生代谢物的生物合成提供理论依据。     

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

采用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类,其中包括多种化合物的代谢途径和次生代谢产物的生物合成途径,其中有较多的次生物质代谢途径,如萜类化合物生物合成途径、黄酮类生物合成代谢途径以及花青素的生物合成途径等。     

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

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

基于转录组测序的滇山茶花叶呈色机理分析

该研究采用Illumina Hi-Seq2500高通量测序技术,对滇山茶的花叶(绿叶区和金叶区)的cDNA进行了转录组测序,分析滇山茶的花叶呈色机理。测序结果得到了228 862条单基因序列,筛选得到了4 851条差异表达基因。2 291条差异表达基因在GO数据库中有功能注释,1 826条差异表达基因被注释到COG数据库,这些差异表达基因被分为23个主要功能类别,其中大部分基因与滇山茶的生长和代谢有关。1 363条差异表达基因注释到了KEGG数据库,进一步筛选出了差异表达基因注释序列最多的15条代谢通路,其中卟啉和叶绿素代谢途径中的GluRS、δ-ALAD、ALAS基因和类黄酮合成途径中ANS、LAR、CHS基因与滇山茶的花叶形成密切相关。研究表明,滇山茶的花叶呈色是由于叶色相关代谢通路中大量的基因表达量发生变化引起的,并不是由单基因突变所致,病毒诱导的多个代谢途径基因沉默可能是滇山茶花叶呈色的主要原因。该研究丰富了滇山茶的转录组信息,并为滇山茶的彩叶品种培育提供了一定的参考。     

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

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

不同鲜食品质橄榄果实转录组测序及酚类代谢途径相关调控基因挖掘

【目的】果实酚类物质类别、含量是与橄榄营养及风味密切相关的重要品质性状,探究橄榄酚类物质生物合成的分子调控机制。【方法】以总酚含量差异显著的橄榄（低酚/高酚）为试材,分别取花后80-160 d的果实进行转录组分析,对莽草酸-水解单宁/苯丙烷-类黄酮生物合成通路差异基因进行表征,并对差异表达基因进行WGCNA分析,从中挖掘与酚类代谢途径相关的转录因子。【结果】转录组测序共获得296 314条Unigene,其中73%的Unigene被注释到数据库;4个品种（系）橄榄成熟果共鉴定到1 628个差异表达基因（DEGs）,KEGG分析显示,DEGs在酚类代谢途径的“类黄酮生物合成”通路上显著富集;进一步对果实成熟过程莽草酸-水解单宁/苯丙烷-类黄酮生物合成途径的DEGs进行表征,结合WGCNA分析中每个模块与性状的R2和P值,筛选到4个关键模块,根据模块内基因的调控关系利用MCC拓扑分析法以度值≥1挖掘模块内关键转录因子,挖掘到137个转录因子Unigene与30个酚类合成结构基因Unigene共表达,转录因子Unigene功能注释来自35个基因家族,最多的为锌指蛋白（C2C2、C3H、C2H...     

陆地棉MYB类转录因子基因GhTT2克隆及功能初步分析

原花青素作为植物重要的次生代谢产物,是植物应对生物和非生物胁迫的一种重要防御手段,也是影响植物发育和品质的重要因素。原花青素作为花青素生物合成的一条末端通路在模式植物中已有研究,但是具体代谢和调控机制尚不明确;原花青素作为棕色棉纤维呈色的主要物质,其棉纤维呈色的生化与分子机制仍未完全阐明。本研究从陆地棉(Gossypium hirsutum)中克隆了一个MYB类转录因子基因GhTT2 (transparent testa 2),并对其基因结构、表达模式、亚细胞定位及功能进行了分析。结果表明：GhTT2转录因子具有典型的MYB结构域,在纤维中优势表达,其转录水平随花青素含量增加而降低;该基因可被原核诱导表达;与GFP融合的重组蛋白定位在细胞核;酵母转化结果表明GhTT2具有转录激活功能;在棉花中沉默GhTT2基因的表达,导致原花青素含量显著降低,表明其可能参与调控陆地棉原花青素的生物合成。本研究结果为深入阐明MYB类转录因子参与调控植物原花青素生物合成途径的分子机制提供参考。     

紫化茶树生理生化及其花青素调控机理研究进展

紫化茶树是一种特异的茶树种质资源,富含具有多种生理活性的花青素。本文综述了紫化茶树生理生化、花青素生物合成途径及影响因素,旨在为进一步探明茶树花青素合成和代谢的调控机理奠定基础。     

黄秋葵花和果荚转录组测序及类黄酮代谢差异表达分析

采用Illumina Hi seq 2500转录组高通量测序,构建黄秋葵花和果荚转录组文库,并利用测序评估、差异基因功能注释等生物信息学方法进行分析。研究结果表明：(1)分别获得黄秋葵花和果荚有效数据7.12 Gb和8.14 Gb,碱基百分比(Q30)均达到91.0%以上。(2)获得差异表达基因(DEGs)1 336个,其中上调基因319个,下调基因1 017个。(3)获得功能注释的基因有1 131个,GO将455注释基因归成41个功能小类,主要涉及代谢过程、催化活性、单一生物过程和细胞过程等过程;KOG注释了472个DEGs,涉及23个功能分类,其中与次生代谢直接相关过程O和Q类别获得111个注释结果;KEGG将372个DEGs注释到80条代谢通路中,获得F3H、F3′5′H、DFR、ANR、ANS、GT、LAR共10个关键差异基因,其中F3H、DFR在黄秋葵花中表现上调效应,F3′5′H、ANR、LAR在黄秋葵果荚中表现显著上调效应,ANS、GT则分别在花和果荚中均有上调或下调效应。(4)实时定量PCR分析显示,其中7个差异表达基因得到的相对表达量与转录组表达谱分析趋势完全一致。(5)类黄酮代谢途径分析表明,黄秋葵花通过F3H、DFR、ANS、GT途径将NAR催化为生成花青素苷;果荚则将NAR通过F3′5′H将催化为DHM,后在FLS催化下生成黄酮醇类物质等;部分NAR在F3H、 DFR催化下,生成无色飞燕草素苷元,再分别在ANS、LAR作用下,进入花青素苷元和原花青素合成途径。该研究结果丰富了黄秋葵转录组信息,为黄秋葵类黄酮物质纯化和功能利用提供参考依据。     

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

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

高丛蓝莓品种“瑞蓝”的幼叶转录组分析

蓝莓营养丰富,富含花青素等抗氧化性活性物质,被誉为"小浆果之王"。为了加深对蓝莓优良品种"瑞蓝"生长发育分子机制的了解,利用第二代测序技术对其2年生苗的幼叶转录组进行测序。获得总共约1 G的有效测序数据,共拼装成41 987条unigenes,其中23 555条unigenes能被GO数据库成功注释,分别涉及到生物学过程、细胞成分及分子功能相关的共35种生理代谢功能;20 379条unigenes能被KOG数据库注视到,共涉及25条代谢通路;9 771条unigenes能被KEGG数据库成功注释,共涉及到5个大类、20个中类、124条代谢通路;27 293条unigenes被Nr数据库成功注释,这些unigenes共对应32 949条蛋白质的编码区。同时,生物信息学分析还显示:全部unigenes中,有5 221条unigenes编码转录因子,涉及到54个家族;2 487条unigenes编码抗性基因,涉及19个家族;共检测到6 261个SSR序列,其中2碱基重复的SSR占68.2%。转录组测序分析结果阐明了"瑞蓝"的生长发育的部分分子信息,对利用其进一步开展分子标记辅助育种等具有较为重要的意义。     

MBW复合体在植物花青素合成途径中的研究进展

花青素是植物呈色的一种关键物质,它既能赋予植物丰富的色彩,又具有广泛的生物学功能,如抗氧化、抗紫外线、抗病虫害等。植物花青素生物合成途径由一系列结构基因和调节基因协同完成。结构基因表达受MYB、bHLH和WD40类转录因子以其组成的MBW复合体调控。综述了近年来MYB、bHLH和WD40类转录因子及MBW复合物在植物花青素合成中的功能研究进展,总结了花青素合成的各种生物过程和调控网络。     

藏茵陈基源植物皱边喉毛花的全长转录组信息分析

皱边喉毛花为藏药藏茵陈基源植物之一,其包含丰富的药用成分。为进一步了解皱边喉毛花转录组,丰富其基因注释、代谢通路等遗传信息,该研究利用PacBio测序平台对皱边喉毛花叶片进行全长转录组测序。结果表明:(1)全长转录组测序共获得17 Gb的高质量数据,对795 698 个CCS序列进行聚类和去冗余,最终获得87 814 条高质量的全长转录本。(2)与7个数据库比对后,共有277 451 条转录本注释成功,其中注释到NR数据库的转录本最多,有39 214 条。26 396 条转录本成功注释到KOG数据库中,共有26 个子类。39 104 条转录本注释到KEGG数据库中,涉及6 个主要通路和40 个子通路。39 102 条转录本注释到GO数据库中,按分子功能、生物学过程和细胞成分3大类对注释成功的转录本进行分类。(3)SSR分析共鉴定到22 861 个SSR,其中单碱基重复最为丰富; 共检测到1 874 个转录因子和15 166 个长非编码RNA(LncRNA),而注释到转录本最多的转录因子家族是C3H。(4)筛选出55 条与单萜类及黄酮类化合物合成相关的转录本。该研究结果丰富了皱边喉毛花的转录组信息,为进一步筛选皱边喉毛花药用成分合成相关的关键基因提供了重要的遗传资源。     

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

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

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

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

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

榆瘿蚜取食侵染榆树叶片形成了榆树虫瘿,本研究采用新一代的高通量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通路。本研究通过二代高通量转录组测序技术研究榆瘿蚜侵染下榆树虫瘿的相关基因,为今后研究榆瘿蚜侵染榆树叶片形成虫瘿的分子机理提供了基础资料。     

微藻(Chlorella sorokiniana)的转录组分析:油脂生物合成相关的途径解析和基因挖掘

【目的】为了构建和改造产油微藻(Chloralla sorokiniana),需要从基因组水平解析脂肪酸、三酰甘油、淀粉等生物大分子的合成和降解途径,并分析和确定其中的重要基因。【方法】本研究选取了初始氮浓度分别为KNO3∶8g/L和KNO3∶2 g/L的培养条件培养微藻C.sorokiniana,培养至84h收集藻细胞进行Illumina Hiseq2000双端测序,利用Trinity进行de novo拼接,得到的转录本通过Nr数据库、UniProtKB/Swiss-Prot数据库、COG数据库进行功能注释及分类,通过KEGG数据库进行相关代谢途径注释。最后利用RSEM计算每个转录本的RPKM值,并对相关代谢途径中的转录本的表达差异进行了初步分析。【结果】Illumina Hiseq2000双端测序共获得49M 100 bp的raw reads,de novo拼接得到49885个转录本,其长度范围在300bp到14100bp之间,N50为1941 bp。其中26479个转录本成功注释到功能,2357个转录本注释到了EC编号,利用这些转录本注释到207条代谢途径。在此基础上构建了C.sorokiniana的脂肪酸、三酰甘油、淀粉等生物大分子的合成及降解途径。并初步确定了代谢途径中基因的上调及下调水平。【结论】通过RNA-seq分析实现了对非模式藻株C.sorokiniana的基因组解析,在此基础上构建了脂肪酸、三酰甘油、淀粉的生物合成和降解途径及重要基因,所构建的代谢途径与模式藻株莱氏衣藻是一致的,并比较了相关途径中的基因的表达差异,这些信息有助于将来对微藻进行遗传改造提高其产油能力。     

文冠果果实转录组测序及分析

为了研究文冠果果实转录组中功能基因表达情况,采样RNA-seq技术对文冠果不同发育时期果实进行测序分析。结果显示:两个时期样品材料测序组装后共获得68 298个Unigene序列,有24 691个Unigene得到注释,占36.15%;GO数据库注释显示14 766条Unigene分为细胞组分、分子功能及生物学过程3大类54个功能组;KOG数据库中注释到的13 994条Unigene功能系统分为25类;以KEGG代谢途径数据库为依据,可将8 284个文冠果果实转录组Unigene分为127个代谢通路;在文冠果果实转录组中发现11 732个SSR位点,其中最多的为单核苷酸SSR,占60.85%。本研究为文冠果果实分子生物学研究提供了一定的基础和参考。