湿地松转录组SSR分析及EST-SSR标记开发

【目的】湿地松是南方地区优先推广的优质产脂树种。但湿地松分子遗传基础薄弱,基因组序列信息匮乏,影响了湿地松基因组学的深入研究。目前,湿地松分子研究所用的SSR标记主要来自其他近缘种或利用公共数据库中有限的基因序列资源开发的SSR标记,其多态性和通用性较差。为解决这一问题, 笔者根据湿地松转录组测序数据开发EST-SSR位点,并揭示其在转录组序列中的分布类型及特征,为湿地松分子标记辅助育种奠定基础。【方法】利用MISA软件对转录组序列进行SSRs查找和分布特征分析。查找标准参数设置为: 单核苷酸重复>10次,二核苷酸重复>6次,三、四、五、六核苷酸重复>5次。根据SSR位点两端的保守区域,利用Primer3.0设计并随机挑选120对SSR引物,通过琼脂糖电泳和毛细管电泳对来自美国和吉安的113份家系个体进行遗传多样性分析,确定引物多态性。【结果】79 574条unigenes序列中搜索到3 818个SSR位点,出现频率为4.80%,平均18.27 kb出现1个SSR位点,3 373个unigenes含有SSR位点,SSR发生频率(含SSR位点的序列数/搜索序列总数)为4.24%,其中2 980条序列含1个SSR位点,含1个以上SSR位点的序列有393条。在检测到的3 818个SSR标记中,单核苷酸分布最多,其次是二核苷酸和三核苷酸,SSR数量分别占总数的63.54%、19.15%和16.27%,而四、五、六核苷酸重复类型所占比例较小,分别为0.52%、0.13%和0.31%。SSR重复单元的重复次数分布在5~22次之间,除单核苷酸重复外的1 391个SSR中,重复5次的数量最多,为498个(35.80%);重复6次和7次的次之,分别为417个(29.98%)和198个(14.23%);重复10次以上的仅有38个(2.73%)。在检测到的731个二核苷酸重复SSR中,最常出现的重复单元为AT/AT,数量为491个(12.86%),AG/CT和AC/GT两种类型的重复单元出现的次数次之,分别为156(4.09%)和81个(2.12%)。在检测到的621个三核苷酸重复中,AAT/ATT是出现频率最高的单元,共139个(3.64%),其次是AAG/CTT,共122个(3.20%)。3 818个SSR中有24.59%的位置未知,其余的SSR则分布在非编码区域(untranslated region,UTR)或者编码区(coding sequence,CDS)上,分布数量表现为3'UTR>5'UTR>CDS。参试的120对SSR引物,有92对扩增成功(76.78%),其中24对呈现多态(20%)。24对引物(13个二核苷酸重复、7个三核苷酸重复和4个四核苷酸重复)共检测出81个等位基因,每个位点的等位基因数为2~9,平均为3.38个。多态性信息含量(polymorphism information content,PIC)为0.103~0.726,平均为0.349。【结论】通过对湿地松转录组数据的挖掘,共获得 3 818个SSR位点,主要重复单元为AT/AT和AAT/ATT,可扩增出多态性位点的引物重复单元以二、三核苷酸重复为主。基于湿地松转录组序列的SSR标记开发是可行的。