北京地区郁金香常见品种的数量分类

分 子植物育种， 2017 年，第 15 卷，第 7 期，第 2863- 2873 页Molecular Plant Breeding, 2017, Vol.15, No.7, 2863- 2873研究报告Research Report北京地区郁金香常见品种的数量分类高星1吕彤2臧凤岐1吕英民1*1 北京林业大学园林学院 , 北京 , 100083; 2 北京市植物园花卉研究所 , 北京 , 100093* 通讯作者 , Luyingminbjfu.edu.cn摘 要 对于北京园林中常见的 69 个郁金香品种进行调查 。根据 中华人民共和国郁金香属 DUS测试指南 ，选定 28 个观赏性状并对其进行编码 。运用 EXCEL2010 软件计算平均值 、标准差 、变异系数 、极差 、百分位数 、多样性指数等；运用 SPSS23.0 软件对选定的 28 个形态性状进行 R 聚类分析 。根据分析选定 24 个形态性状进行 Q 聚类分析，运用 SPSS23.0 软件完成 。结果显示郁金香 69 个品种数量分类结果与形态特征基本一致 。中国的郁金香育种工作还处于起步的阶段，正在大力的推进 。现有品种的整理和准确的品种分类是推动品种选育工作的前提条件和基础 。利用数量分类的方式对品种进行聚类分析，为进一步研究的开展提供相关依据和参考 。关键词 郁金香 , 品种 , 数量分类 , 聚类分析Classification of the Tulipa gesneriana L. Common Varieties in Beijing byNumerical TaxonomyGao Xing1Lv Tong2Zang Fengqi1Lv Yingmin1*1 College of Landscape Architecture, Beijing Forestry University, Beijing, 100083; 2 Institute of Flowers, Beijing Botanical Garden, Beijing, 100093* Corresponding author, luyingminbjfu.edu.cnDOI: 10.13271/j.mpb.015.002863Abstract The paper made the investigation of the tulip (Tulipa gesneriana L.) common varieties in Beijing. Onthe basis of Guidelines for the Conduct of Tests for Distinctness, Uniformity and Stability tulipa (Tulipa L.),selected and encoded 28 morphological traits of 69 varieties. Using EXCEL 2010 to calculate average, standarddeviation, coefficient of variation, range, Shannon-Wiener diversity index etc. And using SPSS23.0 softwaremade R cluster analysis with selected 28 morphological traits. Using SPSS23.0 software made Q cluster analysiswith analyzed and selected 24 morphological traits. The results showed that the result of the numerical taxonomybasically consistent with morphological traits in 69 varieties tulip. Tulip breeding still in the start stage in Chinaand making a big push to improve it. Sorting out varieties and varieties' classification accurately are preconditionsand foundations to promote tulip breeding. We classified tulip varieties by numerical taxonomy and provided thebasis for the further research.Keywords Tulip, Varieties, Numerical taxonomy, Cluster analysis引用格式： Gao X., Lv T., Zang F.Q., and Lv Y.M., 2017, Classification of the Tulipa gesneriana L. common varieties in Beijing bynumerical taxonomy,FenziZhiwuYuzhong(MolecularPlantBreeding),15(7):2863-2873(高星 , 吕彤 , 臧凤岐 , 吕英民 , 2017, 北京地区郁金香常见品种的数量分类 , 分子植物育种 , 15(7): 2863-2873)郁金香 (Tulipa gesneriana L.)又名洋荷花 、草麝香，是百合科 (Liliaceae)郁金香属 (Tulipa L.)的球根花卉 。世界范围内，郁金香属的植物包括约 150 个种 (费砚良等 , 2008, 中国农业出版社 , pp.459-477)。郁金香栽培品种十分丰富，目前世界范围内大约有 8000 多个栽培品种，并且 6 000 多个郁金香栽培品种已经进行了注册，呈逐年增加的趋势 。郁金香花大色艳，栽培品种众多，在世界各地的园林绿地中得到了广泛的推广和应用，深受全世界人民的喜爱 。世界范围内，荷兰 Keukenhof 公园的郁分 子植物育种Molecular Plant Breeding表 1 数量性状基本分析Table 1 Quantitative traits fundamental analysis数量性状Quantitative traits茎长Stem length茎直径Stem diameter花瓣长度Petal length极小值Minimum16.30.33.0极大值Maximum56.60.98.3极差Range40.30.65.3平均值Average32.2960.5505.393标准差Standard deviation9.1360.1291.188变异系数Coefficient of variation0.2830.2350.220金香球根植物展最负盛名 。中国的不少地区也有郁金香花展，如西安植物园 、武汉植物园 、杭州太子湾公园 、上海鲜花港等 (宋碧琰等 , 2015)。1977 年，荷兰女王曾赠送给中国 39 箱郁金香种球，次年北京首届郁金香花展向市民开放参观 (赖娜娜等 , 2010)。早期荷兰皇家园艺协会 、英国皇家园艺协会等权威机构都曾为郁金香进行品种归类 (张金政 , 1999)，先前所用的郁金香品种分类方法是荷兰皇家种球协会按照花期为第一分类等级分出的十四大类 。然而郁金香育种的飞速发展，使得 1981 年出版的分类方法已不能适用于现今的众多品种，很多品种的特性也与原始分类出现不符的现象 。形态学指标上的分类方法是最基础的一种分类方式，数量分类学的迅速发展使得分类学提高了精确性和数量性 。数量分类学诞生于上世纪 50 年代，这个方法使得植物的形态分类不再只是调查和数据罗列方面，而是融入数学思想和计算机技术，使得最终分析结果更具科学性 、准确性 、目标性 (张广进等 , 2006; 骆菁菁等 , 2012)。查阅文献，数量分类学研究以大量应用于海棠 (刘志强和汤庚国 , 2006)、梅花 (毛汉书等 , 1992)、芍药 (刘春迎和王莲英 , 1995)、桂花 (唐东芹 , 1998)、桃花 (张春英等 , 1999)等花卉分类研究中，为形态分类学提供了准确的数学计算依据 。张震 (2003)对中国郁金香属的分布情况进行野外调查和分析 。除此以外没有看到相关郁金香种类区分及品种分类的相关文献，有关郁金香数量分类相关论文更是没有人涉及 。现将调查的 69 个郁金香品种选用 24 个形态性状进行分类，并进行 Q 聚类分析，以期为郁金香品种分类体系建立 、梳理品种亲缘关系 、郁金香新品种选育等提供一定依据 。1 结果与分析1.1 数量性状在品种间的差异性分析郁金香的数量性状相比于其他物种而言数量较少，主要有 3 个：茎长 、茎直径 、花瓣长度 。3 个数量性状的极小值 、极大值 、极差 、平均值 、标准差 、变异系数等已列于表中 (表 1)。通过表中的数值可知， 3 个数量性状的变异系数都在 0.20.3 之间，说明品种间在茎干长度及粗度 、花瓣长度方面的差异比较显著 。茎长的变异范围为 16.356.6，标准差 9.136，品种间变异比较大 。茎干直径因数值较小，相对变化范围也不算大，同时其因存在环境及营养等方面因素的影响，可靠性较另外两种较低 。1.2 质量性状遗传多样性分析对所调查的郁金香质量性状，也就是二元及多态性状进行初步整理分析，发现不同品种间的差异性较大，形态和遗传多样性较高 。其花枝 、叶片 、茎干等描述性状在所有调查品种所占比例 、各性状遗传多样性指数 (H') (表 2)。由此表可以看出，郁金香由于属于观花类草花，花朵性状信息所占比例较大， 25 个多元性状中有 18 个有关花的性状 。叶子形状 、外花被片次色分布 、花朵颜色等性状形态多样性较高 。H值也均在 1 以上 。郁金香主要以单瓣花为主，复瓣品种所占比例较小，其中花色多样性最为丰富，复色品种所占比例大，复色分布情况及颜色多样，其中复色分布在花瓣边缘的情况最多 。叶子出现斑纹和花朵边缘带须边的情况比例比较小，在具体特征的描述中缺失项也就比较多 。一些性状的不同特征分布比例比较平均，如茎干花青苷显色 、叶缘波状等 。郁金香的花期主要以中花期品种为主，这也是园林中郁金香花田布置的必然选择，保证了大面积郁金香的最佳效果，也用少量早花期 、晚花期品种延长了观赏时期 。1.3 R 聚类结果分析根据 R 分析结果可知，对于郁金香品种分类性状的选择基本上是正确的，各形状之间遗传距离较大，也就是性状间可以做到彼此独立影响品种的变异和2864表 2 质量性状多样性Table 2 Variations of qualitative traits质量性状Qualitative trait花枝数量Stem: number of flowers茎花青苷显色Position of anthocyanin coloration叶子形状Leaf shape叶子斑纹Leaf variegation斑纹分布Leaf distribution of variegation斑纹颜色Leaf color of variegation叶边缘波状Leaf undulation of margin花类型Flower type花形状Flower shape花外部颜色数量Flower: number of color outside外花被片次色分布Flower distribution of secondarycolor on the outsider花须边Flower fringe花须边程度Conspicuousness of fringe花须边位置Flower position of fringe on tepals外花被片先端形状Shape of tip of outer tepal是否有色斑Flower color of macule on pental色斑边缘有否颜色Different color of border of maculeon pental雄蕊颜色Stamen color of filament花粉颜色Stamen color of pollen花期Plant beginning of floweringH' 值H' value0.345 60.997 91.477 90.302 50.356 70.379 20.692 30.325 10.770 60.707 81.565 20.366 90.418 60.450 10.331 80.302 50.673 10.950 90.818 10.978 51 个 89One 89不显色 52Absent 52椭圆形窄 1Narrow elliptic 1卵圆形窄 33Narrow ovate 33无 91Absent 91边缘 1On margin 1白色 3White 3无 52Absent 52单瓣 90Single 90重瓣 12Double 121 种 24One 24无 28Absent 28火焰状 9Flamed 9无 88Absent 88弱 7Weak 7仅在顶部 6Top only 6渐尖 33Acuminate 33无 9Absent 9无 59Absent 59黄色 40Yellow 40黄色 70Yellow 70早 24Early 24>1 个 11More than one 11顶端显色 32Distal part only32椭圆形 23Medium elliptic 23卵圆形 28Medium ovate 28有 9Present 9边缘环状 4Margin zone 4黄绿色 1Yellow green 1有 48Present 48复瓣 10Double 10椭圆形 77Ellipsoid 772 种 72Two 72边缘 10Margin 10晕染状 13Flushed 13有 12Present 12中 4Intermediate 4整个边缘 6All margin 6急尖 15Acute 15有 91Present 91有 41Present 41紫色 38Purple 38黄 +紫色 16Yellow+purple16中 57Medium 57整个显色 16Whole stem 16椭圆形宽 11Broad elloptic 11卵圆形宽 4Broad ovate 4条状 3Stripes 3红色 1Red 1卵圆形 4Ovoid 43 种及以上 4Three and more 4边缘环状 36Margin zone 36底色状 4Ar base 4圆形 19Rounded 19蓝色 9Blue 9紫 , 黑色 14Purple, black 14晚 19Late 19紫色 3Purple 3百合花形 7Lily flower 7微凹 33Emarginated 33多样性特征 (%)Variation (%)北京地区郁金香常见品种的数量分类Classification of the Tulipa gesneriana L. Common Varieties in Beijing by Numerical Taxonomy2865分 子植物育种Molecular Plant Breeding图 1 28 个表型性状 R 聚类分析Figure 1 R cluster analysis of 28 traits图 2 69 个品种 Q 聚类分析注 : 169 见表 3Figure 2 Q cluster analysis of 69 varietiesNote: 169 in table 3进化，不会相互间造成不必要的干扰 。但是在部分性状间存在关联性大的现象 (图 1)。其中，花须边 、花须边程度 、花须边位置这 3 个性状以及叶子斑纹 、斑纹分布 、斑纹颜色这 3 个性状之间的关联性很大，并且根据前文所述，两类性状存在很多缺失项，因此可以仅保留花须边有无及叶子斑纹有无这 2 个性状，其余 4 个子性状在 Q 聚类分析中进行去除 。1.4 Q 聚类结果分析利用 69 个品种 24 个性状进行品种之间的聚类分析 。根据聚类结果并结合所做直线将品种分类 (图 2)。结合线 L1 将 69 个品种按照花瓣长度分为两大类，2866第 一类包括从 44 号品种到 25 号品种在内的 48 个品种，它们的共同特征是花瓣长度较长，基本上都大于平均值以上；第二类包括从 21 号品种到 39 号品种在内的 21 个品种，它们的共同特征是花瓣长度较短，基本上都大于平均值以下 。结合线 L2 又可以将分出的两大类别进行细化分类 。首先将第一大类分为两小类 。第一类包括从44号品种到 37 号品种在内的 23 个品种，它们的共同特征是花粉的颜色基本上全部带有紫色成分；第二类包括从 13 号品种到 25 号品种在内的 25 个品种，它们的共同特征是花粉的颜色全部为黄色，没有紫色成分 。第二大类根据线 L2 可以分出四小类 。第一类主要是 干杯 、皇室维尔京 、南极洲 、狐步舞 这些白色花色的品种；第二类主要是 阿巴 、黑色英雄 、维多娜 、黄绣球 等重瓣品种；第三类主要是玲曲 、查托 、粉色博瑞 、积累 等红紫色系品种；品种 绿色田野 被单独划为一类，花色以绿色为主 。这样第二大类就按照花色及花类型进行聚类 。不难发现，这些主要分类依据都属于花的性状 。郁金香作为观花球根植物，主要的品种性状差异存在于花朵上是必然的，这也说明在区分郁金香的形态 、培育郁金香的新品种时等，可以将重点放在花朵的性状上进行更多考虑 。2 讨论2.1 数量分类学应用于观赏植物品种鉴定与分类的有效性人们在进行植物的分类和鉴定时，最基础也最直接的方式就是进行形态上的调查和研究 。数量分类学也越来越多的被应用于品种资源分类的分析上 。骆菁菁等 (2012)对北京地区的海棠品种进行调查，以 25 个品种为材料，选取了 20 个性状以数量分类的方法进行了聚类分析 。聚类结果将海棠品种按照花色与花的类型分为粉色重瓣 、红色重瓣 、深色单瓣 、浅色单瓣这 4 种品种群，同时得出结论，数量分类聚类的结果与形态分类结果基本一致 。洪艳等 (2012)以 811 个菊花品种为试验材料，利用色差仪测色的方法对这些品种的花色表型值进行测定，同时进行数量分类研究 。结果得出，只通过花色表型得到的数值来进行聚类分析的结果无法将菊花花色的特征完全体现出来，而通过 ISCC-NBS 方法来定义花色优于前者 。笔者认为其原因有可能是提供的分析性状较少无法全面涵盖花色信息 。我们对北京常见的郁金香品种进行形态上的调查 。并通过数量性状分析 、质量形状遗传多样性分析 、R 聚类分析三种主要方式首先从 28 个性状中筛选出 24 个主要性状进行 Q 聚类分析 。通过对 69 种郁金香品种 24 个观赏性状进行调查并进行聚类分析发现，可以按照花色 、花瓣长度 、花粉颜色 、花类型等重要性状分成五大类，基本上符合形态上的分类，对郁金香的分类也是有效的，同时也对其他观赏植物品种鉴定 、分类等方面具有借鉴意义 。2.2 解决的问题及提出的新问题本研究通过统计分析和数量分类，将郁金香按照 5 个花性状将 69 个品种进行分类 。荷兰种球协会曾在 1981 年提出按照花期为一级分类标准将郁金香分为三大类 14 个小类 。通过实际调查发现，很多品种与实际情况不符 。例如分类中提到的重瓣早花型 (DET)中所包含的 导火线 、重瓣晚花型 (DLT)中所包含的 查托 、须边型 (FT)并没有包含明显带有须边的 查托 、绿色田野 、幸福一代 等品种等 。笔者认为可以得出，近些年郁金香育出了千余个新品种，而依旧使用几十年前的分类方法已无法适应众多的新品种需要的基本结论 。本研究提出的分类方式也将给新的郁金香分类方法制定提供依据 。在实际调查的过程中，笔者发现部分质量性状存在较大的主观性，使得在判断测量时存在一定的误差 。例如花形状 、花次色分布情况等性状很可能存在不同人测量不同标准的情况，使得推广性较低 。再者由于栽植郁金香的小环境不同，导致相同品种间的花期 、株高等性状受到外在环境 、营养条件等方面的影响，存在较大的误差 。然而通过结果可以看出，总体来讲没有对分类造成很大的影响 。由此，完善形态性状鉴定标准和方法 、同时使用多种分类方法进行综合分析以及从品种演化的角度进行分类和鉴定品种间关系是十分有必要发展的方向 。近年来，包括孢粉学分类 、同工酶标记 、分子标记等手段已大量应用于观赏植物品种分类，也将是郁金香品种分类工作下一步的目标 。3 材料与方法3.1 试验材料共调查北京植物园常用品种 69 种 (表 3)。对其进行了完整生长阶段的观察 、形态纪录和拍照 (图 3)。以这 69 个品种作为分类运算单位 (OTU)，每个样本北京地区郁金香常见品种的数量分类Classification of the Tulipa gesneriana L. Common Varieties in Beijing by Numerical Taxonomy2867分 子植物育种Molecular Plant Breeding表 3 调查郁金香品种Table 3 Name of tulip varieties编号No.123456789101112131415161718品种Varieties阿巴Abba阿伯拉Abra王子Ad rem美国梦American dream南极洲Antarctica杏色印记Apricot impression阿玛尼Armani班雅Banja luka玲曲Bell song领袖Big chief黑色英雄Black hero蓝钻石Blue diamond羞美人Blushing beauty粉色博瑞Bolroyal pink花边红葡萄酒Burgundy lace卡马乔Camargue新烛光Candela坎迪王子Candy prince编号No.3637383940414243444546474849505152品种Varieties雄鹅狂想曲Gander rhapsody金色牛津Golden oxford金色检阅Golden parade绿色田野Green land幸福一代Happy generation荷兰美人Holland beauty积累Jackpot硕大粉红Jumbo pink功夫Kung-fu利欧Leo visser玛里琳Marilyn蒙戴尔Madame lefeber小黑人Mistress新构思New design武士Ninja奥利奥斯Ollioules牛津精华Oxford's elite编号No.5354555657585960616263646566676869品种Varieties检阅Parade粉色印记Pink impression普瑞斯玛Purissima紫旗Purple flag红色印象Red impresion皇室维尔京Royal virgin撒哈拉赛车Sahara rally肉色艾奇Salmon van eijk花花女郎Shirley/playgirl纯金Strong gold甜蜜的爱Sugar love西内德阿莫Synaeda amor维多娜Verona白王朝White dynasty世界的喜爱World's favourite黄绣球Yellow pomponnet阿波罗精华Apeldorn's elite编号No.1920212223242526272829303132333435品种Varieties卡罗拉Carola查托Chato干杯Cheers圣诞颂歌Christmas carol净水Clearwater征服者Conqueror破碎的甜蜜Crispion sweet水晶星Crystal star道琼斯Dow Jones端庄小姐Elegant lady逃亡Escape一级品First class网络旗Flaming flag热情鹦鹉Flaming parrot导火线Flash point狐步舞Foxtrot饰边优雅Fringed elegance每项性状指标测取 5 个重复 。3.2 分类性状的选取和编码根据农业部制定的 中华人民共和国郁金香属DUS 测试指南 ，选取 28 个性状进行调查 。其中，二元性状 6 个，多元性状 19 个，数量性状 3 个 。根据品种性状调查的实际情况和 DUS 测试指南，对这些二元性状和多元性状进行编码，分别采用三种编码方法： (1)二元性状，如叶片是否有斑纹等性状，以 1， 2表示； (2)多元性状，如茎干花青苷显色情况等性状，以 1， 2， 3 表示； (3)数量性状，如株高等性状，用实际测量值表示 。所选性状观测方法及所选性状的详细编码方法列于表中 (表 4; 表 5)。将所有数据测定后录入Excel，计算各性状平均值并汇总，称为原始数据矩阵 。2868表 4 性状观测方法及性状类型Table 4 Type of trait and the observation methods编号No.12345678910111213141516171819202122性状Trait茎长Stem length茎直径Stem diam花枝数量Number of stem茎花青苷显色Position of anthocyanin coloration叶子形状Leaf shape叶子斑纹Leaf variegation斑纹分布Leaf distribution of variegation斑纹颜色Leaf color of variegation叶边缘波状Leaf: undulation of margin花类型Flower type花长度Flower length花形状Flower shape外花被片外主色Flower main color花外部颜色数量Flower number of color outside外花被片次色分布Flower distribution of secondarycolor on the outsider外花被片内次色Flower secondary color outside花须边Flower fringe花须边程度Conspicuousness of fringe花须边位置Flower position of fringe on tepals外花被片先端形状Shape of tip of outer tepal外花被片内主色Flower main color inside外花被片内次色Flower secondary color inside观测方法Observation method用直尺测量茎干的长度 (cm), 随机重复 5 次取均值Measure the distant between top to bottom by ruler, repeat 5 times用直尺测量茎干最大直径 (cm), 随机重复 5 次取均值Measure maximum diameter by ruler, repeat 5 times目测Visual observation目测Visual observation目测Visual observation目测Visual observation根据郁金香属新品种测试指南内容进行记录Record data according to UPOV guidelines Tulipa根据郁金香属新品种测试指南内容进行记录Record data according to UPOV guidelines Tulipa目测Visual observation目测Visual observation用直尺测量花瓣长度 (cm), 随机重复 5 次取均值Measure length of petal by ruler, repeat 5 times根据郁金香属新品种测试指南内容进行记录Record data according to UPOV guidelines Tulipa用英国皇家比色卡 (RHS)测出数值Measure by RHS card目测Visual observation根据郁金香属测试指南内容进行记录Record data according to UPOV guidelines Tulipa用英国皇家比色卡 (RHS)测出数值Measure by RHS card目测Visual observation根据郁金香属测试指南内容进行记录Record data according to UPOV guidelines Tulipa根据郁金香属测试指南内容进行记录Record data according to UPOV guidelines Tulipa根据郁金香属测试指南内容进行记录Record data according to UPOV guidelines Tulipa用英国皇家比色卡 (RHS)测出数值Measure by RHS card用英国皇家比色卡 (RHS)测出数值Measure by RHS card性状类型Type of trait数量性状Quantitative数量性状Quantitative二元性状Binary trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait二元性状Binary trait数量性状Quantitative多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait二元性状Binary trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait北京地区郁金香常见品种的数量分类Classification of the Tulipa gesneriana L. Common Varieties in Beijing by Numerical Taxonomy2869分 子植物育种Molecular Plant Breeding图 3 部分郁金香品种注 : 1: 粉色博瑞 ; 2: 坎迪王子 3: 狐步舞 ; 4: 玲曲 ; 5: 道琼斯 ; 6:Menton unique; 7: 网络旗 ; 8:紫旗 ; 9: 阿波罗精华Figure 3 Parts of Tulip varietiesNote: 1: Bolroyal pink; 2: Candela; 3: Foxtort; 4: Bell song; 5:Dow jones; 6: Menton unique; 7: Flaming parrot; 8:Purple flag;9: Apeldorn's elite编号No.232425262728性状Trait花内有否色斑Flower color of macule on petal色斑边缘有否颜色Different color of border of macule on petal色斑边缘颜色Color of border of macule on petal雄蕊颜色Stamen color of filament花粉颜色Stamen color of pollen花期Plant beginning of flowering观测方法Observation method目测Visual observation目测Visual observation用英国皇家比色卡 (RHS)测出数值Measure by RHS card目测Visual observation目测Visual observation目测Visual observation性状类型Type of trait二元性状Binary trait二元性状Binary trait多元性状Multi trait多元性状Multi trait多元性状Multi trait多元性状Multi trait续表 4Continuing table 43.3 数量性状分析方法按照表 1 将所观测的数量性状数据化，对于数量性状，计算其平均值 、标准差 、极值 、极差 、变异系数，其中变异系数计算公式为： CV=/，其中 是标准差， 是均值 (许莹修 , 2005)。3.4 遗传多样性指数分析遗传多样性指数即 Shannon-Wiener 多样性指数(Shannon-Wiener diversity index, H)可以利用各性状的相对出现频率描述质量性状的总体多样性，同时对此进行评价 (张辕 , 2014):H'=- ni=1Piln Pi(i=1, 2, 3,., 10)上式中， Pi为某性状第 n 级别内其品种数量占品种总数的百分比， Ln 是自然对数 。遗传多样性指数的计算由 Microsoft Excel 2010 完成 。3.5 形态性状 R 分析本研究利用 SPSS23.0 软件对经过前两步的分析过后所选出的性状进行 R 分析 。R 分析是以性状为分类基本单位进行的分类运算，结果表现的是不同性状之间地内在关联 。3.6 Q 聚类分析以所研究的 69 个郁金香品种作为分类运算单位，根据之前步骤中分析的结果，选择出 24 个基本性状进行聚类分析 。聚类分析应用 SPSS23.0 软件首先进行 Z 得分的标准化，之后选用欧式距离平方 、最远临近法进行 Q 聚类分析 。作者贡献高星是本研究的设计和实验研究的执行人，并完成数据分析和处理 、论文初稿的写作和修改；吕彤和臧凤岐参与部分实验操作；吕英民是总构思者及负责人，指导实验设计，数据分析，论文写作与修改 。2870表 5 性状及编码Table 5 Traits and codes编号No.123456789101112131415161718192021性状Trait茎长Stem length茎直径Stem diameter花枝数量Number of squids茎花青苷显色Position of anthocyanin coloration叶子形状Leaf shape叶子斑纹Leaf variegation斑纹分布Leaf distribution of variegation斑纹颜色Leaf color of variegation叶边缘波状Leaf undulation of margin花类型Flower type花长度Flower length花形状Flower shape外花被片内主色Flower main color花外部颜色数量Flower number of color outside外花被片次色分布Flower distribution of secondary color on the outsider外花被片内次色Flower secondary color outside花须边Flower fringe花须边程度Conspicuousness of fringe花须边位置Flower position of fringe on tepals外花被片先端形状Shape of tip of outer tepal外花被片内主色Flower main color inside编码Code实际数值Actual numerical value实际数值Actual numerical value(1)1 个 ; (2)大于 1 个(1) One; (2) More than one(1)不显色 ; (2)仅在顶端显色 ; (3)整个茎干显色(1)