蝴蝶兰种质资源主要数量性状变异与概率分级
2020-08-04钟海丰陈剑锋陈宇华邱思鑫黄敏玲
钟海丰 陈剑锋 陈宇华 邱思鑫 黄敏玲
摘 要：以53份蝴蝶兰种质资源为研究材料，对24个主要数量性状的变异情况、正态性检验以及概率分级进行研究。结果表明：24个主要数量性状的种间变异幅度为8.22%～45.31%，其中花序花朵数量变异系数最大，合蕊柱长度变异系数最小。种内变异幅度为5.07%～22.50%，变异系数最大的为花序长度，变异系数最小的为花朵宽度。经K-S检验植株大小、叶片长度等23个数量性状Sig值>0.05，符合正态分布，花序花朵数量经c2检验，符合卡方分布，呈现偏态分布，所有考察的数量性状认为均符合正态分布。24个数量性状均分为5级，统一采用（X?1.2818S）、（X?0.5246S）、（X+0.5246S）、（X+1.2818S） 4个点，从而使性状值落入1～5级的概率分别为10%、20%、40%、20%和10%。
关键词：蝴蝶兰;数量性状;概率分级;DUS测试中图分类号：S31 文献标识码：A
Variation and Probability Grading of Main Quantitative Traits ofPhalaenopsisGermplasm Resources
ZHONG Haifeng1， 2， CHEN Jianfeng1， 2， CHEN Yuhua1， 2， QIU Sixin1， 2， HUANG Minling1， 2*
1. Crop Research Institute， Fujian Academy of Agricultural Sciences， Fuzhou， Fujian 350013， China; 2. Fuzhou Sub-center for New Plant Variety Tests， Ministry of Agriculture and Rural Affairs， Fuzhou， Fujian 350013， China
Abstract： The variation， normality test and probability grading analysis of 24 main quantitative traits were carried out， using 53Phalaenopsis germplasmresources as the research materials. The coefficient of interspecific variation of the traits varied from 8.22% to 45.31%. Among them， the coefficient of variation of the flower number was the biggest， while the length of gynostemium was the least. Meanwhile， the coefficient of intraspecific variation of the traits varied from 5.07% to 22.50%. The coefficient of variation of inflorescence length was the biggest， and the flower width was the least. The Sig value of 23 quantitative traits， such as plant size and leaf length， was more than 0.05 by K-S test， which was in accordance with normal distribution， and the number of inflorescence showed skewed distribution according to chi-square test. All the quantitative traits examined were in accordance with normal distribution. The 24 quantitative traits could be divided into five grades. Four points （X-1.2818S）， （X-0.5246S）， （X+0.5246S） and （X+1.2818S） are used uniformly， so that the probability of the trait value falling into 1-5 grades is 10%， 20%， 40%， 20% and 10%， respectively.
Keywords： Phalaenopsis; quantitative traits; probability grading; DUS testing
DOI： 10.3969/j.issn.1000-2561.2020.06.007
品種特异性（Distinctness）、一致性（Unifo?rm?ity）和稳定性（Stability）（简称DUS）测试是定义品种的科学方法，也是目前我国作物品种管理、新品种保护、育种评价的重要技术支撑[1-3]。数量性状是植物形态描述的主要指标体系，数量性状的评价标准是开展种质评价和DUS测试的重要内容[4-5]。同时，数量性状容易受到环境因素的影响，稳定性较差[6]。因此，有必要对数量性状的变异情况、分布类型以及分级标准进行研究分析，建立科学的评价与分级体系。
蝴蝶兰（Phalaenopsis）是兰科（Orchidaceae）蝴蝶兰属植物，其花形奇特，犹如翩翩起舞的蝴蝶，色彩丰富艳丽且花期长，素有“洋兰皇后”的美誉，观赏价值极高。蝴蝶兰于2010年被列入我国第8批农业植物品种保护名录，并于2012年12月7日公告发布了《植物新品种特异性、一致性和稳定性测试指南 蝴蝶兰》[7-8]，但并未明确其数量性状的分级标准。目前，DUS测试中所采用的数量性状分级标准一般基于传统等差分级以及最小显著差（LSD）的等距分级方法[9]，国内对于蝴蝶兰数量性状分级研究主要是基于最小显著差（LSD）的等距分级方法[10-11]，而基于概率分级的蝴蝶兰数量性状分级方法鲜有报道。
数量性状的概率分级是基于数量性状取值的概率分布特征的分级方法，刘孟军[12]于1992年提出，现已广泛应用于大花蕙兰[13]、草莓[14]、多花黑麦草[4]、梯棱羊肚菌[15]、酸枣[16]和桑树[17]等作物的数量性状分级，其能较好地反映性状变异的中值、离散程度以及不同品种性状值在总体变异中的系统位置[18]，可有效服务于种质资源的科学描述、评价、创新利用等。
本研究通过对53份蝴蝶兰种质资源的24个主要数量性状的变异情况与分布规律进行分析，并应用概率分布理论，探索建立蝴蝶兰数量性状等级划分标准。目前，关于蝴蝶兰数量性状的分级研究均是基于等距分级，本文主要利用概率分级方法，对蝴蝶兰数量性状进行非等距分级，使其分布频率符合概率分级的要求，为蝴蝶兰种质资源的综合评价、创新利用以及蝴蝶兰种质资源描述规范和数据标准化提供参考。
1 材料与方法
1.1 材料
试验材料为收集保存于农业农村部植物新品种测试福州分中心的53份蝴蝶兰种质资源。
1.2方法
本研究主要调查《植物新品种特异性、一致性和稳定性测试指南 蝴蝶兰》测试指南中列出的19个数量性状以及测试指南之外新增的5个数量性状（表1）[7-8]。所有数量性状数量均以单株为单位测量，每个品种观测20株，每个植株取样数量为1个。花器官数量性状的观测在花序上有50%的花已开之时，选择最近开花的最大花朵并在其花色没有褪去之前观测，花的长和宽以及花的各部分的观测在花保持自然的状态下进行。叶片的观测应在开花植株的最长叶上进行[8]。
1.3 数据处理
1.3.1 数量性状稳定性与变异情况分析 采用Excel 2017对各供试蝴蝶兰品种数量性状的最大值、最小值、极差、平均值、中值、标准差和变异系数等进行统计分析。
计算单个品种每个性状观测数据的平均值、标准差以及变异系数，再求该性状在所有蝴蝶兰种质资源中变异系数的平均值，即为该性状的品种内变异系数。
根据单个品种每个性状观测数据的平均值，计算出该性状在所有品种中的平均值、标准差以及变异系数，该变异系数即为品种间变异系数。
1.3.2 数量性状正态性检验 采用单样本K-S（Kolmogorov-Smirnov）检验法[19]和卡方（c2）检验法[20]分析所采集的数量性状数据是否符合正态分布或卡方分布。
1.3.3 数量性状概率分级 对符合正态分布的数量性状，参照杨俊霞等[17-18， 21]的概率分级方法，分为5级或3级，这2种分级方法统一采用（X?1.2818S）、（X?0.5246S）、（X+0.5246S）和（X+1.2818S）4个分点或（X?0.5246S）、（X+0.5246S）2个分点。5级中1～5类别出现概率分别为10%、20%、40%、20%、10%，3级中1～3类出现的概率分别为30%、40%、30%，其中X为该性状的平均值，S为该性状的标准差。
对于不符合正态分布的性状，每个数量性状以其0.5～1倍于其标准差的级差进行调整，使其分布于各组的频数符合卡方分布。
2 结果与分析
2.1 蝴蝶兰品种资源数量性状变异情况分析
对53份蝴蝶兰品种资源的数量性状数据进
行统计分析，结果显示，蝴蝶兰各数量性状在品种内的变异范围为5.07%～22.50%，变异系数最大为花序长度（C4），为22.50%;变异系数最小的为花朵宽度（C9），为5.07%（表2）。另外，种内变异系数在10%以上的有花序长度（C4）、花序花朵数量（C5）、花序梗长度（C6），反映这3个性状种内稳定性较差，容易受环境因素影响。在种间变异情况方面，各数量性状的变异系数为8.22%～45.31%，品种间变异丰富，其中花序花朵数量（C5）的变异幅度最大，变异系数为45.31%，合蕊柱长度（C18）的变异幅度最小，变异系数为8.22%。在24个数量性状中，种间变异系数在15%以下的性状3个，分别为唇瓣中裂片长度（C14）、叶片厚度（C17）、合蕊柱长度（C18）。
2.2蝴蝶兰品种资源数量性状的正态性检验
对蝴蝶兰24个数量性状进行K-S正态性檢验，检验结果表明，除花序花朵数量（C5）指标Sig值小于0.05，差异不显著，不符合正态分布外，其他23个数量性状均符合正态分布。对不符合正态分布的数量性状再进行卡方检验，其c2-P值为1.000，符合卡方分布（表3）。因此，本研究所观测的24个数量性状均符合正态分布。
2.3 蝴蝶兰品种资源数量性状分级
对符合正态分布的各数量性状均采用（X? 1.2818S）、（X?0.5246S）、（X+0.5246S）和（X+ 1.2818S）4个分点，分为5级，分别以1、3、5、7、9表示，代表极小、小、中、大、极大或极少、少、中、多、极多等。根据分点计算结果，得到蝴蝶兰各数量性状分级的分值点（表4），并确定每个数量性状分级的取值范围，并计算出其分布频率（表5）。
3 讨论
性状变异是物种进化及新品种和新物种形成的前提，历来被植物分类和育种学所重视[22-23]，变异系数的差异反映了性状在进化过程中遗传可塑性和保守性方面的不同[14]。在种间变异方面，各数量性状的变异系数在8.22%～45.31%之间，表明蝴蝶兰种质资源在植株大小、叶片大小、花部相关数量性状均存在丰富的遗传多样性，这与陈和明等[10]、张鹏等[11]的研究结果相一致。丰富的遗传多样性可为种质资源创新利用提供更大的选择潜力。从种内性状稳定性来看，本研究的花序长度、花序花朵数量、花序梗长度3个数量性状种内变异系数较大，表明3个性状的种内稳定性较差，容易受环境影响，在开展性状数据采集时需确保植株长势、花芽分化条件以及栽培管理条件尽可能一致。
通常认为，在自然状态下，生物现象的连续性变量或间断性变量是符合正态分布的[24]，本研究观测的数量性状在K-S正态性检验中，仅花序花朵数量Sig值小于0.05，差异不显著，不符合正态分布，通过卡方检验，花序花朵数量c2-P值为1.000，符合卡方分布，呈现偏态分布。在进行数量性状的测定时，这一分布形态也可作为正态分布性来看。这和陈和明等[10，25]研究结果相一致。偏态分布在大部分植物上都普遍存在[26-27]。这可能与长期的人工选择压力下，群体向某一特定方向进化，从而使性状分布发生偏移的结果。在蝴蝶兰育种和栽培过程中，对多花性状的长期选择，可能是导致花序花朵数量呈现偏态分布的一个原因。
数量性状的合理分级对种质资源的利用具有重要指导价值[17]。传统的等距离分级或经验分级存在一定的局限性，往往无法准确反映性状变异的中值和离散程度以及各级取值在总体变异中的系统位置，也缺乏统一的分级标准。本研究采用概率分级方法对蝴蝶兰24个数量性状进行分级分析，24个数量性状均可划分为5个等级。与陈和明等[10]采用等距分级方法的蝴蝶兰数量性状分级研究相比，本研究在植株大小、叶片长度、叶片宽度、花朵长度、花朵宽度等性状的等级1取值均较大。通过比较分布频率大小，本研究在植株大小、叶片宽度2个性状的等级1取值范围更为合理。另外，在等级9的取值上，二者研究结果整体相对一致。在花朵数量性状分级方面，与陈和明等[10]的研究结果基本一致，等级1的分布频率均极低，因此，该性状的分级标准有必要在此基础上进行进一步调整。本研究通过分析蝴蝶兰24个主要数量性状的分布类型和分级标准，初步构建了蝴蝶兰种质资源数量性状概率分级指标体系，下一步将继续探讨概率分级方法应用于蝴蝶兰DUS测试中数量性状分级的可行性与科学性。
参考文献
[1] Zhuang X Y， Qin L， Leng P S，et al. Comparison and analysis on DUS screening traits of two laveder species[J]. Agricultural Science and Technology， 2015， 16（5）： 912-917， 925.
[2] 王永行， 白立华， 单飞彪， 等. 基于DUS测试性状的普通小麦测试品种的遗传多样性分析[J]. 黑龙江农业科学， 2017（1）： 4-8.
[3] 钟海丰， 黄敏玲， 钟淮钦， 等. 中国农业植物新品种保护与DUS测试技术发展现状[J]. 热带作物学报， 2017， 38（6）： 1155-1162.
[4] 刘 欢， 马 啸， 张新全， 等. 多花黑麦草品种DUS测试数量性状分级分析[J]. 植物遗传资源学报， 2016， 17（5）： 846-853， 860.
[5] 褚云霞， 邓 姗， 黄志城， 等. 非洲菊新品种DUS测试数量性状分级及形态性状多样性研究[J]. 植物遗传资源学报， 2015， 16（4）： 920-926.
[6] 陈海荣， 杨 华， 王加红， 等. 基于DUS测试性狀的玉米标准品种形态多样性分析[J]. 玉米科学， 2015， 23（2）： 46-51.
[7] 陈和明， 朱根发， 吕复兵， 等. 蝴蝶兰新品种DUS测试指南的研制[J]. 中国农学通报， 2014， 30（10）： 182-185.
[8] 中华人民共和国农业部. 植物新品种特异性、一致性和稳定性测试指南 蝴蝶兰： NY/T 2230—2012 [S]. 北京： 中国农业出版社， 2013.
[9] 王凤华， 郝彩环， 周海涛， 等. 玉米DUS测试主要数量性状分级方法的研究[J]. 玉米科学， 2011， 19（2）： 144-147.
[10] 陈和明， 吕复兵， 李 佐， 等. 蝴蝶兰品质性状综合评价体系的构建[J]. 中国农业大学学报， 2017， 22（8）： 83-94.
[11] 张 鹏， 王江民， 管俊娇， 等. 蝴蝶兰品种数量性状与分组性状的DUS判定[J]. 西北农林科技大学学报（自然科学版）， 2018， 46（11）： 81-88.
[12] 刘孟军. 桃树部分经济性状的种内变异及其分级标准研究[J]. 北京农学院学报， 1992（2）： 98-104.
[13] 张韶伊， 宁惠娟， 范义荣. 大花蕙兰主要数量性状的变异及概率分级[J]. 浙江农业科学， 2013（9）： 1118-1121.
[14] 杨 雷， 杨 莉， 李 莉， 等. 草莓种质资源果实主要数量性状变异及概率分级[J]. 西南农业学报， 2007（5）： 1067-1069.
[15] 陈 影， 唐 杰， 彭卫红， 等. 梯棱羊肚菌种质资源数量性状变异及概率分级[J]. 食用菌学报， 2017， 24（1）： 27-32.
[16] 杨 雷， 周俊义， 刘 平， 等. 酸枣种质资源果实主要数量性状变异及概率分级[J]. 河北农业大学学报， 2006（1）： 34-37.
[17] 杨俊霞， 郭宝林. 桑树数量性状的概率分级研究[J]. 蚕业科学， 1998（2）： 122-124， 130.
[18] 刘孟军. 枣树数量性状的概率分级研究[J]. 园艺学报， 1996（2）： 105-109.
[19] 马小河， 赵旗峰， 董志刚， 等. 鲜食葡萄品种资源果实数量性状变异及概率分级[J]. 植物遗传资源学报， 2013， 14（6）： 1185-1189.
[20] 张新跃， 李元华， 苟文龙， 等. 多花黑麦草研究进展[J]. 草业科学， 2009， 26（1）： 55-60.
[21] 刘 平， 刘孟军， 周俊义， 等. 枣树数量性状的分布类型及其概率分级指标体系[J]. 林业科学， 2003（6）： 77-82.
[22] 王力荣. 我国果树种质资源科技基础性工作30年回顾与发展建议[J]. 植物遗传资源学报， 2012， 13（3）： 343-349.
[23] 王述民， 李立会， 黎 裕， 等. 中国粮食和农业植物遗传资源状况报告（II）[J]. 植物遺传资源学报， 2011， 12（2）： 167-177.
[24] 杜庆鑫， 刘攀峰， 庆 军， 等. 杜仲种质资源雄花主要数量性状变异及概率分级[J]. 北京林业大学学报， 2016， 38（11）： 42-49.
[25] 陈和明， 吕复兵， 朱根发， 等. 蝴蝶兰DUS测试主要数量性状分布特点研究[C]//中国观赏园艺研究进展， 北京： 中国林业出版社， 2013： 94-97.
[26] 胡宏霞， 穆国俊， 侯名语， 等. 河北省花生地方品种基于EST-SSR的遗传多样性及性状标记相关分析[J]. 植物遗传资源学报， 2013， 14（6）： 1118-1123， 1129.
[27] 林 茂， 李正强， 郑治洪， 等. 贵州省花生地方品种的遗传多样性[J]. 作物学报， 2012， 38（8）： 1387-1396.
