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

    利用抗旱選擇導(dǎo)入系定位向日葵產(chǎn)量性狀QTL

    2018-03-13 06:54:26于海峰侯建華
    作物學(xué)報(bào) 2018年3期
    關(guān)鍵詞:供體親本向日葵

    呂 品 于海峰 侯建華,*

    ?

    利用抗旱選擇導(dǎo)入系定位向日葵產(chǎn)量性狀QTL

    呂 品1于海峰2侯建華1,*

    1內(nèi)蒙古農(nóng)業(yè)大學(xué)農(nóng)學(xué)院;2內(nèi)蒙古農(nóng)牧業(yè)科學(xué)院作物育種與栽培研究所, 內(nèi)蒙古呼和浩特 010010

    干旱是造成向日葵減產(chǎn)的最主要因素之一。利用綜合性狀優(yōu)良的自交系K55作為輪回親本與抗旱自交系K58雜交構(gòu)建回交導(dǎo)入系, 在干旱條件下進(jìn)行單株產(chǎn)量篩選, 得到45個(gè)BC3F2抗旱定向選擇導(dǎo)入系。通過(guò)全基因組SSR及SNP標(biāo)記掃描, 以方差分析和基于遺傳搭車(chē)原理的卡方檢驗(yàn)對(duì)呼和浩特市及武川縣兩點(diǎn)、兩種水分條件下的5個(gè)產(chǎn)量性狀進(jìn)行QTL檢測(cè)。方差分析檢測(cè)到的QTL根據(jù)不同環(huán)境下的表達(dá)情況分為三類(lèi), 第一類(lèi)在兩種水分條件下穩(wěn)定表達(dá), 包括武川的4個(gè)百粒重QTL及呼和浩特的2個(gè)單株產(chǎn)量QTL、3個(gè)單株實(shí)粒數(shù)QTL, 這些QTL可能對(duì)向日葵抗旱性有直接貢獻(xiàn); 第二類(lèi)受干旱脅迫表達(dá), 包括呼和浩特的30個(gè)和武川的27個(gè); 第三類(lèi)僅在正常供水條件下被檢測(cè)到, 包括呼和浩特的38個(gè)和武川的64個(gè)??ǚ綑z驗(yàn)檢測(cè)到極顯著位點(diǎn)274個(gè)。用兩種方法共檢測(cè)到一致性位點(diǎn)14個(gè), 可能是與向日葵抗旱性相關(guān)的關(guān)鍵位點(diǎn)。本研究結(jié)果可為向日葵高效抗旱分子育種奠定基礎(chǔ)并提供相關(guān)材料。

    向日葵; 抗旱性; 選擇回交導(dǎo)入系; 產(chǎn)量性狀; QTL定位

    我國(guó)是向日葵種植大國(guó), 西北、華北及東北的干旱、半干旱地區(qū)是主要的向日葵種植區(qū), 這些地區(qū)降水量不足且時(shí)空分布不均, 干旱是向日葵減產(chǎn)的最主要因素之一。干旱對(duì)作物的影響最終體現(xiàn)在產(chǎn)量上, 直接以產(chǎn)量進(jìn)行選擇可以有效提高作物的抗旱性[1-2]。產(chǎn)量及抗旱性均屬于多基因控制的數(shù)量性狀, 隨著分子標(biāo)記技術(shù)的發(fā)展, 可定位目標(biāo)數(shù)量性狀位點(diǎn)(QTL), 通過(guò)QTL的準(zhǔn)確追蹤, 挖掘與耐旱基因緊密關(guān)聯(lián)的分子標(biāo)記, 促進(jìn)分子標(biāo)記輔助選育(MAS), 可以加快培育耐旱品種的進(jìn)程[3]。

    向日葵抗旱QTL定位研究多利用自交系群體進(jìn)行[4-7]。自交系定位易受背景干擾, 研究群體脫離育種實(shí)踐, 再加上產(chǎn)量及抗旱性本身的復(fù)雜性, 存在定位精度不高的問(wèn)題。高代回交導(dǎo)入系群體可有效克服此類(lèi)缺陷, 尤其是Li等[8]提出的選擇導(dǎo)入系概念, 推動(dòng)了作物QTL準(zhǔn)確定位與基因型分析的發(fā)展[9]。利用選擇導(dǎo)入系進(jìn)行QTL定位的研究在水稻、大豆中報(bào)道較多[10-18], 卓壯[11]利用黃華占為輪回親本構(gòu)建的兩個(gè)抗旱選擇群體分別定位到20個(gè)和11個(gè)產(chǎn)量性狀位點(diǎn); 趙秀琴等[12]利用Lemont導(dǎo)入到特青背景的高代回交導(dǎo)入系在2種水分條件下檢測(cè)到籽粒產(chǎn)量及其穗部相關(guān)性狀QTL 32個(gè), 其中對(duì)耐旱性可能有直接貢獻(xiàn)的為在兩種條件下均檢測(cè)到的3個(gè)QTL (、和)及影響兩種條件下性狀差值的9個(gè)QTL; 李燦東等[17]對(duì)36個(gè)耐旱定向選擇導(dǎo)入系進(jìn)行方差分析共檢測(cè)到14個(gè)控制相對(duì)發(fā)芽率的QTL, 其中Satt156、Satt423、Sattt167等位點(diǎn)與卡方測(cè)驗(yàn)檢測(cè)的結(jié)果具有一致性。但利用抗旱選擇群體進(jìn)行QTL的定位在向日葵上還未見(jiàn)報(bào)道。本研究利用油用向日葵高代回交抗旱選擇導(dǎo)入系群體, 以SSR及SNP標(biāo)記進(jìn)行基因型分析, 在干旱脅迫和正常供水條件下定位產(chǎn)量性狀的QTL, 研究結(jié)果可為向日葵抗旱及產(chǎn)量性狀QTL的精細(xì)定位、克隆及向日葵高效抗旱分子育種奠定基礎(chǔ)。

    1 材料與方法

    1.1 抗旱選擇回交導(dǎo)入群體的構(gòu)建

    選取油用向日葵抗旱自交系K58為供體親本、綜合性狀優(yōu)良的非抗旱自交系K55為輪回親本雜交后再回交2代, 產(chǎn)生BC2F1群體。2013年春季, 將BC2F1群體種植于內(nèi)蒙古農(nóng)業(yè)大學(xué)試驗(yàn)農(nóng)場(chǎng), 與輪回親本K55繼續(xù)回交并進(jìn)行抗旱初次篩選, 田間干旱脅迫處理為作水播種后全生育期雨養(yǎng), 成熟期收獲單株產(chǎn)量高于輪回親本的單株, 共得到80個(gè)BC3F1單株。2013年冬季, 將80個(gè)BC3F1單株材料種植于海南, 進(jìn)一步抗旱篩選, 田間處理同2013年春季, 自交套袋, 成熟期收獲產(chǎn)量高于輪回親本的單株, 共45株作為試驗(yàn)材料。

    1.2 田間試驗(yàn)設(shè)計(jì)

    2014年和2015年春季, 將篩選到的45個(gè)BC3F2株系及其親本分別種植于內(nèi)蒙古呼和浩特市內(nèi)蒙古農(nóng)業(yè)大學(xué)農(nóng)場(chǎng)(兩年間全生育期降雨量分別為364.9 mm和274.9 mm)和內(nèi)蒙古武川縣農(nóng)業(yè)部?jī)?nèi)蒙古耕地保育科學(xué)觀測(cè)實(shí)驗(yàn)站(兩年間全生育期降雨量分別為366.1 mm和307.5 mm), 兩點(diǎn)均采用隨機(jī)區(qū)組設(shè)計(jì), 3次重復(fù), 2個(gè)處理(干旱脅迫和正常供水條件, 兩處理間設(shè)置隔離區(qū)), 行長(zhǎng)4 m, 行株距0.50 m× 0.35 m, 干旱脅迫處理為作水播種后全生育期完全雨養(yǎng), 正常供水條件為播種期、現(xiàn)蕾期和花末期灌水3次。

    1.3 產(chǎn)量測(cè)定

    考察2點(diǎn)(呼和浩特和武川)、2種水分條件下(干旱脅迫和正常供水)的5個(gè)產(chǎn)量性狀, 成熟后每小區(qū)選取5株測(cè)量花盤(pán)直徑(DS, cm), 單株收獲。收獲后單株脫粒曬干至種子含水量8%以下, 測(cè)量單株產(chǎn)量(SY, g plant-1)、百粒重(HSW, g), 數(shù)取單株實(shí)粒數(shù)(FSP), 計(jì)算結(jié)實(shí)率(SSR, %), 結(jié)實(shí)率=單株實(shí)粒數(shù)/ (單株實(shí)粒數(shù)+空粒數(shù))×100%。

    1.4 基因型分析

    選取向日葵苗期鮮嫩真葉, 用TIANGEN試劑盒提取雙親及選擇導(dǎo)入系群體的基因組DNA。利用已有在親本K55、K58間具有多態(tài)性的177對(duì)SSR引物[7]對(duì)抗旱選擇群體進(jìn)行標(biāo)記分析, 參照房冬梅改進(jìn)的PCR體系及程序 [19], 由北京百邁克生物科技有限公司對(duì)雙親及選擇導(dǎo)入系群體進(jìn)行高通量測(cè)序開(kāi)發(fā)SNP標(biāo)記。

    SSR及SNP標(biāo)記條帶統(tǒng)計(jì)以親本為依據(jù), 將與輪回親本K55相同的帶型記為“A”, 與供體親本K58相同的帶型記做“B”, 雜合型條帶記為“H”, 模糊不清或者缺失的帶型記為“–”。使用引物名稱(chēng)命名SSR標(biāo)記, 同一引物擴(kuò)增出多條多態(tài)性位點(diǎn)的情況以“–1、–2、–3”等區(qū)分, 順序?yàn)槠斡纱蟮叫? 使用M+數(shù)字命名SNP標(biāo)記。

    1.5 QTL定位

    根據(jù)遺傳搭車(chē)?yán)碚揫20]及基因型數(shù)據(jù), 將選擇群體的供體等位基因?qū)腩l率與BC3代的預(yù)期值(0.0625)的偏離情況進(jìn)行卡方檢測(cè), 顯著水平為<0.005, 當(dāng)某一片段可能攜帶有與目標(biāo)性狀有關(guān)的QTL時(shí), 其供體片段的基因型頻率將顯著偏離其期望頻率。

    以5株平均值代表該株系性狀值, 2年間3次重復(fù)的平均值進(jìn)行數(shù)據(jù)分析和QTL定位。對(duì)抗旱選擇群體及親本的目標(biāo)性狀計(jì)算平均值及標(biāo)準(zhǔn)差, 利用兩點(diǎn)、兩種水分條件下分別測(cè)得的產(chǎn)量性狀數(shù)據(jù), 采用SAS PROC GLM 軟件進(jìn)行性狀-標(biāo)記間單向方差分析檢測(cè)QTL, 顯著水平為<0.001。在QTL定位過(guò)程中, 為了避免采用同一個(gè)固定的閾值時(shí)出現(xiàn)在一個(gè)環(huán)境下檢測(cè)到QTL在另一個(gè)環(huán)境下卻檢測(cè)不到的情況即統(tǒng)計(jì)學(xué)上第II類(lèi)錯(cuò)誤, 可將受環(huán)境影響較大的性狀的閾值進(jìn)行調(diào)整[21], 用<0.01或<0.05作為該性狀QTL定位的顯著水平, 將結(jié)果一并列出。采用地點(diǎn)+QTL+性狀命名QTL, 性狀以英文縮寫(xiě)表示, 地點(diǎn)以大寫(xiě)首字母代表(呼和浩特H, 武川W)。同一性狀的多個(gè)QTL以“–1、–2、–3”等區(qū)分。

    2 結(jié)果與分析

    2.1 選擇導(dǎo)入系群體及親本的性狀表現(xiàn)

    從輪回親本及供體親本兩點(diǎn)的產(chǎn)量表現(xiàn)(表1)看出, 干旱脅迫條件下兩點(diǎn)輪回親本單株產(chǎn)量降低8.61 g (30.16%)、9.61 g (33.60%), 供體親本單株產(chǎn)量降低6.30 g (19.02%)、7.02 g (21.89%)。與正常供水條件相比, 干旱脅迫下兩親本的其余各性狀值在兩點(diǎn)間亦表現(xiàn)降低, 且不同性狀對(duì)水分脅迫響應(yīng)有所差異。兩點(diǎn)間兩種水分條件下輪回親本百粒重均優(yōu)于供體親本, 盤(pán)徑略?xún)?yōu)于供體親本, 其他性狀則供體親本優(yōu)于輪回親本。

    選擇群體各產(chǎn)量性狀在兩點(diǎn)間均表現(xiàn)為正常供水條件下優(yōu)于干旱條件。與輪回親本相比, 選擇群體干旱脅迫條件下除呼和浩特的盤(pán)徑及兩點(diǎn)間百粒重低于輪回親本外, 其余各性狀均優(yōu)于輪回親本, 其中2點(diǎn)間單株產(chǎn)量高于輪回親本1.65 g (8.27%)、0.70 g (3.69%), 結(jié)實(shí)率高于輪回親本4.70%、10.00%, 單株實(shí)粒數(shù)高于輪回親本85.7個(gè)(18.12%)、62.8個(gè)(13.39%); 正常供水條件下除武川的盤(pán)徑及呼和浩特的百粒重低于輪回親本外, 其余各性狀均優(yōu)于輪回親本, 其中兩點(diǎn)間單株產(chǎn)量高于輪回親本2.74 g (9.60%)、9.78 g (34.20%), 結(jié)實(shí)率高于輪回親本6.16%、10.11%, 單株實(shí)粒數(shù)高于輪回親本100.8個(gè)(17.11%)、210.0個(gè)(36.78%)。說(shuō)明回交育種結(jié)合外部選擇的策略能夠使輪回親本優(yōu)良性狀得到保持的同時(shí)彌補(bǔ)了輪回親本的不足, 使輪回親本得到改良。

    2.2 卡方檢驗(yàn)

    通過(guò)178個(gè)SSR標(biāo)記和2433個(gè)高質(zhì)量SNP標(biāo)記分析得到的供體基因?qū)腩l率與BC3代預(yù)期值(0.0625)進(jìn)行卡方檢驗(yàn), 檢測(cè)到0.005水平下的極顯著位點(diǎn)274個(gè)(表2), 卡方值范圍8.01~70.86。這些位點(diǎn)可能與抗旱性密切相關(guān)。

    表1 選擇群體及親本的性狀表現(xiàn)

    選擇群體大小為45個(gè)株系。P1: 輪回親本; P2: 供體親本; 環(huán)境E1和E2為2個(gè)試驗(yàn)點(diǎn), 分別在呼和浩特和武川。群體的性狀表型值為2年3次重復(fù)的平均值±標(biāo)準(zhǔn)差。

    Selected backcrossing introgression lines including 45 lines. P1: recurrentparent; P2: donor parent; E1 and E2 are two sites Hohhot and Wuchuan. The phenotypic values of selected backcrossing introgression lines are mean of three replicates in two years ± standard deviation.

    表2 選擇群體卡方檢驗(yàn)結(jié)果

    (續(xù)表2)

    2.3 產(chǎn)量及其相關(guān)性狀的QTL定位

    在呼和浩特市和武川縣兩點(diǎn)的干旱脅迫及正常供水條件下, 通過(guò)性狀-標(biāo)記間單項(xiàng)方差分析檢測(cè)到154個(gè)QTL (表3和表4) , 可分為三類(lèi), 第一類(lèi)為兩種水分條件下穩(wěn)定表達(dá)的QTL; 第二類(lèi)為受干旱脅迫誘導(dǎo)表達(dá)的QTL; 第三類(lèi)為僅在正常供水條件下表達(dá)的QTL。

    檢測(cè)到影響單株產(chǎn)量的QTL 44個(gè): 第一類(lèi)QTL僅在呼和浩特檢被測(cè)到2個(gè)(和), 增加目標(biāo)性狀的等位基因均來(lái)自供體親本K58。第二類(lèi)QTL在呼和浩特被檢測(cè)到7個(gè)(), 武川被檢測(cè)到12個(gè)(), 其中、在呼和浩特被同時(shí)檢測(cè)到(和), 增加目標(biāo)性狀的等位基因均來(lái)自供體親本K58。第三類(lèi)QTL在呼和浩特被檢測(cè)到15個(gè)(), 武川被檢測(cè)到10個(gè)(), 除和由K55提供等位基因外, 其余增加目標(biāo)性狀的等位基因均來(lái)自供體親本K58。

    檢測(cè)到影響結(jié)實(shí)率的QTL 28個(gè):兩點(diǎn)均未檢測(cè)到第一類(lèi)QTL。第二類(lèi)QTL僅在呼和浩特被檢測(cè)到1個(gè)(), 增加目標(biāo)性狀的等位基因來(lái)自供體親本K58。第三類(lèi)QTL在呼和浩特被檢測(cè)到5個(gè)(), 均由K55提供等位基因; 在武川被檢測(cè)到22個(gè)(), 其中共6個(gè)QTL增加目標(biāo)性狀的等位基因來(lái)自供體親本K58, 其余QTL增加目標(biāo)性狀的等位基因均來(lái)自輪回親本K55。

    檢測(cè)到影響百粒重的QTL 30個(gè): 第一類(lèi)QTL僅在武川被檢測(cè)到4個(gè)(), 增加目標(biāo)性狀的等位基因均來(lái)自輪回親本K55。第二類(lèi)QTL僅在呼和浩特被檢測(cè)到1個(gè)(), 增加目標(biāo)性狀的等位基因來(lái)自供體親本K58。第三類(lèi)QTL在呼和浩特被檢測(cè)到8個(gè)(), 除由K55提供等位基因外, 其余7個(gè)QTL增加目標(biāo)性狀的等位基因均來(lái)自供體親本K58; 在武川被檢測(cè)到17個(gè)(), 共有5個(gè)QTL (、)增加目標(biāo)性狀的等位基因來(lái)自供體親本K58, 其余QTL增加目標(biāo)性狀的等位基因均由輪回親本K55提供。

    檢測(cè)到影響單株實(shí)粒數(shù)的QTL 17個(gè), 第一類(lèi)QTL在呼和浩特被檢測(cè)到3個(gè)(), 而且這3個(gè)QTL在武川干旱條件下被檢測(cè)到()。第二類(lèi)QTL在呼和浩特被檢測(cè)到13個(gè)(、), 在武川被檢測(cè)到11個(gè)(~), 其中在呼和浩特被同時(shí)檢測(cè)到()。第三類(lèi)QTL在呼和浩特被檢測(cè)到1個(gè)(), 其在武川正常供水條件下被同時(shí)檢測(cè)到(), 兩點(diǎn)、兩種水分條件下檢測(cè)到的單株實(shí)粒數(shù)QTL增加目標(biāo)性狀的等位基因均來(lái)自供體親本K58。

    到影響盤(pán)徑的QTL 35個(gè), 兩點(diǎn)均未檢測(cè)到第一類(lèi)QTL。第二類(lèi)QTL在呼和浩特被檢測(cè)到8個(gè)(), 武川被檢測(cè)到4個(gè)(~), 增加目標(biāo)性狀的等位基因均來(lái)自輪回親本K55。第三類(lèi)QTL在呼和浩特被檢測(cè)到9個(gè)(), 增加目標(biāo)性狀的等位基因均來(lái)自供體親本K58; 在武川被檢測(cè)到14個(gè)(~), 其中、、由K55提供等位基因, 其余QTL增加目標(biāo)性狀的等位基因來(lái)自供體親本K58。

    2.4 QTL的“一因多效”

    方差分析檢測(cè)到“一因多效”位點(diǎn)14個(gè), ORS883-1、M26331、M65672和M48345處同時(shí)檢測(cè)到影響單株產(chǎn)量和百粒重的QTL; M40130、M45907、M45908、M45909、M45910、M159248、M159250、M53398、M53399、M53340處同時(shí)定位到影響單株產(chǎn)量和單株實(shí)粒數(shù)的QTL。

    2.5 兩種方法定位一致性位點(diǎn)

    方差分析檢測(cè)到的QTL中, ORS883-1、M70648、M38509、ORS525-1、M67939、M67940、M91206、M24719、M158109、ORS710-2、M48359、M35489、M73285、M108409共14個(gè)位點(diǎn)與卡方檢驗(yàn)位點(diǎn)一致, 這些可能是控制向日葵抗旱性的關(guān)鍵QTL。

    3 討論

    3.1 抗旱QTL分析

    多數(shù)抗(耐)旱QTL定位分析認(rèn)為, 水、旱環(huán)境下表現(xiàn)相似的QTL可能對(duì)抗旱性有直接貢獻(xiàn), 而受干旱誘導(dǎo)的QTL可能與耐旱適應(yīng)性反應(yīng)有關(guān), 未必對(duì)耐旱性有實(shí)質(zhì)性貢獻(xiàn)[3,22-23]。本研究中檢測(cè)到的QTL受干旱脅迫表達(dá)水平不同,方差分析檢測(cè)到的154個(gè)QTL中, 102個(gè)在單一試驗(yàn)點(diǎn)的正常供水條件下被檢測(cè)到, 占比66.23%; 57個(gè)在單一試驗(yàn)點(diǎn)的干旱脅迫條件下被檢測(cè)到, 占比37.01%, 其中同一水分條件下2個(gè)試驗(yàn)點(diǎn)重復(fù)檢測(cè)到影響單株產(chǎn)量、單株實(shí)粒數(shù)的QTL共10個(gè), 占比6.49%, 這些QTL受環(huán)境誘導(dǎo)在不同點(diǎn)穩(wěn)定表達(dá), 其與抗旱性關(guān)系有待研究。同一試驗(yàn)點(diǎn)兩種水分條件下檢測(cè)到影響單株產(chǎn)量、百粒重及單株實(shí)粒數(shù)的9個(gè)QTL, 占比5.84%, 這些QTL可能對(duì)向日葵抗旱性有直接貢獻(xiàn), 其中影響實(shí)粒數(shù)的3個(gè)QTL在兩點(diǎn)、多種水分條件下均被檢測(cè)到, 且加性效應(yīng)方向一致, 可靠性較高, 可能是向日葵抗旱的重要QTL, 在未來(lái)的工作中可以進(jìn)一步驗(yàn)證, 對(duì)向日葵抗旱性分子育種研究具有重要意義。

    3.2 統(tǒng)計(jì)分析方法比較

    本文利用方差分析和卡方檢驗(yàn)兩種方法對(duì)向日葵定向選擇群體進(jìn)行QTL定位。經(jīng)干旱脅迫, 通過(guò)單株產(chǎn)量篩選獲得了明顯優(yōu)于輪回親本的抗旱回交選擇導(dǎo)入系。選擇群體經(jīng)連續(xù)定向選擇, 理論上應(yīng)保留有大部分的耐旱相關(guān)基因, 與抗(耐)旱性有關(guān)的等位基因?qū)腩l率偏離嚴(yán)重, 基于遺傳搭車(chē)原理的卡方檢驗(yàn)使抗旱相關(guān)位點(diǎn)的定位得以實(shí)現(xiàn), 方差分析則能夠直接定位與性狀表型數(shù)據(jù)對(duì)應(yīng)的位點(diǎn), 兩者結(jié)合可提高定位結(jié)果的可靠性[16]。本試驗(yàn)用卡方檢驗(yàn)檢測(cè)到274個(gè)極顯著偏分離位點(diǎn), 其中有14個(gè)與方差分析檢測(cè)到的QTL一致, 這些位點(diǎn)是向日葵抗旱性相關(guān)的關(guān)鍵位點(diǎn)。兩種方法檢測(cè)到的一致性位點(diǎn)較少, 與鄭天清等人的研究結(jié)果類(lèi)似, 一些在方差分析中較穩(wěn)定表達(dá)的QTL在導(dǎo)入位點(diǎn)上未表現(xiàn)為等位基因位點(diǎn)的顯著偏分離, 而一些極顯著偏分離位點(diǎn)則未在方差分析中定位到相關(guān)性狀的QTL, 其研究表明這些定位到的位點(diǎn)無(wú)法簡(jiǎn)單地以表型選擇或改變?nèi)后w的基因頻率來(lái)加以利用[24]。張金巍的研究也發(fā)現(xiàn)一些在單向方差分析中檢測(cè)到的解釋表型貢獻(xiàn)率較大的位點(diǎn)并未在卡方測(cè)驗(yàn)法中被檢測(cè)到, 其推測(cè)可能是由于選擇的力度抵消了位點(diǎn)原本較高的效應(yīng)[18]。

    3.3 SNP分子標(biāo)記的應(yīng)用

    本試驗(yàn)的研究結(jié)果表明, 利用高代回交結(jié)合抗旱篩選得到的選擇群體進(jìn)行向日葵抗旱QTL的發(fā)掘定位是可行且有效的。盡可能多地鑒定出影響目標(biāo)性狀的供體導(dǎo)入片段是采用導(dǎo)入系定位QTL的目的與前提, 采用更飽和的分子標(biāo)記可以防止小的供體導(dǎo)入片段被漏檢, 從而鑒定出盡可能多的供體導(dǎo)入片段[25]。SNP標(biāo)記因其具有高通量、遺傳穩(wěn)定性好、密度高等優(yōu)點(diǎn), 在QTL研究中應(yīng)用廣泛。本試驗(yàn)首次將SNP標(biāo)記結(jié)合SSR標(biāo)記對(duì)選擇群體基因組進(jìn)行掃描, 比僅用SSR標(biāo)記掃描得到了更多的位點(diǎn)(2611個(gè)), 在此基礎(chǔ)上進(jìn)行QTL檢測(cè), 卡方檢驗(yàn)得到274個(gè)偏分離位點(diǎn), 占總位點(diǎn)數(shù)的10.49%; 方差分析每個(gè)性狀在兩點(diǎn)、兩種水分條件下檢測(cè)到的QTL數(shù)目為17~44個(gè), 平均為30.8個(gè); 每個(gè)性狀在單個(gè)試驗(yàn)點(diǎn)、單種水分條件下單獨(dú)統(tǒng)計(jì)檢測(cè)到的QTL數(shù)目為0~22個(gè), 平均為8.85個(gè)。與以往利用SSR標(biāo)記檢測(cè)QTL相比, 本試驗(yàn)檢測(cè)到了更多的QTL, 提高了QTL檢測(cè)的效率, 未來(lái)定位到圖譜上則體現(xiàn)為區(qū)間更小, 定位更加精準(zhǔn), 同時(shí)在目標(biāo)QTL區(qū)域進(jìn)一步開(kāi)發(fā)SNP標(biāo)記, 可實(shí)現(xiàn)QTL的精細(xì)定位。

    3.4 QTL應(yīng)用探討

    選擇導(dǎo)入系是經(jīng)過(guò)目標(biāo)性狀定向選擇的群體, 高強(qiáng)度的選擇導(dǎo)致群體的規(guī)模較小、等位基因頻率偏離嚴(yán)重, 因而無(wú)法實(shí)施傳統(tǒng)的連鎖作圖分析。借助已構(gòu)建完成的物理圖譜可以明確QTL在圖譜上的相對(duì)位置[26-27], 但向日葵的參考圖譜目前尚未明確[28],因而本試驗(yàn)中僅檢測(cè)到與QTL緊密連鎖的標(biāo)記而無(wú)法明確其在連鎖群的位置, 也無(wú)法排除檢測(cè)到的QTL存在重復(fù)檢測(cè)的可能, 但是仍具有參考價(jià)值。隨著向日葵圖譜的發(fā)展, 借助參考圖譜可明確本研究定位到的QTL在連鎖群的位置, 并進(jìn)行QTL定位結(jié)果間的比較分析, 從而應(yīng)用于分子標(biāo)記輔助育種。

    4 結(jié)論

    利用45個(gè)選擇導(dǎo)入系方差分析檢測(cè)到154個(gè)QTL, 其中2種水分條件下穩(wěn)定表達(dá)的9個(gè), 包括武川影響百粒重的4個(gè)()、呼和浩特影響單株產(chǎn)量的2個(gè)(和)及影響單株實(shí)粒數(shù)的3個(gè)(), 這些QTL可能對(duì)抗旱性有直接貢獻(xiàn)??ǚ綑z驗(yàn)檢測(cè)到極顯著偏分離位點(diǎn)274個(gè), 其中14個(gè)與方差分析的結(jié)果一致, 可能是與抗旱性相關(guān)的關(guān)鍵位點(diǎn)。

    [1] 王韻. 利用雙向?qū)胂到馕鏊井a(chǎn)量相關(guān)性狀的遺傳背景效應(yīng)及環(huán)境互作效應(yīng). 沈陽(yáng)農(nóng)業(yè)大學(xué)博士學(xué)位論文, 遼寧沈陽(yáng), 2009 Wang Y. Genetic Background Effect and Environment Interaction Effect on QTL Expression of Yield Related Traits Revealed by Reciprocal Introgression Lines in Rice. PhD Dissertation of Shenyang Agricultural University, Shenyang, China, 2009 (in Chinese with English abstract)

    [2] 王陽(yáng), 劉成, 王天宇, 石云素, 宋燕春, 黎裕. 干旱脅迫和正常灌溉條件下玉米產(chǎn)量性狀的QTL分析. 植物遺傳資源學(xué)報(bào), 2007, 8: 179–183Wang Y, Liu C, Wang T Y, Shi Y S, Song Y C, Li Y. QTL analysis of yield components in maize under different water regimes., 2007, 8: 179–183 (in Chinese with English abstract)

    [3] 趙秀琴, 徐建龍, 朱苓華, 黎志康. 利用高代回交導(dǎo)入系定位水、旱條件下影響水稻根系及產(chǎn)量的QTL. 中國(guó)農(nóng)業(yè)科學(xué), 2008, 41: 1887–1893 Zhao X Q, Xu J L, Zhu L H, Li Z K. QTL Mapping of yield and root traits under irrigation and drought conditions using advanced backcrossing introgression lines in rice., 2008, 41: 1887–1893 (in Chinese with English abstract)

    [4] Adiredjo A L, Navaud O, Munos S, Langlade N B, Lamaze T, Grieu P. Genetic control of water use efficiency and leaf carbon isotope discrimination in sunflower (L.) subjected to two drought scenarios., 2014, 9: e101218

    [5] Abdi1 N, Darvishzadeh R, Jafaril M, Pirzad A, Haddadi P. Genetic analysis and QTL mapping of agro-morphological traits in sunflower (L.) under two contrasting water treatment conditions., 2012, 5: 149–158

    [6] Haddadi P, Yazdi-samadi B, Naghavi M R, Kalantari A, Maury P, Sarrafi A. QTL analysis of agronomic traits in recombinant inbred lines of sunflower under partial irrigation., 2011, 5: 135–146

    [7] 呂品, 于海峰, 于志賢, 張永虎, 張艷芳, 王婷婷, 侯建華. 向日葵高密度遺傳連鎖圖譜構(gòu)建及兩種水分條件下芽期性狀的QTL分析. 作物學(xué)報(bào), 2017, 43: 19–30 Lyu P, Yu H F, Yu Z X, Zhang Y H, Zhang Y F, Wang T T, Hu J H. Construction of high-density genetic map and qtl mapping for seed germination traits in sunflower under two water conditions., 2017, 43: 19–30 (in Chinese with English abstract)

    [8] Li Z K, Fu Z K, Gao Y M, Xu J L, Ali J, Lafittle H R, Jiang Y Z, Domingo R J, Vijayakumar C H M, Maghirang R, Zheng T Q, Zhu L H. Genome–wide introgression lines and their use in genetic and molecular dissection of complex phenotypes in rice (L. )., 2005, 59: 33–52

    [9] 劉海燕. 水稻高代回交導(dǎo)入系抗旱相關(guān)性狀QTL定位研究. 中國(guó)農(nóng)業(yè)科學(xué)院碩士學(xué)位論文, 北京, 2009 Liu H Y. QTL Mapping of Traits Related to Drought Tolerance of Advanced Backcrossing Introgression Lines in Rice. MS Thesis of Chinese Academy of Agricultural Sciences, Beijing, China, 2009 (in Chinese with English abstract)

    [10] 黎志康. 我國(guó)水稻分子育種計(jì)劃的策略. 分子植物育種, 2005, 3: 603–608 Li Z K. Strategies for molecular rice breeding in China., 2005, 3: 603–608

    [11] 卓壯. 黃華占回交導(dǎo)入系抗旱篩選及QTL分析. 四川農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 四川雅安, 2013 Zhuo Z. Screening and QTL Mapping of Advanced Backcrossing Introgression Lines for Drought Tolerance in Huang Huazhan. MS Thesis of Sichuan Agricultural University, Sichuan Ya’an, China, 2013 (in Chinese with English abstract)

    [12] 趙秀琴, 朱苓華, 徐建龍, 黎志康. 灌溉與自然降雨條件下水稻高代回交導(dǎo)入系產(chǎn)量QTL的定位. 作物學(xué)報(bào), 2007, 33: 1536–1542 Zhao X Q, Zhu L H , Xu J L, Li Z K. QTL mapping of yield under irrigation and rainfed field conditions for advanced backcrossing introgression lines in rice.. 2007, 33: 1536–1542 (in Chinese with English abstract)

    [13] 王利鋒, 郝憲彬, 何云霞, 華澤田, 高用明, 黎志康. 利用選擇回交導(dǎo)入群體進(jìn)行粳稻優(yōu)良恢復(fù)系抗旱性改良與QTL定位.安徽農(nóng)業(yè)大學(xué)學(xué)報(bào), 2007, 34: 311–319Wang L F, Hao X B, He Y X, Hua Z T, Gao Y M, Li Z K. Improvement of drought tolerance for elite restorer ofrice (L.) and identification of drought tolerant QTLs using selected introgression populations by backcross.2007, 34: 311–319 (in Chinese with English abstract)

    [14] 王英. 利用回交導(dǎo)入系篩選水稻高產(chǎn)、抗旱和耐鹽株系及選擇導(dǎo)入系相關(guān)性狀的QTL定位. 中國(guó)農(nóng)業(yè)科學(xué)院博士學(xué)位論文, 北京, 2013 Wang Y. Screening of High Yield, Drought and Salt Tolerant Plants from Backcross Introgression Lines and QTL Detection for Related Traits in Rice. PhD Dissertation of Chinese Academy of Agricultural Sciences, Beijing, China, 2013 (in Chinese with English abstract)

    [15] 任潔, 趙秀琴, 丁在松, 項(xiàng)超, 張晶, 王超, 張俊巍, Joseph C A, 張強(qiáng), 龐昀龍, 高用明, 石英堯. 利用選擇導(dǎo)入系進(jìn)行水稻耐低磷鑒定與QTL定位分析.中國(guó)水稻科學(xué), 2015, 29: 1–13Ren J, Zhao X Q , Ding Z S, Xiang C, Zhang J, Wang C, Zhang J W, Joseph C A, Zhang Q, PongY L, Gao Y M, Shi Y Y. Dissection and QTL mapping of low-phosphorus tolerance using selected introgression lines in rice.,2015, 29: 1–13 (in Chinese with English abstract)

    [16] 李芳, 程立銳, 許美容, 周政, 張帆, 孫勇, 周永力, 朱苓華, 徐建龍, 黎志康. 利用品質(zhì)性狀的回交選擇導(dǎo)入系挖掘水稻抗紋枯病 QTL. 作物學(xué)報(bào), 2009, 35: 1729?1737 Li F, Cheng L R, Xu M R, Zhou Z, Zhang F, Sun Y, Zhou Y L, Zhu L H, Xu J L, Li Z K. QTL mining for sheath blight resistance using the backcross selected introgression lines for grain quality in rice., 2009, 35: 1729–1737 (in Chinese with English abstract)

    [17] 李燦東, 蔣洪蔚, 劉春燕, 邱鵬程, 張聞博, 李文福, 高運(yùn)來(lái), 陳慶山, 胡國(guó)華. 大豆定向選擇群體耐旱性位點(diǎn)基因型分析及QTL定位. 中國(guó)油料作物學(xué)報(bào), 2009, 31: 285–292 Li C D, Jiang H W, Liu C Y, Qiu P C, Zhang W B, Li T F, Gao Y L, Chen Q S. Genotype and QTL analysis of drought tolerance loci for directional population in soybean.2009, 31: 285–292 (in Chinese with English abstract)

    [18] 張金巍, 韓粉霞, 陳明陽(yáng), 孫君明, 韓廣振, 閆淑榮, 楊華. 利用選擇回交導(dǎo)入群體定位大豆蛋白質(zhì)含量 QTL. 中國(guó)油料作物學(xué)報(bào), 2015, 37: 433–442 Zhang J W, Han F X, Chen M Y, Sun J M, Han G Z, Yan S R, Yang H. Detection of protein content QTL with random and extremely selected BC2F2populations in soybean., 2015, 37: 433–442 (in Chinese with English abstract)

    [19] 房冬梅, 呂品, 侯建華. 油葵SSR-PCR反應(yīng)體系的優(yōu)化及引物篩選. 中國(guó)農(nóng)學(xué)通報(bào), 2015, 31(12): 205–209 Fang D M, Lyu P, Hou J H. Optimization of SSR-PCR reaction system and primer screening in oil sunflower., 2015, 31(12): 205–209 (in Chinese with English abstract)

    [20] Schlotterer C. Hitchhiking mapping-functional genomics from the population genetics perspective., 2003, 19: 32–38

    [21] 邱鵬程, 張聞博, 李燦東, 蔣洪蔚, 劉春燕, 范冬梅, 曾慶力, 胡國(guó)華, 陳慶山. 利用選擇導(dǎo)入系分析大豆芽期和苗期耐旱性的遺傳重疊. 作物學(xué)報(bào), 2011, 37: 477–483 Qiu P C, Zhang W B, Li C D, Jiang H W, Liu C Y, Fan D M, Zeng Q L, Hu G H, Chen Q S. Genetic overlap of drought- tolerance loci between germination stage and seedling stage analyzed using introgression lines in soybean.2011, 37: 477–483 (in Chinese with English abstract)

    [22] 彭勃, 王陽(yáng), 李永祥, 劉成, 劉志齋, 王迪, 譚巍巍, 張巖, 孫寶成, 石云素, 宋燕春, 王天宇, 黎裕. 不同水分環(huán)境下玉米產(chǎn)量構(gòu)成因子及籽粒相關(guān)性狀的QTL分析. 作物學(xué)報(bào), 2010, 36: 1832–1842 Peng B, Wang Y, Li Y X, Liu C, Liu Z Z, Wang D, Tan W W, Zhang Y, Sun B C, Shi Y S, Song Y C, Wang T Y, Li Y. QTL analysis for yield components and kernel-related traits in maize under different water regimes., 2010, 36: 1832–1842 (in Chinese with English abstract)

    [23] 譚巍巍, 李永祥, 王陽(yáng), 劉成, 劉志齋, 彭勃, 王迪, 張巖, 孫寶成, 石云素, 宋燕春, 楊德光, 王天宇, 黎裕. 在干旱和正常水分條件下玉米穗部性狀QTL分析. 作物學(xué)報(bào), 2011, 37: 235–248Tan W W, Li Y X, Wang Y, Liu C, Liu Z Z, Peng B, Wang D, Zhang Y, Sun B C, Shi Y S, Song Y C, Yang D G, Wang T Y, Li Y. QTL mapping of ear traits of maize under different water regimes., 2011, 37: 235–248 (in Chinese with English abstract)

    [24] 鄭天清, 徐建龍, 傅彬英, 高用明, Veruka S, Lafitte R, 翟虎渠, 萬(wàn)建民, 朱苓華, 黎志康. 遺傳搭車(chē)與方差分析在水稻定向選擇群體的抗旱性位點(diǎn)分析中的初步應(yīng)用. 作物學(xué)報(bào), 2007, 33: 799–804 Zheng T Q, Xu J L, Fu B Y, Gao Y M, Veruka S, Lafitte R, Zhai H Q, Wan J M, Zhu L H, Li Z K. application of genetic hitch-hiking and ANOVA in identification of loci for drought tolerance in populations of rice from directional selection., 2007, 33: 779–804 (in Chinese with English abstract)

    [25] 何云霞. 應(yīng)用回交育種和QTL聚合改良粳稻抗旱性. 沈陽(yáng)農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 遼寧沈陽(yáng), 2009 He Y X. Improving Drought Tolerance ofRice by BC Breeding and QTL Pyramiding. MS Thesis of Shenyang Agricultural University, Shenyang, China, 2009 (in Chinese with English abstract)

    [26] 石英堯. 水稻選擇導(dǎo)入系產(chǎn)量和抗倒伏性狀改良及QTL定位.安徽農(nóng)業(yè)大學(xué)碩士學(xué)位論文, 安徽合肥, 2013 Shi Y R. Improvement and QTL Mapping for Yield and Loding Resistance Using Selected Backcrossing Introgression Lines in Rice. MS Thesis of Anhui Agricultural University, Hefei, China, 2013 (in Chinese with English abstract)

    [27] 李燦東, 蔣洪蔚, 劉春燕, 郭泰, 王志新, 吳秀紅, 鄭偉, 邱鵬程, 張聞博, 宋英博, 欒奕娜, 陳慶山, 胡國(guó)華. 大豆耐旱選擇群體QTL定位. 作物學(xué)報(bào), 2011, 37: 603–611 Li C D, Jiang H W, Liu C Y, Guo T, Wang Z X, Wu X H, Zheng W, Qiu P C, Zhang W B, Song Y B, Luan Y N, Chen Q S, Hu G H. QTL identification of drought tolerance to soybean in selection population.,2011, 37: 603–611 (in Chinese with English abstract)

    [28] Sharma R, Xia X, Cano L M, Evangelisti E, Kemen E, Judelson H, Oome S, Sambles C, van den Hoogen D J, Kitner M, Klein J, Meijer H J G, Spring O, Win J, Zipper R, Bode H B, Govers F, Kamoun S, Schornack S, Studholme D J, Van den Ackerveken G, Thines M. Genome analyses of the sunflower pathogenprovide insights into effector evolution in downy mildews and phytophthora., 2015, 16: 2–23

    QTL Mapping of Yield Traits Using Drought Tolerance Selected Backcrossing Introgression Lines in Sunflower

    LYU Pin1, YU Hai-Feng2, and HOU Jian-Hua1,*

    1College of Agriculture, Inner Mongolia Agricultural University;2Institute of Crop Research, Inner Mongolia Academy of Agricultural and Animal Husbandry Sciences, Hohhot 010010, Inner Mongolia, China

    Drought is one of the most important factors to decrease the yield of sunflower. The BC3F2selected backcrossing introgression lines of oil sunflower population including 45 lines were developed by drought-tolerance screening for yield usingK55 with excellent comprehensive characters and drought sensitivity as recurrent parent and K58 with drought-tolerance as donor parent. The genotypes of selected backcrossing introgression lines were obtained with the whole genome SSR and SNP markers.QTLs affecting five yield traits were detected under both drought stress and well watered conditions in Hohhot and Wuchuan respectively by one-way ANOVA and Chi-square test based on the Genetic Hitchhiking Effect.The QTLs detected by one-way ANOVA were grouped into three types based on their behaviors: type I was the QTLs detected in both watered conditions, including four QTLs affecting hundred-seed weight (HSW) in Wuchuan and two QTLs affecting seed yield (SY) and three QTLs affecting filled seeds per plant (FSP) in Hohhot, which were considered to be able to directly contribute to drought tolerance; type II was the QTLs detected only under drought stress, including 30 QTLs in Hohhot and 27 QTLs in Wuchuan; type III was the QTLs detected only in well watered condition, including 38 QTLs in Hohhot and 64 QTLs in Wuchuan. There 274 loci were detected by chi-square test, among them 14 loci could be detected by ANOVA and chi-square test simultaneously, which might be the key loci for drought tolerance of sunflower. The results lay a foundation of efficient drought tolerance molecular breeding and provide useful drought tolerance materials for sunflower.

    sunflower; droughttolerance; selected backcrossing introgression lines; yield traits; QTL

    2017-05-15;

    2017-11-21;

    2017-12-25.

    10.3724/SP.J.1006.2018.00385

    本研究由國(guó)家自然科學(xué)基金項(xiàng)目(31160288)資助。

    This study was supported by the National Natural Science Foundation of China (31160288).

    Corresponding author侯建華, E-mail: houjh68@163.com

    E-mail: 18847123096@163.com

    http://kns.cnki.net/kcms/detail/11.1809.S.20171225.0843.002.html

    猜你喜歡
    供體親本向日葵
    甘蔗親本農(nóng)藝性狀評(píng)價(jià)與分析
    向日葵
    幾種蘋(píng)果砧木實(shí)生后代與親本性狀的相關(guān)性
    向日葵
    云瑞10系列生產(chǎn)性創(chuàng)新親本2種方法評(píng)價(jià)
    不同氫供體對(duì)碳納米管負(fù)載鈀催化劑催化溴苯脫溴加氫反應(yīng)的影響
    向日葵
    高齡供體的肝移植受者生存分析
    終末熱灌注對(duì)心臟移植術(shù)中豬供體心臟的保護(hù)作用
    油用向日葵親本繁殖及雜交制種技術(shù)
    国产视频一区二区在线看| 97人妻精品一区二区三区麻豆| 久久久久久久亚洲中文字幕 | 午夜免费激情av| 日本a在线网址| 天堂动漫精品| 国产乱人视频| 好男人电影高清在线观看| 国产爱豆传媒在线观看| 国产高潮美女av| 欧美区成人在线视频| 又紧又爽又黄一区二区| 女同久久另类99精品国产91| 白带黄色成豆腐渣| 99国产精品一区二区三区| 国产精品99久久久久久久久| 免费av观看视频| 国产午夜精品论理片| 制服人妻中文乱码| 性欧美人与动物交配| x7x7x7水蜜桃| 夜夜夜夜夜久久久久| 精品熟女少妇八av免费久了| 午夜激情福利司机影院| 久99久视频精品免费| 久99久视频精品免费| 国产真实乱freesex| 美女cb高潮喷水在线观看| 日韩高清综合在线| 美女cb高潮喷水在线观看| 欧美日韩一级在线毛片| 午夜激情欧美在线| 国产淫片久久久久久久久 | 少妇丰满av| eeuss影院久久| 五月伊人婷婷丁香| 啦啦啦免费观看视频1| 亚洲成人中文字幕在线播放| 村上凉子中文字幕在线| 神马国产精品三级电影在线观看| 最好的美女福利视频网| 亚洲成人精品中文字幕电影| 18禁黄网站禁片免费观看直播| 国内揄拍国产精品人妻在线| 久久精品夜夜夜夜夜久久蜜豆| 他把我摸到了高潮在线观看| av欧美777| 精品乱码久久久久久99久播| 免费高清视频大片| 国产高清视频在线播放一区| 国产在视频线在精品| 免费人成视频x8x8入口观看| 亚洲国产精品sss在线观看| 啦啦啦观看免费观看视频高清| 久久久色成人| 一二三四社区在线视频社区8| 1000部很黄的大片| 丁香六月欧美| 国产v大片淫在线免费观看| 搡老妇女老女人老熟妇| 国产精品一区二区三区四区久久| 人人妻,人人澡人人爽秒播| 看黄色毛片网站| x7x7x7水蜜桃| 国内毛片毛片毛片毛片毛片| 亚洲专区中文字幕在线| 在线观看免费午夜福利视频| 国产精品美女特级片免费视频播放器| 黄片大片在线免费观看| 久久午夜亚洲精品久久| 国产高清三级在线| 亚洲精品影视一区二区三区av| 岛国在线免费视频观看| 亚洲精华国产精华精| 最近最新免费中文字幕在线| 亚洲精品久久国产高清桃花| 亚洲五月婷婷丁香| 成年人黄色毛片网站| 欧美成狂野欧美在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲欧美日韩高清专用| 很黄的视频免费| 成人鲁丝片一二三区免费| 又爽又黄无遮挡网站| 亚洲男人的天堂狠狠| www.熟女人妻精品国产| 国产精品电影一区二区三区| a级一级毛片免费在线观看| 人人妻人人看人人澡| 欧美日韩瑟瑟在线播放| 女同久久另类99精品国产91| 动漫黄色视频在线观看| 成人av一区二区三区在线看| 国内少妇人妻偷人精品xxx网站| 亚洲美女黄片视频| 少妇丰满av| 深夜精品福利| 宅男免费午夜| 国产一区二区激情短视频| 最近视频中文字幕2019在线8| 香蕉丝袜av| 国产精品 国内视频| 久久久久国内视频| 欧美成狂野欧美在线观看| 一进一出抽搐动态| x7x7x7水蜜桃| 成熟少妇高潮喷水视频| 国产aⅴ精品一区二区三区波| 免费在线观看日本一区| 可以在线观看毛片的网站| 搞女人的毛片| 小说图片视频综合网站| 999久久久精品免费观看国产| 成人国产一区最新在线观看| 欧美av亚洲av综合av国产av| 草草在线视频免费看| 一二三四社区在线视频社区8| 69人妻影院| 亚洲五月天丁香| 亚洲黑人精品在线| 成年人黄色毛片网站| 亚洲成a人片在线一区二区| 欧美国产日韩亚洲一区| 精品国内亚洲2022精品成人| 国产亚洲精品av在线| 两性午夜刺激爽爽歪歪视频在线观看| 小说图片视频综合网站| 精品电影一区二区在线| 久久国产乱子伦精品免费另类| 久久久久久久久大av| 国产高清视频在线观看网站| 综合色av麻豆| 成人午夜高清在线视频| 久久久久九九精品影院| 美女 人体艺术 gogo| 欧美另类亚洲清纯唯美| 午夜精品一区二区三区免费看| 国产亚洲精品一区二区www| 日韩精品中文字幕看吧| 小说图片视频综合网站| 人人妻,人人澡人人爽秒播| 黄片小视频在线播放| 天堂动漫精品| 99国产综合亚洲精品| 最后的刺客免费高清国语| 非洲黑人性xxxx精品又粗又长| 婷婷丁香在线五月| 最新在线观看一区二区三区| 天堂影院成人在线观看| 亚洲在线观看片| 免费看光身美女| 亚洲人成网站在线播放欧美日韩| 12—13女人毛片做爰片一| 无遮挡黄片免费观看| 国产精品久久久久久精品电影| 国产三级黄色录像| 偷拍熟女少妇极品色| 成人无遮挡网站| 99热6这里只有精品| 熟女少妇亚洲综合色aaa.| 国产在线精品亚洲第一网站| 在线观看舔阴道视频| 在线十欧美十亚洲十日本专区| 亚洲中文字幕一区二区三区有码在线看| 1000部很黄的大片| 少妇人妻一区二区三区视频| 国产精品久久久久久精品电影| 亚洲av美国av| 免费看美女性在线毛片视频| 欧美成人a在线观看| 在线观看免费午夜福利视频| 天堂影院成人在线观看| 亚洲成av人片在线播放无| 亚洲人成网站在线播放欧美日韩| 男女午夜视频在线观看| 国产精品久久久久久精品电影| 国产亚洲欧美在线一区二区| 亚洲国产精品成人综合色| 97超视频在线观看视频| 欧美黑人欧美精品刺激| 久久久久久国产a免费观看| 最近最新免费中文字幕在线| 狂野欧美白嫩少妇大欣赏| 蜜桃久久精品国产亚洲av| 免费人成在线观看视频色| 国产精品99久久99久久久不卡| 国产三级中文精品| 久久久成人免费电影| 露出奶头的视频| 久久精品91无色码中文字幕| 在线看三级毛片| 日本一本二区三区精品| 日韩欧美在线乱码| 成人精品一区二区免费| 亚洲av成人精品一区久久| 熟女少妇亚洲综合色aaa.| 国产亚洲欧美98| 亚洲av免费高清在线观看| 日本a在线网址| 亚洲国产欧美人成| 国产私拍福利视频在线观看| 日本与韩国留学比较| 12—13女人毛片做爰片一| 波野结衣二区三区在线 | 国产av麻豆久久久久久久| 午夜精品一区二区三区免费看| 九色国产91popny在线| 一本精品99久久精品77| 丰满人妻熟妇乱又伦精品不卡| 蜜桃亚洲精品一区二区三区| 最新在线观看一区二区三区| 99热6这里只有精品| 黄色片一级片一级黄色片| 欧洲精品卡2卡3卡4卡5卡区| 一本综合久久免费| 12—13女人毛片做爰片一| 久久中文看片网| 岛国在线免费视频观看| 少妇人妻精品综合一区二区 | 午夜福利18| 久久精品国产清高在天天线| 一个人免费在线观看的高清视频| 青草久久国产| 99久久成人亚洲精品观看| 成人高潮视频无遮挡免费网站| 亚洲精品在线观看二区| 亚洲国产欧洲综合997久久,| 18禁裸乳无遮挡免费网站照片| 久久久久精品国产欧美久久久| 欧美不卡视频在线免费观看| 国产免费一级a男人的天堂| 久久精品国产清高在天天线| 啪啪无遮挡十八禁网站| 免费看日本二区| 亚洲成av人片在线播放无| 波多野结衣高清无吗| 99精品欧美一区二区三区四区| 成人无遮挡网站| 校园春色视频在线观看| 国产伦一二天堂av在线观看| 波多野结衣高清无吗| 国产成人福利小说| 久久久国产成人免费| 内地一区二区视频在线| 国产免费男女视频| 午夜免费成人在线视频| 色综合站精品国产| 久久久久久久精品吃奶| 亚洲av不卡在线观看| 免费人成在线观看视频色| 日本黄色视频三级网站网址| 99久久精品热视频| 免费观看人在逋| 亚洲成a人片在线一区二区| 亚洲国产欧美网| 少妇人妻一区二区三区视频| 搡老妇女老女人老熟妇| 国产精品亚洲av一区麻豆| 亚洲av中文字字幕乱码综合| 精品一区二区三区人妻视频| 久久久成人免费电影| or卡值多少钱| 免费人成视频x8x8入口观看| 高清毛片免费观看视频网站| 欧美乱妇无乱码| 一级黄片播放器| 91麻豆av在线| 亚洲av电影不卡..在线观看| 亚洲avbb在线观看| svipshipincom国产片| 国内精品久久久久精免费| 制服人妻中文乱码| 啦啦啦韩国在线观看视频| 免费在线观看日本一区| 国产在线精品亚洲第一网站| 国内精品久久久久精免费| 欧美最黄视频在线播放免费| 中文字幕av在线有码专区| 国产亚洲精品av在线| 一个人看的www免费观看视频| 欧美三级亚洲精品| 国内少妇人妻偷人精品xxx网站| 亚洲精品在线观看二区| 午夜福利18| 黄色女人牲交| 欧美成人一区二区免费高清观看| 亚洲av不卡在线观看| 精品久久久久久久久久免费视频| 亚洲中文字幕日韩| 国产成年人精品一区二区| 久久久久九九精品影院| 亚洲在线观看片| 国内精品美女久久久久久| 十八禁网站免费在线| 午夜精品久久久久久毛片777| 精品国内亚洲2022精品成人| 国产乱人视频| 久久久久久国产a免费观看| 精品无人区乱码1区二区| 亚洲午夜理论影院| 国产精品久久久久久亚洲av鲁大| 国产国拍精品亚洲av在线观看 | 欧美日韩亚洲国产一区二区在线观看| 成人永久免费在线观看视频| 精华霜和精华液先用哪个| 久久亚洲真实| 亚洲国产中文字幕在线视频| 国产精品一区二区三区四区免费观看 | 欧美另类亚洲清纯唯美| 白带黄色成豆腐渣| 香蕉av资源在线| 欧美色视频一区免费| 国产av麻豆久久久久久久| 91av网一区二区| 18美女黄网站色大片免费观看| 精品久久久久久久久久免费视频| 我的老师免费观看完整版| 亚洲精品国产精品久久久不卡| 又爽又黄无遮挡网站| h日本视频在线播放| 欧美午夜高清在线| 国产精华一区二区三区| 国产精品电影一区二区三区| 久久精品人妻少妇| 听说在线观看完整版免费高清| 99久国产av精品| 一区二区三区激情视频| 90打野战视频偷拍视频| 亚洲精品456在线播放app | 高清毛片免费观看视频网站| 国产高清三级在线| 两个人的视频大全免费| 亚洲无线在线观看| 亚洲av不卡在线观看| 麻豆成人av在线观看| 国产av不卡久久| 成年免费大片在线观看| 91在线观看av| 97人妻精品一区二区三区麻豆| 久久精品91无色码中文字幕| 久久人妻av系列| 中文字幕av成人在线电影| 午夜亚洲福利在线播放| 亚洲国产精品久久男人天堂| 国产午夜福利久久久久久| 国产成人福利小说| 首页视频小说图片口味搜索| h日本视频在线播放| 变态另类丝袜制服| 欧美黄色淫秽网站| 网址你懂的国产日韩在线| 国内少妇人妻偷人精品xxx网站| 午夜免费成人在线视频| 在线视频色国产色| 亚洲av五月六月丁香网| 少妇的逼水好多| 亚洲av成人精品一区久久| 91麻豆av在线| 99热精品在线国产| 午夜福利高清视频| 国产不卡一卡二| 桃红色精品国产亚洲av| 成人精品一区二区免费| 一级a爱片免费观看的视频| 亚洲精品久久国产高清桃花| 国产伦精品一区二区三区四那| 最好的美女福利视频网| 国产精品久久视频播放| 窝窝影院91人妻| 美女免费视频网站| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 啪啪无遮挡十八禁网站| xxx96com| 一进一出抽搐gif免费好疼| 91在线精品国自产拍蜜月 | 免费av毛片视频| 日韩有码中文字幕| 日韩精品中文字幕看吧| 欧美黄色片欧美黄色片| av片东京热男人的天堂| av天堂中文字幕网| 美女高潮喷水抽搐中文字幕| 国产麻豆成人av免费视频| 免费无遮挡裸体视频| 床上黄色一级片| 91麻豆av在线| 国产一区二区三区视频了| 香蕉丝袜av| 国产色爽女视频免费观看| 神马国产精品三级电影在线观看| 久久久久久久精品吃奶| 看免费av毛片| 精华霜和精华液先用哪个| 免费观看的影片在线观看| 国产午夜精品论理片| 香蕉av资源在线| aaaaa片日本免费| 免费av不卡在线播放| 国产黄a三级三级三级人| 少妇的丰满在线观看| 亚洲中文字幕一区二区三区有码在线看| 香蕉丝袜av| 国产色婷婷99| 中文字幕人成人乱码亚洲影| 免费观看的影片在线观看| 国产伦在线观看视频一区| 免费电影在线观看免费观看| 最新在线观看一区二区三区| 日本成人三级电影网站| 成人特级av手机在线观看| 尤物成人国产欧美一区二区三区| 亚洲av成人不卡在线观看播放网| 亚洲黑人精品在线| 99久久精品一区二区三区| 中文字幕人成人乱码亚洲影| 国产淫片久久久久久久久 | 日本在线视频免费播放| 国产黄片美女视频| 十八禁网站免费在线| 欧美一级毛片孕妇| 老汉色∧v一级毛片| 国产高清有码在线观看视频| 中文字幕久久专区| 久久精品国产自在天天线| 国产精品美女特级片免费视频播放器| 国产黄片美女视频| 中文字幕熟女人妻在线| 九色国产91popny在线| 美女免费视频网站| 亚洲人成电影免费在线| 国产精品98久久久久久宅男小说| 亚洲无线在线观看| av黄色大香蕉| 日本与韩国留学比较| 国产欧美日韩一区二区精品| 亚洲欧美日韩东京热| 国产高清视频在线观看网站| 欧美日韩亚洲国产一区二区在线观看| 三级男女做爰猛烈吃奶摸视频| 热99在线观看视频| 母亲3免费完整高清在线观看| svipshipincom国产片| 国产精品99久久99久久久不卡| 男女下面进入的视频免费午夜| 51午夜福利影视在线观看| www国产在线视频色| 国产精品99久久久久久久久| 99国产极品粉嫩在线观看| 哪里可以看免费的av片| 久久久久国内视频| 级片在线观看| 国产视频一区二区在线看| 色吧在线观看| 桃红色精品国产亚洲av| 天堂影院成人在线观看| 亚洲真实伦在线观看| 亚洲欧美日韩无卡精品| 国产成人影院久久av| 欧美日韩精品网址| 黄色视频,在线免费观看| 欧美一级a爱片免费观看看| 国内揄拍国产精品人妻在线| 又爽又黄无遮挡网站| 日韩人妻高清精品专区| www.www免费av| 免费看光身美女| 亚洲性夜色夜夜综合| 精品国内亚洲2022精品成人| 香蕉av资源在线| 欧美成人免费av一区二区三区| 神马国产精品三级电影在线观看| 一本久久中文字幕| 老熟妇乱子伦视频在线观看| 日韩欧美免费精品| 国产激情偷乱视频一区二区| 国产精品久久久人人做人人爽| 人妻久久中文字幕网| 午夜免费成人在线视频| 在线观看日韩欧美| 成年女人看的毛片在线观看| 亚洲精华国产精华精| 婷婷亚洲欧美| a级毛片a级免费在线| 99热精品在线国产| 色吧在线观看| 亚洲国产精品久久男人天堂| 搡女人真爽免费视频火全软件 | 夜夜夜夜夜久久久久| 日本免费一区二区三区高清不卡| 午夜福利视频1000在线观看| 黄色成人免费大全| 狠狠狠狠99中文字幕| 国产熟女xx| 噜噜噜噜噜久久久久久91| 高潮久久久久久久久久久不卡| 免费av毛片视频| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 亚洲av免费高清在线观看| 久久香蕉精品热| 色精品久久人妻99蜜桃| 内射极品少妇av片p| 久久精品综合一区二区三区| 国产色婷婷99| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 成人三级黄色视频| 久久精品国产99精品国产亚洲性色| 神马国产精品三级电影在线观看| 精品一区二区三区视频在线 | 国产精品爽爽va在线观看网站| 欧美国产日韩亚洲一区| 一级a爱片免费观看的视频| 久久性视频一级片| 欧美日本视频| 久久这里只有精品中国| 亚洲国产色片| 好看av亚洲va欧美ⅴa在| 亚洲片人在线观看| 麻豆成人午夜福利视频| 在线播放无遮挡| 亚洲美女视频黄频| 亚洲人成电影免费在线| 欧美乱码精品一区二区三区| 国产一区二区在线观看日韩 | 高清毛片免费观看视频网站| 丁香六月欧美| 少妇人妻一区二区三区视频| www.www免费av| 少妇的逼水好多| 婷婷六月久久综合丁香| bbb黄色大片| 国内精品一区二区在线观看| 99热只有精品国产| 成年免费大片在线观看| 波野结衣二区三区在线 | 18禁在线播放成人免费| 国产精品 欧美亚洲| ponron亚洲| 岛国在线免费视频观看| 免费观看精品视频网站| 精华霜和精华液先用哪个| 国产亚洲av嫩草精品影院| 国产精品亚洲美女久久久| 午夜免费观看网址| 非洲黑人性xxxx精品又粗又长| 国产精品国产高清国产av| 波多野结衣高清作品| 国内精品久久久久精免费| 国产高清视频在线观看网站| 亚洲第一电影网av| 搡女人真爽免费视频火全软件 | 性色avwww在线观看| 日本成人三级电影网站| 午夜老司机福利剧场| 熟女电影av网| 男人的好看免费观看在线视频| 亚洲成人免费电影在线观看| 99riav亚洲国产免费| 亚洲中文字幕日韩| 高清毛片免费观看视频网站| 欧美区成人在线视频| 久久久久九九精品影院| 久久精品91蜜桃| 真人一进一出gif抽搐免费| 国内少妇人妻偷人精品xxx网站| 中亚洲国语对白在线视频| 在线天堂最新版资源| 桃红色精品国产亚洲av| 又紧又爽又黄一区二区| 性色avwww在线观看| 三级国产精品欧美在线观看| 久久久国产成人免费| 日韩人妻高清精品专区| 久久久久性生活片| 美女大奶头视频| 蜜桃久久精品国产亚洲av| 久久精品国产亚洲av香蕉五月| 日本三级黄在线观看| 久久亚洲精品不卡| 日韩欧美国产在线观看| 中文字幕人妻丝袜一区二区| 国产av麻豆久久久久久久| 国产精品爽爽va在线观看网站| 香蕉久久夜色| 成人午夜高清在线视频| av天堂中文字幕网| 午夜福利视频1000在线观看| 国产免费一级a男人的天堂| 午夜激情欧美在线| 色在线成人网| 亚洲avbb在线观看| 色综合亚洲欧美另类图片| 中文在线观看免费www的网站| 十八禁网站免费在线| 看片在线看免费视频| 亚洲狠狠婷婷综合久久图片| 欧美乱码精品一区二区三区| 一夜夜www| 99国产精品一区二区三区| 两个人视频免费观看高清| 麻豆成人av在线观看| 亚洲av免费高清在线观看| 亚洲精品亚洲一区二区| 欧美中文日本在线观看视频| av欧美777| 国内精品久久久久精免费| 欧美一区二区国产精品久久精品| 久久精品国产综合久久久| 51午夜福利影视在线观看| 又黄又爽又免费观看的视频| 久久天躁狠狠躁夜夜2o2o| av在线蜜桃| 99热6这里只有精品|