邢 媛 ，賈馥翠 ，陳 凌 ，王海崗 ，王計(jì)平 ，喬治軍 ，王瑞云

（1.山西農(nóng)業(yè)大學(xué) 農(nóng)學(xué)院，山西 太谷 030801；2.山西農(nóng)業(yè)大學(xué) 農(nóng)業(yè)基因資源研究中心/農(nóng)業(yè)農(nóng)村部黃土高原作物基因資源與種質(zhì)創(chuàng)制重點(diǎn)實(shí)驗(yàn)室/雜糧種質(zhì)資源發(fā)掘與遺傳改良山西省重點(diǎn)實(shí)驗(yàn)室，山西太原 030031）

干旱是農(nóng)作物減產(chǎn)的重要原因之一。全球的干旱半干旱地區(qū)相對(duì)于陸地面積的占比高達(dá)35%，涉及到全球60 多個(gè)國(guó)家和地區(qū)[1]。干旱影響植物的正常生長(zhǎng)和發(fā)育，進(jìn)而導(dǎo)致農(nóng)作物減產(chǎn)嚴(yán)重[2]，在所有非生物對(duì)農(nóng)作物造成的損失中，干旱脅迫占首位[3]。農(nóng)作物具有抵御和適應(yīng)干旱環(huán)境的能力，但是不同農(nóng)作物間及作物的不同品種之間的抗旱能力和抗旱機(jī)制存在著顯著差異[4]。隨著全球溫室效應(yīng)的不斷影響，溫室氣體濃度的增加促使地球表面溫度加快上升，進(jìn)而導(dǎo)致土地干旱程度加重、沙漠化面積增大，這對(duì)農(nóng)作物的生產(chǎn)將造成很大的影響，而選育抗旱品種是緩解干旱的有效策略[5]。

黍稷是起源于我國(guó)的古老作物，抗旱耐瘠，營(yíng)養(yǎng)豐富[6]。黍稷適應(yīng)土壤能力強(qiáng)，在多數(shù)土壤中均能正常生長(zhǎng)[7]。黍稷種子吸收25%自身干質(zhì)量的水分就可發(fā)芽，而小麥和玉米分別需要44% 和45%。當(dāng)黍稷苗期遭遇高溫干旱，葉片會(huì)卷縮成筒狀，氣孔閉合，水分蒸發(fā)量減少，推遲或停止生長(zhǎng)，呈萎蔫狀以適應(yīng)干旱和減少失水造成的危害。不同農(nóng)作物間及作物的不同品種之間的抗旱能力和抗旱機(jī)制具有很大的差異。此前，人們關(guān)于黍稷抗旱性的研究層次遠(yuǎn)不如其他主要糧食作物如玉米、水稻、大豆、小麥等。山侖等[8]比較了黍稷、谷子、高粱和玉米在各個(gè)時(shí)期的抗旱性，結(jié)果顯示，幾種作物在芽的生長(zhǎng)階段對(duì)水分脅迫的敏感度最高；王綸等[9]研究表明，黍稷資源的抗旱性在每個(gè)個(gè)體中都表現(xiàn)出比較大的差異[10]；張盼盼等[11]用模擬干旱脅迫的方法，對(duì)黍稷芽期抗旱性的指標(biāo)進(jìn)行鑒定，篩選出芽期抗旱性指標(biāo)為相對(duì)根長(zhǎng)、相對(duì)芽長(zhǎng)、相對(duì)發(fā)芽率、相對(duì)芽鞘長(zhǎng)和相對(duì)芽干質(zhì)量。張美俊等[12]研究發(fā)現(xiàn)，干旱脅迫抑制糜子種子萌發(fā)活力，降低了萌發(fā)速度和發(fā)芽率。王海茹等[13]研究認(rèn)為，在干旱脅迫的情況下，適當(dāng)?shù)脑鍪┑士梢蕴岣呙缙谑蝠⒌娜~綠素含量和根系活力，增加根系總表面積、總體積與總根長(zhǎng)，降低根系丙二醛含量，在一定程度上緩解干旱脅迫的影響。李鵬等[14]研究發(fā)現(xiàn)，低磷脅迫下，耐低磷品種、磷中間型品種植株相對(duì)吸磷量均顯著高于磷敏感型品種。李占成等[15]研究發(fā)現(xiàn)，鹽脅迫條件下糜子種子的相對(duì)發(fā)芽勢(shì)和相對(duì)發(fā)芽指數(shù)均顯著低于對(duì)照，各種鹽對(duì)糜子種子發(fā)芽抑制作用由大到小為NaCl、Na2SO4、MgCl2。閆江艷等[16]研究發(fā)現(xiàn)，干旱脅迫下引起2個(gè)黍稷品種各項(xiàng)生理指標(biāo)及復(fù)水條件下變化幅度不同，2 個(gè)品種根系活力均明顯下降，但抗旱性強(qiáng)的隴糜4號(hào)根系活力下降幅度明顯低于抗旱性弱的晉黍7號(hào)。不同植物的萌發(fā)特性及干旱脅迫下的反應(yīng)都有自身的特點(diǎn)，尚無(wú)統(tǒng)一的評(píng)價(jià)方法和體系。在實(shí)際生產(chǎn)過(guò)程中，輕度干旱與重度干旱對(duì)不同作物種子萌發(fā)具有不同的影響，表明種子的活力及綜合性的抗旱評(píng)價(jià)至關(guān)重要。目前，干旱對(duì)黍稷種子活力影響、抗旱性指標(biāo)篩選及抗旱性評(píng)價(jià)研究較少。

本研究以東北地區(qū)35 份黍稷資源為試驗(yàn)材料，在PEG 脅迫的干旱條件下，綜合分析黍稷在干旱脅迫下2 個(gè)重要的生長(zhǎng)發(fā)育時(shí)期的發(fā)芽勢(shì)、發(fā)芽率、根長(zhǎng)、根鮮質(zhì)量、根干質(zhì)量、苗長(zhǎng)、苗鮮質(zhì)量、苗干質(zhì)量、葉綠素含量，進(jìn)而對(duì)比它們之間的抗旱性差異，建立抗旱性評(píng)價(jià)體系，旨在篩選出黍稷抗旱資源。

1 材料和方法

1.1 試驗(yàn)材料

本研究所用的35份北方黍稷種質(zhì)資源材料包含黑龍江19份、吉林10份、遼寧6份（表1）。

表1 參試品種及來(lái)源Tab.1 The test varieties and their origin

1.2 試驗(yàn)方法

從35 份黍稷品種中分別選取180 粒飽滿種子，用0.3%的高錳酸鉀溶液消毒15 min，再用無(wú)菌水沖洗5 遍。試驗(yàn)設(shè)置濃度25%的PEG-6000（D）和正常供水（CK）2 個(gè)處理，每個(gè)處理3 次重復(fù)，每個(gè)重復(fù)30 粒種子，每個(gè)培養(yǎng)皿為1 次重復(fù)。在每個(gè)無(wú)菌培養(yǎng)皿中放入2 層濾紙，放入30 粒消毒充分的黍稷種子，其中3 個(gè)培養(yǎng)皿中加入5 mL 配制好的PEG 溶液進(jìn)行干旱脅迫，剩余3 個(gè)培養(yǎng)皿正常供水作為對(duì)照組，最后將不同處理的培養(yǎng)皿隨機(jī)放入25 ℃恒溫光照培養(yǎng)箱（光照∶黑暗=12∶12）中進(jìn)行種子萌發(fā)試驗(yàn)。為保持脅迫濃度相對(duì)穩(wěn)定，發(fā)芽期間每4 d 更換濾紙和PEG 脅迫溶液。

1.3 測(cè)定指標(biāo)

在種子萌發(fā)階段分別統(tǒng)計(jì)不同處理下各處理組的發(fā)芽勢(shì)（I1）、發(fā)芽率（I2），計(jì)算種子的相對(duì)發(fā)芽勢(shì)（I1′）、相對(duì)發(fā)芽率（I2′）；在苗期試驗(yàn)結(jié)束后，每個(gè)育苗盤中挑選生長(zhǎng)一致幼苗，測(cè)量其根長(zhǎng)（C1）、苗長(zhǎng)（C2），并測(cè)定根鮮質(zhì)量（C3）、苗鮮質(zhì)量（C4）、根干質(zhì)量（C5）、苗干質(zhì)量（C6）、葉綠素含量（C7）。

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

采用Microsoft Excel 2010 和SPSS 25.0 進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析及作圖。

2 結(jié)果與分析

2.1 黍稷芽、苗期對(duì)干旱脅迫的響應(yīng)

PEG-6000 脅迫對(duì)不同品種黍稷種子萌發(fā)的影響如表2 所示。

表2 不同品種黍稷芽、苗期PEG 脅迫較對(duì)照組的差值Tab.2 The difference of PEG stress between different varieties of broomcorn millet at bud and seedling stages compared with the control group

續(xù)表2 不同品種黍稷芽、苗期PEG 脅迫較對(duì)照組的差值Tab.2（Continued） The difference of PEG stress between different varieties of broomcorn millet at bud and seedling stages compared with the control group

從表2 可以看出，與對(duì)照組相比，PEG 脅迫下9 個(gè)指標(biāo)數(shù)值均有所下降，其中，發(fā)芽勢(shì)和發(fā)芽率下降的尤為明顯，分別下降93.007%和94.449%，可見(jiàn)各指標(biāo)均受PEG 脅迫的抑制，其中，發(fā)芽勢(shì)和發(fā)芽率受抑制較明顯。從表3 可以看出，通過(guò)成對(duì)樣本t檢驗(yàn)發(fā)現(xiàn)，CK 與PEG 脅迫間t值變化范圍為［-0.234，-76.479］，說(shuō)明35 份黍稷種質(zhì)資源中參與測(cè)試的品種彼此之間差異較大，性狀變化明顯，有助于抗旱性檢測(cè)。其中，發(fā)芽勢(shì)和發(fā)芽率的t值分別為-51.423 和-76.479，t值的絕對(duì)值越大表明材料間差異越顯著，所以，PEG 干旱脅迫對(duì)種子芽期的生長(zhǎng)形成明顯抑制。P值多數(shù)小于0.05，表明CK與處理間差異顯著。測(cè)定PEG 脅迫材料苗期各性狀衡量指標(biāo)，得出其變異系數(shù)介于0.171～0.664，CK 組 變異 系 數(shù) 為0.154～0.612，試驗(yàn) 組 和CK 變異系數(shù)相差0.017～0.051，跨度較大，表明試驗(yàn)中35 份材料的每項(xiàng)指標(biāo)對(duì)PEG 脅迫較為敏感，具有代表性。

表3 25%PEG-6000 脅迫下35 份黍稷品種各抗旱指標(biāo)統(tǒng)計(jì)分析Tab.3 Statistical analysis of drought resistance indexes of 35 broomcorn millet varieties under 25% PEG-6000 stress

2.2 PEG 脅迫下黍稷各指標(biāo)的相關(guān)性分析

對(duì)PEG 脅迫下的芽期與苗期的9 項(xiàng)相對(duì)值指標(biāo)進(jìn)行Pearson 相關(guān)性分析可知（表4），芽期的相對(duì)發(fā)芽勢(shì)（I1′）和相對(duì)發(fā)芽率（I2′）與苗期的7 個(gè)指標(biāo)之間相關(guān)性均不顯著，與相對(duì)葉綠素含量呈正相關(guān)，與其余6 個(gè)指標(biāo)呈負(fù)相關(guān)；苗期7 個(gè)指標(biāo)兩兩之間均呈極顯著正相關(guān)，相對(duì)苗鮮質(zhì)量（C4′）與相對(duì)苗干質(zhì)量（C6′）相關(guān)系數(shù)最大，為0.934，二者相關(guān)性最高；相對(duì)苗長(zhǎng)（C2′）和相對(duì)苗鮮質(zhì)量（C4′）間的相關(guān)系數(shù)為0.897，相關(guān)性次之；相對(duì)根干質(zhì)量（C5′）與相對(duì)根長(zhǎng)（C1′）和相對(duì)根鮮質(zhì)量（C3′）間的相關(guān)系數(shù)分別為0.474 和0.823。通過(guò)數(shù)據(jù)可以得出，苗期各項(xiàng)指標(biāo)結(jié)果與芽期各項(xiàng)指標(biāo)結(jié)果間無(wú)相關(guān)性，相對(duì)葉綠素含量與苗期及芽期階段的相關(guān)性均較低。綜上所述，多數(shù)指標(biāo)間都呈現(xiàn)極顯著相關(guān)，說(shuō)明在PEG 脅迫下對(duì)黍稷的相對(duì)值作用效果明顯，但是直接利用這些指標(biāo)對(duì)黍稷做出抗旱性評(píng)價(jià)太過(guò)單一不夠準(zhǔn)確，為了得到更為準(zhǔn)確的結(jié)果，采用主成分分析法對(duì)黍稷的抗旱性進(jìn)行綜合性評(píng)價(jià)。

表4 干旱脅迫下各指標(biāo)的相關(guān)系數(shù)Tab.4 Correlation coefficients of each index under drought stress

2.3 黍稷苗期抗旱性分析

用主成分分析法對(duì)苗期相關(guān)抗旱性指標(biāo)的相對(duì)值進(jìn)行分析，篩選出可用于黍稷苗期抗旱性鑒定的主要指標(biāo)。由表5 可知，將特征值大于1 的指標(biāo)納入主成分，獲得2 個(gè)主成分。前2 個(gè)主成分貢獻(xiàn)率分別為66.776%和14.907%，累計(jì)貢獻(xiàn)率為81.683%（符合主成分分析累計(jì)貢獻(xiàn)率大于80%的要求），說(shuō)明前2個(gè)主因子所包含的要素信息量可以反映出9 個(gè)抗旱性指標(biāo)原始特征參數(shù)的大部分信息，因此，9 個(gè)抗旱性指標(biāo)可用前2 個(gè)主成分因子表示。

表5 主成分分析中各主成分的特征值及貢獻(xiàn)率Tab.5 Eigenvalues and contribution rates of each principal component in principal component analysis

由于因子載荷陣是不唯一的，所以對(duì)因子載荷矩陣進(jìn)行旋轉(zhuǎn)，使因子載荷陣結(jié)構(gòu)簡(jiǎn)化，運(yùn)用方差最大法得到各主成分的因子載荷矩陣（表6）：決定第1 主成分主要為相對(duì)苗長(zhǎng)（0.939）和相對(duì)苗鮮質(zhì)量（0.950）2 個(gè)指標(biāo)，說(shuō)明黍稷苗期生長(zhǎng)情況可以作為評(píng)價(jià)抗旱性的主要指標(biāo)，決定第2 主成分的主要為相對(duì)葉綠素含量（0.946），說(shuō)明黍稷的抗旱性強(qiáng)弱能夠影響植株光合作用。根據(jù)主成分分析結(jié)果，相對(duì)苗鮮質(zhì)量、相對(duì)苗長(zhǎng)和相對(duì)葉綠素含量3 個(gè)指標(biāo)與主成分相關(guān)且在主成分中相關(guān)系數(shù)較大，可以作為評(píng)價(jià)黍稷抗旱性的重要指標(biāo)。

表6 各主成分的因子載荷矩陣Tab.6 Factor load matrix of each principal component

參照寶力格等[17]的方法計(jì)算成分得分系數(shù)，并進(jìn)一步得到綜合得分計(jì)算公式：F=0.594 93F1+0.143 24F2。

根據(jù)公式計(jì)算得出所有黍稷資源苗期抗旱性綜合得分（F）進(jìn)行排序（表7），并使用沃德聯(lián)接對(duì)35 份黍稷材料進(jìn)行聚類分析（圖1）。

表7 黍稷種質(zhì)資源苗期抗旱性排序Tab.7 Ranking of drought resistance of broomcorn millet germplasm resources at seedling stages

圖1 35 份黍稷品種的苗期抗旱性聚類Fig.1 Cluster analysis of drought resistance of 35 broomcorn millet varieties at seedling stage

由圖1 可知，供試材料可以劃歸5 個(gè)類群。第1類群包括2 個(gè)品種（5.7%），F(xiàn) 值在3.495～3.647，基本不受PEG 脅迫影響，可以認(rèn)定為苗期極強(qiáng)抗旱材 料。第2 類 群 包 括3 個(gè) 品 種（8.6%），F(xiàn) 值 在1.832～2.217，PEG 干旱脅迫下受到的影響很小，可以認(rèn)定為苗期強(qiáng)抗旱材料。第3 類群一共包括8 個(gè)品種（22.9%），苗期綜合得分F 值在1.739～2.195，在PEG 干旱脅迫下受到的影響較小，可以認(rèn)定為苗期中度抗旱材料。第4 類群共包括11 個(gè)品種（31.4%），苗期綜合得分F 值為1.408～1.696，在PEG 干旱脅迫下受到的影響比較嚴(yán)重，可以認(rèn)定為苗期極不抗旱材料。第5 類群共包括11 個(gè)品種（31.4%），F(xiàn) 值在0.747～1.375，PEG 脅迫下幾乎失去抗旱能力，可以認(rèn)定為苗期極不抗旱材料。根據(jù)上述方法及分析結(jié)果，得到苗期強(qiáng)抗旱材料5 種，分別是千溝黃糜子（00005708）、大白黍（00000440）、Jan-55（00000281）、西北天糜子（黍）（00004230）、熟谷（00001987）。

3 結(jié)論與討論

作物種類決定干旱脅迫帶來(lái)的影響[18]，干旱可能對(duì)某一作物的某一指標(biāo)影響較明顯，但另一作物對(duì)對(duì)應(yīng)指標(biāo)卻影響甚微[19]。近年來(lái)，關(guān)于黍稷基因元件解析、種質(zhì)資源遺傳多樣性評(píng)價(jià)[20]、營(yíng)養(yǎng)價(jià)值評(píng)估[21]、生態(tài)栽培區(qū)劃[22]、生物和非生物脅迫等方面[23]已有許多研究。此前研究中不乏關(guān)于黍稷抗旱性鑒定，主要集中在萌芽和幼苗期，如劉天鵬等[24]采用反復(fù)干旱法對(duì)全國(guó)56 份黍稷品種進(jìn)行芽期抗旱性鑒定，結(jié)果發(fā)現(xiàn)，相對(duì)根干質(zhì)量、相對(duì)芽干質(zhì)量、相對(duì)發(fā)芽率、相對(duì)貯藏物質(zhì)轉(zhuǎn)運(yùn)率可用來(lái)評(píng)估芽期抗旱性；張盼盼[25]以小宗糧豆研究中心40 份黍稷資源為材料，采用PEG-6000 模擬干旱脅迫，篩選出芽期抗旱性指標(biāo)為相對(duì)根長(zhǎng)、相對(duì)芽長(zhǎng)、相對(duì)芽干質(zhì)量、相對(duì)發(fā)芽率和相對(duì)芽鞘長(zhǎng)，用以鑒定抗旱程度；邵歡歡等[26]以230 份黍稷資源為材料，采用相關(guān)性分析，篩選出相對(duì)葉綠素含量、相對(duì)苗鮮質(zhì)量和相對(duì)苗長(zhǎng)用以鑒定苗期抗旱性。本試驗(yàn)選取黍稷芽期和苗期共9 個(gè)指標(biāo)，覆蓋了能反映生長(zhǎng)狀況的指標(biāo)，采用相關(guān)性分析找到生育時(shí)期呈現(xiàn)顯著水平的抗旱指標(biāo)，其中，相對(duì)苗鮮質(zhì)量、相對(duì)苗長(zhǎng)和相對(duì)葉綠素含量在主成分中相關(guān)系數(shù)較大，可以作為評(píng)價(jià)黍稷抗旱性的重要指標(biāo)，這與前人研究結(jié)果基本相符。

PEG 引發(fā)是近20 a 來(lái)興起的一項(xiàng)種子處理技術(shù)，其主要功能是減緩萌發(fā)初期水分進(jìn)入種子的速率，減輕種子吸脹過(guò)程中膜系統(tǒng)的損傷，有利于受損膜系統(tǒng)的修復(fù)[27]。PEG 被廣泛應(yīng)用于種子引發(fā)和模擬干旱試驗(yàn)，本試驗(yàn)通過(guò)對(duì)PGE 干旱脅迫與正常供水情況下進(jìn)行相關(guān)性分析,得出全生育時(shí)期的9 個(gè)指標(biāo)間絕大多數(shù)呈顯著或極顯著水平，其中，相對(duì)苗鮮質(zhì)量與相對(duì)苗干質(zhì)量的相關(guān)性最高，芽期各指標(biāo)與苗期各種指標(biāo)間均相關(guān)性不顯著，相對(duì)葉綠素含量與苗期及芽期相關(guān)性均較低。本試驗(yàn)將黍稷抗旱性這一主觀的、模糊的評(píng)價(jià)方法，在數(shù)理統(tǒng)計(jì)的基礎(chǔ)上作了定量表達(dá)，這在黍稷的抗旱性研究中具有重要的意義[28]。此前已有研究表明，通過(guò)PEG 模擬干旱脅迫法可測(cè)定和比較不同種子的耐旱能力[29]。但也有結(jié)果表明，PEG 脅迫會(huì)延緩種子的萌發(fā)[30]，這可能是由于種子類型不同所致，有待進(jìn)一步的研究。通過(guò)計(jì)算得出所有黍稷資源苗期抗旱性綜合得分（F）并采用聚類分析對(duì)抗旱性進(jìn)行綜合分析，得出苗期強(qiáng)抗旱以上的品種有5 種，分別是千溝黃糜子（00005708）、大白黍（00000440）、Jan-55（00000281）、西北天糜子（黍）（00004230）、熟谷（00001987），可用于后續(xù)抗旱基因發(fā)掘和抗旱品種選育。