Dataset Information

A statistical model for describing and simulating microbial community profiles.

ABSTRACT: Many methods have been developed for statistical analysis of microbial community profiles, but due to the complex nature of typical microbiome measurements (e.g. sparsity, zero-inflation, non-independence, and compositionality) and of the associated underlying biology, it is difficult to compare or evaluate such methods within a single systematic framework. To address this challenge, we developed SparseDOSSA (Sparse Data Observations for the Simulation of Synthetic Abundances): a statistical model of microbial ecological population structure, which can be used to parameterize real-world microbial community profiles and to simulate new, realistic profiles of known structure for methods evaluation. Specifically, SparseDOSSA's model captures marginal microbial feature abundances as a zero-inflated log-normal distribution, with additional model components for absolute cell counts and the sequence read generation process, microbe-microbe, and microbe-environment interactions. Together, these allow fully known covariance structure between synthetic features (i.e. "taxa") or between features and "phenotypes" to be simulated for method benchmarking. Here, we demonstrate SparseDOSSA's performance for 1) accurately modeling human-associated microbial population profiles; 2) generating synthetic communities with controlled population and ecological structures; 3) spiking-in true positive synthetic associations to benchmark analysis methods; and 4) recapitulating an end-to-end mouse microbiome feeding experiment. Together, these represent the most common analysis types in assessment of real microbial community environmental and epidemiological statistics, thus demonstrating SparseDOSSA's utility as a general-purpose aid for modeling communities and evaluating quantitative methods. An open-source implementation is available at http://huttenhower.sph.harvard.edu/sparsedossa2.

SUBMITTER: Ma S

PROVIDER: S-EPMC8491899 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:BackgroundSocial-stress mouse model, based on the resident-intruder paradigm was used to simulate features of human post-traumatic stress disorder (PTSD). The model involved exposure of an intruder (subject) mouse to a resident aggressor mouse followed by exposure to trauma reminders with rest periods. C57BL/6 mice exposed to SJL aggressor mice exhibited behaviors suggested as PTSD-in-mouse phenotypes: intermittent freezing, reduced locomotion, avoidance of the aggressor-associated cue and apparent startled jumping. Brain tissues (amygdala, hippocampus, medial prefrontal cortex, septal region, corpus striatum and ventral striatum) from subject (aggressor exposed: Agg-E) and control C57BL/6 mice were collected at one, 10 and 42 days post aggressor exposure sessions. Transcripts in these brain regions were assayed using Agilent's mouse genome-wide arrays.ResultsPathways and biological processes associated with differentially regulated genes were mainly those thought to be involved in fear-related behavioral responses and neuronal signaling. Expression-based assessments of activation patterns showed increased activations of pathways related to anxiety disorders (hyperactivity and fear responses), impaired cognition, mood disorders, circadian rhythm disruption, and impaired territorial and aggressive behaviors. In amygdala, activations of these pathways were more pronounced at earlier time-points, with some attenuation after longer rest periods. In hippocampus and medial prefrontal cortex, activation patterns were observed at later time points. Signaling pathways associated with PTSD-comorbid conditions, such as diabetes, metabolic disorder, inflammation and cardiac infarction, were also significantly enriched. In contrast, signaling processes related to neurogenesis and synaptic plasticity were inhibited.ConclusionsOur data suggests activations of behavioral responses associated with anxiety disorders as well as inhibition of neuronal signaling pathways important for neurogenesis, cognition and extinction of fear memory. These pathways along with comorbid-related signaling pathways indicate the pervasive and multisystem effects of aggressor exposure in mice, potentially mirroring the pathologic conditions of PTSD patients.

Project description:Soil microbial community assembly is crucial for understanding the mechanisms of microbial communities that regulate ecosystem-level functioning. The relative contributions of stochastic and deterministic processes to microbial community assembly remain poorly defined, and major questions exist concerning the soil organic carbon (SOC) dynamics of microbial community assembly in deep soil. Here, the bacterial community assembly processes were explored across five soil profile depths (up to 80 cm) during a 15-year field experiment involving four fertilization regimes. We found that the bacterial community assembly was initially governed by deterministic selection in topsoil but was progressively structured by increasing stochastic dispersal with depth. The migration rate (m) and β-null deviation pattern supported the hypothesis of a relatively greater influence of dispersal in deep soil, which was correlated with bacterial community assembly by stochastic processes. These changes in the entire community assembly reflected consistent assembly processes of the two most dominant phyla, Acidobacteria and Chloroflexi Structural equation modeling showed that soil features (pH and total phosphorus) and bacterial interactions (competition and network complexity) were significantly related to bacterial community assembly in the 0-to-10-cm and 10-to-20-cm layers. Partial Mantel tests, structural equation modeling, and random forest modeling consistently indicated a strong and significant correlation between bacterial community assemblages and SOC dynamics, implying that bacterial assembly processes would potentially suppress SOC metabolism and mineralization when the contributions of stochastic dispersal to communities increased in deeper layers. Our results have important implications for integrating bacterial community assembly processes into the predictions of SOC dynamics.IMPORTANCE We have provided a framework to better understand the mechanisms governing the balance between stochastic and deterministic processes and to integrate the shifts in community assembly processes with microbial carbon metabolism. Our study reinforced that environmental filtering and bacterial cooccurrence patterns influence the stochastic/deterministic continuum of soil bacterial community assembly and that stochasticity may act through deeper soil layers to influence carbon metabolism. Delineating theoretically the potential linkages between community assembly and SOC dynamics across a broad range of microbial systems represents an interesting topic for future research.

Project description:Uncertainty analysis is the process of identifying limitations in scientific knowledge and evaluating their implications for scientific conclusions. In the context of microbial risk assessment, the uncertainty in the predicted microbial behavior can be an important component of the overall uncertainty. Conventional deterministic modeling approaches which provide point estimates of the pathogen's levels cannot quantify the uncertainty around the predictions. The objective of this study was to use Bayesian statistical modeling for describing uncertainty in predicted microbial thermal inactivation of Salmonella enterica Typhimurium DT104. A set of thermal inactivation data in broth with water activity adjusted to 0.75 at 9 different temperature conditions obtained from the ComBase database (www.combase.cc) was used. A log-linear microbial inactivation was used as a primary model while for secondary modeling, a linear relation between the logarithm of inactivation rate and temperature was assumed. For comparison, data were fitted with a two-step and a global Bayesian regression. Posterior distributions of model's parameters were used to predict Salmonella thermal inactivation. The combination of the joint posterior distributions of model's parameters allowed the prediction of cell density over time, total reduction time and inactivation rate as probability distributions at different time and temperature conditions. For example, for the time required to eliminate a Salmonella population of about 107 CFU/ml at 65°C, the model predicted a time distribution with a median of 0.40 min and 5th and 95th percentiles of 0.24 and 0.60 min, respectively. The validation of the model showed that it can describe successfully uncertainty in predicted thermal inactivation with most observed data being within the 95% prediction intervals of the model. The global regression approach resulted in less uncertain predictions compared to the two-step regression. The developed model could be used to quantify uncertainty in thermal inactivation in risk-based processing design as well as in risk assessment studies.

Dataset Information

A statistical model for describing and simulating microbial community profiles.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets