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

    基于SRAP和SSR標(biāo)記的小麥粒長和千粒質(zhì)量QTL定位及效應(yīng)分析

    2017-09-15 03:14:58郭利建王竹林汪世娟劉香利趙惠賢
    關(guān)鍵詞:粒長染色體遺傳

    郭利建,王竹林,汪世娟,劉香利,趙惠賢

    (1. 西北農(nóng)林科技大學(xué) 生命科學(xué)學(xué)院,陜西楊凌 712100;2. 西北農(nóng)林科技大學(xué) 農(nóng)學(xué)院,陜西楊凌 712100)

    基于SRAP和SSR標(biāo)記的小麥粒長和千粒質(zhì)量QTL定位及效應(yīng)分析

    郭利建1,王竹林2,汪世娟1,劉香利1,趙惠賢1

    (1. 西北農(nóng)林科技大學(xué) 生命科學(xué)學(xué)院,陜西楊凌 712100;2. 西北農(nóng)林科技大學(xué) 農(nóng)學(xué)院,陜西楊凌 712100)

    為探究小麥粒長和千粒質(zhì)量性狀的QTL,以千粒質(zhì)量差異較大的小麥品種‘西農(nóng)981’和‘陜麥159’雜交構(gòu)建的169株F2群體和F2:3家系為研究材料,利用SRAP標(biāo)記和SSR標(biāo)記進(jìn)行遺傳圖譜的構(gòu)建,通過統(tǒng)計(jì)分析楊凌及三原2個(gè)環(huán)境下的F2:3家系的表型數(shù)據(jù),利用完備區(qū)間作圖法對(duì)粒長和千粒質(zhì)量進(jìn)行QTL定位。結(jié)果表明,2個(gè)環(huán)境條件下,親本‘西農(nóng)981’和‘陜麥159’在粒長和千粒質(zhì)量上均表現(xiàn)出極顯著差異,F(xiàn)2:3家系表現(xiàn)出明顯的超親分離現(xiàn)象。QTL定位共檢測到31個(gè)粒長和千粒質(zhì)量相關(guān)的QTL,其中粒長相關(guān)QTL檢測到7個(gè),分布于5A、6A、7A、2B、3B、4B染色體上,可解釋表型變異的4.36%~14.80%,在3B染色體上檢測到1個(gè)能在2個(gè)環(huán)境點(diǎn)穩(wěn)定表達(dá)的粒長QTL;千粒質(zhì)量相關(guān)QTL檢測到24個(gè),分布于2A、3A、5A、6A、7A、2B、3B、4B、5B、7B染色體上,可解釋表型變異的0.25%~8.23%。另外,在2B染色體上檢測到5個(gè)控制千粒質(zhì)量的QTL,表明2B染色體與千粒質(zhì)量關(guān)系密切。

    小麥;粒長;千粒質(zhì)量;SRAP;QTL

    小麥?zhǔn)鞘澜缛笾骷Z之一,其產(chǎn)量直接關(guān)系著糧食安全問題,因此,尋找控制小麥產(chǎn)量性狀的基因,提高小麥產(chǎn)量一直是科研工作者的目標(biāo)。千粒質(zhì)量是小麥產(chǎn)量的三大主要組成因素之一,對(duì)小麥的產(chǎn)量具有極其重要的影響,而籽粒的大小及形態(tài)則直接影響著千粒質(zhì)量的大小[1],粒長(kernel length,KL)作為籽粒大小的重要標(biāo)志,與千粒質(zhì)量存在重要的相關(guān)性[2],因此,通過研究粒長及千粒質(zhì)量相關(guān)的QTL可為后續(xù)的產(chǎn)量研究提供基礎(chǔ)。有關(guān)小麥產(chǎn)量相關(guān)性狀的QTL前人已有較多研究,檢測到的控制小麥千粒質(zhì)量的QTL幾乎遍布于21條染色體上[3-7], Campbell等[3]將控制粒長的QTL定位到1B、2B、2D、3B、7D染色體上。Tyagi等[8]利用‘中國春’和‘Rye Selection 111’所構(gòu)建的92個(gè)重組自交系為材料,利用SSR(基因間重復(fù)序列多態(tài)性)標(biāo)記和AFLP標(biāo)記共檢測到45個(gè)與粒長、粒寬、容積相關(guān)的QTL,這些QTL分布在除2D和3D外的19條染色體上,可解釋表型變異的6.97%~29.87%。Li等[9]利用594個(gè)SNP標(biāo)記和404個(gè)SSR標(biāo)記構(gòu)建一張高密度遺傳連鎖圖譜,檢測到14個(gè)與千粒質(zhì)量相關(guān)的QTL,2個(gè)與穗數(shù)相關(guān)的QTL,5個(gè)與穗粒數(shù)相關(guān)QTL,這些加性QTL分布于13條染色體上。姚琴等[10]利用SSR和DArT分子標(biāo)記,通過2個(gè)環(huán)境點(diǎn)數(shù)據(jù)共檢測到4個(gè)千粒質(zhì)量QTL,分布于3B、4D、6D染色體上,可解釋表型變異的4.79%~31.37%。

    前人雖已檢測到許多控制小麥粒長和千粒質(zhì)量相關(guān)的QTL,但這些性狀遺傳基礎(chǔ)復(fù)雜,易受環(huán)境及材料等影響,因此用不同研究材料及不同分子標(biāo)記揭示能在不同環(huán)境下穩(wěn)定表達(dá)的QTL及其效應(yīng)仍有重要價(jià)值。SRAP(基因表達(dá)區(qū)序列)標(biāo)記是由Li等[11]首先于蕓薹屬植物中建立,是一種基于基因表達(dá)區(qū)序列的功能性分子標(biāo)記。目前多被應(yīng)用于植物遺傳多樣性研究[12-13],目前小麥產(chǎn)量性狀QTL研究主要采用SSR、AFLP等分子標(biāo)記。而在小麥數(shù)量性狀QTL定位研究應(yīng)用SRAP標(biāo)記的報(bào)道較少。本研究利用近年來審定的農(nóng)藝性狀優(yōu)良的小麥品種‘西農(nóng)981’和‘陜麥159’構(gòu)建的F2群體和F2:3家系,通過SRAP標(biāo)記和SSR標(biāo)記相結(jié)合的方法進(jìn)行遺傳圖譜的構(gòu)建,在2個(gè)環(huán)境條件下對(duì)小麥千粒質(zhì)量和粒長性狀進(jìn)行QTL定位及效應(yīng)分析。旨在尋找能在不同環(huán)境條件下穩(wěn)定表達(dá)的QTL位點(diǎn),為分子標(biāo)記輔助小麥的遺傳改良提供參考信息。

    1 材料與方法

    1.1 植物材料及種植

    用近年來審定的農(nóng)藝性狀優(yōu)良的2個(gè)小麥品種‘西農(nóng)981’和‘陜麥159’進(jìn)行雜交,構(gòu)建含有169個(gè)單株的F2群體和F2:3家系;分別用于作圖群體和表型定位。

    2012年4月以‘西農(nóng)981’為母本,以‘陜麥159’為父本,進(jìn)行雜交。6月收獲雜交種子,10月將雜交種子種植于西北農(nóng)林科技大學(xué)楊凌小麥試驗(yàn)田。2013年6月收獲F1種子,同年10月,將F1種子種植于楊凌小麥試驗(yàn)田,同時(shí)種植2個(gè)親本,田間常規(guī)管理。2014年4月分單株掛牌標(biāo)記169個(gè)F2單株,同年6月按單株收獲并脫粒。2014年10月將收獲的169個(gè)F2單株的種子一分為二,分別種植于陜西楊凌試驗(yàn)站(環(huán)境E1)(34°17′N,108°4′E,海拔 525 m),和陜西三原試驗(yàn)站(環(huán)境E2)(34°36′N,108°52′E,海拔429 m),形成2個(gè)環(huán)境點(diǎn)的F2:3家系群體;采用隨機(jī)區(qū)組分布,試驗(yàn)設(shè)3次重復(fù);單行區(qū),行長2 m,行距0.24 m;同時(shí)種植親本‘西農(nóng)981’和‘陜麥159’,‘西農(nóng)979’作對(duì)照,田間常規(guī)管理。

    1.2 SRAP引物與SSR引物

    SRAP引物信息來源于Li等[11]以及西北農(nóng)林科技大學(xué)農(nóng)學(xué)院胡勝武實(shí)驗(yàn)室,其中正向引物34條,反向引物18條,共能組成612對(duì)引物,引物由英濰捷基(上海)貿(mào)易有限公司合成。

    SSR引物包括Xbarc系列、Xwmc系列、Xgdm系列、Xcfa系列、Xcfd系列等共計(jì)1 550對(duì),由西北農(nóng)林科技大學(xué)農(nóng)學(xué)院王竹林老師提供。

    1.3 DNA提取及分子標(biāo)記檢測

    2個(gè)親本、F1代植株以及F2分離群體中掛牌各單株分別采取幼苗葉片,采用CTAB法[14]進(jìn)行小麥基因組DNA的提取,提取后通過瓊脂糖凝膠電泳及超微量微孔板分光光度計(jì)分析DNA純度及測量DNA濃度。

    SRAP引物PCR反應(yīng)體系為10 μL:ddH2O 2.4 μL,2×EsTaqMasterMix(康為世紀(jì))5 μL,SRAP引物各0.3 μL(10 μmol·L-1),模版DNA 2 μL(20 ng· μL-1)。SRAP引物擴(kuò)增程序?yàn)椋?5 ℃預(yù)變性60 s,94 ℃變性60 s,35 ℃退火60 s,72 ℃延伸60 s,5個(gè)循環(huán);94 ℃變性60 s,50 ℃退火60 s,72 ℃延伸60 s,35個(gè)循環(huán);72 ℃延伸7 min;4 ℃保存。

    SSR引物PCR反應(yīng)體系為15 μL,ddH2O 5.5 μL,2×EsTaqMasterMix(康為世紀(jì))7.5 μL,SSR引物各0.5 μL(10 μmol·L-1),模版DNA 1 μL(20ng· μL-1)。SSR引物擴(kuò)增程序?yàn)椋?5 ℃預(yù)變性4 min,94 ℃變性45 s,60 ℃退火60 s(每個(gè)循環(huán)降低0.5 ℃),72 ℃延伸55 s,10個(gè)循環(huán);94 ℃變性40 s,55 ℃退火40 s,72 ℃延伸60 s,35個(gè)循環(huán);72 ℃延伸10 min;4 ℃保存。

    擴(kuò)增產(chǎn)物用80 g/L的聚丙烯酰胺凝膠電泳分離,銀染法檢測并拍照統(tǒng)計(jì)擴(kuò)增條帶。

    1.4 性狀調(diào)查

    2015年6月分別收獲2個(gè)環(huán)境點(diǎn)的169個(gè)F2:3家系;每株系隨機(jī)選取6個(gè)單株分別收獲、脫粒及進(jìn)行性狀分析。粒長采用游標(biāo)卡尺單粒測定,每個(gè)單株重復(fù)測定30次,取平均值;千粒質(zhì)量測定采用200粒種子稱量,重復(fù)3次,計(jì)算平均值并換算成千粒質(zhì)量。上述性狀調(diào)查結(jié)果取各株系平均值,采用SPSS 19.0軟件進(jìn)行統(tǒng)計(jì)分析。

    1.5 連鎖圖譜構(gòu)建及QTL定位

    遺傳連鎖圖譜構(gòu)建采用軟件Icimapping4.1(www.isbreeding.net)。帶型記錄參照軟件要求;各參數(shù)值依照軟件初始設(shè)置,LOD≥2.5。用Kosambi函數(shù)和完備區(qū)間作圖法作圖[15]。

    2 結(jié)果與分析

    2.1 親本及F2:3家系產(chǎn)量性狀表型變異

    親本及F2:3家系產(chǎn)量性狀表型分析見表1。親本‘西農(nóng)981’在楊凌環(huán)境點(diǎn)和三原環(huán)境點(diǎn)粒長分別為7.32 mm和7.45 mm,而親本‘陜麥159’在這2個(gè)環(huán)境點(diǎn)的粒長分別為6.25 mm和6.22 mm;親本‘西農(nóng)981’在2環(huán)境點(diǎn)的千粒質(zhì)量分別為50.95 g和50.05 g,親本‘陜麥159’在2個(gè)環(huán)境點(diǎn)的千粒質(zhì)量分別為42.3 g和42.51 g。

    雙親中,粒長、千粒質(zhì)量親本‘西農(nóng)981’均顯著高于‘陜麥159’。對(duì)親本間進(jìn)行單因素方差分析表明,粒長和千粒質(zhì)量在2環(huán)境條件下均達(dá)到極顯著差異。F2:3家系中,粒長性狀在2環(huán)境點(diǎn)下變異系數(shù)分別為2.45%和2.48%;千粒質(zhì)量性狀在2環(huán)境點(diǎn)下變異系數(shù)分別為8.19%和6.65%,且表現(xiàn)出超親分離現(xiàn)象,說明雙親及F2:3家系存在真實(shí)的遺傳差異,符合數(shù)量性狀遺傳分析的要求。

    表1 2個(gè)環(huán)境條件下親本及F2:3群體表型及分析Table 1 Phenotypic analysis of F2:3 population and parents in two environment sites

    注:**極顯著差異。

    Note:** level of highly significant difference.

    2.2 分子標(biāo)記連鎖圖譜

    利用SRAP引物和SSR引物在‘西農(nóng)981’和‘陜麥159’2個(gè)親本之間進(jìn)行單重復(fù)篩選及雙重復(fù)驗(yàn)證,最終篩選獲得81對(duì)有多態(tài)性的SRAP引物和150對(duì)有多態(tài)性SSR引物。將這些在親本間有多態(tài)性的引物再用F2群體進(jìn)行檢測,SRAP引物有66對(duì)出現(xiàn)多態(tài)性,SSR引物有63對(duì)具有多態(tài)性。其中,SSR引物和SRAP引物在F2群體的部分個(gè)體中擴(kuò)增多態(tài)性如圖1和圖2所示。最終在F2群體存在115個(gè)多態(tài)性標(biāo)記(55個(gè)SRAP標(biāo)記,60個(gè)SSR標(biāo)記),這些標(biāo)記被正確定位到遺傳連鎖圖譜上。圖譜全長3 203.67 cM,標(biāo)記主要定位于A染色體組和B染色體組,其中2B染色體標(biāo)記數(shù)量最多,共有21個(gè)標(biāo)記,D基因組1D染色體最少,定位到3個(gè)標(biāo)記。

    P1.西農(nóng)981 Xinong 981;P2.陜麥159 Shaanmai 159;1~88.部分F2群體 Part of F2 population;下圖同 The same below 圖1 SSR引物Xgwm66在F2群體的部分個(gè)體之間的擴(kuò)增多態(tài)性檢測圖譜Fig.1 Polymorphism results of SSR primer pair Xgwm66 in partial F2 population

    圖2 SRAP引物對(duì)Me11Em1在F2群體的部分個(gè)體之間的擴(kuò)增多態(tài)性檢測圖譜Fig.2 Polymorphism results of SRAP primer pair Me11Em1 in partial F2 population

    2.3 小麥粒長和千粒質(zhì)量性狀的QTL定位及效應(yīng)分析

    在楊凌和三原2個(gè)環(huán)境點(diǎn)共檢測到31個(gè)粒長和千粒質(zhì)量相關(guān)QTL(圖3及表2);其中粒長相關(guān)QTL檢測到7個(gè),分布于5A、6A、7A、2B、3B、4B染色體上;千粒質(zhì)量相關(guān)QTL檢測到24個(gè),分布于2A、3A、5A、6A、7A、2B、3B、4B、5B、7B染色體上。

    粒長QTL中,楊凌環(huán)境點(diǎn)只檢測到1個(gè)QTL,位于3B染色體Xgwm376-Em1Me10-3,可解釋表型變異的8.18%。三原環(huán)境點(diǎn)檢測到6個(gè)粒長相關(guān)QTL,分別分布于2B、3B、4B、5A、6A和7A染色體上;其中,在三原環(huán)境點(diǎn)3B染色體Xgwm376-Em1Me10-3同樣檢測到1個(gè)粒長QTL,可解釋表型變異的5.47%,表明該控制粒長的QTL能在2個(gè)不同環(huán)境條件下穩(wěn)定表達(dá);在5A染色體Xgwm126-Xgwm291檢測到1個(gè)粒長QTL,可解釋表型變異的5.81%;在6A染色體上檢測到1個(gè)粒長QTL,位于Em8Me18-2-Em5Me28-2,可解釋表型變異的4.36%。在7A染色體上檢測到1個(gè)粒長QTL,位于WMC273-Xgwm282-1,可解釋表型變異的4.69%。在2B染色體上檢測到1個(gè)粒長QTL,位于Em8Me32-Em9Me12,可解釋表型變異的14.8%??赡転橐恢餍TL。

    ●.楊凌粒長 Kernel length measured in Yangling;▲.三原粒長 Kernel length measured in Sanyuan;◆.楊凌千粒質(zhì)量 1 000-grain mass measured in Yangling▼.三原千粒質(zhì)量 1 000-grain mass measured in Sanyuan

    圖3 小麥粒長和千粒質(zhì)量相關(guān)QTL位點(diǎn)定位Fig.3 Location of QTLs related kernel length and 1 000-grain mass in wheat表2 完備區(qū)間作圖法檢測粒長和千粒質(zhì)量的QTLTable 2 QTLs determined by inclusive composite interval mapping

    千粒質(zhì)量QTL中,在楊凌環(huán)境點(diǎn)共檢測到21個(gè)QTL。其中3A染色體上1個(gè);5A染色體上3個(gè);6A染色體上1個(gè),7A染色體上1個(gè);2B染色體上4個(gè),3B染色體上3個(gè);4B染色體上2個(gè),5B染色體上2個(gè);7B染色體上4個(gè)。加性效應(yīng)值表明,千粒質(zhì)量QTL的增效等位基因來源于親本‘西農(nóng)981’。三原環(huán)境點(diǎn)檢測到3個(gè)千粒質(zhì)量QTL,分別位于2A、2B、4B染色體上。

    3 討 論

    3.1 SRAP標(biāo)記能有效增加基于SSR標(biāo)記的小麥遺傳圖譜中分子標(biāo)記的密度

    SSR標(biāo)記是基于高等植物基因組中基因間短的中度或高度重復(fù)序列在物種間的多態(tài)性而建立的分子標(biāo)記,其在作物遺傳多樣性、遺傳連鎖圖的構(gòu)建及數(shù)量性狀QTL定位、目標(biāo)基因定位及分子標(biāo)記輔助育種選擇等方面已得到廣泛應(yīng)用[16]。SRAP標(biāo)記技術(shù)是首先在蕓薹屬植物中建立的,李媛媛等[17]利用SRAP標(biāo)記、SSR標(biāo)記和AFLP標(biāo)記對(duì)甘藍(lán)型油菜進(jìn)行遺傳圖譜的構(gòu)建,發(fā)現(xiàn)SRAP多態(tài)性檢出率最高。本研究構(gòu)建小麥遺傳連鎖圖譜時(shí),多態(tài)性篩選結(jié)果表明SRAP標(biāo)記的多態(tài)性檢出率高于SSR標(biāo)記。SRAP引物不受物種限制,不同引物可以隨意組合,大大提高引物的使用效率,進(jìn)而降低引物合成的成本。因此,SRAP標(biāo)記在作物遺傳圖譜構(gòu)建中具有很大的利用空間。在小麥遺傳連鎖圖譜構(gòu)建中,SRAP標(biāo)記能夠彌補(bǔ)SSR標(biāo)記的不足,使基于SSR標(biāo)記的遺傳連鎖圖的標(biāo)記密度能夠進(jìn)一步增加。

    此外,由于本試驗(yàn)所用SSR標(biāo)記多位于A染色體組和B染色體組,因此所得遺傳圖譜中,分子標(biāo)記多位于A染色體組和B染色體組,D染色體組較少,賈繼增等[18]對(duì)來自世界11個(gè)國家的小麥品種進(jìn)行遺傳多樣性分析,結(jié)果發(fā)現(xiàn)A染色體組和B染色體組遺傳多樣性較高,而D基因組較差,這可能是導(dǎo)致D染色體組標(biāo)記較少的原因之一。

    3.2 小麥粒長和千粒質(zhì)量QTL定位結(jié)果比較

    由于小麥基因組龐大且復(fù)雜,因此利用不同的材料和不同的方法對(duì)小麥進(jìn)行QTL定位得到的結(jié)果往往有所差異并可能得到新的結(jié)果,這與不同的材料攜帶的基因不同有關(guān),同時(shí)也受環(huán)境差異和不同標(biāo)記方法的共同影響。本研究利用近年來審定的農(nóng)藝性狀優(yōu)良的小麥品種‘西農(nóng)981’和‘陜麥159’構(gòu)建的F2群體和F2:3家系為研究材料,在三原和楊凌2個(gè)環(huán)境條件下對(duì)小麥千粒質(zhì)量和粒長性狀進(jìn)行QTL定位及效應(yīng)分析,結(jié)果共檢測到31個(gè)小麥粒長和千粒質(zhì)量QTL位點(diǎn);其中,其7個(gè)粒長QTL;24個(gè)千粒質(zhì)量QTL。前人已研究得到的粒長QTL基本遍布于小麥21條染色體上[19-21]。余曼麗等[20]將粒長QTL定位到5A、7A、1B、2B、3B、1D、6D、7D染色體上,其中5A、1B、2B、3B染色體上的粒長QTL可在2環(huán)境條件下穩(wěn)定表達(dá);本研究同樣在5A、7A、2B、3B染色體上檢測到粒長QTL,但是只在3B染色體上檢測到1個(gè)能在2環(huán)境條件下穩(wěn)定表達(dá)的粒長QTL;本研究在3B染色體所得的粒長QTL與余曼麗等[20]在3B染色體所得的粒長QTL所處位置不同,鄰近標(biāo)記也不同,但是否為新的QTL仍需要進(jìn)一步增加標(biāo)記密度才能確定。此外,本研究在2B染色體上的粒長QTL可解釋表型變異的14.8%,可能為一主效QTL;Ramya等[21]利用SSR標(biāo)記在小麥2B染色相近位置同樣發(fā)現(xiàn)了一個(gè)控制粒長的QTL;但是由于本研究在該位置臨近標(biāo)記為SRAP標(biāo)記,所以,無法確定二者是否為同一QTL位點(diǎn),還需進(jìn)一步研究分析才能明確。本研究除在5A、7A、2B、3B染色體上檢測到粒長QTL外,還在6A、4B染色體上檢測到控制粒長的QTL,其中位于6A染色體上的QTL前人未檢測到,可能為一新的粒長QTL。

    已報(bào)道的小麥千粒質(zhì)量QTL分布于小麥21條染色體上[3,22-23]。Ramya等[21]將小麥千粒質(zhì)量QTL定位于1A、2B、4B、5B、6B、1D、2D染色體上,其中位于2B染色體上的千粒質(zhì)量QTL能夠在4個(gè)環(huán)境條件下穩(wěn)定表達(dá);周淼平等[24]和丁安明等[25]也在2B染色體上檢測到千粒質(zhì)量相關(guān)QTL;本研究在2B染色體上未檢測到2環(huán)境下穩(wěn)定表達(dá)的千粒質(zhì)量QTL,但是在2B染色體共檢測到5個(gè)千粒質(zhì)量QTL,這5個(gè)千粒質(zhì)量QTL與Ramya等[21]以及周淼平等[24]和丁安明等[25]所檢測得到的千粒質(zhì)量QTL所處位置均不相同,這5個(gè)千粒質(zhì)量QTL多數(shù)與SRAP標(biāo)記相關(guān),由于SRAP標(biāo)記與其他標(biāo)記如SSR標(biāo)記的擴(kuò)增區(qū)間不同,因此推測可能為新的千粒質(zhì)量QTL,但是其貢獻(xiàn)率都較低,屬于微效QTL。王瑞霞等[2]曾在3B染色體上檢測到3個(gè)千粒質(zhì)量QTL;本研究也在3B染色體上檢測到3個(gè)千粒質(zhì)量QTL,這表明2B及3B染色體上可能存在控制千粒質(zhì)量相關(guān)的基因,但是3B染色體上檢測得到的千粒質(zhì)量QTL貢獻(xiàn)率也較低。此外,本研究還在5A和7B染色體上分別檢測到3個(gè)和4個(gè)千粒質(zhì)量QTL,這與前人在相應(yīng)染色體上得到的千粒質(zhì)量QTL的位置有所不同,是否為新的千粒質(zhì)量QTL位點(diǎn)尚需進(jìn)一步分析。本研究只在粒長性狀上檢測到1個(gè)能在2環(huán)境條件下穩(wěn)定表達(dá)的QTL;2個(gè)環(huán)境條件下檢測得到的粒長QTL及千粒質(zhì)量QTL結(jié)果差異較大,這可能是受2地環(huán)境條件影響所導(dǎo)致;李巧云等[26]研究曾指出小麥產(chǎn)量性狀易受地理環(huán)境和生產(chǎn)條件的影響,同一品種在不同地區(qū)或同一地區(qū)不同種植地點(diǎn)都可能存在明顯差異。

    本研究所采用的2個(gè)親本‘西農(nóng)981’和‘陜麥159’是近年來審定的農(nóng)藝性狀優(yōu)良的小麥品種,本研究結(jié)果對(duì)利用這2個(gè)親本進(jìn)行新品種的選育具有一定的借鑒;在后續(xù)的研究中將重點(diǎn)對(duì)2B及3B染色體區(qū)段進(jìn)行加密研究,為利用分子標(biāo)記輔助育種及培育具有高產(chǎn)、優(yōu)質(zhì)的品種提供有用信息。

    Reference:

    [1] GIURA A,SAULESCU N N.Chromosomal location of genes controlling grain size in a large grained selection of wheat (TriticumaestivumL.)[J].Euphytica,1996,89(1):77-80.

    [2] 王瑞霞,張秀英,伍 玲,等.不同生態(tài)環(huán)境下冬小麥籽粒大小相關(guān)性狀的QTL分析[J].中國農(nóng)業(yè)科學(xué),2009,42(2):398-407.

    WANG R X,ZHANG X Y,WU L,etal.QTL analysis of grain size and related traits in winter wheat under different ecological environments.[J].ScientiaAgriculturaSinica,2009,42(2):398-407(in Chinese with English abstract)

    [3] CAMPBELL K G,BERGMAN C J,GUALBERTO D G,etal.Quantitative trait loci associated with kernel traits in a soft × hard wheat cross[J].CropScience,1999,39(4):1184-1195.

    [4] VARSHNEY R K,PRASAD M,ROY J K,etal.Identification of eight chromosomes and a microsatellite marker on 1AS associated with QTL for grain weight in bread wheat[J].Theoretical&AppliedGenetics,2000,100(8):1290-1294.

    [5] AMMIRAJU J S S,DHOLAKIA B B,SANTRA D K,etal.Identification of inter simple sequence repeat (ISSR) markers associated with seed size in wheat[J].Theoretical&AppliedGenetics,2001,102(5):726-732.

    [6] 李斯深.小麥產(chǎn)量性狀QTL的分子標(biāo)記定位[D].濟(jì)南:山東農(nóng)業(yè)大學(xué),2002.

    LI S SH.QTL mapping for yield traits using molecular markers in wheat[D].Jinan:Shandong Agricultural University,2002(in Chinese).

    [7] 宋彥霞,景蕊蓮,霍納新,等.普通小麥(T.aestivumL.)不同作圖群體抽穗期QTL分析[J].中國農(nóng)業(yè)科學(xué),2006,39(11):2186-2193.

    SONG Y X,JING R L,HUO N X,etal.Detection of QTLs for heading in common wheat (T.aestivumL.) using different populations[J].ScientiaAgriculturaSinica,2006,39(11):2186-2193 (in Chinese with English abstract).

    [8] TYAGI S,MIR R R,BALYAN H S,etal.Interval mapping and meta-QTL analysis of grain traits in common wheat (TriticumaestivumL.)[J].Euphytica,2015,201(3):367-380.

    [9] LI C,BAI G,CARVER B F,etal.Single nucleotide polymorphism markers linked to QTL for wheat yield traits[J].Euphytica,2015,206(1):89-101.

    [10] 姚 琴,周榮華,潘昱名,等.小麥品種偃展1號(hào)與品系早穗30重組自交系群體遺傳連鎖圖譜構(gòu)建及重要農(nóng)藝性狀的QTL分析[J].中國農(nóng)業(yè)科學(xué),2010,43(20):4130-4139.

    YAO Q,ZHOU R H,PAN Y M,etal.Construction of genetic linkage map and QTL analysis of agronomic important traits based on a RIL population derived from common wheat variety Yanzhan 1 and Zaosui 30.[J].ScientiaAgriculturaSinica,2010,43(20):4130-4139(in Chinese with English abstract).

    [11] LI G,QUIROS C F.Sequence-related amplified polymorphism (SRAP),a new marker system based on a simple PCR reaction:its application to mapping and gene tagging in Brassica[J].Theoretical&AppliedGenetics,2001,103:455-461.

    [12] 祁建民,梁景霞,陳美霞,等.應(yīng)用ISSR與SRAP分析煙草種質(zhì)資源遺傳多樣性及遺傳演化關(guān)系[J].作物學(xué)報(bào),2012,38(8):1425-1434.

    QI J M,LIANG J X,CHEN M X,etal.Genetic diversity and evolutionary analysis of tobacco (NicotianatabacumL.)germplasm resources based on ISSR and SRAP markers[J].ActaAgronomicaSinica,2012,38(8):1425-1434(in Chinese with English abstract).

    [13] 韓 芳,亓佳佳,馬守才,等.黃淮麥區(qū)部分小麥品種(系)遺傳多樣性的SRAP分析[J].西北農(nóng)業(yè)學(xué)報(bào),2014,23(12):60-67.

    HAN F,QI J J,MA SH C,etal.Genetic diversity of some wheat varieties in Huang-Huai wheat area revealed by SRAP marker[J].ActaAgriculturaeBoreali-occidentalisSinica,2014,23(12):60-67(in Chinese with English abstract).

    [14] SAGHAI-MAROOF M A,SOLIMAN K M,JORGENSEN R A,etal.Ribosomal DNA spacer-length polymorphisms in barley:mendelian inheritance,chromosomal location,and population dynamics[J].ProcedingoftheNationalAcademyofSciencesoftheUnitedStatesofAmerica,1984,81(24):8014-8018.

    [15] JIANG Y F,LAN X J,LUO W,etal.Genome-wide quantitative trait locus mapping identifies multiple major loci for brittle rachis and threshability in Tibetan semi-wild wheat (Triticumaestivumssp.tibetanumShao)[J].PlosOne,2013,9(12):114066-114066.

    [16] 程小毛,黃曉霞.SSR標(biāo)記開發(fā)及其在植物中的應(yīng)用[J].中國農(nóng)學(xué)通報(bào),2011,27(5):304-307.

    CHENG X M,HUANG X X.Development and application of SSR markers in plants[J].ChineseAgriculturalScienceBulletin,2011,27(5):304-307(in Chinese with English abstract).

    [17] 李媛媛,沈金雄,王同華,等.利用SRAP、SSR和AFLP標(biāo)記構(gòu)建甘藍(lán)型油菜遺傳連鎖圖譜[J].中國農(nóng)業(yè)科學(xué),2007,40(6):1118-1126.

    LI Y Y,SHEN J X,WANG T H,etal.Construction of a linkage map using SRAP,SSR and AFLP markers inBrassicanapusL.[J].ScientiaAgriculturaSinica,2007,40(6):1118-1126(in Chinese with English abstract).

    [18] 賈繼增,張正斌.小麥21條染色體RFLP作圖位點(diǎn)遺傳多樣性分析[J].中國科學(xué)(C輯:生命科學(xué)),2001,31(1):13-21.

    JIA J Z,ZHANG ZH B.Genetic diversity of 21 chromosome for wheat by RFLP marker[J].ScienceinChina(SeriesC),2001,31(1):13-21(in Chinese).

    [19] LI Q,ZHANG Y,LIU T,etal.Genetic analysis of kernel weight and kernel size in wheat (TriticumaestivumL.) using unconditional and conditional QTL mapping[J].MolecularBreeding,2015,35(10):1-15.

    [20] 余曼麗,趙林姝,郭會(huì)君,等.小麥籽粒性狀的 QTL 定位[J].麥類作物學(xué)報(bào),2014,34(8):1029-1035.

    YU M L,ZHAO L SH,GUO H J,etal.QTL mapping for kernel traits in wheat[J].JournalofTriticeaeCrops,2014,34(8):1029-1035(in Chinese with English abstract).

    [21] RAMYA P,CHAUBAL A,KULKARNI K,etal.QTL mapping of 1 000-kernel weight,kernel length,and kernel width in bread wheat (TriticumaestivumL.)[J].JournalofAppliedGenetics,2010,51(4):421-429.

    [22] HUANG X Q,KEMPF H,GANAL M W,etal.Advanced backcross QTL analysis in progenies derived from a cross between a German elite winter wheat variety and a synthetic wheat (TriticumaestivumL.)[J].TheoreticalandAppliedGenetics,2004,109(5):933-943.

    [23] HUANG X Q,CLOUTIER S,LYCAR L,etal.Molecular detection of QTLs for agronomic and quality traits in a doubled haploid population derived from two Canadian wheats (TriticumaestivumL.)[J].Theoretical&AppliedGenetics,2006,113(4):753-766.

    [24] 周淼平,任麗娟,張 旭,等.小麥產(chǎn)量性狀的 QTL 分析[J].麥類作物學(xué)報(bào),2006,26(4):35-40.

    ZHOU M P,RENG L J,ZHANG X,etal.Analysis of QTLs for yield traits of wheat [J].JournalofTriticeaeCrops,2006,26(4):35-40(in Chinese with English abstract).

    [25] 丁安明,崔 法,李 君,等.小麥單株產(chǎn)量與株高的 QTL 分析[J].中國農(nóng)業(yè)科學(xué),2011,44(14):2857-2867.

    DING A M,CUI F,LI J,etal.QTL analysis on grain yield per plant and plant height in wheat[J].ScientiaAgriculturaSinica,2011,44(14):2857-2867(in Chinese with English abstract).

    [26] 李巧云,李 磊,劉萬代,等.河南省小麥產(chǎn)量及其構(gòu)成因素變化規(guī)律分析[J].河南農(nóng)業(yè)科學(xué),2011,40(4):38-40.

    LI Q Y,LI L,LIU W D,etal.Analyses of yield and its composition property of wheat in henan province [J].JournalofHenanAgriculturalSciences,2011,40(4):38-40(in Chinese with English abstract).

    (責(zé)任編輯:史亞歌 Responsible editor:SHI Yage)

    The QTL Mapping and Effect Analysis of Wheat Kernel Length and 1 000-grain Mass Based on SRAP and SSR Markers

    GUO Lijian1, WANG Zhulin2, WANG Shijuan1, LIU Xiangli1and ZHAO Huixian1

    (1.College of Life Sciences,Northwest A&F University, Yangling Shaanxi 712100, China;2. College of Agronomy, Northwest A&F University, Yangling Shaanxi 712100, China)

    In order to explore QTLs related to kernel length and 1 000-grain mass, wheat variety ‘Xinong 981’ and ‘Shaanmai 159’ were used to develop 169 F2population and F2:3lines. The F2population was exploited to construct genetic map and the F2:3lines grown in two different environments of Yangling and Sanyuan for testing the phenotypic data. The results showed that highly significant difference of environment variance was detected for kernel length and 1 000-grain mass.Thirty-one related QTLs were detected in two environments. Of all, seven QTLs account for 4.36%-14.80% phenotypic variance of wheat kernel length, scattered in 5A, 6A, 7A, 2B, 3B, and 4B chromosomes, and on 3B chromosome, one kernel length-related QTL was detected in both locations; 24 QTLs for 1 000-grain mass, distributed on 2A, 3A, 5A, 6A, 7A, 1B, 2B, 3B, 4B,5B, and 7B chromosome and can explain 0.25%-8.23% phenotypic variance. Five QTLs related to 1 000-grain mass were detected on 2B chromosome, suggesting that 2B chromosome might closely related to kernel mass.

    Wheat; Kernel length; Grain mass; SRAP; QTL

    2016-04-14 Returned 2016-06-02

    The National Natural Science Foundation of China (No.31471482).

    GUO Lijian, male,master student. Research area:wheat QTL reseach. E-mail:baiheichaye@163.com

    ZHAO Huixian, female, doctoral supervisor. Research area:biochemistry and molecular biology. E-mail:hxzhao212@nwsuaf.edu.cn

    日期:2017-08-18

    2016-04-14

    2016-06-02

    國家自然科學(xué)基金(31471482)。 第一作者:郭利建,男,碩士研究生,從事小麥Q(jìng)TL定位研究。E-mail:baiheichaye@163.com

    趙惠賢,女,博士,教授,主要從事生化與分子生物學(xué)研究。E-mail:hxzhao212@nwsuaf.edu.cn

    S512.1

    A

    1004-1389(2017)08-1165-08

    網(wǎng)絡(luò)出版地址:http://kns.cnki.net/kcms/detail/61.1220.S.20170818.0938.020.html

    猜你喜歡
    粒長染色體遺傳
    非遺傳承
    水稻粒長遺傳及其功能基因研究進(jìn)展
    還有什么會(huì)遺傳?
    還有什么會(huì)遺傳
    還有什么會(huì)遺傳?
    多一條X染色體,壽命會(huì)更長
    為什么男性要有一條X染色體?
    能忍的人壽命長
    秈稻粒長與稻米品質(zhì)的相關(guān)性及其育種應(yīng)用
    日本晴/R1126水稻重組自交系群體粒形性狀QTL定位
    一级a做视频免费观看| 网址你懂的国产日韩在线| 精品久久久久久成人av| 国产黄色视频一区二区在线观看| 亚洲欧美日韩卡通动漫| 热99在线观看视频| 日韩欧美 国产精品| 人妻制服诱惑在线中文字幕| 久热久热在线精品观看| 亚洲国产欧美在线一区| av国产免费在线观看| 永久网站在线| 久久鲁丝午夜福利片| 午夜精品一区二区三区免费看| 久久国内精品自在自线图片| 91久久精品国产一区二区三区| 五月玫瑰六月丁香| 日韩人妻高清精品专区| 成年女人看的毛片在线观看| 久久久久久久久久人人人人人人| 久久久久久久久中文| 老女人水多毛片| 亚洲久久久久久中文字幕| 少妇猛男粗大的猛烈进出视频 | 精品少妇黑人巨大在线播放| 国产成人精品一,二区| 简卡轻食公司| 日本免费a在线| 精品一区在线观看国产| 国产av不卡久久| 日韩不卡一区二区三区视频在线| 精品不卡国产一区二区三区| 国产午夜精品一二区理论片| 成人二区视频| 久久久久九九精品影院| 一级二级三级毛片免费看| 日韩三级伦理在线观看| 最近中文字幕2019免费版| 大又大粗又爽又黄少妇毛片口| 国产精品一区二区三区四区久久| 日韩人妻高清精品专区| 久久久久久久久久久免费av| 午夜精品一区二区三区免费看| 国产69精品久久久久777片| 一区二区三区乱码不卡18| or卡值多少钱| 日韩制服骚丝袜av| 少妇人妻一区二区三区视频| 亚洲在久久综合| 久久久久久久亚洲中文字幕| 十八禁国产超污无遮挡网站| 男人和女人高潮做爰伦理| 乱人视频在线观看| 欧美另类一区| 免费av毛片视频| 亚洲国产精品sss在线观看| 国产真实伦视频高清在线观看| 免费高清在线观看视频在线观看| 人妻夜夜爽99麻豆av| 久久久久国产网址| av在线亚洲专区| 欧美人与善性xxx| av黄色大香蕉| 99热全是精品| 中国国产av一级| 精品午夜福利在线看| 亚洲内射少妇av| 有码 亚洲区| 麻豆国产97在线/欧美| 少妇人妻一区二区三区视频| 午夜福利高清视频| 亚洲av二区三区四区| 国产伦精品一区二区三区四那| 免费人成在线观看视频色| 秋霞伦理黄片| 亚洲性久久影院| 精品久久久久久久人妻蜜臀av| 蜜桃久久精品国产亚洲av| 亚洲精品中文字幕在线视频 | 午夜福利在线在线| 人妻夜夜爽99麻豆av| 成人一区二区视频在线观看| 国产片特级美女逼逼视频| 18禁裸乳无遮挡免费网站照片| 能在线免费看毛片的网站| 国产成人免费观看mmmm| 精品熟女少妇av免费看| 亚洲最大成人中文| 免费看不卡的av| 久久久久久久亚洲中文字幕| 免费看日本二区| 精品人妻熟女av久视频| 国产成人aa在线观看| 一级毛片久久久久久久久女| 欧美bdsm另类| 99视频精品全部免费 在线| 尾随美女入室| 80岁老熟妇乱子伦牲交| 亚洲在线观看片| 久99久视频精品免费| 成人无遮挡网站| 亚洲av一区综合| 成人漫画全彩无遮挡| 久久99热这里只有精品18| 深夜a级毛片| 97人妻精品一区二区三区麻豆| 少妇丰满av| 国产永久视频网站| 三级经典国产精品| 亚洲经典国产精华液单| 国产熟女欧美一区二区| 成人毛片a级毛片在线播放| 久久久午夜欧美精品| 两个人的视频大全免费| 国产免费福利视频在线观看| 2022亚洲国产成人精品| 美女主播在线视频| 最近最新中文字幕免费大全7| 3wmmmm亚洲av在线观看| 亚洲欧美成人精品一区二区| 国产综合懂色| 老女人水多毛片| 久久精品熟女亚洲av麻豆精品 | 啦啦啦韩国在线观看视频| 最近最新中文字幕大全电影3| 99热这里只有精品一区| 亚洲在线观看片| 秋霞伦理黄片| 18+在线观看网站| .国产精品久久| 亚洲精品视频女| 纵有疾风起免费观看全集完整版 | 少妇丰满av| 亚洲精品日韩av片在线观看| 国产精品久久久久久久电影| 中文字幕av成人在线电影| 草草在线视频免费看| 国产精品久久久久久精品电影| 精品人妻偷拍中文字幕| 波野结衣二区三区在线| 亚洲电影在线观看av| 国产黄色免费在线视频| 男女啪啪激烈高潮av片| 国产一级毛片七仙女欲春2| 精品不卡国产一区二区三区| 亚洲av国产av综合av卡| 国产精品av视频在线免费观看| 热99在线观看视频| 又大又黄又爽视频免费| 在线观看一区二区三区| 国产伦精品一区二区三区视频9| 亚洲国产精品专区欧美| 精品国产三级普通话版| 在线免费观看的www视频| 久久99热这里只有精品18| 国产午夜精品久久久久久一区二区三区| 日日啪夜夜撸| av播播在线观看一区| 亚洲自拍偷在线| 亚洲乱码一区二区免费版| 日韩欧美精品v在线| 国产一区二区三区av在线| 91aial.com中文字幕在线观看| 国产三级在线视频| 欧美日韩在线观看h| 精品人妻视频免费看| 日本午夜av视频| 少妇被粗大猛烈的视频| 男人舔女人下体高潮全视频| 男人舔女人下体高潮全视频| 国产精品av视频在线免费观看| 精品久久久久久电影网| 免费电影在线观看免费观看| 亚洲av男天堂| 亚洲精品乱久久久久久| 一级片'在线观看视频| a级毛片免费高清观看在线播放| 久久久国产一区二区| 国产亚洲最大av| 一个人免费在线观看电影| 春色校园在线视频观看| 又爽又黄无遮挡网站| 日韩 亚洲 欧美在线| 国产精品嫩草影院av在线观看| 日本欧美国产在线视频| 18禁动态无遮挡网站| 大又大粗又爽又黄少妇毛片口| 日本熟妇午夜| 99re6热这里在线精品视频| 国产黄片美女视频| ponron亚洲| 精品国产三级普通话版| 免费av毛片视频| 久久这里只有精品中国| 久久99热这里只有精品18| 欧美另类一区| 久久97久久精品| 亚洲欧美中文字幕日韩二区| 51国产日韩欧美| 插阴视频在线观看视频| 国产精品一区www在线观看| 国产精品一二三区在线看| av在线老鸭窝| 国产亚洲精品av在线| 国产免费福利视频在线观看| 国产一区二区三区av在线| 九色成人免费人妻av| 色播亚洲综合网| 欧美xxxx性猛交bbbb| 免费在线观看成人毛片| 内地一区二区视频在线| 午夜免费男女啪啪视频观看| 中国美白少妇内射xxxbb| 久久久久久久久久久免费av| 人人妻人人澡欧美一区二区| 国产av不卡久久| 免费观看无遮挡的男女| 午夜福利视频1000在线观看| 国产精品久久视频播放| 午夜精品在线福利| 一级av片app| 中文字幕av在线有码专区| 一个人看的www免费观看视频| 人体艺术视频欧美日本| av.在线天堂| 国产三级在线视频| 最近最新中文字幕免费大全7| 搡老乐熟女国产| 六月丁香七月| 精品人妻熟女av久视频| 国产在线一区二区三区精| 欧美高清成人免费视频www| 高清午夜精品一区二区三区| 男的添女的下面高潮视频| 日产精品乱码卡一卡2卡三| 久久精品夜夜夜夜夜久久蜜豆| 国产男人的电影天堂91| 欧美一区二区亚洲| 午夜福利在线观看吧| 人妻少妇偷人精品九色| 国产精品伦人一区二区| 国产色爽女视频免费观看| 日韩在线高清观看一区二区三区| 九草在线视频观看| 两个人的视频大全免费| av免费观看日本| 国产精品一区二区三区四区久久| 18禁在线播放成人免费| 一级毛片我不卡| 国产视频内射| 真实男女啪啪啪动态图| 亚洲无线观看免费| 大陆偷拍与自拍| 久久久亚洲精品成人影院| 一级毛片 在线播放| 一级二级三级毛片免费看| 日日干狠狠操夜夜爽| 在线 av 中文字幕| 中文字幕免费在线视频6| 麻豆成人av视频| 国产精品国产三级专区第一集| 国产精品国产三级专区第一集| 亚洲自拍偷在线| 国产成人91sexporn| 国产精品精品国产色婷婷| 国产成人aa在线观看| 一级爰片在线观看| 欧美激情久久久久久爽电影| 一级毛片黄色毛片免费观看视频| 亚洲欧美精品专区久久| 舔av片在线| 久久精品夜色国产| 青春草视频在线免费观看| 国产伦在线观看视频一区| 免费无遮挡裸体视频| 人人妻人人澡人人爽人人夜夜 | 午夜福利在线观看免费完整高清在| 国产精品av视频在线免费观看| 卡戴珊不雅视频在线播放| 成人二区视频| 丝袜喷水一区| 高清视频免费观看一区二区 | 成人特级av手机在线观看| 国产亚洲av嫩草精品影院| 青春草国产在线视频| 中文字幕av在线有码专区| 在线观看av片永久免费下载| 精品久久国产蜜桃| 69人妻影院| 午夜激情欧美在线| 亚洲丝袜综合中文字幕| 午夜福利在线观看免费完整高清在| 黄片无遮挡物在线观看| 国产成人福利小说| 国产成人aa在线观看| 成人午夜精彩视频在线观看| 国产永久视频网站| 美女国产视频在线观看| 日日摸夜夜添夜夜爱| 日本一二三区视频观看| 色视频www国产| 精品久久久久久久人妻蜜臀av| 寂寞人妻少妇视频99o| 日韩三级伦理在线观看| 嫩草影院新地址| 精品欧美国产一区二区三| 极品教师在线视频| 成人午夜精彩视频在线观看| 不卡视频在线观看欧美| 永久免费av网站大全| 美女黄网站色视频| 久久99热这里只频精品6学生| 三级国产精品片| 少妇熟女欧美另类| 久久久久久久亚洲中文字幕| 国内精品一区二区在线观看| 国产免费一级a男人的天堂| 菩萨蛮人人尽说江南好唐韦庄| 国产高清国产精品国产三级 | 久久鲁丝午夜福利片| 人人妻人人澡欧美一区二区| 麻豆国产97在线/欧美| 高清欧美精品videossex| 熟妇人妻不卡中文字幕| 熟妇人妻久久中文字幕3abv| 亚洲精品自拍成人| 777米奇影视久久| 欧美xxxx黑人xx丫x性爽| 男人狂女人下面高潮的视频| 国产又色又爽无遮挡免| 边亲边吃奶的免费视频| 亚洲av免费在线观看| 亚洲精品国产av成人精品| 免费黄网站久久成人精品| 欧美日本视频| 亚洲精品,欧美精品| 亚洲精品久久午夜乱码| 午夜免费激情av| 国产单亲对白刺激| 91精品国产九色| 国产久久久一区二区三区| 黄色一级大片看看| 久久精品国产自在天天线| 丰满少妇做爰视频| 观看免费一级毛片| 日韩欧美 国产精品| 久久久久久久久久黄片| 性插视频无遮挡在线免费观看| 伊人久久国产一区二区| 国产午夜福利久久久久久| 国产一区二区亚洲精品在线观看| 国产亚洲一区二区精品| 日韩国内少妇激情av| freevideosex欧美| 国产精品一及| 免费观看无遮挡的男女| 亚洲精品中文字幕在线视频 | 国产高清国产精品国产三级 | 九草在线视频观看| 永久免费av网站大全| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品久久午夜乱码| 99久久人妻综合| 大话2 男鬼变身卡| 搡老乐熟女国产| 国产精品久久久久久久久免| 18+在线观看网站| 波野结衣二区三区在线| 精品久久久久久久久av| 色视频www国产| 成人亚洲精品一区在线观看 | 亚洲国产欧美人成| 狂野欧美白嫩少妇大欣赏| 伦精品一区二区三区| 搡女人真爽免费视频火全软件| videossex国产| 尤物成人国产欧美一区二区三区| 亚洲,欧美,日韩| 亚洲乱码一区二区免费版| 日本爱情动作片www.在线观看| 色尼玛亚洲综合影院| 亚洲欧美精品自产自拍| 偷拍熟女少妇极品色| 最近的中文字幕免费完整| 一级毛片aaaaaa免费看小| 国产一级毛片七仙女欲春2| 日本wwww免费看| 一本一本综合久久| 最近最新中文字幕大全电影3| 大陆偷拍与自拍| 国产亚洲午夜精品一区二区久久 | 在线观看免费高清a一片| 欧美极品一区二区三区四区| 三级经典国产精品| 久久久久久久久大av| 99久久精品国产国产毛片| 精品久久久噜噜| 22中文网久久字幕| 亚洲最大成人av| 99热全是精品| 亚洲av中文字字幕乱码综合| 日产精品乱码卡一卡2卡三| 毛片一级片免费看久久久久| 又黄又爽又刺激的免费视频.| 午夜日本视频在线| 国产高清国产精品国产三级 | 亚洲经典国产精华液单| 国产成人精品婷婷| 国产精品久久久久久久久免| 国产免费福利视频在线观看| 国产精品女同一区二区软件| 熟妇人妻不卡中文字幕| 成年版毛片免费区| 爱豆传媒免费全集在线观看| 久久99热6这里只有精品| 在线a可以看的网站| 舔av片在线| 在线观看人妻少妇| 日日啪夜夜撸| 国产亚洲最大av| 国产精品久久久久久精品电影小说 | 日本三级黄在线观看| av在线蜜桃| 久久综合国产亚洲精品| 欧美 日韩 精品 国产| 搞女人的毛片| 国产一区二区在线观看日韩| 久久97久久精品| 国产一级毛片在线| 国产精品一区二区在线观看99 | 久久99精品国语久久久| 久久久久九九精品影院| 久久亚洲国产成人精品v| 日日撸夜夜添| 国产亚洲一区二区精品| 久久亚洲国产成人精品v| 少妇丰满av| 亚洲av免费高清在线观看| 久久精品国产亚洲av涩爱| 最新中文字幕久久久久| 99久国产av精品国产电影| 国产亚洲精品av在线| 爱豆传媒免费全集在线观看| 国产v大片淫在线免费观看| 国产精品一区二区性色av| 国产免费视频播放在线视频 | 亚洲欧洲国产日韩| 一级毛片久久久久久久久女| 国产一级毛片在线| 夜夜看夜夜爽夜夜摸| 一级二级三级毛片免费看| 精品久久久久久久久久久久久| 日本与韩国留学比较| 啦啦啦中文免费视频观看日本| 国产久久久一区二区三区| 高清午夜精品一区二区三区| 中国国产av一级| 亚洲天堂国产精品一区在线| freevideosex欧美| 丰满少妇做爰视频| 在线播放无遮挡| www.色视频.com| 久久精品国产亚洲网站| 国精品久久久久久国模美| 久久草成人影院| 精品熟女少妇av免费看| 欧美性感艳星| 中文字幕人妻熟人妻熟丝袜美| 天天躁夜夜躁狠狠久久av| av在线播放精品| 夜夜爽夜夜爽视频| 国产 一区精品| 大香蕉97超碰在线| 黄色配什么色好看| 少妇的逼好多水| 别揉我奶头 嗯啊视频| 蜜臀久久99精品久久宅男| 国产在视频线精品| 亚洲,欧美,日韩| 国产精品三级大全| videos熟女内射| 五月伊人婷婷丁香| 99热全是精品| 永久免费av网站大全| 久久6这里有精品| 国产激情偷乱视频一区二区| 国产精品av视频在线免费观看| 少妇人妻精品综合一区二区| 成人亚洲欧美一区二区av| 一夜夜www| 国产精品不卡视频一区二区| 亚洲av福利一区| 伊人久久国产一区二区| 亚洲成人精品中文字幕电影| 五月天丁香电影| 麻豆久久精品国产亚洲av| 欧美成人一区二区免费高清观看| 久久久久国产网址| 色综合站精品国产| 国产色爽女视频免费观看| 国内揄拍国产精品人妻在线| 亚洲一区高清亚洲精品| 日韩一本色道免费dvd| 日本欧美国产在线视频| .国产精品久久| 视频中文字幕在线观看| 小蜜桃在线观看免费完整版高清| 最近最新中文字幕免费大全7| 国产69精品久久久久777片| 亚洲电影在线观看av| 91精品伊人久久大香线蕉| 高清视频免费观看一区二区 | 国产在视频线在精品| 国产不卡一卡二| 欧美丝袜亚洲另类| 成人毛片60女人毛片免费| 午夜福利在线观看吧| 免费大片黄手机在线观看| 听说在线观看完整版免费高清| 久久精品熟女亚洲av麻豆精品 | 国产高清有码在线观看视频| av在线亚洲专区| av一本久久久久| 中文字幕久久专区| 亚州av有码| av在线老鸭窝| 我的女老师完整版在线观看| 熟妇人妻不卡中文字幕| av国产免费在线观看| 日本黄色片子视频| 久久久午夜欧美精品| 成人av在线播放网站| 精品熟女少妇av免费看| 国语对白做爰xxxⅹ性视频网站| 国产精品一区二区在线观看99 | 人妻少妇偷人精品九色| 亚洲第一区二区三区不卡| videossex国产| av天堂中文字幕网| 亚洲精品影视一区二区三区av| 国内精品一区二区在线观看| 最近中文字幕高清免费大全6| 亚洲国产日韩欧美精品在线观看| 一级a做视频免费观看| 国产一区有黄有色的免费视频 | 日本免费a在线| 伦精品一区二区三区| 久久久精品欧美日韩精品| 国产午夜精品一二区理论片| 国产国拍精品亚洲av在线观看| 黄色欧美视频在线观看| 午夜激情久久久久久久| 中文在线观看免费www的网站| 久久精品久久久久久噜噜老黄| 国产精品伦人一区二区| 亚洲自拍偷在线| 欧美3d第一页| 亚洲成人精品中文字幕电影| 国产视频内射| 免费高清在线观看视频在线观看| 在线免费观看不下载黄p国产| 亚洲精品aⅴ在线观看| 国产亚洲精品久久久com| 你懂的网址亚洲精品在线观看| 校园人妻丝袜中文字幕| 日韩国内少妇激情av| 日本爱情动作片www.在线观看| 国产伦在线观看视频一区| 狂野欧美激情性xxxx在线观看| av专区在线播放| 日韩av免费高清视频| 欧美一级a爱片免费观看看| 亚洲一级一片aⅴ在线观看| 亚洲av在线观看美女高潮| 成人国产麻豆网| 成年av动漫网址| 又大又黄又爽视频免费| 成人亚洲精品av一区二区| eeuss影院久久| 中文字幕亚洲精品专区| 91av网一区二区| 晚上一个人看的免费电影| 婷婷色综合大香蕉| 高清毛片免费看| 日韩三级伦理在线观看| 真实男女啪啪啪动态图| av又黄又爽大尺度在线免费看| 亚洲三级黄色毛片| 99热网站在线观看| 中文字幕人妻熟人妻熟丝袜美| 久久久久久伊人网av| av黄色大香蕉| 国产 亚洲一区二区三区 | 男人狂女人下面高潮的视频| 国产一区二区在线观看日韩| 内地一区二区视频在线| 国产成人91sexporn| 日本一本二区三区精品| 2018国产大陆天天弄谢| 午夜视频国产福利| 搡老乐熟女国产| 国产精品99久久久久久久久| 久久精品久久精品一区二区三区| 大香蕉久久网| 精品久久久精品久久久| 亚洲欧美精品自产自拍| 免费黄色在线免费观看| 欧美日韩在线观看h| 狠狠精品人妻久久久久久综合| 菩萨蛮人人尽说江南好唐韦庄| 天堂影院成人在线观看| 爱豆传媒免费全集在线观看| 一级爰片在线观看| 床上黄色一级片| 小蜜桃在线观看免费完整版高清|