Ecological modeling from time-series inference: insight into dynamics and stability of intestinal microbiota.
Ontology highlight
ABSTRACT:
SUBMITTER: Stein RR
PROVIDER: S-EPMC3861043 | biostudies-literature | 2013
REPOSITORIES: biostudies-literature
Publications
Ecological modeling from time-series inference: insight into dynamics and stability of intestinal microbiota.
Stein Richard R RR Bucci Vanni V Toussaint Nora C NC Buffie Charlie G CG Rätsch Gunnar G Pamer Eric G EG Sander Chris C Xavier João B JB
PLoS computational biology 20131212 12
The intestinal microbiota is a microbial ecosystem of crucial importance to human health. Understanding how the microbiota confers resistance against enteric pathogens and how antibiotics disrupt that resistance is key to the prevention and cure of intestinal infections. We present a novel method to infer microbial community ecology directly from time-resolved metagenomics. This method extends generalized Lotka-Volterra dynamics to account for external perturbations. Data from recent experiments ...[more]
Similar Datasets
Project description:The gastrointestinal (GI) microbiome is a densely populated ecosystem where dynamics are determined by interactions between microbial community members, as well as host factors. The spatial organization of this system is thought to be important in human health, yet this aspect of our resident microbiome is still poorly understood. In this study, we report significant spatial structure of the GI microbiota, and we identify general categories of spatial patterning in the distribution of microbial taxa along a healthy human GI tract. We further estimate the biotic interaction structure in the GI microbiota, both through time series and cooccurrence modeling of microbial community data derived from a large number of sequentially collected fecal samples. Comparison of these two approaches showed that species pairs involved in significant negative interactions had strong positive contemporaneous correlations and vice versa, while for species pairs without significant interactions, contemporaneous correlations were distributed around zero. We observed similar patterns when comparing these models to the spatial correlations between taxa identified in the adherent microbiota. This suggests that colocalization of microbial taxon pairs, and thus the spatial organization of the GI microbiota, is driven, at least in part, by direct or indirect biotic interactions. Thus, our study can provide a basis for an ecological interpretation of the biogeography of the human gut. IMPORTANCE The human gut microbiome is the subject of intense study due to its importance in health and disease. The majority of these studies have been based on the analysis of feces. However, little is known about how the microbial composition in fecal samples relates to the spatial distribution of microbial taxa along the gastrointestinal tract. By characterizing the microbial content both in intestinal tissue samples and in fecal samples obtained daily, we provide a conceptual framework for how the spatial structure relates to biotic interactions on the community level. We further describe general categories of spatial distribution patterns and identify taxa conforming to these categories. To our knowledge, this is the first study combining spatial and temporal analyses of the human gut microbiome. This type of analysis can be used for identifying candidate probiotics and designing strategies for clinical intervention.
| S-EPMC5585691 | biostudies-literature
Project description:Estimation of the period length of time-course data from cyclical biological processes, such as those driven by the circadian pacemaker, is crucial for inferring the properties of the biological clock found in many living organisms. We propose a methodology for period estimation based on spectrum resampling (SR) techniques. Simulation studies show that SR is superior and more robust to non-sinusoidal and noisy cycles than a currently used routine based on Fourier approximations. In addition, a simple fit to the oscillations using linear least squares is available, together with a non-parametric test for detecting changes in period length which allows for period estimates with different variances, as frequently encountered in practice. The proposed methods are motivated by and applied to various data examples from chronobiology.
| S-EPMC3988453 | biostudies-literature
Project description:Recent and ongoing revolutions in measurement technologies imply completely new possibilities for genome research: today, time-resolved, quantitative, and systems-level data are available. Nevertheless, without a corresponding revolution in methods for data analysis, these new data tend to drown researchers and doctors, rather than provide clear and useful insights. Such new methods are developed within the field of systems biology. Systems biology has two main approaches: mechanistically detailed and well-determined simulation models for small subsystems, and more approximative statistical models for the entire genome. However, there are few, if any, methods that combine the strengths of these two approaches. Herein, we present LASSIM, a new simulation-based approach, which can be applied to systems of the size of the entire genome. The superior performance of LASSIM is demonstrated in three examples: i) an example with simulated data shows that unlike traditional large-scale methods, LASSIM correctly identifies the true behavior between measured data-points, ii) LASSIM outperforms the winner of a previous DREAM challenge, the most competitive benchmarking approach available, iii) based on new data from TH2 differentiation, LASSIM identifies a first mechanistic model for the entire genome. The key predictions of this model are typically enriched for DNA bindings, which suggests that most predicted interactions are direct. Moreover, in silico knockdowns were experimentally validated. In summary, LASSIM opens the door to a new type of model-based data analysis: to models that combine the strengths of reliable mechanistic models with truly systems-level data.
2017-06-12 | GSE60681 | GEO
Project description:Live-cell imaging can be used to capture spatio-temporal aspects of cellular responses that are not accessible to fixed-cell imaging. As the use of live-cell imaging continues to increase, new computational procedures are needed to characterize and classify the temporal dynamics of individual cells. For this purpose, here we present the general experimental-computational framework SAPHIRE (Stochastic Annotation of Phenotypic Individual-cell Responses) to characterize phenotypic cellular responses from time series imaging datasets. Hidden Markov modeling is used to infer and annotate morphological state and state-switching properties from image-derived cell shape measurements. Time series modeling is performed on each cell individually, making the approach broadly useful for analyzing asynchronous cell populations. Two-color fluorescent cells simultaneously expressing actin and nuclear reporters enabled us to profile temporal changes in cell shape following pharmacological inhibition of cytoskeleton-regulatory signaling pathways. Results are compared with existing approaches conventionally applied to fixed-cell imaging datasets, and indicate that time series modeling captures heterogeneous dynamic cellular responses that can improve drug classification and offer additional important insight into mechanisms of drug action. The software is available at http://saphire-hcs.org.
| S-EPMC5058786 | biostudies-literature
Project description:The rates of cell growth, division, and carbon loss of microbial populations are key parameters for understanding how organisms interact with their environment and how they contribute to the carbon cycle. However, the invasive nature of current analytical methods has hindered efforts to reliably quantify these parameters. In recent years, size-structured matrix population models (MPMs) have gained popularity for estimating division rates of microbial populations by mechanistically describing changes in microbial cell size distributions over time. Motivated by the mechanistic structure of these models, we employ a Bayesian approach to extend size-structured MPMs to capture additional biological processes describing the dynamics of a marine phytoplankton population over the day-night cycle. Our Bayesian framework is able to take prior scientific knowledge into account and generate biologically interpretable results. Using data from an exponentially growing laboratory culture of the cyanobacterium Prochlorococcus, we isolate respiratory and exudative carbon losses as critical parameters for the modeling of their population dynamics. The results suggest that this modeling framework can provide deeper insights into microbial population dynamics provided by size distribution time-series data.
| S-EPMC8794270 | biostudies-literature
Project description:Background. Establishing health-related causal relationships is a central pursuit in biomedical research. Yet, the interdependent non-linearity of biological systems renders causal dynamics laborious and at times impractical to disentangle. This pursuit is further impeded by the dearth of time series that are sufficiently long to observe and understand recurrent patterns of flux. However, as data generation costs plummet and technologies like wearable devices democratize data collection, we anticipate a coming surge in the availability of biomedically-relevant time series data. Given the life-saving potential of these burgeoning resources, it is critical to invest in the development of open source software tools that are capable of drawing meaningful insight from vast amounts of time series data. Results. Here we present CauseMap, the first open source implementation of convergent cross mapping (CCM), a method for establishing causality from long time series data (?25 observations). Compared to existing time series methods, CCM has the advantage of being model-free and robust to unmeasured confounding that could otherwise induce spurious associations. CCM builds on Takens' Theorem, a well-established result from dynamical systems theory that requires only mild assumptions. This theorem allows us to reconstruct high dimensional system dynamics using a time series of only a single variable. These reconstructions can be thought of as shadows of the true causal system. If reconstructed shadows can predict points from opposing time series, we can infer that the corresponding variables are providing views of the same causal system, and so are causally related. Unlike traditional metrics, this test can establish the directionality of causation, even in the presence of feedback loops. Furthermore, since CCM can extract causal relationships from times series of, e.g., a single individual, it may be a valuable tool to personalized medicine. We implement CCM in Julia, a high-performance programming language designed for facile technical computing. Our software package, CauseMap, is platform-independent and freely available as an official Julia package. Conclusions. CauseMap is an efficient implementation of a state-of-the-art algorithm for detecting causality from time series data. We believe this tool will be a valuable resource for biomedical research and personalized medicine.
| S-EPMC4359046 | biostudies-literature
Project description:The ability to measure the properties of proteins at the single-molecule level offers an unparalleled glimpse into biological systems at the molecular scale. The interpretation of single-molecule time series has often been rooted in statistical mechanics and the theory of Markov processes. While existing analysis methods have been useful, they are not without significant limitations including problems of model selection and parameter nonidentifiability. To address these challenges, we introduce the use of nonparametric Bayesian inference for the analysis of single-molecule time series. These methods provide a flexible way to extract structure from data instead of assuming models beforehand. We demonstrate these methods with applications to several diverse settings in single-molecule biophysics. This approach provides a well-constrained and rigorously grounded method for determining the number of biophysical states underlying single-molecule data.
| S-EPMC4317543 | biostudies-literature
Project description:BackgroundComparing metabolic profiles under different biological perturbations has become a powerful approach to investigating the functioning of cells. The profiles can be taken as single snapshots of a system, but more information is gained if they are measured longitudinally over time. The results are short time series consisting of relatively sparse data that cannot be analyzed effectively with standard time series techniques, such as autocorrelation and frequency domain methods. In this work, we study longitudinal time series profiles of glucose consumption in the yeast Saccharomyces cerevisiae under different temperatures and preconditioning regimens, which we obtained with methods of in vivo nuclear magnetic resonance (NMR) spectroscopy. For the statistical analysis we first fit several nonlinear mixed effect regression models to the longitudinal profiles and then used an ANOVA likelihood ratio method in order to test for significant differences between the profiles.ResultsThe proposed methods are capable of distinguishing metabolic time trends resulting from different treatments and associate significance levels to these differences. Among several nonlinear mixed-effects regression models tested, a three-parameter logistic function represents the data with highest accuracy. ANOVA and likelihood ratio tests suggest that there are significant differences between the glucose consumption rate profiles for cells that had been--or had not been--preconditioned by heat during growth. Furthermore, pair-wise t-tests reveal significant differences in the longitudinal profiles for glucose consumption rates between optimal conditions and heat stress, optimal and recovery conditions, and heat stress and recovery conditions (p-values <0.0001).ConclusionWe have developed a nonlinear mixed effects model that is appropriate for the analysis of sparse metabolic and physiological time profiles. The model permits sound statistical inference procedures, based on ANOVA likelihood ratio tests, for testing the significance of differences between short time course data under different biological perturbations.
| S-EPMC3114728 | biostudies-literature
Project description:Understanding the nature of interpopulation interactions in host-associated microbial communities is critical to understanding gut colonization, responses to perturbations, and transitions between health and disease. Characterizing these interactions is complicated by the complexity of these communities and the observation that even if populations can be cultured, their in vitro and in vivo phenotypes differ significantly. Dynamic models are the cornerstone of computational systems biology and a key objective of computational systems biologists is the reconstruction of biological networks (i.e., network inference) from high-throughput data. When such computational models reflect biology, they provide an opportunity to generate testable hypotheses as well as to perform experiments that are impractical or not feasible in vivo or in vitro. We modeled time-series data for murine microbial communities using statistical approaches and systems of ordinary differential equations. To obtain the dense time-series data, we sequenced the 16S ribosomal RNA (rRNA) gene from DNA isolated from the fecal material of germfree mice colonized with cecal contents of conventionally raised animals. The modeling results suggested a lack of mutualistic interactions within the community. Among the members of the Bacteroidetes, there was evidence for closely related pairs of populations to exhibit parasitic interactions. Among the Firmicutes, the interactions were all competitive. These results suggest future animal and in silico experiments. Our modeling approach can be applied to other systems to provide a greater understanding of the dynamics of communities associated with health and disease.
| S-EPMC3890833 | biostudies-other
Project description:BackgroundGrowth rates, interactions between community members, stochasticity, and immigration are important drivers of microbial community dynamics. In sequencing data analysis, such as network construction and community model parameterization, we make implicit assumptions about the nature of these drivers and thereby restrict model outcome. Despite apparent risk of methodological bias, the validity of the assumptions is rarely tested, as comprehensive procedures are lacking. Here, we propose a classification scheme to determine the processes that gave rise to the observed time series and to enable better model selection.ResultsWe implemented a three-step classification scheme in R that first determines whether dependence between successive time steps (temporal structure) is present in the time series and then assesses with a recently developed neutrality test whether interactions between species are required for the dynamics. If the first and second tests confirm the presence of temporal structure and interactions, then parameters for interaction models are estimated. To quantify the importance of temporal structure, we compute the noise-type profile of the community, which ranges from black in case of strong dependency to white in the absence of any dependency. We applied this scheme to simulated time series generated with the Dirichlet-multinomial (DM) distribution, Hubbell's neutral model, the generalized Lotka-Volterra model and its discrete variant (the Ricker model), and a self-organized instability model, as well as to human stool microbiota time series. The noise-type profiles for all but DM data clearly indicated distinctive structures. The neutrality test correctly classified all but DM and neutral time series as non-neutral. The procedure reliably identified time series for which interaction inference was suitable. Both tests were required, as we demonstrated that all structured time series, including those generated with the neutral model, achieved a moderate to high goodness of fit to the Ricker model.ConclusionsWe present a fast and robust scheme to classify community structure and to assess the prevalence of interactions directly from microbial time series data. The procedure not only serves to determine ecological drivers of microbial dynamics, but also to guide selection of appropriate community models for prediction and follow-up analysis.
| S-EPMC6022718 | biostudies-literature