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

    不同水分脅迫處理對(duì)陸地棉核心種質(zhì)資源的影響

    2024-06-16 15:51:55李江博高文舉運(yùn)曉東趙杰銀耿世偉韓春斌陳全家陳琴
    關(guān)鍵詞:隸屬函數(shù)水分脅迫棉花

    李江博 高文舉 運(yùn)曉東 趙杰銀 耿世偉 韓春斌 陳全家 陳琴

    摘要:為研究不同水分脅迫下陸地棉對(duì)干旱的響應(yīng)能力并篩選棉花抗旱關(guān)鍵指標(biāo)和優(yōu)異抗旱種質(zhì)資源,選取30份陸地棉核心種質(zhì)為試驗(yàn)材料,以全生育期正常灌水處理為對(duì)照,設(shè)置播種后和花鈴期各澆水1次和花鈴期斷水2次2種脅迫處理,在蕾期、花鈴期和吐絮期測(cè)定相應(yīng)指標(biāo),通過(guò)描述性統(tǒng)計(jì)分析、差異分析、主成分分析和相關(guān)性分析,采用抗旱綜合度量值(D)進(jìn)行各材料的抗旱性評(píng)價(jià)。結(jié)果表明,干旱脅迫對(duì)不同棉花材料的生長(zhǎng)發(fā)育均有不同程度的影響,播種后和花鈴期各澆一水處理的D 值離散較大(0.285~0.774),能更清晰地區(qū)分不同品種的抗旱性。利用D 值可將30份材料分為4類(lèi):第Ⅰ類(lèi)為抗旱材料,包括‘中棉所41‘新陸早7號(hào)等6個(gè)品種;第Ⅱ類(lèi)為中抗材料,包括‘晉棉46‘新陸早31號(hào)等11個(gè)品種;第Ⅲ類(lèi)為敏感材料,包括‘中棉所17‘魯1138等5個(gè)品種;第Ⅳ類(lèi)為高敏材料,包括‘新陸中8號(hào)‘酒棉8號(hào)等8個(gè)品種。果枝始節(jié)數(shù)、果枝數(shù)、株高、皮面產(chǎn)量和衣分5個(gè)指標(biāo)對(duì)干旱較為敏感,可作為棉花抗旱評(píng)價(jià)的關(guān)鍵指標(biāo)。以上研究結(jié)果可為棉花抗旱育種提供參考。

    關(guān)鍵詞:棉花;水分脅迫;隸屬函數(shù);抗旱指標(biāo);抗旱性評(píng)價(jià)

    doi:10.13304/j.nykjdb.2022.0732

    中圖分類(lèi)號(hào):S562 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):1008‐0864(2024)03‐0026‐14

    中國(guó)是全球最大的棉花生產(chǎn)國(guó)和消費(fèi)國(guó),棉花產(chǎn)業(yè)對(duì)中國(guó)國(guó)民經(jīng)濟(jì)的發(fā)展具有重要意義。新疆是我國(guó)最大的棉花種植區(qū),棉花是新疆主要的經(jīng)濟(jì)作物,棉花種植面積占新疆農(nóng)業(yè)總面積的1/3以上[1‐2]。隨著全球氣候變暖和環(huán)境日益惡化,持續(xù)干旱、土地荒漠化、大氣污染、水污染等問(wèn)題極大影響了我國(guó)農(nóng)業(yè)生產(chǎn)。新疆農(nóng)業(yè)生產(chǎn)的主要限制因子是干旱缺水,干旱對(duì)農(nóng)業(yè)造成的損失相當(dāng)于其他非生物自然災(zāi)害造成損失的總和,雖然棉花是相對(duì)比較抗旱的作物,但干旱仍會(huì)對(duì)棉花的品質(zhì)和產(chǎn)量產(chǎn)生一定影響[3]。為了降低干旱缺水對(duì)棉花生產(chǎn)造成的損失,必需加快培育高產(chǎn)抗旱品種,這是降低干旱影響的一種有效方法[4]。增強(qiáng)棉花抗旱性有助于提高棉花植株的水分利用效率,從而提高棉花產(chǎn)量,減少因干旱造成的直接經(jīng)濟(jì)損失。因此,篩選和鑒定抗旱性強(qiáng)的親本資源材料對(duì)棉花抗旱品種培育有重要意義。

    棉花種質(zhì)資源收集與保存數(shù)量日益增加,使得高效利用這些資源變得十分困難,因此對(duì)種質(zhì)資源的評(píng)價(jià)顯得越來(lái)越重要。Frankel等[5]于1984年首次提出了核心種質(zhì)的概念,即以最少數(shù)量的遺傳資源最大限度地代表整個(gè)資源群體的遺傳多樣性及整個(gè)群體的地理分布,以此來(lái)提高種質(zhì)資源的利用效率。核心種質(zhì)遴選了不同地理位置的材料、表型性狀變異廣泛的材料、品種之間具有差異(野生種/馴化種/地方品種等)的材料。材料的選擇直接決定了核心種質(zhì)的表型和遺傳變異組成,因此,核心種質(zhì)的選擇應(yīng)盡可能地保證種質(zhì)資源的遺傳多樣性。

    單一的指標(biāo)和單一的評(píng)價(jià)方法只能反映脅迫條件下某一性狀對(duì)干旱的敏感性,而不能有效地反映干旱脅迫下作物的綜合表現(xiàn)??购敌跃C合度量值(drought resistance comprehensive evaluation value,D)是使用多個(gè)性狀的加權(quán)隸屬函數(shù)值,其整合了不同性狀的抗旱系數(shù)[6‐7],能有效反映干旱脅迫下農(nóng)作物的綜合性能。這種綜合評(píng)價(jià)方法已在多種作物中應(yīng)用[8-15],已有30個(gè)性狀被作為耐旱性評(píng)價(jià)的重要指標(biāo)。這些指標(biāo)被廣泛用于小麥、棉花和玉米的抗旱性研究[16-18]。本研究設(shè)置2種水分脅迫處理,基于12個(gè)性狀對(duì)棉花的抗旱性進(jìn)行綜合評(píng)價(jià),同時(shí)對(duì)2種脅迫方式進(jìn)行評(píng)價(jià),以探究2種脅迫之間的差異,通過(guò)綜合分析法篩選抗旱種質(zhì)和關(guān)鍵性狀。

    1 材料與方法

    1.1 試驗(yàn)材料

    在課題組前期研究[19-21]基礎(chǔ)上,根據(jù)Frankel等[5]對(duì)核心種質(zhì)的定義,以不同地理來(lái)源、表型變異廣泛為原則,本研究選取了來(lái)自國(guó)內(nèi)3大棉區(qū)(西北內(nèi)陸4 份、長(zhǎng)江流域3 份、黃河流域13 份)的20份種質(zhì)資源和來(lái)自國(guó)外5個(gè)地區(qū)(美國(guó)6份、澳大利亞1份、非洲1份、法國(guó)1份、保加利亞1份)的10份種質(zhì)資源。不同材料的12個(gè)表型性狀在正常處理下變異系數(shù)在8.47%~30.64%。本研究的30份棉花種質(zhì)資源(表1)均由新疆農(nóng)業(yè)大學(xué)農(nóng)學(xué)院作物遺傳改良與種質(zhì)創(chuàng)新重點(diǎn)實(shí)驗(yàn)室收集并提供。

    1.2 試驗(yàn)方法

    試驗(yàn)于2020 年4—10 月在新疆沙灣市(44°31′N(xiāo)、85°41′E)進(jìn)行。沙灣市屬于大陸性中溫帶干旱氣候區(qū),年有效積溫較大。無(wú)霜期164~173 d,最長(zhǎng)190 d。年均降雨量174.3 mm,年均降雪量74.3 mm。沙灣市獨(dú)特的光熱資源和晝夜溫差為棉花的生長(zhǎng)發(fā)育提供了有利條件。

    田間試驗(yàn)采用隨機(jī)區(qū)組設(shè)計(jì),小區(qū)為1膜3行種植,膜幅1.80 m,種植行距76 cm,行長(zhǎng)3.14 m。設(shè)置正常灌水處理(CK)、干旱脅迫處理Ⅰ(W1)和干旱脅迫處理Ⅱ(W2),每個(gè)處理2次重復(fù)。正常灌水處理:在棉花播種后1 d(4月中旬)灌水1次,灌水時(shí)間和灌水量同干旱脅迫處理,其后灌地9次。干旱脅迫處理Ⅰ:在棉花播種后1 d(4月中旬)灌水1次,花鈴期再灌水1次。干旱脅迫處理Ⅱ:在棉花的花鈴期(當(dāng)?shù)貢r(shí)間7月5日)停水2次(停水周期10 d),其余灌溉時(shí)間(7月25日后)和灌水量同正常處理。

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

    于蕾期測(cè)定果枝始節(jié)數(shù)(the frist node fruitbranch,F(xiàn)NFB)和始節(jié)高度(height of the frist nodefruit branch,HFNFB)。始節(jié)高度(HFNFB)為棉花現(xiàn)蕾后從子葉節(jié)至第1果枝節(jié)位的高度。

    花期測(cè)定果枝數(shù)(fruit branch,F(xiàn)B)、株高(plant height,PH)、上五果枝葉片數(shù)(leaves on topfive fruit branches number,LTFFBN)、主莖葉片數(shù)(main stem leaves number,MSLN)、開(kāi)花數(shù)(flowernumber,F(xiàn)N)、棉鈴數(shù)(boll number,BN)。果枝數(shù)(FB)為棉株主莖果枝數(shù)量,株高(PH)為棉株子葉節(jié)到主莖頂端的主莖高度,上五果枝葉片數(shù)(LTFFBN)為棉株上五臺(tái)果枝葉片數(shù)量,主莖葉片數(shù)(MSLN)為棉株主莖葉片數(shù)量,盛花期記錄每日開(kāi)花數(shù)(FN),吐絮期記錄棉株結(jié)鈴數(shù)(BN)。

    吐絮期測(cè)定籽棉產(chǎn)量(cotton seed yield,CSY)、皮棉產(chǎn)量(cotton lint yield,CLY)、衣分(lint percent,LP)、單株產(chǎn)量(yield per plant,YPP)。籽棉產(chǎn)量(CSY)為收獲20鈴直接稱重為籽棉重,皮棉產(chǎn)量(CLY)為軋花后稱取20鈴棉絮質(zhì)量為皮棉產(chǎn)量,完整收獲5株植株的所有棉鈴,其籽棉產(chǎn)量的均值記作該品種的單株產(chǎn)量(YPP)。

    從每個(gè)處理的2個(gè)重復(fù)中連續(xù)選擇3株長(zhǎng)勢(shì)均勻的植株進(jìn)行測(cè)量,按照《棉花種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)》[22]中的方法對(duì)30份種質(zhì)材料的相關(guān)指標(biāo)進(jìn)行調(diào)查記錄。

    1.4 數(shù)據(jù)處理

    用Excel 2013進(jìn)行基礎(chǔ)數(shù)據(jù)統(tǒng)計(jì)和計(jì)算;用SPSS 25.0 進(jìn)行主成分分析(principal componentanalysis,PCA);用R-4.0.5 進(jìn)行差異分析、聚類(lèi)分析、相關(guān)性分析??购翟u(píng)價(jià)包括干旱變異指數(shù)(drought variability index,VId)、抗旱系數(shù)(droughtresistance coefficient,DC)、抗旱性綜合度量值(D)。

    2 結(jié)果與分析

    2.1 單一指標(biāo)在不同處理下的統(tǒng)計(jì)分析

    30份棉花種質(zhì)資源在不同干旱脅迫處理下的描述性統(tǒng)計(jì)結(jié)果如表2所示,不同指標(biāo)對(duì)干旱的響應(yīng)程度不同,其中有9個(gè)性狀的均值在2種干旱脅迫處理下均低于正常灌水處理。不同處理下各組性狀的變異系數(shù)表現(xiàn)為:正常灌水處理(CK)在8.47%~30.64%,干旱脅迫處理Ⅰ(W1)在9.24%~30.06%,干旱脅迫處理Ⅱ(W2)在6.42%~35.43%,3種處理的單株產(chǎn)量變異系數(shù)均最大。

    變異系數(shù)顯示了不同材料間各性狀存在的差異,其變化幅度反映了某一指標(biāo)對(duì)干旱的響應(yīng)程度,變化幅度越大說(shuō)明其對(duì)干旱脅迫越敏感。為了更準(zhǔn)確地描述這種變化幅度,引入干旱變異指數(shù)(VId)。不同程度的干旱脅迫使棉花性狀產(chǎn)生了不同幅度的變化。由表3可知,干旱脅迫處理Ⅰ(W1)中只有果枝始節(jié)數(shù)、上五果枝葉片數(shù)、單株產(chǎn)量3個(gè)性狀的干旱變異指數(shù)小于10%,其余9個(gè)指標(biāo)的干旱變異指數(shù)在10.12%~57.55%,干旱變異指數(shù)最大的是果枝始節(jié)高(57.55%),最小的是單株產(chǎn)量(1.91%)。

    干旱脅迫處理Ⅱ(W2)只有果枝始節(jié)數(shù)和皮棉產(chǎn)量2個(gè)性狀的干旱變異指數(shù)小于10%,其余10個(gè)指標(biāo)的干旱變異指數(shù)在10.79%~56.13%,干旱變異指數(shù)最大的是主莖葉片數(shù)(56.13%),最小的是皮棉產(chǎn)量(7.40%)。在2種干旱脅迫下共有8項(xiàng)指標(biāo)的變異指數(shù)大于10%,說(shuō)明選取的指標(biāo)具有一定的代表性。果枝始節(jié)高是2種脅迫下綜合干旱變異指數(shù)最大的,說(shuō)明此指標(biāo)受環(huán)境影響較大,對(duì)干旱較敏感。

    2.2 棉花不同性狀對(duì)干旱的響應(yīng)分析

    為了進(jìn)一步探究棉花種質(zhì)資源在不同干旱脅迫處理下各性狀的表現(xiàn),對(duì)12種指標(biāo)進(jìn)行差異分析。由圖1可知,干旱脅迫處理Ⅰ(W1)下,12個(gè)性狀與對(duì)照(CK)相比均有顯著差異;干旱脅迫處理Ⅱ(W2)下,12個(gè)性狀中有果枝數(shù)、株高、開(kāi)花數(shù)、棉鈴數(shù)、上五果枝葉片數(shù)、主莖葉片數(shù)、籽棉產(chǎn)量、皮棉產(chǎn)量、單株產(chǎn)量9個(gè)性狀與對(duì)照(CK)有顯著差異。2種脅迫處理相比,果枝始節(jié)高、株高、籽棉產(chǎn)量、皮棉產(chǎn)量、單株產(chǎn)量5個(gè)性狀差異顯著。

    2種脅迫處理與對(duì)照相比,有差異性的性狀不同,說(shuō)明不同性狀指標(biāo)對(duì)不同程度干旱有不同的響應(yīng),相較于其他性狀,W2處理下的始節(jié)數(shù)、始節(jié)高、衣分3個(gè)性狀與對(duì)照均沒(méi)有顯著差異。2種脅迫處理相比,營(yíng)養(yǎng)生長(zhǎng)階段的一些指標(biāo)如果枝始節(jié)數(shù)、果枝數(shù)、上五果枝葉片數(shù)、主莖葉片數(shù)沒(méi)有明顯的差異,但生殖生長(zhǎng)階段的相關(guān)指標(biāo)如籽棉產(chǎn)量、皮棉產(chǎn)量、單株產(chǎn)量卻差異顯著,可能是棉花在受到干旱脅迫影響時(shí)營(yíng)養(yǎng)生長(zhǎng)和生殖生長(zhǎng)同時(shí)進(jìn)行的平衡會(huì)被打破,而優(yōu)先進(jìn)行營(yíng)養(yǎng)生長(zhǎng)來(lái)滿足植株生長(zhǎng)的基本所需。

    2.3 主成分分析篩選關(guān)鍵抗旱指標(biāo)分析

    主成分分析(PCA)是將多個(gè)指標(biāo)降維后轉(zhuǎn)換為幾個(gè)關(guān)鍵指標(biāo),降低數(shù)據(jù)維度,而不會(huì)丟失原始信息,本研究利用12個(gè)指標(biāo)進(jìn)行主成分分析。由表4可知,干旱脅迫處理Ⅰ(W1)下通過(guò)對(duì)抗旱系數(shù)進(jìn)行主成分分析,其中特征值>1的5個(gè)主成分解釋的總方差達(dá)到77.307%。主成分1(principalcomponent 1,PC1)的特征值為3.015,貢獻(xiàn)率為25.128%,其中主莖葉片數(shù)、果枝始節(jié)高,主莖葉片數(shù)的荷載最大,代表植株生長(zhǎng)因子; 主成分2(principal component 2,PC2)的特征值為2.564,貢獻(xiàn)率為21.365%,其中皮棉產(chǎn)量、衣分,皮棉產(chǎn)量荷載最大,代表棉花質(zhì)量因子;主成分3(principalcomponent 3,PC3)的特征值為1.381,貢獻(xiàn)率為11.506%,其中單株產(chǎn)量荷載最大,代表產(chǎn)量因子;主成分4(principal component 4,PC4)特征值為1.303,貢獻(xiàn)率為 10.862%,上五果枝葉片數(shù)荷載最大,代表光和作用因子;在主成分5(principalcomponent 5,PC5)中,特征值為 1.014,貢獻(xiàn)率為8.447%,其中株高荷載最大,代表株高因子。主成分分析結(jié)果表明,在干旱脅迫處理Ⅰ(W1)后,主莖葉片數(shù)、皮棉產(chǎn)量、單株產(chǎn)量、上五果枝葉片數(shù)和株高變化較明顯。

    由表5可知,干旱脅迫處理Ⅱ(W2)下通過(guò)對(duì)抗旱系數(shù)進(jìn)行主成分分析,其中特征值>1的4個(gè)主成分解釋總方差的65.778%。PC1的特征值為2.938,貢獻(xiàn)率為24.480%,其中果枝始節(jié)高、主莖葉片數(shù),果枝始節(jié)高的荷載最大,代表植株發(fā)育因子;PC2的特征值為2.048,貢獻(xiàn)率為17.064%,其中衣分荷載最大,代表棉花質(zhì)量因子;PC3的特征值為1.670,貢獻(xiàn)率為13.916%,其中皮棉產(chǎn)量荷載最大,代表產(chǎn)量構(gòu)成因子;PC4的特征值為1.238,貢獻(xiàn)率為 10.319%,其中上五果枝葉片數(shù)荷載最大,代表光和作用因子。主成分分析結(jié)果表明,在干旱脅迫處理Ⅱ(W2)后,果枝始節(jié)高、衣分、皮棉產(chǎn)量、上五果枝葉片數(shù)變化較明顯。

    2.4 棉花的抗旱性評(píng)價(jià)分析

    棉花的抗旱性受多種因素的影響,靠單一指標(biāo)不能準(zhǔn)確評(píng)價(jià)各材料間抗旱性的差異,有必要對(duì)各種質(zhì)資源進(jìn)行綜合評(píng)價(jià)。根據(jù)指標(biāo)權(quán)重和隸屬函數(shù)計(jì)算抗旱性綜合度量值(D),D 值越大表明其抗旱能力越強(qiáng),反之其對(duì)干旱越敏感。由表6可知,W1處理下‘中棉所41‘晉棉46‘新陸早7號(hào)等材料具有較好的抗旱性;‘魯棉2號(hào)‘新陸中8號(hào)‘中棉所50較為敏旱。在W2處理下,‘新陸早7號(hào)‘徐州6號(hào)‘泗168具有較好的抗旱性;‘保2367‘新陸中8號(hào)‘酒棉8號(hào)較為敏旱。

    通過(guò)試驗(yàn)發(fā)現(xiàn),干旱脅迫處理Ⅰ(W1)的D 值離散較大(0.285~0.774),產(chǎn)量損失較大(表2),只有正常產(chǎn)量的58%。干旱脅迫處理Ⅱ(W2)的D值離散較小(0.231~0.697),產(chǎn)量損失較少(表2),占正常產(chǎn)量的77%。干旱脅迫處理Ⅰ(W1)的D值離散程度較大,說(shuō)明該處理能更清晰地凸顯棉花的抗旱性,能較好地區(qū)分抗旱性材料和敏旱性材料,但以產(chǎn)量損失嚴(yán)重為代價(jià)。干旱脅迫處理Ⅱ(W2)雖然產(chǎn)量損失較少,但判斷抗旱性相近的資源材料時(shí)沒(méi)有更明顯的效果。

    2.5 不同處理下各性狀與D 值的相關(guān)性分析

    為進(jìn)一步探究各性狀與抗旱性綜合度量值(D)之間的關(guān)系,利用12個(gè)性狀在2種脅迫處理下的抗旱系數(shù)與D 值進(jìn)行相關(guān)性分析。由圖2可知,在干旱脅迫處理Ⅰ(W1)下,果枝始節(jié)數(shù)(FNFB)、果枝始節(jié)高(HFNFB)、果枝數(shù)(FN)、株高(PH)、主莖葉片數(shù)(MSLN)、皮棉產(chǎn)量(CLY)和衣分(LP)7個(gè)指標(biāo)與D 值存在極顯著正相關(guān);開(kāi)花數(shù)(FN)、棉鈴數(shù)(BN)、上五果枝葉片數(shù)(LTFFBN)、籽棉產(chǎn)量(CSY)、單株產(chǎn)量(YPP)與D 值相關(guān)不顯著。

    由圖3可知,在干旱脅迫處理Ⅱ(W2)下,果枝始節(jié)數(shù)(FNFB)、果枝數(shù)(FN)、株高(PH)、皮棉產(chǎn)量(CLY)和衣分(LP)5個(gè)指標(biāo)與D值存在極顯著正相關(guān);單株產(chǎn)量(YPP)與D 值存在極顯著負(fù)相關(guān);果枝始節(jié)高(HFNFB)、開(kāi)花數(shù)(FN)、棉鈴數(shù)(BN)、上五果枝葉片數(shù)(LTFFBN)、主莖葉片數(shù)(MSLN)、籽棉產(chǎn)量(CSY)與D值相關(guān)性不顯著。

    相關(guān)性結(jié)果分析發(fā)現(xiàn),2 種脅迫處理下與D值存在極顯著正相關(guān)的共有指標(biāo)有5個(gè),分別為果枝始節(jié)數(shù)、果枝數(shù)、株高、皮棉產(chǎn)量和衣分,說(shuō)明這5個(gè)指標(biāo)對(duì)干旱較為敏感,能較早對(duì)干旱脅迫做出響應(yīng),因此,這5個(gè)指標(biāo)可以作為抗旱性評(píng)價(jià)的優(yōu)選指標(biāo)。

    2.6 棉花耐旱指標(biāo)綜合篩選分析

    不同性狀對(duì)干旱的響應(yīng)顯示,2種脅迫與對(duì)照相比,果枝始節(jié)數(shù)、株高、開(kāi)花數(shù)、棉鈴數(shù)、上五果枝葉片數(shù)、主莖葉片數(shù)、籽棉產(chǎn)量、皮棉產(chǎn)量、單株產(chǎn)量9個(gè)性狀有顯著性差異。PCA降維法通過(guò)簡(jiǎn)化多個(gè)變量選取關(guān)鍵主成分作為代表性變量發(fā)現(xiàn),在2種處理下皮棉產(chǎn)量和上五果枝葉片數(shù)都有較高的貢獻(xiàn)率,被選為關(guān)鍵指標(biāo)。不同處理下,各性狀與D 值的相關(guān)性結(jié)果顯示,2種脅迫下與D值都顯著相關(guān)的指標(biāo)共有5個(gè),分別為果枝始節(jié)數(shù)、果枝數(shù)、株高、皮棉產(chǎn)量和衣分。綜合上述結(jié)果發(fā)現(xiàn),在所有的分析中皮棉產(chǎn)量與棉花的抗旱性都高度相關(guān),因此認(rèn)為皮棉產(chǎn)量可以作為關(guān)鍵的棉花耐旱指標(biāo)。

    2.7 2 種脅迫處理綜合聚類(lèi)分析

    根據(jù)指標(biāo)權(quán)重和隸屬函數(shù)計(jì)算D 值,并根據(jù)2種脅迫處理的D 值綜合劃分不同材料的類(lèi)別。根據(jù)聚類(lèi)結(jié)果(圖4),30份棉花資源材料可聚為4類(lèi),第Ⅰ類(lèi)為抗旱材料,第Ⅱ類(lèi)為中抗材料,第Ⅲ類(lèi)為敏感材料,第Ⅳ類(lèi)為高敏材料。其中,第Ⅰ類(lèi)主要包括‘中棉所41‘新陸早7號(hào)等6份材料;第Ⅱ類(lèi)主要包括‘晉棉46‘新陸早31號(hào)等11份材料;第Ⅲ類(lèi)主要包括‘中棉所17‘魯1138等5份材料;第Ⅳ類(lèi)主要包括‘新陸中8號(hào)‘酒棉8號(hào)等8份材料。6份抗旱材料中5份為我國(guó)3大棉區(qū)推廣種植的材料,8份高敏材料中4份為國(guó)外引進(jìn)材料,這可能是由于地理環(huán)境的差異所引起的。

    3 討論

    俞希根等[23]、李少昆等[24]研究發(fā)現(xiàn),花鈴期干旱對(duì)棉花產(chǎn)量和品質(zhì)的影響最嚴(yán)重。因此,在花鈴期研究脅迫條件下植株的變化情況十分重要,本研究的2種脅迫處理選擇在花鈴期進(jìn)行。已有研究對(duì)作物的抗旱性評(píng)價(jià)多數(shù)只設(shè)置1種脅迫處理,本研究運(yùn)用了2種脅迫處理,鑒定出的抗旱種質(zhì)和關(guān)鍵指標(biāo)更有說(shuō)服力,同時(shí)本研究也鑒別了2種脅迫處理間的差異,為棉花抗旱育種提供理論依據(jù)。

    作物抗旱性主要體現(xiàn)在產(chǎn)量方面,因此出現(xiàn)了以抗旱系數(shù)、抗旱指數(shù)等產(chǎn)量指標(biāo)為依據(jù)的單一指標(biāo)直接評(píng)價(jià)法[25‐26]。隨著對(duì)作物抗旱問(wèn)題的深入研究,不同的抗旱評(píng)價(jià)方法也被提出。黎裕等[27]認(rèn)為,抗旱性和產(chǎn)量都具有復(fù)雜的遺傳機(jī)制,應(yīng)該從多方面進(jìn)行抗旱性評(píng)價(jià)。近年來(lái),普遍認(rèn)為結(jié)合多種指標(biāo)、多種方法的綜合抗旱性評(píng)價(jià)比較可靠,也研究和提出了基于抗旱系數(shù)、主成分分析和權(quán)重分析的加權(quán)隸屬函數(shù)法。徐銀萍等[28]利用綜合度量值對(duì)50份大麥資源進(jìn)行抗旱性鑒定發(fā)現(xiàn),基于D值鑒定的種質(zhì)與聚類(lèi)結(jié)果與作物田間實(shí)際抗旱性更為接近。劉光輝等[29]通過(guò)對(duì)90份陸地棉的抗旱鑒定發(fā)現(xiàn),D 值遺傳率最高,是適合進(jìn)行評(píng)價(jià)的方法。Sun等[30]通過(guò)3年的水分脅迫試驗(yàn)發(fā)現(xiàn),D 值與棉花的減產(chǎn)值存在線性關(guān)系,且在3年內(nèi)呈顯著相關(guān),表明D 值可以表示抗旱強(qiáng)度。

    本研究通過(guò)2種水分脅迫處理并利用多種分析方法相結(jié)合的方式對(duì)棉花的抗旱性進(jìn)行評(píng)價(jià)。主成分分析通過(guò)降維簡(jiǎn)化原有信息,可以用來(lái)解釋和描述種質(zhì)抗旱性和耐鹽性的重要指標(biāo)[31-33]。本研究通過(guò)主成分分析分別得到了5(W1處理)和4(W2處理)個(gè)主成分,這些主成分分別解釋了77.307% 和65.778% 的總變異,接著計(jì)算各個(gè)材料的隸屬函數(shù)值,最后計(jì)算D 值,再根據(jù)D 值大小對(duì)各個(gè)材料進(jìn)行抗旱性排序。最后基于D 值對(duì)2種處理下的30份棉花資源材料進(jìn)行綜合性聚類(lèi)分析,鑒定出了相對(duì)抗旱與敏旱的種質(zhì)資源。聚類(lèi)分析將30份棉花種質(zhì)資源分為4類(lèi):抗旱性材料6份、中抗材料11份、敏感材料5份、高敏材料8份,分組類(lèi)型與D 值的鑒定結(jié)果一致。綜合評(píng)價(jià)的結(jié)果與劉鵬鵬等[34]、閆成川等[19]、鄭巨云等[35]鑒定結(jié)果相比,中間材料存在較多交集,極端材料存在較少交集,這可能是由于評(píng)價(jià)指標(biāo)的選擇、試驗(yàn)環(huán)境、試驗(yàn)設(shè)計(jì)不同造成的。

    植物的耐旱機(jī)制非常復(fù)雜,對(duì)此許多研究者在指標(biāo)選擇上從各個(gè)方面做了大量研究[36]。本研究運(yùn)用的抗旱性綜合度量值(D)是基于12個(gè)指標(biāo)綜合得出的結(jié)果,通過(guò)D 值的評(píng)價(jià)可以篩選出綜合抗旱能力較強(qiáng)的種質(zhì),但為了能更準(zhǔn)確可靠地描述指標(biāo)與D 值之間的關(guān)系,應(yīng)該選擇關(guān)鍵性的指標(biāo)。本研究通過(guò)不同性狀對(duì)干旱的響應(yīng),主成分分析不同處理下各性狀與D 值的相關(guān)性分析逐步篩選出了1個(gè)關(guān)鍵抗旱指標(biāo),即皮棉產(chǎn)量。研究結(jié)果為棉花抗旱育種提供了資源。

    參 考 文 獻(xiàn)

    [1] WANG C, ISODA A, WANG P. Growth and yield performance

    of some cotton cultivars in Xinjiang China an arid area with

    short growing period [J]. J. Agron. Crop Sci., 2004, 190(3):

    177-183.

    [2] KANG Y H, WANG R S, WAN S Q, et al .. Effects of different

    water levels on cotton growth and water use through drip

    irrigation in an arid region with saline ground water of

    Northwest China [J]. Agric. Water Manage., 2012, 109(6):

    117-126.

    [3] ULLAH A, SUN H, YANG X Y, et al .. Drought coping

    strategies in cotton: increased crop per drop [J]. Plant Biotech.

    J., 2017,15(3):271-284.

    [4] CATTIVELLI L, RIZZA F, BADECK F W, et al .. Drought

    tolerance improvement in crop plants: an integrated view from

    breeding to genomics [J]. Field Crops Res., 2008,105(1):1-14.

    [5] FRANKEL O H, BROWN A H D. Current plant genetic

    resources a critical appraisal [J]. Genetics: New Front., 1984,

    4: 3-13.

    [6] ALIREZA N G, GADIR N G, DAVOUD H, et al .. Effects of

    drought stress condition on the yield and yield components of

    advanced wheat genotypes in Ardabil Iran [J]. J. Food Agric.

    Environ., 2009,7(3):228-234.

    [7] LIU C Y, YANG Z Y, HU Y G. Drought resistance of wheat

    alien chromosome addition lines evaluated by membership

    function value based on multiple traits and drought resistance

    index of grain yield [J]. Field Crops Res., 2015,179(8):103-112.

    [8] 王士強(qiáng),胡銀崗,佘奎軍,等.小麥抗旱相關(guān)農(nóng)藝性狀和生理

    生化性狀的灰色關(guān)聯(lián)度分析[J].中國(guó)農(nóng)業(yè)科學(xué),2007,40(11):

    2452-2459.

    WANG S Q, HU Y G, SHE K J, et al .. Gray relational grade

    analysis of agronomical and physi-biochemical traits related to

    drought tolerance in wheat [J]. Sci. Agric. Sin., 2007, 40(11):

    2452-2459.

    [9] 李貴全,張海燕,季蘭,等.不同大豆品種抗旱性綜合評(píng)價(jià)[J].

    應(yīng)用生態(tài)學(xué)報(bào),2006,17(12):2408-2412.

    LI G Q, ZHANG H Y, JI L, et al .. Comprehensive evaluation

    on drought-resistance of different soybean varieties [J]. Chin. J.

    Appl. Ecol., 2006,17(12):2408-2412.

    [10] 謝小玉,張霞,張兵.油菜苗期抗旱性評(píng)價(jià)及抗旱相關(guān)指標(biāo)變

    化分析[J].中國(guó)農(nóng)業(yè)科學(xué),2013,46(3):476-485.

    XIE X Y, ZHANG X, ZHANG B. Evaluation of drought

    resistance and analysis of variation of relevant parameters at

    seedling stage of rapeseed (Brassica napus L.) [J]. Sci.Agric.

    Sin., 2013,46(3):476-485.

    [11] 柯貞進(jìn),尹美強(qiáng),溫銀元,等.干旱脅迫下聚丙烯酰胺浸種對(duì)

    谷子種子萌發(fā)及幼苗期抗旱性的影響[J].核農(nóng)學(xué)報(bào),2015,

    29(3):563-570.

    KE Z J, YIN M Q, WEN Y Y, et al .. Effects of polyacrylamide

    seed soaking on seed germination and drought resistance of

    millet (Setaria italica) seedlings under drought stress [J]. J.

    Nucl. Agric. Sci., 2015,29(3):563-570.

    [12] 徐蕊,王啟柏,張春慶,等.玉米自交系抗旱性評(píng)價(jià)指標(biāo)體系

    的建立[J].中國(guó)農(nóng)業(yè)科學(xué),2009,42(1):72-84.

    XU R, WANG Q B, ZHANG C Q, et al .. Drought-resistance

    evaluation system of maize inbred [J]. Sci. Agric. Sin., 2009,

    42(1): 72-84.

    [13] 祁旭升,王興榮,張彥軍,等.胡麻成株期抗旱指標(biāo)篩選與種

    質(zhì)抗性鑒定[J].中國(guó)農(nóng)業(yè)科學(xué),2010,43(15):3076-3087.

    QI X S,WANG X R, ZAHNG Y J, et al .. Drought resistance

    indexes selection for oil flax at the adult stage and drought

    resistance identification of oil flax germplasm [J]. Sci. Agric.

    Sin., 2010,43(15):3076-3087.

    [14] 姚玉波.不同品種亞麻種子萌發(fā)期抗旱性鑒定[J].核農(nóng)學(xué)

    報(bào),2015,29(10):2033-2039.

    YAO Y B. Drought resistance identification of different flax

    varieties in seed germination stage [J]. J. Nucl. Agric. Sci.,

    2015,29(10):2033-2039.

    [15] 石運(yùn)慶,苗華榮,胡曉輝,等.花生耐鹽堿性鑒定指標(biāo)的研究

    及應(yīng)用[J].核農(nóng)學(xué)報(bào),2015,29(3):442-447.

    SHI Y Q, MIAO H R, HU X H, et al .. Research and

    application of saline-alkaline tolerance indexes of peanut [J]. J.

    Nucl. Agric. Sci., 2015,29(3): 442-447.

    [16] KRAMER D M, EVANS J R. The importance of energy

    balance in improving photosynthetic productivity [J]. Plant

    Physiol., 2011,155(1):70-78.

    [17] SOOMROO M H, MARKHAND G S. Screening Pakistani

    cotton for drought tolerance [J]. Pak. J. Bot., 2011, 44(1):

    383-388.

    [18] WIGGINS M S, LEIB B G, MUELLER T C, et al .. Investigation

    of physiological growth, fiber quality, yield, and yield stability

    of upland cotton varieties in differing environments [J]. J. Cott.

    Sci., 2013,17(3):140-148.

    [19] 閆成川,曾慶濤,陳琴,等.陸地棉花鈴期抗旱指標(biāo)篩選及評(píng)

    價(jià)[J].中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào),2022,24(7):46-57.

    YAN C C, ZENG Q T, CHEN Q, et al .. Screening and

    evaluation of drought resistance indicators at flowering and

    boll stage of upland cotton [J]. J. Agric. Sci. Technol., 2022,

    24(7):46-57.

    [20] 孫豐磊,張玻,曲延英,等.花鈴期干旱脅迫對(duì)不同棉花品種

    光合特性影響及抗旱性評(píng)價(jià)[J].干旱地區(qū)農(nóng)業(yè)研究,2018,

    36(5):7-13.

    SUN F L, ZHANG B, QU Y Y, et al .. Effects of drought stress

    during the blooming period on photosynthetic characteristics

    and assessment drought resistance of different cotton varieties [J].

    Agric. Res. Arid Areas, 2018,36(5):7-13.

    [21] 孫豐磊,曲延英,陳全家,等. 棉花抗旱相關(guān)指標(biāo)綜合評(píng)價(jià)及

    灰色關(guān)聯(lián)分析[J].干旱地區(qū)農(nóng)業(yè)研究,2019,37(1):233-239.

    SUN F L, QU Y Y, CHEN Q J, et al .. Comprehensive

    evaluation of cotton drought tolerance indexes and gray relational

    analysis [J]. Agric. Res. Arid Areas, 2019,37(1):233-239.

    [22] 杜雄明,周忠麗.棉花種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)[M].北

    京:中國(guó)農(nóng)業(yè)出版社,2005:1-89.

    [23] 俞希根,孫景生,肖俊夫,等.棉花適宜土壤水分下限和干旱

    指標(biāo)研究[J].棉花學(xué)報(bào),1999,27(1):36-39.

    YU X G, SUN J S, XIAO J F, et al .. A study on drought indices

    and lower limit of suitable soil moisture of cotton [J]. Cott. Sci.,

    1999,27(1):36-39.

    [24] 李少昆,肖璐,黃文華.不同時(shí)期干旱脅迫對(duì)棉花生長(zhǎng)和產(chǎn)量

    的影響Ⅱ棉花生長(zhǎng)發(fā)育及生理特性的變化[J].石河子大學(xué)

    學(xué)報(bào)(自然科學(xué)版),1999,3(4):259-264.

    LI S K, XIAO L, HUANG W H. Effect of drought stress on

    cotton growth and lint yield at different growing stage [J]. J.

    Shihezi. Univ. (Nat. Sci), 1999,3(4):259-264.

    [25] 蘭巨生,胡福順,張景瑞,等.作物抗旱指數(shù)的概念和統(tǒng)計(jì)方

    法[J].華北農(nóng)學(xué)報(bào),1990,5(2):20-25.

    LAN J S, HU F S, ZHANG J R, et al .. The concept and

    statistical method of drought resistance index in crops [J]. Acta

    Agric. Boreali-Sin., 1990,5(2):20-25.

    [26] ABRAHAM B, WAYNE R J. Breeding crop varieties for stress

    environments [J]. Crit. Rev. Plant Sci., 1985,2(3):199-238.

    [27] 黎裕,王天宇,劉成,等.玉米抗旱品種的篩選指標(biāo)研究[J].植

    物遺傳資源學(xué)報(bào),2004,5(3):210-215.

    LI Y, WANG T Y, LIU C, et al .. Analysis on criteria for

    screening drought tolerant maize hybrids [J]. J. Plant Genet.

    Resour., 2004,5(3):210-215.

    [28] 徐銀萍,潘永東,劉強(qiáng)德,等.大麥種質(zhì)資源成株期抗旱性鑒

    定及抗旱指標(biāo)篩選[J].作物學(xué)報(bào),2020,46(3):448-461.

    XU Y P, PAN Y D, LIU Q D, et al .. Drought resistance

    identification and drought resistance indexes screening of

    barley resources at mature period [J]. Acta Agron. Sin., 2020,

    46(3):448-461.

    [29] 劉光輝,陳全家,吳鵬昊,等.棉花花鈴期抗旱性綜合評(píng)價(jià)及

    指標(biāo)篩選[J].植物遺傳資源學(xué)報(bào),2016,17(1): 53-62.

    LIU G H, CHEN Q J, WU P H, et al .. Screening and

    comprehensive evaluation of drought resistance indices of

    cotton at blossing and boll-forming stages [J]. J. Plant Genet.

    Resour., 2016,17(1):53-62,

    [30] SUN F L, CHEN Q, CHEN Q J, et al .. Screening of key drought

    tolerance indices for cotton at the flowering and boll setting

    stage using the dimension reduction method [J/OL]. Front.

    Plant Sci., 2021, 12:619926 [2022-08-01]. https://doi.org/

    10.3389/fpls.2021.619926.

    [31] BO W, FU B C, QIN G J, et al .. Evaluation of drought

    resistance in Iris germanica L. based on subordination function

    and principal component analysis [J]. Emirates J. Food Agric.,

    2017,29(10):770-778.

    [32] NEGRAO S, SCHMOCKEL S M, TESTER M. Evaluating

    physiological responses of plants to salinity stress [J]. Ann.

    Bot., 2017,119(1): 1-11.

    [33] KAKAR N, JUMAA S H, REDONA E D, et al .. Evaluating rice

    for salinity using potculture provides a systematic tolerance

    assessment at the seedling stage [J]. Rice, 2019,57(12):1-14.

    [34] 劉鵬鵬,陳全家,曲延英,等.棉花種質(zhì)資源抗旱性評(píng)價(jià)[J].新

    疆農(nóng)業(yè)科學(xué),2014,51(11):1961-1969.

    LIU P P, CHEN Q J, QU Y Y, et al .. The drought resistance

    evaluation of cotton germplasm resources [J]. Xinjiang Agric.

    Sci. Technol., 2014,51(11):1961-1969.

    [35] 鄭巨云,桑志偉,王俊鐸,等.棉花品種抗旱性相關(guān)指標(biāo)分析

    與綜合評(píng)價(jià)[J].中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào), 2022,24(10):23-34.

    ZHENG J Y, SANG Z W, WANG J D, et al .. Indexes analysis

    and comprehensive evaluation of drought resistance of cotton

    varieties [J]. J. Agric. Sci. Technol., 2022,24(10):23-34.

    [36] 宋鳳斌,徐世昌.玉米抗旱性鑒定指標(biāo)的研究[J].中國(guó)生態(tài)

    農(nóng)業(yè)學(xué)報(bào),2004,12(1):127-129.

    SONG F B, XU S C. Study on the drought-resistant

    identification indexes in maize [J]. Chin. J. Eco-Agric., 2004,

    12(1):127-129.

    (責(zé)任編輯:胡立霞)

    猜你喜歡
    隸屬函數(shù)水分脅迫棉花
    棉花是花嗎?
    棉花
    小讀者(2020年4期)2020-06-16 03:33:54
    不同玉米品種萌芽期抗旱篩選
    混合鹽堿脅迫對(duì)苜蓿種子萌發(fā)特性的影響
    水分脅迫對(duì)蘭花“曙光”生理生化指標(biāo)的影響研究
    基于鉆進(jìn)過(guò)程的油井安全綜合評(píng)估
    水分和氮素脅迫對(duì)辣椒產(chǎn)量形成影響及抗逆性評(píng)價(jià)
    心中的“棉花糖”
    第三講 棉花肥害診斷及其防治
    久久人妻熟女aⅴ| 国产一区二区激情短视频 | 亚洲精品久久成人aⅴ小说| 在线av久久热| 在线天堂中文资源库| 精品亚洲乱码少妇综合久久| 性色av乱码一区二区三区2| 1024香蕉在线观看| 一本一本久久a久久精品综合妖精| 欧美变态另类bdsm刘玥| 少妇被粗大的猛进出69影院| 性色av乱码一区二区三区2| 色婷婷av一区二区三区视频| 人妻人人澡人人爽人人| 亚洲国产毛片av蜜桃av| 大话2 男鬼变身卡| 日韩av免费高清视频| 国产在视频线精品| 亚洲欧美色中文字幕在线| 亚洲av电影在线进入| 在线观看免费午夜福利视频| 久久精品久久精品一区二区三区| 亚洲精品国产av蜜桃| 亚洲,欧美,日韩| 国产精品国产三级专区第一集| 国产激情久久老熟女| 无限看片的www在线观看| 乱人伦中国视频| 久久鲁丝午夜福利片| 只有这里有精品99| 欧美日韩视频高清一区二区三区二| 亚洲欧美日韩另类电影网站| 欧美在线黄色| 大型av网站在线播放| 91老司机精品| 一级黄片播放器| 纯流量卡能插随身wifi吗| 国产成人系列免费观看| 可以免费在线观看a视频的电影网站| 中文字幕人妻丝袜制服| 一级毛片黄色毛片免费观看视频| 亚洲五月婷婷丁香| 国产成人免费无遮挡视频| 亚洲精品第二区| 久久亚洲精品不卡| 母亲3免费完整高清在线观看| 婷婷丁香在线五月| 电影成人av| 丝瓜视频免费看黄片| 久久久久久久大尺度免费视频| 涩涩av久久男人的天堂| 真人做人爱边吃奶动态| 日本午夜av视频| 在线天堂中文资源库| 19禁男女啪啪无遮挡网站| 国产又爽黄色视频| 热re99久久精品国产66热6| 另类精品久久| 亚洲av在线观看美女高潮| 成年av动漫网址| 国产精品久久久av美女十八| 蜜桃国产av成人99| 久久精品国产亚洲av高清一级| 91精品三级在线观看| 免费看不卡的av| 视频区欧美日本亚洲| 国产精品成人在线| 亚洲情色 制服丝袜| 一级a爱视频在线免费观看| 男人添女人高潮全过程视频| 美女大奶头黄色视频| 欧美日韩亚洲国产一区二区在线观看 | 国产熟女午夜一区二区三区| 久久青草综合色| 国产无遮挡羞羞视频在线观看| 久久av网站| 免费高清在线观看视频在线观看| 看免费av毛片| 亚洲精品一区蜜桃| 99国产精品一区二区蜜桃av | 久久精品熟女亚洲av麻豆精品| 国产免费视频播放在线视频| 国产片内射在线| xxx大片免费视频| 久久这里只有精品19| 亚洲av片天天在线观看| 少妇粗大呻吟视频| 精品免费久久久久久久清纯 | 老司机影院成人| 又黄又粗又硬又大视频| 国产精品久久久久成人av| 国产在视频线精品| 在线观看人妻少妇| 中文字幕人妻丝袜制服| 在线观看国产h片| 国产亚洲一区二区精品| 久久精品人人爽人人爽视色| av片东京热男人的天堂| 老司机在亚洲福利影院| 韩国精品一区二区三区| 成年动漫av网址| 国产又爽黄色视频| a级片在线免费高清观看视频| 国产欧美日韩精品亚洲av| 在线看a的网站| 国产97色在线日韩免费| 狂野欧美激情性bbbbbb| a级片在线免费高清观看视频| 欧美国产精品va在线观看不卡| 一本—道久久a久久精品蜜桃钙片| e午夜精品久久久久久久| 18禁裸乳无遮挡动漫免费视频| 国产成人av激情在线播放| 2021少妇久久久久久久久久久| avwww免费| 日韩一本色道免费dvd| 久久国产精品人妻蜜桃| 一边摸一边做爽爽视频免费| 欧美av亚洲av综合av国产av| 在线天堂中文资源库| 黄色一级大片看看| 久久精品国产亚洲av涩爱| 亚洲中文av在线| 男女边摸边吃奶| 国产97色在线日韩免费| 国产在线一区二区三区精| 欧美精品啪啪一区二区三区 | 另类亚洲欧美激情| 欧美日韩视频高清一区二区三区二| 精品一区在线观看国产| 高潮久久久久久久久久久不卡| 啦啦啦视频在线资源免费观看| 五月天丁香电影| 亚洲成人免费av在线播放| 中文字幕人妻丝袜制服| 久久人人爽av亚洲精品天堂| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲一码二码三码区别大吗| 18禁黄网站禁片午夜丰满| 性少妇av在线| 老司机午夜十八禁免费视频| 少妇的丰满在线观看| 国产视频首页在线观看| 精品少妇内射三级| 免费不卡黄色视频| 五月开心婷婷网| 久久鲁丝午夜福利片| 欧美激情高清一区二区三区| 男女之事视频高清在线观看 | 五月天丁香电影| 国产片特级美女逼逼视频| 蜜桃国产av成人99| 日本av手机在线免费观看| 老司机亚洲免费影院| 自线自在国产av| 大香蕉久久网| 久久中文字幕一级| 又大又黄又爽视频免费| 欧美大码av| 又黄又粗又硬又大视频| 午夜影院在线不卡| 国产av精品麻豆| 视频在线观看一区二区三区| 人人妻人人澡人人看| cao死你这个sao货| 熟女少妇亚洲综合色aaa.| 男人爽女人下面视频在线观看| 国产精品一二三区在线看| 日韩视频在线欧美| 看十八女毛片水多多多| 成人手机av| 成人国产av品久久久| 性色av一级| 美女大奶头黄色视频| 欧美xxⅹ黑人| av片东京热男人的天堂| 少妇 在线观看| 性高湖久久久久久久久免费观看| 日韩中文字幕视频在线看片| 飞空精品影院首页| 男的添女的下面高潮视频| av天堂久久9| 狂野欧美激情性bbbbbb| 亚洲男人天堂网一区| 一区二区av电影网| 精品少妇一区二区三区视频日本电影| 久久精品人人爽人人爽视色| 久久热在线av| 久久国产精品男人的天堂亚洲| 国产成人影院久久av| 丁香六月欧美| 宅男免费午夜| 精品欧美一区二区三区在线| 三上悠亚av全集在线观看| 51午夜福利影视在线观看| av天堂久久9| 人人妻,人人澡人人爽秒播 | 啦啦啦中文免费视频观看日本| 欧美av亚洲av综合av国产av| 欧美精品一区二区免费开放| 国产精品99久久99久久久不卡| 99国产精品免费福利视频| 各种免费的搞黄视频| 国产视频一区二区在线看| 成人免费观看视频高清| 国产精品一区二区免费欧美 | 久久久欧美国产精品| 18禁黄网站禁片午夜丰满| a 毛片基地| 成年人免费黄色播放视频| 精品久久蜜臀av无| 国产在线一区二区三区精| 中国美女看黄片| 亚洲,一卡二卡三卡| 另类精品久久| 青春草视频在线免费观看| 在现免费观看毛片| 不卡av一区二区三区| 人人妻,人人澡人人爽秒播 | 国产精品久久久久久人妻精品电影 | www.自偷自拍.com| 欧美黄色片欧美黄色片| 国产精品 国内视频| 18禁黄网站禁片午夜丰满| 天堂8中文在线网| 十八禁高潮呻吟视频| 视频在线观看一区二区三区| 欧美精品啪啪一区二区三区 | 亚洲欧美一区二区三区黑人| 久久 成人 亚洲| 国产成人系列免费观看| 国产熟女欧美一区二区| 欧美黄色淫秽网站| 香蕉丝袜av| 国产免费现黄频在线看| 一区二区三区精品91| 国产精品亚洲av一区麻豆| 欧美乱码精品一区二区三区| 熟女av电影| 中文字幕色久视频| 伊人久久大香线蕉亚洲五| 母亲3免费完整高清在线观看| 丁香六月欧美| 91精品三级在线观看| 欧美日韩亚洲国产一区二区在线观看 | 日本五十路高清| 一区二区三区精品91| 欧美国产精品va在线观看不卡| 国产在线免费精品| svipshipincom国产片| 99久久人妻综合| 日韩av不卡免费在线播放| 久久久久久久久免费视频了| 亚洲视频免费观看视频| 亚洲国产欧美在线一区| 在线亚洲精品国产二区图片欧美| 伊人亚洲综合成人网| 一区二区三区四区激情视频| 男女床上黄色一级片免费看| 免费av中文字幕在线| 人人妻,人人澡人人爽秒播 | 精品国产超薄肉色丝袜足j| 免费高清在线观看日韩| 国产99久久九九免费精品| 国产av一区二区精品久久| 操出白浆在线播放| 超碰97精品在线观看| 伊人亚洲综合成人网| av电影中文网址| 十八禁网站网址无遮挡| 青青草视频在线视频观看| 午夜福利乱码中文字幕| 精品久久久久久电影网| 秋霞在线观看毛片| 国产精品久久久久久人妻精品电影 | 亚洲欧洲国产日韩| 中文字幕亚洲精品专区| 如日韩欧美国产精品一区二区三区| 老司机影院成人| 老司机影院毛片| 国产精品一区二区在线不卡| 一区二区日韩欧美中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 夫妻性生交免费视频一级片| 国产精品一区二区免费欧美 | 中国美女看黄片| 2021少妇久久久久久久久久久| 亚洲av日韩精品久久久久久密 | www日本在线高清视频| 日韩熟女老妇一区二区性免费视频| 日韩一卡2卡3卡4卡2021年| 免费日韩欧美在线观看| 欧美激情高清一区二区三区| 美国免费a级毛片| 亚洲久久久国产精品| 国产激情久久老熟女| 波多野结衣av一区二区av| 午夜老司机福利片| 老熟女久久久| 久久精品久久久久久噜噜老黄| 啦啦啦啦在线视频资源| 免费高清在线观看视频在线观看| 国产视频一区二区在线看| cao死你这个sao货| 亚洲欧美中文字幕日韩二区| 一级毛片黄色毛片免费观看视频| 最近手机中文字幕大全| 女性生殖器流出的白浆| 美女福利国产在线| 女人久久www免费人成看片| 久久影院123| 国产精品香港三级国产av潘金莲 | 少妇被粗大的猛进出69影院| 欧美少妇被猛烈插入视频| 亚洲国产毛片av蜜桃av| 成年女人毛片免费观看观看9 | 又大又黄又爽视频免费| 91字幕亚洲| 成人黄色视频免费在线看| 亚洲精品第二区| 久久99一区二区三区| 咕卡用的链子| 成在线人永久免费视频| av电影中文网址| 老司机影院毛片| 国产成人a∨麻豆精品| 在线精品无人区一区二区三| 一级毛片女人18水好多 | 精品少妇内射三级| 波野结衣二区三区在线| 亚洲国产欧美网| 在线天堂中文资源库| 黄色a级毛片大全视频| 黄片小视频在线播放| √禁漫天堂资源中文www| 在线观看一区二区三区激情| 爱豆传媒免费全集在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 777米奇影视久久| 一个人免费看片子| 一本色道久久久久久精品综合| 老熟女久久久| 最新在线观看一区二区三区 | 99香蕉大伊视频| 性色av一级| 免费高清在线观看视频在线观看| 亚洲精品在线美女| 69精品国产乱码久久久| 成年av动漫网址| 久9热在线精品视频| 韩国高清视频一区二区三区| 婷婷色麻豆天堂久久| 久久久久国产精品人妻一区二区| 十八禁人妻一区二区| 亚洲欧美精品综合一区二区三区| 日韩一本色道免费dvd| 777久久人妻少妇嫩草av网站| 交换朋友夫妻互换小说| 97精品久久久久久久久久精品| 精品福利永久在线观看| 精品福利观看| 国产成人影院久久av| 99热全是精品| 制服诱惑二区| 天天躁夜夜躁狠狠躁躁| 日本vs欧美在线观看视频| 欧美+亚洲+日韩+国产| 中文字幕人妻熟女乱码| 好男人电影高清在线观看| 婷婷色综合大香蕉| 男人舔女人的私密视频| 99热国产这里只有精品6| 欧美 日韩 精品 国产| 婷婷色综合大香蕉| 亚洲欧洲日产国产| 亚洲精品一区蜜桃| 成年女人毛片免费观看观看9 | 国产成人av教育| 丰满饥渴人妻一区二区三| 最近最新中文字幕大全免费视频 | 青春草亚洲视频在线观看| 国产欧美日韩一区二区三区在线| 日韩av在线免费看完整版不卡| 一本大道久久a久久精品| 99国产精品免费福利视频| 国产精品三级大全| 国产精品久久久久久精品古装| 女性生殖器流出的白浆| 久久久久视频综合| 美国免费a级毛片| 日本黄色日本黄色录像| 美女扒开内裤让男人捅视频| 亚洲伊人色综图| 中国美女看黄片| av国产久精品久网站免费入址| 精品国产一区二区三区四区第35| 日本色播在线视频| 午夜两性在线视频| 日本a在线网址| 欧美日韩视频高清一区二区三区二| 老司机靠b影院| 黄色视频不卡| 久久影院123| 老汉色∧v一级毛片| 久久中文字幕一级| 欧美日韩亚洲国产一区二区在线观看 | 黄片播放在线免费| tube8黄色片| 精品熟女少妇八av免费久了| 亚洲精品国产区一区二| 亚洲精品日韩在线中文字幕| 热99国产精品久久久久久7| 夫妻性生交免费视频一级片| 日日爽夜夜爽网站| 青春草视频在线免费观看| 国产欧美日韩综合在线一区二区| 成年人午夜在线观看视频| 高清视频免费观看一区二区| 人人妻人人澡人人看| 韩国高清视频一区二区三区| 国产精品久久久久久人妻精品电影 | 悠悠久久av| 精品福利永久在线观看| 男男h啪啪无遮挡| 成人午夜精彩视频在线观看| 七月丁香在线播放| 国产又爽黄色视频| 久久ye,这里只有精品| 中文字幕人妻丝袜一区二区| 亚洲av片天天在线观看| 国产成人精品久久二区二区91| 9热在线视频观看99| a级毛片在线看网站| 性色av一级| 99热国产这里只有精品6| 十八禁高潮呻吟视频| 777米奇影视久久| 亚洲人成网站在线观看播放| 亚洲欧美日韩高清在线视频 | 亚洲成人国产一区在线观看 | 欧美日韩亚洲综合一区二区三区_| 老司机亚洲免费影院| 日韩av在线免费看完整版不卡| 超碰97精品在线观看| 亚洲欧洲国产日韩| 成人国语在线视频| 久久人人爽人人片av| 日韩一本色道免费dvd| 热re99久久精品国产66热6| 免费观看a级毛片全部| 乱人伦中国视频| 中文精品一卡2卡3卡4更新| 一区二区日韩欧美中文字幕| 久久人人爽av亚洲精品天堂| 午夜福利在线免费观看网站| 一边摸一边抽搐一进一出视频| 亚洲精品自拍成人| 美女视频免费永久观看网站| 亚洲免费av在线视频| 国产xxxxx性猛交| 制服人妻中文乱码| 黄色一级大片看看| 热99国产精品久久久久久7| 性高湖久久久久久久久免费观看| 亚洲精品国产av成人精品| a级毛片在线看网站| 欧美亚洲日本最大视频资源| 老司机午夜十八禁免费视频| 国产主播在线观看一区二区 | 啦啦啦 在线观看视频| 一区二区三区乱码不卡18| 国产福利在线免费观看视频| 国产精品国产三级国产专区5o| 欧美人与性动交α欧美精品济南到| www.自偷自拍.com| 婷婷丁香在线五月| 亚洲成人国产一区在线观看 | 大型av网站在线播放| 国产精品一二三区在线看| 99精品久久久久人妻精品| 久久久国产欧美日韩av| 亚洲欧洲精品一区二区精品久久久| www.999成人在线观看| 女人爽到高潮嗷嗷叫在线视频| 男女下面插进去视频免费观看| 国产视频一区二区在线看| 日本一区二区免费在线视频| 999久久久国产精品视频| av又黄又爽大尺度在线免费看| 国产91精品成人一区二区三区 | 亚洲av电影在线观看一区二区三区| 女性生殖器流出的白浆| 亚洲午夜精品一区,二区,三区| 欧美精品一区二区免费开放| 91成人精品电影| 国产高清videossex| 日韩,欧美,国产一区二区三区| 中文字幕高清在线视频| 51午夜福利影视在线观看| 一区二区三区激情视频| 国产一区二区在线观看av| 欧美在线黄色| 我的亚洲天堂| 国产有黄有色有爽视频| 国产精品国产三级国产专区5o| 国产麻豆69| 久久久精品免费免费高清| 国产精品一区二区在线不卡| 午夜福利免费观看在线| 中文乱码字字幕精品一区二区三区| 这个男人来自地球电影免费观看| 999精品在线视频| 一区在线观看完整版| 欧美 日韩 精品 国产| 亚洲美女黄色视频免费看| 91精品国产国语对白视频| 中文字幕亚洲精品专区| 丝袜人妻中文字幕| 97人妻天天添夜夜摸| 亚洲国产欧美一区二区综合| 久久精品久久久久久久性| 日本色播在线视频| 精品人妻一区二区三区麻豆| 久久久久久亚洲精品国产蜜桃av| 国产xxxxx性猛交| 一级毛片电影观看| 亚洲久久久国产精品| 中文字幕高清在线视频| 51午夜福利影视在线观看| 久久久久国产精品人妻一区二区| 日本wwww免费看| 天天躁夜夜躁狠狠躁躁| 亚洲av成人精品一二三区| 亚洲成人免费av在线播放| 中文精品一卡2卡3卡4更新| 欧美变态另类bdsm刘玥| 国产高清国产精品国产三级| 久久毛片免费看一区二区三区| 在线观看国产h片| 欧美人与善性xxx| 国产免费一区二区三区四区乱码| 久久人妻熟女aⅴ| 亚洲欧美精品综合一区二区三区| 成年人黄色毛片网站| 自线自在国产av| 成人国产av品久久久| 1024香蕉在线观看| 黄色a级毛片大全视频| www.av在线官网国产| av不卡在线播放| 这个男人来自地球电影免费观看| 999精品在线视频| av福利片在线| 搡老岳熟女国产| 国产成人啪精品午夜网站| 91麻豆精品激情在线观看国产 | 脱女人内裤的视频| 黄色视频不卡| 美国免费a级毛片| 人成视频在线观看免费观看| 99热全是精品| 黄色片一级片一级黄色片| 少妇精品久久久久久久| 午夜激情久久久久久久| 黄片小视频在线播放| 亚洲国产欧美网| 国产成人精品久久二区二区91| 国产日韩欧美视频二区| 成人手机av| 一级黄色大片毛片| 成年人黄色毛片网站| 精品亚洲乱码少妇综合久久| 大码成人一级视频| 在线观看免费视频网站a站| 国产精品一区二区免费欧美 | 国产女主播在线喷水免费视频网站| 韩国高清视频一区二区三区| 国产成人av教育| 日韩电影二区| 黄色怎么调成土黄色| 人人妻人人澡人人爽人人夜夜| 美女福利国产在线| 高清黄色对白视频在线免费看| 麻豆乱淫一区二区| 香蕉国产在线看| 黑人巨大精品欧美一区二区蜜桃| 国产一区亚洲一区在线观看| 80岁老熟妇乱子伦牲交| 午夜激情久久久久久久| 国产精品久久久人人做人人爽| 50天的宝宝边吃奶边哭怎么回事| 91精品三级在线观看| netflix在线观看网站| 欧美在线一区亚洲| 黄网站色视频无遮挡免费观看| 18禁黄网站禁片午夜丰满| 国产日韩欧美亚洲二区| 中文精品一卡2卡3卡4更新| 久久亚洲国产成人精品v| 青春草亚洲视频在线观看| 久久精品人人爽人人爽视色| 一本综合久久免费| 人成视频在线观看免费观看| 热re99久久精品国产66热6| 国产麻豆69| 国产主播在线观看一区二区 | 欧美另类一区| 久久久国产欧美日韩av| 成年女人毛片免费观看观看9 | 成人午夜精彩视频在线观看| xxxhd国产人妻xxx| 欧美日韩亚洲综合一区二区三区_| 日韩 亚洲 欧美在线| 别揉我奶头~嗯~啊~动态视频 | 1024香蕉在线观看|