基于转录组测序的波纹唇鱼SSR和SNP多态特征分析

波纹唇鱼是一种极度濒危的珍贵观赏鱼和高价值食用鱼,全面了解其转录组中SSR和SNP位点的分布及序列特征,可有助于开展波纹唇鱼遗传资源的保存和合理开发利用。利用二代高通量RNA-seq测序技术对波纹唇鱼进行转录组测序,通过MISA和Samtools对所得Unigene进行SSR与SNP位点的发掘与分析。结果显示,在150 218条Unigene序列中共发现22 428个SSR,出现频率为14.93%,平均约5.35 kb出现1个SSR。波纹唇鱼SSR的主要重复单元类型为单碱基和二碱基重复,分别占SSR总数的61.32%和19.12%,除复合类型以外,所有重复基元共65种,其中(A/T)n所占比例最高(55.64%),(AG/CT)n和(AC/GT)n分别占8.31%和6.80%。22 427个SSR处于CDS中,其中1 773个位于编码区。SSR重复次数集中在5～15次,序列平均长度为13.3 bp。3 438个SSR位点共获得669对候选引物。随机选取50对验证,发现21对可扩增出与预期产物长度大小一致的特异性条带,其中5对在4个波纹唇鱼个体间具有多态性。在Unigene中共发现245 373个SNP位点,发生频率为1/490 bp,其中转换类型发生频率显著高于颠换类型,转换类型中A/G (33.61%)和C/T (32.51%)发生频率最高,颠换类型中则是A/T的最高(11.12%)。波纹唇鱼转录组中SSR和SNP位点非常丰富,可为波纹唇鱼遗传多样性分析、亲缘关系鉴定与遗传资源开发利用等方面提供丰富的基础数据信息。     

重要食药同源植物余甘子转录组微卫星特征分析

为全面了解余甘子转录组SSR位点的分布特征和变异规律,本研究利用Illumina Hiseq 4000平台对余甘子叶片转录组进行测序,通过MISA软件对获得的Unigenes进行SSR位点搜索和统计分析。结果发现9 538条包含SSR位点的Unigenes,共检测到9 991个SSR位点,平均每5.49 kB出现1个SSR。单碱基和二碱基为余甘子转录组SSR主要重复类型,分别占SSR总数的42.3%和30.79%。位于基因编码区的SSR位点共有1 731个,出现频率为0.039 SSRs/kB,优势重复类型为三碱基重复。余甘子转录组SSR中共有169种重复基元,其中所占比例最高的是A/T（42.10%）,其次是AG/CT（22.91%）和AAG/CTT（5.02%）。SSR各基元的重复次数波动于4~75次,且多数集中于4~20次。重复片段长度≥ 20 bp的SSR占21.20%,且SSR发生频率与片段长度呈显著负相关（P<0.01）,相关系数为-0.561。本研究获得的余甘子转录组SSR位点出现频率较高、分布密度较大、低级重复基元较多,重复次数较高、长片段较多,大多数SSR位点的多态性潜能较高,用于余甘子遗传多样性分析的潜力较大,为下一步余甘子转录组SSR标记的大规模开发和群体遗传学研究提供了重要的数据信息,进而为余甘子野生资源的保护和合理开发利用提供了参考依据。     

山地虎耳草和棒腺虎耳草转录组SSR和SNP分析

利用Illumina HiSeqTM2000对山地虎耳草和棒腺虎耳草进行转录组测序,分析和比较其SSR和SNP特征。结果表明:山地虎耳草63 763条Unigene序列中含有4 622个SSR,发生频率为7.25%,有110种重复基元,平均每10.00kb出现一个SSR位点;棒腺虎耳草60 972条Unigene序列中含有4 542个SSR,发生频率为7.45%,有85种重复基元,平均每10.40kb出现一个SSR位点,略低于山地虎耳草。山地虎耳草和棒腺虎耳草转录组序列的SSR优势基元均为三核苷酸重复。2个物种的转录组SSR以5~10次的较低重复次数为主,长度主要集中于12~30bp,具有较高的多态性。山地虎耳草和棒腺虎耳草中分别获得118 424个和112 006个SNP位点,编码区的SNP位点分别占30.40%和28.59%,且在编码SNP中同义突变所占比例(30.27%、28.48%)远高于非同义突变(0.13%、0.11%)。比较发现,2个物种的各项检索结果基本一致,推测与选取的组织部位、组织的发育阶段以及物种的亲缘关系有关。     

沟眶象转录组微卫星特征分析

本研究基于高通量测序获得的沟眶象Eucryptorrhynchus chinensis(Olivier)转录组数据,采用MISA(Micro Satellite)软件对其简单序列重复(simple sequence repeat,SSR)位点进行了高通量发掘。对筛选得到的1 kb以上的Unigene(17590条,占Unigene总数的26.98%)进行了SSR分析,结果发现含SSR的序列有1665条,共发现1823个SSR,平均32.52 kb出现一个SSR。沟眶象微卫星主要以三碱基重复类型为主(637,34.94%),其次为单碱基重复(576,31.60%)。共发现104种碱基重复基元,所占比例最高的为(A/T)n(30.39%),其次为(AT/AT)n(6.91%)。沟眶象SSR的平均长度为15.74 bp,长度大于20 bp的SSR仅占总数的12.70%。研究还发现沟眶象微卫星出现的频率和其长度呈极显著负相关(P0.01),相关系数为-0.572。研究结果为开发高多态性微卫星引物来进行沟眶象功能基因组学、种群遗传结构、种群遗传多样性等研究奠定了基础。     

中间锦鸡儿转录组EST-SSR标记系统性识别与引物筛选

旨在对中间锦鸡儿转录组数据库EST信息进行SSR系统性识别和初步验证,为进一步SSR分子标记开发提供依据。对Hi Seq2000测序技术获得的中间锦鸡儿转录组Unigenes进行SSR位点搜索,共获得45 706个SSR位点,出现频率为10.38%,平均4.30kb出现一个SSR位点。SSR重复类型以单核苷酸重复序列基元为主,所占比例为56.47%;二核苷酸、三核苷酸重复序列基元的数量所占比例分别是20.56%和21.04%;其他数量的基元所占比例仅为1.9%。多核苷酸重复类型中最多的为2核苷酸重复AG/CT;其次为3核苷酸重复AAG/CTT。针对EST-SSR位点随机挑选了150对引物,通过琼脂糖凝胶电泳进行PCR验证,其中有79对能获得扩增条带,21对引物扩增出单一条带,比例为14.0%。     

基于转录组数据的密花香薷SSR位点特征分析

利用MISA软件对密花香薷转录组42 362条Unigene进行SSR位点搜索,并对其SSR序列结构及分布特征进行了分析。结果表明:(1)密花香薷转录组Unigene序列中共检测到17 564个SSR重复序列,分布于11 903条Unigene上,出现频率为28.10%,平均每3 200 bp出现一个SSR位点。(2)单、二、三核苷酸重复类型为密花香薷转录组SSR位点的主导基序类型,占总SSR位点的97.27%,3种主导基序类型中,单核苷酸所形成基元类型数量最多,共检测到169个基元类型(51.22%),单核苷酸(A/T)n基元类型占明显优势,二核苷酸重复类型(AG/CT)n基元类型占优,分别占总SSR位点的50.60%和12.17%。(3)单核苷酸SSR位点所包含重复次数最多(49),重复次数介于10～66,同一基序类型不同重复次数所形成的SSR位点数量差异较大,随重复次数的增加,SSR位点数呈下降趋势。(4)密花香薷转录组二至六核苷酸基序SSR序列长度集中在12～30 bp区间,共包含有8 190个SSR位点,占所统计SSR位点的95.60%,1 589 (≥20 bp)个SSR序列具有极高的多态性,占所统计SSR位点的18.54%。综合出现频率、分布密度、基元重复次数和长度变异等多个研究结果发现,密花香薷转录组检索到的SSR序列表现出较高的多态性潜能,具有较大的开发价值。该研究为后续密花香薷SSR分子标记引物开发奠定了理论基础。     

基于转录组的沙棘木蠹蛾简单重复序列特征分析

为开发沙棘木蠹蛾微卫星信息,利用已获得的转录组数据,对其EST-SSR位点进行发掘,进而分析其特征。结果发现含SSR的序列5126条,识别的SSR总数为7499个,SSR出现频率为51.41%。微卫星序列主要以单碱基重复为主,发生频率为39.52%。研究共发现77种碱基重复基元,所占比例最高的为(A/T)n(73.74%),其次是(AT/AT)n(3.37%)。微卫星多为重复次数为10且长度为10 bp的短序列。研究结果为沙棘木蠹蛾的SSR分子标记研究,遗传多样性分析,种群遗传结构以及关键性状基因的发掘等研究奠定基础。     

龙眼顶芽转录组简单重复序列(SSR)标记信息分析及分子标记开发

随着新一代测序技术的发展,大量的转录组数据和表达序列标签(EST)成为开发简单重复序列(SSR)标记的可利用资源。本研究利用MISA软件筛选龙眼(Dimocarpus longan)顶芽转录组数据库序列,从114 445条龙眼转录组unigene序列中发现11 546个SSR位点,SSR出现频率为10.09%。其中1 975条unigene含有两个或两个以上EST-SSR位点,占所有SSR位点的比例为17.10%,SSR出现的平均距离为7.52 kb。从龙眼转录组SSR核苷酸基序类型来看,二核苷酸(52.11%)和三核苷酸(46.15%)出现频率最高,占所有核苷酸出现频率的99.26%。在龙眼转录组SSR中二核苷酸重复基元出现频率最高的是AG/CT(4 250个,占36.81%),三核苷酸重复基元出现频率最高的是AAG/CTT(1 109个,占9.61%)。对含SSR位点的9 571条unigene序列进行引物设计,共设计出了8 347对SSR位点特异引物。随机挑选合成50对EST-SSR引物,以‘石硖’、‘储良’、‘古山2号’、‘立冬本’等四份龙眼材料的基因组DNA为模板对这批引物进行PCR扩增、筛选,结果表明,其中21对引物能产生理想的PCR产物,有效扩增率为42%;16对引物扩增条带具有多态性,占有效引物的76.2%;16对多态性引物共扩增获得50个条带,其中多态性片段21个,每对引物平均产生1.31个多态性片段。     

基于转录组数据的荔枝蒂蛀虫SSR位点信息分析

荔枝蒂蛀虫Conopomorpha sinensis Bradley是专一性危害我国荔枝和龙眼的重要害虫。简单重复序列标记(Simple sequence repeat,SSR)为短串联重复序列或微卫星标记,其在荔枝蒂蛀虫偏嗜选择寄主的遗传进化机制研究和荔枝蒂蛀虫综合治理中具有重要意义。本研究基于高通量测序获得的荔枝蒂蛀虫转录组数据,利用MISA软件从68996条转录组unigenes结果中发掘出10521个SSR位点,出现频率为15.25%。荔枝蒂蛀虫转录组中SSR的主要重复类型为单碱基重复,占SSR总数的66.22%。其次是三碱基重复,占SSR总数的24.94%。在发现的33种重复基元中共筛选获得8种优势重复基元,其中A/T在单碱基重复基元中所占的比例达98.55%。基于筛选的SSR设计的9对引物中,有4对引物得到扩增预期大小的条带。荔枝蒂蛀虫SSR位点的信息分析将为探究荔枝蒂蛀虫种群遗传结构、遗传多样性和进化关系、害虫综合治理等研究提供重要科学依据。     

基于高通量测序的怒江红山茶微卫星位点特征分析

以怒江红山茶叶片为材料,采用Illumina Hiseq 2000平台测序,共获得140 996条无冗余的序列,进行SSR位点搜索后,得到32 696个SSR位点,出现频率为23.2%。所搜索的SSR以二核苷酸重复类型最多,三核苷酸和单核苷酸次之,四、五、六核苷酸重复类型较少（<1%）。单核苷酸重复类型中以A/T基元较丰富（10.92%）;二核苷酸中AG/CT基元出现频率最大,达到49.72%,AT/AT基元和AC/GT基元所占比例相差不多,而CG/CG基元所占比例最少,为0.07%;三核苷酸重复类型中AAG/CTT最多,ACC/GGT、ATC/ATG和AGG/CCT基元次之,CCG/GGC、ACT/AGT和ACG/CGT基元较低,都小于1%;四、五、六核苷酸类型中各重复基元均较少。在怒江红山茶转录组中,微卫星的数量随着对应的重复类型、重复次数的增加而降低,也随重复区段碱基长度的增加而降低。     

基于转录组数据的黑腹胃蝇微卫星位点信息分析

黑腹胃蝇是中国新疆荒漠地区马科动物马胃蝇蛆病最主要的病原体。本研究基于已获得的黑腹胃蝇的转录组数据,利用MISA软件(1. 0版)对黑腹胃蝇转录组中1 kb以上的Unigene进行SSR位点分析。在黑腹胃蝇转录组25 847条Unigenes中筛选得到SSR总数为12 187个,存在于8 037条Unigenes当中,出现频率为31. 09%。微卫星的序列中单碱基(38. 39%)和三碱基(33. 45%)为优势重复类型。研究共发现103种重复基元,出现最多的为A/T (37. 52%),其次是AC/GT (13. 10%)。由三碱基构成的SSR基元占比较多,分别为15 bp (20. 94%),12 bp (14. 48%)和18 bp (14. 01%)。本研究是黑腹胃蝇开发EST-SSR的基础研究,为黑腹胃蝇的遗传多样性分析、变异水平分析和功能基因发掘奠定基础。     

基于转录组数据的桔小实蝇微卫星位点信息分析

以桔小实蝇为材料,从其转录组数据库中筛选功能微卫星(EST-SSR)序列,并进行SSR位点的信息分析.共获得1890个EST-SSR位点,可用于引物设计的SSR数为1296个.EST-SSR平均分布频率为1/10.21 kb,但这种分布频率在不同重复类型SSR之间相差很大.其中,三碱基重复SSR在该种昆虫的EST-SSR中出现的频率最高,结合其他文献推断三碱基重复可能是所有昆虫EST-SSR的优势类型.本文共设计了42对桔小实蝇EST-SSR引物,其中有18对引物可以扩增得到预期大小的条带.最后,探讨了基于转录组数据发掘昆虫SSR的前景和挑战,以及进行转录组EST-SSR筛选时应注意的问题.
     

扇贝“渤海红”与墨西哥湾扇贝转录组SSR信息分析及unigene功能注释

墨西哥湾扇贝是中国南方海域养殖的重要经济贝类,由于连续累代养殖,种质退化严重,扇贝“渤海红”为山东省近几年培育的国家级新品种,性状优势明显。课题组从北方引进扇贝“渤海红”拟与南方墨西哥湾扇贝进行育成杂交,创新墨西哥湾扇贝种质资源。本研究利用MISA软件对两种扇贝转录组测序获得的unigene序列进行检索,开发SSR标记,同时对SSR位点的生物学信息进行分析,对含SSR的unigene进行功能和KEGG代谢通路注释,并利用primer5软件设计SSR引物。结果显示:扇贝“渤海红”转录组中平均8.27 kb出现1个SSR,出现频率为0.11,检索到的11 512个SSR分布于9 681条unigene上,涉及到137种重复基元,二核苷酸SSR数最多(5 469),占比为46.66%,其次是单核苷酸(4 039),占比为34.46%。墨西哥湾扇贝转录组中,共发现12 241个SSR位点,分布于10 362条unigene序列上,平均8.78 kb出现1个SSR,出现频率为0.10,共包含140种重复基元,SSR重复类型以二核苷酸为主,占SSR总数的45.61%。两种扇贝转录组SSR的重复次数主要集中在5～10次。A/T、AT和ATC/ATG分别是单核苷酸、二核苷酸和三核苷酸的优势重复基元。两种扇贝转录组中共筛选出33个含SSR的生长相关unigene,且被注释到KEGG代谢通路的含特异SSR的unigene全部富集在新陈代谢和遗传信息处理两大一级通路中。两种扇贝SSR标记的开发,为杂交双亲遗传差异分析、杂交子代鉴定和杂种优势预测等奠定了分子辅助育种基础。     

星天牛转录组SSR位点特征分析

【目的】为了获得星天牛Anoplophora chinensis的SSR位点信息并开发其SSR分子标记技术,进一步为其遗传多样性以及综合治理提供理论依据。【方法】利用MISA软件,对星天牛转录组数据进行简单重复序列(SSR)位点筛选与分析;使用Primer3软件设计引物,采用PCR扩增以及电泳检测,筛选SSR引物,开发星天牛SSR分子标记技术。【结果】在9 325条unigene序列中共挖掘到2 360个SSR位点,出现频率为25.31%,涉及SSR位点序列1 758条,发生频率为18.85%。星天牛转录组中SSR的主要重复类型为单碱基重复,其次是三碱基重复,分别占总数的79.03%、12.54%。在核苷酸重复类型中,A/T基元种类数目最多,所占比例高达99.30%。SSR长度为10-11 bp的占比最高,为56.10%;重复次数为10次的数量最多,SSR位点数为1 188(50.34%)。重复次数和长度的分析结果对SSR位点的多态性获得了初步验证。在随机挑选序列设计的60对引物中,53对扩增产物达到预期大小,候选引物可用率高达88%,可在今后的研究中利用。【结论】本文对星天牛SSR位点的信息分析以及引物的设计与验证将有助于星天牛基因挖掘、种群遗传结构、遗传多样性、进化关系和综合治理的研究。     

基于转录组测序的缺须盆唇鱼微卫星标记特征分析

缺须盆唇鱼(Placocheilus cryptonemus)是云南省怒江流域特有的一种高原经济鱼类.通过在Illumi-naHiSeq 2500平台上对缺须盆唇鱼进行转录组测序(RNA-seq)分析,利用MicroSAtellite (MISA)微卫星识别工具对1～6个碱基重复、碱基数在10 bp以上的微卫星进行鉴定.结果 表明,在缺须盆唇鱼转录组中共发现分布在36764条Unigene上的30500个SSR位点,包含序列的SSR数量为18078,发生频率为49.17％.获得的单纯SSR中共有215种重复基元,其中单核苷酸重复基元类型数量最多,为88个,六碱基重复的SSR最少,其他的重复从二、三、四、五、六碱基开始,SSR的形成呈现递减规律.SSR位点的重复次数主要发生在6～10次和11～15次重复之间,发生频率分别为38.59％和36.02％.本研究对缺须盆唇鱼的转录组数据进行测序鉴定,得到了数量较多的SSR标记,同时对得到的缺须盆唇鱼SSR数据特点进行归纳总结.本研究不仅提供了一个有价值的转录组资源,而且建立了一套SSR标记,可以为缺须盆唇鱼的基础研究、应用研究提供有效的分子标记.     

陆地棉EST长度多态性与其SSR分布特征相关性分析

目的:分析陆地棉EST长度多态性与其SSR分布特征的相关性。方法:从NCBI公共数据库下载陆地棉EST序列,应用SSRIT搜索SSR,分析20 000条无冗余的EST序列。结果:在剔除低质量和冗余的序列后,得到全长为7 363.878kb的无冗余EST序列7 322条,其中含有SSR位点的EST序列数520条,占被分析EST比例的2.60%。长度在400bp以下的EST序列含SSR的比例为1.46%;长度在400bp以上的EST序列含SSR的比例为8.94%。在1～6bp的重复基元中,二核苷酸重复基元的SSR重复频率最高,占总数的63.46%,其次是三核苷酸,占总数的34.04%。二核苷酸类型(AG)n、(AT)n和三核苷酸类型(AAG)n、(ACC)n、(ACT)n、(AAT)n是SSR的主要重复基元。结论:棉花EST-SSR可用于棉花分子标记,为有针对性设计陆地棉EST-SSR引物奠定基础。     

基于陆地棉茎尖转录组序列genic-SSR标记的开发与评价

利用软件MISA和SSR Locator对陆地棉(Gossypium hirsutum L.)Coker312茎尖转录组测序得到的73515条Unigene序列进行分析,查找到4507个SSR位点,分布于4039条转录本序列中,SSR位点出现频率为5.5%,非编码区的SSR位点要显著高于编码区(2853/1654)。SSR重复基元中,三核苷酸重复基元占主导地位(51.03%),且AAG/CTT(20.08%)类型最多;其次是二核苷酸重复基元(28.76%),主要为AG/CT(55.47%)。利用引物批量设计软件共开发1569对SSR引物,在陆地棉TM-1、海岛棉3-79(G.barbadense L.)、阿非利加棉(G.herbaceum L.)、雷蒙德氏棉(G.raimondii Ulbrich.)4个棉种中进行初步评价,其中1117对引物能在4个棉种间扩增出稳定的条带,扩增率为71.2%。选择650对能稳定扩增出条带的引物,在涉及5个棉种(增加了亚洲棉(G.arboreum L.))的13份材料中进一步筛选,83对引物能在13份材料中扩增出特异性条带,PIC变幅在0.121~0.648之间,平均值为0.422,其中80对引物能在7份陆地棉材料中扩增出特异性条带,平均PIC值为0.336;挑选5对在本实验室作图亲本鲁棉研22和鲁原343间有明显多态的引物进行连锁分析,其中4对成功连锁到本实验室构建的遗传图谱上。本研究丰富了棉属genic-SSR的数量,为棉属遗传作图、性状关联分析提供了更多引物选择。     

基于转录组测序的茄子SSR标记开发

应用Trinity软件对茄子转录组测序数据组装,得到总长度为47919660 bp的45404条Unigenes,MISA从中检索到8316个SSR位点,发生频率为18.32%,平均5.63 kb一个位点。SSR位点中单碱基重复类型最多,为5372个,占到64.60%;其次为三碱基重复1628个,占到19.58%。三碱基重复中AAG/CTT是优势重复单元,占三碱基重复数的31.6%;二碱基重复中AG/CT是优势重复单元,占二碱基重复数的42.3%。利用Primer 3设计引物,共得到858对SSR引物,随机选取100对引物对17份茄子材料进行扩增,结果表明:有84对可以扩增出条带清晰的片段,有47对引物扩增片段为多态性片段。对47对多态性引物进行分析,多态性信息含量范围为0.10~0.64,平均多态性信息含量为0.32,UPGMA聚类分析可将17份材料分为3类。以上结果表明,基于茄子转录组测序开发的SSR标记可以为茄子的遗传多样性分析和遗传图谱构建提供更加丰富的标记来源。     

山地虎耳草转录组SSR信息分析

利用MISA（MicroSatellite）软件对山地虎耳草转录组拼接序列进行微卫星位点信息分析,为后期SSR标记的开发和物种遗传多样性检测提供候选序列。结果发现,在拼接得到的63 763条Unigene序列中含有4 622个SSR,发生频率为7.25%,有110种重复基元,平均每10.00 kB出现一个SSR位点。山地虎耳草转录组序列的SSR主要集中在三核苷酸重复（55.50%）,其次为二核苷酸重复（30.23%）。二核苷酸重复和三核苷酸重复中的优势重复基元分别为AG/TC和AAG/TTC。二核苷酸重复基元的重复次数类型最多,跨度最大,具有更高的多态性,三核苷酸次之,而四、五、六核苷酸重复类型很少。山地虎耳草转录组SSR以5~9次重复为主,且SSR数量随着重复次数的增加逐渐减少,基序长度主要集中于12~30 bp,多态性均在中等以上。     

两种黄连转录组中SSR位点信息分析与多态性引物开发

本研究通过对两种黄连,即黄连(Coptis chinensis Franch.)和云南黄连(C.teeta Wall.)转录组中SSR位点信息的分析,设计SSR引物,为开发新的SSR标记提供理论基础。利用MISA工具筛选黄连及云南黄连转录组测序获得的55 903和49 741条Unigenes,对其SSR位点信息进行了分析;Primer3设计SSR引物,随机选取60对引物对两种黄连,每种20株,共40株进行多态性扩增分析。黄连和云南黄连转录组分别搜索到4 071和4 041个SSR位点。两种黄连中二核苷酸和三核苷酸重复皆是主要的类型,分别占总SSR的87.71%和87.58%。黄连中,二核苷酸重复基元出现最多的为AG/CT,共953个,占总SSR的23.40%,三核苷酸重复基元出现最多的为AAG/CTT,共746个,占总SSR的18.3%;在云黄连中,AG/CT和AAG/CTT分别为934(23.10%)和773(19.10%)。随机选取60对引物进行PCR扩增,其中12对(20.00%)表现出多态性差异。利用Populations软件作图,将40个材料分为2类。两种黄连转录组中SSR位点出现频率高,类型丰富;大量的SSR为其遗传多样性分析和遗传图谱构建提供了丰富的候选分子标记。