邵千順 楊琳

摘要 [目的]對(duì)鷹嘴豆種質(zhì)資源進(jìn)行篩選及多樣性分析。[方法]采用主成分分析、聚類分析等方法，通過5個(gè)產(chǎn)量相關(guān)因子研究鷹嘴豆遺傳多樣性并篩選優(yōu)質(zhì)種質(zhì)資源。[結(jié)果]多樣性指數(shù)最高的是百粒重，其次是粒型；性狀變異系數(shù)最大的是株型，其次是單株粒數(shù)。主成分分析表明，主要信息集中在6個(gè)主成分，其累計(jì)貢獻(xiàn)率達(dá)84.48%；通過聚類分析將50份鷹嘴豆在歐氏距離1.059 2處劃分為4類，第一類百粒重最大，即籽粒較大；第二類株高最高，產(chǎn)量居中；第三類各性狀都相對(duì)居中；第四類株高最低，單株莢數(shù)和單株粒數(shù)最大。[結(jié)論]該研究為鷹嘴豆的開發(fā)利用提供理論支撐。

關(guān)鍵詞 鷹嘴豆；種質(zhì)資源；多樣性

中圖分類號(hào) S529 文獻(xiàn)標(biāo)識(shí)碼 A 文章編號(hào) 0517-6611（2017）34-0016-04

Abstract [Objective] To screen the Cicer arietinum germplasm resources，and to carry out their diversity analysis.[Method] Five yield correlation factors were used to research the genetic diversity of C.arietinum and to screen the highquality germplasm resources by Principal Component Analysis，Cluster Analysis and so on.[Result] 100grain weight had the highest diversity index，followed with grain shape.Plant type had the highest coefficient of character variation，followed with grains per plant.Results of Principal Component Analysis showed that the major information was concentrated in 6 principal components，and thier accumulative contribution rate reached 84.48%.50 samples of C.arietinum were divided into four types at 1.059 2 Euclidean distance by Cluster Analysis.The first type had the greatest 100grain weight，and the grains were relatively big.The second type had the greatest plant height and middle yield.The characters of C.arietinum in the third type were in the middle.The fourth type had the lowest plant height，but the greatest pods per plant and grains per plant.[Conclusion] This research provides theoretical support for the development and utilization of C.arietinum.

Key words C.arietinum；Germplasm resources；Diversity

鷹嘴豆（Cicer arietinum）起源于亞洲西部和近東地區(qū)，是繼大豆之后第二大重要的谷物豆類[1]，鷹嘴豆耐旱[1-2]、耐鹽堿[3]、抗病和抗蟲，適宜在冷涼地區(qū)種植[4]，目前分布于全球 89.89%的地區(qū)，產(chǎn)量占這些地區(qū)的 86.73%。在我國(guó)，鷹嘴豆種植面積最大的區(qū)域是新疆[5]，其次在甘肅青海等地也有分布。

鷹嘴豆育種始于印度（1905年），之后全球多個(gè)國(guó)家先后開始鷹嘴豆育種，并在抗病、抗旱育種等方面取得了一定成果[6]。但我國(guó)在鷹嘴豆育種方面進(jìn)展緩慢。種質(zhì)資源是育種工作的基礎(chǔ)，了解資源遺傳多樣性、篩選優(yōu)質(zhì)種質(zhì)資源可為鷹嘴豆新品種選育提供親本選擇的基礎(chǔ)。韓文革等[7]利用聚類分析方法研究8個(gè)鷹嘴豆品種的親緣關(guān)系，并將8個(gè)品種聚為3類。張金波等[8]利用主成分分析、聚類分析等分析方法對(duì)新疆引進(jìn)的部分鷹嘴豆資源農(nóng)藝性狀進(jìn)行遺傳多樣性分析，發(fā)現(xiàn)鷹嘴豆資源各性狀遺傳多樣性指數(shù)變異較大，單株生物學(xué)產(chǎn)量的多樣性指數(shù)最高，但是不同材料間變異系數(shù)存在很大差別。這是由于農(nóng)藝性狀受環(huán)境影響較大，同一品種同一性狀在不同地區(qū)表現(xiàn)可能不同。

鷹嘴豆系寧夏地區(qū)首次引進(jìn)種植，在栽培、育種方面的研究還處于空白狀態(tài)。該研究以寧夏農(nóng)林科學(xué)院引進(jìn)保存的50份鷹嘴豆種質(zhì)資源為研究對(duì)象，利用主成分分析、相關(guān)性分析、聚類分析等方法對(duì)鷹嘴豆主要農(nóng)藝性狀間的遺傳關(guān)系及多樣性進(jìn)行分析，為深入開展鷹嘴豆資源新基因的挖掘與創(chuàng)新利用，以及在寧夏開展鷹嘴豆育種研究提供理論基礎(chǔ)。

1 材料與方法

1.1 試驗(yàn)地概況

2016年在寧夏農(nóng)林科學(xué)院固原分院頭營(yíng)試驗(yàn)基地進(jìn)行，試驗(yàn)基地海拔1 550 m，年降雨量420 mm，≥10 ℃積溫為2 691.9 ℃；年日照時(shí)數(shù)2 200～3 000 h，生育期日照時(shí)數(shù)1 030～1 120 h。

1.2 試驗(yàn)材料

50份鷹嘴豆種質(zhì)資源均引自中國(guó)農(nóng)業(yè)科學(xué)院國(guó)家種質(zhì)資源庫，供試材料的名稱及編號(hào)詳見表1。

1.3 試驗(yàn)方法 試驗(yàn)采用隨機(jī)區(qū)組設(shè)計(jì)，3次重復(fù)；每品種為1個(gè)小區(qū)，每小區(qū)3行，行長(zhǎng)6 m，行距50 cm，株距8～10 cm。

田間性狀調(diào)查參照《鷹嘴豆種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)》[9]進(jìn)行，主要考查株高、百粒重、單株莢數(shù)、單株粒數(shù)、單株粒質(zhì)5個(gè)農(nóng)藝性狀。待植株成熟時(shí)每小區(qū)取10株考種樣于室內(nèi)調(diào)查株高、百粒重、單株莢數(shù)、單株粒數(shù)、單株粒重5個(gè)性狀。

1.4 數(shù)據(jù)統(tǒng)計(jì)分析

采用Excel 2010統(tǒng)計(jì)處理基本數(shù)據(jù)，求性狀的平均值、最大值、最小值、標(biāo)準(zhǔn)差（SD）、變異系數(shù)（CV）；不同品種間性狀的差異用變異系數(shù)表示，遺傳多樣性指數(shù)的計(jì)算采用Shannon -Weaver 信息指數(shù)表示。計(jì)算公式：H′=-ΣPilnPi，其中Pi為某一性狀第i個(gè)級(jí)別出現(xiàn)的概率。計(jì)算多樣性指數(shù)時(shí)的劃級(jí)方法如下：先計(jì)算參試材料總體平均數(shù)（X）和標(biāo)準(zhǔn)差（σ），然后劃分為10級(jí)，從第1級(jí)[Xi<（X-2σ）]到第10級(jí)[Xi>（X+2σ）]，每0.5σ為1級(jí)。每一級(jí)的相對(duì)頻率用于計(jì)算多樣性指數(shù)。多樣性指數(shù)H′=-PilnPi，式中Pi為某性狀第i級(jí)別內(nèi)材料份數(shù)占總份數(shù)的百分比，ln為自然對(duì)數(shù)。利用DPS 9.50數(shù)據(jù)處理系統(tǒng)進(jìn)行鷹嘴豆品種的株高、百粒重、單株莢數(shù)、單株粒數(shù)、單株粒重等5個(gè)性狀的主成分分析；利用SAS 8.2軟件對(duì)50份鷹嘴豆種質(zhì)資源進(jìn)行聚類分析。

2 結(jié)果與分析

2.1 鷹嘴豆資源農(nóng)藝性狀指標(biāo)多樣性分析

參試資源5個(gè)產(chǎn)量相關(guān)因子多樣性指數(shù)最高的是百粒重（7.750 4 g），其次是單株莢數(shù)（7.387 1）；變異系數(shù)最大的是單株粒重（63.96%），其次是單株粒數(shù)（60.13%）。株高變異系數(shù)最小。單株粒重平均9.93 g，變幅2.28～24.50 g，主要分布在10.0～15.0 g，單株粒重超過20 g的資源數(shù)占供試資源總數(shù)的4.00%。單株粒數(shù)平均40.35個(gè)，變幅12.50～88.60個(gè)。參試資源單株粒數(shù)主要分布在30～45個(gè)，占25.20%，單株粒數(shù)超過75個(gè)的品種僅占參試資源的6.70%。百粒重平均24.06 g，變幅16.4～32.1 g。百粒重主要分布在20～30 g范圍內(nèi)，百粒重超過30 g的種質(zhì)占供試品種的11.92%。單株莢數(shù)平均56.43個(gè)，變幅2.3～113.3。品種單株莢數(shù)主要分布在45～75個(gè)，占參試資源總數(shù)的61.00%，單株莢數(shù)超過90個(gè)的品種占供試品種的4.60%。

2.2 鷹嘴豆資源農(nóng)藝性狀間的相關(guān)分析

將參試資源5個(gè)性狀兩兩之間的相關(guān)系數(shù)值列于表3中。對(duì)鷹嘴豆5個(gè)性狀的相關(guān)分析表明，植株株高與株型呈負(fù)相關(guān)，與單株莢數(shù)呈極顯著正相關(guān)。這表明植株越高，植株越呈現(xiàn)披散型，所結(jié)莢數(shù)也越多。株高、百粒重、單株莢數(shù)、單株粒數(shù)和單株粒重5個(gè)產(chǎn)量相關(guān)性狀間存在明顯的相關(guān)關(guān)系。百粒重與單株莢數(shù)和單株粒數(shù)存在極顯著負(fù)相關(guān)，與單株粒重呈極顯著正相關(guān)。這表明單株鷹嘴豆所結(jié)莢數(shù)越多，或單株上所產(chǎn)粒數(shù)越多，鷹嘴豆顆粒越小，百粒重越低。反之，籽粒越大，百粒重越高，單株所產(chǎn)粒質(zhì)量越大。單株莢數(shù)與單株粒數(shù)呈顯著正相關(guān)。單株粒數(shù)與單株粒重呈極顯著正相關(guān)，表明單株粒數(shù)越多，單株產(chǎn)量越高。

2.3 鷹嘴豆種質(zhì)主要農(nóng)藝性狀的主成分分析

2.3.1 主成分特征值。利用DPS軟件計(jì)算出5個(gè)主要農(nóng)藝性狀的特征向量及貢獻(xiàn)率（表4）。根據(jù)各向量的絕對(duì)值將不同性狀指標(biāo)劃分到不同的主成分中，同一指標(biāo)在各因子中的最大絕對(duì)值所在位置即為其所屬主成分。從表4可以看出，在所有主成分構(gòu)成中，主要信息集中在3個(gè)主成分，其累計(jì)貢獻(xiàn)率達(dá)90.62%。第1主成分特征值為2.26，貢獻(xiàn)率為45.13%；第2主成分特征值為1.20，貢獻(xiàn)率為23.97%；第3主成分特征值為1.08，貢獻(xiàn)率為21.53%。

2.3.2 主成分分析。3個(gè)主要成分因子對(duì)5個(gè)農(nóng)藝性狀的影響見表5。第1主成分主要影響產(chǎn)量性狀，單株粒數(shù)和單株粒重是主要指標(biāo)，向量值分別為0.60和0.61。第2主成分主要影響植株高度，向量值為0.82。第3主成分影響籽粒重量，百粒重是主要指標(biāo)，向量值為0.94。

2.4 鷹嘴豆種質(zhì)資源的聚類分析

利用 SAS 8.2軟件，采用參試資源形態(tài)性狀間的歐氏距離，對(duì)50份鷹嘴豆種質(zhì)資源進(jìn)行聚類分析。由表7、圖1可知，50份鷹嘴豆種質(zhì)資源在歐氏距離0.9655處聚為4類。第一類包含6份材料，單株粒數(shù)和單株莢數(shù)平均值最高，相應(yīng)的變異系數(shù)也較低，分別為12.05%、10.28%。第二類包含32份材料，百粒重最高，單株莢數(shù)和單株粒數(shù)次之，但變異系數(shù)較高，為15.60%。第三類包含10份材料，株高最高，百粒重和單株莢數(shù)次之，相應(yīng)變異系數(shù)最高。第四類包含2份材料，株高、百粒重和單株莢數(shù)等都相對(duì)較低。

3 結(jié)論與討論

種質(zhì)資源是育種工作的基礎(chǔ)，作物種子產(chǎn)業(yè)的每一次更新?lián)Q代，都離不開優(yōu)質(zhì)種質(zhì)資源的開發(fā)利用。因此，鷹嘴豆產(chǎn)業(yè)的發(fā)展，品種改良的關(guān)鍵依然是鷹嘴豆種質(zhì)資源的有效開發(fā)及利用。盡管各種同功酶標(biāo)記和DNA 分子標(biāo)記已經(jīng)被廣泛地應(yīng)用于植物種質(zhì)資源的鑒定和分類研究，但是農(nóng)藝性狀的鑒定和描述仍然是種質(zhì)資源研究的最基本的方法和途徑，農(nóng)藝性狀數(shù)據(jù)是種以上或種內(nèi)分類不可缺少的重要依據(jù)之一[10-12]。

該研究通過5個(gè)農(nóng)藝性狀的多樣性分析結(jié)果看出，參試50份鷹嘴豆資源多樣性豐富，遺傳背景廣泛、改良潛力較大，可為寧夏在鷹嘴豆新品種的選育，親本選配等方面提供優(yōu)異的種質(zhì)基礎(chǔ)。通過聚類分析將50份資源聚為4類，而這4類都有其特點(diǎn)：第一類百粒重最大，即籽粒較大；第二類株高最高，產(chǎn)量居中；第三類各性狀居中；第四類株高最低，單株莢數(shù)和單株粒數(shù)最大。從平均值和變異系數(shù)綜合考慮，第一類適合選育大粒型品種，第二類適合選育高桿兼高產(chǎn)型，第四類適合選育矮桿高產(chǎn)型。這一研究結(jié)果將為這50份鷹嘴豆資源的雜交組合選配、新品種的定向選育等提供理論依據(jù)。可根據(jù)不同的栽培方式或機(jī)械化程度，利用不同的類型，構(gòu)建不同的雜交組合，定向選育不同類型的鷹嘴豆品種。

遺傳多樣性反映的是生物種類遺傳信息的變化，可以表現(xiàn)在分子、細(xì)胞和個(gè)體等多個(gè)方面[13]。該研究?jī)H從5個(gè)農(nóng)藝性狀方面入手研究鷹嘴豆多樣性，結(jié)果顯示鷹嘴豆種質(zhì)資源遺傳多樣性豐富，但這一結(jié)果只是說明鷹嘴豆在農(nóng)藝性狀層面存在豐富的多樣性，因此有一定的局限性。為更準(zhǔn)確的揭示鷹嘴豆遺傳多樣性，還需從細(xì)胞學(xué)水平、生理生化水平及分子水平進(jìn)行多角度的深入研究，才能充分鑒定鷹嘴豆的遺傳多樣性。

參考文獻(xiàn)

[1] VARSHMEY R K，THUDI M，NAGAK S N，et al.Genetic dissection of drought tolerance in chickpea（Cicer arietinum L.）[J].

Theor Appl Genet，2014，127（11）：445-462.

[2] BHATTARAI T，F(xiàn)ETTIG S.Isolation and characterization of a dehydrin gene from Cicer pinnatifiduma，a droughtresistant wild relative of chickpea[J].Physiologia plantarum，2005，123（4）：452-458.

[3] SHUKLA R K，RAHA S，TRIPATHI V，et al.Expression of CAP2，an APETALA2-family transcription factor from chickpea，

enhances growth and tolerance to dehydration and salt stress in transgenic tobacco[J].Plant Physiol，2006，142（1）：113-123.

[4] 靳曉麗.不同鷹嘴豆材料ISSR遺傳多樣性分析及在蘭州地區(qū)的適應(yīng)性研究[D].蘭州：甘肅農(nóng)業(yè)大學(xué)，2013.

[5] 聶石輝，彭琳，王仙，等.鷹嘴豆種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析[J].植物遺傳資源學(xué)報(bào)，2015，16（1）：64-70.

[6] 黎太爵.食用鷹嘴豆及其優(yōu)良品種[J].食品工業(yè)科技，2001，22（4）：117-118.

[7] 韓文革，于曉春.聚類分析在鷹嘴豆農(nóng)藝性狀分類上的應(yīng)用[J].內(nèi)蒙古農(nóng)業(yè)科技，2006（6）：40-41，45.

[8] 張金波，李利民，苗昊翠，等.鷹嘴豆種質(zhì)資源主要農(nóng)藝性狀遺傳多樣性研究[J].新疆農(nóng)業(yè)科學(xué)，2014，51（1）：110-117.

[9] 宗緒曉，關(guān)建平，李玲，等.鷹嘴豆種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[M].北京：中國(guó)農(nóng)業(yè)科學(xué)技術(shù)出版社，2012.

[10] HE C，POYSA V，YU K.Development and characterization of simple sequence repeat（SSR）marker s and their use in determining relationships among Lycopersicon esculentum cultivars[J].Theor Appl Genet ，2003 ，106（2）：363-373 .

[11] ARESHCHENKOVA T，GANAL M W.Comparative analysis of polymorphism and chromosomal location of tomato microsatellite mark ersisolated from different sources[J].Theor Appl Genet，2002，104（2/3）：229-235.

[12] 王貴元，夏仁學(xué)，曾祥國(guó)，等.DMSO 處理對(duì)紅肉臍橙果肉主要色素和糖含量變化的影響及其相關(guān)性[J].西北農(nóng)業(yè)學(xué)報(bào)，2008，17（2）：165-168 .

[13] 楊忠芳.鷹嘴豆優(yōu)質(zhì)高產(chǎn)栽培技術(shù)[J].新疆農(nóng)業(yè)科技，2008（1）：12.