Dataset Information

Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.

ABSTRACT: The taxonomic composition of microbial communities can be assessed using universal marker amplicon sequencing. The most common taxonomic markers are the 16S rDNA for bacterial communities and the internal transcribed spacer (ITS) region for fungal communities, but various other markers are used for barcoding eukaryotes. A crucial step in the bioinformatic analysis of amplicon sequences is the identification of representative sequences. This can be achieved using a clustering approach or by denoising raw sequencing reads. DADA2 is a widely adopted algorithm, released as an R library, that denoises marker-specific amplicons from next-generation sequencing and produces a set of representative sequences referred to as 'Amplicon Sequence Variants' (ASV). Here, we present Dadaist2, a modular pipeline, providing a complete suite for the analysis that ranges from raw sequencing reads to the statistics of numerical ecology. Dadaist2 implements a new approach that is specifically optimised for amplicons with variable lengths, such as the fungal ITS. The pipeline focuses on streamlining the data flow from the command line to R, with multiple options for statistical analysis and plotting, both interactive and automatic.

SUBMITTER: Ansorge R

PROVIDER: S-EPMC8157834 | biostudies-literature | 2021 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.

Ansorge Rebecca R Birolo Giovanni G James Stephen A SA Telatin Andrea A

International journal of molecular sciences 20210518 10

The taxonomic composition of microbial communities can be assessed using universal marker amplicon sequencing. The most common taxonomic markers are the 16S rDNA for bacterial communities and the internal transcribed spacer (ITS) region for fungal communities, but various other markers are used for barcoding eukaryotes. A crucial step in the bioinformatic analysis of amplicon sequences is the identification of representative sequences. This can be achieved using a clustering approach or by denoi ...[more]

PMID: 34069990

Similar Datasets

Project description:Low-cost phenological experiments with cut twigs are increasingly used to study bud development in response to spring warming and photoperiod. However, a broader variety of species needs to be tackled and in particular the influence of insufficient winter chilling deserves more attention. Therefore, we investigated if and how chilling requirements can be efficiently investigated by cut twigs and how this low-tech approach could be successfully implemented as a citizen science or school project. We conducted an experiment on bud burst and leaf development of Corylus avellana L. twigs, with natural chilling outdoors on a shrub (S) and another chilling treatment as cut twigs in containers (C), and subsequent forcing indoors. Subsampling of the number of cutting dates and number of twigs was used to infer minimum required sample sizes. Apart from insufficiently chilled twigs, ~80% of the twigs (both S and C) reached leaf out. For multiple definitions of chilling and forcing, a negative exponential relationship was revealed between chilling and amount of forcing needed to reach certain developmental stages. At least 5 out of 15 cutting dates or alternatively half of the 10 twig repetitions, but especially those mirroring low chilling conditions, were needed to describe the chilling-forcing relationship with some degree of robustness. In addition, for cutting dates with long chilling, i.e., from January onwards, freshly cut twigs (S) required significantly more forcing to reach bud burst than twigs from containers (C), although the effect was small. In general, chilling conditions of mature shrubs were well captured by cut twigs, therefore opening the possibility of chilling through refrigeration. We conclude that experimental protocols as outlined here are feasible for citizen scientists, school projects, and science education, and would have the potential to advance the research field if carried out on a large scale. We provide an easy-to-use Shiny simulation app to enable citizen scientists to build up a bud development model based on their own experimental data and then simulate future phenological development with winter and/or spring warming. This may encourage them to further study other aspects of climate change and the impacts of climate change.

Project description:BACKGROUND: Examining the distribution of variation has proven an extremely profitable technique in the effort to identify sequences of biological significance. Most approaches in the field, however, evaluate only the conserved portions of sequences - ignoring the biological significance of sequence differences. A suite of sophisticated likelihood based statistical models from the field of molecular evolution provides the basis for extracting the information from the full distribution of sequence variation. The number of different problems to which phylogeny-based maximum likelihood calculations can be applied is extensive. Available software packages that can perform likelihood calculations suffer from a lack of flexibility and scalability, or employ error-prone approaches to model parameterisation. RESULTS: Here we describe the implementation of PyEvolve, a toolkit for the application of existing, and development of new, statistical methods for molecular evolution. We present the object architecture and design schema of PyEvolve, which includes an adaptable multi-level parallelisation schema. The approach for defining new methods is illustrated by implementing a novel dinucleotide model of substitution that includes a parameter for mutation of methylated CpG's, which required 8 lines of standard Python code to define. Benchmarking was performed using either a dinucleotide or codon substitution model applied to an alignment of BRCA1 sequences from 20 mammals, or a 10 species subset. Up to five-fold parallel performance gains over serial were recorded. Compared to leading alternative software, PyEvolve exhibited significantly better real world performance for parameter rich models with a large data set, reducing the time required for optimisation from approximately 10 days to approximately 6 hours. CONCLUSION: PyEvolve provides flexible functionality that can be used either for statistical modelling of molecular evolution, or the development of new methods in the field. The toolkit can be used interactively or by writing and executing scripts. The toolkit uses efficient processes for specifying the parameterisation of statistical models, and implements numerous optimisations that make highly parameter rich likelihood functions solvable within hours on multi-cpu hardware. PyEvolve can be readily adapted in response to changing computational demands and hardware configurations to maximise performance. PyEvolve is released under the GPL and can be downloaded from http://cbis.anu.edu.au/software.

Dataset Information

Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.

Publications

Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets