四球茶转录组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%。     

基于茶树芽叶转录组序列的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分布特性进行了分析比较,以期为后续的分子标记开发提供参考。     

基于银鲳RNA-seq数据中SSR标记的信息分析

[目的]开发银鲳分子标记技术。[方法]通过对银鲳(Pampus argenteus)进行高通量转录组测序(RNA-seq)获得银鲳转录组原始试验数据,经过拼接后,获得3 715 603条unigene序列。采用生物信息学分析软件MISA对所有unigene进行简单重复序列(simple sequence repeat,SSR)位点鉴定。同时,利用软件对银鲳转录组SSR的多态性进行评价。[结果]银鲳转录组水平上,共鉴定出107 007个SSR位点,分布在97 289条unigene中,发生频率为2.62%,SSR平均密度为476个/Mbp。在银鲳转录组SSRs中,单核苷酸与二核苷酸重复序列为主要重复类型,分别占总SSRs的48.14%和34.10%。银鲳转录组数据中SSR序列共包括424种重复基元类型,单核苷酸重复基元A占较高比例,占同一重复类型SSRs的49.57%,二核苷酸重复基元TG/AC和三核苷酸重复基元GAG/AAC是优势重复基元,分别占同一重复类型SSRs的42.43%和9.61%。重复序列长度在12 bp以上的SSR标记位点数占总SSR的76.95%,具有丰富的多态性。[结论]银鲳SSRs位点具有极大的可开发性,可为银鲳遗传多样性研究、遗传图谱构建及分子辅助育种提供有效工具。     

柴胡转录组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分子标记的潜力,将其用于柴胡的遗传多样性分析、种质资源评价及分子标记辅助育种等研究。     

基于银鲳RNA－seq数据中SR标记的信息分析

[目的]开发银鲳分子标记技术。[方法]通过对银鲳（Pampus argenteus）进行高通量转录组测序（RNA－seq ）获得银鲳转录组原始试验数据,经过拼接后,获得3715603条unigene序列。采用生物信息学分析软件MISA对所有unigene进行简单重复序列（ simple se-quence repeat,SSR）位点鉴定。同时,利用软件对银鲳转录组SSR的多态性进行评价。[结果]银鲳转录组水平上,共鉴定出107007个SSR位点,分布在97289条unigene中,发生频率为2．62％,SSR平均密度为476个／Mbp。在银鲳转录组SSRs中,单核苷酸与二核苷酸重复序列为主要重复类型,分别占总SSRs的48．14％和34．10％。银鲳转录组数据中SSR序列共包括424种重复基元类型,单核苷酸重复基元A占较高比例,占同一重复类型SSRs的49．57％,二核苷酸重复基元TG／AC和三核苷酸重复基元GAG／AAC是优势重复基元,分别占同一重复类型SSRs的42．43％和9．61％。重复序列长度在12 bp以上的SSR标记位点数占总SSR的76．95％,具有丰富的多态性。[结论]银鲳SSRs位点具有极大的可开发性,可为银鲳遗传多样性研究、遗传图谱构建及分子辅助育种提供有效工具。     

基于棉花转录组测序的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分子标记筛选

[目的]筛选出可用于黄唇鱼(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分子标记,可为黄唇鱼的遗传多样性分析、亲缘鉴定和遗传图谱构建等后续研究提供基础资料.     

牙鲆EST资源的SSR信息分析

应用生物信息学方法,对牙鲆EST-SSRs的分布频率及碱基重复特点进行了归纳与分析。结果表明,5 927条牙鲆非冗余EST序列通过在线搜索共检索出471个SSRs,分布于390条EST中,出现频率为7.95%,平均分布距离为7.9 kb。包括了112种重复基元,其中二核苷酸重复基元占主导地位,占总SSR的59.02%。AC是二核苷酸中的优势基元,占二核苷酸重复类型的16.91%。三核苷酸重复基元、四核苷酸重复基元和六核苷酸重复基元分布比较分散,三核苷酸重复基元中GAG数量最多,但仅占三核苷酸重复类型的7.26%。研究结果为牙鲆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引物开发提供候选序列。     

生物土壤结皮土生对齿藓转录组SSR位点分析

[目的]分析生物土壤结皮土生对齿藓转录组SSR位点。[方法]以生物土壤结皮土生对齿藓为材料,分析土生对齿藓转录组数据,开发土生对齿藓SSR分子标记技术。[结果]通过MISA软件对土生对齿藓高通量测序获得的73 084条unigenes进行SSR筛选。从搜索序列中发现4 103个符合条件的SSRs,分布在3 420条unigenes上,出现频率为5.62%。三核苷酸和二核苷酸为主要重复类型,分别占SSR总数的47.9%和26.1%,AGC/CTG(31.9%)和AG/CT(75.5%)分别是其主要重复基元。69.0%的基序长度为12~20 bp。[结论]土生对齿藓转录组SSR位点出现频率高、类型丰富、多态性潜能大,为土生对齿藓分子生物学机制研究等方面提供理论依据。     

美国大灵芝（Ganoderma oregonense）转录组SSR位点的生物信息学分析

为了解美国大灵芝（Ganoderma oregonense）转录组中的简单重复序列（simple sequence repeat,SSR）位点信息,并为开发灵芝分子标记辅助育种技术奠定基础,将高通量转录组测序获得的美国大灵芝转录组数据通过Trinity软件拼接组装,再利用MISA软件分析其SSR位点信息。随后采用Primer 3根据SSR位点设计特异性鉴别引物,并以灵芝DNA为模板,对其有效性进行验证。共获得65 064条单基因簇（unigene）,在其中2 981条unigene中检测到3 380个符合软件参数的SSR位点,发生频率为5.19%,平均分布距离为12.60 kb。美国大灵芝转录组SSR位点总长度为52 168 bp,平均长度为15 bp。单核苷酸、三核苷酸是美国大灵芝转录组SSR的主要重复类型,分别占总SSR的35.33%和38.67%。SSR基元重复次数范围为5~73次,其中5次重复的发生次数最高。SSR重复基元共有98个类型,其中C/G是最常见的重复基元,比例为17.72%,其次是A/T,比例为17.60%。针对3 380个SSR位点设计出2 324对SSR引物,随机选择15对引物对3种灵芝栽培品种进行PCR扩增,10对引物可稳定扩增出条带,其中8对在3个灵芝栽培品种中表现出多态性。说明美国大灵芝转录组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分子标记是一种切实可行的方法，能开发出通用性较高、数量较多、覆盖性较广的分子标记。棘胸蛙具有中度偏高的遗传多样性，可作为种质材料进一步开发利用。     

油松胚珠转录组微卫星特征分析

为利用微卫星分子标记技术进一步研究油松胚珠发育分子调控机制以及微卫星对油松胚珠发育的影响,本文以前期高通量测序结果为基础,用MISA软件对所获得的转录组数据进行SSR位点的发掘和分析,为后续研究提供重要信息资源和数据保障。得到的1412个SSR位点存在于1274条序列中,平均每2.15kb出现一个SSR。油松胚珠转录组微卫星中存在87种重复基元,其中(A/T)n比例最高,约占45.89%;出现频率较高的重复类型是单核苷酸(46.88%)和三核苷酸(34.99%)重复。所得微卫星多为长度<20bp的短重复序列,≥20bp仅占总数的5.67%,重复序列出现频率与长度呈负相关。所得的SSR位点大部分位于非编码区,位于编码区的仅232个,有18个跨越了编码区和非编码区;编码区中占比例最高的是三核苷酸重复序列(105个,45.26%)。      

基于RNA-seq数据大规模挖掘蜜蜂球囊菌的SSR分子标记

基于已获得的球囊菌转录组数据预测微卫星标记(SSR),并进行SSR位点的信息分析和SSR引物的挖掘.利用软件MISA对球囊菌转录组中42610条unigenes进行搜索,共预测出7968个SSRs,分布于5233条unigenes中,其中最主要的重复类型为三核苷酸重复(53.15%),其次为二核苷酸重复(32.32%)和四核苷酸重复(8.46%).二核苷酸重复中的基序主要是AG/CT(占总量的15.8%).针对所有的SSR位点,利用Primer Premier 5软件设计出6956对SSR特异性引物,随机选取20对引物对两个不同来源的球囊菌样品进行SSR位点扩增,共有6对引物成功扩增出符合预期的目的片段.研究结果表明,利用球囊菌转录组数据开发SSR引物是可行的.     

番茄枯萎病菌基因组微卫星位点分析与引物设计

通过对前期测定的番茄枯萎病菌Ｆｕｓａｒｉｕｍｏｘｙｓｐｏｒｕｍｆ．ｓｐ．ｌｙｃｏｐｅｒｓｉｃ基因组１４５０１条ｇｅｎｅ序列进行搜索,发现４６６个简单重复序列（Ｓｉｍｐｌｅｓｅｑｕｅｎｃｅｒｅｐｅａｔ,ＳＳＲ）分布在４２５条ｕｎｉｇｅｎｅｓ中,发生频率为２.９３％．在番茄枯萎病基因组的ＳＳＲ中,三核苷酸重复占主导地位,占ＳＳＲ总数的６５.８８％．优势重复单元为ＡＧＣ／ＣＴＧ和ＡＡＧ／ＣＴＴ,分别占１６.５２％和１２.６６％．基于筛选的ＳＳＲｓ,运用Ｐｒｉｍｅｒ３软件进行引物的批量设计,共设计出４５１对引物．结果说明番茄枯萎病菌基因组中含有丰富的ＳＳＲ位点．