Dadaist2: A Toolkit to Automate and Simplify Statistical Analysis and Plotting of Metabarcoding Experiments.
Ontology highlight
ABSTRACT:
SUBMITTER: Ansorge R
PROVIDER: S-EPMC8157834 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Biological small-angle X-ray scattering (SAXS) is an increasingly popular technique used to obtain nanoscale structural information on macromolecules in solution. However, radiation damage to the samples limits the amount of useful data that can be collected from a single sample. In contrast to the extensive analytical resources available for macromolecular crystallography (MX), there are relatively few tools to quantitate radiation damage for SAXS, some of which require a significant level of manual characterization, with the potential of leading to conflicting results from different studies. Here, computational tools have been developed to automate and standardize radiation damage analysis for SAXS data. RADDOSE-3D, a dose calculation software utility originally written for MX experiments, has been extended to account for the cylindrical geometry of the capillary tube, the liquid composition of the sample and the attenuation of the beam by the capillary material to allow doses to be calculated for many SAXS experiments. Furthermore, a library has been written to visualize and explore the pairwise similarity of frames. The calculated dose for the frame at which three subsequent frames are determined to be dissimilar is defined as the radiation damage onset threshold (RDOT). Analysis of RDOTs has been used to compare the efficacy of radioprotectant compounds to extend the useful lifetime of SAXS samples. Comparison of the RDOTs shows that, for radioprotectant compounds at 5 and 10 mM concentration, glycerol is the most effective compound. However, at 1 and 2 mM concentrations, dithiothreitol (DTT) appears to be most effective. Our newly developed visualization library contains methods that highlight the unusual radiation damage results given by SAXS data collected using higher concentrations of DTT: these observations should pave the way to the development of more sophisticated frame merging strategies.
| S-EPMC5182020 | biostudies-literature
Project description:Validating statistical analysis methods for RNA sequencing (RNA-seq) experiments is a complex task. Researchers often find themselves having to decide between competing models or assessing the reliability of results obtained with a designated analysis program. Computer simulation has been the most frequently used procedure to verify the adequacy of a model. However, datasets generated by simulations depend on the parameterization and the assumptions of the selected model. Moreover, such datasets may constitute a partial representation of reality as the complexity or RNA-seq data is hard to mimic. We present the use of plasmode datasets to complement the evaluation of statistical models for RNA-seq data. A plasmode is a dataset obtained from experimental data but for which come truth is known. Using a set of simulated scenarios of technical and biological replicates, and public available datasets, we illustrate how to design algorithms to construct plasmodes under different experimental conditions. We contrast results from two types of methods for RNA-seq: (1) models based on negative binomial distribution (edgeR and DESeq), and (2) Gaussian models applied after transformation of data (MAANOVA). Results emphasize the fact that deciding what method to use may be experiment-specific due to the unknown distributions of expression levels. Plasmodes may contribute to choose which method to apply by using a similar pre-existing dataset. The promising results obtained from this approach, emphasize the need of promoting and improving systematic data sharing across the research community to facilitate plasmode building. Although we illustrate the use of plasmode for comparing differential expression analysis models, the flexibility of plasmode construction allows comparing upstream analysis, as normalization procedures or alignment pipelines, as well.
| S-EPMC3775431 | biostudies-literature
Project description:Ultrafast diffraction at X-ray free-electron lasers (XFELs) has the potential to yield new insights into important biological systems that produce radiation-sensitive crystals. An unavoidable feature of the `diffraction before destruction' nature of these experiments is that images are obtained from many distinct crystals and/or different regions of the same crystal. Combined with other sources of XFEL shot-to-shot variation, this introduces significant heterogeneity into the diffraction data, complicating processing and interpretation. To enable researchers to get the most from their collected data, a toolkit is presented that provides insights into the quality of, and the variation present in, serial crystallography data sets. These tools operate on the unmerged, partial intensity integration results from many individual crystals, and can be used on two levels: firstly to guide the experimental strategy during data collection, and secondly to help users make informed choices during data processing.
| S-EPMC4321488 | biostudies-literature
Project description:Octopus-toolkit is a stand-alone application for retrieving and processing large sets of next-generation sequencing (NGS) data with a single step. Octopus-toolkit is an automated set-up-and-analysis pipeline utilizing the Aspera, SRA Toolkit, FastQC, Trimmomatic, HISAT2, STAR, Samtools, and HOMER applications. All the applications are installed on the user's computer when the program starts. Upon the installation, it can automatically retrieve original files of various epigenomic and transcriptomic data sets, including ChIP-seq, ATAC-seq, DNase-seq, MeDIP-seq, MNase-seq and RNA-seq, from the gene expression omnibus data repository. The downloaded files can then be sequentially processed to generate BAM and BigWig files, which are used for advanced analyses and visualization. Currently, it can process NGS data from popular model genomes such as, human (Homo sapiens), mouse (Mus musculus), dog (Canis lupus familiaris), plant (Arabidopsis thaliana), zebrafish (Danio rerio), fruit fly (Drosophila melanogaster), worm (Caenorhabditis elegans), and budding yeast (Saccharomyces cerevisiae) genomes. With the processed files from Octopus-toolkit, the meta-analysis of various data sets, motif searches for DNA-binding proteins, and the identification of differentially expressed genes and/or protein-binding sites can be easily conducted with few commands by users. Overall, Octopus-toolkit facilitates the systematic and integrative analysis of available epigenomic and transcriptomic NGS big data.
| S-EPMC5961211 | biostudies-literature
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.
| S-EPMC7498619 | biostudies-literature
Project description:ChIP-exo/nexus experiments present modifications on the commonly used ChIP-seq protocol for high resolution mapping of transcription factor binding sites. Although many aspects of the ChIP-exo data analysis are similar to those of ChIP-seq, these high throughput experiments pose a number of unique quality control and analysis challenges. We develop a statistical quality control pipeline and accompanying R package, ChIPexoQual, to enable exploration and analysis of ChIP-exo and related experiments. ChIPexoQual evaluates a number of key issues including strand imbalance, library complexity, and signal enrichment of data. Assessment of these features are facilitated through diagnostic plots and summary statistics calculated over regions of the genome with varying levels of coverage. We evaluated our QC pipeline with both large collections of public ChIP-exo/nexus data and multiple, new ChIP-exo datasets from E. coli. ChIPexoQual analysis of these datasets resulted in guidelines for using these QC metrics across a wide range of sequencing depths and provided further insights for modelling ChIP-exo data. Finally, although ChIP-exo experiments have been compared to ChIP-seq experiments with single-end (SE) sequencing, we provide, for the first time, comparisons with paired-end (PE) ChIP-seq experiments. We illustrate that, at fixed sequencing depths, ChIP-exo provides higher sensitivity, specificity, and spatial resolution than PE ChIP-seq and both significantly outperform their SE ChIP-seq counterpart.
2017-03-01 | GSE84830 | GEO
Project description:ChIP-exo/nexus experiments rely on innovative modifications of the commonly used ChIP-seq protocol for high resolution mapping of transcription factor binding sites. Although many aspects of the ChIP-exo data analysis are similar to those of ChIP-seq, these high throughput experiments pose a number of unique quality control and analysis challenges. We develop a novel statistical quality control pipeline and accompanying R/Bioconductor package, ChIPexoQual, to enable exploration and analysis of ChIP-exo and related experiments. ChIPexoQual evaluates a number of key issues including strand imbalance, library complexity, and signal enrichment of data. Assessment of these features are facilitated through diagnostic plots and summary statistics computed over regions of the genome with varying levels of coverage. We evaluated our QC pipeline with both large collections of public ChIP-exo/nexus data and multiple, new ChIP-exo datasets from Escherichia coli. ChIPexoQual analysis of these datasets resulted in guidelines for using these QC metrics across a wide range of sequencing depths and provided further insights for modelling ChIP-exo data.
| S-EPMC5587812 | biostudies-literature
Project description:We have studied the primary nucleation of adipic acid from aqueous solutions in thousands of microdroplets generated in a fully automated microfluidic setup. By varying supersaturation in solution and residence time, we were able to estimate nucleation rates and growth times, while accounting for the stochastic nature of nucleation, the variability in microdroplet volumes (which is kept below 2%, thanks to a carefully designed experimental protocol), and the uncertainty in the automated image analysis procedure. Through a thorough statistical analysis we have obtained exact expressions for the expected values and the variances of all the random variables involved, all the way to the nucleation rate and the growth time associated with each supersaturation level explored and to the model parameters appearing in the corresponding constitutive equations. We have analyzed what controls the overall uncertainty in the estimation of the physical quantities above. We have shown that the distribution of droplet volumes at the level observed here is not limiting, whereas the detection technique and the image analysis algorithm play a critical role, together with the fact that the supersaturation levels and residence times that can be reasonably explored are limited. The tools and methods presented and made available to the scientific community will help in making microfluidics-based studies of nucleation more effective.
| S-EPMC6961308 | biostudies-literature
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.
| S-EPMC317364 | biostudies-literature
Project description:BACKGROUND: Systems biology experiments studying different topics and organisms produce thousands of data values across different types of genomic data. Further, data mining analyses are yielding ranked and heterogeneous results and association networks distributed over the entire genome. The visualization of these results is often difficult and standalone web tools allowing for custom inputs and dynamic filtering are limited. RESULTS: We have developed POMO (http://pomo.cs.tut.fi), an interactive web-based application to visually explore omics data analysis results and associations in circular, network and grid views. The circular graph represents the chromosome lengths as perimeter segments, as a reference outer ring, such as cytoband for human. The inner arcs between nodes represent the uploaded network. Further, multiple annotation rings, for example depiction of gene copy number changes, can be uploaded as text files and represented as bar, histogram or heatmap rings. POMO has built-in references for human, mouse, nematode, fly, yeast, zebrafish, rice, tomato, Arabidopsis, and Escherichia coli. In addition, POMO provides custom options that allow integrated plotting of unsupported strains or closely related species associations, such as human and mouse orthologs or two yeast wild types, studied together within a single analysis. The web application also supports interactive label and weight filtering. Every iterative filtered result in POMO can be exported as image file and text file for sharing or direct future input. CONCLUSIONS: The POMO web application is a unique tool for omics data analysis, which can be used to visualize and filter the genome-wide networks in the context of chromosomal locations as well as multiple network layouts. With the several illustration and filtering options the tool supports the analysis and visualization of any heterogeneous omics data analysis association results for many organisms. POMO is freely available and does not require any installation or registration.
| S-EPMC3880012 | biostudies-literature