杜鹃花叶片转录组测序数据组装及功能注释

采用2代 Illumina Hi-Seq 测序技术对2年生杜鹃花白凤4号(Rhododendron pulchurum cv. BaiFeng 4)扦插苗叶片的转录组进行测序,获得了213 723 424条 Clean reads数据,经过质控和De novo。组装共获得平均长度为930 nt的53 568个Unigenes。与Nr、Swiss-Prot、KOG、KEGG四大数据库进行BLAST信息比对(E-value≤10^-5),共获得28 877个注释基因。与Nr数据库序列同源性比较发现,杜鹃花与葡萄(Vitis vinifera)具有较高的同源性,而与其他物种的同源性较低。杜鹃花转录组的Unigenes在KOG数据库比对上30 508个注释,分为25类。在GO数据库比对上17 218个Unigenes,其中与抗逆有关的Unigenes有11 928个。在KEGG数据库的132条代谢通路中富集了6 475个杜鹃花Unigenes,其中注释到植物激素生物合成代谢途径和信号转导途径的有384个。同时,有1 062个Unigenes被注释为转录因子,有8 738条SSR标记被挖掘。     

茶树转录组中SSR位点的信息分析

利用茶树全转录组的高通量测序获得的127 094条Unigenes来发掘茶树转录组SSR功能性标记。在这些序列中共搜索出12 242个SSRs,分布于10 325条Unigenes中,出现频率为9.63%。茶树转录组SSRs的平均长度为16 bp,平均分布频率是1/3.68 kb。在茶树转录组的SSRs中,二核苷酸重复是主要的类型,占总SSRs的63.78%。茶树转录组SSRs共包含181种重复基元,二核苷酸重复基元CT/AG和TC/GA是优势重复基元类型,分别占总SSRs的23.84%和23.58%。同时对这些SSR的可用性进行了评价。     

基于棉花转录组测序的SSR分子标记的开发

通过对棉花转录组进行测序获得的57 695条Unigenes,利用MISA软件进行搜索,得出简单重复序列(simple sequence repeat,简称SSR)的分布情况,共挖掘检索得到1 886个SSR位点,SSR的出现频率为3.27%。所有SSR位点分布于1 759条Unigenes上,发生频率为3.05%,平均每20.8 kb会出现1个SSR位点。有117条Unigenes含有1个以上SSR位点,含有复合型SSR的Unigenes数目为56条,其中三核苷酸基元类型分布最多,SSR数量为1 582个,占挖掘出总SSR数量的83.88%,而单、二、四、五、六核苷酸重复类型所占的比例较小。最后通过Primer 3程序进行SSR引物设计,得到部分引物序列,可用于棉花遗传多样性分析、分子标记辅助育种以及种质资源的保存等。     

罗氏沼虾转录组SSR标记信息分析

对罗氏沼虾肝胰腺组织转录组测序获得的33 450条Unigenes,进行简单重复序列(simple sequence repeat,简称SSR)检测与分析。结果共检测出15 356个SSR位点,分布在12 930条Unigenes上,发生频率为38.65%,平均距离为1.91 kb。其中,优势重复基序为单核苷酸、二核苷酸、三核苷酸。SSR序列重复长度主要为12～16 bp和18 bp,占罗氏沼虾总SSR的86.33%。单核苷酸重复基序中多数为A/T,二核苷酸重复基序以GA/TC、AG/CT及AT/TA为主,三核苷酸重复基序主要为AAT/TTA、AAG/TTC、ATT/TAA、CCT/GGA、CTC/GAG和CTT/GAA。最后,对SSR所在的Unigene进行了KOG分类注释。通过对罗氏沼虾转录组SSR位点信息分析,获得多态性较为丰富的SSR位点,具有极大的开发价值,可为罗氏沼虾遗传多样性的分析、遗传图谱的构建及分子标记辅助育种提供有效工具。     

基于茶树芽叶转录组序列的EST-SSR分布特征研究

采用生物信息学方法对茶树(Camellia Sinensis)叶部不同发育阶段(单芽/三叶)转录组中的EST-SSR位点信息进行分析,以期为茶树分子标记辅助育种提供有效的参考。利用SSRFINDER软件,对茶树芽叶转录组的高通量测序获得的97454条Unigenes序列(100.91Mb)进行大通量SSR位点的筛选,共获得36249个SSR位点,分布于27913条Unigenes序列中,其发生频率为28.64%,平均分布密度为1/2.78 kb。在茶树芽叶的转录组序列中共发现962种碱基重复基元,其中占主导的是以(AG/CT)n为主(占总SSRs的23.78%)的二核苷酸重复,占到总SSRs的59.19%,其次是三核苷酸和单核苷酸重复,其占比分别为20.92%和13.13%。茶树芽叶转录组所含不同重复基元SSRs的平均长度,相对全器官转录组较短,其中占比最高的是重复长度为18 bp的中长重复序列,占到总SSRs的23.37%,而大于25 bp的较长序列重复极少。同时对茶树不同器官转录组测序的SSR分布特性进行了分析比较,以期为后续的分子标记开发提供参考。     

基于甘薯耐盐转录组测序的SSR和SNP特征分析

旨在对甘薯耐盐转录组中潜在的SSR和SNP位点进行挖掘及特征分析，以完善甘薯分子标记。试验以不同耐盐甘薯品种的转录组序列为基础，利用MISA、GATK软件分析了转录本中的SSR和SNP信息。结果显示，研究共获得SSR位点33192个，分布在24323条Unigenes中，出现频率为21.11%。其中，6271条Unigenes含有超过1个以上的SSR位点。SSR类型以单核苷酸重复SSR和双核苷酸重复SSR为主。在重复基元中，A/T和AG/CT所占的比例最高，为优势基元。在157252条Unigenes中共挖掘到SNP位点7691906个，分布密度为0.08个/bp。其中转换类型有4729922个，颠换类型有2961984个。在所有突变类型中，C/T和A/G含量最高，分别为占总数的32.34%和29.15%。综上所述，甘薯转录组具有较丰富的SSR和SNP标记，多态性较高，可为甘薯抗逆品种筛选、种质培育等研究奠定理论基础。     

四球茶转录组SSR位点信息分析

对基于高通量测序拼接获得的99 741条四球茶嫩叶Unigenes进行了SSR(simple sequence repeats)位点分析,共揭示了23 732个SSR位点,分布于18 552条Unigenes中,SSR发生频率为23.79%;含有SSR位点的Unigenes的总长度为226 188 bp,SSR位点间的平均距离为1/2.07 kb,重复长度平均长度为9.53 bp。在四球茶转录组SSR中,二核苷酸重复和三核苷酸重复是主要的重复类型,分别占总SSRs的47.53%和25.92%;在181种重复基元类型中,优势重复基元分别是AG/CT、A/T、ACC/GGT,分别占总SSRs的41.38%、22.10%和6.49%,共占总SSRs的比例达69.97%。     

基于转录组测序技术的黄唇鱼SSR分子标记筛选

[目的]筛选出可用于黄唇鱼(Bahaba flavolabiata)遗传资源评价及亲缘鉴定的SSR分子标记,为开展其保护遗传学研究打下基础.[方法]利用Illumina HiSeq 4000测序平台对野生黄唇鱼样本进行转录组测序,采用TGICL对转录组数据进行聚类去冗余以获得Unigene,并通过生物信息学方法进行功能注释分类、代谢通路和SSR特征分析等.[结果]黄唇鱼混合组织的转录组测序共获得13.43 Gb数据,经组装并去冗余后获得65047条Unigenes.采用BLAST对组装获得的Unigenes进行七大功能数据库(NR、NT、GO、COG、KEGG、Swissprot和Interpro)注释,结果共有55583条Unigenes被注释(占85.45%).通过MISA对黄唇鱼的Unigene进行SSR搜索,共检测到29797个SSR位点分布于19664条Unigenes中.SSR位点分布类型及其特征分析结果显示,最主要的重复类型为二核苷酸重复基元(11855个,占39.79%),其次是单核苷酸重复基元(9695个,占32.54%)和三核苷酸重复基元(6663个,22.36%),而四核苷酸、五核苷酸和六核苷酸重复基元的数量相对较少.二核苷酸重复基元类型中重复最多的基序是AC/GT(8355条),在三核苷酸重复基元类型中重复最多的基序是AGG/CCT(1866条).从随机选取的186对SSR引物中,最终筛选出47对扩增稳定性好且具有多态性的SSR引物,以二核苷酸和四核苷酸重复基元的扩增结果居多.[结论]通过转录组测序能有效筛选出具有多样性的黄唇鱼SSR分子标记,可为黄唇鱼的遗传多样性分析、亲缘鉴定和遗传图谱构建等后续研究提供基础资料.     

茶树生长素合成酶基因YUCCA10克隆与表达分析

【目的】探究YUCCA10基因在正常花和不育花中的表达规律,为深入探索YUCCA10基因和生长素合成对茶树花器官发育的调控机制提供理论基础。【方法】前期研究茶树正常花与不育花转录组差异时,获得一条与YUCCA10基因高度同源的基因序列,利用‘云茶1号'茶树全长转录组数据库以及PCR技术克隆验证茶树YUCCA10基因,并对其生物信息学特征进行分析,同时运用qRT-PCR技术研究其表达特征。【结果】经验证YUCCA10基因含有1个1137 bp的开放阅读框,编码379个氨基酸,分子量为42.63 kDa,等电点为8.88,茶树YUCCA10蛋白具有NADPH结合位点和FAD结合位点,与多种植物YUCCA10蛋白同源。qRT-PCR分析CsYUCCA10在正常花中随着花的生长发育表达量增加,而在不育花的花瓣、雄蕊、雌蕊以及不同发育期中的表达均低于正常花。【结论】茶树CsYUCCA10基因在不育花中不能正常合成,从而影响了花的发育导致不育。     

卵形鲳鲹卵巢不同发育时期的转录组测序分析

【目的】鉴定筛选出与卵形鲳鲹卵巢发育相关的候选基因及信号通路,为揭示其卵巢性成熟过程的分子机制打下基础。【方法】挑选卵巢发育处于?期和Ш期的雌性卵形鲳鲹,分别构建卵形鲳鲹卵巢?期和Ш期的cDNA文库,采用Illumina HiSeqTM 2500进行转录组测序,经过滤、质量控制及拼接组装后获得的Unigenes在七大数据库（Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO）中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行功能注释及信号通路富集分析,并采用MISA和GATK3进行SSR鉴定及SNP分析。【结果】卵形鲳鲹卵巢组织转录组测序获得的325156432条Raw reads,经过滤筛选得到317206752条Clean reads,拼接组装后得到59554条Unigenes;69.65%的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等七大数据库中注释成功,其中有24599条Unigenes被注释到GO数据库,15997条Unigenes被注释到KEGG数据库。在卵形鲳鲹卵巢组织的2个发育时期共鉴定获得56115个基因,经差异表达分析后获得17737个差异基因,其中8169个基因在卵巢Ш期上调表达、9568个基因在卵巢Ш期下调表达。GO功能注释分析发现,卵形鲳鲹卵巢差异表达基因主要注释在细胞过程、氮化合物代谢过程、初级代谢过程、核、核部分、离子结合及水解酶活性等条目上;而KEGG信号通路富集分析结果显示,17737个差异表达基因显著富集在318条代谢途径上,其中前20条KEGG信号通路包括2-氧代羧酸代谢、PI3K-Akt信号通路、甲状腺激素信号通路、磷脂酶D信号通路、Fc εRI信号通路和细胞周期等。卵形鲳鲹卵巢转录组（59554条Unigenes）中共存在30133个SSRs和82490个SNPs。【结论】GnRHR、FSHR、FSHβ、CYP11A、SIRT3和PEG3等差异表达基因及PI3K-Akt信号通路和VEGF信号通路等与卵形鲳鲹卵巢的发育密切相关,共同调节卵巢的发育与成熟,在卵巢性成熟过程中发挥重要作用。     

粘虫转录组中SSR位点的信息分析

粘虫Mythimna separata(Walker)是一种迁飞性害虫,严重危害玉米、水稻、小麦等粮食作物。SSR是指以1~6个核苷酸为基本重复单位的串联重复DNA序列。SSR位点的信息分析为粘虫扩散、迁飞和交配等行为分子机制的研究以及粘虫的综合防治奠定理论基础。本研究基于高通量测序获得的粘虫转录组数据,利用软件msatcommander发掘粘虫SSR位点。结果从20 776条转录组Unigenes中共搜索出400个SSR,分布于372条Unigenes中。在粘虫转录组SSR中,三核苷酸重复的数量最为丰富,有271个;其次是二核苷酸和单核苷酸重复,分别是70个和49个;四至六核苷酸重复的数量都很少,共10个。粘虫转录组SSR共包含24种重复基元,其中CCG/CGG是优势重复基元类型,有69个;其次是AAG/CTT,有57个。CG/CG有18个,在二核苷酸重复基元中所占的比例达到25.7%。此研究发掘到的SSR位点将为粘虫遗传图谱的构建、遗传多样性分析、亲缘关系分析等提供丰富的分子标记。     

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

[目的]分析藤茶高通量转录组序列,从中挖掘出黄酮类化合物合成相关基因,为进一步揭示藤茶黄酮类化合物生物合成调控机制提供理论参考.[方法]分别采集藤茶的幼叶和成熟叶,提取其总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序列.[结论]藤茶在细胞过程、代谢过程、单有机体过程、细胞和细胞部分、结合和催化活性能力分布的基因较丰富,在一般功能、翻译、翻译后修饰、蛋白质翻转及分子伴侣的基因表达量较高,具有较强的碳水化合物代谢能力.多种关键酶基因参与藤茶黄酮类化合物的生物合成,推测其生物合成途径存在多条分支,调控机制也较复杂.     

2种不同生长规格墨瑞鳕肌肉组织转录组测序分析

【目的】掌握不同生长规格墨瑞鳕个体间差异表达基因的表达特点,为其功能相关基因深度挖掘及分子遗传育种提供科学依据。【方法】挑选同一养殖条件下极大个体和极小个体的墨瑞鳕,构建肌肉组织cDNA文库后,采用Illumina HiSeqTM 4000测序平台对存在生长差异的墨瑞鳕肌肉组织进行转录组测序分析,获得的Unigenes在Nr、Nt、Pfam、KOG/COG、Swiss-Prot、KEGG和GO等数据库中进行比对;通过FPKM及DEGseq筛选出差异表达基因,以GOseq和KOBAS对差异表达基因分别进行GO功能注释及KEGG信号通路富集分析,并采用MISA进行SSR鉴定分析。【结果】从墨瑞鳕肌肉组织中共测序获得39749条Unigenes,其长度范围在301~55230 bp,平均长度为1705 bp。注释到Nt、Nr、Swiss-Prot、Pfam数据库的Unigenes分别有27046、20824、18268和17772条,在7个数据库中均得到注释的Unigenes共计6742条,占Unigenes总数的16.96%。根据差异表达基因筛选条件P<0.05且|log₂ Fold Change|>1,共筛选出722个差异表达基因,其中上调基因308个、下调基因414个。差异表达基因GO功能注释分析结果表明,注释基因数目较多的GO功能条目包括细胞过程、代谢过程、膜、细胞器及结合等;KEGG信号通路富集分析发现,差异表达基因被成功富集到234条信号通路上,主要涉及磷脂酰肌醇3激酶/蛋白激酶信号通路、MAPK信号通路、胰岛素信号通路及FoxO信号通路等。在39749条Unigenes中鉴定筛选出22120个SSRs,占Unigene总数的55.65%,SSR的平均间距为3063 bp。【结论】基于转录组测序分析获得的墨瑞鳕肌肉组织差异表达基因以发挥结合、细胞过程及代谢过程等功能为主,且主要富集在PI3K-Akt信号通路、核糖体信号通路、FoxO信号通路及细胞凋亡等能量代谢相关通路上,通过共同协调而对墨瑞鳕的生长发育起调控作用。     

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

为进一步丰富鲫鱼基因组数据，对肝脏、血液和卵巢组织进行高通量转录组测序。采用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位点。结果可为进一步克隆和挖掘鲫鱼功能基因、多态性检测及群体遗传多样性分析等方面研究奠定基础。     

大竹蛏高通量转录组测序数据组装和分析

为进一步开发大竹蛏Solen grandis的基因资源,采用2代Illumina Hi-seq测序技术对大竹蛏的鳃组织进行了转录组测序,构建了转录组数据库,获得338 483 476条Clean Reads数据;拼接组装后获得190 856条Unigene数据,平均长度为1147 bp;与NR、NT、KO、Swiss Prot、PFAM、GO、KOG等数据库进行Blast信息比对(E-value为10-5),共获得63 337个注释基因;与NR数据库比对发现,大竹蛏转录组基因序列与长牡蛎Crassostrea gigas具有较高的同源性,为53.3%;将大竹蛏转录组的Unigene的功能通过与KOG数据库进行注释比对划分为25类;GO数据库注释可分为三类,即细胞组分、生物过程和分子功能,共包括65个分支;KEGG分析发现,大竹蛏转录组数据中按照代谢通路可分为92类,利用Blast蛋白库比对和Estscan软件进行ORF预测,获得长度大于300 nt的ORF共50 681个;通过SSR分析,共获得73 089个SSR标记。本研究中获得的转录组信息可为今后进行大竹蛏分子标记的开发和关键基因的克隆及功能分析等研究提供基础数据。     

柴胡转录组SSR的分布及序列特征分析

【目的】分析柴胡转录组中SSR的分布及序列特征,为开发多态性良好的、功能基因相关的SSR标记提供理论依据。【方法】以不同温度处理的柴胡种子为材料,经高通量转录组测序后,使用Trinity对测序结果进行de novo组装,并利用MISA对组装得到的Unigenes进行SSR位点搜索,最后统计分析SSR的分布及序列特征。【结果】从转录组数据组装获得244194条Unigenes,其N50值为1036 bp,平均长度791 bp,总长度193138105 bp。从Unigenes序列中共检测到50303个SSR位点,去除20405个复合型SSR位点,以29898个单一型SSR位点为后续分析对象。转录组SSR的出现频率为12.24%,平均分布距离6.46 kb,主要重复基元类型为二核苷酸,共17105个,占SSR总数的57.21%,其次为三核苷酸和单核苷酸,分别占SSR总数的20.75%和20.46%,四核苷酸~六核苷酸重复基元数量均较少。转录组SSR中,共有97种重复基元,其中二核苷酸和三核苷酸重复基元分别以AT/AT和ATC/GAT为主,分别占SSR总数的33.33%和3.98%;重复次数5~10次的SSR位点数量最多,共27942个,占SSR总数的93.46%。转录组SSR序列长度存在明显差异（12~76 bp）,平均长度15.28 bp。【结论】柴胡转录组的SSR位点出现频率较高,类型较丰富,具有开发出高多态性SSR分子标记的潜力,将其用于柴胡的遗传多样性分析、种质资源评价及分子标记辅助育种等研究。     

基于转录组测序的棘胸蛙SSR和SNP分子标记开发

【目的】基于转录组开发SSR和SNP分子标记用于评价棘胸蛙（Quasipaa spinosa）遗传多样性,为其种质资源的创新利用提供理论支撑。【方法】采用TRIzol试剂盒提取棘胸蛙肝脏、肌肉和肾脏组织总RNA,构建cDNA文库后利用Illumina HiSeq 2500测序平台进行高通量测序,通过MISA对棘胸蛙转录组测序数据进行SSR检索,并以SAMtools和VarScan v.2.2.7进行SNP查找。【结果】棘胸蛙转录组测序共获得93887条非冗余基因（Unigenes）,序列总长度达91352712 bp,且所有转录组Q30均超过95.00%。在93887条Unigenes中发现33019个SSRs,其中21966条Unigenes含有SSR,6688条Unigenes含有超过1个SSR;以单核苷酸重复型SSR数最多,达25788个,且出现频率最高（27.47%）。SSR的平均长度以四核苷酸重复型最长,达35.47 bp;棘胸蛙SSR以（A/T）n为绝对优势重复基元,占总SSR的65.51%,然后依次为（C/G）_n、（AT/AT）_n、（AC/GT）n、（AG/CT）_n、（AAT/ATT）_n、（AGG/CCT）_n,分别占总SSR的12.59%、5.66%、5.55%、3.31%、1.55%和1.30%。在33019个SSRs中,核苷酸重复次数主要集中在5~25次,占总SSR的99.91%,且大部分SSR位于非编码区,仅有1633个SSRs位于编码区;长度≥12 bp的SSR共计17244个,占总SSR的58.53%。挑选120对SSR引物进行引物有效性验证,发现有57对引物扩增出单一条带,且条带大小与预期结果一致。对棘胸蛙转录组序列进行SNP检索,共发现87634个SNPs,其中,56300个SNPs属于转换位点、31334个SNPs属于颠换位点,转化/颠换比达1.80,碱基的转换频率明显高于颠换频率。【结论】利用高通量转录组测序开发棘胸蛙SSR和SNP分子标记是一种切实可行的方法,能开发出通用性较高、数量较多、覆盖性较广的分子标记。棘胸蛙具有中度偏高的遗传多样性,可作为种质材料进一步开发利用。     

基于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个)。     

茶树EST SSRs特点及应用

利用NCBI公共数据库下载获得的3 306条茶树EST,去除其中低质量和冗余的序列,得到全长为635.46 kb的1 335条无冗余EST,在这些序列中共发现216个SSR,分布于185个EST中,出现的频率为16.18%,平均每2.94 kb的表达序列中有1个SSR。茶树的EST-SSR类型丰富,其中二核苷酸重复出现的频率最高,占主导地位。利用Primer 5.0软件从185条含有SSR的EST中设计引物100对,从中随机选取60对引物检测EST-SSR在云南茶树资源中的多态性。结果表明,有30对引物在供试的63份茶树资源中得到有效扩增,并具有多态性;共检测到等位基因数71个,每对引物检测出的等位基因数2~8个,平均4.9个;多态信息量(PIC)变幅为0~0.731,平均为0.499;观察杂合度变幅为0~0.970,平均0.392。EST-SSR标记分析表明云南茶树资源具有较高的遗传多样性。