基于SSR标记的桃种质资源遗传多样性研究

江西农业学报㊀２０１８，３０（１１）：１４ １８
ＡｃｔａＡｇｒｉｃｕｌｔｕｒａｅＪｉａｎｇｘｉ
㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀
ｈｔｔｐ：／／ｗｗｗ．ｊｘｎｙｘｂ．ｃｏｍ
ＤＯＩ：１０．１９３８６／ｊ．ｃｎｋｉ．ｊｘｎｙｘｂ．２０１８．１１．０４
基于ＳＳＲ标记的桃种质资源遗传多样性研究
凌士鹏，孙萍，林贤锐，沈建生∗
㊀
㊀收稿日期：２０１８－０６－２１
基金项目：浙江省农业（果品）新品种选育重大科技专项（２０１６Ｃ０２０５２－５）㊂
作者简介：凌士鹏（１９６３─），男，浙江兰溪人，高级农艺师，主要从事果树种质资源收集及栽培推广㊂∗通讯作者沈建生㊂
（浙江省金华市农业科学研究院，浙江金华３２１０００）
摘㊀要：评价了桃种质资源的遗传多样性水平及遗传结构，为其保存提供理论依据㊂利用１２对ＳＳＲ引物对５３份桃种质资源进行了遗传多样性分析，利用Ｎｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇ聚类对其遗传结构进行了分析比较㊂结果表明：（１）１２对引物在５３个样品中的平均等位基因数（Ｎａ）为７．９１７，Ｓｈａｎｎｏｎ多样性指数（Ｉ）为１．４０９，平均有效等位基因数（Ｎｅ）为３．３９９，平均观察杂合度（Ｈｏ）为０．６５７，期望杂合度（Ｈｅ）平均值为０．６７８㊂（２）基于Ｎｅｉ距离的Ｎｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇ聚类分析表明：５３个桃样品被分为３个组，同时能清楚地鉴别出同物异名的种质资源㊂由此可知，供试的桃种质资源表现出较高的遗传多样性水平；大部分遗传关系较近的品种以及部分形态学性状相近的品种聚在一类，说明聚类分析结果与系谱及生物学特征具有一定的关联性，同时聚类结果与品种的地理分布具有一定相关性㊂
关键词：桃；ＳＳＲ；遗传多样性；遗传结构
中图分类号：Ｓ６６２．１㊀文献标志码：Ａ㊀文章编号：１００１－８５８１（２０１８）１１－００１４－０５
ＡＳｔｕｄｙｏｎＧｅｎｅｔｉｃＤｉｖｅｒｓｉｔｙｏｆＰｅａｃｈＧｅｒｍｐｌａｓｍ
ＲｅｓｏｕｒｃｅｓＢａｓｅｄｏｎＳＳＲＭａｒｋｅｒｓ
ＬＩＮＧＳｈｉ－ｐｅｎｇ，ＳＵＮＰｉｎｇ，ＬＩＮＸｉａｎ－ｒｕｉ，ＳＨＥＮＪｉａｎ－ｓｈｅｎｇ∗
（ＪｉｎｈｕａＡｃａｄｅｍｙｏｆＡｇｒｉｃｕｌｔｕｒａｌＳｃｉｅｎｃｅｓｉｎＺｈｅｊｉａｎｇ，Ｊｉｎｈｕａ３２１０００，Ｃｈｉｎａ）
Ａｂｓｔｒａｃｔ：Ｔｈｅｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙａｎｄｇｅｎｅｔｉｃｓｔｒｕｃｔｕｒｅｏｆｐｅａｃｈｇｅｒｍｐｌａｓｍｒｅｓｏｕｒｃｅｓｗｅｒｅａｓｓｅｓｓｅｄｔｏｐｒｏｖｉｄｅｔｈｅｏｒｅｔｉｃａｌｂａｓｉｓｆｏｒｔｈｅｉｒｃｏｎｓｅｒｖａｔｉｏｎ．ＴｗｅｌｖｅｐａｉｒｓｏｆＳＳＲ（ＳｉｍｐｌｅＳｅｑｕｅｎｃｅＲｅｐｅａｔｓ）ｐｒｉｍｅｒｓｗｅｒｅｕｓｅｄｔｏｅｖａｌｕａｔｅｔｈｅｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙｏｆ５３ｐｅａｃｈａｃｃｅｓｓｉｏｎｓ，ａｎｄＮｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇｃｌｕｓｔｅｒｉｎｇｗａｓｕｓｅｄｔｏａｎａｌｙｚｅａｎｄｃｏｍｐａｒｅｔｈｅｉｒｇｅｎｅｔｉｃｓｔｒｕｃｔｕｒｅ．Ｔｈｅｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔ：（１）Ｔｈｅｍｅａｎｎｕｍｂｅｒｏｆｐｏｌｙｍｏｒｐｈｉｃａｌｌｅｌｅｓ（Ｎａ）ｏｆ１２ｐａｉｒｓｏｆｐｒｉｍｅｒｓｉｎ５３ｐｅａｃｈａｃｃｅｓｓｉｏｎｓｗａｓ７．９１７，ａｎｄｔｈｅａｖｅｒ⁃ａｇｅｎｕｍｂｅｒｏｆｅｆｆｅｃｔｉｖｅａｌｌｅｌｅｓ（Ｎｅ）ｐｅｒｌｏｃｕｓｗａｓ３．３９９；ｔｈｅＳｈａｎｎｏｎｄｉｖｅｒｓｉｔｙｉｎｄｅｘ（Ｉ）ｗａｓ１．４０９；ｔｈｅａｖｅｒａｇｅｏｂｓｅｒｖｅｄｈｅｔ⁃ｅｒｏｚｙｇｏｓｉｔｙ（Ｈｏ）ｗａｓ０．６５７，ａｎｄｔｈｅａｖｅｒａｇｅｅｘｐｅｃｔｅｄｈｅｔｅｒｏｚｙｇｏｓｉｔｙ（Ｈｅ）ｗａｓ０．６７８．（２）ＴｈｅＮｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇｃｌｕｓｔｅｒｉｎｇａ⁃ｎａｌｙｓｉｓｂａｓｅｄｏｎＮｅｉ ｓｄｉｓｔａｎｃｅｄｉｖｉｄｅｄ５３ｐｅａｃｈｓａｍｐｌｅｓｉｎｔｏ３ｇｒｏｕｐｓ，ａｎｄｉｔｃｏｕｌｄｃｌｅａｒｌｙｉｄｅｎｔｉｆｙｔｈｅｓｙｎｏｎｙｍｓｏｆｔｈｅｓｅｇｅｒｍ⁃ｐｌａｓｍｒｅｓｏｕｒｃｅｓ．Ｔｈｅｔｅｓｔｅｄｐｅａｃｈｇｅｒｍｐｌａｓｍｒｅｓｏｕｒｃｅｓｅｘｈｉｂｉｔｅｄａｈｉｇｈｌｅｖｅｌｏｆｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙ，ａｎｄｍｏｓｔｏｆｔｈｅｖａｒｉｅｔｉｅｓｗｉｔｈｃｌｏｓｅｒｇｅｎｅｔｉｃｒｅｌａｔｉｏｎｓｈｉｐｏｒｓｏｍｅｓｉｍｉｌａｒｍｏｒｐｈｏｌｏｇｉｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓｗｅｒｅｃｌｕｓｔｅｒｅｄｉｎｔｏｔｈｅｓａｍｅｇｒｏｕｐ，ｉｎｄｉｃａｔｉｎｇｔｈａｔｔｈｅｒｅ⁃ｓｕｌｔｓｏｆｃｌｕｓｔｅｒｉｎｇａｎａｌｙｓｉｓｗｅｒｅｒｅｌａｔｅｄｔｏｔｈｅｇｅｎｅａｌｏｇｙａｎｄｂｉｏｌｏｇｉｃａｌｃｈａｒａｃｔｅｒｉｓｔｉｃｓ．Ａｔｔｈｅｓａｍｅｔｉｍｅ，ｔｈｅｃｌｕｓｔｅｒｉｎｇｒｅｓｕｌｔｓｗｅｒｅａｌｓｏｃｏｒｒｅｌａｔｅｄｗｉｔｈｔｈｅｇｅｏｇｒａｐｈｉｃａｌｄｉｓｔｒｉｂｕｔｉｏｎｏｆｔｈｅｖａｒｉｅｔｉｅｓ．
Ｋｅｙｗｏｒｄｓ：Ｐｅａｃｈ；ＳＳＲ；Ｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙ；Ｇｅｎｅｔｉｃｓｔｒｕｃｔｕｒｅ
㊀㊀桃是蔷薇科（Ｒｏｓａｃｅａｅ）李属（Ｐｒｕｎｕｓ）植物，起源地是中国的西部地区，分布点主要是欧㊁亚㊁美㊁澳㊁
非等五大洲，其种植范围特别广㊁适应性比较强㊁种质资源极为丰富［１］㊂桃经过多年的自然选择和人工驯化种植，形成了许多的类型和品系，这些资源为进一步改良桃的相关性状提供了基础，了解这些资源是利用的前提［２］㊂目前，相关研究表明：我国拥有北京㊁南京㊁郑州３个国家级桃种质资源圃以及很多地方桃种质资源圃，收集并保存了６个种，共２０００多份种质资源［３］，桃作为世界上栽培最为广泛的落叶果
树之一，现已遍布全球，居核果类果树之首［４］㊂自２０
世纪８０年代以来，国内很多学者对桃亚属种间分类进行了研究，目前已经在形态学［１］㊁酶学［５］㊁孢粉学［６］㊁细胞学［７］㊁分子标记［８－１１］等方面获得了一些桃资源分类及亲缘关系等的结论和观点㊂但是研究发现，从形态学方面对品种进行分类需要看经验，因此对于一些种间的杂种进行归类十分困难，更无法判断其起源；另外，利用同工酶㊁染色体核型等试验手段也同样面临着更多的问题㊂因此开展桃分子标记辅助分类具有重要意义㊂
简单序列重复（ＳｉｍｐｌｅＳｅｑｕｅｎｃｅＲｅｐｅａｔ，ＳＳＲ）标
记的特点是等位基因变异数多㊁信息含量及多态性
高，稳定性㊁重复性好，属共显性遗传，特异性强，并
且在种属间有良好的通用性［１４］，是目前果树上鉴定种质资源㊁分析遗传多样性及亲缘关系应用中最为广
泛的分子标记㊂同时，Ｖｅｒｄｅ等［１２］基于桃不同器官转录组数据研究后上传至ＮＣＢＩ共８０７９７条ＥＳＴｓ序列，为挑选更好的桃ＳＳＲ引物提供了丰富的资源㊂笔者所在课题组已于前期对部分桃种质资源进行了调查㊁收集保存，因此在此基础上，本研究利用ＳＳＲ分子标记技术对收集保存的桃种质资源进行遗传多样性的评价，以期为后续桃种质资源的综合利用提供借鉴㊂１㊀材料与方法
１．１㊀植物材料
于各地收集桃地方品种（品系）共５３份，嫁接于
浙江省金华国家高科技农业园区的桃李品种园，详细
内容见表１㊂试验材料为健康幼嫩叶片，采集并分别
置于采集袋内，同时利用变色硅胶迅速干燥，带回实
验室，挑取干燥的叶片进行ＤＮＡ提取㊂
表１㊀供试的桃品种（品系）
编号品种（品系）编号品种（品系）１５１８油桃３４天鲜红
２超丽春３５香梅露
３千年红３６湖景蜜露
４双喜红３７硬肉桃
５京东巨油３８日本桃
６中油１４号３９新玉
７超早油桃４０雪雨露
９红芸果油桃４１晚熟１号桃１１金山早红４２２号桃
１２艳光４３特早４号
１３中油５号４４八仙蜜露
１４大观４５仙岛红２号１５金华大白桃４６东溪小仙桃１６早霞露４７Ｘ２－５
１７江山早４８Ｘ
１８沪油０１８４９早露蟠
１９沙红５０玉露蟠
２３夏香姬５１红冠
２４加纳桃５２叶１０月桃２５白玉桃５３天台１０月桃２６新白花５４春蜜
２７夏之梦５５春瑞
２８南国红５６Ｔ突围
３０红清水５７锦秀
３２上汗红桃５８新川岛
３３红燕露５９春雪
６０锦圆１．２㊀方法
１．２．１㊀ＤＮＡ的提取㊀采用ＣＴＡＢ法［１３］从采集的叶片中提取总ＤＮＡ，用双蒸水稀释到１０ ３０ｎｇ／μＬ，储存于－２０ħ冰箱中备用㊂
１．２．２㊀ＳＳＲ标记㊀随机抽取８个桃样品的ＤＮＡ模板，选用５６对ＳＳＲ引物进行ＰＣＲ扩增，将扩增产物利用６％的聚丙烯酰胺凝胶电泳约２．５ｈ，后银染显色㊂从中筛选出１２对扩增条带清晰㊁稳定的引物进行ＳＳＲ扩增，详见表２㊂试验前对筛选好的１２对引物进行荧光修饰，即在单向引物的５ᶄ端或３ᶄ端加上荧光修饰基团，便可实现微卫星多态性的荧光半自动检测㊂荧光引物由上海英骏生物技术有限公司合成㊂㊀㊀ＰＣＲ反应体系及ＳＳＲ扩增程序详见孙萍等［１７］２０１７年发表的文章㊂扩增结束后对ＰＣＲ产物用灭菌的双蒸水稀释１００ ２００倍，利用ＡＢＩ３１３０测序仪进行毛细管电泳检测；数据采集采用ＧｅｎｅＭａｐｐｅｒ４．０仪器在电泳结果６０ ３５０ｂｐ片段大小范围内进行㊂１．３㊀数据分析
将测序的结果导入ＧｅｎｅｔｉｃＰｒｏｆｉｌｅｖ１．０软件中，参照峰值的强度及分子内标大小，确定微卫星位点扩增片段的区间，对片段大小进行统计，如果出现多峰型的样品，则需重新进行ＰＣＲ试验后重新检测，以Ｅｘｃｅｌ格式保存统计结果，以备后续分析［１７］㊂
利用遗传数据分析软件ＧｅｎＡｌＥｘ６．３［１８］计算多态性等位基因数（Ｎａ）㊁ＳＳＲ位点的有效等位基因数（Ｎｅ）㊁Ｓｈａｎｎｏｎ多样性指数（Ｉ）㊁观察杂合度（Ｈｏ）㊁期望杂合度（Ｈｅ）等遗传多样性指标㊂
用Ｐｏｐｕｌａｔｉｏｎ１．２［１９］软件对５３个桃种质资源进行Ｎｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇ聚类分析，并描绘和修改聚类图［２０］㊂
２㊀结果与分析
２．１㊀ＳＳＲ扩增产物的多态性
从表３可以看出，１２对ＳＳＲ引物在所有桃样品中的多态性等位基因数（Ｎａ）从６到１１不等，等位基因数平均值为７．９１７㊂有效等位基因数（Ｎｅ）从１．７７３（ＣＰＰＣＴ０４２）到４．６５７（ＢＰＰＣＴ０２６），平均值为３．３９９㊂其中，引物ＣＰＰＣＴ０２２扩增出的等位基因数为１１个，可以区分５３个品种中的２２个品种；多态性较好的引物为ＣＰＰＣＴ０２２㊁ＢＰＰＣＴ０３８㊁ＣＰＰＣＴ０２６㊁ＢＰＰＣＴ０２５㊁ＵＤＰ９６－００５，能区分几乎所有５３个桃地方品种㊂在本研究中，观察杂合度（Ｈｏ）为０．４３６（ＣＰＰＣＴ０４２） ０．９２２（ＢＰＰＣＴ０３８），平均值为０．６５７㊂期望杂合度（Ｈｅ）为０．４３６（ＣＰＰＣＴ０４２） ０．７８５（ＢＰＰＣＴ０２６），平均期望杂合度为０．６７８㊂Ｓｈａｎｎｏｎ多样性指数（Ｉ）为０．９３９（ＣＰＰＣＴ０４２） １．６９８（ＢＰＰＣＴ０３８），平均值为１．４０９㊂反映了所收集保存的桃种质资源的遗传变异程度较高，遗传多样性较高㊂
５１
㊀１１期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀凌士鹏等：基于ＳＳＲ标记的桃种质资源遗传多样性研究
表２㊀ＳＳＲ引物及扩增的长度范围
引物名称引物序列（５ᶄ－３ᶄ）等位基因范围／ｂｐ参考文献ＢＰＰＣＴ００１Ｆ：ＡＡＴＴＣＣＣＡＡＡＧＧＡＴＧＴＧＴＡＴＧＡＧ，Ｒ：ＣＡＧＧＴＧＡＡＴＧＡＧＣＣＡＡＡＧＣ１２８ １７０［１４］ＢＰＰＣＴ０２５Ｆ：ＴＣＣＴＧＣＧＴＡＧＡＡＧＡＡＧＧＴＡＧＣ，Ｒ：ＣＧＡＣＡＴＡＡＡＧＴＣＣＡＡＡＴＧＧＣ１４７ １９９［１４］ＢＰＰＣＴ０１７Ｆ：ＴＴＡＡＧＡＧＴＴＴＧＴＧＡＴＧＧＧＡＡＣＣ，Ｒ：ＡＡＧＣＡＴＡＡＴＴＴＡＧＣＡＴＡＡＣＣＡＡＧＣ１３９ １８２［１４］ＢＰＰＣＴ００７Ｆ：ＴＣＡＴＴＧＣＴＣＧＴＣＡＴＣＡＧＣ，Ｒ：ＣＡＧＡＴＴＴＣＴＧＡＡＧＴＴＡＧＣＧＧＴＡ１３４ １６１［１４］ＢＰＰＣＴ０３８Ｆ：ＴＡＴＡＴＴＧＴＴＧＧＣＴＴＣＴＴＧＣＡＴＧ，Ｒ：ＴＧＡＡＡＧＴＧＡＡＡＣＡＡＴＧＧＡＡＧＣ１０６ １５３［１４］Ｐｃｈｇｍｓ６Ｆ：ＣＡＴＴＧＴＴＣＡＴＧＧＧＡＧＧＡＡＴＴ，Ｒ：ＡＧＡＡＣＡＴＴＣＣＴＡＡＡＧＧＡＧＣＡ２１６ ２７５［１５］ＣＰＰＣＴ０２６Ｆ：ＡＧＡＣＧＣＡＧＣＡＣＣＣＡＡＡＣＴＡＣ，Ｒ：ＣＡＴＴＡＣＡＴＣＡＣＣＧＣＣＡＡＣＡＡ１５０ ２２０［１５］ＣＰＰＣＴ０２２Ｆ：ＣＡＡＴＴＡＧＣＴＡＧＡＧＡＧＡＡＴＴＡＴＴＧ，Ｒ：ＧＡＣＡＡＧＡＡＧＣＡＡＧＴＡＧＴＴＴＧ２１７ ３０５［１５］ＣＰＰＣＴ０４２Ｆ：ＡＣＴＣＧＴＧＧＡＧＣＡＧＡＣＣＡＣＴＴ，Ｒ：ＡＴＣＣＡＡＧＣＡＴＧＣＣＡＡＧＡＡＡＧ１４９ ２００［１５］ＣＰＰＣＴ０２９Ｆ：ＣＣＡＡＡＴＴＣＣＡＡＡＴＣＴＣＣＴＡＡＣＡ，Ｒ：ＧＡＴＣＡＡＣＴＴＴＧＡＧＡＴＴＴＧＴＴＧＡＡ１７３ １９８［１５］ＵＤＰ９８－０２５Ｆ：ＧＧＧＡＧＧＴＴＡＣＴＡＴＧＣＣＡＴＧＡＡＴ，Ｒ：ＣＧＣＡＧＡＣＡＴＧＴＡＧＧＡＣＣＴＣ１２９ １３７［１６］ＵＤＰ９６－００５Ｆ：ＧＴＡＡＣＧＣＴＣＧＣＴＡＣＣＡＣＡＡＡ，Ｒ：ＣＣＴＧＣＡＴＡＴＣＡＣＣＡＣＣＣＡＧ１００ １７７［１６］
表３㊀５３个栽培桃品种１２对引物的遗传参数引物ＮａＮｅＩＨｏＨｅＢＰＰＣＴ０３８１０４．４８５１．６９８０．９２２０．７７７ＣＰＰＣＴ０４２７１．７７３０．９３９０．４３６０．４３６ＣＰＰＣＴ０２６９４．６５７１．６９１０．８５５０．７８５Ｐｃｈｇｍｓ６６３．１８３１．２９００．５２７０．６８６ＢＰＰＣＴ０１７７３．２６１１．３３１０．７１７０．６９３ＢＰＰＣＴ００７８３．５５７１．５３３０．６０００．７１９ＣＰＰＣＴ０２２１１４．００４１．６８１０．７２７０．７５０ＢＰＰＣＴ００１６２．９３５１．３３８０．５０００．６５９ＣＰＰＣＴ０２９６１．９７２０．９７３０．４５５０．４９３ＢＰＰＣＴ０２５９４．１１８１．５７６０．６７３０．７５７ＵＤＰ９８－０２５７２．６６２１．２５５０．５０００．６２４ＵＤＰ９６－００５９４．１７８１．６０７０．８８９０．７６１平均７．９１７３．３９９１．４０９０．６５７０．６７８２．２㊀聚类分析
通过ＳＳＲ数据对５３个地方桃品种进行聚类分析，构建亲缘关系图（图１）㊂从聚类分析图可以得出：５３个地方桃品种分成了３组，第一组中包括大观㊁江山早㊁加纳桃㊁红清水等１１个品种，其中大观和江山早的相似系数为０．９６，晚熟１号桃和２号桃的相似系数为０．９８，说明它们之间遗传关系很近；第二组包括了３５个品种，其中５１８油桃和千年红的相似系数为１．００，即为同物异名的现象；第三组包括了其余的９个品种，其中新玉和东溪小仙桃的相似系数为
０．９４，也表现出了很近的亲缘关系㊂
图１㊀基于Ｎｅｉ的遗传距离的Ｎｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇ聚类分析图
３㊀小结与讨论３．１㊀遗传多样性分析
杂合度的高低及等位基因数的多少是反映所选用ＳＳＲ位点鉴别植物基因型能力的重要指标［１７］㊂本
６１江㊀西㊀农㊀业㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀３０卷
研究结果表明：所选用的１２个ＳＳＲ引物其多态性较高，可以用于地方桃品种的鉴定和遗传多样性分析，另外也表明了笔者前期所收集保存的地方桃种质资源遗传变异较为丰富㊂
１２对引物在所有样品中等位基因的平均值为７．９１７，平均有效等位基因数（Ｎｅ）为３．３９９，平均观察杂合度（Ｈｏ）为０．６５７，期望杂合度（Ｈｅ）的平均值为０．６７８，与史红丽等［２１］的研究结果相似（Ｎｅ＝３．４），比Ｇｈａｆｆａｒｉ等［２２］在扁桃上的研究结果稍低（Ｎａ＝８．８，Ｈｏ＝０．７，Ｈｅ＝０．７９）㊁但显著高于Ａｒａｎｚａｎａ等［２３］在桃上的研究结果（Ｎａ＝６．６ ７．３，Ｈｏ＝０．３５，Ｈｅ＝０．５０ ０．５５）㊂另外，本研究中桃种质资源的遗传多样性水平与杏［２４］（Ｎａ＝３．５，Ｈｏ＝０．５８）㊁樱桃［２５］（Ｎａ＝３．５ ６．０，Ｈｅ＝０．４９ ０．６６）等其他李属植物的研究结果相比，有明显差异，表现出较高的遗传多样性水平㊂这为桃种质资源圃的建立奠定了物质基础㊂
３．２　遗传结构分析
基于Ｎｅｉ距离的Ｎｅｉｇｈｂｏｒ－Ｊｏｉｎｉｎｇ聚类分析清楚检测到桃供试样品的遗传结构：５３个地方桃品种分成了３组，其中第一组中的大观和江山早均为早熟桃品种，形态学性状相近；晚熟１号桃和２号桃均为晚熟桃品种，形态学性状也极其相近；另外，第一组中的加纳㊁红清水㊁上汗这３个桃品种均为日本选育的品种㊂第二组的５１８油桃和千年红的相似系数为１．００，为同物异名的现象；春瑞㊁春雪㊁春蜜㊁突围这４个桃均为早熟桃品种，形态学特征相似；红芸果㊁艳光㊁中油５号㊁超丽春㊁千年红㊁京东巨油㊁超早油桃这几个品种均为油桃品种，聚在一起；红冠㊁夏香姬㊁沙红㊁雪雨露这４个品种形态学特征相似，聚在一起㊂第三组中新玉和东溪小仙桃的相似系数为０．９４，也表现出了很近的遗传关系㊂
在不同的桃类群中，大部分遗传关系较近的品种以及部分形态学性状相近的品种聚在一起，说明其聚类分析结果与桃系谱以及生物学特征具有一定的关联性，但也不是完全一致的㊂部分形态学性状相差较大的桃品种，如锦绣和红燕露㊁天鲜红在聚类图中聚集在一起，说明这３个桃品种之间的基因型存在相互渗透的现象，亲缘关系较近㊂造成这种现象的原因可能是控制桃肉色性状的基因发生了几个碱基的变化，而通过ＳＳＲ分析检测不到碱基的微小变异，也可能与桃本身遗传结构相对复杂有关㊂
本研究结果也显示出聚类结果与桃品种的地理分布之间存在一定相关性，如日本桃品种㊁部分地方品种分别较为集中的聚在一起㊂但也有部分来自不同地区甚至地理距离较远的桃品种，其遗传距离表现非常相近，造成这种现象的原因可能与品种所处的环境条件相似，或者人工干预较为相近有一定的关系㊂因此在选择杂交亲本时，应注重扩大地域选择范围，这对于提高遗传多样性具有明显意义㊂
另外，研究还存在着普通桃和油桃聚为一类的交叉现象，说明各桃品种之间通过多年的自然选择，以及杂交育种过程中基因之间相互渗透，从而形成了较为丰富的遗传多样性，同时通过生态条件的变化以及人工干预，形成了不同的品种（品系）类型㊂
本研究采用ＳＳＲ标记来研究桃种质资源的遗传多样性，结果表明：ＳＳＲ标记是一种经济㊁可靠的分子标记技术㊂从本研究的结论可以看出，所收集到的桃种质资源遗传变异较为丰富，遗传多样性水平较高，对今后桃种质资源的保存和利用具有一定的参考价值㊂但本研究发现部分性状差别较大的地方桃品种在聚类图上聚为一类，这有待于进一步研究和论证㊂
参考文献：
［１］汪祖华，庄恩及．中国果树志 桃卷［Ｍ］．北京：中国林业出版社，２００１．
［２］叶宇芸，王清明，马建伟，等．基于ＳＳＲ标记的桃品种遗传多样性分析及遗传相似系数比较［Ｊ］．华南农业大学学报，２０１７，３８（３）：３９－４５．
［３］李雄伟，孟宪桥，贾惠娟，等．桃品种特异性荧光ＳＳＲ分子标记数据库构建［Ｊ］．果树学报，２０１３，３０（６）：９２４－９３２．［４］俞明亮，马瑞娟，沈志军，等．中国桃种质资源研究进展［Ｊ］．江苏农业学报，２０１０，２６（６）：１４１８－１４２３．
［５］高锁柱，马德伟，刘景芬，等．几种桃的过氧化物酶同工酶酶谱分析比较［Ｊ］．河北农业大学学报，１９８７（１）：２３－２７．［６］张秀英，王雁，王桂萍．桃花种质资源花粉形态的观察与比较［Ｊ］．北京林业大学学报，１９９７，１９（２）：５７－６２．［７］郭振怀，葛会波，王秀伶，等．桃属植物染色体核型及种间亲缘关系分析［Ｊ］．园艺学报，１９９６，２３（３）：２２３－２２６．［８］陈巍，王力荣，张绍铃，等．利用ＳＳＲ研究不同国家桃育成品种的遗传多样性［Ｊ］．果树学报，２００７，２４（５）：５８０－５８４．［９］杨英军，张开春，林珂．常见桃属植物ＲＡＰＤ多态性及亲缘关系分析［Ｊ］．河南农业大学学报，２００２，３６（２）：１８７－１９０．［１０］王富荣，佟兆国，赵剑波，等．桃野生种和地方品种种质资源亲缘关系的ＡＦＬＰ分析［Ｊ］．果树学报，２００８，２５（３）：３０５－３１１．
［１１］ＢｏｕｈａｄｉｄａＭ，ＭｏｒｅｎｏＭＡ，ＧｏｎｚａｌｏＭＪ，ｅｔａｌ．Ｇｅｎｅｔｉｃｖａｒ⁃ｉａｂｉｌｉｔｙｏｆｉｎｔｒｏｄｕｃｅｄａｎｄｌｏｃａｌＳｐａｎｉｓｈｐｅａｃｈｃｕｌｔｉｖａｒｓｄｅｔｅｒ⁃ｍｉｎｅｄｂｙＳＳＲｍａｒｋｅｒｓ［Ｊ］．ＴｒｅｅＧｅｎｅｔＧｅｎｏｍｅｓ，２０１１，７（２）：２５７－２７０．
［１２］ＶｅｒｄｅＩ，ＡｂｂｏｔｔＡＧ，ＳｃａｌａｂｒｉｎＳ，ｅｔａｌ．Ｔｈｅｈｉｇｈ－ｑｕａｌｉｔｙｄｒａｆｔｇｅｎｏｍｅｏｆｐｅａｃｈ（Ｐｒｕｎｕｓｐｅｒｓｉｃａ）ｉｄｅｎｔｉｆｉｅｓｕｎｉｑｕｅｐａｔｔｅｒｎｓｏｆｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙ，ｄｏｍｅｓｔｉｃａｔｉｏｎａｎｄｇｅｎｏｍｅｅｖｏｌｕ⁃ｔｉｏｎ［Ｊ］．ＮａｔＧｅｎｅｔ，２０１３，４５（５）：４８７－４９４．
７１
㊀１１期㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀凌士鹏等：基于ＳＳＲ标记的桃种质资源遗传多样性研究
［１３］ＤｏｙｌｅＪＪ，ＤｏｙｌｅＪＬ．ＡｒａｐｉｄＤＮＡｉｓｏｌａｔｉｏｎｐｒｏｃｅｄｕｒｅｆｏｒｓｍａｌｌｑｕａｎｔｉｔｉｅｓｏｆｆｒｅｓｈｌｅａｆｔｉｓｓｕｅ［Ｊ］．ＰｈｙｔｏｃｈｅｍＢｕｌｌ，１９８７（１９）：１１－１５．
［１４］ＤｉｒｌｅｗａｎｇｅｒＥ，ＣｏｓｓｏｎＰ，ＴａｖａｕｄＭ，ｅｔａｌ．Ｄｅｖｅｌｏｐｍｅｎｔｏｆｍｉｃｒｏｓａｔｅｌｌｉｔｅｍａｒｋｅｒｓｉｎｐｅａｃｈ［Ｐｒｕｎｕｓｐｅｒｓｉｃａ（Ｌ．）Ｂａｔｓｃｈ］ａｎｄｔｈｅｉｒｕｓｅｉｎｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙａｎａｌｙｓｉｓｉｎｐｅａｃｈａｎｄｓｗｅｅｔｃｈｅｒｒｙ（ＰｒｕｎｕｓａｖｉｕｍＬ．）［Ｊ］．ＴｈｅｏｒＡｐｐｌＧｅｎｅｔ，２００２，１０５（１）：１２７－１３８．
［１５］ＡｒａｎｚａｎａＭＪ，Ｇａｒｃｉａ－ＭａｓＪ，ＣａｒｂóＪ，ｅｔａｌ．Ｄｅｖｅｌｏｐｍｅｎｔａｎｄｖａｒｉａｂｉｌｉｔｙａｎａｌｙｓｉｓｏｆｍｉｃｒｏｓａｔｅｌｌｉｔｅｍａｒｋｅｒｓｉｎｐｅａｃｈ［Ｊ］．ＰｌａｎｔＢｒｅｅｄｉｎｇ，２００２，１２１（１）：８７－９２．
［１６］ＣｉｐｒｉａｎｉＧ，ＬｏｔＧ，ＨｕａｎｇＷＧ，ｅｔａｌ．ＡＣ／ＧＴａｎｄＡＧ／ＣＴｍｉｃｒｏｓａｔｅｌｌｉｔｅｒｅｐｅａｔｓｉｎｐｅａｃｈ［Ｐｒｕｎｕｓｐｅｒｓｉｃａ（Ｌ）Ｂａｔｓｃｈ］：ｉｓｏｌａｔｉｏｎ，ｃｈａｒａｃｔｅｒｉｓａｔｉｏｎａｎｄｃｒｏｓｓ－ｓｐｅｃｉｅｓａｍｐｌｉfiｃａｔｉｏｎｉｎＰｒｕｎｕｓ［Ｊ］．ＴｈｅｏｒＡｐｐｌＧｅｎｅｔ，１９９９，９９（１／２）：６５－７２．［１７］孙萍，林贤锐，沈建生．基于ＳＳＲ标记的李种质资源遗传多样性研究［Ｊ］．江西农业学报，２０１７，２９（２）：３３－３９．［１８］ＰｅａｋａｌｌＲ，ＳｍｏｕｓｅＰＥ．ＧＥＮＡＬＥＸ６：ｇｅｎｅｔｉｃａｎａｌｙｓｉｓｉｎＥｘ⁃ｃｅｌ．Ｐｏｐｕｌａｔｉｏｎｇｅｎｅｔｉｃｓｏｆｔｗａｒｅｆｏｒｔｅａｃｈｉｎｇａｎｄｒｅｓｅａｒｃｈ［Ｊ］．ＭｏｌｅｃｕｌａｒＥｃｏｌｏｇｙＮｏｔｅｓ，２００６，６（１）：２８８－２９５．［１９］ＬａｎｇｅｌｌａＯ．ＰＯＰＵＬＡＴＩＯＮＳ１．２：ｐｏｐｕｌａｔｉｏｎｇｅｎｅｔｉｃｓｏｆｔｗａｒｅ，ｉｎｄｉｖｉｄｕａｌｓｏｒｐｏｐｕｌａｔｉｏｎｄｉｓｔａｎｃｅ，ｐｈｙｌｏｇｅｎｅｔｉｃ
ｔｒｅｅｓ［ＥＢ／ＯＬ］．ｈｔｔｐ：／／ｂｉｏｉｎｆｏｒｍａｔｉｃｓ．ｏｒｇ／－ｔｒｙｐｈｏｎ／ｐｏｐｕ⁃
ｌａｔｉｏｎｓ／，２０００．［２０］ＨｕｓｏｎＤＨ，ＳｃｏｒｎａｖａｃｃａＣ．Ｄｅｎｄｒｏｓｃｏｐｅ３：ａｎｉｎｔｅｒａｃｔｉｖｅｔｏｏｌｆｏｒｒｏｏｔｅｄｐｈｙｌｏｇｅｎｅｔｉｃｔｒｅｅｓａｎｄｎｅｔｗｏｒｋｓ［Ｊ］．ＳｙｓｔｅｍａｔｉｃＢｉｏｌｏｇｙ，２０１２，６１（６）：１０６１－１０６７．
［２１］史红丽，韩明玉，赵彩平．桃遗传多样性的ＳＲＡＰ和ＳＳＲ标记分析［Ｊ］．华北农学报，２００９，２４（６）：１８７－１９２．［２２］ＧｈａｆｆａｒｉＭＲ．Ａｓｓｅｓｓｍｅｎｔｏｆｔｈｅｇｅｎｅｔｉｃｄｉｖｅｒｓｉｔｙｏｆａｌｍｏｎｄ（Ｐｒｕｎｕｓｄｕｌｃｉｓ）ｕｓｉｎｇｍｉｃｒｏｓａｔｅｌｌｉｔｅｍａｒｋｅｒｓａｎｄｍｏｒｐｈｏｌｏｇｉｃａｌｔｒａｉｔｓ［Ｊ］．ＩｒａｎｉａｎＪｏｕｒｎａｌｏｆＢｉｏｔｅｃｈｎｏｌｏｇｙ，２００８（２）：９８－１０６．
［２３］ＡｒａｎｚａｎａＭＪ，ＣａｒｂóＪ，ＡｒúｓＰ．Ｍｉｃｒｏｓａｔｅｌｌｉｔｅｖａｒｉａｂｉｌｉｔｙｉｎｐｅａｃｈ［Ｐｒｕｎｕｓｐｅｒｓｉｃａ（Ｌ．）Ｂａｔｓｃｈ］：ｃｕｌｔｉｖａｒｉｄｅｎｔｉｆｉｃａｔｉｏｎ，ｍａｒｋｅｒｍｕｔａｔｉｏｎ，ｐｅｄｉｇｒｅｅｉｎｆｅｒｅｎｃｅｓａｎｄｐｏｐｕｌａｔｉｏｎｓｔｒｕｃｔｕｒｅ［Ｊ］．Ｔｈｅｏｒｅｔｉｃａｌ＆ＡｐｐｌｉｅｄＧｅｎｅｔｉｃｓ，２００３，１０６（８）：１３４１－１３５２．
［２４］ＲｕｔｈｎｅｒＳＺ，ＰｅｄｒｙｃＡ，ＫｒｉｓｋａＢ，ｅｔａｌ．Ｍｏｌｅｃｕｌａｒｃｈａｒａｃ⁃ｔｅｒｉｚａｔｉｏｎｏｆａｐｒｉｃｏｔ（ＰｒｕｎｕｓａｒｍｅｎｉａｃａＬ．）ｃｕｌｔｉｖａｒｓｕｓｉｎｇｃｒｏｓｓｓｐｅｃｉｅｓＳＳＲａｍｐｌｉｆｉｃａｔｉｏｎｗｉｔｈｐｅａｃｈｐｒｉｍｅｒｓ［Ｊ］．ＩｎｔＪＨｏｒｔＳｃｉ，２００６（１２）：５３－５７．
［２５］ＷüｎｓｃｈＡ，ＨｏｒｍａｚａＪＩ．Ｍｏｌｅｃｕｌａｒｃｈａｒａｃｔｅｒｉｚａｔｉｏｎｏｆｓｗｅｅｔｃｈｅｒｒｙ（ＰｒｕｎｕｓａｖｉｕｍＬ．）ｇｅｎｏｔｙｐｅｓｕｓｉｎｇｐｅａｃｈ［Ｐｒｕｎｕｓｐｅｒｓｉｃａ（Ｌ．）Ｂａｔｓｃｈ］ＳＳＲｓｅｑｕｅｎｃｅｓ［Ｊ］．Ｈｅｒｅｄｉｔｙ，２００２，８９（１）：５６－６３．
（责任编辑：曾小军）
８１江㊀西㊀农㊀业㊀学㊀报㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀㊀３０卷。