董 揚

(黑龍江省農(nóng)業(yè)科學(xué)院齊齊哈爾分院，黑龍江齊齊哈爾 161006)

糜子(PanicummiliaceumL.)是C4作物，能忍受鹽堿、干旱、高溫等極端條件[1]，是干旱、半干旱地區(qū)的重要特色雜糧作物[2]。我國具有豐富的糜子種質(zhì)資源，目前共收集10 200余份，國家已存儲 8 850 余份，尤其是地方品種資源較多[3]，具有廣泛的遺傳變異性和豐富的抗逆基因[4-5]。糜子的生育期短，是東北等冷涼地區(qū)的救災(zāi)抗逆作物，在北方有較大的生產(chǎn)優(yōu)勢和地區(qū)優(yōu)勢。

植物的耐冷性是指作物個體對低溫冷害所表現(xiàn)出的抵抗性或耐受性[6]，是植物對低溫逆境而表現(xiàn)出的一種生理應(yīng)激反應(yīng)[7]。農(nóng)業(yè)生產(chǎn)中，減輕低溫冷害所造成的損失最有效的手段就是選用耐冷品種，而農(nóng)作物耐冷性評價體系的建立是鑒定和篩選耐冷種質(zhì)資源的基礎(chǔ)與前提[8]。種子在萌發(fā)階段對低溫反應(yīng)最為敏感，低溫破壞了種子細(xì)胞的生物活性，從而大幅度地降低了種子的發(fā)芽勢、發(fā)芽率、發(fā)芽指數(shù)等[9]。近些年，在耐冷性鑒定分析評價方法上，前人已經(jīng)做了大量工作，種子萌發(fā)的時間和活力最先受到低溫影響[10]。環(huán)境溫度越低，種子萌發(fā)所需的時間越長[9]。長時間的萌發(fā)會使種子受到土壤中病原菌的侵害，使種子腐爛，降低發(fā)芽率[11-12]。在種子萌發(fā)期，發(fā)芽率、發(fā)芽勢、發(fā)芽指數(shù)、根長、芽長等指標(biāo)會發(fā)生不同程度的變化，是植物受到低溫脅迫最直接、具體的表現(xiàn)[13-14]。徐偉豪將冷脅迫14 ℃條件下水稻芽期的發(fā)芽率設(shè)置5個耐冷評價標(biāo)準(zhǔn)，發(fā)芽率>90%為極強，80%<發(fā)芽率≤90%為強，70%<發(fā)芽率≤80%為中，60%<發(fā)芽率≤70%為弱，發(fā)芽率≤60%為極弱[15]。Marocco等的研究表明，玉米的最適生長溫度為 25～28 ℃，在溫度低于10 ℃時植株生長會變得明顯緩慢，6～8 ℃ 時就不再生長發(fā)育[16]。王俊娟等的研究表明，棉花最低萌發(fā)溫度為10～12 ℃，最適溫度為 20～30 ℃[17-18]。辣椒的適宜發(fā)芽生長溫度在20～30 ℃ 之間，長時間在12 ℃以下生長會受到冷害[19]。唐雙勤等設(shè)計 5、8、10 ℃等3 個低溫處理直播秈稻，結(jié)果表明 8 ℃處理最適鑒別品種間耐冷性差異[20]。楊光等對296份玉米自交系進行了耐冷性鑒定研究，結(jié)果表明在4 ℃低溫下，yw9706自交系全部發(fā)芽，而自交系龍抗11完全沒有發(fā)芽[21]，說明同一作物不同品種間耐冷性存在差異。糜子種子發(fā)芽最適溫度為20～30 ℃[22]，當(dāng)溫度低于 12 ℃ 時容易對糜子造成冷害脅迫[23]。

目前，對糜子的研究主要集中在抗旱、抗倒伏和抗鹽堿等抗逆方面[24-26]，對耐冷能力的研究鮮見報道。本研究以240份糜子種質(zhì)資源為參試材料，通過逐級低溫篩選，確定糜子的最低發(fā)芽溫度，通過萌發(fā)期不同低溫處理，調(diào)查耐冷性相關(guān)指標(biāo)，并利用模糊數(shù)學(xué)隸屬函數(shù)、正態(tài)分布檢驗、相關(guān)性分析、主成分分析等統(tǒng)計方法進行綜合分析，以期建立糜子耐冷性綜合評價體系，為耐冷糜子品種資源的鑒定提供方法，并對資源加以利用。

1 材料與方法

1.1 供試材料

試驗于2021年在黑龍江省農(nóng)業(yè)科學(xué)院齊齊哈爾分院實驗室內(nèi)進行，該科研基地(123°45′E，47°15′N)位于黑龍江省西部干旱半干旱地區(qū)，海拔141 m，年平均氣溫為3.37 ℃，積溫為2 600～2 700 ℃，生長季總降水量為415 mm。試驗材料是從適宜在該科研基地種植的820份糜子品種資源中，根據(jù)2017—2019年測量的株高、落粒性、穗質(zhì)量、千粒質(zhì)量等農(nóng)藝性狀和產(chǎn)量性狀指標(biāo)初步篩選出的240份優(yōu)質(zhì)糜子資源，選種以來源廣泛、適宜當(dāng)?shù)胤N植、成熟度好、農(nóng)藝性狀和產(chǎn)量性狀突出為原則。

240份糜子資源來自21個國家，其中我國資源209份，國外資源31份，我國資源中包括新疆11份、山西27份、寧夏16份、陜西16份、青海10份、內(nèi)蒙古33份、黑龍江31份、遼寧8份、西藏1份、吉林8份、云南2份、山東7份、江蘇3份、甘肅17份、湖北1份、海南2份、北京3份、河北13份。

1.2 試驗設(shè)計

將糜子種子在1%的NaClO溶液中消毒15 min后，用自來水沖洗5 min，將種子置于陰涼通風(fēng)處陰干后置于鋪有3層濾紙的大號培養(yǎng)皿中，每個培養(yǎng)皿中均等擺放籽粒飽滿、大小顏色一致的種子30粒。將培養(yǎng)皿置于人工氣候箱內(nèi)，每個試驗均設(shè)置3次重復(fù)。設(shè)25 ℃(對照)及12、10、8、6、4 ℃ 6個溫度處理試驗，黑暗24 h，設(shè)置濕度為70%。

1.2.1 耐冷種子最低發(fā)芽溫度試驗方法 將240份糜子資源逐級放入25、12、10、8、6、4 ℃中篩選，25、12、10 ℃處理的發(fā)芽率以80%為標(biāo)準(zhǔn)，發(fā)芽率≥80%的資源進入下一水平低溫處理，8、6、4 ℃處理的發(fā)芽率以60%為標(biāo)準(zhǔn)，發(fā)芽率≥60%的資源進入下一水平低溫處理。

1.2.2 最適低溫篩選試驗方法 為確定糜子萌發(fā)期耐冷性鑒定的最適溫度，根據(jù)不同地域來源從供試糜子資源中隨機選出25份材料，分別于25、12、10、8、6 ℃(試驗“1.2.1”節(jié)確定出的最低發(fā)芽溫度)5個溫度條件下處理10 d，之后轉(zhuǎn)入25 ℃條件下恢復(fù)生長6 d。每個試驗均設(shè)置3次重復(fù)。低溫脅迫處理期間無光照，生長恢復(fù)階段07：00—19：00光照，每天12 h，設(shè)置2 500 lx的光照強度和70%的相對濕度。每天記錄發(fā)芽過程中的發(fā)芽粒數(shù)，計算發(fā)芽指數(shù)，生長恢復(fù)第3天測定發(fā)芽勢，第6天測發(fā)芽率，并計算活力指數(shù)。

1.2.3 萌發(fā)期耐冷性鑒定試驗 將240份糜子資源分別放在10 ℃(最適低溫篩選溫度)和25 ℃(對照)2個溫度條件下萌發(fā)，培養(yǎng)方法同“1.2.2”節(jié)。發(fā)芽過程中每天記錄發(fā)芽粒數(shù)，計算出發(fā)芽指數(shù)，恢復(fù)生長第3天測定發(fā)芽勢，第6天測定發(fā)芽率、根長、根鮮質(zhì)量、芽長、芽鮮質(zhì)量，計算活力指數(shù)和冷害率等。根長為從種胚到最長根尖的長度；芽長為從種胚到最長葉尖的長度。萌發(fā)期生長指標(biāo)計算公式如下：發(fā)芽勢=第3天發(fā)芽種子數(shù)/處理種子數(shù)×100%；發(fā)芽率=第6天發(fā)芽種子數(shù)/處理種子數(shù)×100%；發(fā)芽指數(shù)=∑(第n天的發(fā)芽種子數(shù)/相應(yīng)發(fā)芽天數(shù))；活力指數(shù)=發(fā)芽指數(shù)×胚根鮮質(zhì)量；發(fā)芽相對冷害率=(對照發(fā)芽率-處理發(fā)芽率)/對照發(fā)芽率×100%；相對值=各指標(biāo)處理值/對照值。

1.3 數(shù)據(jù)處理與分析

用Excel 2019軟件進行數(shù)據(jù)統(tǒng)計分析，用SPSS 19.0軟件進行相關(guān)性分析、主成分分析，用R studio軟件進行聚類分析。

采用標(biāo)準(zhǔn)差系數(shù)賦予權(quán)重法[27]進行耐冷性綜合評價：采用公式(1)計算隸屬函數(shù)值，用公式(2)計算標(biāo)準(zhǔn)差系數(shù)，用公式(3)計算各指標(biāo)的權(quán)重，用公式(4)計算各材料的綜合評價值。

(1)

式中：U(Xj)表示第j個指標(biāo)的隸屬函數(shù)值；Xj表示第j個綜合指標(biāo)；Xmin表示第j個綜合指標(biāo)的最小值；Xmax表示第j個綜合指標(biāo)的最大值。

(2)

(3)

(4)

式中：D表示各品種用綜合指標(biāo)評價所得的綜合評價值。

2 結(jié)果與分析

2.1 低溫萌發(fā)溫度評價

2.1.1 糜子萌發(fā)期最低發(fā)芽溫度的確定 對240份糜子資源進行逐級降溫篩選，結(jié)果(圖1)表明，在25 ℃萌發(fā)時出芽率均大于98%，說明參試糜子資源種子質(zhì)量良好，因此全部進行逐級低溫發(fā)芽試驗。在12 ℃溫度下發(fā)芽率大于80%的糜子資源有195份，占總數(shù)的81.25%，195份資源在10 ℃培養(yǎng)箱中發(fā)芽率大于80%的糜子資源有164份，占總數(shù)的68.33%，164份資源在8 ℃培養(yǎng)箱中發(fā)芽率大于60%的資源有84份，占總數(shù)的35.00%，將篩選出的84份糜子放入6 ℃培養(yǎng)箱中發(fā)芽，結(jié)果顯示發(fā)芽率大于60%的糜子資源有22份，占總數(shù)的9.16%，22份糜子資源在4 ℃培養(yǎng)箱中均沒有發(fā)芽。因此，糜子萌發(fā)期最低發(fā)芽溫度為6 ℃。

2.1.2 不同低溫處理下各指標(biāo)的耐冷響應(yīng) 由表1可知，在12、10、8、6 ℃低溫脅迫下，25份糜子資源的4個發(fā)芽性狀均表現(xiàn)出不同程度的下降趨勢，且各指標(biāo)的相對值大部分都小于1，平均值均小于1，說明4個低溫處理對25份糜子的萌發(fā)均造成一定程度的冷害，并且相對發(fā)芽勢的變異系數(shù)范圍為9.12%～113.14%，相對發(fā)芽率的變異系數(shù)范圍為7.10%～89.42%，相對發(fā)芽指數(shù)的變異系數(shù)范圍為8.01%～94.33%，相對活力指數(shù)的變異系數(shù)范圍為17.12%～105.00%，表明測定的指標(biāo)對4個低溫反應(yīng)敏感，通過對供試材料萌發(fā)期相對指標(biāo)的統(tǒng)計分析可知，25份抽選樣本的遺傳多樣性較高，能夠代表240份參試糜子材料，可以進一步分析以確定糜子耐冷性評價的最適溫度。

2.1.3 確定最適耐冷性評價溫度 確定篩選耐冷糜子資源的最適低溫，用隸屬函數(shù)法計算出25個糜子資源在4個低溫下的芽期指標(biāo)綜合得分(D)(表2)，對D進行正態(tài)分布檢驗和QQ-Plot分析(圖2)。

表1 不同低溫處理下25份糜子材料的萌發(fā)表現(xiàn)

從圖2-A可以看出，12、10、8 ℃正態(tài)分布檢驗的峰度值分別為-0.90、-0.72、-1.46，均小于0，呈平頂峰，6 ℃處理下的峰度值為1.46，大于0，呈尖頂峰，10 ℃的峰度值最接近于0，6 ℃正態(tài)檢驗偏度值為1.36，呈較大程度的正向偏離，12、10、8 ℃的偏度值分別為 -0.62、-0.12、-0.43，均小于0，呈負(fù)向偏離。結(jié)果表明10 ℃處理下的偏度絕對值僅為0.12，是所有處理中偏度最小的，正態(tài)分布擬合度最高。隸屬函數(shù)的QQ-Plot分析結(jié)果(圖2-B)表明，25份糜子材料在10 ℃處理下相對發(fā)芽勢、相對發(fā)芽率、相對發(fā)芽指數(shù)和相對活力指數(shù)的隸屬函數(shù)值與趨勢線擬合度最高(R2=0.985)，最符合正態(tài)分布。所以，10 ℃處理為篩選耐冷糜子資源的最適溫度。

2.2 萌發(fā)期耐冷性評價

2.2.1 萌發(fā)期糜子資源對冷脅迫的響應(yīng) 從表3可以看出，在10 ℃條件下240份糜子品種資源萌發(fā)期的各項指標(biāo)與對照比的相對值大部分都小于1，平均值均小于1，說明10 ℃低溫處理對參試材料的萌發(fā)造成了一定程度的冷害。相對發(fā)芽勢的范圍為0.00～1.00，相對發(fā)芽率的范圍為0.08～1.00，相對發(fā)芽指數(shù)的范圍為0.02～0.38，相對芽長的范圍為0.09～0.93，相對根長的范圍為0.21～2.98，相對芽鮮質(zhì)量的范圍為0.03～1.33，相對根鮮質(zhì)量的范圍為0.18～2.13，相對活力指數(shù)的范圍為0.02～0.69，相對發(fā)芽冷害率的范圍為0.00～0.92。其中，相對根鮮質(zhì)量、相對活力指數(shù)、相對發(fā)芽冷害率的差異性較大，變異系數(shù)分別是43.33%、58.45%、100.03%。表明不同品種資源的耐冷性存在差異，具有豐富的遺傳多樣性，參試材料的測定指標(biāo)對 10 ℃ 低溫脅迫反應(yīng)敏感，測試材料可作為在冷脅迫下鑒定耐冷性狀的品種資源。

2.2.2 萌發(fā)期冷脅迫下各性狀的相關(guān)性分析 從相關(guān)性分析結(jié)果(表4)可以看出，相對發(fā)芽勢、相對發(fā)芽率、相對發(fā)芽指數(shù)分別與相對芽鮮質(zhì)量、相對根鮮質(zhì)量、相對活力指數(shù)、相對冷害率等指標(biāo)間呈極顯著負(fù)相關(guān)(P<0.01)；相對根長與相對發(fā)芽率呈顯著正相關(guān)(P<0.05)，與相對芽鮮質(zhì)量、相對根鮮質(zhì)量、相對活力指數(shù)、相對冷害率呈顯著負(fù)相關(guān)(P<0.05)，與相對芽長、相對發(fā)芽指數(shù)、相對發(fā)芽勢呈極顯著正相關(guān)(P<0.01)；相對發(fā)芽勢、相對發(fā)芽率、相對芽長、相對發(fā)芽指數(shù)這4個指標(biāo)相互間均呈極顯著正相關(guān)(P<0.01)；相對根鮮質(zhì)量、相對芽鮮質(zhì)量、相對活力指數(shù)、相對冷害率這4個指標(biāo)相互間均呈極顯著正相關(guān)(P<0.01)。相對發(fā)芽勢、相對發(fā)芽率和相對發(fā)芽指數(shù)相互間的相關(guān)系數(shù)均大于0.80，相對芽鮮質(zhì)量、相對活力指數(shù)和相對根鮮質(zhì)量相互間的相關(guān)系數(shù)均大于0.90，說明各性狀間受到冷脅迫后的響應(yīng)是相互關(guān)聯(lián)的。

表3 萌發(fā)期10 ℃冷脅迫下各性狀相對值的統(tǒng)計分析結(jié)果

表4 萌發(fā)期冷脅迫下各性狀相對值的相關(guān)系數(shù)

2.2.3 萌發(fā)期冷脅迫下糜子測量性狀相對值的主成分分析 對240份糜子冷脅迫下9個相對性狀進行主成分分析，選取3個特征值大于1的主成分。由表5可知，第1主成分的貢獻(xiàn)率為52.57%，第2主成分的貢獻(xiàn)率為22.87%，第3主成分的貢獻(xiàn)率為11.78%，累計貢獻(xiàn)率為87.22%，表明3個主成分涵蓋了9個相對指標(biāo)87%以上的數(shù)據(jù)信息。由表6可知，9個性狀相對值的變化趨勢可以通過前3個主成分代表，它們之間的相關(guān)性可通過因子載荷量反映出來。 第1主成分與相對發(fā)芽勢、相對發(fā)芽率和相對發(fā)芽指數(shù)呈較大的正相關(guān)，表明冷脅迫下糜子發(fā)芽性狀可作為耐冷性評價的重要指標(biāo)。第2主成分與相對芽鮮質(zhì)量和相對根鮮質(zhì)量呈較大的正相關(guān)，表明冷脅迫與種子發(fā)芽后的鮮質(zhì)量有較大關(guān)系。第3主成分與相對根長呈較明顯的正相關(guān)，表明根長對其起主要作用。

表5 3個主成分的特征值以及貢獻(xiàn)率

表6 各因子載荷矩陣

3 討論

3.1 糜子耐冷溫度評價

近年來，各國學(xué)者對不同農(nóng)作物種子萌發(fā)期耐冷性方面做了大量研究。張鶴等采用模糊數(shù)學(xué)隸屬函數(shù)方法研究認(rèn)為，花生萌發(fā)期進行耐冷性評價最佳耐冷溫度為6 ℃[28]；孫世臣等研究發(fā)現(xiàn)，寒地水稻耐冷性評價最適溫度為10 ℃[29]。雖然前人對多種作物種子萌發(fā)期的耐冷性做了大量研究，但在處理糜子的耐冷低溫和評價指標(biāo)上還未見報道，因此，確定糜子耐冷性評價的適宜溫度條件是進行糜子耐冷性鑒定工作的前提。本研究以240份糜子種質(zhì)資源為參試材料，分別在 25、12、10、8、6、4 ℃的低溫條件下處理，逐級篩選后發(fā)現(xiàn)，部分糜子可在6 ℃發(fā)芽，而在4 ℃條件下未見發(fā)芽，因此在25、12、10、8、6 ℃的溫度下，對根據(jù)不同地域來源隨機選出25份材料進行處理，分別測定低溫條件下的相對發(fā)芽勢、相對發(fā)芽率、相對發(fā)芽指數(shù)和相對活力指數(shù)作為評價指標(biāo)，利用標(biāo)準(zhǔn)差系數(shù)賦予權(quán)重法、標(biāo)準(zhǔn)正態(tài)分布檢驗和QQ-Plot對數(shù)據(jù)進行分析，結(jié)果顯示在6、8 ℃低溫條件下萌發(fā)糜子對種子的活力影響較大，無法對耐冷能力進行準(zhǔn)確判斷；在12 ℃低溫條件下種子萌發(fā)對活力的影響較小， 品種間的差異并不顯著；10 ℃處理后的各項活力指標(biāo)均呈明顯的正態(tài)分布，且品種間差異顯著。因此，在低溫10 ℃時對糜子種子做耐冷發(fā)芽試驗可作為糜子種質(zhì)資源的耐冷性鑒定條件。糜子耐冷性篩選的最適溫度與寒地水稻[29]、玉米雜交種[30]等相同，但與花生[28]和高粱[31]等作物相比有所提高。說明不同作物的最適耐冷性評價溫度存在差異。

3.2 糜子萌發(fā)期耐冷性狀評價

目前還未見有糜子資源萌發(fā)期耐冷性狀的評價。袁雨豪等對100份糜子資源通過鹽脅迫測定萌發(fā)期的發(fā)芽率、發(fā)芽勢、發(fā)芽指數(shù)等指標(biāo)，通過主成分分析，確定糜子萌發(fā)期耐鹽性鑒定性狀指標(biāo)[26]。本研究以240個糜子品種資源為試驗材料，測定 10 ℃ 低溫脅迫下糜子萌發(fā)期的相對發(fā)芽率、相對發(fā)芽勢、相對冷害率等指標(biāo)(以25 ℃為對照)，通過相關(guān)性分析，發(fā)現(xiàn)耐冷性狀之間受到冷脅迫的程度是相互關(guān)聯(lián)的。主成分分析后將發(fā)芽率、發(fā)芽勢和發(fā)芽指數(shù)作為糜子萌發(fā)期耐冷性的鑒定指標(biāo)，與油菜[32]、水稻[33]等作物的耐冷性鑒定指標(biāo)有一定相似性。

4 結(jié)論

通過逐級降溫法發(fā)現(xiàn)部分糜子的最低發(fā)芽溫度為6 ℃，利用標(biāo)準(zhǔn)差系數(shù)賦予權(quán)重法計算綜合隸屬函數(shù)值，通過正態(tài)分布等確定最適耐冷性評價溫度，結(jié)合主成分分析建立糜子資源萌芽期耐冷綜合評價體系：在萌發(fā)期，10 ℃處理，以發(fā)芽率、發(fā)芽勢和發(fā)芽指數(shù)作為耐冷性鑒定指標(biāo)。