Longitudinal Prediction of the Infant Gut Microbiome with Dynamic Bayesian Networks

Methodology for microbiome data analysis: An overview

It is known that microbiome and health are related, in addition, recent research has found that microbiome has potential clinical uses. These facts highlight the importance of the microbiome in actual science. However, microbiome data has some characteristics that makes its statistical study challenging. In recent years, longitudinal and non-longitudinal methods have been designed to analyze the microbiota and knowing more about the bacterial behavior. In this article in the form of a review we summarize the characteristics of microbiome data and the statistical methods most widespread to analyze it. We have taken into account if the strategies are longitudinal or not. We also classify the methods based on their specific analytical objectives and based on their mathematical characteristics. The methods are structured according to their biological goals and mathematical features, ensuring that the insights provided are both relevant and accessible to professionals in biology and statistics. We present this review as a reference for the most widely used methods in microbiome data analysis and as a foundation for identifying potential areas for future research. We want to point out that this review can be particularly useful to remark the importance of the methodology designed in order to study microbiome longitudinal datasets.

Bayesian networks for network inference in biology

Bayesian networks (BNs) have been used for reconstructing interactions from biological data, in disciplines ranging from molecular biology to ecology and neuroscience. BNs learn conditional dependencies between variables, which best 'explain' the data, represented as a directed graph which approximates the relationships between variables. In the 2000s, BNs were a popular method that promised an approach capable of inferring biological networks from data. Here, we review the use of BNs applied to biological data over the past two decades and evaluate their efficacy. We find that BNs are successful in inferring biological networks, frequently identifying novel interactions or network components missed by previous analyses. We suggest that as false positive results are underreported, it is difficult to assess the accuracy of BNs in inferring biological networks. BN learning appears most successful for small numbers of variables with high-quality datasets that either discretize the data into few states or include perturbative data. We suggest that BNs have failed to live up to the promise of the 2000s but that this is most likely due to experimental constraints on datasets, and the success of BNs at inferring networks in a variety of biological contexts suggests they are a powerful tool for biologists.

Exploring omics approaches in probiotics: Contemporary developments and prospective pathways

The application of omics technologies in combination with bench investigations has brought about a significant transformation in the field of probiotics, enabling a thorough investigation of the basic elements contributing to the probiotic activity. Genomics studies have decoded the complete set of genes of probiotic organisms, shedding light on beneficial traits and mechanisms of probiotic action. Transcriptomics analyses focus on gene expression patterns and investigate probiotic adaptation and functionality. Proteomic studies have revealed the intricate connections between proteins in probiotic cells and their relationship with the host environment. Metabolomic profiling has provided a comprehensive perspective on the metabolic pathways related to probiotic metabolism and the production of bioactive substances. The ongoing development of omics technology presents exciting opportunities for probiotic research, as it allows for a deeper exploration of probiotic-host interactions and the creation of advanced and tailored probiotics that offer specific health advantages. A comprehensive analysis of recent progress in genomics, transcriptomics, proteomics, and metabolomics related to probiotics is presented in this review.

Exploring the female genital tract mycobiome in young South African women using metaproteomics

BACKGROUND

Female genital tract (FGT) diseases such as bacterial vaginosis (BV) and sexually transmitted infections are prevalent in South Africa, with young women being at an increased risk. Since imbalances in the FGT microbiome are associated with FGT diseases, it is vital to investigate the factors that influence FGT health. The mycobiome plays an important role in regulating mucosal health, especially when the bacterial component is disturbed. However, we have a limited understanding of the FGT mycobiome since many studies have focused on bacterial communities and have neglected low-abundance taxonomic groups, such as fungi. To reduce this knowledge deficit, we present the first large-scale metaproteomic study to define the taxonomic composition and potential functional processes of the FGT mycobiome in South African reproductive-age women.

RESULTS

We examined FGT fungal communities present in 123 women by collecting lateral vaginal wall swabs for liquid chromatography-tandem mass spectrometry. From this, 39 different fungal genera were identified, with Candida dominating the mycobiome (53.2% relative abundance). We observed changes in relative abundance at the protein, genus, and functional (gene ontology biological processes) level between BV states. In women with BV, Malassezia and Conidiobolus proteins were more abundant, while Candida proteins were less abundant compared to BV-negative women. Correspondingly, Nugent scores were negatively associated with total fungal protein abundance. The clinical variables, Nugent score, pro-inflammatory cytokines, chemokines, vaginal pH, Chlamydia trachomatis, and the presence of clue cells were associated with fungal community composition.

CONCLUSIONS

The results of this study revealed the diversity of FGT fungal communities, setting the groundwork for understanding the FGT mycobiome. Video Abstract.

Wise Roles and Future Visionary Endeavors of Current Emperor: Advancing Dynamic Methods for Longitudinal Microbiome Meta-Omics Data in Personalized and Precision Medicine

Understanding the etiological complexity of diseases requires identifying biomarkers longitudinally associated with specific phenotypes. Advanced sequencing tools generate dynamic microbiome data, providing insights into microbial community functions and their impact on health. This review aims to explore the current roles and future visionary endeavors of dynamic methods for integrating longitudinal microbiome multi-omics data in personalized and precision medicine. This work seeks to synthesize existing research, propose best practices, and highlight innovative techniques. The development and application of advanced dynamic methods, including the unified analytical frameworks and deep learning tools in artificial intelligence, are critically examined. Aggregating data on microbes, metabolites, genes, and other entities offers profound insights into the interactions among microorganisms, host physiology, and external stimuli. Despite progress, the absence of gold standards for validating analytical protocols and data resources of various longitudinal multi-omics studies remains a significant challenge. The interdependence of workflow steps critically affects overall outcomes. This work provides a comprehensive roadmap for best practices, addressing current challenges with advanced dynamic methods. The review underscores the biological effects of clinical, experimental, and analytical protocol settings on outcomes. Establishing consensus on dynamic microbiome inter-studies and advancing reliable analytical protocols are pivotal for the future of personalized and precision medicine.

Interpolation of microbiome composition in longitudinal data sets

The human gut microbiome significantly impacts health, prompting a rise in longitudinal studies that capture microbiome samples at multiple time points. Such studies allow researchers to characterize microbiome changes over time, but importantly, also present major analytical challenges due to incomplete or irregular sampling. To address this challenge, longitudinal microbiome studies often employ various interpolation methods, aiming to infer missing microbiome data. However, to date, a comprehensive assessment of such microbiome interpolation techniques, as well as best practice guidelines for interpolating microbiome data, is still lacking. This work aims to fill this gap, rigorously implementing and systematically evaluating a large array of interpolation methods, spanning several different categories, for longitudinal microbiome interpolation. To assess each method and its ability to accurately infer microbiome composition at missing time points, we used three longitudinal microbiome data sets that follow individuals over a long period of time and a leave-one-out approach. Overall, our analysis demonstrated that the K-nearest neighbors algorithm consistently outperforms other methods in interpolation accuracy, yet, accuracy varied widely across data sets, individuals, and time. Factors such as microbiome stability, sample size, and the time gap between interpolated and adjacent samples significantly influenced accuracy, allowing us to develop a model for predicting the expected interpolation accuracy at a missing time point. Our findings, combined, suggest that accurate interpolation in longitudinal microbiome data is feasible, especially in dense cohorts. Furthermore, using our predictive model, future studies can interpolate data only in time points where the expected interpolation accuracy is high.

IMPORTANCE

Since missing samples are common in longitudinal microbiome dataset due to inconsistent collection practices, it is important to evaluate and benchmark different interpolation methods for predicting microbiome composition in such samples and facilitate downstream analysis. Our study rigorously evaluated several such methods and identified the K-nearest neighbors approach as particularly effective for this task. The study also notes significant variability in interpolation accuracy among individuals, influenced by factors such as age, sample size, and sampling frequency. Furthermore, we developed a predictive model for estimating interpolation accuracy at a specific time point, enhancing the reliability of such analyses in future studies. Combined, our study, thus, provides critical insights and tools that enhance the accuracy and reliability of data interpolation methods in the growing field of longitudinal microbiome research.

Prediction of Overall Disease Burden in (y)pN1 Breast Cancer Using Knowledge-Based Machine Learning Model

BACKGROUND

We aimed to construct an expert knowledge-based Bayesian network (BN) model for assessing the overall disease burden (ODB) in (y)pN1 breast cancer patients and compare ODB across arms of ongoing trials.

METHODS

Utilizing institutional data and expert surveys, we developed a BN model for (y)pN1 breast cancer. Expert-derived probabilities and disability weights for radiotherapy-related benefit (e.g., 7-year disease-free survival [DFS]) and toxicities were integrated into the model. ODB was defined as the sum of disability weights multiplied by probabilities. In silico predictions were conducted for Alliance A011202, PORT-N1, RAPCHEM, and RT-CHARM trials, comparing ODB, 7-year DFS, and side effects.

RESULTS

In the Alliance A011202 trial, 7-year DFS was 80.1% in both arms. Axillary lymph node dissection led to higher clinical lymphedema and ODB compared to sentinel lymph node biopsy with full regional nodal irradiation (RNI). In the PORT-N1 trial, the control arm (whole-breast irradiation [WBI] with RNI or post-mastectomy radiotherapy [PMRT]) had an ODB of 0.254, while the experimental arm (WBI alone or no PMRT) had an ODB of 0.255. In the RAPCHEM trial, the radiotherapy field did not impact the 7-year DFS in ypN1 patients. However, there was a mild ODB increase with a larger irradiation field. In the RT-CHARM trial, we identified factors associated with the major complication rate, which ranged from 18.3% to 22.1%.

CONCLUSIONS

The expert knowledge-based BN model predicted ongoing trial outcomes, validating reported results and assumptions. In addition, the model demonstrated the ODB in different arms, with an emphasis on quality of life.

A digital twin of the infant microbiome to predict neurodevelopmental deficits

Despite the recognized gut-brain axis link, natural variations in microbial profiles between patients hinder definition of normal abundance ranges, confounding the impact of dysbiosis on infant neurodevelopment. We infer a digital twin of the infant microbiome, forecasting ecosystem trajectories from a few initial observations. Using 16S ribosomal RNA profiles from 88 preterm infants (398 fecal samples and 32,942 abundance estimates for 91 microbial classes), the model (Q-net) predicts abundance dynamics with R2 = 0.69. Contrasting the fit to Q-nets of typical versus suboptimal development, we can reliably estimate individual deficit risk (Mδ) and identify infants achieving poor future head circumference growth with ≈76% area under the receiver operator characteristic curve, 95% ± 1.8% positive predictive value at 98% specificity at 30 weeks postmenstrual age. We find that early transplantation might mitigate risk for ≈45.2% of the cohort, with potentially negative effects from incorrect supplementation. Q-nets are generative artificial intelligence models for ecosystem dynamics, with broad potential applications.

Modeling Microbial Community Networks: Methods and Tools for Studying Microbial Interactions

Microbial interactions function as a fundamental unit in complex ecosystems. By characterizing the type of interaction (positive, negative, neutral) occurring in these dynamic systems, one can begin to unravel the role played by the microbial species. Towards this, various methods have been developed to decipher the function of the microbial communities. The current review focuses on the various qualitative and quantitative methods that currently exist to study microbial interactions. Qualitative methods such as co-culturing experiments are visualized using microscopy-based techniques and are combined with data obtained from multi-omics technologies (metagenomics, metabolomics, metatranscriptomics). Quantitative methods include the construction of networks and network inference, computational models, and development of synthetic microbial consortia. These methods provide a valuable clue on various roles played by interacting partners, as well as possible solutions to overcome pathogenic microbes that can cause life-threatening infections in susceptible hosts. Studying the microbial interactions will further our understanding of complex less-studied ecosystems and enable design of effective frameworks for treatment of infectious diseases.

Progress on network modeling and analysis of gut microecology: a review

The gut microecological network is a complex microbial community within the human body that plays a key role in linking dietary nutrition and host physiology. To understand the complex relationships among microbes and their functions within this community, network analysis has emerged as a powerful tool. By representing the interactions between microbes and their associated omics data as a network, we can gain a comprehensive understanding of the ecological mechanisms that drive the human gut microbiota. In addition, the network-based approach provides a more intuitive analysis of the gut microbiota, simplifying the study of its complex dynamics and interdependencies. This review provides a comprehensive overview of the methods used to construct and analyze networks in the context of gut microecological background. We discuss various types of network modeling approaches, including co-occurrence networks, causal networks, dynamic networks, and multi-omics networks, and describe the analytical techniques used to identify important network properties. We also highlight the challenges and limitations of network modeling in this area, such as data scarcity and heterogeneity, and provide future research directions to overcome these limitations. By exploring these network-based methods, researchers can gain valuable insights into the intricate relationships and functional roles of microbial communities within the gut, ultimately advancing our understanding of the gut microbiota's impact on human health.

Estimating the risk and benefit of radiation therapy in (y)pN1 stage breast cancer patients: A Bayesian network model incorporating expert knowledge (KROG 22-13)

BACKGROUND

We aimed to evaluate the risk and benefit of (y)pN1 breast cancer patients in a Bayesian network model.

METHOD

We developed a Bayesian network (BN) model comprising three parts: pretreatment, intervention, and risk/benefit. The pretreatment part consisted of clinical information from a tertiary medical center. The intervention part regarded the field of radiotherapy. The risk/benefit component encompasses radiotherapy (RT)-related side effects and effectiveness, including factors such as recurrence, cardiac toxicity, lymphedema, and radiation pneumonitis. These factors were evaluated in terms of disability weights and probabilities from a nationwide expert survey. The overall disease burden (ODB) was calculated as the sum of the probability multiplied by the disability weight. A higher value of ODB indicates a greater disease burden for the patient.

RESULTS

Among the 58 participants, a BN model utilizing discretization and clustering techniques revealed five distinct clusters. Overall, factors associated with breast reconstruction and RT exhibited high discrepancies (24-34 %), while RT-related side effects demonstrated low discrepancies (3-11 %) among the experts. When incorporating recurrence and RT-related side effects, the mean ODB of (y)pN1 patients was 0.258 (range, 0.244-0.337), with a higher tendency observed in triple-negative breast cancer (TNBC) or mastectomy cases. The ODB for TNBC patients undergoing mastectomy without postmastectomy radiotherapy was 0.327, whereas for non-TNBC patients undergoing breast conserving surgery with RT, the disease burden was 0.251. There was an increasing trend in ODB as the field of RT increased.

CONCLUSION

We developed a Bayesian network model based on an expert survey, which helps to understand treatment patterns and enables precise estimations of RT-related risk and benefit in (y)pN1 patients.

Attention Deficit Hyperactivity Disorder (ADHD) and the gut microbiome: An ecological perspective

Attention Deficit Hyperactivity Disorder (ADHD) is an increasingly prevalent neuropsychiatric disorder characterized by hyperactivity, inattention, and impulsivity. Symptoms emerge from underlying deficiencies in neurocircuitry, and recent research has suggested a role played by the gut microbiome. The gut microbiome is an ecosystem of interdependent taxa involved in an exponentially complex web of interactions, plus host gene and reaction pathways, some of which involve neurotransmitters with roles in ADHD neurocircuitry. Studies have analyzed the ADHD gut microbiome using macroscale metrics such as diversity and differential abundance, and have proposed several taxa as elevated or reduced in ADHD compared to Control. Few studies have delved into the complex underlying dynamics ultimately responsible for the emergence of such metrics, leaving a largely incomplete, sometimes contradictory, and ultimately inconclusive picture. We aim to help complete this picture by venturing beyond taxa abundances and into taxa relationships (i.e. cooperation and competition), using a publicly available gut microbiome dataset (targeted 16S, v3-4 region, qPCR) from an observational, case-control study of 30 Control (15 female, 15 male) and 28 ADHD (15 female, 13 male) undergraduate students. We first perform the same macroscale analyses prevalent in ADHD gut microbiome literature (diversity, differential abundance, and composition) to observe the degree of correspondence, or any new trends. We then estimate two-way ecological relationships by producing Control and ADHD Microbial Co-occurrence Networks (MCNs), using SparCC correlations (p ≤ 0.01). We perform community detection to find clusters of taxa estimated to mutually cooperate along with their centroids, and centrality calculations to estimate taxa most vital to overall gut ecology. We finally summarize our results, providing conjectures on how they can guide future experiments, some methods for improving our experiments, and general implications for the field.

Capturing the dynamics of microbial interactions through individual-specific networks

Longitudinal analysis of multivariate individual-specific microbiome profiles over time or across conditions remains dauntin. Most statistical tools and methods that are available to study microbiomes are based on cross-sectional data. Over the past few years, several attempts have been made to model the dynamics of bacterial species over time or across conditions. However, the field needs novel views on handling microbial interactions in temporal analyses. This study proposes a novel data analysis framework, MNDA, that combines representation learning and individual-specific microbial co-occurrence networks to uncover taxon neighborhood dynamics. As a use case, we consider a cohort of newborns with microbiomes available at 6 and 9 months after birth, and extraneous data available on the mode of delivery and diet changes between the considered time points. Our results show that prediction models for these extraneous outcomes based on an MNDA measure of local neighborhood dynamics for each taxon outperform traditional prediction models solely based on individual-specific microbial abundances. Furthermore, our results show that unsupervised similarity analysis of newborns in the study, again using the notion of a taxon's dynamic neighborhood derived from time-matched individual-specific microbial networks, can reveal different subpopulations of individuals, compared to standard microbiome-based clustering, with potential relevance to clinical practice. This study highlights the complementarity of microbial interactions and abundances in downstream analyses and opens new avenues to personalized prediction or stratified medicine with temporal microbiome data.

Univariate and Multivariate Statistical Analysis of Microbiome Data: An Overview

Microbiome data is high dimensional, sparse, compositional, and over-dispersed. Therefore, modeling microbiome data is very challenging and it is an active research area. Microbiome analysis has become a progressing area of research as microorganisms constitute a large part of life. Since many methods of microbiome data analysis have been presented, this review summarizes the challenges, methods used, and the advantages and disadvantages of those methods, to serve as an updated guide for those in the field. This review also compared different methods of analysis to progress the development of newer methods.

Discovering gene regulatory networks of multiple phenotypic groups using dynamic Bayesian networks

Dynamic Bayesian networks (DBNs) can be used for the discovery of gene regulatory networks (GRNs) from time series gene expression data. Here, we suggest a strategy for learning DBNs from gene expression data by employing a Bayesian approach that is scalable to large networks and is targeted at learning models with high predictive accuracy. Our framework can be used to learn DBNs for multiple groups of samples and highlight differences and similarities in their GRNs. We learn these DBN models based on different structural and parametric assumptions and select the optimal model based on the cross-validated predictive accuracy. We show in simulation studies that our approach is better equipped to prevent overfitting than techniques used in previous studies. We applied the proposed DBN-based approach to two time series transcriptomic datasets from the Gene Expression Omnibus database, each comprising data from distinct phenotypic groups of the same tissue type. In the first case, we used DBNs to characterize responders and non-responders to anti-cancer therapy. In the second case, we compared normal to tumor cells of colorectal tissue. The classification accuracy reached by the DBN-based classifier for both datasets was higher than reported previously. For the colorectal cancer dataset, our analysis suggested that GRNs for cancer and normal tissues have a lot of differences, which are most pronounced in the neighborhoods of oncogenes and known cancer tissue markers. The identified differences in gene networks of cancer and normal cells may be used for the discovery of targeted therapies.

Temporal Alignment of Longitudinal Microbiome Data

A major challenge in working with longitudinal data when studying some temporal process is the fact that differences in pace and dynamics might overshadow similarities between processes. In the case of longitudinal microbiome data, this may hinder efforts to characterize common temporal trends across individuals or to harness temporal information to better understand the link between the microbiome and the host. One possible solution to this challenge lies in the field of “temporal alignment” – an approach for optimally aligning longitudinal samples obtained from processes that may vary in pace. In this work we investigate the use of alignment-based analysis in the microbiome domain, focusing on microbiome data from infants in their first years of life. Our analyses center around two main use-cases: First, using the overall alignment score as a measure of the similarity between microbiome developmental trajectories, and showing that this measure can capture biological differences between individuals. Second, using the specific matching obtained between pairs of samples in the alignment to highlight changes in pace and temporal dynamics, showing that it can be utilized to predict the age of infants based on their microbiome and to uncover developmental delays. Combined, our findings serve as a proof-of-concept for the use of temporal alignment as an important and beneficial tool in future longitudinal microbiome studies.

Dynamic networks of psychological symptoms, impairment, substance use, and social support: The evolution of psychopathology among emerging adults

BACKGROUND

Subthreshold/attenuated syndromes are established precursors of full-threshold mood and psychotic disorders. Less is known about the individual symptoms that may precede the development of subthreshold syndromes and associated social/functional outcomes among emerging adults.

METHODS

We modeled two dynamic Bayesian networks (DBN) to investigate associations among self-rated phenomenology and personal/lifestyle factors (role impairment, low social support, and alcohol and substance use) across the 19Up and 25Up waves of the Brisbane Longitudinal Twin Study. We examined whether symptoms and personal/lifestyle factors at 19Up were associated with (a) themselves or different items at 25Up, and (b) onset of a depression-like, hypo-manic-like, or psychotic-like subthreshold syndrome (STS) at 25Up.

RESULTS

The first DBN identified 11 items that when endorsed at 19Up were more likely to be reendorsed at 25Up (e.g., hypersomnia, impaired concentration, impaired sleep quality) and seven items that when endorsed at 19Up were associated with different items being endorsed at 25Up (e.g., earlier fatigue and later role impairment; earlier anergia and later somatic pain). In the second DBN, no arcs met our a priori threshold for inclusion. In an exploratory model with no threshold, >20 items at 19Up were associated with progression to an STS at 25Up (with lower statistical confidence); the top five arcs were: feeling threatened by others and a later psychotic-like STS; increased activity and a later hypo-manic-like STS; and anergia, impaired sleep quality, and/or hypersomnia and a later depression-like STS.

CONCLUSIONS

These probabilistic models identify symptoms and personal/lifestyle factors that might prove useful targets for indicated preventative strategies.

Explainable Machine Learning for Longitudinal Multi-Omic Microbiome

Over the years, research studies have shown there is a key connection between the microbial community in the gut, genes, and immune system. Understanding this association may help discover the cause of complex chronic idiopathic disorders such as inflammatory bowel disease. Even though important efforts have been put into the field, the functions, dynamics, and causation of dysbiosis state performed by the microbial community remains unclear. Machine learning models can help elucidate important connections and relationships between microbes in the human host. Our study aims to extend the current knowledge of associations between the human microbiome and health and disease through the application of dynamic Bayesian networks to describe the temporal variation of the gut microbiota and dynamic relationships between taxonomic entities and clinical variables. We develop a set of preprocessing steps to clean, filter, select, integrate, and model informative metagenomics, metatranscriptomics, and metabolomics longitudinal data from the Human Microbiome Project. This study accomplishes novel network models with satisfactory predictive performance (accuracy = 0.648) for each inflammatory bowel disease state, validating Bayesian networks as a framework for developing interpretable models to help understand the basic ways the different biological entities (taxa, genes, metabolites) interact with each other in a given environment (human gut) over time. These findings can serve as a starting point to advance the discovery of novel therapeutic approaches and new biomarkers for precision medicine.

Quantifying the impact of ecological memory on the dynamics of interacting communities

Ecological memory refers to the influence of past events on the response of an ecosystem to exogenous or endogenous changes. Memory has been widely recognized as a key contributor to the dynamics of ecosystems and other complex systems, yet quantitative community models often ignore memory and its implications.

Recent modeling studies have shown how interactions between community members can lead to the emergence of resilience and multistability under environmental perturbations. We demonstrate how memory can be introduced in such models using the framework of fractional calculus. We study how the dynamics of a well-characterized interaction model is affected by gradual increases in ecological memory under varying initial conditions, perturbations, and stochasticity.

Our results highlight the implications of memory on several key aspects of community dynamics. In general, memory introduces inertia into the dynamics. This favors species coexistence under perturbation, enhances system resistance to state shifts, mitigates hysteresis, and can affect system resilience both ways depending on the time scale considered. Memory also promotes long transient dynamics, such as long-standing oscillations and delayed regime shifts, and contributes to the emergence and persistence of alternative stable states. Our study highlights the fundamental role of memory in communities, and provides quantitative tools to introduce it in ecological models and analyse its impact under varying conditions.

An ecosystem is said to exhibit ecological memory when its future states do not only depend on its current state but also on its initial state and trajectory. Memory may arise through various mechanisms as organisms adapt to their environment, modify it, and accumulate biotic and abiotic material. It may also emerge from phenotypic heterogeneity at the population level. Despite its commonness in nature, ecological memory and its potential influence on ecosystem dynamics have been so far overlooked in many applied contexts. Here, we use modeling to investigate how memory can influence the dynamics, composition, and stability landscape of communities. We incorporate long-term memory effects into a multi-species model recently introduced to investigate alternative stable states in microbial communities. We assess the impact of memory on key aspects of model behavior and further examine our findings using a model parameterized by empirical data from the human gut microbiota. Our approach for modeling long-term memory and studying its implications has the potential to improve our understanding of microbial community dynamics and ultimately our ability to predict, manipulate, and experimentally design microbial ecosystems. It could also be applied more broadly in the study of systems composed of interacting components.

Longitudinal 16S rRNA gut microbiota data of infant triplets show partial susceptibility to host genetics

The question of whether host genetics plays a role in the development of the infant gut microbiota does not, as yet, have a clear answer. In order to throw additional light on this question, we have analyzed 16S rRNA amplicon sequences from 99 valid fecal samples of five sets of dichorionic triplet babies born by C-section from 1 to 36 months of age. Beta diversity analysis showed that monozygotic twins were more similar to each other than their dizygotic siblings. Monozygotic twins also tended to share more amplicon sequence variants between them. Heritability analysis showed that the genera Bacteroides and Veillonella are particularly susceptible to host genetics. We conclude that infant gut microbiota development is influenced by host genetics, but this effect is subtle and may affect only certain bacterial taxa during a limited time period early in life.

ARZIMM: A Novel Analytic Platform for the Inference of Microbial Interactions and Community Stability from Longitudinal Microbiome Study

Dynamic changes of microbiome communities may play important roles in human health and diseases. The recent rise in longitudinal microbiome studies calls for statistical methods that can model the temporal dynamic patterns and simultaneously quantify the microbial interactions and community stability. Here, we propose a novel autoregressive zero-inflated mixed-effects model (ARZIMM) to capture the sparse microbial interactions and estimate the community stability. ARZIMM employs a zero-inflated Poisson autoregressive model to model the excessive zero abundances and the non-zero abundances separately, a random effect to investigate the underlining dynamic pattern shared within the group, and a Lasso-type penalty to capture and estimate the sparse microbial interactions. Based on the estimated microbial interaction matrix, we further derive the estimate of community stability, and identify the core dynamic patterns through network inference. Through extensive simulation studies and real data analyses we evaluate ARZIMM in comparison with the other methods.

Conceptual strategies for characterizing interactions in microbial communities

Understanding the sets of inter- and intraspecies interactions in microbial communities is a fundamental goal of microbial ecology. However, the study and quantification of microbial interactions pose several challenges owing to their complexity, dynamic nature, and the sheer number of unique interactions within a typical community. To overcome such challenges, microbial ecologists must rely on various approaches to distill the system of study to a functional and conceptualizable level, allowing for a practical understanding of microbial interactions in both simplified and complex systems. This review broadly addresses the role of several conceptual approaches available for the microbial ecologist’s arsenal, examines specific tools used to accomplish such approaches, and describes how the assumptions, expectations, and philosophies underlying these tools change across scales of complexity.

Resource Competition and Host Feedbacks Underlie Regime Shifts in Gut Microbiota

AbstractThe spread of an enteric pathogen in the human gut depends on many interacting factors, including pathogen exposure, diet, host gut environment, and host microbiota, but how these factors jointly influence infection outcomes remains poorly characterized. Here we develop a model of host-mediated resource competition between mutualistic and pathogenic taxa in the gut that aims to explain why similar hosts, exposed to the same pathogen, can have such different infection outcomes. Our model successfully reproduces several empirically observed phenomena related to transitions between healthy and infected states, including (1) the nonlinear relationship between pathogen inoculum size and infection persistence, (2) the elevated risk of chronic infection during or after treatment with broad-spectrum antibiotics, (3) the resolution of gut dysbiosis with fecal microbiota transplants, and (4) the potential protection from infection conferred by probiotics. We then use the model to explore how host-mediated interventions-namely, shifts in the supply rates of electron donors (e.g., dietary fiber) and respiratory electron acceptors (e.g., oxygen)-can potentially be used to direct gut community assembly. Our study demonstrates how resource competition and ecological feedbacks between the host and the gut microbiota can be critical determinants of human health outcomes. We identify several testable model predictions ready for experimental validation.

Network analysis methods for studying microbial communities: A mini review

Microorganisms including bacteria, fungi, viruses, protists and archaea live as communities in complex and contiguous environments. They engage in numerous inter- and intra- kingdom interactions which can be inferred from microbiome profiling data. In particular, network-based approaches have proven helpful in deciphering complex microbial interaction patterns. Here we give an overview of state-of-the-art methods to infer intra-kingdom interactions ranging from simple correlation- to complex conditional dependence-based methods. We highlight common biases encountered in microbial profiles and discuss mitigation strategies employed by different tools and their trade-off with increased computational complexity. Finally, we discuss current limitations that motivate further method development to infer inter-kingdom interactions and to robustly and comprehensively characterize microbial environments in the future.

Dynamic Bayesian Networks for Integrating Multi-omics Time Series Microbiome Data

A key challenge in the analysis of longitudinal microbiome data is the inference of temporal interactions between microbial taxa, their genes, the metabolites that they consume and produce, and host genes. To address these challenges, we developed a computational pipeline, a pipeline for the analysis of longitudinal multi-omics data (PALM), that first aligns multi-omics data and then uses dynamic Bayesian networks (DBNs) to reconstruct a unified model. Our approach overcomes differences in sampling and progression rates, utilizes a biologically inspired multi-omic framework, reduces the large number of entities and parameters in the DBNs, and validates the learned network. Applying PALM to data collected from inflammatory bowel disease patients, we show that it accurately identifies known and novel interactions. Targeted experimental validations further support a number of the predicted novel metabolite-taxon interactions.

IMPORTANCE While a number of large consortia collect and profile several different types of microbiome and genomic time series data, very few methods exist for joint modeling of multi-omics data sets. We developed a new computational pipeline, PALM, which uses dynamic Bayesian networks (DBNs) and is designed to integrate multi-omics data from longitudinal microbiome studies. When used to integrate sequence, expression, and metabolomics data from microbiome samples along with host expression data, the resulting models identify interactions between taxa, their genes, and the metabolites that they produce and consume, as well as their impact on host expression. We tested the models both by using them to predict future changes in microbiome levels and by comparing the learned interactions to known interactions in the literature. Finally, we performed experimental validations for a few of the predicted interactions to demonstrate the ability of the method to identify novel relationships and their impact.

CODY enables quantitatively spatiotemporal predictions on in vivo gut microbial variability induced by diet intervention

Quantitatively understanding and predicting spatiotemporal dynamics of microbiota is imperative for development of tailored microbiome-directed therapeutics treatments. However, the complexity of microbial variations, due to interactions with the host, other microbes, and environmental factors, makes it challenging to identify how microbiota colonize in the human gut. Here, we describe a novel multiscale framework for COmputing the DYnamics of the gut microbiota (CODY), which enables the quantification of spatiotemporal-specific variations of gut microbiome abundance profiles, without a prior knowledge of microbiome interactions. Importantly, the predictive power of CODY is demonstrated using cross-sectional data from two longitudinal metagenomics studies—the microbiota development during early infancy and during short-term diet intervention of obese adults.

Microbial variations in the human gut are harbored in temporal and spatial heterogeneity, and quantitative prediction of spatiotemporal dynamic changes in the gut microbiota is imperative for development of tailored microbiome-directed therapeutics treatments, e.g. precision nutrition. Given the high-degree complexity of microbial variations, subject to the dynamic interactions among host, microbial, and environmental factors, identifying how microbiota colonize in the gut represents an important challenge. Here we present COmputing the DYnamics of microbiota (CODY), a multiscale framework that integrates species-level modeling of microbial dynamics and ecosystem-level interactions into a mathematical model that characterizes spatial-specific in vivo microbial residence in the colon as impacted by host physiology. The framework quantifies spatiotemporal resolution of microbial variations on species-level abundance profiles across site-specific colon regions and in feces, independent of a priori knowledge. We demonstrated the effectiveness of CODY using cross-sectional data from two longitudinal metagenomics studies—the microbiota development during early infancy and during short-term diet intervention of obese adults. For each cohort, CODY correctly predicts the microbial variations in response to diet intervention, as validated by available metagenomics and metabolomics data. Model simulations provide insight into the biogeographical heterogeneity among lumen, mucus, and feces, which provides insight into how host physical forces and spatial structure are shaping microbial structure and functionality.

Inferring directional relationships in microbial communities using signed Bayesian networks

BACKGROUND

Microbe-microbe and host-microbe interactions in a microbiome play a vital role in both health and disease. However, the structure of the microbial community and the colonization patterns are highly complex to infer even under controlled wet laboratory conditions. In this study, we investigate what information, if any, can be provided by a Bayesian Network (BN) about a microbial community. Unlike the previously proposed Co-occurrence Networks (CoNs), BNs are based on conditional dependencies and can help in revealing complex associations.

RESULTS

In this paper, we propose a way of combining a BN and a CoN to construct a signed Bayesian Network (sBN). We report a surprising association between directed edges in signed BNs and known colonization orders.

CONCLUSIONS

BNs are powerful tools for community analysis and extracting influences and colonization patterns, even though the analysis only uses an abundance matrix with no temporal information. We conclude that directed edges in sBNs when combined with negative correlations are consistent with and strongly suggestive of colonization order.

Inter-site and interpersonal diversity of salivary and tongue microbiomes, and the effect of oral care tablets

Background: Oral microbiota has been linked to both health and disease. Specifically, tongue-coating microbiota has been implicated in aspiration pneumonia and halitosis. Approaches altering one's oral microbiota have the potential to improve oral health and prevent diseases. Methods: Here, we designed a study that allows simultaneous monitoring of the salivary and tongue microbiomes during an intervention on the oral microbiota. We applied this study design to evaluate the effect of single-day use of oral care tablets on the oral microbiome of 10 healthy individuals. Tablets with or without actinidin, a protease that reduces biofilm formation in vitro, were tested. Results: Alpha diversity in the saliva was higher than that on the tongue without the intervention. The core operational taxonomic units (OTUs) common to both sites were identified. The salivary and tongue microbiomes of one individual tended to be more similar to one another than to those of other individuals. The tablets did not affect the alpha or beta diversity of the oral microbiome, nor the abundance of specific bacterial species. Conclusions: While the salivary and tongue microbiomes differ significantly in terms of bacterial composition, they show inter- rather than intra-individual diversity. A one-day usage of oral care tablets did not alter the salivary or tongue microbiomes of healthy adults. Whether the use of oral tablets for a longer period on healthy people or people with greater tongue coating accumulation shifts their oral microbiome needs to be investigated.

Inter-site and interpersonal diversity of salivary and tongue microbiomes, and the effect of oral care tablets

Background: Oral microbiota has been linked to both health and diseases. Specifically, tongue-coating microbiota has been implicated in aspiration pneumonia and halitosis. Approaches altering one's oral microbiota have the potential to improve oral health and prevent diseases. Methods: Here, we designed a study that allows simultaneous monitoring of the salivary and tongue microbiomes during an intervention on the oral microbiota. We applied this study design to evaluate the effect of single-day use of oral care tablets on the oral microbiome of 10 healthy individuals. Tablets with or without actinidin, a protease that reduces biofilm formation in vitro, were tested. Results: Alpha diversity of the tongue microbiome was significantly lower than that of the salivary microbiome, using both the number of observed amplicon sequence variants (254 ± 53 in saliva and 175 ± 37 in tongue; P = 8.9e-7, Kruskal-Wallis test) and Shannon index (6.0 ± 0.4 in saliva and 5.4 ± 0.3 in tongue; P = 2.0e-7, Kruskal-Wallis test). Fusobacterium periodonticum, Saccharibacteria sp. 352, Streptococcus oralis subsp . dentisani, Prevotella melaninogenica, Granulicatella adiacens, Campylobacter concisus, and Haemophilus parainfluenzae were the core operational taxonomic units (OTUs) common to both sites. The salivary and tongue microbiomes of one individual tended to be more similar to one another than to those of other individuals. The tablets did not affect the alpha or beta diversity of the oral microbiome, nor the abundance of specific bacterial species. Conclusions: While the salivary and tongue microbiomes differed significantly in terms of bacterial composition, they showed inter- rather than intra-individual diversity. A one-day usage of oral care tablets did not alter the salivary or tongue microbiomes of healthy adults. Whether the use of oral tablets for a longer period on healthy people or people with greater tongue coating accumulation shifts their oral microbiome needs to be investigated.

Zero-Inflated gaussian mixed models for analyzing longitudinal microbiome data

Motivation

The human microbiome is variable and dynamic in nature. Longitudinal studies could explain the mechanisms in maintaining the microbiome in health or causing dysbiosis in disease. However, it remains challenging to properly analyze the longitudinal microbiome data from either 16S rRNA or metagenome shotgun sequencing studies, output as proportions or counts. Most microbiome data are sparse, requiring statistical models to handle zero-inflation. Moreover, longitudinal design induces correlation among the samples and thus further complicates the analysis and interpretation of the microbiome data.

Results

In this article, we propose zero-inflated Gaussian mixed models (ZIGMMs) to analyze longitudinal microbiome data. ZIGMMs is a robust and flexible method which can be applicable for longitudinal microbiome proportion data or count data generated with either 16S rRNA or shotgun sequencing technologies. It can include various types of fixed effects and random effects and account for various within-subject correlation structures, and can effectively handle zero-inflation. We developed an efficient Expectation-Maximization (EM) algorithm to fit the ZIGMMs by taking advantage of the standard procedure for fitting linear mixed models. We demonstrate the computational efficiency of our EM algorithm by comparing with two other zero-inflated methods. We show that ZIGMMs outperform the previously used linear mixed models (LMMs), negative binomial mixed models (NBMMs) and zero-inflated Beta regression mixed model (ZIBR) in detecting associated effects in longitudinal microbiome data through extensive simulations. We also apply our method to two public longitudinal microbiome datasets and compare with LMMs and NBMMs in detecting dynamic effects of associated taxa.

Mobilization of Antibiotic Resistance: Are Current Approaches for Colocalizing Resistomes and Mobilomes Useful?

Antimicrobial resistance (AMR) poses a global human and animal health threat, and predicting AMR persistence and transmission remains an intractable challenge. Shotgun metagenomic sequencing can help overcome this by enabling characterization of AMR genes within all bacterial taxa, most of which are uncultivatable in laboratory settings. Shotgun sequencing, therefore, provides a more comprehensive glance at AMR "potential" within samples, i.e., the "resistome." However, the risk inherent within a given resistome is predicated on the genomic context of various AMR genes, including their presence within mobile genetic elements (MGEs). Therefore, resistome risk stratification can be advanced if AMR profiles are considered in light of the flanking mobilizable genomic milieu (e.g., plasmids, integrative conjugative elements (ICEs), phages, and other MGEs). Because such mediators of horizontal gene transfer (HGT) are involved in uptake by pathogens, investigators are increasingly interested in characterizing that resistome fraction in genomic proximity to HGT mediators, i.e., the "mobilome"; we term this "colocalization." We explored the utility of common colocalization approaches using alignment- and assembly-based techniques, on clinical (human) and agricultural (cattle) fecal metagenomes, obtained from antimicrobial use trials. Ordination revealed that tulathromycin-treated cattle experienced a shift in ICE and plasmid composition versus untreated animals, though the resistome was unaffected during the monitoring period. Contrarily, the human resistome and mobilome composition both shifted shortly after antimicrobial administration, though this rebounded to pre-treatment status. Bayesian networks revealed statistical AMR-MGE co-occurrence in 19 and 2% of edges from the cattle and human networks, respectively, suggesting a putatively greater mobility potential of AMR in cattle feces. Conversely, using Mobility Index (MI) and overlap analysis, abundance of de novo-assembled contigs supporting resistomes flanked by MGE increased shortly post-exposure within human metagenomes, though > 40 days after peak dose such contigs were rare (∼2%). MI was not substantially altered by antimicrobial exposure across all cattle metagenomes, ranging 0.5-4.0%. We highlight that current alignment- and assembly-based methods estimating resistome mobility yield contradictory and incomplete results, likely constrained by approach-specific data inputs, and bioinformatic limitations. We discuss recent laboratory and computational advancements that may enhance resistome risk analysis in clinical, regulatory, and commercial applications.

Compositional Lotka-Volterra describes microbial dynamics in the simplex

Dynamic changes in microbial communities play an important role in human health and disease. Specifically, deciphering how microbial species in a community interact with each other and their environment can elucidate mechanisms of disease, a problem typically investigated using tools from community ecology. Yet, such methods require measurements of absolute densities, whereas typical datasets only provide estimates of relative abundances. Here, we systematically investigate models of microbial dynamics in the simplex of relative abundances. We derive a new nonlinear dynamical system for microbial dynamics, termed “compositional” Lotka-Volterra (cLV), unifying approaches using generalized Lotka-Volterra (gLV) equations from community ecology and compositional data analysis. On three real datasets, we demonstrate that cLV recapitulates interactions between relative abundances implied by gLV. Moreover, we show that cLV is as accurate as gLV in forecasting microbial trajectories in terms of relative abundances. We further compare cLV to two other models of relative abundance dynamics motivated by common assumptions in the literature—a linear model in a log-ratio transformed space, and a linear model in the space of relative abundances—and provide evidence that cLV more accurately describes community trajectories over time. Finally, we investigate when information about direct effects can be recovered from relative data that naively provide information about only indirect effects. Our results suggest that strong effects may be recoverable from relative data, but more subtle effects are challenging to identify.

Dynamic changes in microbial communities play an important role in human health and disease. Specifically, deciphering how microbial species in a community interact with each other and their environment can elucidate mechanisms of disease, a problem typically investigated using tools from community ecology. Yet, such methods require measurements of absolute densities, whereas typical only provide estimates of relative abundances. We investigate methods for describing microbial dynamics in terms of relative abundances using approaches from machine learning and dynamical systems. Across three real datasets, we show that relative abundances are sufficient to describe compositional dynamics. Additionally, we show that models trained on relative abundances alone predict future compositions as well models trained on absolute abundances. Finally, we provide criteria for when direct effects, which typically can only be learned from absolute abundances, are recoverable for relative data. As a proof of concept, we recapitulate a previously proposed interaction network for C. difficile colonization.

Applications of weighted association networks applied to compositional data in biology

Next-generation sequencing technologies have generated, and continue to produce, an increasingly large corpus of biological data. The data generated are inherently compositional as they convey only relative information dependent upon the capacity of the instrument, experimental design and technical bias. There is considerable information to be gained through network analysis by studying the interactions between components within a system. Network theory methods using compositional data are powerful approaches for quantifying relationships between biological components and their relevance to phenotype, environmental conditions or other external variables. However, many of the statistical assumptions used for network analysis are not designed for compositional data and can bias downstream results. In this mini-review, we illustrate the utility of network theory in biological systems and investigate modern techniques while introducing researchers to frameworks for implementation. We overview (1) compositional data analysis, (2) data transformations and (3) network theory along with insight on a battery of network types including static-, temporal-, sample-specific- and differential-networks. The intention of this mini-review is not to provide a comprehensive overview of network methods, rather to introduce microbiology researchers to (semi)-unsupervised data-driven approaches for inferring latent structures that may give insight into biological phenomena or abstract mechanics of complex systems.

Molecular networks in Network Medicine: Development and applications

Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.

Gut Microbiome Composition Is Associated with a Pathologic Response After Preoperative Chemoradiation in Patients with Rectal Cancer

PURPOSE

There are ongoing investigations to find promising biomarkers for predicting a complete response (CR) after concurrent chemoradiation (CCRT) in rectal cancer. We aimed to find the predictive value in the gut microbiome in terms of response after preoperative CCRT.

METHODS AND MATERIALS

We collected a total of 45 fecal samples from patients with rectal cancer before CCRT. Tumor response after CCRT was assessed according to the American Joint Committee on Cancer tumor regression grading system. Analysis of linear discriminant analysis effect size and MetaCyc pathway abundance predictions were performed to compare composition and metabolic function of microbiome between patients with and without CR. We also established a Bayesian network model to identify microbial networks and species to be related with CCRT response.

RESULTS

Seven patients (15.6%) demonstrated pathologically CR, and 38 patients (84.4%) showed non-CR after preoperative CCRT. Between CR and non-CR patients, there was a significant difference in terms of β-diversity (P = .028), but no difference in α-diversity was found. Bacteroidales (Bacteroidaceae, Rikenellaceae, Bacteroides) were relatively more abundant in patients with non-CR than those with CR. Pathways related to anabolic function predominated in CR patients. According to Bayesian network analysis, Duodenibacillus massiliensis was linked with the improved CR rate.

CONCLUSIONS

From the fecal microbiome using samples obtained before preoperative CCRT, differences in microbial community composition and functions were observed between patients with and without CR in rectal cancer. However, the finding that a specific taxon may be linked with the improved therapeutic response should be verified in a prospective setting.

Resources and tools for the high-throughput, multi-omic study of intestinal microbiota

The human gut microbiome impacts several aspects of human health and disease, including digestion, drug metabolism and the propensity to develop various inflammatory, autoimmune and metabolic diseases. Many of the molecular processes that play a role in the activity and dynamics of the microbiota go beyond species and genic composition and thus, their understanding requires advanced bioinformatics support. This article aims to provide an up-to-date view of the resources and software tools that are being developed and used in human gut microbiome research, in particular data integration and systems-level analysis efforts. These efforts demonstrate the power of standardized and reproducible computational workflows for integrating and analysing varied omics data and gaining deeper insights into microbe community structure and function as well as host-microbe interactions.

Gut microbial metabolites associated with HIV infection

HIV infection has been associated with alterations in gut microbiota and related microbial metabolite production. However, the mechanisms of how these functional microbial metabolites may affect HIV immunopathogenesis and comorbidities, such as cardiovascular disease and other metabolic diseases, remain largely unknown. Here we review the current understanding of gut microbiota and related metabolites in the context of HIV infection. We focus on several bacteria-produced metabolites, including tryptophan catabolites, short-chain fatty acids and trimethylamine-N-oxide (TMAO), and discuss their implications in HIV infection and comorbidities. We also prospect future studies using integrative multiomics approaches to better understand host-microbiota-metabolites interactions in HIV infection, and facilitate integrative medicine utilizing the microbiota in HIV infection.

Dynamic interaction network inference from longitudinal microbiome data

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.

Host and microbiome multi-omics integration: applications and methodologies

The study of the microbial community-the microbiome-associated with a human host is a maturing research field. It is increasingly clear that the composition of the human's microbiome is associated with various diseases such as gastrointestinal diseases, liver diseases and metabolic diseases. Using high-throughput technologies such as next-generation sequencing and mass spectrometry-based metabolomics, we are able to comprehensively sequence the microbiome-the metagenome-and associate these data with the genomic, epigenomics, transcriptomic and metabolic profile of the host. Our review summarises the application of integrating host omics with microbiome as well as the analytical methods and related tools applied in these studies. In addition, potential future directions are discussed.

Addressing Global Ruminant Agricultural Challenges Through Understanding the Rumen Microbiome: Past, Present, and Future

The rumen is a complex ecosystem composed of anaerobic bacteria, protozoa, fungi, methanogenic archaea and phages. These microbes interact closely to breakdown plant material that cannot be digested by humans, whilst providing metabolic energy to the host and, in the case of archaea, producing methane. Consequently, ruminants produce meat and milk, which are rich in high-quality protein, vitamins and minerals, and therefore contribute to food security. As the world population is predicted to reach approximately 9.7 billion by 2050, an increase in ruminant production to satisfy global protein demand is necessary, despite limited land availability, and whilst ensuring environmental impact is minimized. Although challenging, these goals can be met, but depend on our understanding of the rumen microbiome. Attempts to manipulate the rumen microbiome to benefit global agricultural challenges have been ongoing for decades with limited success, mostly due to the lack of a detailed understanding of this microbiome and our limited ability to culture most of these microbes outside the rumen. The potential to manipulate the rumen microbiome and meet global livestock challenges through animal breeding and introduction of dietary interventions during early life have recently emerged as promising new technologies. Our inability to phenotype ruminants in a high-throughput manner has also hampered progress, although the recent increase in “omic” data may allow further development of mathematical models and rumen microbial gene biomarkers as proxies. Advances in computational tools, high-throughput sequencing technologies and cultivation-independent “omics” approaches continue to revolutionize our understanding of the rumen microbiome. This will ultimately provide the knowledge framework needed to solve current and future ruminant livestock challenges.

Microbiome Data Accurately Predicts the Postmortem Interval Using Random Forest Regression Models

Death investigations often include an effort to establish the postmortem interval (PMI) in cases in which the time of death is uncertain. The postmortem interval can lead to the identification of the deceased and the validation of witness statements and suspect alibis. Recent research has demonstrated that microbes provide an accurate clock that starts at death and relies on ecological change in the microbial communities that normally inhabit a body and its surrounding environment. Here, we explore how to build the most robust Random Forest regression models for prediction of PMI by testing models built on different sample types (gravesoil, skin of the torso, skin of the head), gene markers (16S ribosomal RNA (rRNA), 18S rRNA, internal transcribed spacer regions (ITS)), and taxonomic levels (sequence variants, species, genus, etc.). We also tested whether particular suites of indicator microbes were informative across different datasets. Generally, results indicate that the most accurate models for predicting PMI were built using gravesoil and skin data using the 16S rRNA genetic marker at the taxonomic level of phyla. Additionally, several phyla consistently contributed highly to model accuracy and may be candidate indicators of PMI.

Systems biology of the human microbiome

Recent research has shown that the microbiome-a collection of microorganisms, including bacteria, fungi, and viruses, living on and in a host-are of extraordinary importance in human health, even from conception and development in the uterus. Therefore, to further our ability to diagnose disease, to predict treatment outcomes, and to identify novel therapeutics, it is essential to include microbiome and microbial metabolic biomarkers in Systems Biology investigations. In clinical studies or, more precisely, Systems Medicine approaches, we can use the diversity and individual characteristics of the personal microbiome to enhance our resolution for patient stratification. In this review, we explore several Systems Medicine approaches, including Microbiome Wide Association Studies to understand the role of the human microbiome in health and disease, with a focus on 'preventive medicine' or P4 (i.e., personalized, predictive, preventive, participatory) medicine.

Experimental design and quantitative analysis of microbial community multiomics

Studies of the microbiome have become increasingly sophisticated, and multiple sequence-based, molecular methods as well as culture-based methods exist for population-scale microbiome profiles. To link the resulting host and microbial data types to human health, several experimental design considerations, data analysis challenges, and statistical epidemiological approaches must be addressed. Here, we survey current best practices for experimental design in microbiome molecular epidemiology, including technologies for generating, analyzing, and integrating microbiome multiomics data. We highlight studies that have identified molecular bioactives that influence human health, and we suggest steps for scaling translational microbiome research to high-throughput target discovery across large populations.

Computational Approaches for Integrative Analysis of the Metabolome and Microbiome

The study of the microbiome, the totality of all microbes inhabiting the host or an environmental niche, has experienced exponential growth over the past few years. The microbiome contributes functional genes and metabolites, and is an important factor for maintaining health. In this context, metabolomics is increasingly applied to complement sequencing-based approaches (marker genes or shotgun metagenomics) to enable resolution of microbiome-conferred functionalities associated with health. However, analyzing the resulting multi-omics data remains a significant challenge in current microbiome studies. In this review, we provide an overview of different computational approaches that have been used in recent years for integrative analysis of metabolome and microbiome data, ranging from statistical correlation analysis to metabolic network-based modeling approaches. Throughout the process, we strive to present a unified conceptual framework for multi-omics integration and interpretation, as well as point out potential future directions.

Modeling time-series data from microbial communities

As sequencing technologies have advanced, the amount of information regarding the composition of bacterial communities from various environments (for example, skin or soil) has grown exponentially. To date, most work has focused on cataloging taxa present in samples and determining whether the distribution of taxa shifts with exogenous covariates. However, important questions regarding how taxa interact with each other and their environment remain open thus preventing in-depth ecological understanding of microbiomes. Time-series data from 16S rDNA amplicon sequencing are becoming more common within microbial ecology, but methods to infer ecological interactions from these longitudinal data are limited. We address this gap by presenting a method of analysis using Poisson regression fit with an elastic-net penalty that (1) takes advantage of the fact that the data are time series; (2) constrains estimates to allow for the possibility of many more interactions than data; and (3) is scalable enough to handle data consisting of thousands of taxa. We test the method on gut microbiome data from white-throated woodrats (Neotoma albigula) that were fed varying amounts of the plant secondary compound oxalate over a period of 22 days to estimate interactions between OTUs and their environment.

Disentangling Interactions in the Microbiome: A Network Perspective

Microbiota are now widely recognized as being central players in the health of all organisms and ecosystems, and subsequently have been the subject of intense study. However, analyzing and converting microbiome data into meaningful biological insights remain very challenging. In this review, we highlight recent advances in network theory and their applicability to microbiome research. We discuss emerging graph theoretical concepts and approaches used in other research disciplines and demonstrate how they are well suited for enhancing our understanding of the higher-order interactions that occur within microbiomes. Network-based analytical approaches have the potential to help disentangle complex polymicrobial and microbe-host interactions, and thereby further the applicability of microbiome research to personalized medicine, public health, environmental and industrial applications, and agriculture.

Host-Microbiome Interaction and Cancer: Potential Application in Precision Medicine

It has been experimentally shown that host-microbial interaction plays a major role in shaping the wellness or disease of the human body. Microorganisms coexisting in human tissues provide a variety of benefits that contribute to proper functional activity in the host through the modulation of fundamental processes such as signal transduction, immunity and metabolism. The unbalance of this microbial profile, or dysbiosis, has been correlated with the genesis and evolution of complex diseases such as cancer. Although this latter disease has been thoroughly studied using different high-throughput (HT) technologies, its heterogeneous nature makes its understanding and proper treatment in patients a remaining challenge in clinical settings. Notably, given the outstanding role of host-microbiome interactions, the ecological interactions with microorganisms have become a new significant aspect in the systems that can contribute to the diagnosis and potential treatment of solid cancers. As a part of expanding precision medicine in the area of cancer research, efforts aimed at effective treatments for various kinds of cancer based on the knowledge of genetics, biology of the disease and host-microbiome interactions might improve the prediction of disease risk and implement potential microbiota-directed therapeutics. In this review, we present the state of the art of sequencing and metabolome technologies, computational methods and schemes in systems biology that have addressed recent breakthroughs of uncovering relationships or associations between microorganisms and cancer. Together, microbiome studies extend the horizon of new personalized treatments against cancer from the perspective of precision medicine through a synergistic strategy integrating clinical knowledge, HT data, bioinformatics, and systems biology.