利用转录组数据开发意大利蜜蜂的SSR分子标记

意大利蜜蜂(简称意蜂)是自然界重要的授粉昆虫,广泛用于世界各国的养蜂生产,具有很高的生态和经济价值。为了开发意蜂的SSR分子标记,利用MISA软件对基于前期获得的意蜂幼虫肠道转录组数据组装得到的12 115条unigenes进行搜索,共预测出分布于2 149条unigenes的6 312个SSR位点,其中主要的重复类型为二核苷酸重复(54.42%)和三核苷酸重复(32.49%),主要基元为AT/AT(31.2%)和AG/CT(18.6%)。进一步利用软件设计出18 444对特异性引物,随机选取24对引物对国内3个不同来源的意蜂幼虫样品进行SSR位点扩增,有22对成功扩增出目的片段。研究开发出的SSR分子标记可用于意蜂的种群遗传和分子进化等研究,结果表明利用转录组数据鉴定非模式生物SSR位点的方法可行且高效。     

大蒜转录组简单重复序列标记分析与分子标记开发

为探明大蒜转录组简单重复序列标记(SSR)的特征,开发适合大蒜育种的SSR引物,利用MISA软件对大蒜转录组序列进行SSR位点检测与信息分析,并通过PCR扩增检测引物的有效性和多态性。结果表明:大蒜转录组测序获得444 865条unigenes,共鉴定出141 132个SSR位点,发生频率为23.02%,平均每3.45 kb出现一个SSR位点。单核苷酸和二核苷酸为SSR位点中优势类型,分别占总SSR的68.02%和24.12%;SSR位点共有82种重复基序,以AT和ATAT数量较多,占总SSR位点的65.02%和12.37%。利用Primer 3.0共设计出125 616对SSR引物,在随机选取的14对引物中有12对引物能够有效扩增,其中6对引物在35份大蒜资源中表现出多态性,扩增共得到26条多态性条带,每对引物平均产生4.33个条带。UPGMA分析显示,在遗传相似系数0.756 9处,35份大蒜资源可被分成5大类群。综上所述,利用大蒜转录组数据进行SSR分子标记开发能够获得较高频率的SSR位点,且开发的SSR标记可用性强。     

大蒜转录组简单重复序列标记分析与分子标记开发

为探明大蒜转录组简单重复序列标记(SSR)的特征,开发适合大蒜育种的SSR引物,利用MISA软件对大蒜转录组序列进行SSR位点检测与信息分析,并通过PCR扩增检测引物的有效性和多态性。结果表明:大蒜转录组测序获得444 865条unigenes,共鉴定出141 132个SSR位点,出现频率为31.72%,平均每3.45 kb出现1个SSR位点。单核苷酸和二核苷酸为SSR位点中优势类型,分别占总SSR的68.02%和24.12%;SSR位点共有82种重复基序,以A/T和AT/AT数量较多,占总SSR位点的65.02%和12.37%。利用Primer 3.0共设计出125 616对SSR引物,在随机选取的14对引物中有12对引物能够有效扩增,其中6对引物在35份大蒜资源中表现出多态性,扩增共得到26条多态性条带,每对引物平均产生4.33条条带。UPGMA分析显示,在遗传相似系数0.756 9处,35份大蒜资源可被分成5大类群。综上所述,利用大蒜转录组数据进行SSR分子标记开发能够获得较高频率的SSR位点,且开发的SSR标记可用性强。     

基于RNA-seq数据的小麦条锈菌SSR标记开发

为开发小麦条锈菌基因组编码区内多态性高的SSR标记,应用MISA软件从4个小麦条锈菌分离物转录组测序数据中搜索SSR位点并设计相应引物。结果表明,在小麦条锈菌转录组2 941条unigenes中发现6 083个SSR位点,SSR发生频率为14.13%,平均每9.84 kb出现1个SSR位点。获得的SSR有144种重复基序,单核苷酸和三核苷酸是主要重复类型,分别占SSR总数的42.33%和39.90%。T和AAC、ATC、CAT分别是单核苷酸和三核苷酸的优势重复基序。转录组中SSR重复次数以5~12次为主,基序长度主要集中在12~20 bp。随机合成100对SSR引物对3个小麦条锈菌分离物基因组DNA进行PCR扩增,有87对引物可以扩增出条带。总之,这些基于转录组发掘的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分子标记,可为黄唇鱼的遗传多样性分析、亲缘鉴定和遗传图谱构建等后续研究提供基础资料。     

利用芥菜转录组信息挖掘SSR标记及用于种质分析

【目的】为开发高效的芥菜Brassicajuncea.分子标记。【方法】本研究通过对芥菜转录组测序的信息搜索SSR位点并分析其分布特点,应用Primer 3.0软件设计SSR引物,随机选择50对引物扩增44份芥菜种质,检测其多态性。【结果】芥菜转录组测序共获得55 636条unigene,全部序列有48 193 376 bp;有7 834条unigene包含SSR的序列,从中鉴定出9 526个SSR位点,其中有1 371条序列包含1个以上SSR,复合SSR有572个,SSR发生频率为14.08%。优势重复基序为三核苷酸和二核苷酸,分别占总SSR的51.12%和41.91%。二核苷酸重复基元中以AG/CT为优势重复基元,占总位点的34.74%,三核苷酸重复基元以AAG/CTT为主,占总位点的18.30%。共设计出21 282对SSR引物,随机选择50对引物进行PCR扩增,其中41对扩增出清晰可重复的预期条带,17对(占34%)在44份芥菜种质中表现出多态性。应用UPGMA得到将44份供试材料分为4大类聚类图,可准确地体现了芥菜种质材料的关系。【结论】根据芥菜转录组数据能开发出类型丰富、效率较高的SSR标记,为芥菜亲缘关系分析和遗传图谱构建等提供更可靠的标记。     

银杏转录组数据中EST-SSR位点的生物信息学分析

银杏(Ginkgo biloba L.)是雌雄异株植物,其植株价值因性别不同而异。银杏转录组数据中EST-SSR位点的生物信息学分析将为银杏遗传学研究开展提供重要的理论与方法支持。首先通过高通量测序技术获得银杏大、小孢子叶球转录组数据,然后开展数据拼接组装与EST-SSR位点挖掘及相应的生物信息学分析。转录组数据处理及拼接组装后共获得108 307条unigenes,然后利用MISA软件发掘银杏转录组数据中的SSR位点,最终从8 178条unigenes中检索出9 668个SSR位点。其中,单核苷酸重复的数量最多,有5 663个;其次是二核苷酸和三核苷酸,重复数量分别为2 471、1 438个;四核苷酸至六核苷酸重复的数量相对较少,共有96个。银杏转录组EST-SSR位点共包含147种重复基元。在单核苷酸重复中,A和T是优势重复基元类型,分别有2 808、2 685个;在二核苷酸重复基元中,AT与TA数量较多,分别为469、383个,所占比例为34.48%。此外,设计得到6 809对银杏EST-SSR位点特异引物。银杏转录组EST-SSR位点的发掘将为银杏遗传图谱构建、遗传性状分析、幼年期性别鉴定方法的建立等提供有力的理论与方法支持。     

牡丹转录组SSR信息分析及其分子标记开发

利用Misa软件对牡丹转录组数据进行SSR信息分析和分子标记开发,并通过PCR扩增进行引物多态性分析。结果显示,对牡丹转录组测序获得的131 265条Unigene进行筛选,共获得44 181个SSR位点,分布于38 316条Unigene中,SSR的发生频率和平均分布距离分别为29.19%和2.23 kb。单核苷酸、二核苷酸和三核苷酸为优势重复类型,分别占SSR总数的75.19%、15.02%和6.88%。依据转录组中的SSR信息序列,随机合成60对引物在16份牡丹品种中进行多态性扩增,其中15对引物表现出多态性,各个位点的等位基因数为4～13个,平均等位基因数8.13个,观测杂合度(H_o)、期望杂合度(H_e)、香农多样性指数(I)、多态信息含量(PIC)的平均值分别是0.757 8、0.790 3、1.718 6和0.730 8。     

美国大灵芝(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信息分析及分子标记开发

分析黑果枸杞Lycium ruthenicum转录组中简单重复序列(SSR)位点信息,开发SSR分子标记。采用MISA软件在黑果枸杞转录组中共检测到73 896个SSR位点,分布于56 170条单基因簇(unigenes)中,出现频率为26.36%,平均分布距离为3.55 kb。优势重复基序为单核苷酸、二核苷酸和三核苷酸,分别占总SSR位点的74.33%, 13.30%和11.81%;共发现84种重复基序,重复次数为5～33次,其中出现频率最高的基序为A/T, AG/CT, AT/AT, C/G和AAC/GTT。针对SSR位点设计了12 674对引物,随机挑选128对引物进行多态性验证,其中74对(57.8%)得到清晰的扩增产物;用其中的28对高多态性引物对24份枸杞种质进行扩增,共检测到等位基因256个,平均为9.1个。平均主效等位基因频率、观察杂合度、期望杂合度和多态信息量分别为0.432, 0.439, 0.712和0.678。以上结果表明:基于黑果枸杞转录组测序得到的Unigene信息可作为开发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位点出现频率高、类型丰富、多态性潜能大,为土生对齿藓分子生物学机制研究等方面提供理论依据。     

基于茶树芽叶转录组序列的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分子标记用于评价棘胸蛙（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分子标记是一种切实可行的方法，能开发出通用性较高、数量较多、覆盖性较广的分子标记。棘胸蛙具有中度偏高的遗传多样性，可作为种质材料进一步开发利用。     

四球茶转录组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的分布及序列特征分析

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

中华蜜蜂幼虫肠道参考转录组的de novo组装及SSR分子标记鉴定

【目的】利用RNA seq技术对中华蜜蜂(Apis cerana cerana,简称中蜂)幼虫肠道参考转录组进行de novo组装,并进行功能及代谢通路注释,进而利用该转录组数据进行中蜂幼虫的SSR分子标记鉴定。【方法】实验室条件下饲养中蜂幼虫,将纯化的蜜蜂球囊菌(Ascosphaera apis,简称球囊菌)孢子饲喂3日龄幼虫,剖取4、5和6日龄幼虫肠道,液氮速冻。将健康幼虫肠道与感染球囊菌的幼虫肠道同时进行Illumina测序。通过对raw reads的过滤得到clean reads,利用Trinity软件组装得到unigenes。通过BLASTx(E-value10-5)比对NCBI Nr、Swiss-Prot、KOG和KEGG数据库,对unigenes进行功能和代谢通路注释。利用MISA软件对所有unigenes进行SSR搜索,并利用Primer Premier 5软件设计特异性SSR引物,通过常规PCR对来源于北京、辽宁兴城和四川成都的中蜂幼虫肠道样本进行SSR位点鉴定。【结果】中蜂幼虫肠道的RNA seq共得到35 670 000条reads,de novo组装得到43 557个unigenes,平均长度为898 nt。共有18 225个unigenes可被注释到上述公共蛋白数据库,单独注释到NCBI Nr、Swiss-Prot、KOG和KEGG数据库的unigenes数分别为3 899、443、37和10个。KOG注释结果显示,11 442条unigenes分布于25个基因家族,其中注释到RNA加工和修饰家族的基因数最多,达1 249个。9 679个unigenes的GO分类结果显示,在生物学进程分类中,注释到细胞进程的基因最多,达4 201个,在分子功能和细胞组分类中,注释到结合与细胞的基因数最多,分别为4 935和2 900个。4 517个unigenes可注释到KEGG数据库中的216个代谢通路,注释到核糖体的基因数最多,达385个。利用MISA软件,可在7 763个unigenes搜索到13 448个SSR位点,随机选取20对SSR引物对国内3个不同来源的中蜂幼虫肠道样本的SSR位点进行扩增,有6对引物可鉴定出SSR分子标记。【结论】成功组装并注释了中蜂幼虫肠道参考转录组,可为中蜂及其幼虫的分子生物学及组学研究提供重要的参考信息,也可用于补充、丰富和检验东方蜜蜂的参考基因组,基于此转录组数据开发出6个中蜂的SSR分子标记,可应用于中蜂的基因图谱构建、基因多样性分析、基因定位等研究,也说明利用转录组数据开发非模式生物SSRs的方法可行。     

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

【目的】 以6份不同地理来源的糜子资源为试验材料,基于前期转录组测序获得的1 000对SSR引物,选出200对进行多态性检测,以期构建一批可以准确评估糜子种质遗传差异的分子标记。【方法】 用Primer Premier 5.0软件设计引物,改良CTAB法提取DNA,PCR扩增DNA和聚丙烯酰胺凝胶电泳筛选引物多态性;用PowerMarker 3.25和PopGen 1.32计算遗传多样性参数。【结果】 200对引物中97对呈单态性,80对呈多态性。单碱基序列重复引物有20对,10对具多态性,其重复基元是A（50%）和T（50%）。二碱基序列重复引物有36对,15对具多态性,其碱基重复类型有7种（AG最多,TC、GC和GA次之,CA、TA和AC最少）。三碱基序列重复引物有144对,55对具多态性,其碱基重复类型有24种（GGC、GCG和GCC最多,GAA、GCT和CGC等次之,ACC、AGG、CAG、CGT、AAG、AAC、TCG、CGA、ATT、CAA和CCA最少）。就引物分辨率(Rp)值而言,1—3碱基序列重复引物分别为0.67—4.67（平均2.07）、1.33—4.33（平均2.73）和0.67—4.00（平均1.83）。基于Rp值分析SSR分布频次,发现80个标记分布在5个区间：0—1、1—2、2—3、3—4和4—5,分别包含17（21.25%）、36（45.00%）、11（13.75%）、14（17.50%）和2个（2.50%）。就$\overline{\text{Rp}}$值而言,1—3碱基序列重复引物分别为0.33—0.67（平均0.51）、0.40—0.78（平均0.59）和0.33—0.83（平均0.59）。单碱基序列重复标记共检测到22个等位变异,每个位点为2—3个（平均2.2000）,其中8个和2个位点分别具2和3个变异;二碱基序列重复标记共检测到38个等位变异,每个位点为2—3个（平均2.5333）,其中7个和8个位点分别具2和3个变异;三碱基重复标记共检测到136个等位变异,每个位点为2—3个（平均2.4727）,其中29个和26个位点分别具2和3个变异。就多态性信息含量而言,1—3碱基序列重复引物分别为0.3750—0.5355（平均0.4293）、0.2392—0.7438（平均0.4293）和0.2392—0.7438（平均0.3989）。就多样性指数而言,1—3碱基序列重复分别为0.6365—1.0776（平均0.7497）、0.5623—1.0986（平均0.8339）和0.5623—1.0889（平均0.8312）。【结论】 基于转录组测序结果,检测200对SSR引物的多态性,发现177对（88.5%）可以扩增出完整条带,其中80对具多态性,多态率为40%。     

蒙药冷蒿转录组SSR信息分析

利用MISA(Micro SAtelite)软件对测序得到的蒙药冷蒿转录组序列143 700条跨叠群(contigs)进行简单重复序列(SSR)位点的挖掘,发现3 614条序列含有3 753个SSR位点,发生频率为2.51%,共有122种重复基元,平均每18.46 kb出现1个SSR位点。冷蒿转录组序列的SSR主要集中在三核苷酸重复(56.12%),其次是二核苷酸重复(31.60%)。AC/TG、AT/TA、CA/GT、AAT/TTA和AAC/TTG是二核苷酸、三核苷酸中的优势重复基元。冷蒿转录组SSR以5~12次重复为主,基序长度主要集中于12~36 bp。冷蒿转录组共注释43 415个contigs,其中578个SSRs位于编码区,主要以三核苷酸重复为主(397,68.69%)。从分子水平和生物信息学角度介绍了蒙药冷蒿转录组SSR信息的开发利用,其出现频率高、重复类型丰富,将为冷蒿的分子标记辅助育种、遗传多样性分析、遗传图谱构建和功能基因挖掘提供了候选序列。