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

    粗山羊草種質(zhì)遺傳多樣性及群體結(jié)構(gòu)的ISSR分析

    2023-08-20 07:59:07魏颯張?;?/span>王玉泉吳曉軍胡喜貴茹振鋼
    廣西植物 2023年7期
    關(guān)鍵詞:遺傳多樣性聚類分析

    魏颯 張?;? 王玉泉 吳曉軍 胡喜貴 茹振鋼

    摘 要:? 粗山羊草分布范圍廣,遺傳變異豐富,被認(rèn)為是改良普通小麥的重要基因源。為深入了解不同來源粗山羊草種質(zhì)的遺傳多樣性和群體結(jié)構(gòu),該研究利用ISSR分子標(biāo)記對56份粗山羊草種質(zhì)進(jìn)行了遺傳多樣性和群體結(jié)構(gòu)分析。結(jié)果表明:(1)16個ISSR引物共檢測170條多態(tài)性位點,每個ISSR引物多態(tài)性位點為3~18條,平均為10.63條;多態(tài)性信息(PIC)變異范圍為0.17~0.85,平均為0.67。(2)粗山羊草4個群體的遺傳多樣性比較顯示,中亞粗山羊草的群體遺傳多樣性水平最高(He=0.225 4,I=0.355 7),群體間的基因流較低(Nm=1.638 6)。(3)聚類結(jié)果在遺傳相似系數(shù)約0.67處,來源于塔吉克斯坦6份和土庫曼斯坦2份粗山羊草種質(zhì)材料聚成一類(Group 2);其他48份種質(zhì)材料形成一大類(Group 1),其中Group 1可進(jìn)一步分成3個Sub亞類,呈現(xiàn)出來源相同的粗山羊草種質(zhì)材料傾向聚在一起。(4)群體結(jié)構(gòu)分析將56份粗山羊草種質(zhì)分為5個群體,其中,來源于西亞伊朗V群體種質(zhì)材料遺傳背景比較一致,混雜程度相對較低;進(jìn)一步分析各群體Q值,發(fā)現(xiàn)IV群體種質(zhì)材料親緣關(guān)系的來源相對復(fù)雜,遺傳多樣最為豐富。該研究結(jié)果可為粗山羊草種質(zhì)親緣關(guān)系解析、種質(zhì)多樣性保護(hù)提供重要參考依據(jù),為其科學(xué)利用以及進(jìn)化研究奠定基礎(chǔ)。

    關(guān)鍵詞: 粗山羊草, 遺傳多樣性, 聚類分析, 群體結(jié)構(gòu), ISSR分析

    中圖分類號:? Q948

    文獻(xiàn)標(biāo)識碼:? A

    文章編號:? 1000-3142(2023)07-1317-09

    收稿日期:? 2022-08-26

    基金項目:? 河南省重大公益專項(201300110800); 河南省高校重點科研項目(21A210007)。

    第一作者: 魏颯(1998-),碩士研究生,研究方向為種質(zhì)資源創(chuàng)新與利用,(E-mail)1478441271@qq.com。

    通信作者:? 胡喜貴,博士,副教授,研究方向為作物種質(zhì)資源挖掘與利用,(E-mail)xiguihu@126.com。

    Genetic diversity and population structure of Aegilops

    tauschii accessions based on ISSR method

    WEI Sa1, ZHANG Haihui2, WANG Yuquan1, WU Xiaojun1, HU Xigui1*, RU Zhengang1

    ( 1. Henan Institute of Science and Technology, Xinxiang 453003, Henan, China; 2. Xinxiang

    Institute of Engineering, Xinxiang 453700, Henan, China )

    Abstract:? Aegilops tauschii is considered as an important gene source for improving common wheat, which has wide distribution and rich genetic variation. In order to understand genetic diverstity and population structure of A. tauschii from different origins, ISSR markers were used to evaluate genetic diversity and population structure of 56 A. tauschii accessions. The results were as follows:(1) The 170 polymorphic bands were detected by 16 ISSR primers, and polymorphic bands of each ISSR primer ranged from 3 to 18, with an average of 10.63. The variation of polymorphism information content(PIC)ranged from 0.17 to 0.85, with an averaged of 0.67. (2) The comparison among four populations indicated that the population genetic diversity of Central Asia was the highest (He=0.225 4, I=0.355 7) and gene flow between populations was relatively low (Nm=1.638 6). (3)The clustering results showed that 56 A. tauschii accessions were divided into two groups at the genetic similarity coefficient 0.67, of which eight accessions from Tajikistan and Turkmenistan were clustered in Group 2. And, the Group 1 including 48 accessions could be further divided into three sub-groups, which indicated that A. tauschii accessions with clustering together have the same origin. (4) Based on population structure analysis, 56 A. tauschii accessions were divided into five populations, of which the V population from Iran in West Asia had relatively consistent genetic background and relatively low degree of hybridization. Furthermore, the Q value analysis of populations showed that the genetic relationship of IV population were relatively complex, producing the most abundant genetic diversity. The results of this study can provide an important reference for analysis of genetic relationship, protection of biodiversity, and lay a foundation for the scientific utilization and evolution research of A. tauschii.

    Key words: Aegilops tauschii, genetic diversity, clustering analysis, population structure, ISSR analysis

    粗山羊草(Aegilops tauschii, DD, 2n=2x=14)別稱為節(jié)節(jié)麥,屬于禾本科小麥族(Triticeae)山羊草屬(Aegilops),被認(rèn)為是普通小麥(Triticum aestivum, AABBDD, 2n=6x=42)D基因組祖先的供體種(Lagudah et al., 1991; Dvorak et al., 2012)。粗山羊草主要分布在歐亞大陸中部,其起源中心位于里海南部海岸和阿塞拜疆,而后分別向東越過土庫曼斯坦的科彼特山脈擴(kuò)散到中國的新疆和黃河中部地區(qū),向西穿越土耳其東南的山谷傳播至敘利亞中部地區(qū)(Lubbers et al., 1991)。已有研究表明,粗山羊草D基因組具有豐富遺傳變異性,如抗逆性(Abbas et al., 2021)、抗病性(Olson et al., 2013)、優(yōu)異品質(zhì)基因(Hsam et al., 2001)等,已被作為改良普通小麥的重要種質(zhì)資源之一(Lagudah et al., 1991; Dvorak et al., 2012; 趙昕鵬等, 2019; 郜曉峰等, 2021)。

    核心種質(zhì)的收集和評價是種質(zhì)的應(yīng)用、創(chuàng)新及作物新品種選育的基礎(chǔ),尤其對供體物種遺傳多樣性和群體結(jié)構(gòu)的研究不僅能了解物種親緣關(guān)系,還可以為其合理高效的保護(hù)和利用提供理論基礎(chǔ)(Mourad et al., 2020)。粗山羊草的遺傳多樣性已從形態(tài)學(xué)、生理學(xué)、種子儲藏蛋白、同工酶、分子標(biāo)記等進(jìn)行廣泛研究(William et al., 1993; Ghasemzade et al., 2008; Mahjoob et al., 2021)。其中,因為分子標(biāo)記不受表型條件局限,可早期完成,所以被公認(rèn)為是研究遺傳多樣性和群體結(jié)構(gòu)高效、靈活的方法(Abouzied et al., 2013)。Lubbers等(1991)對不同來源102份粗山羊草25個位點的RFLP分析,發(fā)現(xiàn)來自里海的粗山羊草變異最大,阿富汗次之,而土耳其和巴基斯坦的變異最小,并進(jìn)一步支持粗山羊草起源于里海南部海岸??琢钭尩龋?998)利用RAPD標(biāo)記分析粗山羊草兩個亞種的基因組DNA多態(tài)性,表明Aegilops tauschii ssp. tauschii 多態(tài)性明顯高于A. tauschii ssp. strangulata。Pestsova等(2000)利用SSR標(biāo)記分析113份粗山羊草遺傳多樣性,發(fā)現(xiàn)來自高加索地區(qū)的粗山羊草多樣性比中亞地區(qū)的豐富,同時,將其劃分為兩大類且與地理分布一致。

    基于內(nèi)部簡單重復(fù)序列(inter-simple sequence repeats,ISSR)是一種顯性分子標(biāo)記,具有較高的多態(tài)性、重復(fù)性、穩(wěn)定性且操作簡便等諸多優(yōu)點,已廣泛應(yīng)用在埃及小麥(Egyptian wheat)(Abdel-Lateif & Hewedy, 2018)、硬粒小麥(Triticum durum)(Aslan-Parviz et al., 2020)、鉤刺山羊(Aegilops triuncialis)(Khodaee et al., 2021)、大麥(Hordeum vulgare)(張超等,2020)等小麥及野生近緣屬種的遺傳多樣性研究。目前,利用ISSR標(biāo)記研究粗山羊草遺傳多樣性的研究報道較少,僅李玉閣等(2017)利用9個ISSR標(biāo)記研究了75份中國粗山羊草的遺傳多樣性,結(jié)果將中國粗山羊草劃分為黃河流域和新疆兩大群體,并篩選到5份具有獨特變異的黃河流域類型,進(jìn)一步研究認(rèn)為黃河流域種群是從伊朗或土庫曼斯坦南部地區(qū)直接傳播而來(Wei et al., 2008; Su et al., 2020)。普通小麥的D基因組來源于有限粗山羊草類群,但多倍化和進(jìn)化“瓶頸”導(dǎo)致其遺傳基礎(chǔ)日益狹窄,與A、B基因組相比,小麥D基因組的遺傳多樣性尤其匱乏。為了利用不同來源粗山羊草類群拓展小麥D基因組,首先要解析粗山羊草種質(zhì)遺傳多樣性。因此,本研究利用16個ISSR標(biāo)記對不同來源的56份粗山羊草種質(zhì)開展遺傳多樣性和群體結(jié)構(gòu)研究,旨在了解它們的遺傳多樣性和群體結(jié)構(gòu),解析它們的親緣關(guān)系,以期為種質(zhì)多樣性保護(hù)和科學(xué)利用提供理論基礎(chǔ)。

    1 材料與方法

    1.1 試驗材料

    56份粗山羊草(Aegilops tauschii)種質(zhì)材料(表1),其中來源于西亞23份(伊朗18份、土耳其和阿塞拜疆各2份、格魯吉亞1份),中亞26份(阿富汗11份、土庫曼斯坦9份、塔吉克斯坦6份),南亞3份(巴基斯坦3份),東亞1份(中國1份),未知3份,均由美國種質(zhì)資源庫(National Genetic Resource Program, NPGS)饋贈,在河南科技學(xué)院小麥中心繁殖并保存。

    1.2 DNA提取

    按照Yan等(2002)2×CTAB的方法,從56份粗山羊草種質(zhì)材料的幼嫩葉片中提取基因組DNA,并經(jīng)Nanodrop 2000(美國賽默飛世爾)檢測DNA濃度和純度。

    1.3 ISSR-PCR擴(kuò)增分析

    參照Khodaee等(2021)研究,選用擴(kuò)增清晰且重復(fù)性好的16個ISSR引物(上海生物工程技術(shù)服務(wù)有限公司合成)對56份粗山羊草種質(zhì)材料進(jìn)行遺傳多樣性分析(表2)。

    ISSR-PCR擴(kuò)增體系20 μL:10 μL 2×Taq PCR StarMix(北京康潤誠業(yè)生物科技有限公司)、6 μL ddH2O、2 μL ISSR引物(10 pmol·μL-1)、2 μL DNA(50~60 ng·μL-1)。其PCR擴(kuò)增程序如下:94 ℃預(yù)變性5 min;94 ℃變性30 s,53.7~60.5 ℃(視不同引物而定)退化45 s,72 ℃延伸2 min,共35個循環(huán);72 ℃再延伸7 min,4 ℃繼續(xù)保存。PCR擴(kuò)增產(chǎn)物經(jīng)1.2%變性瓊脂糖凝膠電泳分離,D2000 bpDNA分子量標(biāo)準(zhǔn)(天根生化科技有限公司)作為對照,觀察結(jié)果,照相并保存。

    1.4 數(shù)據(jù)分析

    依據(jù)ISSR擴(kuò)增電泳圖,分別讀取數(shù)據(jù),相同遷移率處有帶為1,無帶為0,模糊不清楚為999,并建立二元數(shù)據(jù)矩陣。利用PowerMarker 3.2.5(Liu & Muse, 2005)軟件分析多態(tài)性信息含量(polymorphism information content,PIC);應(yīng)用Popgene 1.32(Yeh et al., 1999)軟件計算群體觀察到的等位基因(Na)和有效等位基因數(shù)(Ne)、Neis基因多樣性指數(shù)(He)和Shannons信息指數(shù)(I)、不同群體間的總基因多樣性(Ht)、群體內(nèi)的基因多樣性(Hs)、群體間的基因多樣性(Dst)、遺傳分化系數(shù)(Gst)和基因流(Nm);采用Ntsys 2.1(Rohlf, 2000)軟件,基于簡單匹配系數(shù)(simple matching coefficient,SM)進(jìn)行非加權(quán)平均法(UPGMA)聚類分析。

    將數(shù)據(jù)轉(zhuǎn)換成Structure 2.3.4(Pritchard et al., 2000)軟件的格式,進(jìn)行群體遺傳結(jié)構(gòu)分析。將K值設(shè)為2~20,基于模型每個K值進(jìn)行10次分析,再通過Structure Harvester(Earl & Vonholdt, 2012)確定最合適的K值,構(gòu)建遺傳結(jié)構(gòu)圖。同時,計算群體中種質(zhì)材料變異源于該群體中的概率Q值,用于分析其遺傳成分。

    2 結(jié)果與分析

    2.1 56份粗山羊草種質(zhì)ISSR引物多態(tài)性

    16個ISSR引物在56份粗山羊草種質(zhì)材料中擴(kuò)增出175條清晰穩(wěn)定的條帶,其中,ISSR13引物擴(kuò)增條帶最多(19條),ISSR23和ISSR28引物擴(kuò)增帶最少(3條),平均為10.94條。各ISSR引物的多態(tài)性條帶數(shù)為3~18條,平均為10.63條,多態(tài)性位點百分率為97.14%。ISSR引物PIC在0.17(ISSR28)~0.85(ISSR5和ISSR27)之間,平均為0.67。其中12個ISSR引物的PIC值高于平均值,占75%(表2)。Botstein等(1980)提出了衡量基因變異程度高低的PIC,并將PIC分成3個程度:高度(PIC>0.50)、中度(0.25≤PIC≤0.50)和低度(PIC<0.25)。本研究中16個ISSR引物PIC中,除引物ISSR3(PIC=0.49)、ISSR7(PIC=0.19)、ISSR23(PIC=0.39)和ISSR28(PIC=0.17)外,其他均表現(xiàn)為高度多態(tài)性信息。

    2.2 56份粗山羊草種質(zhì)群體遺傳多樣性和分化程度分析

    56份粗山羊草種質(zhì)材料按中亞、西亞、南亞和其他4個群體進(jìn)行遺傳多樣性分析,由表3可知,4個群體的觀測等位基因數(shù)(Na)、有效等位基因數(shù)(Ne)、Neis基因多樣性指數(shù)(He)與Shannons信息指數(shù)(I)由高到低排序為中亞>西亞>其他>南亞。這說明中亞粗山羊草的群體遺傳多樣性水平最高,西亞次之,其他和南亞的相對較低。進(jìn)一步分析群體基因的分化程度,由表4可知,56份粗山羊草種質(zhì)材料的總基因多樣性(Ht)為0.236 0,4個群體間發(fā)生一定的遺傳分化(Gst=0.233 8),而群體內(nèi)基因多樣性比值為0.766 2(Hs/Ht), 并存在中等程度的基因交流(Nm=1.638 6)。這說明粗山羊草種質(zhì)群體間的遺傳分化水平較低,大部分變異以群體內(nèi)為主。

    2.3 56份粗山羊草種質(zhì)的聚類分析

    利用Ntsys 2.1軟件對56份粗山羊草種質(zhì)材料的ISSR遺傳多樣性結(jié)果進(jìn)行聚類分析,可將52份粗山羊草種質(zhì)材料分開,占92.86%,其種質(zhì)材料間的配對遺傳相似系數(shù)(genetic similarity coefficient,GS)變化范圍為0.64~0.94(圖1)。由聚類結(jié)果(圖1)可知,在遺傳相似系數(shù)約0.67處,來源于塔吉克斯坦6份(PI662084、PI662095、PI662105、PI662106、PI662112和PI662116) 和土庫曼斯坦2份(PI662076和PI662078)種質(zhì)材料聚成一類(Group 2),其他48份粗山羊草種質(zhì)材料形成一大類(Group 1)。而在遺傳相似系數(shù)約0.70處,可進(jìn)一步將Group 1分成3個亞類:來源于8個國家的26份種質(zhì)材料組成Sub 1亞類,其中具有相同來源仍聚在一起,如土庫曼斯坦(PI662065、PI662066、PI662067、PI662068、PI662069、PI662070、PI662072)、阿富汗(CIae3、CIae4、CIae5、CIae6、PI220331、PI220641、PI220642、PI317392、PI317394、PI317398、PI511366)等;除3份種質(zhì)材料(土耳其PI486269和未知PI330489、PI369627)外,Sub 2亞類中種質(zhì)材料主要來源于伊朗(CIae8、CIae9、CIae10、CIae13、CIae15、CIae16、CIae17、CIae18、CIae19、CIae20、CIae21、PI268210、PI276985、PI511368、PI511378、PI511379、PI511382);Sub 3亞類種質(zhì)材料來源于阿塞拜疆(PI349037和PI428564)。同時,發(fā)現(xiàn)未知CIae51與阿富汗PI220642,未知PI330489和PI369627與伊朗PI276985聚在一起,可推測未知CIae51來源于阿富汗,而未知PI330489和PI369627來源于伊朗。本研究聚類分析結(jié)果將大部分粗山羊草種質(zhì)材料可聚類分開,說明供試粗山羊草種質(zhì)之間具有一定的遺傳差異,也證明了ISSR標(biāo)記在粗山羊草種質(zhì)分析方面具有較好的效果。

    2.4 56份粗山羊草種質(zhì)的遺傳結(jié)構(gòu)分析

    基于ISSR數(shù)據(jù),利用Structure 2.3.4軟件分析56份粗山羊草種質(zhì)群體的遺傳結(jié)構(gòu),繪制K與ΔK的關(guān)系圖(圖2),K=5時,ΔK最大。因此,可將56份粗山羊草種質(zhì)材料分成5個群體。由圖3可知,V群體的種質(zhì)材料主要來源于西亞的伊朗,遺傳背景比較一致,其混雜程度相對較低,而其他群體間混合程度相對偏高,尤其IV群體。從遺傳結(jié)構(gòu)整體分析也發(fā)現(xiàn),相同來源的種質(zhì)材料傾向聚在一起,這與聚類分析(圖1)的結(jié)果相似。

    由56份粗山羊草種質(zhì)材料在各群體中Q值分布(表5)可知,15份種質(zhì)材料Q值小于0.6,占26.8%,Q值大于0.8和0.9的種質(zhì)材料分別占50.0%和42.9%,這說明群體中大部分種質(zhì)材料的親緣關(guān)系比較單一。對各群體Q值大于0.8的進(jìn)行對比,發(fā)現(xiàn)V群體的Q值最大,為83.3%,而IV最小,僅為26.7%。這揭示了IV群體種質(zhì)材料親緣關(guān)系的來源相對復(fù)雜,遺傳多樣最為豐富。

    3 討論與結(jié)論

    3.1 ISSR標(biāo)記分析粗山羊草種質(zhì)遺傳多樣性

    近年來,分子標(biāo)記在生物多樣性和進(jìn)化關(guān)系方面研究中發(fā)揮著重要的作用,ISSR分子標(biāo)記是基于PCR技術(shù)發(fā)展起來的一種DNA多態(tài)性檢測技術(shù)。本研究利用16個ISSR標(biāo)記對56份粗山羊草種質(zhì)的遺傳多樣性進(jìn)行分析,結(jié)果表明多態(tài)性位點百分率為97.14%,觀測等位基因數(shù)為1.977 1,有效等位基因數(shù)為1.387 5,Neis基因多樣性指數(shù)為0.242 8,Shannons信息指數(shù)為0.382 6, 群體內(nèi)基因遺傳多樣性為0.180 9; 李玉閣等(2017)用9個ISSR標(biāo)記研究75份中國粗山羊草(新疆、河南和陜西)的遺傳多樣性,多態(tài)性位點百分率、觀測等位基因數(shù)、有效等位基因數(shù)、Neis基因多樣性指數(shù)、Shannons信息指數(shù)、群體內(nèi)基因遺傳多樣性分別為89.31%、1.993 4、1.471 5、0.273 4、0.415 2、0.236 7。經(jīng)比較同物種的兩個群體,本研究粗山羊草群體(主要來源西亞和中亞)除多態(tài)性位點百分率外,各項指標(biāo)均略低于中國粗山羊草群體,表明粗山羊草種質(zhì)從西亞和中亞向東傳入中國的新疆和黃河中部地區(qū)后,中國區(qū)域內(nèi)群體種質(zhì)發(fā)生部分變異,進(jìn)而提升了遺傳多樣性水平。另外,利用ISSR研究其他物種遺傳多樣性, 如鉤刺山羊引物多態(tài)性條帶平均為9.50條、PIC為0.30(Khodaee et al., 2021);硬粒小麥引物多態(tài)性條帶平均為10.00條、PIC為0.42(Aslan-Parviz et al., 2020);大麥引物多態(tài)性條帶平均為8.313條(張超等,2020)。與以上不同物種比較,粗山羊草ISSR引物多態(tài)性及PIC較高,這有可能與本研究粗山羊草種質(zhì)材料的地理來源廣有關(guān)。

    3.2 粗山羊草群體遺傳結(jié)構(gòu)

    Mizuno等(2010)和Sohail等(2012)對不同地理分布的粗山羊草種質(zhì)遺傳結(jié)構(gòu)進(jìn)行分析,將其分成兩個親緣譜系:L1(隸屬于Aegilops tauschii ssp. tauschii)和L2(隸屬于A. tauschii ssp. strangulata),并表明少數(shù)L2譜系的粗山羊草參與了普通小麥D基因組的起源;而未參與起源的L1譜系的粗山羊草具有豐富遺傳變異。Wang等(2013)進(jìn)一步按地理分布,將兩個親緣譜系細(xì)劃分成西(W)亞系和東(E)亞系,其中L1W亞系位于土耳其東部、 亞美尼亞、 阿塞拜疆和伊朗西部,L1E亞系位于伊朗中部到中國西部;L2W和L2E亞系分別位于亞美尼亞到阿塞拜疆區(qū)域和阿塞拜疆南部到伊朗北部的里海沿岸區(qū)域。Pestsova等(2000)利用SSR標(biāo)記分析不同地理來源的113份粗山羊草遺傳多樣性,發(fā)現(xiàn)高加索地區(qū)(格魯吉亞、亞美尼亞和俄羅斯的達(dá)吉斯坦地區(qū))的粗山羊草多樣性最為豐富,而中亞(烏茲別克斯坦和土庫曼斯坦)粗山羊草多樣性最低。本研究56份粗山羊草種質(zhì)材料按4個群體的遺傳多樣性對比,表明來自中亞(阿富汗、土庫曼斯坦和塔吉克斯坦)的粗山羊草群體遺傳多樣性水平最高(He=0.225 4,I=0.355 7),這與Pestsova等(2000)研究結(jié)果相反,這可能歸因于來自中亞阿富汗和塔吉克斯坦的粗山羊草種質(zhì)變異較大。王亞娟等(2010)基于SSR標(biāo)記對78份粗山羊草進(jìn)行聚類分析,在遺傳相似系數(shù)0.77處劃分成6個主要類群;而Sohail等(2012)基于DarT標(biāo)記對81份粗山羊草進(jìn)行聚類分析,將其分成3大類(Group A/B/C),Group C類主要來源伊朗。本研究聚類結(jié)果在遺傳相似系數(shù)約0.67處,56份粗山羊草種質(zhì)材料主要聚成2大類(Group 1和Group 2),其中Group 1進(jìn)一步分成3個Sub亞類;群體結(jié)構(gòu)將其分成5個群體。兩種分析方法在分類數(shù)量上略有差異,但均呈現(xiàn)出來源相同的粗山羊草種質(zhì)材料傾向聚在一起,本研究結(jié)果進(jìn)一步支持種質(zhì)材料間的親緣關(guān)系與其地理來源有一定的相關(guān)性(王亞娟等,2010; Sohail et al., 2012)。同時,本研究56份粗山羊草種質(zhì)材料中Q值小于0.6的15份,僅占26.8%,也進(jìn)一步說明大部分種質(zhì)材料的親緣關(guān)系比較單一,以群體內(nèi)或地理區(qū)域變異為主。因此,粗山羊草種質(zhì)在傳播過程中,群體內(nèi)種質(zhì)的基因交流或重組是造成不同群體粗山羊草差異的主要原因。

    3.3 粗山羊草保護(hù)生物學(xué)意義

    自然遺傳變異是植物生物學(xué)最重要的基礎(chǔ)資源之一(Koornneef et al., 2004)。普通小麥?zhǔn)撬谋扼w小麥和粗山羊草天然雜交后經(jīng)過染色體加倍而形成的,其起源、馴化及現(xiàn)代育種過程,排斥了供體物種的大量遺傳變異(Dubcovsky & Dvorak, 2007; Dvorak et al., 2012)。模擬普通小麥起源過程而創(chuàng)制的人工合成小麥,是重新引入供體物種遺傳變異的重要橋梁,在小麥育種中的應(yīng)用潛力很大(Hao et al., 2019)。但Waines(1998)研究發(fā)現(xiàn),許多粗山羊草物種的棲息地已經(jīng)遭受破壞甚至面臨危機(jī),尤其是在地中海東部沿岸地區(qū)。因此,從遺傳多樣性保護(hù)來看,了解粗山羊草種群結(jié)構(gòu)十分重要。本研究對不同來源的56份粗山羊草種質(zhì)進(jìn)行遺傳多樣性和群體結(jié)構(gòu)分析,了解它們的親緣關(guān)系,對每個種群或地區(qū)選取代表性的種質(zhì)保護(hù),可為今后粗山羊草自然遺傳變異研究與利用提供基礎(chǔ)。

    參考文獻(xiàn):

    ABBAS A, YU HY, CUI HL, et al., 2021. Genetic diversity and synergistic modulation of salinity tolerance genes in Aegilops tauschii Coss [J]. Plants, 10(7): 1393.

    ABDEL-LATEIF KS, HEWEDY OA, 2018. Genetic diversity among Egyptian wheat cultivars using SCoT and ISSR markers [J]. SABRAO J Breed Genet, 50(1): 36-45.

    ABOUZIED HM, ELDEMERY SMM, ABDELLATIF KF, 2013. SSR-based genetic diversity assessment in tetraploid and hexaploid wheat populations [J]. British Biotechnol J, 3(3): 390-404.

    ASLAN-PARVIZ M, OMIDI M, RASHIDI V, et al., 2020. Evaluation of genetic diversity of durum wheat (Triticum durum Desf.) genotypes using inter-simple sequence repeat (ISSR) and caat box-derived polymorphism (CBDP) markers [J]. Genetika, 52(3): 895-909.

    BOTSTEIN D, WHITE RL, SKOLNICK M, et al., 1980. Construction of a genetic linkage map in man using restriction fragment length polymorphisms [J]. Am J Hum Genet, 32(3): 314-331.

    DUBCOVSKY J, DVORAK J, 2007. Genome plasticity a key factor in the success of polyploid wheat under domestication [J]. Science, 316: 1862-1866.

    DVORAK J, DEAL KR, LUO MC, et al., 2012. The origin of spelt and free-threshing hexaploid wheat [J]. J Hered, 103(3): 426-441.

    EARL DA, VONHOLDT BM, 2012. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method [J]. Conser Genet Resour, 4(2): 359-361.

    GAO XF, ZHOU Y, SONG DY, et al., 2021. Construction of core collection of Chinese Aegilops tauschii Coss. germplasm resource based on spike morphological traits and molecular markers [J]. J Plant Genet Resour, 22(2): 361-370.? [郜曉峰, 周儀, 宋丹陽, 等, 2021. 基于穗形特征與分子標(biāo)記進(jìn)行中國節(jié)節(jié)麥核心種質(zhì)的創(chuàng)建 [J]. 植物遺傳資源學(xué)報, 22(2): 361-370.]

    GHASEMZADE R, BEHAMTA MR, AGHAII MJ, et al., 2008. Intra- and inter-population diversity of Iranian Aegilops tauschii based on seed storage protein electrophoresis [J]. Int J Agric Biol, 10(4): 463-465.

    HAO M, ZHANG LQ, ZHAO LB, et al., 2019. A breeding strategy targeting the secondary gene pool of bread wheat: introgression from a synthetic hexaploid wheat [J]. Theor Appl Genet, 132(8): 2285-2294.

    HSAM SLK, KIEFFER R, ZELLER FJ, 2001. Significance of Aegilops tauschii glutenin genes on breadmaking properties of wheat [J]. Cereal Chem, 78(5): 521-525.

    KHODAEE L, AZIZINEZHAD R, ETMINAN AR, et al., 2021. Assessment of genetic diversity among Iranian Aegilops triuncialis accessions using ISSR, SCoT, and CBDP markers [J]. J Genet Eng Biotechnol, 19: 5.

    KONG LR, DONG YC, JIA JZ, 1998. Random amplified polymorphism of DNA analysis in Aegilops tauschii [J]. Acta Bot Sin, 40(3): 223-227.? [孔令讓, 董玉琛, 賈繼增, 1998. 粗山羊草隨機(jī)擴(kuò)增多態(tài)性DNA研究 [J]. 植物學(xué)報, 40(3): 223-227.]

    KOORNNEEF M, ALONSO-BLANCO C, VREUGDENHIL D, 2004. Naturally occurring genetic variation in Arabidopsis thaliana [J]. Ann Rev Plant Biol, 55(1): 141-172.

    LAGUDAH ES, APPELS R, BROWN AHD, 1991. The molecular-genetic analysis of Triticum tauschii, the D-genome donor to hexaploid wheat [J]. Genome, 34(3): 375-386.

    LI YG, SU YZ, ZHANG DL, et al., 2017. Genetic diversity of Aegilops tauschii accessions native to China revealed by ISSR markers [J]. J Triticeae Crops, 37(1): 30-39.? [李玉閣, 蘇亞中, 張大樂, 等, 2017. 中國節(jié)節(jié)麥基于ISSR標(biāo)記的遺傳多樣性分析 [J]. 麥類作物學(xué)報, 37(1): 30-39.]

    LIU K, MUSE SV, 2005. PowerMarker: an integrated analysis environment for genetic marker analysis [J]. Bioinformatics, 21(9): 2128-2129.

    LUBBERS EL, GILL KS, COX TS, et al., 1991. Variation of molecular markers among geographically diverse accessions of Triticum tauschii [J]. Genome, 34(3): 354-361.

    MAHJOOB MMM, GORAFI YSA, KAMAL NM, et al., 2021. Genome-wide association study of morpho-physiological traits in Aegilops tauschii to broaden wheat genetic diversity [J]. Plants, 10(2): 211.

    MIZUNO N, YAMASAKI M, MATSUOKA Y, et al., 2010. Population structure of wild wheat D-genome progenitor Aegilops tauschii Coss.: implications for intraspecific lineage diversification and evolution of common wheat [J]. Mol Ecol, 19(5): 999-1013.

    MOURAD AMI, BELAMKAR V, BAENZIGER PS, 2020. Molecular genetic analysis of spring wheat core collection using genetic diversity, population structure, and linkage disequilibrium [J]. BMC Genom, 21(1): 1-12.

    OLSON EL, ROUSE MN, PUMPHREY MO, et al., 2013. Introgression of stem rust resistance genes SrTA10187 and SrTA10171 from Aegilops tauschii to wheat [J]. Theor Appl Genet, 126(10): 2477-2484.

    PESTSOVA E, KORZUN V, GONCHAROV NP, et al., 2000. Microsatellite analysis of Aegilops tauschii germplasm [J]. Theor Appl Genet, 101: 100-106.

    PRITCHARD JK, STEPHENS M, DONNELLY P, 2000. Inference of population structure using multilocus genotype data [J]. Genetics, 155(2): 945-959.

    ROHLF FJ, 2000. NTSYS-pc: numerical taxonomy and multivariate analysis system Version 2.1 [Z]. New York: Exeter Software.

    SOHAIL Q, SHEHZAD T, KILIAN A, et al., 2012. Development of diversity array technology (DArT) markers for assessment of population structure and diversity in Aegilops tauschii [J]. Breed Sci, 62(1): 38-45.

    SU YZ, ZOU MW, ZHU YM, et al., 2020.Analysis of population structure and origin in Aegilops tauschii Coss. from China through SNP markers [J]. Genet Resour Crop Evol, 67(4): 923-934.

    WAINES JG, 1998. In situ conservation of wild relatives of crop plants in relation to their history [M]//DAMANIA AB, VALKOUN J, WILCOX G, et al. The Origins of Agriculture and Crop Domesticatio. ICARDA, Apeppo, Syria: 300-306.

    WANG JR, LUO MC, CHEN ZX, et al., 2013. Aegilops tauschii single nucleotide polymorphisms shed light on the origins of wheat D-genome genetic diversity and pinpoint the geographic origin of hexaploid wheat [J]. New Phytol, 198(3): 925-937.

    WANG YJ, WANG CY, LIU XL, et al., 2010. Genetic diversity of Aegliops tauschii based on SSR marker [J]. J Agric Biotechnol, 18(3): 493-500.? [王亞娟, 王長有, 劉新倫, 等, 2010. 基于SSR標(biāo)記的粗山羊草遺傳多樣性分析 [J]. 農(nóng)業(yè)生物技術(shù)學(xué)報, 18(3): 493-500.]

    WEI HT, LI J, PENG ZS, et al., 2008. Relationships of Aegilops tauschii revealed by DNA fingerprints: the evidence for agriculture exchange between China and the West [J]. Prog Nat Sci, 18(12): 1525-1531.

    WILLIAM MDHM, PEA RJ, MUJEEB-KAZI A, 1993. Seed protein and isozyme variations in Triticum tauschii (Aegilops squarrosa) [J]. Theor Appl Genet, 87: 257-263.

    YAN ZH, WAN YF, LIU KF, et al., 2002. Identification of a novel HMW glutenin subunit and comparison of its amino acid sequence with those of homologous subunits [J]. Chin Sci Bull, 47(3): 222-226.

    YEH FC, YANG RC, BOYLE T, 1999. POPGENE version 1.32 Microsoft windows-based freeware for populations genetic analysis [Z]. Edmonton: University of Alberta.

    ZHANG C, Cairangzhuoma, KA MX, et al., 2020. Genetic diversity of naked barley analyzed by ISSR molecular markers [J]. Mol Plant Breed, 18(4): 1193-1201.? [張超, 才讓卓瑪, 卡毛先, 等, 2020. 應(yīng)用ISSR分子標(biāo)記分析裸大麥的遺傳多樣性 [J]. 分子植物育種, 18(4): 1193-1201.]

    ZHAO XP, ZHOU Y, L LL, et al., 2019. Genetic diversity of Aegilops tauschii Coss. and its utilization in improving common wheat [J]. Biotechnol Bull, 35(7): 181-189.? [趙昕鵬, 周云, 呂琳琳, 等, 2019. 節(jié)節(jié)麥遺傳多樣性及在改良普通小麥中的應(yīng)用 [J]. 生物技術(shù)通報, 35(7): 181-189.]

    (責(zé)任編輯 周翠鳴)

    猜你喜歡
    遺傳多樣性聚類分析
    茄子種質(zhì)資源農(nóng)藝性狀遺傳多樣性分析
    淺析田間水稻紋枯病抗性鑒定體系的確立與完善
    西藏野核桃的表型特征及其保育措施
    基于聚類分析研究貴州省各地區(qū)經(jīng)濟(jì)發(fā)展綜合評價
    商情(2016年39期)2016-11-21 08:45:54
    新媒體用戶行為模式分析
    農(nóng)村居民家庭人均生活消費(fèi)支出分析
    基于省會城市經(jīng)濟(jì)發(fā)展程度的實證分析
    中國市場(2016年33期)2016-10-18 12:16:58
    基于聚類分析的互聯(lián)網(wǎng)廣告投放研究
    科技視界(2016年20期)2016-09-29 12:32:48
    “縣級供電企業(yè)生產(chǎn)經(jīng)營統(tǒng)計一套”表輔助決策模式研究
    水稻紋枯病抗性鑒定體系的確立與遺傳多樣性研究
    亚洲男人的天堂狠狠| 一二三四社区在线视频社区8| 午夜福利高清视频| 999久久久精品免费观看国产| 又紧又爽又黄一区二区| 亚洲va日本ⅴa欧美va伊人久久| 91在线观看av| 99久久99久久久精品蜜桃| 怎么达到女性高潮| 九色国产91popny在线| 亚洲av电影不卡..在线观看| 又爽又黄无遮挡网站| 人妻久久中文字幕网| 精品日产1卡2卡| 亚洲第一电影网av| 91大片在线观看| 最近最新中文字幕大全电影3| svipshipincom国产片| 一二三四社区在线视频社区8| 国产三级在线视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品在线美女| 女警被强在线播放| 老鸭窝网址在线观看| 欧美成人一区二区免费高清观看 | 国产一区在线观看成人免费| 美女午夜性视频免费| 欧美黑人精品巨大| 五月伊人婷婷丁香| 久久这里只有精品中国| 国产精品av视频在线免费观看| 少妇被粗大的猛进出69影院| 久久精品国产亚洲av香蕉五月| 久久精品国产亚洲av香蕉五月| 丰满的人妻完整版| 亚洲精品色激情综合| 久久久水蜜桃国产精品网| 岛国视频午夜一区免费看| 日韩欧美一区二区三区在线观看| 男人舔女人下体高潮全视频| 搡老岳熟女国产| 免费看a级黄色片| ponron亚洲| 欧美日韩黄片免| ponron亚洲| 久久婷婷成人综合色麻豆| 色哟哟哟哟哟哟| 国产一级毛片七仙女欲春2| 国产精品综合久久久久久久免费| 欧美又色又爽又黄视频| 国产精品久久久久久亚洲av鲁大| 一级a爱片免费观看的视频| 亚洲成人中文字幕在线播放| 亚洲美女视频黄频| 岛国视频午夜一区免费看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品自产拍在线观看55亚洲| 亚洲av第一区精品v没综合| 高清毛片免费观看视频网站| 毛片女人毛片| 精品无人区乱码1区二区| 欧美成人免费av一区二区三区| 美女扒开内裤让男人捅视频| 伦理电影免费视频| 91老司机精品| 一级毛片高清免费大全| www日本在线高清视频| 精品少妇一区二区三区视频日本电影| 一进一出抽搐gif免费好疼| 久久热在线av| 精品一区二区三区av网在线观看| 久久精品国产综合久久久| 757午夜福利合集在线观看| 婷婷精品国产亚洲av| 日本黄大片高清| 久久久久久九九精品二区国产 | 俺也久久电影网| 久久久久国内视频| 制服诱惑二区| 又爽又黄无遮挡网站| 亚洲人成网站高清观看| 女人爽到高潮嗷嗷叫在线视频| 成人国产综合亚洲| 精品福利观看| 熟女电影av网| 日韩欧美国产一区二区入口| 欧美在线黄色| 黄色女人牲交| 19禁男女啪啪无遮挡网站| 一个人免费在线观看的高清视频| 亚洲av第一区精品v没综合| 国产成人精品久久二区二区91| 亚洲国产精品999在线| 欧美黑人精品巨大| 久久香蕉国产精品| 我要搜黄色片| 免费人成视频x8x8入口观看| tocl精华| x7x7x7水蜜桃| 少妇的丰满在线观看| 国内久久婷婷六月综合欲色啪| 最新美女视频免费是黄的| 中文在线观看免费www的网站 | 久久精品aⅴ一区二区三区四区| 丰满的人妻完整版| 国产午夜精品久久久久久| 欧美日韩乱码在线| 99国产精品一区二区蜜桃av| 一本一本综合久久| 国产视频内射| 成人18禁在线播放| 久久久久精品国产欧美久久久| 午夜精品一区二区三区免费看| 亚洲国产日韩欧美精品在线观看 | 亚洲五月天丁香| 一本一本综合久久| 久久草成人影院| 亚洲人成网站高清观看| 成人欧美大片| 757午夜福利合集在线观看| 亚洲精品在线美女| 亚洲成人精品中文字幕电影| 久久婷婷人人爽人人干人人爱| videosex国产| 母亲3免费完整高清在线观看| 亚洲人成伊人成综合网2020| 精品久久久久久久久久久久久| 国产激情偷乱视频一区二区| 日本熟妇午夜| 国产一区二区三区视频了| 午夜影院日韩av| 色噜噜av男人的天堂激情| 免费高清视频大片| cao死你这个sao货| 欧美日韩一级在线毛片| 特大巨黑吊av在线直播| 中文在线观看免费www的网站 | 十八禁网站免费在线| 叶爱在线成人免费视频播放| 亚洲精品久久国产高清桃花| 国产成人一区二区三区免费视频网站| 一卡2卡三卡四卡精品乱码亚洲| 国产激情欧美一区二区| 手机成人av网站| 日韩免费av在线播放| 国产精品av视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 这个男人来自地球电影免费观看| aaaaa片日本免费| 午夜久久久久精精品| 超碰成人久久| 老汉色∧v一级毛片| 国产真人三级小视频在线观看| 久久精品影院6| 亚洲aⅴ乱码一区二区在线播放 | 亚洲电影在线观看av| 亚洲真实伦在线观看| 国产成年人精品一区二区| 国产成人精品无人区| 亚洲成人久久爱视频| 欧美国产日韩亚洲一区| 精品高清国产在线一区| 女同久久另类99精品国产91| 亚洲美女黄片视频| 深夜精品福利| 欧美成人性av电影在线观看| 久99久视频精品免费| 精品久久蜜臀av无| 国产成人av激情在线播放| 少妇熟女aⅴ在线视频| 日韩免费av在线播放| 国产主播在线观看一区二区| 中文字幕熟女人妻在线| 黄片大片在线免费观看| 精品午夜福利视频在线观看一区| 成人精品一区二区免费| 亚洲午夜精品一区,二区,三区| 欧美又色又爽又黄视频| 亚洲,欧美精品.| 久久天躁狠狠躁夜夜2o2o| 久久这里只有精品19| 一本精品99久久精品77| 亚洲午夜精品一区,二区,三区| 国产精品影院久久| 亚洲av美国av| 亚洲av第一区精品v没综合| 可以在线观看的亚洲视频| 制服人妻中文乱码| av免费在线观看网站| 男女午夜视频在线观看| 亚洲国产欧洲综合997久久,| 岛国在线免费视频观看| 俺也久久电影网| 欧美3d第一页| 一级毛片精品| av免费在线观看网站| 精品人妻1区二区| 岛国视频午夜一区免费看| 色播亚洲综合网| 夜夜夜夜夜久久久久| 精品第一国产精品| 国产免费男女视频| 亚洲中文字幕一区二区三区有码在线看 | 亚洲精品美女久久av网站| 亚洲精华国产精华精| 美女 人体艺术 gogo| 正在播放国产对白刺激| 国产伦在线观看视频一区| 99久久久亚洲精品蜜臀av| 国产又色又爽无遮挡免费看| 久久精品综合一区二区三区| 国产av又大| 欧美丝袜亚洲另类 | 欧美一区二区国产精品久久精品 | 桃色一区二区三区在线观看| 啦啦啦韩国在线观看视频| 老司机深夜福利视频在线观看| 国产一区二区三区视频了| 一进一出好大好爽视频| 欧美日韩亚洲综合一区二区三区_| 亚洲欧美日韩高清专用| 中出人妻视频一区二区| 男男h啪啪无遮挡| 我的老师免费观看完整版| 久久精品国产亚洲av香蕉五月| 亚洲一码二码三码区别大吗| 国内精品久久久久久久电影| 国内精品久久久久精免费| 变态另类成人亚洲欧美熟女| 村上凉子中文字幕在线| 久久久久久国产a免费观看| 十八禁人妻一区二区| 午夜成年电影在线免费观看| 国产免费男女视频| 亚洲全国av大片| 婷婷六月久久综合丁香| 午夜亚洲福利在线播放| 日韩有码中文字幕| 国产又黄又爽又无遮挡在线| 婷婷精品国产亚洲av| 午夜福利欧美成人| 免费看美女性在线毛片视频| 国产熟女午夜一区二区三区| 国产亚洲欧美98| 久久久久国产一级毛片高清牌| 最近最新中文字幕大全电影3| 精品久久久久久久久久免费视频| 一级黄色大片毛片| 久久亚洲精品不卡| 一夜夜www| 精品一区二区三区四区五区乱码| 精品国内亚洲2022精品成人| 两个人的视频大全免费| 88av欧美| 亚洲五月婷婷丁香| 成人18禁在线播放| 777久久人妻少妇嫩草av网站| 亚洲美女黄片视频| 男女之事视频高清在线观看| e午夜精品久久久久久久| 日韩大尺度精品在线看网址| 久久婷婷成人综合色麻豆| 十八禁网站免费在线| 成年人黄色毛片网站| 50天的宝宝边吃奶边哭怎么回事| 国产黄片美女视频| 91麻豆精品激情在线观看国产| 十八禁人妻一区二区| 精品久久久久久久久久免费视频| 在线观看66精品国产| 久久精品国产亚洲av香蕉五月| 免费人成视频x8x8入口观看| 国产精品久久电影中文字幕| 日韩欧美在线乱码| 成人国产一区最新在线观看| videosex国产| 国产麻豆成人av免费视频| 国产1区2区3区精品| 每晚都被弄得嗷嗷叫到高潮| 国产午夜精品久久久久久| 在线观看日韩欧美| 一级毛片高清免费大全| 午夜福利在线在线| 90打野战视频偷拍视频| 麻豆成人午夜福利视频| 啦啦啦观看免费观看视频高清| 日本黄色视频三级网站网址| 亚洲18禁久久av| 欧美人与性动交α欧美精品济南到| 亚洲人成网站高清观看| 每晚都被弄得嗷嗷叫到高潮| 日本精品一区二区三区蜜桃| 91国产中文字幕| 精品第一国产精品| 欧美三级亚洲精品| 日韩欧美一区二区三区在线观看| 最近最新中文字幕大全免费视频| 久久久久国产一级毛片高清牌| 黄色成人免费大全| 国产av一区二区精品久久| 精品国产乱码久久久久久男人| 欧美中文日本在线观看视频| 老司机在亚洲福利影院| 两个人免费观看高清视频| 国模一区二区三区四区视频 | 亚洲男人的天堂狠狠| 日韩精品中文字幕看吧| 欧美又色又爽又黄视频| 国产亚洲精品久久久久5区| 久久久国产成人精品二区| 欧美乱码精品一区二区三区| 日本五十路高清| 在线观看美女被高潮喷水网站 | 黄频高清免费视频| 久久精品国产99精品国产亚洲性色| 国产午夜精品论理片| 国产真人三级小视频在线观看| 亚洲aⅴ乱码一区二区在线播放 | 亚洲自拍偷在线| 国产精品久久视频播放| 久久久久精品国产欧美久久久| 三级国产精品欧美在线观看 | 宅男免费午夜| 91国产中文字幕| 可以在线观看毛片的网站| 午夜福利高清视频| av在线天堂中文字幕| 亚洲av电影不卡..在线观看| 国产黄片美女视频| 欧美黑人欧美精品刺激| 男女视频在线观看网站免费 | 国产精品久久久久久精品电影| 亚洲精品久久国产高清桃花| 国产成人av激情在线播放| 国产一级毛片七仙女欲春2| 久久精品综合一区二区三区| 国产97色在线日韩免费| 少妇熟女aⅴ在线视频| 在线永久观看黄色视频| 国产精品久久视频播放| 99热这里只有是精品50| 老汉色∧v一级毛片| 国产一区二区三区在线臀色熟女| www国产在线视频色| 国产乱人伦免费视频| 1024香蕉在线观看| 麻豆国产av国片精品| 黄色毛片三级朝国网站| 亚洲专区中文字幕在线| 一夜夜www| 国产av一区在线观看免费| 日本黄色视频三级网站网址| 亚洲精品美女久久久久99蜜臀| xxxwww97欧美| 白带黄色成豆腐渣| 亚洲18禁久久av| xxxwww97欧美| 三级男女做爰猛烈吃奶摸视频| 日韩三级视频一区二区三区| 国内精品久久久久久久电影| 岛国在线观看网站| 1024香蕉在线观看| 一级作爱视频免费观看| 人人妻人人看人人澡| 亚洲乱码一区二区免费版| 最近最新免费中文字幕在线| 久久久精品大字幕| 亚洲自拍偷在线| 欧美成人一区二区免费高清观看 | 国产精品永久免费网站| 在线观看66精品国产| 黄色视频不卡| 两人在一起打扑克的视频| 白带黄色成豆腐渣| 国产一区二区激情短视频| 性欧美人与动物交配| 熟女电影av网| 99精品欧美一区二区三区四区| 高清毛片免费观看视频网站| 成人18禁高潮啪啪吃奶动态图| 亚洲专区国产一区二区| 国产精品一区二区三区四区久久| 后天国语完整版免费观看| 亚洲av片天天在线观看| 一区福利在线观看| 免费搜索国产男女视频| 特大巨黑吊av在线直播| 国内精品一区二区在线观看| 男女下面进入的视频免费午夜| www.www免费av| 国产午夜精品论理片| 国产精品一及| 女警被强在线播放| 免费无遮挡裸体视频| 久久久久久大精品| 国产av一区在线观看免费| 97碰自拍视频| 又大又爽又粗| 国产熟女xx| 少妇被粗大的猛进出69影院| 午夜老司机福利片| 麻豆成人av在线观看| 91成年电影在线观看| 日本在线视频免费播放| 黄色 视频免费看| av在线天堂中文字幕| 中文在线观看免费www的网站 | 国产精品美女特级片免费视频播放器 | 午夜福利视频1000在线观看| 久久久精品大字幕| 在线视频色国产色| 亚洲国产欧美一区二区综合| 国产1区2区3区精品| 美女免费视频网站| 桃红色精品国产亚洲av| 成人欧美大片| 国产av一区二区精品久久| 日韩有码中文字幕| 亚洲国产精品合色在线| 免费电影在线观看免费观看| 夜夜夜夜夜久久久久| 婷婷精品国产亚洲av| 老汉色av国产亚洲站长工具| 99久久综合精品五月天人人| 亚洲av电影不卡..在线观看| 麻豆久久精品国产亚洲av| av在线天堂中文字幕| 国产成人av激情在线播放| 在线视频色国产色| 亚洲一卡2卡3卡4卡5卡精品中文| 曰老女人黄片| 丰满人妻一区二区三区视频av | 国产乱人伦免费视频| 九色成人免费人妻av| 在线观看美女被高潮喷水网站 | 免费看a级黄色片| 中亚洲国语对白在线视频| 亚洲av成人一区二区三| 国产主播在线观看一区二区| 听说在线观看完整版免费高清| 色尼玛亚洲综合影院| 嫩草影视91久久| 日本 欧美在线| 国产亚洲精品综合一区在线观看 | 特大巨黑吊av在线直播| 三级国产精品欧美在线观看 | 丝袜美腿诱惑在线| 色综合欧美亚洲国产小说| 日韩欧美 国产精品| 男女视频在线观看网站免费 | 女生性感内裤真人,穿戴方法视频| 国产成人啪精品午夜网站| 午夜福利在线在线| 久久性视频一级片| 免费在线观看日本一区| 老司机靠b影院| 波多野结衣高清作品| 99久久无色码亚洲精品果冻| 国模一区二区三区四区视频 | 搡老岳熟女国产| 亚洲熟女毛片儿| 国产精品免费一区二区三区在线| 精品一区二区三区四区五区乱码| 一区二区三区激情视频| 精品免费久久久久久久清纯| 亚洲av中文字字幕乱码综合| 他把我摸到了高潮在线观看| 非洲黑人性xxxx精品又粗又长| 亚洲国产高清在线一区二区三| 成人午夜高清在线视频| 精品高清国产在线一区| 亚洲男人天堂网一区| 国产激情久久老熟女| 精品人妻1区二区| 99在线视频只有这里精品首页| 老熟妇乱子伦视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久久热在线av| avwww免费| 欧美中文日本在线观看视频| 少妇粗大呻吟视频| 亚洲精品久久国产高清桃花| 欧美激情久久久久久爽电影| 亚洲人成77777在线视频| 国内精品久久久久精免费| 男女做爰动态图高潮gif福利片| 免费在线观看成人毛片| 熟女电影av网| a级毛片在线看网站| 女人爽到高潮嗷嗷叫在线视频| 免费高清视频大片| 一本精品99久久精品77| 搡老岳熟女国产| 欧美3d第一页| 嫩草影视91久久| 欧美色视频一区免费| 小说图片视频综合网站| 日日干狠狠操夜夜爽| 欧美乱码精品一区二区三区| 国产亚洲欧美在线一区二区| 人妻丰满熟妇av一区二区三区| 久久精品国产亚洲av香蕉五月| 老司机午夜福利在线观看视频| 成人国产一区最新在线观看| 1024手机看黄色片| 999久久久精品免费观看国产| 亚洲一码二码三码区别大吗| 欧美乱码精品一区二区三区| 国产精品日韩av在线免费观看| 九色国产91popny在线| 亚洲美女视频黄频| 熟女少妇亚洲综合色aaa.| 国产一区在线观看成人免费| 99久久精品国产亚洲精品| 一进一出抽搐gif免费好疼| 国产成人精品久久二区二区91| 亚洲午夜理论影院| 国产精品一区二区三区四区久久| 正在播放国产对白刺激| av欧美777| 99久久综合精品五月天人人| 国产成人系列免费观看| 免费看日本二区| 欧美av亚洲av综合av国产av| 99热6这里只有精品| 999久久久国产精品视频| 国产av又大| 黄色 视频免费看| 桃红色精品国产亚洲av| 国语自产精品视频在线第100页| 日日干狠狠操夜夜爽| www日本在线高清视频| а√天堂www在线а√下载| 亚洲最大成人中文| 成人午夜高清在线视频| 久久久精品大字幕| 亚洲男人天堂网一区| 国产aⅴ精品一区二区三区波| 草草在线视频免费看| 精品日产1卡2卡| 美女大奶头视频| 久久久久久亚洲精品国产蜜桃av| 一区二区三区高清视频在线| 欧美日韩一级在线毛片| 亚洲中文字幕一区二区三区有码在线看 | 国产av在哪里看| 99热6这里只有精品| 亚洲全国av大片| 1024手机看黄色片| 两个人的视频大全免费| 亚洲一卡2卡3卡4卡5卡精品中文| 日韩免费av在线播放| 国产熟女午夜一区二区三区| 老司机靠b影院| av福利片在线| 久久精品人妻少妇| 18禁黄网站禁片午夜丰满| 我的老师免费观看完整版| 禁无遮挡网站| 999久久久精品免费观看国产| 麻豆国产av国片精品| 不卡一级毛片| 国产真实乱freesex| 精品电影一区二区在线| 可以在线观看的亚洲视频| 久久欧美精品欧美久久欧美| 亚洲欧美激情综合另类| 欧美又色又爽又黄视频| 夜夜夜夜夜久久久久| 桃色一区二区三区在线观看| 成人三级黄色视频| 精品欧美一区二区三区在线| 巨乳人妻的诱惑在线观看| 性欧美人与动物交配| 国产一区在线观看成人免费| 日韩欧美在线二视频| 免费电影在线观看免费观看| 精品第一国产精品| 日本精品一区二区三区蜜桃| 99精品在免费线老司机午夜| 亚洲av日韩精品久久久久久密| netflix在线观看网站| 国产伦人伦偷精品视频| 亚洲欧美日韩东京热| 男女视频在线观看网站免费 | 国产伦人伦偷精品视频| 欧美激情久久久久久爽电影| 中亚洲国语对白在线视频| 成人特级黄色片久久久久久久| 真人一进一出gif抽搐免费| 999久久久国产精品视频| x7x7x7水蜜桃| 国产成人系列免费观看| 99热这里只有精品一区 | 我的老师免费观看完整版| 日韩有码中文字幕| 最好的美女福利视频网| 亚洲第一欧美日韩一区二区三区| av在线播放免费不卡| 日韩欧美国产在线观看| 久久久久久久久久黄片| 免费一级毛片在线播放高清视频| 国内毛片毛片毛片毛片毛片| 不卡av一区二区三区| 男女视频在线观看网站免费 | 中国美女看黄片| 全区人妻精品视频| 好男人在线观看高清免费视频| 精品人妻1区二区| 色av中文字幕| 精品乱码久久久久久99久播| 成人一区二区视频在线观看| 国产一区二区三区在线臀色熟女| 超碰成人久久| 丝袜人妻中文字幕|