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

    優(yōu)質(zhì)抗旱小麥種質(zhì)的篩選及功能基因檢測

    2025-01-27 00:00:00薛振宇張康康張?jiān)?/span>閆強(qiáng)強(qiáng)姚立蓉張宏孟亞雄司二靜李葆春馬小樂王化俊汪軍成

    摘要:為有效篩選優(yōu)異抗旱小麥種質(zhì),以163份小麥種質(zhì)資源為材料,對其高分子量麥谷蛋白亞基組成、萌發(fā)期和苗期的抗旱性進(jìn)行分析,篩選優(yōu)質(zhì)抗旱種質(zhì)資源,并對優(yōu)質(zhì)抗旱種質(zhì)資源的已知功能基因/位點(diǎn)進(jìn)行檢測。結(jié)果表明,供試材料中鑒定出包含優(yōu)質(zhì)面條亞基組合1/7+8/5+10和N/14+15/2+12的材料分別有11和1份,包含優(yōu)質(zhì)饅頭亞基組合1/7+9/5+10和N/7+8/2+12的材料分別有5和53份。將萌發(fā)期和苗期抗旱性綜合評價(jià)結(jié)果聯(lián)合后,發(fā)現(xiàn)有14份抗旱性強(qiáng)的小麥種質(zhì)。此外,結(jié)合14份抗旱種質(zhì)的HMW-GS鑒定結(jié)果最終篩選出9份優(yōu)質(zhì)抗旱種質(zhì)資源,分別為中41-42、隴育5號、隴麥079、蘭天31號、92-47、綿麥367、西科麥531號、銅麥3號和臨旱5322。經(jīng)檢測,9份優(yōu)質(zhì)抗旱材料均含有提高粒重基因以及抗赤霉病、條銹病和葉銹病的基因/位點(diǎn),6份材料含有提高多酚氧化酶活性的基因Ppo2-D1,6份材料含有半矮化基因RHT-8,僅隴麥079含有抗旱基因TaSINA-2。以上結(jié)果為小麥育種提供了優(yōu)異種質(zhì)資源與必要的理論支撐。

    關(guān)鍵詞:小麥;HMW-GS;干旱脅迫;功能基因/位點(diǎn);優(yōu)質(zhì)

    doi:10.13304/j.nykjdb.2024.0200

    中圖分類號:S512

    文獻(xiàn)標(biāo)志碼:A

    文章編號:1008?0864(2025)01?0035?15

    小麥?zhǔn)侵匾募Z食作物之一[1],隨著人們生活水平的不斷提高,對小麥品質(zhì)有了更高的要求。小麥品質(zhì)包含加工品質(zhì)、外觀品質(zhì)和營養(yǎng)品質(zhì),其中高分子量麥谷蛋白亞基(high molecular weightglutenin subunit,HMW-GS)組成對小麥加工品質(zhì)發(fā)揮著關(guān)鍵作用[2]。朱保磊等[3]對182份不同品質(zhì)小麥品種的HMW-GS 組成進(jìn)行分析發(fā)現(xiàn),不同HMW-GS 組合對小麥品質(zhì)具有顯著影響,含有1/7+8/5+10亞基組合的小麥品種表現(xiàn)出較優(yōu)的品質(zhì)指標(biāo)[4],含有1/7+9/2+12 和N/14+15/2+12 亞基組合的小麥品種品質(zhì)指標(biāo)較差。張金乾等[5]對隴東旱塬小麥資源的HMW-GS組成進(jìn)行研究,發(fā)現(xiàn)這些種質(zhì)資源亞基變異類型豐富,存在多個優(yōu)質(zhì)亞基可供育種工作者開發(fā)利用。因此,對小麥高分子量麥谷蛋白亞基進(jìn)行分析對小麥品質(zhì)育種中至關(guān)重要。

    據(jù)統(tǒng)計(jì),我國每年有超過2 000萬hm2的耕地受到干旱的影響,淡水資源的短缺使得這一形勢更加嚴(yán)峻[6]。甘肅地處我國西北地區(qū),干旱氣候頻發(fā),這對甘肅省糧食生產(chǎn)造成了巨大威脅,為應(yīng)對干旱對小麥生產(chǎn)的威脅,人們應(yīng)用多種方法對小麥抗旱性進(jìn)行了評價(jià),篩選出了抗旱性強(qiáng)的小麥品種[7?8]。郭瑩等[9]對引自CIMMYT的小麥種質(zhì)進(jìn)行抗旱性篩選,最終篩選出9份抗旱種質(zhì)材料用于今后開展小麥育種工作。但是目前關(guān)于小麥多時期抗旱性篩選的研究較少。

    分子標(biāo)記是以生物基因組上分布的每個標(biāo)記為靶標(biāo),分析特定位點(diǎn)和區(qū)域的基因型變異,即基因型檢測,以全面反映整個基因組的遺傳變化[10]?;蛐蜋z測具有精準(zhǔn)性、高效性、可預(yù)測性、指導(dǎo)農(nóng)業(yè)生產(chǎn)等優(yōu)勢,對小麥育種、遺傳改良和病蟲害防治等方面都具有重要的應(yīng)用價(jià)值。任立世等[11]利用與小麥穗發(fā)芽相關(guān)的標(biāo)記對我國107份歷史品種的穗發(fā)芽抗性進(jìn)行了鑒定,發(fā)現(xiàn)5份抗穗發(fā)芽基因型材料。此外,甘肅地區(qū)小麥病害頻發(fā),對該地區(qū)的小麥生產(chǎn)構(gòu)成了巨大威脅,抗病品種的選育刻不容緩。

    基于此,本研究對163份小麥種質(zhì)材料進(jìn)行萌發(fā)期和苗期的抗旱性綜合評價(jià),結(jié)合這些材料的HMW-GS組分分析結(jié)果,篩選出優(yōu)質(zhì)抗旱的小麥種質(zhì),最后利用靶向測序基因分型技術(shù)對篩選的優(yōu)質(zhì)抗旱小麥種質(zhì)進(jìn)行已知功能基因/位點(diǎn)的檢測,明確這些優(yōu)質(zhì)抗旱種質(zhì)攜帶的已知功能基因/位點(diǎn),為甘肅省優(yōu)質(zhì)抗旱小麥育種提供優(yōu)異種質(zhì)資源和數(shù)據(jù)支撐。

    1 材料與方法

    1.1 試驗(yàn)材料

    本研究所使用的163份小麥種質(zhì)資源由甘肅農(nóng)業(yè)大學(xué)省部共建干旱生境作物學(xué)國家重點(diǎn)實(shí)驗(yàn)室麥類作物育種課題組保存并提供,其中來自甘肅省106份、四川省17份、山西省7份、陜西省8份、河南省2份、江蘇省3份、北京市1份、山東省1份、河北省1份、17份來源未知(表1)。

    1.2 高分子量麥谷蛋白亞基SDS-PAGE 電泳

    小麥高分子量麥谷蛋白亞基(HMW-GS)組分分析為蛋白質(zhì)提取、蛋白質(zhì)電泳、凝膠染色、脫色和成像等過程,具體步驟參考高正等[12]所使用的方法。以已知HMW-GS帶型的小麥材料中國春(N/7+8/2+12)和馬奎斯(1/7+8/5+10)為參考,鑒定供試材料的亞基類型,根據(jù)Payne 等[13]的方法記錄,并進(jìn)行品質(zhì)評價(jià)。

    1.3 種質(zhì)萌發(fā)期抗旱性鑒定

    分別設(shè)置了質(zhì)量體積分?jǐn)?shù)10%、15%、20%和25%的PEG-6000處理進(jìn)行預(yù)試驗(yàn),篩選出用于抗旱性鑒定最適的用量為15%PEG-6000。選取大小一致、飽滿均勻的種子30粒,置于鋪有2層濾紙的發(fā)芽盒中(長12 cm、寬12 cm、高5 cm)澆灌10 mL15%的PEG-6000,以澆灌等量蒸餾水處理為CK,分別設(shè)置3次重復(fù)。然后置于20 ℃的光照培養(yǎng)箱中培養(yǎng),連續(xù)培養(yǎng)7 d,統(tǒng)計(jì)發(fā)芽勢(germinationpotential,GP)、發(fā)芽率(germination rate,GR)。每個重復(fù)選5株幼苗,測量苗長(seedling length,SL)和主根長(root length,RL)、根鮮重(root fresh weight,RFW)、苗鮮重(seedling fresh weight,SFW)、根干重(root dry weight,RDW)、苗干重(seedling dryweight,SDW)等指標(biāo),以上指標(biāo)的測定均參考朱麗麗等[14]的方法。

    1.4 種質(zhì)苗期抗旱性鑒定

    選取大小一致、飽滿均勻的種子40粒置于培養(yǎng)皿中進(jìn)行培養(yǎng),7 d后選擇長勢一致的幼苗,在水培箱中培養(yǎng)至三葉一心期后進(jìn)行處理。用Hoagland營養(yǎng)液配制15% PEG-6000溶液模擬干旱條件,對照組使用Hoagland 營養(yǎng)液,每個處理設(shè)3個重復(fù),每個重復(fù)選3株幼苗。培養(yǎng)10 d后對株高(plant height,PH)和主根長(taproot length,TL)進(jìn)行測量,利用根系掃描儀(Epson, LongBeach, CA, USA)獲取總根長(total root length,TRL)、根體積(root volume,RV)、根表面積(rootsurface area,RSA)等數(shù)據(jù)。并測定根冠比(rootcrown ratio,RCR)、根鮮重和苗鮮重,幼苗烘干至恒重后,稱量苗干重和根干重等指標(biāo)[15],根據(jù)以下公式計(jì)算苗失水率(seedling water loss rate,SWLR)、根失水率(rate of root water loss,RWLR)。

    SWLR=1 - (地上部干重/地上部鮮重)× 100% (1)

    RWLR=1 - (地下部干重/地下部鮮重)× 100% (2)

    1.5 靶向測序基因型分型技術(shù)

    采用石家莊博瑞迪生物技術(shù)有限公司GenoBaits技術(shù),對篩選出的9份優(yōu)質(zhì)抗旱小麥種質(zhì)資源的功能基因/位點(diǎn)進(jìn)行檢測。具體方法如下:首先采用高通量DNA提取試劑盒提取各樣品DNA,并使用1%瓊脂糖凝膠電泳方法分析DNA的純度和完整性,使用Qubit對DNA含量進(jìn)行準(zhǔn)確定量。取質(zhì)檢合格的DNA利用超聲波破碎儀進(jìn)行隨機(jī)物理破碎,對破碎后的小片段進(jìn)行末端修復(fù)、連接Ploy-A、連接測序接頭、純化、PCR擴(kuò)增等過程完成文庫構(gòu)建,文庫檢測合格后上機(jī)測序。待測序完成后使用軟件fastp(version 0.20)對raw reads進(jìn)行過濾;使用軟件BWA(mem 比對方式)將進(jìn)行質(zhì)控后的clean reads與參考基因組(中國春V1.0)進(jìn)行比對,通過比對可以定位cleanreads 在參考基因組上的位置;使用軟件GATK(Broad Institute,4.0)檢測SNP變異。使用基因組注釋工具ANNOVAR[16]對檢測出的基因變異進(jìn)行功能注釋。

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

    采用Excel 2016 統(tǒng)計(jì)數(shù)據(jù),參考宋璐杏等[17]的方法計(jì)算抗旱系數(shù);利用SPSS 25.0 進(jìn)行萌發(fā)期、萌發(fā)期各指標(biāo)相關(guān)性、主成分分析和t 檢驗(yàn);利用R語言進(jìn)行隸屬函數(shù)計(jì)算、聚類分析。

    2 結(jié)果與分析

    2.1 種質(zhì)資源HMW-GS 組成分析

    由表1和表2可知,供試材料的Glu-A1/B1/D1位點(diǎn)上共檢測到10 種高分子量麥谷蛋白亞基。在Glu-A1位點(diǎn)檢測到N、1亞基類型,出現(xiàn)頻率分別為58.28% 和41.72%;Glu-B1 位點(diǎn)檢測到7+8(67.48%)、7+9(30.06%)、13+16(0.61%)、17+18(1.23%)、14+15(0.61%)五種亞基類型;Glu-D1位點(diǎn)檢測到5+10(22.70%)、2+12(73.01%)、2+10(4.29%)3種亞基類型。亞基組成類型較多說明材料來源廣泛,具有良好的統(tǒng)計(jì)意義。供試材料中共有15種亞基組合,鑒定出存在優(yōu)質(zhì)饅頭亞基組合1/7+9/5+10(3.07%)和N/7+8/2+12(32.52%)各有5份和53份。存在優(yōu)質(zhì)面條亞基組合1/7+8/5+10(6.75%)和N/14+15/2+12(0.61%)的材料分別為中41-42、隴麥079、蘭天48 號、蘭10-85、92-47、綿麥367、西科麥510、西科麥531 號、銅麥3號、川麥1247、川麥1546和蘭天37。由以上結(jié)果可知,供試材料中高分子量麥谷蛋白亞基類型豐富,可為小麥優(yōu)質(zhì)品種選育提供數(shù)據(jù)支撐。

    2.2 種質(zhì)資源萌發(fā)期和苗期抗旱性分析

    2.2.1 小麥種質(zhì)性狀分析

    在對照和PEG-6000模擬中度干旱脅迫條件下調(diào)查163 份小麥種質(zhì)資源萌發(fā)期和苗期抗旱相關(guān)指標(biāo)(表3,圖1),發(fā)現(xiàn)各萌發(fā)指標(biāo)和苗期相關(guān)指標(biāo)在模擬中度干旱脅迫下均受到抑制,且干旱脅迫對各萌發(fā)指標(biāo)和苗期指標(biāo)的抑制程度均有差異。此外,發(fā)現(xiàn)在小麥萌發(fā)期苗鮮重受到干旱脅迫的抑制最明顯,在小麥苗期總根長受到干旱脅迫的抑制最嚴(yán)重,因此在評價(jià)小麥抗旱性時可將這2 個指標(biāo)作為評價(jià)依據(jù)。

    2.2.2 小麥抗旱指標(biāo)主成分分析

    對小麥萌發(fā)期和苗期的抗旱系數(shù)進(jìn)行主成分分析(表4),對萌發(fā)期抗旱指標(biāo)分析后共抽取4個主成分,其貢獻(xiàn)率分別為47.1%、22.1%、13.3%和6.9%,累計(jì)貢獻(xiàn)率為89.4%。將苗期抗旱指標(biāo)綜合為3個主成分,貢獻(xiàn)率分別為51.800%、31.900%和3.700%,累計(jì)貢獻(xiàn)率為95.400%。

    基于綜合評價(jià)D 值對供試材料進(jìn)行聚類分析,在萌發(fā)期供試材料被分為3類,第Ⅰ類為干旱敏感型小麥種質(zhì),D 值介于0.277~0.393,占供試材料的41.1%;第Ⅱ類為中等抗旱型小麥種質(zhì),D 值介于0.397~0.518,占供試材料的46.62%;第Ⅲ類為抗旱型小麥種質(zhì),D 值介于0.521~0.550,占供試材料的12.26%。

    基于綜合評價(jià)D 值對供試材料進(jìn)行聚類分析,在苗期供試材料被分為3 類,第Ⅰ類為中等抗旱種質(zhì)共106 份,D 值介于0.540~0.689,占供試材料的65.030%;第Ⅱ類為強(qiáng)抗旱種質(zhì)共30份,D 值介于0.699~0.820,占供試材料的18.400%;第Ⅲ類為弱抗旱種質(zhì)共27份,D 值介于0.166~0.538,占供試材料的16.56%,干旱敏感型材料較少。

    2.3 優(yōu)質(zhì)抗旱小麥種質(zhì)篩選

    結(jié)合各材料萌發(fā)期、苗期抗旱性鑒定結(jié)果以及HMW-GS組分鑒定結(jié)果篩選優(yōu)質(zhì)抗旱小麥材料,如圖3韋恩圖所示,共篩出9份種質(zhì)資源,分別為中41-42、隴育5號、隴麥079、蘭天31號、92-47、綿麥367、西科麥531號、銅麥3號和臨旱5322。9份材料在萌發(fā)期和苗期都表現(xiàn)強(qiáng)抗旱性,含有的亞基類型為1、7+8、7+9、5+10、2+12,其中隴育5號和蘭天31號亞基組合均為7+8/2+12、臨旱5322亞基組合為1/7+9/5+10、中41-42、隴麥079、92-47、綿麥367、西科麥531號和銅麥3號亞基組合均為1/7+8/5+10,這些亞基組合對饅頭和面條品質(zhì)影響顯著。

    2.4 優(yōu)質(zhì)抗旱小麥功能基因/位點(diǎn)檢測

    2.4.1 靶向測序基因型分型數(shù)據(jù)質(zhì)控和SNP位點(diǎn)統(tǒng)計(jì)

    對9份優(yōu)質(zhì)抗旱小麥材料進(jìn)行基因型分型,結(jié)果表明,9份材料過濾后read介于54 203 470~27 028 686個,Q20(質(zhì)控后質(zhì)量大于20)在98.32~98.81,Q30(質(zhì)控后質(zhì)量大于30)在95.74~96.51。進(jìn)一步與參考基因組進(jìn)行比對,比對read數(shù)最多的為臨旱5322,有45 281 639 個,比對率為84.71%,最少的為蘭天31號,有22 371 675個,比對率為82.77%,其中Q20、Q30的樣本堿基數(shù)目占質(zhì)控后樣本總堿基數(shù)目的比值。

    對9份優(yōu)質(zhì)抗旱材料SNP位點(diǎn)進(jìn)行統(tǒng)計(jì),結(jié)果(圖4)如下:位于基因間區(qū)的SNP位點(diǎn)個數(shù)最多,有53 048個,占比48.62%;其次為位于外顯子編碼區(qū)域的核心SNP 位點(diǎn),有21 956 個,占比20.12%;SNP位點(diǎn)位于基因5’端UTR區(qū)域的個數(shù)最少,有869個,占比0.80%。

    2.4.2 優(yōu)質(zhì)抗旱種質(zhì)已知功能基因/位點(diǎn)檢測

    由9份優(yōu)質(zhì)抗旱種質(zhì)材料基因鑒定結(jié)果(表6)可知,所有材料都含有抗赤霉病位點(diǎn)、抗條銹病基因及位點(diǎn)、抗葉銹病基因、抗白粉病基因、大粒籽粒基因、提高籽粒重量的基因、矮桿基因或位點(diǎn),提高多酚氧化酶活性的基因只有隴育5號與蘭天31號沒有鑒定到,隴麥079 中鑒定到了抗旱基因TaSINA-2。檢測出的這些與抗逆、品質(zhì)、抗病相關(guān)的功能基因/位點(diǎn),可為甘肅省優(yōu)良小麥新品種的選育提供優(yōu)異種質(zhì)信息。

    3 討論

    前人對作物的抗旱性已采用多種方法進(jìn)行了研究,但大多都采用的是單一時期,如陳春舟等[18]對16份小麥種質(zhì)材料使用PEG-6000進(jìn)行了苗期抗旱性試驗(yàn),篩選出抗旱性品種40IBWSN931。張樹林等[19]利用主成分分析和聚類分析對萌發(fā)期小麥抗旱性進(jìn)行評價(jià),篩選出可用于抗旱育種的親本材料。不同生育時期小麥抗旱性有所差異[20],因此只進(jìn)行單一時期抗旱性鑒定具有片面性。本研究則采用多指標(biāo)結(jié)合多種分析方法對163份小麥種質(zhì)的抗旱性從萌發(fā)期和苗期兩個時期進(jìn)行了評價(jià),以避免研究單一時期造成的片面性,使鑒定結(jié)果更加可靠。此外,對比2個時期抗旱性結(jié)果,大部分種質(zhì)在苗期和萌發(fā)期所表現(xiàn)的抗旱性有所差異,這與孟雨等[21]研究結(jié)果一致。

    關(guān)于HMW-GS 與小麥品質(zhì)的關(guān)系前人已進(jìn)行了報(bào)道[22]。范家霖等[23]對黃淮麥區(qū)的127份小麥品種(系)進(jìn)行分析后發(fā)現(xiàn)1、7+9、5+10 亞基在Glu-1 位點(diǎn)上出現(xiàn)的頻率最高,本研究的結(jié)果得出7+8和2+12出現(xiàn)的頻率最高,這與前人的研究結(jié)果有所區(qū)別,可能與供試材料不同有關(guān)。此外,本研究共鑒定出15種亞基組合,1/7+9/5+10、N/7+8/2+12、1/7+8/5+10、N/14+15/2+12 亞基組合比為42.95%,可見供試材料中具有優(yōu)良加工品質(zhì)的小麥種質(zhì)占比較高。

    目前,已鑒定出許多與小麥抗病、抗逆、品質(zhì)相關(guān)的基因/功能位點(diǎn)。本研究篩選出的9份優(yōu)異種質(zhì)均含有赤霉病、條銹病、葉銹病和白粉病抗性基因和粒重基因,其中以抗葉銹病基因Lr67和抗白粉病基因Pm12 分布最廣;有7份材料含有多酚氧化酶活性基因/位點(diǎn)、粒長和株高基因/位點(diǎn);含有抗旱基因的種質(zhì)只有1份,其余抗旱材料中未檢測出TaSINA_2 基因,這一問題的出現(xiàn)可能是由于品種的抗旱位點(diǎn)不存在、抗旱位點(diǎn)變異、環(huán)境條件影響等原因所導(dǎo)致,因此后續(xù)需對其進(jìn)行進(jìn)一步探討。隨著小麥生產(chǎn)要求的進(jìn)一步改變,科學(xué)家已經(jīng)應(yīng)用了多種功能基因/位點(diǎn),QFhb.caas-4BS、QFhb.caas-4DS[24]、Yr29[25]、QYr、nwafu-3BS、Yr30、Yr78[26]、Lr46、Lr67、Lr68[27]、Pm12[28]、TaCwi-A1[29]、TaGS5-A1[30]、RHT-8[31]等已廣泛應(yīng)用于小麥育種實(shí)踐中。因此,本研究中篩選出的9份優(yōu)異種質(zhì)可為甘肅省優(yōu)良小麥品種選育提供支持。

    參考文獻(xiàn)

    [1] 蔣赟,張麗麗,薛平,等.我國小麥產(chǎn)業(yè)發(fā)展情況及國際經(jīng)驗(yàn)借鑒[J].中國農(nóng)業(yè)科技導(dǎo)報(bào),2021,23(7):1-10.

    JIANG Y, ZHANG L L, XUE P, et al .. Development status ofwheat industry in China andinternational experience forreference [J]. J. Agric. Sci. Technol., 2021,23(7):1-10.

    [2] 康志鈺,王建軍,尚勛武.關(guān)于拉面品質(zhì)的小麥高分子量谷蛋白亞基評分系統(tǒng)研究[J].中國農(nóng)業(yè)科學(xué),2006,39(12):2415-2421.

    KANG Z Y, WANG J J, SHANG X W. Score system study forhand-extended noodle quality based on HMW-GS index inwheat flour [J]. Sci. Agric. Sin., 2006,39(12):2415-2421.

    [3] 朱保磊,陳宏,尹志剛,等.高分子量麥谷蛋白亞基組成及其對品質(zhì)性狀的影響[J/OL].分子植物育種,2024,1-8 [2024-06-07].https://link.cnki.net/urlid/46.1068.S.20230202.1107.005.html.

    ZHU B L, CHEN H, YIN Z G, et al .. Composition of highmolecular weight glutenin subunits and their effects on qualitytraits [J]. Molecular Plant Breeding, 2024: 1-8 [2024-06-07].https://link.cnki.net/urlid/46.1068.S.20230202.1107.005.html.

    [4] 陳昱利,楊平,李華偉.種植密度和播期對小麥籽粒品質(zhì)影響的模擬研究[J].中國農(nóng)業(yè)科技導(dǎo)報(bào),2022,24(10):143-153.

    CHEN Y L, YANG P, LI H W. Simulation study on the effectsof planting density and sowing date on grain quality of winterwheat [J]. J. Agric. Sci. Technol., 2022,24(10):143-153.

    [5] 張金乾,張成,劉自成,等.隴東旱塬小麥資源的HMW-GS組成及品質(zhì)比較[J].干旱地區(qū)農(nóng)業(yè)研究,2016,34(2):146-151.

    ZHANG J Q, ZHANG C, LIU Z C, et al .. Compositions ofHMW-GS and quality properties of winter wheat varieties inLongdong drought area [J]. Agric. Res. Arid Areas, 2016,34(2):146-151.

    [6] 傅曉藝,王紅光,劉志連,等.水分脅迫對不同小麥幼苗期生長的影響及抗旱品種篩選[J].作物雜志,2023,1(4):224-229.

    FU X Y, WANG H G, LIU Z L, et al .. Effect of water stress onthe growt of different wheat seedlings and screening of droughtresistantvarieties [J]. Crops, 2023,1(4): 224-229.

    [7] 劉海,姜亮亮,劉冰,等.近40年中國干旱特征及其對植被變化的影響[J].生態(tài)學(xué)報(bào),2023,43(19):7936-7949.

    LIU H, JIANG L L, LIU B, et al .. Characteristics of drought inChina and its effect on vegetation chang in recent 40 years [J].Acta Ecol. Sinica, 2023,43(19):7936-7949.

    [8] 吳曉軍,姜豪,宋林通,等.黃淮海麥區(qū)部分小麥品種的抗旱性評價(jià)及相關(guān)分子標(biāo)記的有效性驗(yàn)證[J]. 麥類作物學(xué)報(bào),2023,43(3):270-278.

    WU X J, JIANG H, SONG L T, et al .. Comprehensiveevaluation on drought resistance of some wheat cultivars inHuang-Huai-Hai region and validation of the related molecularmarkers [J]. J. Triticeae Crops, 2023,43(3):270-278.

    [9] 郭瑩,王振平,呂迎春,等.不同麥類種質(zhì)抗旱性評價(jià)及優(yōu)異種質(zhì)篩選[J].麥類作物學(xué)報(bào),2022,42(10):1208-1219.

    GUO Y, WANG Z P, LYU Y C, et al .. Evaluation of droughtresistance and excellent germplasm screening in differentTriticeae crops [J]. J. Triticeae Crops, 2022,42(10):1208-1219.

    [10] 徐云碧,楊泉女,鄭洪建,等.靶向測序基因型檢測(GBTS)技術(shù)及其應(yīng)用[J].中國農(nóng)業(yè)科學(xué),2020,53(15):2983-3004.

    XU Y B, YANG Q N, ZHENG H J, et al .. Genotyping by targetsequencing (GBTS) and its applications [J]. Sci. Agric. Sin.,2019,53(15):2983-3004.

    [11] 任立世,劉進(jìn)英,楊燕.107份小麥歷史品種抗穗發(fā)芽基因型的檢測[J].麥類作物學(xué)報(bào),2015,35(6):752-758.

    REN L S, LIU J Y, YANG Y. Detection of Pre-harvestsprouting resistance genotypes with Vp1A3 and Vp1B3 in 107chinese historical wheat cultivars [J]. J. Triticeae Crops, 2015,35(6):752-758.

    [12] 高正,王曉龍,張曉科,等.小麥HMW-GS檢測方法的優(yōu)化及應(yīng)用[J].麥類作物學(xué)報(bào),2013,33(5):929-934.

    GAO Z, WANG X L, ZHANG X K, et al .. Optimization andapplication of the detection method of HMW-GS in commonwheat [J]. J. Triticeae Crop, 2013,33(5):929-934.

    [13] PAYNE P I, NIGHTINGALE M A, KRATTIGER A F, et al .. Therelationship between HMW glutenin subunit composition and thebread-making quality of British-grown wheat varieties [J]. J. Sci.Food Agric., 1987, 40(1): 51-65.

    [14] 朱麗麗,張發(fā)玉,安寧,等.116份藜麥種質(zhì)資源萌發(fā)期抗旱性綜合評價(jià)[J].干旱地區(qū)農(nóng)業(yè)研究,2024,42(1):23-31.

    ZHU L L, ZHANG F Y, AN N, et al .. Comprehensiveevaluation of drought resistance of 116 quinoa germplasmresources during germination [J]. Agric. Res. Arid Areas, 2024,42(1):23-31.

    [15] 王亞茹,楊向東,趙寒冬,等.油莎豆種質(zhì)資源苗期抗旱性鑒定與評價(jià)[J].干旱地區(qū)農(nóng)業(yè)研究,2022,40(6):12-22.

    WANG Y R, YANG X D, ZHAO H D, et al.. Characterization andassessment of drought resistance in tigernut germplasm resources atseedling stage [J]. Agric. Res. Arid Areas, 2022,40(6):12-22.

    [16] WANG K, LI M Y, HAKONARSON H. ANNOVAR: functionalannotation of genetic variants from high-throughput sequencingdata [J]. Nucleic Acids Res., 2010, 38 (16):164-164.

    [17] 宋璐杏,張閃閃,李玹,等.小麥高代系的抗旱性狀篩選與抗旱性評價(jià)[J].干旱地區(qū)農(nóng)業(yè)研究,2024,42(1):14-22.

    SONG L X, ZHANG S S, LI X, et al . Screening of droughtresistance trait and drought resistance evaluation in wheatadvanced lines [J]. Agric. Res. Arid Areas, 2024,42(1):14-22.

    [18] 陳春舟,馬占軍,孟亞雄,等.小麥種質(zhì)資源抗旱耐鹽性評價(jià)及種質(zhì)篩選[J].分子植物育種,2021,19(14):4820-4835.

    CHEN C Z, MA Z J, MENG Y X, et al .. Evaluation andscreening of wheat germplasm resources for drought and salttolerance [J]. Mol. Plant Breeding, 2021,19(14):4820-4835.

    [19] 張樹林,劉玉玲,田麗,等.不同小麥新品系萌發(fā)期抗旱性的篩選與鑒定[J].分子植物育種, 2018, 16(21):7138-7147.

    ZHANG S L, LIU Y L, TIAN L, et al .. Screening andidentification of drought resistance of different new wheat linesduring germination period [J]. Mol. Plant Breeding, 2018, 16(21):7138-7147.

    [20] 張軍,魏國,彭彥珉,等.8份強(qiáng)筋小麥品種抗旱性評價(jià)[J].麥類作物學(xué)報(bào), 2024, 44(4):442-452.

    ZHANG J, WEI G, PENG Y M, et al .. Drought resistanceEvaluation of eight strong gluten wheat varieties [J]. J.Triticeae Crop, 2024, 44(4):442-452.

    [21] 孟雨,田文仲,溫鵬飛,等.基于不同發(fā)育階段協(xié)同的小麥品種抗旱性綜合評判[J].作物學(xué)報(bào), 2023, 49(2):570-582.

    MENG Y, TIAN W Z, WEN P F, et al .. Comprehensiveevaluation of drought resistance of wheat varieties based onsynergr of different developmental stages [J]. Acta Agron. Sin.,2023, 49(2):570-582.

    [22] 張影全,張曉科,魏益民,等.高分子量麥谷蛋白亞基對小麥蛋白質(zhì)品質(zhì)特性的影響[J].西北農(nóng)業(yè)學(xué)報(bào),2013,22(1):48-53.

    ZHANG Y Q, ZHANG X k, WEI Y M, et al .. Effect of highmolecular weight glutenin subunits on kernel protein qualityproperty of wheat [J]. Acta Agric. Boreali-occident. Sin., 2013,22(1): 48-53.

    [23] 范家霖, 陳曉杰, 張建偉, 等. 高分子量麥谷蛋白亞基組成及其與小麥品質(zhì)性狀的關(guān)系分析[J]. 麥類作物學(xué)報(bào), 2021,41(5): 544-552.

    FAN J L, CHEN X J, ZHANG J W, et al .. Composition of highmolecular weight glutenin subunit and theirs relationship withwheat quality traits [J]. J. Triticeae Crop, 2018, 41(5): 544-552.

    [24] ZHU Z, XU X, FU L, et al .. Molecular mapping of quantitativetrait loci for fusarium head blight resistance in a doubledhaploid population of chinese bread wheat [J]. Plant Dis.,2021, 105(5): 1339-1345.

    [25] LIU S J, LIU D, ZHANG C L, et al .. Slow stripe rusting inChinese wheat Jimai 44 conferred by Yr29 in combination witha major QTL on chromosome arm 6AL [J]. Theor. Appl. Genet.,2023, 136(8): 175-175.

    [26] HUANG S, WU J H, WANG X T, et al.. Utilization of thegenomewide wheat 55K SNP array for genetic analysis of stripe rustresistance in common wheat line P9936 [J]. Phytopathology, 2019,109(5): 819-827.

    [27] SOCHALOVA P L, BOYKO I N, POTESHKINA A A, et al ..Effective leaf rust resistance genes of wheat in NovosibirskProvince in connection with the variability of the Pucciniatriticina population [J]. Труды по Прикладной Ботанике,Генетике и Селекции, 2023, 184(2): 235-244.

    [28] TOMASZ O, SYLWIA O. Chromosomal location of Pm12—Anovel powdery mildew resistance gene from avena sterilis [J].Genes, 2022, 13(12): 2409-2409.

    [29] 相吉山,穆培源,桑偉,等.小麥粒重基因TaCwi-A1 功能標(biāo)記CWI22、CWI21的驗(yàn)證及應(yīng)用[J]. 中國農(nóng)業(yè)科學(xué), 2014, 47(13):2671-2709.

    XIANG J S, MU P Y, SANG W, et al .. Validation andapplication of function markers CWI22 and CWI21 of TaCwi-A1 gene related to wheat kernel weight [J]. Sci. Agric. Sin.,2014, 47(13): 2671-2709.

    [30] WANG S S, YAN X F, WANG Y Y, et al .. Haplotypes of theTaGS5-A1 gene are associated with thousand-kernel weight inChinese bread wheat [J/OL]. Front. Plant Sci., 2016,7: 7783[2024-06-07]. https://doi.org/10.3389/fpls.2016.00783.

    [31] 張凱, 蘭素缺, 金京京, 等. 秋播冬麥區(qū)矮稈基因RhtB1b、RhtD1b 和Rht-8 的分布頻率及其與產(chǎn)量性狀的關(guān)系[J]. 麥類作物學(xué)報(bào), 2023,43(4): 399-408.

    ZHANG K, LAN S Q, JIN J J, et al .. Distribution frequency ofdwarf genes RhtB1b, RhtD1b and Rht-8 in autumn-sownwinter wheat areas and their relationship with yield traits [J]. J.Triticeae Crop, 2023, 43(4): 399-408.

    基金項(xiàng)目:2023年甘肅省級重點(diǎn)人才項(xiàng)目(GSAU-2023-03);甘肅省隴原青年英才項(xiàng)目(GSAU-2023-005);國家自然科學(xué)基金項(xiàng)目(31860377);甘肅省教育廳產(chǎn)業(yè)支撐計(jì)劃項(xiàng)目(2021CYZC-12);甘肅省農(nóng)業(yè)科技支撐項(xiàng)目(KJZC-2023-2);甘肅農(nóng)業(yè)大學(xué)伏羲青年英才計(jì)劃項(xiàng)目(Ganfx-03Y06, GAUfx-04Y011)。

    9191精品国产免费久久| 三级毛片av免费| 亚洲久久久久久中文字幕| 18+在线观看网站| 成人av一区二区三区在线看| 精品久久国产蜜桃| 狂野欧美白嫩少妇大欣赏| 蜜桃亚洲精品一区二区三区| 在线十欧美十亚洲十日本专区| 91久久精品电影网| 免费高清视频大片| 日韩亚洲欧美综合| 国产精品野战在线观看| 神马国产精品三级电影在线观看| 国产精品三级大全| 亚洲成av人片免费观看| 精品久久久久久久久久久久久| 欧美日韩国产亚洲二区| 午夜免费成人在线视频| 一夜夜www| 美女高潮的动态| 国内精品久久久久精免费| 精品不卡国产一区二区三区| 亚洲av美国av| 久久九九热精品免费| 精品人妻一区二区三区麻豆 | 夜夜爽天天搞| 级片在线观看| 国产一区二区三区视频了| 国产蜜桃级精品一区二区三区| 欧美zozozo另类| 99久国产av精品| 少妇的逼水好多| 国产欧美日韩精品一区二区| 男人舔女人下体高潮全视频| 精品人妻视频免费看| 搡女人真爽免费视频火全软件 | 18禁裸乳无遮挡免费网站照片| 在现免费观看毛片| 99在线视频只有这里精品首页| 精品熟女少妇八av免费久了| 日本撒尿小便嘘嘘汇集6| 精品日产1卡2卡| 直男gayav资源| 九色国产91popny在线| av专区在线播放| 欧美+日韩+精品| 一个人看的www免费观看视频| 国产一级毛片七仙女欲春2| 波多野结衣高清无吗| 欧美bdsm另类| 99在线人妻在线中文字幕| 69av精品久久久久久| 97超级碰碰碰精品色视频在线观看| 怎么达到女性高潮| 神马国产精品三级电影在线观看| 91在线精品国自产拍蜜月| 亚洲精品456在线播放app | 无遮挡黄片免费观看| 超碰av人人做人人爽久久| 99热6这里只有精品| 脱女人内裤的视频| 99国产精品一区二区三区| 中国美女看黄片| 亚洲18禁久久av| 精品乱码久久久久久99久播| 国产精品,欧美在线| 精品乱码久久久久久99久播| 久久久久久大精品| 99久久无色码亚洲精品果冻| 久久午夜亚洲精品久久| 首页视频小说图片口味搜索| 一卡2卡三卡四卡精品乱码亚洲| 一本精品99久久精品77| 国产激情偷乱视频一区二区| 人妻制服诱惑在线中文字幕| 久久久久免费精品人妻一区二区| 欧美乱色亚洲激情| 国产精品av视频在线免费观看| 男女视频在线观看网站免费| 精品久久久久久久久久免费视频| 神马国产精品三级电影在线观看| 成人午夜高清在线视频| 国产色爽女视频免费观看| 精品国产亚洲在线| 九九热线精品视视频播放| 久久久久性生活片| 日韩国内少妇激情av| 国内少妇人妻偷人精品xxx网站| 国产一区二区激情短视频| 别揉我奶头~嗯~啊~动态视频| 欧美潮喷喷水| 精品人妻视频免费看| 亚洲欧美日韩卡通动漫| av在线老鸭窝| 波多野结衣高清无吗| 超碰av人人做人人爽久久| 精品日产1卡2卡| 国产精品1区2区在线观看.| 他把我摸到了高潮在线观看| 国产精华一区二区三区| 亚洲中文日韩欧美视频| av国产免费在线观看| 99久久无色码亚洲精品果冻| 真实男女啪啪啪动态图| 久久久久亚洲av毛片大全| 偷拍熟女少妇极品色| xxxwww97欧美| 国产老妇女一区| 国产精品,欧美在线| 久久久国产成人免费| 日日摸夜夜添夜夜添小说| 中文字幕av成人在线电影| 黄色视频,在线免费观看| 亚洲欧美精品综合久久99| 一个人看的www免费观看视频| 久久中文看片网| 麻豆成人午夜福利视频| 国产免费男女视频| 美女黄网站色视频| 桃色一区二区三区在线观看| 夜夜躁狠狠躁天天躁| 国产精品伦人一区二区| 我要看日韩黄色一级片| 色哟哟·www| 亚洲aⅴ乱码一区二区在线播放| 亚洲精华国产精华精| 免费观看的影片在线观看| 精品熟女少妇八av免费久了| 国产久久久一区二区三区| 小蜜桃在线观看免费完整版高清| 国产精品久久久久久久电影| 午夜影院日韩av| 天堂动漫精品| 精品久久久久久久久亚洲 | 97超视频在线观看视频| 亚洲五月天丁香| 美女 人体艺术 gogo| 别揉我奶头 嗯啊视频| 欧美三级亚洲精品| 国产精品精品国产色婷婷| 校园春色视频在线观看| 好男人电影高清在线观看| 亚洲av一区综合| 欧洲精品卡2卡3卡4卡5卡区| 不卡一级毛片| 91在线精品国自产拍蜜月| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 69av精品久久久久久| 亚洲久久久久久中文字幕| 亚洲中文字幕一区二区三区有码在线看| 欧美一区二区国产精品久久精品| 久久伊人香网站| 三级国产精品欧美在线观看| 国产高清有码在线观看视频| 啦啦啦观看免费观看视频高清| 日本五十路高清| 久久久久久久精品吃奶| 免费人成在线观看视频色| 精品久久久久久,| 久久午夜福利片| 国产精品一区二区免费欧美| 亚洲专区国产一区二区| 亚洲精品亚洲一区二区| 精品人妻视频免费看| 亚洲激情在线av| 亚洲熟妇中文字幕五十中出| 久久久久久久精品吃奶| 欧美中文日本在线观看视频| 99久久99久久久精品蜜桃| 成人av一区二区三区在线看| 欧美黑人欧美精品刺激| 亚洲精品影视一区二区三区av| 中出人妻视频一区二区| 亚洲午夜理论影院| 日韩人妻高清精品专区| 国产欧美日韩精品亚洲av| 美女cb高潮喷水在线观看| 国产大屁股一区二区在线视频| 网址你懂的国产日韩在线| 美女免费视频网站| av在线天堂中文字幕| 精品久久久久久成人av| 亚洲成人免费电影在线观看| 性色av乱码一区二区三区2| 在线看三级毛片| 真实男女啪啪啪动态图| 深夜a级毛片| 成年女人永久免费观看视频| 欧美丝袜亚洲另类 | 亚洲精品成人久久久久久| 国产精华一区二区三区| 欧美不卡视频在线免费观看| 露出奶头的视频| 精品久久久久久成人av| 精品人妻1区二区| 色精品久久人妻99蜜桃| 午夜视频国产福利| 久久久久久久久久黄片| 欧美极品一区二区三区四区| 在线观看66精品国产| 亚洲专区中文字幕在线| 亚洲国产精品久久男人天堂| 国产亚洲精品av在线| 波多野结衣高清无吗| 国产一区二区三区在线臀色熟女| 成人av在线播放网站| 天堂影院成人在线观看| 国产精品98久久久久久宅男小说| 男女下面进入的视频免费午夜| 一个人免费在线观看电影| xxxwww97欧美| 亚洲,欧美,日韩| 美女免费视频网站| 嫩草影院新地址| 极品教师在线视频| 757午夜福利合集在线观看| 欧美在线黄色| 69人妻影院| 欧美日韩亚洲国产一区二区在线观看| 久久久久国产精品人妻aⅴ院| eeuss影院久久| 欧美zozozo另类| 一进一出抽搐动态| 国产欧美日韩精品一区二区| 免费在线观看影片大全网站| 欧美午夜高清在线| 69av精品久久久久久| 精品午夜福利视频在线观看一区| 国产探花在线观看一区二区| 国产精品乱码一区二三区的特点| 男人舔奶头视频| 欧美性猛交╳xxx乱大交人| 两人在一起打扑克的视频| 老司机深夜福利视频在线观看| 2021天堂中文幕一二区在线观| 欧美一区二区亚洲| 如何舔出高潮| h日本视频在线播放| 久久精品国产亚洲av香蕉五月| 18禁裸乳无遮挡免费网站照片| 夜夜躁狠狠躁天天躁| 久久精品人妻少妇| 综合色av麻豆| 国产色婷婷99| 精品乱码久久久久久99久播| 99久国产av精品| 午夜福利高清视频| 1024手机看黄色片| 亚洲成人免费电影在线观看| 最近最新免费中文字幕在线| 内射极品少妇av片p| 国产av不卡久久| 久久精品国产亚洲av天美| 国产精品久久久久久亚洲av鲁大| 国产真实乱freesex| 欧美最新免费一区二区三区 | 首页视频小说图片口味搜索| 神马国产精品三级电影在线观看| 国产欧美日韩一区二区精品| 精品国产亚洲在线| 99国产精品一区二区蜜桃av| 久久精品国产亚洲av天美| 黄色一级大片看看| 美女免费视频网站| 亚洲无线在线观看| 精品人妻熟女av久视频| 老司机深夜福利视频在线观看| 亚洲片人在线观看| 久久精品国产亚洲av香蕉五月| 91在线精品国自产拍蜜月| 别揉我奶头 嗯啊视频| 国产视频一区二区在线看| 97超视频在线观看视频| 亚洲国产色片| 国产在线男女| 日韩成人在线观看一区二区三区| 天天躁日日操中文字幕| 别揉我奶头 嗯啊视频| 亚洲一区二区三区色噜噜| 高清日韩中文字幕在线| 最后的刺客免费高清国语| 怎么达到女性高潮| 久久欧美精品欧美久久欧美| 超碰av人人做人人爽久久| 精品久久久久久成人av| 成人特级黄色片久久久久久久| 十八禁网站免费在线| 一进一出抽搐动态| 搞女人的毛片| 国产免费av片在线观看野外av| 男人和女人高潮做爰伦理| 99热精品在线国产| av天堂中文字幕网| 我的女老师完整版在线观看| 欧美色视频一区免费| 色5月婷婷丁香| 精品国产三级普通话版| 午夜精品久久久久久毛片777| 国产在线男女| 美女大奶头视频| 欧美成人一区二区免费高清观看| 久久久久久久亚洲中文字幕 | 亚洲最大成人手机在线| 亚洲午夜理论影院| 色噜噜av男人的天堂激情| 国产亚洲精品久久久com| 日韩欧美在线乱码| 五月伊人婷婷丁香| 18美女黄网站色大片免费观看| 赤兔流量卡办理| 在线免费观看的www视频| 欧洲精品卡2卡3卡4卡5卡区| 热99在线观看视频| 亚洲国产精品久久男人天堂| 久久婷婷人人爽人人干人人爱| 特大巨黑吊av在线直播| ponron亚洲| 日本免费a在线| 免费看日本二区| 又黄又爽又免费观看的视频| 久久中文看片网| 99久久精品国产亚洲精品| 婷婷亚洲欧美| 亚洲av成人av| 国产精品伦人一区二区| 全区人妻精品视频| 国产精品自产拍在线观看55亚洲| 精品日产1卡2卡| 1000部很黄的大片| 欧美日韩乱码在线| 最近最新中文字幕大全电影3| 日韩精品青青久久久久久| 亚洲七黄色美女视频| 国产在视频线在精品| 国产一区二区三区视频了| 男女床上黄色一级片免费看| 一级黄片播放器| 精品福利观看| 99在线人妻在线中文字幕| 老司机福利观看| 久久久久久九九精品二区国产| av视频在线观看入口| 亚洲三级黄色毛片| 97人妻精品一区二区三区麻豆| 男人的好看免费观看在线视频| 国产精品亚洲美女久久久| 色综合亚洲欧美另类图片| 国产人妻一区二区三区在| 亚洲国产精品合色在线| 久久中文看片网| 女人被狂操c到高潮| 亚洲第一区二区三区不卡| 51午夜福利影视在线观看| 啦啦啦观看免费观看视频高清| www日本黄色视频网| 精品一区二区三区视频在线| 午夜福利高清视频| 很黄的视频免费| 亚洲精品亚洲一区二区| 久久99热6这里只有精品| 免费看a级黄色片| 久久99热6这里只有精品| 亚洲精品一区av在线观看| 国产日本99.免费观看| 欧美+亚洲+日韩+国产| 18禁黄网站禁片午夜丰满| 好男人电影高清在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲 国产 在线| bbb黄色大片| 国产精品综合久久久久久久免费| 热99re8久久精品国产| 亚洲av成人不卡在线观看播放网| 日本撒尿小便嘘嘘汇集6| 亚洲精品一区av在线观看| 日本黄色片子视频| 中文字幕人成人乱码亚洲影| 国产久久久一区二区三区| 久久99热这里只有精品18| 丰满人妻熟妇乱又伦精品不卡| 久久99热这里只有精品18| 精品久久久久久久久久久久久| 亚洲综合色惰| 亚洲欧美日韩高清在线视频| 精品久久久久久久久av| 亚洲国产日韩欧美精品在线观看| 日本a在线网址| 久久久久国产精品人妻aⅴ院| 久久久久久久精品吃奶| 欧美区成人在线视频| 国产精华一区二区三区| 久久热精品热| 麻豆av噜噜一区二区三区| 成人三级黄色视频| 国产成人啪精品午夜网站| 久久精品国产亚洲av香蕉五月| 99久久成人亚洲精品观看| 高清日韩中文字幕在线| 国产探花在线观看一区二区| 欧美精品国产亚洲| 亚洲精品乱码久久久v下载方式| 国产精品久久久久久久电影| 少妇的逼好多水| 欧美高清性xxxxhd video| 亚洲精品在线观看二区| 日韩高清综合在线| 国产精品亚洲美女久久久| АⅤ资源中文在线天堂| 国产成人啪精品午夜网站| 日韩欧美精品v在线| 亚洲人成电影免费在线| 网址你懂的国产日韩在线| 成年女人永久免费观看视频| 精品人妻一区二区三区麻豆 | 成年人黄色毛片网站| 一个人看的www免费观看视频| 国产精品精品国产色婷婷| 亚洲三级黄色毛片| 久久精品国产清高在天天线| 精品日产1卡2卡| 最好的美女福利视频网| 精品熟女少妇八av免费久了| 亚洲中文字幕日韩| 91久久精品国产一区二区成人| 日韩欧美在线乱码| 波多野结衣高清作品| 午夜久久久久精精品| 热99在线观看视频| 国产精品久久久久久久久免 | 精品久久久久久久末码| 少妇高潮的动态图| 免费av不卡在线播放| 亚洲欧美清纯卡通| 91午夜精品亚洲一区二区三区 | 国产伦精品一区二区三区四那| 午夜福利免费观看在线| 小蜜桃在线观看免费完整版高清| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 日本成人三级电影网站| 亚洲国产欧美人成| 日本一本二区三区精品| 久久这里只有精品中国| 欧美黄色淫秽网站| 波多野结衣高清无吗| 精品一区二区三区视频在线| 日本黄色视频三级网站网址| 啪啪无遮挡十八禁网站| 亚洲va日本ⅴa欧美va伊人久久| 久久久久久大精品| 免费在线观看日本一区| 精品一区二区免费观看| 国产熟女xx| 国产精品,欧美在线| 国内久久婷婷六月综合欲色啪| 久久国产乱子伦精品免费另类| bbb黄色大片| 麻豆久久精品国产亚洲av| 日韩高清综合在线| 久久久久国产精品人妻aⅴ院| 欧美性猛交╳xxx乱大交人| 天堂av国产一区二区熟女人妻| 久久中文看片网| 色尼玛亚洲综合影院| av欧美777| 国产v大片淫在线免费观看| 日本五十路高清| 狂野欧美白嫩少妇大欣赏| av在线蜜桃| 国产亚洲精品久久久com| 麻豆成人av在线观看| 女人十人毛片免费观看3o分钟| 国产免费男女视频| 欧美性感艳星| 国产亚洲精品综合一区在线观看| 亚洲国产欧美人成| 免费在线观看日本一区| 亚洲熟妇中文字幕五十中出| 观看免费一级毛片| 国产精品一及| 99久久成人亚洲精品观看| 婷婷六月久久综合丁香| 亚洲人成网站在线播| 丁香欧美五月| 首页视频小说图片口味搜索| 亚洲精品一区av在线观看| 一级黄片播放器| 91午夜精品亚洲一区二区三区 | 国产视频内射| 最好的美女福利视频网| 国产中年淑女户外野战色| 成人av一区二区三区在线看| 午夜精品在线福利| 精品国产亚洲在线| 午夜福利视频1000在线观看| 国产黄片美女视频| 女人十人毛片免费观看3o分钟| 伊人久久精品亚洲午夜| 亚洲av一区综合| 久久久久久久久久成人| 国内毛片毛片毛片毛片毛片| 色综合婷婷激情| 国产精品精品国产色婷婷| 亚洲午夜理论影院| 我要看日韩黄色一级片| 91在线观看av| 国产黄a三级三级三级人| 成人一区二区视频在线观看| 国产精品一区二区三区四区久久| 国产高清三级在线| 亚洲av美国av| 精品久久久久久久久久免费视频| 99热只有精品国产| 精品一区二区三区av网在线观看| 亚洲自偷自拍三级| 人妻丰满熟妇av一区二区三区| 亚洲国产欧洲综合997久久,| 人妻丰满熟妇av一区二区三区| 国产av不卡久久| 日日摸夜夜添夜夜添小说| 亚洲国产欧美人成| 午夜免费男女啪啪视频观看 | 日本黄大片高清| 欧美色视频一区免费| 一区二区三区高清视频在线| 熟女人妻精品中文字幕| 1000部很黄的大片| 日韩有码中文字幕| 国产爱豆传媒在线观看| 日韩有码中文字幕| 亚洲av成人不卡在线观看播放网| 在线看三级毛片| 久久国产乱子伦精品免费另类| 免费看日本二区| 久久精品国产亚洲av香蕉五月| 久久人妻av系列| 亚洲av免费在线观看| 日韩大尺度精品在线看网址| 精品久久久久久久久av| 赤兔流量卡办理| 18禁在线播放成人免费| 欧美日韩国产亚洲二区| 国产黄色小视频在线观看| 国语自产精品视频在线第100页| 午夜老司机福利剧场| 九九在线视频观看精品| 午夜免费激情av| 欧美日韩黄片免| 成人av在线播放网站| 亚洲av成人不卡在线观看播放网| 99热这里只有是精品在线观看 | 老司机午夜福利在线观看视频| 欧美在线黄色| 国产av一区在线观看免费| 97超视频在线观看视频| 日韩大尺度精品在线看网址| 国产精品日韩av在线免费观看| 9191精品国产免费久久| 国产黄色小视频在线观看| 国产精品98久久久久久宅男小说| 亚洲精品久久国产高清桃花| 动漫黄色视频在线观看| 麻豆成人午夜福利视频| 男人舔奶头视频| 国产精品伦人一区二区| 搡女人真爽免费视频火全软件 | 国产免费av片在线观看野外av| 最近最新中文字幕大全电影3| АⅤ资源中文在线天堂| 人人妻人人看人人澡| 51午夜福利影视在线观看| 又黄又爽又免费观看的视频| 精品久久久久久久末码| 一区二区三区激情视频| 麻豆国产av国片精品| 国产精品久久久久久亚洲av鲁大| 男人狂女人下面高潮的视频| 男人舔女人下体高潮全视频| 中文在线观看免费www的网站| 欧美精品国产亚洲| 91在线观看av| 69av精品久久久久久| 国产精品久久久久久久电影| 精品人妻熟女av久视频| 亚洲自拍偷在线| 直男gayav资源| 亚洲熟妇熟女久久| 99久久99久久久精品蜜桃| 91在线精品国自产拍蜜月| 在线观看免费视频日本深夜| 国产真实乱freesex| 一a级毛片在线观看| 日韩欧美国产在线观看| 国产成人福利小说| 欧美+日韩+精品| 色精品久久人妻99蜜桃| 亚洲一区二区三区色噜噜| x7x7x7水蜜桃| 久久九九热精品免费| av欧美777| 亚洲无线观看免费| 午夜精品久久久久久毛片777| www.色视频.com| 看片在线看免费视频| 日日干狠狠操夜夜爽| 午夜两性在线视频| 最近在线观看免费完整版| av天堂在线播放| 蜜桃久久精品国产亚洲av| 18禁在线播放成人免费| 一进一出抽搐gif免费好疼| 天堂影院成人在线观看| av视频在线观看入口| 99国产极品粉嫩在线观看|