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

    西北地區(qū)和尚頭小麥遺傳多樣性及農(nóng)藝性狀的關(guān)聯(lián)分析

    2019-02-25 01:57:12,,,,
    草業(yè)學(xué)報(bào) 2019年2期
    關(guān)鍵詞:表型種質(zhì)遺傳

    ,,,,

    (甘肅省農(nóng)業(yè)科學(xué)院作物研究所,甘肅 蘭州 730070)

    小麥(Triticumaestivum)是我國主要的糧食作物之一,尤其是北方地區(qū)重要的食物原材料,在我國已有幾千年的種植歷史。隨著全球氣候變暖、災(zāi)害頻發(fā)、耕作栽培制度變革和人民生活水平的提高,小麥品種的多樣化得到世界各國的高度重視。但是,在小麥育種快速發(fā)展的今天,我國小麥優(yōu)良親本材料遺傳基礎(chǔ)變得越來越狹窄,嚴(yán)重制約了小麥品種的遺傳改良,導(dǎo)致品種的同質(zhì)化問題比較突出。古老地方品種經(jīng)過長期栽培馴化和自然選擇,往往具有抗逆性強(qiáng)、品質(zhì)優(yōu)等特點(diǎn),蘊(yùn)藏著豐富的優(yōu)異基因,是小麥遺傳改良的重要資源。和尚頭小麥?zhǔn)俏鞅钡貐^(qū)的地方品種,且種植歷史悠久,經(jīng)調(diào)查發(fā)現(xiàn)各地被稱作“和尚頭”的小麥有多種多樣,同名異質(zhì)、同質(zhì)異名現(xiàn)象普遍存在,各自的特征特性不是十分清楚,而有關(guān)其遺傳多樣性的研究又鮮見報(bào)道,不利于育種家選擇利用[1]。因此,分析西北地區(qū)和尚頭小麥種質(zhì)資源的遺傳多樣性,挖掘與主要農(nóng)藝性狀相關(guān)的優(yōu)異等位位點(diǎn),區(qū)分同名不同來源的和尚頭小麥,明確其性狀特點(diǎn),有利于提高育種利用效率。

    遺傳資源多樣性的研究是作物育種的前提[2]。早期,作物遺傳多樣性分析主要依靠形態(tài)學(xué)標(biāo)記。不同學(xué)者利用形態(tài)學(xué)指標(biāo)對小麥遺傳多樣性分析發(fā)現(xiàn),育成品種的遺傳多樣性較地方品種有所下降,品種遺傳基礎(chǔ)日益狹窄[3-7]。因此,小麥地方品種資源越來越受到廣大育種家的重視。形態(tài)學(xué)指標(biāo)具有直觀、易于識別,便于掌握等特點(diǎn),但受環(huán)境影響大,準(zhǔn)確性較差。隨著人們對基因的深入研究,分子標(biāo)記得到迅速發(fā)展,尤其SSR分子標(biāo)記由于多態(tài)性好、穩(wěn)定性可靠、經(jīng)濟(jì)方便使其被廣泛應(yīng)用[8-11],利用SSR分子標(biāo)記分析小麥和小麥野生近緣種遺傳多樣性研究報(bào)道較多[12-17],技術(shù)已非常成熟。關(guān)聯(lián)分析是發(fā)掘基因和等位基因的有效方法,對分子標(biāo)記輔助選擇育種具有非常重要的意義,其已被廣泛應(yīng)用于作物種質(zhì)資源分子評價(jià),發(fā)掘了一批優(yōu)異的等位基因[18-19]。根據(jù)本研究前期對甘肅和尚頭小麥的實(shí)地調(diào)查發(fā)現(xiàn)[1],和尚頭具有耐瘠薄、抗逆性強(qiáng)等優(yōu)良性狀,其面粉質(zhì)量好,尤其是蛋白質(zhì)含量高,具有滑潤爽口、味感純正、面筋強(qiáng)、食用方便等特點(diǎn),但抗倒伏性差、產(chǎn)量相對較低。由于小麥1B染色體的短臂被黑麥(Secalecereale)的1R染色體短臂所取代即形成小麥-黑麥1BL/1RS易位系,從而對小麥的品質(zhì)產(chǎn)生影響,馬小樂等[20]在和尚頭小麥中未檢測出影響小麥品質(zhì)的1BL/1RS易位系,其品質(zhì)達(dá)到高筋品種標(biāo)準(zhǔn);王世紅[21]在Glu-B1(Glutenin-B1)位點(diǎn)檢測到小麥品質(zhì)優(yōu)質(zhì)亞基7+8,說明和尚頭小麥品質(zhì)較好。和尚頭小麥有較強(qiáng)的環(huán)境變化耐受性,在改變水分和氮素含量時(shí)氣孔導(dǎo)度和蒸騰速率變化幅度較小,水分利用效率高,能夠維持產(chǎn)量構(gòu)成因素的恒定,提高產(chǎn)量[22-24];高天鵬等[25]研究表明,增強(qiáng)UV-B(Ultraviolet-B)輻射和干旱脅迫下,和尚頭小麥生物量較其他品種降低幅度較小。可見,前人多以某一地區(qū)的和尚頭小麥為研究對象,并未對西北地區(qū)的同名和尚頭小麥進(jìn)行全面系統(tǒng)分析。本研究以來自西北地區(qū)的43份和尚頭小麥為研究對象,采用表型和分子標(biāo)記相結(jié)合的方法,分析供試材料的遺傳多樣性,利用關(guān)聯(lián)分析挖掘與主要農(nóng)藝性狀相關(guān)聯(lián)的優(yōu)異等位位點(diǎn)。闡明同名不同來源和尚頭小麥的特征特性和應(yīng)用價(jià)值,為育種家選擇利用和分子標(biāo)記輔助育種提供參考依據(jù)。

    1 材料與方法

    1.1 試驗(yàn)材料

    選用西北地區(qū)和尚頭小麥種質(zhì)資源共計(jì)43份 (表1),分別由國家農(nóng)作物種質(zhì)資源共享服務(wù)平臺、國家農(nóng)作物種質(zhì)資源共享服務(wù)平臺(甘肅)、“西北干旱區(qū)抗逆農(nóng)作物種質(zhì)資源調(diào)查”項(xiàng)目組提供。

    1.2 表型鑒定

    2016-2017年在甘肅蘭州和張掖進(jìn)行和尚頭小麥的田間鑒定,試驗(yàn)設(shè)3 m行長、0.2 m行距、2行區(qū),3次重復(fù)。按照《小麥種質(zhì)資源描述規(guī)范和數(shù)據(jù)標(biāo)準(zhǔn)》[26]記載芽鞘色、幼苗色、幼苗習(xí)性、株形和葉姿;成熟后每小區(qū)取樣20株進(jìn)行室內(nèi)考種,考取株高、穗長、小穗數(shù)、小穗粒數(shù)、有效分蘗數(shù)、植株整齊度、穗形、殼色、芒形、芒色、粒色、粒質(zhì)、飽滿度和千粒重。

    表1 43份和尚頭種質(zhì)名稱及來源Table 1 The name and origin of 43 Heshangtou germplasms

    1.3 SSR標(biāo)記分析

    1.3.1DNA提取與檢測 將材料放置在培養(yǎng)箱中暗培養(yǎng),在幼苗期取0.2 g新鮮葉片,參照改良CTAB法[27]提取和尚頭小麥DNA,用1%瓊脂糖凝膠電泳檢測質(zhì)量,用NanoDrop 2000c測定其濃度,最后加入100 μL TE緩沖液置于-20 ℃保存?zhèn)溆谩?/p>

    1.3.2SSR擴(kuò)增及電泳檢測 選擇均勻分布于小麥21條染色體上150對已公布的SSR引物,同時(shí)選取甘肅、陜西、寧夏的5份供試材料進(jìn)行PCR擴(kuò)增。根據(jù)擴(kuò)增產(chǎn)物的電泳結(jié)果,選擇擴(kuò)增條帶清晰、多態(tài)性好的SSR引物對43份和尚頭小麥基因組DNA進(jìn)行遺傳多樣性檢測。公共引物序列均來自Graingenes(http://wheat.pw.usda.gov),由上海生工生物工程公司合成。PCR反應(yīng)體系為15 μL,包含:10×buffer (+MgCl2) 1.5 μL;dNTP (2.5 mmol·L-1) 1.2 μL;Taq酶(5 U·μL-1) 0.3 μL;Forward primer (10 μmol·L-1) 0.4 μL;Reverse primer(10 μmol·L-1) 0.4 μL;DNA(20~50 ng·μL-1)2.0 μL;ddH2O 9.2 μL。熱循環(huán)程序?yàn)椋?4 ℃預(yù)變性5 min,94 ℃變性50 s,55 ℃復(fù)性40 s,72 ℃延伸50 s,35次循環(huán),72 ℃延伸10 min。擴(kuò)增產(chǎn)物經(jīng)6%非變性聚丙烯酰胺凝膠120 V穩(wěn)壓電泳,利用銀染法顯色,拍照保存后統(tǒng)計(jì)條帶。

    1.4 數(shù)據(jù)統(tǒng)計(jì)分析

    表2 和尚頭種質(zhì)資源表型性狀賦值和分級標(biāo)準(zhǔn)Table 2 Coden designed for phenotypic traits and grading standard in Heshangtou germplasms

    BSC:芽鞘色;SC:幼苗色;ST:幼苗習(xí)性;PT:株形;LP:葉姿;PU:植株整齊度;SS:穗形;GC:殼色;AS:芒形;AC:芒色;KC:粒色;KT:粒質(zhì);PP:飽滿度;PH:株高;SL:穗長;SN:小穗數(shù);ETN:有效分蘗數(shù);KNS:小穗粒數(shù);KW:千粒重。下同。

    BSC: Bud-sheath color; SC: Seedling color; ST: Seedling type; PT: Plant type; LP: Leaf posture; PU: Plant uniformity; SS: Spike shape; GC: Glume color; AS: Awn shape; AC: Awn color; KC: Kernel color; KT: Kernel texture; PP: Plumpness; PH: Plant height; SL: Spike length; SN: Spikelet number; ETN: Effective tiller number; KNS: Kernel number per spikelet; KW: 1000-kernel weight. The same below.

    1.4.2基因型數(shù)據(jù)統(tǒng)計(jì)分析 在Microsoft Excel 上建立數(shù)據(jù)庫,將SSR擴(kuò)增帶型在相同遷移率位置上有帶記為“1”,無帶記為“0”,缺失記為“9”,利用PIC_Calc 0.6軟件對基因型數(shù)據(jù)進(jìn)行統(tǒng)計(jì)分析,計(jì)算多態(tài)性信息含量(polymorphism information content,PIC值)。利用NTSYSpc(2.10e版)軟件進(jìn)行數(shù)據(jù)處理,應(yīng)用SIMQUAL法計(jì)算各種質(zhì)間遺傳相似系數(shù)(F)和遺傳距離(D)。用SAHN功能,按照非加權(quán)成組配對法(UPGMA)進(jìn)行聚類分析。

    用Structure 2.3.4軟件分析群體的遺傳結(jié)構(gòu),估計(jì)最佳群體組群數(shù)K,假定材料的群體數(shù)為2~10,且假定各位點(diǎn)相互獨(dú)立進(jìn)行分析,將MCMC(markov chain monte carlo)開始時(shí)的不作數(shù)迭代(length of burnin period)設(shè)為10000次,再將不作數(shù)迭代后的MCMC設(shè)為100000,運(yùn)行20次,根據(jù)似然值最大原則選擇合適群體K值,或利用ΔK的變化規(guī)律來確定最適的群體數(shù)目,構(gòu)建遺傳結(jié)構(gòu)圖。ΔK的計(jì)算公式為:ΔK=m(|L(K+1)-2L(K)+L(K-1)|)/s[L(K)]。L(K)為每個(gè)K對應(yīng)的對數(shù)值,s為標(biāo)準(zhǔn)差,m為平均值。

    1.5 關(guān)聯(lián)分析

    利用Tassel 2.1軟件對和尚頭小麥主要農(nóng)藝性狀與SSR標(biāo)記進(jìn)行相關(guān)性分析,采用一般線性模型(general linear model,GLM),以Q作為協(xié)變量進(jìn)行回歸分析,并計(jì)算標(biāo)記對表型變異的解釋率。

    2 結(jié)果與分析

    2.1 和尚頭小麥種質(zhì)表型性狀的遺傳多樣性分析

    2.1.1質(zhì)量性狀多樣性分析 通過對和尚頭小麥種質(zhì)的13個(gè)質(zhì)量性狀統(tǒng)計(jì)分析,43份種質(zhì)芽鞘色全部為綠色,幼苗色中綠色占95.35%,深綠色占4.65%,幼苗習(xí)性以直立為主,半匍匐和匍匐各占4.65%,97.67%的種質(zhì)粒色為紅色,僅有1份材料為白色,上述4個(gè)性狀的遺傳多樣性指數(shù)較小,說明43份和尚頭種質(zhì)的芽鞘色、幼苗色、幼苗習(xí)性和粒色變異度??;在質(zhì)量性狀中,葉姿、植株整齊度、芒形和芒色4個(gè)性狀的遺傳多樣性指數(shù)在1左右,豐富度高,葉姿中挺直、平展和下披分別占20.93%、27.91%、51.16%,58.14%的種質(zhì)無芒,23.26%為短芒材料,甘肅的1份種質(zhì)芒色為黑色,表現(xiàn)較特殊;株形、穗形、殼色、粒質(zhì)和飽滿度5個(gè)性狀的遺傳多樣性指數(shù)介于上述8個(gè)性狀之間(表3),變異度較高。總體來看,43份和尚頭種質(zhì)在質(zhì)量性狀中差異較明顯,表現(xiàn)出了一定的變異度。

    表3 和尚頭種質(zhì)資源13個(gè)質(zhì)量性狀的多樣性分析Table 3 Diversity analysis for 13 qualitative traits of Heshangtou germplasms

    2.1.2數(shù)量性狀多樣性分析 對供試材料的6個(gè)數(shù)量性狀進(jìn)行遺傳多樣性分析,從表4中可以看出,6個(gè)性狀的變異系數(shù)介于11.02%~49.03%,由大到小依次為有效分蘗數(shù)、穗長、小穗粒數(shù)、千粒重、株高和小穗數(shù)。有效分蘗數(shù)的變異系數(shù)最大,且極差大于均值,陜西省編號為W17和W21的和尚頭小麥有效分蘗數(shù)最少,為3個(gè),甘肅省編號為W37的為21.5個(gè),說明43份和尚頭種質(zhì)的有效分蘗數(shù)更為分散,遺傳改良潛力較大。小穗數(shù)的變異范圍最小,為16.00~25.25個(gè)。

    與質(zhì)量性狀相比,6個(gè)數(shù)量性狀的遺傳多樣性指數(shù)普遍較高,均在1.5632以上,其中千粒重的遺傳多樣性指數(shù)最高,為1.9065,小穗粒數(shù)的最低,為1.5632,株高、穗長、小穗數(shù)和有效分蘗數(shù)的遺傳多樣性指數(shù)比較接近,在1.8500左右,說明各性狀比較均勻,豐富度高。對比質(zhì)量性狀和數(shù)量性狀,43份和尚頭小麥種質(zhì)在數(shù)量性狀方面改良潛力大,可為小麥品質(zhì)育種提供材料基礎(chǔ)。

    表4 和尚頭種質(zhì)資源6個(gè)數(shù)量性狀統(tǒng)計(jì)分析Table 4 Statistical analysis of 6 numerical traits of Heshangtou germplasms

    2.2 和尚頭小麥種質(zhì)表型性狀的相關(guān)性分析

    圖1 和尚頭種質(zhì)資源基于18個(gè)形態(tài)性狀的聚類圖Fig.1 Cluster tree based on 18 morphological traits of Heshangtou germplasms

    由于43份材料芽鞘色全部為綠色,因此對表型性狀進(jìn)行相關(guān)性和主成分分析時(shí)剔除了芽鞘色。從表5可以看出,殼色與所有性狀間均不顯著相關(guān),其他不同表型性狀間呈不同程度的相關(guān)性。在質(zhì)量性狀中,幼苗色與芒色和粒色呈極顯著正相關(guān),株形與植株整齊度呈顯著正相關(guān),幼苗習(xí)性與株形、芒形與穗形呈顯著負(fù)相關(guān);數(shù)量性狀中,穗長與株高呈極顯著負(fù)相關(guān),穗長與小穗數(shù)和千粒重呈極顯著正相關(guān),小穗數(shù)與千粒重呈極顯著正相關(guān)。

    2.3 基于表型特征的聚類分析

    利用DPS軟件,對表型數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化轉(zhuǎn)換,采用歐氏距離為遺傳距離,UPGMA方法對43份和尚頭小麥種質(zhì)進(jìn)行表型聚類分析(圖1),在遺傳相似系數(shù)為6.0處將43份和尚頭小麥種質(zhì)分為6類,第Ⅰ類只有陜西的1份材料,其粒色為白色與其他42份材料粒色均不同;第Ⅱ類同樣為1份材料,來自甘肅,芒色為黑色;第Ⅲ類包含2份材料,1份來自甘肅,1份為陜西材料,幼苗色均為深綠色;第Ⅳ類僅有1份材料,來自陜西,為無芒,穗形為棍棒型;第Ⅴ類包括4份材料,1份來自甘肅,3份來自陜西,4份材料的芒色和殼色均為白色;第Ⅵ類包含34份材料,大部分材料被劃分在此類,其中甘肅材料22份,陜西材料11份,寧夏1份。

    2.4 SSR多態(tài)性分析

    從均勻分布于小麥21條染色體上的150對SSR引物中篩選出45對條帶清晰、多態(tài)性好的引物,多態(tài)性比率為30%。45對引物共檢測出151個(gè)等位位點(diǎn),引物等位點(diǎn)數(shù)在2~8個(gè),平均每個(gè)引物等位數(shù)為3.36個(gè)。引物Barc78和gwm174的等位位點(diǎn)數(shù)最多(圖2)。引物等位位點(diǎn)多態(tài)性信息含量PIC變幅為0.044~0.771,平均0.345。以引物gwm174值最高(0.771),引物gwm135值最低(0.044)。

    圖2 SSR引物Barc78對43份和尚頭小麥種質(zhì)資源的PCR擴(kuò)增結(jié)果Fig.2 Amplified result of SSR markers Barc78 of 43 Heshangtou germplasms

    編號Number引物Primer連鎖群Linkage group等位變異數(shù)Alleles number多態(tài)性信息含量PIC編號Number引物Primer連鎖群Linkage group等位變異數(shù)Alleles number多態(tài)性信息含量PICM2gwm1351A20.0444 M47wmc7275A20.2113 M3barc1581A20.3511 M48barc1415A20.1214 M4wmc3671B60.3516 M51gwm1545A20.0444 M7wmc4291D40.2743 M52wmc5375B40.4809 M8wmc1471D20.0444 M53gwm1745D80.7713 M9cfd2821D30.3936 M54gwm4596A20.2932 M11gwm3722A40.6865 M56gwm6176A40.5785 M12gwm3112A50.5481 M58wmc2016A20.1214 M13wmc1542B20.3083 M60gwm2196B20.3716 M15gwm3742B30.1631 M61gwm4696D30.2010 M16wmc5032D50.6577 M62barc1756D40.1688 M17gwm3492D40.2354 M64barc1966D20.1545 M20gwm4803A30.4027 M37gwm6357A30.4980 M22gwm3893B20.3214 M38wmc3357B20.3083 M23gwm2993B40.2897 M39gwm4007B30.2665 M24Barc1643B30.1273 M40wmc1217D40.3942 M25wmc5293D20.0444 M65gwm6357D40.6497 M28wmc4204A20.1844 M66Xwmc7907A30.4496 M29barc1704A60.7407 M33Xwmc6337A50.6093 M30barc784A80.5691 M34BF1459356A30.3989 M32gwm1134B20.0444 M69Xbarc1521D40.5221 M45gwm1494B30.4244 M70Xcfd151A40.4216 M46wmc2854D20.2932

    2.5 基于SSR分子標(biāo)記的聚類分析

    通過NTSYS 2.10e分析軟件,應(yīng)用SIMQUAL法計(jì)算各種質(zhì)間遺傳相似系數(shù)(F)和遺傳距離(D),用SAHN功能,按照非加權(quán)成組配對法(UPGMA)進(jìn)行聚類分析,在距離0.705處可將43份種質(zhì)分為五類(圖3)。第Ⅰ類只包含甘肅的黑芒和尚頭1份材料;第Ⅱ類包括9份材料,全部為陜西材料;第Ⅲ類包含陜西的2份材料;第Ⅳ類只有編號為W7的1份甘肅材料;第Ⅴ類包括30份材料,23份來自甘肅,6份來自陜西,1份為寧夏材料。從分類情況看,整體分布與地理來源存在一定的關(guān)系,大多數(shù)不同來源的材料被劃分在了不同的群組,少數(shù)材料存在相互間的滲透。

    圖3 和尚頭種質(zhì)資源基于SSR分子標(biāo)記的聚類圖Fig.3 Cluster tree based on SSR molecular markers of Heshangtou germplasms

    2.6 基于數(shù)學(xué)模型的遺傳結(jié)構(gòu)分析

    利用篩選出的45對條帶清晰、多態(tài)性好的引物用于群體遺傳結(jié)構(gòu)分析,從圖4A中可以看出,隨著K值的增大,lnP(D)值呈現(xiàn)上升的趨勢,并沒有出現(xiàn)明顯的拐點(diǎn),無法劃分群體的亞群數(shù),因此,參照Evanno等[30]的方法通過ΔK來確定K值。在K=8時(shí),ΔK達(dá)到最大值并出現(xiàn)明顯的峰值(圖4B),因此,可將43份和尚頭小麥種質(zhì)劃分為8個(gè)群組,8個(gè)亞群分別包含了2、2、7、1、17、4、1和2份材料,并繪制了群體遺傳結(jié)構(gòu)圖(圖5)。通過對Q值分析,有36份材料被分在某一群組中的Q值大于0.6,說明其遺傳組分相對單一,被劃分在8個(gè)群組中的某一組,7份材料在8個(gè)群組中Q值均小于0.6,形成混合群組。第Ⅱ、Ⅳ、Ⅴ、Ⅵ、Ⅶ組中材料全部來源于甘肅?。坏冖?、Ⅲ、Ⅷ組中材料全部來自陜西?。换旌辖M群中包含陜西省6份材料和寧夏的1份材料。通過遺傳結(jié)構(gòu)分析表明,不同來源的和尚頭小麥種質(zhì)遺傳背景差異較大,同一省份的材料間遺傳多樣性比較豐富,可為小麥品種的遺傳改良提供優(yōu)異的基礎(chǔ)材料。

    2.7 SSR分子標(biāo)記與主要農(nóng)藝性狀關(guān)聯(lián)分析

    2.7.1和尚頭小麥種質(zhì)資源SSR位點(diǎn)間連鎖不平衡分析 基因間的連鎖不平衡(linkage disequilibrium,LD)是關(guān)聯(lián)分析的前提和基礎(chǔ),分析45對SSR引物位點(diǎn)990種組合中連鎖不平衡有利于了解小麥基因組的連鎖不平衡,為和尚頭小麥種質(zhì)資源的關(guān)聯(lián)分析提供基礎(chǔ)。從圖6可以看出,在和尚頭小麥基因組中存在著不同程度的連鎖不平衡,包括共線性的(同一染色體)和非共線性的(不同染色體),說明不同染色體的不同區(qū)段發(fā)生過重組或者突變,其中R2>0.5的LD有10個(gè)。在P<0.01的概率支持下成對存在的不平衡位點(diǎn)有60個(gè)(圖6中為黑色斜線下方紅、綠、藍(lán)色小格),占整個(gè)組合的6.06%。從分析結(jié)果來看,本研究所選用的45對SSR引物在43份和尚頭小麥種質(zhì)資源中存在連鎖不平衡,可與主要農(nóng)藝性狀進(jìn)行關(guān)聯(lián)分析。

    圖4 K值與LnP(D)、ΔK值折線圖Fig.4 Lines chart of K with lnP(D) and ΔK A:K值與lnP(D)值的折線圖;B:K值與ΔK的折線圖。A: Lines chart of K with lnP(D); B: Lines chart of K with ΔK.

    圖5 43份和尚頭小麥種質(zhì)群體遺傳結(jié)構(gòu)Fig.5 Population genetic structure diagram of 43 Heshangtou germplasms

    圖6 45對SSR多態(tài)性位點(diǎn)間的連鎖不平衡Fig.6 Linkage disequilibrium among 45 SSR polymorphic sites

    2.7.2和尚頭小麥種質(zhì)資源主要農(nóng)藝性狀相關(guān)聯(lián)的SSR位點(diǎn) 將45對有多態(tài)性的SSR標(biāo)記與株高、穗長、小穗數(shù)、有效分蘗數(shù)、小穗粒數(shù)和千粒重等6個(gè)小麥主要農(nóng)藝性狀利用Tassel 2.1軟件,采用GLM模型進(jìn)行關(guān)聯(lián)分析,在P<0.01的水平上, 張掖和蘭州兩個(gè)試驗(yàn)點(diǎn)共獲得與株高、穗長、小穗數(shù)和小穗粒數(shù)相關(guān)聯(lián)的SSR標(biāo)記11個(gè)(表7),而有效分蘗數(shù)和千粒重間沒有發(fā)現(xiàn)與之顯著關(guān)聯(lián)的標(biāo)記。在11個(gè)相關(guān)聯(lián)的標(biāo)記中,與株高相關(guān)聯(lián)的標(biāo)記最多,為6個(gè),穗長的最少為1個(gè),標(biāo)記wmc201與穗長和小穗粒數(shù)均顯著相關(guān)。11個(gè)標(biāo)記表型變異解釋率為8.89%~24.74%,其中與株高相關(guān)聯(lián)的標(biāo)記表型變異解釋率最高,小穗粒數(shù)次之,小穗數(shù)最小。

    3 討論

    3.1 和尚頭小麥種質(zhì)資源遺傳多樣性分析

    遺傳多樣性是指地球上所有生物所攜帶的遺傳信息的總和。一般而言,遺傳多樣性是物種內(nèi)的遺傳多樣性,即種內(nèi)不同群體或一個(gè)群體內(nèi)不同個(gè)體的遺傳變異總和。分析小麥遺傳多樣性,有利于小麥種質(zhì)資源的收集、鑒定、保存和利用,對小麥品種的遺傳改良具有重要意義。小麥遺傳多樣性研究經(jīng)歷了從最早的形態(tài)學(xué)水平[31],到細(xì)胞水平[32]、生理生化水平[33],再到目前主要采用的DNA分子水平[34]4個(gè)階段。根據(jù)本課題組的實(shí)地調(diào)查和前人的研究,和尚頭小麥具有耐瘠薄、抗逆性強(qiáng)等優(yōu)良性狀,其面粉質(zhì)量好,尤其是蛋白質(zhì)含量高,具有滑潤爽口、味感純正、面筋強(qiáng)、食用方便等特點(diǎn),在西北地區(qū)享有很高的名譽(yù),但產(chǎn)量相對較低,而且前人并沒有對不同省份和地區(qū)的和尚頭小麥種質(zhì)資源進(jìn)行遺傳多樣性分析。因此,本研究利用形態(tài)學(xué)標(biāo)記和SSR分子標(biāo)記分析和尚頭小麥種質(zhì)資源的遺傳多樣性,以期為特色小麥品種的遺傳改良提供參考依據(jù)。本研究利用19個(gè)表型性狀和45對多態(tài)性SSR標(biāo)記對來自西北地區(qū)的43份和尚頭小麥種質(zhì)進(jìn)行遺傳多樣性研究,對比基于形態(tài)性狀、SSR標(biāo)記的聚類和群體遺傳結(jié)構(gòu)分析,基于分子標(biāo)記的兩種群體劃分結(jié)果一致性比較高,基于形態(tài)性狀的群體劃分結(jié)果與另外兩種方法稍有差別,但大部分和尚頭小麥種質(zhì)資源可歸于一類。通過綜合分析,可將43份和尚頭小麥種質(zhì)劃分為8類,發(fā)現(xiàn)和尚頭小麥種質(zhì)資源與地理來源存在一定的關(guān)系,大多數(shù)不同來源的材料被劃分在了不同的群組,少數(shù)材料存在相互間的滲透,其中編號為W7,W8的甘肅和尚頭小麥和W16,W17,W20,W29,W31,W33和W35的陜西和尚頭小麥種質(zhì)較其他種質(zhì)遺傳背景差異大;編號為W36的寧夏和尚頭小麥不能夠單獨(dú)被劃分在一個(gè)組群,說明與甘肅和陜西的遺傳背景較近。通過對表型特征進(jìn)行分析,W8為黑芒,千粒重達(dá)到了53.78 g;W16和W17小穗數(shù)較多,為22和23個(gè);W20是從長芒W18中分離出來的一個(gè)頂芒種質(zhì);W29的生長類型為匍匐型;W31的株高較低,為62.8 cm;W33和W35千粒重均較大,分別為41.41和46.36 g,且W35的株高為86.9 cm,相對較低,這些遺傳背景差異較大的和尚頭小麥種質(zhì)資源可供育種家選擇利用。

    表7 與農(nóng)藝性狀顯著相關(guān)的位點(diǎn)及其對表型變異的解釋率Table 7 Loci associated with agronomic traits and percentage of phenotypic variation explained

    本研究發(fā)現(xiàn),不同來源的和尚頭小麥種質(zhì)遺傳多樣性比較豐富,且同一省份的和尚頭小麥種質(zhì)間也存在較大差異,這種同名異質(zhì)現(xiàn)象是由多方面原因引起的,第一種可能是由于老百姓根據(jù)和尚頭小麥種質(zhì)的形態(tài)特征賦予其名;第二種可能原因是在和尚頭小麥長期種植的過程中發(fā)生了基因突變,也有可能是育種家對和尚頭小麥進(jìn)行了選擇或雜交改良,導(dǎo)致這些同名和尚頭小麥種質(zhì)的遺傳背景有所差異,但這需要進(jìn)一步考證。這些遺傳多樣性豐富的和尚頭小麥種質(zhì)資源可為小麥品種的遺傳改良提供優(yōu)異的基礎(chǔ)材料。

    3.2 和尚頭小麥種質(zhì)資源產(chǎn)量相關(guān)性狀關(guān)聯(lián)分析

    發(fā)掘優(yōu)異基因資源是作物種質(zhì)資源分子評價(jià)的重要內(nèi)容,對作物分子標(biāo)記輔助選擇育種具有非常重要的實(shí)踐意義,基于連鎖不平衡(LD)的關(guān)聯(lián)分析是基因發(fā)掘也是等位基因發(fā)掘的有效方法[35]。前人基于全基因組掃描和候選基因進(jìn)行關(guān)聯(lián)分析研究,找到了許多目標(biāo)性狀相關(guān)聯(lián)的標(biāo)記或者數(shù)量性狀遺傳位點(diǎn)(quantitative trait locus,QTL)。2006年Breseghello等[36]首次利用基于全基因組掃描的關(guān)聯(lián)分析發(fā)現(xiàn)了14個(gè)與粒重和粒長相關(guān)的標(biāo)記。Yao等[37]對108份小麥品種的6個(gè)農(nóng)藝性狀進(jìn)行關(guān)聯(lián)分析,發(fā)現(xiàn)14個(gè)SSR標(biāo)記與農(nóng)藝性狀顯著相關(guān);Dodig等[38]利用SSR標(biāo)記對小麥品種農(nóng)藝性狀和產(chǎn)量性狀進(jìn)行關(guān)聯(lián)分析,在不同環(huán)境條件下共檢測到76個(gè)與性狀相關(guān)聯(lián)的標(biāo)記位點(diǎn);Neumann等[39]利用兩種不同的模型方法對96份小麥種質(zhì)的20個(gè)農(nóng)藝性狀進(jìn)行關(guān)聯(lián)分析,共發(fā)現(xiàn)了385個(gè)與目標(biāo)性狀顯著相關(guān)的標(biāo)記位點(diǎn);Zhang等[40]對小麥籽粒性狀進(jìn)行關(guān)聯(lián)分析,發(fā)現(xiàn)27個(gè)SSR標(biāo)記位點(diǎn)與之顯著關(guān)聯(lián),其中23個(gè)位點(diǎn)極顯著相關(guān)。

    本研究采用GLM模型對43份和尚頭小麥種質(zhì)資源的株高、穗長、小穗數(shù)、有效分蘗數(shù)、小穗粒數(shù)和千粒重等6個(gè)主要農(nóng)藝性狀進(jìn)行關(guān)聯(lián)分析,獲得了與株高、穗長、小穗數(shù)和小穗粒數(shù)相關(guān)聯(lián)的SSR標(biāo)記11個(gè),與株高相關(guān)聯(lián)的標(biāo)記最多為6個(gè)。通過分析前人的研究結(jié)果,發(fā)現(xiàn)標(biāo)記wmc537與余侃[41]的研究結(jié)果相一致,標(biāo)記cfd282與趙京嵐[42]的研究結(jié)果相一致,兩個(gè)標(biāo)記均與和尚頭小麥種質(zhì)的株高相關(guān)聯(lián),分別解釋了22.19%和24.74%的表型變異;標(biāo)記gwm372與李小軍等[43]的研究結(jié)果一致,與小穗數(shù)相關(guān)聯(lián),解釋了9.6%的表型變異。同時(shí)發(fā)現(xiàn)標(biāo)記wmc201與穗長和小穗粒數(shù)均顯著相關(guān)聯(lián),而有效分蘗數(shù)和千粒重間沒有發(fā)現(xiàn)與之顯著關(guān)聯(lián)的標(biāo)記。因此,本研究結(jié)果與前人的研究具有相似之處,說明本研究結(jié)果具有較高的參考價(jià)值。

    4 結(jié)論

    利用表型特征和分子標(biāo)記分析了43份和尚頭小麥種質(zhì)資源的遺傳多樣性,通過群體結(jié)構(gòu)分析將其劃分為8個(gè)組群,發(fā)現(xiàn)西北地區(qū)和尚頭小麥種質(zhì)資源具有豐富的遺傳多樣性,其中編號為W7,W8的甘肅和尚頭小麥和W16,W17,W20,W29,W31,W33和W35的陜西和尚頭小麥種質(zhì)較其他種質(zhì)遺傳背景差異大,可為特色小麥品種遺傳改良提供材料基礎(chǔ);通過關(guān)聯(lián)分析找到了11個(gè)與株高、穗長、小穗數(shù)和小穗粒數(shù)相關(guān)聯(lián)的SSR標(biāo)記,為特色小麥分子標(biāo)記輔助選擇提供參考依據(jù)。

    猜你喜歡
    表型種質(zhì)遺傳
    非遺傳承
    華南地區(qū)最大農(nóng)作物種質(zhì)資源保護(hù)庫建成
    還有什么會(huì)遺傳?
    還有什么會(huì)遺傳
    還有什么會(huì)遺傳?
    亞麻抗白粉病種質(zhì)資源的鑒定與篩選
    建蘭、寒蘭花表型分析
    貴州玉米種質(zhì)資源遺傳多樣性及核心種質(zhì)庫構(gòu)建
    紅錐種質(zhì)早期生長表現(xiàn)
    GABABR2基因遺傳變異與肥胖及代謝相關(guān)表型的關(guān)系
    久久人妻福利社区极品人妻图片| 亚洲熟妇熟女久久| 一区二区三区国产精品乱码| 久久精品成人免费网站| 欧美av亚洲av综合av国产av| 国产1区2区3区精品| 欧美亚洲 丝袜 人妻 在线| videos熟女内射| xxxhd国产人妻xxx| av有码第一页| 成年版毛片免费区| 久久久精品国产亚洲av高清涩受| 日韩三级视频一区二区三区| 日韩免费高清中文字幕av| 久久亚洲真实| www.自偷自拍.com| 法律面前人人平等表现在哪些方面| 久久精品aⅴ一区二区三区四区| 久久国产精品人妻蜜桃| 1024视频免费在线观看| 国产一卡二卡三卡精品| 热99国产精品久久久久久7| 最新在线观看一区二区三区| 亚洲久久久国产精品| 亚洲国产欧美网| 免费观看人在逋| 天天躁夜夜躁狠狠躁躁| 黑人猛操日本美女一级片| videosex国产| 女人爽到高潮嗷嗷叫在线视频| 天堂中文最新版在线下载| 国产一区二区三区视频了| 中文字幕色久视频| 欧美激情极品国产一区二区三区| 嫩草影视91久久| 精品高清国产在线一区| 精品一区二区三区四区五区乱码| 成年人免费黄色播放视频| 国产一区二区三区视频了| 国产亚洲精品一区二区www | 一级毛片女人18水好多| 久久久久久亚洲精品国产蜜桃av| 亚洲av欧美aⅴ国产| 美女高潮到喷水免费观看| 久久久精品免费免费高清| 女人被躁到高潮嗷嗷叫费观| 天天躁日日躁夜夜躁夜夜| 精品国产乱子伦一区二区三区| 中文字幕人妻熟女乱码| 亚洲av欧美aⅴ国产| 精品国产一区二区三区久久久樱花| 亚洲精品一二三| 国产不卡一卡二| 精品一区二区三卡| 国产无遮挡羞羞视频在线观看| 真人做人爱边吃奶动态| 欧美日韩亚洲国产一区二区在线观看 | 亚洲五月色婷婷综合| 久久青草综合色| 国产熟女午夜一区二区三区| 亚洲国产欧美日韩在线播放| 韩国av一区二区三区四区| 欧美精品av麻豆av| 精品久久久久久电影网| 18在线观看网站| 大陆偷拍与自拍| 国产欧美日韩一区二区三| 在线视频色国产色| 亚洲熟女毛片儿| 99国产精品99久久久久| aaaaa片日本免费| 亚洲av欧美aⅴ国产| 亚洲精品国产精品久久久不卡| 亚洲国产精品sss在线观看 | 午夜亚洲福利在线播放| 国产又色又爽无遮挡免费看| 好看av亚洲va欧美ⅴa在| 国产成人av教育| 女人久久www免费人成看片| 成人av一区二区三区在线看| 搡老乐熟女国产| 亚洲精品自拍成人| 老司机靠b影院| 建设人人有责人人尽责人人享有的| 国产精品一区二区在线不卡| 99riav亚洲国产免费| 亚洲性夜色夜夜综合| 亚洲成人国产一区在线观看| 国产又爽黄色视频| 久久婷婷成人综合色麻豆| 手机成人av网站| 亚洲精品久久午夜乱码| 国产精品久久电影中文字幕 | 久久久国产成人精品二区 | 一个人免费在线观看的高清视频| videos熟女内射| 飞空精品影院首页| www日本在线高清视频| 精品久久久久久久久久免费视频 | 亚洲少妇的诱惑av| 国产精品亚洲一级av第二区| 欧美日韩精品网址| 又黄又粗又硬又大视频| 丝袜美足系列| 人妻丰满熟妇av一区二区三区 | 亚洲综合色网址| 女警被强在线播放| 高清毛片免费观看视频网站 | 天堂俺去俺来也www色官网| 中文字幕av电影在线播放| 高清av免费在线| 精品一区二区三区四区五区乱码| 久久性视频一级片| 日韩欧美免费精品| 成人亚洲精品一区在线观看| 一二三四社区在线视频社区8| 在线永久观看黄色视频| 亚洲全国av大片| 亚洲 欧美一区二区三区| 国产欧美日韩综合在线一区二区| 亚洲五月天丁香| 电影成人av| 18禁美女被吸乳视频| 国产成人精品久久二区二区91| 国内毛片毛片毛片毛片毛片| 宅男免费午夜| 国产不卡一卡二| 天堂俺去俺来也www色官网| 国产精品乱码一区二三区的特点 | 91字幕亚洲| 欧美午夜高清在线| 精品人妻1区二区| 国产视频一区二区在线看| a级毛片黄视频| 国产精品乱码一区二三区的特点 | 国产精品久久久av美女十八| av天堂在线播放| 国产高清视频在线播放一区| 在线观看www视频免费| 建设人人有责人人尽责人人享有的| 国产在视频线精品| 我的亚洲天堂| 欧美另类亚洲清纯唯美| 欧美乱妇无乱码| 视频在线观看一区二区三区| 亚洲七黄色美女视频| 天堂√8在线中文| 一级黄色大片毛片| videosex国产| 日本一区二区免费在线视频| 在线十欧美十亚洲十日本专区| 欧美日韩视频精品一区| 俄罗斯特黄特色一大片| 亚洲精品粉嫩美女一区| 欧美黑人欧美精品刺激| 久热这里只有精品99| 国产精品久久久av美女十八| 国产精品久久久av美女十八| 村上凉子中文字幕在线| 国产成人系列免费观看| videos熟女内射| 国产亚洲精品第一综合不卡| ponron亚洲| 亚洲精品中文字幕在线视频| 久久久精品区二区三区| 建设人人有责人人尽责人人享有的| 午夜两性在线视频| 麻豆成人av在线观看| 亚洲熟女毛片儿| 国产欧美日韩一区二区三区在线| 日韩免费av在线播放| 亚洲精品一卡2卡三卡4卡5卡| 变态另类成人亚洲欧美熟女 | 一进一出好大好爽视频| 亚洲国产毛片av蜜桃av| videos熟女内射| 国产成人系列免费观看| www.自偷自拍.com| 亚洲成人免费av在线播放| 亚洲全国av大片| 亚洲国产欧美日韩在线播放| 日本黄色视频三级网站网址 | 精品国产美女av久久久久小说| 久久国产精品男人的天堂亚洲| 国产真人三级小视频在线观看| 亚洲第一av免费看| 午夜福利欧美成人| 国产蜜桃级精品一区二区三区 | cao死你这个sao货| 亚洲精品粉嫩美女一区| 国产av精品麻豆| 国产av精品麻豆| 国产av精品麻豆| 亚洲一码二码三码区别大吗| 午夜亚洲福利在线播放| 国产伦人伦偷精品视频| 满18在线观看网站| 叶爱在线成人免费视频播放| 国产一区二区三区在线臀色熟女 | 国产精品国产av在线观看| 久久精品熟女亚洲av麻豆精品| 1024香蕉在线观看| 欧美黑人欧美精品刺激| 大香蕉久久网| 大陆偷拍与自拍| 女人久久www免费人成看片| 国产亚洲精品一区二区www | 99热网站在线观看| 不卡av一区二区三区| 99re在线观看精品视频| 母亲3免费完整高清在线观看| 国产人伦9x9x在线观看| 日本黄色视频三级网站网址 | 国产成人免费无遮挡视频| 男女午夜视频在线观看| 日本a在线网址| 国产成人啪精品午夜网站| 少妇 在线观看| 宅男免费午夜| 免费高清在线观看日韩| 丝袜人妻中文字幕| 操出白浆在线播放| 欧美丝袜亚洲另类 | 久久久精品区二区三区| av免费在线观看网站| 美女国产高潮福利片在线看| 超碰成人久久| 777久久人妻少妇嫩草av网站| 久热爱精品视频在线9| 老熟妇仑乱视频hdxx| 黑人巨大精品欧美一区二区蜜桃| ponron亚洲| 亚洲av成人不卡在线观看播放网| 国产在线一区二区三区精| tube8黄色片| 国产xxxxx性猛交| av在线播放免费不卡| 人人妻人人爽人人添夜夜欢视频| 亚洲精品粉嫩美女一区| 一夜夜www| 成年人黄色毛片网站| 久久中文字幕一级| 性色av乱码一区二区三区2| 亚洲aⅴ乱码一区二区在线播放 | 精品熟女少妇八av免费久了| 一边摸一边抽搐一进一小说 | 久久影院123| 精品第一国产精品| 黄片播放在线免费| 两个人看的免费小视频| 91精品国产国语对白视频| 中文字幕高清在线视频| 不卡av一区二区三区| 午夜激情av网站| 天堂俺去俺来也www色官网| 女人高潮潮喷娇喘18禁视频| 国产精品久久久av美女十八| 两个人免费观看高清视频| 欧美午夜高清在线| 国产三级黄色录像| 51午夜福利影视在线观看| 亚洲精品自拍成人| 亚洲欧美激情综合另类| av天堂久久9| 一边摸一边抽搐一进一出视频| 国产高清视频在线播放一区| 欧美在线黄色| 在线免费观看的www视频| 黄频高清免费视频| 国产熟女午夜一区二区三区| 女人被狂操c到高潮| 观看免费一级毛片| 日韩精品中文字幕看吧| 99久久久亚洲精品蜜臀av| 国产野战对白在线观看| 美女黄网站色视频| 男女之事视频高清在线观看| 久久午夜亚洲精品久久| 欧美成人免费av一区二区三区| 一级毛片女人18水好多| 欧美激情久久久久久爽电影| 嫩草影视91久久| 美女高潮喷水抽搐中文字幕| 日韩欧美免费精品| 精品欧美国产一区二区三| 亚洲男人的天堂狠狠| 亚洲一区二区三区色噜噜| 国产一区二区亚洲精品在线观看| 国产69精品久久久久777片| 欧美乱色亚洲激情| 国产精品久久久久久久电影 | 母亲3免费完整高清在线观看| 亚洲最大成人手机在线| 狂野欧美激情性xxxx| 在线国产一区二区在线| 免费在线观看影片大全网站| 黄片小视频在线播放| 亚洲精华国产精华精| 久久草成人影院| netflix在线观看网站| 麻豆久久精品国产亚洲av| www.色视频.com| 天天躁日日操中文字幕| 桃色一区二区三区在线观看| 美女被艹到高潮喷水动态| 精品福利观看| 又粗又爽又猛毛片免费看| av欧美777| 国产精品一区二区三区四区久久| 搡老妇女老女人老熟妇| 色尼玛亚洲综合影院| 亚洲在线自拍视频| 欧美最黄视频在线播放免费| h日本视频在线播放| 日本免费一区二区三区高清不卡| 精品国产三级普通话版| 国产亚洲欧美在线一区二区| 嫩草影院入口| 国产69精品久久久久777片| 精品一区二区三区av网在线观看| 日韩精品中文字幕看吧| 日韩免费av在线播放| 99久久精品热视频| 长腿黑丝高跟| 午夜老司机福利剧场| 亚洲欧美精品综合久久99| 国内精品久久久久精免费| av在线蜜桃| 久久人人精品亚洲av| 国产精品久久久久久精品电影| 丰满乱子伦码专区| 亚洲欧美一区二区三区黑人| 毛片女人毛片| 一卡2卡三卡四卡精品乱码亚洲| 怎么达到女性高潮| 久久久久免费精品人妻一区二区| 欧美日韩精品网址| 亚洲真实伦在线观看| 97超视频在线观看视频| 我要搜黄色片| 欧美成狂野欧美在线观看| 性欧美人与动物交配| 看免费av毛片| 久久国产精品影院| 亚洲精品一区av在线观看| 国产精品免费一区二区三区在线| 亚洲狠狠婷婷综合久久图片| 最后的刺客免费高清国语| 少妇的丰满在线观看| 国产激情偷乱视频一区二区| 国产午夜福利久久久久久| 高清日韩中文字幕在线| 一个人观看的视频www高清免费观看| 国产美女午夜福利| 一级毛片高清免费大全| 国产午夜福利久久久久久| 欧美午夜高清在线| 亚洲欧美日韩卡通动漫| xxxwww97欧美| 窝窝影院91人妻| 色哟哟哟哟哟哟| 亚洲熟妇熟女久久| 熟女人妻精品中文字幕| 可以在线观看毛片的网站| 老熟妇仑乱视频hdxx| 麻豆国产97在线/欧美| 国产真实乱freesex| 麻豆一二三区av精品| 国产乱人视频| 久久久国产精品麻豆| 久久久久性生活片| 窝窝影院91人妻| 国产极品精品免费视频能看的| 婷婷六月久久综合丁香| 亚洲黑人精品在线| 亚洲不卡免费看| 全区人妻精品视频| 日韩免费av在线播放| 亚洲av免费在线观看| 日本一二三区视频观看| 99热精品在线国产| 久久欧美精品欧美久久欧美| 久久久精品欧美日韩精品| 欧美3d第一页| 一个人看的www免费观看视频| 色综合亚洲欧美另类图片| 久久久久九九精品影院| 免费人成在线观看视频色| 中文字幕av成人在线电影| 乱人视频在线观看| 中出人妻视频一区二区| 国产探花极品一区二区| 搡女人真爽免费视频火全软件 | 久久精品国产自在天天线| 免费看十八禁软件| svipshipincom国产片| 欧美高清成人免费视频www| 色老头精品视频在线观看| 18+在线观看网站| 美女被艹到高潮喷水动态| 亚洲第一欧美日韩一区二区三区| 国产亚洲精品久久久com| 给我免费播放毛片高清在线观看| 99久久综合精品五月天人人| 国产精品1区2区在线观看.| 日韩欧美三级三区| 精品无人区乱码1区二区| 成人国产综合亚洲| 亚洲av免费在线观看| 午夜激情欧美在线| 757午夜福利合集在线观看| 黄片小视频在线播放| 黄色片一级片一级黄色片| 久久精品亚洲精品国产色婷小说| 久久久成人免费电影| 国产精品嫩草影院av在线观看 | 男人舔奶头视频| 久久6这里有精品| 亚洲欧美日韩高清在线视频| 婷婷精品国产亚洲av在线| av视频在线观看入口| 国产高清三级在线| 日本五十路高清| 午夜福利在线观看吧| 啦啦啦韩国在线观看视频| 99国产极品粉嫩在线观看| 欧美zozozo另类| 中文字幕久久专区| 欧美又色又爽又黄视频| 母亲3免费完整高清在线观看| 亚洲在线观看片| 色在线成人网| 少妇丰满av| 可以在线观看的亚洲视频| 日本成人三级电影网站| 观看美女的网站| 人妻丰满熟妇av一区二区三区| 一级a爱片免费观看的视频| 色精品久久人妻99蜜桃| 真人做人爱边吃奶动态| 最新中文字幕久久久久| 成人av在线播放网站| 日韩精品中文字幕看吧| 亚洲av成人av| 又爽又黄无遮挡网站| 亚洲精品在线美女| 搡老妇女老女人老熟妇| 久久久精品大字幕| 啦啦啦免费观看视频1| 久9热在线精品视频| www日本在线高清视频| 十八禁网站免费在线| 亚洲欧美精品综合久久99| 久久精品人妻少妇| 久久精品91无色码中文字幕| 一区二区三区国产精品乱码| 九九久久精品国产亚洲av麻豆| svipshipincom国产片| 在线十欧美十亚洲十日本专区| 国产高清视频在线观看网站| 亚洲欧美激情综合另类| 一个人看的www免费观看视频| 久久久久久久午夜电影| 校园春色视频在线观看| 欧美日韩黄片免| 最近最新中文字幕大全免费视频| 啦啦啦观看免费观看视频高清| 搡老岳熟女国产| 国产精品嫩草影院av在线观看 | 亚洲国产精品sss在线观看| 欧美乱妇无乱码| 99久久九九国产精品国产免费| 久久香蕉国产精品| 日韩欧美在线二视频| 99国产极品粉嫩在线观看| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 日韩欧美三级三区| 国产成人欧美在线观看| 两个人的视频大全免费| 色播亚洲综合网| 一区二区三区国产精品乱码| 精品国产美女av久久久久小说| 桃色一区二区三区在线观看| 91久久精品电影网| 亚洲黑人精品在线| 99久久99久久久精品蜜桃| 免费电影在线观看免费观看| 成年女人毛片免费观看观看9| 两个人的视频大全免费| 一个人观看的视频www高清免费观看| 真人做人爱边吃奶动态| 免费av毛片视频| 久久草成人影院| 嫩草影视91久久| 免费观看的影片在线观看| 亚洲国产精品合色在线| 午夜福利免费观看在线| 两个人看的免费小视频| 国产一区在线观看成人免费| 国产中年淑女户外野战色| 亚洲 欧美 日韩 在线 免费| 丝袜美腿在线中文| 在线观看舔阴道视频| 变态另类丝袜制服| 久久伊人香网站| 搡老岳熟女国产| 亚洲人成电影免费在线| 一个人看的www免费观看视频| 亚洲性夜色夜夜综合| 欧美最黄视频在线播放免费| 精品国产美女av久久久久小说| 国产精品久久久久久久久免 | 啦啦啦观看免费观看视频高清| 欧美在线黄色| 免费在线观看亚洲国产| 我的老师免费观看完整版| 国模一区二区三区四区视频| 亚洲最大成人中文| 久久精品综合一区二区三区| 久久久久久久亚洲中文字幕 | 日本五十路高清| 国产真实乱freesex| 熟女电影av网| 啪啪无遮挡十八禁网站| 欧美乱码精品一区二区三区| 老汉色av国产亚洲站长工具| 黄色丝袜av网址大全| 久久人人精品亚洲av| 午夜激情福利司机影院| 国产视频一区二区在线看| 男女之事视频高清在线观看| 国产伦精品一区二区三区四那| 久久精品国产亚洲av香蕉五月| 国产成人欧美在线观看| 中文字幕av成人在线电影| 久久久久国内视频| 两性午夜刺激爽爽歪歪视频在线观看| 19禁男女啪啪无遮挡网站| 一本精品99久久精品77| 亚洲片人在线观看| www国产在线视频色| 国产精品综合久久久久久久免费| a在线观看视频网站| 成年版毛片免费区| 亚洲av第一区精品v没综合| 日韩欧美在线乱码| 国产成人福利小说| 国产高清三级在线| 亚洲激情在线av| 熟女少妇亚洲综合色aaa.| 亚洲国产精品999在线| 51国产日韩欧美| 亚洲成人久久性| 久久午夜亚洲精品久久| 久久精品91蜜桃| 欧美成人一区二区免费高清观看| 亚洲最大成人中文| 亚洲狠狠婷婷综合久久图片| 可以在线观看毛片的网站| 国产欧美日韩精品亚洲av| 精品久久久久久久久久免费视频| 国产免费av片在线观看野外av| 亚洲男人的天堂狠狠| 无遮挡黄片免费观看| av片东京热男人的天堂| 亚洲专区中文字幕在线| av女优亚洲男人天堂| 成人特级黄色片久久久久久久| 国产精品久久视频播放| 亚洲av二区三区四区| 色播亚洲综合网| 国产午夜精品论理片| 极品教师在线免费播放| 一级毛片高清免费大全| 久久天躁狠狠躁夜夜2o2o| 成人高潮视频无遮挡免费网站| 亚洲av一区综合| 中文资源天堂在线| av在线天堂中文字幕| 国产熟女xx| 午夜福利免费观看在线| 国产激情偷乱视频一区二区| 久久久精品欧美日韩精品| 国产精品野战在线观看| АⅤ资源中文在线天堂| 亚洲在线观看片| 亚洲熟妇熟女久久| 国产欧美日韩精品一区二区| 成年女人永久免费观看视频| 一级毛片女人18水好多| 亚洲熟妇熟女久久| 成人欧美大片| 国产一区二区在线av高清观看| 99热精品在线国产| 有码 亚洲区| 午夜免费男女啪啪视频观看 | 九九热线精品视视频播放| 亚洲中文日韩欧美视频| 最近视频中文字幕2019在线8| 91九色精品人成在线观看| 国产精品爽爽va在线观看网站| 国产熟女xx| 日韩欧美免费精品| 又黄又爽又免费观看的视频| 国产精品久久久久久人妻精品电影| 少妇高潮的动态图| 中亚洲国语对白在线视频| 全区人妻精品视频| 无人区码免费观看不卡| 一个人看的www免费观看视频| 中亚洲国语对白在线视频| 小说图片视频综合网站| 国产精华一区二区三区| 国产精品综合久久久久久久免费|