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

    施氮與不施氮條件下玉米開(kāi)花期相關(guān)性狀的 QTL 定位

    2017-04-14 05:28:11郭向陽(yáng)陳建軍衛(wèi)曉軼祝云芳王安貴劉鵬飛湯繼華陳澤輝
    關(guān)鍵詞:吐絲開(kāi)花期氮素

    郭向陽(yáng),陳建軍,衛(wèi)曉軼,祝云芳,王安貴,劉鵬飛,湯繼華,陳澤輝*

    (1 貴州省旱糧研究所,貴陽(yáng) 550006;2 貴州省遵義市輝煌種業(yè)有限責(zé)任公司,遵義 563100;3 河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,鄭州 450002;4 河南省新鄉(xiāng)市農(nóng)業(yè)科學(xué)院,新鄉(xiāng) 453000)

    施氮與不施氮條件下玉米開(kāi)花期相關(guān)性狀的 QTL 定位

    郭向陽(yáng)1,陳建軍2,衛(wèi)曉軼3,4,祝云芳1,王安貴1,劉鵬飛1,湯繼華3*,陳澤輝1*

    (1 貴州省旱糧研究所,貴陽(yáng) 550006;2 貴州省遵義市輝煌種業(yè)有限責(zé)任公司,遵義 563100;3 河南農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,鄭州 450002;4 河南省新鄉(xiāng)市農(nóng)業(yè)科學(xué)院,新鄉(xiāng) 453000)

    【目的】探究施氮與不施氮對(duì)玉米開(kāi)花期相關(guān)性狀變異的影響,并定位相關(guān) QTLs。 【方法】以玉米的骨干自交系綜 3 為供體親本,許 178 為受體親本,通過(guò)雜交、回交和分子標(biāo)記輔助選擇的方法,構(gòu)建了一套以許178 為遺傳背景的綜 3 染色體單片段代換系 (SSSLs) 群體,其中包含 160 個(gè)單片段代換系。以這套 SSSLs 以及許 178 為材料,采用裂區(qū)試驗(yàn)設(shè)計(jì),在施氮 (N+) 和不施氮 (N-) 條件下,通過(guò)一年三點(diǎn) (貴州貴陽(yáng)、德江和云南羅平) 的表型評(píng)價(jià),利用復(fù)合區(qū)間作圖法對(duì)玉米散粉期、吐絲期和散粉-吐絲間隔期 (ASI) 3 個(gè)開(kāi)花期相關(guān)性狀進(jìn)行 QTL 定位。 【結(jié)果】在基因組范圍內(nèi),2 種氮處理?xiàng)l件下共定位到 54 個(gè)開(kāi)花期相關(guān)性狀 QTLs,主要定位在第 1、3、6、9 和 10 染色體上。其中 5 個(gè) QTLs 在 3 個(gè)環(huán)境中均被檢測(cè)到。施氮條件下,吐絲期 (DTS) 相關(guān)位點(diǎn) qDTS9a,位于第 9 染色體,可解釋表型變異的 3.05%;吐絲散粉間隔期 (ASI) 相關(guān)位點(diǎn) qASI10a,位于第10 染色體,可解釋表型變異的 30.74%。不施氮條件下,散粉期 (DTP) 相關(guān)位點(diǎn) qDTP9,可解釋表型變異的3.43%;吐絲期相關(guān)位點(diǎn) qDTS9a,位于第 9 染色體,可解釋表型變異的 4.08%;ASI 相關(guān)位點(diǎn) qASI10,位于第10 染色體,可解釋表型變異的 50.28%。對(duì)不同處理?xiàng)l件下的定位結(jié)果比較發(fā)現(xiàn),不施氮條件下檢測(cè)到的 QTLs數(shù)目顯著高于施氮條件下的檢測(cè)數(shù)量。 【結(jié)論】不同氮素處理下存在一些共有的控制玉米開(kāi)花期相關(guān)的遺傳區(qū)段,分別位于 Bin9.02 (umc1170-umc1636-bnlg1401-umc1271) 和 Bin10.04 (umc1077-umc1053-umc2350),這些區(qū)段可能在玉米氮素吸收、轉(zhuǎn)運(yùn)和利用過(guò)程中起重要作用,可作為下一步精細(xì)定位和圖位克隆玉米開(kāi)花期相關(guān)基因的重要候選區(qū)域。

    玉米;SSSL;不同氮處理;開(kāi)花期相關(guān)性狀;數(shù)量性狀基因座 (QTLs)

    氮作為玉米必需的大量營(yíng)養(yǎng)元素之一,在其生長(zhǎng)和發(fā)育過(guò)程中起著重要作用。長(zhǎng)期以來(lái),為追求玉米產(chǎn)量,氮肥被大量使用在玉米的栽培過(guò)程中。有報(bào)道顯示,發(fā)達(dá)國(guó)家的農(nóng)業(yè)氮肥使用量超過(guò) 500萬(wàn)噸,發(fā)展中國(guó)家則遠(yuǎn)遠(yuǎn)超過(guò)該數(shù)字。過(guò)量使用氮肥 造 成 土 壤 、 淡 水 以 及 大 氣 等 嚴(yán) 重 污 染[1]。 長(zhǎng) 期 以來(lái),為探索玉米產(chǎn)業(yè)可持續(xù)發(fā)展途徑,研究者從土壤供氮能力、施肥量、作物本身的氮素吸收和利用能力等方面對(duì)氮素利用效率進(jìn)行了大量的探索,揭示了不同地區(qū)的施肥量與玉米高產(chǎn)潛力的關(guān)系,說(shuō)明了不同土壤的供氮能力對(duì)玉米產(chǎn)量提高的作用,并在此基礎(chǔ)上培育出一系列的高產(chǎn)、高效優(yōu)勢(shì)玉米自交系或雜交種。如 Mo17,掖 478 等在低氮條件下往往表現(xiàn)出嚴(yán)重的缺氮癥狀,而大多數(shù)自交系表現(xiàn)為中間類(lèi)型。沈 125 表現(xiàn)出一定程度的低氮耐受性[2]。劉建安等[3]研究表明,不同雜交種對(duì)氮水平的反應(yīng)存在顯著差異。劉宗華等[4]利用許 178 和黃 C 構(gòu)建的F2∶3群體研究了不同氮水平條件下玉米株高變異的遺傳差異,結(jié)果表明許 178 具有較黃 C 更高的低氮耐受性,并定位到了一些與氮素利用效率相關(guān)的重要QTLs。鄭祖平等[5]利用 Mo17 和黃早四構(gòu)建的 RIL 群體 定 位 到42 個(gè) 可 能 與 玉 米 對(duì) 低 氮 脅 迫 相 關(guān) 的QTLs。然而,氮素吸收利用是一個(gè)復(fù)雜的網(wǎng)絡(luò)調(diào)控過(guò)程,涉及到很多微效位點(diǎn)或基因,且不同性狀在不同氮處理?xiàng)l件下,其遺傳差異也不盡相同。因此,針對(duì)不同性狀,構(gòu)建相應(yīng)的分離群體并以不同的表型性狀作為目標(biāo),探究不同氮處理水平對(duì)目標(biāo)表型變異的影響并定位出相應(yīng)的遺傳位點(diǎn),將會(huì)為深入解析玉米氮素利用效率的遺傳基礎(chǔ)提供進(jìn)一步的支持。

    開(kāi)花期作為玉米重要農(nóng)藝性狀,對(duì)玉米籽粒產(chǎn)量、地域與季節(jié)的適應(yīng)性都具有十分重要的決定作用[6]。不同來(lái)源種質(zhì)之間的開(kāi)花期相關(guān)性狀存在豐富的變異,在溫帶材料與熱帶種質(zhì)之間的差異尤其突出。為研究不同種質(zhì)開(kāi)花期變異的遺傳基礎(chǔ),研究者基于不同群體,定位了大量與開(kāi)花期相關(guān)遺傳位點(diǎn)。目前,在 MaizeGDB 上已經(jīng)公布的玉米花期相關(guān) QTLs 達(dá) 201 個(gè)(http://www.maizegdb.org/qtl.php),其中 156 個(gè) QTLs 與散粉期、吐絲期和開(kāi)花-吐絲間隔期 (ASI) 直接相關(guān)[7]。前人利用不同群體對(duì)玉米開(kāi)花期相關(guān)性狀進(jìn)行了定位研究;Veldboom 等[8]利用Mol7 × H99 的 F2∶3家系,將影響散粉期的 QTL 定位到 1、5、6、7、8 和 9 染色體上;影響吐絲期的 5個(gè) QTLs 定位在 1、2、6、7 和 8 染色體上,影響ASI 的 QTLs 定位在第 6 和 8 染色體上。Sari-Gorla等[9]利用 B73 × H99 的 RILs 群體和 Ribaut 等[10]利用Ac7643 × Ac7729/TZSRW 的 F2∶3家系,定位到的與散粉期、吐絲期、ASI 相關(guān)的 QTLs 與前者有所差異。此外,Agrama 等[11]、Khairallah 等[12]、Austin 等[13]也對(duì)玉米花期相關(guān)性狀進(jìn)行了 QTL 定位。蘭進(jìn)好[6]研究表明,花期性狀 QTL 的基因作用方式以加性和部分顯性為主,上位性效應(yīng)對(duì)玉米花期性狀具有十分 重 要 的 作 用 , 韓 婭 楠 等[14]研 究 也 得 出 同 樣 結(jié) 果 。Li 等[15]利用黃早四 × 掖 107 的 F2∶3群體,在灌溉和干旱脅迫下,分別檢測(cè)到 3 個(gè)和 2 個(gè)控制 ASI 的QTLs,分別位于第 1、2、3 和第 2、5 染色體;同時(shí),發(fā)現(xiàn) ASI 與產(chǎn)量性狀 QTL 存在顯著關(guān)系。胡彥民等[16]利用瓊 68 和 K12 的 F2∶3群體定位到 22 個(gè)開(kāi)花期相關(guān)的 QTLs,散粉期 5 個(gè)、吐絲期 4 個(gè)、散粉至吐絲間隔期 5 個(gè)、散粉周期 4 個(gè)、吐絲周期 2個(gè)、花絲活力 2 個(gè),結(jié)果發(fā)現(xiàn)多數(shù) QTL 存在環(huán)境表達(dá)的特異性。Gonzalo 等[17]利用 B73 和 Mo17 構(gòu)建的RILs 群體,對(duì)不同密度下的開(kāi)花期相關(guān)性狀進(jìn)行定位,分別在第 2、8 染色體上檢測(cè)到開(kāi)花期相關(guān)QTLs。但是利用開(kāi)花期相關(guān)性狀作為目標(biāo),探究不同氮素處理對(duì)其表型變異的影響并定位相關(guān) QTLs 的研究則少有報(bào)道。以單片段代換系為材料,采用施氮和不施氮兩種處理,鑒定不同施氮水平下的開(kāi)花期相關(guān)性狀變異情況,并利用符合區(qū)間作圖法定位控制花期相關(guān)性狀的 QTLs,借助不同處理?xiàng)l件下的定位結(jié)果比較,鑒定出兩種條件下的共有位點(diǎn)。通過(guò)對(duì)兩種施氮條件下的玉米開(kāi)花期相關(guān)性狀進(jìn)行QTL 定位,發(fā)掘控制玉米開(kāi)花期相關(guān)的遺傳區(qū)段,借助不同條件下的定位結(jié)果一致性分析,鑒定出與氮素吸收利用顯著相關(guān)的候選遺傳區(qū)段,為揭示玉米氮素效率利用差異以及分子輔助的提高玉米氮利用效率育種提供分子遺傳學(xué)依據(jù)。

    1 材料與方法

    1.1 供試材料

    本研究所用的 160 份玉米單片段代換系由河南農(nóng)業(yè)大學(xué)湯繼華教授提供。該群體是以優(yōu)良玉米自交系綜 3 為供體親本、許 178 為受體親本,通過(guò)雜交、回交、自交,結(jié)合 SSR 分子標(biāo)記輔助選擇的方法構(gòu)建而成。綜 3 和許 178 作為中國(guó)玉米育種中的骨干材料,在生產(chǎn)上發(fā)揮了重要作用,以這兩個(gè)材料為親本之一的重要玉米雜交種如豫玉 22 號(hào) (綜 3 × 87 - 1) 和農(nóng)大 108 (黃 C × 許 178) 等在生產(chǎn)上被大面積推廣應(yīng)用。

    1.2 試驗(yàn)區(qū)概況與試驗(yàn)設(shè)計(jì)

    2015 年分別在貴州省德江縣 (北緯 28°15′51″,東經(jīng) 108°06′39″)、貴陽(yáng)市 (北緯 26°34′13″,東經(jīng)106°42′39″) 和云南省羅平縣 (北緯 24°31′,東經(jīng)103°57′) 進(jìn)行田間試驗(yàn)。試驗(yàn)采用裂區(qū)設(shè)計(jì),不同氮處理水平作為主區(qū),設(shè)施氮 (N+) 和不施氮 (N-) 兩個(gè)水平,三個(gè)親本群體 (SSSLs) 作為副區(qū)。單行區(qū),小區(qū)行長(zhǎng) 4.0 m,寬行 0.7 m,每行 17 株,株距0.235 m,密度為 52500 plant/hm2,每個(gè)處理 2 次重復(fù)。對(duì)代表花期相關(guān)性狀的散粉期、吐絲期和散粉-吐絲間隔期進(jìn)行鑒定和檢測(cè)。散粉期是指播種至全區(qū) 50% 以上植株雄花開(kāi)始散粉的天數(shù),吐絲期是指播種至全區(qū) 50% 以上植株雌穗花絲吐出的天數(shù),散粉-吐絲間隔是指從散粉期到吐絲期之間的天數(shù)。試驗(yàn)點(diǎn)各小區(qū)供試土壤耕層養(yǎng)分含量見(jiàn)表 1。

    表1 0—20 cm 土層基礎(chǔ)肥力狀況Table 1 Soil fertility at the 0-20 cm layer

    施氮處理在苗期和大喇叭口期各施尿素375 kg/hm2(N+) 。 每 個(gè) 小 區(qū) 均 基 施 過(guò) 磷 酸 鈣375 kg/hm2、硫酸鉀 187.5 kg/hm2,田間管理措施同大田。

    1.3 表型數(shù)據(jù)分析與 QTL 定位

    采用 SPSS17.0 統(tǒng)計(jì)軟件,對(duì) 160 份 SSSLs 的開(kāi)花期進(jìn)行單因素的方差分析和 Duncan’s 多重比較。本研究中涉及到的基因型鑒定數(shù)據(jù)的相關(guān)研究結(jié)果[18]已發(fā)表。利用完備區(qū)間作圖法對(duì)生育期相關(guān) QTL 進(jìn)行定位,若某一性狀與許 178 在 P < 0.05 水平上差異顯著,則認(rèn)為該單片段代換系代換片段上存在一個(gè) QTL,并同時(shí)估算各個(gè) QTL 的加性效應(yīng)值及貢獻(xiàn)率[19]。

    加性效應(yīng)值 (A) = (單片段代換系的表型值 - 對(duì)照的表型值)/2;

    貢獻(xiàn)率 (R) = 加性效應(yīng)值/對(duì)照的表型值 × 100%。

    2 結(jié)果與分析

    2.1 開(kāi)花期相關(guān)性狀

    不同 SSSLs 開(kāi)花期相關(guān)性狀的表型調(diào)查結(jié)果列于表 2。施氮條件下,散粉期不同親本群體在德江的變異范圍為 84.5~94.5 d,平均值為 88.3 ± 1.9 d;在貴 陽(yáng) 和 羅 平 的 變 異 范 圍 分 別 為74.0~86.0 d 、68.0~80.5 d。吐絲期不同 SSSLs 在德江、貴陽(yáng)和羅平的變異范圍分別為 90.0~99.5 d、76.5~87.0 d 和71.0~84.0 d。散粉-吐絲間隔期 (ASI) 不同 SSSLs 在德江、貴陽(yáng)和羅平的變異范圍分別為 1.5~6.5 d、1.0~5.0 d 和 1.0~6.5 d。

    表2 SSSLs 開(kāi)花期相關(guān)性狀表型Table 2 Performance of SSSLs population in the three experimental sites

    不施氮條件下,不同 SSSLs 在德江、貴陽(yáng)和羅平 的 散 粉 期 的 變 化 范 圍 分 別 為 85.5~ 97.0 d、76.5~85.5 d 和 70.0~82.0 d;吐絲期的變化范圍分別為 90.0~102.0 d、80.0~88.0 d 和 73.0~85.5 d;ASI 變 化 范 圍 分 別 為2.5~6.0 d、1.5~5.5 d 和1.0~7.5 d。開(kāi)花期相關(guān)性狀在不同處理?xiàng)l件下均表現(xiàn)出較為豐富的表型變異,為后續(xù)的 QTL 定位提供了較為豐富的表型數(shù)據(jù)。

    2.2 控制玉米開(kāi)花期相關(guān)性狀的 QTLs

    對(duì)不同氮素處理?xiàng)l件下的開(kāi)花期性狀 QTL 定位結(jié)果列于表 3。結(jié)果顯示,施氮 (N+) 和不施氮 (N-)兩 種 處 理 條 件 下 , 共 檢 測(cè) 到 54 個(gè) 開(kāi) 花 期 性 狀QTLs,分布在除第 2 和 4 染色體外的其它染色體上。其中在德江試驗(yàn)點(diǎn),兩種氮處理水平下共檢測(cè)到 18 個(gè) QTLs;在貴陽(yáng)試驗(yàn)點(diǎn),兩種氮處理水平下共檢測(cè)到 17 個(gè) QTLs;在羅平試驗(yàn)點(diǎn),兩種氮處理水平下共檢測(cè)到 19 個(gè) QTLs;其中單個(gè) QTL 可解釋目標(biāo)表型變異的變幅為 2.2%~80.0%。

    在施氮條件下,3 個(gè)地點(diǎn)共鑒定出 6 個(gè)控制散粉期的一致性 QTLs,分布于第 1、6、7、10 染色體上,單個(gè) QTL 可解釋的表型變異介于 2.5%~8.8%之間,其中 1 個(gè) QTLqDTP6 在兩個(gè)地點(diǎn) (德江、羅平) 被同時(shí)檢測(cè)到,位于 Bin6.04,在德江和羅平點(diǎn)的貢獻(xiàn)率分別為 8.8% 和 3.4%,且來(lái)自于許 178 的等位基因起增效作用。對(duì)于吐絲期,在 3 個(gè)環(huán)境條件下共鑒定到 7 個(gè) QTLs,分布在第 1、3、5、6、7、9 染色體上,單個(gè) QTL 所解釋的表型變異介于2.5%~7.7% 之間。其中 1 個(gè) QTLqDTS9a 在三個(gè)地點(diǎn)被同時(shí)檢測(cè)到,位于 Bin9.0 區(qū)域,其貢獻(xiàn)率分別為 2.5% (德江)、3.1% (貴陽(yáng)) 和 3.6% (羅平),且來(lái)自許 178 的等位基因起增效作用。對(duì)于 ASI,在 3 個(gè)地點(diǎn)共鑒定到 4 個(gè)共有 QTLs,分布在第 5、9 和 10 染色體上,單個(gè) QTL 所解釋的表型變異介于 13.6%~62.5% 之間。1 個(gè) QTLqASI10a 在三個(gè)地點(diǎn)被同時(shí)檢測(cè)到,位于 Bin10.04 區(qū)域;其貢獻(xiàn)率分別為 13.6% (德江)、50.0% (貴陽(yáng)) 和 28.6% (羅平),來(lái)自于許 178的等位基因起增效作用。

    不施氮條件下,3 個(gè)地點(diǎn)共鑒定到 7 個(gè)散粉期相關(guān) QTLs,分布在第 1、5、6、8、9 和 10 染色體上,單個(gè) QTL 所解釋的表型變異介于 2.2%~6.7%之間。1 個(gè) QTLqDTP9 在三個(gè)地點(diǎn)被同時(shí)檢測(cè)到,位于 Bin9.02,可解釋的表型變異分別為 2.2% (德江)、2.5% (貴陽(yáng)) 和 5.5% (羅平),來(lái)自于許 178 的等位基因起增效作用。在 3 個(gè)地點(diǎn)共鑒定到 8 個(gè)吐絲期 QTL,分布在第 1、3、5、8、9、10 染色體上,單個(gè) QTL 所解釋的表型變異介于 1.8%~7.0% 之間。1 個(gè) QTLqDTS9a 在三個(gè)地點(diǎn)被同時(shí)檢測(cè)到,位Bin9.02,在 3 種環(huán)境的貢獻(xiàn)率分別為 3.2% (德江)、3.0% (貴陽(yáng)) 和 6.0% (羅平),來(lái)自于許 178 的等位基因起增效作用。在 3 個(gè)地點(diǎn)共鑒定到 4 個(gè)散粉-吐絲間隔期 QTL,分布在第 3、9 和 10 染色體上,單個(gè)QTL 所解釋的表型變異介于 12.5%~80.0% 之間,1個(gè) QTLqASI10a 在三個(gè)地點(diǎn)被同時(shí)檢測(cè)到;qASI10a位于 Bin10.04,在 3 種環(huán)境的貢獻(xiàn)率分別為 12.5% (德江)、80.0% (貴陽(yáng)) 和 58.3% (羅平),來(lái)自于許 178的等位基因起增效作用。這些染色體區(qū)段將作為下一步精細(xì)定位和圖位克隆玉米生育期相關(guān)性狀候選基因的重點(diǎn)區(qū)域。

    3 討論

    關(guān)于玉米氮高效利用的 QTL 定位研究前人已有報(bào)道[20-24]。Gallais 等[25]在高氮水平下和低氮水平下檢測(cè)到的 QTL 有很大差別;Agrama 等[21]的研究發(fā)現(xiàn)在低氮和高氮水平下既有相同的 QTL 也有特異表達(dá)的QTL。本研究中,施氮條件下,在染色體 Bin1.03 區(qū)域,檢測(cè)散粉期 QTLqDTP1b 和吐絲期 QTLqDTS1,代換片段均為 phi001-phi095-bnlg182,該研究結(jié)果與前人的定位結(jié)果存在一定的一致性,只是在他們的研究中,該區(qū)域存在一個(gè)分別與株型性狀和絲黑穗病抗性相關(guān)的 QTL-qDH12[26-27],一定程度上反映出生育期與生物量以及抗病性之間是存在一定關(guān)系的。在 Bin6.04 區(qū)域,存在一個(gè)與散粉期和吐絲期顯著相關(guān)的 QTLs (qDTP6 和 qDTS6),代換片段均為umc1014。在 Bin7.05 區(qū)域,檢測(cè)到控制散粉期變異的 QTLs,qDTP7a 和 qDTP7b,吐絲期 QTLqDTS7,而前期報(bào)道則在該區(qū)域檢測(cè)到一個(gè)與株高相關(guān) QTL和葉面積相關(guān) QTL qLA-6[26]。說(shuō)明了在不同氮素處理?xiàng)l件下,一些共有的和特異性遺傳位點(diǎn)共同起用,從而影響玉米氮素利用的效率。

    除此之外,兩種氮處理?xiàng)l件下,共檢測(cè)到 31 個(gè)控制玉米開(kāi)花期相關(guān)性狀的 QTLs。不施氮條件下,在 染 色 體bin8.05 區(qū) 域 , 檢 測(cè) 到2 個(gè) 散 粉 期QTLqDTP8 和 QTLqDTS8,其代換片斷為 bnlg2046-bnlg162-umc1889-umc2401-umc2356-umc1724,在該區(qū) 域 存 在 穗 行 數(shù) 相 關(guān) QTL qKRN8-2QTL[28]。 在Bin9.06 區(qū) 域 , 檢 測(cè) 到1 個(gè) 吐 絲 期 相 關(guān) QTL qDTS9b,代換片段為 umc2345-umc2207-bnlg292;在染色體 Bin10.04 區(qū)域,檢測(cè)到 1 個(gè)吐絲期 QTL qDTS10,代換片段為 umc1077-umc1053-umc2350;與Khairallah 等[29]基于 F2 群體的定位結(jié)果較為一致,他們?cè)谠搮^(qū)域定位了一個(gè)控制 ASI 的 QTL;另外,Ducrocq 等[30]基于不同日照長(zhǎng)短自交系構(gòu)建的近等基因系群體,通過(guò)關(guān)聯(lián)分析的方法也在該區(qū)域定位到一個(gè)與玉米散粉期相關(guān)的主效 QTLs,其貢獻(xiàn)率為40%。進(jìn)一步說(shuō)明了本研究結(jié)果的可靠性。更為重要的是,本研究在兩種氮處理?xiàng)l件下,共檢測(cè)出在 3個(gè)試驗(yàn)環(huán)境中均穩(wěn)定表達(dá)且效應(yīng)較大的一致性QTLs,分別位于 Bin9.02 區(qū)域 (qDTS9a)(umc1170~umc1636~bnlg1401~umc1271) 和 Bin10.04 區(qū)域(qASI10a)(umc1077~umc1053~umc2350)。盡管這些區(qū)域在前人的研究中與穗粗、穗行數(shù)等顯著相關(guān)[27],但是關(guān)于該群段在玉米氮素吸收、轉(zhuǎn)運(yùn)以及利用過(guò)程中的重要性的研究則沒(méi)有。本研究結(jié)果證明了這些區(qū)段在玉米不同氮素處理過(guò)程中均起著重要作用,對(duì)于深入挖掘玉米氮素利用效率將提供更為可靠的分子遺傳學(xué)信息。

    4 結(jié)論

    玉米在不同氮素處理?xiàng)l件下既存在一些共有的主效 QTL,如 qDTS9a、qASI10 在施氮和不施氮條件下均被檢測(cè)到;同時(shí)還存在一些特異性的 QTL 區(qū)段,如 qDTP9 僅在不施氮條件下被檢測(cè)到。這些染色體區(qū)段將為進(jìn)一步改良玉米氮素利用效率以及相關(guān)候選基因的圖位克隆等提供更多的分子遺傳學(xué)支持。

    致謝:貴州省旱糧研究所楊明倫老師在田間管理過(guò)程中給予幫助;貴州省旱糧研究所吳迅博士在文章修改中提出寶貴意見(jiàn)與建議;在此表示感謝。

    [1]Ranmos C. Effect of agricultural practices on the nitrogen losses in the environment [J]. Nutrient Cycling in Agroecosystems, 1995, 43(1): 183-189.

    [2]曹敏建, 衣瑩, 佟占昌, 等. 耐低氮脅迫玉米的篩選 與評(píng)價(jià) [J]. 玉米科學(xué), 2000, 8 (4): 64-69. Cao M J, Yi Y, Tong Z C, et al. Indexes and comprehensive evaluation of low nitrogen tolerance of maize [J]. Journal of Maize Sciences, 2000, 8(4): 64-69.

    [3]劉 建安, 米國(guó)華, 張福鎖. 不 同基因型玉米 氮 效率差異的比 較 研究[J]. 農(nóng)業(yè)生物技術(shù)學(xué)報(bào), 1999, 7 (3): 248-254. Liu J A, Mi G H, Zhang F S. Difference in nitrogen efficiency among maize genotypes[J]. Journal of Agricultural Biotechnology, 1999, 7(3): 248-254.

    [4]劉宗華, 湯繼華, 王春麗, 等. 氮脅迫與非脅迫條件下玉米不同時(shí)期株高的動(dòng)態(tài)QTL定位 [J]. 作物學(xué)報(bào), 2007, 33 (5): 782-789. Liu Z H, Tang J H, Wang C L, et al. QTL analysis of plant height under N-stress and N-input at different stages in maize [J]. Acta Agronomica Sinica, 2007, 33 (5): 782-789.

    [5]鄭祖平, 黃玉碧, 田孟良, 等. 不同供氮水平下玉米株型相關(guān)性狀的QTLs定位和上位性效應(yīng)分析 [J]. 玉米科學(xué), 2007, 15 (2): 16-18. Zheng Z P, Huang Y B, Tian M L, et al. Mapping QTLs and epistasis for plant type traits in maize under two nitrogen levels[J]. Journal of Maize Sciences, 2007, 15(2): 16-18.

    [6]蘭進(jìn)好. 玉米開(kāi)花期相關(guān)性狀的QTL分析 [J]. 西北植物學(xué)報(bào), 2010, 30(3): 471-480. Lan J H. QTL analysis on the flowering related traits in maize [J]. Acta Botanica Boreali-Occidentalia Sinica, 2010, 30 (3): 471-480.

    [7]李賢唐, 丁俊強(qiáng), 王瑞霞, 等. 玉米株型相關(guān)性狀的QTL定位與分析[J]. 江蘇農(nóng)業(yè)科學(xué), 2011, 39 (2): 21-25. Li X T, Ding J Q, Wang R X, et al. Maize plant type traits QTL mapping and analysis[J]. Jiangsu Agricultural Sciences, 2011, 39 (2):21-25.

    [8]Velboom L R, Lee M. Molecular-marker-facilitated studies of morphological traits in maize II. Determination of QTLs for grain yield and yield components [J]. Theoretical and Applied Genetics, 1994, 89 (4): 451-458.

    [9]Sari-Goria M, Krajewski P, Di Fonzo N, et al. Genetic analysis of drough tolerance in maize by molecular markers. II. Plant heigh and flowering [J]. Theoretical and Applied Genetics, 1999, 99: 289-295.

    [10]Ribaut J, Hoisingtond, Deutseh J A. Identification of quantitative trait loci under drought conditions in tropical maize.2. Flowering parameters and the anthesis-silking interval [J]. Theoretical and Applied Genetics, 1996, 92: 905-914.

    [11]Agrama H A, Moussa M E. Mapping QTLs in breeding for drought tolerance in maize (Zea mays L.) [J]. Euphytica, 1996, 91: 89-97.

    [12]Khairallah M M, Bohn M, Jiang C, et al. Molecular mapping of QTL for southwestern corn borer resistance, plant height and flowering in tropical maize[J]. Plant Breeding, 1998, 117(4): 309-318.

    [13]Austin D F, Lee M. Genetic resolution and verification of quantitative trait loci for flowering and plant height with recombinant inbred lines of maize [J]. Genome, 1996, 39: 957-968.

    [14]韓婭楠, 劉福建, 王瑞霞, 等.玉米生育期QTL定位及上位性互作效應(yīng)的遺傳研究[J]. 華北農(nóng)學(xué)報(bào), 2010, 25 (2): 84-87. Han Y N, Liu F J, Wang R X, et al. QTL mapping and epistasis analysis of flowering related traits in maize [J]. Acta Agriculturae Boreali-Sinica, 2010, 25(2): 84-87.

    [15]Li X H, Liu X D, Li M S, et al. Identification of quantitativer trait loci for anthesis-silking interval and yield components under drought stress in maize [J]. Acta Botanica Sinica, 2003, 45(7): 852-857.

    [16]胡彥民, 吳欣, 李翠香, 等. 玉米制種花期相關(guān)性狀的QTL分析 [J].南京農(nóng)業(yè)大學(xué)學(xué)報(bào), 2008, 31(1): 11-16. Hu Y M, Wu X, Li C X, et al. Genetic analysis on the related traits of florescence for hybrid seed production in maize [J]. Journal of Nanjing Agricultural University, 2008, 31 (1): 11-16.

    [17]Gonzalo M, Holland J B, Vyn T J, et al. Direct mapping of density response in a population of B73×Mo17 recombinant inbred lines of maize (Zea Mays L.) [J]. Heredity, 2010, 104: 583-599.

    [18]曹浩飛, 王彬, 毛克舉, 等. 利用單片段代換系對(duì)玉米開(kāi)花期相關(guān)性狀的QTL定位 [J]. 河南農(nóng)業(yè)大學(xué)學(xué)報(bào), 2014, 48 (1): 6-10. Cao H F, Wang B, Mao K J, et al. Mapping of the QTL for flowering related traits in maize using a series of single segment substitution lines [J]. Tournal of Henan Agricultural University, 2014, 48 (1):6-10.

    [19]Guo X, Guo Y P, Ma J, et al. Mapping heterotic loci for yield and agronomic traits using chromosome segment introgression lines in cotton[J]. Journal of Integrative Plant Biology, 2013, 55 (8):759-774.

    [20]Lance R V, Lee M. Genetic mapping of quantitative trait loci in maize in stress and nonstress environments: Ⅰ. Grain yield and yield components[J]. Crop Science, 1996, 36: 1310-1319.

    [21]Agrama H A S, Zakaria A G, Said F B, Tuinstra M. Identification of quantitative trait loci for nitrogen use efficiency in maize [J]. Molecular Breeding, 1999, 5: 187-195.

    [22]Hirel B, Pascal B, Isabelle Q, et al. Towards a better understanding of the genetic and physiological basis for nitrogen use efficiency in maize [J]. Plant Physiology, 2001, 125(3): 1258-1270.

    [23]Olivier L, Sylvain C, Patricia M, et al. Quantitative trait Loci analysis of nitrogen use efficiency in Arabidopsis [J]. Plant Physiology, 2003, 131: 345-358.

    [24]Coque M, Gallais A. Genomic regions involved in response to grain yield selection at high and low nitrogen fertilization in maize [J]. Theoretical and Applied Genetics, 2006, 112: 1205-1220.

    [25]Gallais A, Hirel B. An approach to the genetics of nitrogen use efficiency in maize [J]. Journal of Experimental Botany, 2004, 55(396): 295-306.

    [26]許誠(chéng), 王彬, 毛克舉, 等. 利用單片段代換系群體定位玉米株型性狀QTL [J]. 玉米科學(xué), 2014, 22 (2): 28-34. Xu C, Wang B, Mao K J, et al. QTL mapping for plant-type related traits using single segment substitution lines in maize [J]. Journal of Maize Sciences, 2014, 22(2): 28-34.

    [27]張書(shū)紅, 張世煌, 李新海, 等. 玉米抗病基因一致性圖譜的構(gòu)建 [J].中國(guó)農(nóng)學(xué)通報(bào), 2007, 23 (6): 601-606. Zhang S H, Zhang S H, Li X H, et al. Construction of consensus map of gene in maize [J]. Chinese Agricultural Science Bulletin, 2007, 23(6): 601-606.

    [28]焦付超, 李永祥, 陳林, 等. 特異玉米種質(zhì)四路糯的穗行數(shù)遺傳解析[J]. 中國(guó)農(nóng)業(yè)科學(xué), 2014, 47(7): 1256-1264. Jiao F C, Li Y X , Chen L, et al. Genetic dissection for kernel row number in the specific maize germplasm four-rowed waxy corn [J]. Scientia Agricultura Sinica, 2014, 47(7): 1256-1264.

    [29]Khairallah M M, Bohn M, Jiang C, et al. Molecular mapping of QTL for southwestern corn borer resistance, plant height and flowering in tropical maize [J]. Plant Breeding, 1998, 117(4): 309-318.

    [30]Ducrocq S, Giauffret C, Madur D, et al. Fine mapping and haplotype structure analysis of a major flowering time quantitative trait locus on maize chromosome 10 [J]. Genetics, 2009, 183: 1555-1563.

    QTL mapping of flowering related traits of maize with and without nitrogen application

    GUO Xiang-yang1, CHEN Jian-jun2, WEI Xiao-yi3,4, ZHU Yun-fang1, WANG An-gui1,
    LIU Peng-fei1, TANG Ji-hua3*, CHEN Ze-hui1*
    ( 1 Guizhou Institute of Upland Crops, Guiyang 550006, China; 2 Huihuang Seed Industry Co., LTD, Zunyi Guizhou 563100, China; 3 College of Agronomy, Henan Agricultural University, Zhengzhou 450002, China; 4 Xinxiang Academy of Agricultural Sciences of Henan Province, Xinxiang Henan 453000, China )

    【Objectives】The purpose of this study was to dissect variation of maize flowering-related traits under two conditions of nitrogen applied normally or none, and to map the related QTLs. 【Methods】 One panel included 160 Single Segment Substitution Lines (SSSLs) was constructed by using foundation parental lines of Zong3 and Xu178 as donor and receptor parents, respectively. The panel was treated with two conditions of nitrogen applied normally or none by using the split-plot experiment design, and the nitrogen treatment was set as the main plot and genotypes were set as the sub-plots under three environments, including Guiyang, Dejiang in Guizhou Province and Luoping in Yunnan Province. Three flowering-related traits, including day to pollen (DTP),day to silk (DTS) and anthesis-silking interval (ASI), were investigated in one treatment. 【Results】The results showed that across the whole genome, about 54 flowering-related traits QTLs were identified under the two N treatments, which were mainly located on the chromosome 1, 3, 6, 9 and 10, among them, 3 QTLs were common under three environments. Under N+treatment, QTL qDTS9a controlling day to silling (DTS) was found which was located on chromosome 9, and could explain 3.05% genetic variation; another QTL qASI10a controlling the ASI was also found, which located on chromosome 10, and could explain 30.74% genetic variation. Under N-conditions, QTL qDTP9 was identified, which could explain 3.43% genetic variation of DTP; qDTS9a related with DTS was found, which located on chromosome 9, and could explain 4.08% genetic variation; QTL qASI10 was also identified, which located on chromosome 10, and could explain 50.28% genetic variation. Comparing between N+ and N- conditions showed that, the QTLs identified under N- was much more than those under N+ treatment. 【Conclusions】Two commonly QTLs are existed in the two nitrogen treatments, one is qDTS9a, which is located on Bin9.02, with flank markers of umc1170-umc1636-bnlg1401-umc1271, and the other is qASI10a, which is located on Bin 10.04, with flank markers of umc1077-umc1053-umc2350. These loci may play important role in nitrogen absorbing, transporting, and utilization during maize development, and served as candidate loci in further map-based cloning of maize flowering-related traits.

    maize; single segment substitution lines (SSSLs); nitrogen levels; flowering related traits; quantitative trait locus (QTLs)

    2016-07-18 接受日期:2016-11-22

    貴州省農(nóng)業(yè)科學(xué)院自主創(chuàng)新項(xiàng)目(黔農(nóng)科院自主創(chuàng)新科研專(zhuān)項(xiàng)字[2014]006號(hào));貴州省科技支撐計(jì)劃項(xiàng)目(黔科合支撐[2016]2605號(hào));國(guó)家“七大作物育種”專(zhuān)項(xiàng)(2016YFD0101206-4)資助。

    郭向陽(yáng)(1982—),男,河南滑縣人,博士研究生,副研究員,主要從事玉米遺傳育種研究。E-mail:xyguo0372@163.com

    * 通信作者 E-mail:Tangjihua1@163.com;Tel:0851-83760096,E-mail:chenzh907@sohu.com

    猜你喜歡
    吐絲開(kāi)花期氮素
    高寒草原針茅牧草花期物候變化特征及其影響因子分析
    高速線材廠吐絲管固定方式的改進(jìn)
    吐絲的蜘蛛與吐絲的蠶
    意林(2018年18期)2018-09-30 02:44:00
    小蠶兒吐絲
    會(huì)吐絲的蠶寶寶
    大灰狼(2017年1期)2017-03-01 12:34:08
    SOLVABILITY OF A PARABOLIC-HYPERBOLIC TYPE CHEMOTAXIS SYSTEM IN 1-DIMENSIONAL DOMAIN?
    初春氣象條件對(duì)蘋(píng)果開(kāi)花期的影響分析
    楸樹(shù)無(wú)性系苗期氮素分配和氮素效率差異
    施氮量對(duì)不同開(kāi)花期棉(Gossypium hirsutum L.)鈴纖維細(xì)度和成熟度形成的影響
    基于光譜分析的玉米氮素營(yíng)養(yǎng)診斷
    可以在线观看的亚洲视频| 中文字幕av在线有码专区| 亚洲最大成人中文| av国产免费在线观看| 国产av不卡久久| 久久天躁狠狠躁夜夜2o2o| 三级毛片av免费| 精品久久久久久久久久免费视频| 欧美最新免费一区二区三区| 搡老岳熟女国产| 黄色视频,在线免费观看| 禁无遮挡网站| 亚洲av成人精品一区久久| 一个人看的www免费观看视频| 亚洲国产精品sss在线观看| 久久久久久久午夜电影| 成人性生交大片免费视频hd| 亚洲三级黄色毛片| 欧美一区二区精品小视频在线| 国产精品免费一区二区三区在线| 亚洲欧美清纯卡通| 国产av一区在线观看免费| a级毛片a级免费在线| 午夜免费激情av| 国产久久久一区二区三区| 亚洲欧美日韩东京热| 波多野结衣高清无吗| 亚洲熟妇熟女久久| 男插女下体视频免费在线播放| 蜜桃久久精品国产亚洲av| 日韩国内少妇激情av| 亚洲乱码一区二区免费版| 国产精品日韩av在线免费观看| 国产亚洲精品久久久com| 亚洲精品乱码久久久v下载方式| 国产探花极品一区二区| 男人舔女人下体高潮全视频| 99国产精品一区二区蜜桃av| 老司机福利观看| 伊人久久精品亚洲午夜| 国内揄拍国产精品人妻在线| 中国美女看黄片| 校园春色视频在线观看| 欧美日韩瑟瑟在线播放| 国产免费av片在线观看野外av| 69av精品久久久久久| 久久精品国产99精品国产亚洲性色| 中文资源天堂在线| 美女免费视频网站| 男女之事视频高清在线观看| 亚洲美女搞黄在线观看 | 国产乱人伦免费视频| 成人三级黄色视频| 99久久无色码亚洲精品果冻| 男女那种视频在线观看| 一区二区三区高清视频在线| 九九热线精品视视频播放| 狂野欧美白嫩少妇大欣赏| 精品国内亚洲2022精品成人| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 免费看日本二区| 两性午夜刺激爽爽歪歪视频在线观看| 美女大奶头视频| 精品免费久久久久久久清纯| 精品久久久久久久久久久久久| 波野结衣二区三区在线| 欧美黑人欧美精品刺激| 能在线免费观看的黄片| 国产成人一区二区在线| eeuss影院久久| 岛国在线免费视频观看| 日本在线视频免费播放| 亚洲人成网站在线播放欧美日韩| 长腿黑丝高跟| 国产麻豆成人av免费视频| 精品久久久久久久人妻蜜臀av| 成人毛片a级毛片在线播放| 久久人人爽人人爽人人片va| 精品无人区乱码1区二区| 色综合婷婷激情| 亚洲精品在线观看二区| 91在线精品国自产拍蜜月| 老司机福利观看| 国产白丝娇喘喷水9色精品| 成人国产麻豆网| 国产精品久久久久久av不卡| 亚洲人成网站高清观看| 在线观看av片永久免费下载| 我要搜黄色片| 久久精品国产亚洲av天美| 亚洲精品一卡2卡三卡4卡5卡| 国产免费av片在线观看野外av| 深爱激情五月婷婷| 欧美成人免费av一区二区三区| 久久久久国产精品人妻aⅴ院| 久久精品久久久久久噜噜老黄 | 少妇熟女aⅴ在线视频| 午夜福利18| 校园人妻丝袜中文字幕| 又爽又黄a免费视频| 免费观看在线日韩| 成年女人永久免费观看视频| 免费在线观看成人毛片| 国产精品国产高清国产av| 热99在线观看视频| 在线免费十八禁| 国产激情偷乱视频一区二区| 人人妻人人澡欧美一区二区| 成人无遮挡网站| av黄色大香蕉| 在线观看av片永久免费下载| 亚洲欧美日韩卡通动漫| 久久久久久久亚洲中文字幕| 精品人妻1区二区| 一个人免费在线观看电影| 精品99又大又爽又粗少妇毛片 | 日韩亚洲欧美综合| 国产探花在线观看一区二区| 国产真实伦视频高清在线观看 | 国产极品精品免费视频能看的| 国产91精品成人一区二区三区| 给我免费播放毛片高清在线观看| 中文字幕高清在线视频| 国产高清激情床上av| 午夜老司机福利剧场| 国产在线精品亚洲第一网站| 欧美又色又爽又黄视频| 97热精品久久久久久| 最近最新中文字幕大全电影3| 日本黄大片高清| 国产精品日韩av在线免费观看| 久久99热6这里只有精品| 日韩,欧美,国产一区二区三区 | 十八禁国产超污无遮挡网站| 久久久午夜欧美精品| 国产黄色小视频在线观看| 99热这里只有精品一区| 国产高潮美女av| 偷拍熟女少妇极品色| 欧美成人免费av一区二区三区| 三级男女做爰猛烈吃奶摸视频| 99热只有精品国产| 色视频www国产| 欧美日韩瑟瑟在线播放| 春色校园在线视频观看| 国产色婷婷99| 嫩草影视91久久| 日本免费a在线| 亚洲精品乱码久久久v下载方式| 成年女人看的毛片在线观看| 天美传媒精品一区二区| 色av中文字幕| 在线观看av片永久免费下载| 天天躁日日操中文字幕| 看黄色毛片网站| 国产午夜精品久久久久久一区二区三区 | 3wmmmm亚洲av在线观看| 国产成人影院久久av| 国产欧美日韩精品一区二区| 久久精品国产亚洲av香蕉五月| 日本一二三区视频观看| 国产伦一二天堂av在线观看| 99热只有精品国产| 久久久国产成人精品二区| 国产黄片美女视频| 亚洲av日韩精品久久久久久密| 国产精品久久久久久久久免| av在线蜜桃| 中文字幕熟女人妻在线| 午夜久久久久精精品| 一进一出抽搐gif免费好疼| 深爱激情五月婷婷| 黄色欧美视频在线观看| 一区福利在线观看| 热99re8久久精品国产| 免费看光身美女| 美女黄网站色视频| 桃色一区二区三区在线观看| 国产精品1区2区在线观看.| 十八禁网站免费在线| 日本黄色片子视频| 免费在线观看影片大全网站| 日韩av在线大香蕉| netflix在线观看网站| 国产精品av视频在线免费观看| 男人和女人高潮做爰伦理| 中文资源天堂在线| 色综合亚洲欧美另类图片| 欧美xxxx性猛交bbbb| 伦精品一区二区三区| 久久久成人免费电影| 一进一出抽搐动态| 亚洲人成网站在线播放欧美日韩| 午夜亚洲福利在线播放| 成人一区二区视频在线观看| 美女免费视频网站| 少妇人妻精品综合一区二区 | 国产单亲对白刺激| 久久久久久九九精品二区国产| 一夜夜www| 国产 一区 欧美 日韩| 日本五十路高清| 欧美xxxx黑人xx丫x性爽| 美女被艹到高潮喷水动态| 3wmmmm亚洲av在线观看| 午夜精品久久久久久毛片777| 国产精品乱码一区二三区的特点| 亚洲内射少妇av| 精品午夜福利视频在线观看一区| 久久欧美精品欧美久久欧美| 99久久精品国产国产毛片| 国产亚洲精品av在线| 高清日韩中文字幕在线| 如何舔出高潮| 国产av在哪里看| 日韩欧美精品免费久久| 国产精品人妻久久久久久| 久久婷婷人人爽人人干人人爱| 乱系列少妇在线播放| 色综合色国产| 制服丝袜大香蕉在线| 又黄又爽又刺激的免费视频.| 国产一区二区三区av在线 | 69人妻影院| 亚洲欧美日韩高清在线视频| 日韩欧美一区二区三区在线观看| 99久国产av精品| 身体一侧抽搐| 成人无遮挡网站| 日日摸夜夜添夜夜添av毛片 | 午夜视频国产福利| 国产精品国产三级国产av玫瑰| 干丝袜人妻中文字幕| 免费观看在线日韩| 天天躁日日操中文字幕| 日韩国内少妇激情av| 99在线人妻在线中文字幕| 亚洲男人的天堂狠狠| 欧美+亚洲+日韩+国产| 一边摸一边抽搐一进一小说| 免费看光身美女| 有码 亚洲区| 精品午夜福利视频在线观看一区| 国产精品久久久久久亚洲av鲁大| 偷拍熟女少妇极品色| 欧美另类亚洲清纯唯美| 免费看a级黄色片| 亚洲最大成人av| 精品日产1卡2卡| 22中文网久久字幕| 久久久久九九精品影院| 亚洲最大成人中文| 免费大片18禁| netflix在线观看网站| 丰满人妻一区二区三区视频av| 女人被狂操c到高潮| 麻豆久久精品国产亚洲av| 亚洲欧美日韩卡通动漫| 国产精品99久久久久久久久| 哪里可以看免费的av片| 亚洲av成人精品一区久久| 午夜福利欧美成人| 国产成人av教育| av.在线天堂| 日本a在线网址| av专区在线播放| 久久婷婷人人爽人人干人人爱| 夜夜夜夜夜久久久久| 亚洲中文字幕一区二区三区有码在线看| av在线天堂中文字幕| 国产精品一区二区免费欧美| 日日撸夜夜添| 亚洲不卡免费看| av中文乱码字幕在线| 日韩亚洲欧美综合| av天堂中文字幕网| 欧美zozozo另类| 免费av观看视频| 欧美精品国产亚洲| 中文在线观看免费www的网站| 国产午夜精品久久久久久一区二区三区 | 69人妻影院| 人人妻人人看人人澡| 97超视频在线观看视频| 中文字幕免费在线视频6| 日韩精品青青久久久久久| 桃红色精品国产亚洲av| 国产精品亚洲一级av第二区| 久久亚洲真实| 亚洲av五月六月丁香网| 91在线精品国自产拍蜜月| 国产亚洲精品综合一区在线观看| 国国产精品蜜臀av免费| 日日摸夜夜添夜夜添小说| 国产精品免费一区二区三区在线| 精品人妻一区二区三区麻豆 | 亚洲av.av天堂| 精品日产1卡2卡| 国产男靠女视频免费网站| av女优亚洲男人天堂| 久久精品久久久久久噜噜老黄 | 日本色播在线视频| 啦啦啦啦在线视频资源| 午夜免费成人在线视频| 中文字幕精品亚洲无线码一区| 午夜福利在线观看免费完整高清在 | 午夜福利18| 国产麻豆成人av免费视频| 999久久久精品免费观看国产| 久久亚洲真实| 亚洲色图av天堂| 一夜夜www| 中文资源天堂在线| 极品教师在线免费播放| 日韩欧美在线乱码| 久久久久久大精品| 人人妻人人看人人澡| 亚洲av成人av| 联通29元200g的流量卡| 免费高清视频大片| 亚洲欧美日韩高清在线视频| 香蕉av资源在线| 亚洲avbb在线观看| 久久亚洲精品不卡| 欧美国产日韩亚洲一区| 久久久成人免费电影| 校园春色视频在线观看| 成人国产麻豆网| 特大巨黑吊av在线直播| 日本色播在线视频| 精品日产1卡2卡| 深夜a级毛片| aaaaa片日本免费| 成人毛片a级毛片在线播放| 国内精品宾馆在线| av在线蜜桃| 久久九九热精品免费| 九九在线视频观看精品| 春色校园在线视频观看| 永久网站在线| 久久久国产成人免费| 亚洲中文字幕一区二区三区有码在线看| 最后的刺客免费高清国语| 波野结衣二区三区在线| 中国美女看黄片| 国产精品1区2区在线观看.| 69av精品久久久久久| 天天一区二区日本电影三级| 日本a在线网址| 欧美人与善性xxx| 久久久久久大精品| 欧美日本视频| 国产精品女同一区二区软件 | 久久6这里有精品| 日韩国内少妇激情av| 一级毛片久久久久久久久女| 搞女人的毛片| 日本a在线网址| 两性午夜刺激爽爽歪歪视频在线观看| 欧美成人性av电影在线观看| 国产毛片a区久久久久| 一夜夜www| 久久精品国产亚洲av涩爱 | 麻豆国产97在线/欧美| 免费在线观看日本一区| 亚洲内射少妇av| 成人一区二区视频在线观看| 91久久精品国产一区二区三区| 亚洲av中文字字幕乱码综合| 91久久精品国产一区二区三区| 国产在线男女| 国产一区二区三区在线臀色熟女| av在线亚洲专区| 国内毛片毛片毛片毛片毛片| 欧美精品啪啪一区二区三区| 久久6这里有精品| 成人国产综合亚洲| 毛片一级片免费看久久久久 | 亚洲狠狠婷婷综合久久图片| 国产精品99久久久久久久久| 69av精品久久久久久| 老熟妇乱子伦视频在线观看| 97热精品久久久久久| 国产精品亚洲一级av第二区| 日韩亚洲欧美综合| 久久精品人妻少妇| 国产视频内射| 精品久久久久久,| 夜夜夜夜夜久久久久| 亚洲自拍偷在线| 日本黄色片子视频| 色噜噜av男人的天堂激情| 精品久久久久久成人av| 他把我摸到了高潮在线观看| 国产精品一区二区免费欧美| 有码 亚洲区| 成人精品一区二区免费| 乱人视频在线观看| 亚洲va日本ⅴa欧美va伊人久久| 搡老妇女老女人老熟妇| 色综合色国产| 亚洲va在线va天堂va国产| 亚洲国产高清在线一区二区三| 国内精品宾馆在线| 亚洲 国产 在线| x7x7x7水蜜桃| 日韩欧美精品v在线| 又爽又黄无遮挡网站| eeuss影院久久| 亚洲国产精品久久男人天堂| 97超视频在线观看视频| 国产伦精品一区二区三区四那| 亚洲国产欧美人成| 国产伦在线观看视频一区| 搡老妇女老女人老熟妇| 欧美最新免费一区二区三区| 午夜精品在线福利| 波多野结衣高清无吗| 麻豆av噜噜一区二区三区| 国产精品无大码| 亚洲人成网站在线播| 久久精品国产亚洲av香蕉五月| 乱码一卡2卡4卡精品| 丰满人妻一区二区三区视频av| 国产色婷婷99| 国产精品自产拍在线观看55亚洲| 国产精品乱码一区二三区的特点| 午夜福利在线观看吧| 欧美三级亚洲精品| 热99re8久久精品国产| 色在线成人网| 别揉我奶头 嗯啊视频| 在线免费观看的www视频| 中文字幕久久专区| 国产精品爽爽va在线观看网站| 欧美+日韩+精品| 亚洲黑人精品在线| 成年女人看的毛片在线观看| 少妇人妻精品综合一区二区 | 欧美日韩国产亚洲二区| 91久久精品国产一区二区成人| h日本视频在线播放| 高清毛片免费观看视频网站| 国产乱人视频| 成年女人毛片免费观看观看9| 国产主播在线观看一区二区| 又黄又爽又刺激的免费视频.| 国产一区二区在线av高清观看| 日韩在线高清观看一区二区三区 | 免费一级毛片在线播放高清视频| av在线蜜桃| 乱系列少妇在线播放| 午夜影院日韩av| 国产精品久久电影中文字幕| 久久久久久久久久黄片| 小蜜桃在线观看免费完整版高清| 一区二区三区四区激情视频 | 一个人看视频在线观看www免费| 少妇猛男粗大的猛烈进出视频 | 日本熟妇午夜| 久久久久国内视频| 国产乱人视频| 女的被弄到高潮叫床怎么办 | 国产欧美日韩精品一区二区| 一个人免费在线观看电影| 婷婷精品国产亚洲av| 毛片女人毛片| 男人舔奶头视频| 久久精品国产99精品国产亚洲性色| 九九久久精品国产亚洲av麻豆| 十八禁国产超污无遮挡网站| 国产一区二区三区在线臀色熟女| 白带黄色成豆腐渣| 最近在线观看免费完整版| 日日啪夜夜撸| 蜜桃久久精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 亚洲午夜理论影院| 天天躁日日操中文字幕| 日韩欧美 国产精品| 成人鲁丝片一二三区免费| 男女做爰动态图高潮gif福利片| 99热这里只有是精品50| 欧美性猛交╳xxx乱大交人| 欧美+亚洲+日韩+国产| 久久精品国产99精品国产亚洲性色| 九九久久精品国产亚洲av麻豆| 欧美最黄视频在线播放免费| 制服丝袜大香蕉在线| 欧美激情久久久久久爽电影| 国产亚洲精品久久久com| 午夜福利在线在线| 蜜桃久久精品国产亚洲av| 国产高清视频在线观看网站| 韩国av一区二区三区四区| 91久久精品国产一区二区成人| 欧美日韩黄片免| 成年免费大片在线观看| www.色视频.com| 欧美最黄视频在线播放免费| 最近中文字幕高清免费大全6 | 欧美激情在线99| 亚洲国产高清在线一区二区三| 人妻少妇偷人精品九色| 欧美日韩精品成人综合77777| av黄色大香蕉| 国产精品爽爽va在线观看网站| 看十八女毛片水多多多| 精品人妻熟女av久视频| 男人舔女人下体高潮全视频| 成人一区二区视频在线观看| 岛国在线免费视频观看| 麻豆久久精品国产亚洲av| 又爽又黄无遮挡网站| 国产精品免费一区二区三区在线| 在线观看舔阴道视频| 成人无遮挡网站| 国产三级中文精品| 色尼玛亚洲综合影院| 麻豆国产97在线/欧美| 人妻制服诱惑在线中文字幕| 真人一进一出gif抽搐免费| 国产真实伦视频高清在线观看 | 嫩草影院精品99| av天堂中文字幕网| 国产中年淑女户外野战色| 亚洲美女视频黄频| 成人精品一区二区免费| 日本 欧美在线| 又粗又爽又猛毛片免费看| 精品久久久久久成人av| 久久久久久国产a免费观看| 欧美日韩瑟瑟在线播放| 成人av在线播放网站| 午夜福利欧美成人| 可以在线观看毛片的网站| av视频在线观看入口| 国产女主播在线喷水免费视频网站 | 欧美极品一区二区三区四区| 亚洲内射少妇av| 婷婷亚洲欧美| 国产人妻一区二区三区在| av国产免费在线观看| 免费观看精品视频网站| 丝袜美腿在线中文| 日韩人妻高清精品专区| 久久久久久久久久久丰满 | 亚洲人成网站高清观看| 国产黄色小视频在线观看| 身体一侧抽搐| 亚洲av电影不卡..在线观看| 国产一区二区亚洲精品在线观看| 国产精品98久久久久久宅男小说| 少妇人妻精品综合一区二区 | av福利片在线观看| 美女 人体艺术 gogo| 亚洲美女视频黄频| 日韩在线高清观看一区二区三区 | 黄色配什么色好看| 99精品久久久久人妻精品| 亚洲精品成人久久久久久| 午夜精品一区二区三区免费看| 国产高清有码在线观看视频| 老熟妇仑乱视频hdxx| 亚洲av电影不卡..在线观看| 色综合站精品国产| 在线a可以看的网站| 一级黄片播放器| 国产精品,欧美在线| 久久精品国产亚洲av涩爱 | 久久天躁狠狠躁夜夜2o2o| 婷婷精品国产亚洲av在线| 中文字幕人妻熟人妻熟丝袜美| 乱人视频在线观看| 日本-黄色视频高清免费观看| 蜜桃亚洲精品一区二区三区| 久久精品影院6| 国产乱人伦免费视频| 免费在线观看成人毛片| 深夜精品福利| 色综合色国产| 国内精品久久久久久久电影| 久久久久久久久久黄片| 欧美日韩瑟瑟在线播放| 国产成人av教育| 久久中文看片网| 内射极品少妇av片p| 亚洲欧美日韩高清专用| 热99在线观看视频| 中国美白少妇内射xxxbb| 男人和女人高潮做爰伦理| 日韩精品有码人妻一区| netflix在线观看网站| 亚洲四区av| 欧美日韩精品成人综合77777| 午夜视频国产福利| 久久久久九九精品影院| www日本黄色视频网| 少妇的逼水好多| 午夜老司机福利剧场| 91av网一区二区| 午夜视频国产福利| 一本久久中文字幕| eeuss影院久久| 国产亚洲精品综合一区在线观看| 免费观看的影片在线观看| 五月伊人婷婷丁香| 中文字幕高清在线视频| 国产亚洲精品久久久com| 午夜a级毛片| 亚洲无线观看免费| 麻豆成人午夜福利视频| 免费av毛片视频| 老司机深夜福利视频在线观看|