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

    Web Resources for Metagenomics Studies

    2015-12-21 03:13:16PravinDudhagaraSunilBhavsarbChintanBhagatcAnjanaGhelaniShreyasBhattRajeshPatel
    Genomics,Proteomics & Bioinformatics 2015年5期

    Pravin Dudhagara*,Sunil Bhavsarb,Chintan Bhagatc,Anjana Ghelani, Shreyas Bhatt,Rajesh Patel

    1Department of Biosciences,Veer Narmad South Gujarat University,Surat 395007,India

    2Department of Microbiology,Shree Ramkrishna Institute of Computer Education and Applied Sciences,Surat 395001,India

    3Department of Life Sciences,Hemchandracharya North Gujarat University,Patan 384265,India

    RESOURCE REVIEW

    Web Resources for Metagenomics Studies

    Pravin Dudhagara1,*,a,Sunil Bhavsar1,b,Chintan Bhagat1,c,Anjana Ghelani2,d, Shreyas Bhatt3,e,Rajesh Patel3,f

    1Department of Biosciences,Veer Narmad South Gujarat University,Surat 395007,India

    2Department of Microbiology,Shree Ramkrishna Institute of Computer Education and Applied Sciences,Surat 395001,India

    3Department of Life Sciences,Hemchandracharya North Gujarat University,Patan 384265,India

    Available online 18 November 2015

    Handled by Vladimir B.Bajic

    Metagenomics; Metagenomes; Web resources; Software tools; Synthetic metagenome

    The developmentofnext-generation sequencing(NGS)platforms spawned an enormous volume of data.This explosion in data has unearthed new scalability challenges for existing bioinformatics tools.The analysis of metagenomic sequences using bioinformatics pipelines is complicated by the substantial complexity of these data.In this article,we review several commonly-used online tools for metagenomics data analysis with respect to their quality and detail of analysis using simulated metagenomics data.There are at least a dozen such software tools presently available in the public domain.Among them,MGRAST,IMG/M,and METAVIR are the most well-known tools according to the number of citations by peer-reviewed scientific media up to mid-2015.Here,we describe 12 online tools with respect to their web link,annotation pipelines,clustering methods,online user support,and availability of data storage.We have also done the rating for each toolto screen more potentialand preferentialtools and evaluated five best tools using synthetic metagenome.The article comprehensively deals with the contemporary problems and the prospects of metagenomics from a bioinformatics viewpoint.

    Introduction

    Metagenomics has become an emerging and indispensable field of microbial ecology since the first use of the concept in 1998 [1].It deals with the study of genetic material recovered directly from environmental samples.Originally,metagenomics studies depended on the cultivation of clonalcultures followed by functional expression screening[1],which are not able to represent the total community profile and may overlook the vast majority of the microbial biodiversity [2].The advent of next-generation sequencing(NGS)technologies has made it possible to perform sequencing-basedmetagenomic analyses.NGS platforms allow researchers to retrieve high-throughput data at substantially lower cost,fueling the development of sequence-based metagenomics analyses aimed at decoding the entire microbial community[3].

    Two distinctmetagenomics approaches are commonly used: the first is referred to as marker-gene metagenomics or targeted metagenomics,and the second is shotgun metagenomics.The former involves the use of marker genes(including conserved markers such as 16S rRNA[4,5],18S rRNA,and ITS[6])and an amplicon sequencing method to reassemble the taxonomic structure of a microbial community[7,8].For the latter approach,the random nature of shotgun sequencing ensures adequate coverage to assess the entire community’s structure [9,10]and opens new avenues for discovering structural and functional novelties[11].At present,thousands of markergene assisted and shotgun metagenomics projects have been undertaken,comprising millions of samples available in the public domain.After sequencing of a metagenomics sample, the primary task is to analyze the vast amount of data in order to determine the identities and roles ofthe microbialspeciespresentin the sample[12].Severaloffline toolsare available forclassifying metagenomic reads against known reference datasets [13].However,the assembly of metagenomic sequence data and the identification of the operational taxonomic units (OTUs)are major technical challenges that persist in metagenomic data analysis pipelines[14].Moreover,the high complexity of metagenomics data is a criticalbarrier in analyzing these metadata using stand-alone tools.On the other hand,various web-based workflows may be able to offer in-depth exploration without the need for highly configured computers.

    In this article,we review the most useful online software tools for metagenomics data analysis,including a brief description of how they work in the course of the analysis. More importantly,we have rated the tools on the basis of overall users’experience in five vital criteria:(1)ease of data uploading,(2)availability of online user support,(3)spectrum of data analyses available,(4)number of citations in scientific literature,and(5)size of the stored data(Figure 1, Table 1).While such an evaluation can never be completely thorough,we nevertheless hope that our evaluation of these web resources for metagenomics studies will provide vital information for budding scientists,assisting them to make the best choice for their preferred applications.

    MG-RAST

    The software package Metagenomics Rapid Annotation using Subsystem Technology(MG-RAST)is a user-friendly opensource server initially constructed on the SEED project (http://www.theseed.org/)[15]framework for metagenomics studies[16–18].Debuting in 2007,MG-RAST was one of the earliest online metagenomics analysis tools[16].A new version (GitHub embedded version 3.6)is currently available,which does not entirely depend on SEED technology but uses the SEED subsystem instead as a preferred data source,to enable taxonomic and functional classification of metagenomes.The new version is encapsulated and separated from the data store, allowing far greater scalability and much faster computation than previous versions.Currently,MG-RAST contains 215,773 metagenomic datasets,which have been accessed by more than 12,000 registered users.Other than data access,registered users can also submit their own raw metagenomic data in FASTA,FASTQ,and SFF formats along with detailed sample information.The uploaded data are processed by a multi-step workflow that includes quality control,automated annotation,and analysis of prokaryotic metagenomic shotgun samples as well as targeted samples.Automated annotation is done using numerous reference datasets.The server supports a variety of user-driven analyses,including phylogenetic,functional,metabolic,and comparative analyses of two or more metagenomes.MG-RAST also provides the facility to select a database for quantitative insights into microbialpopulations based on sequencing data.It also offers the ability to output data in multiple clustering forms that are downloadable in flattext format or to export data as FASTA and QIIME files. Registered users have the right to publish their data openly or keep data private and can share data among multiple users with protected confidentiality.In these ways,MG-RAST offers great flexibility in analysis,privacy,and sharing of data. Therefore,it is the most reviewed and explored tool in use today[8,13,19].So far,more than 2000 metagenome projects and over 13,000 metagenomic datasets have been uploaded to the MG-RAST servers each month of 2015,indicating the great popularity of the tool.

    Figure 1 Relative citation of the metagenomics software tools from articles published in peer-reviewed scientific journals by October 2015 Year in the bracket in the legend box indicates the year of original release of the respective tool.Citation of each tool was tracked from Google Scholar.Totalcitation of alltools is considered as 100%in each year and relative percentage of citations of each toolper year was calculated in relative to the total citations of all tools in the respective year.

    ?

    IMG/M

    The software referred to as the Integrated Microbial Genomes and Metagenomes(IMG/M)is a data storage,management, and analysis system for metagenomes hosted by the Joint Genome Institute(JGI)of U.S.Department of Energy (DOE)[20].IMG/M integrates metagenome datasets with isolated microbialgenomes from the IMG system.It is a domainspecific tool of IMG compatible with sequencing data from microbialcommunities.IMG/Mallows for combined analyses with all available draft and complete genomes,plasmids,and viruses in the public domain[21,22].It provides free support for genomic and metagenomic data annotation,integration, and comparative analyses of integrated genomic and metagenomics data.The data content and analyticaltools are updated continually.As of January 2015,IMG/M contained a total of 5234 metagenome datasets with billions of genes,out of which 3193 metagenome datasets were publicly available.More than 13,945 registered users from 95 countries had gained access as of October 2015.Registered users can access IMG/Mdata content and analytical tools.Registered users can submit data online and maintain data’s privacy for two years before data move to the public domain.Furthermore,analysis and downloading ofmetagenomics datasets via JGI’s portals is also possible for registered users.Pre-processing,quality control, and annotation of input data are performed by JGI’s metagenome annotation system.The stored IMG/Mdata are annotated using multiple reference datasets to conduct analyses in three tiers:(i)phylogenetic composition,(ii)functionalor metabolic potentialwithin individualmicrobiomes,and(iii)comparisons across microbiomes.IMG/M provides support for such analyses by integrating metagenome datasets with isolated microbial genomes from the IMG system[23].IMG/M provides the outputdata in multiple-cluster forms,which facilitates its use as an online tool for environmental and organismal metagenomics analysis.In addition to providing integrated microbial metagenomes,IMG/M has developed the customized IMG/ M-HMP(Human Microbiome Project;http://www.hmpdacc. org),which provides analysis tools for the annotation of HMP-generated metagenome datasets in the reference of all publicly-available genomes and metagenomics samples in IMG[24].IMG/M attracts increasing number of researchers due to its in-depth descriptive analyses of metagenomes and genome extraction through binning of metagenomes.

    METAREP

    Metagenomics Reports(METAREP)is openly available software for comparative metagenomics analysis.It is a well-known open browser for annotation,analysis,and comparison of short reads or assemblies of metagenomics data [25].To upload input data,it requires only basic programming knowledge and the skill to operate the stand-alone version of METAREP.Users can simultaneously link a minimum of two datasets at various functional and taxonomic levels for comparative analyses through applying multiple clustering methods.For each of these clustering methods,METAREP provides a download option to export tab-delimited files for downstream analysis.The software integrates highly-scalable platforms Solr/Lucene,R,and CAKEPHP to compare huge volumes of metagenomic annotations.METAREP has been developed and maintained under utmost privacy by the J.C. Venter Institute(JCVI),and it is hosted at GitHub.Users are encouraged to suggest new features of interest for improving and updating the application.Recently,JCVI has updated the software with a better annotation process,clustering procedure,and input format.These improvements are practically demonstrated by the comparative analysis of shotgun metagenomics data from the HMP(http://www.hmpdacc.org/)[26].

    CoMet

    CoMet is a freely accessible tool for comparative metagenomics analysis.It is the right pipeline for rapid analysis of metagenomic short-read data[27].CoMet can easily handle metagenomics data generated by the 454 and Ion Torrent platforms.For the comparative analysis of metagenomic data,it accepts up to 20 metagenomic DNA sequencing files in multiple FASTA format as input;however zip files may also be submitted for faster uploading.Data sized up to 500 MB can be uploaded at a time for analysis through the easy-to-use submission interface.CoMet is specifically designed to compare the functional or metabolic profiles of multiple metagenomes in terms of Pfam protein domains and the associated Gene Ontology(GO)terms.CoMet combines ORF-finding and subsequent assignment of protein sequences to Pfam domain families with a comparative statistical analysis,providing a quick overview of putative functional differences in a set of metagenomic samples.Registration is not mandatory for job submission.The output comprises downloadable text files and also provides figures and a matrix of the statistical results.It has been the most used tool for comparative metagenomics study since 2011.

    METAGENassist

    METAGENassist is a freely available webserver for comparative metagenomic analysis.It can provide a bacterial census from diverse biomes and performs comprehensive multivariate statisticalanalyses on the data[28].The server accepts the CSV file format along with a wide variety of common formats generated by online and standalone metagenomic analysis tools such as mothur,MG-RAST,QIIME,and MEGAN.It is able to analyze marker-gene assisted as well as shotgun metagenomic data.METAGENassist is the only webserver to perform taxonomic-to-phenotypic mapping of the submitted data.Phenotypic information covering nearly 20 significant functional categories is automatically generated from the input data.The server uses univariate and multivariate dataprocessing and follows a multi-step workflow to cluster data into different layouts.Output is available in annotated colorful tables and graphs,which can be downloaded in PNG or PDF formats.

    MetaABC

    MetaABC is an integrated metagenomics platform for data adjustment,binning,and clustering.Its primary focus is data filtration and normalization techniques for improving taxonomic assignment in metagenomic analysis[29].MetaABC was developed with compatibility for the metagenomic data generated by both NGS and Sanger platforms.Currently, more than 50 metagenomic datasets have been stored in MetaABC.Users can simultaneously upload 2–20 SFF, FASTA,or FASTQ files for analysis in My MetaABC.After the uploading of data files,the user can select and execute all the steps for abundance profiling of microbes.The final outputs are delivered to a user-friendly interface and can be retrieved for downloading in different visualization formats. Although MetaABC has advantages,such as integrating a hierarchical clustering program for comparative analysis of metagenomes,it has limitations as it is not applicable for functional and metabolic profiling of metagenomes.Moreover,its data privacy is limited because unregistered use can also access the data,and sharing feature is not available.Furthermore,it does not accept large(>2 MB)input files.Instead,MetaABC provides a stand-alone version to deal with larger datasets.It has been freely available in the public domain for the last four years,though its citations are low and it does not seem to be gaining enough attention in the scientific field.

    MyTaxa

    MyTaxa is the latest novelplatform for the taxonomic assignment of metagenomic sequences with high accuracy.It has homology-based bioinformatics framework to classify metagenomic and genomic data with extraordinary precision[30]. MyTaxa employs all genes present in an unidentified sequence as classifiers,weighting each gene on the basis of its precomputed classifying power at a specified taxonomic level and frequency of horizontal gene transfer[30].MyTaxa is specifically applicable to consigning and classifying unknown genomic and metagenomic sequences.It consists of a standalone indexed database that is freely accessible and downloadable from MyTaxa’s website.The database contains parameters for the gene clusters that are employed in the online analysis for taxonomic classification.The database is regularly updated to achieve high sensitivity and specificity in taxonomic analyses by incorporating additional reference gene sequences that are available through isolated genome or sequencing projects.The toolis scalable to analyze large dataset as well as compatible with any new,faster algorithms besides homology search algorithms.In the web server,users can start MyTaxa analyses by supplying two files,a standard GFF file and a tabular output file.MyTaxa is equipped with a unique classification system based on the genome-aggregate average amino acid identity(AAI)concept.This classification system substantiates MyTaxa as an appropriate toolfor determining the degree of novelty of unknown metagenomic sequences.However,due to the absence of representative sequences for diverse species in the public domain,taxonomic assignment of the lowest taxa(i.e.,species)is a challenging task for MyTaxa.Moreover,its utility is hampered by the inability to resolve the functional or metabolic profiling of metagenomes.Finally,the BLASTx-translated metagenomic protein file is a prerequisite for operating the analytical tools, which turns out to be a key barrier for the new learners.

    metaMicrobesOnline

    The metaMicrobesOnline software is a free web portal for the phylogenetic study of genes from microbialgenomes and metagenomes.The metaMicrobesOnline analysis pipeline starts with contig assembly,gene calling,and translation into protein sequences,followed by scanning the sequence against recognized gene and protein domain families[31].Outputs of the scanning are used to add genes from the metagenomics data to a multiple sequence alignment for each gene family,followed by construction of phylogenetic trees for each gene/domain family using FastTree-2.Users can select the optimal approach for their data,from assembly and gene calling to gene-tree analysis.Genome-context browsing is available, combining metagenomic contigs and isolated genomes.Therefore,users can determine the contiguously-sequenced genome and search functionally-associated genes.metaMicrobesOnline contains a secured database of 155 metagenomes including 123 ecological and 32 host-associated metagenomes.It has been integrated with the MicrobesOnline(http://www.microbesonline.org)isolated genomes since July 2010,which holds thousands of genome sequences from all three domains of life [32].However,currently this tool is not being used widely and has not been frequently cited due to a lack of regular updates,the limited range of analyses,and restricted submission of new metagenomes.

    EBI Metagenomics

    EBI Metagenomics was the first resource developed in Europe by EMBL-EBI in 2011 for taxonomic,functional,and comparative metagenomics analysis[33].It contains QIIME-embedded tools for data management,analysis,storage,and sharing of metagenomes.Users must register to upload raw nucleotide reads that can remain private up to 2 years upon request,and reads are eventually stored with unique accession numbers in the European Nucleotide Archive(ENA).The uploaded datasets,along with essential details,comply with the Genomic Standards Consortium(GSC)and Minimum Information about any(X)Sequence(MIxS)for validation, so that the unambiguous data can be achieved and reused by the scientific community.Furthermore,users can submit raw sequencing reads generated from any NGS platforms through ENA’s Webin tool or ISAcreator.The EBI Metagenomics analysis pipeline includes trimming and quality checks to remove artifacts and streamline of the data.The processed reads are analyzed by rRNAselector and QIIME to predict taxa,followed by functional assignment by a composite protein database called InterProScan.EBI Metagenomics is well suited for the analysis of shotgun and marker-gene metagenomes.At present,it is mainly focused on shotgunmetagenomes;however,the use of rRNASelector for the extraction of rRNA from shotgun metagenomic datasets is a key feature for analyzing marker-gene metagenomes.The results of analyses are easily accessible via the EBI Metagenomics web interface and are downloadable in a variety of forms compatible with further analyses using online or standalone tools.Rapid comparative taxonomic and functional profiling of metadata is also made possible using EBI Metagenomics.The comparative analysis allows extraction and identification of common and unique information among various metagenomes.Due to the multiple advantages offered by EBI Metagenomics,it has quickly become a favorite tool among beginners.EBI Metagenomics is constantly updating the analysis pipeline with the development of additional analysis and visualization tools.Currently,it endeavors to establish a common collaborative platform for NGS computational setups.In near future,multi-omics data analysis should also be possible through the EBI Metagenomics platform.

    CAMERA

    The Community Cyberinfrastructure for Advanced Marine Microbial Ecology Research and Analysis(CAMERA)was developed in 2007 with the aim of monitoring microbial communities of the ocean and their response to environmental changes[34].In 2010,CAMERA expanded its mission to include all metagenomics data derived from different biomes by dropping the word‘‘marine”from its original name[35]. Currently,the Data Distribution Center(DDC)of CAMERA is no longer available at the original domain(http:// camera.calit2.net),but has been relocated to the iMicrobe Project(http://data.imicrobe.us),which holds 128 projects and 2660 samples as of October 2015.Like MG-RAST and IMG/M,CAMERA is a prominent large-scale database that processes,shares,and secures metagenomic datasets.The sharing of data and computed results is currently under the auspices of the GSC.CAMERA permitted the dataset’s publication and was the first web tool to support the GSC’s Minimal Information checklists for metadata.To submit data, users are required to register in the iPlant cyberinfrastructure through the Discovery Environment(DE)web interface. Registered users normally are allotted 100 GB of virtualspace, and an additional allocation of up to 1 TB can be granted upon request.The uploading and importing of huge datasets (>1.9 GB)within the DE is carried out either by Cyberduck for Mac and Windows users or by iDrop Desktop(iCommands)for LINUX users.Alternatively,small data files can be uploaded and imported within the DE using a simple URL.CAMERA’s fundamental task is to create a rich,distinctive data repository and bioinformatics tool resource that overcomes the unique challenges of metagenomics such as complexity,heterogeneity,and fragmented data.The database includes environmental metagenomic and genomic sequence data,associated environmental parameters,pre-computed search results,and software tools to support powerful crossanalysis of environmental samples.It holds a large number of samples from across the world;however,since 2011,its use and resulting citations have consistently decreased because ofthe technically lengthy and complicated data uploading process.The descriptive analysis of metagenomes is also not userfriendly.However:recently,it is incorporated into QIIME and offers rapid online cloud computing service with greater user involvement and tractability during the analysis process.

    METAVIR

    METAVIR is a freely available webserver exclusively designed to annotate the raw and assembled viral metagenomic sequences[36].Recently,a new framework,METAVIR 2, has been released that allows for comprehensive virome analysis[37].The server is useful for exploring viraldiversity from metagenomes and undertakes comparative analyses through pre-computed phylogenies of viral marker genes,rarefaction curves,and multivariate analyses constructed on k-mer signatures and BLAST-based comparisons.METAVIR 2 can handle assembled viromes made up of thousands of large contigs. As of October 2015,METAVIR holds 335 virome projects that comprise more than 64 million sequences from diverse samples.Registered users can upload their sequence datasets only in FASTA file in a private space.However,users can also submit large data with compressed files in zip,gzip,or tar.gz format.The output can be visualized in different interactive and dynamic forms,i.e.,tables,phylogenetic trees,recruitment plots,and maps.METAVIR has an easy-to-use submission interface and offers rapid in-depth virome comparison along with data privacy,which are key features that make it the most cited virome analysis pipeline.

    VIROME

    Viral Informatics Resource for Metagenome Exploration (VIROME)is a web-based application specifically designed for the analysisofmetagenomic data collected from viralassemblages occurring within different environmental contexts.It emphasizes the classification of metagenome viral sequences and putative ORFs based on homology search results against both known and environmental sequences[38].Submission of sequence data to VIROME is initiated through the web-application interface.After filling out a form,the user is contacted by Email correspondence with further instructions regarding sequence data transfer.It accepts multiple file formats,including FASTA,QUAL,FASTQ,and SFF. Metagenome libraries submitted for analysis at VIROME are blasted against the MetagenomesOnline database(MgOl),with a graphical interface and links to various protein metadata enabling exploration of the results.The VIROME workflow consists of two consecutive steps,sequence quality screening and sequence analysis against the UniVec database,followed by three parallel steps that include ORF-calling using MetaGeneAnnotator,known protein identification,and environmental protein characterization by subject databases such as MGOL and UniRef 100.Sequence quality trimming is run using a workflow managementsystem called Ergatis.Data from the sequence processing and BLAST analysis components are stored in a back-end MySQL database.An Adobe ColdFusion server handles communication between the MySQL database and the VIROME web application.VIROMEoffers downloadable outputthatcan be used as inputdata for multivariate analyses using third-party tools.The main hurdle to its widespread use is the length of time needed to analyze even a small-sized metagenome,usually severalweeks to months.

    Table 2 Result of synthetic metagenome analysis using the five selected tools

    Analysis of synthetic metagenome

    Synthetic metagenome was analyzed using five selected tools to demonstrate the efficiency and analysis time course.Synthetic metagenome(Run accession:ERR340319)was downloaded from EMBL-EBI.It contains 32,782,167 bp incorporated in 90,385 raw sequences generated from 454 GS FLX Titanium sequencing.This synthetic metagenome was uploaded into MG-RAST,IMG/M,EBI Metagenomics,CoMetand METAVIR that were selected based on rating criteria for taxonomic assessment using online support provided by the software developers.The result showed the numbers of different microbial families reported from each tool are not exactly similar due to the different annotation pipelines used by these tools. Moreover,the analysis time required by the each toolwas also varied,because analysis depends on the number of running job in server and length of input sequences(Table 2).Though the remaining tools provided online help for analysis,it was not feasible to complete the data submission and taxonomic analysis of such a large synthetic metagenome within a short period of time by beginning users.

    Concluding remarks

    Structural and functional microbial profiling has rapidly become commonplace,thanks to the advent of NGS platforms.Additionally,the comparative metagenomics of various biomes quickly enlarged the field of microbial ecology.To decode the ecosystem comprehensively,new tools,frameworks,and hypotheses will be required for analysis,storage, visualization,and sharing of the dataset.A single platform is therefore not adequate for conducting holistic metagenomics analyses.Longer-read sequences,accurate assembly,and annotation pipelines are anticipated developments from impending metagenomics study.The present review has given an overview of the prevailing metagenomics analysis tools,and we expect more advanced software tools from their original developers in the near future.This review may serve as an atlas for rising researchers to choose the most appropriate webbased platforms for metagenomic data analysis.It is specifically intended to aid those organizations that have limited computational resources with which to analyze their data.

    Competing interests

    The authors have declared no competing interests.

    Acknowledgments

    The authors are thankful to BIT Virtual Centre,Patan Node, Hemchandracharya North Gujarat University,India for providing servers and computational facilities.

    References

    [1]Handelsman J,Rondon MR,Brady SF,Clardy J,Goodman RM. Molecular biological access to the chemistry of unknown soil microbes:a new frontier for natural products.Chem Biol 1998;5: R245–9.

    [2]Daniel R.The metagenomics of soil.Nat Rev Microbiol 2005;3:470–8.

    [3]Thomas T,Gilbert J,Meyer F.Metagenomics–a guide from sampling to data analysis.Microb Inform Exp 2012;2:3.

    [4]Dudhagara P,Ghelani A,Patel R,Chaudhari R,Bhatt S. Bacterial tag encoded FLX titanium amplicon pyrosequencing (bTEFAP)based assessment of prokaryotic diversity in metagenome of Lonar soda lake,India.Genom Data 2015;4:8–11.

    [5]Ghelani A,Patel R,Mangrola A,Dudhagara P.Cultivationindependent comprehensive survey of bacterial diversity in Tulsi Shyam Hot Springs,India.Genom Data 2015;4:54–6.

    [6]Dudhagara P,Ghelani A,Bhavsar S,Bhatt S.Metagenomic data of fungal internal transcribed Spacer and 18S rRNA gene sequences from Lonar lake sediment,India.Data Brief 2015;4:266–8.

    [7]Carlos N,Tang YW,Pei Z.Pearls and pitfalls of genomics-based microbiome analysis.Emerg Microbes Infect 2012;1:e45.

    [8]Oulas A,Pavloudi C,Polymenakou P,Pavlopoulos GA, Papanikolaou N,Kotoulas G,et al.Metagenomics:tools and insights for analyzing next-generation sequencing data derived from biodiversity studies.Bioinform Biol Insights 2015;5:75–88.

    [9]Patel R,Mevada V,Prajapati D,Dudhagara P,Koringa P,Joshi CG.Metagenomic sequence of saline desert microbiota from wild ass sanctuary,Little Rann of Kutch,Gujarat,India.Genom Data 2015;3:137–9.

    [10]Mangrola AV,Dudhagara P,Koringa P,Joshi CG,Patel RK. Shotgun metagenomic sequencing based microbial diversity assessment of Lasundra hot spring,India.Genom Data 2015;4:73–5.

    [11]Singh AH,Doerks T,Letunic I,Raes J,Bork P.Discovering functionalnovelty in metagenomes:examples from light-mediated processes.J Bacteriol2009;191:32–41.

    [12]Ounit R,Wanamaker S,Close TJ,Lonardi S.CLARK:fast and accurate classification of metagenomic and genomic sequences using discriminative k-mers.BMC Genomics 2015;16:236.

    [13]Mandal RS,Saha S,Das S.Metagenomic surveys of gut microbiota.Genomics Proteomics Bioinformatics 2015;13:148–58.

    [14]Behnam E,Smith AD.The Amordad database engine for metagenomics.Bioinformatics 2014;30:2949–55.

    [15]Overbeek R,Begley T,Butler RM,Choudhuri JV,Chuang HY, Cohoon M,etal.The subsystems approach to genome annotation and its use in the projectto annotate 1000 genomes.Nucleic Acids Res 2005;33:5691–702.

    [16]Meyer F,Paarmann D,D’Souza M,Olson R,Glass EM,Kubal M,et al.The metagenomics RAST server–a public resource for the automatic phylogenetic and functional analysis of metagenomes.BMC Bioinformatics 2008;9:386.

    [17]Overbeek R,Olson R,Pusch GD,Olsen GJ,Davis JJ,Disz T, et al.The SEEDand the Rapid Annotation of microbialgenomes using Subsystems Technology(RAST).Nucleic Acids Res 2014;42:D206–14.

    [18]Wilke A,Bischof J,Harrison T,Brettin T,D’Souza M,Gerlach W,et al.A RESTful API for accessing microbialcommunity data for MG-RAST.PLoS Comput Biol 2015;11:e1004008.

    [19]Sun Q,Liu L,Wu L,Li W,Liu Q,Zhang J,et al.Web resources for microbial data.Genomics Proteomics Bioinformatics 2015;13:69–72.

    [20]Markowitz VM,Ivanova NN,Szeto E,Palaniappan K,Chu K, Dalevi D,et al.IMG/M:a data management and analysis system for metagenomes.Nucleic Acids Res 2008;36:D534–8.

    [21]Markowitz VM,Chen IM,Chu K,Szeto E,Palaniappan K, Grechkin Y,et al.IMG/M:the integrated metagenome data management and comparative analysis system.Nucleic Acids Res 2012;40:D123–9.

    [22]Markowitz VM,Chen IM,Chu K,Szeto E,Palaniappan K,Pillay M,et al.IMG/M 4 version of the integrated metagenome comparative analysis system.Nucleic Acids Res 2014;42:D568–73.

    [23]Markowitz VM,Chen IM,Palaniappan K,Chu K,Szeto E, Grechkin Y,et al.IMG:the Integrated Microbial Genomes database and comparative analysis system.Nucleic Acids Res 2012;40:D115–22.

    [24]Markowitz VM,Chen IM,Chu K,Szeto E,Palaniappan K,Jacob B,et al.IMG/M-HMP:a metagenome comparative analysis system for the human microbiome project.PLoS One 2012;7: e40151.

    [25]Goll J,Rusch DB,Tanenbaum DM,Thiagarajan M,LiK,Methe BA,et al.METAREP:JCVI metagenomics reports–an open source tool for high-performance comparative metagenomics. Bioinformatics 2010;26:2631–2.

    [26]GollJ,Thiagarajan M,Abubucker S,Huttenhower C,Yooseph S, Methe′BA.A case study for large-scale human microbiome analysis using JCVI’s metagenomics reports(METAREP).PLoS One 2012;7:e29044.

    [27]Lingner T,Asshauer KP,Schreiber F,Meinicke P.CoMet–a web server for comparative functional profiling of metagenomes. Nucleic Acids Res 2011;39:W518–23.

    [28]Arndt D,Xia J,Liu Y,Zhou Y,Guo AC,Cruz JA,et al. METAGENassist:a comprehensive web server for comparative metagenomics.Nucleic Acids Res 2012;40:W88–95.

    [29]Su CH,Hsu MT,Wang TY,Chiang S,Cheng JH,Weng FC,et al. MetaABC–an integrated metagenomics platform for data adjustment,binning and clustering.Bioinformatics 2011;27:2298–9.

    [30]Luo C,Rodriguez-R LM,Konstantinidis KT.MyTaxa:an advanced taxonomic classifier for genomic and metagenomic sequences.Nucleic Acids Res 2014;42:e73.

    [31]Chivian D,Dehal PS,Keller K,Arkin AP.MetaMicrobesOnline: phylogenomic analysis of microbial communities.Nucleic Acids Res 2013;41:D648–54.

    [32]Dehal PS,Joachimiak MP,Price MN,Bates JT,Baumohl JK, Chivian D,et al.MicrobesOnline:an integrated portal for comparative and functional genomics.Nucleic Acids Res 2010;38:D396–400.

    [33]Hunter S,Corbett M,Denise H,Fraser M,Gonzalez-Beltran A, Hunter C,et al.EBI metagenomics—a new resource for the analysis and archiving of metagenomic data.Nucleic Acids Res 2014;42:D600–6.

    [34]SeshadriR,Kravitz SA,Smarr L,Gilna P,Frazier M.CAMERA: a community resource for metagenomics.PLoS Biol 2007;5:e75.

    [35]Sun S,Chen J,Li W,Altintas I,Lin A,Peltier S,et al.Community cyberinfrastructure for Advanced Microbial Ecology Research and Analysis:the CAMERA resource.Nucleic Acids Res 2011;39: D546–51.

    [36]Roux S,Faubladier M,Mahul A,Paulhe N,Bernard A,Debroas D,et al.METAVIR:a web server dedicated to virome analysis. Bioinformatics 2011;27:3074–5.

    [37]Roux S,Tournayre J,Mahul A,Debroas D,Enault F. METAVIR 2:new tools for viral metagenome comparison and assembled virome analysis.BMC Bioinformatics 2014;15:76.

    [38]Wommack KE,Bhavsar J,Polson SW,Chen J,Dumas M, Srinivasiah S,et al.VIROME:a standard operating procedure for analysis of viral metagenome sequences.Stand Genomic Sci 2012;6:421–33.

    6 August 2015;revised 9 October 2015;accepted 13 October 2015

    .

    E-mail:dudhagarapr@gmail.com(Dudhagara P).aORCID:0000-0002-4041-6751.bORCID:0000-0002-7525-0624.cORCID:0000-0002-2318-5998.dORCID:0000-0003-4638-6925.eORCID:0000-0002-4364-2244.fORCID:0000-0002-4712-1921.

    Peer review under responsibility of Beijing Institute of Genomics, Chinese Academy of Sciences and Genetics Society of China.

    18美女黄网站色大片免费观看| 首页视频小说图片口味搜索| 在线观看免费日韩欧美大片| 麻豆av在线久日| 18禁裸乳无遮挡免费网站照片 | 麻豆国产av国片精品| 可以在线观看的亚洲视频| 欧美色欧美亚洲另类二区| 国产日本99.免费观看| 欧美另类亚洲清纯唯美| 亚洲国产日韩欧美精品在线观看 | 免费观看人在逋| 嫩草影院精品99| 久久中文字幕一级| 十分钟在线观看高清视频www| 黄色毛片三级朝国网站| 久久欧美精品欧美久久欧美| 亚洲无线在线观看| 女警被强在线播放| 99久久99久久久精品蜜桃| 久久婷婷成人综合色麻豆| 午夜免费观看网址| 大香蕉久久成人网| 久久久久久久久免费视频了| 少妇粗大呻吟视频| 日韩三级视频一区二区三区| 真人一进一出gif抽搐免费| 久久人妻av系列| 成人永久免费在线观看视频| 巨乳人妻的诱惑在线观看| 少妇裸体淫交视频免费看高清 | 国产久久久一区二区三区| 成人永久免费在线观看视频| 久久狼人影院| 18美女黄网站色大片免费观看| 国产精品久久电影中文字幕| 欧美日韩瑟瑟在线播放| 两个人免费观看高清视频| 最近最新中文字幕大全免费视频| 日韩欧美在线二视频| 女性被躁到高潮视频| 亚洲欧美日韩无卡精品| 人成视频在线观看免费观看| 亚洲天堂国产精品一区在线| 老汉色∧v一级毛片| 亚洲精品色激情综合| 18禁裸乳无遮挡免费网站照片 | 啦啦啦免费观看视频1| 日韩视频一区二区在线观看| 免费看日本二区| 国产成+人综合+亚洲专区| av中文乱码字幕在线| 亚洲狠狠婷婷综合久久图片| 亚洲欧美激情综合另类| 亚洲最大成人中文| 久久精品国产亚洲av香蕉五月| 色播亚洲综合网| 欧美成人免费av一区二区三区| 欧美乱色亚洲激情| 美女大奶头视频| 性色av乱码一区二区三区2| 欧美黄色淫秽网站| 久久精品亚洲精品国产色婷小说| 午夜a级毛片| 久久久久久久精品吃奶| 成人三级做爰电影| 久热爱精品视频在线9| 婷婷六月久久综合丁香| 日韩大码丰满熟妇| 精品一区二区三区视频在线观看免费| 波多野结衣高清作品| 国产精华一区二区三区| 99国产精品99久久久久| 村上凉子中文字幕在线| 亚洲国产欧美一区二区综合| 国内久久婷婷六月综合欲色啪| av免费在线观看网站| 在线永久观看黄色视频| 欧美色视频一区免费| 日韩大尺度精品在线看网址| 亚洲精品中文字幕在线视频| 午夜福利免费观看在线| 桃色一区二区三区在线观看| 国产精品国产高清国产av| 日韩 欧美 亚洲 中文字幕| 亚洲av熟女| 美女扒开内裤让男人捅视频| 精品一区二区三区视频在线观看免费| 欧美日韩中文字幕国产精品一区二区三区| 国产精品亚洲av一区麻豆| 欧美绝顶高潮抽搐喷水| 久久99热这里只有精品18| 一级片免费观看大全| www.熟女人妻精品国产| 日韩欧美 国产精品| 美女高潮到喷水免费观看| 免费在线观看黄色视频的| 中文字幕av电影在线播放| 又黄又爽又免费观看的视频| 我的亚洲天堂| 国产一级毛片七仙女欲春2 | 国产免费男女视频| 日本 av在线| www.www免费av| 男女那种视频在线观看| 男人舔女人的私密视频| 日韩欧美免费精品| 国产激情久久老熟女| 变态另类成人亚洲欧美熟女| 哪里可以看免费的av片| 久久久久久亚洲精品国产蜜桃av| 日韩国内少妇激情av| 一区福利在线观看| 真人一进一出gif抽搐免费| 久久 成人 亚洲| 欧美又色又爽又黄视频| 色播亚洲综合网| 最近最新中文字幕大全电影3 | 十八禁人妻一区二区| 我的亚洲天堂| 欧美av亚洲av综合av国产av| 18禁裸乳无遮挡免费网站照片 | 久久99热这里只有精品18| 给我免费播放毛片高清在线观看| 欧美日本视频| 国产av一区二区精品久久| 日韩高清综合在线| 色播在线永久视频| 国产精品香港三级国产av潘金莲| 一本久久中文字幕| netflix在线观看网站| 亚洲国产精品999在线| 一个人观看的视频www高清免费观看 | 久久久久亚洲av毛片大全| 国产伦人伦偷精品视频| 精品第一国产精品| 亚洲 国产 在线| 丝袜美腿诱惑在线| 成人18禁高潮啪啪吃奶动态图| 久久精品91蜜桃| 国产精品九九99| 97超级碰碰碰精品色视频在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 深夜精品福利| 最新在线观看一区二区三区| 欧美精品亚洲一区二区| 黄色视频不卡| 国产成人啪精品午夜网站| 黄网站色视频无遮挡免费观看| АⅤ资源中文在线天堂| 国产97色在线日韩免费| 亚洲三区欧美一区| 国产精品自产拍在线观看55亚洲| 欧美成人免费av一区二区三区| 日本成人三级电影网站| 99国产综合亚洲精品| 亚洲精品国产一区二区精华液| 亚洲中文字幕日韩| 丰满的人妻完整版| 日本免费a在线| 一本久久中文字幕| 久久热在线av| 一边摸一边做爽爽视频免费| 国产精品久久久人人做人人爽| 欧美日韩中文字幕国产精品一区二区三区| av福利片在线| 91老司机精品| 国产午夜福利久久久久久| 巨乳人妻的诱惑在线观看| 视频区欧美日本亚洲| 99riav亚洲国产免费| 欧美黄色淫秽网站| 9191精品国产免费久久| 人人妻人人澡欧美一区二区| 极品教师在线免费播放| aaaaa片日本免费| 高潮久久久久久久久久久不卡| 色老头精品视频在线观看| 亚洲专区中文字幕在线| 色综合站精品国产| 可以在线观看的亚洲视频| 久久热在线av| 亚洲中文字幕一区二区三区有码在线看 | 真人一进一出gif抽搐免费| 美女扒开内裤让男人捅视频| 又黄又粗又硬又大视频| 亚洲国产毛片av蜜桃av| 一夜夜www| 免费在线观看亚洲国产| 黄网站色视频无遮挡免费观看| or卡值多少钱| 亚洲国产中文字幕在线视频| 波多野结衣高清无吗| 国产成人精品无人区| 精品久久久久久久久久免费视频| av在线播放免费不卡| 午夜精品久久久久久毛片777| 美国免费a级毛片| 国产精品久久久久久精品电影 | 中文字幕人成人乱码亚洲影| 亚洲欧美日韩无卡精品| 亚洲精品久久国产高清桃花| 国产成人精品久久二区二区免费| 久久精品aⅴ一区二区三区四区| 欧美日韩福利视频一区二区| 麻豆成人av在线观看| 亚洲av日韩精品久久久久久密| 亚洲精品国产一区二区精华液| 12—13女人毛片做爰片一| 亚洲av第一区精品v没综合| 日本一本二区三区精品| 精华霜和精华液先用哪个| 制服诱惑二区| 不卡av一区二区三区| 欧美黄色片欧美黄色片| 国产亚洲欧美精品永久| 999精品在线视频| 国产精品精品国产色婷婷| 亚洲男人的天堂狠狠| 久久久久久人人人人人| cao死你这个sao货| 美国免费a级毛片| 免费在线观看完整版高清| e午夜精品久久久久久久| 一二三四在线观看免费中文在| 国产亚洲欧美在线一区二区| 国产久久久一区二区三区| www.自偷自拍.com| 丝袜人妻中文字幕| 少妇熟女aⅴ在线视频| 免费高清视频大片| √禁漫天堂资源中文www| 日韩免费av在线播放| 母亲3免费完整高清在线观看| 亚洲在线自拍视频| 婷婷六月久久综合丁香| 成年女人毛片免费观看观看9| 老司机靠b影院| 51午夜福利影视在线观看| 黄色a级毛片大全视频| 亚洲国产精品sss在线观看| 国产av又大| 国产精品国产高清国产av| 琪琪午夜伦伦电影理论片6080| 国产亚洲欧美精品永久| 久久久国产欧美日韩av| 久99久视频精品免费| 婷婷精品国产亚洲av| 90打野战视频偷拍视频| 人妻丰满熟妇av一区二区三区| 亚洲国产毛片av蜜桃av| 久久久国产成人免费| 国产黄片美女视频| 一夜夜www| 亚洲精品av麻豆狂野| 久久99热这里只有精品18| 桃色一区二区三区在线观看| 久久精品人妻少妇| 免费看十八禁软件| 国产免费男女视频| 19禁男女啪啪无遮挡网站| svipshipincom国产片| av中文乱码字幕在线| 青草久久国产| 侵犯人妻中文字幕一二三四区| 免费看美女性在线毛片视频| 在线国产一区二区在线| 国产一区二区在线av高清观看| www.精华液| 久久久久亚洲av毛片大全| 欧美成狂野欧美在线观看| 亚洲国产精品合色在线| 免费女性裸体啪啪无遮挡网站| 亚洲欧美日韩高清在线视频| 精品国产美女av久久久久小说| 国产区一区二久久| 伊人久久大香线蕉亚洲五| 一区二区三区精品91| 国产精品国产高清国产av| 国产一区二区三区在线臀色熟女| 国产在线精品亚洲第一网站| 国产又色又爽无遮挡免费看| www国产在线视频色| 午夜福利高清视频| 欧美色视频一区免费| 在线观看66精品国产| 黑丝袜美女国产一区| 可以免费在线观看a视频的电影网站| 天天一区二区日本电影三级| 在线免费观看的www视频| 麻豆国产av国片精品| 男男h啪啪无遮挡| 国产三级在线视频| 午夜福利在线观看吧| 制服人妻中文乱码| 欧美+亚洲+日韩+国产| 久久中文字幕人妻熟女| 日本三级黄在线观看| 日本成人三级电影网站| 欧美性猛交黑人性爽| 亚洲成国产人片在线观看| 欧美日韩精品网址| 国产三级黄色录像| 校园春色视频在线观看| 12—13女人毛片做爰片一| 久久久精品欧美日韩精品| 精品日产1卡2卡| 婷婷精品国产亚洲av| www.999成人在线观看| 午夜老司机福利片| 在线观看66精品国产| 久久精品影院6| 老司机午夜福利在线观看视频| 亚洲av成人不卡在线观看播放网| 亚洲成人久久爱视频| 欧美亚洲日本最大视频资源| 99在线视频只有这里精品首页| 精品少妇一区二区三区视频日本电影| a级毛片在线看网站| 亚洲最大成人中文| 成人午夜高清在线视频 | a级毛片在线看网站| 757午夜福利合集在线观看| 久久人妻福利社区极品人妻图片| 亚洲一区二区三区不卡视频| 嫩草影院精品99| 一级毛片高清免费大全| 久久久久久国产a免费观看| 高潮久久久久久久久久久不卡| 18禁国产床啪视频网站| 亚洲天堂国产精品一区在线| 久久中文看片网| 国产高清videossex| avwww免费| 亚洲欧美精品综合久久99| 国产亚洲精品第一综合不卡| 久久久久久亚洲精品国产蜜桃av| 在线视频色国产色| 精品欧美一区二区三区在线| 国产真人三级小视频在线观看| 国产日本99.免费观看| www.自偷自拍.com| 在线看三级毛片| 亚洲av成人不卡在线观看播放网| 国产日本99.免费观看| 韩国精品一区二区三区| 美女 人体艺术 gogo| 悠悠久久av| 真人一进一出gif抽搐免费| 精品国产超薄肉色丝袜足j| 高清毛片免费观看视频网站| 一区二区三区激情视频| 后天国语完整版免费观看| 妹子高潮喷水视频| tocl精华| 久久久久久九九精品二区国产 | 国产精品野战在线观看| 视频区欧美日本亚洲| 亚洲av中文字字幕乱码综合 | 色老头精品视频在线观看| 亚洲七黄色美女视频| 欧美绝顶高潮抽搐喷水| 69av精品久久久久久| 久久久久国产精品人妻aⅴ院| 性色av乱码一区二区三区2| 熟妇人妻久久中文字幕3abv| 大型av网站在线播放| 欧美精品啪啪一区二区三区| 免费搜索国产男女视频| av免费在线观看网站| 久久天躁狠狠躁夜夜2o2o| 成年版毛片免费区| 国产日本99.免费观看| 日本免费一区二区三区高清不卡| 热99re8久久精品国产| 自线自在国产av| 免费在线观看黄色视频的| 18禁观看日本| 桃色一区二区三区在线观看| 啦啦啦免费观看视频1| 成熟少妇高潮喷水视频| 国产高清激情床上av| 国产精品永久免费网站| 亚洲av美国av| 可以在线观看毛片的网站| 日韩三级视频一区二区三区| 亚洲国产欧美网| 露出奶头的视频| 亚洲午夜理论影院| 国内毛片毛片毛片毛片毛片| 美女国产高潮福利片在线看| 国产精品久久视频播放| 妹子高潮喷水视频| 18禁黄网站禁片午夜丰满| 97碰自拍视频| 色尼玛亚洲综合影院| 午夜免费成人在线视频| 国产精品亚洲av一区麻豆| 精华霜和精华液先用哪个| 亚洲avbb在线观看| 亚洲第一青青草原| 免费观看精品视频网站| 久久久久久大精品| 欧美 亚洲 国产 日韩一| 久久精品国产亚洲av高清一级| 色老头精品视频在线观看| 欧美黑人巨大hd| 国产精品亚洲一级av第二区| 亚洲一卡2卡3卡4卡5卡精品中文| 韩国精品一区二区三区| www国产在线视频色| 国产不卡一卡二| 久久伊人香网站| 精品国产美女av久久久久小说| xxx96com| 欧美性猛交黑人性爽| 搡老熟女国产l中国老女人| 18禁美女被吸乳视频| 欧美激情久久久久久爽电影| 国产成人影院久久av| 妹子高潮喷水视频| 三级毛片av免费| 久久人人精品亚洲av| 99久久综合精品五月天人人| 香蕉av资源在线| 男女那种视频在线观看| 视频区欧美日本亚洲| 香蕉久久夜色| av超薄肉色丝袜交足视频| 夜夜看夜夜爽夜夜摸| 99在线视频只有这里精品首页| 好男人电影高清在线观看| 国产区一区二久久| av在线播放免费不卡| 日韩欧美国产在线观看| 日韩三级视频一区二区三区| 国产精品亚洲美女久久久| 亚洲一码二码三码区别大吗| 亚洲专区国产一区二区| 麻豆成人午夜福利视频| 精品人妻1区二区| 精品欧美一区二区三区在线| 91麻豆精品激情在线观看国产| 国产精品爽爽va在线观看网站 | 婷婷六月久久综合丁香| 精品第一国产精品| 露出奶头的视频| 超碰成人久久| 一本久久中文字幕| 男人舔奶头视频| 老司机午夜十八禁免费视频| 国产不卡一卡二| 成人18禁高潮啪啪吃奶动态图| 在线天堂中文资源库| 中文资源天堂在线| 免费人成视频x8x8入口观看| 精品一区二区三区视频在线观看免费| 99热6这里只有精品| 欧美在线黄色| 午夜久久久久精精品| 国产av不卡久久| 91九色精品人成在线观看| 亚洲精品中文字幕一二三四区| 在线观看66精品国产| 欧美人与性动交α欧美精品济南到| 女人高潮潮喷娇喘18禁视频| 日本在线视频免费播放| 19禁男女啪啪无遮挡网站| 好看av亚洲va欧美ⅴa在| 一个人观看的视频www高清免费观看 | 亚洲欧美激情综合另类| 久久精品夜夜夜夜夜久久蜜豆 | 一本综合久久免费| 18美女黄网站色大片免费观看| 亚洲国产欧美日韩在线播放| 亚洲一区高清亚洲精品| 黄色女人牲交| 欧美最黄视频在线播放免费| 精品免费久久久久久久清纯| 欧美成人午夜精品| 国产成人av教育| 中文字幕人成人乱码亚洲影| 亚洲欧美日韩高清在线视频| 少妇裸体淫交视频免费看高清 | 久久欧美精品欧美久久欧美| www.精华液| 久久精品影院6| 成人永久免费在线观看视频| 俺也久久电影网| 日韩精品青青久久久久久| 久久草成人影院| 1024香蕉在线观看| 精品卡一卡二卡四卡免费| 满18在线观看网站| 老汉色∧v一级毛片| 欧美丝袜亚洲另类 | 天天躁狠狠躁夜夜躁狠狠躁| 中文在线观看免费www的网站 | 夜夜夜夜夜久久久久| 男女下面进入的视频免费午夜 | 他把我摸到了高潮在线观看| 男男h啪啪无遮挡| 欧美在线一区亚洲| 丝袜在线中文字幕| 精品午夜福利视频在线观看一区| 亚洲久久久国产精品| 国产亚洲精品一区二区www| 中国美女看黄片| 欧美+亚洲+日韩+国产| 亚洲avbb在线观看| 国产精品九九99| 国产精品国产高清国产av| 麻豆国产av国片精品| 搡老岳熟女国产| 国产精品一区二区免费欧美| 免费观看人在逋| 日本黄色视频三级网站网址| 十分钟在线观看高清视频www| 久久精品夜夜夜夜夜久久蜜豆 | 欧美不卡视频在线免费观看 | 热re99久久国产66热| 久久亚洲真实| 国产爱豆传媒在线观看 | 一二三四在线观看免费中文在| 免费在线观看黄色视频的| av在线播放免费不卡| 国产精品电影一区二区三区| 日韩av在线大香蕉| 91字幕亚洲| 国产精品久久电影中文字幕| 欧美激情极品国产一区二区三区| 手机成人av网站| 麻豆成人午夜福利视频| 男人舔女人下体高潮全视频| 久久久久精品国产欧美久久久| 91麻豆精品激情在线观看国产| 亚洲人成网站高清观看| 亚洲精品国产精品久久久不卡| 国产激情偷乱视频一区二区| av在线天堂中文字幕| 不卡av一区二区三区| 亚洲aⅴ乱码一区二区在线播放 | 亚洲三区欧美一区| 国产熟女午夜一区二区三区| 国内少妇人妻偷人精品xxx网站 | 一区二区日韩欧美中文字幕| www.熟女人妻精品国产| 99热只有精品国产| 国产成人av激情在线播放| 国产黄a三级三级三级人| 一区二区三区国产精品乱码| 久久久精品欧美日韩精品| 无限看片的www在线观看| 久久精品国产亚洲av高清一级| 婷婷亚洲欧美| 国产精品一区二区免费欧美| 黄片小视频在线播放| 免费在线观看日本一区| 亚洲av日韩精品久久久久久密| 91九色精品人成在线观看| 久久精品国产亚洲av高清一级| 国产乱人伦免费视频| 男人的好看免费观看在线视频 | 欧美av亚洲av综合av国产av| 亚洲avbb在线观看| 久久婷婷人人爽人人干人人爱| 丰满人妻熟妇乱又伦精品不卡| 国产精品免费一区二区三区在线| 日本一本二区三区精品| 国产精品日韩av在线免费观看| 亚洲国产精品合色在线| 18禁黄网站禁片午夜丰满| 在线永久观看黄色视频| 国产成+人综合+亚洲专区| 欧美色视频一区免费| 日韩精品免费视频一区二区三区| 久久久久久久久中文| 嫁个100分男人电影在线观看| 美女大奶头视频| 亚洲精品av麻豆狂野| 中亚洲国语对白在线视频| 亚洲真实伦在线观看| 国产精品一区二区免费欧美| 亚洲av熟女| 老司机靠b影院| √禁漫天堂资源中文www| 黄色视频,在线免费观看| 欧美日韩黄片免| 精品久久久久久成人av| 国产三级黄色录像| 在线观看免费视频日本深夜| www国产在线视频色| 国产一区二区激情短视频| 在线观看免费视频日本深夜| 可以免费在线观看a视频的电影网站| 国产黄a三级三级三级人| 变态另类成人亚洲欧美熟女| 亚洲成国产人片在线观看| 19禁男女啪啪无遮挡网站| 又紧又爽又黄一区二区| 成人国语在线视频| 啦啦啦 在线观看视频| 精品不卡国产一区二区三区| www国产在线视频色| 琪琪午夜伦伦电影理论片6080| 亚洲成人免费电影在线观看| 老司机福利观看| 91成年电影在线观看| 动漫黄色视频在线观看| 国产激情久久老熟女| 国产精品av久久久久免费| 黄色 视频免费看| 高清在线国产一区| 搡老妇女老女人老熟妇| ponron亚洲|