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

    ‘六月早’蜜柚葉綠體基因組及其特征分析

    2021-07-20 14:00:52徐世榮陳燕瓊潘東明潘鶴立
    熱帶作物學報 2021年5期

    徐世榮 陳燕瓊 潘東明 潘鶴立

    摘 ?要:‘六月早蜜柚(Citrus maxima ‘Liuyuezao)為國內(nèi)重要的柚品種琯溪蜜柚早熟芽變品種。以‘六月早蜜柚為材料,利用二代測序進行全基因組重測序,從中組裝獲取葉綠體基因組并對其進行注釋。結(jié)果表明:組裝獲得的‘六月早蜜柚的葉綠體基因組全長160 186 bp,四分體結(jié)構(gòu)由大單拷貝區(qū)(large single copy, LSC)、小單拷貝區(qū)(small single copy, SSC)和反向重復區(qū)(inverted repeat, IR)組成,3個分區(qū)的長度分別為87 939、18 395、26 926 bp。注釋得到133個基因,其中包含89個編碼基因,37個tRNA和8個rRNA。共識別到31個短串聯(lián)重復序列,101個SSR位點。將已公開發(fā)表的29個蕓香科葉綠體基因組使用最大似然法進行系統(tǒng)發(fā)育關(guān)系的構(gòu)建,結(jié)果表明,‘六月早蜜柚與甜橙(C. sinensis)、檸檬(C. limon)和C. platymamma的親緣關(guān)系較近。

    關(guān)鍵詞:琯溪蜜柚;葉綠體基因組;SSR;系統(tǒng)發(fā)育

    中圖分類號:S961.6 ? ? ?文獻標識碼:A

    Abstract: Citrus maxima ‘Liuyuezao is a precocious bud mutation of Guanximiyou, an important pomelo variety of China. In this study, the whole genome resequencing of C. maxima ‘Liuyuezao was conducted using next-generation sequencing. The chloroplast (cp) genome was assembled and annotated. The cp genome had a total length of 160 186 bp, the tetrad structed was consisted with large single copy (LSC), small single copy (SSC) and inverted repeat (IR), with the length in 87 939 bp, 18 395 bp and 26 926 bp respectively. The cp genome contained a total of 133 genes, including 89 protein-coding genes, 37 tRNA genes and 8 tRNA genes. We identified 31 tandem repeat and 101 SSRs. Combined with 29 published cp genome sequences of Rutaceae, we used the maximum likelihood method to reconstruct the phylogenetic relationship of Rutaceae. The result showed that C. maxima ‘Liuyuezao had a close relationship with C. sinensis, C. limon and C. platymamma.

    Keywords: Guanximiyou; chloroplast genome; SSR; phylogeny

    DOI: 10.3969/j.issn.1000-2561.2021.05.004

    葉綠體在維持地球生命中起著至關(guān)重要的作用,通過光合作用和氧氣再釋放將太陽能轉(zhuǎn)化為碳水化合物。葉綠體基因組編碼許多參與光合作用和其他代謝過程的關(guān)鍵蛋白質(zhì)[1]。在被子植物中,葉綠體為母系遺傳,且具有結(jié)構(gòu)穩(wěn)定、序列高度保守和遺傳重組率低等特點,被廣泛應用于重建系統(tǒng)發(fā)育關(guān)系,揭示物種內(nèi)部和物種之間的巨大差異[2–4]。對葉綠體基因組數(shù)據(jù)的了解可提高我們對植物生物學和多樣性的認識。

    柚[Citrus maxima (Burm.) Merr.]是我國特產(chǎn)果樹,為蕓香科(Rutaceae)柑橘屬(Citrus)常綠喬木果樹,而柑橘屬水果是世界上種植最廣泛的水果之一,與蘋果、葡萄、香蕉并稱為世界四大水果。柚作為一個古老的品種,在歷史的繁衍進化過程中通過自身芽變或者與其它柑橘雜交產(chǎn)生了許多果實大小和酸度不同的新品種,如葡萄柚、酸橙等,極大的豐富了柑橘家族[5]。國內(nèi)的柚類品種主要為琯溪蜜柚、沙田柚、四季柚、度尾文旦、強德勒柚、坪山柚和暹羅柚等。其中栽培面積最大的品種為福建平和的琯溪蜜柚,僅福建地區(qū)產(chǎn)量就約占全國柚類總產(chǎn)量的2/3、世界的1/3。

    琯溪蜜柚因其果大皮薄、肉嫩汁多無核、營養(yǎng)物質(zhì)豐富及具有獨特的香氣而廣受消費者的喜愛,目前產(chǎn)品不僅占據(jù)了國內(nèi)的大小商場,還遠銷至歐盟、美國、加拿大等國家。近10多年已逐漸發(fā)展成為柑橘類果樹中最具市場競爭力和發(fā)展?jié)摿Φ乃?。但其成熟期在每年?0—12月,上市的時間過于集中,導致貨架長時間空缺,因此培育不同成熟的琯溪蜜柚是市場的需求。福建平和果農(nóng)于2007年發(fā)現(xiàn)了1株琯溪蜜柚發(fā)生了早熟芽變,經(jīng)過對母樹及突變枝嫁接子代的連續(xù)多年觀察,發(fā)現(xiàn)其性狀穩(wěn)定,成熟期比琯溪蜜柚早約45 d,由于在農(nóng)歷6月就能成熟,故命名為‘六月早。

    本研究利用高通量測序技術(shù)對‘六月早蜜柚進行測序,以已發(fā)表的柚的葉綠體基因組為參考,獲取‘六月早蜜柚的葉綠體基因組序列,進行注釋和序列特征分析,并與其他的蕓香科的物種進行比較分析和系統(tǒng)發(fā)育關(guān)系推斷,以期為蕓香科的分子標記開發(fā)和系統(tǒng)發(fā)育研究提供參考。

    1 ?材料與方法

    1.1 ?材料

    以‘六月早琯溪蜜柚為材料,其幼嫩的葉片組織采于福建平和縣小溪鎮(zhèn)果園。

    1.2 ?方法

    1.2.1 ?總DNA的提取和測序 ?利用TIANGEN試劑盒提取‘六月早琯溪蜜柚的總DNA,用瓊脂糖凝膠電泳確定總DNA具有較好的完整性,并用分光光度計檢測DNA的質(zhì)量。檢測的DNA濃度>50 ng/μL,無色素污染、無明顯RNA、無蛋白質(zhì)等雜質(zhì)污染的樣品視為合格樣品。將檢測合格的DNA樣品送至深圳華大基因科技服務有限公司進行文庫的構(gòu)建,最后在Illumina HiSeq平臺上機測序,獲得質(zhì)控過濾后的Clean data 12 Gb的數(shù)據(jù)量,覆蓋該物種測序深度37×。

    1.2.2 ?葉綠體基因組的組裝與注釋 ?對下機的數(shù)據(jù)進行過濾,包括去除低質(zhì)量的序列和接頭。從GenBank中下載已發(fā)表的柚葉綠體基因組序列(GenBank登錄號:KY055833.1)作為參考序列,使用GetOrganelle toolkit[6]初步組裝出葉綠體基因組的序列,再使用Bandage[7]根據(jù)葉綠體基因組的機構(gòu)特征手動拼接出完整的質(zhì)體基因組。

    使用Geneirous primer v2019.03[8]對組裝好的序列進行注釋,注釋使用已發(fā)表的柚作為注釋參考,再人工校正基因的起始和終止密碼的位置。序列的注釋可視化在OrganellarGenomeDRAW(OGDRAW)v1.3.1(https://chlorobox.mpimp-golm. mpg.de/OGDraw.html)進行[9]。將組裝注釋好的‘六月早蜜柚的葉綠體基因組上傳至GenBank數(shù)據(jù)庫,獲得登錄號:MT527726。

    1.2.3 ?葉綠體基因組的特征分析 ?串聯(lián)重復序列的識別用Tandem repeat finder軟件[10],使用默認參數(shù)。簡單重復序列(simple sequence repeat,SSR)的識別使用Perl腳本MISA(microSAtellite identification tool),參數(shù)設置為單核苷酸(mononucleotide)、二核苷酸(dinucleotide)、三核苷酸(trinucleotide)、四核苷酸(tetranucleotide)、五核苷酸(pentanucleotide)和六核苷酸(hexanucleotide)的重復數(shù)閾值分別為10、5、4、3、3和3。使用在線分析軟件CodonW 1.4.4(https:// galaxy.pasteur.fr/?form=codonw)對葉綠體基因組內(nèi)的編碼基因進行密碼子的使用偏好性分析。

    1.2.4 ?蕓香科物種葉綠體基因組比較分析 ?從GenBank數(shù)據(jù)庫中下載檸檬(Citrus limon)、山小橘(Glycosmis pentaphylla)、假黃皮(Clausena excavata)、蕓香(Ruta graveolens)、花椒(Zanthoxylum bungeanum)、黃檗(Phellodendron amurense)、香肉果(Casimiroa edulis)等29個蕓香科隸屬于11個屬的其他物種的葉綠體基因組。每個屬選1個代表種共11個種進行邊界的擴張與收縮分析,使用在線程序Irscope(https:// irscope.shinyapps.io/irapp/)對IR、SSC和LSC區(qū)域進行邊界變化分析[11]。

    以無患子科(Sapindaceae)欒樹(Koelreuteria paniculata)為外類群,與其他29個蕓香科物種進行系統(tǒng)發(fā)育樹的構(gòu)建。將30個葉綠體基因組序列使用Homblock軟件進行同源共線性的比對[12],識別出的共線性區(qū)域用于系統(tǒng)發(fā)育樹的構(gòu)建。使用IQ-tree對數(shù)據(jù)集進行系統(tǒng)發(fā)育樹的構(gòu)建[13]。

    2 ?結(jié)果與分析

    2.1 ?‘六月早蜜柚葉綠體基因組的基本特征

    ‘六月早蜜柚的葉綠體基因組結(jié)構(gòu)為經(jīng)典四段式結(jié)構(gòu),包含大單拷貝區(qū)(large single copy,LSC),2個反向重復區(qū)(inverted repeat,IR)和小單拷貝區(qū)(small single copy,SSC)(圖1)。葉綠體基因組的總長度為160 186 bp,GC含量為38.5%,其中LSC區(qū)長87 939 bp,IR區(qū)長26 926 bp,SSC區(qū)長18 395 bp。共有基因133個,包括編碼基因89個,tRNA 37個,rRNA8個。其中20個基因含有2個拷貝,包括8個蛋白編碼基因(rpl2、rpl23、rps7、rps19、ndhB、ycf1、ycf2和ycf15),8個tRNA(trnA-UGC、trnG-GCC、trnI-CAU、trnI-GAU、trnL-CAA、trnN-GUU、trnR-ACG和trnV-GAC)和4個rRNA(rrn4.5、rrn5、rrn16和rrn23)。

    2.2 ?‘六月早蜜柚葉綠體基因組的重復序列

    識別葉綠體基因組的重復序列,特別是SSR對于分子標記的開發(fā),有助于為物種的遺傳發(fā)育及分子育種研究提供參考。對‘六月早的重復序列進行檢測,共識別出31個短串聯(lián)重復序列的區(qū)域,絕大部分分布在基因間隔區(qū)(IGS),長度范圍在25~75 bp之間,占基因組序列的0.81%(表1)。共檢測到101個SSR,其中包括5個復合型SSR、68個單堿基重復、3個二堿基重復、13個三堿基重復、11個四堿基重復和1個五堿基重復(表2,圖2)。SSR序列長度為10~35 bp,大部分SSR均分布在基因組的大單拷貝區(qū)(占74.26%)和基因間隔區(qū)(占76.34%,圖2)。

    2.3 ?‘六月早蜜柚葉綠體基因組的密碼子偏好性

    密碼子是基因編碼區(qū)蛋白翻譯的核心元件,分析編碼區(qū)密碼子的使用特征,對基因功能和系統(tǒng)進化研究具有重要意義,為后續(xù)研究提供基礎資料。將‘六月早蜜柚葉綠體基因組中得到的88個編碼基因(coding DNA sequence,CDS)用于分析,為了保證結(jié)果的準確性,去除了重復的基因序列和長度小于300 bp序列,將剩下的52條CDS用于密碼子的偏好性分析。RSCU值大于1.00的密碼子為30個,其中29個以A或U結(jié)尾,1個以G或C結(jié)尾(表3),表明這些密碼子偏向A/U結(jié)尾,‘六月早蜜柚與大多數(shù)被子植物的密碼子偏好性類似。

    2.4 ?IR邊界的收縮與擴張

    選擇了11個物種代表蕓香科的11個屬進行IR邊界的比較分析。從圖3可見,蕓香科植物具有相似的基因結(jié)構(gòu)和基因順序,差異主要發(fā)生在JLB(LSC與IRb邊界)和JSB(IRb與SSC邊界)邊界。在JLB邊界表現(xiàn)出2種類型,第1種類型為無基因橫跨邊界,僅見黃皮屬的假黃皮Clausena excavata和小蕓木屬的Micromelum minutum,位于LSC的rps3為離邊界最近的基因分別距JLB 38 bp和115 bp;第2種類型是rpl22基因橫跨了JLB邊界,其余9個種均為這種類型,其中跨到LSC邊界的長度為67~295 bp。JSB邊界則表現(xiàn)出了3種類型,第1種類型為無基因處于邊界上,僅表現(xiàn)在酒餅簕屬廣東酒餅簕(Atalantia kwangtungensis)中;第2種類型為ycf1基因橫跨JSB邊界,在黃檗屬黃檗(Phellodendron amurense)、花椒屬(Zanthoxylum tragodes)、蕓香屬蕓香(Ruta graveolens)和‘六月早蜜柚4個種中表現(xiàn)有1~35 bp跨到SSC區(qū);第3種類型為ndhF基因橫跨JSB邊界,其余的6個種、蕓香和‘六月早蜜柚有1~31 bp跨越IRa區(qū)域。JLA邊界也表現(xiàn)出3種類型,表現(xiàn)為不同的種在邊界處trnH和rpl22基因的位置不同,第1種為邊界無基因;第2種為trnH跨越邊界;第3種為rpl22跨越邊界。JSA則非常保守,均為ycf1基因跨越了邊界,有1080~1130 bp位于IRa。

    2.5 ?系統(tǒng)發(fā)育分析

    為了確定‘六月早蜜柚在蕓香科中的進化地位,基于葉綠體基因組序列構(gòu)建29個的蕓香科物種系統(tǒng)發(fā)育關(guān)系,其中柑橘屬10種,酒餅簕屬1種,Merrillia屬1種,山小橘屬2種,小蕓木屬1種,咖喱樹屬1種,黃皮屬1種,蕓香屬1種,花椒屬9種,黃檗屬1種和香肉果屬1種,以無患子科欒樹為外類群。由圖4可知,蕓香科由兩大進化支組成,各分支均獲得很高的自展值(BSML均為100);花椒屬、黃檗屬和香肉果屬聚為一支,為最先分化的分支,且均屬于Amyridoideae亞科;柑橘屬、酒餅簕屬、Merrillia屬、山小橘屬、小蕓木屬、咖喱樹屬、黃皮屬和蕓香屬聚為一支,其中蕓香屬為最先分化的屬,其他的7個屬聚為一支;‘六月早蜜柚與其他柑橘屬聚為一支,與Citrus platymamma、檸檬和甜橙分支成姐妹關(guān)系。

    3 ?討論

    ‘六月早蜜柚同其他被子植物一樣,具有四分結(jié)構(gòu)環(huán)狀的葉綠體基因組,其基因順序、GC含量也同其他柑橘屬的植物基本一樣[14–16],這也說明了在柑橘屬內(nèi)葉綠體基因組非常保守。植物的葉綠體基因組長度在100~220 kb之間,蕓香科內(nèi)葉綠體基因組長度在157 434 bp(蕓香)~161 172 bp(假黃皮)之間,而差異主要表現(xiàn)在反向重復(IR)區(qū)的收縮、擴張或缺失,因此,IR區(qū)對于葉綠體基因組的結(jié)構(gòu)穩(wěn)定性和決定葉綠體基因組大小方面的作用不可忽視[17-18]。重復序列的存在從原核生物到真核生物的進化過程中呈遞增趨勢,具有保護編碼基因和作為進化的動力促進形成新基因的作用,是進化的必然結(jié)果[19]。‘六月早蜜柚的重復序列大部分分布于IGS區(qū),且超過一半(16/31)分布于LSC區(qū),這表明重復序列對該物種不具有特定功能。SSR廣泛存在于基因組中,長度在1~6 bp短的重復基序(motif),具有高多態(tài)性、多等位性、共顯性等特點,同時因其具有良好的種間甚至種內(nèi)的遺傳變異區(qū)分精度,一直被廣泛用于植物資源研究[20–23],本研究從‘六月早蜜柚葉綠體基因組中識別到的SSR可提供該品種的種間多態(tài)性,為柑橘屬品種間的親緣關(guān)系和分子鑒定提供可靠的遺傳標記。

    本研究將已公開發(fā)表的所有蕓香科植物的葉綠體基因組從GenBank中獲取,包含29個種,隸屬于11個屬,3個亞科,對其進行蕓香科的系統(tǒng)發(fā)育關(guān)系的重建。在前人的研究中,僅使用葉綠體基因片段進行蕓香科的構(gòu)建,如使用rps16內(nèi)含子和trnL-trnF基因間隔區(qū)對65個蕓香科物種進行分析,發(fā)現(xiàn)亞科與族的關(guān)系都需要重新評估[24]。本研究結(jié)果獲得了穩(wěn)固的系統(tǒng)發(fā)育關(guān)系,各亞科和各屬的單系性均得到了很好的支持,其中柑橘亞科與蕓香亞科成姐妹關(guān)系,構(gòu)成的分支與最早分化的脂檀亞科成姐妹關(guān)系。但由于本研究使用的物種序列有限,構(gòu)建的系統(tǒng)發(fā)育關(guān)系仍不夠全面,有待進一步完善。

    參考文獻

    [1] Daniell H, Lin C S, Yu M, et al. Chloroplast genomes: diversity, evolution, and applications in genetic engineering[J]. Genome Biology, 2016, 17(1): 134.

    [2] Jansen R K, Cai Z, Raubeson L A, et al. Analysis of 81 genes from 64 plastid genomes resolves relationships in angiosperms and identifies genome-scale evolutionary patterns[J]. Proceedings of the National Academy of Sciences, 2007, 104(49): 19 369-19 374.

    [3] Song Y, Yu W B, Tan Y, et al. Evolutionary comparisons of the chloroplast genome in Lauraceae and insights into loss events in the Magnoliids[J]. Genome Biology and Evo-lution, 2017, 9(9): 2354-2364.

    [4] Gao C, Deng Y, Wang J. The complete chloroplast genomes of Echinacanthus species (Acanthaceae): phylogenetic relationships, adaptive evolution, and screening of mo-lecular markers[J]. Frontiers in Plant Science, 2018, 9: 1989.

    [5] Wu G A, Terol J, Ibanez V, et al. Genomics of the origin and evolution of Citrus[J]. Nature, 2018, 554(7692): 311-316.

    [6] Jin J J, Yu W B, Yang J B, et al. GetOrganelle: a fast and versatile toolkit for accurate de novo assembly of organelle genomes[J]. Genome Biology, 2020, 21(1): 1-31.

    [7] Wick R R, Schultz M B, Zobel J, et al. Bandage: interactive visualization of de novo genome assemblies[J]. Bioinformatics, 2015, 31(20): 3350-3352.

    [8] Kearse M, Moir R, Wilson A, et al. Geneious basic: an integrated and extendable desktop software platform for the organization and analysis of sequence data[J]. Bioinformatics, 2012, 28(12): 1647-1649.

    [9] Greiner S, Lehwark P, Bock R. OrganellarGenomeDRAW (OGDRAW) version 1.3.1: expanded toolkit for the graphical visualization of organellar genomes[J]. Nucleic Acids Research, 2019, 47(W1): 59-64.

    [10] Benson G. Tandem repeats finder: a program to analyze DNA sequences[J]. Nucleic Acids Research, 1999, 27(2): 573-580.

    [11] Amiryousefi A, Hyv?nen J, Poczai P. Irscope: an online program to visualize the junction sites of chloroplast genomes[J]. Bioinformatics, 2018, 34(17): 3030-3031.

    [12] Bi G, Mao Y, Xing Q, et al. Homblocks: a multiple-alignment construction pipeline for organelle phylogenomics based on locally collinear block searching[J]. Genomics, 2018, 110(1): 18-22.

    [13] Nguyen L T, Schmidt H A, Von Haeseler A, et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies[J]. Molecular Biology and Evolution, 2015, 32(1): 268-274.

    [14] Bausher M G, Singh N D, Lee S B, et al. The complete chloroplast genome sequence of Citrus sinensis (L.) Osbeck var ‘Ridge Pineapple: organization and phylogenetic relationships to other angiosperms[J]. BMC Plant Biology, 2006, 6(1): 725-730.

    [15] Su H J, Hogenhout S A, Al-Sadi A M, et al. Complete chloroplast genome sequence of Omani lime (Citrus aurantiifolia) and comparative analysis within the rosids[J]. PLoS One, 2014, 9(11): 113049.

    [16] Xu S R, Huang C Y, Deng Y T, et al. The complete chloroplast genome sequence of Citrus maxima (Burm.) Merr. ‘Guanximiyou[J]. Mitochondrial DNA Part B, 2020, 5(1): 482-483.

    [17] Palmer J D, Thompson W F. Chloroplast DNA rearrangements are more frequent when a large inverted repeat sequence is lost[J]. Cell, 1982, 29(2): 537-550.

    [18] Wang R J, Cheng C L, Chang C C, et al. Dynamics and evolution of the inverted repeat-large single copy junctions in the chloroplast genomes of monocots[J]. BMC Evolutionary Biology, 2008, 8(1): 1-14.

    [19] 艾對元. 基因組中重復序列的意義[J]. 生命的化學, 2008(3): 343-345.

    [20] Provan J, Powell W, Hollingsworth P M. Chloroplast microsatellites: new tools for studies in plant ecology and evolution[J]. Trends in Ecology & Evolution, 2001, 16(3): 142-147.

    [21] Kantety R V, La Rota M, Matthews D E, et al. Data mining for simple sequence repeats in expressed sequence tags from barley, maize, rice, sorghum and wheat[J]. Plant Molecular Biology, 2002, 48(5-6): 501-510.

    [22] Lassois L, Denancé C, Ravon E, et al. Genetic diversity, population structure, parentage analysis, and construction of core collections in the French apple germplasm based on SSR markers[J]. Plant Molecular Biology Reporter, 2016, 34(4): 827-844.

    [23] Singh G, Aulakh P S, Sarao N K, et al. Molecular verification of putative zygotic seedlings in different intra-specific crosses in Mandarins (Citrus reticulata) by SSR markers[J]. Agricultural Eesearch, 2019, 8(1): 21-26.

    [24] Groppo M, Pirani J R, Salatino M L F, et al. Phylogeny of Rutaceae based on twononcoding regions from cpDNA[J]. American Journal of Botany, 2008, 95(8): 985-1005.

    責任編輯:黃東杰

    一级毛片女人18水好多 | www.999成人在线观看| 日本欧美国产在线视频| 日本猛色少妇xxxxx猛交久久| 狠狠婷婷综合久久久久久88av| 免费黄频网站在线观看国产| 超碰成人久久| 久久久久国产一级毛片高清牌| 永久免费av网站大全| 亚洲精品美女久久av网站| 丁香六月天网| 咕卡用的链子| 欧美大码av| 亚洲欧洲国产日韩| 亚洲七黄色美女视频| 久久久久久免费高清国产稀缺| 国产男女内射视频| 一区二区三区精品91| 高清视频免费观看一区二区| 久久久国产欧美日韩av| 蜜桃国产av成人99| 欧美乱码精品一区二区三区| 国产日韩一区二区三区精品不卡| 高清av免费在线| 99精品久久久久人妻精品| 这个男人来自地球电影免费观看| 亚洲七黄色美女视频| 看十八女毛片水多多多| 国产又色又爽无遮挡免| 国产欧美日韩一区二区三区在线| 亚洲欧美日韩另类电影网站| 我要看黄色一级片免费的| 久久99热这里只频精品6学生| 国产人伦9x9x在线观看| 国产一区亚洲一区在线观看| 欧美成狂野欧美在线观看| 人成视频在线观看免费观看| 99国产精品免费福利视频| 日本欧美视频一区| 国产野战对白在线观看| 精品国产超薄肉色丝袜足j| 日本欧美国产在线视频| av国产久精品久网站免费入址| 国产片内射在线| 国产成人啪精品午夜网站| av又黄又爽大尺度在线免费看| av在线播放精品| 中文字幕精品免费在线观看视频| 下体分泌物呈黄色| 性色av一级| 亚洲图色成人| 91麻豆精品激情在线观看国产 | 久久人人97超碰香蕉20202| 97在线人人人人妻| av线在线观看网站| 999久久久国产精品视频| av在线app专区| 中文字幕亚洲精品专区| 国产欧美日韩一区二区三区在线| 美女中出高潮动态图| 久久精品熟女亚洲av麻豆精品| 久久精品aⅴ一区二区三区四区| 日韩精品免费视频一区二区三区| 精品国产一区二区三区四区第35| 在线观看免费视频网站a站| 成人午夜精彩视频在线观看| 午夜福利,免费看| 看十八女毛片水多多多| 人人妻人人爽人人添夜夜欢视频| 成年人免费黄色播放视频| 久久女婷五月综合色啪小说| 亚洲av电影在线进入| 亚洲久久久国产精品| 欧美日本中文国产一区发布| 久久精品人人爽人人爽视色| 欧美精品一区二区免费开放| 国产极品粉嫩免费观看在线| 99久久精品国产亚洲精品| 日韩大片免费观看网站| 欧美另类一区| 精品高清国产在线一区| 日本欧美国产在线视频| 黄网站色视频无遮挡免费观看| 欧美+亚洲+日韩+国产| h视频一区二区三区| 成年人免费黄色播放视频| 满18在线观看网站| 制服诱惑二区| 久久久久视频综合| 久久热在线av| 人妻人人澡人人爽人人| 亚洲精品成人av观看孕妇| 亚洲伊人色综图| 两个人免费观看高清视频| 日本黄色日本黄色录像| 亚洲一码二码三码区别大吗| 欧美日韩视频精品一区| 少妇的丰满在线观看| 成人18禁高潮啪啪吃奶动态图| 亚洲中文日韩欧美视频| 亚洲精品第二区| 亚洲专区中文字幕在线| 精品亚洲成国产av| av视频免费观看在线观看| 看十八女毛片水多多多| 久久免费观看电影| 波多野结衣一区麻豆| 一区二区av电影网| 精品欧美一区二区三区在线| 搡老乐熟女国产| 99九九在线精品视频| 老司机影院成人| 国产一级毛片在线| 新久久久久国产一级毛片| 美女中出高潮动态图| 深夜精品福利| 波野结衣二区三区在线| 国产成人精品久久二区二区91| 两性夫妻黄色片| av网站在线播放免费| 免费一级毛片在线播放高清视频 | 悠悠久久av| 国产爽快片一区二区三区| 国产又爽黄色视频| 久久毛片免费看一区二区三区| 精品少妇黑人巨大在线播放| 久久精品人人爽人人爽视色| 色网站视频免费| 日韩人妻精品一区2区三区| 国产极品粉嫩免费观看在线| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜爱| 欧美日韩精品网址| 亚洲伊人色综图| 黄片小视频在线播放| 精品免费久久久久久久清纯 | 咕卡用的链子| 国产一区二区在线观看av| 国产成人欧美在线观看 | 久久中文字幕一级| 亚洲一码二码三码区别大吗| 母亲3免费完整高清在线观看| 男女之事视频高清在线观看 | 国产熟女午夜一区二区三区| e午夜精品久久久久久久| 校园人妻丝袜中文字幕| 国产av国产精品国产| 看免费av毛片| 咕卡用的链子| 亚洲精品一卡2卡三卡4卡5卡 | 伊人久久大香线蕉亚洲五| 青春草视频在线免费观看| 深夜精品福利| 狠狠婷婷综合久久久久久88av| 国产色视频综合| 久久性视频一级片| 99国产综合亚洲精品| 亚洲激情五月婷婷啪啪| 欧美精品一区二区免费开放| 亚洲五月婷婷丁香| 久久九九热精品免费| 久久久久久亚洲精品国产蜜桃av| 久久人人爽av亚洲精品天堂| 亚洲欧美日韩高清在线视频 | 手机成人av网站| 精品亚洲成a人片在线观看| 亚洲第一av免费看| 成在线人永久免费视频| 久久久久久久久久久久大奶| 视频区图区小说| 日本wwww免费看| 亚洲精品久久午夜乱码| 午夜日韩欧美国产| 咕卡用的链子| 免费在线观看影片大全网站 | 亚洲国产精品国产精品| 亚洲精品美女久久av网站| 女性生殖器流出的白浆| 亚洲欧美精品自产自拍| a级毛片在线看网站| 日本色播在线视频| 久久精品久久久久久噜噜老黄| 成人国产av品久久久| 少妇人妻久久综合中文| 欧美亚洲 丝袜 人妻 在线| 大香蕉久久成人网| 亚洲av男天堂| xxx大片免费视频| 久久久久国产一级毛片高清牌| 日韩av免费高清视频| 亚洲欧洲精品一区二区精品久久久| 国产精品免费大片| 美女福利国产在线| 婷婷成人精品国产| 久热这里只有精品99| 精品国产一区二区三区久久久樱花| 丝袜美腿诱惑在线| 中文字幕另类日韩欧美亚洲嫩草| 国产又色又爽无遮挡免| 久久久欧美国产精品| 亚洲久久久国产精品| 久久久久久久久久久久大奶| 欧美人与性动交α欧美精品济南到| 亚洲欧洲日产国产| 91老司机精品| 国产精品熟女久久久久浪| 在线 av 中文字幕| 777久久人妻少妇嫩草av网站| 老司机影院成人| 久久人人97超碰香蕉20202| 一区福利在线观看| 深夜精品福利| 看免费成人av毛片| 满18在线观看网站| 日本av免费视频播放| 国产欧美日韩综合在线一区二区| 久久久国产一区二区| 久久久久久久久免费视频了| 亚洲,一卡二卡三卡| 9热在线视频观看99| 日韩一区二区三区影片| 欧美 日韩 精品 国产| 久久久久久人人人人人| 考比视频在线观看| 女性生殖器流出的白浆| 欧美激情极品国产一区二区三区| 蜜桃在线观看..| 免费在线观看完整版高清| 精品人妻一区二区三区麻豆| 丁香六月天网| 考比视频在线观看| 久久国产精品大桥未久av| 久久免费观看电影| av福利片在线| 一级毛片 在线播放| av在线播放精品| 夜夜骑夜夜射夜夜干| 久久国产精品影院| 欧美日韩亚洲高清精品| 国产欧美日韩一区二区三 | 纵有疾风起免费观看全集完整版| cao死你这个sao货| 首页视频小说图片口味搜索 | 国产欧美日韩一区二区三区在线| tube8黄色片| 亚洲精品久久成人aⅴ小说| 午夜免费观看性视频| 亚洲,欧美精品.| 亚洲av欧美aⅴ国产| 婷婷色麻豆天堂久久| 国产麻豆69| 2018国产大陆天天弄谢| 在线 av 中文字幕| 亚洲av综合色区一区| 男人舔女人的私密视频| 高清视频免费观看一区二区| 丝袜美腿诱惑在线| 99国产精品一区二区蜜桃av | 日韩av在线免费看完整版不卡| 亚洲激情五月婷婷啪啪| 国产亚洲欧美精品永久| 9色porny在线观看| 精品一区二区三区av网在线观看 | 久久99热这里只频精品6学生| 亚洲色图综合在线观看| 久热爱精品视频在线9| 亚洲国产欧美日韩在线播放| tube8黄色片| 免费高清在线观看日韩| 纵有疾风起免费观看全集完整版| 中文字幕另类日韩欧美亚洲嫩草| 亚洲人成电影免费在线| 亚洲精品国产一区二区精华液| 久久久久久久久免费视频了| 秋霞在线观看毛片| 男人舔女人的私密视频| 国产亚洲av片在线观看秒播厂| 亚洲国产成人一精品久久久| 另类亚洲欧美激情| 男女床上黄色一级片免费看| 久久女婷五月综合色啪小说| 一级毛片黄色毛片免费观看视频| 一级片免费观看大全| 桃花免费在线播放| 国语对白做爰xxxⅹ性视频网站| 曰老女人黄片| 亚洲国产最新在线播放| 欧美性长视频在线观看| 亚洲av综合色区一区| 国产欧美日韩一区二区三区在线| 欧美黑人精品巨大| 国产亚洲欧美在线一区二区| 国产精品av久久久久免费| 欧美精品人与动牲交sv欧美| 国产熟女午夜一区二区三区| 伦理电影免费视频| 亚洲精品久久午夜乱码| 国产精品久久久久久精品电影小说| 亚洲精品一区蜜桃| 香蕉丝袜av| 别揉我奶头~嗯~啊~动态视频 | 亚洲男人天堂网一区| 日韩欧美一区视频在线观看| 岛国毛片在线播放| 亚洲精品日韩在线中文字幕| bbb黄色大片| 又粗又硬又长又爽又黄的视频| 老司机亚洲免费影院| 亚洲专区中文字幕在线| 成年人免费黄色播放视频| 男的添女的下面高潮视频| 在线观看免费日韩欧美大片| 久久九九热精品免费| 两个人看的免费小视频| 色婷婷久久久亚洲欧美| www.自偷自拍.com| 国产精品一区二区在线观看99| 日韩av不卡免费在线播放| 在线av久久热| 一区二区三区乱码不卡18| 免费在线观看影片大全网站 | 久久精品国产综合久久久| 欧美亚洲日本最大视频资源| 国产精品久久久av美女十八| 电影成人av| 亚洲成国产人片在线观看| 精品高清国产在线一区| 丝袜喷水一区| 大香蕉久久网| 亚洲精品自拍成人| 精品第一国产精品| 2021少妇久久久久久久久久久| 悠悠久久av| 叶爱在线成人免费视频播放| 后天国语完整版免费观看| 多毛熟女@视频| 91精品伊人久久大香线蕉| 亚洲国产欧美日韩在线播放| 成人午夜精彩视频在线观看| 后天国语完整版免费观看| 精品一区在线观看国产| 国产一级毛片在线| 亚洲,欧美,日韩| 久久99精品国语久久久| 日日摸夜夜添夜夜爱| 亚洲精品中文字幕在线视频| 日本黄色日本黄色录像| 制服诱惑二区| 亚洲av成人精品一二三区| 老司机亚洲免费影院| 国产日韩欧美在线精品| 久久99精品国语久久久| 国产精品亚洲av一区麻豆| 亚洲欧美中文字幕日韩二区| videos熟女内射| 母亲3免费完整高清在线观看| 成年女人毛片免费观看观看9 | av天堂在线播放| 精品免费久久久久久久清纯 | 亚洲av综合色区一区| h视频一区二区三区| 国产免费福利视频在线观看| 性少妇av在线| 我要看黄色一级片免费的| 男女午夜视频在线观看| 亚洲成国产人片在线观看| 日本欧美国产在线视频| 精品人妻1区二区| 老汉色av国产亚洲站长工具| 国产主播在线观看一区二区 | 国产日韩欧美亚洲二区| 99精国产麻豆久久婷婷| 老司机午夜十八禁免费视频| 久久久精品94久久精品| 激情视频va一区二区三区| 热99久久久久精品小说推荐| 国产在线视频一区二区| 熟女av电影| 欧美在线一区亚洲| 一级毛片 在线播放| 国产片特级美女逼逼视频| 国产精品香港三级国产av潘金莲 | 一本—道久久a久久精品蜜桃钙片| 91九色精品人成在线观看| 超碰97精品在线观看| 亚洲精品一二三| 精品国产乱码久久久久久男人| 中文精品一卡2卡3卡4更新| 99国产精品一区二区三区| 18禁观看日本| 国产成人精品久久二区二区91| 国产亚洲午夜精品一区二区久久| 亚洲国产最新在线播放| 亚洲欧美日韩另类电影网站| 精品视频人人做人人爽| 久久精品国产亚洲av高清一级| 可以免费在线观看a视频的电影网站| 久久久久久久精品精品| 婷婷色综合www| 亚洲精品一区蜜桃| 免费黄频网站在线观看国产| 18在线观看网站| 精品一区二区三区av网在线观看 | 久久99热这里只频精品6学生| 一级黄片播放器| 性色av一级| 另类亚洲欧美激情| 交换朋友夫妻互换小说| 91麻豆av在线| 91精品伊人久久大香线蕉| 欧美乱码精品一区二区三区| 飞空精品影院首页| 精品人妻在线不人妻| 国产精品.久久久| 国产成人影院久久av| 亚洲成国产人片在线观看| 久久国产亚洲av麻豆专区| 青春草视频在线免费观看| 高清av免费在线| 自拍欧美九色日韩亚洲蝌蚪91| 91老司机精品| 久久精品国产a三级三级三级| 性高湖久久久久久久久免费观看| 久久久国产一区二区| 2021少妇久久久久久久久久久| 亚洲国产毛片av蜜桃av| 五月开心婷婷网| 啦啦啦在线观看免费高清www| 色婷婷久久久亚洲欧美| 巨乳人妻的诱惑在线观看| 亚洲欧美一区二区三区黑人| 欧美日韩福利视频一区二区| 2021少妇久久久久久久久久久| av在线老鸭窝| 成年美女黄网站色视频大全免费| 一区二区三区精品91| 男女边吃奶边做爰视频| 一级黄片播放器| 亚洲精品美女久久av网站| 精品一区二区三区四区五区乱码 | 狠狠精品人妻久久久久久综合| 男女之事视频高清在线观看 | 久久99热这里只频精品6学生| 国产精品二区激情视频| 黄频高清免费视频| 成年女人毛片免费观看观看9 | 久久狼人影院| 欧美激情极品国产一区二区三区| 可以免费在线观看a视频的电影网站| 五月天丁香电影| 在线观看免费视频网站a站| 亚洲av成人精品一二三区| 丰满少妇做爰视频| 国产精品麻豆人妻色哟哟久久| 久久99精品国语久久久| 色视频在线一区二区三区| 欧美日韩黄片免| 久久久亚洲精品成人影院| 欧美精品av麻豆av| 国产精品国产三级国产专区5o| 亚洲,一卡二卡三卡| 考比视频在线观看| 大型av网站在线播放| 9191精品国产免费久久| av天堂在线播放| 大陆偷拍与自拍| 久久久久久亚洲精品国产蜜桃av| 一区二区三区乱码不卡18| 日本欧美视频一区| 只有这里有精品99| 一边摸一边抽搐一进一出视频| 亚洲第一青青草原| 精品熟女少妇八av免费久了| 久久久久网色| 亚洲人成网站在线观看播放| 亚洲七黄色美女视频| 国产高清videossex| 日本91视频免费播放| 国产女主播在线喷水免费视频网站| 欧美国产精品一级二级三级| 晚上一个人看的免费电影| 黄色一级大片看看| 成人国语在线视频| 一本大道久久a久久精品| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| av在线播放精品| 欧美黑人精品巨大| 只有这里有精品99| 亚洲自偷自拍图片 自拍| 久久精品aⅴ一区二区三区四区| 两性夫妻黄色片| 亚洲专区国产一区二区| 中文字幕人妻丝袜一区二区| 久久久精品国产亚洲av高清涩受| 99九九在线精品视频| 91精品三级在线观看| 色婷婷av一区二区三区视频| 国产精品久久久久成人av| 亚洲av日韩精品久久久久久密 | 只有这里有精品99| 午夜福利视频精品| 国产又爽黄色视频| 日韩大片免费观看网站| 日韩一区二区三区影片| 大陆偷拍与自拍| 精品卡一卡二卡四卡免费| 女人高潮潮喷娇喘18禁视频| 国产福利在线免费观看视频| 中文字幕精品免费在线观看视频| 美女中出高潮动态图| 日韩,欧美,国产一区二区三区| 美女扒开内裤让男人捅视频| 亚洲成人免费电影在线观看 | 大话2 男鬼变身卡| 真人做人爱边吃奶动态| 波多野结衣av一区二区av| 男女边摸边吃奶| 一级黄色大片毛片| 国产精品国产三级国产专区5o| 美女脱内裤让男人舔精品视频| 黄片播放在线免费| 国产真人三级小视频在线观看| av片东京热男人的天堂| 在线观看www视频免费| 精品免费久久久久久久清纯 | 国产精品熟女久久久久浪| 精品国产一区二区三区四区第35| 亚洲一码二码三码区别大吗| 亚洲自偷自拍图片 自拍| 欧美老熟妇乱子伦牲交| 国产精品秋霞免费鲁丝片| 建设人人有责人人尽责人人享有的| 久久国产精品大桥未久av| 亚洲综合色网址| 天堂俺去俺来也www色官网| 大香蕉久久成人网| 亚洲精品国产av成人精品| av在线老鸭窝| 精品欧美一区二区三区在线| 国产欧美日韩一区二区三区在线| 后天国语完整版免费观看| 亚洲一区中文字幕在线| 亚洲精品美女久久久久99蜜臀 | 80岁老熟妇乱子伦牲交| av网站在线播放免费| 1024香蕉在线观看| 欧美日韩国产mv在线观看视频| 精品久久久久久久毛片微露脸 | 精品第一国产精品| 美女脱内裤让男人舔精品视频| 国产成人av激情在线播放| 人体艺术视频欧美日本| 久久精品亚洲熟妇少妇任你| 成人黄色视频免费在线看| 999精品在线视频| av又黄又爽大尺度在线免费看| 在线观看一区二区三区激情| 精品少妇久久久久久888优播| 亚洲av美国av| 国产精品久久久久成人av| 成年美女黄网站色视频大全免费| 蜜桃国产av成人99| 亚洲人成电影免费在线| 在线观看免费视频网站a站| 男女边吃奶边做爰视频| 午夜福利免费观看在线| 黄色毛片三级朝国网站| av网站免费在线观看视频| 一边摸一边抽搐一进一出视频| 亚洲人成电影免费在线| 97在线人人人人妻| 国产一区二区 视频在线| 国产一区二区三区综合在线观看| 久久精品亚洲熟妇少妇任你| 电影成人av| 秋霞在线观看毛片| 久久久久久久久久久久大奶| 亚洲av成人不卡在线观看播放网 | 亚洲一区二区三区欧美精品| 丝瓜视频免费看黄片| 69精品国产乱码久久久| 老司机靠b影院| 叶爱在线成人免费视频播放| 精品国产超薄肉色丝袜足j| 少妇人妻久久综合中文| 亚洲国产精品国产精品| 国产片内射在线| 欧美黄色片欧美黄色片| 高清欧美精品videossex| 亚洲人成77777在线视频| 高清av免费在线| 嫩草影视91久久| e午夜精品久久久久久久| 国产又色又爽无遮挡免| xxxhd国产人妻xxx| 99国产精品一区二区蜜桃av | 后天国语完整版免费观看| 日韩 欧美 亚洲 中文字幕| 男女边吃奶边做爰视频| 熟女av电影| 十八禁人妻一区二区| 99香蕉大伊视频| 亚洲欧美一区二区三区久久| 两个人免费观看高清视频| 自拍欧美九色日韩亚洲蝌蚪91| 美女视频免费永久观看网站| 你懂的网址亚洲精品在线观看| 一本大道久久a久久精品| 性高湖久久久久久久久免费观看| 91精品伊人久久大香线蕉| 一二三四社区在线视频社区8| 两个人看的免费小视频| 久久精品熟女亚洲av麻豆精品| 亚洲精品av麻豆狂野|