微卫星DNA又叫简短串联重复或简单序列重复(simple sequence repeats，SSRs)，一般是以1~6个核苷酸为重复单位的串联重复序列，在种群内不同个体之间的重复数量不同(Tóth et al., 2000)。随着PCR技术的出现，微卫星DNA的这一特性被转化为用途广泛的遗传标记。同时，微卫星DNA在病毒到真核生物的基因组中均有分布，且具有高度的长度多态性(Zane et al., 2002)，加之微卫星这一共显性分子标记有突变率快、多态性高、易于扩增等优点，被广泛应用于种群遗传、谱系地理、个体识别和亲子鉴定等领域(Selkoe & Toonen，2006)。开发微卫星标记的方法虽然很多，但步骤繁琐、耗时耗力的问题在二代测序系统推出前一直存在。2007年，罗氏454公司推出了基于焦磷酸测序法的第二代基因组测序系统——Genome Sequencer FLX System (GS FLX)。454高通量测序除具有二代测序高效、快捷的普遍特点外，还有比其他二代测序的片段读长更大的优点，因此更适合于微卫星标记的开发(Allentoft et al., 2009；Xia et al., 2018)。

南疆沙蜥Phrynocephalus forsythii隶属于鬣蜥科Agamidae沙蜥属Phrynocephalus，是我国特有的一种小型卵胎生爬行动物(头体长36~50 mm，尾长48~62 mm)，分布于新疆维吾尔自治区天山山脉以南海拔1 400~ 3 100 m的广大地区。南疆沙蜥的研究集中在生态适应(戴昆，马鸣，1991)、核型(曾晓茂等，1997)、组织学(吾玛尔·阿布力孜，杨立中，1999；吾玛尔·阿布力孜等，2000)、基于ND4基因的谱系地理格局(Zhang et al., 2010)以及线粒体基因组特征(Chen et al., 2016, 2019；Shao et al., 2016)等。微卫星标记的研究仅见于Nie等(2015)运用转录组测序对青海沙蜥P. vlangalii开展微卫星位点筛选，跨物种在南疆沙蜥中检测得到22个可扩增的微卫星位点。

本研究通过Roche 454 GS FLX高通量测序对南疆沙蜥进行低覆盖度基因组测序，并利用Krait(Du et al., 2018)首次在南疆沙蜥基因组中进行微卫星的查找与统计，旨在对其基因组水平上的微卫星重复序列的种类、数量等分布特征进行初步探索。考虑到微卫星标记多态性，三碱基和四碱基重复微卫星的不易产生由于滑链错配形成影子带(O'reilly & Wright，1995)等优越性(O'connell & Wright，1997)，因此，选取部分三、四碱基重复的微卫星进行引物设计筛选，得到可用于微卫星分析的部分候选引物，以期为利用微卫星标记研究南疆沙蜥种群遗传结构奠定基础。

用于基因组测序的南疆沙蜥标本(标本号：WGXG08351)于2008年6月采自新疆维吾尔自治区若羌县(86.25804°E，41.33690°N，海拔827 m)，标本经75%乙醇固定后，保存于中国科学院成都生物研究所两栖爬行动物标本馆。其肝脏组织样品固定于95%乙醇，-20 ℃保存。将提取的基因组DNA送至上海美吉生物信息科技有限公司进行Roche 454 GS FLX基因组测序。采用GS FLX系统对测序数据进行过滤、整理，使用Newbler 2.6(Roche，2011)对优化后的原始数据从头拼装，并将组装出的contigs和未拼装进去的single read序列合并。

采用Krait(Du et al., 2018)对合并的总序列进行完美型微卫星位点的查找，搜索标准如下：单碱基重复≥12个拷贝，二碱基重复≥7个拷贝、三碱基重复≥5个拷贝，其他碱基重复(四碱基、五碱基、六碱基)≥4个拷贝；重复序列两端的侧翼序列长100 bp；其余参数采用默认设置。同时，对部分三碱基和四碱基重复类型在Krait中基于primer3和primer3-py进行引物设计，参数为：PCR产物长100~450 bp；引物长20~27 bp；引物熔解温度55~65 ℃；GC含量30%~70%；其余采用默认参数设置。对设计出的引物再依据以下原则二次筛选：①引物序列碱基尽量随机分布，3’端不超过3个连续的G或C，避免引物在GC富集序列区的错误引发；②引物自身及引物之间尽量无互补序列，避免引物自身折叠成发卡结构使引物本身复性，引物自身无连续4个碱基或4个碱基的互补；③引物序列3’端尽量无连续的3个相同碱基相连的情况，避免引起错配。

采用Roche 454 GS FLX高通量测序对原始测序数据的每条read质量过滤后，对下机数据进行统计，得到78 927条reads，共39 807 639 bp，最短的22 bp，最长的1 132 bp，平均长度为504.4 bp，本次测序获得的reads长度主要为581~760 bp。利用Newbler对reads拼接得到670个contigs，共341 384 bp，最短的100 bp，最长的15 140 bp，平均为509.5 bp。未拼接上的reads有55 239个，共27 984 639 bp。

将组装出的contigs和未组装进去的single read序列合并，得到合并序列55 909个，共28 326 023 bp。对合并的总序列进行微卫星位点的查找，获得南疆沙蜥全基因组中6种完美型微卫星12 109个，总长度371 160 bp，占基因组的1.32%，相对丰度为427.62个/Mb。其中，数量最多的为四碱基重复类型，占重复类型序列总数的33.34%；其次是二碱基重复类型，占28.09%；其余几种类型占比均小于20%，最少的为六碱基重复类型，只占1.46%，相对丰度为6.25个/Mb(表 1)。

同一种类型的重复序列中，各重复拷贝类别所占的比例也不相同(表 2)。单碱基重复类型中，C最多，有1 206个；二碱基重复类型中，AC最多，有2 101个。AAC(392个)、AAAT(1 592个)、AAAAT(181个)和AACCCT(40个)分别是三碱基、四碱基、五碱基和六碱基重复类型中最多的重复拷贝类别。所有搜索到的完美型微卫星中，AC(17.35%)、AAAT(13.15%)、C(9.96%)、AG(9.67%)、A(3.95%)、ATAG(3.63%)、AAC(3.24%)、AGG(2.96%)、AAT(2.64%)和AAAC(2.58%)是数量最多的前10种重复拷贝类别。其他重复拷贝类别数量均不超过300个，相对丰度均不超过10个/Mb(图 1)。

图 1 南疆沙蜥基因组中相对丰度最高的微卫星基序分布 Fig. 1 Relative abundance of the most relative abundant microsatellite motifs in the genome of Phrynocephalus forsythii

图选项

单碱基重复分布范围为12~32次，主要分布在12~16次，有1 464个，占单碱基重复类型总数的86.94%；二碱基重复分布范围为7~204次，主要分布在7~11次，有1 402个，占二碱基重复类型总数的41.21%；三碱基重复分布范围为5~35次，主要分布在5~13次，有1 832个，占三碱基重复类型总数的80.81%；四碱基重复分布范围为4~121次，主要分布在4~15次，共3 616个，占四碱基重复类型总数的89.57%；五碱基重复分布范围主要为4~7次，共492个，占五碱基重复类型总数的90.77%；六碱基重复分布范围为4~10次，主要分布在4~7次，共169个，占六碱基重复类型总数的95.48%(图 2，表 3)。

图 2 南疆沙蜥不同重复拷贝类型微卫星的重复次数分布 Fig. 2 Distribution of the copy number in different microsatellite motifs for Phrynocephalus forsythii

图选项

挑选了部分三碱基和四碱基重复的微卫星进行引物的设计，对引物二次筛选后，最后保留100对候选引物(表 4)。

本研究利用Roche 454 GS FLX高通量测序平台对南疆沙蜥进行了低覆盖度基因组测序，并利用Krait(Du et al., 2018)进行微卫星的查找与统计。结果表明，不同类型的微卫星在南疆沙蜥基因组中丰度差异很大，以四碱基重复类型最多。与以往对物种基因组中进行微卫星查找的研究相比，这种情况很少见。植物基因组中，已有研究表明占优势的微卫星是二碱基重复(Tóth et al., 2000)或六碱基重复(马秋月等，2013；廖卓毅等，2014)。对动物的研究，无脊椎动物中占主导的微卫星重复类别没有特别明显的规律，最丰富的类型既有单碱基重复(汪自立等，2013)，又有二碱基重复(高焕等，2004)和三碱基重复(魏朝明等，2007；汪自立等，2013)；脊椎动物中单碱基重复占主导的较多(黄杰等, 2012, 2015；戚文华等，2013；李午佼等，2014；聂虎等，2017；崔凯，岳碧松，2018)。有研究推测，多数物种基因组微卫星中单碱基重复丰度最高，可能是由于微卫星序列越长，突变率就越高，因此稳定性就越差(Wierdl et al., 1997)；在很多真核生物中，重复长度和重复频率呈负相关(Katti et al., 2001)。然而，我们对南疆沙蜥的研究结果与之相差甚远，与南疆沙蜥同属于爬行类的红尾蚺Boa constrictor和原矛头蝮Protobothrops mucrosquamatus也并未表现出一致的微卫星重复类型(聂虎等，2017)。总体而言，不同物种基因组中的微卫星重复类型丰度表现不尽相同，重复长度和重复频率的相关关系可能也不宜一概而论。

从南疆沙蜥重复基序的拷贝数分布来看，重复次数少的类型占多数。谱系(Weber & Wong，1993；Schlötterer et al., 1998)、种群(Goldstein & Clark，1995)和克隆株(Wierdl et al., 1997)的微卫星分析结果均支持SSR位点的变异频率与拷贝数存在一定相关性，即重复次数越多，SSR产生变异的可能性越大，但Schlötterer(2000)认为这种相关性是否呈线性关系还有待进一步研究。

据Guichoux等(2011)统计，自1990年以来，有关微卫星的论文数量增长迅猛，一定程度上也促进了二代测序技术的发展。而二代测序出现之前，微卫星重复类型的获得对于非模式生物来说均是一个瓶颈(Guichoux et al., 2011)。随着高通量测序技术的发展，测序成本逐渐降低，可获得的数据越来越多，这为进一步增进对非模式物种的深度了解提供了很好的契机。Nie等(2015)运用Illumina高通量RNA-seq测序技术对青海沙蜥进行转录组测序，鉴定筛选出25个微卫星位点，其中22个位点可在南疆沙蜥中扩增。本研究使用454 GS FLX高通量测序技术对南疆沙蜥进行基因组测序，SSR分子标记开发及分析，利用软件预测首次对南疆沙蜥基因组微卫星进行了搜索统计并设计筛选了100对基于三碱基和四碱基重复的候选引物，为利用微卫星标记研究南疆沙蜥种群遗传结构奠定了基础。