Dataset Information

Evaluation of viral genome assembly and diversity estimation in deep metagenomes.

ABSTRACT:

Background

Viruses have unique properties, small genome and regions of high similarity, whose effects on metagenomic assemblies have not been characterized so far. This study uses diverse in silico simulated viromes to evaluate how extensively genomes can be assembled using different sequencing platforms and assemblers. Further, it investigates the suitability of different methods to estimate viral diversity in metagenomes.

Results

We created in silico metagenomes mimicking various platforms at different sequencing depths. The CLC assembler revealed subpar compared to IDBA_UD and CAMERA , which are metagenomic-specific. Up to a saturation point, Illumina platforms proved more capable of reconstructing large portions of viral genomes compared to 454. Read length was an important factor for limiting chimericity, while scaffolding marginally improved contig length and accuracy. The genome length of the various viruses in the metagenomes did not significantly affect genome reconstruction, but the co-existence of highly similar genomes was detrimental. When evaluating diversity estimation tools, we found that PHACCS results were more accurate than those from CatchAll and clustering, which were both orders of magnitude above expected.

Conclusions

Assemblers designed specifically for the analysis of metagenomes should be used to facilitate the creation of high-quality long contigs. Despite the high coverage possible, scientists should not expect to always obtain complete genomes, because their reconstruction may be hindered by co-existing species bearing highly similar genomic regions. Further development of metagenomics-oriented assemblers may help bypass these limitations in future studies. Meanwhile, the lack of fully reconstructed communities keeps methods to estimate viral diversity relevant. While none of the three methods tested had absolute precision, only PHACCS was deemed suitable for comparative studies.

SUBMITTER: Aguirre de Carcer D

PROVIDER: S-EPMC4247695 | biostudies-literature | 2014 Nov

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Evaluation of viral genome assembly and diversity estimation in deep metagenomes.

Aguirre de Cárcer Daniel D Angly Florent E FE Alcamí Antonio A

BMC genomics 20141118

<h4>Background</h4>Viruses have unique properties, small genome and regions of high similarity, whose effects on metagenomic assemblies have not been characterized so far. This study uses diverse in silico simulated viromes to evaluate how extensively genomes can be assembled using different sequencing platforms and assemblers. Further, it investigates the suitability of different methods to estimate viral diversity in metagenomes.<h4>Results</h4>We created in silico metagenomes mimicking variou ...[more]

PMID: 25407630

Similar Datasets

Project description:BackgroundMetagenomics has transformed our understanding of microbial diversity across ecosystems, with recent advances enabling de novo assembly of genomes from metagenomes. These metagenome-assembled genomes are critical to provide ecological, evolutionary, and metabolic context for all the microbes and viruses yet to be cultivated. Metagenomes can now be generated from nanogram to subnanogram amounts of DNA. However, these libraries require several rounds of PCR amplification before sequencing, and recent data suggest these typically yield smaller and more fragmented assemblies than regular metagenomes.MethodsHere we evaluate de novo assembly methods of 169 PCR-amplified metagenomes, including 25 for which an unamplified counterpart is available, to optimize specific assembly approaches for PCR-amplified libraries. We first evaluated coverage bias by mapping reads from PCR-amplified metagenomes onto reference contigs obtained from unamplified metagenomes of the same samples. Then, we compared different assembly pipelines in terms of assembly size (number of bp in contigs ≥ 10 kb) and error rates to evaluate which are the best suited for PCR-amplified metagenomes.ResultsRead mapping analyses revealed that the depth of coverage within individual genomes is significantly more uneven in PCR-amplified datasets versus unamplified metagenomes, with regions of high depth of coverage enriched in short inserts. This enrichment scales with the number of PCR cycles performed, and is presumably due to preferential amplification of short inserts. Standard assembly pipelines are confounded by this type of coverage unevenness, so we evaluated other assembly options to mitigate these issues. We found that a pipeline combining read deduplication and an assembly algorithm originally designed to recover genomes from libraries generated after whole genome amplification (single-cell SPAdes) frequently improved assembly of contigs ≥10 kb by 10 to 100-fold for low input metagenomes.ConclusionsPCR-amplified metagenomes have enabled scientists to explore communities traditionally challenging to describe, including some with extremely low biomass or from which DNA is particularly difficult to extract. Here we show that a modified assembly pipeline can lead to an improved de novo genome assembly from PCR-amplified datasets, and enables a better genome recovery from low input metagenomes.

Project description:Background & aimsWe aimed to investigate how viral quasispecies of the HBV whole genome evolves and diversifies in response to HBeAg seroconversion and viral control utilising next-generation sequencing (NGS).MethodsFifty HBeAg-positive chronic hepatitis B patients, including 18 treatment-naïve and 32 interferon (IFN)-treated individuals, were recruited. Serial HBV whole genomes in serum were analysed by NGS to determine sequence characteristics and viral quasispecies.ResultsHBV quasispecies diversity, measured by nucleotide diversity, was negatively correlated with viral load and hepatitis activity. Spontaneous HBeAg seroconverters exhibited significantly greater viral quasispecies diversity than treatment-naïve non-seroconverters from >1 year before seroconversion (0.0112 vs. 0.0060, p <0.01) to >1 year after seroconversion (0.0103 vs. 0.0068, p <0.01). IFN-induced HBeAg seroconverters tended to have higher viral genetic diversity than non-seroconverters along with treatment. Particularly, the IFN responders, defined as IFN-induced HBeAg seroconversion with low viraemia, exhibited significantly greater genetic diversity of whole HBV genome at 6 months post-IFN treatment than IFN non-responders (0.0148 vs. 0.0106, p = 0.048). Moreover, spontaneous HBeAg seroconverters and IFN responders exhibited significantly higher evolutionary rates and more intra-host single-nucleotide variants. Interestingly, in spontaneous HBeAg seroconverters and IFN responders, there were distinct evolutionary patterns in the HBV genome.ConclusionsHigher HBV quasispecies diversity is associated with spontaneous HBeAg seroconversion and IFN-induced HBeAg seroconversion with low viraemia, conferring a favourable clinical outcome.Lay summaryHBeAg seroconversion is a landmark in the natural history of chronic HBV infection. Using next-generation sequencing, we found that the nucleotide diversity of HBV was negatively correlated with viral load and hepatitis activity. Patients undergoing HBeAg seroconversion had more diverse HBV genomes and a faster viral evolution rate. Our findings suggest HBeAg seroconversion is driven by host selection pressure, likely immune selection pressure.

Project description:Shotgun metagenomics has greatly advanced our understanding of microbial communities over the last decade. Metagenomic analyses often include assembly and genome binning, computationally daunting tasks especially for big data from complex environments such as soil and sediments. In many studies, however, only a subset of genes and pathways involved in specific functions are of interest; thus, it is not necessary to attempt global assembly. In addition, methods that target genes can be computationally more efficient and produce more accurate assembly by leveraging rich databases, especially for those genes that are of broad interest such as those involved in biogeochemical cycles, biodegradation, and antibiotic resistance or used as phylogenetic markers. Here, we review six gene-targeted assemblers with unique algorithms for extracting and/or assembling targeted genes: Xander, MegaGTA, SAT-Assembler, HMM-GRASPx, GenSeed-HMM, and MEGAN. We tested these tools using two datasets with known genomes, a synthetic community of artificial reads derived from the genomes of 17 bacteria, shotgun sequence data from a mock community with 48 bacteria and 16 archaea genomes, and a large soil shotgun metagenomic dataset. We compared assemblies of a universal single copy gene (rplB) and two N cycle genes (nifH and nirK). We measured their computational efficiency, sensitivity, specificity, and chimera rate and found Xander and MegaGTA, which both use a probabilistic graph structure to model the genes, have the best overall performance with all three datasets, although MEGAN, a reference matching assembler, had better sensitivity with synthetic and mock community members chosen from its reference collection. Also, Xander and MegaGTA are the only tools that include post-assembly scripts tuned for common molecular ecology and diversity analyses. Additionally, we provide a mathematical model for estimating the probability of assembling targeted genes in a metagenome for estimating required sequencing depth.

Dataset Information

Evaluation of viral genome assembly and diversity estimation in deep metagenomes.

Background

Results

Conclusions

Publications

Evaluation of viral genome assembly and diversity estimation in deep metagenomes.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets