中国樱桃14个自然居群遗传多样性和遗传结构的SSR评价

园艺学报，2016，43(6)：1148–1156.ActaHorticulturaeSinica1148doi：10.16420/j.issn.0513-353x.2015-0815；http：//www.ahs.ac.cn收稿日期：2016–04–06；修回日期：2016–06–07基金项目：‘十二五’国家科技 计划 项目进度计划表范例计划下载计划下载计划下载课程教学计划下载 项目（2013BAD02B00）；唐仲英果树育种基金项目（2009YZ033）*通信作者Authorforcorrespondence（E-mail：cylxlcz0673@sina.com；Tel：13186045875）中国樱桃14个自然居群遗传多样性和遗传结构的SSR评价高天翔1，蔡宇良1,*，冯瑛1，赵晓军2（1西北农林科技大学园艺学院，农业部西北园艺种质资源与遗传改良重点开放实验室，陕西杨凌712100；2陕西省铜川市果业局，陕西铜川727100）摘要：以中国樱桃14个自然居群280个个体为材料，采用SSR分子标记技术分析其遗传多样性和遗传结构。结果显示：11对SSR引物共检测到80个等位基因，各引物扩增条带在4~13条之间。基因多样性指数（h）为0.5431~0.7151，Shannon’s信息指数（I）为0.9057~1.4684。分子方差分析（AMOVA） 表 关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf 明，遗传变异主要来自居群内（70.00%），Mantel检验显示总群体的遗传距离和地理距离显著相关（r=0.472，P=0.011)。因此，中国樱桃在居群水平（PPL=100%，h=0.643，I=1.207）和物种水平（PPL=100%，h=0.743，I=1.591）上均具有较高的遗传多样性。关键词：中国樱桃；SSR；遗传多样性；遗传分化中图分类号：S662.5文献标志码：A文章编号：0513-353X（2016）06-1148-09GeneticDiversityandGeneticStructureofPrunuspseudocerasusPopulationsfromChinaasRevealedbySSRMarkersGAOTian-xiang1，CAIYu-liang1,*，FENGYing1，andZHAOXiao-jun2（1CollegeofHorticulture，NorthwestHorticulturalPlantsGeneticandBreedingKeyLaboratoryofMinistryofAgriculture，NorthwestA&FUniversity，Yangling，Shaanxi712100，China；2FruitIndustryBureauofTongchuanCity，Tongchuan，Shaanxi727100，China）Abstract：Simplesequencerepeat（SSR）markerswereemployedtoinvestigatethegeneticdiversityandgeneticstructureof280individualssampledfrom14naturalpopulationsofPrunuspseudocerasus.Atotalof80allelesof11lociweredetected，andthenumberofallelesperlocusrangedfrom4to13.Therelativelyhighlevelsofgenediversity（h，0.5431–0.7151）andShannon’sdiversity（I，0.9057–1.4684）revealedrelativelyrichgeneticdiversityamongthe14Prunuspseudocerasuspopulations.Thehierarchicalanalysisofmolecularvariance（AMOVA）revealedthegeneticdifferentiationmainlywithinpopulations（70.00%）.Manteltestrevealedasignificantcorrelationbetweengeographicandgeneticdistances（r=0.472，P=0.011）.Therefore，thegeneticdiversitywashighbothatthepopulationlevel（PPL=100%，h=0.643，I=1.207）andthespecieslevel（PPL=100%，h=0.743，I=1.591）.Keywords：Prunuspseudocerasus；SSR；geneticdiversity；geneticdifferentiation高天翔，蔡宇良，冯瑛，赵晓军.中国樱桃14个自然居群遗传多样性和遗传结构的SSR评价.园艺学报，2016，43(6)：1148–1156.1149由于长期的定向选择，中国樱桃（PrunuspseudocerasusL.）已表现出遗传基础狭窄，果粒小，肉质薄，不耐贮运等问题（黄晓姣等，2011）。开展中国樱桃居群遗传多样性与遗传结构的研究，对野生优良樱桃资源的收集利用，以及中国樱桃的遗传改良都具有重要意义。对中国樱桃的研究主要以栽培品种为主，有基于ITS序列（何文等，2014）和ISSR（宋常美等，2011）分子标记的种质资源分析，也有遗传资源多样性的研究（黄晓姣等，2013），以及基于经济性状的研究（黄晓姣等，2011）等。而对野生中国樱桃的研究相对较少，仅见基于RAPD（蔡宇良，2006）、PCR-RFLP（曹东伟等，2007）、SSR（陈娇等，2013）分子标记对野生中国樱桃遗传多样性的研究，但样本采集涉及地域有一定局限性。本研究中利用SSR（Simplesequencerepeat）分子标记，对中国樱桃14个自然居群280个个体进行了遗传多样性和遗传结构分析，为充分了解中国樱桃种质资源的遗传信息，以及中国樱桃资源的有效保存和生产利用提供参考。1材料与方法1.1植物材料2014年6—9月，根据中国数字植物标本馆、中国植物志、相关文献记载及当地群众指引，在中国10个省（市）选取14个具代表性的地点对中国樱桃资源的生存现状进行野外调查，根据自然居群现有分布规模对每个居群以株为单位随机取样，每个居群采集20个样本，样本间距离大于50m，所采植株直径不小于10cm，采集幼叶立即放入装有干燥变色硅胶的自封塑料袋中干燥备用。共计14个居群（表1），280个个体。表1野生中国樱桃居群采样点概况Table1Geographiclocalities，samplesizesofPrunuspseudocerasuspopulationsinthisstudy居群代码Populationcode采样地Location海拔/mAltitude经度Longitude纬度LatitudeXZP陕西西安灞桥区西张坡村Xizhangpo，Baqiao，Xi’an，Shaanxi639E109°08.055′N34°13.510′NZ陕西南郑县团山Tuanshan，Nanzheng，Shaanxi871E106°56.638′N32°57.732′TS山东泰安岱岳区小常庄Xiaochangzhuang，Daiyuequ，Tai’an，Shandong197E117°03.290′N36°12.916′LS山东青岛崂山区卧龙村Wolong，Laoshanqu，Qingdao，Shandong121E120°34.446′N36°14.291′BTM河南内乡县宝天曼Baotianman，Neixiang，He’nan1032E111°55.154′N33°29.863′TB河南桐柏县南山Nanshan，Tongbai，He’nan195E113°24.586′N32°20.233′WB甘肃康县王坝乡Wangba，Kangxian，Gansu1202E105°41.348′N33°20.334′YX四川汶川县映秀镇Yingxiu，Wenchuan，Sichuan1153E103°29.051′N31°03.459′JYS重庆缙云山Jinyunshan，Chongqing817E106°23.232′N29°49.959′FHS贵州遵义凤凰山Fenghuangshan，Zunyi，Guizhou948E106°55.926′N27°42.028′YJ贵州铜仁市印江栗子园Liziyuan，Yinjiang，Tongren，Guizhou732E108°36.272′N27°57.070′HLT云南昆明黑龙潭Heilongtan，Kunming，Yunnan1922E102°44.848′N25°07.910′HPS湖南石门县壶瓶山镇黄虎村Huanghu，Hupingshan，Shimen，Hunan412E110°48.107′N29°55.612′SNJ湖北神农架林区木鱼镇Muyu，Shennongjia，Hubei1174E110°24.488′N31°27.122′1.2DNA提取和SSR-PCRDNA提取采用3×CTAB法（Zhangetal.，2008），利用紫外分光光度计和1.0%琼脂糖凝胶电泳检测DNA的浓度和质量，稀释至30~60ng·μL-1，–20℃保存备用。SSR-PCR反应体系：包含0.4μL引物（0.2mmol·L-1），1μLDNA 模板 个人简介word模板免费下载关于员工迟到处罚通告模板康奈尔office模板下载康奈尔 笔记本 模板 下载软件方案模板免费下载 （30~60ng·μL-1），3.6μLddH2O，5μLPCRMIX混合液（康为公司），共10μL。PCR扩增程序为：94℃预变性4min；94℃变性30s，63~58℃退火30s，72℃延伸30s，38个循环后，72℃延伸10min，然后8℃GaoTian-xiang，CaiYu-liang，FengYing，ZhaoXiao-jun.GeneticdiversityandgeneticstructureofPrunuspseudocerasuspopulationsfromChinaasrevealedbySSRmarkers.1150ActaHorticulturaeSinica，2016，43(6)：1148–1156.保存。选择来自李属的100对SSR引物（Dirlewangeretal.，2002；Strussetal.，2003；Mnejjaetal.，2004；Aietal.，2008；Guarinoetal.，2009；Illaetal.，2009；Wünschetal.，2009），筛选出条带清晰、多态性丰富、重复性好、稳定性高的11对引物（表2），对280个个体进行分析。扩增产物（2μL）在1×TBE缓冲液180V恒定电压下，用9%聚丙烯酰胺凝胶分离1.5h，终止电泳，银染显色，照相保存。表2SSR引物信息Table2SSRprimersandtheirinformationsusedforgeneticsdiversityanalysis来源Origin引物Primer引物序列（5′–3′）Primersequence引物获得文献Referenceofprimer桃PrunuspersicaEPDCU5060F：ACCAAATTGGACATGCAACCR：CGGTCGAGAAGACTGAGGAGIllaetal.，2009BPPCT002F：TCGACAGCTTGATCTTGACCR：CAATGCCTACGGAGATAAAAGACWünschetal.，2009UDP98-022F：CTAGTTGTGCACACTCACGCR：GTCGCAGGAACAGTAAGCCTWünschetal.，2009BPPCT014F：TTGTCTGCCTCTCATCTTAACCR：CATCGCAGAGAACTGAGAGCDirlewangeretal.，2002BPPCT034F：CTACCTGAAATAAGCAGAGCCATR：CAATGGAGAATGGGGTGCDirlewangeretal.，2002甜樱桃PrunusaviumUCD-CH14F：GTACACGGACCCAATCCTGR：TCTAACATCATGTTAAACATCGStrussetal.，2003EMPaS12F：TGTGCTAATGCCAAAAATACCR：ACATGCATTTCAACCCACTCGuarinoetal.，2009BPPCT005F：GCTAGCAGGGCACTTGATCR：ACGCGTGTACGGTGGATDirlewangeretal.，2002BPPCT037F：CATGGAAGAGGATCAAGTGCR：CTTGAAGGTAGTGCCAAAGCDirlewangeretal.，2002李PrunussalicinaCPSCT029F：ATGGGCTAGAAGTGGTGGTGR：ATTCCGACTCGAAACGAAGAMnejjaetal.，2004中国樱桃PrunuspseudocerasusPTCR18F：GCTAATATCAAATCCCAGCTCTCR：TGAAGAAGTATGGCTTCTGTGGAietal.，20081.3数据分析根据分子量大小对扩增结果进行条带分析，从大到小依次记为A、B、C……（万宣伍等，2010）。利用GenAlExv6.501（Peakall&Smouse，2006）软件获得用于评价居群遗传多样性的指标，包括多态性位点比率（PPL）、观测等位基因数（Na）、有效等位基因数（Ne）、莱特固定指数（Fis）、Nei’s基因多样性指数（h）、Shannon’s信息指数（I）（Lewontin，1972），及种群间分化系数（Fst）、Nei’s标准遗传距离（GD）和基因流（Nm）（McDermott&McDonald，1993）等遗传结构指标。同时使用Structurev2.3.4软件的贝叶斯聚类法，根据每个个体的遗传背景进行分组（Pritchardetal.，2000）。设置分组数量1~15，每个K值运行20次，使用马尔可夫链法（Markov’schainMonteCarlo，MCMC），设burn-in为100000次、run-length为1000000次迭代，计算不同的K值对应的lnP（D）值，依据拟然值最大的原则确定最佳K值（Evannoetal.，2005）。利用MEGA6软件采用非加权组平均法（UPGMA）进行聚类分析，并且在Nei’s遗传距离的基础上完成居群聚类图。利用GenAlExv6.501软件的分子变异方差分析（Thehierarchicalanalysisofmolecularvariance，AMOVA）（Excoffieretal.，1992）评价基因多样性在群体间、群体内、区域间的差异和贡献。并对遗传距离和地理距离之间的相关性进行Manteltest检验，推断群体是否存在地理隔离（Smouseetal.，1986；Smouse&Long，1992）。高天翔，蔡宇良，冯瑛，赵晓军.中国樱桃14个自然居群遗传多样性和遗传结构的SSR评价.园艺学报，2016，43(6)：1148–1156.11512结果与分析2.1SSR引物扩增结果11对SSR引物中各引物扩增所得等位基因数在4~13之间（表3），平均每个位点7.3个，各引物在各居群内均呈现多态性。其中BPPCT002和BPPCT005扩增等位基因数最少，为4个；BPPCT037最多，为13个。来源自甜樱桃的SSR引物多态性高于来源自桃的SSR引物（平均扩增等位基因数分别为7.75和6.4）。图1为引物BPPCT014在壶瓶山（HPS）居群的扩增图谱。居群在所有位点均偏离Hardy-Weinberg平衡（P<0.05），从表3可以看出莱特固定指数（Fis）都为负值，即说明所有位点均表现出显著的杂合子过剩。图1引物BPPCT014在壶瓶山（HPS）居群上的扩增图谱M：DNA标准分子量；1~20：HPS居群1~20号个体。Fig.1AmplificationpatternofHPSpopulationusingprimerBPPCT014M：50bpDNAmarker；1–20：ThesamplesfromHPSpopulation.表311个SSR引物位点的种群分化和杂合度Table3Estimatesofpopulationdifferentiationandheterozygosityat11SSRprimers来源Origin基因座LocusFis莱特固定指数Wright'sfixationindexFst基因分化系数CoefficientofgenedifferentiationNa观测等位基因TheobservednumberofallelesNe有效等位基因EffecitvenumberofallelesNm基因流GeneflowHo可观察杂合度Observedheterozy-gosityHe期望杂合度ExpectedheterozygosityAn无效等位基因频率NullallelesfrequencyEPDCU5060–0.35730.07829.03.06442.94670.84290.67490.10BPPCT002–0.45590.10364.01.94502.16290.97140.74570.12UDP98-022–0.25340.11447.05.90361.93510.53930.48670.03BPPCT014–0.51430.16366.04.45291.27820.98210.77680.11桃P.persicaBPPCT034–0.42910.10366.03.99852.16370.96070.75120.11UCD-CH14–0.17080.17615.03.79951.16960.71070.73810.01EMPaS12–0.29750.18469.05.02291.10460.91790.86910.02BPPCT005–0.19010.10694.02.90032.08900.69640.65640.02甜樱桃P.aviumBPPCT037–0.21510.145413.06.96521.46900.88930.85800.01李P.salicinaCPSCT029–0.37020.168412.07.54791.23430.94640.83210.06中国樱桃P.pseudocerasusPTCR18–0.33420.11765.03.91171.87520.94290.80230.07平均Mean–0.32990.13567.34.50111.59410.85450.74472.2遗传多样性11对引物在居群间扩增丰富，各居群多态位点比率为100%。在居群水平上，可观察等位基因数（Na）、有效等位基因数（Ne）、香侬信息指数（I）、基因多样性指数（h）均值分别为7.2727、4.5011、1.2069和0.6440（表3，表4）。FHS（PPL=100%，I=1.4684，h=0.7151）遗传多样性最高，TSGaoTian-xiang，CaiYu-liang，FengYing，ZhaoXiao-jun.GeneticdiversityandgeneticstructureofPrunuspseudocerasuspopulationsfromChinaasrevealedbySSRmarkers.1152ActaHorticulturaeSinica，2016，43(6)：1148–1156.（PPL=100%，I=0.9057，h=0.5584）和LS（PPL=100%，I=0.9854，h=0.5431）遗传多样性相对较低（表4）。基因多样性指数（h）是衡量居群遗传多样性最常用的指标，根据基因多样性指数得到中国樱桃14个自然居群遗传多样性水平：FHS>HPS>SNJ>TB>YJ>WB>BTM>JYS>NZ>YX>HLT>XZP>TS>LS（表4）。表4基于11对SSR引物的中国樱桃14个自然居群的遗传多样性水平Table4Geneticdiversityforthe14populationsofPrunuspseudocerasusbasedonSSRmakersfrom11primers居群PopulationPPL多态性比率/%PercentageofpolymorphiclociHo观察杂合度ObservedheterozygosityHe期望杂合度Expectedheterozygosityh基因遗传多样性Nei’sexpectedheterozygostyI香侬信息指数Shannon'sinformationindexBTM1000.84090.68740.67021.2926FHS1000.88640.73340.71511.4684HLT1000.79550.58960.57491.0453HPS1000.85000.70980.69201.3315JYS1000.81360.66080.66431.2489LS1000.77730.55700.54310.9845NZ1000.77270.67770.66081.2762SNJ1000.86820.70520.68761.3097TB1000.90000.69850.68101.3062TS1000.88180.57270.55840.9057WB1000.92730.68920.67191.2585XZP1000.98640.58150.56690.9307YJ1000.81360.69550.67811.3090YX1000.85000.66820.65151.2297平均Mean1000.85460.65900.64401.2069物种水平Specieslevel1000.85450.74470.74331.59122.3居群遗传结构14个自然居群的Nei’s遗传距离介于0.1217（YJ和JYS）到0.8057（XZP和HLT）之间（表5）。采用UPGMA法进行聚类分析，进一步直观分析居群间的遗传关系（图2），可将14个自然居群归为7类：其中宝天曼（BTM）、王坝（WB）、印江（YJ）、缙云山（JYS）、南郑（NZ）、神农架（SNJ）、壶瓶山（HPS）归为一类，崂山（LS）和泰山（TS）归为一类，凤凰山（FHS）、映秀（YX）、黑龙潭（HLT）、西张坡（XZP）、桐柏（TB）5个居群分别独自归为一类。用Structure进行群体遗传结构分析，在K=1~15中，K=7时拟然值最大（Evannoetal.，2005），推断14个自然居群来自7个理论组群（图3），与UPGMA聚类结果一致。Mantel检验显示，总群体的遗传距离和地理距离显著相关（r=0.472，P=0.011）。表5居群间遗传一致度和遗传距离的Nei’s无偏估计值表Table5Nei’sgeneticidentityandgeneticdistancebetweenpopulations居群PopulationBTMFHSHLTHPSJYSLSNZSNJTBTSWBXZPYJYXBTM0.72480.62000.79980.80000.78350.84310.79640.71130.74460.80330.69700.79220.6984FHS0.32190.61700.77210.71840.65070.75430.73000.62410.59940.74420.63200.72120.7104HLT0.47800.48290.65660.52890.61170.56680.63040.50840.58170.58790.44680.61300.5680HPS0.22340.25860.42070.77530.74870.83730.82470.69670.69150.76740.72810.81330.7666JYS0.22320.33070.63690.25450.66040.87730.77680.66880.59270.77910.66720.88540.7189LS0.24400.42970.49150.28940.41480.69070.78650.62360.88360.76100.62230.66230.7372NZ0.17060.28190.56780.17750.13090.37000.79940.68260.63050.81090.70000.83630.7402SNJ0.22760.31470.46150.19270.25260.24010.22390.69690.80100.80510.62320.79910.7516TB0.34070.47150.67650.36140.40220.47230.38180.36110.64330.71100.62600.66210.6540TS0.29490.51180.54180.36890.52310.12370.46120.22190.44120.70090.61070.61500.6872WB0.21910.29550.53120.26470.24960.27310.20960.21680.34110.35540.67350.74320.7240XZP0.36090.45890.80570.31730.40460.47440.35670.47300.46840.49320.39530.70000.5566YJ0.23290.32680.48940.20670.12170.41200.17880.22430.41240.48620.29680.35670.7250YX0.35900.34190.56560.26580.33000.30480.30090.28550.42460.37510.32290.58590.3216注：右上为遗传一致度；左下为遗传距离。Notes：Toprightdiagonal：Geneticidentity；Leftbottomdiagonal：Geneticdistance.高天翔，蔡宇良，冯瑛，赵晓军.中国樱桃14个自然居群遗传多样性和遗传结构的SSR评价.园艺学报，2016，43(6)：1148–1156.1153图2基于Nei’s遗传距离的中国樱桃14个自然居群的UPGMA聚类图Fig.2UPGMAdendrogramof14populationsofPrunuspseudocerasusbasedonNei’sgeneticdistance图3基于Structure分析的中国樱桃14个自然居群的遗传结构图7种颜色代表7个不同的遗传组群（Ⅰ~Ⅶ）；细竖线代表个体，不同颜色所占比例越大，划分到相应组群的可能性越大。Fig.3Geneticstructuralmapof14PrunuspseudocerasuspopulationsbasedonstructureanalysisSevencolorsrepresentedsevendifferentgeneticclusters（Ⅰ–Ⅶ）；Individualswererepresentedasthinverticallinespartitionedintosegmentscorrespondingtotheirmembershipingeneticclustersindicatedbythecolors.2.4居群遗传分化遗传分化是反映遗传结构的重要指标。中国樱桃14个自然居群间基因分化系数平均值Fst为0.1356（表3），低于双子叶植物基因分化系数的平均值（0.273）（蔡宇良，2006）。AMOVA结果（表6）表明，14个自然居群内、居群间均存在显著的遗传变异（P<0.001），群体间的遗传变异占总变异的30.00%，居群内占70.00%，说明遗传变异主要分布在居群内。表6野生中国樱桃居群的分子变异方差分析Table6Analysisofmolecularvariance（AMOVA）analysisamongPrunuspseudocerasuspopulations变异来源Sourceofvariation自由度df平方和Sumofsquare方差分量Variancecomponent方差分量百分率/%PercentageofvariancecomponentP值Pvalue居群间Amongpopulations13619.7932.13530.00<0.001居群内Withinpopulations2661325.9004.98570.00<0.001GaoTian-xiang，CaiYu-liang，FengYing，ZhaoXiao-jun.GeneticdiversityandgeneticstructureofPrunuspseudocerasuspopulationsfromChinaasrevealedbySSRmarkers.1154ActaHorticulturaeSinica，2016，43(6)：1148–1156.3讨论根据已经获得的评价遗传多样性的重要指标，在物种水平上Shannon’s信息指数（I）为1.5912，Nei’s基因多样性指数（h）为0.7433，与近缘物种相比，明显高于杏（I=0.458，h=0.287；Heetal.，2007）、马哈利樱桃（I=0.1720，h=0.1144；Jordano&Godoy，2000）及新疆野苹果（I=0.408，h=0.262；张春雨等，2009），略高于四川野生中国樱桃（I=1.5256，h=0.6953；陈娇等，2013）；低于欧洲甜樱桃（I=1.9963，h=0.7695；艾呈祥等，2007；Zhongetal.，2009），因此不能否认，本研究中中国樱桃14个自然居群确实也具有较高的遗传多样性。而且从根据基因多样性指数得到的中国樱桃14个自然居群遗传多样性水平中，可以看出多数分布在秦岭以南的居群遗传多样性要高于分布在秦岭以北的居群。地理分布范围是决定物种遗传多样性的一个主要因素，一个物种的分布范围越广，其遗传多样性水平越高（Karron，1987）。而中国樱桃在中国广泛分布，这也许是其具有较高遗传多样性水平的原因之一。本研究中所有位点（P<0.05）都显著偏离了Hardy-Weinberg平衡，表现为杂合子过剩，不同于陈娇等（2013）杂合度不足的研究结果。遗传分化是反映遗传结构的重要指标。本研究中居群间遗传分化系数为0.1356，稍高于陈娇等（2013）的0.0844。AMOVA分析表明中国樱桃自然居群的遗传变异主要存在于居群内（70.00%），这与蔡宇良（2006）和陈娇等（2013）的研究结果一致，在毛樱桃的研究中也有类似结果（张琪静，2007；宛甜等，2013）。种子散布机制对居群间的遗传分化有显著影响，中国樱桃主要依靠动物消化系统传播种子，其种子流受到动物活动范围的制约（蔡宇良，2006）。利用等位酶标记研究发现，种子的扩散方式对于居群间和居群内的遗传变异有很大的影响，通过动物食果进行种子扩散的植物种与以其他方式扩散种子的植物种相比，具有典型的高水平的居群内遗传变异（Hamricketal.，1993；Hamrick&Godt，1996）。本研究中发现中国樱桃自然居群内存在着较高的遗传变异，与前人结论（Hamricketal.，1993；Hamrick&Godt，1996）一致。基因流是影响居群遗传分化的重要因素，研究结果普遍认为基因流大于1时就能够抵制群体内的遗传漂变作用，也同时防止了群体间分化的发生（Slatkin，1985）。本研究中中国樱桃自然居群的基因流（Nm=1.5941）>1，基于Structure分析的14个自然居群的遗传结构图中各个颜色有所交织，也表明了各个居群间存在一定的基因交流，抑制了中国樱桃自然居群间的遗传变异，所以中国樱桃自然居群的遗传变异主要存在于居群内。居群聚类中，黑龙潭（HLT）、泰山（TS）和崂山（LS）与其他采样地区距离相对较远，进行的基因交流相对较少，与其他居群之间产生了一定程度的遗传分化，因此独自聚为一类。野生中国樱桃属于通过食果动物体内消化系统进行种子扩散的树种（蔡宇良，2006），西张坡（XZP）、凤凰山（FHS）、映秀（YX）、桐柏（TB）4个居群并没有与其附近的居群聚为一类，而是各自聚为一类，推测很可能是由于该地区与其他采样地区相比，人类活动范围相对较大且比较频繁，影响了动物的活动范围，不能形成有效的种子流，导致这4个居群分别与其附近的居群之间产生了一定程度的遗传分化。类似于赵阿曼等（2003）的研究结果：由于环境恶化和人类活动干扰（过度砍伐、放牧）等导致砂生槐生境片断化，影响居群间基因交流。中国樱桃不但具有较好的经济效益，而且还具有很高的观赏价值，但由于受环境变化和人类活动影响，中国樱桃自然资源大量减少，所以很有必要对中国樱桃采取一定的保护措施。因为遗传多样性对种群和物种的长期生存和发展起着关键作用，所以保护工作应该以保持最大遗传变异基因库为目标（Zhouetal.，2010）。本研究中秦岭以南的中国樱桃自然居群遗传多样性水平更高，因此应重点对秦岭以南的野生中国樱桃自然居群进行保护。高天翔，蔡宇良，冯瑛，赵晓军.中国樱桃14个自然居群遗传多样性和遗传结构的SSR评价.园艺学报，2016，43(6)：1148–1156.1155ReferencesAiCX，YuXM，DahlenburgA，ZhaoHJ，ZhaoY，LiuQZ.2008.DevelopmentandcharacterizationofSSRmarkersinChinesecherry（PrunuspseudocerasusLindl.）.EuropeanJournalofHorticulturalScience，73(3)：104–110.AiCheng-xiang，XinLi，YuXian-mei，ZhangLi-si，WeiHai-rong，YuanKe-jun，SunQing-rong，LiuQing-zhong.2007.AnalysisofgeneticdiversityinCarasusbySSRmarkers.ActaHorticulturaeSinica，34(4)：871–876.(inChinese)艾呈祥，辛力，余贤美，张力思，魏海荣，苑克俊，孙清荣，刘庆忠.2007.樱桃主栽品种的遗传多样性分析.园艺学报，34(4)：871–876.CaiYu-liang.2006.GeneticanalysisofthewildcherrygermplasmandidentificationofcultivatedcherryvarietiesusingDNAfingerprints[Ph.D.Dissertation].Xi’an：NorthwestUniversity.(inChinese)蔡宇良.2006.野生樱桃种质资源的遗传分析及其栽培品种的DNA指纹鉴定[博士论文].西安：西北大学.CaoDong-wei，CaiYu-liang，YangJuan，ZhaoGui-fang.2007.PCR-RFLPanalysisofPrunuspseudocerasus.JournalofNorthwestA&FUniversity：NaturalScienceEdition，35(5)：173–178.(inChinese)曹东伟，蔡宇良，杨娟，赵桂仿.2007.中国樱桃的PCR-RFLP分析.西北农林科技大学学报：自然科学版，35(5)：173–178.ChenJiao，WangXiao-rong，TangHao-ru，ChenTao，HuangXiao-jiao，LiangQin-biao.2013.AssessmentofgeneticdiversityandpopulationsgeneticstructureinwildChinesecherryfromSichuanProvinceusingSSRmarkers.ActaHorticulturaeSinica，40(2)：333–340.(inChinese)陈娇，王小蓉，汤浩茹，陈涛，黄晓娇，梁勤彪.2013.基于SSR标记的四川野生中国樱桃遗传多样性和居群遗传结构分析.园艺学报，40(2)：333–340.DirlewangerE，CossonP，TavaudM，AranzanaM，PoizatC，ZanettoA，LaigretF.2002.Developmentofmicrosatellitemarkersinpeach[Prunuspersica（L.）Batsch]andtheiruseingeneticdiversityanalysisinpeachandsweetcherry（PrunusaviumL.）.TheoreticalandAppliedGenetics，105(1)：127–138.EvannoG，RegnautS，GoudetJ.2005.DetectingthenumberofclustersofindividualsusingthesoftwareSTRUCTURE：asimulationstudy.MolecularEcology，14(8)：2611–2620.ExcoffierL，SmousePE，QuattroJM.1992.AnalysisofmolecularvarianceinferredfrommetricdistancesamongDNAhaplotypes：applicationtohumanmitochondrialDNArestrictiondata.Genetics，131(2)：479–491.GuarinoC，SantoroS，DeSimoneL，CiprianiG.2009.Prunusavium：nuclearDNAstudyinwildpopulationsandsweetcherrycultivars.Genome，52(4)：320–337.HamrickJL，GodtMJW.1996.Effectsoflifehistorytraitsongeneticdiversityinplantspecies.PhilosophicalTransactionsoftheRoyalSocietyofLondonB：BiologicalSciences，351(1345)：1291–1298.HamrickJL，MurawskiDA，NasonJD.1993.Theinfluenceofseeddispersalmechanismsonthegeneticstructureoftropicaltreepopulations.Frugivoryandseeddispersal：ecologicalandevolutionaryaspects.SpringerNetherlands，15：281–297.HeTian-ming，ChenXue-sen，XuZheng，GaoJiang-sheng，LinPei-jun，LiuWen，LiangQing，WuYan.2007.UsingSSRmarkerstodeterminethepopulationgeneticstructureofwildapricot（PrunusarmeniacaL.）intheIlyValleyofWestChina.GeneticResourcesandCropEvolution，54(3)：563–572.HeWen，ZhangJing，HuangZhi-lin，ChenQing，TangHao-ru，TangFu-yi，WangXiao-rong.2014.GeneticdiversityandpopulationgeneticstructureamonglocalChinesecherryvarieties[Cerasuspseudocerasus（Lindl.）G.Don]basedonITSsequence.ActaBotanicaBoreali-OccidentaliaSinica，34(3)：463–472.(inChinese)何文，张静，黄智林，陈清，汤浩茹，汤福义，王小蓉.2014.基于ITS序列对栽培中国樱桃遗传多样性及其群体遗传结构的分析.西北植物学报，34(3)：463–472.HuangXiao-jiao，ChenTao，LiangQin-biao，ChenJiao，WangXiao-rong，TangHao-ru.2011.Preliminaryeconomiccharactersevaluationto14cherrygermplasms（Cerasuspseudocerasus）.SouthChinaFruits，40(5)：9–12.(inChinese)黄晓姣，陈涛，梁勤彪，陈娇，王小蓉，汤浩茹.2011.14份中国樱桃种质经济性状的初步评价.中国南方果树，40(5)：9–12.HuangXiao-jiao，WangXiao-rong，ChenTao，ChenJiao，TangHao-ru.2013.ResearchprogressofgermplasmdiversityinChinesecherry（Cerasuspseudocerasus）.JournalofFruitScience，30(3)：470–479.(inChinese)黄晓姣，王小蓉，陈涛，陈娇，汤浩茹.2013.中国樱桃遗传资源多样性研究进展.果树学报，30(3)：470–479.GaoTian-xiang，CaiYu-liang，FengYing，ZhaoXiao-jun.GeneticdiversityandgeneticstructureofPrunuspseudocerasuspopulationsfromChinaasrevealedbySSRmarkers.1156ActaHorticulturaeSinica，2016，43(6)：1148–1156.IllaE，LambertP，QuilotB，AudergonJM，DirlewangerE，HowadW，BassiD.2009.Linkagemapsaturation，construction，andcomparisoninfourpopulationsofPrunus.JHorticulturalSciBiotechnolISAFRUITSpecIssue，84(6)：168–175.JordanoP，GodoyJA.2000.RAPDvariationandpopulationgeneticstructureinPrunusmahaleb（Rosaceae），ananinmal-dispersedtree.MolecularEcology，9(9)：1293–1305.KarronJD.1987.Acomparisonoflevelsofgeneticpolymorphismandself-compatibilityingeographicallyrestrictedandwidespreadplantcongeners.EvolutionaryEcology，1(1)：47–58.LewontinRC.1972.Theapportionmentofhumandiversity.EvolBiol，6：381–398.McDermottJM，McDonaldBA.1993.Geneflowinplantpathosystems.AnnuRevPhytopathol，31：353–373.MnejjaM，Garcia-MasJ，HowadW，BadenesML，ArúsP.2004.Simple-sequencerepeat（SSR）markersofJapaneseplum（PrunussalicinaLindl.）arehighlypolymorphicandtransferabletopeachandalmond.MolecularEcologyNotes，4(2)：163–166.PeakallR，SmousePE.2006.GENALEX6：geneticanalysisinExcel.Populationgeneticsoftwareforteachingandresearch.MolEcolNotes，6：288–295.PritchardJK，StephensM，DonnellyP.2000.Inferenceofpopulationstructureusingmultilocusgenotypedata.Genetics，155(2)：945–959.SlatkinM.1985.Geneflowinnaturalpopulations.AnnualReviewofEcologyandSystematics，16：393–430.SmousePE，LongJC.1992.Matrixcorrelationanalysisinanthropologyandgenetics.YearbookPhysAnthropol，35：187–213.SmousePE，LongJC，SokalRR.1986.MultipleregressionandcorrelationextensionsoftheManteltestofmatrixcorrespondence.SystematicZoology，35：627–632.SongChang-mei，WenXiao-peng，YangEr-tai.2011.C