DADA2: High-resolution sample inference from Illumina amplicon data.
Ontology highlight
ABSTRACT:
SUBMITTER: Callahan BJ
PROVIDER: S-EPMC4927377 | biostudies-literature | 2016 Jul
REPOSITORIES: biostudies-literature
Publications
DADA2: High-resolution sample inference from Illumina amplicon data.
Callahan Benjamin J BJ McMurdie Paul J PJ Rosen Michael J MJ Han Andrew W AW Johnson Amy Jo A AJ Holmes Susan P SP
Nature methods 20160523 7
We present the open-source software package DADA2 for modeling and correcting Illumina-sequenced amplicon errors (https://github.com/benjjneb/dada2). DADA2 infers sample sequences exactly and resolves differences of as little as 1 nucleotide. In several mock communities, DADA2 identified more real variants and output fewer spurious sequences than other methods. We applied DADA2 to vaginal samples from a cohort of pregnant women, revealing a diversity of previously undetected Lactobacillus crispa ...[more]
Similar Datasets
Project description:MotivationNext-generation amplicon sequencing is a powerful tool for investigating microbial communities. A main challenge is to distinguish true biological variants from errors caused by amplification and sequencing. In traditional analyses, such errors are eliminated by clustering reads within a sequence similarity threshold, usually 97%, and constructing operational taxonomic units, but the arbitrary threshold leads to low resolution and high false positive rates. Recently developed "denoising" methods have proven able to resolve single-nucleotide amplicon variants, but they still miss low frequency sequences, especially those near more frequent sequences, because they ignore the sequencing quality information.ResultsWe introduce AmpliCI, a reference-free, model-based method for rapidly resolving the number, abundance and identity of error-free sequences in massive Illumina amplicon datasets. AmpliCI takes into account quality information and allows the data, not an arbitrary threshold or an external database, to drive conclusions. AmpliCI estimates a finite mixture model, using a greedy strategy to gradually select error-free sequences and approximately maximize the likelihood. AmpliCI has better performance than three popular denoising methods, with acceptable computation time and memory usage.AvailabilitySource code is available at https://github.com/DormanLab/AmpliCI.Supplementary informationSupplementary material are available at Bioinformatics online.
| S-EPMC7850112 | biostudies-literature
Project description:BackgroundAmplicon sequencing of phylogenetic marker genes, e.g., 16S, 18S, or ITS ribosomal RNA sequences, is still the most commonly used method to determine the composition of microbial communities. Microbial ecologists often have expert knowledge on their biological question and data analysis in general, and most research institutes have computational infrastructures to use the bioinformatics command line tools and workflows for amplicon sequencing analysis, but requirements of bioinformatics skills often limit the efficient and up-to-date use of computational resources.ResultsWe present dadasnake, a user-friendly, 1-command Snakemake pipeline that wraps the preprocessing of sequencing reads and the delineation of exact sequence variants by using the favorably benchmarked and widely used DADA2 algorithm with a taxonomic classification and the post-processing of the resultant tables, including hand-off in standard formats. The suitability of the provided default configurations is demonstrated using mock community data from bacteria and archaea, as well as fungi.ConclusionsBy use of Snakemake, dadasnake makes efficient use of high-performance computing infrastructures. Easy user configuration guarantees flexibility of all steps, including the processing of data from multiple sequencing platforms. It is easy to install dadasnake via conda environments. dadasnake is available at https://github.com/a-h-b/dadasnake.
| S-EPMC7702218 | biostudies-literature
Project description:Identification and analysis of fungal communities commonly rely on internal transcribed spacer-based (ITS-based) amplicon sequencing. There is no gold standard used to infer and classify fungal constituents since methodologies have been adapted from analyses of bacterial communities. To achieve high-resolution inference of fungal constituents, we customized a DADA2-based pipeline using a mix of 11 medically relevant fungi. While DADA2 allowed the discrimination of ITS1 sequences differing by single nucleotides, quality filtering, sequencing bias, and database selection were identified as key variables determining the accuracy of sample inference. Due to species-specific differences in sequencing quality, default filtering settings removed most reads that originated from Aspergillus species, Saccharomyces cerevisiae, and Candida glabrata. By fine-tuning the quality filtering process, we achieved an improved representation of the fungal communities. By adapting a wobble nucleotide in the ITS1 forward primer region, we further increased the yield of S. cerevisiae and C. glabrata sequences. Finally, we showed that a BLAST-based algorithm based on the UNITE+INSD or the NCBI NT database achieved a higher reliability in species-level taxonomic annotation compared with the naive Bayesian classifier implemented in DADA2. These steps optimized a robust fungal ITS1 sequencing pipeline that, in most instances, enabled species-level assignment of community members.
| S-EPMC8765055 | biostudies-literature
Project description:Recent experiments have provided Hi-C data at resolution as high as 1 kbp. However, 3D structural inference from high-resolution Hi-C datasets is often computationally unfeasible using existing methods.We have developed miniMDS, an approximation of multidimensional scaling (MDS) that partitions a Hi-C dataset, performs high-resolution MDS separately on each partition, and then reassembles the partitions using low-resolution MDS. miniMDS is faster, more accurate, and uses less memory than existing methods for inferring the human genome at high resolution (10 kbp).A Python implementation of miniMDS is available on GitHub: https://github.com/seqcode/miniMDS .mahony@psu.edu.Supplementary data are available at Bioinformatics online.
| S-EPMC5870652 | biostudies-literature
Project description:Previously we presented Swarm v1, a novel and open source amplicon clustering program that produced fine-scale molecular operational taxonomic units (OTUs), free of arbitrary global clustering thresholds and input-order dependency. Swarm v1 worked with an initial phase that used iterative single-linkage with a local clustering threshold (d), followed by a phase that used the internal abundance structures of clusters to break chained OTUs. Here we present Swarm v2, which has two important novel features: (1) a new algorithm for d = 1 that allows the computation time of the program to scale linearly with increasing amounts of data; and (2) the new fastidious option that reduces under-grouping by grafting low abundant OTUs (e.g., singletons and doubletons) onto larger ones. Swarm v2 also directly integrates the clustering and breaking phases, dereplicates sequencing reads with d = 0, outputs OTU representatives in fasta format, and plots individual OTUs as two-dimensional networks.
| S-EPMC4690345 | biostudies-literature
Project description:BackgroundSequencing of the 16S rRNA gene has been the standard for studying the composition of microbial communities. While it allows identification of bacteria at the level of species, this method does not usually provide sufficient information to resolve communities at the sub-species level. Species-level resolution is not adequate for studies of transmission or stability or for exploring subspecies variation in disease association. Strain level analysis using whole metagenome shotgun sequencing has significant limitations that can make it unsuitable for large-scale studies. Achieving sufficient depth of sequencing can be cost-prohibitive, and even with adequate coverage, deconvoluting complex communities such as the oral microbiota is computationally very challenging. Thus, there is a need for high-resolution, yet cost-effective, high-throughput methods for characterizing microbial communities.ResultsSignificant improvement in resolution for amplicon-based bacterial community analysis was achieved by combining amplicon sequencing of a high-diversity marker gene, the ribosomal 16-23S intergenic spacer region (ISR), with a probabilistic error modeling based denoising algorithm, DADA2. The resolving power of this new approach was compared to that of both standard and high-resolution 16S-based approaches using a set of longitudinal subgingival plaque samples. The ISR strategy resulted in a 5.2-fold increase in community resolution compared to reference-based 16S rRNA gene analysis and showed 100% accuracy in predicting the correct source of a clinical sample. Individuals' microbial communities were highly personalized, and although they exhibited some drift in membership and levels over time, that difference was always smaller than the differences between any two subjects, even after 1 year. The construction of an ISR database from publicly available genomic sequences allowed us to explore genomic variation within species, resulting in the identification of multiple variants of the ISR for most species.ConclusionsThe ISR approach resulted in significantly improved resolution of communities and revealed a highly personalized human oral microbiota that was stable over 1 year. Multiple ISR types were observed for all species examined, demonstrating a high level of subspecies variation in the oral microbiota. The approach is high-throughput, high-resolution yet cost-effective, allowing subspecies-level community fingerprinting at a cost comparable to that of 16S rRNA gene amplicon sequencing. It will be useful for a range of applications that require high-resolution identification of organisms, including microbial tracking, community fingerprinting, and potentially for identification of virulence-associated strains.
| S-EPMC6126016 | biostudies-literature
Project description:Trajectory inference (TI) methods infer cell developmental trajectory from single-cell RNA sequencing data. Current TI methods can be categorized into those using RNA velocity information and those using only single-cell gene expression data. The latter type of methods are restricted to certain trajectory structures, and cannot determine cell developmental direction. Recently proposed TI methods using RNA velocity information have limited accuracy. We present CellPath, a method that infers cell trajectories by integrating single-cell gene expression and RNA velocity information. CellPath overcomes the restrictions of TI methods that do not use RNA velocity information: it can find multiple high-resolution trajectories without constraints on the trajectory structure, and can automatically detect the direction of each trajectory path. We evaluate CellPath on both real and simulated datasets and show that CellPath finds more accurate and detailed trajectories than the state-of-the-art TI methods using or not using RNA velocity information.
| S-EPMC9017235 | biostudies-literature
Project description:In Saccharomyces cerevisiae, specific genes physically relocate from the nucleoplasm to the nuclear periphery concomitant with transcriptional activation, where they associate with the nuclear pore complex (NPC). We took a genomics approach in order to gain insight into the universality and physiological relevance of the interaction between active genes and the NPC. Using synthetic genetic array (SGA) approach, we identified interactions between components of the SAGA histone acetyltransferase complex and the Mlp and Nup60 subunits of the NPC. Cells lacking these SAGA and NPC components display growth defects under optimal growth conditions, in which cells are grown on rich medium containing the preferred sugar glucose as a carbon source and incubated at their optimal temperature. That growth defects were observed under these non-stress conditions suggests that these interactions are indicative of defects in normal cell physiology. These results are consistent with a model where physical interactions between the NPC and SAGA are important for transcription of constitutively expressed genes. To test this hypothesis, we used microarray analysis to assess changes global transcript levels in the absence of Nup60, Ada2, or Nup60 and Ada2. Microarray analysis reveals that the growth defect in these double mutants is correlated with a synthetic reduction in steady-state transcript levels for numerous genes that are strongly expressed in wildtype cells. These results suggest that SAGA and Nup60 cooperate in transcriptional regulation, and are consistent with a model for global regulation of SAGA-dependent transcription at the NPC.
2013-05-13 | GSE21870 | GEO
Project description:Determining the composition of bacterial communities beyond the level of a genus or species is challenging because of the considerable overlap between genomes representing close relatives. Here, we present the mSWEEP pipeline for identifying and estimating the relative sequence abundances of bacterial lineages from plate sweeps of enrichment cultures. mSWEEP leverages biologically grouped sequence assembly databases, applying probabilistic modelling, and provides controls for false positive results. Using sequencing data from major pathogens, we demonstrate significant improvements in lineage quantification and detection accuracy. Our pipeline facilitates investigating cultures comprising mixtures of bacteria, and opens up a new field of plate sweep metagenomics.
| S-EPMC8543175 | biostudies-literature
Project description:In high-throughput sequencing data, performance comparisons between computational tools are essential for making informed decisions at each step of a project. Simulations are a critical part of method comparisons, but for standard Illumina sequencing of genomic DNA, they are often oversimplified, which leads to optimistic results for most tools. ReSeq improves the authenticity of synthetic data by extracting and reproducing key components from real data. Major advancements are the inclusion of systematic errors, a fragment-based coverage model and sampling-matrix estimates based on two-dimensional margins. These improvements lead to more faithful performance evaluations. ReSeq is available at https://github.com/schmeing/ReSeq .
| S-EPMC7896392 | biostudies-literature