• <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ù)分析
    耐鹽保水劑的合成及其性能
    耐鹽高降解蛋白菌株的分離鑒定及其降解條件的研究
    免费在线观看影片大全网站| 日本黄色视频三级网站网址| 国产精品一区二区精品视频观看| 成人国产综合亚洲| 在线观看免费视频日本深夜| 久久婷婷成人综合色麻豆| 亚洲一区二区三区色噜噜| 12—13女人毛片做爰片一| 亚洲 国产 在线| 伊人久久大香线蕉亚洲五| 久久国产精品男人的天堂亚洲| 在线永久观看黄色视频| 亚洲国产看品久久| 日韩精品中文字幕看吧| 久久人人爽av亚洲精品天堂| 久久草成人影院| 色综合欧美亚洲国产小说| 欧美一级毛片孕妇| 丰满人妻熟妇乱又伦精品不卡| 欧美大码av| 精品一区二区三区av网在线观看| 九色国产91popny在线| 日本三级黄在线观看| 69精品国产乱码久久久| 亚洲精品久久国产高清桃花| 天堂动漫精品| 亚洲第一av免费看| 午夜福利欧美成人| 亚洲av美国av| 黄片小视频在线播放| 国产成人系列免费观看| 国产三级在线视频| 国产精品久久久av美女十八| 国产精品影院久久| 成人永久免费在线观看视频| 99国产精品99久久久久| 99久久综合精品五月天人人| 日本三级黄在线观看| 国产片内射在线| 欧美国产日韩亚洲一区| 久久久久久亚洲精品国产蜜桃av| 久久精品国产99精品国产亚洲性色 | 波多野结衣高清无吗| 老司机靠b影院| 最好的美女福利视频网| 一区二区日韩欧美中文字幕| АⅤ资源中文在线天堂| 欧美激情极品国产一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 不卡一级毛片| 欧美+亚洲+日韩+国产| 婷婷六月久久综合丁香| 国产精品一区二区三区四区久久 | 美女高潮喷水抽搐中文字幕| or卡值多少钱| 给我免费播放毛片高清在线观看| 好男人在线观看高清免费视频 | 性欧美人与动物交配| 久久久久国产一级毛片高清牌| 狠狠狠狠99中文字幕| 中亚洲国语对白在线视频| 丝袜人妻中文字幕| 日韩一卡2卡3卡4卡2021年| 亚洲国产中文字幕在线视频| 久久影院123| 99国产精品99久久久久| 日韩精品中文字幕看吧| 午夜两性在线视频| 久热爱精品视频在线9| 亚洲精品一卡2卡三卡4卡5卡| 欧美精品亚洲一区二区| 亚洲最大成人中文| 国产免费男女视频| 一二三四在线观看免费中文在| 大型黄色视频在线免费观看| 国产精品免费一区二区三区在线| 一个人观看的视频www高清免费观看 | 亚洲国产精品999在线| 亚洲精品在线观看二区| 少妇熟女aⅴ在线视频| 国产亚洲av高清不卡| 亚洲国产毛片av蜜桃av| 男女床上黄色一级片免费看| 黄色视频,在线免费观看| 色精品久久人妻99蜜桃| 18禁黄网站禁片午夜丰满| 亚洲国产中文字幕在线视频| 国产日韩一区二区三区精品不卡| 亚洲国产欧美网| 久久久久国产精品人妻aⅴ院| 欧美日韩黄片免| 欧美日本亚洲视频在线播放| 国产在线观看jvid| 国产成年人精品一区二区| 国产亚洲欧美在线一区二区| 免费观看人在逋| 国产av又大| 久久草成人影院| 免费观看精品视频网站| 少妇熟女aⅴ在线视频| 999久久久国产精品视频| 久久亚洲精品不卡| 国产av又大| 热re99久久国产66热| 99国产极品粉嫩在线观看| 欧美日韩亚洲综合一区二区三区_| av片东京热男人的天堂| 欧美 亚洲 国产 日韩一| 国产xxxxx性猛交| 亚洲精华国产精华精| 精品第一国产精品| 日本五十路高清| 少妇裸体淫交视频免费看高清 | 人人妻人人爽人人添夜夜欢视频| 精品无人区乱码1区二区| av电影中文网址| 美女午夜性视频免费| cao死你这个sao货| 可以免费在线观看a视频的电影网站| 亚洲国产中文字幕在线视频| 电影成人av| 亚洲一码二码三码区别大吗| 国产精品 国内视频| 18禁黄网站禁片午夜丰满| 国产精品久久电影中文字幕| 99国产精品99久久久久| 免费女性裸体啪啪无遮挡网站| 黄频高清免费视频| 18禁黄网站禁片午夜丰满| 国内精品久久久久精免费| 精品不卡国产一区二区三区| 国产av在哪里看| 一二三四社区在线视频社区8| 黄频高清免费视频| 大香蕉久久成人网| 精品一区二区三区四区五区乱码| 国产人伦9x9x在线观看| 欧美性长视频在线观看| 99国产极品粉嫩在线观看| 最近最新中文字幕大全免费视频| 亚洲欧美日韩高清在线视频| 神马国产精品三级电影在线观看 | 国产高清激情床上av| 久久性视频一级片| 黄色a级毛片大全视频| 黄色成人免费大全| 97人妻精品一区二区三区麻豆 | 黑人巨大精品欧美一区二区蜜桃| 少妇被粗大的猛进出69影院| 亚洲av第一区精品v没综合| 国产单亲对白刺激| 久久精品亚洲精品国产色婷小说| 国产野战对白在线观看| 成人亚洲精品一区在线观看| 一级作爱视频免费观看| 午夜久久久久精精品| 日韩国内少妇激情av| 日本欧美视频一区| 少妇裸体淫交视频免费看高清 | 悠悠久久av| 一进一出好大好爽视频| 精品免费久久久久久久清纯| 久久香蕉激情| 正在播放国产对白刺激| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣高清无吗| 在线免费观看的www视频| 满18在线观看网站| 亚洲免费av在线视频| 国产高清videossex| 一本大道久久a久久精品| ponron亚洲| 久久狼人影院| 国产精品自产拍在线观看55亚洲| 别揉我奶头~嗯~啊~动态视频| 国产1区2区3区精品| 嫩草影视91久久| 国产av精品麻豆| 可以免费在线观看a视频的电影网站| 亚洲国产精品合色在线| avwww免费| 日韩高清综合在线| 老司机福利观看| 91在线观看av| 成人国产一区最新在线观看| 99国产精品免费福利视频| 一级,二级,三级黄色视频| 色尼玛亚洲综合影院| 美国免费a级毛片| 亚洲伊人色综图| 亚洲中文av在线| 国产熟女午夜一区二区三区| 国产人伦9x9x在线观看| 欧美国产精品va在线观看不卡| 嫩草影视91久久| 国产精品自产拍在线观看55亚洲| 精品国产一区二区久久| 中文字幕av电影在线播放| 亚洲精品中文字幕在线视频| 亚洲成av片中文字幕在线观看| 亚洲成av人片免费观看| 久久影院123| 欧美精品亚洲一区二区| 美女免费视频网站| 国产精品自产拍在线观看55亚洲| 国产亚洲av嫩草精品影院| 午夜福利成人在线免费观看| 老熟妇乱子伦视频在线观看| 中文字幕人成人乱码亚洲影| 亚洲av熟女| 国产伦人伦偷精品视频| 在线观看免费日韩欧美大片| 老司机午夜十八禁免费视频| 高潮久久久久久久久久久不卡| 男人舔女人的私密视频| 十八禁网站免费在线| 亚洲国产欧美一区二区综合| 两人在一起打扑克的视频| 欧美最黄视频在线播放免费| 亚洲 欧美 日韩 在线 免费| 一a级毛片在线观看| 制服丝袜大香蕉在线| 日韩有码中文字幕| 亚洲男人的天堂狠狠| 国产亚洲精品一区二区www| 国产高清视频在线播放一区| 国产免费av片在线观看野外av| 亚洲无线在线观看| 久久亚洲真实| 性欧美人与动物交配| 日日爽夜夜爽网站| 国产99白浆流出| 午夜免费观看网址| 大型av网站在线播放| 欧美另类亚洲清纯唯美| 国产1区2区3区精品| 精品少妇一区二区三区视频日本电影| 一区福利在线观看| 免费观看人在逋| 操美女的视频在线观看| 在线观看舔阴道视频| 国产亚洲av高清不卡| 一个人免费在线观看的高清视频| 18禁美女被吸乳视频| 97超级碰碰碰精品色视频在线观看| 亚洲精品一区av在线观看| 国产蜜桃级精品一区二区三区| 黄频高清免费视频| 18禁黄网站禁片午夜丰满| 中文字幕久久专区| 一进一出抽搐gif免费好疼| 国产主播在线观看一区二区| 亚洲av熟女| 69精品国产乱码久久久| 国内精品久久久久精免费| 老司机深夜福利视频在线观看| 妹子高潮喷水视频| av视频免费观看在线观看| 满18在线观看网站| 婷婷丁香在线五月| 91九色精品人成在线观看| 韩国精品一区二区三区| 99精品欧美一区二区三区四区| av网站免费在线观看视频| 91麻豆av在线| 亚洲国产精品合色在线| 日本 av在线| 桃红色精品国产亚洲av| 欧美日韩亚洲国产一区二区在线观看| 国产精品野战在线观看| 人人妻,人人澡人人爽秒播| 视频区欧美日本亚洲| 大码成人一级视频| 亚洲人成网站在线播放欧美日韩| 一a级毛片在线观看| 久久人妻熟女aⅴ| 国产精品 欧美亚洲| 亚洲欧美激情在线| 91精品三级在线观看| 最新在线观看一区二区三区| 亚洲人成网站在线播放欧美日韩| 天堂√8在线中文| 国产97色在线日韩免费| 无限看片的www在线观看| 成人特级黄色片久久久久久久| 久久精品国产清高在天天线| 狠狠狠狠99中文字幕| 少妇的丰满在线观看| 久久精品国产综合久久久| 国产三级在线视频| or卡值多少钱| 两个人免费观看高清视频| а√天堂www在线а√下载| 成人国语在线视频| netflix在线观看网站| www.自偷自拍.com| 国产欧美日韩一区二区三| 久久人妻福利社区极品人妻图片| 免费久久久久久久精品成人欧美视频| av电影中文网址| 精品久久久久久久人妻蜜臀av | 日韩欧美在线二视频| 亚洲色图综合在线观看| 在线观看www视频免费| 久久精品aⅴ一区二区三区四区| 亚洲熟女毛片儿| 一区二区日韩欧美中文字幕| 亚洲av第一区精品v没综合| 久久久国产成人免费| 国产精品电影一区二区三区| 黑人巨大精品欧美一区二区mp4| 日韩欧美三级三区| www.www免费av| 国产片内射在线| 欧美性长视频在线观看| 美女午夜性视频免费| 视频在线观看一区二区三区| 日韩大尺度精品在线看网址 | 亚洲成a人片在线一区二区| 在线播放国产精品三级| 日韩中文字幕欧美一区二区| 一进一出好大好爽视频| 国产欧美日韩一区二区精品| 两人在一起打扑克的视频| 亚洲色图综合在线观看| 亚洲午夜精品一区,二区,三区| 亚洲少妇的诱惑av| 国产亚洲精品第一综合不卡| 国产一区二区激情短视频| 亚洲性夜色夜夜综合| 满18在线观看网站| 亚洲av电影不卡..在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 精品国产美女av久久久久小说| 国产亚洲精品久久久久5区| 校园春色视频在线观看| 亚洲五月色婷婷综合| 国产欧美日韩一区二区三| 精品免费久久久久久久清纯| 免费在线观看亚洲国产| 黑人操中国人逼视频| 91字幕亚洲| 国产精品av久久久久免费| 欧美日韩一级在线毛片| 黄片播放在线免费| 成人特级黄色片久久久久久久| 欧美午夜高清在线| 日韩av在线大香蕉| 欧美日韩亚洲国产一区二区在线观看| 香蕉久久夜色| 国产亚洲精品第一综合不卡| 免费看美女性在线毛片视频| 老司机午夜十八禁免费视频| 真人做人爱边吃奶动态| 亚洲中文字幕一区二区三区有码在线看 | 色哟哟哟哟哟哟| 十八禁网站免费在线| 一本久久中文字幕| 色综合婷婷激情| 国产精品久久久久久亚洲av鲁大| 少妇裸体淫交视频免费看高清 | 嫩草影视91久久| 欧美激情高清一区二区三区| 色综合欧美亚洲国产小说| 无限看片的www在线观看| 在线观看免费视频网站a站| 久久久久国内视频| 日日爽夜夜爽网站| 在线视频色国产色| 最近最新中文字幕大全免费视频| 国产在线精品亚洲第一网站| 亚洲国产毛片av蜜桃av| 亚洲黑人精品在线| 国产成人啪精品午夜网站| 变态另类成人亚洲欧美熟女 | 自拍欧美九色日韩亚洲蝌蚪91| 曰老女人黄片| 亚洲,欧美精品.| 99在线视频只有这里精品首页| 国产精品乱码一区二三区的特点 | 非洲黑人性xxxx精品又粗又长| 色综合婷婷激情| 无遮挡黄片免费观看| 满18在线观看网站| www.熟女人妻精品国产| 日本黄色视频三级网站网址| 国产一区二区在线av高清观看| 操美女的视频在线观看| 女同久久另类99精品国产91| 国产午夜精品久久久久久| 大码成人一级视频| 免费在线观看影片大全网站| 国产成人精品无人区| 男男h啪啪无遮挡| av欧美777| 国产麻豆69| 757午夜福利合集在线观看| 波多野结衣高清无吗| 久久精品91无色码中文字幕| 超碰成人久久| 国产三级黄色录像| 国产精品一区二区精品视频观看| 999精品在线视频| 亚洲七黄色美女视频| 真人做人爱边吃奶动态| 日本 av在线| 女人被狂操c到高潮| 欧美一区二区精品小视频在线| 亚洲av日韩精品久久久久久密| 国产高清视频在线播放一区| 51午夜福利影视在线观看| 国产精品 欧美亚洲| 中文字幕人成人乱码亚洲影| 黄色片一级片一级黄色片| 一级黄色大片毛片| 极品教师在线免费播放| 久久人人爽av亚洲精品天堂| 国产又色又爽无遮挡免费看| 极品人妻少妇av视频| 国产精品99久久99久久久不卡| 久久久精品欧美日韩精品| 黄色a级毛片大全视频| 国产精品综合久久久久久久免费 | 久久中文字幕一级| 成人18禁高潮啪啪吃奶动态图| 亚洲精品一区av在线观看| 两性夫妻黄色片| 国产乱人伦免费视频| 制服丝袜大香蕉在线| 久久 成人 亚洲| 久久中文看片网| 女人被狂操c到高潮| 999精品在线视频| 日韩av在线大香蕉| 色播在线永久视频| 精品一品国产午夜福利视频| 美女 人体艺术 gogo| 国产一级毛片七仙女欲春2 | 亚洲熟妇中文字幕五十中出| 色在线成人网| 十八禁网站免费在线| 国产aⅴ精品一区二区三区波| 亚洲国产看品久久| 久久草成人影院| 无遮挡黄片免费观看| 亚洲欧美激情综合另类| 老熟妇仑乱视频hdxx| 国产亚洲av嫩草精品影院| 97人妻天天添夜夜摸| 制服丝袜大香蕉在线| 免费一级毛片在线播放高清视频 | 99久久综合精品五月天人人| 999久久久国产精品视频| 亚洲五月婷婷丁香| 欧美激情久久久久久爽电影 | 中文字幕精品免费在线观看视频| 亚洲第一青青草原| 午夜久久久在线观看| 一a级毛片在线观看| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品第一综合不卡| 精品一区二区三区四区五区乱码| 一进一出抽搐gif免费好疼| 午夜精品在线福利| 久久香蕉精品热| 视频在线观看一区二区三区| 国产精品1区2区在线观看.| 久久人妻熟女aⅴ| 中文字幕高清在线视频| 在线免费观看的www视频| 正在播放国产对白刺激| 欧美黑人精品巨大| 成年女人毛片免费观看观看9| 成人国产一区最新在线观看| 国产av又大| 欧美日韩黄片免| 日韩大码丰满熟妇| 精品少妇一区二区三区视频日本电影| 最近最新中文字幕大全免费视频| 在线国产一区二区在线| 亚洲精品美女久久av网站| 国产精品一区二区在线不卡| 亚洲精品久久成人aⅴ小说| 97超级碰碰碰精品色视频在线观看| 亚洲精品国产区一区二| 麻豆av在线久日| 不卡av一区二区三区| 俄罗斯特黄特色一大片| 亚洲av成人一区二区三| 女性生殖器流出的白浆| 欧美最黄视频在线播放免费| 亚洲男人天堂网一区| 久久影院123| 美女高潮喷水抽搐中文字幕| 国产成人av激情在线播放| 久久精品影院6| 国产成人av激情在线播放| 老汉色av国产亚洲站长工具| 可以在线观看毛片的网站| 日本精品一区二区三区蜜桃| 午夜久久久久精精品| 亚洲片人在线观看| 91麻豆av在线| 午夜福利18| 十分钟在线观看高清视频www| 欧美日本亚洲视频在线播放| 日本撒尿小便嘘嘘汇集6| 亚洲专区中文字幕在线| 操出白浆在线播放| 麻豆久久精品国产亚洲av| 91国产中文字幕| 国产午夜福利久久久久久| 成年人黄色毛片网站| 亚洲精品国产区一区二| 欧美色视频一区免费| 成人av一区二区三区在线看| 波多野结衣一区麻豆| 午夜精品国产一区二区电影| 波多野结衣一区麻豆| 国产精品久久久久久亚洲av鲁大| 日本三级黄在线观看| 国产一区二区激情短视频| 国产精品国产高清国产av| 久久香蕉激情| 在线观看一区二区三区| 日日摸夜夜添夜夜添小说| 日韩av在线大香蕉| 午夜精品国产一区二区电影| www.精华液| 国产亚洲欧美在线一区二区| av免费在线观看网站| 无人区码免费观看不卡| 少妇被粗大的猛进出69影院| 99久久99久久久精品蜜桃| 免费看十八禁软件| 午夜免费成人在线视频| 国产日韩一区二区三区精品不卡| 欧美激情高清一区二区三区| 欧美日韩黄片免| 91麻豆av在线| 国产午夜精品久久久久久| 亚洲av五月六月丁香网| 国产极品粉嫩免费观看在线| 国产av一区二区精品久久| 在线国产一区二区在线| 亚洲精品久久国产高清桃花| 亚洲少妇的诱惑av| 9色porny在线观看| 一夜夜www| 悠悠久久av| 欧美另类亚洲清纯唯美| videosex国产| 亚洲欧美激情在线| 久久精品国产99精品国产亚洲性色 | 成在线人永久免费视频| 亚洲三区欧美一区| 成年女人毛片免费观看观看9| 十八禁网站免费在线| 日韩精品中文字幕看吧| 18禁黄网站禁片午夜丰满| 国产精品久久久人人做人人爽| 亚洲中文日韩欧美视频| 日韩高清综合在线| 午夜福利18| 一级毛片女人18水好多| 人成视频在线观看免费观看| 日韩欧美国产一区二区入口| 亚洲天堂国产精品一区在线| 午夜亚洲福利在线播放| 亚洲国产欧美日韩在线播放| 精品久久久久久成人av| 一个人免费在线观看的高清视频| 亚洲熟女毛片儿| 欧美一级毛片孕妇| 深夜精品福利| 色综合婷婷激情| 亚洲一码二码三码区别大吗| 看免费av毛片| 亚洲精品在线观看二区| 国产蜜桃级精品一区二区三区| 亚洲欧美激情综合另类| 亚洲国产看品久久| 亚洲在线自拍视频| 久久中文字幕一级| 在线观看66精品国产| 嫩草影视91久久| aaaaa片日本免费| 在线视频色国产色| 在线播放国产精品三级| 一二三四社区在线视频社区8| 在线国产一区二区在线| 亚洲精品中文字幕一二三四区| 国产精品一区二区精品视频观看| 两性夫妻黄色片| 人人妻,人人澡人人爽秒播| 成熟少妇高潮喷水视频| 亚洲,欧美精品.| 中文字幕高清在线视频| 18禁美女被吸乳视频| 成人特级黄色片久久久久久久| 麻豆国产av国片精品| 精品国产超薄肉色丝袜足j| 黄色片一级片一级黄色片| 国产在线观看jvid| av天堂久久9| 久久久久久大精品| 亚洲午夜精品一区,二区,三区| 国产亚洲av嫩草精品影院| 国产精品久久久久久亚洲av鲁大| 亚洲av电影不卡..在线观看| 深夜精品福利| 女人爽到高潮嗷嗷叫在线视频| a级毛片在线看网站| 午夜视频精品福利| 精品久久蜜臀av无|