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

    粳稻種質(zhì)資源芽期耐鹽性綜合評(píng)價(jià)與篩選

    2017-12-02 03:19:09田蕾陳亞萍劉俊馬曉剛王娜楊兵李瑩郭海東李娟胡慧張銀霞李培富
    中國(guó)水稻科學(xué) 2017年6期
    關(guān)鍵詞:芽長(zhǎng)耐鹽耐鹽性

    田蕾 陳亞萍 劉俊 馬曉剛 王娜 楊兵 李瑩 郭海東 李娟 胡慧 張銀霞李培富

    (寧夏大學(xué) 農(nóng)學(xué)院, 銀川 750021; *通訊聯(lián)系人,E-mail:peifuli@163. com)

    粳稻種質(zhì)資源芽期耐鹽性綜合評(píng)價(jià)與篩選

    田蕾 陳亞萍 劉俊 馬曉剛 王娜 楊兵 李瑩 郭海東 李娟 胡慧 張銀霞李培富*

    (寧夏大學(xué) 農(nóng)學(xué)院, 銀川 750021;*通訊聯(lián)系人,E-mail:peifuli@163. com)

    【目的】土壤鹽漬化是危害水稻生產(chǎn)的重要非生物脅迫之一。鑒定水稻種質(zhì)資源發(fā)芽期耐鹽性,篩選耐鹽指標(biāo),培育耐鹽品種,對(duì)水稻生產(chǎn)的發(fā)展具有重要意義。【方法】利用125 mmol/L NaCl溶液對(duì)64份粳稻種質(zhì)資源進(jìn)行鹽脅迫,于脅迫后3 d測(cè)定發(fā)芽數(shù);脅迫5 d、10 d后,測(cè)定發(fā)芽數(shù)、芽長(zhǎng)和根長(zhǎng),并計(jì)算相對(duì)芽長(zhǎng)、根長(zhǎng)、發(fā)芽勢(shì)、發(fā)芽率、鹽害率,發(fā)芽指數(shù)和活力指數(shù)。運(yùn)用多種統(tǒng)計(jì)學(xué)方法對(duì)各種質(zhì)資源的芽期耐鹽性進(jìn)行綜合評(píng)價(jià),分析典型耐鹽和鹽敏感種質(zhì)鹽脅迫條件下的發(fā)芽特征。【結(jié)果】相對(duì)鹽害率與相對(duì)根長(zhǎng)、相對(duì)發(fā)芽勢(shì)、相對(duì)發(fā)芽率、發(fā)芽指數(shù)和活力指數(shù)均極顯著負(fù)相關(guān);除相對(duì)芽長(zhǎng)外各指標(biāo)間的相關(guān)性均達(dá)到極顯著水平。通過聚類分析將64份粳稻種質(zhì)資源劃分成4個(gè)類群。第Ⅰ、Ⅳ類群分別為典型的鹽敏感和耐鹽類群,第Ⅱ類群為弱耐鹽種質(zhì)為主的混合類群,第Ⅲ類群主要由耐鹽種質(zhì)組成。通過主成分分析將7個(gè)評(píng)價(jià)指標(biāo)轉(zhuǎn)換為3個(gè)主成分,應(yīng)用隸屬函數(shù)和權(quán)重,獲得了客觀評(píng)價(jià)粳稻種質(zhì)資源耐鹽性的綜合評(píng)價(jià)值 D。分別選取D值最高和最低的5份種質(zhì)資源進(jìn)行芽期耐鹽指標(biāo)的差異顯著性分析,結(jié)果表明,兩組種質(zhì)資源鹽脅迫5 d的各評(píng)價(jià)指標(biāo)差異均達(dá)到極顯著水平,10 d的評(píng)價(jià)指標(biāo)除相對(duì)芽長(zhǎng)外,均達(dá)到了顯著差異水平?!窘Y(jié)論】水稻芽期對(duì)鹽脅迫較為敏感,且耐鹽性不同的種質(zhì)間差異顯著。利用逐步回歸和主成分分析獲得發(fā)芽指數(shù)、相對(duì)根長(zhǎng)和相對(duì)鹽害率3個(gè)指標(biāo),可作為快速鑒定粳稻種質(zhì)資源芽期耐鹽性的重要指標(biāo),若采用多元統(tǒng)計(jì)方法評(píng)價(jià)可靠性更高。

    水稻;種質(zhì)資源;芽期耐鹽性;主成分分析;綜合評(píng)價(jià)

    鹽堿土在世界范圍內(nèi)分布廣泛,遍及6大洲的100多個(gè)國(guó)家,總面積已達(dá)10億hm2[1-3]。近年來隨著氣候變暖、化肥的大量施用以及不合理灌溉等農(nóng)業(yè)措施的影響,世界約20 %的農(nóng)業(yè)用地鹽堿化程度表現(xiàn)出不斷加劇的趨勢(shì)[4-5]。中國(guó)鹽堿土面積龐大,種類豐富,總面積約為9913.3萬(wàn) hm2[6]。土壤鹽漬化已成為全球性的資源與生態(tài)問題,對(duì)作物的生長(zhǎng)、發(fā)育造成嚴(yán)重影響,制約了我國(guó)農(nóng)業(yè)生產(chǎn)的發(fā)展。

    水稻(Oryza sativa L.)為一年生禾本科植物,是全世界一半以上人口的主要食物來源[7],是一種鹽敏感作物[8],其芽期耐鹽性是鹽漬土條件下水稻穩(wěn)產(chǎn)的重要影響因素[9]。我國(guó)東北、西北和華北的內(nèi)陸干旱和半干旱地區(qū)以及長(zhǎng)江以北的濱海鹽堿稻作區(qū),水稻秧田和直播田在種子萌發(fā)階段易受到鹽堿的危害,導(dǎo)致發(fā)芽勢(shì)降低、發(fā)芽不齊、芽尖枯黃和彎曲,甚至死亡,嚴(yán)重影響了水稻生產(chǎn)。為確保鹽堿地內(nèi)水稻的安全生產(chǎn),創(chuàng)制耐鹽新種質(zhì)、對(duì)其發(fā)芽期耐鹽性的研究顯得尤為重要。

    鹽脅迫對(duì)水稻發(fā)芽期的影響主要表現(xiàn)在對(duì)發(fā)芽速度、整齊度及生長(zhǎng)狀態(tài)的抑制[10]。郭望模等[11]研究發(fā)現(xiàn),不同濃度的NaCl、品種特性以及NaCl濃度和品種間的互作對(duì)發(fā)芽都有極顯著的影響。賈寶艷等[12]研究發(fā)現(xiàn)隨著NaCl濃度的上升,水稻種子發(fā)芽率降低,發(fā)芽時(shí)間明顯推遲,且材料間差異很大。Wang等在對(duì)1個(gè)水稻重組自交系群體的芽期耐鹽性評(píng)價(jià)中發(fā)現(xiàn),水稻種子萌發(fā)初期(0~5d)可能對(duì)鹽脅迫更加敏感[13],各家系的耐鹽相關(guān)指標(biāo)差異較后期更大。Chen等[14]研究發(fā)現(xiàn),鹽脅迫會(huì)明顯降低水稻種子的發(fā)芽率,使芽長(zhǎng)和根長(zhǎng)與空白對(duì)照相比顯著縮短。Hasthanasombut等[15]的研究結(jié)果也類似。近年來,國(guó)內(nèi)外學(xué)者在水稻發(fā)芽期耐鹽性鑒定方面開展了大量工作[9-16],提出了發(fā)芽指數(shù)法,篩選出一些耐鹽性強(qiáng)的優(yōu)異種質(zhì),并應(yīng)用于育種實(shí)踐。但前人的研究多是對(duì)幾個(gè)水稻芽期耐鹽指標(biāo)的簡(jiǎn)單分析,在利用統(tǒng)計(jì)學(xué)方法評(píng)價(jià)水稻芽期耐鹽性方面,多采用聚類[12]、隸屬函數(shù)[16]等方法進(jìn)行簡(jiǎn)單分析,利用多元統(tǒng)計(jì)方法綜合評(píng)價(jià)水稻芽期耐鹽性的報(bào)道較少。本研究以國(guó)內(nèi)外具有廣泛代表性的64份粳稻種質(zhì)資源為試驗(yàn)材料,采用發(fā)芽指數(shù)法,應(yīng)用主成分、隸屬函數(shù)和逐步回歸等方法進(jìn)行各種質(zhì)資源芽期耐鹽性的綜合評(píng)價(jià),以期篩選芽期耐鹽水稻種質(zhì),初步闡明水稻發(fā)芽期的耐鹽機(jī)制,為水稻耐鹽育種提供優(yōu)異資源和理論支持。

    1 材料與方法

    1.1試驗(yàn)材料

    64份粳稻種質(zhì)資源由寧夏大學(xué)農(nóng)學(xué)院遺傳育種實(shí)驗(yàn)室李培富教授提供,其中,源自亞洲38個(gè)、歐洲19個(gè)、大洋洲4 個(gè)、美洲2 個(gè)、非洲1 個(gè)(表1)。

    1.2試驗(yàn)方法

    試驗(yàn)于 2015年在寧夏大學(xué)農(nóng)學(xué)院遺傳育種實(shí)驗(yàn)室進(jìn)行。將各種質(zhì)資源種子置于42℃恒溫培養(yǎng)箱內(nèi)處理 2 d,以充分打破種子休眠。每份種子分別選取成熟度好、整齊、飽滿的種子 50粒,均勻地置于底部鋪有 2層無(wú)菌濾紙的黑色育苗盤(28 cm×56 cm)中,每份種質(zhì)1穴,每盤32穴,每穴直徑6 cm。用5 %次氯酸鈉溶液浸泡30 min對(duì)種子進(jìn)行消毒,蒸餾水沖洗3~5次。設(shè)置對(duì)照組(蒸餾水)和鹽脅迫組(NaCl溶液)兩個(gè)處理,每處理3個(gè)重復(fù)。黑色育苗盤用塑料盤墊底,每盤加入800 mL 處理液后用塑料保鮮膜密封,放入智能人工氣候箱(HP1500GS型,武漢瑞華儀器設(shè)備有限公司)中進(jìn)行發(fā)芽培養(yǎng),培養(yǎng)條件為28℃下光照12 h,26℃下暗培養(yǎng)12 h。每天進(jìn)行溶液更換以保證脅迫濃度的穩(wěn)定,脅迫處理10d。

    為了選擇最佳的鹽脅迫濃度,根據(jù)各種質(zhì)資源在寧夏石嘴山市平羅縣西大灘鹽堿地改良示范基地鹽堿直播田中的出苗率,選擇了3份出苗率較高的種質(zhì) OB22(漾濞光殼陸稻)、OB23(Bertone)和OB28(Galhardo),出苗率分別為 88.9%、83.3%和81.9%,作為芽期耐鹽種質(zhì);3份出苗率較低的種質(zhì)OB1(合江 21)、OB12(Rizzotto)和 OB18(荒木),出苗率分別為 19.4%、25.0%和 29.2%,作為芽期鹽敏感種質(zhì)。設(shè)置NaCl濃度梯度為0、50、100、125、150、175、200 mmol/L 進(jìn)行預(yù)實(shí)驗(yàn),篩選最佳NaCl濃度用于64份種質(zhì)資源芽期耐鹽性評(píng)價(jià)。

    表1 64份粳稻種質(zhì)資源名稱、來源、綜合評(píng)價(jià)值及排名Table 1. Origin and names of 64 japonica rice germplasm and their D values,comprehensive ranking.

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

    以種子胚根達(dá)到 1 mm 長(zhǎng)為發(fā)芽標(biāo)準(zhǔn)[17],在NaCl脅迫3、5、10 d時(shí)調(diào)查各處理的發(fā)芽數(shù),脅迫5、10 d時(shí)調(diào)查各處理的芽長(zhǎng)和根長(zhǎng)。以相對(duì)芽長(zhǎng)、相對(duì)根長(zhǎng)、相對(duì)發(fā)芽勢(shì)、相對(duì)發(fā)芽率、相對(duì)鹽害率、發(fā)芽指數(shù)和活力指數(shù)作為發(fā)芽期耐鹽性的評(píng)價(jià)指標(biāo),對(duì)各種質(zhì)資源的芽期耐鹽性進(jìn)行評(píng)價(jià),其中發(fā)芽指數(shù)用鹽脅迫3、5、10 d的發(fā)芽數(shù)進(jìn)行計(jì)算。

    發(fā)芽勢(shì)(%)=(培養(yǎng) 5 d種子發(fā)芽數(shù)/供試種子粒數(shù))×100

    相對(duì)發(fā)芽勢(shì)(%) =(NaCl處理種子發(fā)芽勢(shì)/空白對(duì)照種子發(fā)芽勢(shì))×100;

    發(fā)芽率(%)=(培養(yǎng)10 d種子發(fā)芽數(shù)/供試種子粒數(shù))×100;

    相對(duì)發(fā)芽率(%)=(NaCl處理種子發(fā)芽率/空白對(duì)照種子發(fā)芽率)×100;

    相對(duì)芽長(zhǎng)(%)=(NaCl處理種子芽長(zhǎng)/空白對(duì)照種子芽長(zhǎng)) ×100;

    相對(duì)根長(zhǎng)(%)=(NaCl處理種子根長(zhǎng)/空白對(duì)照種子根長(zhǎng)) ×100;

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

    發(fā)芽指數(shù) GI=∑(Gt/Dt),(Gt指第 t天種子發(fā)芽數(shù),Dt指相應(yīng)發(fā)芽天數(shù));[17,18]

    活力指數(shù)VI=GI×S,(S指種苗生長(zhǎng)量,為脅迫10 d的根長(zhǎng))[19]。

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

    使用Excel 2013整理試驗(yàn)數(shù)據(jù),利用SPSS 19.0和SAS 9.2軟件對(duì)7個(gè)水稻芽期耐鹽相關(guān)性狀指標(biāo)進(jìn)行獨(dú)立樣本t檢驗(yàn)、相關(guān)性、差異顯著性、聚類、主成分和逐步線性回歸分析,利用隸屬函數(shù)和標(biāo)準(zhǔn)差系數(shù)賦予權(quán)重法進(jìn)行耐鹽性的綜合評(píng)價(jià),相關(guān)指標(biāo)計(jì)算公式及方法參照文獻(xiàn)[20-21]。

    1.4.1 粳稻種質(zhì)資源各綜合指標(biāo)的隸屬函數(shù)值隸屬函數(shù)值計(jì)算公式:

    式中Xj表示第j個(gè)綜合指標(biāo),Xmin、Xmax為所有參試材料某一指標(biāo)的最小值和最大值。1.4.2 各綜合指標(biāo)的權(quán)重

    式中,Wj表示第j 個(gè)綜合指標(biāo)在所有綜合指標(biāo)中的重要程度即權(quán)重;Pj代表經(jīng)主成分分析所得到各粳稻種質(zhì)資源第j 個(gè)綜合指標(biāo)的貢獻(xiàn)率。

    1.4.3 粳稻種質(zhì)資源綜合耐鹽能力

    式中,D值表示各粳稻種質(zhì)資源在鹽脅迫下由綜合指標(biāo)評(píng)價(jià)所得到的發(fā)芽期耐鹽性綜合評(píng)價(jià)值。

    2 結(jié)果與分析

    2.1粳稻種質(zhì)資源芽期耐鹽性評(píng)價(jià)鹽脅迫最適濃度篩選

    對(duì)3個(gè)芽期耐鹽和3個(gè)芽期鹽敏感種質(zhì)資源進(jìn)行不同濃度NaCl溶液脅迫后發(fā)現(xiàn),隨著NaCl濃度的提高,各種質(zhì)資源的7個(gè)芽期耐鹽相關(guān)指標(biāo)均顯著降低。除100 mmol/L NaCl溶液脅迫10 d 的相對(duì)根長(zhǎng)和活力指數(shù)外,兩類種質(zhì)資源各芽期耐鹽相關(guān)指標(biāo)在較低的NaCl濃度(0、50、100 mmol/L)和較高的濃度(175、200 mmol/L)下,差異均未達(dá)到極顯著水平;而125 mmol/L NaCl濃度下,二者之間各芽期耐鹽相關(guān)指標(biāo)差異均達(dá)到極顯著水平(表2),故選擇該濃度對(duì)64份粳稻種質(zhì)資源的芽期耐鹽性進(jìn)行評(píng)價(jià)。

    2.2粳稻種質(zhì)資源鹽脅迫條件下芽期各評(píng)價(jià)指標(biāo)表現(xiàn)

    t檢驗(yàn)結(jié)果表明,各種質(zhì)資源的7個(gè)芽期耐鹽相關(guān)性狀在鹽處理?xiàng)l件下與空白對(duì)照相比,均表現(xiàn)出極顯著的差異(表3)。t值均為正值,表明與空白條件下相比,鹽脅迫嚴(yán)重影響了水稻種子的萌發(fā),各指標(biāo)的表型值均明顯下降。除鹽脅迫10 d的相對(duì)根長(zhǎng)外,各指標(biāo)均表現(xiàn)出較大的分布范圍,表明這些指標(biāo)在各種質(zhì)資源間存在著較大的變異。其中,鹽脅迫5 d相對(duì)發(fā)芽勢(shì)和相對(duì)鹽害率的分布范圍最廣,分別為0.8%~100.0%和0.0%~99.2%。

    2.3粳稻種質(zhì)資源鹽脅迫條件下芽期各評(píng)價(jià)指標(biāo)與苗期耐鹽級(jí)別的相關(guān)分析

    與正常培養(yǎng)條件相比,鹽脅迫下水稻種子萌發(fā)的形態(tài)學(xué)特征有不同程度的變化,通過形態(tài)生理指標(biāo)間的相關(guān)性分析可揭示指標(biāo)間是否存在依存關(guān)系及相關(guān)關(guān)系的方向與強(qiáng)度[22]。本研究前期已完成64份粳稻種質(zhì)資源苗期耐鹽性評(píng)價(jià)[23,24],利用苗期耐鹽級(jí)別(STS)與芽期各評(píng)價(jià)指標(biāo)進(jìn)行相關(guān)性分析。由表4可以看出,相對(duì)鹽害率與相對(duì)根長(zhǎng)、相對(duì)發(fā)芽勢(shì)、相對(duì)發(fā)芽率、發(fā)芽指數(shù)和活力指數(shù)均呈極顯著負(fù)相關(guān),其中相對(duì)鹽害率與發(fā)芽指數(shù)的相關(guān)系數(shù)最大,為-0.860;相對(duì)芽長(zhǎng)與其他指標(biāo)間的相關(guān)性未達(dá)到顯著水平。苗期耐鹽級(jí)別與相對(duì)芽長(zhǎng)、相對(duì)根長(zhǎng)、相對(duì)發(fā)芽勢(shì)、發(fā)芽指數(shù)均呈顯著正相關(guān),與活力指數(shù)呈極顯著正相關(guān)。

    根據(jù)芽期各指標(biāo)間相關(guān)分析,我們發(fā)現(xiàn)多個(gè)評(píng)價(jià)指標(biāo)間存在極顯著相關(guān)性,會(huì)使其提供的信息發(fā)生重疊,各指標(biāo)變化幅度參差不齊,且發(fā)揮作用的大小也不相同。水稻耐鹽性是一個(gè)復(fù)雜的綜合性狀,直接利用單項(xiàng)指標(biāo)不能準(zhǔn)確、直觀地進(jìn)行耐鹽性評(píng)價(jià),可利用主成分分析對(duì)供試粳稻種質(zhì)資源芽期耐鹽性進(jìn)行深入的綜合評(píng)價(jià)。

    2.4粳稻種質(zhì)資源鹽脅迫條件下芽期耐鹽性的綜合評(píng)價(jià)

    2.4.1 主成分分析

    對(duì)7個(gè)水稻芽期耐鹽單項(xiàng)指標(biāo)進(jìn)行主成分分析(表 5),CI1、CI2兩個(gè)主成分的特征值分別為 4.37和1.20,均大于1,其貢獻(xiàn)率分別為62.5%和17.1%,表明第一主成分能夠提供原變量綜合信息的一半以上。第一、二、三主成分的累計(jì)貢獻(xiàn)率為91.6%,已對(duì)大多數(shù)指標(biāo)作了充分的概括,因此,提取了 3個(gè)主成分,將原有的7個(gè)芽期評(píng)價(jià)指標(biāo)轉(zhuǎn)換為3個(gè)新的相互獨(dú)立的綜合指標(biāo),來代表原始指標(biāo)攜帶的絕大部分信息。

    表2 典型耐鹽和鹽敏感粳稻種質(zhì)資源在不同濃度NaCl脅迫條件下芽期各評(píng)價(jià)指標(biāo)表現(xiàn)Table 2. Phenotypic values of typical salt-tolerant and sensitive rice germplasm accessions under different NaCl concentrations at seed germination stage.

    表3 粳稻種質(zhì)資源鹽脅迫條件下芽期各評(píng)價(jià)指標(biāo)均值及分布范圍Table 3. Performance of salt tolerance-related traits of 64 rice germplasm accessions at seed germination stage

    由表6可知,在第一主成分中,相對(duì)發(fā)芽勢(shì)(X3)、相對(duì)發(fā)芽率(X4)、發(fā)芽指數(shù)(X5)和相對(duì)鹽害率(X7)的載荷絕對(duì)值較高,均在0.4以上,是第一主成分中的主要作用因子;第二主成分中,相對(duì)芽長(zhǎng)(X1)、相對(duì)根長(zhǎng)(X2)的載荷絕對(duì)值較高,分別為 0.707和0.490;相對(duì)芽長(zhǎng)(X1)、相對(duì)根長(zhǎng)(X2)和活力指數(shù)(X6)為第三主成分的主要作用因子。綜合3個(gè)主成分所攜帶的信息,可用相對(duì)發(fā)芽勢(shì)、相對(duì)發(fā)芽率、發(fā)芽指數(shù)及相對(duì)鹽害率4個(gè)指標(biāo)來概括7個(gè)指標(biāo)中的大部分信息。

    2.4.2 隸屬函數(shù)及粳稻種質(zhì)資源芽期耐鹽性綜合評(píng)價(jià)

    利用綜合指標(biāo)值(CI1-CI3)和公式(1)分別計(jì)算64個(gè)粳稻種質(zhì)資源的隸屬函數(shù)值u(Xj)。對(duì)于同一綜合指標(biāo)如CI1而言,鹽脅迫處理下,OB34(Daniela)的u(X1)最大,其值為1,表明此種質(zhì)在主成分一表現(xiàn)為耐鹽性最強(qiáng),而 OB15(寧粳 16)的u(X1)最小,其值為 0,表明該種質(zhì)在這一綜合指標(biāo)上表現(xiàn)為耐鹽性最差。依據(jù)芽期各綜合指標(biāo)貢獻(xiàn)率的大小利用公式(2)求出芽期3個(gè)綜合指標(biāo)的權(quán)重,分別為0.682、0.187和0.131。利用公式(3)以隸屬函數(shù)u(X1)、u(X2)和u(X3) 結(jié)合權(quán)重處理并累加得到綜合評(píng)價(jià) D值(表1),并依此進(jìn)行芽期耐鹽性的排名。

    表4 鹽脅迫處理下水稻芽期各單項(xiàng)指標(biāo)與苗期耐鹽級(jí)別的相關(guān)系數(shù)矩陣Table 4. Correlation matrix of each single index of rice seed germination and salt tolerance score under salt stress.

    表5 各綜合指標(biāo)的特征值及貢獻(xiàn)率Table 5. Eigen values and proportion of comprehensive indexes [Clx].

    表6 各因子載荷矩陣Table 6. Loading matrix of each component.

    綜合評(píng)價(jià)值 D的大小客觀反應(yīng)了各粳稻種質(zhì)資源芽期耐鹽能力的大小,數(shù)值越大說明耐鹽能力越強(qiáng)。分別選取綜合排名前5名和最后5名作為高D值種質(zhì)資源和低D值種質(zhì)資源,即OB34(Daniela)、OB37(Gostima)、 OB45(Hrborio Cyauco)、 OB35(Arborio)和OB41(Rossi),其D值均在0.820以上,排名第一的Daniela芽期綜合耐鹽能力最強(qiáng);低D值種質(zhì)資源為OB30(遼豐8號(hào))、OB64(ARC7042)、OB48(嘉農(nóng) 485)、OB1(合江 21)和 OB15(寧粳 16),D值均在0.370以下,耐鹽性相對(duì)較弱。其中,寧粳16號(hào)芽期綜合耐鹽能力最弱。

    2.5粳稻種質(zhì)資源芽期耐鹽性形態(tài)指標(biāo)聚類分析

    以D值為依據(jù)將64份粳稻種質(zhì)資源進(jìn)行類群劃分:D值范圍在0.211~0.510之間有10份種質(zhì),為鹽敏感種質(zhì);D值范圍在0.511~0.650之間有23份種質(zhì),為弱耐鹽種質(zhì);D 值范圍在 0.651~0.750之間有 21份種質(zhì),為耐鹽種質(zhì);D值范圍在0.751~0.961之間有10份種質(zhì),為高度耐鹽種質(zhì)。為驗(yàn)證用D值評(píng)價(jià)水稻芽期耐鹽性的效果,用除相對(duì)芽長(zhǎng)以外的6個(gè)形態(tài)指標(biāo)數(shù)據(jù),參照戴海芳等[20]的方法進(jìn)行系統(tǒng)聚類分析,發(fā)現(xiàn)在歐氏距離1.08處可劃分為兩大類群,分別為第Ⅰ類群和由Ⅱ、Ⅲ、Ⅳ組成的混合類群;在歐氏距離0.87時(shí)可將混合類群細(xì)劃為Ⅱ、Ⅲ、Ⅳ 3個(gè)類群(圖1)。

    圖1 基于6個(gè)芽期耐鹽指標(biāo)的64份粳稻種質(zhì)資源聚類分析結(jié)果Fig. 1. Cluster analysis of 64 rice germplasm based on the six salt response indexes.

    第Ⅰ類群為鹽敏感類群,包括 OB1、OB30、OB19、OB31、OB64、OB48和 OB15,這 7 份種質(zhì)資源的D值均小于0.50,皆為鹽敏感種質(zhì),原產(chǎn)地均在亞洲,在鹽脅迫條件下,該類群種質(zhì)受到抑制最為明顯。

    第Ⅱ類群為弱耐鹽種質(zhì)為主的混合類群,包括8份弱耐鹽種質(zhì)和3份鹽敏感種質(zhì),除OB27、OB18和OB5以外,均為弱耐鹽種質(zhì),其D值分布范圍為0.484~0.548。第Ⅲ類群為中度耐鹽類群,由34 份種質(zhì)組成,包括2個(gè)高度耐鹽種質(zhì),17個(gè)耐鹽種質(zhì)和15個(gè)弱耐鹽種質(zhì),D值分布范圍在0.570~0.795之間,是耐鹽種質(zhì)相對(duì)集中的1個(gè)類群。第Ⅳ類群為高度耐鹽類群,包括8個(gè)高度耐鹽種質(zhì)和4個(gè)耐鹽種質(zhì),D值分布范圍為0.702~0.961,具有較強(qiáng)的耐鹽性。除鹽敏感類群均來自亞洲外,未發(fā)現(xiàn)芽期耐鹽性與種質(zhì)來源、地理分布有明顯的相關(guān)性。

    2.6低D值種質(zhì)資源與高D值種質(zhì)資源各指標(biāo)的差異顯著性分析

    由于發(fā)芽率與發(fā)芽勢(shì)都可表示為一定條件下第n天的發(fā)芽種子數(shù)占供試種子數(shù)的百分率,因此在對(duì)高(低)D值種質(zhì)資源各指標(biāo)差異顯著性分析時(shí)統(tǒng)一記為相對(duì)發(fā)芽率。對(duì)5份高D值種質(zhì)資源、5份低D值種質(zhì)資源的7個(gè)芽期耐鹽相關(guān)指標(biāo)分別進(jìn)行方差分析,F(xiàn)值的分布范圍為3.35~674.90;除鹽脅迫10 d的相對(duì)芽長(zhǎng)P值為0.078之外,其他各指標(biāo)P值均小于0.001。利用最小顯著差異法(LSD)進(jìn)行多重比較,結(jié)果表明,鹽脅迫5 d后,低D值種質(zhì)資源和高D值種質(zhì)資源的各指標(biāo)均表現(xiàn)出極顯著差異(圖2),尤其是相對(duì)鹽害率二者差異最為明顯。鹽處理10 d后,相對(duì)芽長(zhǎng)的差異未達(dá)到顯著水平,其他各指標(biāo)均表現(xiàn)為極顯著差異。表明兩類種質(zhì)在鹽脅迫5 d時(shí)已有顯著差異,主要表現(xiàn)在發(fā)芽速度的快慢,芽和根系伸長(zhǎng)的程度等方面,隨著發(fā)芽時(shí)間的延長(zhǎng),發(fā)芽率和芽長(zhǎng)的差異逐漸減小,主要表現(xiàn)為根系伸長(zhǎng)的差異(圖2-A~B)。

    圖2 125 mmol/L NaCl 脅迫下典型耐(敏)鹽粳稻種質(zhì)資源芽期形態(tài)指標(biāo)表現(xiàn)Fig. 2. Morphology indexes of typical salt tolerant andsensitive japonica rice germplasm at germination stage under salt stress.

    2.7水稻芽期和苗期高度耐鹽種質(zhì)資源篩選

    芽期耐鹽性綜合排名7、8的漾濞光殼陸稻和湟羅表現(xiàn)出很強(qiáng)的芽期耐鹽性(圖3-A、表1),其鹽脅迫10 d 的相對(duì)根長(zhǎng)分別為7.8%和8.1%,發(fā)芽指數(shù)分別為14.6和14.5,相對(duì)鹽害率分別為1.0%和2.0%;它們的苗期耐鹽性也很出色,耐鹽級(jí)別分別為6.3和6.6。此外,綜合排名14的Bertone也是一個(gè)苗期高度耐鹽種質(zhì),耐鹽級(jí)別7.5[24],125 mmol/L NaCl脅迫6 d,僅第1、2葉萎蔫,有3~4片完整的綠葉(圖3-B),可作為水稻耐鹽育種的優(yōu)良供體親本。

    2.8回歸模型的建立及芽期耐鹽主要評(píng)價(jià)指標(biāo)篩選

    利用已獲得芽期耐鹽綜合評(píng)價(jià)D值和7個(gè)芽期耐鹽相關(guān)指標(biāo),建立可用于水稻芽期耐鹽性評(píng)價(jià)的數(shù)學(xué)模型,把D值作為因變量,把各單項(xiàng)指標(biāo)作為自變量進(jìn)行逐步回歸分析,建立最優(yōu)回歸方程:D=-0.083+0.025X5+0.033X2+0.002X4+0.001X1。方程相關(guān)系數(shù)R=0.993,決定系數(shù)R2=0.985,F(xiàn)值為978.606。由方程可知,7個(gè)指標(biāo)中有4個(gè)指標(biāo)對(duì)水稻芽期耐鹽性有顯著影響,分別是X5(發(fā)芽指數(shù))、X2(相對(duì)根長(zhǎng))、X4(相對(duì)發(fā)芽率)和X1(相對(duì)芽長(zhǎng)),它們與D值均呈極顯著相關(guān),P值除相對(duì)芽長(zhǎng)為0.007外均小于0.001。結(jié)合這4個(gè)指標(biāo)的偏回歸系數(shù),發(fā)芽指數(shù)和相對(duì)根長(zhǎng)可作為水稻芽期耐鹽性評(píng)價(jià)的關(guān)鍵指標(biāo)。

    3 討論

    3.1鹽脅迫對(duì)不同耐鹽性水稻種子萌發(fā)的影響

    植物種子在萌發(fā)期受到鹽脅迫時(shí),發(fā)芽時(shí)間會(huì)延遲,發(fā)芽率降低,生育進(jìn)程受到抑制[11,25]。肖文斐等[10]在對(duì)31個(gè)秈稻品種的芽期耐鹽性評(píng)價(jià)中發(fā)現(xiàn),各品種在9 g/L NaCl溶液處理下發(fā)芽起始時(shí)間明顯延遲,發(fā)芽變緩,發(fā)芽指數(shù)和活力指數(shù)均顯著降低。本研究在對(duì)64份粳稻種質(zhì)資源的芽期鹽脅迫中發(fā)現(xiàn),與空白對(duì)照相比,各種質(zhì)的發(fā)芽起始時(shí)間普遍延遲了0.5~1.5 d,發(fā)芽時(shí)間均比空白對(duì)照的3~4 d延長(zhǎng)了1 d以上,該現(xiàn)象在鹽敏感種質(zhì)中更為明顯。如OB48(嘉農(nóng)485)在鹽脅迫10 d時(shí)相對(duì)發(fā)芽率僅為58.9 %,其種子萌發(fā)受到了很大程度的抑制。崔江慧等在對(duì)10份高粱材料的鹽處理時(shí)發(fā)現(xiàn),側(cè)根鹽害率的變化最大,表明側(cè)根的發(fā)生和伸長(zhǎng)對(duì)鹽脅迫最敏感[26],水稻中也有類似的報(bào)道[13]。本研究中,各種質(zhì)資源10 d的相對(duì)根長(zhǎng)分布范圍為0.6 %~13.4 %(表3),根的伸長(zhǎng)受到了很嚴(yán)重的影響,典型耐鹽種質(zhì)的相對(duì)根長(zhǎng)極顯著大于鹽敏感種質(zhì)(圖2-B);而10 d的相對(duì)芽長(zhǎng)分布范圍為7.3 %~40.7 %,表明水稻種子萌發(fā)時(shí)根的伸長(zhǎng)對(duì)鹽脅迫最為敏感,芽次之,這與前人的研究結(jié)果一致。

    圖3 125 mmol/L NaCl脅迫條件下典型耐鹽和鹽敏感水稻種質(zhì)資源的表現(xiàn)型Fig. 3. Phenotypes of two salt tolerant rice germplasm accessions and two salt sensitive rice germplasm accessions under 125 mmol/L NaCl for five and six days,respectively in germination and seedling stages.

    3.2水稻芽期耐鹽性評(píng)價(jià)關(guān)鍵指標(biāo)的選擇

    作物種子萌發(fā)期進(jìn)行耐鹽評(píng)價(jià)時(shí),發(fā)芽勢(shì)、發(fā)芽率、芽長(zhǎng)、根長(zhǎng)、發(fā)芽指數(shù)、活力指數(shù)等是常用的耐鹽評(píng)價(jià)指標(biāo)[27]。眾多學(xué)者利用上述指標(biāo)及其相對(duì)值對(duì)水稻[9-16]、高粱[26,28]、偃麥草[27]、小麥[29]和苜蓿[18]等作物的芽期耐鹽性開展了大量研究,但對(duì)芽期耐鹽評(píng)價(jià)的關(guān)鍵指標(biāo)未達(dá)成統(tǒng)一的結(jié)論。賈寶艷等[12]認(rèn)為,發(fā)芽指數(shù)、發(fā)芽率和相對(duì)受害率是衡量水稻芽期耐鹽能力的主要指標(biāo)。馮鐘慧等[16]認(rèn)為相對(duì)芽長(zhǎng)和相對(duì)根長(zhǎng)有望成為水稻早期耐鹽性快速鑒定的主要指標(biāo)。本研究通過對(duì)綜合評(píng)價(jià)D值與7個(gè)單項(xiàng)指標(biāo)的逐步線性回歸分析,獲得了1個(gè)最優(yōu)回歸方程,結(jié)合方程中4個(gè)自變量指標(biāo)的偏回歸系數(shù),篩選到了發(fā)芽指數(shù)和相對(duì)根長(zhǎng)兩個(gè)關(guān)鍵指標(biāo)。同時(shí),結(jié)合主成分分析中第一主成分的貢獻(xiàn)率及其各載荷因子的系數(shù),以及典型耐(敏)鹽種質(zhì)各指標(biāo)的差異顯著性分析,認(rèn)為相對(duì)鹽害率也是一個(gè)很重要的水稻芽期耐鹽相關(guān)指標(biāo)。利用上述3個(gè)指標(biāo)結(jié)合相近的NaCl濃度,不同學(xué)者進(jìn)行了粳稻[11,12,16]、秈稻[9-11,13]、爪哇稻[11]的芽期耐鹽性評(píng)價(jià),均取得了較好地鑒定結(jié)果,表明本研究所建立的粳稻種質(zhì)資源芽期耐鹽性評(píng)價(jià)方法和重要指標(biāo)也可用于秈稻等其他類型水稻的芽期耐鹽性評(píng)價(jià)。

    3.3水稻芽期耐鹽性的多元統(tǒng)計(jì)學(xué)分析

    作物耐鹽性是受諸多因素影響的復(fù)雜性狀,是外界環(huán)境與其自身基因型互作的結(jié)果,簡(jiǎn)單地使用某些單項(xiàng)指標(biāo)很難全面的反映其耐鹽本質(zhì)[20]。本研究選用7個(gè)水稻芽期耐鹽指標(biāo)進(jìn)行測(cè)定,并在此基礎(chǔ)上利用多元統(tǒng)計(jì)分析對(duì)粳稻種質(zhì)資源芽期耐鹽能力進(jìn)行了綜合評(píng)價(jià)。利用主成分分析,將7個(gè)指標(biāo)轉(zhuǎn)化為3個(gè)彼此獨(dú)立的新指標(biāo);結(jié)合隸屬函數(shù)及權(quán)重獲得了可以客觀反映各種質(zhì)耐鹽性的綜合評(píng)價(jià)D值,利用該指標(biāo)將64份粳稻種質(zhì)資源劃分為高度耐鹽、耐鹽、弱耐鹽和鹽敏感4種類型。利用逐步回歸構(gòu)建了1個(gè)較為可靠的水稻芽期耐鹽性評(píng)價(jià)模型,利用該模型,可在相同的脅迫條件下預(yù)測(cè)目標(biāo)種質(zhì)耐鹽性的強(qiáng)弱,為水稻耐鹽育種及資源的鑒定提供一定的依據(jù)。

    3.4水稻種質(zhì)資源芽期與苗期耐鹽的相關(guān)性分析及耐鹽種質(zhì)在育種中的應(yīng)用策略

    吳家富等[30]的研究表明,水稻耐鹽性有明顯的階段發(fā)育特異性,不同發(fā)育階段的耐鹽性之間無(wú)明顯相關(guān)性。本研究通過相關(guān)分析發(fā)現(xiàn),苗期耐鹽級(jí)別與包括D值在內(nèi)的多個(gè)芽期耐鹽相關(guān)指標(biāo)表現(xiàn)為極顯著或顯著正相關(guān)(表4),與D值的相關(guān)系數(shù)為0.334,呈極顯著正相關(guān),與前人的研究結(jié)果存在一定的差異,這可能與選用的種質(zhì)資源材料特性不同有關(guān)。

    水稻直播對(duì)實(shí)現(xiàn)稻米生產(chǎn)的專業(yè)化、輕型化、規(guī)?;兄匾饬x。近年來我國(guó)水稻直播面積不斷擴(kuò)大[31],因此較強(qiáng)的芽期耐鹽性是鹽漬地區(qū)水稻直播品種必須優(yōu)先具備的特性。這就要求我們?cè)谟N實(shí)踐中要綜合考慮種質(zhì)資源的特性并合理利用。“設(shè)計(jì)育種”最早由Peleman等提出,是通過有目的地聚合有利基因,培育和創(chuàng)造優(yōu)良品種或資源的育種方法。根據(jù)Li等[32]提出的分子設(shè)計(jì)育種策略,可選擇生產(chǎn)中大面積推廣的品種(如富源4號(hào))作為輪回親本,以耐鹽水稻種質(zhì)資源(如本研究篩選到的3個(gè)芽期和苗期都耐鹽的種質(zhì))為供體,培育高代回交群體,結(jié)合分子標(biāo)記輔助選擇和客觀準(zhǔn)確的耐鹽篩選和鑒定,獲得可直接用于生產(chǎn)的優(yōu)良耐鹽滲入系。同時(shí)通過基因組學(xué)、轉(zhuǎn)錄組學(xué)和代謝組學(xué)等手段可進(jìn)行目標(biāo)基因的克隆及功能標(biāo)記開發(fā),進(jìn)一步改良水稻品種的耐鹽性。寧粳16和合江21是芽期和苗期高度鹽敏感的種質(zhì)(圖3),可用于水稻芽期、苗期鹽敏感基因的分離和克隆。

    4 結(jié)論

    在人工氣候箱中利用125 mmol/L NaCl溶液對(duì)64 份粳稻種質(zhì)資源芽期耐鹽性進(jìn)行鑒定評(píng)價(jià)。利用逐步回歸和主成分分析等統(tǒng)計(jì)學(xué)方法,從測(cè)定的7個(gè)單項(xiàng)指標(biāo)中篩選出發(fā)芽指數(shù)、相對(duì)根長(zhǎng)、相對(duì)鹽害率3個(gè)關(guān)鍵指標(biāo),建立了回歸方程,可用于水稻種質(zhì)資源芽期耐鹽性的快速鑒定與預(yù)測(cè)。利用綜合評(píng)價(jià)D值獲得了5份芽期高度耐鹽種質(zhì),為水稻種質(zhì)資源篩選與鑒定,耐鹽品種的選育、推廣提供了材料和依據(jù)。

    [1] 張翼夫,李洪文,胡紅,陳婉芝,王憲良,牛琪. 鹽堿地后懸掛隨動(dòng)式打孔通氣機(jī)的設(shè)計(jì)與試驗(yàn). 農(nóng)業(yè)工程學(xué)報(bào),2016,32(18): 42-49.Zhang Y F,Li H W,Hu H,Chen W Z,Wang X L,Niu Q.Design and experiment on rear suspended passive aerator in saline-alkali land. Trans Chin Soc Agric Eng,2016,32(18): 42-49. (in Chinese with English abstract)

    [2] Wang X,Wang J G,Liu H L,Zou D T,Zhao H W.Influence of natural saline-alkali stress on chlorophyll content and chloroplast ultrastructure of two contrasting rice (Oryza sativa L. japonica) cultivars. Aust J Crop Sci,2013,7(2): 289-292.

    [3] 王善仙,劉宛,李培軍,吳海燕. 鹽堿土植物改良研究進(jìn)展. 中國(guó)農(nóng)學(xué)通報(bào),2011,27(24): 1-7.Wang S X,Liu W,Li P J,Wu H Y. Advances of researches in plant-improvement of saline-alkaline soil.Chin Agric Sci Bull,2011,27(24): 1-7. (in Chinese with English abstract)

    [4] Yu J B,Wang Z C,Meixner F X,Yang F,Wu H F,Chen X B. Biogeochemical characterizations and reclamation strategies of saline sodic soil in northeastern China.Clean–Soil,Air,Water,2010,38(11): 1010-1016.

    [5] Tian L,Tan L B,Liu F X,Cai H W,Sun C Q.Identification of quantitative trait loci associated with salt tolerance at seedling stage from Oryza rufipogon. J Genet Genom,2011,38(12): 593-601.

    [6] 曾玲玲,季生棟,王俊強(qiáng),張成亮,趙剛. 植物耐鹽機(jī)理的研究進(jìn)展. 黑龍江農(nóng)業(yè)科學(xué),2009(3): 156-159.Zeng L L,Ji S D,Wang J Q,Zhang C L,Zhao G.Advances on the mechanism of plant salt tolerance.Heilongjiang Agric Sci,2009(3): 156-159. (in Chinese with English abstract)

    [7] Wang Z F,Cheng J P,Chen Z W,Huang J,Bao Y M,Wang J F,Zhang H S. Identification of QTLs with main,epistatic and QTL×environment interaction effects for salt tolerance in rice seedlings under different salinity conditions. Theor Appl Genet,2012,125(4): 807-815.

    [8] Ghaffaria A,Gharechahib J,Nakhodaa B,Salekdeh G H.Physiology and proteome responses of two contrasting rice mutants and their wild type parent under salt stress conditions at the vegetative stage. J Plant Physiol,2014,171(1): 31-44.

    [9] Wang Z F,Wang J F,Bao Y M,Wu Y Y,Zhang H S.Quantitative trait loci controlling rice seed germination under salt stress. Euphytica,2011,178 (3): 297-307.

    [10] 肖文斐,馬華升,陳文岳,裘劼人,童建新,鄭桂珍,忻雅,王淑珍,方獻(xiàn)平,阮松林. 秈稻耐鹽性與稻米品質(zhì)性狀的關(guān)聯(lián)分析. 核農(nóng)學(xué)報(bào),2013,27(12):1938-1947.Xiao W F,Ma H S,Chen W Y,Qiu J R,Tong J X,Zheng G Z,Xin Y,Wang S Z,Fang X P,Ruan S L. Correlation analysis of salt tolerance and grain quality traits in Indica rice (Oryza sativa L.). Acta Agric Nucl Sin,2013,27(12):1938-1947. (in Chinese with English abstract)

    [11] 郭望模,傅亞萍,孫宗修,鄭鎮(zhèn)一. 鹽脅迫下不同水稻種質(zhì)形態(tài)指標(biāo)與耐鹽性的相關(guān)分析. 植物遺傳資源學(xué)報(bào),2003,4(3): 245-251.Guo W M,Fu Y P,Sun Z X,Zheng Z Y. The Correlation analysis between the morphological indices and salt tolerance in different rice germplasm under the salt stress.J Plant Genet Resour,2003,4(3): 245-251. (in Chinese with English abstract)

    [12] 賈寶艷,周嬋嬋,孫曉雪,董立強(qiáng),黃元財(cái),王巖,王術(shù). 遼寧省水稻種質(zhì)資源的耐鹽性鑒定評(píng)價(jià). 作物雜志,2013,4: 57-62.Jia B Y,Zhou C C,Sun X X,Dong L Q,Hang Y C,Wang Y,Wang S. The evaluation of salt tolerance ability of rice varieties in Liaoning Province. Crops,2013,4: 57-62. (in Chinese with English abstract)

    [13] Wang Z F,Wang J F,Bao Y M,Wu Y Y,Su X,Zang H S.Inheritance of rice seed germination ability under salt stress. Rice Sci,2010,17(2): 105-110.

    [14] Chen D F,Li Y L,Fang T,Shi X L,Chen X W. Specific roles of tocopherols and tocotrienols in seed longevity and germination tolerance to abiotic stress in transgenic rice. Plant Sci,2016,244: 31-39.

    [15] Hasthanasombut S,Supaibulwatana K,Mii M,Nakamura I. Genetic manipulation of japonica rice using the OsBADH1 gene from Indica rice to improve salinity tolerance. Plant Cell Tissue Organ Cult,2011,104(1):79-89.

    [16] 馮鐘慧,劉曉龍,姜昌杰,梁正偉. 吉林省粳稻種質(zhì)萌發(fā)期耐堿性和耐鹽性綜合評(píng)價(jià). 土壤與作物,2016,5(2):120-127.Feng Z H,Liu X L,Jiang C J,Liang Z W.Comprehensive evaluation of rice(Oryza sativa japonica)germplasm for alkaline /saline tolerance at germination stage from Jilin province,China. Soil Crop,2016,5(2):120-127. (in Chinese with English abstract)

    [17] 朱麗偉,曹棟棟,付玉營(yíng),胡琦娟,利站,關(guān)亞靜,胡偉民,胡晉. 可溶性寡糖和小分子的熱激蛋白與雜交水稻種子成熟過程中發(fā)芽能力及種子活力相關(guān). 作物學(xué)報(bào),2016,42(5): 714-724.Zhu L W,Cao D D,Fu Y Y,Hu Q J,Li Z,Guan Y J,Hu W M,Hu J. Soluble oligosaccharide and small heat shock protein correlated with seed germination and vigor during hybrid rice seed maturation. Acta Agron Sin,2016,42(5):714-724. (in Chinese with English abstract)

    [18] 魚小軍,肖紅,徐長(zhǎng)林,景媛媛,柴成燕. 扁蓿豆和苜蓿種子萌發(fā)期抗旱性和耐鹽性比較. 植物遺傳資源學(xué)報(bào),2015,16(2): 405-410.Yu X J,Xiao H,Xu C L,Jing Y Y,Chai C Y.Comparative study on drought resistance and salt tolerance of Medicago ruthenica and Medicago varia at seed germination period. J Plant Genet Resour,2015,16(2): 405-410. (in Chinese with English abstract)

    [19] 李振,廖同慶,馮青春,張東彥,王強(qiáng),王曉楠,張琳潔,王秀. 基于圖像處理技術(shù)的黃瓜種子活力指數(shù)檢測(cè)系統(tǒng)設(shè)計(jì). 種子,2015,34(6): 111-115.Li Z,Liao T Q,Feng Q C,Zhang D Y,Wang Q,Wang X N,Zhang L J,Wang X. A system design on cucumber seed vigor index detection on based on image processing.Seed,2015,34(6): 111-115. (in Chinese with English abstract)

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

    [21] 齊國(guó)昌,余春磊,張鹍飛,羅小嬌,馮宗云. 青稞苗期耐鹽性鑒定及評(píng)價(jià). 麥類作物學(xué)報(bào),2014,37(7):950-956.Qi G C,Yu C L,Zhang K F,Luo X J,Feng Z Y.Identification and assessment on salt resistance of hulless barley at seedling stage. J Triticeae Crops,2014,37(7):950-956. (in Chinese with English abstract)

    [22] 管志勇,陳發(fā)棣,滕年軍,陳素梅,劉浦生. 5種菊花近緣種屬植物的耐鹽性比較. 中國(guó)農(nóng)業(yè)科學(xué),2010,43(4):784-787.Guan Z Y,Chen F D,Teng N J,Chen S M,Liu P S. Study on the NaCl tolerance in five plant species from Dendranthema and its relatives. Sci Agric Sin,2010,43(4):784-787. (in Chinese with English abstract)

    [23] 田蕾,王娜,張雪艷,楊斌林,孫佳瑩,李敏,李培富.鹽脅迫下不同粳稻品種的形態(tài)和生理特性. 廣東農(nóng)業(yè)科學(xué),2014,41(23): 1-6.Tian L,Wang N,Zhang X Y,Yang B L,Sun J Y,Li M,Li P F. Morphological and physiological characteristics of different japonica rice varieties under salt stress.Guangdong Agric Sci,2014,41(23): 1-6. (in Chinese with English abstract)

    [24] 王娜,陳亞萍,田蕾,張得雯,王瑞智,楊苗,李培富.粳稻種質(zhì)資源苗期根系形態(tài)特征與耐鹽性相關(guān)分析.廣東農(nóng)業(yè)科學(xué),2015,42(10): 1-10.Wang N,Chen Y P,Tian L,Zhang D W,Wang R Z,Yang M,Li P F. Correlation between root morphological characteristics of japonica rice germplasm and salt tolerance at seedling stage. Guangdong Agric Sci,2015,42(10): 1-10. (in Chinese with English abstract)

    [25] Averinaa N G,Gritskevicha E R,Vershilovskayaa I V,Usatov A V,Yaronskaya E B. Mechanisms of salt stress tolerance development in barley plants under the Influence of 5-aminolevulinic acid. Russ J Plant Physiol-Engl Tr,2010,57(6): 792-798.

    [26] 崔江慧,謝登磊,常金華. 高粱材料耐鹽性綜合評(píng)價(jià)方法的初步建立與驗(yàn)證. 植物遺傳資源學(xué)報(bào),2012,13(1):35-41.Cui J H,Xie D L,Chang J H. Establishment and verification of comprehensive evaluation method for salt tolerance of sorghum materials. J Plant Genet Resour,2012,13(1): 35-41. (in Chinese with English abstract)

    [27] 李培英,孫宗玖. 33 份偃麥草種質(zhì)芽期耐鹽性評(píng)價(jià).草業(yè)科學(xué),2015,32(4): 593-600.Li P Y,Sun Z J. Evaluation on the salt resistance of germplasm resources of 33 Elytrigria repens during seed germination period. Pratac Sci,2015,32(4): 593-600. (in Chinese with English abstract)

    [28] 何曉蘭,徐照龍,張大勇,黃益洪,彭陳,邵宏波,王為,郭士偉. 65個(gè)高粱種質(zhì)萌芽期的耐鹽指標(biāo)比較及其耐鹽性綜合評(píng)價(jià). 植物資源與環(huán)境學(xué)報(bào),2015,24(4):52-60.He X L,Xu Z L,Zhang D Y,Huang Y H,Peng C,Shao H B,Wang W,Guo S W. Comparison on salt tolerance indexes of 65 germplasms of Sorghum bicolor at germination stage and comprehensive evaluation on their salt tolerance. J Plant Resour Environ,2015,24(4): 52-60.(in Chinese with English abstract)

    [29] 王萌萌,姜奇彥,胡正,張輝,樊守金,馮瀝,張海玲.小麥品種資源耐鹽性鑒定. 植物遺傳資源學(xué)報(bào),2012,13(2): 189-194.Wang M M,Jiang Q Y,Hu Z,Zhang H,Fan S J,Feng L,Zhang H L. Evaluation for salt tolerance of wheat cultivars. J Plant Genet Resour,2012,13(2): 189-194. (in Chinese with English abstract)

    [30] 吳家富,楊博文,向珣朝,許亮,顏李梅. 不同水稻種質(zhì)在不同生育期耐鹽鑒定的差異. 植物學(xué)報(bào),2017,52(1): 77-88.Wu J F,Yang B W,Xiang X C,Xu L,Yan L M.Identification of salt tolerance in different rice germplasm at different growth stages. Chin Bull Bot,2017,52 (1):77-88. (in Chinese with English abstract)

    [31] 張銀平,杜瑞成,刁培松,楊善東. 山東省水稻免耕旱直播試驗(yàn)及可行性分析. 農(nóng)業(yè)工程學(xué)報(bào),2016,32(12):24-30.Zhang Y P,Du R C,Diao P S,Yang S D. Experiment of no-tillage and drought direct sowing rice and feasibility analysis in Shandong Province. Trans Chin Soc Agric Eng,2016,32(12): 24-30. (in Chinese with English abstract)

    [32] Li Z K,Fu B Y,Gao Y M,Xu J L,Ali J ,Lafitte H R,Jiang Y Z,Rey J D,Vijayakumar C H M,Maghirang R,Zheng T Q,Zhu L H. Genome-wide introgression lines and their use in genetic and molecular dissection of complex phenotypes in rice (Oryza sativa L.). Plant Mol Biol,2005,59(1): 33-52.

    Comprehensive Evaluation and Selection of Rice (Oryza sativa japonica) Germplasm for Saline Tolerance at Germination Stage

    TIAN Lei,CHEN Yaping,LIU Jun,MA Xiaogang,WANG Na,YANG Bing,LI Ying,GUO Haidong,LI Juan,HU Hui,ZHANG Yinxia,LI Peifu*
    (College of Agriculture,Ningxia University,Yinchuan 750021,China; *Corresponding author,E-mail: peifuli@163. com)

    【Objective】Soil salinity is one of the major abiotic stresses affecting rice growth and production. Salt resistance identification at germination stage,indices screening,and cultivar breeding of rice germplasm resources are important for rice production.【Method】Seeds of sixty-four japonica rice germplasm accessions germinated in the 125 mmol L-1NaCl solution. On the third day after salt stress,germination number of the seeds was measured. After salt stress for five and ten days,germination number,shoot length and root length were measured and then the relative shoot length,relative root length,relative germination potential,relative germination rate,relative salt damage rate,germination index and vigor index were calculated. Multiple statistical methods were used to comprehensively evaluate the salt tolerance of 64 rice japonica germplasm and the germination characteristics of typical salt-tolerant and salt-sensitive rice germplasm under salt stress were analyzed. 【Result】The relative salt damage rate was significantly negatively correlated with relative root length,relative germination potential,relative germination rate,germination index and vigor index.Except for relative shoot length,these six indexes were significantly correlated with each other. Through cluster analysis,64 japonica rice germplasm were divided into four groups. The group Ⅰ and Ⅳwere typical salt-sensitive and salt-tolerant germplasm,respectively. The group Ⅱ was mainly composed of low salt-resistant rice and the group Ⅲ was moderately salt-tolerant. With the principal component analysis these seven single indexes could be converted to three principal components. An objective comprehensive evaluation value (D value) of salt tolerance was obtained using membership function with index weight method. The five highest and lowest D value germplasm resources were selected to carry on significance analysis. The results showed that except for relative shoot length,all of these seven indexes measured after stress for five and ten days were significantly different from the two types. 【Conclusion】Rice was sensitive to salt stress at germination stage and there were significant differences between different germplasm. Using stepwise regression and principal component analysis,germination index,relative root length and relative salt damage rate were selected as the key indexes for the rapid identification of salt tolerance of rice germplasm. If multivariate statistical method was used,the reliability will be improved.

    rice; germplasm resource; salt tolerance at germination stage; principal component analysis; comprehensive evaluation

    10.16819/j.1001-7216.2017.6168

    2016-12-20;修改稿收到日期:2017-04-29。

    國(guó)家自然科學(xué)基金資助項(xiàng)目(31401361);寧夏農(nóng)業(yè)育種專項(xiàng)(2013NYYZ0302)。

    Q945.78; S511.2+2

    A

    1001-7216(2017)06-0631-12

    猜你喜歡
    芽長(zhǎng)耐鹽耐鹽性
    有了這種合成酶 水稻可以耐鹽了
    郁金香耐鹽性鑒定方法研究
    玉米種子芽大小對(duì)農(nóng)桿菌轉(zhuǎn)化的影響
    播種不同芽長(zhǎng)種薯芽塊對(duì)早春錯(cuò)期馬鈴薯生育進(jìn)程及產(chǎn)量的影響
    不同芽長(zhǎng)對(duì)馬鈴薯理化指標(biāo)的影響
    源于大麥小孢子突變體的苗期耐鹽性研究
    三個(gè)大豆品種萌發(fā)期和苗期的耐鹽性比較
    大豆發(fā)芽期和苗期耐鹽性的隸屬函數(shù)分析
    耐鹽保水劑的合成及其性能
    耐鹽高降解蛋白菌株的分離鑒定及其降解條件的研究
    亚洲av日韩精品久久久久久密| 中文字幕av成人在线电影| a级一级毛片免费在线观看| 亚洲欧美激情综合另类| 亚洲aⅴ乱码一区二区在线播放| 黄色视频,在线免费观看| 午夜免费观看网址| 老汉色∧v一级毛片| 亚洲国产日韩欧美精品在线观看 | 舔av片在线| 久久亚洲精品不卡| 高潮久久久久久久久久久不卡| 亚洲精品久久国产高清桃花| 欧美性猛交黑人性爽| 精品国产美女av久久久久小说| 国产真人三级小视频在线观看| 在线天堂最新版资源| 国产熟女xx| a级一级毛片免费在线观看| 国产精品一区二区免费欧美| 女人被狂操c到高潮| 国产精品美女特级片免费视频播放器| 桃红色精品国产亚洲av| 国产精品自产拍在线观看55亚洲| 麻豆国产97在线/欧美| 国产真人三级小视频在线观看| av中文乱码字幕在线| 人人妻,人人澡人人爽秒播| 小说图片视频综合网站| 国产免费男女视频| 久久精品国产亚洲av涩爱 | 日韩av在线大香蕉| 禁无遮挡网站| 18美女黄网站色大片免费观看| a级毛片a级免费在线| 在线观看一区二区三区| av中文乱码字幕在线| 桃红色精品国产亚洲av| 亚洲av第一区精品v没综合| 国产一区二区激情短视频| 宅男免费午夜| www日本黄色视频网| 伊人久久大香线蕉亚洲五| 国产极品精品免费视频能看的| 美女被艹到高潮喷水动态| 99热这里只有精品一区| 午夜激情福利司机影院| 亚洲欧美日韩高清专用| 给我免费播放毛片高清在线观看| 51午夜福利影视在线观看| a级毛片a级免费在线| 欧美黑人欧美精品刺激| 国产一区在线观看成人免费| 波多野结衣高清无吗| 久久国产乱子伦精品免费另类| 搡老熟女国产l中国老女人| 久久久久久久亚洲中文字幕 | 一夜夜www| 国产一区二区三区在线臀色熟女| 亚洲精品色激情综合| 91av网一区二区| 啦啦啦免费观看视频1| 久久欧美精品欧美久久欧美| 我的老师免费观看完整版| 黄色成人免费大全| 又黄又粗又硬又大视频| 国产毛片a区久久久久| 少妇高潮的动态图| 成人国产综合亚洲| 无限看片的www在线观看| 丁香六月欧美| 免费在线观看日本一区| 又爽又黄无遮挡网站| 日本 欧美在线| 18美女黄网站色大片免费观看| а√天堂www在线а√下载| 91在线精品国自产拍蜜月 | 欧美中文日本在线观看视频| 欧美日韩一级在线毛片| 国产精品 国内视频| 69人妻影院| 日本免费a在线| 久久精品影院6| 他把我摸到了高潮在线观看| 美女免费视频网站| 国产三级在线视频| 精品久久久久久,| 国产野战对白在线观看| 欧美又色又爽又黄视频| 欧美+亚洲+日韩+国产| 亚洲 国产 在线| 又粗又爽又猛毛片免费看| 又粗又爽又猛毛片免费看| av在线蜜桃| 日韩欧美 国产精品| 亚洲成人精品中文字幕电影| 草草在线视频免费看| 国产精品av视频在线免费观看| 国产三级在线视频| 国产真实伦视频高清在线观看 | 亚洲国产色片| 免费看美女性在线毛片视频| 两个人视频免费观看高清| 超碰av人人做人人爽久久 | 女警被强在线播放| 成人鲁丝片一二三区免费| x7x7x7水蜜桃| 精品熟女少妇八av免费久了| 男人和女人高潮做爰伦理| 免费人成视频x8x8入口观看| 国产精品久久视频播放| 少妇裸体淫交视频免费看高清| 午夜亚洲福利在线播放| 偷拍熟女少妇极品色| 天天添夜夜摸| 亚洲av熟女| 国产精品 国内视频| 亚洲国产精品成人综合色| 狠狠狠狠99中文字幕| 在线观看午夜福利视频| 国产精品嫩草影院av在线观看 | 99热这里只有是精品50| 韩国av一区二区三区四区| 国产成+人综合+亚洲专区| 99久久99久久久精品蜜桃| 成人精品一区二区免费| 亚洲国产日韩欧美精品在线观看 | 成人特级av手机在线观看| 一二三四社区在线视频社区8| 一二三四社区在线视频社区8| 久久久久久久午夜电影| 天天添夜夜摸| 深夜精品福利| 天天添夜夜摸| 搡老妇女老女人老熟妇| 免费搜索国产男女视频| 日本 欧美在线| 一进一出抽搐动态| 精品熟女少妇八av免费久了| 欧美3d第一页| 国产69精品久久久久777片| 国产亚洲欧美98| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 丰满人妻一区二区三区视频av | 色综合欧美亚洲国产小说| 又爽又黄无遮挡网站| 一本精品99久久精品77| 午夜精品一区二区三区免费看| 人妻丰满熟妇av一区二区三区| 午夜免费男女啪啪视频观看 | 中文资源天堂在线| 国语自产精品视频在线第100页| 国产精品1区2区在线观看.| 三级国产精品欧美在线观看| 免费在线观看影片大全网站| 国产高潮美女av| www日本黄色视频网| 欧美bdsm另类| 国产伦精品一区二区三区视频9 | 91久久精品国产一区二区成人 | 亚洲第一电影网av| av专区在线播放| 亚洲中文日韩欧美视频| 午夜日韩欧美国产| 国产精品国产高清国产av| 欧美不卡视频在线免费观看| 午夜福利在线观看免费完整高清在 | 伊人久久大香线蕉亚洲五| 精品国产亚洲在线| 国产免费男女视频| 激情在线观看视频在线高清| 麻豆成人av在线观看| 身体一侧抽搐| 人人妻人人澡欧美一区二区| 国产精品香港三级国产av潘金莲| 欧美午夜高清在线| 久久久成人免费电影| 亚洲,欧美精品.| 国产av不卡久久| 国产一区二区亚洲精品在线观看| 69av精品久久久久久| 亚洲精品久久国产高清桃花| 亚洲一区高清亚洲精品| 成人永久免费在线观看视频| 亚洲av成人不卡在线观看播放网| 99久久精品热视频| 啦啦啦免费观看视频1| 午夜福利在线在线| 动漫黄色视频在线观看| 婷婷精品国产亚洲av| av在线蜜桃| www.色视频.com| 亚洲精品影视一区二区三区av| 真人一进一出gif抽搐免费| 国产日本99.免费观看| 国产色爽女视频免费观看| 精品久久久久久久人妻蜜臀av| 熟妇人妻久久中文字幕3abv| 成人特级黄色片久久久久久久| 日本黄大片高清| 人人妻人人澡欧美一区二区| 国产精品久久久久久精品电影| 男女之事视频高清在线观看| 三级毛片av免费| 久久久久久久亚洲中文字幕 | 国产精品 欧美亚洲| 久久久久久久久中文| 日本成人三级电影网站| 国产一区在线观看成人免费| 综合色av麻豆| 国产精品精品国产色婷婷| 两个人视频免费观看高清| 久久精品91无色码中文字幕| av女优亚洲男人天堂| 少妇的丰满在线观看| 一本精品99久久精品77| 日日夜夜操网爽| 在线十欧美十亚洲十日本专区| 黄色成人免费大全| 久久久久亚洲av毛片大全| 国产免费av片在线观看野外av| 国产高清视频在线观看网站| 嫩草影院精品99| 中国美女看黄片| 欧美日韩一级在线毛片| 国产久久久一区二区三区| 天堂av国产一区二区熟女人妻| 亚洲av五月六月丁香网| 精品电影一区二区在线| 在线观看一区二区三区| 有码 亚洲区| 日本 欧美在线| 日本撒尿小便嘘嘘汇集6| 精品人妻一区二区三区麻豆 | 婷婷亚洲欧美| 午夜福利18| 午夜亚洲福利在线播放| 色av中文字幕| 波野结衣二区三区在线 | 国产淫片久久久久久久久 | 操出白浆在线播放| 成人欧美大片| 国产成人影院久久av| 国产中年淑女户外野战色| 日日夜夜操网爽| 此物有八面人人有两片| 日本a在线网址| 欧美中文综合在线视频| 免费av不卡在线播放| 国产精品香港三级国产av潘金莲| 欧美日韩乱码在线| 老司机午夜十八禁免费视频| www.熟女人妻精品国产| 香蕉久久夜色| 每晚都被弄得嗷嗷叫到高潮| 国产精品免费一区二区三区在线| 免费人成视频x8x8入口观看| 色视频www国产| 制服人妻中文乱码| 国产精品99久久久久久久久| 亚洲成人免费电影在线观看| 国产野战对白在线观看| 亚洲18禁久久av| 国产黄色小视频在线观看| 麻豆国产av国片精品| 老司机福利观看| 国产伦精品一区二区三区四那| 波野结衣二区三区在线 | 老汉色av国产亚洲站长工具| 99国产精品一区二区蜜桃av| 国产精品美女特级片免费视频播放器| 久久伊人香网站| 中文资源天堂在线| 不卡一级毛片| 久久精品夜夜夜夜夜久久蜜豆| 国产精品爽爽va在线观看网站| 久久九九热精品免费| 色av中文字幕| 韩国av一区二区三区四区| 身体一侧抽搐| 99久久精品热视频| 国内毛片毛片毛片毛片毛片| 给我免费播放毛片高清在线观看| 成熟少妇高潮喷水视频| 一级毛片高清免费大全| 九九在线视频观看精品| 成人三级黄色视频| 琪琪午夜伦伦电影理论片6080| 一区二区三区免费毛片| 国产亚洲欧美98| 丁香六月欧美| 亚洲 欧美 日韩 在线 免费| 色噜噜av男人的天堂激情| 18禁在线播放成人免费| 久久性视频一级片| 亚洲国产精品久久男人天堂| 三级国产精品欧美在线观看| 色吧在线观看| 欧美性感艳星| 一进一出抽搐动态| 久久久精品大字幕| 女同久久另类99精品国产91| 五月伊人婷婷丁香| 欧美一区二区亚洲| 午夜激情福利司机影院| 国产精品一区二区免费欧美| 国产精品久久久久久久电影 | 99久久综合精品五月天人人| 欧美区成人在线视频| 有码 亚洲区| 亚洲av不卡在线观看| av在线蜜桃| 亚洲av熟女| 99热这里只有是精品50| 国产精品1区2区在线观看.| eeuss影院久久| 99国产综合亚洲精品| 三级毛片av免费| 99热精品在线国产| 亚洲电影在线观看av| 久久这里只有精品中国| 天美传媒精品一区二区| 又爽又黄无遮挡网站| 可以在线观看毛片的网站| 国产午夜福利久久久久久| 午夜免费观看网址| 一边摸一边抽搐一进一小说| 亚洲av不卡在线观看| 熟妇人妻久久中文字幕3abv| 久久久久亚洲av毛片大全| 有码 亚洲区| 亚洲欧美日韩高清专用| 最新美女视频免费是黄的| 婷婷精品国产亚洲av| 男人舔女人下体高潮全视频| 午夜福利免费观看在线| 好男人在线观看高清免费视频| 国产精品一区二区免费欧美| 国产精品亚洲一级av第二区| 3wmmmm亚洲av在线观看| 18禁美女被吸乳视频| 亚洲一区高清亚洲精品| 久久精品亚洲精品国产色婷小说| 波多野结衣高清无吗| 欧美xxxx黑人xx丫x性爽| 99久国产av精品| 亚洲最大成人中文| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 最新在线观看一区二区三区| 亚洲国产日韩欧美精品在线观看 | 欧美黑人巨大hd| 99热这里只有精品一区| 国产成+人综合+亚洲专区| 午夜福利在线在线| 噜噜噜噜噜久久久久久91| 18禁黄网站禁片午夜丰满| 免费人成在线观看视频色| av天堂在线播放| 国产亚洲av嫩草精品影院| 麻豆一二三区av精品| 91av网一区二区| 欧美日韩一级在线毛片| 内射极品少妇av片p| 美女高潮喷水抽搐中文字幕| 欧美黑人巨大hd| 黄色丝袜av网址大全| 亚洲熟妇中文字幕五十中出| 中文字幕人妻丝袜一区二区| 欧美乱妇无乱码| 国产成人av激情在线播放| 亚洲专区中文字幕在线| 18禁在线播放成人免费| 91麻豆av在线| 亚洲专区国产一区二区| 国产私拍福利视频在线观看| 欧美乱妇无乱码| 啦啦啦观看免费观看视频高清| 亚洲精品一区av在线观看| 成人无遮挡网站| 观看免费一级毛片| 亚洲色图av天堂| 久久亚洲精品不卡| 久久久久国内视频| 精品久久久久久久久久久久久| 国产一区二区三区视频了| 久久久精品欧美日韩精品| 99久久综合精品五月天人人| 国产精品亚洲一级av第二区| 午夜福利免费观看在线| 日本撒尿小便嘘嘘汇集6| 又黄又粗又硬又大视频| 久久精品夜夜夜夜夜久久蜜豆| 毛片女人毛片| 18+在线观看网站| 欧美av亚洲av综合av国产av| 啪啪无遮挡十八禁网站| 天堂影院成人在线观看| 又粗又爽又猛毛片免费看| 母亲3免费完整高清在线观看| 尤物成人国产欧美一区二区三区| 久久香蕉精品热| 淫秽高清视频在线观看| www国产在线视频色| 日本熟妇午夜| 欧美黄色片欧美黄色片| 国产成人aa在线观看| 国产亚洲精品一区二区www| 日韩欧美精品v在线| 18禁美女被吸乳视频| 内地一区二区视频在线| 成年女人看的毛片在线观看| 国产精品精品国产色婷婷| 亚洲男人的天堂狠狠| 我的老师免费观看完整版| 国产高清有码在线观看视频| 99热这里只有是精品50| 日韩免费av在线播放| 免费看美女性在线毛片视频| 国产成人a区在线观看| 亚洲自拍偷在线| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 欧美在线黄色| 国产真实乱freesex| 搡老熟女国产l中国老女人| 九九久久精品国产亚洲av麻豆| 国产成人aa在线观看| 欧美绝顶高潮抽搐喷水| 国产成年人精品一区二区| 日本黄色片子视频| 国产99白浆流出| 男女之事视频高清在线观看| 真人做人爱边吃奶动态| 国产91精品成人一区二区三区| 一进一出抽搐动态| 成人国产综合亚洲| 九色成人免费人妻av| 最新在线观看一区二区三区| 夜夜看夜夜爽夜夜摸| 俄罗斯特黄特色一大片| 夜夜夜夜夜久久久久| 亚洲专区国产一区二区| 亚洲av第一区精品v没综合| 亚洲av电影不卡..在线观看| 国产高潮美女av| 中国美女看黄片| 免费av观看视频| 丰满人妻一区二区三区视频av | 老司机午夜福利在线观看视频| 亚洲 欧美 日韩 在线 免费| 欧美一区二区亚洲| 欧美大码av| 久久人人精品亚洲av| 亚洲av免费高清在线观看| 亚洲色图av天堂| 国产亚洲欧美在线一区二区| 欧美成人免费av一区二区三区| 国产伦人伦偷精品视频| 国产欧美日韩精品亚洲av| 国产v大片淫在线免费观看| 国内精品美女久久久久久| 超碰av人人做人人爽久久 | 亚洲 欧美 日韩 在线 免费| 精品久久久久久久毛片微露脸| 脱女人内裤的视频| 婷婷六月久久综合丁香| 热99re8久久精品国产| 国产视频一区二区在线看| 免费大片18禁| 国产精品野战在线观看| 亚洲激情在线av| 国产精品久久久久久久久免 | 国产欧美日韩精品一区二区| eeuss影院久久| 在线观看日韩欧美| 日韩中文字幕欧美一区二区| 超碰av人人做人人爽久久 | 亚洲色图av天堂| 国产三级黄色录像| 久久久久精品国产欧美久久久| 波多野结衣高清无吗| 日韩欧美精品免费久久 | 国产激情欧美一区二区| 欧美bdsm另类| 国产美女午夜福利| 亚洲专区国产一区二区| 久久久精品欧美日韩精品| 国产aⅴ精品一区二区三区波| 三级男女做爰猛烈吃奶摸视频| 19禁男女啪啪无遮挡网站| 亚洲av不卡在线观看| 亚洲专区中文字幕在线| 亚洲人成电影免费在线| 一区二区三区激情视频| 女人十人毛片免费观看3o分钟| 欧美高清成人免费视频www| 高清日韩中文字幕在线| 老司机午夜福利在线观看视频| 91久久精品电影网| 欧美在线一区亚洲| 少妇的逼水好多| 美女黄网站色视频| 亚洲久久久久久中文字幕| 日本与韩国留学比较| 亚洲最大成人手机在线| www日本在线高清视频| 麻豆国产97在线/欧美| 色尼玛亚洲综合影院| 亚洲精品在线观看二区| 男人的好看免费观看在线视频| 无限看片的www在线观看| 19禁男女啪啪无遮挡网站| 日本三级黄在线观看| 色噜噜av男人的天堂激情| 国产日本99.免费观看| 亚洲一区二区三区色噜噜| 69av精品久久久久久| 日韩欧美在线乱码| 国产高清视频在线观看网站| 天堂影院成人在线观看| 一进一出好大好爽视频| 一卡2卡三卡四卡精品乱码亚洲| 久久欧美精品欧美久久欧美| 成年人黄色毛片网站| 国产精品1区2区在线观看.| 国产精品精品国产色婷婷| 天堂av国产一区二区熟女人妻| 精品国产美女av久久久久小说| 成年女人毛片免费观看观看9| 伊人久久大香线蕉亚洲五| 午夜免费激情av| 成人国产综合亚洲| 在线观看美女被高潮喷水网站 | 俄罗斯特黄特色一大片| 国产av麻豆久久久久久久| 国产私拍福利视频在线观看| 亚洲av电影不卡..在线观看| 午夜免费成人在线视频| av女优亚洲男人天堂| 亚洲 欧美 日韩 在线 免费| 免费人成视频x8x8入口观看| 国产精品亚洲美女久久久| 国产一区二区亚洲精品在线观看| 久久久久精品国产欧美久久久| 精品国产三级普通话版| 国产单亲对白刺激| 午夜福利在线在线| 一级作爱视频免费观看| 亚洲国产高清在线一区二区三| 国产高清有码在线观看视频| 中文资源天堂在线| 亚洲熟妇熟女久久| 啦啦啦观看免费观看视频高清| 波多野结衣高清作品| 男女那种视频在线观看| 午夜福利18| 欧美中文日本在线观看视频| 日本免费一区二区三区高清不卡| 成人av在线播放网站| 午夜福利在线观看免费完整高清在 | 欧美3d第一页| 午夜影院日韩av| 狂野欧美激情性xxxx| 亚洲成人中文字幕在线播放| 午夜福利在线观看吧| 亚洲av中文字字幕乱码综合| 香蕉av资源在线| 国产97色在线日韩免费| 亚洲电影在线观看av| 桃红色精品国产亚洲av| 一区二区三区激情视频| 性色avwww在线观看| 麻豆成人av在线观看| 日本 欧美在线| 麻豆成人av在线观看| 99久久无色码亚洲精品果冻| 一本综合久久免费| 国产av不卡久久| 小蜜桃在线观看免费完整版高清| 亚洲性夜色夜夜综合| 欧美高清成人免费视频www| 欧美性猛交╳xxx乱大交人| netflix在线观看网站| 亚洲性夜色夜夜综合| 欧美日本视频| 国产精品自产拍在线观看55亚洲| 9191精品国产免费久久| 噜噜噜噜噜久久久久久91| 日日摸夜夜添夜夜添小说| 国产三级黄色录像| 色综合婷婷激情| 最近视频中文字幕2019在线8| 亚洲av中文字字幕乱码综合| 国产精品日韩av在线免费观看| 欧美最黄视频在线播放免费| 亚洲av成人av| xxxwww97欧美| 欧美又色又爽又黄视频| 特大巨黑吊av在线直播| 精品乱码久久久久久99久播| 啦啦啦韩国在线观看视频| 国产三级黄色录像| 日本熟妇午夜| 国产亚洲精品久久久com| 亚洲精品在线观看二区| 精品99又大又爽又粗少妇毛片 | 黄色片一级片一级黄色片| 亚洲av成人不卡在线观看播放网| 美女 人体艺术 gogo| 欧美日本亚洲视频在线播放| 国产成人欧美在线观看| 国产精品av视频在线免费观看| www.熟女人妻精品国产| 黄片小视频在线播放|