解慧芳 牛靜 邢璐 王素英 王磊 宋慧 蔣自可

摘要：為研究和創(chuàng)制谷子種質(zhì)資源，加深對種質(zhì)進化和分類的認識，以便進一步為遺傳育種和性狀改良提供依據(jù)，以117份國內(nèi)外谷子種質(zhì)資源為材料，開展13個表型性狀多樣性鑒定，并進行了變異系數(shù)分析、相關(guān)性分析、主成分分析和聚類分析。結(jié)果表明，供試谷子種質(zhì)表型性狀變異較豐富，13個性狀變異系數(shù)范圍為7.00%～74.11%，變異系數(shù)最大的為單株粒質(zhì)量（74.11%），表明谷子單株粒質(zhì)量變異最為豐富。相關(guān)性分析表明，單株穗質(zhì)量、主莖單穗質(zhì)量和主莖單穗粒質(zhì)量與葉面積、穗粗呈極顯著正相關(guān)關(guān)系。主成分分析把13個性狀綜合為5個主成分，累計貢獻率達到76.76%，表明5個主成分反映了谷子種質(zhì)資源大部分的性狀信息。第1主成分主要為產(chǎn)量相關(guān)性狀，表明種質(zhì)間產(chǎn)量性狀差異很豐富，能夠為選育高產(chǎn)品種提供豐富的材料來源。聚類分析將117份谷子種質(zhì)分為5類，第二大類占比最大，占總數(shù)的62%，通過比較發(fā)現(xiàn)其生育期、株高和節(jié)數(shù)等性狀較高，特別是主莖單穗質(zhì)量和主莖單穗粒質(zhì)量，可以作為以生物產(chǎn)量為目標性狀的育種材料。

關(guān)鍵詞：谷子；種質(zhì)資源；表型性狀；遺傳多樣性；變異系數(shù)分析；相關(guān)性分析；主成分分析；聚類分析

中圖分類號：S515.02? 文獻標志碼：A

文章編號：1002-1302（2023）13-0076-06

谷子別稱粟，屬禾本科狗尾草屬，是起源于我國的二倍體自花授粉作物，在燦爛的中華農(nóng)耕文明中起著舉足輕重的作用［1］。谷子不僅具有抗旱耐瘠薄等種植特點，還有豐富的營養(yǎng)成分，具有健脾胃和補益虛損等保健功效［2-5］，符合我國居民飲食結(jié)構(gòu)需求，在人民膳食健康的需求中發(fā)揮著重要的作用。隨著谷子種業(yè)的發(fā)展，谷子品種遺傳改良、種子繁育及品種推廣成為谷子發(fā)展的重點研究方向。目前，我國保存有多種谷子種質(zhì)資源，占世界谷子種質(zhì)資源保有量的80%左右，是占有量最多、研究最深入的國家［6-7］。豐富的種質(zhì)資源是品種改良的物質(zhì)基礎(chǔ)，種質(zhì)資源研究為雜交育種提供了基礎(chǔ)，遺傳多樣性分析和相關(guān)性研究是種質(zhì)資源常用的方法。

谷子種質(zhì)多樣性的研究主要集中在產(chǎn)量相關(guān)性狀，穗長、穗粗、單穗質(zhì)量、穗粒質(zhì)量、千粒質(zhì)量等是影響產(chǎn)量的重要因素［8-9］。趙利蓉等研究發(fā)現(xiàn)不同生態(tài)區(qū)16個谷子品種農(nóng)藝性狀的變異幅度為10.79%～32.61%［10］。魏萌涵等對華北夏谷區(qū)谷子農(nóng)藝性狀的研究表明，產(chǎn)量與株高、穗粒質(zhì)量、單穗質(zhì)量、出谷率呈正相關(guān)，與生育期，穗長、穗粗、千粒質(zhì)量呈負相關(guān)［11］。聚類分析也廣泛應用于谷子多樣性的劃分，相關(guān)表型性狀的分類可以為種質(zhì)資源的利用提供參考［12-15］。

了解種質(zhì)遺傳變異的大小以及性狀之間的關(guān)系，能夠加深種質(zhì)進化和分類的認識，對遺傳育種和性狀改良都有重要意義。本研究以不同來源的117份谷子核心種質(zhì)資源為研究對象，對其表型性狀進行遺傳多樣性分析、相關(guān)性分析和聚類分析，并對其進行全面評價，以期為谷子資源保護和利用奠定基礎(chǔ)，為種質(zhì)資源創(chuàng)新利用和雜交親本選配提供參考。

1 材料與方法

1.1 試驗材料

供試材料117份，其中野生品種5份、農(nóng)家品種56份、育成品種36份、國外品種20份，由中國農(nóng)業(yè)科學院作物科學研究所提供（表1）。

1.2 試驗設(shè)計

試驗于2020年在河南省安陽市北關(guān)區(qū)柏莊鎮(zhèn)安陽市農(nóng)業(yè)科學院柏莊試驗基地進行，5月20日播種，小區(qū)面積為2 m2，種植密度為60萬株/hm2，行長為 1 m，行距為40 cm，2行區(qū)種植，采用隨機區(qū)組排列。該基地地勢開闊平坦，田間排灌設(shè)施齊全，無建筑物、樹木等遮陰，田間管理根據(jù)當?shù)卮筇锷a(chǎn)習慣實施。

1.3 試驗方法

每個小區(qū)選取分布均勻且長勢一致的10株材料進行表型性狀調(diào)查，根據(jù)谷子種質(zhì)資源描述規(guī)范和數(shù)據(jù)標準對生育期、節(jié)數(shù)、株高、主穗長度、葉面積、穗粗、主穗碼數(shù)、草質(zhì)量、單株穗質(zhì)量、主莖單穗質(zhì)量、主莖單穗粒質(zhì)量、單株粒質(zhì)量、千粒質(zhì)量共13個表型性狀進行田間表型性狀調(diào)查。

1.4 數(shù)據(jù)分析

采用Excel 2010進行數(shù)據(jù)整理和變異分析，使用SPSS Statistics 23 軟件進行相關(guān)性分析和主成分分析，利用Origin 2018進行聚類分析。

2 結(jié)果與分析

2.1 117份谷子表型性狀的多樣性比較

由表2可知，117份材料的13個表型性狀的變異系數(shù)范圍為7.00%～74.11%，變異系數(shù)最大的為單株粒質(zhì)量，最小的是生育期，均值為38.5%。其中單株穗質(zhì)量、主莖單穗質(zhì)量、主莖單穗粒質(zhì)量、單株粒質(zhì)量的變異系數(shù)均高于50%，說明穗質(zhì)量、粒質(zhì)量相關(guān)表型性狀表現(xiàn)出更好的變異性；節(jié)數(shù)、株高、穗粗、千粒質(zhì)量、主穗長度、葉面積、主穗碼數(shù)變異系數(shù)在20%～40%之間，草質(zhì)量的變異系數(shù)為48.08%。以上分析表明117份谷子種質(zhì)資源具有豐富的遺傳多樣性，品種間性狀差異較大，適宜進行品種比較分析和品種鑒定評價。

2.2 117份谷子表型性狀的相關(guān)性分析

相關(guān)性分析是檢驗性狀間相關(guān)程度的重要方法，對育種性狀的選擇具有重要意義。由表3可知，谷子生育期、節(jié)數(shù)、株高等13個表型性狀之間存在不同程度相關(guān)性。其中生育期、節(jié)數(shù)、株高、主穗長度、葉面積和草質(zhì)量之間呈極顯著正相關(guān)；產(chǎn)量構(gòu)成因素單株穗質(zhì)量、主莖單穗質(zhì)量和主莖單穗粒質(zhì)量與葉面積、穗粗呈極顯著正相關(guān)關(guān)系；單株穗質(zhì)量和主莖單穗質(zhì)量越高單株粒質(zhì)量和千粒質(zhì)量越高，且達到極顯著水平。與生育期、株高、葉面積相關(guān)性最高的均為節(jié)數(shù)，其中與節(jié)數(shù)相關(guān)性最高的為株高；與主穗長度相關(guān)性最高的為生育期；與穗粗相關(guān)性最高的為主莖單穗粒質(zhì)量；與主穗碼數(shù)相關(guān)性最高的為穗粗；與草質(zhì)量相關(guān)性最高的為單株穗質(zhì)量，與單株穗質(zhì)量相關(guān)性最高的為主莖單穗質(zhì)量；與單株粒質(zhì)量相關(guān)性最高的為主莖單穗粒質(zhì)量。單株穗質(zhì)量與主莖單穗質(zhì)量的相關(guān)性系數(shù)接近于1，這主要是因為谷子植株多為單株，少有分蘗。

2.3 117份谷子表型性狀的主成分分析

本研究對117份谷子品種的13個表型性狀進行主成分分析（表4），提取了5個主要成分，結(jié)果顯示，前5個主成分的累計貢獻率達76.76%，它們的特征值均大于1，表明前5個主成分代表了谷子表型性狀的大部分信息，可以作為谷子性狀評價的綜合指標。第1主成分的特征值為4.536，其貢獻率最大，為34.89%，其中，主莖單穗質(zhì)量、主莖單穗粒質(zhì)量、單株穗質(zhì)量分值較高，說明第1主成分主要是與產(chǎn)量穗質(zhì)量相關(guān)的性狀。第2主成分特征值為2.191，貢獻率為16.85%，節(jié)數(shù)、葉面積、株高和生育期分值較高，說明第2主成分主要是谷子農(nóng)藝性狀的綜合反映。第3主成分特征值為1.190，貢獻率為9.15%，主穗碼數(shù)和穗粗分值較高，反映了谷子穗型相關(guān)性狀。第4主成分特征值為1.062，貢獻率為8.17%，千粒質(zhì)量的分值較高。第5主成分特征值為1.001，貢獻率為7.70%，主穗長度的分值較高?？偟膩碚f， 前2個主成分主要反映的是谷子產(chǎn)量和農(nóng)藝性狀有關(guān)的情況。

2.4 117份谷子表型性狀的聚類分析

根據(jù)性狀變異對117份谷子種質(zhì)進行系統(tǒng)的聚類分析（圖1），分析發(fā)現(xiàn)種質(zhì)聚類與地里位置關(guān)系不明顯，主要由表型性狀進行聚類。根據(jù)表型性狀可將參試種質(zhì)分為五大類。第1類包括37個品種，占總數(shù)的32%左右。第2類包括73個品種，占總數(shù)的62%左右，第3類包括1個品種，第4類包括5個品種，第5類包括1個品種。第2類品種的生育期、株高和節(jié)數(shù)等性狀高于第1類品種，尤其是主莖單穗質(zhì)量和主莖單穗粒質(zhì)量差異較大，結(jié)果說明第2類品種可以作為選育高產(chǎn)的目標親本加以利用。由聚類分析結(jié)果可知，各類種質(zhì)資源品種類型較為混亂，這可能是近年來育種材料的來源豐富導致，育種過程中加強了農(nóng)家種、國外品種及野生種的利用。

3 討論與結(jié)論

作物表型性狀是種質(zhì)遺傳多樣性比較直接有效的性狀指標，為選育品種提供一定的數(shù)據(jù)參考。通過對不同谷子種質(zhì)遺傳多樣性分析發(fā)現(xiàn)，品種間性狀變異廣泛，但是其變異系數(shù)較低，其中單株穗質(zhì)量和單株粒質(zhì)量的變異系數(shù)相對較［8，15-19］。本研究通過對117份谷子材料的13個表型性狀進行分析也得到了類似的研究結(jié)果，不同谷子間的表型性狀變異較為豐富，其中主莖單穗粒質(zhì)量和單株粒質(zhì)量的變異系數(shù)較大，分別為68.44%、74.11%，結(jié)果表明，盡管種子產(chǎn)量形狀受環(huán)境因素的影響較大，但是種子形成和分配主要取決于作物自身的遺傳特性，相對穩(wěn)定。其他表型性狀的變異范圍在7.00%～59.83%之間，較王海崗等對878份谷子的研究變異性較大，遺傳多樣性更為豐富［16］。谷子品種變異較為豐富，不僅能夠較好地代表谷子資源的遺傳多樣性，在優(yōu)質(zhì)品種篩選和培育中也有較好的利用價值。

本研究對117份谷子材料13個表型性狀進行了相關(guān)分析，發(fā)現(xiàn)生育期、節(jié)數(shù)、株高、主穗長度、葉面積和草質(zhì)量之間呈極顯著正相關(guān)，且單株穗質(zhì)量、主莖單穗質(zhì)量和主莖單穗粒質(zhì)量與葉面積、穗粗呈極顯著正相關(guān)關(guān)系。趙芳等研究表明單株穗粒質(zhì)量與穗質(zhì)量、單株草質(zhì)量呈極顯著正相關(guān)關(guān)系［15］，王歡歡等研究也表明穗粒質(zhì)量與單穗質(zhì)量、穗粗呈正相關(guān)關(guān)系［20］，但魏萌涵等研究表明，穗粗與單穗質(zhì)量和穗粒質(zhì)量呈負相關(guān)關(guān)系［11］，這可能與研究所選用的材料來源不同相關(guān)。

主成分分析和聚類分析是評估作物遺傳多樣性和篩選優(yōu)異種質(zhì)資源中普遍的一種分析方法，通過此方法能夠反映作物種質(zhì)間的遺傳差異和相關(guān)性狀特點，在作物中已經(jīng)得到廣泛應用［8，16，21-26］。本研究通過對13個表型性狀進行主成分分析，提取了前5個主成分，反映了谷子主要表型性狀76.76%的遺傳信息，第1主成分主要為產(chǎn)量相關(guān)性狀，表明種質(zhì)間產(chǎn)量性狀差異交豐富，能夠為選育高產(chǎn)品種提供豐富的材料來源，在前人研究中也得到相似的研究結(jié)果［13，15］。根據(jù)谷子的表型性狀進行聚類分析，將種質(zhì)聚類為五大類，第二大類品種，占總數(shù)的62%左右，主要由農(nóng)家品種和育成品種構(gòu)成，通過比較發(fā)現(xiàn)其生育期、株高和節(jié)數(shù)等性狀較高，特別是主莖單穗質(zhì)量和主莖單穗粒質(zhì)量，可以作為以生物產(chǎn)量為目標性狀的育種材料。此外，本研究發(fā)現(xiàn)在第2聚類中河北省的品種較多，可能由于地域位置和環(huán)境相似，有較好的育種和研究價值［27］。

通過對117份谷子的13個表型性狀進行統(tǒng)計分析，表明不同種質(zhì)間存在較大的變異系數(shù)和豐富的遺傳多樣性。通過主成分分析和聚類分析提取了主要主成分和品種聚類，為核心種質(zhì)的應用以及新品種的選育和培育提供了參考。

參考文獻：

［1］Yang X Y，Wan Z W，Perry L，et al. Early millet use in Northern China［J］. Proceedings of the National Academy of Sciences of the United States of America，2012，109（10）：3726-3730.

［2］章潔瓊，劉 輝，鄒 軍，等. 不同小米品種理化品質(zhì)及感官品質(zhì)分析與評價［J］. 中國糧油學報，2022，37（5）：70-78.

［3］王軍鋒，周顯青，張玉榮. 小米的營養(yǎng)特性與保健功能及產(chǎn)品開發(fā)［J］. 糧食加工，2012，37（3）：60-63.

［4］李暮男，蘭鳳英. 小米的營養(yǎng)成分及保健功能研究進展［J］. 河北北方學院學報（自然科學版），2017，33（7）：56-60.

［5］劉宇杰，陳銀煥，楊修仕，等. 小米營養(yǎng)及功能成分研究進展［J］. 糧食與油脂，2020，33（5）：1-3.

［6］劉 旭，黎 裕，曹永生，等. 中國禾谷類作物種質(zhì)資源地理分布及其富集中心研究［J］. 植物遺傳資源學報，2009，10（1）：1-8，171，169.

［7］賈冠清，刁現(xiàn)民. 谷子［Setaria italica （L.） P.Beauv.］作為功能基因組研究模式植物的發(fā)展現(xiàn)狀及趨勢［J］. 生命科學，2017，29（3）：292-301.

［8］田伯紅. 谷子地方品種和育成品種的遺傳多樣性研究［J］. 植物遺傳資源學報，2010，11（2）：224-228.

［9］楊天育，黃相國，何繼紅，等. 谷子遺傳資源多樣性研究進展［J］. 西北農(nóng)業(yè)學報，2003，12（1）：43-47.

［10］趙利蓉，馬 珂，張麗光，等. 不同生態(tài)區(qū)谷子品種農(nóng)藝性狀和品質(zhì)分析［J］. 作物雜志，2022（2）：44-53.

［11］魏萌涵，解慧芳，邢 璐，等. 華北地區(qū)谷子產(chǎn)量與農(nóng)藝性狀的綜合評價分析［J］. 作物雜志，2018（4）：42-47.

［12］相吉山，張恒儒，劉 涵，等. 不同生態(tài)區(qū)谷子種質(zhì)資源表型比較分析［J］. 中國農(nóng)業(yè)科技導報，2020，22（9）：31-41.

［13］任芹勇，樊巧利，李 濤，等. 65份谷子品種農(nóng)藝性狀聚類和相關(guān)性［J］. 分子植物育種，2017，15（12）：5178-5188.

［14］丁銀燈，聶石輝，王 仙，等. 谷子主要育成品種在新疆的遺傳多樣性研究［J］. 植物遺傳資源學報，2018，19（2）：232-242.

［15］趙 芳，魏 瑋，張曉磊，等. 224個谷子品種農(nóng)藝性狀聚類和相關(guān)性分析［J］. 種子，2022，41（1）：74-83.

［16］王海崗，賈冠清，智 慧，等. 谷子核心種質(zhì)表型遺傳多樣性分析及綜合評價［J］. 作物學報，2016，42（1）：19-30.

［17］王曉娟，祁旭升，王興榮，等. 甘肅省谷子地方種質(zhì)資源遺傳多樣性分析［J］. 干旱地區(qū)農(nóng)業(yè)研究，2009，27（6）：129-133，153.

［18］孫加梅，王雪梅，王東建，等. 谷子種質(zhì)資源遺傳多樣性研究［J］. 山東農(nóng)業(yè)科學，2013，45（3）：33-37，60.

［19］楊慧卿，王根全，郝曉芬，等. 山西省谷子地方種質(zhì)資源表型多樣性分析［J］. 江蘇農(nóng)業(yè)科學，2022，50（13）：20-25.

［20］王歡歡，賈亞濤，許鵬飛，等. 谷子主要農(nóng)藝性狀與產(chǎn)量的關(guān)系［J］. 山西農(nóng)業(yè)科學，2014，42（7）：657-659，754.

［21］胡標林，萬 勇，李 霞，等. 水稻核心種質(zhì)表型性狀遺傳多樣性分析及綜合評價［J］. 作物學報，2012，38（5）：829-839.

［22］魏添梅，昌小平，閔東紅，等. 小麥抗旱品種的遺傳多樣性分析及株高優(yōu)異等位變異挖掘［J］. 作物學報，2010，36（6）：895-904.

［23］姚堅強，鮑堅東，朱金慶，等. 中國糯玉米wx基因種質(zhì)資源遺傳多樣性［J］. 作物學報，2013，39（1）：43-49.

［24］王 瑗，郁萬文，周 凱，等. 苦豆子表型性狀多樣性分析及綜合評價［J］. 廣西植物，2019，39（9）：1147-1158.

［25］李洪立，胡文斌，洪青梅，等. 火龍果種質(zhì)資源果實特性的遺傳多樣性分析［J］. 熱帶亞熱帶植物學報，2019，27（4）：432-438.

［26］李 偉，王 攀，其其格，等. 藍莓種質(zhì)資源表型多樣性研究［J］. 北京林業(yè)大學學報，2020，42（2）：124-134.

［27］祁 娟，曹文俠，閆偉紅. 披堿草屬野生居群表型多樣性及其與環(huán)境關(guān)系研究［J］. 西北植物學報，2013，33（5）：1027-1033.