亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        高粱粒重遺傳研究進(jìn)展

        2019-06-04 07:04:42韓立杰才宏偉
        生物技術(shù)通報(bào) 2019年5期
        關(guān)鍵詞:粒重同源高粱

        韓立杰 才宏偉

        (1. 河北廣播電視大學(xué),石家莊 050080;2. 中國(guó)農(nóng)業(yè)大學(xué)農(nóng)學(xué)院,北京 100193)

        高粱[Sorghum bicolor(L.)Moench]為禾本科、高粱屬、一年生草本植物,是人類最早栽培的重要谷類作物之一,其產(chǎn)量?jī)H次于玉米、水稻、小麥和大麥,居世界第五位。高粱具有抗旱、耐澇、耐鹽堿、耐瘠薄等特性,是部分干旱、半干旱、低洼易澇及瘠薄鹽堿地區(qū)人們的主要糧食和飼料來(lái)源[1]。近年來(lái),人口增長(zhǎng)、耕地面積減少、水資源短缺及氣候變化等問(wèn)題日益嚴(yán)峻,嚴(yán)重影響世界糧食安全問(wèn)題[2]。預(yù)計(jì)到2050年,全球糧食產(chǎn)量只有增長(zhǎng)70%才能滿足人們未來(lái)的糧食需求,因此培育高產(chǎn)的高粱品種將對(duì)緩解世界糧食安全問(wèn)題起到重要作用。粒重是影響高粱產(chǎn)量的三要素之一,增加粒重是提高高粱產(chǎn)量的有效途徑[3]。粒重是由數(shù)量性狀基因控制的,遺傳機(jī)理較復(fù)雜。20多年來(lái),眾多研究者利用數(shù)量性狀基因座(Quantitative trait locus,QTL)定位的方法來(lái)挖掘控制高粱粒重的相關(guān)基因,以期利用分子標(biāo)記輔助選擇的方法進(jìn)行育種,并揭示粒重基因的調(diào)控機(jī)理,取得了一定進(jìn)展。本文對(duì)高粱粒重的遺傳特性、高粱粒重的影響因素以及高粱粒重QTL定位研究進(jìn)展進(jìn)行了綜述,對(duì)已定位的高粱粒重QTLs進(jìn)行了比較分析,查找水稻和玉米中已克隆的粒重相關(guān)基因在高粱中的同源基因,并與高粱粒重QTLs定位區(qū)間進(jìn)行比較,從而為高粱高產(chǎn)的分子標(biāo)記輔助育種及高粱粒重主效QTLs的精細(xì)定位和克隆提供參考。

        1 高粱粒重的遺傳分析及影響因素

        高粱粒重與穗粒數(shù)和有效穗數(shù)共同構(gòu)成高粱單位面積產(chǎn)量的三要素,提高粒重是增加高粱產(chǎn)量的重要育種方向[3]。

        高粱粒重的遺傳力高,親本的粒重高低直接影響著F1代粒重的表現(xiàn)。盧慶善[4]認(rèn)為高粱千粒重性狀遺傳力高,是由少數(shù)主效基因控制的。此外,楊偉光等[5]和程寶成等[6]研究發(fā)現(xiàn)高粱粒重的遺傳符合加性-顯性模型,高值親本比低值親本遺傳傳遞力強(qiáng),并且雜交一代多數(shù)接近或超過(guò)高值親本,這對(duì)于高粱雜交育種親本材料的選用具有重要的指導(dǎo)意義。

        除遺傳因素外,環(huán)境也會(huì)對(duì)高粱粒重造成影響,進(jìn)而影響高粱產(chǎn)量。開(kāi)花前發(fā)生干旱會(huì)影響穗的大小、穗粒數(shù)和籽粒產(chǎn)量,開(kāi)花后發(fā)生干旱,則會(huì)使籽粒灌漿受阻,造成減產(chǎn)。研究表明,灌漿初、后期的溫度和光照會(huì)對(duì)高粱粒重造成影響[7]。另外高粱病蟲(chóng)害也會(huì)對(duì)高粱產(chǎn)量造成影響。除上述生物和非生物脅迫以外,Kiniry[8]和 Tao等[9]研究發(fā)現(xiàn)高粱千粒重和穗粒數(shù)呈負(fù)相關(guān)性。Murry等[10]對(duì)高粱籽粒和莖稈含糖量的遺傳平衡進(jìn)行研究發(fā)現(xiàn)莖糖產(chǎn)量和濃度的升高在一定程度上制約了籽粒產(chǎn)量的增加。

        綜上所述,高粱粒重主要是由遺傳因素決定的,并間接受到環(huán)境影響。開(kāi)展高粱粒重基因的定位和克隆研究對(duì)于進(jìn)行高粱分子標(biāo)記輔助育種具有重要意義。

        2 高粱粒重QTL定位研究進(jìn)展

        高粱為二倍體植物(2n=20),基因組小,其單倍體基因組DNA含量約為730 Mbp,這些特點(diǎn)使其成為研究甘蔗族及其他C4植物功能基因組學(xué)的模式植物[11]。隨著高粱全基因組測(cè)序的完成以及大量新型分子標(biāo)記的開(kāi)發(fā),完整覆蓋高粱全基因組的高密度遺傳連鎖圖譜已構(gòu)建完成,一些控制高粱粒重的QTLs也通過(guò)分子標(biāo)記定位于相應(yīng)的遺傳連鎖群上。

        2.1 高粱粒重QTL初步定位研究進(jìn)展

        近幾十年來(lái),越來(lái)越多的科學(xué)家利用QTL定位的方法來(lái)揭示作物復(fù)雜性狀的遺傳控制原理,以用于分子標(biāo)記輔助選擇育種或相關(guān)基因的克隆。目前已有來(lái)自16篇研究報(bào)道的108個(gè)控制高粱粒重的QTLs被定位到(表1),做圖群體有F2群體、F2:3群體、重組自交系群體和回交群體等,群體大小為90-370株,遺傳連鎖群數(shù)量為10-17個(gè),遺傳連鎖圖全長(zhǎng)為878 cM-2225.3 cM,標(biāo)記間的平均距離介于0.5 cM-14.8cM之間,具體定位結(jié)果匯總?cè)缦隆?/p>

        Paterson等[12]利用來(lái)自 BTx623×S.propinquum的F2群體定位到6個(gè)控制高粱粒重的QTLs,分別位于 1、2、3、6、8和10號(hào)染色體上。Brown等[13]用來(lái)自BTx623×IS3620C的119個(gè)重組自交系將3個(gè)高粱百粒重QTLs分別定位到4、8和9號(hào)染色體上,并在9號(hào)染色體上定位到一個(gè)高粱穗粒重QTL。Feltus等[14]用來(lái)自BTx623×IS3620C的重組自交系群體定位到3個(gè)控制高粱千粒重的QTLs(位于4、6和10號(hào)染色體),用來(lái)自BTx623×S. propinquum的F2群體定位到7個(gè)控制高粱千粒重的QTLs(位于1、2、3、4、6、7和8號(hào)染色體),解釋的表型變異分別為35%和49%,其中在不同群體中定位到的兩組QTLs(位于4號(hào)染色體上的QKwt.txs-D與Kwt.uga-F和位于6號(hào)染色體上的QKwt.txs-I與QKwt.uga-D)為相同的QTLs。Srinivas等[15]用來(lái)自296B×IS18551的168個(gè)重組自交系定位到3個(gè)高粱粒重QTLs,分別位于1、4 和6號(hào)染色體上,解釋的表型變異為7%-14.8%,來(lái)自大粒親本296B的等位基因增加了粒重。Phuong等[16]利用來(lái)自IS2449×IS1488的100個(gè)重組自交系分別在干旱脅迫下和充足灌溉的環(huán)境下定位到1個(gè)(位于5號(hào)染色體)和2個(gè)(位于2和5號(hào)染色體)高粱百粒重QTLs,解釋的表型變異為12.7%-18.9%,其中在不同環(huán)境下定位到的位于5號(hào)染色體的QTL為相同的QTL。Rajkumar等[17]利用來(lái)自E36-1×SPV570的184個(gè)重組自交系對(duì)高粱粒重QTL進(jìn)行定位,將4個(gè)粒重QTLs分別定位到2(2個(gè)QTLs)、3和4號(hào)染色體上,解釋的表型變異最高為7.82%。Shehzad和Okuno[18]利用Red Kafir×Takakibi組配的F2和F2:3家系對(duì)高粱產(chǎn)量相關(guān)性狀進(jìn)行QTL定位,結(jié)果在F2群體中定位到9個(gè)百粒重QTLs,分布于1、2、3、4、5和7號(hào)染色體上,其中位于2號(hào)染色體的qtl2GW解釋的表型變異最高,為36.3%,LOD值為33.8,利用F2:3家系能夠重復(fù)檢測(cè)到5個(gè)相同的QTLs,qtl2GW的LOD值最高,為33.8,說(shuō)明qtl2GW是一個(gè)控制高粱百粒重的主效QTL。Han等[19]利用來(lái)自 SA2313×Hiro-1的F2群體對(duì)高粱粒重QTL進(jìn)行定位,用種于北京的138株F2群體定位到3個(gè)QTLs,位于1、2和7號(hào)染色體上,1號(hào)染色體的qGW1解釋了最高的表型變異,為22%,并用種于海南的F2群體中47株極端小粒個(gè)體和47株極端大粒個(gè)體對(duì)定位結(jié)果進(jìn)行驗(yàn)證,定位到6個(gè)QTLs,其中位于1號(hào)和2號(hào)染色體的qGW1和qGW2能夠被重復(fù)檢測(cè)到,解釋的表型變異率分別為40%和27%,是主效的QTLs。此外,在ATx623×SA2313和K-385×SA2313的兩個(gè)F2群體中也檢測(cè)到了qGW1,證明qGW1在不同的遺傳背景中均能被檢測(cè)到。Sukumaran等[20]利用測(cè)序基因分型(Genotyping-by-sequencing,GBS)的方法對(duì)定位群體進(jìn)行基因型鑒定,結(jié)合高粱產(chǎn)量表型共定位到7個(gè)高粱產(chǎn)量QTLs,分布在1、2、3、4和6號(hào)染色體上,其中1號(hào)染色體的qYLD1.1在不同的種植環(huán)境下均能夠被檢測(cè)到,解釋的表型變異最高,為16%。Boyles等[21]利用兩個(gè)不同的定位群體定位到6個(gè)高粱千粒重QTLs,分布于1、2、5、6和8號(hào)染色體上,其中位于5號(hào)染色體上的QTL是最新發(fā)現(xiàn)的,解釋了21%的遺傳變異。Tao等[9]利用一個(gè)優(yōu)質(zhì)粒用高粱品種R931945-2-2和其野生祖先S. bicolorsubsp.Verticilliflorum組配的BC1F5群體定位到17個(gè)控制高粱千粒重的QTLs,其中有5個(gè)QTLs與前人的定位結(jié)果相同。Pereira等[22],Tuinstra等[23],Rami等[24],Murray 等[10]和 Upadhyaya等[25]也對(duì)高粱粒重QTL的定位也進(jìn)行了報(bào)道。近年來(lái),Zhang等[26]和 Boyles等[27]利用全基因組關(guān)聯(lián)分析法(Genome-wide association study,GWAS)分別檢測(cè)到10個(gè)和17個(gè)與高粱粒重相關(guān)聯(lián)的基因區(qū)域。根據(jù)以上研究結(jié)果發(fā)現(xiàn)高粱粒重相關(guān)QTL在10條染色體上均有分布。

        2.2 高粱粒重QTL精細(xì)定位研究進(jìn)展

        雖然對(duì)高粱粒重QTL初步定位研究已有很多報(bào)道,但高粱粒重QTL的精細(xì)定位或基因克隆研究卻寥寥無(wú)幾,遠(yuǎn)遠(yuǎn)落后于水稻、玉米等其他作物粒重相關(guān)基因的研究。

        Han等[19]首次對(duì)高粱粒重QTL精細(xì)定位進(jìn)行了研究,他們選取在不同的環(huán)境和遺傳背景下均能被檢測(cè)到的主效QTLqGW1為研究對(duì)象進(jìn)行精細(xì)定位,在F2群體中選取qGW1處為雜合,其它QTL位點(diǎn)處為隱性的大粒交換單株自交得到的F3群體為精細(xì)定位群體,將qGW1定位到1號(hào)染色體短臂標(biāo)記HLJ-10和HLJ-18之間大約101 kb的范圍內(nèi),在這個(gè)范圍內(nèi)包含13個(gè)候選基因,根據(jù)這13個(gè)基因在不同大粒和小粒品種間的測(cè)序結(jié)果,初步將Sobic.001G038300基因作為qGW1的候選基因,但沒(méi)作進(jìn)一步驗(yàn)證,韓立杰[28]還對(duì)另一個(gè)主效QTLqGW2進(jìn)行了精細(xì)定位,將其定位到2號(hào)染色體標(biāo)記Xtxp100和SB01474之間大約1 674 kb的范圍內(nèi)。

        Boyles等[21]在5號(hào)染色體發(fā)現(xiàn)了一個(gè)新的高粱千粒重QTL,解釋了較高的表型變異率為21%,該位點(diǎn)的優(yōu)勢(shì)等位基因來(lái)自小粒親本BTxARG-1,經(jīng)過(guò)進(jìn)一步檢測(cè),最高峰處的標(biāo)記位于Sobic.005G188400基因上,該基因編碼一個(gè)功能未知的的remorin蛋白,主要在早期花序組織和發(fā)育的胚中表達(dá);而水稻中的Grain setting defect1(GSD1)基因也編碼一個(gè)remorin蛋白,通過(guò)調(diào)控胞間連絲的通透性,進(jìn)而調(diào)控光合同化物的裝載和分配,使光合同化物由光合位點(diǎn)向韌皮部轉(zhuǎn)運(yùn)過(guò)程受阻所致,最終影響水稻種子結(jié)實(shí)率[29]。

        3 高粱粒重QTL的比較分析

        目前已有多個(gè)控制高粱粒重的QTLs被定位到,盡管研究者所用研究材料不盡相同,但有些QTLs可能是相同的,為了對(duì)這些定位結(jié)果進(jìn)行有效的比較和利用,Mace和 Jordan[30]與 Tao等[9]分別將1995-2010年和1995-2018年報(bào)道的高粱粒重相關(guān)QTLs整合到同一遺傳連鎖圖譜上。根據(jù)這些標(biāo)記在高粱遺傳連鎖圖上的距離及Li等[31]和Yonemaru等[32]發(fā)布的相關(guān)標(biāo)記在高粱染色體上的物理距離,將已定位高粱粒重相關(guān)QTLs在染色體上的位置進(jìn)行整理,見(jiàn)表2。

        表1 高粱粒重QTL定位信息

        從表2中可以看出,來(lái)自16篇研究報(bào)道的108個(gè)控制高粱粒重的QTLs被整合為50個(gè)QTLs。1號(hào)、3號(hào)和7號(hào)染色體上富含控制高粱粒重的基因區(qū)域。1號(hào)染色體涵蓋了數(shù)量最多的QTLs,占全部QTLs的 16%,qGW1-2、qGW1-4、qGW1-5和qGW1-64個(gè)QTLs均被兩篇及以上研究報(bào)道共同定位到,其中qGW1-2在5篇報(bào)道中均被定位到,解釋的遺傳變異最高為 22%[19],該 QTL 與 Mace和 Jordan[30]報(bào)道的元QTLQKWT_meta1.1為相同的QTL。3號(hào)和7號(hào)染色體上各有7個(gè)QTLs,各占全部QTLs的14%,但是僅有qGW3-3、qGW3-7和qGW7-5、qGW7-7被兩篇及以上研究報(bào)道共同定位到。2、4、6、8和10號(hào)染色體上各包含4個(gè)QTLs,其中2、4和10號(hào)染色體上各有3個(gè)QTLs,6和8號(hào)染色體上各有2個(gè)QTLs,在兩篇及以上研究報(bào)道中被定位到。5號(hào)染色體上的5個(gè)QTLs均為單獨(dú)定位到的QTL。9號(hào)染色體上的QTLs最少,僅有3個(gè),其中qGW9-2被兩篇研究報(bào)道共同定位到。在整合的50個(gè)QTLs中,有22個(gè)QTLs被兩篇及以上研究報(bào)道共同定位到,這些QTLs均能解釋較高的遺傳變異,并且研究者

        所用的材料均不盡相同,說(shuō)明高粱粒重相關(guān)QTL在不同遺傳背景和環(huán)境條件下一致表達(dá),具有重要的實(shí)際意義,不同報(bào)道中定位到的QTLs相互間也對(duì)定位結(jié)果的準(zhǔn)確性進(jìn)行了驗(yàn)證,這些QTLs和標(biāo)記將幫助育種者對(duì)高粱粒重性狀進(jìn)行遺傳改良,從而培育高產(chǎn)的高粱品種。

        表2 高粱粒重QTL 比較結(jié)果

        續(xù)表

        續(xù)表

        續(xù)表

        4 水稻和玉米已克隆粒重相關(guān)基因在高粱中的同源基因分析

        玉米、水稻、小麥、高粱和谷子等主要作物都是由其野生種馴化而來(lái),在漫長(zhǎng)的馴化過(guò)程中,一些共同的優(yōu)良特性,如果實(shí)或種子的增大、種子落粒性的喪失、種子休眠的打破和生育期的同步等,被人們有意或無(wú)意的選擇并保存下來(lái)[33]。這些相似表型變化的分子基礎(chǔ)有些是不同作物間同源基因?qū)е碌?。Lin等[34]在高粱中克隆的Sh1基因在水稻及玉米中的同源基因被證實(shí)可控制水稻及玉米的落粒,Liu等[35]在高粱中克隆的HD1基因被證實(shí)在高粱、水稻和谷子中受到平行選擇,水稻中克隆的粒重相關(guān)基因GW2和GS3在玉米中的同源基因均被證實(shí)影響著玉米籽粒大小和粒重[36-37]。因此,在谷物馴化過(guò)程中,某些馴化性狀基因的同源性和共線性被很好地保存下來(lái),可能存在功能保守的基因組區(qū)域,是各種谷物種子大小變化的基礎(chǔ)[12]。目前,水稻和玉米中已克隆了大量粒重相關(guān)基因,查找這些基因在高粱中的同源基因,并結(jié)合其在馴化過(guò)程中是否受到選擇和QTL定位結(jié)果,對(duì)于克隆高粱粒重相關(guān)基因并闡明其分子機(jī)理具有重要意義。

        Han等[19]和 Tao等[38]均查找了水稻和玉米中已克隆粒重相關(guān)基因在高粱中的同源基因,Tao等對(duì)這些基因在馴化過(guò)程中是否受到人工選擇也進(jìn)行了研究,現(xiàn)將這些基因信息進(jìn)行整理并與高粱粒重QTL整合信息進(jìn)行比較,見(jiàn)表3。將22個(gè)水稻粒重相關(guān)基因和8個(gè)玉米粒重相關(guān)基因的CDS序列與高粱基因組數(shù)據(jù)庫(kù)Sorghum bicolorv3.1.1(https://phytozome.jgi.doe.gov)進(jìn)行比對(duì),得到32個(gè)同源基因,分布在高粱10條染色體上。將同源基因在基因組上的位置和表2中整合的QTLs側(cè)翼標(biāo)記的物理距離進(jìn)行比較,發(fā)現(xiàn)18個(gè)同源基因位于高粱粒重QTL區(qū)間內(nèi),其中有12個(gè)基因具有馴化選擇信號(hào),即水稻小圓粒基因SRS5[39]在高粱中的同源基因Sobic.001G107100位于qGW1-2所在區(qū)間內(nèi);水稻粒長(zhǎng)基因GL3.1[40]的同源基因Sobic.001G154900位于qGW1-4所在區(qū)間內(nèi);水稻粒長(zhǎng)基因GL7[41]在高粱基因組中存在2個(gè)同源性較高的基因,Sobic.001G335800位于qGW1-7所在區(qū)間內(nèi),Sobic.002G367300位于qGW2-4所在區(qū)間內(nèi);水稻粒長(zhǎng)粒重基因GS3[42]和玉米粒重基因zmGS3[37]的同源基因?yàn)橄嗤幕騍obic.001G341700,位于qGW1-7所在區(qū)間內(nèi);水稻短?;騍G1[43]的同源基因Sobic.002G226500位于qGW2-3所在區(qū)間內(nèi);水稻粒重基因GW8[44]在高粱基因組中存在3個(gè)同源性較高的基因,Sobic.002G257900位于qGW2-3所在區(qū)間內(nèi),Sobic.003G406600位于qGW3-7所在區(qū)間內(nèi),Sobic.007G193500位于qGW7-7所在區(qū)間內(nèi) ;水稻小圓粒基因DEP2[45]/SRS1[46]的同源基因Sobic.002G374400位于qGW2-4所在區(qū)間內(nèi) ;水稻籽粒大小基因GS2[47]/GL2[48]的同源基因Sobic.004G269900位于qGW4-3所在區(qū)間內(nèi);水稻小粒基因D11[49]的同源基因Sobic.006G114600位于qGW6-3所在區(qū)間內(nèi)。

        以上這些同源基因位于前人定位的高粱粒重QTL所在區(qū)間內(nèi),并在高粱馴化過(guò)程中表現(xiàn)出選擇的特征,表明這些QTLs在高粱馴化過(guò)程中受到了影響,并且這些同源基因可能是調(diào)控高粱粒重的相關(guān)基因,在作物馴化過(guò)程中受到了平行馴化,但仍需要進(jìn)一步驗(yàn)證。

        5 高粱粒重基因的應(yīng)用及展望

        進(jìn)一步提高產(chǎn)量是高粱育種的重要目標(biāo),可通過(guò)傳統(tǒng)雜交育種和分子標(biāo)記輔助育種兩種途徑培育高粱高產(chǎn)新品種。由于高粱粒重超親遺傳現(xiàn)象較為普遍,親本粒重的高低影響著雜交一代粒重的表現(xiàn),因此在傳統(tǒng)雜交育種中,親本材料特別是母本(不育系)材料最好選用大粒材料。分子標(biāo)記輔助育種是快速培育優(yōu)質(zhì)作物品種的有效途徑,在水稻、玉米和小麥等作物中的應(yīng)用已取得很大進(jìn)展,通過(guò)QTL定位手段對(duì)高粱粒重關(guān)鍵基因進(jìn)行定位和克隆,可為育種家選育高產(chǎn)的高粱新品種提供分子標(biāo)記,加快育種速度,提高育種效率。

        表3 水稻、玉米已克隆粒重相關(guān)基因在高粱中的同源基因

        目前已定位到大量的高粱粒重相關(guān)QTLs,通過(guò)對(duì)不同研究者定位到的QTLs進(jìn)行比較,可對(duì)QTLs進(jìn)行相互驗(yàn)證,并能發(fā)現(xiàn)一些遺傳變異率高,在不同遺傳背景和環(huán)境下均能被檢測(cè)到的主效QTLs,與之連鎖的分子標(biāo)記可用于高粱高產(chǎn)分子標(biāo)記輔助育種,對(duì)這些主效QTLs可作進(jìn)一步的精細(xì)定位和克隆研究。

        在谷物馴化過(guò)程中,大粒的品種易受到人工選擇,不同谷物間功能保守的大粒性狀基因的同源性和共線性可以被很好的保存下來(lái)。查找水稻和玉米中已克隆粒重相關(guān)基因在高粱中的同源基因,研究其在馴化過(guò)程中是否受到選擇,并與QTL定位結(jié)果進(jìn)行比較,可發(fā)掘高粱粒重候選基因,對(duì)于克隆高粱粒重相關(guān)基因并闡明其分子機(jī)理具有重要意義。

        猜你喜歡
        粒重同源高粱
        藥食同源
        ——紫 蘇
        基于GBS測(cè)序和連鎖分析的藜麥單株粒重QTL定位
        種子(2023年9期)2023-11-22 13:10:56
        我終于認(rèn)識(shí)高粱了
        兩岸年味連根同源
        高粱名稱考釋
        高粱紅了
        以同源詞看《詩(shī)經(jīng)》的訓(xùn)釋三則
        離體穗培養(yǎng)條件下C、N供給對(duì)小麥穗粒數(shù)、粒重及蛋白質(zhì)含量的影響
        玉米自交系京92改良后代單穗粒重的雜種優(yōu)勢(shì)研究
        虔誠(chéng)書(shū)畫(huà)乃同源
        91日本在线精品高清观看| 久久精品夜色国产亚洲av| 国产精品久久久久影院嫩草| 在线观看av手机网址| 国产目拍亚洲精品二区| 国产色婷婷久久又粗又爽| 精品免费久久久久久久| 国产精品6| 免费美女黄网站久久久| 一区二区三区国产高清视频| 国产精品欧美一区二区三区| 日韩AV无码免费二三区| 亚洲处破女av一区二区| 女优一区二区三区在线观看| 精品国产一区二区三区av片| 欧美国产日本精品一区二区三区 | 无码爽视频| 毛片24种姿势无遮无拦| 色播中文字幕在线视频| 五月婷婷开心五月播五月| 性饥渴的农村熟妇| 精品国产一区二区三区久久久狼| 亚洲区精品久久一区二区三区女同| 蜜臀av一区二区三区久久| 国产精品久久久久久久妇| 国产无线乱码一区二三区| 亚洲产在线精品亚洲第一页| 精品人妻系列无码人妻漫画| 无码人妻黑人中文字幕| 国产精品福利久久香蕉中文| 人妻经典中文字幕av| 国产又色又爽又高潮免费视频麻豆| 国产欧美精品一区二区三区–老狼| 久久久精品国产亚洲av网不卡| 成年女人免费v片| 成人精品综合免费视频| 欧美巨大xxxx做受中文字幕 | 免费国产在线精品一区| 久久国产色av| 一区二区三区岛国av毛片| 亚洲成在人线视av|