• <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.

    責任編輯:黃東杰

    国产一区二区在线观看日韩| 秋霞伦理黄片| 亚洲国产色片| 两个人视频免费观看高清| 中文精品一卡2卡3卡4更新| 欧美日本亚洲视频在线播放| 成人毛片60女人毛片免费| 欧美区成人在线视频| 岛国毛片在线播放| 99热精品在线国产| 一级爰片在线观看| 一级爰片在线观看| 日日摸夜夜添夜夜添av毛片| 人妻制服诱惑在线中文字幕| 欧美变态另类bdsm刘玥| 久久99热这里只有精品18| 尾随美女入室| 国产女主播在线喷水免费视频网站 | 水蜜桃什么品种好| 男插女下体视频免费在线播放| 精品免费久久久久久久清纯| 亚洲国产欧美在线一区| 欧美日韩在线观看h| 亚洲欧美日韩无卡精品| 国产高清不卡午夜福利| 欧美三级亚洲精品| 成年女人看的毛片在线观看| 亚洲欧美精品自产自拍| 精品人妻视频免费看| 亚洲av成人av| 嫩草影院新地址| 久久综合国产亚洲精品| 久久精品国产亚洲网站| 亚洲,欧美,日韩| 神马国产精品三级电影在线观看| 18禁在线播放成人免费| 女人被狂操c到高潮| 国产又色又爽无遮挡免| 国产三级中文精品| 麻豆国产97在线/欧美| 男人和女人高潮做爰伦理| 亚洲天堂国产精品一区在线| 亚洲成人av在线免费| 午夜a级毛片| 久久99热这里只频精品6学生 | 亚洲欧美精品自产自拍| 国产精品三级大全| 亚洲欧美一区二区三区国产| 小蜜桃在线观看免费完整版高清| 自拍偷自拍亚洲精品老妇| 一级毛片电影观看 | 在线a可以看的网站| 哪个播放器可以免费观看大片| av在线天堂中文字幕| 午夜激情福利司机影院| 不卡视频在线观看欧美| 热99在线观看视频| 美女cb高潮喷水在线观看| 国产精品永久免费网站| 成年女人看的毛片在线观看| 男女边吃奶边做爰视频| 免费黄色在线免费观看| 国产伦一二天堂av在线观看| 亚洲国产精品sss在线观看| 成人漫画全彩无遮挡| 国产精品美女特级片免费视频播放器| 国产亚洲一区二区精品| 性色avwww在线观看| 久久精品久久精品一区二区三区| 国产一区二区在线av高清观看| 变态另类丝袜制服| 国产高清视频在线观看网站| 少妇丰满av| 免费av毛片视频| 亚洲av免费在线观看| 国产一区亚洲一区在线观看| 亚洲av.av天堂| 3wmmmm亚洲av在线观看| 国产精品一区二区性色av| 中文亚洲av片在线观看爽| 日韩一本色道免费dvd| 国产三级在线视频| 国产国拍精品亚洲av在线观看| 欧美性猛交黑人性爽| 听说在线观看完整版免费高清| 亚洲精品国产成人久久av| 亚洲激情五月婷婷啪啪| 欧美日本亚洲视频在线播放| 一个人观看的视频www高清免费观看| 精品99又大又爽又粗少妇毛片| 岛国毛片在线播放| 亚洲av免费在线观看| 18禁在线无遮挡免费观看视频| 大香蕉97超碰在线| 小蜜桃在线观看免费完整版高清| 国产三级在线视频| 亚洲精品影视一区二区三区av| 毛片女人毛片| 亚洲不卡免费看| 欧美高清成人免费视频www| 国产精品国产三级国产av玫瑰| 一个人免费在线观看电影| 三级国产精品片| 国产欧美日韩精品一区二区| h日本视频在线播放| 在现免费观看毛片| 一级毛片我不卡| av卡一久久| or卡值多少钱| 精品免费久久久久久久清纯| 成人av在线播放网站| 高清在线视频一区二区三区 | 久久精品熟女亚洲av麻豆精品 | 又爽又黄无遮挡网站| 男人狂女人下面高潮的视频| 观看免费一级毛片| 国产精品乱码一区二三区的特点| av在线蜜桃| 91aial.com中文字幕在线观看| 亚洲国产精品合色在线| 真实男女啪啪啪动态图| 精品不卡国产一区二区三区| 蜜桃久久精品国产亚洲av| 直男gayav资源| 日韩视频在线欧美| 国产伦理片在线播放av一区| 搞女人的毛片| 精品一区二区三区人妻视频| 国产免费一级a男人的天堂| 桃色一区二区三区在线观看| 丝袜喷水一区| 91狼人影院| 麻豆久久精品国产亚洲av| 色噜噜av男人的天堂激情| a级毛片免费高清观看在线播放| 亚洲成人中文字幕在线播放| 1000部很黄的大片| 久久精品久久久久久噜噜老黄 | 日本一二三区视频观看| 国产精品久久视频播放| 亚洲aⅴ乱码一区二区在线播放| 青春草亚洲视频在线观看| 少妇的逼好多水| 亚洲一级一片aⅴ在线观看| 久久精品久久久久久噜噜老黄 | 男女边吃奶边做爰视频| 七月丁香在线播放| 午夜福利在线观看免费完整高清在| 国产淫语在线视频| 麻豆成人午夜福利视频| 最近中文字幕高清免费大全6| 色综合亚洲欧美另类图片| 日韩精品青青久久久久久| 一个人看的www免费观看视频| 亚洲三级黄色毛片| 晚上一个人看的免费电影| 岛国在线免费视频观看| 亚洲av一区综合| 我的女老师完整版在线观看| 九色成人免费人妻av| 91精品伊人久久大香线蕉| 青春草国产在线视频| 久久精品熟女亚洲av麻豆精品 | 国语自产精品视频在线第100页| 亚洲人与动物交配视频| 深爱激情五月婷婷| 麻豆精品久久久久久蜜桃| 别揉我奶头 嗯啊视频| 亚洲图色成人| 国产免费男女视频| 国产爱豆传媒在线观看| 亚洲自偷自拍三级| 欧美激情在线99| 国产片特级美女逼逼视频| 成人鲁丝片一二三区免费| 久久久久久久国产电影| 亚洲丝袜综合中文字幕| 91精品一卡2卡3卡4卡| 一级爰片在线观看| 成人亚洲欧美一区二区av| 91精品一卡2卡3卡4卡| 伦精品一区二区三区| 国语自产精品视频在线第100页| 午夜福利网站1000一区二区三区| 男人舔女人下体高潮全视频| 中文精品一卡2卡3卡4更新| 91在线精品国自产拍蜜月| 禁无遮挡网站| 午夜久久久久精精品| 国产精品不卡视频一区二区| 国产不卡一卡二| 亚洲va在线va天堂va国产| 美女被艹到高潮喷水动态| 2021天堂中文幕一二区在线观| 毛片一级片免费看久久久久| 国产在线男女| 亚洲av男天堂| 99热6这里只有精品| 日本熟妇午夜| 精品久久国产蜜桃| 久久精品久久久久久久性| 天堂av国产一区二区熟女人妻| 国产一区二区三区av在线| 亚洲婷婷狠狠爱综合网| 亚洲人成网站在线观看播放| 菩萨蛮人人尽说江南好唐韦庄 | 国产精品电影一区二区三区| 最新中文字幕久久久久| 午夜激情欧美在线| 国产白丝娇喘喷水9色精品| 亚洲精品色激情综合| 建设人人有责人人尽责人人享有的 | .国产精品久久| 亚洲国产欧美在线一区| 啦啦啦啦在线视频资源| 亚洲国产高清在线一区二区三| 色综合站精品国产| 国产精品一二三区在线看| 一级毛片aaaaaa免费看小| 精品不卡国产一区二区三区| 我的女老师完整版在线观看| 99久国产av精品国产电影| 日本三级黄在线观看| 亚洲,欧美,日韩| 日韩精品青青久久久久久| 国产亚洲91精品色在线| 性色avwww在线观看| 欧美精品国产亚洲| 寂寞人妻少妇视频99o| 久久这里只有精品中国| 午夜老司机福利剧场| 99久久人妻综合| 国产精品国产高清国产av| 欧美日韩在线观看h| 中国国产av一级| 午夜激情欧美在线| 欧美日本视频| 91av网一区二区| 美女xxoo啪啪120秒动态图| 国产精品美女特级片免费视频播放器| 欧美日韩一区二区视频在线观看视频在线 | 欧美人与善性xxx| av又黄又爽大尺度在线免费看 | 欧美又色又爽又黄视频| 久久久久网色| 亚洲成人av在线免费| 乱人视频在线观看| 夜夜爽夜夜爽视频| 婷婷色综合大香蕉| 非洲黑人性xxxx精品又粗又长| 波野结衣二区三区在线| 视频中文字幕在线观看| 亚洲av中文字字幕乱码综合| 在线免费观看的www视频| 干丝袜人妻中文字幕| 国产极品精品免费视频能看的| 日本av手机在线免费观看| 大又大粗又爽又黄少妇毛片口| 欧美日韩精品成人综合77777| 精品人妻熟女av久视频| 寂寞人妻少妇视频99o| 99国产精品一区二区蜜桃av| 色5月婷婷丁香| 国产成人精品婷婷| 又爽又黄a免费视频| 人人妻人人澡欧美一区二区| 精品一区二区三区视频在线| 亚洲成色77777| 亚洲av电影在线观看一区二区三区 | 一区二区三区免费毛片| 成人高潮视频无遮挡免费网站| 国产一区二区亚洲精品在线观看| 国产高清三级在线| 美女国产视频在线观看| 国产精品久久久久久久久免| 亚洲精品日韩在线中文字幕| 国内精品宾馆在线| 最近中文字幕2019免费版| 天美传媒精品一区二区| 看黄色毛片网站| 国产精品一及| av线在线观看网站| 久久精品久久久久久噜噜老黄 | 国产精品av视频在线免费观看| 秋霞在线观看毛片| 天天躁夜夜躁狠狠久久av| 亚洲va在线va天堂va国产| 免费看光身美女| 91久久精品国产一区二区三区| videos熟女内射| 婷婷色综合大香蕉| 国产单亲对白刺激| 久久精品国产自在天天线| 亚洲va在线va天堂va国产| 特级一级黄色大片| 亚洲最大成人av| 久久久久久久亚洲中文字幕| 久久久成人免费电影| 欧美3d第一页| 国产三级在线视频| 免费在线观看成人毛片| 91久久精品电影网| 亚洲国产最新在线播放| h日本视频在线播放| 我要搜黄色片| 免费av观看视频| 国产成人精品婷婷| 亚洲av成人精品一区久久| 成人国产麻豆网| 久久精品国产99精品国产亚洲性色| 久久久久久九九精品二区国产| 午夜免费激情av| 成人高潮视频无遮挡免费网站| 亚洲最大成人av| 成人一区二区视频在线观看| 国产av不卡久久| 直男gayav资源| 亚洲精品成人久久久久久| 麻豆成人午夜福利视频| 麻豆久久精品国产亚洲av| 六月丁香七月| 一个人看视频在线观看www免费| a级毛片免费高清观看在线播放| 天堂网av新在线| 色综合亚洲欧美另类图片| 中文字幕制服av| 2021天堂中文幕一二区在线观| 能在线免费观看的黄片| 日韩欧美 国产精品| 最后的刺客免费高清国语| 欧美不卡视频在线免费观看| 色播亚洲综合网| 伊人久久精品亚洲午夜| 免费观看性生交大片5| 国产精品人妻久久久久久| 国产黄a三级三级三级人| 国产在视频线在精品| 国产一区二区在线观看日韩| 亚洲欧美精品综合久久99| 国产午夜福利久久久久久| or卡值多少钱| 国产久久久一区二区三区| 人妻制服诱惑在线中文字幕| 国产久久久一区二区三区| 欧美日本视频| 亚洲国产欧美在线一区| 看非洲黑人一级黄片| 日日干狠狠操夜夜爽| 赤兔流量卡办理| 精品久久久久久电影网 | 岛国毛片在线播放| 51国产日韩欧美| 亚洲不卡免费看| 成年女人看的毛片在线观看| 欧美一区二区国产精品久久精品| 中文字幕免费在线视频6| 九九在线视频观看精品| 丰满乱子伦码专区| 久久草成人影院| 中文亚洲av片在线观看爽| 久久久久久大精品| 亚洲国产最新在线播放| 老女人水多毛片| 欧美又色又爽又黄视频| 六月丁香七月| 免费av观看视频| 日韩制服骚丝袜av| 午夜激情福利司机影院| 大话2 男鬼变身卡| 嫩草影院入口| 可以在线观看毛片的网站| 亚洲av成人精品一二三区| 欧美极品一区二区三区四区| 日本wwww免费看| 精品久久久久久成人av| 久久这里只有精品中国| 国产在线一区二区三区精 | 午夜激情福利司机影院| 亚洲真实伦在线观看| 麻豆乱淫一区二区| 欧美日韩在线观看h| 免费看av在线观看网站| 少妇被粗大猛烈的视频| 色尼玛亚洲综合影院| 久久99热这里只有精品18| 99久久成人亚洲精品观看| 美女国产视频在线观看| 成人美女网站在线观看视频| 国产亚洲午夜精品一区二区久久 | 久久精品国产鲁丝片午夜精品| 亚洲美女视频黄频| 精品人妻一区二区三区麻豆| av又黄又爽大尺度在线免费看 | 久久99热这里只有精品18| 国产在线一区二区三区精 | 蜜桃亚洲精品一区二区三区| 日韩中字成人| 欧美另类亚洲清纯唯美| 久久婷婷人人爽人人干人人爱| 色尼玛亚洲综合影院| 免费av毛片视频| 高清在线视频一区二区三区 | 成人亚洲精品av一区二区| 国产不卡一卡二| 九九爱精品视频在线观看| 久久亚洲精品不卡| 真实男女啪啪啪动态图| 国产成人福利小说| 欧美日本视频| 亚洲欧美日韩卡通动漫| 亚洲av二区三区四区| 好男人视频免费观看在线| 99视频精品全部免费 在线| 男插女下体视频免费在线播放| 国产一级毛片七仙女欲春2| 国产片特级美女逼逼视频| 国产爱豆传媒在线观看| 国产亚洲一区二区精品| 日韩中字成人| 中文在线观看免费www的网站| 国产成年人精品一区二区| 免费观看精品视频网站| 我要看日韩黄色一级片| 一个人看视频在线观看www免费| 久久6这里有精品| 日韩,欧美,国产一区二区三区 | 寂寞人妻少妇视频99o| 久久久国产成人免费| 一级黄色大片毛片| 亚洲欧美日韩卡通动漫| av在线亚洲专区| 国产精品爽爽va在线观看网站| 国产成人午夜福利电影在线观看| 一级毛片久久久久久久久女| 黄色一级大片看看| 最近手机中文字幕大全| 久久热精品热| 国产精品伦人一区二区| 日本wwww免费看| 十八禁国产超污无遮挡网站| 国国产精品蜜臀av免费| 国产午夜精品论理片| 69人妻影院| 夫妻性生交免费视频一级片| 国内精品宾馆在线| 久久久午夜欧美精品| 国产一级毛片在线| 久久婷婷人人爽人人干人人爱| 成人性生交大片免费视频hd| 欧美日韩精品成人综合77777| 亚洲国产色片| 国产精品一区二区在线观看99 | 在线a可以看的网站| 最近中文字幕高清免费大全6| 久久久久免费精品人妻一区二区| 乱人视频在线观看| 少妇熟女aⅴ在线视频| 亚洲婷婷狠狠爱综合网| 国产免费一级a男人的天堂| 亚洲成人久久爱视频| 亚洲乱码一区二区免费版| 午夜福利成人在线免费观看| 丰满人妻一区二区三区视频av| 精品国内亚洲2022精品成人| 18禁在线播放成人免费| 亚洲成色77777| 毛片女人毛片| 波多野结衣高清无吗| 免费看日本二区| 亚洲美女视频黄频| 国产成人精品一,二区| 有码 亚洲区| 在线播放无遮挡| 亚洲av熟女| 精品国产一区二区三区久久久樱花 | 亚洲精品乱码久久久v下载方式| 亚洲综合色惰| 99久久精品一区二区三区| 综合色丁香网| 在线观看66精品国产| 亚洲欧美日韩高清专用| 国产精品国产三级专区第一集| 久热久热在线精品观看| 国产精品,欧美在线| 天堂av国产一区二区熟女人妻| 一级毛片久久久久久久久女| 精品不卡国产一区二区三区| 国产美女午夜福利| 欧美日韩精品成人综合77777| 级片在线观看| 国产精品日韩av在线免费观看| 国产成人一区二区在线| videossex国产| 国产免费一级a男人的天堂| a级毛色黄片| 精品国产三级普通话版| 亚洲天堂国产精品一区在线| 九九热线精品视视频播放| 国产成人一区二区在线| 亚洲国产最新在线播放| 91午夜精品亚洲一区二区三区| 蜜桃亚洲精品一区二区三区| 久久久久久久久久黄片| 中文资源天堂在线| videossex国产| 22中文网久久字幕| 国产一级毛片在线| 国产成年人精品一区二区| kizo精华| 国产精品久久久久久av不卡| 国产爱豆传媒在线观看| 美女高潮的动态| 成人国产麻豆网| 国产免费又黄又爽又色| 免费av毛片视频| 久久久a久久爽久久v久久| 国产极品精品免费视频能看的| 九九久久精品国产亚洲av麻豆| 久久精品国产99精品国产亚洲性色| 人人妻人人澡欧美一区二区| 久久久久性生活片| 国产高潮美女av| 能在线免费观看的黄片| 欧美一区二区精品小视频在线| 成人午夜高清在线视频| 国产免费男女视频| 久久久久久国产a免费观看| 国产一区有黄有色的免费视频 | 最近2019中文字幕mv第一页| 日韩视频在线欧美| 美女国产视频在线观看| 久久精品国产亚洲网站| 岛国毛片在线播放| 一区二区三区四区激情视频| 日韩国内少妇激情av| 日韩成人av中文字幕在线观看| 亚洲性久久影院| 综合色丁香网| 免费看光身美女| 99国产精品一区二区蜜桃av| 欧美激情久久久久久爽电影| 两个人视频免费观看高清| 少妇高潮的动态图| 性插视频无遮挡在线免费观看| 99久久九九国产精品国产免费| 一个人观看的视频www高清免费观看| 色综合色国产| 毛片一级片免费看久久久久| 黄色日韩在线| 最后的刺客免费高清国语| 久久精品91蜜桃| 2021少妇久久久久久久久久久| 美女脱内裤让男人舔精品视频| 丰满少妇做爰视频| 大话2 男鬼变身卡| 亚洲最大成人手机在线| 国产极品精品免费视频能看的| 国产真实乱freesex| 成人三级黄色视频| 亚洲在线自拍视频| 久久99热这里只频精品6学生 | 欧美性猛交╳xxx乱大交人| 网址你懂的国产日韩在线| 1024手机看黄色片| 久久精品久久久久久噜噜老黄 | 欧美zozozo另类| ponron亚洲| 久久精品国产鲁丝片午夜精品| 青春草亚洲视频在线观看| 夜夜看夜夜爽夜夜摸| 嘟嘟电影网在线观看| 一边亲一边摸免费视频| 国产黄片美女视频| 日韩一区二区三区影片| 爱豆传媒免费全集在线观看| 久久鲁丝午夜福利片| 熟女人妻精品中文字幕| 色综合亚洲欧美另类图片| 亚洲性久久影院| 欧美色视频一区免费| 午夜a级毛片| 日韩在线高清观看一区二区三区| 亚洲综合精品二区| 久久99精品国语久久久| 成人美女网站在线观看视频| 两个人视频免费观看高清| 22中文网久久字幕| 国产精品永久免费网站| 成人漫画全彩无遮挡| 国产午夜福利久久久久久| 中文字幕制服av| 欧美一区二区国产精品久久精品| 小说图片视频综合网站| 99在线视频只有这里精品首页| 日韩精品青青久久久久久| 国产一级毛片在线| 大又大粗又爽又黄少妇毛片口| 亚洲美女视频黄频| 直男gayav资源| 九色成人免费人妻av| 中国美白少妇内射xxxbb| 深爱激情五月婷婷| 国产精品国产高清国产av| av国产久精品久网站免费入址| 成年女人看的毛片在线观看| 国产精品久久久久久精品电影小说 | 久久久久免费精品人妻一区二区| 岛国在线免费视频观看| 在线观看66精品国产| 1024手机看黄色片| 日韩精品青青久久久久久| 91精品国产九色| 99热网站在线观看| 成人欧美大片| 舔av片在线|