花椒cpSSR标记开发及在种间、种内的通用性分析

李思巧，韦 伊，刘洪妤，张志东，张 野，王丽华，刘玉林

（1.西北农林科技大学 林学院，陕西 杨凌712100；2.四川省林业科学研究院 生物技术与良种研究所，四川 成都 610081）

简单序列重复（SSR）广泛分布于真核生物的核基因组，叶绿体基因组以及线粒体基因组中，常指一类以1～6个核苷酸为重复基元的串联重复DNA序列［1］。因其分布广泛、多态性高且具共显性，SSR标记被广泛开发和应用于种质资源鉴定与遗传多样性分析［2-3］。在植物中，由于叶绿体基因组表现为母系遗传，以半保留方式进行自我复制，几乎不发生重组，保守性远高于核基因组。因此，叶绿体基因组中的SSR标记（chloroplast SSR，cpSSR）可在更高分类水平上进行种间、种内鉴定和遗传进化分析［4］。随着高通量测序技术的发展，大量植物的叶绿体基因组被测序与公布，相应的叶绿体SSR也得到了广泛的开发与应用［5-7］。花椒Zanthoxylum bungeanum是芸香科Rutaceae花椒属Zanthoxylum植物，其果皮晒干后是中国传统的调味料和中草药。在市场上流通的花椒果皮主要来自于花椒Z.bungeanum（又称红花椒或大红袍花椒）和竹叶花椒Z.armatum（又称青花椒），而在日本主要利用日本花椒Z.piperitum（又称胡椒木或朝仓花椒）进行麻味物质和其他药用活性成分的提取分析［8］。在前人的研究中，基于核基因组内部序列差异开发的 SCAR（sequence-related amplified polymorphism）标记［9］、 ISSR（inter-SSR）标记［10］和 ESTSSR（expressed sequence tags-SSR）标记［11］等已利用到花椒属植物的遗传分析中，而cpSSR标记在花椒属植物的开发与应用还未见报道。鉴于此，本研究对花椒cpSSR位点进行筛选，分析其分布特点并利用花椒、竹叶花椒和日本花椒种质资源开发多态性标记，以期为花椒属植物的种质资源鉴定、遗传进化分析等提供可靠的标记资源。

1 材料与方法

1.1 SSR位点搜索与引物设计

基于本研究前期获得的花椒叶绿体基因组序列［12］，利用TRF（tandem repeats finder）在线软件（http://tandem.bu.edu/trf/trf.html）搜索单核苷酸重复序列，设置重复数≥10次；利用SSRHunter软件［13］搜索核苷酸重复数位为2，3，4，5的SSR位点，对应重复数设置为5，3，3，3次。合并2种软件的检测结果，统计1～5 bp核苷酸重复所在花椒叶绿体基因组中的位置：长单拷贝序列区（LSC）、短单拷贝序列区（SSC）和2个反向重复序列区（IR）。SSR引物利用Primer3软件（http://primer3.ut.ee）设计。原则为：扩增片段为150～250 bp；引物长度为18～22 bp；退火温度为57～63℃；GC质量分数为30%～55%。引物由上海生工生物工程有限公司合成。

1.2 种质资源采集与DNA提取

本研究共选用16份花椒属种质资源进行SSR多态性的筛选，其中包括来自中国四川省和重庆市的5份竹叶花椒种质资源和1份日本花椒种质资源（均采自四川省林业科学研究院生物技术与良种研究所花椒种质圃）以及来自四川、贵州、陕西、山东和山西等5个省的10份花椒种质资源（均采自西北农林科技大学凤县花椒试验示范站）进行多态性的筛选（表1）。所有种质资源选取3株，取鲜嫩叶片混合，硅胶干燥，常温保存备用。

供试种质资源全基因组DNA利用 “植物全基因组DNA提取试剂盒（DP305）”（北京天根生物有限公司）提取，质量分数为1%琼脂糖结合紫外分光光度计检测DNA质量和浓度后稀释至20 mg·L-1，4℃保存备用。

1.3 PCR扩增及数据统计分析

PCR扩增体系为20 μL，包括 10 μL 2×TaqMaster Mix（南京诺唯赞生物科技有限公司），模板DNA（20 mg·L-1）5.0 μL， 正、 反向引物（2 μmol·L-1）各 1.5 μL， 双蒸水（ddH2O）2.0 μL。 PCR 扩增程序为： 95℃预变性5 min；然后95℃变性30 s，55℃退火30 s，72℃延伸30 s，35个循环；72℃延伸8 min后4℃保存。取1.5 μL的PCR扩增产物，质量分数为8%的聚丙烯酰胺凝胶电泳（恒定功率100 W，100 min）分离检测，使用pBR322 DNA/MspI Marker作为标准分子量，银染拍照保存。

表1 供试16份花椒属种质资源来源地Table 1 Origin of 16 germplasms of Zanthoxylum

根据电泳图像，统计扩增条带并构建0～1数据矩阵，即有共迁移条带记为1，无条带则记为0。利用NTSYSpc-2.10软件计算遗传相似系数，采用非加权配对算术平均法（UPGMA）进行遗传相似性聚类。

2 结果与分析

2.1 花椒叶绿体基因组中SSR的分布特点

花椒叶绿体基因组全长为158401 bp，共搜索出了144个SSR位点（图1和表2），即在花椒叶绿体基因组中平均每隔1100.00 bp出现1个SSR位点。

根据SSR位点在花椒叶绿体基因组中的位置（图1）可知：LSC区分布85个SSR位点，占SSR总数的59.03%，平均1010.56 bp出现1个SSR位点；SSC区分布23个SSR位点，占SSR总数的15.97%，平均765.65 bp分布1个SSR位点；在 2个 IR区（IRA和IRB），均分布18个反向重复的SSR位点，合计占SSR总数的25.00%，平均1524.72 bp分布1个SSR位点。

图1 SSR在花椒叶绿体基因组中的分布特点Figure 1 Distribution of SSRs in the chloroplast genome of Zanthoxylum bungeanum

数量上，单核苷酸重复和三核苷酸重复分别有62和70个，两者占SSR总数的91.67%，且在LSC区、SSC区和IR区均有分布。二核苷酸重复和四核苷酸重复分别只有5个（3.47%）和7个（4.86%），但二核苷酸重复仅分布于SSC区，而四核苷酸重复仅分布于LSC区。

重复次数上，单核苷酸重复次数多集中在10～15次，共54个位点，占单核苷酸重复总数的87.10%，而重复数＞15次的共有8个位点；二核苷酸重复数集中在5次，三核苷酸和四核苷酸重复数集中在3次。根据重复基元类型，单核苷酸重复中A/T重复共60个，而C/G重复位点仅有2个；二核苷酸重复基元均为AT/TA；AAG/CTT和AAT/ATT为三核苷酸重复中的优势重复基元，均为25个，其余重复基元总计为20个；此外，AAAT为四核苷酸重复中的优势重复基元（3个），其余重复基元总计为4个（表2）。

2.2 多态性标记开发与通用性分析

根据SSR位点两侧序列，由Primer3软件设计了30对引物（命名为YLT-01～YLT-30）并合成。其中，YLT-01～YLT-20扩增片段包含的是2个及2个以上的核苷酸重复位点（YLT-01～YLT-12位于LSC；YLT-14～YLT-18位于SSC区；YLT-13、19、20位于IR区），YLT-21～YLT-30扩增片段包含的是单核苷酸重复位点（YLT-21～YLT-27位于LSC区；YLT-30位于SSC区；YLT-28、29位于IR区）。

表2 花椒cpSSR的重复类型与次数Table 2 Repeat types and times of cpSSR in Zanthoxylum bungeanum

经PCR扩增筛选，30对引物均能扩增出稳定的目的片段。其中，10对引物在供试种质中呈现多态性，占筛选引物总数的33.33%。如表3所示：10对引物共检测到24个等位基因。其中，YLT-26包含4个等位基因，YLT-21和YLT-29包含3个等位基因，其余的7个位点（YLT-01，YLT-17，YLT-22，YLT-23，YLT-24，YLT-25和YLT-27）均检测出2个等位基因。在包含2个等位基因的7个位点中，YLT-01，YLT-17，YLT-24，YLT-25和YLT-27能区分出日本花椒，花椒与竹叶花椒为一类（扩增片段大小一致），而YLT-22和YLT-23的扩增结果为：10份花椒种质同一带型，竹叶花椒与日本花椒为同一带型；YLT-21和YLT-29能将花椒、竹叶花椒和日本花椒区分开；YLT-26不仅能区分出花椒、竹叶花椒和日本花椒3个花椒属的不同种，且能区分出竹叶花椒种质中的狗椒。此外，未检测到有标记在10份花椒种质呈现多态性（图2）。

在设计的20对包含2个及2个以上的核苷酸重复位点的引物中，仅有2对引物表现出多态性；而包含单核苷酸重复位点的10对引物中，有8对引物呈现多态性，多态性比例远高于包含多核苷酸重复的位点。

表3 多态性SSR引物的基本信息Table 3 Characterizations of polymorphic SSR primers pairs

2.3 聚类分析

基于10个多态性标记的检测结果，利用NTSYSpc-2.10软件对供试的16份花椒属种质进行UPGMA聚类分析。如图2所示：在遗传相似系数阈值为0.22时，可将16份花椒种质聚为2类，即花椒和竹叶花椒为1类，日本花椒为1类；在遗传相似系数阈值为0.59时，可将花椒和竹叶花椒分为2个小类；而在遗传系数阈值为0.91时，可将竹叶花椒中的狗椒与其他4种竹叶花椒种质区分为2个小类。聚类结果与多态性标记的鉴定结果保持一致。

图2 16份花椒属种质的UPGMA聚类分析Figure 2 Cluster diagram for 16 germplasms of Zanthoxylum by UPGMA method

3 讨论

本研究在花椒的叶绿体基因组中共筛选到144个SSR位点，其中三核苷酸重复的数量要高于单核苷酸重复以及二核苷酸和四核苷酸重复的数量，与杨丽等［14］在麻竹Dendrocalamus latiflorus和绿竹Bambusa oldhamii叶绿体基因组中获得的SSR位点的规律一致，但与蒺藜苜蓿Medicao truncatula［15］，棉花Gossypium hirsutum［16］和霍山石斛Dendrobium huoshanense［17］等的叶绿体基因组中均以单核苷酸重复为主不同。原则上，碱基重复基元越多，越难形成SSR位点。本研究在三核苷酸重复次数的筛选上，选择重复数大于或等于3次的位点，而大于3次的重复位点仅占三核苷酸重复总数的11.43%（8个）。因此，本研究的三核苷酸重复数量多于单核苷酸重复是由软件参数设置不同引起。此外，花椒cpSSR在LSC和SSC区的分布密度要高于IR区。在花椒叶绿体基因组中，GC质量分数在LSC区，SSC区和IR区中分别为36.87%，33.51%和42.54%［12］。由于GC间有3个氢键相连，打破GC键所需的能量要高于AT间的2个氢键，因此更难造成碱基数量的增加或减少［18］。由此推断，花椒叶绿体基因组中IR区SSR位点相对稀疏的原因与GC质量分数相关。

叶绿体基因组由于保守性高，不涉及基因重组，因此，其中突变位点也更容易保留，基于叶绿体开发的SSR标记除了保留核基因组SSR标记的特点外，其可重复性也更高，在种质资源鉴定中也更具说服力［19］。在本研究获得的10个多态性引物中，有2对引物可以把本研究所选花椒属的3个种区分开，1对引物可以将竹叶花椒种内的 ‘狗椒’区分出。然而，10对多态性引物均未在10份花椒种质资源中表现出多态性。由此推断，相对于竹叶花椒，花椒种内的种质资源在叶绿体基因组进化上保守性更高。

本研究开发的花椒cpSSR可作为花椒属种间鉴定的有效标记，但要应用于种内鉴定，需进一步筛选种间多态性更高的cpSSR标记或利用EST-SSR、核基因组SSR（Genomic-SSR）等保守性相对较低的标记资源进行区分。