BEEM-Static: Accurate inference of ecological interactions from cross-sectional microbiome data.
Ontology highlight
ABSTRACT: Conclusion
BEEM-Static provides new opportunities for mining ecologically interpretable interactions and systems insights from the growing corpus of microbiome data.
SUBMITTER: Li C
PROVIDER: S-EPMC8452072 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:BackgroundThe dynamics of microbial communities is driven by a range of interactions from symbiosis to predator-prey relationships, the majority of which are poorly understood. With the increasing availability of high-throughput microbiome taxonomic profiling data, it is now conceivable to directly learn the ecological models that explicitly define microbial interactions and explain community dynamics. The applicability of these approaches is severely limited by the lack of accurate absolute cell density measurements (biomass).MethodsWe present a new computational approach that resolves this key limitation in the inference of generalized Lotka-Volterra models (gLVMs) by coupling biomass estimation and model inference with an expectation-maximization algorithm (BEEM).ResultsBEEM outperforms the state-of-the-art methods for inferring gLVMs, while simultaneously eliminating the need for additional experimental biomass data as input. BEEM's application to previously inaccessible public datasets (due to the lack of biomass data) allowed us to construct ecological models of microbial communities in the human gut on a per-individual basis, revealing personalized dynamics and keystone species.ConclusionsBEEM addresses a key bottleneck in "systems analysis" of microbiomes by enabling accurate inference of ecological models from high throughput sequencing data without the need for experimental biomass measurements.
| S-EPMC6706891 | biostudies-literature
Project description:Background: Hypertension is one of the most common metabolic diseases in the elderly and its pathogenesis is associated with microbiota dysbiosis. Recent evidence suggests that oral microbiota dysbiosis is also an important factor in the development of hypertension. However, the relationship between hypertension and oral flora in the elderly has not been adequately investigated. Objective: The aim of this cross-sectional study was to investigate the structure of the oral microbiota and its correlation with hypertension in elderly hypertensive patients. To provide new ideas for the prevention and treatment of hypertension. Methods: 206 subjects aged 60 ~ 89 years were selected and divided into normal (CON) and hypertensive (HTN) groups, according to the 2018 Chinese Guidelines for the Management of Hypertension. The oral microbiome composition of saliva samples was determined by 16S rRNA gene sequencing. Results: Although there was no significant difference in α and β diversity between the two groups, systolic and diastolic blood pressure were the most important factors influencing the structure of the oral microbiota. At the phylum level, the relative abundance of the spirochete phylum and the mutualistic bacterial phylum was higher in the HT group than in the CON group (p < 0.05). Diastolic blood pressure was negatively correlated with Streptococcus. Furthermore, we analyzed HTN patients with 120 mmHg<systolic blood pressure<160 mmHg and systolic blood pressure>160 mmHg separately and found that the abundance of Saccharibacteria_(TM7) was significantly increased in the HTN_2 group. Conclusions: Our study identified specific oral microbiota in elderly hypertensive patients, confirming the relationship between oral microbiota and hypertension. This enhances our understanding of the important role of oral microbiota in the pathogenesis of hypertension and accumulates more evidence for microbial involvement in the development of hypertension.
2024-09-17 | GSE277296 | GEO
Project description:BackgroundOne of the main challenges of microbiome analysis is its compositional nature that if ignored can lead to spurious results. Addressing the compositional structure of microbiome data is particularly critical in longitudinal studies where abundances measured at different times can correspond to different sub-compositions.ResultsWe developed coda4microbiome, a new R package for analyzing microbiome data within the Compositional Data Analysis (CoDA) framework in both, cross-sectional and longitudinal studies. The aim of coda4microbiome is prediction, more specifically, the method is designed to identify a model (microbial signature) containing the minimum number of features with the maximum predictive power. The algorithm relies on the analysis of log-ratios between pairs of components and variable selection is addressed through penalized regression on the "all-pairs log-ratio model", the model containing all possible pairwise log-ratios. For longitudinal data, the algorithm infers dynamic microbial signatures by performing penalized regression over the summary of the log-ratio trajectories (the area under these trajectories). In both, cross-sectional and longitudinal studies, the inferred microbial signature is expressed as the (weighted) balance between two groups of taxa, those that contribute positively to the microbial signature and those that contribute negatively. The package provides several graphical representations that facilitate the interpretation of the analysis and the identified microbial signatures. We illustrate the new method with data from a Crohn's disease study (cross-sectional data) and on the developing microbiome of infants (longitudinal data).Conclusionscoda4microbiome is a new algorithm for identification of microbial signatures in both, cross-sectional and longitudinal studies. The algorithm is implemented as an R package that is available at CRAN ( https://cran.r-project.org/web/packages/coda4microbiome/ ) and is accompanied with a vignette with a detailed description of the functions. The website of the project contains several tutorials: https://malucalle.github.io/coda4microbiome/.
| S-EPMC9990256 | biostudies-literature
Project description:Ecological communities are characterized by complex networks of trophic and nontrophic interactions, which shape the dy-namics of the community. Machine learning and correlational methods are increasingly popular for inferring networks from co-occurrence and time series data, particularly in microbial systems. In this study, we test the suitability of these methods for inferring ecological interactions by constructing networks using Dynamic Bayesian Networks, Lasso regression, and Pear-son's correlation coefficient, then comparing the model networks to empirical trophic and nontrophic webs in two ecological systems. We find that although each model significantly replicates the structure of at least one empirical network, no model significantly predicts network structure in both systems, and no model is clearly superior to the others. We also find that networks inferred for the Tatoosh intertidal match the nontrophic network much more closely than the trophic one, possibly due to the challenges of identifying trophic interactions from presence-absence data. Our findings suggest that although these methods hold some promise for ecological network inference, presence-absence data does not provide enough signal for models to consistently identify interactions, and networks inferred from these data should be interpreted with caution.
| S-EPMC5541006 | biostudies-literature
Project description:Abstract Accelerometers in animal‐attached tags are powerful tools in behavioural ecology, they can be used to determine behaviour and provide proxies for movement‐based energy expenditure. Researchers are collecting and archiving data across systems, seasons and device types. However, using data repositories to draw ecological inference requires a good understanding of the error introduced according to sensor type and position on the study animal and protocols for error assessment and minimisation. Using laboratory trials, we examine the absolute accuracy of tri‐axial accelerometers and determine how inaccuracies impact measurements of dynamic body acceleration (DBA), a proxy for energy expenditure, in human participants. We then examine how tag type and placement affect the acceleration signal in birds, using pigeons Columba livia flying in a wind tunnel, with tags mounted simultaneously in two positions, and back‐ and tail‐mounted tags deployed on wild kittiwakes Rissa tridactyla. Finally, we present a case study where two generations of tag were deployed using different attachment procedures on red‐tailed tropicbirds Phaethon rubricauda foraging in different seasons. Bench tests showed that individual acceleration axes required a two‐level correction to eliminate measurement error. This resulted in DBA differences of up to 5% between calibrated and uncalibrated tags for humans walking at a range of speeds. Device position was associated with greater variation in DBA, with upper and lower back‐mounted tags varying by 9% in pigeons, and tail‐ and back‐mounted tags varying by 13% in kittiwakes. The tropicbird study highlighted the difficulties of attributing changes in signal amplitude to a single factor when confounding influences tend to covary, as DBA varied by 25% between seasons. Accelerometer accuracy, tag placement and attachment critically affect the signal amplitude and thereby the ability of the system to detect biologically meaningful phenomena. We propose a simple method to calibrate accelerometers that can be executed under field conditions. This should be used prior to deployments and archived with resulting data. We also suggest a way that researchers can assess accuracy in previously collected data, and caution that variable tag placement and attachment can increase sensor noise and even generate trends that have no biological meaning.
| S-EPMC9303593 | biostudies-literature
Project description:BACKGROUND:Several studies have focused on the microbiota living in environmental niches including human body sites. In many of these studies, researchers collect longitudinal data with the goal of understanding not only just the composition of the microbiome but also the interactions between the different taxa. However, analysis of such data is challenging and very few methods have been developed to reconstruct dynamic models from time series microbiome data. RESULTS:Here, we present a computational pipeline that enables the integration of data across individuals for the reconstruction of such models. Our pipeline starts by aligning the data collected for all individuals. The aligned profiles are then used to learn a dynamic Bayesian network which represents causal relationships between taxa and clinical variables. Testing our methods on three longitudinal microbiome data sets we show that our pipeline improve upon prior methods developed for this task. We also discuss the biological insights provided by the models which include several known and novel interactions. The extended CGBayesNets package is freely available under the MIT Open Source license agreement. The source code and documentation can be downloaded from https://github.com/jlugomar/longitudinal_microbiome_analysis_public . CONCLUSIONS:We propose a computational pipeline for analyzing longitudinal microbiome data. Our results provide evidence that microbiome alignments coupled with dynamic Bayesian networks improve predictive performance over previous methods and enhance our ability to infer biological relationships within the microbiome and between taxa and clinical factors.
| S-EPMC6446388 | biostudies-literature
Project description:Accurate functional annotation of regulatory elements is essential for understanding global gene regulation. Here, we report a genome-wide map of 827,000 transcription factor binding sites in human lymphoblastoid cell lines, which is comprised of sites correspond-ing to 239 position weight matrices of known transcription factor binding motifs, and 49 novel sequence motifs. To generate this map, we developed a probabilistic framework that integrates cell- or tissue-specific experimental data such as histone modifications and DNa-seI cleavage patterns with genomic information such as gene annotation and evolutionary conservation. Comparison to empirical ChIP-seq data suggests that our method is highly accurate yet has the advantage of targeting many factors in a single assay. We anticipate that this approach will be a valuable tool for genome-wide studies of gene regulation in a wide variety of cell-types or tissues under diverse conditions. DNaseI-Seq on two YRI Hapmap cell lines. Each individual sequenced on 8 lanes of the Illumina Genome Analyzer II
2010-11-15 | E-GEOD-25341 | biostudies-arrayexpress
Project description:BACKGROUND:RNA-Seq based transcriptome assembly has become a fundamental technique for studying expressed mRNAs (i.e., transcripts or isoforms) in a cell using high-throughput sequencing technologies, and is serving as a basis to analyze the structural and quantitative differences of expressed isoforms between samples. However, the current transcriptome assembly algorithms are not specifically designed to handle large amounts of errors that are inherent in real RNA-Seq datasets, especially those involving multiple samples, making downstream differential analysis applications difficult. On the other hand, multiple sample RNA-Seq datasets may provide more information than single sample datasets that can be utilized to improve the performance of transcriptome assembly and abundance estimation, but such information remains overlooked by the existing assembly tools. RESULTS:We formulate a computational framework of transcriptome assembly that is capable of handling noisy RNA-Seq reads and multiple sample RNA-Seq datasets efficiently. We show that finding an optimal solution under this framework is an NP-hard problem. Instead, we develop an efficient heuristic algorithm, called Iterative Shortest Path (ISP), based on linear programming (LP) and integer linear programming (ILP). Our preliminary experimental results on both simulated and real datasets and comparison with the existing assembly tools demonstrate that (i) the ISP algorithm is able to assemble transcriptomes with a greatly increased precision while keeping the same level of sensitivity, especially when many samples are involved, and (ii) its assembly results help improve downstream differential analysis. The source code of ISP is freely available at http://alumni.cs.ucr.edu/~liw/isp.html.
| S-EPMC4331715 | biostudies-literature
Project description:Ecological memory cross-reactivity sequence data
| PRJEB48837 | ENA
Project description:BackgroundEstimating the parameters that describe the ecology of viruses,particularly those that are novel, can be made possible using metagenomic approaches. However, the best-performing existing methods require databases to first estimate an average genome length of a viral community before being able to estimate other parameters, such as viral richness. Although this approach has been widely used, it can adversely skew results since the majority of viruses are yet to be catalogued in databases.ResultsIn this paper, we present ENVirT, a method for estimating the richness of novel viral mixtures, and for the first time we also show that it is possible to simultaneously estimate the average genome length without a priori information. This is shown to be a significant improvement over database-dependent methods, since we can now robustly analyze samples that may include novel viral types under-represented in current databases. We demonstrate that the viral richness estimates produced by ENVirT are several orders of magnitude higher in accuracy than the estimates produced by existing methods named PHACCS and CatchAll when benchmarked against simulated data. We repeated the analysis of 20 metavirome samples using ENVirT, which produced results in close agreement with complementary in virto analyses.ConclusionsThese insights were previously not captured by existing computational methods. As such, ENVirT is shown to be an essential tool for enhancing our understanding of novel viral populations.
| S-EPMC7394321 | biostudies-literature