• <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ù)
    午夜免费男女啪啪视频观看| 少妇高潮的动态图| 麻豆av噜噜一区二区三区| 亚洲国产精品成人久久小说 | 乱码一卡2卡4卡精品| 国产精品永久免费网站| 男人狂女人下面高潮的视频| 在线免费观看不下载黄p国产| 男人和女人高潮做爰伦理| 日韩欧美精品免费久久| 亚洲av熟女| 久久久久久国产a免费观看| 噜噜噜噜噜久久久久久91| 中国国产av一级| 亚洲欧洲日产国产| a级毛片a级免费在线| 中文字幕熟女人妻在线| 久久综合国产亚洲精品| 欧美日韩乱码在线| 亚洲av第一区精品v没综合| 18禁在线播放成人免费| 久久综合国产亚洲精品| 搡老妇女老女人老熟妇| 一个人观看的视频www高清免费观看| 欧美zozozo另类| 欧美成人一区二区免费高清观看| 欧美xxxx黑人xx丫x性爽| 国产老妇伦熟女老妇高清| 亚洲国产精品成人久久小说 | 成人特级黄色片久久久久久久| 日本黄大片高清| 成人毛片a级毛片在线播放| 麻豆久久精品国产亚洲av| 我的老师免费观看完整版| 成人特级黄色片久久久久久久| 麻豆久久精品国产亚洲av| 91午夜精品亚洲一区二区三区| 国产 一区精品| 一区福利在线观看| 久久精品久久久久久噜噜老黄 | 97人妻精品一区二区三区麻豆| 亚洲成人中文字幕在线播放| 亚洲精品成人久久久久久| 菩萨蛮人人尽说江南好唐韦庄 | 久久精品国产亚洲av涩爱 | 一夜夜www| 99国产精品一区二区蜜桃av| 午夜老司机福利剧场| 国产一区二区三区在线臀色熟女| 久久综合国产亚洲精品| 美女xxoo啪啪120秒动态图| 天美传媒精品一区二区| av在线播放精品| 夜夜爽天天搞| 99久久无色码亚洲精品果冻| 精品久久久久久久末码| 久久国内精品自在自线图片| 国产精品美女特级片免费视频播放器| 午夜激情欧美在线| 99热全是精品| 成人二区视频| 国产精品永久免费网站| 精品久久久噜噜| 麻豆成人av视频| h日本视频在线播放| 婷婷精品国产亚洲av| .国产精品久久| 不卡一级毛片| 欧美一区二区精品小视频在线| 亚洲真实伦在线观看| 免费观看a级毛片全部| 午夜a级毛片| 国产精品国产三级国产av玫瑰| 99热只有精品国产| 深夜精品福利| 91久久精品国产一区二区三区| 欧美日韩乱码在线| 日本av手机在线免费观看| 免费人成视频x8x8入口观看| 国产精品久久视频播放| 欧美色视频一区免费| 美女 人体艺术 gogo| 男人狂女人下面高潮的视频| 欧洲精品卡2卡3卡4卡5卡区| 中文在线观看免费www的网站| av在线老鸭窝| 亚洲丝袜综合中文字幕| 少妇高潮的动态图| 成年女人看的毛片在线观看| ponron亚洲| 自拍偷自拍亚洲精品老妇| 久久精品久久久久久噜噜老黄 | 99在线人妻在线中文字幕| 欧美bdsm另类| 日韩欧美三级三区| 国产一区二区在线av高清观看| 天美传媒精品一区二区| 国产成人影院久久av| 乱系列少妇在线播放| 中出人妻视频一区二区| 99久久成人亚洲精品观看| 啦啦啦观看免费观看视频高清| av在线亚洲专区| 小说图片视频综合网站| 亚洲欧美精品自产自拍| 久久精品人妻少妇| 热99在线观看视频| 国产国拍精品亚洲av在线观看| 舔av片在线| 性插视频无遮挡在线免费观看| 亚洲欧美日韩高清专用| 欧美丝袜亚洲另类| 欧美另类亚洲清纯唯美| 国产色爽女视频免费观看| 久久久欧美国产精品| www.av在线官网国产| av.在线天堂| 国产三级在线视频| 美女国产视频在线观看| 国产一区二区激情短视频| 丝袜美腿在线中文| 春色校园在线视频观看| 亚洲人与动物交配视频| 国产精品久久久久久久久免| 国产亚洲av片在线观看秒播厂 | 蜜臀久久99精品久久宅男| www.色视频.com| 亚洲自偷自拍三级| 91在线精品国自产拍蜜月| 男人的好看免费观看在线视频| 麻豆乱淫一区二区| 尤物成人国产欧美一区二区三区| 亚洲人成网站在线观看播放| 日韩,欧美,国产一区二区三区 | 国产熟女欧美一区二区| 又粗又爽又猛毛片免费看| 日韩制服骚丝袜av| 嫩草影院精品99| 中国美女看黄片| 99热全是精品| av黄色大香蕉| 欧美日韩国产亚洲二区| 欧美日韩一区二区视频在线观看视频在线 | 亚洲丝袜综合中文字幕| 欧美高清成人免费视频www| 熟妇人妻久久中文字幕3abv| 看十八女毛片水多多多| 久久精品91蜜桃| 精华霜和精华液先用哪个| 亚洲精品粉嫩美女一区| 亚洲高清免费不卡视频| 99久国产av精品国产电影| 成熟少妇高潮喷水视频| 日韩大尺度精品在线看网址| 黄色一级大片看看| 一区二区三区高清视频在线| 91aial.com中文字幕在线观看| 长腿黑丝高跟| 欧美zozozo另类| 51国产日韩欧美| 尤物成人国产欧美一区二区三区| 欧美高清成人免费视频www| 女同久久另类99精品国产91| 日本与韩国留学比较| 少妇的逼好多水| av视频在线观看入口| 精品人妻熟女av久视频| 高清在线视频一区二区三区 | 青春草国产在线视频 | 熟女电影av网| 久久国产乱子免费精品| 99在线人妻在线中文字幕| 夜夜爽天天搞| 久久精品久久久久久久性| 国产成人精品一,二区 | 日本av手机在线免费观看| 国内精品宾馆在线| 国产v大片淫在线免费观看| 久久这里有精品视频免费| 熟妇人妻久久中文字幕3abv| 国产精品一区二区性色av| 男女做爰动态图高潮gif福利片| 我要看日韩黄色一级片| 国产亚洲av嫩草精品影院| 国产一级毛片在线| 亚洲av中文字字幕乱码综合| 国产黄片美女视频| 免费无遮挡裸体视频| 你懂的网址亚洲精品在线观看 | 一级黄色大片毛片| 网址你懂的国产日韩在线| 高清午夜精品一区二区三区 | 欧美性感艳星| 午夜免费激情av| 人妻久久中文字幕网| 看片在线看免费视频| 欧美日韩综合久久久久久| 男女啪啪激烈高潮av片| 草草在线视频免费看| 青青草视频在线视频观看| 麻豆国产av国片精品| 在线播放国产精品三级| 日日摸夜夜添夜夜添av毛片| 国产熟女欧美一区二区| 深夜a级毛片| 精品久久久久久久久av| 一个人观看的视频www高清免费观看| 三级毛片av免费| 久久久精品大字幕| 国产黄片美女视频| 亚洲国产精品sss在线观看| 一边亲一边摸免费视频| 免费人成视频x8x8入口观看| 久久精品国产亚洲网站| 精品久久久久久久久久免费视频| 天天一区二区日本电影三级| 亚洲精品乱码久久久久久按摩| 我的女老师完整版在线观看| 69人妻影院| 午夜激情欧美在线| 精品国内亚洲2022精品成人| 99视频精品全部免费 在线| 国产精品无大码| 成人特级av手机在线观看| 男人的好看免费观看在线视频| 国产精品久久电影中文字幕| 国产伦理片在线播放av一区 | 久久中文看片网| 国产精品女同一区二区软件| 老司机影院成人| 午夜久久久久精精品| 国产白丝娇喘喷水9色精品| 最好的美女福利视频网| or卡值多少钱| 一进一出抽搐动态| 男女那种视频在线观看| av福利片在线观看| 在线观看一区二区三区| 日日摸夜夜添夜夜爱| 女的被弄到高潮叫床怎么办| 国产伦在线观看视频一区| 国产女主播在线喷水免费视频网站 | 精品人妻偷拍中文字幕| 国产精品人妻久久久影院| 成人无遮挡网站| 在线播放国产精品三级| 亚洲色图av天堂| 哪个播放器可以免费观看大片| 日韩强制内射视频| 能在线免费观看的黄片| 亚洲人成网站高清观看| 大香蕉久久网| av免费观看日本| 欧美激情在线99| 国产亚洲欧美98| 中文资源天堂在线| 此物有八面人人有两片| 亚洲欧洲日产国产| 国产美女午夜福利| 精品久久久久久久人妻蜜臀av| 欧美三级亚洲精品| 高清日韩中文字幕在线| 久久午夜福利片| 又爽又黄a免费视频| 熟女电影av网| 久久久成人免费电影| 国产成人影院久久av| 乱码一卡2卡4卡精品| 国产午夜福利久久久久久| 亚洲精品色激情综合| 99久国产av精品国产电影| 啦啦啦韩国在线观看视频| 色哟哟哟哟哟哟| 九色成人免费人妻av| 精品欧美国产一区二区三| 中文在线观看免费www的网站| 国产高清不卡午夜福利| 91麻豆精品激情在线观看国产| 禁无遮挡网站| 国产成人a区在线观看| 美女脱内裤让男人舔精品视频 | 久久韩国三级中文字幕| 婷婷色av中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 久久精品人妻少妇| 亚洲国产欧美人成| 午夜激情欧美在线| 国产在线男女| 自拍偷自拍亚洲精品老妇| 成人漫画全彩无遮挡| 国产单亲对白刺激| 欧美激情国产日韩精品一区| 精品人妻视频免费看| av卡一久久| av在线老鸭窝| 91久久精品电影网| 有码 亚洲区| 热99re8久久精品国产| 青春草视频在线免费观看| 欧美激情久久久久久爽电影| 国产高潮美女av| 久久久久久久午夜电影| 乱系列少妇在线播放| 91精品一卡2卡3卡4卡| 国产 一区精品| 91久久精品国产一区二区成人| 国产精品永久免费网站| 国产熟女欧美一区二区| 精品久久久久久久久av| 99久久无色码亚洲精品果冻| 日本黄色片子视频| 亚洲最大成人中文| 国产伦精品一区二区三区视频9| 中国美女看黄片| 国产高清三级在线| av在线蜜桃| 菩萨蛮人人尽说江南好唐韦庄 | 国产成人91sexporn| 美女内射精品一级片tv| 天堂影院成人在线观看| 九九在线视频观看精品| 国产精品人妻久久久影院| 久久久国产成人免费| 在线国产一区二区在线| 美女大奶头视频| 日日摸夜夜添夜夜爱| 亚洲国产精品国产精品| 国产极品精品免费视频能看的| 麻豆国产97在线/欧美| 夜夜夜夜夜久久久久| 少妇熟女aⅴ在线视频| 深夜a级毛片| 嫩草影院入口| 国产成人精品一,二区 | 国产人妻一区二区三区在| 国产精品不卡视频一区二区| 亚洲精品自拍成人| 免费大片18禁| 直男gayav资源| 我的老师免费观看完整版| 综合色av麻豆| 午夜福利成人在线免费观看| 你懂的网址亚洲精品在线观看 | 丰满乱子伦码专区| 久久精品久久久久久久性| 变态另类成人亚洲欧美熟女| 国产成人精品婷婷| 亚洲国产精品成人久久小说 | 久久久久久久午夜电影| 中出人妻视频一区二区| av.在线天堂| 久久久精品欧美日韩精品| 特大巨黑吊av在线直播| 高清午夜精品一区二区三区 | 国产欧美日韩精品一区二区| 中文欧美无线码| 六月丁香七月| 久久久久久久亚洲中文字幕| 99国产极品粉嫩在线观看| 一个人观看的视频www高清免费观看| 国产精品麻豆人妻色哟哟久久 | 成熟少妇高潮喷水视频| 欧美日韩综合久久久久久| 精品久久久久久久久av| 一级av片app| 国产精品久久久久久久电影| 精品久久国产蜜桃| 欧美成人精品欧美一级黄| 在线天堂最新版资源| 国产视频首页在线观看| 亚洲人成网站在线播| 亚洲欧美精品自产自拍| 听说在线观看完整版免费高清| 亚洲精品乱码久久久v下载方式| 22中文网久久字幕| 日本爱情动作片www.在线观看| 丝袜喷水一区| 2022亚洲国产成人精品| 亚洲高清免费不卡视频| 看片在线看免费视频| 久久亚洲精品不卡| av女优亚洲男人天堂| 搞女人的毛片| 日韩三级伦理在线观看| 中文字幕精品亚洲无线码一区| 九草在线视频观看| 免费电影在线观看免费观看| 韩国av在线不卡| 日本成人三级电影网站| 色综合站精品国产| 九九爱精品视频在线观看| 欧美精品一区二区大全| 日韩av在线大香蕉| 欧美最黄视频在线播放免费| 国产亚洲5aaaaa淫片| 在线观看午夜福利视频| 波多野结衣巨乳人妻| 久久精品国产清高在天天线| 最近的中文字幕免费完整| 黄色日韩在线| a级毛片免费高清观看在线播放| 亚洲第一电影网av| 精品人妻偷拍中文字幕| 国产精品久久久久久av不卡| 欧美高清成人免费视频www| 97超视频在线观看视频| 成人二区视频| 国内精品美女久久久久久| 国产黄片视频在线免费观看| 国产又黄又爽又无遮挡在线| 黄片无遮挡物在线观看| 男人和女人高潮做爰伦理| 国产毛片a区久久久久| 亚洲精品国产成人久久av| 国内久久婷婷六月综合欲色啪| 2021天堂中文幕一二区在线观| 久久久久网色| 国产精品一区二区三区四区免费观看| 国产成人福利小说| 午夜久久久久精精品| 最好的美女福利视频网| 禁无遮挡网站| 大香蕉久久网| 乱人视频在线观看| 日本一本二区三区精品| www日本黄色视频网| 午夜精品在线福利| 少妇高潮的动态图| a级毛片免费高清观看在线播放| 精品人妻一区二区三区麻豆| 蜜臀久久99精品久久宅男| 久久久色成人| 精品不卡国产一区二区三区| 综合色丁香网| 麻豆一二三区av精品| av免费在线看不卡| 欧美3d第一页| 两个人的视频大全免费| 毛片一级片免费看久久久久| 真实男女啪啪啪动态图| 永久网站在线| 欧美在线一区亚洲| 亚洲欧美成人精品一区二区| 黑人高潮一二区| 欧美丝袜亚洲另类| 91精品一卡2卡3卡4卡| 菩萨蛮人人尽说江南好唐韦庄 | 欧美日本亚洲视频在线播放| 男人的好看免费观看在线视频| 精品少妇黑人巨大在线播放 | 欧美性猛交黑人性爽| 老女人水多毛片| 网址你懂的国产日韩在线| 国国产精品蜜臀av免费| 亚洲国产精品久久男人天堂| 亚洲精品久久久久久婷婷小说 | 长腿黑丝高跟| 成人亚洲精品av一区二区| 亚洲精品亚洲一区二区| 中国美女看黄片| 一个人看视频在线观看www免费| 亚洲无线在线观看| 99热这里只有是精品50| 日韩精品有码人妻一区| 校园人妻丝袜中文字幕| 国产日本99.免费观看| 国产亚洲精品久久久久久毛片| 日本爱情动作片www.在线观看| 国产av麻豆久久久久久久| 国产成人91sexporn| 成人二区视频| 国产女主播在线喷水免费视频网站 | 99热网站在线观看| 桃色一区二区三区在线观看| 精品熟女少妇av免费看| 日韩av在线大香蕉| kizo精华| 99在线人妻在线中文字幕| 内地一区二区视频在线| 只有这里有精品99| 国产真实乱freesex| 国产亚洲精品久久久久久毛片| av在线播放精品| 亚洲国产日韩欧美精品在线观看| 亚洲国产欧洲综合997久久,| 亚洲美女搞黄在线观看| 人人妻人人澡欧美一区二区| 精品无人区乱码1区二区| 91久久精品电影网| 亚洲欧美日韩高清专用| 三级男女做爰猛烈吃奶摸视频| av视频在线观看入口| 久久久久免费精品人妻一区二区| 精品午夜福利在线看| 久久久久久久久久久丰满| 欧美日韩精品成人综合77777| 亚洲成人久久爱视频| 日本黄大片高清| 日韩欧美精品免费久久| 看十八女毛片水多多多| 日本爱情动作片www.在线观看| 亚洲性久久影院| 最近视频中文字幕2019在线8| 亚洲人成网站在线观看播放| 精品久久久久久久久av| 欧美潮喷喷水| 国产一区二区亚洲精品在线观看| 亚洲人成网站高清观看| 久久中文看片网| 老师上课跳d突然被开到最大视频| 午夜a级毛片| 人妻久久中文字幕网| 婷婷色综合大香蕉| 亚洲欧美精品综合久久99| 亚洲成人精品中文字幕电影| 亚洲一区高清亚洲精品| 最好的美女福利视频网| 少妇熟女欧美另类| 波多野结衣高清无吗| 乱系列少妇在线播放| 国产精品久久久久久久电影| 最后的刺客免费高清国语| 白带黄色成豆腐渣| 中国国产av一级| 日韩在线高清观看一区二区三区| 哪里可以看免费的av片| 男女那种视频在线观看| 高清日韩中文字幕在线| 天堂影院成人在线观看| 国产成人91sexporn| 蜜桃久久精品国产亚洲av| 婷婷六月久久综合丁香| 精品久久久久久成人av| eeuss影院久久| 成人亚洲精品av一区二区| 男女下面进入的视频免费午夜| 免费黄网站久久成人精品| 久久精品国产自在天天线| 晚上一个人看的免费电影| 中文欧美无线码| 国产单亲对白刺激| 精华霜和精华液先用哪个| 国产精品久久久久久久久免| 亚洲欧美日韩无卡精品| 色视频www国产| 欧美+亚洲+日韩+国产| 日韩 亚洲 欧美在线| 日本免费一区二区三区高清不卡| 国产精品一区www在线观看| 大型黄色视频在线免费观看| 51国产日韩欧美| 欧美+亚洲+日韩+国产| 韩国av在线不卡| 欧美日韩在线观看h| 嘟嘟电影网在线观看| 在线观看一区二区三区| 国产91av在线免费观看| 婷婷六月久久综合丁香| 精品人妻一区二区三区麻豆| 美女黄网站色视频| 日韩视频在线欧美| 国产一区二区在线av高清观看| av专区在线播放| 激情 狠狠 欧美| 久久久久久久久久成人| av视频在线观看入口| 国产免费男女视频| 成人午夜高清在线视频| 婷婷精品国产亚洲av| 午夜老司机福利剧场| 免费无遮挡裸体视频| 国产精品久久视频播放| 成人特级黄色片久久久久久久| 午夜精品在线福利| 国产乱人视频| 麻豆精品久久久久久蜜桃| 91aial.com中文字幕在线观看| 国产精品蜜桃在线观看 | 国产成人a∨麻豆精品| 在线观看66精品国产| 一本久久中文字幕| 三级国产精品欧美在线观看| 熟女人妻精品中文字幕| 日产精品乱码卡一卡2卡三| 又爽又黄无遮挡网站| 国产不卡一卡二| 午夜久久久久精精品| 亚洲精品色激情综合| 国产欧美日韩精品一区二区| av在线天堂中文字幕| 黄色一级大片看看| 中国美女看黄片| 亚洲国产高清在线一区二区三| 久久久久网色| 男女啪啪激烈高潮av片| 啦啦啦啦在线视频资源| 亚洲欧美日韩无卡精品| 欧洲精品卡2卡3卡4卡5卡区| 美女国产视频在线观看| 真实男女啪啪啪动态图| 99在线人妻在线中文字幕| 久久久欧美国产精品| 嫩草影院新地址| 老熟妇乱子伦视频在线观看| 亚洲高清免费不卡视频| 国产午夜福利久久久久久| 国产美女午夜福利| 国产三级在线视频| 亚洲美女搞黄在线观看| 久久婷婷人人爽人人干人人爱| 亚洲av成人精品一区久久| 久久久久久大精品| 天天一区二区日本电影三级| 国产高清有码在线观看视频| 日韩国内少妇激情av|