海陸回交群體主要農藝性狀與纖維品質關系的探討

李生梅, 楊濤, 黃雅婕, 任丹, 耿世偉,李典鵬, 芮存, 高文偉*

(1.新疆農業大學農學院, 教育部棉花工程研究中心, 烏魯木齊 830052;2.新疆農業大學草業與環境科學學院, 烏魯木齊 830052)

陸地棉和海島棉是兩個最重要的四倍體栽培種經濟作物，前者因產量高而主導著世界棉花的生產；后者的纖維品質較好，但單產較低，由此其種植面積受到了限制[1-2]。快速有效地培育高產優質棉花新品種是當前育種的主攻目標。雜交有利于將優異基因進行轉移和集中，快速創造具有優異農藝性狀和品質的新材料[3]。由于纖維品質與農藝性狀間存在復雜的相互作用，因此，通過簡單的相關分析、偏回歸分析、通徑分析和主成分分析等[4-8]統計方法，探究棉花纖維品質與主要農藝性狀的關系，發掘與棉花纖維品質密切相關的農藝性狀，有助于育種工作者通過對農藝性狀的選擇和改良來提高纖維品質。

只有對各性狀充分了解，才能降低表型性狀鑒定的盲目性，提高優異種質的篩選效率。國內外學者對陸地棉群體的纖維品質開展了多方面研究[9-12]，對海陸種間雜交群體纖維品質的遺傳特性也進行了系統研究[13-20]，發現海島棉在纖維品質和產量性狀方面表現出明顯優勢，以篩選的優異材料作為改良親本在產量和品質上也有較好的配合力。李慧琴等[21]對270份陸地棉種質資源進行遺傳多樣性分析，篩選出9份材料可用于棉花產量、馬克隆值和株型的改良。劉翔宇等[22]對126份新疆陸地棉資源進行綜合性評價，篩選出了7個類別的棉花品種，同時對應類別提出陸地棉的改良方案。唐淑榮等[23]以2005—2014年國家531份參試品種的纖維品質數據為研究對象，對雜交棉和常規棉品種的纖維品質進行了綜合評價。但是，以上研究多針對單一的農藝性狀或纖維品質性狀，有關海陸回交群體主要農藝性狀與纖維品質間關系的研究鮮見報道。

因此，本研究選取各性狀趨于穩定的海陸高代回交群體BC4F2:4為材料，擬通過相關性分析、主成分分析、逐步回歸分析和通徑分析等方法，探究海陸回交群體纖維品質與農藝性狀的關系，為新疆及周邊地區選育高產優質的棉花新品種提供理論參考。

1 材料與方法

1.1 試驗材料

以陸地棉系9、海島棉新海16號以及它們的回交群體(BC4F2∶4)為試驗材料，均由新疆農業大學農學院遺傳育種系提供。其中，輪回親本陸地棉系9是新疆自育材料，豐產性好、抗逆性強，但纖維品質相對較差；供體親本海島棉新海16號具有優良的纖維品質。

1.2 試驗設計

田間試驗于2018年在沙灣144團新疆農業大學棉花基地(E43°20′—45°20′，N84°45′—86°40′)進行，采用隨機區組設計。2018年4月21日播種，采用一膜六行的種植模式，播幅2.35 m。每份材料種植兩行，行長2 m，行距0.7 m，株距0.1 m。田間管理同大田管理模式。

1.3 測定項目與方法

于2018年9月25日，每份材料選取連續的5株調查株高(plant height，PH)、始節高(height of the first branch，BSH)、始節數(number of the first branch，NOB)、果枝數(fruit branch number，FBN)、有效果枝數(effective fruit branch，EFNB)、鈴數(boll number per plant，BN)和有效鈴數(effective boll number per，EBN)。于10月3日收獲單株進行室內考種，測定單株籽棉重和皮棉重，計算衣分(lint percentage，LP)，同時稱取10 g纖維樣品，送往石河子農墾科學院棉花品質監督檢驗中心(HV11000)測定纖維品質，包括纖維長度(fiber length，FL)、纖維強度(fiber strength，FS)、馬克隆值(micronaire，FM)和纖維整齊度(fiber uniformity，FU)。具體測定方法參照《農作物田間試驗記載項目及標準》[24]。

1.4 數據處理與分析

采用Excel 2010對試驗數據進行整理分析，采用SPSS 21.0進行描述性統計、相關分析、主成分分析、多元逐步回歸分析和通徑分析。

2 結果與分析

2.1 親本及群體各性狀的描述性分析

親本及群體各性狀的描述性統計分析結果(表1)表明，輪回親本陸地棉系9較供體親本新海16有效鈴數較多、衣分較高、株型較矮，但纖維品質一般；而供體親本新海16擁有優良的纖維品質，但有效鈴數較少、衣分較低，株型較高。海陸回交群體株高、纖維長度和纖維強度的平均值均高于輪回親本系9；衣分和馬克隆值的平均值高于供體親本新海16；其余性狀的平均值均低于雙親。群體中不同性狀的變異系數為0.11%～48.14%，多數性狀的變異幅度較大，表明群體的遺傳多樣性較高。各性狀的變異系數從大到小依次為：鈴數(48.14%)>有效鈴數(46.61%)>有效果枝數(31.59%)>果枝數(30.43%)>始節高(25.38%)>始節數(19.37%)>馬克隆值(16.95%)>株高(16.41%)>纖維強度(9.76%)>纖維長度(6.22%)>纖維整齊度(1.68%)>衣分(0.11%)。其中，鈴數的變異系數最大，衣分的變異系數最小。

表1 親本及群體各性狀的描述性統計Table 1 Descriptive statistics analysis of all traits in parents and populations

2.2 海陸回交群體主要農藝性狀與纖維品質性狀的相關性分析

海陸回交群體主要農藝性狀與纖維品質性狀的相關分析(表2)表明，纖維長度與株高、有效果枝數、鈴數和有效鈴數呈顯著正相關；纖維強度與始節高呈顯著正相關；馬克隆值與衣分呈極顯著正相關，與株高、有效果枝數、鈴數和有效鈴數呈顯著負相關；纖維整齊度與株高和有效果枝數呈顯著正相關。且纖維品質和農藝性狀的各自性狀之間也存在不同程度的相關性，如纖維長度與纖維強度和纖維整齊度呈極顯著正相關，與馬克隆值呈極顯著負相關；株高與始節高、始節數、有效果枝數和有效鈴數呈極顯著正相關。由此可見，纖維品質與農藝性狀之間存在一定相關性。

表2 海陸回交群體纖維品質與農藝性狀的相關分析Table 2 Correlation analysis of fiber quality and agronomic traits in backcross populations of Gossypium hirsutum L.×Gossypium barbadense L.

2.3 海陸回交群體纖維品質與農藝性狀的主成分分析

為了更加充分的了解纖維品質與農藝性狀的關系，采用主成分分析法獲取綜合變量研究農藝性狀對棉花纖維發育的影響(表3)。提取的4個主成分特征值分別為2.803、2.157、1.330和1.235，貢獻率分別為28.031%、21.567%、13.299%和12.355%，累計貢獻率達75.252%。其中，PC1主要由有效鈴數(0.819)、鈴數(0.816)和纖維長度(0.627)決定，將其命名為“棉花成鈴數相關因子”，但馬克隆值的特征向量為負值，說明群體有效鈴數多、鈴數多、纖維長度長的材料馬克隆值較低；PC2特征向量值最高的為始節高(0.603)，其次是纖維強度(0.590)、長度(0.564)和整齊度(0.558)，將其命名為“纖維品質選擇因子”；PC3載荷較高的是株高和始節高，與纖維品質的特征向量均為負值，故命名為“低品質控制因子”；PC4載荷較高的是衣分和馬克隆值，可作為衣分和馬克隆值選擇的綜合因子。

表3 海陸回交群體纖維品質與農藝性狀的主成分載荷矩陣Table 3 Principal component loading matrix of fiber quality and agronomic traits in backcross populations of Gossypium hirsutum L.×Gossypium barbadense L.

2.4 海陸回交群體纖維品質與農藝性狀的多元逐步回歸分析

為進一步確定纖維品質性狀的決定因子，以4個纖維品質性狀為因變量，8個農藝性狀為自變量進行多元逐步回歸分析，構建最佳回歸方程，結果表明(表4)，每個回歸模型都具有顯著的統計學意義(P<0.05)。其中，鈴數能夠單獨解釋纖維長度變異的4.8%，鈴數和果枝數共同解釋其變異的8.5%，表明鈴數和果枝數與棉花纖維長度關系密切；始節高和鈴數與纖維強度有較強正相關關系；衣分、有效鈴數和始節高與馬克隆值有較強相關性；而株高只與纖維整齊度有較強正相關關系。由此可見，始節數和有效果枝數對纖維品質影響較小，可在進一步的分析中將其剔除。

表4 海陸回交群體纖維品質與農藝性狀的逐步回歸分析Table 4 Stepwise regression analysis of fiber quality and agronomic traits in Gossypium hirsutum L.× Gossypium barbadense L. backcross populations

2.5 海陸回交群體纖維品質與農藝性狀的通徑分析

為進一步明確海陸回交群體中農藝性狀對纖維品質的影響，以纖維品質性狀為因變量，以農藝性狀為自變量進行通徑分析，結果(表5)表明，果枝數(-0.234)、始節高(0.223)和鈴數(0.182)對棉花纖維長度有直接影響；鈴數(0.647)、有效鈴數(-0.421)和始節高(0.152)對棉花纖維強度有直接影響；馬克隆值為因變量時，有效鈴數(-0.362)和衣分(0.286)對其有直接影響，有效鈴數多，馬克隆值較小；纖維整齊度受鈴數(0.204)和株高(0.091)的直接影響較大。

表5 海陸回交群體纖維品質與農藝性狀的通徑分析Table 5 Path analysis of fiber quality and agronomic traits in backcross populations of Gossypium hirsutum L.×Gossypium barbadense L.

除直接影響外，株高通過始節高對棉花纖維長度的間接影響較大(0.102)；果枝數(0.076)和有效鈴數(0.178)通過鈴數對棉花纖維長度產生間接影響；有效鈴數(0.634)、果枝數(0.272)和株高(0.164)通過鈴數對棉花纖維強度產生間接影響；鈴數(-0.354)、果枝數(-0.165)、始節高(0.082)和衣分(0.037)通過有效鈴數對馬克隆值產生間接影響；株高(0.052)通過鈴數對棉花纖維整齊度產生間接影響。綜上所述，株高、鈴數和有效鈴數對纖維品質的間接作用占主導地位。

3 討論

棉花纖維品質受基因型與環境共同影響，屬于數量性狀，因此，對其性狀的操控能力是決定棉花育種效率的關鍵[25]。本研究以新疆陸地棉系9和海島棉新海16號及其回交后代群體為材料，對農藝性狀和纖維品質進行分析，結果表明，12個性狀的變異系數除纖維整齊度和衣分外均大于5%，其中，鈴數的變異系數最大，與尹會會等[26]和李海明等[27]研究結果相一致；而衣分的變異系數較小，可能是由于親本間衣分差異較大，導致回交群體衣分的變異系數很小。總體來看，海陸回交群體的表型變異較為豐富，說明將海島棉引入陸地棉的遺傳改良對于擴大陸地棉遺傳多樣性具有重要意義。

植株形態結構的改變會影響作物的產量和品質[28]，因此棉花農藝性狀與纖維品質性狀之間關系密切。研究這些性狀間的相互關系以及對纖維品質的影響，有助于分清主次，明確主攻方向，對棉花雜交種的選育和豐產栽培有著重要的指導意義[29-31]。石建斌等[32]以陸地棉F5群體為研究對象發現，株高和果枝數對棉花纖維長度和纖維強度有較大影響。李慧琴等[21]以270份陸地棉種質資源進行研究發現，纖維品質性狀與株高均呈顯著正相關，表明株高與纖維品質關系密切。董承光等[7]對429份優異陸地棉種質進行表型性狀鑒定發現，株高與7個農藝性狀和5個纖維品質性狀均呈顯著正相關。羅海華等[33]對常規棉和雜交棉的農藝性狀與品質性狀進行相關性分析發現，在常規棉群體中纖維長度和纖維強度與株高和果枝數均呈顯著負相關，而雜交棉群體中則未檢測到相關性。馮文林等[34]研究發現，海島棉與陸地棉正反雜交群體的農藝性狀與纖維品質間均存在一定相關性，如衣分與絨長呈顯著負相關，馬克隆值與株高呈顯著正相關。由于試驗環境的不同和研究材料在株型、產量和品質等方面差異較大，本試驗相關性結果與前人研究結果并不完全一致，纖維長度與株高、有效果枝數、鈴數和有效鈴數呈顯著正相關；而馬克隆值與這些農藝性狀呈顯著負相關，與衣分呈極顯著正相關；纖維整齊度與株高、有效果枝數呈顯著正相關；纖維強度與始節高呈顯著正相關。由此表明，對農藝性狀提升的同時有利于改善棉花纖維品質。

主成分分析將海陸回交群體的性狀劃分為4個主成分，累計貢獻率達75.252%，分別為棉花成鈴數相關因子、纖維品質選擇因子、低品質控制因子和衣分與馬克隆值選擇的綜合因子。由于棉花農藝性狀之間也存在著不同程度的相關性，僅通過農藝性狀與纖維品質的相關分析和主成分分析，難以反映農藝性狀對棉花纖維品質性狀的作用，故本研究進一步進行了逐步回歸分析和通徑分析。通過多元逐步回歸分析，發現始節數和有效果枝數對纖維品質的影響較小，可剔除。鈴數和果枝數與棉花纖維長度具有較強相關性；始節高和鈴數與纖維強度有較強相關性；衣分、有效鈴數和始節高與馬克隆值有較強相關性；纖維整齊度只與株高有較強相關性。通徑分析進一步確定果枝數和鈴數對纖維長度分別產生負向直接效應和正向直接效應；鈴數和始節高對纖維強度產生正向直接效應；有效鈴數對馬克隆值產生負向直接效應，而衣分對其產生正向直接效應；株高對纖維整齊度產生正向直接效應。總體來看，有效鈴數增多，馬克隆值減小；而株高、鈴數和有效鈴數對纖維品質的間接作用占主導地位，與石建斌等[32]研究結果相一致。因此，在優質棉的選育中應兼顧農藝性狀和纖維品質，可通過對農藝性狀的選擇實現對纖維品質的間接選擇；海陸雜交可以作為棉花纖維品質改良的有效方法。