文章信息
- 赵勇, 赵培方, 赵俊, 赵丽萍, 昝逢刚, 姚丽, 刘家勇
- ZHAO Yong, ZHAO Peifang, ZHAO Jun, ZHAO Liping, ZAN Fenggang, YAO Li, LIU Jiayong
- 基于投影寻踪分类法评价43份澳大利亚甘蔗种质资源
- Evaluation of 43 Australian sugarcane germplasm based on projection pursuit clustering
- 亚热带农业研究, 2019, 15(1): 7-13
- Subtropical Agriculture Research, 2019, 15(1): 7-13.
- DOI: 10.13321/j.cnki.subtrop.agric.res.2019.01.002
-
文章历史
- 收稿日期: 2018-12-02
种质资源评价是甘蔗杂交育种的重要工作, 主要对工农艺性状进行比较,为杂交利用提供参考。我国甘蔗杂交亲本体系的建立主要基于引进种质资源与国内自主创制相结合[1]。国外甘蔗种质资源在我国甘蔗亲本体系和杂交育种中占有重要地位[2],外引亲本对我国甘蔗品种改良发挥了重要作用[3-4]。据不完全统计,自20世纪50年代至2009年,我国利用引进种质育成的甘蔗品种有163个[3]。目前对甘蔗种质资源的研究主要集中在两方面。一是表型评价[5-8],即通过测量、比较表型性状而实现评价和分类,是最直观的评价;二是基于分子标记技术对遗传基础深入解析[9-13]。其中,研究甘蔗表型性状主要采用综合分析法[14-16], 即通过比较单项指标以及指标之间的因果关系或相关性而对样本加以分析。但不同因子对综合评价的影响程度(即权重)主要凭研究者的经验进行评判, 往往受到研究者主观性的影响[17-18]。
1974年Friedman和Turkey提出了投影寻踪分类法(projectionpursuit clustering, PPC)[19]。所谓投影,就是从不同角度观察数据,寻找能够最大程度地反映数据特征和最能充分挖掘数据信息的最优投影方向[20]。其特点是在未知权重系数的情况下,通过把高维数据投影到低维(1~3维)子空间上[21]。该方法省去了专家打分评定的步骤,避免了专家打分的人为干扰,因此更为准确和便捷,在定量评价指标数据的处理上更具优势。目前该方法已应用于水质评价[22]、环境质量综合评价[23]、工业经济[24]等方面。在农作物研究上的应用主要集中在作物补偿效应[25]、水分利用评价[26]等,尚未见甘蔗上的应用报道。基于此, 本研究利用投影寻踪分类法对引进的澳大利亚甘蔗种质资源进行评价,挖掘优异的种质资源,以期为我国甘蔗种质创新提供参考,同时也为甘蔗种质资源评价探索新方法。
1 材料与方法 1.1 试验材料参试的澳大利亚杂交种质资源共计43份,其中,Q型40份、KQ型3份。
1.2 试验设计 1.2.1 试验地与田间种植试验在云南省农业科学院甘蔗研究所育种基地进行。该试验地为水田,海拔1 051.8 m,属亚热带高原季风气候,光照资源充足,年日照时数2 382 h,年均气温20 ℃,年均降雨量771.1 mm,年潜在蒸发量1 987 mm,无霜期341 d左右。2016年12月29日进行田间试验(新植)。采用随机区组设计,3次重复,行长4.0 m、行距1.1 m,每个种质、每个重复种植2行,每小区下芽52个,四周设保护行。试验地管理同大田生产相一致。
1.2.2 农艺性状调查农艺性状调查采用分级方法,主要基于本试验的分级需要并结合育种经验而设定指标参数范围。2017年10月25日至27日,田间调查43份种质的株高、茎径、有效茎、叶部病害严重度和总体长势共5个农艺性状。其中,叶部病害严重度指田间自然发病的花叶病、锈病、褐条病等叶部病害的综合反映。之后由3~4名本研究团队成员组成的分级小组以观察到的性状表现为分级依据,对小区种质群体进行整体评估、讨论分级。每个性状指标均分为5级,分级的数值范围见表 1。
分级 | 株高/cm | 茎径/cm | 有效茎 | 叶部病害 严重度/% |
总体长势 |
株·cm-1 | |||||
1 | ≥280 | ≥3.0 | ≥10 | ≤5 | 根据观察,对种质材料进行综合评价并分级。 |
2 | 220~280 | 2.5~3.0 | 8~9 | 5~10 | |
3 | 160~220 | 2.0~2.5 | 6~7 | 10~20 | |
4 | 100~160 | 1.5~2.0 | 4~5 | 20~30 | |
5 | ≤100 | ≤1.5 | ≤3 | ≥30 |
于2017年11月至2018年3月对43份种质进行锤度、蔗糖分、蔗汁糖分、简纯度和纤维分检测。各个种质每次随机检测6条蔗茎(每重复3条蔗茎),共检测5次。其中,糖分、蔗汁糖分、简纯度和纤维分检测参考Liu et al[27]的方法。
1.3 数据分析利用Excel 2013整理数据。采用DPS统计分析软件对各农艺性状分级数据和主要品质性状指标进行方差分析和基本参数估计;采用投影寻踪分类法对43份种质资源进行投影方向和投影值大小分析,并基于投影值的聚类分析对各种质进行分类评价,以筛选优异种质。其中,农艺性状指标原始分级数据先进行倒数处理,然后再进行相关统计分析。
2 结果与分析 2.1 主要工农艺性状差异分析 2.1.1 农艺性状方差分析表明, 各种质农艺性状差异均极显著(P < 0.01),变异系数在24.30%~35.21%之间,变异广泛(表 2)。
农艺性状 | 均值 | 均方 | F值 | 变异系数 |
% | ||||
总体长势 | 2.91** | 1.41 | 4.36 | 28.89 |
茎径 | 2.98** | 1.05 | 3.18 | 24.30 |
有效茎 | 2.98** | 1.09 | 2.56 | 24.85 |
叶部病害 | 2.94** | 1.41 | 2.46 | 28.54 |
株高 | 2.42** | 1.45 | 2.54 | 35.21 |
1)**表示差异达极显著水平(P < 0.01)。 |
方差分析表明, 各种质品质性状间均存在极显著差异(P < 0.01),变异系数在5.31%~11.99%之间(表 3)。除纤维分变异较大外,其他变异较稳定。
品质性状 | 均值 | 均方 | F值 | 变异系数 |
% | ||||
锤度 | 20.89** | 2.46 | 7.74 | 5.31 |
蔗糖分 | 14.69** | 1.42 | 5.56 | 5.74 |
蔗汁糖分 | 17.96** | 2.68 | 6.68 | 6.45 |
简纯度 | 85.75** | 6.09 | 2.93 | 2.04 |
纤维分 | 13.04** | 4.89 | 6.38 | 11.99 |
1)**表示差异达极显著水平(P < 0.01)。 |
根据投影寻踪分类法原理,投影方向各分量的大小反映了该指标对甘蔗种质综合评价的影响程度。利用DPS统计分析软件对43份参试种质的主要农艺性状分级数据(倒数)和品质性状(11月—翌年3月均值)进行投影寻踪分类法分析,得到各性状的投影方向(图 1)。从图 1可见,叶部病害严重度对甘蔗综合指标的影响最大,其次为甘蔗总体长势;茎径、有效茎、株高和蔗糖分的影响相当,纤维分和简纯度等其他品质性状指标影响较小。
2.2.2 投影值计算确定投影方向后,基于主要工农艺性状的投影寻踪分析,得到各甘蔗种质资源的投影值(图 2)。从图 2可以看出,43份种质投影值在0~2之间。其中,Q179最高(1.83),Q72最低(0.32),且各种质投影值大小差异显著。
2.2.3 基于聚类分析对43份种质进行分类基于参试种质投影值大小,采用系统聚类,以欧氏距离-WPGMA法对43份种质进行分类。在0.16~0.32的距离上可把43份种质分为5类(图 3)。其中,Ⅰ~Ⅴ类群分别有5、8、14、14、2个。进一步对各类群投影值进行分析表明,第Ⅰ~Ⅴ类群投影均值分别为1.41、1.73、1.02、0.72、0.37,且类群间存在显著差异(表 4)。第Ⅰ~Ⅳ类群投影值变异系数在2.78%~17.51%之间,第Ⅳ类群变异最大(17.51%)。第Ⅴ类群个数太少,故未统计变异系数。
类群 | 个数 | 投影值 | 变异系数/% | ||
均值 | 最大值 | 最小值 | |||
Ⅰ | 5 | 1.41 | 1.48 | 1.37 | 2.78 |
Ⅱ | 8 | 1.73 | 1.83 | 1.59 | 4.81 |
Ⅲ | 14 | 1.02 | 1.27 | 0.92 | 12.18 |
Ⅳ | 14 | 0.72 | 0.83 | 0.52 | 17.51 |
Ⅴ | 2 | 0.37 | 0.42 | 0.32 |
为了鉴别聚类结果是否可靠,采用多类逐步判别法对其进行判别(表 5)。结果表明,第Ⅰ、Ⅱ、Ⅴ类群判别准确率均为100.00%,第Ⅲ、Ⅳ类群判别准确率分别为92.86%、78.57%。根据判别分析原理可知,判别准确率越高,表明种质被正确分类的机率越大。从表 5可以看出,每个种质被正确分类的概率都较高,平均判别准确率为90.70%。
类群 | 聚类个数 | 预测个数 | 判别准确率 | ||||||
个 | Ⅰ类群 | Ⅱ类群 | Ⅲ类群 | Ⅳ类群 | Ⅴ类群 | % | |||
Ⅰ | 5 | 5 | 0 | 0 | 0 | 0 | 100.00 | ||
Ⅱ | 8 | 0 | 8 | 0 | 0 | 0 | 100.00 | ||
Ⅲ | 14 | 1 | 0 | 13 | 0 | 0 | 92.86 | ||
Ⅳ | 14 | 0 | 0 | 2 | 11 | 1 | 78.57 | ||
Ⅴ | 2 | 0 | 0 | 0 | 0 | 2 | 100.00 |
从参试种质综合投影值的聚类结果来看,第Ⅰ、Ⅱ类群综合投影值较高,且判别准确率为100%。根据投影寻踪原理,投影值越大,其综合性状表现越好,故筛选出13份优异种质(Ⅰ类群5份,Ⅱ类群8份)。从表 6可以看出,13份种质蔗糖分分级在13.89%~16.02%之间;总体长势分级在1.33~2.33之间;茎径分级在2.33~3.00之间;有效茎分级在1.33~3.00之间;叶片病害分级在1.67~3.33之间;株高分级在1.00~2.67之间(3重复平均值)。综合来看, 主要工农艺性状总体评价较好。
种质资源 | 类群 | 投影值 | 蔗糖分 | 总体长势 | 茎径 | 有效茎 | 叶部病害 | 株高 |
Q183 | Ⅱ | 1.83 | 16.00 | 2.00 | 3.00 | 2.00 | 2.33 | 1.33 |
Q179 | Ⅱ | 1.79 | 15.15 | 1.33 | 2.00 | 2.33 | 2.33 | 1.00 |
Q127 | Ⅱ | 1.79 | 15.81 | 2.00 | 2.33 | 2.00 | 2.00 | 2.00 |
Q197 | Ⅱ | 1.69 | 13.89 | 2.00 | 2.00 | 1.67 | 2.00 | 1.00 |
KQ01-1390 | Ⅱ | 1.68 | 15.27 | 2.00 | 2.00 | 3.00 | 1.67 | 2.67 |
Q124 | Ⅱ | 1.68 | 15.14 | 1.67 | 2.00 | 2.33 | 2.00 | 1.67 |
Q172 | Ⅱ | 1.67 | 14.28 | 2.00 | 2.33 | 1.33 | 2.00 | 1.00 |
Q189 | Ⅱ | 1.59 | 14.60 | 2.00 | 2.33 | 2.00 | 2.00 | 1.00 |
Q188 | Ⅰ | 1.48 | 15.81 | 2.00 | 2.67 | 2.67 | 3.00 | 1.33 |
Q152 | Ⅰ | 1.40 | 14.78 | 2.00 | 3.00 | 2.00 | 2.67 | 1.67 |
KQ01-1261 | Ⅰ | 1.40 | 15.77 | 2.33 | 2.67 | 3.00 | 2.00 | 2.33 |
Q206 | Ⅰ | 1.40 | 15.51 | 2.33 | 2.67 | 2.67 | 2.00 | 2.67 |
Q121 | Ⅰ | 1.37 | 16.02 | 2.00 | 3.00 | 2.00 | 3.33 | 2.67 |
1)以上分级数据为3次重复平均值。 |
近年来,我国不断引进国外甘蔗种质应用于育种[28]。因此,客观评价引进种质对筛选优异种质具有现实意义。蒲光兰等[29]研究了四川核桃种质资源坚果的数量性状变异及概率分级,发现数量性状分级数据1~5级出现的频率符合正态分布;蔡志翔等[30]完善了桃种质资源托叶长度评价方法,同时提出了桃托叶长度分级指标;Zhao et al[31]将甘蔗总体长势分为9级应用于新品种选育,发现在育种第一阶段总体长势遗传变异大于锤度,应当利用总体长势和锤度对亲本进行选择,以优化杂交组合。基于上述研究结果,本研究将重要农艺性状分级法应用于甘蔗种质评价,发现43份澳大利亚甘蔗种质农艺性状分级数据差异达到极显著水平(P < 0.01),变异广泛。
对43份澳大利亚甘蔗种质进行投影寻踪分类法分析表明,甘蔗叶部病害、总体长势的投影方向值最大,有效茎、茎径、蔗糖分、株高等次之;简纯度和纤维分投影方向较小。根据投影寻踪分类法原理,各指标对种质综合评价的影响程度依次为总体长势、有效茎、茎径、蔗糖分、株高;简纯度和纤维分影响最低。育种中,叶部病害发生程度决定了甘蔗材料成为推广品种的潜力,而总体长势代表了甘蔗品种对环境的适应性水平。同时,选种时通常给予甘蔗糖分和产量较大的选择压力,其他品质性状指标选择压力相对较小。投影寻踪分类法较好地体现了各性状指标的权重,客观评价了各种质资源实际情况。43份种质主要农艺性状变异系数在24.30%~35.21%之间,主要品质性状变异系数在5.51%~11.99%之间,表明主要农艺性状变异较广、糖分品质变异较小。进一步说明投影寻踪分类法确定主要农艺性状有较高投影方向的合理性。根据投影值大小,筛选优良种质13份(Ⅰ、Ⅱ类群),占全部种质数的30.23%;中等或以下的种质合计30份,占69.77%。说明在数量众多的甘蔗品种资源中,表现优良的种质数量并不多,因此对数量众多的种质资源进行客观和有效的评价具有重要的意义。
通过对主要农艺性状进行分级和检测蔗糖分等品质性状,本研究建立了以甘蔗叶部病害、总体长势、有效茎、茎径、株高和甘蔗糖分为主要权重的投影寻踪分类法以评价甘蔗种质的方法,筛选了13份主要工农艺性状优良(投影值较高)外引种质用于杂交育种。投影寻踪分类法只要合理选择评价指标, 其评价结果符合实际, 往往能达到预期效果,同时有利于提高种质精准评价和杂交利用的针对性。
[1] | 彭绍光. 甘蔗育种学[M]. 北京: 农业出版社, 1990: 290-342. |
[2] | 吴才文. 甘蔗亲本创新与突破性品种培育的探讨[J]. 西南农业学报, 2005, 18(6): 858–861. DOI: 10.3969/j.issn.1001-4829.2005.06.044 |
[3] | 张琼, 齐永文, 张垂明, 等. 我国大陆甘蔗骨干亲本亲缘关系分析[J]. 广东农业科学, 2009(10): 44–48. DOI: 10.3969/j.issn.1004-874X.2009.10.011 |
[4] | 齐永文, 邓海华, 李奇伟. 我国大陆甘蔗种质资源利用进展[J]. 作物研究, 2012, 26(5): 443–446. DOI: 10.3969/j.issn.1001-5280.2012.05.08 |
[5] | 黄忠兴, 周峰, 王勤南, 等. 国内外割手密资源农艺性状表型遗传多样性分析[J]. 植物遗传资源学报, 2012, 13(5): 825–829. DOI: 10.3969/j.issn.1672-1810.2012.05.020 |
[6] | 吴才文, 范源洪, 刘家勇, 等. 甘蔗引进品种的区试表现及利用前景[J]. 亚热带农业研究, 2005, 1(2): 9–14. |
[7] | 赵俊, 吴才文, 赵培方, 等. 引进甘蔗种质工艺与农艺性状的相关性及聚类分析[J]. 湖南农业大学学报(自然科学版), 2012, 38(5): 476–481. |
[8] | 刘家勇, 陈学宽. 19个云南自育甘蔗品种(系)基础种质分析[J]. 亚热带农业研究, 2004, 11(3): 30–33. |
[9] | SINGH R B, SINGH B, SINGH R K. Study of genetic diversity of sugarcane (Saccharum) species and commercial varieties through TRAP molecular markers[J]. Indian Journal of Plant Physiology, 2017, 22(7): 1–7. |
[10] | HEMAPRABHA G, GOVINDARAJ P, BALASUNDARAM N, et al. Genetic diversity analysis of indian sugarcane breeding pool based on sugarcane specific STMS markers[J]. Sugar Tech, 2005, 7(2/3): 9–14. |
[11] | 刘新龙, 李旭娟, 刘洪博, 等. 云南甘蔗常用亲本资源遗传多样性的SSR分析[J]. 植物遗传资源学报, 2015, 16(6): 1214–1222. |
[12] | SELVI A, MUKUNTHAN N, SHANTHI R M, et al. Assessment of genetic relationships and marker identification in sugarcane cultivars with different levels of top borer resistance[J]. Sugar Tech, 2008, 10(1): 53–59. DOI: 10.1007/s12355-008-0009-x |
[13] | 宋弦弦.甘蔗及其近缘植物遗传多样性的SRAP和TRAP标记分析[D].福州: 福建农林大学, 2009. http://cdmd.cnki.com.cn/Article/CDMD-10389-2009170718.htm |
[14] | 周会, 杨荣仲, 方锋学, 等. 桂糖系列甘蔗种质资源数量性状的主成分分析和聚类分析[J]. 西南农业学报, 2012, 25(2): 390–395. DOI: 10.3969/j.issn.1001-4829.2012.02.008 |
[15] | 王水琦, 杨坤正, 郭陈福, 等. 三种分析法在甘蔗品种综合评价中的应用[J]. 甘蔗, 1996, 3(4): 22–24. |
[16] | 刘新龙, 马丽, 苏火生, 等. 甘蔗杂交品种核心种质重要农艺性状评价及亲缘关系分析[J]. 植物遗传资源学报, 2014, 15(1): 67–73. |
[17] | 刘家勇, 陈学宽, 范源洪, 等. 甘蔗品种模糊综合评判[J]. 亚热带农业研究, 2005, 1(1): 4–7. |
[18] | 经艳芬, 边芯, 郎荣斌, 等. 云瑞创新亲本选配32个甘蔗杂交组合的模糊综合评判[J]. 亚热带农业研究, 2015, 11(4): 217–221. |
[19] | FRIEDMAN J H, TUKEY J W. A projection pursuit algorithm for exploratory data analysis[J]. IEEE Transactions on Computers, 1974, 23(9): 881–890. |
[20] | JONES M C, SIBSON R. What is projection pursuit?[J]. Journal of the Royal Statistical Society Series A, 1987, 150(1): 1–37. DOI: 10.2307/2981662 |
[21] | BOLTON R J, KRZANOWSKI W J. Projection pursuit clustering for exploratory data analysis[J]. Journal of Computational and Graphical Statistics, 2003, 12(1): 121–142. DOI: 10.1198/1061860031374 |
[22] | WANG S J, YANG Z F, DING J. Projection pursuit cluster model and its application in water quality assessment[J]. Journal of Environmental Sciences, 2004, 16(6): 994–995. |
[23] | 王顺久, 杨志峰. 区域农业生态环境质量综合评价投影寻踪模型研究[J]. 中国生态农业学报, 2006, 14(1): 173–175. |
[24] | 姚小芹, 董催军. 中国重化工业循环经济发展研究-基于投影寻踪的评价分析[J]. 科学学研究, 2010, 28(12): 1855–1860. |
[25] | 郝树荣, 郭相平, 张展羽. 投影寻踪分类模型在作物补偿效应评价中的应用[J]. 农业机械学报, 2010, 41(1): 59–62, 33. DOI: 10.3969/j.issn.1000-1298.2010.01.011 |
[26] | 吴新新. 基于改进投影寻踪模型的地表覆盖对作物水分利用的影响评价[J]. 水资源开发与管理, 2016(2): 75–78. |
[27] | LIU J Y, BASNAYAKE J, JACKSONP A, et al. Growth and yield of sugarcane genotypes are strongly correlatedacross irrigated and rainfed environments[J]. Field Crops Research, 2016, 196: 418–425. DOI: 10.1016/j.fcr.2016.07.022 |
[28] | 赵培方, 夏红明, 刘家勇, 等. 6个Q型甘蔗品种杂交利用的遗传力和配合力[J]. 湖南农业大学学报(自然科学版), 2013, 39(4): 348–353. |
[29] | 蒲光兰, 肖千文, 蔡利娟, 等. 四川核桃种质资源坚果的数量性状变异及概率分级[J]. 湖南农业大学学报(自然科学版), 2015, 41(6): 647–650. |
[30] | 蔡志翔, 严娟, 沈志军, 等. 桃种质资源托叶长度评价与分级体系的建立[J]. 植物遗传资源学报, 2016, 17(3): 461–465. |
[31] | ZHAO D L, COMSTOCK J C, GLAZ B, et al. Vigor rating and brix for first clonal selection stage of the canal point sugarcane cultivar development program[J]. Journal of Crop Improvement, 2012, 26(1): 60–75. DOI: 10.1080/15427528.2011.611926 |