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

        ?

        園藝作物基因組測序研究進展

        2016-04-11 01:34:17倪曉鵬高志紅
        江蘇農(nóng)業(yè)科學(xué) 2016年2期
        關(guān)鍵詞:園藝作物基因組學(xué)

        倪曉鵬++高志紅

        摘要:回顧了園藝作物全基因組測序的發(fā)展歷程,介紹3代測序技術(shù)的特點和應(yīng)用現(xiàn)狀??偨Y(jié)了葡萄、番木瓜、草莓等果樹,黃瓜、白菜、番茄等蔬菜以及蓮、康乃馨等花卉在內(nèi)的25種園藝作物基因組測序簡況。重點介紹了全基因組測序在基因注釋、比較基因組學(xué)研究、重測序、全基因組關(guān)聯(lián)分析、轉(zhuǎn)錄組學(xué)、SNP芯片開發(fā)等方面的應(yīng)用。最后討論了全基因測序研究的難點和今后的研究方向。

        關(guān)鍵詞:園藝作物;基因組測序;基因組學(xué)

        中圖分類號: Q78文獻標志碼: A 文章編號:1002-1302(2016)02-0009-04

        園藝業(yè)是農(nóng)業(yè)種植業(yè)的重要組成部分,對豐富人類營養(yǎng)和美化、改善人類生存環(huán)境有重要意義。我國是世界園藝大國,園藝是我國農(nóng)業(yè)的重要組成部分,其中蔬菜播種面積由1990年的600萬hm2增加到2013年的2 100萬hm2,產(chǎn)量由1990年的1.95億t增加到2013年的7.35億t[1]。種植面積和產(chǎn)量分別占世界的43%和49%,居世界第一位。隨著動植物全基因組測序的不斷發(fā)展,2001年由美國能源部推動的人類基因組工作草圖的發(fā)表被認為是人類基因組計劃成功的里程碑[2]。破譯人類遺傳信息,將對醫(yī)學(xué)、生物學(xué)乃至整個生物科學(xué)產(chǎn)生不可估量的影響。這一計劃促成了大批動植物全基因測序的完成,也有望從基因組水平上分析基因的結(jié)構(gòu)、組成、調(diào)控和物種進化,從而避免傳統(tǒng)分子生物學(xué)帶來的種種弊端,大大促進全基因組測序技術(shù)在園藝作物上的應(yīng)用發(fā)展。

        1全基因測序的發(fā)展歷程及其在園藝植物測序中的應(yīng)用

        距世界第一例模式植物擬南芥全基因組測序發(fā)表已有十余年[3],在此期間測序技術(shù)飛速發(fā)展,每年發(fā)表的測序文章都穩(wěn)定增加。擬南芥的測序就是基于第一代測序技術(shù)Sanger測序完成的,運用了鏈終止和斷裂技術(shù)。第一代測序技術(shù)已經(jīng)規(guī)?;?,且具有測序讀長較長、測序準確率高等特點,但是由于其時間長、成本高、通量低等缺點而無法滿足現(xiàn)在實驗需求。

        2004年,美國國家人類基因組研究所(NHGRI)發(fā)起了一項融資計劃[4],目標是在未來10年將人類基因組測序的費用減少到1 000美元,這刺激了商業(yè)化的新一代測序(next generation sequencing,NGS)技術(shù)的快讀萌芽。第二代測序技術(shù)大幅提升基因組測序的輸出與成本比,相較第一代測序技術(shù)涵蓋范圍更加廣泛,可以同時測定多個平行DNA片段,從而輸出更多測序閱讀量,但一般長度較短、質(zhì)量低。第一個NGS技術(shù)是羅氏公司(Roche)公司的454焦磷酸測序法[5],Solexa和Illumina測序平臺在1年后也進行商業(yè)化應(yīng)用。第二代測序技術(shù)具有成本低、單次數(shù)據(jù)量大、消耗時間短等特點,故又被稱為第二代高通量測序技術(shù),并漸漸成為大規(guī)模全基因組測序技術(shù)的主導(dǎo)。

        Roche 454測序平臺是基于焦磷酸測序,由于低數(shù)據(jù)率和相對較短的讀取,最初適用于細菌基因組測序,隨后技術(shù)改進,利用Roche 454與Sanger測序結(jié)合起來運用于更復(fù)雜的基因組測序,完成蘋果基因組測序[6],取代Sanger測序作為主要數(shù)據(jù)源完成可可基因組[7]和甜瓜[8]基因組測序。Illumina 公司提供的Solexa系統(tǒng)比Roche 454系統(tǒng)數(shù)據(jù)產(chǎn)出量多且花費較少的費用。自推出以來,讀取長度獲得明顯改善,并成功運用到黃瓜的測序工作中[9]。

        野生草莓[10]作為第一個采用Roche 454、Illumina、SOLiDTM三大平臺共同完成測序的植物,標志著測序技術(shù)向多平臺合作而不是獨立運行的模式發(fā)展。近幾年來,Illumina測序成為第二代測序平臺的主導(dǎo),截至2014年11月,已為多種園藝作物如中國白菜[11]、馬鈴薯[12]、香蕉[13]、橙[14]和西瓜[15]基因組測序(表1)。

        近年來,新一代測序平臺已經(jīng)出現(xiàn)并被稱為“第三代測序技術(shù)”,與第二代測序方法相比有著進一步模式轉(zhuǎn)變,有2個突出特點,一是測序前不再需要PCR擴增,二是熒光或電流信號都可以在互補鏈加核苷酸的酶反應(yīng)中被檢測到。2種測序平臺現(xiàn)在已經(jīng)投入商業(yè)化運營,分別是Pacific Biosciences公司的PacBio RS平臺(http://www. Pacificbio scie nces.com)和Ion Torrent 公司的Personal Genome平臺 (Life Technologies,http://www.iontorren t.com)。PacBio測序平臺使用實時零模式波導(dǎo)檢測單個DNA聚合酶的活性[31],大規(guī)模并行單分子實時(SMRT)測序保證了高通量測序。其突出特點是序列讀長,經(jīng)報道其單程長讀取對的準確率為84.2%~97.8%[32],相較于第二代測序技術(shù),第三代測序技術(shù)可以讀取更長序列、提高測序通量、節(jié)省試劑成本,但其錯誤率比第二代測序技術(shù)高,3代測序技術(shù)優(yōu)缺點比較見表2。

        2全基因組測序的應(yīng)用

        2.1基因注釋

        基因注釋是利用生物信息學(xué)方法給測序完成的物種序列附上生物學(xué)信息的過程,通過識別不編碼蛋白質(zhì)的基因片段,識別基因上的元素(基因預(yù)測),給元素附上生物學(xué)信息的手段進行重復(fù)序列的識別、非編碼RNA的預(yù)測、基因結(jié)構(gòu)預(yù)測和基因功能注釋,發(fā)現(xiàn)與農(nóng)藝性狀相關(guān)基因,例如花期調(diào)控、植物生長習(xí)性、耐寒性、休眠、果實性狀與品質(zhì)等。Fukuoka等通過鑒定蘋果(Malus×domestica)的146個MADS-box基因,與擬南芥和水稻MADS-box基因聚成6個亞組,預(yù)測MADS-box基因在17條染色體上的密度,在根、莖、葉、花組織和果實發(fā)育的5個階段進行了分析,表明MADS-box基因參與了蘋果的生理和發(fā)育過程的各個方面[33]。Yu等將超過6 000 個茄子組織和處理組的cDNA克隆組成1個含有16 245 個獨特序列的單基因組,對單基因組進行測序并采用SAGE策略對茄子單基因集在轉(zhuǎn)錄組中的功能分析,相當數(shù)量的短序列標簽被成功注釋[34]。

        2.2比較基因組學(xué)研究

        比較基因組學(xué)研究是基于基因組圖譜和測序技術(shù),對已知的基因和基因組結(jié)構(gòu)進行比較以了解相關(guān)基因的功能、表達機制和物種親緣關(guān)系。

        番茄測序小組完成了對栽培番茄及其近緣野生種醋栗番茄全基因組的精細序列分析。在解碼的番茄基因組中共鑒定出 34 727 個基因,其中97.4% 的基因已經(jīng)精確定位到染色體上。通過比較基因組分析發(fā)現(xiàn)了番茄果實進化和發(fā)育的基因組學(xué)基礎(chǔ),番茄基因組經(jīng)歷的2次三倍化使基因家族產(chǎn)生了特異控制果實發(fā)育及營養(yǎng)品質(zhì)的新成員[19]。Zhang等研究發(fā)現(xiàn)在梅花基因組中存在和抗病相關(guān)的PR基因家族,由PR編碼的蛋白質(zhì)可介導(dǎo)植物產(chǎn)生防御反應(yīng),以此來抵抗病害和不良環(huán)境。例如:梅花基因組中存在的[WTBX][STBX]PR10[WTBZ][STBZ]與梅花的葉和根抗鹽度、抗干旱和抗病菌的機制相關(guān)[21]。

        2.3重測序

        通過基因組重測序,可以對栽培種和野生種基因組之間的差異進行比較,從而揭示物種起源以及馴化過程,為鑒定有價值的遺傳資源以及園藝作物育種提供重要參考。

        Qi等對115個黃瓜品系進行了深度重測序,鑒定出112個假定的馴化延伸,其中1個區(qū)域包含著參與黃瓜果實苦味消失過程的1個基因,以及1個重要的馴化特性,也調(diào)查了栽培品種中分離度的基因組基礎(chǔ),在β-胡蘿卜素羥化酶發(fā)現(xiàn)了一個自然遺傳變異,可以用來培育含有更高營養(yǎng)價值的黃瓜,揭露黃瓜進化的遺傳歷程,為未來遺傳育種提供了保證[35]。

        2.4全基因組關(guān)聯(lián)分析

        全基因組關(guān)聯(lián)分析(genome-wide association study,GWAS)是應(yīng)用基因組中數(shù)以百萬計的單核苷酸多態(tài)性(single nucleotide ploymorphism,SNP)作為分子遺傳標記,進行全基因組水平上的對照分析或相關(guān)性分析,通過比較發(fā)現(xiàn)影響復(fù)雜性狀的基因變異的一種新策略。隨著基因組學(xué)研究以及基因芯片技術(shù)的發(fā)展,人們已通過GWAS方法發(fā)現(xiàn)并鑒定了大量與復(fù)雜性狀相關(guān)聯(lián)的遺傳變異。近年來,這種方法在農(nóng)業(yè)動物重要經(jīng)濟性狀主效基因的篩查和鑒定中得到了應(yīng)用。

        Kenta等為了研究栽培番茄遺傳和表型多樣性的關(guān)系,重測序6個栽培番茄品系,基于全基因連鎖不平衡分析表明,連鎖不平衡的延長依賴于染色質(zhì)的性質(zhì),全基因組關(guān)聯(lián)分析表明所鑒定的SNP與農(nóng)藝重要的性狀有著密切的關(guān)系,3個基因(FAS、SP、U)被發(fā)現(xiàn)與性狀變異相關(guān),證明通過重測序得到的大量SNP進行全基因組關(guān)聯(lián)分析可以鑒別特異性位點,表明栽培番茄在基因組學(xué)和遺傳學(xué)發(fā)展上又邁進了一大步[36]。

        2.5與轉(zhuǎn)錄組數(shù)據(jù)比較分析

        轉(zhuǎn)錄組研究是基因功能及結(jié)構(gòu)研究的基礎(chǔ)和出發(fā)點,是全基因組測序完成后首先要面對的問題。最近科學(xué)家們將高通量測序技術(shù)應(yīng)用于轉(zhuǎn)錄組分析開發(fā)出了RNA測序技術(shù)(RNA-SEq),該技術(shù)能夠在全基因組范圍內(nèi)檢測基因表達情況,進行差異基因篩選分析。由于RNA-SEq技術(shù)具有通量高、可重復(fù)性好、檢測范圍寬、定量準等特點,已經(jīng)廣泛應(yīng)用于細菌、擬南芥、水稻和人類等生物轉(zhuǎn)錄組的研究。

        Zenoni等利用RNA-Seq技術(shù)對葡萄栽培種Corvina漿果發(fā)育3個階段的轉(zhuǎn)錄組變化進行了研究,對比Pinot Noir 40024參考基因組,分析測量基因表達水平,共檢測到17 324個基因在果實發(fā)育過程中表達,其中6 695個在不同階段會特異性表達,這表明了基因功能的多樣性與表達特異性以及轉(zhuǎn)錄組的高度復(fù)雜性[37]。

        2.6SNP芯片的開發(fā)

        單核苷酸多態(tài)性(SNP)根據(jù)其在基因中的位置,可以分為基因編碼區(qū)、基因非編碼區(qū)、基因間隔區(qū)(基因之間的區(qū)域),在基因組中是多態(tài)性最豐富的DNA分子標記,具有數(shù)量眾多、分布均勻、分型方便等特點,它可以識別遺傳變異和關(guān)聯(lián)的表型基因分型。

        Chagné等利用生物信息學(xué)工具,從包含蘋果品系中350 000 條序列的EST數(shù)據(jù)庫中開發(fā)SNP標記,結(jié)果識別了71 482個假定SNP分子標記。設(shè)計了464個PCR引物對,對PCR產(chǎn)物進行測序,重新獲得的SNP標記映射到蘋果參考基因(Royal Gala×A689-24雜交系和Malling 9×Robusta 5)SNP基因分型采用高分辨率熔解 (HRM)技術(shù),共包含210個編碼SNP的93個新標記被成功映射到參考圖譜。此方法為使用數(shù)量性狀定位(QTL)技術(shù)了解重要農(nóng)藝性狀和連鎖不平衡分析的基因調(diào)控提供了借鑒,也成為物理和遺傳圖譜結(jié)合的有效標記[38]。

        3全基因組測序存在的挑戰(zhàn)

        測序任何基因組的挑戰(zhàn)包括倍性和雜合度。與動物的基因組相比,植物基因組具有巨大的基因組大小、高度重復(fù)序列和全基因組或者片段基因組復(fù)制。

        第二代測序平臺Roche 454、 Illumina和SoLiD不僅吞吐量顯著增加,還減少了錯誤率,增加了閱讀長度,使大范圍物種測序都能采用測序平臺經(jīng)濟地進行。然而,植物基因組的獨特性給測序帶來了許多挑戰(zhàn),如他們的重復(fù)特性,裝配完整基因組是具有挑戰(zhàn)性的,這是由于植物基因組轉(zhuǎn)座因子的高拷貝數(shù)和擴張?zhí)匦詻Q定的。例如,玉米基因組測序通過BAC-by-BAC方法,基因組中85%含有轉(zhuǎn)座元件測序[39]。植物基因組中全基因組、節(jié)段和串聯(lián)重復(fù)的頻率[40]也給旁系同源種類帶來組裝問題,如果是最近一次復(fù)制,其序列同源性會高。因此,盡管第二代測序平臺能夠測序一大部分基因組,但迄今為止的全基因組測序方案還是缺乏足夠的代表性,質(zhì)量評估在組裝過程中還會丟失部分基因。

        另一個問題是,不是所有植物品種都是純合二倍體自交系,盡管使用WGS方法為雜合組織測序是可行的。然而高雜合度,包括交叉和無性繁殖的物種,如馬鈴薯,會阻礙WGS組裝。為了解決葡萄基因組(高度雜合二倍體)問題,Jaillonet等[JP3]采用了高度近交的品種,減少了雜合程度,進而才完成了測序。使用長閱讀可以提高一個基因組中組裝分離單體型的能力[16]。

        植物倍性也是在植物denovo測序和組裝的一個巨大挑戰(zhàn),其結(jié)果取決于物種是一個同源多倍體是異源多倍體,迄今為止,所有多倍體測序依賴于倍性降低或者染色體的物理分離。例如,絕大多數(shù)馬鈴薯品系是四倍體,最初是取雜合二倍體馬鈴薯基因 (RH89-039-16)進行測序從而避免同源染色體之間的高度雜合性,最后利用1個特殊的馬鈴薯基因型,1個加倍的單倍體 (DM1-3 516 R44),含有純合的12條染色體作為參考馬鈴薯基因組[12]。栽培草莓(Fragaria ananassa)是來自4個不同祖先的異源八倍體,野生草莓(Fragaria vesca)是二倍體,被用來測序是為了避免測序多個基因組帶來的困難[10]。

        4全基因組測序的展望

        為了面對世界基因組測序技術(shù)發(fā)展如此迅猛的挑戰(zhàn),大量全基因組測序數(shù)據(jù)亟待深度挖掘,應(yīng)制定長遠的戰(zhàn)略性基因組學(xué)研究計劃,不僅局限于栽培物種,更應(yīng)該深度開發(fā)我國重要野生近緣物種的測序,促進重要基因資源的挖掘、保護和利用。擺脫先前依賴外觀表現(xiàn)型而轉(zhuǎn)入到基因型依賴型的研究當中去,從單一基因研究深入到全基因組關(guān)聯(lián)分析研究中。大力推動我國基因組學(xué)在基因克隆與分子育種領(lǐng)域的應(yīng)用研究,提高園藝作物育種能力和水平。加強與轉(zhuǎn)錄組學(xué)、代謝組學(xué)、蛋白質(zhì)組學(xué)和降解組學(xué)的相關(guān)研究,促進基因組學(xué)生物信息的共享與利用。加強生物信息學(xué)教育的投入,并將其應(yīng)用于實踐中。隨著測序技術(shù)的不斷發(fā)展,相信在不久的將來園藝作物全基因組測序?qū)⑦M入快速發(fā)展的階段,為世界園藝產(chǎn)業(yè)帶來巨大貢獻。

        [HS2][HT8.5H]參考文獻:[HT8.SS]

        [1][ZK(#]中華人民共和國農(nóng)業(yè)部.中國農(nóng)業(yè)統(tǒng)計資料:2013[M]. 北京:中華人民共和國農(nóng)業(yè)部,2014.

        [2]Collins F S,Morgan M,Patrinos A. The human genome project:lessons from large-scale biology[J]. Science,2003,300(5617):286-290.

        [3]Arabidopsis Genome Initiative. Analysis of the genome sequence of the flowering plant Arabidopsis thaliana[J]. Nature,2000,408(6814):796-815.

        [4]Schloss J A. How to get genomes at one ten-thousandth the cost[J]. Nature Biotechnology,2008,26(10):1113-1115.

        [5]Margulies M,Egholm M,Altman W E,et al. Genome sequencing in microfabricated high-density picolitre reactors[J]. Nature,2005,437(757):376-380.

        [6]Velasco R,Zharkikh A,Affourtit J,et al. The genome of the domesticated apple (Malus × domestica Borkh.)[J]. Nature Genetics,2010,42(10):833-839.

        [7]Argout X,Salse J,Aury J M,et al. The genome of Theobroma cacao[J]. Nature Genetics,2011,43(2):101-108.

        [8]Rodríguez-Moreno L,González V M,Benjak A,et al. Determination of the melon chloroplast and mitochondrial genome sequences reveals that the largest reported mitochondrial genome in plants contains a significant amount of DNA having a nuclear origin[J]. BMC Genomics,2011,12:424.

        [9]Huang S,Li R,Zhang Z,et al. The genome of the cucumber,Cucumis sativus L.[J]. Nature Genetics,2009,41(12):1275-1281.

        [10][ZK(#]Shulaev V,Sargent D J,Crowhurst R N,et al. The genome of woodland strawberry (Fragaria vesca)[J]. Nature Genetics,2011,43(2):109-116.

        [11]Wang X,Wang H,Wang J,et al. The genome of the mesopolyploid crop species Brassica rapa[J]. Nature Genetics,2011,43(10):1035-1039.

        [12]Potato Genome Sequencing Consortium,Xu X,Pan S,et al. Genome sequence and analysis of the tuber crop potato[J]. Nature,2011,475(7355):189-195.

        [13]Dhont A,Denoeud F,Aury J M,et al. The banana (Musa acuminata) genome and the evolution of monocotyledonous plants[J]. Nature,2012,488(7410):213-217.

        [14]Xu Q,Chen L L,Ruan X,et al. The draft genome of sweet orange (Citrus sinensis)[J]. Nature Genetics,2013,45(1):59-66.

        [15]Guo S,Zhang J,Sun H,et al. The draft genome of watermelon (Citrullus lanatus) and resequencing of 20 diverse accessions[J]. Nature Genetics,2013,45(1):51-58.

        [16]Jaillon O,Aury J M,Noel B,et al. The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla[J]. Nature,2007,449(7161):463-467.

        [17]Ming R,Hou S,F(xiàn)eng Y,et al. The draft genome of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus)[J]. Nature,2008,452(7190):991-996.

        [18]Al-Dous E K,George B,Al-Mahmoud M E,et al. De novo genome sequencing and comparative genomics of date palm (Phoenix dactylifera)[J]. Nature Biotechnology,2011,29(6):521-527.

        [19]Sato S,Tabata S,Hirakawa H,et al. The tomato genome sequence provides insights into fleshy fruit evolution[J]. Nature,2012,485(740):635-641.

        [20][JP3]Wu J,Wang Z,Shi Z,et al. The genome of the pear (Pyrus bretschneideri Rehd.)[J]. Genome Research,2013,23(2):396-408.

        [21]Zhang Q,Chen W,Sun L,et al. The genome of Prunus mume[J]. Nature Communications,2012,3:1318.

        [22]International Peach Genome Initiative,Verde I,Abbott A G,et al. The high-quality draft genome of peach (Prunus persica) identifies unique patterns of genetic diversity,domestication and genome evolution[J]. Nature Genetics,2013,45(5):487-494.

        [23]Ming R,Vanburen R,Liu Y,et al. Genome of the long-living sacred lotus (Nelumbo nucifera Gaertn.)[J]. Genome Biology,2013,14(5):R41.

        [24]Cheng S,van den Bergh E,Zeng P,et al. The tarenaya hassleriana genome provides insight into reproductive trait and genome evolution of crucifers[J]. The Plant Cell,2013,25(8):2813-2830.

        [25]Huang S,Ding J,Deng D,et al. Draft genome of the kiwifruit Actinidia chinensis[J]. Nature Communications,2013,4:2640.

        [26]Yagi M,Kosugi S,Hirakawa H,et al. Sequence analysis of the genome of carnation (Dianthus caryophyllus L.)[J]. DNA Research:an International Journal for Rapid Publication of Reports on Genes and Genomes,2014,21(3):231-241.

        [27]Qin C,Yu C,Shen Y,et al. Whole-genome sequencing of cultivated and wild peppers provides insights into Capsicum domestication and specialization[J]. Proceedings of the National Academy of Sciences of the United States of America,2014,111(14):5135-5140.

        [28]Kitashiba H,Li F,Hirakawa H,et al. Draft sequences of the radish (Raphanus sativus L.) genome[J]. DNA Research:an International Journal for Rapid Publication of Reports on Genes and Genomes,2014,21(5):481-490.

        猜你喜歡
        園藝作物基因組學(xué)
        基于基因組學(xué)數(shù)據(jù)分析構(gòu)建腎上腺皮質(zhì)癌預(yù)后模型
        云南省農(nóng)業(yè)科學(xué)院園藝作物研究所
        辣椒雜志(2021年4期)2021-04-14 08:28:18
        新疆和西藏少數(shù)民族的群體基因組學(xué)研究
        科學(xué)(2020年2期)2020-08-24 07:56:48
        系統(tǒng)基因組學(xué)解碼反芻動物的演化
        科學(xué)(2020年2期)2020-08-24 07:56:44
        園藝作物栽培中相關(guān)技術(shù)問題
        花卉(2020年20期)2020-01-09 15:19:03
        園藝作物栽培中相關(guān)的技術(shù)問題
        植物生長調(diào)節(jié)劑在園藝作物上的應(yīng)用
        華法林出血并發(fā)癥相關(guān)藥物基因組學(xué)研究進展
        云斑白條天牛對園藝作物的危害及習(xí)性
        園藝作物藥害問題探討
        一区二区丝袜美腿视频| 成片免费观看视频大全| 色妺妺视频网| 成年视频网站在线观看777| 亚洲精品456在线播放狼人 | 波多野结衣不打码视频| 性动态图av无码专区| 人妻av一区二区三区av免费| 亚洲国产综合精品一区最新| 熟女体下毛荫荫黑森林| 人妻av中文字幕无码专区| 国产一区二区精品久久凹凸| av男人的天堂第三区| 久久99精品久久久久婷婷| 午夜亚洲www湿好爽| 亚洲国产精品久久久性色av| 亚洲精品国产av成拍色拍| 大地资源高清在线视频播放| 国产精品国产三级国产av′| 天堂Av无码Av一区二区三区| 人妻精品一区二区三区蜜桃| av免费不卡国产观看| 这里有精品可以观看| 加勒比东京热综合久久| 精品亚洲一区二区三区四区五| 亚洲熟妇av日韩熟妇在线| 国产成人拍精品免费视频| 亚洲欧美日韩在线一区| 国产成人一区二区三区视频免费蜜 | 亚洲avav天堂av在线网爱情| 亚洲综合国产精品一区二区99 | 国产精品9999久久久久仙踪林| 亚洲一区爱区精品无码| 日本在线中文字幕一区| 国产一区二区三区av天堂| 国产亚洲精品久久久久婷婷瑜伽| 精精国产xxxx视频在线播放器| 青青草免费手机直播视频| 欧美狠狠入鲁的视频777色| 香蕉成人啪国产精品视频综合网 | 国模91九色精品二三四|