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

    A Brief Review of Software Tools for Pangenomics

    2015-02-06 08:00:04JingfaXiaoZhewenZhangJiayanWuJunYu
    Genomics,Proteomics & Bioinformatics 2015年1期

    Jingfa Xiao*,Zhewen Zhang,Jiayan Wu,Jun Yu

    CAS Key Laboratory of Genome Sciences and Information,Beijing Institute of Genomics,Chinese Academy of Sciences, Beijing 100101,China

    A Brief Review of Software Tools for Pangenomics

    Jingfa Xiao*,a,Zhewen Zhangb,Jiayan Wuc,Jun Yud

    CAS Key Laboratory of Genome Sciences and Information,Beijing Institute of Genomics,Chinese Academy of Sciences, Beijing 100101,China

    Pangenomics;

    Pangenome;

    Comparative analysis;

    Genomic dynamics;

    Core genes

    Since the proposal for pangenomic study,there have been a dozen software tools actively in use for pangenomic analysis.By the end of 2014,Panseq and the pan-genomes analysis pipeline (PGAP)ranked as the top two most popular packages according to cumulative citations of peerreviewed scientifc publications.The functions of the software packages and tools,albeit variable among them,include categorizing orthologous genes,calculating pangenomic profles,integrating gene annotations,and constructing phylogenies.As epigenomic elements are being gradually revealed in prokaryotes,it is expected that pangenomic databases and toolkits have to be extended to handle information of detailed functional annotations for genes and non-protein-coding sequences including non-coding RNAs,insertion elements,and conserved structural elements.To develop better bioinformatic tools,user feedback and integration of novel features are both of essence.

    Introduction

    In the past decade or so,the remarkable advancement of DNA sequencing technology and application has led to an astronomical accumulation of genomic data.This is especially true for the prokaryotic genomes as individual of them is only a few megabases in size.It is expected that in the next decade or two, there will be more data collected than what we can actually handle.Therefore,database construction,improvement,and

    *Corresponding author.

    E-mail:xiaojf@big.ac.cn(Xiao J).

    aORCID:0000-0002-2835-4340.

    bORCID:0000-0002-9422-822X.

    cORCID:0000-0001-6048-405X.

    dORCID:0000-0002-2702-055X.

    Peer review under responsibility of Beijing Institute of Genomics, Chinese Academy of Sciences and Genetics Society of China. consolidation,as well as new tool development,are especially welcome.In this way,the sibling felds of genomics,such as pangenomics and metagenomics,can all be ready for curating, sharing,and mining foods of the incoming genomic big data. Coming back to the reality and focusing on pangenomics, there were,as of December 2014,more than 40 bacterial species that have over 20 fully-assembled genomes from different strains and isolates,allowing for comprehensive pangenomic studies.The concept of pangenome was frst proposed in 2005 by Tettelin et al.[1,2],which is defned as the entire genomic repertoire of a given species or phylogenetic clade when multiple species are defned by systematics.According to the defnition,gene profle(content)of a pangenome is divided into three groups:core(shared by all genomes),dispensable, and strain-(or isolate-)specifc genes.A series of pangenomic studies have been performed in genomic dynamics[3–6], pathogenesis and drug resistance[7–9],bacterial toxins[10],and species evolution[11].The concept has also been extended to viral[12],plant[13–15],and fungal genome studies[16].A review on ten-year history and feld achievement of pangenomics has just been published at the beginning of 2014[2], which detailed major projects as well as methodology and technology advancements.

    Here,we provide a brief review on the pangenomic software packages and tools,including their basic function,general utility,and popularity based on their cumulative citation by peer-reviewed scientifc publications.Although such a singlecriterion evaluation may never be adequate and thorough,we hope that it provides a feld guide for students and young scientists to make the right choice for their preferred applications.

    Highlights of the software packages and tools

    Since 2010,we have seen a dozen or so software packages and tools being put forward,which are capable of clustering orthologous genes,identifying single nucleotide polymorphisms(SNPs),constructing phylogenies,and profling core/ shared/isolate-specifc genes.Although they may share similar functions,each has its own characteristics and limitations, leaving rooms for further improvement.

    Among the early-developed packages,Panseq[17]and PanCGHweb[18]were published in 2010,followed by CAMBer[19]and the Prokaryotic-genome Analysis Tool (PGAT)[20]in 2011.PanCGHweb is a web tool for pangenomic microarray analysis based on PanCGH algorithm[21]. It enables users to group genes into orthologs and to construct gene-based phylogenies ofrelated strains and isolates. However,this package is rather specifc for handling microarray data but not RNA-seq data.Panseq,another online pangenomic tool,is able to determine core and accessory regions of genome assemblies based on MUMmer and BLASTn,as well as to identify SNPs among the core genomic regions.In addition,Panseq also has a locus selector module that selects the most discriminatory loci among the accessory loci or core gene SNPs[17].Panseq,however,is not able to provide pangenomic profle and functional enrichment analysis that is important for the biologists to flter out functional relevance of the pangenomic elements.The later released CAMBer is designed to identify multi-gene families from multiple bacterial strains and isolates.These multi-gene families can be used for sequencing errordetection,mutation identifcation,and pangenomic profle computation[19].CAMBer is supreme in refning gene function prediction according to multi-gene family information,but it does not provide tools for comparative or evolutionary analysis among strains and isolates.As a webbased database,PGAT integrates several useful functions, such as plotting the presence and absence of genes among members of a pangenome,identifying SNPs among orthologs and syntenic regions,comparing gene orders among different strains and isolates,providing KEGG pathway analysis tools, and searching for genes through different annotations such as the Cluster of Orthologous Groups of proteins(COG), PSORT,SignalP,the Tied Mixture Hidden Markov Model (TMHMM),and Pfam.However,PGAT is just a database with a limited number of species curated and it cannot perform analysis for new sequencing data from users.

    PGAP is a stand-alone program developed by Zhao et al.in 2012,which contains fve functional models[22].Based on functional gene clustering and analysis,PGAP presents pangenomic profle(partitions of pangenomic elements or gene categories),genetic variation,species evolution,and function enrichment of different strains and isolates of a given pangenome.In addition,all analyses are performed with a single command,and such integration is rather user-friendly and effcient. Nonetheless,PGAP has its limitation as well.For instance,all its output fles of the fve models are text fles,which lacks of intuitiveness.Contreras-Moreira et al.subsequently proposed a program called GET_HOMOLOGUES in 2013,which is also a versatile software package for pangenomics[23].This software package integrated data download,sequence feature extraction,homologous gene identifcation,pangenome profling,graphical display,and phylogenetic tree construction into one powerful toolkit.Several other tools were also available in 2013,such as PanCake[24]and PANNOTATOR[25]. PanCake was developed for identifying singletons and core regions in arbitrary sequence sets,while PANNOTATOR,a web-based automated pipeline,was designed for the annotation of closely-related genomes for pangenomic analysis.However, these two tools only focus on simple functions,such as clustering homologous genes and gene curation.In 2014,a powerful and fexible toolkit,the Integrated Toolkit for Exploration of microbial Pan-genomes(ITEP),was published by Benedict and colleagues[26].ITEP integrates plenty of existing bioinformatics tools for pangenomic analysis,including protein family prediction,ortholog detection,functional domain analysis, pangenomic profling,and metabolic network integration. Moreover,ITEP also integrates some visualization scripts that assist biologists in phylogenetic tree construction,annotation curation,and specifc query for conserved protein domain identifcation.In 2014,another rapid core-genome alignment and visualization pangenomic software package,Harvest,was proposed by Treangen et al.Harvest contains tools,such as Parsnp and Gingr,for core gene alignment,variant calling,recombination detection,and phylogenetic trees construction[27].To analyzepangenomic profle in a larger scale,a software package PanGP was developed with a graphic interface by Zhao et al.in 2014[28].Spine and AGEnt were also developed in 2014,which are capable of profling pangenomes based on both fnished and draft genomic sequences[29].We summarized all the software packages and tools inTable 1,highlighting their platforms and major features.We went one step further and ranked them according to their citations by peer-reviewed scientifc publications(Figure 1),which were collected from ISI Web of Science-Science Citation Index Expanded.Our summary indicates that Panseq and PGAP have been the most popular packages up to the end of 2014.

    A wish list for improving the current software

    Although single-tool solution could not usually satisfy the need for understanding the whole picture,a wish list from the users is always helpful for prioritizing goals for the package developers,hence providing directions for improving each package.

    First,the performance of pangenomic analysis strongly depends on the accuracy of genome assembly and annotation. Therefore,an adequate number of complete sequence assemblies are a prerequisite.Currently,most of the existing bacterial genome sequences are actually incomplete(in most of thecases,contigs are not joined together into single chromosomes),and some only have high-quality and high-coverage raw data available.The inclusion of incomplete genome assemblies for pangenomic analysis may need scaffold building that requires reformatting of the contig data fles.Despite the development of the third-generation sequencing technology,which would certainly help the assembly and fnishing of prokaryotic genomes[30],incomplete prokaryotic genomes are expected to be deposited into public databases in mass.It would be a waste if such data are left unused.

    Table 1 Software tools for pangenomic studies

    Figure 1 Relative citation of the pangenomic software tools from peer-reviewed scientifc publications

    Second,orthologous gene identifcation is a key step in pangenomic analysis.At present,the existing software for ortholog detection is mainly based on sequence similarity, phylogenetic relationship,or other annotation information such as functional information.The development of novel and more effcient ortholog identifcation method for multiple closely-related strains and isolates can greatly improve the accuracy of pangenomic analysis.One possibility is to integrate gene gain-and-loss information for phylogeny building among strains and isolates.

    Third,sampling is also important for pangenomics in a couple of counts.One is how many strains or isolates to choose for a pangenomic analysis.The other is how to implement a flter that differentiates more diverse strains or isolates from the less diverse for pangenomic analysis.For instance,if we choose all genomes of a species for an analysis,which include one or a few divergent genomes,the core genome will be much shorter or reduced.Obviously,individual genomes should be selected and regrouped for better representation of average nucleotide identity(ANI).ANI is one of the most useful measurements for species delineation[31].Therefore,for a better pangenomic analysis,detailed information for the available samples is of essence,which should include their genotypes, phenotypes,and habitats.

    Fourth,the current tools have not incorporated some recent advancements in prokaryotic genomics,such as the so-called genome-organization frameworks(GOFs),which are not only unique to each species but also provide guidance for sequence assembly and fnishing[32].Other annotation information,such as that of non-coding RNAs,pseudogenes, and epigenetic elements,remains to be implemented into the relevantsoftwarepackages.Finally,a never-ending improvementofpangenomictoolsisvisualization that provides not only better displays but also quality graphics for publication.

    Concluding remarks

    We provide an overview on the existing pangenomic analysis tools and hope to see improvements of the software tools from their original developers.We certainly express our enthusiasm for new tools to join the competition,and after all,for a piece of bioinformatic work,a database or a toolkit,the survival or winning game is in its long-term maintenance and constant improvement.

    Competing interests

    The authors declared that there are no competing interests.

    Acknowledgements

    This study was supported by the National High-tech R&D Program(863 Program;Grant No.2012AA020409)from theMinistry of Science and Technology of China,the Key Program of the Chinese Academy of Sciences(Grant No. KSZD-EW-TZ-009-02),and the National Natural Science Foundation of China(Grant Nos.31471248 and 31271386).

    [1]Tettelin H,Masignani V,Cieslewicz MJ,Donati C,Medini D, Ward NL,et al.Genome analysis of multiple pathogenic isolates ofStreptococcus agalactiae:implications for the microbial‘‘pangenome’’.Proc Natl Acad Sci U S A 2005;102:13950–5.

    [2]Vernikos G,Medini D,Riley DR,Tettelin H.Ten years of pangenome analyses.Curr Opin Microbiol 2014;23C:148–54.

    [3]Rasmussen TB,Danielsen M,Valina O,Garrigues C,Johansen E, Pedersen MB.Streptococcus thermophiluscore genome:comparative genome hybridization study of 47 strains.Appl Environ Microbiol 2008;74:4703–10.

    [4]WillenbrockH,HallinPF,WassenaarTM,UsseryDW. Characterization of probioticEscherichia coliisolates with a novel pan-genome microarray.Genome Biol 2007;8:R267.

    [5]Zhang A,Yang M,Hu P,Wu J,Chen B,Hua Y,et al. Comparative genomic analysis ofStreptococcus suisreveals signifcant genomic diversity among different serotypes.BMC Genomics 2011;12:523.

    [6]KaasRS,FriisC,UsseryDW,AarestrupFM.Estimatingvariation within the genes and inferring the phylogeny of 186 sequenced diverseEscherichia coligenomes.BMC Genomics 2012;13:577.

    [7]Park J,Zhang Y,Buboltz AM,Zhang X,Schuster SC,Ahuja U, et al.Comparative genomics of the classicalBordetellasubspecies: the evolution and exchange of virulence-associated diversity amongst closely related pathogens.BMC Genomics 2012;13:545.

    [8]D’Auria G,Jimenez-Hernandez N,Peris-Bondia F,Moya A, Latorre A.Legionella pneumophilapangenome reveals strainspecifc virulence factors.BMC Genomics 2010;11:181.

    [9]Hu P,Yang M,Zhang A,Wu J,Chen B,Hua Y,et al. Comparative genomics study of multi-drug-resistance mechanisms in the antibiotic-resistantStreptococcus suisR61 strain.PLoS One 2011;6:e24988.

    [10]Fang Y,Li Z,Liu J,Shu C,Wang X,Zhang X,et al.A pangenomic study ofBacillus thuringiensis.J Genet Genomics 2011;38:567–76.

    [11]Konstantinidis KT,Ramette A,Tiedje JM.The bacterial species defnition in the genomic era.Philos Trans R Soc Lond B Biol Sci 2006;361:1929–40.

    [12]AherfS,Pagnier I,Fournous G,Raoult D,La Scola B,Colson P. CompletegenomesequenceofCannes8virus,anewmemberofthe proposed family‘‘Marseilleviridae’’.Virus Genes 2013;47:550–5.

    [13]Cao J,Schneeberger K,Ossowski S,Gunther T,Bender S,Fitz J, et al.Whole-genome sequencing of multipleArabidopsis thalianapopulations.Nat Genet 2011;43:956–63.

    [14]Read BA,Kegel J,Klute MJ,Kuo A,Lefebvre SC,Maumus F, et al.Pan genome of the phytoplanktonEmilianiaunderpins its global distribution.Nature 2013;499:209–13.

    [15]Li YH,Zhou G,Ma J,Jiang W,Jin LG,Zhang Z,et al.De novoassembly of soybean wild relatives for pan-genome analysis of diversity and agronomic traits.Nat Biotechnol 2014;32:1045–52.

    [16]Dunn B,Richter C,Kvitek DJ,Pugh T,Sherlock G.Analysis of theSaccharomyces cerevisiaepan-genome reveals a pool of copy number variants distributed in diverse yeast strains from differing industrial environments.Genome Res 2012;22:908–24.

    [17]Laing C,Buchanan C,Taboada EN,Zhang YX,Kropinski A, Villegas A,et al.Pan-genome sequence analysis using Panseq:an online tool for the rapid analysis of core and accessory genomic regions.BMC Bioinformatics 2010;11:461.

    [18]Bayjanov JR,Siezen RJ,van Hijum SAFT.PanCGHweb:a web tool for genotype calling in pangenome CGH data.Bioinformatics 2010;26:1256–7.

    [19]Wozniak M,Wong L,Tiuryn J.CAMBer:an approach to support comparative analysis of multiple bacterial strains.BMC Genomics 2011;12:S6.

    [20]Brittnacher MJ,Fong C,Hayden HS,Jacobs MA,Radey M, Rohmer L.PGAT:a multistrain analysis resource for microbial genomes.Bioinformatics 2011;27:2429–30.

    [21]Bayjanov JR,Wels M,Starrenburg M,van Hylckama Vlieg JE, Siezen RJ,Molenaar D.PanCGH:a genotype-calling algorithm for pangenome CGH data.Bioinformatics 2009;25:309–14.

    [22]Zhao Y,Wu J,Yang J,Sun S,Xiao J,Yu J.PGAP:pan-genomes analysis pipeline.Bioinformatics 2012;28:416–8.

    [23]Contreras-Moreira B,Vinuesa P.GET_HOMOLOGUES,a versatile software package for scalable and robust microbial pangenome analysis.Appl Environ Microbiol 2013;79:7696–701.

    [24]Ernst C,Rahmann S.PanCake:a data structure for pangenomes. In:Bei?barth T,Kollmar M,Leha A,Morgenstern B,Schultz A, editors.German conference on bioinformatics 2013.Dagstuhl, Germany:Schloss Dagstuhl-Leibniz-Zentrum fuer Informatik; 2013.p.35–45.http://dx.doi.org/10.4230/OASIcs.GCB.2013.35.

    [25]Santos AR,Barbosa E,Fiaux K,Zurita-Turk M,Chaitankar V, Kamapantula B,et al.PANNOTATOR:an automated tool for annotation of pan-genomes.Genet Mol Res 2013;12:2982–9.

    [26]Benedict MN,Henriksen JR,Metcalf WW,Whitaker RJ,Price ND.ITEP:an integrated toolkit for exploration of microbial pangenomes.BMC Genomics 2014;15:8.

    [27]Treangen TJ,Ondov BD,Koren S,Phillippy AM.The Harvest suite for rapid core-genome alignment and visualization of thousands of intraspecifc microbial genomes.Genome Biol 2014;15:524.

    [28]Zhao Y,Jia X,Yang J,Ling Y,Zhang Z,Yu J,et al.PanGP:a toolfor quickly analyzing bacterialpan-genome profle. Bioinformatics 2014;30:1297–9.

    [29]Ozer EA,Allen JP,Hauser AR.Characterization of the core and accessory genomes ofPseudomonas aeruginosausing bioinformatic tools Spine and AGEnt.BMC Genomics 2014;15:737.

    [30]Koren S,Harhay GP,Smith TP,Bono JL,Harhay DM,Mcvey SD,et al.Reducing assembly complexity of microbial genomes with single-molecule sequencing.Genome Biol 2013;14:R101.

    [31]Kim M,Oh HS,Park SC,Chun J.Towards a taxonomic coherence between average nucleotide identity and 16S rRNA gene sequence similarity for species demarcation of prokaryotes. Int J Syst Evol Microbiol 2014;64:346–51.

    [32]Kang Y,Gu C,Yuan L,Wang Y,Zhu Y,Li X,et al.Flexibility and symmetry of prokaryotic genome rearrangement reveal lineage-associated core-gene-defned genome organizational frameworks.MBio 2014;5:e01867.

    Received 9 January 2015;revised 22 January 2015;accepted 25 January 2015

    Available online 23 February 2015

    Handled by Vladimir Bajic

    http://dx.doi.org/10.1016/j.gpb.2015.01.007

    1672-0229?2015 The Authors.Production and hosting by Elsevier B.V.on behalf of Beijing Institute of Genomics,Chinese Academy of Sciences and Genetics Society of China.

    This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    久久精品国产综合久久久| 国产aⅴ精品一区二区三区波| 日韩免费高清中文字幕av| 免费观看a级毛片全部| 91精品三级在线观看| 黄网站色视频无遮挡免费观看| 黑人操中国人逼视频| 老熟妇乱子伦视频在线观看| av电影中文网址| √禁漫天堂资源中文www| 国产精品 欧美亚洲| 国产日韩一区二区三区精品不卡| 国产亚洲欧美98| 18禁裸乳无遮挡动漫免费视频| 三级毛片av免费| 后天国语完整版免费观看| 国产精品自产拍在线观看55亚洲 | 日韩有码中文字幕| aaaaa片日本免费| 国产激情久久老熟女| 在线播放国产精品三级| 窝窝影院91人妻| 欧美日韩国产mv在线观看视频| 久久精品国产综合久久久| www日本在线高清视频| 少妇 在线观看| 下体分泌物呈黄色| 国产精品免费视频内射| 国产欧美日韩一区二区三区在线| 99国产极品粉嫩在线观看| 日韩精品免费视频一区二区三区| 人妻 亚洲 视频| 99热只有精品国产| 黑丝袜美女国产一区| 满18在线观看网站| 久久午夜亚洲精品久久| 久久性视频一级片| 99久久精品国产亚洲精品| 国产成人啪精品午夜网站| 国产精品电影一区二区三区 | 成人免费观看视频高清| 久久精品人人爽人人爽视色| 最近最新中文字幕大全免费视频| 交换朋友夫妻互换小说| 老熟女久久久| 一级黄色大片毛片| 自线自在国产av| 午夜福利欧美成人| 亚洲精品美女久久av网站| 一夜夜www| 亚洲色图 男人天堂 中文字幕| 午夜精品在线福利| 三级毛片av免费| 欧美一级毛片孕妇| 我的亚洲天堂| av电影中文网址| 国产成+人综合+亚洲专区| 亚洲精品成人av观看孕妇| 女同久久另类99精品国产91| 国产深夜福利视频在线观看| 美女午夜性视频免费| 99在线人妻在线中文字幕 | 久久热在线av| 欧美最黄视频在线播放免费 | av网站在线播放免费| 人人澡人人妻人| 在线国产一区二区在线| 成人永久免费在线观看视频| 亚洲精品在线美女| 啦啦啦免费观看视频1| 日本黄色日本黄色录像| 久久精品国产亚洲av高清一级| 视频在线观看一区二区三区| 香蕉国产在线看| 成人三级做爰电影| 国精品久久久久久国模美| www日本在线高清视频| 999久久久国产精品视频| 国产主播在线观看一区二区| 亚洲av电影在线进入| 国产精品久久视频播放| 日韩一卡2卡3卡4卡2021年| 最新美女视频免费是黄的| 无限看片的www在线观看| 美女视频免费永久观看网站| 99re在线观看精品视频| 老司机影院毛片| 午夜精品在线福利| 免费久久久久久久精品成人欧美视频| 国精品久久久久久国模美| 老熟妇乱子伦视频在线观看| 免费看a级黄色片| 亚洲久久久国产精品| 久久精品国产亚洲av高清一级| 欧美人与性动交α欧美精品济南到| 女人被狂操c到高潮| 国产99久久九九免费精品| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩亚洲高清精品| 巨乳人妻的诱惑在线观看| 露出奶头的视频| 亚洲色图 男人天堂 中文字幕| 成熟少妇高潮喷水视频| av国产精品久久久久影院| 欧美性长视频在线观看| 日韩免费av在线播放| a级片在线免费高清观看视频| 亚洲性夜色夜夜综合| 欧美乱妇无乱码| 欧美中文综合在线视频| 韩国av一区二区三区四区| 国产精品综合久久久久久久免费 | 国产成人啪精品午夜网站| 女人爽到高潮嗷嗷叫在线视频| 麻豆乱淫一区二区| 欧美在线黄色| 免费高清在线观看日韩| 99精品久久久久人妻精品| 久久影院123| 欧美 亚洲 国产 日韩一| 女人被躁到高潮嗷嗷叫费观| 一进一出好大好爽视频| 狠狠狠狠99中文字幕| 精品无人区乱码1区二区| 最新美女视频免费是黄的| 亚洲国产精品合色在线| 国产欧美日韩一区二区三区在线| 午夜福利一区二区在线看| 一夜夜www| 两人在一起打扑克的视频| 欧美最黄视频在线播放免费 | 国产aⅴ精品一区二区三区波| 欧美 亚洲 国产 日韩一| 18禁观看日本| 18禁国产床啪视频网站| 国产精品久久久久久人妻精品电影| 久久这里只有精品19| xxxhd国产人妻xxx| 久久久精品区二区三区| 99re在线观看精品视频| 国产伦人伦偷精品视频| 美女 人体艺术 gogo| 黄色a级毛片大全视频| 黄色成人免费大全| 亚洲欧美色中文字幕在线| 人人妻人人澡人人爽人人夜夜| 欧美成狂野欧美在线观看| 免费av中文字幕在线| 国产在视频线精品| 国产在线一区二区三区精| 亚洲av美国av| 中国美女看黄片| √禁漫天堂资源中文www| 午夜福利影视在线免费观看| 悠悠久久av| 丁香欧美五月| 欧美激情极品国产一区二区三区| 久久亚洲真实| 操出白浆在线播放| 制服人妻中文乱码| 欧美乱码精品一区二区三区| 18禁观看日本| 99re6热这里在线精品视频| 交换朋友夫妻互换小说| 99香蕉大伊视频| 免费久久久久久久精品成人欧美视频| 一本综合久久免费| 精品电影一区二区在线| 欧美激情 高清一区二区三区| 亚洲七黄色美女视频| 精品欧美一区二区三区在线| 99久久99久久久精品蜜桃| 身体一侧抽搐| 人妻 亚洲 视频| 51午夜福利影视在线观看| 日本wwww免费看| av国产精品久久久久影院| av免费在线观看网站| 日韩 欧美 亚洲 中文字幕| 国产免费现黄频在线看| bbb黄色大片| 99国产精品99久久久久| 亚洲国产欧美日韩在线播放| 亚洲精品一二三| 交换朋友夫妻互换小说| 精品久久蜜臀av无| 香蕉久久夜色| 天天躁夜夜躁狠狠躁躁| 国产午夜精品久久久久久| 久久午夜亚洲精品久久| 精品国产乱子伦一区二区三区| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 国产一卡二卡三卡精品| 母亲3免费完整高清在线观看| av片东京热男人的天堂| 免费在线观看日本一区| 国产亚洲精品久久久久久毛片 | 久久精品国产99精品国产亚洲性色 | 免费在线观看日本一区| 国产亚洲精品久久久久久毛片 | 亚洲精品国产一区二区精华液| 啦啦啦视频在线资源免费观看| 亚洲精品在线美女| 亚洲精品一二三| 亚洲精品自拍成人| 成人三级做爰电影| 欧美性长视频在线观看| 亚洲欧美色中文字幕在线| 亚洲一码二码三码区别大吗| 高清av免费在线| av电影中文网址| 下体分泌物呈黄色| 好男人电影高清在线观看| 欧美精品高潮呻吟av久久| 久久精品国产清高在天天线| 视频在线观看一区二区三区| 成人影院久久| av免费在线观看网站| 国产不卡av网站在线观看| 国产色视频综合| 啦啦啦在线免费观看视频4| 天天躁夜夜躁狠狠躁躁| 老司机影院毛片| 精品无人区乱码1区二区| 91九色精品人成在线观看| 无遮挡黄片免费观看| 精品免费久久久久久久清纯 | 人人妻,人人澡人人爽秒播| 亚洲av成人av| 后天国语完整版免费观看| 免费一级毛片在线播放高清视频 | 国内久久婷婷六月综合欲色啪| 一二三四在线观看免费中文在| 久久国产乱子伦精品免费另类| 中文字幕另类日韩欧美亚洲嫩草| 欧美人与性动交α欧美精品济南到| 两个人看的免费小视频| 亚洲五月色婷婷综合| 国产男靠女视频免费网站| 亚洲精品美女久久久久99蜜臀| 久久亚洲真实| 国产高清videossex| 亚洲欧美精品综合一区二区三区| 亚洲国产毛片av蜜桃av| 国产三级黄色录像| a级毛片在线看网站| 亚洲精品国产色婷婷电影| 亚洲熟女毛片儿| 天堂√8在线中文| 精品久久久久久电影网| 咕卡用的链子| 国产一区在线观看成人免费| 一边摸一边抽搐一进一小说 | 亚洲伊人色综图| 亚洲男人天堂网一区| 国产一卡二卡三卡精品| 欧美激情 高清一区二区三区| 亚洲av日韩精品久久久久久密| 91国产中文字幕| 欧美最黄视频在线播放免费 | 午夜91福利影院| 亚洲伊人色综图| 国产成人精品久久二区二区91| 成人三级做爰电影| 国产精品98久久久久久宅男小说| 午夜福利欧美成人| 亚洲av第一区精品v没综合| 欧美一级毛片孕妇| 激情在线观看视频在线高清 | 国产区一区二久久| 欧美日韩av久久| 久久狼人影院| 夜夜爽天天搞| www.精华液| 精品一区二区三卡| 国产xxxxx性猛交| 亚洲aⅴ乱码一区二区在线播放 | 好看av亚洲va欧美ⅴa在| 身体一侧抽搐| 在线永久观看黄色视频| 丝袜人妻中文字幕| 成人影院久久| 亚洲中文字幕日韩| 99精品欧美一区二区三区四区| 一夜夜www| 91麻豆av在线| 亚洲成人免费电影在线观看| 久热这里只有精品99| 精品免费久久久久久久清纯 | av天堂在线播放| 久久中文字幕人妻熟女| 可以免费在线观看a视频的电影网站| 精品一区二区三区av网在线观看| 午夜精品国产一区二区电影| 亚洲,欧美精品.| 欧美激情极品国产一区二区三区| 女警被强在线播放| 高清欧美精品videossex| 国产精品久久久久成人av| 国产xxxxx性猛交| tocl精华| 99国产精品99久久久久| 在线十欧美十亚洲十日本专区| 999久久久精品免费观看国产| 日韩欧美一区视频在线观看| 久久久国产成人免费| 久久精品国产综合久久久| 国产不卡一卡二| 如日韩欧美国产精品一区二区三区| 丰满迷人的少妇在线观看| 在线国产一区二区在线| 女人高潮潮喷娇喘18禁视频| 在线免费观看的www视频| av福利片在线| 无限看片的www在线观看| 日韩欧美一区视频在线观看| 午夜精品久久久久久毛片777| 久久久国产一区二区| 淫妇啪啪啪对白视频| www.自偷自拍.com| 中文字幕另类日韩欧美亚洲嫩草| 国产成人免费观看mmmm| 欧美日韩瑟瑟在线播放| 久久精品国产亚洲av香蕉五月 | 成人黄色视频免费在线看| 中文字幕制服av| 日本黄色日本黄色录像| 这个男人来自地球电影免费观看| 欧美久久黑人一区二区| av不卡在线播放| 国产在线观看jvid| a级毛片在线看网站| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩另类电影网站| 五月开心婷婷网| 日日夜夜操网爽| 久久久精品区二区三区| 日韩欧美一区二区三区在线观看 | 99精品久久久久人妻精品| 欧美日韩瑟瑟在线播放| 国产欧美日韩一区二区三| 香蕉国产在线看| 国产亚洲精品一区二区www | 国产av又大| 好看av亚洲va欧美ⅴa在| 在线免费观看的www视频| av福利片在线| 亚洲 欧美一区二区三区| 成人国产一区最新在线观看| 亚洲aⅴ乱码一区二区在线播放 | 欧美午夜高清在线| 色精品久久人妻99蜜桃| 两个人看的免费小视频| 午夜免费成人在线视频| 桃红色精品国产亚洲av| 一本综合久久免费| 好男人电影高清在线观看| cao死你这个sao货| 久久午夜综合久久蜜桃| 亚洲专区中文字幕在线| 亚洲性夜色夜夜综合| 极品人妻少妇av视频| 亚洲专区国产一区二区| 欧美黑人精品巨大| 精品无人区乱码1区二区| 日本黄色视频三级网站网址 | 黄色a级毛片大全视频| 国产亚洲精品第一综合不卡| 亚洲九九香蕉| 高潮久久久久久久久久久不卡| 久久久精品区二区三区| 久久人妻福利社区极品人妻图片| 国产精品久久电影中文字幕 | 日本五十路高清| 国产精品久久视频播放| 亚洲欧洲精品一区二区精品久久久| 999久久久国产精品视频| 日韩视频一区二区在线观看| 正在播放国产对白刺激| 免费在线观看影片大全网站| 精品人妻在线不人妻| 久久这里只有精品19| 水蜜桃什么品种好| 在线天堂中文资源库| 男女高潮啪啪啪动态图| 免费看a级黄色片| 亚洲黑人精品在线| 99久久99久久久精品蜜桃| 日韩欧美一区二区三区在线观看 | 伦理电影免费视频| 一级毛片精品| 国产精品国产av在线观看| 人人妻人人澡人人看| 又黄又粗又硬又大视频| 亚洲精品一二三| 欧美日韩福利视频一区二区| 久久国产精品影院| 久久中文看片网| 国产成人av教育| 久久精品亚洲精品国产色婷小说| 欧美黄色片欧美黄色片| 精品欧美一区二区三区在线| 日韩免费av在线播放| 最近最新免费中文字幕在线| 亚洲精品久久午夜乱码| 美女午夜性视频免费| 国产欧美日韩综合在线一区二区| 亚洲av片天天在线观看| 99国产精品免费福利视频| 久久久久视频综合| 中文字幕精品免费在线观看视频| 天堂动漫精品| 免费在线观看黄色视频的| 天堂俺去俺来也www色官网| 国产不卡av网站在线观看| 国产欧美日韩综合在线一区二区| 国产精品电影一区二区三区 | 欧美丝袜亚洲另类 | 免费观看精品视频网站| 午夜视频精品福利| 亚洲精华国产精华精| 成人精品一区二区免费| 久久精品aⅴ一区二区三区四区| 精品一区二区三区四区五区乱码| 色老头精品视频在线观看| 国产精品.久久久| 搡老乐熟女国产| 乱人伦中国视频| 最近最新中文字幕大全免费视频| 中文字幕av电影在线播放| 国产高清视频在线播放一区| 日本a在线网址| 91老司机精品| 午夜福利一区二区在线看| 狂野欧美激情性xxxx| 丝瓜视频免费看黄片| 久久99一区二区三区| 夜夜躁狠狠躁天天躁| 女人高潮潮喷娇喘18禁视频| 最近最新中文字幕大全免费视频| 99国产精品99久久久久| 黄片小视频在线播放| 制服诱惑二区| 亚洲专区国产一区二区| 男女床上黄色一级片免费看| svipshipincom国产片| 国产成人系列免费观看| 国产色视频综合| 好看av亚洲va欧美ⅴa在| 久久天堂一区二区三区四区| 99久久精品国产亚洲精品| 国产无遮挡羞羞视频在线观看| 777久久人妻少妇嫩草av网站| 欧美精品高潮呻吟av久久| 亚洲色图 男人天堂 中文字幕| 午夜老司机福利片| 国产高清视频在线播放一区| 国产免费现黄频在线看| 国内久久婷婷六月综合欲色啪| 久久草成人影院| 久久精品国产a三级三级三级| 久久精品亚洲av国产电影网| 麻豆成人av在线观看| 啪啪无遮挡十八禁网站| 欧美乱码精品一区二区三区| 90打野战视频偷拍视频| 在线观看免费高清a一片| 宅男免费午夜| 法律面前人人平等表现在哪些方面| 亚洲少妇的诱惑av| 日本wwww免费看| 男人舔女人的私密视频| 后天国语完整版免费观看| 日本黄色视频三级网站网址 | 欧美日韩瑟瑟在线播放| 久久精品人人爽人人爽视色| 国产激情欧美一区二区| 一二三四社区在线视频社区8| 伦理电影免费视频| 免费少妇av软件| 久久这里只有精品19| 久久人人97超碰香蕉20202| 国产欧美日韩一区二区精品| 身体一侧抽搐| 欧美成狂野欧美在线观看| 热re99久久精品国产66热6| 久久精品国产亚洲av高清一级| av电影中文网址| 91成年电影在线观看| 丝袜美足系列| 久久久久久久午夜电影 | 丁香六月欧美| 亚洲熟妇中文字幕五十中出 | av一本久久久久| 亚洲国产看品久久| 精品欧美一区二区三区在线| 久久国产精品男人的天堂亚洲| 99久久国产精品久久久| 伊人久久大香线蕉亚洲五| 久9热在线精品视频| 亚洲精品中文字幕在线视频| 国产欧美日韩一区二区精品| 国产无遮挡羞羞视频在线观看| 久久久久国内视频| 亚洲三区欧美一区| 亚洲精品自拍成人| 中文亚洲av片在线观看爽 | 成年版毛片免费区| 日韩人妻精品一区2区三区| 变态另类成人亚洲欧美熟女 | 成人国产一区最新在线观看| 美女 人体艺术 gogo| 精品人妻1区二区| 色婷婷av一区二区三区视频| 亚洲三区欧美一区| 欧美色视频一区免费| 大型黄色视频在线免费观看| 精品久久久久久电影网| 天天躁夜夜躁狠狠躁躁| x7x7x7水蜜桃| 免费看a级黄色片| 中文字幕色久视频| 日韩 欧美 亚洲 中文字幕| www.精华液| 欧美老熟妇乱子伦牲交| 黑人猛操日本美女一级片| 日韩欧美一区视频在线观看| av天堂久久9| 亚洲av美国av| 亚洲精品av麻豆狂野| 成人亚洲精品一区在线观看| 99久久99久久久精品蜜桃| 国产精品一区二区在线观看99| 日日爽夜夜爽网站| 国产一区二区激情短视频| 亚洲国产欧美网| 少妇裸体淫交视频免费看高清 | 午夜福利在线观看吧| 亚洲色图av天堂| 久久99一区二区三区| 王馨瑶露胸无遮挡在线观看| 90打野战视频偷拍视频| 欧美人与性动交α欧美软件| 交换朋友夫妻互换小说| 777久久人妻少妇嫩草av网站| 国产精品 国内视频| 一区二区三区激情视频| 国产欧美日韩一区二区精品| 成人国产一区最新在线观看| 99热网站在线观看| 久久狼人影院| 丰满人妻熟妇乱又伦精品不卡| 一区福利在线观看| 丰满迷人的少妇在线观看| 亚洲视频免费观看视频| xxxhd国产人妻xxx| 黄色视频,在线免费观看| 757午夜福利合集在线观看| 亚洲一区二区三区欧美精品| 在线观看免费视频日本深夜| 波多野结衣一区麻豆| 伊人久久大香线蕉亚洲五| 亚洲第一欧美日韩一区二区三区| 国产xxxxx性猛交| 国产人伦9x9x在线观看| 久久久久精品人妻al黑| 国产日韩一区二区三区精品不卡| 亚洲久久久国产精品| 中文字幕av电影在线播放| 国产成人精品在线电影| www日本在线高清视频| 欧美日韩成人在线一区二区| 999久久久国产精品视频| 可以免费在线观看a视频的电影网站| 国产真人三级小视频在线观看| 大片电影免费在线观看免费| 80岁老熟妇乱子伦牲交| 久久草成人影院| 午夜免费观看网址| 日日爽夜夜爽网站| 99精国产麻豆久久婷婷| 免费日韩欧美在线观看| 久久亚洲真实| 每晚都被弄得嗷嗷叫到高潮| 亚洲免费av在线视频| 久久久精品免费免费高清| 成人影院久久| 黄色视频不卡| 亚洲专区字幕在线| 51午夜福利影视在线观看| 亚洲欧美日韩另类电影网站| 后天国语完整版免费观看| 涩涩av久久男人的天堂| 精品一区二区三区四区五区乱码| 精品久久久久久,| 欧洲精品卡2卡3卡4卡5卡区| 99热只有精品国产| av国产精品久久久久影院| 国产成人精品在线电影| 美女国产高潮福利片在线看| 亚洲色图 男人天堂 中文字幕| 老熟女久久久| 久久草成人影院| 精品乱码久久久久久99久播| 亚洲av成人一区二区三| 人成视频在线观看免费观看| 妹子高潮喷水视频| 黄色片一级片一级黄色片| 精品国产一区二区三区四区第35| 国产精品欧美亚洲77777| 国产精品二区激情视频| 91在线观看av| 亚洲精品自拍成人| 丝瓜视频免费看黄片|