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

    谷子芽期耐鹽堿綜合鑒定及評(píng)價(jià)

    2020-09-14 07:22:52陳二影王潤(rùn)豐楊延兵黎飛飛張華文王海蓮孔清華管延安
    作物學(xué)報(bào) 2020年10期
    關(guān)鍵詞:芽長(zhǎng)鹽堿耐鹽性

    陳二影 王潤(rùn)豐 秦 嶺 楊延兵 黎飛飛 張華文 王海蓮 劉 賓 孔清華 管延安,,*

    谷子芽期耐鹽堿綜合鑒定及評(píng)價(jià)

    陳二影1王潤(rùn)豐1秦 嶺1楊延兵1黎飛飛1張華文1王海蓮1劉 賓1孔清華2管延安1,2,*

    1山東省農(nóng)業(yè)科學(xué)院作物研究所/ 山東省特色作物工程實(shí)驗(yàn)室, 山東濟(jì)南 250100;2山東師范大學(xué), 山東濟(jì)南 250014

    以全國(guó)主推的53個(gè)谷子品種為材料, 在100 mmol L-1混合鹽堿(NaCl∶NaHCO3= 4∶1)脅迫下研究了不同谷子品種的耐鹽堿性。結(jié)果表明, 在鹽堿脅迫下, 53個(gè)谷子品種的發(fā)芽勢(shì)、發(fā)芽率、根長(zhǎng)、芽長(zhǎng)、根鮮重和芽鮮重均受到不同程度的抑制, 以對(duì)根長(zhǎng)的影響最大; 相對(duì)發(fā)芽勢(shì)與相對(duì)發(fā)芽率、相對(duì)根長(zhǎng)與相對(duì)芽長(zhǎng)及相對(duì)根鮮重與相對(duì)芽鮮重均呈顯著或極顯著正相關(guān)。通過(guò)主成分分析將14個(gè)單項(xiàng)性狀指標(biāo)轉(zhuǎn)化為4個(gè)主成分, 累積貢獻(xiàn)率為90.4%; 以4個(gè)主成分的得分值通過(guò)隸屬函數(shù)分析獲得不同品種耐鹽堿的綜合得分值, 并通過(guò)聚類(lèi)分析將53個(gè)谷子品種劃分為6種耐鹽堿類(lèi)型, 其中強(qiáng)耐鹽堿品種2個(gè), 耐鹽堿品種16個(gè), 中間型品種17個(gè), 鹽堿敏感品種6個(gè), 不耐鹽堿品種9個(gè)和極不耐鹽堿的品種3個(gè)。同時(shí)利用回歸分析建立了可用于評(píng)價(jià)谷子耐鹽堿性的回歸方程¢= 0.298 + 0.0372+ 0.1443+ 0.0186+ 0.2097-0.1839+ 0.11511-0.20112+ 0.11213-0.10114+ 0.28415, 相對(duì)發(fā)芽率、根長(zhǎng)鹽害率、芽長(zhǎng)鹽害率和根冠比鹽害率可以作為谷子耐鹽堿性的評(píng)價(jià)指標(biāo)。

    谷子; 鹽堿脅迫; 主成分分析; 隸屬函數(shù); 回歸分析

    土壤鹽堿化是影響作物生長(zhǎng)發(fā)育和產(chǎn)量的重要限制因子[1-3]。全球約有30%左右的耕地受鹽堿的影響[4]。我國(guó)具有大量的鹽堿土地, 總面積約9913萬(wàn)公頃[5], 且呈逐年遞增的趨勢(shì)。芽期和苗期是作物受鹽堿危害的敏感時(shí)期, 提高作物芽期和苗期的耐鹽堿性是應(yīng)對(duì)鹽堿的重要途徑, 篩選具有芽期和苗期強(qiáng)耐鹽堿性的品種是重要的手段。景宇鵬等[6]利用主成分分析在玉米萌發(fā)和幼苗期進(jìn)行了耐鹽性品種的篩選和綜合評(píng)價(jià)。信彩云等[7]在水稻苗期先通過(guò)主成分分析確定綜合性狀指標(biāo), 后通過(guò)回歸方程確定了幼苗鮮重、葉綠素含量和SOD酶活作為單項(xiàng)耐鹽性的鑒定指標(biāo)。于瑩等[8]研究發(fā)現(xiàn), 模糊隸屬函數(shù)在玉米萌發(fā)期可作為耐鹽性分級(jí)評(píng)價(jià)的方法。高春華等[9]的研究指出, 主成分分析和模糊隸屬函數(shù)在高粱芽期耐鹽性的評(píng)價(jià)中存在一致性, 已在大豆[10]上共同用于苗期耐鹽性品種的篩選和評(píng)價(jià)。

    谷子[(L.) Beauv.]是起源于中國(guó)的傳統(tǒng)糧食作物, 具有較強(qiáng)的抗旱、耐瘠薄的能力[11-13], 是典型的環(huán)境友好型作物。同時(shí), 谷子具有較強(qiáng)的抗鹽性[14-16], 培育耐鹽堿性強(qiáng)的品種是應(yīng)對(duì)鹽堿災(zāi)害的重要途徑。在中性鹽NaCl脅迫下, 不同谷子品種發(fā)芽率、根長(zhǎng)和芽長(zhǎng)存在品種間差異[17], 在低鹽脅迫下, 發(fā)芽率和發(fā)芽勢(shì)無(wú)顯著變化[18], 在高鹽脅迫下, 根長(zhǎng)和芽長(zhǎng)均受到一定的抑制[19]; 在堿性鹽脅迫下谷子發(fā)芽率和發(fā)芽勢(shì)顯著降低[20]。通過(guò)主成分分析確定, 發(fā)芽率、發(fā)芽指數(shù)和相對(duì)芽長(zhǎng)可作為谷子芽期耐鹽性評(píng)價(jià)的關(guān)鍵指標(biāo)[21]。但前人研究多集中在單一的中性鹽或堿性鹽條件下, 而鹽堿地實(shí)際生產(chǎn)中存在鹽堿共同脅迫的問(wèn)題, 有關(guān)谷子混合鹽堿脅迫研究尚未見(jiàn)報(bào)道; 且谷子芽期耐鹽性的鑒定方法較為單一, 缺乏評(píng)價(jià)的系統(tǒng)性和準(zhǔn)確性。因此本研究以全國(guó)不同生態(tài)區(qū)主推谷子品種為材料, 在混合鹽堿脅迫下通過(guò)主成分分析、模糊隸屬函數(shù)、聚類(lèi)分析和回歸分析進(jìn)行谷子芽期耐鹽堿的綜合評(píng)價(jià), 以期為篩選芽期強(qiáng)耐鹽堿的谷子品種和評(píng)價(jià)指標(biāo)提供材料和方法。

    1 材料與方法

    1.1 試驗(yàn)材料

    采用全國(guó)生產(chǎn)上主推的53個(gè)谷子品種為材料,見(jiàn)表1。

    1.2 試驗(yàn)設(shè)計(jì)

    2018年在山東省農(nóng)業(yè)科學(xué)院人工氣候室中進(jìn)行培養(yǎng), 每個(gè)品種挑選均勻、大小一致的50粒種子, 先用75%的酒精消毒15 min, 然后用去離子水沖洗3次, 吸水紙吸干后放入鋪有2層濾紙的培養(yǎng)皿中培養(yǎng)。以濟(jì)谷19、濟(jì)谷20、冀谷41、豫谷18、保谷23、秦谷3號(hào)、晉谷46、龍谷31共8個(gè)品種進(jìn)行預(yù)試驗(yàn), 確定篩選的最適鹽堿濃度。本試驗(yàn)根據(jù)濱海鹽堿地的特征[22-23]確定篩選的鹽堿比為4∶1 (NaCl∶NaHCO3), 以此鹽堿比設(shè)定40、60、80、100、120、140、160、180和200 mmol L-1不同的鹽堿濃度, 蒸餾水做對(duì)照, 各處理重復(fù)4次。根據(jù)發(fā)芽勢(shì)、發(fā)芽率、根長(zhǎng)和葉長(zhǎng)確定各指標(biāo)與對(duì)照差異最顯著的100 mmol L-1鹽堿濃度為篩選濃度。人工氣候室設(shè)置晝夜溫度28℃/22℃, 濕度65%, 晝夜時(shí)長(zhǎng)各為12 h, 光照強(qiáng)度為240mmol m-2s-1。

    1.3 測(cè)定項(xiàng)目與方法

    培養(yǎng)第4天和第10天統(tǒng)計(jì)各品種發(fā)芽率, 并在第10天測(cè)定根長(zhǎng)、芽長(zhǎng)和稱(chēng)取鮮重, 每培養(yǎng)皿稱(chēng)量10株, 求平均值。按以下公式進(jìn)行各性狀指標(biāo)的計(jì)算:

    發(fā)芽勢(shì)(%) = 第4天發(fā)芽種子數(shù)/種子總數(shù)′100%;

    發(fā)芽率(%) = 第10天發(fā)芽種子數(shù)/種子總數(shù)′100%;

    相對(duì)發(fā)芽勢(shì)(%) = 處理發(fā)芽勢(shì)/對(duì)照發(fā)芽勢(shì)′100%;

    相對(duì)發(fā)芽率(%) = 處理發(fā)芽率/對(duì)照發(fā)芽率′100%;

    相對(duì)根長(zhǎng)(%) = 處理根長(zhǎng)/對(duì)照根長(zhǎng)′100%;

    相對(duì)芽長(zhǎng)(%) = 處理芽長(zhǎng)/對(duì)照芽長(zhǎng)′100%;

    相對(duì)根鮮重(%) = 處理根鮮重/對(duì)照根鮮重′100%;

    相對(duì)芽鮮重(%) = 處理芽鮮重/對(duì)照芽鮮重′100%;

    相對(duì)根冠比(%) = 相對(duì)根鮮重/相對(duì)芽鮮重′100%;

    發(fā)芽勢(shì)鹽害率(%) = (對(duì)照發(fā)芽勢(shì)-處理發(fā)芽勢(shì))/對(duì)照發(fā)芽勢(shì)′100%;

    表1 供試谷子品種及育成單位

    SDAAS: Shandong Academy of Agricultural Sciences; CAAS: Chinese Academy of Agricultural Sciences; HBAFS: Hebei Academy of Agriculture and Forestry Sciences; LCAAS: Liaocheng Academy of Agricultural Sciences; TAAAS: Taian Academy of Agricultural Sciences; HDAAS: Handan Academy of Agricultural Sciences; HNAAS: Henan Academy of Agricultural Sciences; AYAAS: Anyang Academy of Agricultural Sciences; HSAAS: Hengshui Academy of Agricultural Sciences; BDAAS: Baoding Academy of Agricultural Sciences; CZAFS: Cangzhou Academy of Agriculture and Forestry Sicences; YAIAS: Yanan Institute of Agricultural Sciences; SXAAS: Shanxi Academy of Agricultural Sciences; GAAAS: Gansu Academy of Agricultural Sciences; HLJAAS: Heilongjiang Academy of Agricultural Sciences; JLAAS: Jilin Academy of Agricultural Sciences; CFAAS: Chifeng Academy of Agricultural Sciences.

    發(fā)芽率鹽害率(%) = (對(duì)照發(fā)芽率-處理發(fā)芽率)/對(duì)照發(fā)芽率′100%;

    根長(zhǎng)鹽害率(%) = (對(duì)照根長(zhǎng)-處理根長(zhǎng))/對(duì)照根長(zhǎng)′100%;

    芽長(zhǎng)鹽害率(%) = (對(duì)照芽長(zhǎng)-處理芽長(zhǎng))/對(duì)照芽長(zhǎng)′100%;

    根鮮重鹽害率(%) = (對(duì)照根鮮重-處理根鮮重)/對(duì)照根鮮重′100%;

    芽鮮重鹽害率(%) = (對(duì)照芽鮮重-處理芽鮮重)/對(duì)照芽鮮重′100%;

    根冠比鹽害率(%) = (對(duì)照根冠比-處理根冠比)/對(duì)照根冠比′100%;

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

    采用Microsoft Excel 2007軟件處理數(shù)據(jù), 采用SPSS 18.0進(jìn)行方差分析、主成分分析、相關(guān)性分析和聚類(lèi)分析, 采用SPSS 18.0和R語(yǔ)言作圖。采用隸屬函數(shù)綜合評(píng)價(jià)性狀指標(biāo), 按以下公式進(jìn)行計(jì)算:

    (X) = (Xmin)/(maxmin)= 1, 2, 3, …,(1)

    式中,(X)表示第個(gè)綜合指標(biāo)的隸屬函數(shù)值,X表示第個(gè)綜合指標(biāo)值,max表示第個(gè)綜合指標(biāo)的最大值,min表示第個(gè)綜合指標(biāo)的最小值。

    W=P/∑=1P= 1, 2, 3, …,(2)

    式中,W表示第個(gè)綜合指標(biāo)在所有綜合指標(biāo)中權(quán)重,P為品種第個(gè)綜合指標(biāo)的貢獻(xiàn)率。

    = ∑=1[(X)′W]= 1, 2, 3, …,(3)

    式中,表示各品種耐鹽堿能力的綜合評(píng)價(jià)值。

    2 結(jié)果與分析

    2.1 鹽堿脅迫對(duì)谷子萌發(fā)各性狀的影響

    谷子萌發(fā)期鹽堿脅迫下, 除根冠比外, 各品種的發(fā)芽勢(shì)、發(fā)芽率、根長(zhǎng)、芽長(zhǎng)、根鮮重和芽鮮重的相對(duì)值均小于1, 表明除根冠比外, 各性狀值在鹽堿脅迫下均降低(表2)。各性狀的相對(duì)值以相對(duì)根長(zhǎng)最小, 相對(duì)發(fā)芽勢(shì)和相對(duì)發(fā)芽率較大。不同性狀間以相對(duì)根長(zhǎng)的變異最大, 變幅為8.7%~56.4%, 變異系數(shù)為45.8%; 相對(duì)發(fā)芽勢(shì)的變異最低, 變幅為53.3%~97.9%, 變異系數(shù)為12.2%。在鹽堿脅迫下各性狀的鹽害率與相對(duì)值表現(xiàn)相反的趨勢(shì), 以根長(zhǎng)的鹽害率最大, 均值為77.6%, 發(fā)芽勢(shì)的鹽害率最低, 均值為14.3%。各性狀的鹽害率在品種間的變異以根冠比鹽害率的變異最大, 變幅為-88.9%~25.4%, 變異系數(shù)為-908.4%; 以根長(zhǎng)鹽害率的變異最小, 變幅為43.6%~91.3%, 變異系數(shù)為13.2%。

    表2 不同谷子品種鹽堿脅迫下的性狀指標(biāo)

    RGP: relative germination potential; RGR: relative germination rate; RRL: relative root length; RSL: relative shoot length; RRFW: relative root fresh weight; RSFW: relative shoot fresh weight; RRSR: relative root to shoot rate; SIRGP: salt injury rate of germination potential; SIRGR: salt injury rate of germination rate; SIRRL: salt injury rate of root length; SIRSL: salt injury rate of shoot length; SIRRFW: salt injury rate of root fresh weight; SIRSFW: salt injury rate of shoot fresh weight; SIRRS: salt injury rate of root to shoot; CV: coefficient of variation.

    2.2 鹽堿脅迫下谷子各性狀指標(biāo)的頻率分布

    由圖1 (a, b)可知, 相對(duì)發(fā)芽勢(shì)和相對(duì)發(fā)芽率的分布范圍分別為50%~100%和45%~100%, 兩性狀均主要分布在80%~90%和90%~100%, 頻數(shù)分別為15、25和11、2, 分布頻率分別為28.3%、47.2%和20.8%、49.1%。不同谷子品種相對(duì)根長(zhǎng)和相對(duì)芽長(zhǎng)的分布范圍分別為5.0%~60.0%和35.0%~100.0% (圖2-a, b), 相對(duì)根長(zhǎng)主要集中分布在10%~20%和20%~30%, 分布頻數(shù)和頻率分別為27、16和50.9%、30.2%; 相對(duì)芽長(zhǎng)主要集中分布在70%~80%和80%~90%, 分布頻數(shù)和頻率分別為18、13和34.0%、24.5%。相對(duì)根鮮重和相對(duì)芽鮮重的分布范圍為50.0%~100.0%和40.0%~100.0% (圖2-c, d), 相對(duì)根鮮重主要集中分布在70.0%~80.0%和80.0%~90.0%, 分布頻數(shù)和頻率分別為19、16和35.8%、35.2%; 相對(duì)芽鮮重主要集中分布在80.0%~90.0%和90.0%~100.0%, 分布頻數(shù)和頻率分別為17、12和32.1%、22.6%。

    圖1 不同谷子品種相對(duì)發(fā)芽性狀頻率分布圖

    a: 相對(duì)發(fā)芽勢(shì)頻率分布; b: 相對(duì)發(fā)芽率頻率分布。

    a: frequency distribution diagram of relative germination potential; b: frequency distribution diagram of relative germination rate.

    圖2 不同谷子品種相對(duì)性狀頻率分布圖

    a: 相對(duì)根長(zhǎng)頻率分布; b: 相對(duì)芽長(zhǎng)頻率分布; c: 相對(duì)根鮮重頻率分布; d: 相對(duì)芽鮮重頻率分布。

    a: frequency distribution diagram of relative root length; b: frequency distribution diagram of relative shoot length; c: frequency distribution diagram of relative root fresh weight; d: frequency distribution diagram of relative shoot fresh weight.

    2.3 鹽堿脅迫下各性狀的相關(guān)性分析

    由表3可知, 相對(duì)發(fā)芽勢(shì)與相對(duì)發(fā)芽率呈極顯著正相關(guān), 相關(guān)系數(shù)為0.736 (<0.01)。相對(duì)發(fā)芽勢(shì)和相對(duì)發(fā)芽率與相對(duì)芽長(zhǎng)、相對(duì)芽鮮重和相對(duì)根鮮重均呈顯著或極顯著正相關(guān), 相關(guān)系數(shù)分別為0.281 (<0.05)和0.409 (<0.01)、0.287 (<0.05)和0.470 (<0.01)、0.430 (<0.01)和0.590 (<0.01), 與相對(duì)根冠比呈顯著和極顯著負(fù)相關(guān), 相關(guān)系數(shù)分別為-0.320 (<0.05)和-0.361 (<0.01); 相對(duì)根長(zhǎng)僅與相對(duì)發(fā)芽率呈顯著正相關(guān), 相關(guān)系數(shù)為0.347 (<0.05),與相對(duì)發(fā)芽勢(shì)相關(guān)性不顯著。相對(duì)芽長(zhǎng)與相對(duì)根長(zhǎng)和相對(duì)芽鮮重與相對(duì)根鮮重間均呈極顯著正相關(guān), 相關(guān)系數(shù)分別為0.451和0.458; 相對(duì)芽長(zhǎng)與相對(duì)芽鮮重間呈極限著正相關(guān), 相關(guān)系數(shù)為0.593, 而相對(duì)根長(zhǎng)與相對(duì)根鮮重間相關(guān)性不顯著。相對(duì)根冠比與相對(duì)根長(zhǎng)、相對(duì)芽長(zhǎng)和相對(duì)芽鮮重均呈顯著或極顯著負(fù)相關(guān), 相關(guān)系數(shù)分別為-0.290 (<0.01)、-0.475 (<0.05)和-0.744 (<0.05), 與相對(duì)根鮮重呈正相關(guān),但相關(guān)性不顯著。各性狀鹽害率的相關(guān)性與各性狀相對(duì)值間相關(guān)性表現(xiàn)出相反的趨勢(shì), 且兩者之間呈極顯著負(fù)相關(guān)。

    2.4 鹽堿脅迫下谷子各性狀的主成分分析

    主成分分析結(jié)果顯示(表4), 第I主成分的貢獻(xiàn)率為48.4%, 主要與發(fā)芽率、芽長(zhǎng)和芽鮮重的相對(duì)值和鹽害率有關(guān), 主要反映谷子的發(fā)芽率和地上部的鹽害性狀。第II主成分的貢獻(xiàn)率為18.6%, 主要與根鮮重的相對(duì)值和鹽害率有關(guān), 主要反映根鮮重的鹽害性狀。第III主成分的貢獻(xiàn)率為13.8%, 主要與發(fā)芽勢(shì)有關(guān), 主要反映谷子的萌發(fā)性狀。第IV主成分的貢獻(xiàn)率為9.5%, 主要與根長(zhǎng)的相對(duì)值和鹽害率有關(guān), 主要反映根長(zhǎng)的鹽害性狀。4個(gè)主成分的累積貢獻(xiàn)率為90.4%, 能有效地反應(yīng)數(shù)據(jù)的變化趨勢(shì), 符合主成分的分析要求。

    2.5 不同谷子芽期耐鹽堿的綜合評(píng)價(jià)

    根據(jù)主成分的貢獻(xiàn)率, 利用公式(2)計(jì)算出4個(gè)主成分的權(quán)重分別為0.536、0.206、0.153和0.105。Z1~Z4為各性狀指標(biāo)經(jīng)主成分分析后的4個(gè)主成分的綜合得分值, 利用主成分的綜合性狀得分值通過(guò)公式(1)計(jì)算出綜合性狀指標(biāo)的隸屬函數(shù)值, 根據(jù)各主成分的權(quán)重和綜合性狀的隸屬函數(shù)值利用公式(3)計(jì)算出各品種的綜合得分值并排序(表5)。根據(jù)各品種的綜合得分值進(jìn)行聚類(lèi)分析(圖3), 劃分為6類(lèi)不同類(lèi)型的品種: 強(qiáng)耐鹽堿型品種2個(gè)(編號(hào)2、 10), 耐鹽堿型品種16個(gè)(編號(hào)1、4、5、6、7、9、11、12、14、19、25、27、40、41、50、51), 中間型品種17個(gè)(編號(hào)3、8、15、16、17、18、21、24、26、28、29、33、34、39、44、48、53), 鹽堿敏感型品種6個(gè)(編號(hào)20、37、42、43、47、49), 不耐鹽堿型品種9個(gè)(編號(hào)13、22、23、30、35、36、45、46、52)和極不耐鹽堿型品種3個(gè)(編號(hào)31、32、38)。

    以綜合得分值D值為因變量, 各單項(xiàng)性狀指標(biāo)為自變量進(jìn)行逐步回歸分析, 得到回歸分析方程¢=0.298+0.0372+0.1443+0.0186+0.2097-0.1839+0.11511-0.20112+0.11213-0.10114+0.28415,方程的決定系數(shù)為2=1.000,<0.01。對(duì)綜合得分值和回歸值進(jìn)行檢測(cè), 回歸值和觀測(cè)值的估計(jì)精度均在99%以上(表6), 證明該方程對(duì)谷子的耐鹽堿性的評(píng)價(jià)可靠, 可用于耐鹽堿性的評(píng)價(jià)。在14個(gè)性狀指標(biāo)中, 該方程包含的10個(gè)性狀指標(biāo)中的相對(duì)發(fā)芽率(3)、根長(zhǎng)鹽害率(11)、芽長(zhǎng)鹽害率(12)和根冠比鹽害率(13) 4個(gè)性狀指標(biāo)均達(dá)顯著性的水平(<0.01), 可以作為谷子芽期耐鹽堿的評(píng)價(jià)指標(biāo)。

    3 討論

    3.1 鹽堿脅迫下各性狀的變化及頻率分布

    實(shí)際生產(chǎn)中, 鹽堿地存在鹽堿共同脅迫, 通過(guò)單一中性鹽[17,21]或堿性鹽[20]篩選的品種不能滿(mǎn)足實(shí)際生產(chǎn)的需求, 需要在混合鹽堿脅迫下進(jìn)行谷子耐鹽堿品種的篩選。本研究表明, 在混合鹽堿脅迫下谷子發(fā)芽勢(shì)、發(fā)芽率、根長(zhǎng)、芽長(zhǎng)、根鮮重和芽鮮重均降低, 但存在品種間的差異, 這與前人研究結(jié)果基本一致[24]。另有研究指出鹽堿脅迫下地下部根受抑制的程度大于地上部[8], 這與本研究谷子相對(duì)根長(zhǎng)顯著小于相對(duì)芽長(zhǎng)的研究結(jié)果一致。同時(shí)本研究結(jié)果還表明不同谷子品種相對(duì)發(fā)芽勢(shì)、相對(duì)發(fā)芽率和相對(duì)根長(zhǎng)的分布比較集中, 相對(duì)芽長(zhǎng)、相對(duì)根鮮重和相對(duì)芽鮮重的分布比較分散, 且各性狀相對(duì)值超過(guò)80%的品種數(shù)分別為0、16、40、37、13和29, 表明鹽堿脅迫對(duì)谷子的根長(zhǎng)影響最大, 根鮮重、芽長(zhǎng)和芽鮮重次之, 對(duì)發(fā)芽勢(shì)和發(fā)芽率的影響最小。

    3.2 谷子鹽堿脅迫下的相關(guān)性、主成分分析和綜合評(píng)價(jià)

    韓飛等[21]研究表明, 在鹽脅迫下谷子相對(duì)芽長(zhǎng)和相對(duì)根長(zhǎng)呈顯著正相關(guān), 鹽害率與相對(duì)芽長(zhǎng)和根長(zhǎng)均呈極顯著負(fù)相關(guān)。田伯紅等[17]研究表明, 在鹽脅迫下谷子相對(duì)根長(zhǎng)和相對(duì)芽長(zhǎng)呈顯著正相關(guān), 但與鹽害率無(wú)顯著相關(guān)性。本研究表明, 在混合鹽堿脅迫下, 谷子的相對(duì)根長(zhǎng)與相對(duì)芽長(zhǎng)及相對(duì)根鮮重與相對(duì)芽鮮重均呈顯著正相關(guān), 這與前人研究結(jié)果一致; 且相對(duì)發(fā)芽率與相對(duì)根長(zhǎng)、相對(duì)芽長(zhǎng)、相對(duì)根鮮重和相對(duì)芽鮮重均呈顯著正相關(guān), 表明鹽堿脅迫對(duì)谷子地上部和地下部的影響存在協(xié)同性, 且均對(duì)發(fā)芽率有著顯著的作用。這與于菘等[25]認(rèn)為在鹽堿脅迫下綠豆各性狀均存在協(xié)同性的研究一致。

    表4 主成分的特征根值、貢獻(xiàn)率和載荷矩陣

    縮寫(xiě)同表2。Abbreviations are the same as those given in Table 2.

    作物的耐鹽性是一個(gè)綜合性狀, 通過(guò)單一性狀或少量性狀的評(píng)價(jià)不能真實(shí)有效地反映作物的耐鹽性, 需要進(jìn)行多項(xiàng)指標(biāo)的綜合評(píng)價(jià)[26-27]。主成分分析是一項(xiàng)重要的降維評(píng)價(jià)方法, 可將多個(gè)因素合成為一個(gè)或幾個(gè)因素(主成分)進(jìn)行分析, 在作物的性狀綜合分析中[28-30]得到了廣泛的應(yīng)用。隸屬函數(shù)是一項(xiàng)重要的等級(jí)評(píng)價(jià)方法, 已在作物的耐鹽性評(píng)價(jià)[5,21]中進(jìn)行應(yīng)用。隸屬函數(shù)側(cè)重于評(píng)價(jià)結(jié)果, 對(duì)每個(gè)因素進(jìn)行等級(jí)評(píng)價(jià), 但缺乏考慮對(duì)每個(gè)因素的貢獻(xiàn)率, 所以需要兩者結(jié)合進(jìn)行綜合的分析和評(píng)價(jià)。本研究在鹽堿脅迫下通過(guò)主成分分析選取了特征根和貢獻(xiàn)率較大的4個(gè)主成分, 累計(jì)貢獻(xiàn)率為90.35%, 符合主成分分析的要求。同時(shí)以4個(gè)主成分的得分值及其權(quán)重通過(guò)隸屬函數(shù)獲得綜合得分值, 并進(jìn)行聚類(lèi)分析, 獲得6種不同耐鹽堿類(lèi)型的品種。關(guān)于作物耐鹽性指標(biāo)的鑒定, 前人多利用主成分分析及相關(guān)性進(jìn)行確定[21,27], 但僅獲得綜合性狀指標(biāo), 因此需要通過(guò)建立回歸方程進(jìn)行回歸分析從而獲得單項(xiàng)鑒定指標(biāo)[26,31]。本研究以綜合得分值和各單項(xiàng)指標(biāo)值進(jìn)行回歸分析, 建立了回歸方程, 回歸精度均在99.0%以上, 該方程具有有效性(<0.01)。同時(shí)利用回歸方程確定相對(duì)發(fā)芽率(<0.01)、根長(zhǎng)鹽害率(<0.01)、芽長(zhǎng)鹽害率(<0.01)和根冠比鹽害率(<0.01)可以作為谷子芽期耐鹽堿性的關(guān)鍵鑒定指標(biāo)。

    3.3 谷子耐鹽堿性的預(yù)期及研究方向

    本研究確定了谷子芽期耐鹽堿的鑒定方法和指標(biāo), 過(guò)表達(dá)谷子耐鹽性基因?qū)λ久劝l(fā)期和成株期的耐鹽性是存在差異的[11], 因此谷子芽期和幼苗期的耐鹽性和成株期的耐鹽性是否存在一致性, 仍需進(jìn)一步研究。前人在水稻上建立了耐鹽性的回歸方程, 并指出地上部的生物量、葉片光合速率和蒸騰速率可以作為水稻成株期的耐鹽性的鑒定指標(biāo)[32], 且鑒定指標(biāo)隨著環(huán)境而變化[33]。本研究團(tuán)隊(duì)前期在大田條件下進(jìn)行了谷子成株期耐鹽堿性的篩選和鑒定, 指出地上部生物量、單穗重和地上部含水量可以作為耐鹽堿性的鑒定指標(biāo)[34], 與在水稻上的研究結(jié)果相似。但由于受田間試驗(yàn)品種的數(shù)量的限制, 仍需進(jìn)一步選擇不同代表類(lèi)型的品種在田間進(jìn)行谷子成株期耐鹽堿性的鑒定, 從而獲得更有代表型的鑒定指標(biāo)和方法, 同時(shí)應(yīng)明確成株期是否與芽期或苗期耐鹽堿性存在一致性。近年來(lái), 谷子在分子水平上的研究已取得重要進(jìn)展, 在水稻和擬南芥中通過(guò)轉(zhuǎn)錄或過(guò)表達(dá)谷子的耐鹽性基因[11,14,16], 均能提高轉(zhuǎn)錄作物的耐鹽性。隨著谷子轉(zhuǎn)化水平的提高, 應(yīng)進(jìn)一步加強(qiáng)在谷子中轉(zhuǎn)化和應(yīng)用研究, 在分子水平上提高谷子的耐鹽性, 為谷子抗鹽性的育種提供材料和方法。

    圖3 不同谷子品種D值聚類(lèi)分析圖

    Fig. 3 Cluster analysis of D-values of different millet cultivars

    表6 回歸方程的精度分析

    (續(xù)表6)

    4 結(jié)論

    對(duì)生產(chǎn)中主要推廣的53個(gè)谷子品種進(jìn)行芽期耐鹽堿綜合鑒定, 各性狀均受到不同程度的抑制, 以根長(zhǎng)受到的抑制性最大, 相對(duì)發(fā)芽勢(shì)的影響最小; 相對(duì)發(fā)芽率、根長(zhǎng)鹽害率、芽長(zhǎng)鹽害率和根冠比鹽害率可以作為谷子芽期耐鹽堿鑒定的評(píng)價(jià)指標(biāo)。

    [1] Liang W X, Ma X L, Wang P, Liu L Y. Plant salt-tolerance mechanism: a review., 2018, 495: 286–291.

    [2] Sun H, Meng M H, Yan Z H, Lin Z X, Nie X H, Yang X Y. Genome-wide association mapping of stress tolerance traits in cotton., 2019, 7: 77–88.

    [3] Wu H H. Plant salt tolerance and Na+sensing and tranport., 2018, 6: 215–225.

    [4] Shrivastava P, Kumar R. Soil salinity: a serious environmental issue and plant growth promoting bacteria as one of the tools for its alleviation., 2015, 22: 123–131.

    [5] 李豐先, 周宇飛, 王藝陶, 孫璐, 白薇, 閆彤, 許文娟, 黃瑞冬. 高粱品種萌發(fā)期耐堿性篩選與綜合鑒定. 中國(guó)農(nóng)業(yè)科學(xué), 2013, 46: 25–34. Li F X, Zhou Y F, Wang Y T, Sun L, Bai W, Yan T, Xu W J, Huang R D. Screening and identification of sorghum cultivars for alkali tolerance during germination., 2013, 46: 25–34 (in Chinese with English abstract).

    [6] 景宇鵬, 連海飛, 李煥春, 史培, 杜超, 劉梅, 常新娟. 玉米耐鹽堿能力及評(píng)價(jià)指標(biāo)篩選研究. 中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào), 2018, 20(11): 94–104. Jing Y P, Lian H F, Li H C, Shi P, Du C, Liu M, Chang X J. Study on saline-alkaline tolerance ability and evaluation index selection for corn., 2018, 20(11): 94–104 (in Chinese with English abstract).

    [7] 信彩云, 周學(xué)標(biāo), 劉奇華, 謝先芝, 趙慶雷, 王瑜, 朱其松, 陳博聰, 王佳, 馬佳清. 黃河三角洲鹽堿土條件下水稻苗期耐鹽指標(biāo)篩選. 山東農(nóng)業(yè)科學(xué), 2017, 49 (1): 23–28. Xin C Y, Zhou X B, Liu Q H, Xie X Z, Zhao Q L, Wang Y, Zhu Q S, Chen B C, Wang J, Ma X Q. Screening of salt tolerant indexes for rice seedlings under salinity conditions of the Yellow River delta., 2017, 49(1): 23–28 (in Chinese with English abstract).

    [8] 于瑩, 張樹(shù)權(quán), 郭永利, 姚玉波, 周菲, 趙東升, 劉繼忠, 李文華. 31份玉米自交系萌發(fā)期耐鹽堿性綜合評(píng)價(jià). 東北農(nóng)業(yè)大學(xué)學(xué)報(bào), 2018, 49(9): 9–19. Yu Y, Zhang S Q, Guo Y L, Yao Y B, Zhou F, Zhao D S, Liu J Z, Li W H. Comprehensive evaluation of saline-alkaline tolerance of 31 maize inbred lines at germination stage., 2018, 49(9): 9–19 (in Chinese with English abstract).

    [9] 高春華, 朱金英, 張華文, 田藝心, 高鳳菊. 38 個(gè)粒用高粱品種芽期耐鹽性的綜合鑒定及評(píng)價(jià). 核農(nóng)學(xué)報(bào), 2019, 33: 1841–1855. Gao C H, Zhu J Y, Zhang H W, Tian Y X, Gao F J. Comprehensive identification and evaluation of 38 grain sorghum cultivars for tolerance during germination., 2019, 33: 1841–1855 (in Chinese with English abstract).

    [10] 曹帥, 杜仲陽(yáng), 向殿軍, 李志剛, 孫賀祥, 宮曉旭, 劉鵬. 18 份大豆品種耐鹽堿性篩選與綜合鑒定. 大豆科學(xué), 2019, 38: 344–352. Cao S, Du Z Y, Xiang D J, Li Z G, Sun H X, Gong X Y, Liu P. Salt and alkaline tolerance screening and comprehensive identification of eighteen soybean varieties., 2019, 38: 44–52 (in Chinese with English abstract).

    [11] 寧蕾, 王曙光, 琚鵬舉, 柏星軒, 葛林豪, 齊欣, 姜奇彥, 孫現(xiàn)軍, 陳明, 孫黛珍. 過(guò)表達(dá)谷子對(duì)水稻耐鹽性的影響. 中國(guó)農(nóng)業(yè)科學(xué), 2018, 51: 23–34. Ning L, Wang S G, Ju P J, Bai X X, Ge L H, Qi X, Jiang Q Y, Sun X J, Chen M, Sun D Z. Rice over expression of milletgene increases salt tolerance., 2013, 51: 23–34 (in Chinese with English abstract).

    [12] 陳二影, 楊延兵, 程炳文, 秦嶺, 張華文, 劉賓, 王海蓮, 陳桂玲, 管延安. 不同夏谷品種的產(chǎn)量與氮肥利用效率. 中國(guó)土壤與肥料, 2015, (2): 101–105. Chen E Y, Yang Y B, Cheng B W, Qin L, Zhang H W, Liu B, Wang H L, Chen G L, Guan Y A. Yield and nitrogen use efficiency in summer foxtail millet., 2015, (2): 101–105 (in Chinese with English abstract).

    [13] 陳二影, 秦嶺, 程炳文, 楊延兵, 張華文, 劉賓, 王海蓮, 陳桂玲, 管延安. 夏谷氮、磷、鉀肥的效應(yīng)研究. 山東農(nóng)業(yè)科學(xué), 2015, 47(1): 67–71. Chen E Y, Qin L, Cheng B W, Yang Y B, Zhang H W, Liu B, Wang H L, Chen G L, Guan Y A. Application effect of N, P, K fertilizers on summer foxtail millet., 2017, 47(1): 67–71 (in Chinese with English abstract).

    [14] 馮露, 鐘理, 陳丹丹, 馬有志, 徐兆師, 李連城, 周永斌, 陳明, 張小紅. 過(guò)表達(dá)谷子液泡H+-ATPase E亞基基因在擬南芥中的耐鹽性. 作物學(xué)報(bào), 2015, 41: 66–75. Feng L, Zhong L, Chen D D, Ma Y Z, Xu Z S, Li L C, Zhou Y B, Chen M, Zhang X H. Over expression of vacuole H+-ATPase E subunit genefrom foxtail millet enhances salt resistance in transgenic., 2015, 41: 66–75 (in Chinese with English abstract).

    [15] 黃鎖, 胡利芹, 徐東北, 李微微, 徐兆師, 李連城, 周永斌, 刁現(xiàn)民, 賈冠清, 馬有志, 陳明. 谷子轉(zhuǎn)錄因子SiNF-YA5通過(guò)ABA非依賴(lài)途徑提高轉(zhuǎn)基因擬南芥耐鹽性. 作物學(xué)報(bào), 2016, 42: 1787–1797. Huang S, Hu L Q, Xu D B, Li W W, Xu Z S, Li L C, Zhou Y B, Diao X M, Jia G Q, Ma Y Z, Chen M. Transcription factor SiNF-YA5 from foxtail millet () conferred tolerance to high-salt stress through ABA-independent pathway in transgenic., 2016, 42: 1787–1797 (in Chinese with English abstract).

    [16] 秦玉海, 張小紅, 馮露, 李微微, 徐兆師, 李連城, 周永斌, 馬有志, 刁現(xiàn)民, 賈冠清, 陳明, 閔東紅. 谷子轉(zhuǎn)錄因子基因在擬南芥中對(duì)高鹽和ABA的響應(yīng). 中國(guó)農(nóng)業(yè)科學(xué), 2016, 49: 3276–3286. Qin Y H, Zhang X H, Feng L, Li H H, Xu Z S, Li L C, Zhou Y B, Ma Y Z, Diao X M, Jia G Q, Chen M, Min D H. Response of millet transcription factor geneto high salt and ABA treatment in transgenic., 2016, 49: 3276–3286 (in Chinese with English abstract).

    [17] 田伯紅, 王素英, 李雅靜, 王建廣, 張立新, 梁鳳芹, 翟玉柱, 劉金榮. 谷子地方品種發(fā)芽期和苗期對(duì)NaCl脅迫的反應(yīng)和耐鹽品種篩選. 作物學(xué)報(bào), 2008, 34: 164–168. Tian B H, Wang S Y, Li Y J, Wang J G, Zhang L X, Liang F Q, Zhai G Z, Liu J R. Response to sodium chloride stress at germination and seedling and identification of salinity tolerant genotypes in foxtail millet landraces originated from China., 2008, 34: 164–168 (in Chinese with English abstract).

    [18] 崔興國(guó), 時(shí)麗冉. 衡水地區(qū)14份夏谷品種種子萌發(fā)期耐鹽性研究. 作物雜志, 2011, (4): 117–119. Cui X G, Shi L R. Evaluation of salt tolerance of 14 foxtail millet summer varieties in Hengshui., 2011, (4): 117–119 (in Chinese with English abstract).

    [19] 崔興國(guó), 時(shí)麗冉. 鹽脅迫對(duì)不同品種谷子萌發(fā)及幼苗生長(zhǎng)的影響. 黑龍江農(nóng)業(yè)科學(xué), 2011, (6): 14–16. Cui X G, Shi L R. Effect of salt stress on different seed germination and seedling growth of(L.)., 2011, (6): 14–16 (in Chinese with English abstract).

    [20] 張永芳, 宋喜娥, 王潤(rùn)梅, 湛潤(rùn)生, 趙治海, 范光宇. Na2CO3脅迫對(duì)谷子種子萌發(fā)的影響. 種子, 2015,34(11): 102–105. Zhang Y F, Song X E, Wang R M, Zhan R S, Zhao Z H, Fan G Y. Effect of Na2CO3stress on different seed germination of millet., 2015, 34(11): 102–105 (in Chinese with English abstract).

    [21] 韓飛, 諸葛玉平, 婁燕宏, 王會(huì), 張乃丹, 何偉, 晁贏. 63份谷子種質(zhì)的耐鹽綜合評(píng)價(jià)及耐鹽品種篩選. 植物遺傳資源學(xué)報(bào), 2018, 19: 103–111. Han F, Zhu-Ge Y P, Lou Y H, Wang H, Zhang N D, He W, Chao Y. Evaluation of salt tolerance and screening for salt tolerant accessions of 63 foxtail millet germplasm., 2018, 19: 103–111 (in Chinese with English abstract).

    [22] 張鵬銳, 李旭霖, 崔德杰, 胡景田. 濱海重鹽堿地不同土地利用方式的水鹽特征. 水土保持學(xué)報(bào), 2015, 29(2): 119–123. Zhang P R, Li X L, Cui D J, Hu J T. Characteristics of water and salt under different land use in heavy coastal saline-alkaline land., 2015, 29(2): 119–123 (in Chinese with English abstract).

    [23] 高彥花. 渤海灣鹽堿地土壤水鹽動(dòng)態(tài)及耐鹽植物改良鹽漬土研究. 中國(guó)林業(yè)科學(xué)研究院博士學(xué)位論文, 北京, 2011. pp 18–20.Gao Y H. Study on the Dynamics of Water and Salt in Saline Soil and Ameliortion of Saline Soil by Salt-resistance Plants in Bohai Bay. PhD Dissertation of Graduate School of Chinese Academy of Forestry, Beijing, China, 2011. pp 18–20 (in Chinese with English abstract).

    [24] 郭瑞鋒, 張永福, 任月梅, 楊忠. 混合鹽堿脅迫對(duì)谷子萌發(fā)、幼芽生長(zhǎng)的影響及耐鹽堿品種篩選. 作物雜志, 2017, (4): 63–66. Guo R F, Zhang Y F, Ren Y M, Yang Z. Effects of saline-alkali stress on millet germination and shoots growth and saline-alkali tolerance variety screening., 2017, (4): 63–66 (in Chinese with English abstract).

    [25] 于菘, 郭瀟瀟, 梁海蕓, 付鸞鴻, 史京京, 張翼飛, 闖磊. 不同基因型綠豆萌發(fā)期耐鹽堿性分析及其鑒定指標(biāo)的篩選. 植物生理學(xué)報(bào), 2017, 53: 1629–1639. Yu S, Guo X X, Liang H Y, Fu L H, Shi J J, Zhang Y F, Chuang L. Analysis of saline-alkaline tolerance and screening of identification indicators at the germination stage among different mung bean genotypes., 2017, 53: 1629–1639 (in Chinese with English abstract).

    [26] 段文學(xué), 張海燕, 解備濤, 汪寶卿, 張立明. 甘薯苗期耐鹽性鑒定及其指標(biāo)篩選. 作物學(xué)報(bào), 2018, 44: 137–147. Duan W X, Zhang H Y, Xie B T, Wang B Q, Zhang L M. Identification of salt tolerance and screening for its indicators in sweet potato varieties during seedling stage., 2018, 44: 137–147 (in Chinese with English abstract).

    [27] 苑乂川, 陳小雨, 李明明, 李萍, 賈亞濤, 韓淵懷, 邢國(guó)芳. 谷子苗期耐低磷種質(zhì)篩選及其根系保護(hù)酶系統(tǒng)對(duì)低磷脅迫的響應(yīng). 作物學(xué)報(bào), 2019, 45: 121–132. Fan Y C, Chen X Y, Li M M, Li P, Jia Y T, Han Y H, Xing G F. Screening of germplasm tolerant to low phosphorus of seedling stage and response of root protective enzymes to low phosphorus in foxtail millet., 2019, 45: 121–132 (in Chinese with English abstract).

    [28] 褚能明, 柯劍鴻, 袁亮. 不同鮮食甜糯玉米揮發(fā)性風(fēng)味物質(zhì)主成分分析. 核農(nóng)學(xué)報(bào), 2017, 31: 2175–2185. Chu N M, Ke J H, Yuan L. Principal components analysis for volatility of flavor compositions in different fresh sweet glutinous corn., 2017, 31: 2175–2185 (in Chinese with English abstract).

    [29] 田茂成, 鄧小華, 陸中山, 田峰, 陳治鋒, 張明發(fā), 張黎明. 基于灰色效果測(cè)度和主成分分析的湘西州煙葉物理特性綜合評(píng)價(jià). 核農(nóng)學(xué)報(bào), 2017, 31: 187–193. Tian M C, Deng X H, Lu Z S, Tian F, Chen Z F, Zhang M F, Zhang L M. Gray effect measure and principal component analysis-based comprehensive evaluation for physical properties of flue-cured tobacco leaves from Xiangxi area., 2017, 31: 187–193 (in Chinese with English abstract).

    [30] 孫東雷, 卞能飛, 陳志德, 邢興華, 徐澤俊, 齊玉軍, 王曉軍, 王幸. 花生萌發(fā)期耐鹽性綜合評(píng)價(jià)及耐鹽種質(zhì)篩選. 植物遺傳資源學(xué)報(bào), 2017, 18: 1079–1087. Sun D L, Bian N F, Chen Z D, Xing X H, Xu Z J, Qi Y J, Wang X J, Wang X. Comprehensive evaluation of salt tolerance and screening for salt tolerant accessions of peanut (L) at germination stage., 2017, 18: 1079–1087 (in Chinese with English abstract).

    [31] 戴海芳, 武輝, 阿曼古麗·買(mǎi)買(mǎi)提阿力, 王立紅, 麥麥提·阿皮孜, 張巨松. 不同基因型棉花苗期耐鹽性分析及其鑒定指標(biāo)篩選. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47: 1290–1300. Dai H F, Wu H, Amanguli Maimaitiali, Wang L H, Maimaiti Apizi, Zhang J S. Analysis of salt-tolerance and determination of salt-tolerant evaluation indicators in cotton seedlings of different genotypes., 2014, 47: 1290–1300 (in Chinese with English abstract).

    [32] Radanielson A M, Angeles O, Li T, Ismail A M, Gaydon D S. Describing the physiological responses of different rice genotypes to salts tress using sigmoid and piecewise linear functions., 2018, 220: 46–56.

    [33] Burman D, Maji B, Singh S, Mandal S, Sarangi S K, Bandyopadhyay, B K, Bal A R, Sharma D K, Krishnamurthy S L, Singh H N, Delosreyes A S, Villanueva D, Paris T, Singh U S, Haefele S M, Ismail A M. Participatory evaluation guides the development and selection of farmers’ preferred rice varieties for salt- and flood-affected coastal deltas of south and southeast Asia., 2018, 220: 67–77.

    [34] 陳二影, 秦嶺, 楊延兵, 黎飛飛, 王潤(rùn)豐, 張華文, 王海蓮, 劉賓, 孔清華, 管延安. 生產(chǎn)條件下谷子品種鹽堿耐性的差異及綜合評(píng)價(jià). 中國(guó)農(nóng)業(yè)科學(xué), 2019, 52: 4050–4065. Chen E Y, Qin L, Yang Y B, Li F F, Wang R F, Zhang H W, Wang H L, Liu B, Kong Q H, Guan Y A. Variation and comprehensive evaluation of salt and alkali tolerance of different foxtail millet cultivars under production conditions., 2019, 52: 4050–4065 (in Chinese with English abstract).

    Comprehensive identification and evaluation of foxtail millet for saline-alkaline tolerance during germination

    CHEN Er-Ying1, WANG Run-Feng1, QIN Ling1, YANG Yan-Bing1, LI Fei-Fei1, ZHANG Hua-Wen1, WANGHai-Lian1, LIU Bin1, KONG Qing-Hua2, and GUAN Yan-An1,2,*

    1Crop Research Institute, Shandong Academy of Agricultural Sciences / Featured Crops Engineering Laboratory of Shandong Province, Jinan 250100, Shandong, China;2Shandong Normal University, Jinan 250014, Shandong, China

    In the present study, the evaluation of 53 main foxtail millet cultivars was carried out under saline-alkaline mixed stress (100 mmol L-1, NaCl : NaHCO3= 4 : 1). The results showed that germination potential, germination rate, root length, shoot length, fresh root weight and fresh shoot weight of the 53 cultivars were inhibited, among which with root length most affected by the salt and alkaline condition. Exposed to such a stress condition, significant or extremely significant positive correlations were observed for relative germination potential and relative germination rate, relative root length and relative shoot length, and relative fresh root weight and relative fresh shoot weight. 14 traits were integrated into four principal components with a cumulative contribution rate of 90.4% through principal component analysis (PCA). Composite scores for saline-alkaline tolerance were calculated from membership function with scores of the four principal components. The 53 cultivars were assigned to six different groups of saline-alkaline resistance by using cluster analysis, including 2 highly resistant, 16 moderately resistant, 17 lower resistant, 6 sensitive, 9 susceptible and 3 extremely intolerant foxtail millet cultivars. Meanwhile, a regression equation,¢= 0.298 + 0.0372+ 0.1443+ 0.0186+ 0.2097-0.1839+ 0.11511-0.20112+ 0.11213-0.10114+ 0.28415, was established to estimate the tolerance of foxtail millet cultivars to saline-alkaline stress. Relative germination rate, salt injury rates for root length and shoot length, and root-shoot ratio could be regarded as the indicators of assessing the resistance of foxtail millet to saline-alkaline mixed stress.

    foxtail millet; saline-alkaline stress; principal component analysis; membership function; regression analysis.

    10.3724/SP.J.1006.2020.04064

    本研究由國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2019YFD1001703, 2019YFD1001700, 2019YFD1002703), 山東省自然科學(xué)基金項(xiàng)目(ZR2017YL010), 國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系項(xiàng)目(CARS-06-13.5-A19)和山東省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)科技創(chuàng)新工程項(xiàng)目(CXGC2018E01)資助。

    This study was supported by the National Key Research and Development Program of China (2019YFD1001703, 2019YFD1001700, 2019YFD1002703), the Natural Science Foundation of Shandong Province (ZR2017YL010), the China Agriculture Research System (CARS-06-13.5-A19), and the Agricultural Scientific and Technological Innovation Project of Shandong Academy of Agricultural Sciences (CXGC2018E01).

    管延安, E-mail: Yguan65@163.com, Tel: 0531-66658115

    E-mail: chenerying_001@163.com, Tel: 0531-66658115

    2020-03-12;

    2020-06-02;

    2020-07-01.

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

    猜你喜歡
    芽長(zhǎng)鹽堿耐鹽性
    鹽堿荒灘的“底色之變”
    金橋(2021年11期)2021-11-20 06:37:14
    鹽堿地區(qū)南美白對(duì)蝦、青蝦混養(yǎng)技術(shù)
    郁金香耐鹽性鑒定方法研究
    玉米種子芽大小對(duì)農(nóng)桿菌轉(zhuǎn)化的影響
    播種不同芽長(zhǎng)種薯芽塊對(duì)早春錯(cuò)期馬鈴薯生育進(jìn)程及產(chǎn)量的影響
    不同芽長(zhǎng)對(duì)馬鈴薯理化指標(biāo)的影響
    源于大麥小孢子突變體的苗期耐鹽性研究
    鹽堿脅迫對(duì)三種西藏牧草種子萌發(fā)及幼苗影響
    西藏科技(2016年8期)2016-09-26 09:00:55
    三個(gè)大豆品種萌發(fā)期和苗期的耐鹽性比較
    大豆發(fā)芽期和苗期耐鹽性的隸屬函數(shù)分析
    嫁个100分男人电影在线观看| 亚洲国产精品999在线| 国产精品98久久久久久宅男小说| 国产99白浆流出| 中文亚洲av片在线观看爽| 亚洲精品中文字幕在线视频| 少妇粗大呻吟视频| 人人妻,人人澡人人爽秒播| 国产精品久久久人人做人人爽| 久久人妻福利社区极品人妻图片| 在线看a的网站| 国产精品秋霞免费鲁丝片| 亚洲精品国产色婷婷电影| 国产在线精品亚洲第一网站| 欧美激情久久久久久爽电影 | 亚洲自拍偷在线| 亚洲av美国av| 成在线人永久免费视频| 18禁黄网站禁片午夜丰满| 操出白浆在线播放| 99久久国产精品久久久| 香蕉国产在线看| 亚洲人成电影观看| 91av网站免费观看| 香蕉丝袜av| 国产免费现黄频在线看| 欧美日韩中文字幕国产精品一区二区三区 | 在线天堂中文资源库| 黄色女人牲交| 亚洲aⅴ乱码一区二区在线播放 | 亚洲少妇的诱惑av| 亚洲精品美女久久久久99蜜臀| 啦啦啦 在线观看视频| 午夜免费激情av| 久久草成人影院| av国产精品久久久久影院| 日本黄色视频三级网站网址| 一a级毛片在线观看| 亚洲av熟女| 国产成人一区二区三区免费视频网站| 丰满人妻熟妇乱又伦精品不卡| 老司机在亚洲福利影院| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区在线臀色熟女 | 日韩精品青青久久久久久| 天堂动漫精品| 啪啪无遮挡十八禁网站| 精品久久久久久电影网| 精品国内亚洲2022精品成人| 999久久久国产精品视频| 在线永久观看黄色视频| 最近最新免费中文字幕在线| 色综合欧美亚洲国产小说| 在线观看免费高清a一片| 精品一区二区三区视频在线观看免费 | 国产av一区二区精品久久| 19禁男女啪啪无遮挡网站| 一级,二级,三级黄色视频| 久久草成人影院| 久久精品国产综合久久久| 久久久久久免费高清国产稀缺| 999久久久国产精品视频| 国产精品一区二区在线不卡| 身体一侧抽搐| 精品少妇一区二区三区视频日本电影| 黄色a级毛片大全视频| 性少妇av在线| 精品一区二区三区四区五区乱码| av福利片在线| 国产欧美日韩一区二区精品| tocl精华| 欧美日韩一级在线毛片| 怎么达到女性高潮| 欧美日韩亚洲高清精品| 男女床上黄色一级片免费看| 中国美女看黄片| 久久精品国产99精品国产亚洲性色 | 国产成人精品在线电影| 18美女黄网站色大片免费观看| 国产激情久久老熟女| 国产精品秋霞免费鲁丝片| 国产激情久久老熟女| 久热爱精品视频在线9| 久久久久久人人人人人| 久久精品人人爽人人爽视色| a在线观看视频网站| 日韩欧美在线二视频| 亚洲国产精品合色在线| 国产成人系列免费观看| 亚洲七黄色美女视频| 桃色一区二区三区在线观看| 国产成人精品久久二区二区免费| www.自偷自拍.com| 极品教师在线免费播放| 欧美黄色淫秽网站| 欧美乱色亚洲激情| 久久久国产欧美日韩av| 日韩三级视频一区二区三区| 国产极品粉嫩免费观看在线| 欧美最黄视频在线播放免费 | 欧美日韩中文字幕国产精品一区二区三区 | 在线观看免费视频网站a站| videosex国产| 国产97色在线日韩免费| 欧美性长视频在线观看| av欧美777| 91精品三级在线观看| 亚洲欧美一区二区三区黑人| 老司机午夜十八禁免费视频| 99riav亚洲国产免费| 男人舔女人的私密视频| 亚洲午夜理论影院| 亚洲情色 制服丝袜| 在线观看日韩欧美| 男人的好看免费观看在线视频 | 亚洲av日韩精品久久久久久密| 久久久久久大精品| 高清毛片免费观看视频网站 | 国产成人免费无遮挡视频| 国产99久久九九免费精品| 可以在线观看毛片的网站| 成人18禁在线播放| 亚洲 欧美一区二区三区| 制服人妻中文乱码| 欧美乱色亚洲激情| 亚洲三区欧美一区| av天堂在线播放| 久久精品国产清高在天天线| 黄频高清免费视频| 99精品久久久久人妻精品| 久久久久久久久免费视频了| 亚洲av五月六月丁香网| 日韩免费av在线播放| 两个人免费观看高清视频| 叶爱在线成人免费视频播放| 欧美色视频一区免费| 法律面前人人平等表现在哪些方面| 香蕉久久夜色| 久久久久亚洲av毛片大全| 成人三级黄色视频| 欧美老熟妇乱子伦牲交| 久久久久久久午夜电影 | 9191精品国产免费久久| 久久香蕉精品热| 999精品在线视频| 波多野结衣高清无吗| 男女午夜视频在线观看| 亚洲国产毛片av蜜桃av| 欧美黄色片欧美黄色片| 国产有黄有色有爽视频| 国产精品美女特级片免费视频播放器 | 欧美中文日本在线观看视频| 中文字幕色久视频| 桃色一区二区三区在线观看| 天堂影院成人在线观看| 人妻丰满熟妇av一区二区三区| 在线十欧美十亚洲十日本专区| 亚洲五月色婷婷综合| 日本五十路高清| 国产xxxxx性猛交| 亚洲视频免费观看视频| 国产精品香港三级国产av潘金莲| 免费在线观看完整版高清| 欧美日韩中文字幕国产精品一区二区三区 | 国产熟女午夜一区二区三区| 国产成人一区二区三区免费视频网站| 精品一区二区三区av网在线观看| 丁香六月欧美| 久久精品国产综合久久久| 人妻丰满熟妇av一区二区三区| 成人亚洲精品av一区二区 | 国产欧美日韩一区二区三区在线| 午夜影院日韩av| 久久中文看片网| 叶爱在线成人免费视频播放| 久9热在线精品视频| 性欧美人与动物交配| 在线观看免费午夜福利视频| 精品午夜福利视频在线观看一区| 国产成人免费无遮挡视频| 三上悠亚av全集在线观看| 80岁老熟妇乱子伦牲交| 国产男靠女视频免费网站| 日韩精品青青久久久久久| 国产91精品成人一区二区三区| 啪啪无遮挡十八禁网站| av片东京热男人的天堂| xxx96com| 日本黄色视频三级网站网址| 亚洲精品久久成人aⅴ小说| av在线天堂中文字幕 | 99在线视频只有这里精品首页| 日本三级黄在线观看| 亚洲自偷自拍图片 自拍| 国产在线观看jvid| 日韩视频一区二区在线观看| 波多野结衣一区麻豆| 久久青草综合色| 激情在线观看视频在线高清| 男人舔女人下体高潮全视频| 国产精品av久久久久免费| 青草久久国产| 少妇的丰满在线观看| 多毛熟女@视频| 国产亚洲精品综合一区在线观看 | 精品国产乱子伦一区二区三区| 久久久久久久久中文| 变态另类成人亚洲欧美熟女 | xxx96com| 丝袜美足系列| 国产在线观看jvid| 桃红色精品国产亚洲av| 亚洲一区中文字幕在线| 欧美激情 高清一区二区三区| 一区在线观看完整版| 一级毛片高清免费大全| 亚洲成人免费av在线播放| 夜夜爽天天搞| 美女福利国产在线| av在线播放免费不卡| 丝袜人妻中文字幕| 女性被躁到高潮视频| a级片在线免费高清观看视频| 精品久久久久久电影网| 少妇 在线观看| 男女床上黄色一级片免费看| 国产区一区二久久| 激情视频va一区二区三区| videosex国产| 亚洲欧洲精品一区二区精品久久久| 午夜免费成人在线视频| 黄片大片在线免费观看| 一级毛片精品| 午夜福利一区二区在线看| 黄色片一级片一级黄色片| 岛国在线观看网站| 亚洲av日韩精品久久久久久密| 午夜久久久在线观看| 老汉色∧v一级毛片| 99国产综合亚洲精品| 国产成人啪精品午夜网站| 999久久久精品免费观看国产| 久久精品人人爽人人爽视色| 中文字幕最新亚洲高清| 精品久久久久久,| 亚洲人成电影观看| 久久人人爽av亚洲精品天堂| 亚洲专区中文字幕在线| 日韩人妻精品一区2区三区| 我的亚洲天堂| 国产在线观看jvid| 国产欧美日韩综合在线一区二区| 欧美日韩瑟瑟在线播放| 88av欧美| 国产精品亚洲av一区麻豆| 黑人巨大精品欧美一区二区mp4| 女性生殖器流出的白浆| 美女扒开内裤让男人捅视频| 国产精品综合久久久久久久免费 | 久久精品国产清高在天天线| 国产区一区二久久| 亚洲精品成人av观看孕妇| 久久久久久久午夜电影 | 色婷婷久久久亚洲欧美| 91麻豆精品激情在线观看国产 | 999久久久精品免费观看国产| 国产成年人精品一区二区 | av中文乱码字幕在线| 亚洲熟妇熟女久久| 99riav亚洲国产免费| 丁香六月欧美| 自线自在国产av| 一个人观看的视频www高清免费观看 | 一二三四社区在线视频社区8| 亚洲成国产人片在线观看| 动漫黄色视频在线观看| 精品福利永久在线观看| 中文字幕av电影在线播放| 日本五十路高清| 男女之事视频高清在线观看| 国产成人免费无遮挡视频| 午夜视频精品福利| 啦啦啦 在线观看视频| 可以免费在线观看a视频的电影网站| 亚洲熟女毛片儿| av在线天堂中文字幕 | 999精品在线视频| 叶爱在线成人免费视频播放| 日日夜夜操网爽| 国产一区在线观看成人免费| 久9热在线精品视频| 一边摸一边抽搐一进一出视频| 亚洲第一青青草原| www.熟女人妻精品国产| 亚洲成人免费av在线播放| 操出白浆在线播放| 天堂俺去俺来也www色官网| а√天堂www在线а√下载| 人妻丰满熟妇av一区二区三区| 日本欧美视频一区| av片东京热男人的天堂| 亚洲成人免费电影在线观看| 交换朋友夫妻互换小说| 久久精品国产综合久久久| 亚洲精品美女久久久久99蜜臀| 久久精品国产99精品国产亚洲性色 | 国产1区2区3区精品| 国产xxxxx性猛交| 亚洲欧美一区二区三区黑人| 久久青草综合色| 神马国产精品三级电影在线观看 | 亚洲国产欧美网| 日韩欧美一区二区三区在线观看| 欧美丝袜亚洲另类 | 青草久久国产| 国产成人精品久久二区二区免费| 午夜福利一区二区在线看| 性少妇av在线| ponron亚洲| a级毛片在线看网站| 一本综合久久免费| 亚洲熟妇中文字幕五十中出 | 在线观看一区二区三区| 日韩一卡2卡3卡4卡2021年| 又紧又爽又黄一区二区| 免费av毛片视频| 欧美日韩福利视频一区二区| 天堂俺去俺来也www色官网| 午夜亚洲福利在线播放| 亚洲欧洲精品一区二区精品久久久| 伊人久久大香线蕉亚洲五| 麻豆av在线久日| 成人永久免费在线观看视频| 午夜福利免费观看在线| 无人区码免费观看不卡| 激情视频va一区二区三区| 亚洲av熟女| 色播在线永久视频| 亚洲久久久国产精品| 亚洲男人的天堂狠狠| 男女午夜视频在线观看| 久久久精品国产亚洲av高清涩受| 欧美+亚洲+日韩+国产| 这个男人来自地球电影免费观看| 国产成+人综合+亚洲专区| 日韩三级视频一区二区三区| 黑人猛操日本美女一级片| 精品久久久久久,| 99久久国产精品久久久| 色哟哟哟哟哟哟| 大香蕉久久成人网| 777久久人妻少妇嫩草av网站| 亚洲性夜色夜夜综合| 91精品国产国语对白视频| 欧美老熟妇乱子伦牲交| 18禁黄网站禁片午夜丰满| 久久精品亚洲精品国产色婷小说| 欧美黑人精品巨大| 一级,二级,三级黄色视频| 高清av免费在线| 欧美人与性动交α欧美软件| 精品人妻在线不人妻| 久久久久久久久久久久大奶| 免费av中文字幕在线| 一级黄色大片毛片| 欧美午夜高清在线| 如日韩欧美国产精品一区二区三区| 最近最新中文字幕大全电影3 | 男女下面进入的视频免费午夜 | 极品人妻少妇av视频| 国产人伦9x9x在线观看| 在线播放国产精品三级| 久久精品亚洲精品国产色婷小说| 男女下面进入的视频免费午夜 | 人人妻人人澡人人看| 亚洲视频免费观看视频| 麻豆av在线久日| 久久人妻av系列| 亚洲一区二区三区色噜噜 | 自线自在国产av| 国产精品 国内视频| 亚洲中文日韩欧美视频| 欧美黑人精品巨大| 国产精品美女特级片免费视频播放器 | 成年版毛片免费区| 99精品欧美一区二区三区四区| 操出白浆在线播放| 一本大道久久a久久精品| 国产又色又爽无遮挡免费看| 亚洲一区二区三区欧美精品| 免费一级毛片在线播放高清视频 | 亚洲国产欧美日韩在线播放| 亚洲av电影在线进入| 国产精品久久电影中文字幕| 丝袜在线中文字幕| 精品久久久久久久毛片微露脸| 国产成人影院久久av| 亚洲欧美精品综合久久99| 夫妻午夜视频| 亚洲国产欧美日韩在线播放| 久久精品影院6| 欧美成人免费av一区二区三区| 女人被躁到高潮嗷嗷叫费观| 欧美日本中文国产一区发布| 色婷婷av一区二区三区视频| 国产精品99久久99久久久不卡| 亚洲成av片中文字幕在线观看| 精品久久久久久电影网| 咕卡用的链子| 神马国产精品三级电影在线观看 | 亚洲va日本ⅴa欧美va伊人久久| 欧美日本中文国产一区发布| 国产精品久久久人人做人人爽| 亚洲国产精品一区二区三区在线| 精品午夜福利视频在线观看一区| 最近最新中文字幕大全电影3 | 性色av乱码一区二区三区2| 后天国语完整版免费观看| 一本综合久久免费| 亚洲精品中文字幕在线视频| 欧美成人性av电影在线观看| 久久人人爽av亚洲精品天堂| 99久久久亚洲精品蜜臀av| 女生性感内裤真人,穿戴方法视频| 操美女的视频在线观看| 在线观看免费高清a一片| 在线视频色国产色| 又大又爽又粗| 午夜日韩欧美国产| 日本 av在线| 国产亚洲欧美98| 亚洲欧美激情在线| 日韩精品免费视频一区二区三区| 婷婷精品国产亚洲av在线| 日韩欧美三级三区| 成人国语在线视频| 日韩免费av在线播放| 亚洲欧美一区二区三区黑人| 我的亚洲天堂| 免费少妇av软件| 一级作爱视频免费观看| 黄色成人免费大全| 久久久水蜜桃国产精品网| 99久久精品国产亚洲精品| 亚洲专区国产一区二区| 久久久精品欧美日韩精品| 自拍欧美九色日韩亚洲蝌蚪91| 一区二区三区激情视频| 欧美日韩中文字幕国产精品一区二区三区 | √禁漫天堂资源中文www| 成人特级黄色片久久久久久久| 日本免费a在线| 高清在线国产一区| 99久久人妻综合| 欧美精品一区二区免费开放| 欧美日韩国产mv在线观看视频| 欧美黄色片欧美黄色片| 成年女人毛片免费观看观看9| 曰老女人黄片| 婷婷六月久久综合丁香| 99在线视频只有这里精品首页| 精品卡一卡二卡四卡免费| 国产真人三级小视频在线观看| 久久性视频一级片| 中文字幕精品免费在线观看视频| 亚洲精品粉嫩美女一区| 最近最新中文字幕大全电影3 | 国产精品亚洲av一区麻豆| 三上悠亚av全集在线观看| 真人一进一出gif抽搐免费| 麻豆成人av在线观看| а√天堂www在线а√下载| 日韩国内少妇激情av| 免费女性裸体啪啪无遮挡网站| 欧美日韩国产mv在线观看视频| 亚洲av美国av| svipshipincom国产片| 亚洲精品久久成人aⅴ小说| 免费久久久久久久精品成人欧美视频| 亚洲精品在线观看二区| 国产熟女午夜一区二区三区| 在线观看www视频免费| 精品免费久久久久久久清纯| 国产成人精品久久二区二区免费| www.999成人在线观看| 国产高清国产精品国产三级| 丁香欧美五月| 成人国语在线视频| 天天影视国产精品| 中文字幕色久视频| 精品国产一区二区三区四区第35| 久99久视频精品免费| 丁香六月欧美| 一区福利在线观看| 日本一区二区免费在线视频| 麻豆成人av在线观看| 极品教师在线免费播放| 日日摸夜夜添夜夜添小说| 天堂动漫精品| 欧美大码av| 嫩草影视91久久| 电影成人av| 国产精品一区二区免费欧美| 日韩欧美一区二区三区在线观看| 很黄的视频免费| 精品少妇一区二区三区视频日本电影| 一个人免费在线观看的高清视频| 精品久久久久久久久久免费视频 | 久久国产亚洲av麻豆专区| 婷婷六月久久综合丁香| 久久久久九九精品影院| 午夜91福利影院| 久久天堂一区二区三区四区| 欧美中文日本在线观看视频| 大型黄色视频在线免费观看| 男人的好看免费观看在线视频 | 欧美精品一区二区免费开放| 777久久人妻少妇嫩草av网站| 99国产综合亚洲精品| 高潮久久久久久久久久久不卡| 亚洲国产看品久久| 女人精品久久久久毛片| 在线十欧美十亚洲十日本专区| 亚洲男人天堂网一区| 热99re8久久精品国产| 天堂俺去俺来也www色官网| 视频在线观看一区二区三区| www日本在线高清视频| 亚洲av电影在线进入| 国产精品久久久久成人av| 中文字幕人妻丝袜制服| 色综合婷婷激情| 国产伦人伦偷精品视频| 国产伦一二天堂av在线观看| 国产一区二区在线av高清观看| 高清在线国产一区| 国产真人三级小视频在线观看| 日本撒尿小便嘘嘘汇集6| 男人舔女人的私密视频| 一级作爱视频免费观看| 在线观看www视频免费| 久久久久国内视频| 欧美精品啪啪一区二区三区| 天天添夜夜摸| 欧美乱码精品一区二区三区| 欧美在线一区亚洲| 色老头精品视频在线观看| 欧美日韩亚洲综合一区二区三区_| 啦啦啦在线免费观看视频4| 美女高潮喷水抽搐中文字幕| 69精品国产乱码久久久| 欧美不卡视频在线免费观看 | 一级黄色大片毛片| 18禁美女被吸乳视频| 91在线观看av| 成人亚洲精品一区在线观看| 夜夜躁狠狠躁天天躁| 亚洲第一av免费看| 亚洲成a人片在线一区二区| 国产成人欧美| 日日干狠狠操夜夜爽| 中文字幕高清在线视频| 国产黄a三级三级三级人| 国产高清国产精品国产三级| 国产主播在线观看一区二区| x7x7x7水蜜桃| 国产熟女xx| 91精品三级在线观看| 丰满人妻熟妇乱又伦精品不卡| 在线视频色国产色| 极品人妻少妇av视频| 国产欧美日韩一区二区三区在线| 超色免费av| 热99国产精品久久久久久7| 在线永久观看黄色视频| 97超级碰碰碰精品色视频在线观看| 久久人人97超碰香蕉20202| 香蕉国产在线看| 一区二区三区激情视频| 视频区欧美日本亚洲| 国产单亲对白刺激| 亚洲自拍偷在线| 久久午夜亚洲精品久久| 精品少妇一区二区三区视频日本电影| 电影成人av| 午夜激情av网站| 电影成人av| 亚洲午夜精品一区,二区,三区| 性色av乱码一区二区三区2| 韩国av一区二区三区四区| 黄色视频,在线免费观看| 欧美中文综合在线视频| av在线播放免费不卡| 欧美丝袜亚洲另类 | 黄频高清免费视频| svipshipincom国产片| 久久亚洲精品不卡| 国产精品美女特级片免费视频播放器 | 国产区一区二久久| 国产精品亚洲一级av第二区| 88av欧美| 久久精品国产99精品国产亚洲性色 | 亚洲精品国产精品久久久不卡| 人人妻人人添人人爽欧美一区卜| 99热只有精品国产| 久久国产精品人妻蜜桃| 99热只有精品国产| 看片在线看免费视频| 亚洲精品美女久久av网站| 国内久久婷婷六月综合欲色啪| 欧美黄色片欧美黄色片| 色播在线永久视频| 免费高清在线观看日韩| 在线观看66精品国产| aaaaa片日本免费|