杨梅（Myrica rubra Sieb. et Zucc.）属杨梅科（Myricaceae）杨梅属（Myrica L.）常绿植物，为中国南方特色果树。我国杨梅主要分布于长江流域以南地区，以浙江栽培面积最大，其次为江苏、福建、广东、湖南等地[1-3]。杨梅果实初夏成熟，色泽艳丽，鲜嫩多汁，酸甜适口，富含糖、酸、花青苷、黄酮、维生素和氨基酸等营养物质[4-5]，风味浓郁，既可鲜食，也可加工，深受广大消费者喜爱。
品种资源分析评价是实现植物优良种质创新、品种选育和高效生产的前提和关键[6]。全球杨梅科植物共4 个属约50 余种，主要分布在南美洲、北美洲、欧洲、非洲东部及东亚地区。中国只有杨梅属1属6 个种，包括杨梅（M.rubra）、毛杨梅（M.esculenta）、矮杨梅（M. nana）、青杨梅（M. adenophora）、全缘叶杨梅（M. integrifolia）和大杨梅（M. arboresceus）[1]。我国杨梅种质资源较为丰富，据陈慧等[7]报道我国杨梅有305个品种和105个品系，已定名品种为268 个。但与其他大宗果树相比，我国杨梅育种工作起步晚，种质资源研究相对滞后。关于杨梅品种资源的研究，国内科研工作者主要通过分子标记鉴定[8-9]、基因组学分析[10-11]、芽变选种[12-13]等方法进行杨梅品种的遗传多样性分析和新品种选育研究，但尚未见较为系统全面的表型性状相关研究报道。
杨梅果实、果核、叶片和花序等表型性状不仅是杨梅种质资源描述的重要内容，也是区分杨梅品种资源的主要性状。本研究以浙江省临海市国家杨梅良种繁育基地内收集和引进的95 份杨梅种质为试材，其中包括杨梅种93份、毛杨梅种1份、矮杨梅种和杨梅种杂交1 份。通过对果实、果核、叶片、花序等44 个表型性状进行观测统计，利用相关性分析、主成分分析、聚类分析等方法从形态学水平上研究其遗传多样性，探寻其性状的变异特点，以期为杨梅种质资源研究和育种提供依据和参考。
1 材料和方法
1.1 试验材料
试验于2017—2020年在浙江省临海市小芝镇国家杨梅良种繁育基地进行。试验材料为收集保存的95份杨梅种质（表1），以嫁接方式保存，常规管理。
表1 95 份杨梅种质名称
Table 1 The name of 95 Chinese bayberry germplasms
注：编号1～56 引种自浙江，57～66 引种自福建，67～70 引种自广东，71～82 引种自湖南，83～94 引种自江苏，95 引种自云南。
Note:No.1-No.56 were introduced from Zhejiang province,No.57-No.66 were introduced from Fujian province,No.67-No.70 were introduced from Guangdong province,No.71-No.82 were introduced from Hunan province,No.83-No.94 were introduced from Jiangsu province,and No.95 introduced from Yunnan province.
编号No.1 2 3 4 5 6 7 8 9编号No.25编号No.49种质名称Name of material东魁杨梅Dongkuiyangmei临海早大梅Linhaizaodamei临海早水梅Linhaizaoshuimei临海大黑炭Linhaidaheitan杜桥洋平Duqiaoyangping箬溪变种Ruoxibianzhong临海白梅Linhaibaimei松山早野Songshanzaoye松山糖梅Songshantangmei临海水梅Linhaishuimei松山大野Songshandaye临海大白梅Linhaidabaimei涌泉1号Yongquan1hao涌泉2号Yongquan2hao涌泉3号Yongquan3hao黄岩早梅Huangyanzaomei黄岩大炭梅Huangyandatanmei三门桐子梅Sanmentongzimei黑晶Heijing温岭大梅Wenlingdamei荸荠种杨梅Biqizhongyangmei早荠蜜梅Zaojimimei余姚早酸梅Yuyaozaosuanmei余姚乌大种Yuyaowudazhong编号No.73 26 50 74 27 51 75 28 52 76 29 53 77 30 54 78 31 55 79 32 56 80 33 57 81 10 34 58 82 11 35 59 83 12 36 60 84 13 37 61 85 14 38 62 86 15 39 63 87 16 40 64 88 17 41 65 89 18 42 66 90 19 43 67 91 20 44 68 92 21 45 69 93 22 46 70 94 23 47 71 95种质名称Name of material湖南大叶梅Hunandayemei湖南白杨梅Hunanbaiyangmei上冲梅Shangchongmei墙背梅Qiangbeimei冒顶梅Maodingmei冬瓜梅Dongguamei黑瑞林Heiruilin光贵早Guangguizao王子安海Wangzianhai安海明珠Anhaimingzhu乌梅Wumei桃红Taohong小叶细蒂Xiaoyexidi大叶细蒂Dayexidi甜山Tianshan东方明珠Dongfangmingzhu常熟白杨梅Changshubaiyangmei西山紫条Xishanzitiao张1 Zhang1果1 Guo1王1 Wang1福建晚梅Fujianwanmei毛杨梅Box myrtle 24 48种质名称Name of material余姚水晶Yuyaoshuijing余姚纽扣梅Yuyaoniukoumei余姚粉红种Yuyaofenhongzhong晚荠蜜梅Wanjimimei象山小暑水梅Xiangshanxiaoshushuimei象山本地炭梅Xiangshanbenditanmei乌紫梅Wuzimei奉化野杨梅Fenghuayeyangmei奉化青叶杨梅Fenghuaqingyeyangmei宁海特早香Ninghaitezaoxiang宁海真梅Ninghaizhenmei宁海水晶Ninghaishuijing宁海大麻叶Ninghaidamaye上虞深红种Shangyushenhongzhong上虞水晶Shangyushuijing丁岙梅Dingaomei萧山早色Xiaoshanzaose萧山迟色Xiaoshanchise余杭大炭梅Yuhangdatanmei余杭荔枝梅Yuhanglizhimei金钱炭梅Jinqiantanmei矮杨梅×早色Dwarf bayberry×Zaose兰溪早梅Lanxizaomei兰溪木叶梅Lanximuyemei 72种质名称Name of material兰溪大水梅Lanxidashuimei兰溪杨柳梅Lanxiyangliumei晚稻杨梅Wandaoyangmei定海红杨梅Dinghaihongyangmei大叶岗Dayegang大红袍Dahongpao小夹山Xiaojiashan大夹山Dajiashan硬丝Yingsi安海早Anhaizao胭脂红Yanzhihong安海变种Anhaibianzhong大叶早Dayezao土变Tubian安海中梅Anhaizhongmei福建特早梅Fujiantezaomei安海大粒早Anhaidalizao福建白梅Fujianbaimei广东火炭Guangdonghuotan广东红腊Guangdonghongla广东大虾Guangdongdaxia广东大粒酥Guangdongdalisu湖南小冲梅Hunanxiaochongmei湖南太婆Hunantaipo
1.2 性状观测
每个品种选取3株生长发育良好的结果树作为试验对象。参考《植物新品种特异性、一致性和稳定性测试指南 杨梅》（NY/T 2761—2015）[14]和高志红等[15]的方法，开展果实、果核、叶片、花序等表型性状观测。其中，果实、果核的性状观测主要在5月下旬至6月果实成熟期（即当种质全树约有75%果实的大小、性状和颜色等表现出该种质的固有特性的时间）进行；叶片性状观测主要在11月秋梢生长停止后进行；花序性状观测主要在2月下旬至3月花期进行。
果实、果核：随机选取树冠东南西北中不同部位成熟果实。观测记录果实形状、果实颜色、果实风味、贮藏性、果核形状、果核离核情况、果核缝合线、茸毛颜色和核表颜色，测量果实单果质量、纵径、横径、色泽指标、可溶性固形物含量、可滴定酸含量、果柄长度、果柄粗度、单果核质量、核长度、核宽度、核厚度。为保证测定数据的准确性和一致性，除贮藏性外，其他指标测定均在果实采摘当天完成。
叶片：随机选取老熟春梢中上部的叶片。观测记录叶片形状、叶尖形状、叶缘形状、叶片颜色，测量叶片长度、宽度、厚度。春梢初抽生期，观测记录嫩叶颜色。
花序：随机选取树冠中上部向阳面短果枝中心花枝上中上部花序。选择含苞待放花序测定花序长度、粗度，记录花序形状。开花后观测记录雌花开张度、雌花色泽以及单花序花朵数。
单果质量、单果核质量用准确度为0.01 g 的电子天平测定；果实纵径、果实横径、果柄长度、果柄粗度、果核长度、果核宽度、果核厚度、叶片长度、叶片宽度、花序长度和花序粗度用数显游标卡尺测定；叶片厚度用YH-1叶片厚度计测定；果实可溶性固形物含量用PAL-1 数显糖度计测定；果实可滴定酸含量用PAL-Easy ACID F5数显酸度计测定；果实色泽指标L*、a*、b*用CR-400 色差仪测定。参照Zhang等[16]的方法计算红色葡萄果实颜色指数CIRG 值。果形指数为果实纵径与横径比值。固酸比为果实可溶性固形物含量与可滴定酸含量比值。
1.3 数据分析
观测数据经Excel 2010整理统计。参照曾少敏等[17]方法计算20 个描述型性状的遗传多样性指数Shannon-Wiener指数（H）和Simpson指数（D）。利用SPSS26.0计算24个数量型性状的最小值、最大值、平均值、标准差和变异系数，并进行Pearson相关性分析和主成分分析。通过SPSS26.0以平方欧式距离计算95份种质资源数量型性状遗传距离，利用MEGA X软件按UPGMA遗传距离进行聚类分析并绘制聚类图。
2 结果与分析
2.1 描述型性状的多样性分析
将杨梅种质的20 个描述型性状进行分类（表2），其中果实颜色、果核核表颜色等分组类型较多，具体各性状的分布情况如表3所示。
表2 杨梅种质描述型性状及分级标准
Table 2 Descriptive phenotypic traits and grading standards of Chinese bayberry germplasms
描述型性状Descriptive trait果实形状Fruit shape(FS)果实颜色Fruit color(FC)果肉质地Flesh quality(FQ)果汁量Fruit juice(FJ)果实风味Fruit flavor(FF)果实香气Fruit aroma(FA)果实贮藏性Fruit storability(FST)果核形状Stone shape(SS)果核离核情况Stone adherence to flesh(SA)果核缝合线Stone sutural line(SSL)果核茸毛颜色Stone fur color(SF)果核核表颜色Stone color(SC)叶形Leaf shape(LS)叶尖形状Leaf apex(LA)叶色Leaf color(LC)叶缘Leaf margin(LM)嫩叶颜色Young leaf color(YLC)花序形状Inflorescence shape(IS)雌花开张度Female flower angle(FFA)雌花色泽Female flower color(FFC)分级标准Description of grading 1高圆球形Highly round紫黑Purplish black较硬Firm多Much甜酸Sweet-sour有香味Fragrance较好Good近圆形Suborbicular较离核Free明显Evident淡黄褐色Light tawny淡黄绿色Light yellowish-green窄倒披针形Narrow Oblanceolate急尖Acute深绿色Dark green全缘Entire淡紫红色Light purplish red圆筒Cylinder V形V-shape紫红色Purplish red 2 3 4 5圆球形Round紫红Purplish red中Medium中Medium酸甜Sour-sweet无Absent中Medium卵形或扁卵形Oval or flat oval中Semi-free稍明显Some obvious黄褐色Tawny黄绿色Yellowish-green倒披针形Oblanceolate渐尖Acuminate绿色Green浅锯齿Shallow serration橘黄色Orang长圆筒Long cylinder倒人字形Inverted“人”shape玫瑰红色Rose red扁圆球形Oblate深红Dark red较柔软Softer少Little淡甜酸Light sweet-sour淡松脂味Light turpentine odor较差worse长卵形Long oval较黏核Cling不明显Nothing浅棕色Light brown青绿色Bluish green窄倒卵形Narrow Obovate钝Blunt黄绿色Yellowish green深锯齿Tip entire淡绿色Light green短圆筒Short cylinder M形M-shape红色Red红色Red水红色或白色Water red or white较酸Sour浓松脂味Dense turpentine odor椭圆形Elliptic棕色Brown黄褐色Tawny倒卵形Obovate微缺Intaglio灰褐色Taupe带波状Wave-like褐红色Brown red淡粉色Pink
表3 杨梅种质描述型性状频率分布及多样性
Table 3 Frequency distribution and diversity of phenotype traits in Chinese bayberry germplasms
注：表中数字代表的描述型性状详见表2。Note:The numbers resprsented descriptive traist are shown in Table 2.
描述型性状Descriptive trait果实形状Fruit shape(FS)果实颜色Fruit color(FC)果肉质地Flesh quality(FQ)果汁量Fruit juice(FJ)果实风味Fruit flavor(FF)果实香气Fruit aroma(FA)果实贮藏性Fruit storability(FST)果核形状Stone shape(SS)果核离核情况Stone adherence to flesh(SA)果核缝合线Stone sutural line(SSL)果核茸毛颜色Stone fur color(SF)果核核表颜色Stone color(SC)叶形Leaf shape(LS)叶尖形状Leaf apex(LA)叶色Leaf color(LC)叶缘Leaf margin(LM)嫩叶颜色Young leaf color(YLC)花序形状Inflorescence shape(IS)雌花开张度Female flower angle(FFA)雌花色泽Female flower color(FFC)各级频率Frequency of classification/%1 12.63 18.95 35.79 38.95 65.26 6.32 54.74 10.53 22.10 43.16 49.47 26.31 3.16 8.42 36.84 84.21 26.32 68.42 34.74 55.79 Simpson指数Simpson index,D 0.422 4 0.728 8 0.641 3 0.534 3 0.520 4 0.482 7 0.581 9 0.493 1 0.615 8 0.644 2 0.662 2 0.766 0 0.603 4 0.473 3 0.515 0 0.271 9 0.560 0 0.482 0 0.527 8 0.615 0 2 3 4 5 73.68 42.11 21.05 55.79 21.05 69.47 31.58 68.42 26.32 22.10 16.84 26.32 55.79 70.52 58.95 13.69 60.00 15.79 58.95 7.37 13.69 20.00 43.16 5.26 6.32 14.74 13.68 15.79 51.58 7.37 23.16 27.37 20.00 10.53 4.21 1.05 9.47 15.79 6.31 20.00 9.47 9.47 7.37 9.47 5.26 10.53 11.58 21.05 10.53 8.42 1.05 4.21 16.84 Shannon-Wiener指数Shannon-Wiener index,H 0.758 6 1.447 7 1.058 4 0.847 7 0.973 3 0.933 0 0.966 0 0.943 0 1.026 4 1.063 6 1.224 0 1.515 3 1.084 6 0.928 6 0.812 8 0.512 5 1.014 4 0.842 6 0.853 2 1.139 6
杨梅种质果实形状包括高圆球形、圆球形和扁圆球形，以圆球形为主（73.68%）。杨梅果实颜色可分为紫黑色、紫红色、深红色、红色、水红色或白色，以紫红色为主（42.11%），其次为深红色（20.00%）和紫黑色（18.95%），而红色占9.47%，水红色或白色只占9.47%。杨梅果实肉质可分为较硬、中、较柔软，以较柔软为多。果实果汁量较多，只有5份种质果汁少。果实风味以甜酸风味为主（65.26%）。杨梅果实大部分没有香味（69.47%），少数有清香味（6.32%），部分有松脂味；其中9份有较浓的松脂味，14份有淡松脂味。不同种质贮藏性存在差异，按贮藏性分为较好、中和较差，比例分别为54.74%、31.58%和13.68%。
杨梅种质果核形状可分为近圆形、卵形或扁卵形、长卵形和椭圆形，以卵形或扁卵形为主（68.42%）。杨梅果核离核情况可分为较离核、中和较黏核，以较黏核为主。大部分果核缝合线较明显。果核茸毛颜色包括淡黄褐色、黄褐色、浅棕色和棕色，以淡黄褐色比例最高（49.47%）。果核核表颜色包括淡黄绿色、黄绿色、青绿色、黄褐色和灰褐色。
杨梅种质叶片形状包括窄倒披针形、倒披针形、窄倒卵形和倒卵形，倒披针形比例最高（55.79%）。叶尖包括急尖、渐尖、钝和微缺，以渐尖为主（70.52%）。叶色包括深绿色、绿色和黄绿色，以绿色为主（58.95%）。叶缘有全缘、浅锯齿、深锯齿和带波状，以全缘为主（84.21%）。嫩叶颜色包括淡紫红色、橘黄色、淡绿色和褐红色，橘黄色比例最高（60.00%）。
杨梅种质雌花花序形状可分为圆筒形、长圆筒形和短圆筒形3类，以圆筒形为主（68.42%）。雌花开张度可分为V形、倒人字形、M形3类，以倒人字形比例最高（58.95%）。雌花色泽包括紫红色、玫瑰红色、红色和淡粉色，以紫红色比例最高（55.79%）。
由表3 可知，20 个描述型性状的Shannon-Wiener 指数（H）和Simpson 指数（D）变化范围分别为0.512 5～1.515 3 和0.271 9～0.766 0，其中果实颜色、果核茸毛颜色和果核核表颜色的Shannon-Wiener指数（H）分别为1.447 7、1.224 0 和1.515 3，数值较高，说明多样性较为丰富。
2.2 数量型性状多样性分析
从表4可知，95份杨梅种质资源的24个数量型性状变异系数范围为1.28%～63.99%。仅果实纵径、横径、果形指数、可溶性固形物含量、可食率和果核厚度6 个性状的变异系数小于10%，其他18 个性状的变异系数均大于10%，说明这些性状个体间存在较大差异。其中，变异系数以果实色泽指标b*值最高（63.99%），其次果柄长度为(33.84%)，变异系数排序为b*值＞果柄长度＞L*＞a*＞单果核质量。
表4 杨梅种质24 个数量型性状的变异情况
Table 4 Coefficient of varation on 24 quantitative traits of Chinese bayberry germplasms
性状Trait果实性状Fruit trait果核性状Stone trait叶片性状Leaf trait花序性状Inflorescence trait单果质量Single fruit mass(FM)/g果实纵径Fruit vertical diameter(FVD)/mm果实横径Fruit transverse diameter(FTD)/mm果形指数Fruit shape index(FSI)果柄长度Fruit stalk length(FSL)/mm果柄粗度Fruit stalk thickness(FST)/mm L*a*b*CIRG w（可溶性固形物）Total soluble solids content(TSS)/%w（可滴定酸）Titratable acid content(TA)/%固酸比Solid acid ratio(SAR)可食率Edible rate(ER)/%单果核质量Single stone mass(SM)/g果核长度Stone length(SL)/mm果核宽度Stone width(SW)/mm果核厚度Stone thickness(ST)/mm叶片长度Leaf length(LL)/cm叶片宽度Leaf width(LW)/cm叶片厚度Leaf thickness(LT)/mm花序长度Inflorescence length(IL)/mm花序粗度Inflorescence thickness(IT)/mm单花序花朵数Number of flowers per inflorescence(NFI)最小值Min 5.01 21.51 21.17 0.86 2.92 0.62 16.22 6.30 1.05 1.78 8.62 0.74 4.74 90.92 0.34 10.29 7.70 6.27 6.23 2.09 0.28 4.36 1.62 3.20最大值Max 25.34 36.15 36.42 1.04 18.87 2.17 42.99 24.68 21.14 6.35 13.40 2.11 16.21 96.50 1.51 18.46 13.34 10.97 12.37 4.13 0.58 19.60 2.96 10.40平均值Mean 12.47 27.68 28.54 0.97 7.80 1.55 20.75 15.03 7.08 4.40 11.37 1.14 10.37 94.32 0.69 13.08 9.96 8.20 9.49 2.91 0.44 8.69 2.24 6.78标准差s 3.30 2.63 2.62 0.03 2.64 0.26 6.38 4.33 4.53 1.11 0.68 0.23 1.94 1.20 0.20 1.64 1.02 0.78 1.22 0.47 0.05 2.32 0.29 1.47变异系数CV/%26.46 9.50 9.58 3.12 33.84 16.64 30.76 28.82 63.99 25.23 6.01 20.52 18.70 1.28 28.34 12.50 10.27 9.55 12.88 16.26 11.06 26.71 12.88 21.76
果实性状：果实单果质量、纵径、横径的均值分别为12.47 g、27.68 mm、28.54 mm，其中单果质量的变异系数为26.46%。种质箬溪变种单果质量最大（25.34 g）、纵径最长（36.15 mm）、横径最长（36.42 mm）。毛杨梅单果质量最小（5.01 g）、纵径最短（21.51 mm）。单果质量14 g以上的种质比例为24.21%，10 g以上种质的比例为85.26%。果形指数最大为广东大虾（1.04），最小为湖南白杨梅（0.86）。果柄最长为丁岙梅（18.87 mm），最短为甜山（2.92 mm），最粗为东方明珠（2.17 mm），最细为奉化野杨梅（0.62 mm）。
果实色泽指标L*、a*、b*、CIRG 均值分别为20.75、15.03、7.08、4.40，变异系数最大的为b*。明亮度L*值最大为宁海水晶（42.99），最小为安海中梅（16.22）。红绿色度a*值最大为余姚早酸梅（24.68），最小为临海白梅（6.30）。黄蓝色度b*值最大为临海大白梅（21.14），最小为西山紫条（1.05）。色泽指标CIRG值最大为西山紫条（6.35），最小为临海大白梅（1.78）。
果实可溶性固形物含量、可滴定酸含量、固酸比和可食率均值分别为11.37%、1.14%、10.37%和94.32%，变异系数最大为可滴定酸含量（20.52%）。果实可溶性固形物含量最高为晚稻杨梅（13.40%），最低为毛杨梅（8.62%）。可滴定酸含量最高为广东红腊（2.11%），最低为安海变种（0.74%）。固酸比最高为上冲梅（16.21），最低为毛杨梅（4.74）。可食率最高为东魁杨梅（96.50%），最低为胭脂红（90.92%）。
果核性状：杨梅种质单果核质量、果核长度、果核宽度和果核厚度均值分别为0.69 g、13.08 mm、9.96 mm和8.20 mm，变异系数最大的为单果核质量（28.34%）。单果核质量最大为箬溪变种（1.51 g），最小为余姚纽扣梅（0.34 g）。果核长度最长为箬溪变种（18.46 mm），最短为安海早（10.29 mm）。果核宽度最宽为乌紫梅（13.34 mm），最窄为余姚纽扣梅（7.70 mm）。果核厚度最厚为乌紫梅（10.97 mm），最薄为桃红（6.27 mm）。
叶片性状：杨梅种质叶片长度、宽度、厚度均值分别为9.49 cm、2.91 cm 和0.44 mm，其中叶片宽度变异系数较大（16.26%）。叶片长度最长为甜山（12.37 cm），最短为福建晚梅（6.23 cm）。叶片宽度最宽为临海水梅（4.13 cm），最窄为广东火炭（2.09 cm）。叶片厚度最厚为东方明珠（0.58 mm），最薄为毛杨梅雌株（0.28 mm）。
花序性状：95 份杨梅种质雌性花序长度、粗度和单花序花朵数均值分别为8.69 mm、2.24 mm 和6.78 朵。其中花序长度最长为丁岙梅（19.60 mm），最短为晚稻杨梅（4.36 mm）。花序粗度最粗为涌泉2 号（2.96 mm），最细为晚稻杨梅（1.62 mm）。单花序花朵数最多为湖南大叶梅（10.4 朵），花朵数最少为晚稻杨梅（3.20朵）。
2.3 数量型性状的相关性分析
对95份杨梅种质的24个数量型性状进行Pearson 相关性分析，结果见表5，其中72 对相关系数达极显著水平（p＜0.01），23 对相关系数达显著水平（p＜0.05）。单果质量、纵径和横径三者之间互相呈极显著正相关，相关系数分别为0.962、0.941、0.949。单果质量、纵径、横径均与果柄粗度、可食率、单果核质量、果核长度、果核宽度、果核厚度、叶片长度、叶片宽度、叶片厚度、花序粗度呈极显著正相关。单果质量与果实可滴定酸含量呈极显著负相关，与固酸比呈显著正相关。果实纵径与果实可滴定酸含量呈显著负相关。横径与果形指数、固酸比呈显著正相关，与果实可滴定酸含量呈极显著负相关。果形指数与果实可滴定酸含量、果核长度呈极显著正相关，与可食率呈极显著负相关，与L*呈显著负相关。果柄长度与果核宽度、花序长度、单花序花朵数呈极显著正相关。果柄粗度与果核长度、叶片厚度呈极显著正相关，与单果核质量呈显著正相关。
表5 杨梅质资源数量型性状相关性分析
Table 5 Correlation analysis of quantitative traits of Chinese bayberry germplasms
NFI 1 IT 10.077 IL 10.156 0.441**LT 10.142 0.225*0.130 LW 10.407**LL 10.638**0.241*0.445**0.176 0.207*0.213*0.361**0.471**0.157 ST 10.162 0.186 0.185 0.404**SW 10.675**0.139 0.236*0.373**SL 10.649**0.262*0.613**0.344**0.171 0.290**0.248*0.267**SM 10.772**0.359**0.632**0.386**0.342**0.418**0.212*0.219*0.201 0.180 0.124 0.126 ER 1-0.408**-0.199 0.717**-0.184-0.149 0.036 0.048 0.185 0.061 0.258*-0.145 SAR 10.375**-0.073-0.077-0.082 0.025 0.046-0.020 0.154-0.040-0.216*0.048 TA 1-0.918**-0.420**0.040 0.068 0.148-0.014-0.010 0.025-0.188 0.060 0.173 0.046 0.221*0.050 0.044 0.042 0.067 0.037 0.043 0.147 0.011 0.187 TSS 10.045-0.104-0.053 0.222*0.181 0.048 0.051 CIRG 1-0.097-0.186-0.152-0.062-0.036-0.015-0.134-0.068-0.042-0.007 b* 1-0.970**0.040 0.186-0.230*-0.165 0.157 0.064 0.077-0.021 0.005 0.036 0.137 0.088 0.035 0.024 a* 10.249*-0.349**0.214*0.135-0.093-0.119 0.048 0.135-0.038-0.065-0.100-0.079 0.054 0.062 0.075 0.044 L* 1-0.175 0.876**-0.839**-0.027 0.074-0.144-0.117 0.089-0.061 0.008-0.048-0.016 0.001 0.085-0.009-0.052-0.044 FST 1-0.015 0.033 0.028-0.006 0.023 0.072-0.083 0.118 0.252*0.299**0.008 0.117 0.195 0.081 0.329**-0.094 0.178 0.006 FSL 10.069 0.025 0.135 0.138-0.120 0.171 0.080 0.008-0.011 0.133 0.167 0.291**0.089 0.151 0.130 0.140 0.722**0.114 0.309**FSI 10.087 0.261*-0.249*0.189-0.109 0.116 0.164 0.284**-0.187-0.305**0.185 0.418**0.114 0.045 0.100-0.032-0.073-0.057-0.018 0.182 FTD 1-0.225*0.117 0.275**0.024 0.040 0.243*0.388**0.585**0.546**0.474**0.595**0.329**0.418**0.531**0.201 0.319**0.037 FVD 10.949**-0.166-0.018-0.125-0.319**0.092 0.144 0.366**0.083 0.191 0.297**0.658**0.696**0.521**0.623**0.366**0.421**0.515**0.186 0.320**0.101 FM 10.962**0.941**-0.064-0.104-0.062-0.070-0.237*0.006 0.131 0.374**-0.013-0.085 0.001 0.028-0.087-0.278**0.220*0.309**0.684**0.674**0.512**0.617**0.360**0.424**0.534**0.169 0.293**0.068性Trait状FMFVD FTD FSI FSL FST L*a*b*CIRG TSS TASAR ERSMSLSWSTLLLWLTILITNFI。关相著显平0.01 水p＜0.05 和p＜在示表别**分*和4。表见详状性型量数的表代母字中：表注Note:The character respresented quantitative trait is shown in Table 4.*and**indicated significant correlation at p＜0.05 and p＜0.01,respectively.
L*与b*呈极显著正相关，与CIRG呈极显著负相关。a*与b*、可溶性固形物含量呈显著正相关，与CIRG呈极显著负相关。b*与CIRG呈极显著负相关，与固酸比呈显著负相关。CIRG与固酸比呈显著正相关。果实可溶性固形物含量与a*、固酸比呈显著正相关。果实可滴定酸含量与固酸比、可食率呈极显著负相关。固酸比与可食率呈极显著正相关。可食率和单果核质量呈极显著负相关，与花序粗度呈显著正相关。
单果核质量、果核长度、果核宽度、果核厚度之间均呈极显著正相关。单果核质量、果核长度均与叶片长度、叶片宽度、叶片厚度呈极显著正相关。果核长度与花序长度、花序粗度呈显著正相关。果核宽度与叶片宽度、花序长度呈极显著正相关，与叶片长度、叶片厚度、单花序花朵数呈显著正相关。果核厚度与叶片宽度、叶片厚度呈极显著正相关。叶片长度与叶片宽度、单花序花朵数呈极显著正相关，与叶片厚度、花序粗度呈显著正相关。叶片宽度与叶片厚度、花序粗度呈极显著正相关，与花序长度呈显著正相关。叶片厚度与花序粗度呈显著正相关。花序长度与单花序花朵数呈极显著正相关。
2.4 主成分分析
对杨梅种质资源的24 个数量性状进行主成分分析，结果见表6。根据特征值大于1的原则提取前8 个主成分，累计贡献率达到81.827%，基本可以反映原始因子所代表的大部分信息。其中，主成分1贡献率为26.187%，根据特征向量绝对值大小得出，起决定作用的有单果质量、果实纵径和果实横径等性状，主要反映果实大小因子；主成分2 贡献率为14.167%，起决定作用的有果实色泽指标CIRG、b*、L*、可滴定酸含量和固酸比等性状，主要反映果实色泽因子和糖酸因子；主成分3 贡献率为10.651%，起决定作用的有L*、b*、CIRG 和果形指数等性状，主要反映果实色泽因子；主成分4 贡献率为8.223%，起决定作用的有花序长度和果柄长度等性状，主要反映果实果柄长短因子；主成分5 贡献率为6.963%，起决定作用的有花序粗度、固酸比和可溶性固形物含量等性状，主要反映花序粗度和果实糖酸因子；主成分6 贡献率为6.169%，起决定作用的有果梗粗度、a*和可溶性固形物含量等性状，主要反映果梗粗度、果实色泽和果实糖度因子等性状；主成分7 贡献率为5.12%，起决定作用的有叶片长度和叶片宽度等性状，主要反映叶片因子；主成分8 贡献率为4.349%，起决定作用的有果柄粗度和a*等性状。
表6 杨梅种质数量型性状的主成分分析
Table 6 Principal component analysis of quantitative traits of Chinese bayberry germplasms
注：下划线指示特征值较大的性状。
Note:The underline indicates the traits with higher eigenvalue.
性状Trait单果质量Single fruit mass(FM)果实纵径Fruit vertical diameter(FVD)果实横径Fruit transverse diameter(FTD)果形指数Fruit shape index(FSI)果柄长度Fruit stalk length(FSL)果柄粗度Fruit stalk thickness(FST)L*a*b*CIRG(CIRG)可溶性固形物含量Total soluble solids content(TSS)可滴定酸含量Titratable acid content(TA)固酸比Solid acid ratio(SAR)可食率Edible rate(ER)单果核质量Single stone mass(SM)果核长度Stone length(SL)果核宽度Stone width(SW)果核厚度Stone thickness(ST)叶片长度Leaf length(LL)叶片宽度Leaf width(LW)叶片厚度Leaf thickness(LT)花序长度Inflorescence length(IL)花序粗度Inflorescence thickness(IT)单花序花朵数Number of flowers per inflorescence(NFI)特征值Eigenvalue方差贡献率Contributive percentage/%累计方差贡献率Total percentage/%主成分Principle component PC1 0.907 0.911 0.860 0.098 0.284 0.353-0.005-0.035 0.055-0.026 0.029-0.116 0.101 0.092 0.793 PC2-0.211-0.224-0.279 0.204 0.249 0.006 0.602 0.339 0.776 PC3-0.173-0.088-0.292 0.634 PC4-0.150-0.139-0.094-0.134 0.692 PC6 0.056 0.064-0.054 0.364 0.023 0.604 PC8 0.097 0.113 0.021 0.304 0.203 0.439 0.100 0.039-0.701-0.157 0.542 0.111-0.590-0.256-0.076 0.125 0.006-0.006 0.299-0.171 0.297 0.198-0.266-0.189 0.027-0.184 0.246 0.158 0.034 0.737 PC5 0.007-0.015-0.084 0.209-0.001-0.102 0.009 0.243 0.060-0.103 0.467 PC7-0.112-0.116-0.135 0.063-0.338 0.013 0.159-0.297 0.023 0.009 0.261-0.064 0.140-0.194 0.025-0.092-0.127-0.054 0.630 0.209-0.391-0.780 0.578 0.134 0.661-0.637-0.417 0.535 0.792 0.698 0.758 0.485 0.608 0.622 0.338 0.383 0.266-0.588 0.263 0.194 0.220 0.037-0.002 0.044-0.023 0.188 0.053 0.178 0.121 0.461-0.337-0.439 0.128 0.270 0.213 0.121 0.092 0.017-0.252 0.087 0.004 0.266 0.058 0.564-0.336 0.269 0.185 0.094 0.179-0.258-0.348 0.072-0.145-0.611 0.072-0.130 0.344 0.010 0.049 0.307-0.155 0.041-0.368-0.266 0.065-0.034 0.256-0.193 0.274-0.041 0.462 0.098 0.069-0.323-0.056 0.269 0.057 0.028-0.327-0.056-0.002 0.014-0.012 0.072-0.106-0.249 0.156-0.290-0.178 0.118-0.338 0.193 1.044 4.349 81.827 6.285 26.187 26.187 3.400 14.167 40.354 2.556 10.651 51.005 1.973 8.223 59.227 1.671 6.963 66.190 1.480 6.169 72.359 1.229 5.120 77.479
2.5 聚类分析
如图1 所示，基于数量型性状的UPGMA 聚类，在遗传距离为20 处可将95 份杨梅种质资源划分为5 大类群，各类群性状数值分布情况参见表7。其中，类群Ⅰ只有毛杨梅1个品种，该种质在所有参试种质中，单果质量最小，果实纵径最小，可溶性固形物含量最低，固酸比最小，叶片厚度最薄，可作为特殊种质资源利用。类群Ⅱ包括乌紫梅、箬溪变种、东方明珠、东魁杨梅4 份种质，主要特征表现为大果型，平均单果质量均大于20 g，果实纵横径均大于33 mm，果实色泽为紫红色或深红色，L*值为17.78～19.68，a*值为15.10～22.33，b*值为4.34～9.33，CIRG值为3.69～4.90，可溶性固形物含量为11.06%～12.25%，可滴定酸含量为0.84%～1.19%，固酸比为9.81～13.45，该类可用于选育优质大果型品种。类群Ⅲ包括余姚纽扣梅、余姚水晶和奉化野杨梅3 份种质，主要特征表现为小果型种质，单果质量为6.47～8.67 g，果实颜色为红色或水红色，L*值为21.43～33.29，a*值 为20.10～22.42，b*值 为10.36～17.50，CIRG值为2.32～3.34，可滴定酸含量为1.16%～1.51%，固酸比为6.76～9.65。类群Ⅳ包括福建白梅、宁海水晶、临海白梅、常熟白杨梅、上虞水晶、胭脂红、临海大白梅和湖南白杨梅8份种质，主要特征表现为白梅类种质，果实色泽为水红色或白色，色泽指标L*值为35.41～42.99，a*值为6.30～17.39，b*值为16.27～21.14，CIRG 值为1.78～2.21，可溶性固形物含量为10.61%～11.62%，可滴定酸含量为1.00%～1.28%，固酸比为8.93～11.59，单果质量为10.40～18.40 g，该类可用于选育优质白杨梅品种。类群Ⅴ包含79 份种质，占供试材料的83.16%，单果质量为7.38～17.91 g，果实色泽包括紫黑、紫红、深红或红色，L*值为16.22～27.14，a*值为8.26～24.68，b*值为1.05～14.96，CIRG 值为2.66～6.35，可溶性固形物含量为10.10%～13.40%，可滴定酸含量为0.74%～2.11%。
图1 95 份杨梅种质基于24 个数量型性状的UPGMA 聚类分析
Fig.1 UPGMA cluster analysis of 95 Chinese bayberry germplasms based on 24 quantitative traits
表7 杨梅种质不同类群数量型性状分布情况
Table 7 The variation on 24 quantitative traits in different groups of Chinese bayberry germplasms
性状Trait单果质量Single fruit mass(FM)/g果实纵径Fruit vertical diameter(FVD)/mm果实横径Fruit transverse diameter(FTD)/mm果形指数Fruit shape index(FSI)果柄长度Fruit stalk length(FSL)/mm果柄粗度Fruit stalk thickness(FST)/mm L*a*b*CIRG w（可溶性固形物）Total soluble solids content(TSS)/%w（可滴定酸）Titratable acid content(TA)/%固酸比Solid acid ratio(SAR)可食率Edible rate(ER)/%单果核质量Single stone mass(SM)/g果核长度Stone length(SL)/mm果核宽度Stone width(SW)/mm果核厚度Stone thickness(ST)/mm叶片长度Leaf length(LL)/cm叶片宽度Leaf width(LW)/cm叶片厚度Leaf thickness(LT)/mm花序长度Inflorescence length(IL)/mm花序粗度Inflorescence thickness(IT)/mm单花序花朵数Number of flowers per inflorescence(NFI)Ⅰ类GroupⅠ(n=1)5.01 21.51 21.44 1.00 4.61 1.24 17.22 12.61 4.41 5.26 8.62 1.82 4.74 91.60 0.42 12.58 10.73 7.03 9.85 2.38 0.28 6.80 2.10 6.50Ⅱ类Group Ⅱ(n=4)20.24～25.34 33.47～36.15 34.50～36.42 0.96～1.00 6.59～11.09 1.71～2.17 17.78～19.68 15.10～22.33 4.34～9.33 3.69～4.90 11.06～12.25 0.84～1.19 9.81～13.45 93.77～96.50 0.86～1.51 15.28～18.46 10.94～13.34 9.06～10.97 9.60～11.50 3.20～3.78 0.45～0.58 6.71～9.73 2.40～2.82 5.20～8.70Ⅲ类Group Ⅲ(n=3)6.47～8.67 22.40～23.60 22.43～23.91 0.97～1.00 3.29～9.39 0.62～2.04 21.43～33.29 20.10～22.42 10.36～17.50 2.32～3.34 10.21～11.68 1.16～1.51 6.76～9.65 91.90～94.70 0.34～0.58 11.23～11.46 7.70～9.45 6.66～8.26 7.11～9.42 2.13～2.49 0.38～0.40 5.36～8.32 1.68～2.18 5.50～5.80Ⅳ类Group Ⅳ(n=8)10.40～18.40 25.10～31.87 27.76～33.14 0.86～1.00 4.85～11.13 1.23～1.80 35.41～42.99 6.30～17.39 16.27～21.14 1.78～2.21 10.61～11.62 1.00～1.28 8.93～11.59 90.92～95.40 0.59～1.11 11.53～15.20 9.10～12.26 7.57～9.73 8.85～10.43 2.64～3.32 0.43～0.47 5.82～13.02 1.82～2.57 4.50～9.60Ⅴ类Group Ⅴ(n=79)7.38～17.91 22.01～31.47 21.17～32.57 0.90～1.04 2.92～18.87 1.05～2.09 16.22～27.14 8.26～24.68 1.05～14.96 2.66～6.35 10.10～13.40 0.74～2.11 5.55～16.21 90.94～96.45 0.39～1.12 10.29～16.14 8.20～11.89 6.27～10.46 6.23～12.37 2.09～4.13 0.34～0.56 4.36～19.60 1.62～2.96 3.20～10.40
3 讨 论
表型性状变异是生物遗传多样性在形态水平上的表现，反映了植物自身遗传因素与外部生活环境相互作用的结果[18-19]。表型性状观测是研究遗传多样性简便易行的方法[20]。通过表型性状遗传多样性研究，不仅能从整体了解表型遗传多样性的丰富程度，更为分子生物学研究提供可靠的表型数据[21]。在本研究中通过对95 份杨梅种质资源44 个表型性状的多样性分析，发现杨梅种质资源的表型性状具有丰富的遗传变异。20 个描述型性状的遗传多样性指数Shannon-Wiener 指数（H）和Simpson 指数（D）变化范围分别为0.512 5～1.515 3 和0.271 9～0.766 0，其中核表颜色、果实颜色、果肉质地、果核离核情况等多样性指数较高。24 个数量型性状的变异系数为1.28%～63.99%，其中果实色泽指标、单果质量、单果核质量、花序长度、果柄长度和可滴定酸含量等数量型性状表现出较丰富的变异性。
杨梅种质描述型性状的分布频率差异较大，其中占比较大的种质表现为果实圆球形，紫红色，果实甜酸，果核卵形，较黏核，茸毛淡黄褐色，叶片倒披针形，叶尖渐尖，叶色绿色，叶缘全缘，嫩叶橘黄色，圆筒形花序，雌花紫红色，柱头开张度为倒人字形。
在种质筛选过程中，需注重各农艺性状之间的相关性，综合各性状特点兼具多种优良性状，可为品种选育奠定基础[22]。笔者在本试验中通过数量性状的相关性分析，发现杨梅种质单果质量、果实纵径、果实横径与单果核质量、果核长度、果核宽度、叶片长度、叶片宽度、叶片厚度之间存在极显著正相关性，说明杨梅果实的质量、大小与果核的质量、大小有关，与叶片的大小也存在相关性。果柄长度与花序长度呈极显著正相关，说明花序较长的品种，其果实果柄较长。果实色泽指标a*值与果实可溶性固形物含量呈显著正相关，说明在一定范围内，杨梅果实颜色红色程度越深，果实可溶性固形物含量越高，这与邱立军等[23]研究结果一致。通过数量性状的主成分分析，根据贡献率的大小，从主成分中筛选出单果质量、果实纵横径、果实色泽指标CIRG、b*、L*及可滴定酸含量、固酸比等性状是形成杨梅种质资源表型差异的主要因素。通过数量性状的聚类分析，杨梅种质被划分为5 个类群，其中类群Ⅰ只有毛杨梅1 个品种。毛杨梅属于杨梅属毛杨梅种，而其他种质属于杨梅种及矮杨梅种和杨梅种的杂交，因此在性状表现上差异较大。早荠蜜梅和晚荠蜜梅是荸荠种中芽变选育出的2个品种[24-25]，在本文聚类分析图中，3个品种聚类在同一小组，说明表型性状的聚类一定程度上也能反映出种质间的亲缘关系。
4 结 论
杨梅种质资源遗传变异丰富，表型性状多样性高。通过主成分分析，根据贡献率的大小，从主成分中筛选出果实大小因子、外观色泽因子、糖酸含量因子等可用于作为对杨梅品质综合评定的重要指标。聚类分析在欧式距离20处可将杨梅种质资源分为5类，其中类群Ⅱ可用于选育大果型优质品种，类群Ⅳ可用于选育优质白梅类品种。
[1] 何新华，陈力耕，陈怡，郭长禄.中国杨梅资源及利用研究评述[J].果树学报，2004，21(5)：467-471.HE Xinhua，CHEN Ligeng，CHEN Yi，GUO Changlu. Review on germplasm resources of Myrica and their exploitation in China[J].Journal of Fruit Science，2004，21(5)：467-471.
[2] CHEN K S，XU C J，ZHANG B，FERGUSON I B. Red bayberry：botany and horticulture[M]//JANICK J. Horticultural Reviews.Oxford，UK：John Wiley&Sons，Inc.，2010：83-114.
[3] 周超超，陈竹韵，汪国云，朱奕凡，巨鹏举，赵岚，陈金辉，焦云，高中山.余姚杨梅种质资源研究和开发利用[J].果树资源学报，2022，3(4)：1-6.ZHOU Chaochao，CHEN Zhuyun，WANG Guoyun，ZHU Yifan，JU Pengju，ZHAO Lan，CHEN Jinhui，JIAO Yun，GAO Zhongshan.Research and utilization of the red bayberry germplasm resources in Yuyao city[J].Journal of Fruit Resources，2022，3(4)：1-6.
[4] 张泽煌，钟秋珍，林旗华.杨梅果实氨基酸组成及营养评价[J].热带作物学报，2012，33(12)：2279-2283.ZHANG Zehuang，ZHONG Qiuzhen，LIN Qihua. The content of amino acid of Chinese bayberry fruit and its nutritive evaluation[J]. Chinese Journal of Tropical Crops，2012，33(12)：2279-2283.
[5] 李伟，郜海燕，陈杭君，吴伟杰，房祥军.基于主成分分析的不同品种杨梅果实综合品质评价[J]. 中国食品学报，2017，17(6)：161-171.LI Wei，GAO Haiyan，CHEN Hangjun，WU Weijie，FANG Xiangjun.Evaluation of comprehensive quality of different varieties of bayberry based on principal components analysis[J].Journal of Chinese Institute of Food Science and Technology，2017，17(6)：161-171.
[6] 钟声远，罗宇婷，赵勇，王振兴，管志勇，房伟民，陈发棣，王海滨.切花菊品种资源表型多样性分析[J].植物资源与环境学报，2021，30(5)：22-33.ZHONG Shengyuan，LUO Yuting，ZHAO Yong，WANG Zhenxing，GUAN Zhiyong，FANG Weimin，CHEN Fadi，WANG Haibin.Analysis on phenotypic diversity of cut chrysanthemum cultivar resources[J]. Journal of Plant Resources and Environment，2021，30(5)：22-33.
[7] 陈慧，唐威，费艳.杨梅种质资源遗传多样性研究进展[J].现代园艺，2016(3)：5-7.CHEN Hui，TANG Wei，FEI Yan. Research progress on genetic diversity of Myrica rubra germplasm resources[J]. Contemporary Horticulture，2016(3)：5-7.
[8] 张水明.基于AFLP 和SSR 分子标记的中国杨梅遗传多样性分析[D].杭州：浙江大学，2009.ZHANG Shuiming. Genetic diversity analysis of Chinese bayberry (Myrica rubra Sieb. & Zucc.) based on AFLP and SSR markers[D].Hangzhou：Zhejiang University，2009.
[9] 焦云，柴春燕，舒巧云.基于SSR 分子标记的早熟杨梅优株遗传多样性分析[J].中国南方果树，2019，48(4)：42-45.JIAO Yun，CHAI Chunyan，SHU Qiaoyun. Analysis of genetic diversity of early- maturing elite plants of Chinese bayberry(Myrica rubra)with SSR markers[J].South China Fruits，2019，48(4)：42-45.
[10] 戚行江，任海英，梁森苗，郑锡良，吴阳春.杨梅全基因组测序结果初报[J].浙江农业科学，2015，56(10)：1564-1566.QI Xingjiang，REN Haiying，LIANG Senmiao，ZHENG Xiliang，WU Yangchun. Preliminary report on sequencing results of whole genome of Myrica rubra[J]. Journal of Zhejiang Agricultural Sciences，2015，56(10)：1564-1566.
[11] JIA H M，JIA H J，CAI Q L，WANG Y，ZHAO H B，YANG W F，WANG G Y，LI Y H，ZHAN D L，SHEN Y T，NIU Q F，CHANG L，QIU J，ZHAO L，XIE H B，FU W Y，JIN J，LI X W，JIAO Y，ZHOU C C，TU T，CHAI C Y，GAO J L，FAN L J，VAN DE WEG E，WANG J Y，GAO Z S.The red bayberry genome and genetic basis of sex determination[J]. Plant Biotechnology Journal，2019，17(2)：397-409.
[12] 陈方永，倪海枝，王引，颜帮国.大果杨梅新品种‘永冠’[J].园艺学报，2018，45(6)：1213-1214.CHEN Fangyong，NI Haizhi，WANG Yin，YAN Bangguo. A new big fruit Chinese bayberry cultivar‘Yongguan’[J]. Acta Horticulturae Sinica，2018，45(6)：1213-1214.
[13] 郭秀珠，宋洋，刘冬峰，胡丹，陈巍.杨梅新品种早炭的选育[J].果树学报，2021，38(10)：1821-1823.GUO Xiuzhu，SONG Yang，LIU Dongfeng，HU Dan，CHEN Wei.A new Chinese bayberry cultivar zaotan[J].Journal of Fruit Science，2021，38(10)：1821-1823.
[14] 中华人民共和国农业部.植物新品种特异性、一致性和稳定性测试指南杨梅：NY/T 2761—2015[S].北京：中国农业出版社，2015.Ministry of Agriculture of the People’s Republic of China.Guidelines for the conduct of tests for distinctness，uniformity and stability. Chinese Bayberry (Myrica Linn.)：NY/T 2761—2015[S].Beijing：China Agriculture Press，2015.
[15] 高志红，黄颖宏，倪照君，梁森苗.杨梅种质资源描述规范和数据标准[M].北京：中国农业科学技术出版社，2020.GAO Zhihong，HUANG Yinghong，NI Zhaojun，LIANG Senmiao. Descriptors and data standards for Chinese bayberry(Myrica rubra Sieb. et Zucc.)[M]. Beijing：China Agricultural Science and Technology Press，2020.
[16] ZHANG X N，HUANG H Z，ZHANG Q L，FAN F J，XU C J，SUN C D，LI X，CHEN K S. Phytochemical characterization of Chinese bayberry (Myrica rubra Sieb. et Zucc.) of 17 cultivars and their antioxidant properties[J]. International Journal of Molecular Sciences，2015，16(6)：12467-12481.
[17] 曾少敏，陈小明，黄新忠.福建地方梨资源果实性状多样性分析及其数量分类研究[J].园艺学报，2019，46(2)：237-251.ZENG Shaomin，CHEN Xiaoming，HUANG Xinzhong. Fruit character diversity analysis and numerical classification of local pear germplasm resources in Fujian[J].Acta Horticulturae Sinica，2019，46(2)：237-251.
[18] 吴如健，万继锋，韦晓霞，陈瑾，胡菡青，潘少霖.橄榄种质资源果实表型性状多样性分析及其数量分类研究[J].果树学报，2015，32(5)：797-805.WU Rujian，WAN Jifeng，WEI Xiaoxia，CHEN Jin，HU Hanqing，PAN Shaolin. Fruit character diversity analysis and numerical classification of Chinese olive germplasm resources[J]. Journal of Fruit Science，2015，32(5)：797-805.
[19] 林存学，杨晓华，刘海荣.东北寒地96 份李种质资源表型性状遗传多样性分析[J].园艺学报，2020，47(10)：1917-1929.LIN Cunxue，YANG Xiaohua，LIU Hairong. Genetic diversity analysis of 96 plum germplasm resources by phenotypic traits in northeast cold area[J]. Acta Horticulturae Sinica，2020，47(10)：1917-1929.
[20] 颉刚刚，欧阳丽婷，谢军，刘贝贝，耿文娟.新疆地区欧洲李叶片表型性状多样性及亲缘关系分析[J].植物资源与环境学报，2018，27(3)：72-78.XIE Ganggang，OUYANG Liting，XIE Jun，LIU Beibei，GENG Wenjuan.Analyses on diversity of leaf phenotypic traits and genetic relationships of Prunus domestica in Xinjiang region[J].Journal of Plant Resources and Environment，2018，27(3)：72-78.
[21] 孙珍珠，李秋月，王小柯，赵婉彤，薛杨，冯锦英，刘小丰，刘梦雨，江东.宽皮柑橘种质资源表型多样性分析及综合评价[J].中国农业科学，2017，50(22)：4362-4383.SUN Zhenzhu，LI Qiuyue，WANG Xiaoke，ZHAO Wantong，XUE Yang，FENG Jinying，LIU Xiaofeng，LIU Mengyu，JIANG Dong. Comprehensive evaluation and phenotypic diversity analysis of germplasm resources in mandarin[J]. Scientia Agricultura Sinica，2017，50(22)：4362-4383.
[22] 赫卫，张慧.基于表型性状和SRAP 标记的观赏用辣椒种质资源遗传多样性分析[J].中国瓜菜，2022，35(1)：16-23.HE Wei，ZHANG Hui. Analysis of genetic diversity of pepper germplasm resources for ornamental based on phenotypic traits and SRAP markers[J].China Cucurbits and Vegetables，2022，35(1)：16-23.
[23] 邱立军，张真真，刘洁琪.晚稻杨梅果实外观指标与可溶性固形物含量的关系[J].浙江农业科学，2021，62(4)：768-770.QIU Lijun，ZHANG Zhenzhen，LIU Jieqi. Correlation analysis between fruit appearance and soluble solid content of red bayberry var. Wandao[J]. Journal of Zhejiang Agricultural Sciences，2021，62(4)：768-770.
[24] 戚行江，梁森苗，郑锡良，陈宗良，柴春燕.特早熟杨梅新品种‘早荠蜜梅’[J].园艺学报，2003，30(6)：759.QI Xingjiang，LIANG Senmiao，ZHENG Xiliang，CHEN Zongliang，CHAI Chunyan.A new bayberry cultivar‘zaoqi mimei’[J].Acta Horticulturae Sinica，2003，30(6)：759.
[25] 戚行江，梁森苗，郑锡良，汪国云.优质晚熟杨梅新品种‘晚荠蜜梅’[J].园艺学报，2004，31(1)：136.QI Xingjiang，LIANG Senmiao，ZHENG Xiliang，WANG Guoyun.A new bayberry variety‘Wanqi Mimei’[J].Acta Horticulturae Sinica，2004，31(1)：136.