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

    小麥品種揚(yáng)麥16赤霉病抗擴(kuò)展QTL定位及分析

    2020-01-02 01:16:04胡文靜陸成彬王鳳菊劉金棟蔣正寧王金平朱展望徐小婷郝元峰何中虎高德榮
    作物學(xué)報(bào) 2020年2期
    關(guān)鍵詞:檢測(cè)

    胡文靜 張 勇 陸成彬 王鳳菊 劉金棟 蔣正寧 王金平朱展望 徐小婷 郝元峰 何中虎,3 高德榮,*

    小麥品種揚(yáng)麥16赤霉病抗擴(kuò)展QTL定位及分析

    胡文靜1張 勇1陸成彬1王鳳菊2劉金棟2蔣正寧1王金平2朱展望2徐小婷2郝元峰2何中虎2,3高德榮1,*

    1江蘇里下河地區(qū)農(nóng)業(yè)科學(xué)研究所/ 農(nóng)業(yè)農(nóng)村部長(zhǎng)江中下游小麥生物學(xué)與遺傳育種重點(diǎn)實(shí)驗(yàn)室, 江蘇揚(yáng)州 225007;2國(guó)家小麥改良中心 / 中國(guó)農(nóng)業(yè)科學(xué)院作物科學(xué)研究所, 北京 100081;3國(guó)際玉米小麥改良中心(CIMMYT)中國(guó)辦事處, 北京 100081

    揚(yáng)麥系列品種赤霉病抗性在世界范圍內(nèi)得到重視, 但其抗性遺傳機(jī)制尚不清楚。揚(yáng)麥16是近年來(lái)大面積推廣的抗赤霉病品種, 本研究以揚(yáng)麥16與中麥895雜交構(gòu)建的174個(gè)雙單倍體(double haploid lines, DH)系為材料, 于2017—2019年連續(xù)3年對(duì)該群體采用單花滴注進(jìn)行赤霉病抗擴(kuò)展鑒定。利用660K SNP芯片構(gòu)建高密度遺傳圖譜, 共檢測(cè)到6個(gè)抗性QTL, 分別位于2DL、3BL、4BS、4DS、5BL和6AS染色體上。除4BS位點(diǎn)外, 其他5個(gè)抗性等位基因均來(lái)源于揚(yáng)麥16。和均在多年被檢測(cè)到, 可解釋8.8%~15.0%的表型變異;、僅在1年被檢測(cè)到, 分別解釋10.5%和14.7%的表型變異;和來(lái)源于中麥895的僅在1年被檢測(cè)到且效應(yīng)僅為6.4%和8.3%。QTL效應(yīng)分析結(jié)果表明, 相較于單個(gè)位點(diǎn), 多個(gè)抗性QTL的聚合可顯著降低赤霉病嚴(yán)重度。揚(yáng)麥16抗赤霉病QTL將為揭示揚(yáng)麥品種抗性遺傳機(jī)制及開(kāi)發(fā)相應(yīng)分子標(biāo)記奠定基礎(chǔ)。

    小麥; 赤霉病; QTL; 標(biāo)記輔助育種

    小麥赤霉病(Wheat scab, Fusarium head blight, FHB)是由禾谷鐮刀菌()等引起的一種世界性真菌病害[1]。我國(guó)長(zhǎng)江中下游麥區(qū)是赤霉病的常發(fā)區(qū)和重發(fā)區(qū), 2012—2015年, 江蘇省年均赤霉病發(fā)生面積約120萬(wàn)公頃[2]。近年來(lái), 隨著氣候變暖和玉米秸稈還田, 小麥赤霉病發(fā)生愈來(lái)愈重, 正迅速擴(kuò)展到黃淮麥區(qū)。據(jù)統(tǒng)計(jì), 近10年, 河南省赤霉病年均發(fā)生面積達(dá)110萬(wàn)公頃左右, 其中2012年達(dá)333萬(wàn)公頃, 2016年達(dá)117萬(wàn)公頃[3]。

    赤霉病的侵染時(shí)期主要是小麥開(kāi)花期, 赤霉菌侵染小花后迅速在穗部擴(kuò)展, 在小麥籽粒灌漿成熟過(guò)程中不斷繁殖, 進(jìn)而積累各種毒素, 例如脫氧雪腐鐮孢菌烯醇(deoxynivalenol, DON)、雪腐鐮孢菌烯醇(nivalenol, NIV)和玉米赤霉烯酮(zearalenol, ZEN), 嚴(yán)重影響小麥產(chǎn)量, 并對(duì)人畜健康造成巨大傷害, 成為糧食安全的主要威脅[4]。選育和種植抗病品種是應(yīng)對(duì)赤霉病危害最經(jīng)濟(jì)、安全的途徑。小麥赤霉病的抗性主要可以分為抗侵染(type I)、抗擴(kuò)展(type II)、抗DON積累(type III)和籽??剐?type IV) 4種類型, 其中type II類型研究最深入[5]。多年來(lái), 已報(bào)道超過(guò)250個(gè)抗赤霉病QTL, 覆蓋小麥全部21條染色體[6], 但目前明確的小麥抗赤霉病基因只有7個(gè), 其中和是抗擴(kuò)展類型,和是抗侵染類型, 僅已經(jīng)被克隆[7-8]。蘇麥3號(hào)是國(guó)際上研究和利用最廣泛的赤霉病抗源, 望水白是已知的攜帶抗赤霉病位點(diǎn)最多的小麥品種之一[3]。它們均攜帶主效抗病基因, 但兩者綜合農(nóng)藝性狀較差, 稈高易倒伏、穗數(shù)和穗粒數(shù)少、粒重低, 雖利用它們選育出一批抗性好的材料, 例如寧7840等, 但因其綜合豐產(chǎn)性差未得到大面積應(yīng)用[9]。據(jù)報(bào)道, 聚合和其他抗赤霉病位點(diǎn), 可累加抗性效應(yīng), 提高抗性水平[10-11]。揚(yáng)麥系列品種一直是長(zhǎng)江中下游麥區(qū)主體品種, 多數(shù)品種表現(xiàn)中抗赤霉病。揚(yáng)麥16是近年來(lái)長(zhǎng)江中下游推廣面積最大的小麥品種, 利用其為親本選育出優(yōu)質(zhì)豐產(chǎn)抗赤霉病的小麥品種如揚(yáng)麥23和揚(yáng)麥28等, 表明揚(yáng)麥16在小麥遺傳改良中可以作為穩(wěn)定的抗赤親本使用。經(jīng)中國(guó)農(nóng)業(yè)科學(xué)院開(kāi)發(fā)的小麥功能基因標(biāo)記檢測(cè), 揚(yáng)麥16不含[12-13], 可能存在其他主效抗赤霉病基因。因此, 挖掘揚(yáng)麥16的抗赤霉病基因, 拓寬小麥赤霉病抗源, 對(duì)小麥抗赤霉病育種具有重要意義。本研究利用來(lái)自揚(yáng)麥16/中麥895 的174份雙單倍體家系, 經(jīng)660K SNP芯片檢測(cè), 構(gòu)建高密度遺傳圖譜, 結(jié)合多年表型鑒定結(jié)果, 挖掘揚(yáng)麥16抗赤霉病QTL, 為開(kāi)展標(biāo)記輔助抗赤霉病育種提供技術(shù)支撐。

    1 材料與方法

    1.1 試驗(yàn)材料及田間試驗(yàn)

    以揚(yáng)麥16 (Yangmai 16, 簡(jiǎn)稱YM16, 系譜: 揚(yáng)91F138 (揚(yáng)麥158選系)/揚(yáng)90-30)為母本, 中麥895 (Zhongmai 895, 簡(jiǎn)稱ZM895, 系譜: 周麥16/荔墾4號(hào))為父本雜交, 利用玉米花粉誘導(dǎo)雙單倍體培育174份DH家系。以蘇麥3號(hào)(Sumai 3, 簡(jiǎn)稱SM3)和周麥18 (Zhoumai 18, 簡(jiǎn)稱ZM18)為抗赤霉病和感赤霉病對(duì)照。2016—2017年度(簡(jiǎn)稱為2017年)于國(guó)際玉米小麥改良中心(International Maize and Wheat Improvement Center, CIMMYT)溫室種植, 花盆直徑20 cm, 每盆播種8粒, 最后定苗至5株, 2次重復(fù), 參照CIMMYT的正常管理方法管理溫室。2017—2018年度和2018—2019年度(簡(jiǎn)稱為2018年和2019年)于江蘇里下河地區(qū)農(nóng)業(yè)科學(xué)研究所基地種植, 采用完全隨機(jī)設(shè)計(jì)試驗(yàn), 2行區(qū), 每行30粒, 2次重復(fù), 行長(zhǎng)1.3 m, 行距0.3 m, 參照當(dāng)?shù)氐恼T耘喾椒ü芾泶筇铩?/p>

    1.2 赤霉病鑒定和表型數(shù)據(jù)處理

    參照 Yu 等[14]方法制備赤霉菌孢子懸浮液(5000孢子 mL–1), 采用單花滴注法接種。于小麥開(kāi)花初期, 用注射器吸取10 μL孢子液注射入麥穗中間小穗的一個(gè)小花中(頂部小穗開(kāi)始自上而下第6個(gè)小穗), 接種每系15個(gè)穗子。溫室單花滴注接種后套袋保濕3 d, 田間單花接種后立即采用人工彌霧保濕(每半小時(shí)田間彌霧噴水5 min, 6:00–18:00)。接種21 d后調(diào)查病小穗數(shù), 計(jì)算病小穗率。病小穗率(percentage of scabbed spikelets, 簡(jiǎn)稱PSS)作為赤霉病嚴(yán)重度的衡量指標(biāo)[9]。病小穗率(%) = (病小穗數(shù)/總小穗數(shù)) ×100%。用 DPS 軟件和IciMapping的ANOVA功能進(jìn)行群體表型數(shù)據(jù)的相關(guān)分析、方差分析及遺傳力的估算。

    1.3 分子標(biāo)記的檢測(cè)

    取小麥幼苗, 采用CTAB法提取基因組DNA[15], 利用Affymetrix SNP技術(shù)檢測(cè)平臺(tái)(北京博奧生物有限公司)小麥660K SNP芯片分析親本和群體, 利用Genomestudio v1.0軟件進(jìn)行多態(tài)性分析和整合(Illumina, http://www.illumina.com/)。將13個(gè)包括和[12]等重要農(nóng)藝性狀功能基因的相應(yīng)KASP標(biāo)記也整合入多態(tài)性標(biāo)記中。

    1.4 遺傳圖譜構(gòu)建和QTL定位

    篩選出缺失率<10%, 最小等位基因頻率>30%的多態(tài)性標(biāo)記, 使用IciMapping v4.1 (http://www. isbreeding.net/)的BIN功能進(jìn)行去冗余。根據(jù)SNP側(cè)翼序列比對(duì)到中國(guó)春參考基因組的物理位置(IWGSC RefSeq v1.0, https://wheat-urgi.versailles. inra.fr/)和SNP在小麥660K整合圖譜(趙光耀, 私人通訊) KJ-RIL遺傳圖譜中的位置[16], 挑選461個(gè)分布在小麥21條染色體的SNP作為Anchor標(biāo)記, 利用IciMapping v4.1的MAP功能將標(biāo)記分組, 利用MSTmap(http://mstmap.org/)的Kosambi功能對(duì)各組標(biāo)記排序和計(jì)算遺傳距離。利用JoinMap v4.0[17]校正遺傳長(zhǎng)度, 得到最終遺傳圖譜。利用MapChart 2.3 (https://www.wur.nl/en/show/Mapchart.htm)繪制遺傳圖譜。利用IciMapping v4.1的完備區(qū)間作圖法(Inclusive composite interval mapping, ICIM)檢測(cè)抗赤霉病QTL, LOD閾值設(shè)為2.5[18]。將同一染色體上檢測(cè)到的峰值所在遺傳位置之間距離小于10 cM的QTL視為同一個(gè)位點(diǎn)。QTL的命名方法為“Q”加性狀縮寫, 加單位名稱縮寫, 加QTL所在的染色體, 同一染色體上多個(gè)QTL用.1、.2、.3、……來(lái)表示。例如代表 5A 染色體上檢測(cè)出來(lái)的控制赤霉病抗擴(kuò)展的第一個(gè)QTL。將在2個(gè)或者超過(guò)2個(gè)年份下能夠檢測(cè)到的QTL定義為“穩(wěn)定QTL”, 將效應(yīng)值大于10%的QTL定義為“具有較大表型貢獻(xiàn)率的QTL”。為了與前人結(jié)果比較, 將連鎖標(biāo)記或者基因的序列與中國(guó)春參考基因組序列的EnsemblPlants數(shù)據(jù)庫(kù)(http://plants.ensembl.org/)比對(duì), 獲得標(biāo)記或者基因的物理位置[19]。

    2 結(jié)果與分析

    2.1 親本及DH群體的抗性表現(xiàn)

    由表1和圖1可以看出, 3年試驗(yàn)的PSS范圍, 揚(yáng)麥16是14.2%~23.3%, 中麥895是58.7%~69.7%, 抗病對(duì)照蘇麥3號(hào)是5.2%~6.1%, 感病對(duì)照周麥18是60.2%~70.1%。群體的偏度和峰度絕對(duì)值均<1, 分布呈連續(xù)的正態(tài)分布, 可用于QTL定位研究[20]。

    圖1 揚(yáng)麥16/中麥895 DH群體赤霉病嚴(yán)重度頻次分布(2017–2019)

    橫坐標(biāo)代表病小穗率, 縱坐標(biāo)代表相關(guān)病小穗率范圍內(nèi)的家系數(shù)目。

    The abscissa indicates the percentage of scabbed spikelets, and the ordinate indicates the line number.

    3年試驗(yàn)結(jié)果均表明雙親赤霉病的病小穗率存在極顯著的差異, DH群體赤霉病嚴(yán)重度的最小值和最大值之間差異明顯, 存在超親分離, 可以初步判斷雙親都攜帶抗性基因且位點(diǎn)互補(bǔ), 通過(guò)雜交可以得到超親分離的株系。抗病和感病對(duì)照的病小穗率亦存在極顯著差異。

    表1 揚(yáng)麥16/中麥895 DH群體、親本和對(duì)照的赤霉病嚴(yán)重度

    數(shù)據(jù)后不同字母表示不同基因型間赤霉病嚴(yán)重度差異顯著(< 0.01)。

    Values followed by different letters are significant by different in FHB severity among genotypes (< 0.01).

    2.2 遺傳連鎖圖譜

    利用660K SNP芯片, 過(guò)濾和篩選得到152,310個(gè)高質(zhì)量多態(tài)性SNP標(biāo)記, 經(jīng)過(guò)去冗余分析后, 最終“揚(yáng)麥16/中麥895”高密度遺傳圖譜的上圖標(biāo)記為14,480個(gè)。圖譜覆蓋小麥21條染色體, 長(zhǎng)度為3681.7 cM, 密度為3.9個(gè)標(biāo)記 cM–1。分布于小麥A、B和D染色體組的標(biāo)記數(shù)分別為5036、7634和1810個(gè), 連鎖長(zhǎng)度分別為1223.2、1035.6和1423.0 cM。其中, 標(biāo)記密度最高的是1B染色體, 平均11.7個(gè)標(biāo)記 cM–1, 標(biāo)記密度最低的是7D染色體, 平均0.8個(gè)標(biāo)記 cM–1。

    2.3 抗赤霉病QTL定位

    共檢測(cè)到6個(gè)抗赤霉病QTL, 分別位于2DL、3BL、4BS、4DS、5BL和6AS上。除外, 其他5個(gè)QTL均來(lái)自揚(yáng)麥16。在3年中均被檢測(cè)到, 抗性貢獻(xiàn)率8.8%~13.7%,可在2年同時(shí)被檢測(cè)到, 抗性貢獻(xiàn)率9.5%~15.0%。這2個(gè)QTL是本研究檢測(cè)到的穩(wěn)定的抗赤霉病位點(diǎn)。和分別在2017和2019年被檢測(cè)到, 抗性貢獻(xiàn)率分別為10.5%和14.7%。和均僅在2018年被檢測(cè)到, 抗性貢獻(xiàn)率分別為8.3%和6.4% (表2和圖2)。

    表2 揚(yáng)麥16/中麥895 DH 群體抗赤霉病QTL定位

    #加性效應(yīng)為負(fù)說(shuō)明赤霉病抗性來(lái)源于揚(yáng)麥16, 加性效應(yīng)為正說(shuō)明赤霉病抗性來(lái)源于中麥895。

    #The resistance alleles at loci with negative additive effects are from YM16 and that with a positive additive effect is from ZM895.

    圖2 揚(yáng)麥16/中麥895 DH群體抗赤霉病QTL定位情況

    連鎖群右邊是標(biāo)記名稱, 左邊是遺傳位置(CM)。紅色、藍(lán)色和綠色分別代表2017、2018和2019年定位到的QTL, QTL右邊對(duì)應(yīng)LOD值。

    Markers’ names are shown to the right of vertical axis, and their genetic positions are shown in cM to the left. Red, blue, and green represent 2017, 2018, and 2019, respectively. LOD values of QTL are shown on the right side.

    2.4 QTL效應(yīng)分析

    將在2個(gè)或者超過(guò)2個(gè)環(huán)境下能夠檢測(cè)到的QTL定義為“穩(wěn)定QTL”, 本研究分別在3年和2年環(huán)境下均檢測(cè)到和, 根據(jù)它們的基因型將174個(gè)家系分為4種類型, 用3年赤霉病嚴(yán)重度平均值進(jìn)行類型間的差異分析。38個(gè)家系同時(shí)含有和, 赤霉病嚴(yán)重度范圍9.9%~57.4% (平均32.3%), 多數(shù)家系PSS小于40%, 僅有8個(gè)家系PSS在40%~60%范圍內(nèi); 28個(gè)家系只含有, 赤霉病嚴(yán)重度范圍18.9%~78.8% (平均44.8%), 18個(gè)家系PSS小于50%, 10個(gè)家系PSS在50%~80%范圍內(nèi); 有62個(gè)家系只含有, 赤霉病嚴(yán)重度范圍18.6%~71.3% (平均43.2%), 41個(gè)家系PSS小于50%, 21個(gè)家系PSS在50%~80%范圍內(nèi); 46個(gè)家系不含上述2個(gè)QTL, 赤霉病嚴(yán)重度范圍26.1%~88.0% (平均54.3%), 多數(shù)家系PSS大于40%, 僅有6個(gè)家系PSS在20%~40%范圍。聚合2個(gè)抗性位點(diǎn)的家系比含有1個(gè)抗性位點(diǎn)和不含有抗性位點(diǎn)的家系赤霉病嚴(yán)重度顯著降低(表3)。

    表3 揚(yáng)麥16/中麥895 DH群體中存在不同穩(wěn)定赤霉病抗擴(kuò)展QTL家系的赤霉病嚴(yán)重度情況

    數(shù)據(jù)后不同字母表示不同基因型間赤霉病嚴(yán)重度差異顯著(< 0.01)。

    Values followed by different letters are significant by different in FHB severity among genotypes (< 0.01).

    根據(jù)本研究中定位到的效應(yīng)值大于10%的抗性QTL、、和將174個(gè)DH家系分為5種類型, 用3年P(guān)SS平均值進(jìn)行類型間的差異分析。10個(gè)家系同時(shí)含有4個(gè)抗性QTL, 發(fā)病較輕, 平均PSS是30.5%; 46個(gè)家系含有3個(gè)抗性QTL, 平均PSS是36.9%; 65個(gè)家系含有2個(gè)抗性QTL, 平均PSS是41.8%; 46個(gè)家系只含有1個(gè)抗性QTL, 平均PSS是53.84%; 7個(gè)家系不含有任何一個(gè)上述4個(gè)抗性QTL, 平均PSS是67.1% (表4和圖3)。

    表4 揚(yáng)麥16/中麥895 DH群體中存在不同赤霉病抗擴(kuò)展QTL家系的赤霉病嚴(yán)重度情況

    PSS (percentage of scabbed spikelets)代表赤霉病嚴(yán)重度。#數(shù)據(jù)后不同字母表示不同基因型間赤霉病嚴(yán)重度差異極顯著(< 0.01)。

    PSS (percentage of scabbed spikelets) indicates the FHB severity.#Values followed by different letters are significant by different in FHB severity among genotypes (< 0.01).

    圖3 DH家系所含QTL數(shù)目與其平均赤霉病嚴(yán)重度的關(guān)系

    PSS (percentage of scabbed spikelets)代表赤霉病嚴(yán)重度。

    PSS (percentage of scabbed spikelets) indicates the FHB severity.

    3 討論

    3.1 QTL定位結(jié)果與已知抗赤霉病QTL/基因比較

    在揚(yáng)麥16/中麥895 DH群體中共檢測(cè)到6個(gè)赤霉病抗擴(kuò)展QTL, 其中和在多年均可被檢測(cè)到, 是來(lái)源于揚(yáng)麥16的穩(wěn)定QTL位點(diǎn)。根據(jù)比對(duì)結(jié)果,與緊密連鎖, Srinivasachary等[21]報(bào)道附近確實(shí)存在赤霉病抗侵染QTL, 來(lái)源于Spark,可顯著降低小麥對(duì)赤霉病抗侵染能力, 但其與赤霉病抗擴(kuò)展性不相關(guān)[22], Buerstmayr等[23]研究表明攜帶位點(diǎn)的小麥品種易感染赤霉病, 表明較高的株高營(yíng)造的穗部局部低濕度環(huán)境可以減少赤霉病孢子侵染小麥。本研究采用一年CIMMYT溫室單花滴注然后套袋保濕, 兩年揚(yáng)州田間單花滴注試驗(yàn)采取彌霧保濕, 株高對(duì)穗部的微環(huán)境存在一定的影響, 從而影響發(fā)病的嚴(yán)重度。將2017—2019年DH群體株高數(shù)據(jù)平均數(shù)(未發(fā)表資料)與赤霉病嚴(yán)重度平均值進(jìn)行相關(guān)性分析, 發(fā)現(xiàn)株高均值與赤霉病嚴(yán)重度均值相關(guān)系數(shù)–0.33 (<0.01), 表明該群體赤霉病抗性與株高存在一定相關(guān)性。群體中攜帶抗病等位基因家系的平均株高顯著高于攜帶感病等位基因家系的平均株高(< 0.01), 表明該位點(diǎn)確實(shí)與株高顯著相關(guān)。該位點(diǎn)為“一因多效”或者“緊密連鎖”, 可通過(guò)后續(xù)增加溫室環(huán)境下單花滴注試驗(yàn)進(jìn)行驗(yàn)證。另?yè)?jù)報(bào)道, 小麥粒重基因啟動(dòng)子區(qū)域的單倍型變異除了與粒重相關(guān)外, 還與株高顯著相關(guān)[24], 本研究定位到了位于6AS染色體上的抗赤霉病QTL, DH群體2017—2019年的株高平均值(未發(fā)表資料)與赤霉病嚴(yán)重度平均值相關(guān)系數(shù)為–0.33 (<0.01), 表明該位點(diǎn)的赤霉病抗性與株高也存在一定的關(guān)系。

    與來(lái)源于揚(yáng)麥13的位置相近[6]。揚(yáng)麥13的系譜是揚(yáng)88-84//Maris Dove/揚(yáng)麥3號(hào), 揚(yáng)麥3號(hào)是阿夫的紅殼選系; 揚(yáng)麥16的系譜是揚(yáng)麥158優(yōu)系(91F138)/揚(yáng)90-30, 揚(yáng)麥158的系譜是揚(yáng)麥4號(hào)//St1472/506, 揚(yáng)麥4號(hào)的系譜是南大2419/勝利麥1-3-2//阿夫選系。由于揚(yáng)麥品種的親本中存在共同的親本阿夫, 且親緣關(guān)系較近, 推測(cè)與可能是相同抗病位點(diǎn), 后續(xù)我們會(huì)繼續(xù)研究上述材料, 以確定的來(lái)源。

    其他QTL僅在1年被檢測(cè)到, 其中和的效應(yīng)值較大, 來(lái)自揚(yáng)麥16??山忉尡硇妥儺?0.5%, 經(jīng)兩側(cè)標(biāo)記序列與中國(guó)春參考基因組序列比對(duì)發(fā)現(xiàn), 其物理位置區(qū)間是512.7~532.0 Mb, 與已報(bào)道的Wuhan-1中檢測(cè)到的2DL上抗赤霉病位點(diǎn)位置相近(513.1 Mb)[25], Wuhan-1與揚(yáng)麥16均來(lái)自長(zhǎng)江中下游麥區(qū), 推測(cè)可能是相同的QTL。能解釋表型變異14.7%, 經(jīng)比對(duì)發(fā)現(xiàn)無(wú)同位置抗赤霉病QTL報(bào)道, 推測(cè)可能為新的抗赤霉病位點(diǎn)。小麥赤霉病抗性屬于數(shù)量性狀, 發(fā)病受環(huán)境影響很大,還需增加多環(huán)境試驗(yàn)驗(yàn)證上述QTL的有效性和穩(wěn)定性。

    3.2 抗性改良與QTL/基因聚合

    DH家系中僅存在單個(gè)抗赤霉病位點(diǎn)時(shí), 多數(shù)家系的PSS在50%以上, 但當(dāng)存在2位點(diǎn)或者2位點(diǎn)以上的QTL聚合時(shí), 平均PSS明顯降低。聚合包含4個(gè)QTL的家系具有最好的赤霉病抗性, 其中2個(gè)家系的PSS在0~25%區(qū)間, 達(dá)到高抗水平, 其余8個(gè)家系的PSS在25%~50%區(qū)間, 達(dá)到中抗水平, 表明多個(gè)抗病基因聚合能表現(xiàn)出較好的抗性。Jia等[10]報(bào)道與其他赤霉病抗性QTL/基因的聚合能有效提高小麥的赤霉病抗擴(kuò)展和抗侵染水平。Arruda等[11]研究表明, 聚合多個(gè)3B上抗性位點(diǎn)能使赤霉病抗性水平提高。因此, 將多個(gè)抗病基因聚合到同一品種(系)中, 對(duì)于提高赤霉病抗性具有重要意義。

    3.3 抗赤霉病基因的利用

    我國(guó)長(zhǎng)江中下游麥區(qū)是赤霉病的常發(fā)區(qū)和重發(fā)區(qū), 國(guó)內(nèi)外廣泛使用的主效抗病基因區(qū)段與不利農(nóng)藝性狀存在連鎖[26], 因此單純利用蘇麥3號(hào)和望水白背景的不足以解決該麥區(qū)赤霉病的危害。關(guān)于這一點(diǎn)程順和等[27]早有闡述, 并提出抗赤霉病育種的第2條路線, 即選用豐產(chǎn)、中感-中抗材料雜交, 后代注重豐產(chǎn), 兼顧抗赤性等選擇, 選育出揚(yáng)麥158、揚(yáng)麥11、揚(yáng)麥16等一批抗赤性與豐產(chǎn)性相結(jié)合的大面積種植品種。揚(yáng)麥系列多數(shù)品種的赤霉病抗性好且遺傳力較強(qiáng)正在被大范圍用作育種親本, 更有價(jià)值的是其抗病基因不是, 需要進(jìn)一步發(fā)掘其中的主效赤霉病抗性QTL并開(kāi)發(fā)相應(yīng)的標(biāo)記, 為揚(yáng)麥抗赤霉病基因的高效利用奠定基礎(chǔ)。由于多數(shù)抗病揚(yáng)麥品種不具, 將導(dǎo)入揚(yáng)麥背景, 提高抗赤霉病等級(jí), 不失為解決長(zhǎng)江中下游麥區(qū)赤霉病危害的有效途徑, 如揚(yáng)16-157 (含), 其在2017—2019年度國(guó)家小麥良種重大聯(lián)合攻關(guān)試驗(yàn)中表現(xiàn)高抗赤霉病, 與蘇麥3號(hào)相當(dāng), 2017—2018年品比試驗(yàn)中產(chǎn)量為6325.5 kg hm–2, 比對(duì)照揚(yáng)麥20增產(chǎn)7.0%, 2018—2019年品比試驗(yàn)中產(chǎn)量為6886.5 kg hm–2, 比對(duì)照增產(chǎn)4.5%。

    4 結(jié)論

    揚(yáng)麥16中抗赤霉病, 具有主效QTL, 其中4個(gè)QTL具有較大的表型貢獻(xiàn)率, 2個(gè)QTL多年被檢測(cè)到, QTL聚合能顯著降低小麥赤霉病嚴(yán)重度。揚(yáng)麥系列小麥品種QTL的發(fā)掘可以為小麥抗赤霉病改良提供優(yōu)異基因資源。

    [1] Buerstmayr H, Ban T, Anderson J A. QTL mapping and marker-assisted selection forhead blight resistance in wheat: a review., 2009, 128: 1–26.

    [2] 吳佳文, 楊榮明, 朱鳳, 田子華. 2015 年江蘇省小麥赤霉病發(fā)生特點(diǎn)與防控對(duì)策探討. 中國(guó)植保導(dǎo)刊, 2016, 36(10): 31–34. Wu J W, Yang R M, Zhu F, Tian Z H. The epidemic characteristics of wheathead blight in Jiangsu province in 2015 and discussion of its control measures., 2016, 36(10): 31–34 (in Chinese).

    [3] 張愛(ài)民, 陽(yáng)文龍, 李欣, 孫家柱. 小麥抗赤霉病研究現(xiàn)狀與展望. 遺傳, 2018, 40: 858–873. Zhang A M, Yang W L, Li X, Sun J Z. Current status and perspective on research againsthead blight in wheat., 2018, 40: 858–873 (in Chinese with English abstract).

    [4] 史建榮, 劉馨, 仇劍波, 祭芳, 徐劍宏, 董飛, 殷憲超, 冉軍艦. 小麥中鐮刀菌毒素脫氧雪腐鐮刀菌烯醇污染現(xiàn)狀與防控研究進(jìn)展. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47: 3641–3654. Shi J R, Liu X, Qiu B, Ji F, Xu J H, Dong F, Yin X C, Ran J J. Deoxynivalenol contamination in wheat and its management., 2014, 47: 3641–3654 (in Chinese with English abstract).

    [5] Ren J D, Wang Z, Du Z Y, Chen M Z, Zhang Y B, Quan W, Wang Y J, Jiang X, Zhang Z J. Detection and validation of a novel major QTL for resistance tohead blight fromin the terminal region of chromosome 7DL., 2019, 132: 241–255.

    [6] Yi X, Cheng J, Jiang Z, Hu W, Bie T, Gao D, Li D, Wu R, Li Y, Chen S, Cheng X, Liu J, Zhang Y, Cheng S. Genetic analysis ofhead blight resistance in CIMMYT bread wheat line C615 using traditional and conditional QTL mapping., 2018, 9: 573. https://doi.org/10.3389/ fpls.2018.00573.

    [7] Su Z Q, Bernardo A, Tian B, Chen H, Wang S, Ma H X, Cai S B, Liu D T, Zhang D D, Li T, Trick H, St Amand P, Yu J M, Zhang Z Y, Bai G H. A deletion mutation inconfersresistance tohead blight in wheat., 2019, https://doi.org/10.1038/s41588-019-0425-8.

    [8] Li G Q, Zhou J Y, Jia H Y, Gao Z X, Fan M, Luo Y J, Zhao P T, Xue S L, Li N, Yuan Y, Ma S W, Kong Z X, Jia L, An X, Jiang G, Liu W X, Cao W J, Zhang R R, Fan J C, Xu X W, Liu Y F, Kong Q Q, Zheng S H, Wang Y, Qin B, Cao S Y, Ding Y X, Shi J X, Yan H S, Wang X, Ran C F, Ma Z Q. Mutation of a histidine-rich calcium-binding-protein gene in wheat confers resistance tohead blight., 2019, https://doi.org/10.1038/s41588-019-0426-7.

    [9] Li T, Luo M, Zhang D, Wu D, Li L, Bai G. Effective marker alleles associated with type 2 resistance tohead blight infection in fields., 2016, 66: 350–357.

    [10] Jia H Y, Zhou J Y, Xue S L, Li G Q, Yan H S, Ran C F, Zhang Y D, Shi J X, Jia L, Wang X, Luo J, Ma Z Q. A journey to understand wheathead blight resistance in the Chinese wheat landrace Wangshuibai., 2018, 48–59.

    [11] Arruda M P, Brown P, Brown-Guedira G, Krill A M, Thurber C, Merrill K R, Foresman B J, Kolb F L. Genome-wide association mapping ofhead blight resistance in wheat using Genotyping-by-sequencing., 2016,9. https://doi: 10.3835/plantgenome2015.04.0028.

    [12] Rasheed A, Wen W E, Gao F M, Zhai S N, Jin H, Liu J D, Guo Q, Zhang Y J, Dreisigacker S, Xia X C, He Z H. Deve-lopment and validation of KASP assays for genes underpinning key economic traits in bread wheat., 2016, 129: 1843–1860.

    [13] 張宏軍, 宿振起, 柏貴華, 張旭, 馬鴻翔, 李騰, 鄧云, 買春艷, 于立強(qiáng), 劉宏偉, 楊麗, 李洪杰, 周陽(yáng). 利用基因功能標(biāo)記選擇提高黃淮冬麥區(qū)小麥品種對(duì)赤霉病的抗性. 作物學(xué)報(bào), 2018, 44: 505–511. Zhang H J, Su Z Q, Bai G H, Zhang X, Ma H X, Li T, Deng Y, Mai C Y, Yu L Q, Liu H W, Yang L, Li H J, Zhou Y. Improvement of resistance of wheat cultivars tohead blight in the Yellow-Huai rivers valley winter wheat zone with functional marker selection ofgene., 2018, 44: 505–511 (in Chinese with English abstract).

    [14] Yu J B, Bai G H, Cai S B, Ban T. Marker-assisted characterization of Asian wheat lines for resistance tohead blight., 2006, 113: 308–320

    [15] Stacey J, Isaac P G. Isolation of DNA from plants., 1994, 28: 9–15.

    [16] Cui F, Zhang N, Fan Xiao, Zhang W, Zhao C, Yang L, Pan R, Chen M, Han J, Zhao X, Ji J, Tong Y, Zhang H, Jia J, Zhao G, Li J. Utilization of a wheat 660K SNP array-derived high-density genetic map for high-resolution mapping of a major QTL for kernel number., 2017, https://doi:10.1038/s41598-017-04028-6.

    [17] Stam P. Construction of integrated genetic linkage maps by means of a new computer package: JoinMap., 1993, 3: 739–744.

    [18] 王建康. 數(shù)量性狀基因的完備區(qū)間作圖方法. 作物學(xué)報(bào), 2009, 35: 239–245. Wang J K. Inclusive composite interval mapping of quantitative trait genes., 2009, 35: 239–245 (in Chinese with English abstract).

    [19] Zhai S N, He Z H, Wen W E, Jin H, Liu J D, Zhang Y, Liu Z Y, Xia X C. Genome-wide linkage mapping of flour color-related traits and polyphenol oxidase activity in common wheat., 2016, 129: 377–394.

    [20] 李慧慧, 張魯燕, 王建康. 數(shù)量性狀基因定位研究中若干常見(jiàn)問(wèn)題的分析與解答. 作物學(xué)報(bào), 2010, 36: 918–931. Li H H, Zhang L Y, Wang J K. Analysis and answers to frequently asked questions in quantitative trait locus mapping., 2010, 36: 918–931 (in Chinese with English abstract).

    [21] Srinivasachary, Gosman N, Steed A, Simmonds J, Leverington- Waite M, Wang Y, Snape J, Nicholson P. Susceptibility tohead blight is associated with thesemi-dwarfing allele in wheat., 2008, 116: 1145–1153.

    [22] Srinivasachary, Gosman N, Steed A, Hollins T, Bayles R, Jennings P, Nicholson P. Semi-dwarfingandloci of wheat differ significantly in their influence on resistance tohead blight., 2008, 118: 695–702.

    [23] Buerstmaye H, Ban T, Anderson J A. QTL mapping and marker-assisted selection forhead blight resistance in wheat: a review., 2009, 128: 1–26.

    [24] Jaiswal V, Gahlaut V, Mathur S, Agarwal P, Khandelwal M K, Khurana J P, Tyagi A K, Balyan H S, Gupta P K. Identification of novel SNP in promoter sequence ofassociated with grain weight and other agronomic traits in wheat (L.)., 2015, 10: e0129400.

    [25] Somers D J, Fedak G, Savard M. Molecular mapping of novel genes controllinghead blight resistance and deoxynivalenol accumulation in spring wheat., 2003, 46: 555–564.

    [26] 朱展望, 徐登安, 程順和, 高春保, 夏先春, 郝元峰, 何中虎. 中國(guó)小麥品種抗赤霉病基因的鑒定與溯源. 作物學(xué)報(bào), 2018, 44: 473–482. Zhu Z W, Xu D A, Cheng S H, Gao C B, Xia X C, Hao Y F, He Z H. Characterization of Fusarium head blight resistance geneand its putative ancestor in Chinese wheat germplasm., 2018, 44: 473–482 (in Chinese with English abstract).

    [27] 程順和, 張勇, 張伯橋, 高德榮, 吳宏亞, 陸成彬, 呂國(guó)鋒, 王朝順. 小麥抗赤霉病育種2條技術(shù)路線的探討. 揚(yáng)州大學(xué)學(xué)報(bào)(農(nóng)業(yè)與生命科學(xué)版), 2003, 24(1): 59–62. Cheng S H, Zhang Y, Zhang B Q, Gao D R, Wu H Y, Lu C B, Lyu G F, Wang C S. Discussion of two ways of breeding scab resistance in wheat.(Agric Life Sci Edn), 2003, 24(1): 59–62 (in Chinese with English abstract).

    Mapping and genetic analysis of QTLs forhead blight resistance to disease spread in Yangmai 16

    HU Wen-Jing1, ZHANG Yong1, LU Cheng-Bin1, WANG Feng-Ju2, LIU Jin-Dong2, JIANG Zheng-Ning1, WANG Jin-Ping2, ZHU Zhan-Wang2, XU Xiao-Ting2, HAO Yuan-Feng2, HE Zhong-Hu2,3, and GAO De-Rong1,*

    1Lixiahe Institute of Agriculture Sciences / Key Laboratory of Wheat Biology and Genetic Improvement for Low & Middle Yangtze Valley, Ministry of Agriculture and Rural Affairs, Yangzhou 225007, Jiangsu, China;2National Wheat Improvement Center / Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China;3CIMMYT-China Office, Chinese Academy of Agricultural Sciences, Beijing 100081, China

    head blight (FHB) resistance of Yangmai wheat cultivars has been paid much attention, but the underlying genetic mechanism is unclear. In recent years, Yangmai 16 is a predominant wheat cultivar durably resistant to FHB in production. A population of 174 double haploid lines (DH) produced by crossing Yangmai 16 (YM16) with the susceptible cultivar Zhongmai 895 (ZM895) was evaluated for FHB response using point inoculation from 2017 to 2019. The DH population was genotyped with wheat 660K SNP array and a high-density genetic map was constructed. Six resistance QTLs were detected, and among them, five were from the resistant parent Yangmai 16 and one from Zhongmai 895.andwere detected at least in two years, explaining 8.8% to 15.0% of the phenotypic variances, respectively.andwere detected only in one year, accounting for 10.5% and 14.7% of the phenotypic variances.andwere detected in one year, too, accounting for 6.4% and 8.3% of the phenotypic variances, respectively. Pyramiding ofmultiple resistant loci with large effects (>10%) is an effective approach to increase FHB resistance. The QTLs identified from Yangmai 16 in the present study will provide a starting point for genetic studies of other Yangmai cultivars, and the QTLs closely linked to markers will be useful for marker-assisted selection in wheat FHB improvement.

    ;head blight; QTL; marker-assisted breeding

    10.3724/SP.J.1006.2020.91048

    本研究由國(guó)家自然科學(xué)基金項(xiàng)目(31901544), 國(guó)家現(xiàn)代農(nóng)業(yè)產(chǎn)業(yè)技術(shù)體系建設(shè)專項(xiàng)(CARS-03-03B, CARS-3-2-11), 國(guó)家重點(diǎn)研發(fā)計(jì)劃項(xiàng)目(2017YFD0100801, 2017YFD0101802)和江蘇省自然科學(xué)基金項(xiàng)目(BK20171279)資助。

    This study was supported by the National Natural Science Foundation of China, China Agriculture Research System (CARS-03-03B, CARS-3-2-11), the National Key Research and Development Program of China, and the Natural Science Foundation of Jiangsu Province (BK20171279).

    高德榮, E-mail: gdr@wheat.org.cn

    E-mail: huren2008@126.com

    2019-07-22;

    2019-09-26;

    2019-10-09.

    URL: http://kns.cnki.net/kcms/detail/11.1809.s.20191008.1730.006.html

    猜你喜歡
    檢測(cè)
    QC 檢測(cè)
    “不等式”檢測(cè)題
    “一元一次不等式”檢測(cè)題
    “一元一次不等式組”檢測(cè)題
    “幾何圖形”檢測(cè)題
    “角”檢測(cè)題
    “有理數(shù)的乘除法”檢測(cè)題
    “有理數(shù)”檢測(cè)題
    “角”檢測(cè)題
    “幾何圖形”檢測(cè)題
    母亲3免费完整高清在线观看| 久久午夜综合久久蜜桃| 91九色精品人成在线观看| 亚洲中文字幕日韩| 午夜免费观看网址| 一边摸一边抽搐一进一小说| 成人三级做爰电影| 91九色精品人成在线观看| 亚洲专区中文字幕在线| 午夜福利免费观看在线| 成人手机av| 一本精品99久久精品77| 精品国产亚洲在线| 国产精品,欧美在线| 久久亚洲真实| 特大巨黑吊av在线直播 | 久久久久久久久久黄片| 久久午夜亚洲精品久久| 可以在线观看毛片的网站| 变态另类成人亚洲欧美熟女| 午夜免费观看网址| 久久天躁狠狠躁夜夜2o2o| 精品欧美一区二区三区在线| 亚洲va日本ⅴa欧美va伊人久久| 午夜免费成人在线视频| 日韩视频一区二区在线观看| 欧美中文日本在线观看视频| 99国产精品一区二区三区| 2021天堂中文幕一二区在线观 | 美国免费a级毛片| 韩国精品一区二区三区| 免费看美女性在线毛片视频| 人妻久久中文字幕网| 免费在线观看视频国产中文字幕亚洲| 国产又色又爽无遮挡免费看| 麻豆av在线久日| svipshipincom国产片| 亚洲最大成人中文| 99re在线观看精品视频| 免费人成视频x8x8入口观看| 欧美国产日韩亚洲一区| 亚洲中文字幕日韩| 变态另类丝袜制服| 国产日本99.免费观看| 免费观看精品视频网站| 身体一侧抽搐| 亚洲中文av在线| 欧美日韩一级在线毛片| 人人妻人人看人人澡| 亚洲国产欧美日韩在线播放| 国产欧美日韩精品亚洲av| 视频在线观看一区二区三区| 精品久久久久久久久久免费视频| 久久久久久人人人人人| 欧美日韩亚洲综合一区二区三区_| 三级毛片av免费| 中文字幕高清在线视频| 国产在线观看jvid| 中国美女看黄片| 日本在线视频免费播放| 黄色成人免费大全| 日韩欧美一区二区三区在线观看| 色播亚洲综合网| 成熟少妇高潮喷水视频| av免费在线观看网站| 男女下面进入的视频免费午夜 | 真人做人爱边吃奶动态| 女人高潮潮喷娇喘18禁视频| 亚洲人成77777在线视频| 午夜精品在线福利| 精品久久久久久久人妻蜜臀av| 亚洲专区中文字幕在线| 一个人免费在线观看的高清视频| 国产精品美女特级片免费视频播放器 | 制服丝袜大香蕉在线| 欧美色欧美亚洲另类二区| 亚洲午夜理论影院| 99re在线观看精品视频| 久久性视频一级片| 波多野结衣巨乳人妻| 亚洲av成人不卡在线观看播放网| 两个人视频免费观看高清| 色av中文字幕| 久热爱精品视频在线9| 久久99热这里只有精品18| 一边摸一边抽搐一进一小说| 亚洲精品av麻豆狂野| 久久久久久九九精品二区国产 | 久久久精品国产亚洲av高清涩受| 美国免费a级毛片| 成人三级黄色视频| 这个男人来自地球电影免费观看| 久久 成人 亚洲| 日韩精品青青久久久久久| 日韩视频一区二区在线观看| www.精华液| 又紧又爽又黄一区二区| 久久久国产精品麻豆| 在线看三级毛片| av电影中文网址| 亚洲第一av免费看| 日韩大码丰满熟妇| 国产1区2区3区精品| 国产精品久久久久久亚洲av鲁大| 色播亚洲综合网| 国产人伦9x9x在线观看| 亚洲全国av大片| 黑丝袜美女国产一区| 国产单亲对白刺激| 黑人巨大精品欧美一区二区mp4| 免费电影在线观看免费观看| 国产成人av教育| 国产伦一二天堂av在线观看| 香蕉av资源在线| 一本一本综合久久| 中文字幕精品免费在线观看视频| 巨乳人妻的诱惑在线观看| 老鸭窝网址在线观看| 久久人妻av系列| 每晚都被弄得嗷嗷叫到高潮| 少妇裸体淫交视频免费看高清 | 夜夜躁狠狠躁天天躁| 日韩免费av在线播放| av在线天堂中文字幕| 99久久综合精品五月天人人| √禁漫天堂资源中文www| 国产精品99久久99久久久不卡| 欧美激情 高清一区二区三区| 真人一进一出gif抽搐免费| 久久国产亚洲av麻豆专区| 国产av在哪里看| 俄罗斯特黄特色一大片| 中文在线观看免费www的网站 | 国产蜜桃级精品一区二区三区| a在线观看视频网站| 国产99久久九九免费精品| 国产精品影院久久| 亚洲精品久久国产高清桃花| 真人做人爱边吃奶动态| 草草在线视频免费看| 成人手机av| 中文字幕高清在线视频| 少妇被粗大的猛进出69影院| 男人操女人黄网站| 99国产精品一区二区三区| 色综合亚洲欧美另类图片| 久久久久九九精品影院| 香蕉久久夜色| 欧美黄色片欧美黄色片| 亚洲人成伊人成综合网2020| 成年版毛片免费区| 亚洲av成人av| 91麻豆av在线| 国产高清视频在线播放一区| 一夜夜www| 国产片内射在线| 国产精品综合久久久久久久免费| 69av精品久久久久久| 久久人妻av系列| 自线自在国产av| 亚洲七黄色美女视频| 天堂影院成人在线观看| av中文乱码字幕在线| 亚洲国产精品成人综合色| 国产成人av激情在线播放| 91字幕亚洲| 少妇裸体淫交视频免费看高清 | 久久 成人 亚洲| 精华霜和精华液先用哪个| www.www免费av| 麻豆一二三区av精品| 一级a爱片免费观看的视频| 日韩大码丰满熟妇| 欧美性长视频在线观看| 亚洲,欧美精品.| 一进一出抽搐动态| 国产亚洲精品久久久久久毛片| 国产精品九九99| 久久久久国产一级毛片高清牌| 又紧又爽又黄一区二区| 露出奶头的视频| 精品日产1卡2卡| 亚洲精品在线美女| 免费高清在线观看日韩| 成人18禁高潮啪啪吃奶动态图| 男人舔女人的私密视频| 久久久久国内视频| 国产伦一二天堂av在线观看| 深夜精品福利| 免费在线观看成人毛片| 亚洲人成网站在线播放欧美日韩| 亚洲九九香蕉| 免费观看精品视频网站| 日本 av在线| 亚洲色图 男人天堂 中文字幕| 精品高清国产在线一区| 亚洲精品久久成人aⅴ小说| 后天国语完整版免费观看| 亚洲av片天天在线观看| 午夜福利一区二区在线看| 在线观看日韩欧美| 少妇 在线观看| 亚洲中文av在线| 欧美中文综合在线视频| 999久久久国产精品视频| 国产成人精品久久二区二区91| 人妻丰满熟妇av一区二区三区| 欧美日本亚洲视频在线播放| 精品熟女少妇八av免费久了| 日本成人三级电影网站| 18禁美女被吸乳视频| 成人18禁高潮啪啪吃奶动态图| 男女午夜视频在线观看| 国产精华一区二区三区| 国产精品国产高清国产av| 19禁男女啪啪无遮挡网站| 色综合欧美亚洲国产小说| 麻豆久久精品国产亚洲av| 性色av乱码一区二区三区2| 91老司机精品| 国产av一区二区精品久久| 99久久国产精品久久久| 国产99白浆流出| 制服丝袜大香蕉在线| 50天的宝宝边吃奶边哭怎么回事| 成人三级黄色视频| 免费在线观看亚洲国产| 国产午夜福利久久久久久| 中文字幕另类日韩欧美亚洲嫩草| 黄色a级毛片大全视频| 精华霜和精华液先用哪个| 两性夫妻黄色片| 搞女人的毛片| 在线免费观看的www视频| 欧美中文日本在线观看视频| 97碰自拍视频| 香蕉国产在线看| 99国产精品一区二区三区| 老汉色av国产亚洲站长工具| 国产aⅴ精品一区二区三区波| 最好的美女福利视频网| 欧美人与性动交α欧美精品济南到| 天天一区二区日本电影三级| 亚洲性夜色夜夜综合| 两个人看的免费小视频| 一区二区日韩欧美中文字幕| 亚洲第一欧美日韩一区二区三区| 97碰自拍视频| 中文字幕人妻熟女乱码| 黄网站色视频无遮挡免费观看| 日本免费a在线| 日韩大尺度精品在线看网址| 欧美激情 高清一区二区三区| 亚洲第一欧美日韩一区二区三区| 国产精品久久久久久亚洲av鲁大| av福利片在线| 亚洲精品中文字幕在线视频| 日韩高清综合在线| 又黄又爽又免费观看的视频| 欧美激情高清一区二区三区| 精品国产国语对白av| 精品乱码久久久久久99久播| svipshipincom国产片| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 一区二区三区激情视频| 夜夜爽天天搞| 国产精品国产高清国产av| 国产亚洲欧美在线一区二区| 制服人妻中文乱码| 波多野结衣av一区二区av| 午夜精品久久久久久毛片777| 女性被躁到高潮视频| 成人亚洲精品一区在线观看| 亚洲精品av麻豆狂野| 久久天堂一区二区三区四区| 欧美国产日韩亚洲一区| 午夜a级毛片| 国产日本99.免费观看| 色哟哟哟哟哟哟| 久久亚洲真实| 成年版毛片免费区| 亚洲av熟女| 中文字幕高清在线视频| 国产一区二区激情短视频| 国产人伦9x9x在线观看| 人人妻人人看人人澡| 一区二区三区高清视频在线| 1024手机看黄色片| 久久久久久免费高清国产稀缺| 99国产精品一区二区蜜桃av| 欧美乱色亚洲激情| 黑丝袜美女国产一区| 亚洲 欧美一区二区三区| 国产欧美日韩一区二区三| 国产单亲对白刺激| 精品久久久久久,| 国产国语露脸激情在线看| 夜夜看夜夜爽夜夜摸| 99久久无色码亚洲精品果冻| 色婷婷久久久亚洲欧美| 狠狠狠狠99中文字幕| 人人妻人人看人人澡| 一二三四社区在线视频社区8| 伦理电影免费视频| 精品国产美女av久久久久小说| 成年女人毛片免费观看观看9| 亚洲专区中文字幕在线| 日本一区二区免费在线视频| 久久久久久人人人人人| 给我免费播放毛片高清在线观看| 日韩国内少妇激情av| 琪琪午夜伦伦电影理论片6080| 91老司机精品| 天堂动漫精品| 中文字幕久久专区| 宅男免费午夜| 日韩欧美国产在线观看| 91麻豆精品激情在线观看国产| 一卡2卡三卡四卡精品乱码亚洲| 91成人精品电影| 国产成人av教育| 12—13女人毛片做爰片一| 久久九九热精品免费| 一级a爱视频在线免费观看| 大型av网站在线播放| 欧美绝顶高潮抽搐喷水| 免费看美女性在线毛片视频| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 亚洲 欧美一区二区三区| 久久久久久久久久黄片| 国产成人av激情在线播放| 国产精品久久久人人做人人爽| 校园春色视频在线观看| 国产又色又爽无遮挡免费看| www.精华液| 男人舔女人下体高潮全视频| 国产精品久久电影中文字幕| 久热这里只有精品99| 亚洲av电影不卡..在线观看| 久久狼人影院| 欧美中文综合在线视频| 丰满的人妻完整版| 久久人妻av系列| 91麻豆精品激情在线观看国产| 一区二区三区高清视频在线| 欧美乱妇无乱码| 亚洲三区欧美一区| 亚洲成av人片免费观看| 丝袜在线中文字幕| 久久久国产成人免费| 国产一区二区三区视频了| 高清毛片免费观看视频网站| 动漫黄色视频在线观看| 亚洲色图 男人天堂 中文字幕| 精品国产超薄肉色丝袜足j| 黄片大片在线免费观看| 一区二区三区激情视频| 1024手机看黄色片| 亚洲全国av大片| 亚洲成人久久爱视频| 国产高清视频在线播放一区| 三级毛片av免费| 自线自在国产av| 亚洲人成电影免费在线| 自线自在国产av| 热99re8久久精品国产| ponron亚洲| 国产极品粉嫩免费观看在线| 男人舔奶头视频| 精品无人区乱码1区二区| 伊人久久大香线蕉亚洲五| 国产成人av教育| 亚洲专区中文字幕在线| 免费在线观看亚洲国产| 在线视频色国产色| 国产激情久久老熟女| 亚洲国产看品久久| 欧美成人免费av一区二区三区| 久久国产乱子伦精品免费另类| 最新在线观看一区二区三区| 国产精品影院久久| www国产在线视频色| 国产又爽黄色视频| 精品久久久久久成人av| 搡老岳熟女国产| 国产视频内射| 国产一区在线观看成人免费| 午夜久久久在线观看| 国产又色又爽无遮挡免费看| 伦理电影免费视频| www.999成人在线观看| 50天的宝宝边吃奶边哭怎么回事| 他把我摸到了高潮在线观看| 99久久综合精品五月天人人| 久久99热这里只有精品18| 中文字幕最新亚洲高清| 天堂√8在线中文| 亚洲人成网站高清观看| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 视频区欧美日本亚洲| 男男h啪啪无遮挡| 成在线人永久免费视频| 熟妇人妻久久中文字幕3abv| 欧美zozozo另类| 白带黄色成豆腐渣| 一本久久中文字幕| 欧美一区二区精品小视频在线| 18禁黄网站禁片午夜丰满| 婷婷精品国产亚洲av在线| 国产97色在线日韩免费| avwww免费| 国内久久婷婷六月综合欲色啪| 国产精品久久久久久精品电影 | 久久久久久久午夜电影| 久久精品亚洲精品国产色婷小说| 免费搜索国产男女视频| 在线观看www视频免费| 精品欧美国产一区二区三| 国产高清视频在线播放一区| 亚洲成a人片在线一区二区| 97人妻精品一区二区三区麻豆 | 亚洲精品一区av在线观看| 成年免费大片在线观看| 国产主播在线观看一区二区| 国产又黄又爽又无遮挡在线| 19禁男女啪啪无遮挡网站| 国产亚洲精品综合一区在线观看 | 老司机午夜福利在线观看视频| 亚洲国产看品久久| 久久久久精品国产欧美久久久| 黄片小视频在线播放| 精品国内亚洲2022精品成人| 亚洲久久久国产精品| 欧美最黄视频在线播放免费| 俄罗斯特黄特色一大片| 亚洲色图av天堂| 这个男人来自地球电影免费观看| 一区二区三区高清视频在线| 十分钟在线观看高清视频www| 美女午夜性视频免费| 亚洲成国产人片在线观看| 一a级毛片在线观看| 伦理电影免费视频| 国产精品九九99| 香蕉久久夜色| 日本 欧美在线| 妹子高潮喷水视频| xxxwww97欧美| 久久精品aⅴ一区二区三区四区| 午夜久久久久精精品| 啦啦啦 在线观看视频| 中文资源天堂在线| 色精品久久人妻99蜜桃| 午夜视频精品福利| 久久精品国产99精品国产亚洲性色| 久久精品aⅴ一区二区三区四区| 日本一本二区三区精品| 亚洲av日韩精品久久久久久密| 熟妇人妻久久中文字幕3abv| 免费看十八禁软件| 欧美一区二区精品小视频在线| 成年免费大片在线观看| 丰满的人妻完整版| 男人舔女人的私密视频| 久久草成人影院| 午夜免费观看网址| 热re99久久国产66热| 久久精品国产亚洲av高清一级| 黑人巨大精品欧美一区二区mp4| 麻豆成人午夜福利视频| 波多野结衣av一区二区av| 亚洲第一电影网av| 国产午夜精品久久久久久| 91九色精品人成在线观看| 一个人观看的视频www高清免费观看 | 成年人黄色毛片网站| 日韩欧美国产一区二区入口| 精品久久久久久成人av| 天堂√8在线中文| 美国免费a级毛片| av超薄肉色丝袜交足视频| 久久婷婷成人综合色麻豆| xxxwww97欧美| 国产欧美日韩一区二区三| 丁香六月欧美| 满18在线观看网站| 久久久久久免费高清国产稀缺| 亚洲人成伊人成综合网2020| 久热这里只有精品99| 精品国内亚洲2022精品成人| 搡老妇女老女人老熟妇| 亚洲中文字幕日韩| 99久久精品国产亚洲精品| 免费看a级黄色片| 国产精品 欧美亚洲| 日韩三级视频一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 一本久久中文字幕| 神马国产精品三级电影在线观看 | 欧美精品亚洲一区二区| 国产真实乱freesex| 欧美av亚洲av综合av国产av| 观看免费一级毛片| 丝袜美腿诱惑在线| 一边摸一边做爽爽视频免费| 国产日本99.免费观看| 国产免费av片在线观看野外av| 黄色毛片三级朝国网站| 99久久无色码亚洲精品果冻| 日韩免费av在线播放| 精品日产1卡2卡| 十八禁网站免费在线| 老司机深夜福利视频在线观看| 国产成人影院久久av| 亚洲av成人不卡在线观看播放网| svipshipincom国产片| 国产亚洲欧美在线一区二区| 午夜激情av网站| 亚洲成人久久爱视频| 在线观看一区二区三区| 国产精品电影一区二区三区| 国产伦人伦偷精品视频| 国产午夜福利久久久久久| 国产精品亚洲av一区麻豆| 在线免费观看的www视频| 久久久久免费精品人妻一区二区 | 757午夜福利合集在线观看| tocl精华| 欧美在线黄色| 免费在线观看日本一区| 久久久久久久精品吃奶| 少妇粗大呻吟视频| 欧美另类亚洲清纯唯美| 国产午夜精品久久久久久| 欧美不卡视频在线免费观看 | 一区二区三区国产精品乱码| 中亚洲国语对白在线视频| www.自偷自拍.com| 久久久久精品国产欧美久久久| 一个人免费在线观看的高清视频| 亚洲五月色婷婷综合| 老司机午夜十八禁免费视频| 久久天堂一区二区三区四区| 女人爽到高潮嗷嗷叫在线视频| 十八禁人妻一区二区| 91成人精品电影| 黄片播放在线免费| 51午夜福利影视在线观看| 免费在线观看影片大全网站| 亚洲精品国产一区二区精华液| 亚洲中文av在线| 变态另类成人亚洲欧美熟女| 国产视频内射| 亚洲成av片中文字幕在线观看| 无限看片的www在线观看| 大香蕉久久成人网| 久久精品人妻少妇| 欧美黄色淫秽网站| 国产又色又爽无遮挡免费看| 国产精品野战在线观看| 久久天躁狠狠躁夜夜2o2o| 999久久久精品免费观看国产| 黄色片一级片一级黄色片| 国产精品久久久av美女十八| 脱女人内裤的视频| 亚洲国产欧美一区二区综合| 日本五十路高清| 亚洲人成网站高清观看| 满18在线观看网站| 久久精品人妻少妇| 中文在线观看免费www的网站 | 可以免费在线观看a视频的电影网站| 亚洲精品美女久久久久99蜜臀| 激情在线观看视频在线高清| 午夜老司机福利片| 国产精品免费一区二区三区在线| 国产成+人综合+亚洲专区| 国产亚洲av嫩草精品影院| www日本黄色视频网| 黄色毛片三级朝国网站| 久久精品成人免费网站| 欧美av亚洲av综合av国产av| 可以在线观看毛片的网站| 麻豆av在线久日| 久久国产乱子伦精品免费另类| 国内精品久久久久精免费| 母亲3免费完整高清在线观看| 日本五十路高清| 怎么达到女性高潮| 精品久久久久久久末码| 国产野战对白在线观看| 夜夜躁狠狠躁天天躁| 男女视频在线观看网站免费 | 免费观看人在逋| 亚洲成人久久性| 国内久久婷婷六月综合欲色啪| 国产精品av久久久久免费| 免费在线观看影片大全网站| 欧美国产日韩亚洲一区| 免费观看人在逋| 精品电影一区二区在线| 制服人妻中文乱码| 国产成人精品无人区| 韩国精品一区二区三区| 黄片播放在线免费| 欧美黑人巨大hd| 国产野战对白在线观看| 亚洲 欧美一区二区三区| 这个男人来自地球电影免费观看| 97人妻精品一区二区三区麻豆 | 亚洲熟妇中文字幕五十中出| 国产精品1区2区在线观看.|