• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    燕麥種質(zhì)資源耐鹽堿性鑒定評價及耐鹽堿種質(zhì)篩選

    2023-04-12 01:38:22高文博張宗文周海濤
    作物學(xué)報 2023年6期
    關(guān)鍵詞:鹽堿發(fā)芽勢燕麥

    張 靜 高文博 晏 林 張宗文,2 周海濤 吳 斌,*

    燕麥種質(zhì)資源耐鹽堿性鑒定評價及耐鹽堿種質(zhì)篩選

    張 靜1高文博1晏 林1張宗文1,2周海濤3吳 斌1,*

    1中國農(nóng)業(yè)科學(xué)院作物科學(xué)研究所, 北京 100081;2國際生物多樣性中心, 北京 100081;3張家口市農(nóng)業(yè)科學(xué)院, 河北張家口 075000

    受氣候變化及人類生產(chǎn)活動影響, 世界鹽堿地范圍不斷擴(kuò)大, 土壤鹽堿化現(xiàn)已成為限制農(nóng)業(yè)生產(chǎn)發(fā)展的重要因素。燕麥?zhǔn)且环N耐鹽堿性較強(qiáng)的作物, 為了評價燕麥種質(zhì)資源耐鹽堿性及篩選耐鹽堿的燕麥種質(zhì), 本研究利用125 mmol L–1NaCl、Na2SO4和NaHCO3(物質(zhì)的量比為1∶1∶1)的混合鹽堿溶液對485份燕麥核心種質(zhì)材料進(jìn)行發(fā)芽脅迫處理, 測定了發(fā)芽期發(fā)芽勢、發(fā)芽率、根長、芽長、根鮮重、芽鮮重、根干重和芽干重8項指標(biāo), 利用相關(guān)性分析、主成分分析、隸屬函數(shù)分析及聚類分析等方法對參試燕麥種質(zhì)進(jìn)行了耐鹽堿綜合評價及篩選。結(jié)果表明: 鹽堿脅迫對測定的8項指標(biāo)均表現(xiàn)出抑制作用; 鹽堿脅迫下各指標(biāo)的隸屬函數(shù)值之間, 及各指標(biāo)隸屬函數(shù)值與綜合評價值之間均呈現(xiàn)出極顯著的正相關(guān)性; 通過主成分分析將8項測定指標(biāo)轉(zhuǎn)換成了2個主成分, 累計方差貢獻(xiàn)率為76.926%; 結(jié)合隸屬函數(shù)分析與聚類分析共篩選出2份高耐鹽堿型的燕麥材料: 燕1606和海闊夫斯基596號, 將485份燕麥種質(zhì)的耐鹽堿性劃分成了5個等級, 其中高耐鹽堿材料2份、耐鹽堿材料49份、中等耐鹽堿材料147份、鹽堿敏感材料129份、鹽堿極敏感材料158份; 綜合相關(guān)性分析、主成分分析和逐步回歸方程結(jié)果, 確定芽長作為發(fā)芽期燕麥耐鹽堿篩選的首選指標(biāo), 其次根鮮重、芽鮮重、發(fā)芽勢和發(fā)芽率也是進(jìn)行發(fā)芽期燕麥耐鹽堿綜合評價及篩選的重要指標(biāo)。

    燕麥; 種質(zhì)資源; 發(fā)芽期; 耐鹽堿; 綜合評價

    土壤鹽堿化是源于氣候條件、水文地質(zhì)、不合理灌溉及長期過量施用化肥等因素, 使得鹽分在土壤表層不斷的積累形成大量的可溶性鹽堿的現(xiàn)象[1]。由于氣候變化和人為因素的影響, 世界范圍內(nèi)鹽堿土面積平均每年以1.0×106~1.5×106hm2的速度增長[2-3]。根據(jù)聯(lián)合國教科文組織和糧農(nóng)組織的不完全統(tǒng)計, 全世界鹽堿地的面積約為9.54億公頃, 約占整個陸地面積的7%[4]。土壤鹽堿化現(xiàn)已成為全球所面臨的重大生態(tài)危機(jī), 是限制農(nóng)業(yè)可持續(xù)發(fā)展的重要因素[1,5]。中國鹽堿土面積大約有9913萬公頃, 占全國耕地總面積的10%, 主要分布在西北、東北西部、華北以及長江以北濱海地區(qū), 遍布全國的17個省區(qū), 是土壤鹽堿化危害最嚴(yán)重的國家之一[6]。鹽堿化使土壤理化性質(zhì)惡化, 土壤板結(jié)、肥力下降, 不利于農(nóng)作物吸收養(yǎng)分, 對農(nóng)作物生長發(fā)育產(chǎn)生危害, 阻礙作物生長[7]。當(dāng)土壤中的可溶性鹽堿含量超過0.2%時, 就會對作物產(chǎn)生不同程度的脅迫作用[8]。土地的有限性決定了人們需要對鹽堿土壤進(jìn)行利用與修復(fù), 而通過種植耐鹽堿作物對鹽堿地進(jìn)行生物改良, 可以在有效利用鹽堿地的同時, 吸收土壤中的礦質(zhì)元素, 增加空氣濕度, 減少地表蒸發(fā), 改善農(nóng)田小氣候, 抑制返鹽, 是一種改良鹽堿化土壤的重要方法[9-10]。

    燕麥(L.)是禾本科燕麥屬糧飼兼用型作物, 根據(jù)籽粒類型可分為皮燕麥和裸燕麥, 具有抗旱、抗寒、耐貧瘠等優(yōu)點, 在中國華北和西北農(nóng)牧交混地帶廣泛栽培種植[11], 而這些地區(qū)多有鹽堿地分布。目前, 國內(nèi)外在小麥[12-13]、玉米[14]、水稻[15]、大豆[16-17]等作物種質(zhì)資源耐鹽堿鑒定方面已有大量的研究報道。在燕麥耐鹽堿性鑒定研究方面, 前人分別對燕麥的耐鹽性和耐堿性開展了鑒定研究[18-19], 并且對少量燕麥品種進(jìn)行了耐鹽堿性鑒定[20-21], 這些研究為耐鹽堿性燕麥篩選奠定了良好的基礎(chǔ), 但與其它作物相比, 燕麥種質(zhì)資源的耐鹽堿性研究仍然存在一些不足之處: 一是相關(guān)燕麥耐鹽堿鑒定評價的研究較多側(cè)重于單一的耐鹽性或者耐堿性脅迫上, 而鹽堿地常存在鹽堿共同脅迫的問題, 針對實際生產(chǎn)中存在的混合鹽堿脅迫的相關(guān)研究較少; 二是缺乏燕麥種質(zhì)資源耐鹽堿性鑒定評價體系, 鑒定評價指標(biāo)較多且不盡相同, 不同研究結(jié)果間難以綜合比較; 三是相關(guān)研究選取的材料少, 代表性低, 缺乏針對我國現(xiàn)有燕麥種質(zhì)資源的系統(tǒng)性、大規(guī)模耐鹽堿性鑒定評價研究。因而開展燕麥種質(zhì)資源耐鹽堿性鑒定評價及耐鹽堿種質(zhì)篩選, 為培育耐鹽堿的燕麥品種提供優(yōu)異親本, 對于改良修復(fù)鹽堿土壤, 改善生態(tài)環(huán)境具有重要意義。本研究以國家作物種質(zhì)庫保存的燕麥種質(zhì)資源為基礎(chǔ), 選擇地理來源廣泛、遺傳多樣性豐富的485份燕麥種質(zhì)資源為材料, 以NaCl、Na2SO4和NaHCO3混合鹽堿溶液進(jìn)行耐鹽堿性鑒定, 測定了燕麥種子的發(fā)芽勢、發(fā)芽率、根長、芽長、根鮮重、根干重、芽鮮重和芽干重8項指標(biāo), 通過方差分析、相關(guān)性分析、主成分分析、聚類分析和隸屬函數(shù)分析對燕麥種質(zhì)資源耐鹽堿性進(jìn)行綜合評價鑒定, 以明晰中國現(xiàn)有燕麥種質(zhì)資源耐鹽堿性的現(xiàn)狀, 篩選出耐鹽堿的燕麥種質(zhì), 為相關(guān)育種研究提供優(yōu)異親本材料。

    1 材料與方法

    1.1 供試材料

    供試材料選擇地理來源廣泛、遺傳多樣性豐富的485份燕麥核心種質(zhì)資源, 均由中國農(nóng)業(yè)科學(xué)院國家作物種質(zhì)庫提供, 其中國內(nèi)種質(zhì)345份, 包括中國華北地區(qū)219份, 中國東北地區(qū)20份, 中國西南地區(qū)26份和中國西北地區(qū)80份; 國外種質(zhì)140份, 包括東歐地區(qū)43份, 西歐地區(qū)43份, 大洋洲地區(qū)11份, 東北亞地區(qū)9份, 美洲地區(qū)34份。

    1.2 試驗方法

    1.2.1 耐鹽堿鑒定適宜鹽堿濃度的篩選 隨機(jī)選取本試驗材料中的4個裸燕麥品種和4個皮燕麥品種進(jìn)行預(yù)試驗。選用NaCl、Na2SO4和NaHCO3(物質(zhì)的量比1∶1∶1)混合鹽堿溶液進(jìn)行鹽堿脅迫處理, 鹽堿濃度設(shè)置50、75、100、125、150、175和200 mmol L–17個水平, 用蒸餾水(CK)作為對照處理。選用籽粒飽滿、大小均勻一致的燕麥種子, 浸泡在1%的次氯酸鈉溶液中, 并放在每分鐘180轉(zhuǎn)的搖床上振蕩15 min, 再用流水沖洗3~5 min進(jìn)行種子消毒處理。消毒完成的種子浸泡在蒸餾水中3~4 h, 接種前將種子放置在紗布上吸干種子表面水分。發(fā)芽采用卷紙發(fā)芽法, 每個處理3次重復(fù), 釆用人工氣候箱對供試材料進(jìn)行培養(yǎng), 條件為光照/黑暗時間16 h/8 h, 溫度25℃/20℃。每天添加蒸餾水以彌補(bǔ)水分損失, 調(diào)查處理組和對照組的發(fā)芽率和發(fā)芽勢。

    1.2.2 耐鹽堿性鑒定 燕麥種質(zhì)資源耐鹽堿性鑒定評價參考《燕麥種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)》進(jìn)行, 略有修改[22]。以試驗1.2.1篩出的最佳鹽堿濃度125 mmol L–1為耐鹽堿鑒定處理濃度, 蒸餾水(CK)為對照, 對485份燕麥核心種質(zhì)進(jìn)行發(fā)芽期耐鹽堿性鑒定, 重復(fù)3次, 每個重復(fù)選取50粒籽粒飽滿、大小均勻一致且完整的燕麥種子, 種子消毒方式、發(fā)芽方法及培養(yǎng)條件均與預(yù)備試驗相同。

    1.2.3 高耐鹽堿型燕麥品種的極端耐鹽堿濃度確定 對試驗1.2.2篩選出的高耐鹽堿型材料進(jìn)一步進(jìn)行極端耐鹽堿濃度確定試驗, 共設(shè)置5個處理, 125、150、175和200mmol L–1的NaCl、Na2SO4和NaHCO3(物質(zhì)的量比1∶1∶1)混合鹽堿溶液為鹽堿脅迫處理, 蒸餾水(CK)作為對照處理, 每個處理3次重復(fù), 每個重復(fù)選取50粒燕麥種子, 消毒方式、發(fā)芽方法及培養(yǎng)條件均與預(yù)備試驗相同。

    1.3 指標(biāo)測定

    通過測定燕麥種子發(fā)芽期的發(fā)芽勢(GP)、發(fā)芽率(GR)、根長(RL)、芽長(BL)、根鮮重(RFW)、根干重(RDW)、芽鮮重(BFW)和芽干重(BDW) 8項指標(biāo), 綜合評價燕麥耐鹽堿性。

    以根長至少與種子等長, 芽長至少達(dá)到種子1/2的長度為發(fā)芽標(biāo)準(zhǔn)。培養(yǎng)第4天統(tǒng)計發(fā)芽勢, 第7天統(tǒng)計發(fā)芽率, 發(fā)芽勢和發(fā)芽率參照陳新等[19]的方法按照公式(1)和(2)計算。培養(yǎng)的第7天從每份材料的每個重復(fù)中隨機(jī)選取具有代表性的10株燕麥幼苗, 用直尺測量最長根長和芽長[19], 將選取的幼苗吸干表面水分后分成地上和地下兩部分稱量根鮮重和芽鮮重, 再烘干至恒重稱量根干重和芽干重。

    1.4 統(tǒng)計分析

    利用Microsoft Excel進(jìn)行數(shù)據(jù)的基礎(chǔ)整理, SPSS26.0、R4.1.3和Origin 2022軟件對數(shù)據(jù)進(jìn)行相關(guān)性分析、主成分分析、隸屬函數(shù)分析以及聚類分析, 具體公式如下[23-24]:

    采用煙梗自動開麻包系統(tǒng)后,由于梗開包過程是在一個密閉的負(fù)壓箱體內(nèi)完成,煙梗里的麻毛和粉塵在倒料過程中被充分揚(yáng)起,然后被除塵風(fēng)抽走,改善了環(huán)境。此外,新型梗投料對麻袋不產(chǎn)生破壞,麻袋可以重復(fù)利用,節(jié)約生產(chǎn)成本。

    運(yùn)用公式(3)以每一個測定指標(biāo)3次重復(fù)的平均值為基礎(chǔ)數(shù)據(jù)。計算每個測定指標(biāo)的耐鹽堿系數(shù)。

    采用隸屬函數(shù)分析對485份燕麥材料的耐鹽堿系數(shù)按照公式(4)計算隸屬函數(shù)值, 式中X表示第個指標(biāo)的耐鹽堿系數(shù),max、min表示第個指標(biāo)中耐鹽堿系數(shù)的最大值和最小值,(X)表示第個指標(biāo)的隸屬函數(shù)值。

    對各指標(biāo)的隸屬函數(shù)值進(jìn)行主成分分析, 確定主成分的個數(shù), 按照供公式(5)計算權(quán)重, 式中W表示第個指標(biāo)的權(quán)重,P表示第個指標(biāo)的方差貢獻(xiàn)率。

    運(yùn)用公式(6)計算各材料的綜合評價值, 進(jìn)行燕麥種質(zhì)資源耐鹽堿綜合評價, 綜合評價值越大說明該材料的耐鹽堿性越強(qiáng)。

    2 結(jié)果與分析

    2.1 適宜耐鹽堿鑒定的鹽堿脅迫濃度確定

    為了確定適宜耐鹽堿鑒定的鹽堿脅迫濃度, 從485份燕麥核心種質(zhì)中隨機(jī)選出4份皮燕麥品種和4份裸燕麥品種進(jìn)行不同鹽堿濃度的試驗處理。試驗結(jié)果如圖1所示, 隨著鹽堿濃度的不斷增加, 4份皮燕麥品種和4份裸燕麥品種的發(fā)芽勢和發(fā)芽率均呈現(xiàn)出不同程度的下降, 并且不同鹽堿濃度下燕麥品種間的發(fā)芽勢和發(fā)芽率存在差異。4份皮燕麥品種和4份裸燕麥品種在低濃度的鹽堿濃度處理下與對照的發(fā)芽勢和發(fā)芽率差異較小。當(dāng)鹽堿脅迫濃度大于100 mmol L–1時, 4份皮燕麥品種和4份裸燕麥品種間的發(fā)芽勢和發(fā)芽率與對照處理相比均存在顯著或極顯著差異。當(dāng)鹽堿脅迫濃度大于150 mmol L–1時, 4份皮燕麥品種和4份裸燕麥品種的發(fā)芽勢和發(fā)芽率均大幅度下降, 50%燕麥品種的發(fā)芽勢和發(fā)芽率降為0, 鹽堿濃度大于150 mmol L–1時燕麥不發(fā)芽,此時濃度已不適合進(jìn)行燕麥種質(zhì)的耐鹽堿篩選, 而當(dāng)鹽堿濃度為125 mmol L–1時, 4份皮燕麥品種和4份裸燕麥品種的發(fā)芽勢和發(fā)芽率與對照相比都表現(xiàn)出顯著性或極顯著性的差異, 在該鹽堿濃度下無論是皮燕麥還是裸燕麥均能夠有效的區(qū)分出各燕麥材料之間的耐鹽堿性, 因此確定125 mmol L–1的鹽堿濃度為適宜的燕麥種質(zhì)耐鹽堿鑒定濃度。

    圖1 不同鹽堿濃度下8份燕麥資源發(fā)芽指標(biāo)變化

    A: 裸燕麥發(fā)芽勢變化趨勢; B: 皮燕麥發(fā)芽勢變化趨勢; C: 裸燕麥發(fā)芽率變化趨勢; D: 皮燕麥發(fā)芽率變化趨勢; 同一品種不同字母表示處理間差異顯著(< 0.05)。

    A: the trend of germination potential of Naked oat; B: the trend of germination potential of hulled oat; C: the trend of germination rate of naked oat; D: the trend of germination rate of hulled oat. Different letters of the same variety indicate significant difference between treatments at< 0.05.

    2.2 鹽堿脅迫對燕麥發(fā)芽性狀的影響

    在125 mmol L–1混合鹽堿脅迫下, 485份燕麥材料發(fā)芽期各性狀變化如圖2和表1所示。8項測定指標(biāo)發(fā)芽勢、發(fā)芽率、根鮮重、芽鮮重、根干重、芽干重、根長和芽長的平均值在鹽堿脅迫下與對照相比都呈現(xiàn)出下降趨勢, 說明鹽堿脅迫對燕麥種子萌發(fā)和幼苗的生長均有不同程度的抑制作用。由表2可知鹽堿脅迫對根長的抑制作用最明顯, 均值為1.217 cm, 與對照條件下相比下降了88.85%, 說明根長對鹽堿脅迫反應(yīng)最敏感; 鹽堿脅迫對種子的發(fā)芽勢和發(fā)芽率的抑制作用最小, 均值分別為23.615%和26.245%, 與對照相比分別下降了49.83%和46.85%。

    由表1中的變異系數(shù)可知, 485份燕麥材料的8項測定指標(biāo)在對照和鹽堿脅迫條件下均表現(xiàn)出不同程度的變異。鹽堿脅迫處理下各指標(biāo)間的變異系數(shù)均大于對照處理下各個測定指標(biāo)的變異系數(shù)。對照處理下變異系數(shù)的變化范圍是12.88%~60.26%, 鹽堿脅迫處理下變異系數(shù)變化范圍是35.19%~71.68%, 說明本研究所選用的檢測指標(biāo)顯著受到鹽堿脅迫的影響, 因此這些指標(biāo)可以用于耐鹽堿性鑒定評價。而這些指標(biāo)的變異系數(shù)變化范圍較大, 說明鹽堿脅迫對各指標(biāo)的影響程度不同, 因此不能用單項指標(biāo)來評價燕麥的耐鹽堿性。

    圖2 鹽堿脅迫對燕麥發(fā)芽期性狀的影響

    CK: 對照處理; T: 125 mmol L–1鹽堿脅迫處理; GP: 發(fā)芽勢; GR: 發(fā)芽率; RL: 根長; BL: 芽長; RFW: 根鮮重; BFW: 芽鮮重; RDW: 根干重; BDW: 芽干重; ***: 表示在0.001概率水平顯著相關(guān)。

    CK: control treatment; T: 125 mmol L–1saline-alkali stress treatment; GP: germination potential; GR: germination rate; RL: root length; BL: bud length; RFW: root fresh weight; BFW: bud fresh weight; RDW: root dry weight; BDW: bud dry weight; *** indicates significant correlation at the 0.001 probability level.

    表1 485份燕麥種質(zhì)在對照和鹽堿脅迫條件下各性狀描述統(tǒng)計

    縮寫同圖2。Abbreviations are the same as those given in Fig. 2.

    2.3 燕麥耐鹽堿性綜合評價

    2.3.1 燕麥發(fā)芽期各指標(biāo)隸屬函數(shù)值及綜合評價值的相關(guān)性分析 如圖3所示, 相關(guān)性分析結(jié)果表明, 燕麥發(fā)芽期發(fā)芽勢、發(fā)芽率、根長、芽長、根鮮重、芽鮮重、根干重和芽干重之間存在著極顯著正相關(guān), 并且與值之間也存著極顯著正相關(guān)。其中芽長的隸屬函數(shù)值與值之間的相關(guān)系數(shù)最大為0.81, 發(fā)芽勢、發(fā)芽率、根鮮重的隸屬函數(shù)值與值之間也有較大的相關(guān)系數(shù)(均為0.80), 表明發(fā)芽勢、發(fā)芽率、根鮮重和根長與值之間都存著不少于80%的重疊信息??梢哉f這4個指標(biāo)與燕麥的耐鹽堿性關(guān)系最為密切。相關(guān)性分析結(jié)果顯示, 燕麥發(fā)芽期各性狀之間聯(lián)系密切且存在著信息的重疊, 因此要對其進(jìn)行降維分析, 提取主要的影響因子。

    圖3 燕麥各性狀隸屬函數(shù)值及D值的相關(guān)性分析

    縮寫同圖2;-value: 綜合評價值; **表示在0.01水平上顯著相關(guān)。

    Abbreviations are the same as those given in Fig. 2;-value: comprehensive evaluation value; ** indicates significant correlation at the 0.01 probability level.

    由相關(guān)性分析結(jié)果表明各性狀之間存在著極顯著正相關(guān), 信息重疊性較高, 因此對8項測定指標(biāo)進(jìn)行了主成分分析。通過對燕麥發(fā)芽期測定的8項指標(biāo)的隸屬函數(shù)值進(jìn)行主成分分析適應(yīng)性檢驗和Bartlett球形檢驗, 得到的KMO值為0.778, 大于0.7,=0.000, 小于0.01, 說明各性狀之間的關(guān)聯(lián)程度較高, 主成分分析結(jié)果可靠。對燕麥性狀的隸屬函數(shù)值進(jìn)行主成分分析, 根據(jù)特征值大于1以及累積貢獻(xiàn)率大于75%的原則, 提取主成分。主成分分析結(jié)果如表2所示, 根據(jù)主成分提取原則, 選取了前2個主成分, 第一主成分的特征值為4.629, 方差貢獻(xiàn)率為57.867%; 第二主成分的特征值為1.525, 方差貢獻(xiàn)率為19.060%。2個主成分的累計方差貢獻(xiàn)率達(dá)到了76.926%, 可有效反映出8項測定指標(biāo)絕大部分信息。因此, 將測定的8項指標(biāo)轉(zhuǎn)換成了2個綜合指標(biāo), 用于下一步燕麥耐鹽堿性綜合評價值計算。

    主成分載荷矩陣能反映出各個評價指標(biāo)對主成分的影響程度。由表3可知, 根鮮重、芽鮮重和芽長是第一主成分的主要指標(biāo), 相關(guān)系數(shù)分別為0.847、0.851和0.897, 該主成分與幼苗的生長密切相關(guān)。第二主成分與發(fā)芽勢和發(fā)芽率的相關(guān)性較大, 相關(guān)系數(shù)分別為0.745、0.747, 說明該主成分與燕麥的發(fā)芽指標(biāo)較為密切, 反映了燕麥種子的發(fā)芽情況。

    2.3.3 燕麥耐鹽堿性聚類分析 根據(jù)計算出的值, 對485份燕麥種質(zhì)進(jìn)行耐鹽堿排序, 利用系統(tǒng)聚類方法, 對其進(jìn)行聚類分析, 結(jié)果如圖4所示。485份燕麥種質(zhì)可劃分為5大類。第1類是高耐鹽堿型材料有2份, 占供試材料總數(shù)的0.41%, 第2類是耐鹽堿型材料有49份, 占供試材料總數(shù)的10.10%;第3類是中等耐鹽堿型材料有147份, 占供試材料總數(shù)的30.31%; 第4類是鹽堿敏感型材料有129份, 占供試材料總數(shù)的26.60%, 第5類是鹽堿極敏感型材料有158份, 占供試材料總數(shù)的32.58%。高耐鹽堿型的2份燕麥材料編號為411 (海闊夫斯基596號)和413 (燕1606)。

    表2 燕麥發(fā)芽期各指標(biāo)的主成分方差貢獻(xiàn)率

    表3 主成分載荷矩陣

    縮寫同圖2。Abbreviations are the same as those given in Fig. 2.

    圖4 485份燕麥種質(zhì)資源耐鹽堿性聚類分析

    I: 高耐鹽堿型; II: 耐鹽堿型; III: 中等耐鹽堿型; IV: 鹽堿敏感型; V: 鹽堿極敏感型。

    I: the highly salt-alkaline-tolerant accessions; II: the salt-alkaline-tolerant accessions; III: the medium salt-alkaline-tolerant accessions; IV: the salt-alkaline-sensitive accessions; V: the high salt-alkaline-sensitive accessions.

    2.3.4 高耐鹽堿型材料鹽堿極限濃度鑒定 通過對485份燕麥核心種質(zhì)資源的耐鹽堿性鑒定評價, 篩選出了2份高耐鹽堿型的燕麥材料分別為411 (海闊夫斯596號)和413 (燕1606)。如表4所示, 411 (海闊夫斯596號)發(fā)芽勢和發(fā)芽率在鹽堿脅迫為200 mmol L–1時降到了10%以下, 與對照相比表現(xiàn)出顯著性差異。413 (燕1606)的發(fā)芽勢在鹽堿脅迫為175 mmol L–1時降低到了10%以下, 發(fā)芽率在200 mmol L–1時降到了10%以下, 與對照相比呈現(xiàn)出顯著性差異, 說明這2個品種在鹽堿脅迫濃度為175 mmol L–1時種子發(fā)芽才受到嚴(yán)重抑制。

    表4 不同鹽堿濃度對燕麥發(fā)芽指標(biāo)的影響

    縮寫同圖2。Abbreviations are the same as those given in Fig. 2

    2.3.5 燕麥發(fā)芽期耐鹽堿回歸模型建立及評價指標(biāo)篩選 為了便于燕麥種質(zhì)資源的鑒定與篩選, 建立了燕麥耐鹽堿多元線性回歸數(shù)學(xué)模型。將485份燕麥種質(zhì)按照2∶1的比例隨機(jī)分成了回歸組和檢驗組2組, 回歸組323份材料8個性狀的耐鹽堿系數(shù)作為自變量, 耐鹽堿性D值作為因變量, 進(jìn)行多元逐步回歸分析。結(jié)合8個性狀指標(biāo)的偏回歸系數(shù), 最終建立了多元線性回歸模型:=0.245X+ 0.0522+0.1713+0.2954–0.01, 方程中代表耐鹽堿性值、1代表芽長、2代表發(fā)芽勢、3代表根鮮重、4代表根長。方程決定系數(shù)2=0.99,= 7741.329,<0.01, 該回歸方程對觀測數(shù)據(jù)的真實反映程度為99%。運(yùn)用該回歸方程計算檢驗組162份材料耐鹽堿的預(yù)測值, 并與對應(yīng)的值進(jìn)行相關(guān)性分析, 結(jié)果顯示相關(guān)性系數(shù)為0.996, 兩者之間呈極顯著正相關(guān), 表明該回歸方程能夠很好的預(yù)測燕麥種質(zhì)資源的耐鹽堿性。在進(jìn)行燕麥種質(zhì)資源耐鹽堿性鑒定評價時, 試驗條件相同的情況下, 只需要測定燕麥發(fā)芽期的芽長、發(fā)芽勢、根鮮重和根干重這4個指標(biāo), 利用該回歸方程計算耐鹽堿值, 從而能便捷、高效地評估不同燕麥種質(zhì)資源的耐鹽堿性強(qiáng)弱。

    2.4 不同地理來源燕麥種質(zhì)耐鹽堿性評價

    根據(jù)地理來源的不同將485份燕麥種質(zhì)資源分為9組, 對各組種質(zhì)的耐鹽堿性D值進(jìn)行統(tǒng)計分析。結(jié)果如表5所示, 不同地理來源的燕麥種質(zhì)資源之間耐鹽堿性存在差異, 來自于東北亞地區(qū)的燕麥種質(zhì)資源耐鹽堿性最強(qiáng), 與東歐地區(qū)、美洲地區(qū)和中國西北、西南和華北地區(qū)資源均存在著顯著差異; 來自于中國西南地區(qū)的燕麥種質(zhì)資源的耐鹽堿性最弱, 與來源于其他地區(qū)的燕麥種質(zhì)資源均呈顯著差異。中國東北和西北地區(qū)、歐洲地區(qū)、大洋洲地區(qū)和美洲地區(qū)的燕麥種質(zhì)資源的耐鹽堿性之間差異不顯著, 這幾個地區(qū)的燕麥種質(zhì)資源耐鹽堿性較為一致且耐鹽堿程度為中等水平。

    表5 不同地區(qū)燕麥D值描述統(tǒng)計

    3 討論

    土壤鹽堿化是世界性的生態(tài)環(huán)境問題, 鹽堿地??煞譃辂}地、堿地或含多種鹽堿成分的混合鹽堿地,中國大面積的鹽堿地主要是含有Na+、Cl–、CO32–、HCO3–、SO42–等多種離子的復(fù)合鹽堿地, 組分復(fù)雜, 并且鹽化與堿化往往相伴發(fā)生[25]。研究表明, 混合鹽堿脅迫對作物的危害并不是鹽、堿脅迫危害的簡單疊加, 其危害程度要遠(yuǎn)大于單一的鹽或堿脅迫, 通過單鹽篩選的耐鹽性作物并不都適合在鹽堿土上生長[26], 因此只采用單一的中性鹽或堿性鹽進(jìn)行脅迫試驗, 得到的結(jié)果對于燕麥耐鹽堿機(jī)制分析以及育種改良均具有一定的局限性。本研究以NaCl、Na2SO4和NaHCO3混合溶液來模擬鹽堿脅迫環(huán)境, 進(jìn)行燕麥耐鹽堿性鑒定, 與單一的中性鹽或堿性鹽相比, 在此條件下篩選出的品種更便于生產(chǎn)利用。

    鹽堿脅迫對作物不同生長發(fā)育時期的影響不同, Hosseini等[27]、徐玲秀等[28]研究發(fā)現(xiàn)大豆在芽期、苗期及花莢期對鹽堿的耐受性并不相同, 而且不同生育階段的耐鹽堿性并沒有直接的相關(guān)性, 水稻[29]中也有類似研究結(jié)果, 這說明耐鹽堿作物的篩選有賴于不同生育期的系統(tǒng)鑒定。種子發(fā)芽期是植株生長發(fā)育的初始階段, 發(fā)芽階段的耐鹽堿性對其能否在鹽堿地中存活起到至關(guān)重要的作用, 而且處于發(fā)芽期的作物對鹽堿脅迫更為敏感, 影響著作物后期的生長發(fā)育[30-31]。在實際生產(chǎn)中, 作物發(fā)芽期的耐鹽堿性直接決定能否在鹽堿脅迫條件下保證較高水平的出苗率, 對作物耐鹽堿性評價具有重要參考價值, 同時發(fā)芽期鑒定具有耗時短、容量大、操控性強(qiáng)和重復(fù)性好等特點, 方便不同研究結(jié)果進(jìn)行比較[32], 因此本研究選擇在發(fā)芽期對燕麥進(jìn)行耐鹽堿性鑒定評價。

    作物的耐鹽堿性是一個受多個數(shù)量性狀控制的復(fù)雜性狀, 單一的指標(biāo)并不能夠準(zhǔn)確地評價作物的耐鹽堿性[33-34]。但是對多個指標(biāo)進(jìn)行測定, 不僅需要耗費(fèi)大量的時間、人力、物力和財力, 而且測定的各個指標(biāo)之間還存在著相關(guān)性, 代表信息會出現(xiàn)重疊交叉現(xiàn)象, 也不能很好的評價作物耐鹽堿性[35]。因此, 要選取適宜的篩選指標(biāo)才能避免單一指標(biāo)或是多個指標(biāo)鑒定的局限性。本研究選取了8個性狀, 涵蓋了發(fā)芽期耐鹽堿篩選的所有指標(biāo), 綜合相關(guān)性分析、主成分分析和耐鹽堿多元線性回歸方程的結(jié)果進(jìn)行了最適宜鑒定指標(biāo)的篩選。相關(guān)性分析顯示根長與值的相關(guān)系數(shù)最大; 其次發(fā)芽勢、發(fā)芽率、根鮮重和芽長。因此, 從相關(guān)性分析結(jié)果, 可將芽長、發(fā)芽勢、發(fā)芽率和根鮮重作為評價燕麥發(fā)芽期耐鹽堿性的指標(biāo)。主成分分析顯示第一主成分的方差貢獻(xiàn)率為57.867%, 其中各項性狀指標(biāo)中芽長的載荷數(shù)最大; 其次是芽鮮重和根鮮重。第二主成分的方差貢獻(xiàn)率為19.060%, 其中各性狀指標(biāo)中發(fā)芽勢和發(fā)芽率的載荷數(shù)最大。因此, 主成分分析結(jié)果表明芽長、根鮮重、芽鮮重、發(fā)芽勢和發(fā)芽率可以作為評價燕麥發(fā)芽期耐鹽堿性的指標(biāo)。綜合相關(guān)性分析和主成分分析結(jié)果, 本研究表明芽長可作為燕麥發(fā)芽期耐鹽堿篩選的首選指標(biāo), 其次是根鮮重、發(fā)芽勢和發(fā)芽率也是進(jìn)行燕麥發(fā)芽期耐鹽堿綜合評價及篩選的重要指標(biāo)。本研究結(jié)果與前人篩選出的耐鹽堿指標(biāo)具有相似性。陳新等[19]認(rèn)為與裸燕麥萌發(fā)期耐鹽性關(guān)系最密切的指標(biāo)是發(fā)芽勢和發(fā)芽率。付鸞鴻等[21]認(rèn)為要鑒定燕麥萌發(fā)期耐鹽堿性的主要指標(biāo)是發(fā)芽勢、胚根長和胚根鮮重。斯日古楞等[36]認(rèn)為與燕麥萌發(fā)期耐鹽性鑒定最理想的指標(biāo)是胚芽長和種子萌發(fā)指數(shù)。同時本研究構(gòu)建了燕麥耐鹽堿性鑒定評價的多元線性回歸方程, 篩選出芽長、發(fā)芽勢、根鮮重和根干重4個重要指標(biāo)可快速評價燕麥種質(zhì)發(fā)芽期耐鹽堿性。在相同逆境條件下, 也可測定這4項指標(biāo), 利用該回歸方程快速預(yù)測燕麥種質(zhì)的耐鹽堿性, 為燕麥種質(zhì)資源耐鹽堿性鑒定及育種研究提供重要的信息。

    評價作物耐鹽堿性的指標(biāo)多種多樣, 且各項指標(biāo)之間存在著一定的線性關(guān)系, 使得代表的信息之間出現(xiàn)交叉或重疊。本研究結(jié)果也顯示了多個指標(biāo)的耐鹽堿系數(shù)間存在極顯著的相關(guān)性, 且不同單項指標(biāo)的耐鹽堿系數(shù)在不同種質(zhì)中的大小排序也有所不同。因此有必要利用綜合評價的方法來鑒定耐鹽堿這一個復(fù)雜的綜合性狀, 才能較全面、準(zhǔn)確及有效地反映出不同種質(zhì)資源之間的差異性?,F(xiàn)階段運(yùn)用綜合評價的方法鑒定作物的抗逆性已廣泛應(yīng)用與小麥[37]、玉米[38]和大豆[39]等作物中。在小麥抗逆性鑒定方面, 張婷婷等[40]通過主成分分析、隸屬函數(shù)分析和聚類分析綜合評價了30份春小麥萌發(fā)期和苗期的耐鹽堿性, 并成功篩選出2份高度耐鹽堿的小麥品種。段雅娟等[41]鑒定了玉米萌發(fā)期的10項生理指標(biāo), 運(yùn)用主成分分析、隸屬函數(shù)值分析和聚類分析綜合評價了390份玉米自交系的耐鹽性, 將其劃分為了6個耐鹽等級。在數(shù)據(jù)處理和耐鹽堿性評價過程中, 運(yùn)用綜合評價的方法既減少了指標(biāo)間代表信息的交叉重疊, 又考慮了綜合指標(biāo)重要性的差異, 避免了單一指標(biāo)的片面性, 較為科學(xué)和系統(tǒng)的反映了種質(zhì)資源的耐鹽堿性差異[42]。本研究基于燕麥發(fā)芽期不同性狀指標(biāo), 采用相關(guān)性分析、主成分分析、隸屬函數(shù)分析和聚類分析對485份燕麥種質(zhì)資源的耐鹽堿能力進(jìn)行了綜合鑒定評價, 篩選出了發(fā)芽期高耐鹽堿型的2份燕麥材料分別是海闊夫斯596號基和燕1606, 2份材料的值分別為0.712和0.858, 且當(dāng)鹽堿濃度提高到175 mmol L–1時兩者的發(fā)芽率才降低到10%以下。但還需在大田中進(jìn)一步研究這2份材料苗期和成株期的耐鹽堿性。這些耐鹽堿性材料的鑒定, 為耐鹽堿品種的選育和耐鹽堿基因的發(fā)掘提供了重要資源。

    4 結(jié)論

    在混合鹽堿脅迫下燕麥種子的發(fā)芽指標(biāo)和生長指標(biāo)都受到了抑制。從485份燕麥種質(zhì)資源中篩選出了海闊夫斯基596號和燕1606兩份高耐鹽型的種質(zhì)資源, 為耐鹽堿品種的選育和耐鹽堿基因的挖掘提供了資源。篩選出了芽長可作為燕麥發(fā)芽期耐鹽堿篩選的首選指標(biāo), 其次是根鮮重、芽鮮重、發(fā)芽勢和發(fā)芽率也是進(jìn)行燕麥發(fā)芽期耐鹽堿鑒定評價及篩選的重要指標(biāo), 在同等試驗條件下可作為燕麥發(fā)芽期耐鹽堿性鑒定評價及篩選的指標(biāo)。

    [1] Allbed A, Kumar L. Soil salinity mapping and monitoring in arid and semi-arid regions using remote sensing technology: a review., 2013, 2: 373–385.

    [2] Wang S, Sun L, Ling N, Zhu C, Chi F Q, Li W Q, Hao X Y, Bian J Y, Chen L, Wei D. Exploring soil factors determining composition and structure of the bacterial communities in saline-alkali soils of Songnen Plain., 2020, 10, 2902–2902.

    [3] Kovda Victor A. Loss of productive land due to salinization., 1983: 91–93.

    [4] 李彬, 王志春, 孫志高, 陳淵, 楊福. 中國鹽堿地資源與可持續(xù)利用研究. 干旱地區(qū)農(nóng)業(yè)研究, 2005, (2): 154–158.Li B, Wang Z C, Sun Z G, Chen Y, Yang F. Resources and sustainable resource exploitation of salinized land in China., 2005, (2): 154–158 (in Chinese with Englis abstract).

    [5] Amundson R, Berhe A A, Hopmans J W, Olson C, Sztein A E, Sparks D L. Soil and human security in the 21st century., 2015, 348: 1261071.

    [6] 王遵親. 中國鹽堿土. 北京: 科學(xué)出版社, 1993. Wang Z Q. Saline-alkaline Soils of China. Beijing: Science Press, 1993 (in Chinese).

    [7] Zhang W, Chong W, Rui X U E, Wang L J. Effects of salinity on the soil microbial community and soil fertility., 2019, 18: 1360 –1368.

    [8] Doran J C. Australian trees and shrubs: species for land rehabilitation and farm planting in the tropics., 1997, 24: 314–319.

    [9] 馬晨, 馬履一, 劉太祥, 左海軍, 張博, 劉寅. 鹽堿地改良利用技術(shù)研究進(jìn)展. 世界林業(yè)研究, 2010, 23(2): 28–32. Ma C, Ma L Y, Liu T X, Zuo H J, Zhang B, Liu Y. Research progress on saline land improvement technology., 2010, 23(2): 28–32 (in Chinese with English abstract).

    [10] 丁海榮, 洪立洲, 楊智青, 王茂文, 朱小梅, 劉沖. 鹽堿地及其生物措施改良研究現(xiàn)狀. 現(xiàn)代農(nóng)業(yè)科技, 2010, (6): 299–300. Ding H R, Hong L Z, Yang Z Q, Wang M W, Zhu X M, Liu C. Research status of saline-alkali land and its biological measures improvement., 2010, (6): 299–300 (in Chinese).

    [11] 任長忠, 胡躍高. 中國燕麥學(xué). 北京: 中國農(nóng)業(yè)出版社, 2013.Ren C Z, Hu Y G. Chinese Oatology. Beijing: China Agriculture Press, 2013 (in Chinese).

    [12] Quamruzzaman M, Manik S M N, Livermore M, Johnson P, Zhou M, Shabala S. Multidimensional screening and evaluation of morpho-physiological indices for salinity stress tolerance in wheat., 2022,208: 454–471.

    [13] Kingsbury R W, Epstein E. Selection for salt-resistant spring wheat., 1984, 24: 310–315.

    [14] Giaveno C D, Ribeiro R V, Souza G M, de Oliveira R F. Screening of tropical maize for salt stress tolerance., 2007, 7: 304–313.

    [15] Rasel M, Tahjib-Ul-Arif M, Hossain M A, Hassan L, Farzana S, Brestic M. Screening of salt-tolerant rice landraces by seedling stage phenotyping and dissecting biochemical determinants of tolerance mechanism., 2021, 40: 1853–1868.

    [16] 韓毅強(qiáng), 高亞梅, 杜艷麗, 張玉先, 杜吉到, 張文慧, 潘紹玉. 大豆耐鹽堿種質(zhì)資源鑒定. 中國油料作物學(xué)報, 2021, 43: 1016–1024. Han Y Q, Gao Y M, Du Y L, Zhang Y X, Du J D, Zhang W H, Pan S Y. Identification of saline-alkali tolerant germplasm resources of soybean during the whole growth stage., 2021, 43: 1016–1024 (in Chinese with English abstract).

    [17] 李春紅, 姚興東, 鞠寶韜, 朱明月, 王海英, 張惠君, 敖雪, 于翠梅, 謝甫綈, 宋書宏. 不同基因型大豆耐蔭性分析及其鑒定指標(biāo)的篩選. 中國農(nóng)業(yè)科學(xué), 2014, 47: 2927–2939. Li C H, Yao X D, Ju B T, Zhu M Y, Wang H Y, Zhang H J, Ao X, Yu C M, Xie F D, Song S H. Analysis of shade-tolerance and determination of shade-tolerance evaluation indicators in different soybean genotypes., 2014, 47: 2927–2939 (in Chinese with English abstract).

    [18] Bai J, Yan W, Wang Y, Yin Q, Liu J, Wight C, Ma B. Screening oat genotypes for tolerance to salinity and alkalinity., 2018, 9: 1302.

    [19] 陳新, 張宗文, 吳斌. 裸燕麥萌發(fā)期耐鹽性綜合評價與耐鹽種質(zhì)篩選. 中國農(nóng)業(yè)科學(xué), 2014, 47: 2038–2046.Chen X, Zhang Z W, Wu B. Comprehensive evaluation of salt tolerance and screening for salt tolerant accessions of naked oat (L.) at germination stage., 2014, 47: 2038–2046 (in Chinese with English abstract).

    [20] 楊科, 張保軍, 胡銀崗, 王淑華, 薛曉峰. 混合鹽堿脅迫對燕麥種子萌發(fā)及幼苗生理生化特性的影響. 干旱地區(qū)農(nóng)業(yè)研究, 2009, 27(3): 188–192. Yang K, Zhang B J, Hu Y G, Wang S H, Xue X F. Effects of complex saline-alkali stress on seeds germination and physiological and biochemical characteristics of oat seedlings., 2009, 27(3): 188–192 (in Chinese with English abstract).

    [21] 付鸞鴻, 于崧, 于立河, 薛盈文, 郭偉. 不同基因型燕麥萌發(fā)期耐鹽堿性分析及其鑒定指標(biāo)的篩選. 作物雜志, 2018, (6): 27–35. Fu L H, Yu S, Yu L H, Xue Y W, Guo W. Analysis of salt-alkali tolerance and screening of identification indexes of different Oat at germination stage., 2018, (6): 27–35 (in Chinese with English abstract).

    [22] 鄭殿生, 王曉鳴, 張京. 燕麥種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn). 北京: 中國農(nóng)業(yè)出版社, 2006. Zheng D S, Wang X M, Zhang J. Description Specification and Data Standard for Oat Germplasm Resources. Beijing: China Agriculture Press, 2006 (in Chinese).

    [23] 梁瀟, 侯向陽, 王艷榮, 任衛(wèi)波, 趙青山, 王照蘭, 李怡, 武自念. 羊草種質(zhì)資源耐鹽堿性綜合評價. 中國草地學(xué)報, 2019, 41(3): 1–9. Liang X, Hou X Y, Wang Y R, Ren W B, Zhao Q S, Wang Z L, Li Y, Wu Z N. Comprehensive evaluation on saline-alkali tolerance ofgermplasm resources., 2019, 41(3): 1–9 (in Chinese with English abstract).

    [24] 陳二影, 王潤豐, 秦嶺, 楊延兵, 黎飛飛, 張華文, 王海蓮, 劉賓, 孔清華, 管延安. 谷子芽期耐鹽堿綜合鑒定及評價. 作物學(xué)報, 2020, 46: 1591–1604. Chen E Y, Wang R F, Qin L, Yang Y B, Li F F, Zhang H W, Wang H L, Liu B, Kong Q H, Guan Y A. Comprehensive identification and evaluation of saline-alkali tolerance of Foxtail millet at bud stage., 2020, 46: 1591–1604 (in Chinese with English abstract).

    [25] 趙可夫. 植物抗鹽生理. 北京: 中國科學(xué)技術(shù)出版社, 1993. pp 3–4. Zhao K F. Plant Salt Resistance Physiology. Beijing: Science and Technology of China Press, 1993. pp 3–4 (in Chinese).

    [26] 王佺珍, 劉倩, 高婭妮, 柳旭. 植物對鹽堿脅迫的響應(yīng)機(jī)制研究進(jìn)展. 生態(tài)學(xué)報, 2017, 37: 5565–5577. Wang Q Z, Liu Q, Gao Y N, Liu X. Review on the mechanism of the responses to saline-alkali stress in plants., 2017, 37: 5565–5577 (in Chinese with English abstract).

    [27] Hosseini M K, Powell A A, Bingham I J. Comparison of the seed germination and early seedling growth of soybean in saline conditions., 2002, 12: 165–172.

    [28] 徐玲秀. 大豆種質(zhì)資源不同生育時期耐堿性評價. 東北農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 黑龍江哈爾濱, 2018. Xu L X.Evaluation of Alkali Tolerance of Soybean Germplasm Resources at Different Growth Stages.MS Thesis of Northeast Agricultural University, Harbin, Heilongjiang, China, 2018 (in Chinese with English abstract).

    [29] Zhu G Y, Kinet J M, Lutts S. Characterization of rice (L.) F3populations selected for salt resistance: I. Physiological behavior during vegetative growth., 2001, 121: 251–263.

    [30] Khan A A, Rao S A, McNeilly T. Assessment of salinity tolerance based upon seedling root growth response functions in maize (L.)., 2003, 131: 81–89.

    [31] Verma O P S, Yadava R B R. Salt tolerance of some oats (L.) varieties at germination and seedling stage., 1986, 156: 123–127.

    [32] Huang R. Research progress on plant tolerance to soil salinity and alkalinity in sorghum., 2018, 17: 739–746.

    [33] 馬曉軍, 金峰學(xué), 楊姍, 楊德光, 李曉輝. 作物耐鹽堿數(shù)量性狀基因座(QTL)定位. 分子植物育種, 2015, 13: 221–227. Ma X J, Jin F X, Yang S, Yang D G, Li X H. Mapping QTLs for salt and alkaline tolerance in crops., 2015, 13: 221–227 (in Chinese with English abstract).

    [34] Trachsel S, Messmer R, Stamp P, Ruta N, Hund A. QTLs for early vigor of tropical maize., 2010, 25: 91–103.

    [35] Tuyen D D, Zhang H M, Xu D H. Validation and high-resolution mapping of a major quantitative trait locus for alkaline salt tolerance in soybean using residual heterozygous line., 2013, 31: 79–86.

    [36] 斯日古楞, 張玉霞, 馬群, 鄧麗媛, 董航源. 8個飼用燕麥品種萌發(fā)期耐鹽性綜合評價. 草原與草坪, 2020, 40(6): 118–123. Siriguleng, Zhang Y X, Ma Q, Deng L Y, Dong H Y. Comprehensive evaluation of salt tolerance of 8 forage oat varieties., 2020, 40(6): 118–123 (in Chinese with English abstract).

    [37] Khosravi S, Azizinezhad R, Baghizadeh A, Maleki M. Evaluation and comparison of drought tolerance in some wild diploid populations, tetraploid and hexaploid cultivars of wheat using stress tolerance indices., 2020, 115: 105–112.

    [38] Zhang W, Zhao Z, Bai G, Fu F. Study and evaluation of drought resistance of different genotype maize inbred lines., 2008, 2: 428–434.

    [39] Li G, Zhang H Y, Ji L, Zhao E, Liu J, Li L, Zhang J. Comprehensive evaluation on drought-resistance of different soybean varieties., 2006, 17: 2408–2412.

    [40] 張婷婷, 于崧, 于立河, 李琳, 金珊珊, 郭建華, 張靜. 松嫩平原春小麥耐鹽堿性鑒定及品種(系)篩選. 麥類作物學(xué)報, 2016, 36: 1008–1019. Zhang T T, Yu S, Yu L H, Li L, Jin S S, Guo J H, Zhang J. Saline-alkali tolerance identification and varieties (lines) screening of spring wheat in Songnen Plain., 2016, 36: 1008–1019 (in Chinese with English abstract).

    [41] 段雅娟, 曹士亮, 于滔, 李文躍, 楊耿斌, 王成波, 劉寶民, 劉長華. 玉米自交系萌發(fā)期耐鹽性鑒定. 作物雜志, 2022, (1): 213–219. Duan Y J, Cao S L, Yu T, Li W Y, Yang G B, Liu B M, Liu C H. Identification of salt tolerance during germination of maize in-bred lines., 2022, (1): 213–219 (in Chinese with English abstract).

    [42] 虞曉芬, 傅玳. 多指標(biāo)綜合評價方法綜述. 統(tǒng)計與決策, 2004, (11): 119–121.Yu X F, Fu D. Summary of multi-index comprehensive evaluation methods., 2004, (11): 119–121 (in Chinese).

    Identification and evaluation of salt-alkali tolerance and screening of salt-alkali tolerant germplasm of oat (L.)

    ZHANG Jing1, GAO Wen-Bo1, YAN Lin1, ZHANG Zong-Wen1,2, ZHOU Hai-Tao3, and WU Bin1,*

    1Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China;2Biodiversity International, Beijing 100081, China;3Zhangjiakou Academy of Agricultural Sciences, Zhangjiakou 075000, Hebei, China

    Affected by climate change and human production activities, the world’s salt-alkali land is expanding, and soil salinization has become an important factor limiting the development of agricultural production. Oat is a crop with the strong saline-alkali tolerance. To evaluate the salinity tolerance of oat germplasm resources in China, 485 accessions were stress treated with 125mmol L–1NaCl, Na2SO4,and NaHCO3(1:1:1 molar concentration) solution at germination stage. Eight growth indexes, including germination potential, germination rate, root length, bud length, root fresh weight, bud fresh weight, root dry weight, and bud dry weight, were identified at the germination stage of oats. A comprehensive evaluation and screening of oat germplasm resources for salt-alkali tolerance was performed by the correlation analysis, principal component analysis, membership function analysis, and cluster analysis. The results showed that salt-alkali inhibited all eight indexes identified, and a significant positive correlation was observed between the affiliation function values of the indexes under salt-alkali tolerance and with the comprehensive evaluation value. Eight evaluation indexes were converted into two comprehensive indexes by principal component analysis, with a cumulative variance contribution of 76.926%.The membership function analysis combined with cluster analysis screened a total of two oat accessions with high salt-alkaline tolerance (Oat 1606 and Heikowski 596) and classified 485 accessions into five classes, including two highly saline-alkaline tolerance, 49 salt-alkaline-tolerant accessions, 147 medium saline-alkaline tolerance accessions, 129 sensitive saline-alkaline accessions, and 158 high saline-alkaline sensitive accessions.The results of correlation analysis, principal component analysis, and stepwise multiple regression analysis, and the bud length was screened as the preferred indicator for the screening of salt-alkali tolerance in germinating oats, followed by root fresh weight, bud fresh weight, germination potential, and germination rate, which were also important indexes for the comprehensive evaluation and screening of salt-alkali tolerance in germinating oats.

    L.; germplasm resources; germination stage; salt-alkali tolerance; comprehensive evaluation

    10.3724/SP.J.1006.2023.21032

    本研究由財政部和農(nóng)業(yè)農(nóng)村部國家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(CARS-07-A-1), 河北省重點研發(fā)計劃項目(21326305D), 中國農(nóng)業(yè)科學(xué)院科技創(chuàng)新工程項目和農(nóng)作物種質(zhì)資源保護(hù)與利用專項(2020NWB036-06)資助。

    This study was supported by the China Agriculture Research System of MOF and MARA (CARS-07-A-1), the Key Research and Development Program of Hebei Province (21326305D), the Agricultural Science and Technology Innovation Program (ASTIP) in CAAS, and the Crop Germplasm Resources Protection (2020NWB036-06).

    吳斌, E-mail: wubin03@caas.cn

    E-mail: zhangjing12290218@163.com

    2022-04-28;

    2022-10-10;

    2022-11-02.

    URL: https://kns.cnki.net/kcms/detail/11.1809.S.20221101.1656.004.html

    This is an open access article under the CC BY-NC-ND license (http://creativecommons.org/licenses/by-nc-nd/4.0/).

    猜你喜歡
    鹽堿發(fā)芽勢燕麥
    可嚼燕麥,營養(yǎng)打折
    中老年保健(2022年5期)2022-11-25 14:16:14
    鹽堿荒灘的“底色之變”
    金橋(2021年11期)2021-11-20 06:37:14
    燕麥的栽培技術(shù)
    晉??笛帑?守護(hù)您的健康
    鹽堿地區(qū)南美白對蝦、青蝦混養(yǎng)技術(shù)
    不同濃度NaCl溶液對5種草坪草種子發(fā)芽勢的影響
    5種植物激素對黑豆種子萌發(fā)的影響
    種子(2019年6期)2019-07-13 07:45:16
    耐旱耐鹽棉花種子萌發(fā)的影響因子分析初報
    種子(2019年4期)2019-05-28 02:04:14
    種子發(fā)芽勢對作物田間出苗率的重要性分析
    行了,我像所有的他們一樣
    延河(2017年7期)2017-07-19 21:01:10
    纯流量卡能插随身wifi吗| 国产精品偷伦视频观看了| 精品99又大又爽又粗少妇毛片| 久久午夜综合久久蜜桃| 国产极品天堂在线| 80岁老熟妇乱子伦牲交| 精品少妇一区二区三区视频日本电影 | 人人妻人人添人人爽欧美一区卜| 国产极品天堂在线| 亚洲国产日韩一区二区| 伊人久久大香线蕉亚洲五| 日日爽夜夜爽网站| 亚洲美女视频黄频| 在线天堂最新版资源| 日韩精品免费视频一区二区三区| 亚洲精品成人av观看孕妇| 国语对白做爰xxxⅹ性视频网站| av免费在线看不卡| 纵有疾风起免费观看全集完整版| 少妇的逼水好多| 亚洲国产欧美日韩在线播放| 色播在线永久视频| 国产精品人妻久久久影院| 一级,二级,三级黄色视频| 不卡视频在线观看欧美| 精品亚洲乱码少妇综合久久| 亚洲精品日本国产第一区| 在线观看国产h片| 亚洲欧美成人综合另类久久久| 香蕉丝袜av| 国产精品国产av在线观看| 香蕉精品网在线| 中文字幕精品免费在线观看视频| 在线观看三级黄色| 久久精品亚洲av国产电影网| 久久久国产精品麻豆| 亚洲欧美色中文字幕在线| 亚洲欧美一区二区三区国产| 国产片特级美女逼逼视频| 亚洲综合色网址| 人成视频在线观看免费观看| 精品亚洲成国产av| 久久久久国产一级毛片高清牌| 亚洲精品一二三| 啦啦啦在线观看免费高清www| 亚洲精品久久久久久婷婷小说| 亚洲av综合色区一区| 最近中文字幕2019免费版| 狠狠精品人妻久久久久久综合| 中文字幕最新亚洲高清| 亚洲av欧美aⅴ国产| tube8黄色片| 午夜老司机福利剧场| 男女啪啪激烈高潮av片| 久久精品熟女亚洲av麻豆精品| 久久 成人 亚洲| 香蕉国产在线看| 亚洲成人av在线免费| 欧美日韩亚洲高清精品| 综合色丁香网| 国产精品成人在线| 国产精品不卡视频一区二区| 欧美精品亚洲一区二区| 亚洲av福利一区| 在线天堂最新版资源| 精品亚洲成国产av| 欧美国产精品va在线观看不卡| 制服诱惑二区| 亚洲男人天堂网一区| 韩国高清视频一区二区三区| www.熟女人妻精品国产| 亚洲av欧美aⅴ国产| 在线观看一区二区三区激情| 久久 成人 亚洲| 在线天堂最新版资源| 三上悠亚av全集在线观看| 久久影院123| 国产免费福利视频在线观看| 丁香六月天网| 精品久久久精品久久久| 看非洲黑人一级黄片| 日韩在线高清观看一区二区三区| 久久久久国产精品人妻一区二区| 日韩三级伦理在线观看| 欧美日本中文国产一区发布| av电影中文网址| 最新的欧美精品一区二区| 999精品在线视频| 成人18禁高潮啪啪吃奶动态图| 最新的欧美精品一区二区| 热99久久久久精品小说推荐| 秋霞在线观看毛片| 搡老乐熟女国产| 欧美亚洲 丝袜 人妻 在线| 亚洲视频免费观看视频| 精品少妇一区二区三区视频日本电影 | 亚洲欧美精品自产自拍| 亚洲少妇的诱惑av| 91精品国产国语对白视频| 99香蕉大伊视频| 亚洲精品美女久久久久99蜜臀 | 人人妻人人添人人爽欧美一区卜| 欧美亚洲 丝袜 人妻 在线| 叶爱在线成人免费视频播放| 韩国精品一区二区三区| 久久久久国产精品人妻一区二区| 亚洲精品自拍成人| 免费少妇av软件| 飞空精品影院首页| 国产成人精品婷婷| 亚洲av电影在线进入| 国产免费又黄又爽又色| 国产片特级美女逼逼视频| 精品午夜福利在线看| 欧美xxⅹ黑人| 在线观看人妻少妇| 午夜免费鲁丝| 免费观看a级毛片全部| 久久久国产精品麻豆| 91成人精品电影| 国产精品熟女久久久久浪| xxx大片免费视频| 寂寞人妻少妇视频99o| 男人添女人高潮全过程视频| 水蜜桃什么品种好| 久久久久人妻精品一区果冻| 大香蕉久久成人网| 国产免费现黄频在线看| 日韩制服丝袜自拍偷拍| 久久午夜福利片| 少妇的丰满在线观看| 亚洲欧美精品自产自拍| 久久精品久久精品一区二区三区| 免费高清在线观看视频在线观看| 高清视频免费观看一区二区| 女人被躁到高潮嗷嗷叫费观| 天天躁日日躁夜夜躁夜夜| 七月丁香在线播放| 欧美精品人与动牲交sv欧美| 午夜福利一区二区在线看| 亚洲欧美中文字幕日韩二区| 有码 亚洲区| 制服丝袜香蕉在线| 欧美97在线视频| 9热在线视频观看99| 日韩大片免费观看网站| 午夜福利网站1000一区二区三区| 熟女av电影| 亚洲精品国产av蜜桃| 大片免费播放器 马上看| 少妇的丰满在线观看| 在线观看免费日韩欧美大片| 大香蕉久久网| av在线老鸭窝| 久久精品夜色国产| 啦啦啦视频在线资源免费观看| 日日撸夜夜添| 亚洲av中文av极速乱| 日韩,欧美,国产一区二区三区| 久久久久久人妻| 国产在视频线精品| 一级毛片电影观看| 欧美人与性动交α欧美精品济南到 | 亚洲成人手机| 色哟哟·www| 欧美 亚洲 国产 日韩一| 国产成人精品久久久久久| 宅男免费午夜| 欧美成人午夜免费资源| 亚洲精品国产av成人精品| 亚洲欧洲精品一区二区精品久久久 | 亚洲美女视频黄频| 美女大奶头黄色视频| 亚洲人成网站在线观看播放| 亚洲精品国产av成人精品| 自拍欧美九色日韩亚洲蝌蚪91| 大香蕉久久成人网| 亚洲人成77777在线视频| 日本欧美视频一区| 麻豆乱淫一区二区| 热re99久久国产66热| 国产精品亚洲av一区麻豆 | 天天躁夜夜躁狠狠躁躁| 老司机亚洲免费影院| 黑人欧美特级aaaaaa片| 最黄视频免费看| 中文天堂在线官网| 亚洲成av片中文字幕在线观看 | 久久99热这里只频精品6学生| 久久久久久伊人网av| 国产女主播在线喷水免费视频网站| 宅男免费午夜| 国产在线一区二区三区精| 亚洲精品日本国产第一区| 久久ye,这里只有精品| 欧美bdsm另类| 男女免费视频国产| 亚洲第一青青草原| 国产精品久久久久久精品古装| 婷婷色综合大香蕉| 天堂中文最新版在线下载| 国产精品香港三级国产av潘金莲| 变态另类成人亚洲欧美熟女 | 久久婷婷成人综合色麻豆| 亚洲欧洲精品一区二区精品久久久| 日本免费a在线| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区激情短视频| 一a级毛片在线观看| 自线自在国产av| 国产亚洲精品久久久久久毛片| 韩国av一区二区三区四区| 精品一区二区三区av网在线观看| 可以免费在线观看a视频的电影网站| 一个人观看的视频www高清免费观看 | 乱人伦中国视频| 国产乱人伦免费视频| 久久国产亚洲av麻豆专区| 日韩欧美在线二视频| 精品一区二区三区四区五区乱码| 午夜福利在线观看吧| 少妇裸体淫交视频免费看高清 | 久久久久久大精品| 国产av一区二区精品久久| 亚洲伊人色综图| 欧美亚洲日本最大视频资源| 精品福利观看| 一本大道久久a久久精品| 亚洲欧美日韩高清在线视频| 免费av中文字幕在线| 亚洲九九香蕉| 欧美老熟妇乱子伦牲交| 国产不卡一卡二| 99久久人妻综合| 亚洲欧美日韩另类电影网站| 99国产精品一区二区三区| 成人精品一区二区免费| 国内久久婷婷六月综合欲色啪| 亚洲成国产人片在线观看| 丝袜人妻中文字幕| а√天堂www在线а√下载| 国产三级黄色录像| 亚洲全国av大片| 91精品国产国语对白视频| 首页视频小说图片口味搜索| 国产成人影院久久av| 人成视频在线观看免费观看| 久久中文字幕一级| 一夜夜www| 天天添夜夜摸| 国产欧美日韩一区二区三| 国产xxxxx性猛交| 一进一出抽搐gif免费好疼 | 丝袜在线中文字幕| 日韩国内少妇激情av| av国产精品久久久久影院| 9色porny在线观看| 亚洲欧洲精品一区二区精品久久久| 亚洲成国产人片在线观看| 国产极品粉嫩免费观看在线| 一级片免费观看大全| 少妇裸体淫交视频免费看高清 | 欧美另类亚洲清纯唯美| 免费av中文字幕在线| 亚洲男人的天堂狠狠| 97超级碰碰碰精品色视频在线观看| 国产色视频综合| 亚洲五月天丁香| 午夜亚洲福利在线播放| 久久久精品欧美日韩精品| 女生性感内裤真人,穿戴方法视频| 亚洲第一青青草原| 久久精品国产清高在天天线| 91麻豆精品激情在线观看国产 | 黄频高清免费视频| 国产精品爽爽va在线观看网站 | 久久久久久久午夜电影 | 这个男人来自地球电影免费观看| 日本免费a在线| 香蕉丝袜av| 国产精品亚洲av一区麻豆| 99在线视频只有这里精品首页| 中文欧美无线码| a在线观看视频网站| 一级片免费观看大全| 天堂中文最新版在线下载| av在线天堂中文字幕 | 99热只有精品国产| 他把我摸到了高潮在线观看| 男女下面进入的视频免费午夜 | 成人18禁高潮啪啪吃奶动态图| 欧洲精品卡2卡3卡4卡5卡区| 欧洲精品卡2卡3卡4卡5卡区| 欧美人与性动交α欧美精品济南到| 90打野战视频偷拍视频| 亚洲av日韩精品久久久久久密| 在线天堂中文资源库| 亚洲久久久国产精品| 欧美黄色片欧美黄色片| 国产野战对白在线观看| ponron亚洲| 激情在线观看视频在线高清| 亚洲色图av天堂| 黄色成人免费大全| 水蜜桃什么品种好| 精品少妇一区二区三区视频日本电影| 国产乱人伦免费视频| 一二三四在线观看免费中文在| 国产精品一区二区精品视频观看| 亚洲专区国产一区二区| 精品少妇一区二区三区视频日本电影| 精品国内亚洲2022精品成人| 国产精品野战在线观看 | 亚洲成人久久性| 久久久久亚洲av毛片大全| 女人高潮潮喷娇喘18禁视频| 久久人人精品亚洲av| 国产又爽黄色视频| 亚洲精品久久成人aⅴ小说| 国产欧美日韩一区二区三| 日韩欧美三级三区| 91精品三级在线观看| 啦啦啦免费观看视频1| 97碰自拍视频| 色在线成人网| 国产精品久久久人人做人人爽| 丰满的人妻完整版| 欧美激情 高清一区二区三区| 免费在线观看完整版高清| 99国产精品免费福利视频| 88av欧美| 人人妻人人爽人人添夜夜欢视频| 岛国视频午夜一区免费看| 在线免费观看的www视频| 国产三级黄色录像| 欧美黄色片欧美黄色片| 琪琪午夜伦伦电影理论片6080| 日本黄色日本黄色录像| 日本黄色日本黄色录像| 五月开心婷婷网| 久久久国产精品麻豆| 国产av精品麻豆| 嫩草影院精品99| 露出奶头的视频| 叶爱在线成人免费视频播放| 女人被躁到高潮嗷嗷叫费观| ponron亚洲| 涩涩av久久男人的天堂| 无遮挡黄片免费观看| 国产精品久久久av美女十八| av在线播放免费不卡| 日韩欧美一区视频在线观看| 黄色视频不卡| 琪琪午夜伦伦电影理论片6080| 欧美日本中文国产一区发布| 十八禁网站免费在线| 黑丝袜美女国产一区| 久久精品成人免费网站| 亚洲国产欧美一区二区综合| 国产精品偷伦视频观看了| 一级作爱视频免费观看| 欧美黄色淫秽网站| 性少妇av在线| 亚洲第一欧美日韩一区二区三区| 国产人伦9x9x在线观看| 怎么达到女性高潮| 日韩国内少妇激情av| 亚洲人成网站在线播放欧美日韩| 动漫黄色视频在线观看| 久久人人97超碰香蕉20202| 亚洲 国产 在线| 亚洲国产毛片av蜜桃av| 日本 av在线| 亚洲av美国av| 久久中文字幕一级| 亚洲欧美一区二区三区黑人| 午夜免费激情av| 国产成人av教育| 亚洲国产中文字幕在线视频| 欧美乱色亚洲激情| 欧美精品啪啪一区二区三区| 精品一区二区三区视频在线观看免费 | 在线观看午夜福利视频| 久久香蕉激情| 久久久国产成人精品二区 | 身体一侧抽搐| 亚洲伊人色综图| 久久久国产欧美日韩av| 一a级毛片在线观看| 成人三级黄色视频| 两个人免费观看高清视频| 久久久国产精品麻豆| 亚洲五月天丁香| 热99re8久久精品国产| 成人三级做爰电影| 丰满人妻熟妇乱又伦精品不卡| 免费在线观看日本一区| 免费在线观看视频国产中文字幕亚洲| 国产精品99久久99久久久不卡| 黄网站色视频无遮挡免费观看| 午夜精品国产一区二区电影| 亚洲第一av免费看| 一进一出好大好爽视频| 国产av精品麻豆| 男人的好看免费观看在线视频 | 两性午夜刺激爽爽歪歪视频在线观看 | av电影中文网址| 久久香蕉精品热| 久久欧美精品欧美久久欧美| 久久久久久人人人人人| 日韩成人在线观看一区二区三区| av视频免费观看在线观看| av在线播放免费不卡| 国产三级黄色录像| 亚洲欧美一区二区三区黑人| 国产熟女午夜一区二区三区| 狂野欧美激情性xxxx| 国产99久久九九免费精品| 欧美不卡视频在线免费观看 | 在线观看一区二区三区| 国产精品一区二区三区四区久久 | 91精品国产国语对白视频| 亚洲国产精品合色在线| 精品国产超薄肉色丝袜足j| 天天影视国产精品| 亚洲精品一二三| 欧美激情高清一区二区三区| 亚洲一区二区三区欧美精品| 亚洲成国产人片在线观看| 可以在线观看毛片的网站| 热99国产精品久久久久久7| 精品国产超薄肉色丝袜足j| 老熟妇仑乱视频hdxx| 美女大奶头视频| 久久精品影院6| 久久天躁狠狠躁夜夜2o2o| 国产精品亚洲av一区麻豆| 制服人妻中文乱码| 亚洲情色 制服丝袜| 真人做人爱边吃奶动态| 一级黄色大片毛片| 亚洲色图 男人天堂 中文字幕| 女警被强在线播放| 高清欧美精品videossex| 免费少妇av软件| 高清毛片免费观看视频网站 | 中文字幕精品免费在线观看视频| 在线观看免费视频网站a站| 久久天躁狠狠躁夜夜2o2o| xxxhd国产人妻xxx| 国产精品一区二区精品视频观看| 99久久人妻综合| 老汉色av国产亚洲站长工具| 视频在线观看一区二区三区| 99久久精品国产亚洲精品| 色尼玛亚洲综合影院| av电影中文网址| 麻豆国产av国片精品| 亚洲精品国产色婷婷电影| 精品一区二区三卡| 国产精品久久久久久人妻精品电影| 香蕉丝袜av| 日韩视频一区二区在线观看| 欧美精品一区二区免费开放| www.熟女人妻精品国产| 手机成人av网站| 一级a爱片免费观看的视频| 色综合站精品国产| 十八禁人妻一区二区| 狂野欧美激情性xxxx| 精品卡一卡二卡四卡免费| 国产精品99久久99久久久不卡| 国产亚洲av高清不卡| 老熟妇乱子伦视频在线观看| 激情视频va一区二区三区| 一级毛片女人18水好多| 最近最新中文字幕大全电影3 | 国产精品 国内视频| 亚洲aⅴ乱码一区二区在线播放 | 亚洲av五月六月丁香网| 纯流量卡能插随身wifi吗| 免费在线观看完整版高清| 亚洲精品中文字幕一二三四区| 国产区一区二久久| 免费一级毛片在线播放高清视频 | 日本免费a在线| 视频区图区小说| 亚洲精品av麻豆狂野| 欧美不卡视频在线免费观看 | 在线观看一区二区三区激情| 久久精品影院6| 亚洲五月婷婷丁香| 两个人看的免费小视频| 中文字幕精品免费在线观看视频| 久久久久亚洲av毛片大全| 久久久久久久午夜电影 | 制服人妻中文乱码| 黄色成人免费大全| 两个人免费观看高清视频| 国产成人精品久久二区二区91| 欧美黄色淫秽网站| 亚洲色图综合在线观看| 在线观看舔阴道视频| 自线自在国产av| 国产免费现黄频在线看| 久久久精品国产亚洲av高清涩受| 在线观看免费高清a一片| 亚洲 国产 在线| av在线播放免费不卡| 国产av精品麻豆| 午夜免费激情av| 精品一区二区三区av网在线观看| 一边摸一边做爽爽视频免费| 亚洲精品一二三| 成年人黄色毛片网站| 免费在线观看日本一区| 女人爽到高潮嗷嗷叫在线视频| 亚洲自拍偷在线| 丝袜美腿诱惑在线| 成熟少妇高潮喷水视频| 在线观看免费高清a一片| 97人妻天天添夜夜摸| 一区二区日韩欧美中文字幕| 色精品久久人妻99蜜桃| 法律面前人人平等表现在哪些方面| 国产精品一区二区精品视频观看| 少妇被粗大的猛进出69影院| www.www免费av| 级片在线观看| 亚洲欧美日韩另类电影网站| 丰满的人妻完整版| 18禁裸乳无遮挡免费网站照片 | 另类亚洲欧美激情| 这个男人来自地球电影免费观看| 亚洲情色 制服丝袜| 12—13女人毛片做爰片一| 国产一区二区激情短视频| 99久久国产精品久久久| 精品卡一卡二卡四卡免费| 又紧又爽又黄一区二区| www.999成人在线观看| 99久久国产精品久久久| 日韩国内少妇激情av| 亚洲精品久久成人aⅴ小说| a级毛片在线看网站| 两个人看的免费小视频| 国产高清视频在线播放一区| 久久久水蜜桃国产精品网| 在线看a的网站| 国产高清国产精品国产三级| 身体一侧抽搐| 国产一区二区三区在线臀色熟女 | 性欧美人与动物交配| 91在线观看av| 精品国产超薄肉色丝袜足j| 国产亚洲欧美在线一区二区| 老司机午夜十八禁免费视频| 久久国产精品人妻蜜桃| 国产亚洲精品久久久久久毛片| 宅男免费午夜| 日本一区二区免费在线视频| 精品卡一卡二卡四卡免费| 国产又爽黄色视频| 十八禁人妻一区二区| 制服人妻中文乱码| 久热这里只有精品99| 欧美乱妇无乱码| 美女国产高潮福利片在线看| 亚洲精品一二三| 国产欧美日韩一区二区精品| 精品国产国语对白av| 在线视频色国产色| 国产精品 国内视频| x7x7x7水蜜桃| 久久午夜综合久久蜜桃| 老熟妇仑乱视频hdxx| 午夜福利免费观看在线| 桃红色精品国产亚洲av| 看黄色毛片网站| 成人三级黄色视频| 嫁个100分男人电影在线观看| 啦啦啦在线免费观看视频4| 国产在线观看jvid| a级毛片在线看网站| 国产在线精品亚洲第一网站| 黄色 视频免费看| 久久精品国产亚洲av高清一级| 国产区一区二久久| 一级黄色大片毛片| 久久久国产精品麻豆| 美女扒开内裤让男人捅视频| 夜夜看夜夜爽夜夜摸 | 日本wwww免费看| 日韩欧美在线二视频| 大型黄色视频在线免费观看| 满18在线观看网站| 国产日韩一区二区三区精品不卡| 99久久久亚洲精品蜜臀av| videosex国产| 校园春色视频在线观看| 亚洲三区欧美一区| 美女高潮到喷水免费观看| 91国产中文字幕| 欧美性长视频在线观看| 精品国产超薄肉色丝袜足j| 在线视频色国产色| 成人三级黄色视频| 亚洲黑人精品在线| 少妇的丰满在线观看| 女人爽到高潮嗷嗷叫在线视频| 日韩成人在线观看一区二区三区| 亚洲av日韩精品久久久久久密| 国产单亲对白刺激| 色老头精品视频在线观看| 亚洲视频免费观看视频|