DOI：10.19462/j.cnki.zgzy.20240830001

摘要：為了對柳州地區(qū)種植的甘蔗新品系進行評價，為甘蔗和相關(guān)農(nóng)作物選育提供科學的決策支持，采用主成分分析（PCA）、K-均值聚類、相關(guān)性分析等方法對23個甘蔗新品系的出苗率、分蘗率、枯心率、梢腐病等9個性狀進行統(tǒng)計分析。結(jié)果表明，分蘗率的變異系數(shù)最大，達62.49%，有效莖數(shù)、出苗率、產(chǎn)量變異系數(shù)也在10.00%以上；相關(guān)性分析表明，出苗率和有效莖數(shù)高的品系往往具有更高的產(chǎn)量，分蘗率高的品系具有更高的有效莖數(shù)但糖分會受到影響；主成分分析將9個性狀降維為特征值大于1.0的4個主成分，累計方差貢獻率達76.242%；23個新品系通過K-均值聚類被分為5個類群，類群1～5分別包含1、3、3、5、11個品系，類群1、類群2、類群3品系綜合表現(xiàn)較差，不建議在柳州地區(qū)進行推廣種植，類群3的桂糖17-847、桂南蔗17-416、桂糖18-220整體上糖分含量較高，可以作為高糖親本加以利用，類群4、類群5的桂糖18-229、柳糖16-302、河糖16-179、桂糖18-268、桂糖18-277、桂糖13-113等16個品系可以進入下一步的區(qū)域試驗和大田示范。研究證明，多元統(tǒng)計方法的綜合利用可以為柳州選種和引種提供更加全面的品系綜合評價依據(jù)。

關(guān)鍵詞：柳州；甘蔗育種；多元統(tǒng)計分析；主成分分析；K-均值聚類

Evaluation of New Sugarcane Varieties in Liuzhou Based on Multivariate Statistical Analysis

YANG Yebin，SU Xiaohui，DENG Si，WEI Huiming，QIN Yaoguan，BEI Jianshe，

ZHOU Jianping，HE Xiaomei，MING Fanqiang，WEI Kaijun

（Liuzhou Branch of Guangxi Academy of Agricultural Sciences/Liuzhou Agricultural

Sciences Research Center，Liuzhou 545000，Guangxi）

甘蔗在全球范圍內(nèi)被視為重要的糖源作物，對于食品工業(yè)和能源產(chǎn)業(yè)具有不可替代的作用。廣西作為全國最大的糖生產(chǎn)地區(qū)之一，2023-2024年榨季廣西預計糖料蔗入榨量5000萬t，食糖產(chǎn)量600萬t，因此甘蔗的種植和加工對于廣西經(jīng)濟的發(fā)展具有顯著貢獻[1]。全球氣候變化和市場需求的不斷演變對甘蔗產(chǎn)業(yè)提出了諸多要求，包括提高糖分產(chǎn)量、增強作物的抗逆性和適應性，以及滿足可持續(xù)發(fā)展的要求等。甘蔗育種是應對這些挑戰(zhàn)的關(guān)鍵手段，傳統(tǒng)的育種方法依賴于田間選擇和雜交，而現(xiàn)代分子生物學技術(shù)的應用，如基因工程和基因組選擇，為甘蔗的遺傳改良提供了新的途徑[2-4]。無論采取何種育種技術(shù)，核心目標始終是培育出高產(chǎn)、高糖、抗病性強和適應性強的甘蔗新品種[5-6]。由于廣西不同地區(qū)的氣候和土壤條件各異，所以選育出適宜當?shù)胤N植的品種對于促進當?shù)馗收岙a(chǎn)業(yè)的持續(xù)健康發(fā)展具有重要意義[7]。

本研究采用多元統(tǒng)計分析，包括主成分分析（PCA）、K-均值聚類等，深入分析了23個甘蔗新品系的9個性狀，這些分析方法系統(tǒng)揭示了性狀之間的內(nèi)在聯(lián)系，全面評價了品系表現(xiàn)，旨在篩選出具有優(yōu)良遺傳潛力的甘蔗品系。前人通過此類方法對菜心、甘藍、甘蔗等作物育種提供了科學的決策支

持[8-10]。此外，本研究還探討了主成分分析和聚類分析在甘蔗育種中的應用，以期為培育符合未來市場需求的新品種提供新的視角和方法。

1 材料與方法

1.1 試驗材料 本研究所用23個甘蔗新品系來源于廣西壯族自治區(qū)農(nóng)業(yè)科學院甘蔗研究所、廣西南亞熱帶農(nóng)業(yè)科學研究所、柳州市農(nóng)業(yè)科學研究中心、百色市農(nóng)業(yè)科學研究所、河池市農(nóng)業(yè)科學研究所，柳城05-136、桂糖42號為試驗對照（CK1、CK2）

（表1）。

1.2 試驗地概況 試驗地點位于廣西壯族自治區(qū)柳州市柳北區(qū)沙塘鎮(zhèn)柳州市農(nóng)業(yè)科學研究中心試驗基地，地理位置24°26′N，109°22′E，海拔123m，地形以丘陵、沖積平原為主，土壤質(zhì)地紅壤土，肥力中等，前茬作物為甘蔗。

1.3 試驗設計 試驗采用隨機區(qū)組設計，每小區(qū)4行，行長8.0m，行距1.2m，面積38.4m2。采用雙芽段品字形種植，下種量為10萬芽/hm2，四周設保護行。田間管理按照當?shù)馗收嵘a(chǎn)常規(guī)管理進行。于4月中旬調(diào)查出苗率，5月中旬調(diào)查總苗數(shù)，7月調(diào)查枯心率、梢腐病，11月至翌年1月進行糖分測定及甘蔗株高、莖徑、產(chǎn)量的調(diào)查。

1.4 統(tǒng)計分析 采用主成分分析（PCA）、K-均值聚類和相關(guān)性分析等多元統(tǒng)計方法對試驗數(shù)據(jù)進行處理和分析。

2 結(jié)果與分析

2.1 農(nóng)藝性狀變異分析 變異系數(shù)是衡量作物遺傳多樣性的重要指標，變異系數(shù)越大越能反映品種豐富的遺傳多樣性，從中選出優(yōu)良個體的概率越大[11-13]。由表2可知，23份材料的性狀差異明顯，各性狀變異系數(shù)在7.56%～62.49%之間。其中，分蘗率的變異系數(shù)最大，這表明在參與試驗的甘蔗品系中，分蘗能力的遺傳差異性最大；分蘗率變幅在6.8%～161.8%之間，最高的是桂糖18-276，最低的是桂糖17-847。有效莖數(shù)變異系數(shù)為20.14%，變幅在48983～99529株/hm2之間，最多的是桂糖17-932，最少的是桂糖17-847。出苗率變異系數(shù)為19.78%，變幅在47.3%～91.5%之間，最高的是河糖16-179，最低的是河糖16-546。產(chǎn)量變異系數(shù)為15.58%，變幅在38363～72547kg/hm2之間，最高的是桂糖18-229，最低的是桂糖17-847。株高、莖徑、糖分的變異系數(shù)均小于10.00%，說明這3個性狀變異程度較低，遺傳性狀相對穩(wěn)定。

2.2 性狀相關(guān)性分析 相關(guān)性分析揭示了不同性狀間的相互作用，對于指導育種工作至關(guān)重要。由表3可以看出，出苗率與有效莖數(shù)呈極顯著正相關(guān)，與產(chǎn)量呈顯著正相關(guān)，與分蘗率、梢腐病、枯心率、莖徑、株高、糖分呈負相關(guān)，這表明出苗率高的品系往往能夠帶來更高的產(chǎn)量，但糖分、抗病性會受到一定影響。分蘗率與有效莖數(shù)呈極顯著正相關(guān)，但與糖分含量呈顯著負相關(guān)，與梢腐病、莖徑、株高呈負相關(guān)，說明分蘗率越高的品種糖分含量可能越低。除此之外，有效莖數(shù)與產(chǎn)量呈極顯著正相關(guān)，說明甘蔗有效莖數(shù)是產(chǎn)量提高的關(guān)鍵因子。因此，在育種過程中需要對各性狀間的相關(guān)關(guān)系進行平衡，以確保既能提高有效莖數(shù)和產(chǎn)量，又不會降低糖分含量及抗病性。

2.3 主成分分析 主成分分析揭示了不同性狀對甘蔗總體表現(xiàn)的貢獻度。針對9個農(nóng)藝性狀進行分析，以特征值大于1.0為標準提取主成分。結(jié)果表明，前4個主成分的特征值累計方差貢獻率達76.242%，表明這4個主成分可以反映出9個農(nóng)藝性狀的大部分信息（表4）。第1主成分特征值為3.0931，方差貢獻率為30.931%，有效莖數(shù)和產(chǎn)量載荷值較大；第2主成分特征值為1.7075，方差貢獻率為17.075%，梢腐病抗性和莖徑載荷值較大；第3主成分特征值為1.6043，方差貢獻率為16.043%，莖徑、出苗率和產(chǎn)量載荷值較大；第4主成分特征值為1.2193，方差貢獻率為12.193%，株高和梢腐病抗性載荷值較大。

2.4 K-均值聚類分析 對25份參試材料進行K-均值聚類分析，可將其分為5個類群。由表5和表6可知，類群1包含1個品系（百蔗17-344），表現(xiàn)為抗枯心病強、抗梢腐病相對稍弱、莖徑大、產(chǎn)量中等。類群2包含3個品系（桂糖14-817、桂糖18-226、河糖16-546），出苗率中等，抗枯心病和梢腐病強，莖徑大，有效莖數(shù)、糖分含量、產(chǎn)量均中等。類群3包含3個品系（桂糖17-847、桂南蔗17-416、桂糖18-220），抗枯心病和梢腐病強，莖徑小，糖分含量高，產(chǎn)量低。類群4包含2個對照和5個品系（桂糖18-229、柳糖16-302、河糖16-179、桂糖17-605、桂糖18-496），出苗率高，抗枯心病和梢腐病強，莖徑、有效莖數(shù)、糖分含量均中等，產(chǎn)量高。類群5包含11個品系（桂糖18-277、桂糖18-268、桂糖13-113、桂糖17-932、桂糖17-1187、桂南蔗14-2234、桂糖18-313、桂糖17-862、桂糖18-276、桂糖16-1129、桂糖17-877），出苗率中等，分蘗率高，抗枯心病和梢腐病強，莖徑中等，有效莖數(shù)高，糖分含量一般，產(chǎn)量較高。

表5 參試材料K-均值聚類

類群 品種（系）

1 百蔗17-344

2 桂糖14-817、桂糖18-226、河糖16-546

3 桂糖17-847、桂南蔗17-416、桂糖18-220

4 柳城05-136（CK1）、桂糖42號（CK2）、桂糖18-229、柳糖16-302、河糖16-179、桂糖17-605、桂糖18-496

5 桂糖18-277、桂糖18-268、桂糖13-113、桂糖17-932、桂糖17-1187、桂南蔗14-2234、桂糖18-313、桂糖17-862、

桂糖18-276、桂糖16-1129、桂糖17-877

3 討論與結(jié)論

本研究采用多元統(tǒng)計分析方法對23個甘蔗新品系進行適應性綜合評價分析。在9個性狀中株高、莖徑變異系數(shù)小于10.00%，說明遺傳性狀相對穩(wěn)定。糖分為甘蔗品種選育中的一個重要指標，通過對比新品系中變異系數(shù)大于10.00%的分蘗率、有效莖數(shù)、出苗率、產(chǎn)量4個性狀發(fā)現(xiàn)，分蘗率低于平均值的新品系糖分絕大部分高于平均值，出苗率及分蘗率高的品系有效莖數(shù)均高于平均值，出苗率高的品系往往分蘗率較低。研究表明，分蘗率提高不僅能增加有效莖數(shù)，還能直接影響甘蔗的產(chǎn)量和糖分含量，這與前人[14-17]研究結(jié)果一致，因此甘蔗分蘗的增加對甘蔗生產(chǎn)力的提高有重要作用，如何平衡出苗率和分蘗率是品種選育過程中需要著重考慮的問題。

K-均值聚類將23個甘蔗品系分為5個類群，類群4（桂糖18-229、柳糖16-302、河糖16-179等）和類群5（桂糖18-277、桂糖18-268、桂糖13-113等）表現(xiàn)出高產(chǎn)和高糖分的潛力，可以作為下一步區(qū)域試驗和大田示范的候選品系。類群3的桂糖17-847、桂南蔗17-416、桂糖18-220產(chǎn)量較低，但糖分含量高，可以將其作為高糖親本進行改良。

本文結(jié)合了主成分分析、聚類分析、相關(guān)性分析等多元統(tǒng)計方法，研究再次證明利用變異系數(shù)分析可評估育成品系的遺傳多樣性，以此篩選出遺傳背景豐富、性狀表現(xiàn)多樣的品系作為候選育種材料；利用相關(guān)性分析可篩選出與目標性狀呈正相關(guān)或負相關(guān)的性狀，作為未來品種改良和維系性狀平衡的依據(jù)；主成分分析能夠確定對目標性狀貢獻最大的主成分，從各主成分代表性狀中篩選出表現(xiàn)優(yōu)異的品系；K-均值聚類分析可將種質(zhì)資源分為不同類群，識別具有特定優(yōu)勢和改良潛力的品系，便于其進一步試驗和示范。通過綜合評價體系的建立，可使新品系得到更全面的評價，可為柳州甘蔗品種選育和引種提供科學依據(jù)。本文僅對甘蔗9項生理指標進行了多元評價，更合理的選育和引種評價體系的構(gòu)建有待進一步深入研究。

參考文獻

[1] 李楊瑞．關(guān)于廣西的甘蔗育種．廣西糖業(yè)，2019（3）：3-7

[2] 劉新龍，李旭娟，劉洪博，馬麗，徐超華，范源洪．云南甘蔗常用親本資源遺傳多樣性的SSR分析．植物遺傳資源學報，2015，16（6）：1214-1222

[3] 王英，陳守俊，朱相成，莊南生．80份甘蔗種質(zhì)RAMP標記遺傳多樣性分析．植物遺傳資源學報，2011，12（4）：525-532

[4] Atkin F C，Dieters M J，Stringeer J K．系譜深度和歷史資料對甘蔗育種方案加性方差和育種值估計的影響．理論與應用遺傳學，2009，119（3）：555-565

[5] 楊榮仲，周會，唐仕云，劉昔輝，梁強，黃海榮．廣西甘蔗育種思考．中國糖料，2021，43（3）：18-27

[6] 胡朝輝，潘永保，Werapon P．蔗糖主產(chǎn)國甘蔗單產(chǎn)及產(chǎn)糖率改良進展．中國糖料，2021，43（4）：75-80

[7] 蘇振，鄧宇馳，鐘德發(fā)，王春玲，林飛劍，陳家福．甘蔗新品種的篩選試驗及評價．中國種業(yè)，2023（5）：58-61

[8] 郭強，江清梅，何洪良，黎正英，梁永檢，秦昌鮮，唐利球．基于主成分分析和聚類分析評價141份中蔗優(yōu)良品系資源．熱帶作物學報，2024，45（1）：49-59

[9] 趙勇，趙俊，昝逢剛，趙麗萍，趙培方，朱建榮，Burner D M，吳才文，劉家勇．86份甘蔗種質(zhì)資源工藝性狀的評價．湖南農(nóng)業(yè)大學學報：自然科學版，2019，45（5）：466-471

[10] 史衛(wèi)東，梁勁，張力，琚茜茜，羅海玲，康紅衛(wèi)．菜心和芥藍表型多樣性的主成分分析與聚類分析．西南農(nóng)業(yè)學報，2020，33（12）：2726-2735

[11] 王恩軍，陳垣，韓多紅，蔡子平，張芬琴，張勇．菘藍農(nóng)藝性狀與藥材產(chǎn)量的相關(guān)和通徑分析．核農(nóng)學報，2018，32（2）：399-406

[12] 王慧敏，彭振英，李新國，萬書波，張智猛，丁紅，高文偉．67個花生品種主要農(nóng)藝性狀的變異及相關(guān)性分析．山東農(nóng)業(yè)科學，2019，51（9）：91-96

[13] 梁森苗，張淑文，鄭錫良，任海英，朱婷婷，戚行江．楊梅生長指標與果實品質(zhì)間的相關(guān)性分析．核農(nóng)學報，2019，33（4）：751-758

[14] 范業(yè)賡，丘立杭，陳榮發(fā)，周慧文，李燕嬌，閆海鋒，羅含敏，周忠鳳，鄧宇馳，吳建明．不同分蘗特性甘蔗品種生產(chǎn)力差異和相關(guān)性研究．中國農(nóng)業(yè)科技導報，2023，25（12）：35-43

[15] 陸鑫，毛鈞，應雄美，劉新龍，蘇火生，馬麗，蔡青．甘蔗創(chuàng)新種質(zhì)的因子分析與聚類分析．西南農(nóng)業(yè)學報，2011，24（6）：2072-2076

[16] 張革民，楊榮仲，劉海斌，方位寬．割手密主要數(shù)量性狀的主成分及聚類分析．西南農(nóng)業(yè)學報，2006，19（6）：1127-1131

[17] 俞華先，經(jīng)艷芬，安汝東，郎榮斌，田春艷，董立華，桃聯(lián)安，孫有芳，楊李和，邊芯，周清明．基于主成分與聚類分析的大莖野生種血緣后代育種潛力評價．江西農(nóng)業(yè)學報，2019，31（10）：16-22

（收稿日期：2024-08-30）

基金項目：廣西科技重大專項（桂科AA22117002）；科技先鋒隊“強農(nóng)富民”“六個一”專項行動（桂農(nóng)科盟202403-9）；國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系廣西甘蔗創(chuàng)新團隊建設項目（nycytxgxcxtd-2021-03）

通信作者：韋開軍