香榧转录组测序及生物信息学基础分析

香榧具有重要的经济价值,但其基因组信息相对匮乏,限制了其分子生物学和基因功能的研究。本文以不同组织的香榧作为研究对象,采用新一代高通量测序技术平台Illumina Hi Seq?2000对香榧转录组进行测序和数据分析,共得到37,349,086个reads片段,总碱基数为4.35 G。利用组装软件,对获得的高质量序列进行组装,共得到104,636个Unigene,平均长度为784 nt,N50为1,702。将Unigene序列与公共数据库进行比对,28,766个Unigenes获得了注释。其中26,856个Unigene在NR蛋白数据库中获得注释,24,003个Unigenes在NT数据库中获得注释,21,401个Unigene在Swiss-Prot蛋白数据库中获得注释,16,137个Unigene在COG数据库中获得注释,11,410个Unigene在GO数据库中获得注释。根据KEGG注释信息,18,564个Unigene被划分到256个代谢途径中。SSR位点搜索发现,在4,217个Unigene中含有4,706个SSR位点。分析所获得的转录组数据,将为香榧功能基因的克隆,基因的表达,指纹图谱构建和分子标记辅助选育奠定基础。     

红脚艾蒿的转录组解析

[目的]利用高通量测序技术解析红脚艾(Artemisia rubripes Nakai)的转录组信息特征.[方法]通过高通量测序平台Illumina HiSeq 2500对红脚艾进行转录组测序,通过Trinity软件de novo组装获得Unigene,并基于序列同源性对Unigene进行功能注释,得到红脚艾的转录组信息.[结果]测序数据经过质控后共获得24126043条高质量的reads,通过de novo组装获得173093个转录本,对组装的转录本去冗余后共获得85991个Unigene,平均长度为616.87 bp,N50为925 bp.共有47216个Unigene在NR、KEGG、COG、KOG、GO数据库获得功能注释,40802个Unigene在NR数据库注释,显示红脚艾与向日葵(Helianthus annuus)的单基因匹配率最高,16846个Unigene被KEGG数据库注释到130条代谢途径中,26171个Unigene被注释到25个KOG功能分类中,23203个Unigene被GO注释到生物过程、细胞组成和分子功能三大类51个功能分类,12810个Unigene被注释到25个COG功能分类中.[结论]利用高通量测序技术获得了红脚艾转录组信息特征,这些数据将为后期开展功能基因鉴定、解析化合物次生代谢途径及其调控机制奠定研究基础.     

基于RNA-Seq的杜仲转录组微卫星特征分析

对杜仲(Eucommia ulmoides)国审良种‘华仲6号’和‘华仲10号’花后70和160d的种仁共4个样本进行转录组测序,对测序数据进行组装和功能注释分类,并对转录组获得的单基因簇(unigene)进行微卫星特征分析。利用新一代高通量测序技术Illumina HiSeq~(TM)2000对杜仲样品进行转录组测序,采用软件Trinity进行组装;利用BLAST软件将unigene序列分别与Nr、GO、COG和KEGG等数据库比对分析;利用MISA软件对转录组的96 469条unigenes进行SSR搜索。结果表明:转录组测序分析,共得到72 791 399个高质量的序列读取片段(Clean reads),包含了14 702 548 161个的碱基序列(bp)信息。对reads进行序列组装,共获得96 469个平均长度为690bp的unigene,序列信息量达到了66.56 Mb。同源性分析结果显示,有49 856个与其它物种同源的unigenes得到注释,占All-unigene的51.68%。将杜仲转录组中的unigene与GO数据库进行比对分析,根据其功能可将注释到的38 983条unigene分成3大类(细胞组分、分子功能和生物学过程)56个分支;根据COG功能可将注释的14 796条unigene基因划分成25个类别;KEGG数据库作为参照,可将注释到的11 260条unigene定位到117个代谢途径分支;SSR位点搜索结果显示,96 469条unigenes中共包含9 621个完整型SSR位点,占总SSR位点的84.14%。完整型SSR位点共包含55种重复基元,其中出现频率最高的重复基序类型为单核苷酸重复中的A/T(4 597个),其次是AG/CT(2 597个)、AT/AT(439个)。     

荸荠优良品种‘桂蹄3号’球茎发育过程转录组测序分析

【目的】开展荸荠球茎发育过程转录组测序研究,为研究荸荠球茎发育过程相关基因表达信息提供参考。【方法】运用高通量测序技术,对荸荠球茎不同发育时期进行转录组测序研究。【结果】组装共得到223 182条转录本和90 542条Unigene,平均长度为809 bp,N50为1119。组装完整性较高,效果较好。对所得Unigene进行不同数据库注释,共有50 583条Unigene成功注释到7个数据库(NR、GO、NT、Pfam、KEGG、KOG及Swiss-prot)。对差异表达基因分析,发现球茎膨大初期的差异表达基因数目最多,为最活跃阶段。GO功能富集分析结果表明,33 205个基因获得功能注释,分为分子功能、细胞组分和生物学过程等3大类和54个亚类。COG功能分类结果表明,17 743个基因分布于25个功能区域,其中碳水化合物代谢占重要地位。KEGG代谢通路注释结果表明有20 667个基因获得功能注释,共有116条代谢途径,其中淀粉-蔗糖代谢占主要作用。【结论】利用高通量转录组测序技术首次建立了荸荠优良品种‘桂蹄3号’球茎的转录组数据库,为进一步研究荸荠球茎淀粉生物合成相关基因的功能及形成的分子机制提供了数据基础。     

基于转录组测序的紫薇SSR特征分析

以紫薇叶片为材料，通过Illumina Hiseq X Ten测序平台进行高通量测序，利用MISA软件分析转录组的SSR相关信息。结果表明：测序共获得31 461条Unigene序列，总长度34 408 463 bp，其中有5 617条Unigene序列含有SSR位点，共发现SSR位点6 761个，出现频率为21.49%，平均分布距离为5.09 kb；主要的优势重复类型为二核苷酸重复和三核苷酸重复，分别占总SSR的44.95%和36.41%;AG/CT和A/T为SSR优势重复基元，出现频率分别为40.23%(2 720个)和16.57%(1120个)。分析表明，紫薇转录组中SSR类型丰富、位点出现频率高、多态性潜力高，在紫薇种质遗传多样性研究、品种鉴定、杂种鉴定等方面具有应用价值。     

基于荧光标记的紫娟茶树转录组EST-SSR标记开发

为开发紫娟茶树转录组EST-SSR标记,基于前期对紫娟茶树芽、第2叶、开面叶、成熟叶转录组高通量测序所得到的242 757条Unigene进行多态性分析与评价,再利用荧光标记PCR技术,规模化开发茶树EST-SSR标记。结果表明:从紫娟茶树转录组中搜索得到46 041条Unigene含有57 976个SSR位点,出现频率为23.88%。EST-SSR类型以一核苷酸、二核苷酸、三核苷酸重复为主,占总SSR的98.54%。设计合成138对荧光SSR引物,采用荧光标记PCR技术对4个差异较大的茶树品种进行引物筛选,其中44对引物得到高质量的EST-SSR标记位点。这些EST-SSR可用于茶树遗传多样性分析、分子育种等。     

基于高通量测序的北柴胡根转录组SSR位点信息分析

以北柴胡根为试材,采用高通量转录组测序方法,获得uingene序列信息,研究分析SSR分布、基元特征,设计、验证EST-SSR引物的有效性,为开发、丰富EST-SSR分子标记奠定基础。结果表明：从北柴胡根转录组中共筛选到31 176个SSR位点,分布于21 732条unigenes,出现的频率为20.08%,分布密度为3.20个/10 kb,主要重复基元类型以二核苷酸最为丰富,共21 056个,占SSR位点总数的67.54%;其次是单核苷酸和三核苷酸,分别占总SSR的19.05%和12.07%;四核苷酸～六核苷酸重复基元数量均较少。不同核苷酸基序类型共有131种基元,重复基元次数在5～11次的数量最多,占SSR总数的91.7%,且长度基本小于15 bp,其中AT/AT和AC/GT这两种基元在二核苷酸中出现频率最高。2 064条和1 004条含有SSR位点的unigenes分别注释到GO和KEGG通路,获得注释的基因功能主要与基础代谢相关。从22 090对具有潜在多态性的EST-SSR引物中,随机挑选20对引物对3个北柴胡种质DNA进行PCR扩增,扩增成功率最高达85%。北柴胡根转录组S...     

杜仲雌雄株转录组测序数据组装及基因功能注释

以10年以上树龄的杜仲雌株当年新发枝条上的幼果、嫩芽、叶片和树皮和雄株新发枝条上嫩芽、叶片和树皮为材料,采用Illumina Hi SeqTM2000高通量测序技术进行转录组测序,获得雌株51,574,000条、雄株52,430,502条Clean Reads数据,分别包含总长度为4,641,660,000nt和4,718,745,180nt核苷酸序列数据信息;经拼接组装,获得雌株基因信息长达69,461,730nt的423,339个Contig片段,获得雄株基因信息长达94,814,201nt的542,383个Contig片段;经进一步拼接,分别获得平均长度为288nt的雌株159,434个Unigene片段和平均长度为231nt的雄株257,288个Unigene片段,共有48,761个表达序列标签(EST)。以BLAST(E-value1.0E-5)将Unigene对NR、NT、KEGG和COG数据库进行比对,获得CDS序列35,541条,再通过ESTscan分析获得CDS片段13,220条,共获得48,761条CDS片段。与NR数据库比对发现杜仲雌、雄株转录组Unigene与葡萄相似序列最多(33.8%),其次是蓖麻(11.4%)和杨树(11.2%),与拟南芥的相似序列仅2.3%;根据Unigene与COG数据库比对结果,可将有COG功能的7,571条Unigene分为24类,而根据GO数据库注释,杜仲转录组有GO功能注释的23,314条Unigene可分为生物过程、细胞组分和分子功能3大类55分支。与KEGG数据库比对,杜仲雌、雄株转录组17,468条Ungenes分属128类代谢通路,其中有2,399条属于次生物质代谢途径,314条参与萜类化合物生物合成途径。     

尾叶紫薇全长转录组测序及微卫星标记开发

[目的]尾叶紫薇是培育香花紫薇的重要种质资源,目前花香功能基因的研究及已开发的SSR分子标记较少,无法满足分子育种研究需求.[方法]采用PacBio sequel平台对尾叶紫薇3个开花时期混样进行全长转录组测序,通过公共数据库对Unigene进行功能注释,利用MISA软件检索全长转录组的SSR位点,Primer3批量设计SSR引物并检测.[结果]测序获得39087条非冗余、高质量全长转录本,平均长度为2109 bp,N50长度为2427 bp,长度主要分布在500～4000 bp.有38559条序列(98.64％)获得功能注释,功能主要集中于细胞进程及代谢进程.其中参与花香相关萜类物质合成的有56个,参与脂肪酸衍生物合成的有922个.获得的16425个SSR位点分布在12060条序列中,平均每5.02 kb出现1个SSR序列,3141条序列含有多个SSR位点.SSR位点主要重复类型是二核苷酸和三核苷酸,分别占所有SSR位点的41.95％和31.58％,重复基序类型有96种,AG/CT重复类型的出现频率最高.SSR的重复次数以5～15次为主,长度为10～86 bp,平均为16.84 bp.随机选取52对SSR引物进行PCR检测,有40对引物能产生单一、清晰的条带,有效扩增率为76.92％.[结论]获得与花香物质合成相关基因序列和大量有效SSR标记,将为深入了解尾叶紫薇花香物质合成与代谢、紫薇属植物种质资源评价与保护、分子标记辅助育种等研究奠定基础.     

冬瓜高通量转录组测序及分析

【目的】获得冬瓜转录组序列、遗传变异等信息,从中挖掘冬瓜基因数据及SSR分子标记,为冬瓜后续研究提供数据支撑。【方法】以冬瓜嫩叶为材料,利用Illumina HiSeq~(TM)2000技术对冬瓜进行转录组测序,构建数据库从中获得干净序列。经De novo拼接组装后,将获得的单基因簇(Unigene)数据在非冗余蛋白数据库(nonredundant protein database,Nr)、蛋白质序列数据库(Swiss Prot protein database,Swiss Prot)、基因本体论数据库(gene ontology,GO)、蛋白质真核同源数据库(eukaryotic orthologous groups,KOG)、东京基因与基金组百科全书(Kyoto encyclopedia of genes and genomes,KEGG)、蛋白质家族域数据库(protein families database,Pfam)6个公共数据库中进行比对,最终得到冬瓜单基因簇注释信息。利用MISA软件对转录组单基因簇进行搜索,获得单基因簇中的SSR位点。【结果】从冬瓜嫩叶中得到62 021 032条高品质序列,组装后获得40 611条单基因簇,平均长度955 bp。将所有单基因簇在Nr和Swiss Prot数据库中进行比对,结果分别比对到27 474及19 573条单基因簇;在GO数据库中,所注释到的10 659条单基因簇分别匹配到生物功能、分子功能和细胞组分3个本体的47个功能组中;与KOG数据库进行注释比对,根据其功能将注释到的单基因簇划分为25类;KEGG数据库比对注释到10 799条冬瓜的单基因簇,可分为5个大类、19个亚类、125条代谢途径;在Pfam数据库中比对到17 990条单基因簇,分属于369个类群。SSR位点搜索发现,有5 086条单基因簇包含SSR序列,获得5 474个SSR位点。【结论】利用高通量测序获得大量冬瓜转录组信息,有助于从分子水平对冬瓜进行深入研究。     

基于转录组测序的花椒属物种EST-SSR标记开发

【目的】开发适用于花椒(Zanthoxylum bungeanum)和竹叶花椒(Zanthoxylum armatum)的EST-SSR引物,为花椒属物种的鉴定及遗传多样性探究提供分子标记。【方法】对花椒和竹叶花椒分别进行转录组测序(RNA-seq),基于得到的EST序列开发SSR引物。在2种花椒的EST-SSR引物中分别随机挑选60对引物,用12份材料(4份花椒和8份竹叶花椒)为样本,利用琼脂糖凝胶电泳验证其特异性。从有特异性条带的引物中,再随机选取花椒和竹叶花椒SSR引物各15对,选取6个(2份花椒和4份竹叶花椒)样本,利用聚丙烯酰胺银染电泳进行多态性检测。【结果】花椒共有64 944条Unigene,碱基对长度共54 073 890bp,含有12 746个SSR位点,分布在10 595条Unigene上,SSR位点出现频率为19.63%,平均分布距离4.24kb。竹叶花椒Unigene共75 669条,碱基对长度58 975 053bp,共15 096个SSR位点,分布在12 612条Unigene上。SSR位点出现频率为19.95%,平均分布距离为3.91kb。60对花椒引物中,有46对成功扩增出特异性条带,扩增效率76.67%;60对竹叶花椒引物中,有41对扩增出特异性条带,扩增效率68.33%。花椒和竹叶花椒的特异性条带大小在100～300bp,主要集中在150～250bp。多态性验证试验中,15对花椒引物中有12对产生了多态性条带,多态率80.00%;15对竹叶花椒引物中则有14对产生了多态性条带,多态率93.33%。【结论】基于花椒和竹叶花椒转录组高通量测序结果开发的EST-SSR分子标记引物,具有较明显的特异性和多态性。     

基于转录组测序的花椒属物种EST-SSR标记开发

【目的】开发适用于花椒（Zanthoxylum bungeanum）和竹叶花椒（Zanthoxylum armatum）的EST-SSR引物，为花椒属物种的鉴定及遗传多样性探究提供分子标记。【方法】对花椒和竹叶花椒分别进行转录组测序(RNA-seq)，基于得到的EST序列开发SSR引物。在2种花椒的EST-SSR引物中分别随机挑选60对引物，用12份材料（4份花椒和8份竹叶花椒）为样本，利用琼脂糖凝胶电泳验证其特异性。从有特异性条带的引物中，再随机选取花椒和竹叶花椒SSR引物各15对，选取6个（2份花椒和4份竹叶花椒）样本，利用聚丙烯酰胺银染电泳进行多态性检测。【结果】花椒共有64 944条Unigene，碱基对长度共54 073 890 bp，含有12 746个SSR位点，分布在10 595条Unigene上，SSR位点出现频率为19.63%，平均分布距离4.24 kb。竹叶花椒Unigene共75 669条，碱基对长度58 975 053 bp，共15 096个SSR位点，分布在12 612条Unigene上。SSR位点出现频率为19.95%，平均分布距离为3.91 kb。60对花椒引物中，有46对成功扩增出特异性条带，扩增效率76.67%；60对竹叶花椒引物中，有41对扩增出特异性条带，扩增效率68.33%。花椒和竹叶花椒的特异性条带大小在100~300 bp，主要集中在150~250 bp。多态性验证试验中，15对花椒引物中有12对产生了多态性条带，多态率80.00%；15对竹叶花椒引物中则有14对产生了多态性条带，多态率93.33%。【结论】基于花椒和竹叶花椒转录组高通量测序结果开发的EST-SSR分子标记引物，具有较明显的特异性和多态性。     

辣木转录组SSR位点信息分析

采用MISA对15 824条长度大于1kb的辣木Unigene进行SSR位点分布频率和基本特征的分析,共检测到8 272个SSR位点,分布于6 003条Unigene,SSR位点出现频率为52.28%,共包含114种重复单元。辣木转录组的SSR以单核苷酸重复为主,占总SSR的55.69%。其次是二、三核苷酸重复,分别占总SSR的32.34%和11.24%。A/T、AG/CT和AAG/CTT是单核苷酸、二核苷酸和三核苷酸中的优势重复基元。辣木转录组SSR以10~20次重复为主,基序长度主要集中于12~19bp。结果表明辣木转录组中含有大量SSR,且类型丰富,可为辣木SSR引物开发提供候选序列。     

c-GMP诱导对盐胁迫下番茄的转录组分析

以市场常见的番茄为材料,采用Illumina高通量测序技术对NaCl胁迫、c-GMP诱导以及两者组合处理下番茄幼苗进行转录组测序。结果表明:共获得83.35 Gb的原始数据,Q30碱基百分比在90.91%及以上,分别将各样品的原始读数与指定的参考基因组进行序列比对,比对效率从88.03%到92.74%不等。同时将得到的31 734条序列注释到不同的数据库(GO、KEGG、KOG、NR、Pfam、Swiss-Prot和egg NOG),通过GO数据库,将其划分到细胞组分、分子功能以及生物学过程三大注释类别下的51个亚功能组中,通过COG数据库与KEGG数据库将其分别注释到25个类别与128个代谢通路。在NaCl胁迫、c-GMP诱导以及两者组合处理下共发现SNP位点127 944个,其中仅有22.6%的SNP位点位于非编码区。差异基因结果显示,c-GMP能诱导盐胁迫下更多基因的表达。该转录组测序数据可靠,结果覆盖面广,为番茄c-GMP诱导盐胁迫下基因挖掘与功能研究提供了理论依据。     

基于Illumina高通量测序的岩原鲤转录组分析

为获得岩原鲤转录组信息,发掘功能基因,本研究采用Illumina高通量测序技术对岩原鲤全组织转录组进行测序。结果获得64 257 918个EST序列,经拼接和组装得到83 252条单基因序列(unigene),平均长度787 bp,长度范围201~16 572 bp。利用NCBI的蛋白质非冗余数据库(Nr)对所有unigene进行相似性搜索,共有37 157条unigene(44.63%)与数据库中的已知序列同源。利用Blast2GO v2.5软件对unigene进行注释,共得到29 919条(35.93%)注释基因,根据GO功能分类将其分为生物过程、细胞组分和分子功能3大类56亚类。经KOG注释及分类,共有17 869条(21.49%)unigene成功注释到真核直系同源组中,并将其分为26个功能组分。经KEGG代谢通路分析可分为5大类(细胞过程、环境信息处理、遗传信息处理、代谢和有机系统)32小类共267个代谢通路。本研究通过高通量测序技术,对岩原鲤转录组进行测序,获得了大量的转录组信息,为岩原鲤功能基因克隆及基因组学研究提供了基础。     

鲫鱼血液、肝脏、卵巢组织转录组比较分析

为进一步丰富鲫鱼基因组数据，对肝脏、血液和卵巢组织进行高通量转录组测序。采用Illumina Hiseq高通量测序平台进行转录组测序，对所获得序列进行比对、基因注释、差异表达和SNP、SSR分析。结果提示，组装得到127 801条unigene,平均长度为735 bp。与相应数据库比对，最终获得117 414(91.87%)个有注释信息的unigene,有105条unigene获得GO注释，并分类到生物学过程、细胞组分和分子功能三大类42个功能中；有20 725条unigene被归纳到KOG的26个直系同源功能类型；获得KEGG注释的unigene共参与五大类32小类代谢通路。在鲫鱼肝脏、血液和卵巢组织中分别检测到307 740、343 516、348 795个SNP位点，共检测到48 808个SSR位点。提示获得有注释信息的unigene共117 414条，此外，获得一些SNP和SSR位点。结果可为进一步克隆和挖掘鲫鱼功能基因、多态性检测及群体遗传多样性分析等方面研究奠定基础。     

椰子织蛾转录组分析

为获取入侵性检疫害虫椰子织蛾的基因信息,利用高通量测序技术对该害虫进行了转录组测序。结果表明:经拼接共获得了37 416条unigenes,平均长度为803 bp,N50为1 415 bp;通过BlastX比对,有19 154条unigenes (51.19%)成功在Nr中被注释。在Nr中成功注释的unigene,与黑脉金斑蝶同源性最高,达到68.45%;在GO注释中,共有14 131 (37.76%)条unigenes被成功注释到细胞组分、分子功能和生物过程3大类的54个分支;在COG数据库中,成功注释的8 606条unigenes被注释到26个分类;在KEGG中,共有6 171条(16.49%) unigenes被注释到5大类通路的260条途径上。基于构建的转录组数据库,鉴定得到3 248个微卫星。     

山羊5个组织的全长转录组构建与分析

【目的】构建山羊（Capra hircus）全长转录组,旨在改进参考基因组的功能注释。【方法】利用PacBio长读段测序技术对南江黄羊肝、脾、肾周脂肪、背最长肌和睾丸5个组织的混合样本进行高通量测序;对原始数据进行校正、聚类去冗余,最后用KEGG数据库对预测转录本序列进行功能注释。【结果】质控后获得28 432条高质量一致性全长序列,其中99.49%的序列被成功比对到山羊基因组;得到源自9 879个基因的19 115个转录本。进一步分析得到表达基因总数的41.39%发生了可变剪接;其中起始外显子可变剪接占比最高（36.47%）,外显子互斥占比最低（1.78%）。另外,144条序列未被比对到参考基因组上,这些序列的平均GC含量高达56.50%;利用KEGG数据库对其进行功能注释,揭示28个转录本与免疫系统相关。【结论】获得的测序数据全面涵盖混合样本中表达的转录本,证明了山羊基因组转录的复杂性。研究结果总体上提高了每个基因平均表达的转录本数量,证明可变剪接的广泛存在;未成功比对的转录本序列GC含量较高且许多转录本来自免疫基因的表达。     

基于RNA-Seq技术的鲮转录组分析

为满足标记辅助育种的要求,通过454测序平台首次开展了鲮Cirrhina molitorella全鱼转录组深度测序,并用 Newbler 等软件进行数据精细分析。结果表明：共获得了1297479条 reads,总碱基数为486586191 bp,组装后得到19962条contigs,平均长度为1269 bp, N50为1509 bp。基因功能注释研究共获取了10577个特异蛋白,根据特异蛋白注释结果进行GO分析,有7314条contigs有GO注释,包含5381个特异蛋白;采用GO 功能分类工具可将已注释转录物序列划分为分子功能、生物途径和细胞成分3类,为下一步开展生长等性状相关基因功能验证研究提供丰富的序列资源;共鉴定出5931个具有完整的ORF的全长cDNA序列,并且鉴定出2438个微卫星和5014个SNP位点。本研究中,还建立了鲮转录组数据库和网站,方便同行随时调取数据,这为深入开展鲮分子标记辅助的遗传育种、种群遗传学和资源评估等研究提供了丰富的标记资源。     

藤茶高通量转录组分析及黄酮类化合物合成相关基因挖掘

[目的]分析藤茶高通量转录组序列,从中挖掘出黄酮类化合物合成相关基因,为进一步揭示藤茶黄酮类化合物生物合成调控机制提供理论参考.[方法]分别采集藤茶的幼叶和成熟叶,提取其总RNA构建cDNA文库,采用Il-lumina HiSeqTM 4000高通量测序平台对藤茶叶片进行转录组测序,经过滤处理后运用Trinity组装,将获得的Unigene与Nr、Nt、Pfam、Swiss-Prot、GO、KO和KOG 7个数据库进行比对注释,并预测Unigenes的编码区序列(CDS);基于KEGG信号通路富集分析,发掘藤茶黄酮类化合物合成相关基因.[结果]藤茶叶片转录组测序获得82126236条原始测序序列(Raw reads),过滤处理后得到80156972条高质量序列(Clean reads),进一步组装拼接得到92472条Unige-nes,平均长度为1208 bp,N50长度为1780 bp,其中,至少在1个数据库注释的Unigenes有84217条,占Unigenes总数的91.07%,有8944条Unigenes在7个数据库均被注释,占Unigenes总数的9.67%.在GO数据库成功注释的41116条Unige-nes可分为生物学过程、细胞组分和分子功能三大类,共56个小类;在KOG数据库注释的14553条Unigenes可分成25类,其中,一般功能预测注释成功的Unigenes最多(1946条);其次是翻译后修饰、蛋白质翻转、分子伴侣(1776条),参与次生代谢物质的生物合成、转运和降解的Unigenes较少,仅有319条;KEGG信号通路富集分析发现,共有15262条Unigenes注释到128条KEGG信号通路,以注释为代谢的Unigenes最多,为8694条,其中筛选获得有98个黄酮类化合物合成相关基因,分别编码苯丙烷代谢通路的3种关键酶和类黄酮代谢通路的14种关键酶.藤茶叶片转录组Unigenes与Swiss-Prot和Nr数据库比对,获得52582条CDS序列,ESTScan 3.0.3预测获得35535条CDS序列.[结论]藤茶在细胞过程、代谢过程、单有机体过程、细胞和细胞部分、结合和催化活性能力分布的基因较丰富,在一般功能、翻译、翻译后修饰、蛋白质翻转及分子伴侣的基因表达量较高,具有较强的碳水化合物代谢能力.多种关键酶基因参与藤茶黄酮类化合物的生物合成,推测其生物合成途径存在多条分支,调控机制也较复杂.