Identifying keystone species in the human gut microbiome from metagenomic timeseries using sparse linear regression

Biobased short chain fatty acid production - Exploring microbial community dynamics and metabolic networks through kinetic and microbial modeling approaches

In recent years, there has been growing interest in harnessing anaerobic digestion technology for resource recovery from waste streams. This approach has evolved beyond its traditional role in energy generation to encompass the production of valuable carboxylic acids, especially volatile fatty acids (VFAs) like acetic acid, propionic acid, and butyric acid. VFAs hold great potential for various industries and biobased applications due to their versatile properties. Despite increasing global demand, over 90% of VFAs are currently produced synthetically from petrochemicals. Realizing the potential of large-scale biobased VFA production from waste streams offers significant eco-friendly opportunities but comes with several key challenges. These include low VFA production yields, unstable acid compositions, complex and expensive purification methods, and post-processing needs. Among these, production yield and acid composition stand out as the most critical obstacles impacting economic viability and competitiveness. This paper seeks to offer a comprehensive view of combining complementary modeling approaches, including kinetic and microbial modeling, to understand the workings of microbial communities and metabolic pathways in VFA production, enhance production efficiency, and regulate acid profiles through the integration of omics and bioreactor data.

How aging influences the gut-bone marrow axis and alters hematopoietic stem cell regulation

The gut microbiome has come to prominence across research disciplines, due to its influence on major biological systems within humans. Recently, a relationship between the gut microbiome and hematopoietic system has been identified and coined the gut-bone marrow axis. It is well established that the hematopoietic system and gut microbiome separately alter with age; however, the relationship between these changes and how these systems influence each other demands investigation. Since the hematopoietic system produces immune cells that help govern commensal bacteria, it is important to identify how the microbiome interacts with hematopoietic stem cells (HSCs). The gut microbiota has been shown to influence the development and outcomes of hematologic disorders, suggesting dysbiosis may influence the maintenance of HSCs with age. Short chain fatty acids (SCFAs), lactate, iron availability, tryptophan metabolites, bacterial extracellular vesicles, microbe associated molecular patterns (MAMPs), and toll-like receptor (TLR) signalling have been proposed as key mediators of communication across the gut-bone marrow axis and will be reviewed in this article within the context of aging.

Kinetics-based inference of environment-dependent microbial interactions and their dynamic variation

Microbial communities in nature are dynamically evolving as member species change their interactions subject to environmental variations. Accounting for such context-dependent dynamic variations in interspecies interactions is critical for predictive ecological modeling. In the absence of generalizable theoretical foundations, we lack a fundamental understanding of how microbial interactions are driven by environmental factors, significantly limiting our capability to predict and engineer community dynamics and function. To address this issue, we propose a novel theoretical framework that allows us to represent interspecies interactions as an explicit function of environmental variables (such as substrate concentrations) by combining growth kinetics and a generalized Lotka-Volterra model. A synergistic integration of these two complementary models leads to the prediction of alterations in interspecies interactions as the outcome of dynamic balances between positive and negative influences of microbial species in mixed relationships. The effectiveness of our method was experimentally demonstrated using a synthetic consortium of two Escherichia coli mutants that are metabolically dependent (due to an inability to synthesize essential amino acids) but competitively grow on a shared substrate. The analysis of the E. coli binary consortium using our model not only showed how interactions between the two amino acid auxotrophic mutants are controlled by the dynamic shifts in limiting substrates but also enabled quantifying previously uncharacterizable complex aspects of microbial interactions, such as asymmetry in interactions. Our approach can be extended to other ecological systems to model their environment-dependent interspecies interactions from growth kinetics.IMPORTANCEModeling environment-controlled interspecies interactions through separate identification of positive and negative influences of microbes in mixed relationships is a new capability that can significantly improve our ability to understand, predict, and engineer the complex dynamics of microbial communities. Moreover, the prediction of microbial interactions as a function of environmental variables can serve as valuable benchmark data to validate modeling and network inference tools in microbial ecology, the development of which has often been impeded due to the lack of ground truth information on interactions. While demonstrated against microbial data, the theory developed in this work is readily applicable to general community ecology to predict interactions among macroorganisms, such as plants and animals, as well as microorganisms.

The MTIST platform: a microbiome time series inference standardized test

The human gut microbiome is a promising therapeutic target, but interventions are hampered by our limited understanding of microbial ecosystems. Here, we present a platform to develop, evaluate, and score approaches to learn ecological interactions from microbiome time series data. The microbiome time series inference standardized test (MTIST) comprises: a simulation framework for the in silico generation of microbiome study data akin to what is obtained with quantitative next-generation sequencing approaches, a compilation of a large curated data set generated by the simulation framework representing 648 simulated microbiome studies containing 18,360 time series, with a total of 2,182,800 species abundance measurements, and a scoring method to rank ecological inference algorithms. We use the MTIST platform to rank five implementations of microbiome inference approaches, revealing that while all algorithms performed well on ecosystems with few species (3 and 10), all algorithms failed to infer most interaction in a large ecosystem with 100 member species. However, we do find that the strongest interactions within a large ecosystem are inferred with higher success by all algorithms. Finally, we use the MTIST platform to compare different microbiome study designs, characterizing tradeoffs between samples per subject and number of subjects. Interestingly, we find that when only few samples can be collected per subject, ecological inference is most successful when these samples are collected with highest feasible temporal frequency. Taken together, we provide a computational tool to aid the development of better microbiome ecosystem inference approaches, which will be crucial towards the development of reliable and predictable therapeutic approaches that target the microbiome ecosystem.

Designing function-specific minimal microbiomes from large microbial communities

Microorganisms exist in large communities of diverse species, exhibiting various functionalities. The mammalian gut microbiome, for instance, has the functionality of digesting dietary fibre and producing different short-chain fatty acids. Not all microbes present in a community contribute to a given functionality; it is possible to find a minimal microbiome, which is a subset of the large microbiome, that is capable of performing the functionality while maintaining other community properties such as growth rate and metabolite production. Such a minimal microbiome will also contain keystone species for SCFA production in that community. In this work, we present a systematic constraint-based approach to identify a minimal microbiome from a large community for a user-proposed function. We employ a top-down approach with sequential deletion followed by solving a mixed-integer linear programming problem with the objective of minimising the L1-norm of the membership vector. Notably, we consider quantitative measures of community growth rate and metabolite production rates. We demonstrate the utility of our algorithm by identifying the minimal microbiomes corresponding to three model communities of the gut, and discuss their validity based on the presence of the keystone species in the community. Our approach is generic, flexible and finds application in studying a variety of microbial communities. The algorithm is available from https://github.com/RamanLab/minMicrobiome .

Gut microbial interactions based on network construction and bacterial pairwise cultivation

Association networks are widely applied for the prediction of bacterial interactions in studies of human gut microbiomes. However, the experimental validation of the predicted interactions is challenging due to the complexity of gut microbiomes and the limited number of cultivated bacteria. In this study, we addressed this challenge by integrating in vitro time series network (TSN) associations and co-cultivation of TSN taxon pairs. Fecal samples were collected and used for cultivation and enrichment of gut microbiome on YCFA agar plates for 13 days. Enriched cells were harvested for DNA extraction and metagenomic sequencing. A total of 198 metagenome-assembled genomes (MAGs) were recovered. Temporal dynamics of bacteria growing on the YCFA agar were used to infer microbial association networks. To experimentally validate the interactions of taxon pairs in networks, we selected 24 and 19 bacterial strains from this study and from the previously established human gut microbial biobank, respectively, for pairwise co-cultures. The co-culture experiments revealed that most of the interactions between taxa in networks were identified as neutralism (51.67%), followed by commensalism (21.67%), amensalism (18.33%), competition (5%) and exploitation (3.33%). Genome-centric analysis further revealed that the commensal gut bacteria (helpers and beneficiaries) might interact with each other via the exchanges of amino acids with high biosynthetic costs, short-chain fatty acids, and/or vitamins. We also validated 12 beneficiaries by adding 16 additives into the basic YCFA medium and found that the growth of 66.7% of these strains was significantly promoted. This approach provides new insights into the gut microbiome complexity and microbial interactions in association networks. Our work highlights that the positive relationships in gut microbial communities tend to be overestimated, and that amino acids, short-chain fatty acids, and vitamins are contributed to the positive relationships.

Resource competition predicts assembly of gut bacterial communities in vitro

Microbial community dynamics arise through interspecies interactions, including resource competition, cross-feeding and pH modulation. The individual contributions of these mechanisms to community structure are challenging to untangle. Here we develop a framework to estimate multispecies niche overlaps by combining metabolomics data of individual species, growth measurements in spent media and mathematical models. We applied our framework to an in vitro model system comprising 15 human gut commensals in complex media and showed that a simple model of resource competition accounted for most pairwise interactions. Next, we built a coarse-grained consumer-resource model by grouping metabolomic features depleted by the same set of species and showed that this model predicted the composition of 2-member to 15-member communities with reasonable accuracy. Furthermore, we found that incorporation of cross-feeding and pH-mediated interactions improved model predictions of species coexistence. Our theoretical model and experimental framework can be applied to characterize interspecies interactions in bacterial communities in vitro.

Metabolic model-based ecological modeling for probiotic design

The microbial community composition in the human gut has a profound effect on human health. This observation has lead to extensive use of microbiome therapies, including over-the-counter 'probiotic' treatments intended to alter the composition of the microbiome. Despite so much promise and commercial interest, the factors that contribute to the success or failure of microbiome-targeted treatments remain unclear. We investigate the biotic interactions that lead to successful engraftment of a novel bacterial strain introduced to the microbiome as in probiotic treatments. We use pairwise genome-scale metabolic modeling with a generalized resource allocation constraint to build a network of interactions between taxa that appear in an experimental engraftment study. We create induced sub-graphs using the taxa present in individual samples and assess the likelihood of invader engraftment based on network structure. To do so, we use a generalized Lotka-Volterra model, which we show has strong ability to predict if a particular invader or probiotic will successfully engraft into an individual's microbiome. Furthermore, we show that the mechanistic nature of the model is useful for revealing which microbe-microbe interactions potentially drive engraftment.

Bacteroides and related species: The keystone taxa of the human gut microbiota

Microbial communities play a significant role in maintaining ecosystems in a healthy homeostasis. Presently, in the human gastrointestinal tract, there are certain taxonomic groups of importance, though there is no single species that plays a keystone role. Bacteroides spp. are known to be major players in the maintenance of eubiosis in the human gastrointestinal tract. Here we review the critical role that Bacteroides play in the human gut, their potential pathogenic role outside of the gut, and their various methods of adapting to the environment, with a focus on data for B. fragilis and B. thetaiotaomicron. Bacteroides are anaerobic non-sporing Gram negative organisms that are also resistant to bile acids, generally thriving in the gut and having a beneficial relationship with the host. While they are generally commensal organisms, some Bacteroides spp. can be opportunistic pathogens in scenarios of GI disease, trauma, cancer, or GI surgery, and cause infection, most commonly intra-abdominal infection. B. fragilis can develop antimicrobial resistance through multiple mechanisms in large part due to its plasticity and fluid genome. Bacteroidota (formerly, Bacteroidetes) have a very broad metabolic potential in the GI microbiota and can rapidly adapt their carbohydrate metabolism to the available nutrients. Gastrointestinal Bacteroidota species produce short-chain fatty acids such as succinate, acetate, butyrate, and occasionally propionate, as the major end-products, which have wide-ranging and many beneficial influences on the host. Bacteroidota, via bile acid metabolism, also play a role in in colonization-resistance of other organisms, including Clostridioides difficile, and maintenance of gut integrity.

Study of amino acids absorption and gut microbiome on consumption of pea protein blended with enzymes-probiotics supplement

The current randomized, double-blind, crossover clinical trial was conducted to evaluate changes in the amino acid absorption and gut microbiota on consumption of pea protein supplemented with an enzymes-probiotics blend (Pepzyme Pro). A total of 15 healthy subjects were instructed to take test (pea protein + Pepzyme Pro) or placebo (pea protein + maltodextrin) for 15 days with a 30-day washout period. Blood samples were analyzed for plasma-free amino acids, insulin, and C-reactive protein (CRP). Additionally, nitrogen levels in urine and feces, along with the composition of gut microbiota, were evaluated. On day 15, the test arm showed a tendency to increase the rate of absorption and total absorption (AUC) of amino acids compared with the placebo arm, though the increase was statistically insignificant. In addition, 15-day test supplementation showed a tendency to reduce Tmax of all the amino acids (statistically insignificant except alanine, p = 0.021 and glycine, p = 0.023) in comparison with the placebo supplementation. There were no changes in urine and fecal nitrogen levels as well as serum CRP levels in the test and placebo arm. The increase in serum insulin level after 4 h was statistically significant in both arms, whereas the insulin level of the placebo and test arm at 4 h was not statistically different. Supplementation showed changes with respect to Archaea and few uncharacterized species but did not show statistically significant variations in microbiome profile at the higher taxonomic levels. A study with large sample size and detailed gut microbiome analysis is warranted to confirm the results statistically as well as to characterize altered species. However, the current study could provide an inkling of a positive alteration in protein digestibility, amino acid absorption, and gut microbiome with regular consumption of protein and enzymes-probiotics blend. Clinical Trial Registration:clinicaltrials.gov/; identifier [CTRI/2021/10/037072].

Identifying keystone species in microbial communities using deep learning

Previous studies suggested that microbial communities can harbour keystone species whose removal can cause a dramatic shift in microbiome structure and functioning. Yet, an efficient method to systematically identify keystone species in microbial communities is still lacking. Here we propose a data-driven keystone species identification (DKI) framework based on deep learning to resolve this challenge. Our key idea is to implicitly learn the assembly rules of microbial communities from a particular habitat by training a deep-learning model using microbiome samples collected from this habitat. The well-trained deep-learning model enables us to quantify the community-specific keystoneness of each species in any microbiome sample from this habitat by conducting a thought experiment on species removal. We systematically validated this DKI framework using synthetic data and applied DKI to analyse real data. We found that those taxa with high median keystoneness across different communities display strong community specificity. The presented DKI framework demonstrates the power of machine learning in tackling a fundamental problem in community ecology, paving the way for the data-driven management of complex microbial communities.

Stable States of a Microbial Community Are Formed by Dynamic Metabolic Networks with Members Functioning to Achieve Both Robustness and Plasticity

A more detailed understanding of the mechanisms underlying the formation of microbial communities is essential for the efficient management of microbial ecosystems. The stable states of microbial communities are commonly perceived as static and, thus, have not been extensively examined. The present study investigated stabilizing mechanisms, minority functions, and the reliability of quantitative ana-lyses, emphasizing a metabolic network perspective. A bacterial community, formed by batch transferred cultures supplied with phenol as the sole carbon and energy source and paddy soil as the inoculum, was analyzed using a principal coordinate ana-lysis (PCoA), mathematical models, and quantitative parameters defined as growth activity, community-changing activity, community-forming activity, vulnerable force, and resilience force depending on changes in the abundance of operational taxonomic units (OTUs) using 16S rRNA gene amplicon sequences. PCoA showed succession states until the 3rd transferred cultures and stable states from the 5th to 10th transferred cultures. Quantitative parameters indicated that the bacterial community was dynamic irrespective of the succession and stable states. Three activities fluctuated under stable states. Vulnerable and resilience forces were detected under the succession and stable states, respectively. Mathematical models indicated the construction of metabolic networks, suggesting the stabilizing mechanism of the community structure. Thirteen OTUs coexisted during stable states, and were recognized as core OTUs consisting of majorities, middle-class, and minorities. The abundance of the middle-class changed, whereas that of the others did not, which indicated that core OTUs maintained metabolic networks. Some extremely low abundance OTUs were consistently exchanged, suggesting a role for scavengers. These results indicate that stable states were formed by dynamic metabolic networks with members functioning to achieve robustness and plasticity.

Methodological Considerations in Longitudinal Analyses of Microbiome Data: A Comprehensive Review

Biological processes underlying health and disease are inherently dynamic and are best understood when characterized in a time-informed manner. In this comprehensive review, we discuss challenges inherent in time-series microbiome data analyses and compare available approaches and methods to overcome them. Appropriate handling of longitudinal microbiome data can shed light on important roles, functions, patterns, and potential interactions between large numbers of microbial taxa or genes in the context of health, disease, or interventions. We present a comprehensive review and comparison of existing microbiome time-series analysis methods, for both preprocessing and downstream analyses, including differential analysis, clustering, network inference, and trait classification. We posit that the careful selection and appropriate utilization of computational tools for longitudinal microbiome analyses can help advance our understanding of the dynamic host-microbiome relationships that underlie health-maintaining homeostases, progressions to disease-promoting dysbioses, as well as phases of physiologic development like those encountered in childhood.

Sparsity of higher-order landscape interactions enables learning and prediction for microbiomes

Microbiome engineering offers the potential to leverage microbial communities to improve outcomes in human health, agriculture, and climate. To translate this potential into reality, it is crucial to reliably predict community composition and function. But a brute force approach to cataloging community function is hindered by the combinatorial explosion in the number of ways we can combine microbial species. An alternative is to parameterize microbial community outcomes using simplified, mechanistic models, and then extrapolate these models beyond where we have sampled. But these approaches remain data-hungry, as well as requiring an a priori specification of what kinds of mechanisms are included and which are omitted. Here, we resolve both issues by introducing a mechanism-agnostic approach to predicting microbial community compositions and functions using limited data. The critical step is the identification of a sparse representation of the community landscape. We then leverage this sparsity to predict community compositions and functions, drawing from techniques in compressive sensing. We validate this approach on in silico community data, generated from a theoretical model. By sampling just [Formula: see text]1% of all possible communities, we accurately predict community compositions out of sample. We then demonstrate the real-world application of our approach by applying it to four experimental datasets and showing that we can recover interpretable, accurate predictions on composition and community function from highly limited data.

Gut commensals and their metabolites in health and disease

Purpose of review

This review comprehensively discusses the role of the gut microbiome and its metabolites in health and disease and sheds light on the importance of a holistic approach in assessing the gut.

Recent findings

The gut microbiome consisting of the bacteriome, mycobiome, archaeome, and virome has a profound effect on human health. Gut dysbiosis which is characterized by perturbations in the microbial population not only results in gastrointestinal (GI) symptoms or conditions but can also give rise to extra-GI manifestations. Gut microorganisms also produce metabolites (short-chain fatty acids, trimethylamine, hydrogen sulfide, methane, and so on) that are important for several interkingdom microbial interactions and functions. They also participate in various host metabolic processes. An alteration in the microbial species can affect their respective metabolite concentrations which can have serious health implications. Effective assessment of the gut microbiome and its metabolites is crucial as it can provide insights into one's overall health.

Summary

Emerging evidence highlights the role of the gut microbiome and its metabolites in health and disease. As it is implicated in GI as well as extra-GI symptoms, the gut microbiome plays a crucial role in the overall well-being of the host. Effective assessment of the gut microbiome may provide insights into one's health status leading to more holistic care.

Inferring microbial interactions with their environment from genomic and metagenomic data

Microbial communities assemble through a complex set of interactions between microbes and their environment, and the resulting metabolic impact on the host ecosystem can be profound. Microbial activity is known to impact human health, plant growth, water quality, and soil carbon storage which has lead to the development of many approaches and products meant to manipulate the microbiome. In order to understand, predict, and improve microbial community engineering, genome-scale modeling techniques have been developed to translate genomic data into inferred microbial dynamics. However, these techniques rely heavily on simulation to draw conclusions which may vary with unknown parameters or initial conditions, rather than more robust qualitative analysis. To better understand microbial community dynamics using genome-scale modeling, we provide a tool to investigate the network of interactions between microbes and environmental metabolites over time. Using our previously developed algorithm for simulating microbial communities from genome-scale metabolic models (GSMs), we infer the set of microbe-metabolite interactions within a microbial community in a particular environment. Because these interactions depend on the available environmental metabolites, we refer to the networks that we infer as metabolically contextualized, and so name our tool MetConSIN: Metabolically Contextualized Species Interaction Networks.

Phylosymbiosis: The Eco-Evolutionary Pattern of Insect-Symbiont Interactions

Insects harbor diverse assemblages of bacterial and fungal symbionts, which play crucial roles in host life history. Insects and their various symbionts represent a good model for studying host-microbe interactions. Phylosymbiosis is used to describe an eco-evolutionary pattern, providing a new cross-system trend in the research of host-associated microbiota. The phylosymbiosis pattern is characterized by a significant positive correlation between the host phylogeny and microbial community dissimilarities. Although host-symbiont interactions have been demonstrated in many insect groups, our knowledge of the prevalence and mechanisms of phylosymbiosis in insects is still limited. Here, we provide an order-by-order summary of the phylosymbiosis patterns in insects, including Blattodea, Coleoptera, Diptera, Hemiptera, Hymenoptera, and Lepidoptera. Then, we highlight the potential contributions of stochastic effects, evolutionary processes, and ecological filtering in shaping phylosymbiotic microbiota. Phylosymbiosis in insects can arise from a combination of stochastic and deterministic mechanisms, such as the dispersal limitations of microbes, codiversification between symbionts and hosts, and the filtering of phylogenetically conserved host traits (incl., host immune system, diet, and physiological characteristics).

Universal abundance fluctuations across microbial communities, tropical forests, and urban populations

The growth of complex populations, such as microbial communities, forests, and cities, occurs over vastly different spatial and temporal scales. Although research in different fields has developed detailed, system-specific models to understand each individual system, a unified analysis of different complex populations is lacking; such an analysis could deepen our understanding of each system and facilitate cross-pollination of tools and insights across fields. Here, we use a shared framework to analyze time-series data of the human gut microbiome, tropical forest, and urban employment. We demonstrate that a single, three-parameter model of stochastic population dynamics can reproduce the empirical distributions of population abundances and fluctuations in all three datasets. The three parameters characterizing a species measure its mean abundance, deterministic stability, and stochasticity. Our analysis reveals that, despite the vast differences in scale, all three systems occupy a similar region of parameter space when time is measured in generations. In other words, although the fluctuations observed in these systems may appear different, this difference is primarily due to the different physical timescales associated with each system. Further, we show that the distribution of temporal abundance fluctuations is described by just two parameters and derive a two-parameter functional form for abundance fluctuations to improve risk estimation and forecasting.

Analysis of strain, sex, and diet-dependent modulation of gut microbiota reveals candidate keystone organisms driving microbial diversity in response to American and ketogenic diets

BACKGROUND

The gut microbiota is modulated by a combination of diet, host genetics, and sex effects. The magnitude of these effects and interactions among them is important to understanding inter-individual variability in gut microbiota. In a previous study, mouse strain-specific responses to American and ketogenic diets were observed along with several QTLs for metabolic traits. In the current study, we searched for genetic variants underlying differences in the gut microbiota in response to American and ketogenic diets, which are high in fat and vary in carbohydrate composition, between C57BL/6 J (B6) and FVB/NJ (FVB) mouse strains.

RESULTS

Genetic mapping of microbial features revealed 18 loci under the QTL model (i.e., marginal effects that are not specific to diet or sex), 12 loci under the QTL by diet model, and 1 locus under the QTL by sex model. Multiple metabolic and microbial features map to the distal part of Chr 1 and Chr 16 along with eigenvectors extracted from principal coordinate analysis of measures of β-diversity. Bilophila, Ruminiclostridium 9, and Rikenella (Chr 1) were identified as sex- and diet-independent QTL candidate keystone organisms, and Parabacteroides (Chr 16) was identified as a diet-specific, candidate keystone organism in confirmatory factor analyses of traits mapping to these regions. For many microbial features, irrespective of which QTL model was used, diet or the interaction between diet and a genotype were the strongest predictors of the abundance of each microbial trait. Sex, while important to the analyses, was not as strong of a predictor for microbial abundances.

CONCLUSIONS

These results demonstrate that sex, diet, and genetic background have different magnitudes of effects on inter-individual differences in gut microbiota. Therefore, Precision Nutrition through the integration of genetic variation, microbiota, and sex affecting microbiota variation will be important to predict response to diets varying in carbohydrate composition. Video Abstract.

Microbial interactions in theory and practice: when are measurements compatible with models?

Most predictive models of ecosystem dynamics are based on interactions between organisms: their influence on each other's growth and death. We review here how theoretical approaches are used to extract interaction measurements from experimental data in microbiology, particularly focusing on the generalised Lotka-Volterra (gLV) framework. Though widely used, we argue that the gLV model should be avoided for estimating interactions in batch culture - the most common, simplest and cheapest in vitro approach to culturing microbes. Fortunately, alternative approaches offer a way out of this conundrum. Firstly, on the experimental side, alternatives such as the serial-transfer and chemostat systems more closely match the theoretical assumptions of the gLV model. Secondly, on the theoretical side, explicit organism-environment interaction models can be used to study the dynamics of batch-culture systems. We hope that our recommendations will increase the tractability of microbial model systems for experimentalists and theoreticians alike.

A computational method to dissect colonization resistance of the gut microbiota against pathogens

The mammalian gut microbiome protects the host through colonization resistance (CR) against the incursion of exogenous and often harmful microorganisms, but identifying the exact microbes responsible for the gut microbiota-mediated CR against a particular pathogen remains a challenge. To address this limitation, we developed a computational method: generalized microbe-phenotype triangulation (GMPT). We first systematically validated GMPT using a classical population dynamics model in community ecology and demonstrated its superiority over baseline methods. We then tested GMPT on simulated data generated from the ecological network inferred from a real community (GnotoComplex microflora) and real microbiome data on two mouse studies on Clostridioides difficile infection. We demonstrated GMPT's ability to streamline the discovery of microbes that are potentially responsible for microbiota-mediated CR against pathogens. GMPT holds promise to advance our understanding of CR mechanisms and facilitate the rational design of microbiome-based therapies for preventing and treating enteric infections.

Characteristics of rhizosphere and bulk soil microbial community of Chinese cabbage (Brassica campestris) grown in Karst area

Understanding the rhizosphere soil microbial community and its relationship with the bulk soil microbial community is critical for maintaining soil health and fertility and improving crop yields in Karst regions. The microbial communities in the rhizosphere and bulk soils of a Chinese cabbage (Brassica campestris) plantation in a Karst region, as well as their relationships with soil nutrients, were examined in this study using high-throughput sequencing technologies of 16S and ITS amplicons. The aim was to provide theoretical insights into the healthy cultivation of Chinese cabbage in a Karst area. The findings revealed that the rhizosphere soil showed higher contents of organic matter (OM), alkaline hydrolyzable nitrogen (AN), available phosphorus (AP), total phosphorus (TP), available potassium (AK), total potassium (TK), total nitrogen (TN), catalase (CA), urease (UR), sucrase (SU), and phosphatase (PHO), in comparison with bulk soil, while the pH value showed the opposite trend. The diversity of bacterial and fungal communities in the bulk soil was higher than that in the rhizosphere soil, and their compositions differed between the two types of soil. In the rhizosphere soil, Proteobacteria, Acidobacteriota, Actinobacteriota, and Bacteroidota were the dominant bacterial phyla, while Olpidiomycota, Ascomycota, Mortierellomycota, and Basidiomycota were the predominant fungal phyla. In contrast, the bulk soil was characterized by bacterial dominance of Proteobacteria, Acidobacteriota, Chloroflexi, and Actinobacteriota and fungal dominance of Ascomycota, Olpidiomycota, Mortierellomycota, and Basidiomycota. The fungal network was simpler than the bacterial network, and both networks exhibited less complexity in the rhizosphere soil compared with the bulk soil. Moreover, the rhizosphere soil harbored a higher proportion of beneficial Rhizobiales. The rhizosphere soil network was less complicated than the network in bulk soil by building a bacterial-fungal co-occurrence network. Furthermore, a network of relationships between soil properties and network keystone taxa revealed that the rhizosphere soil keystone taxa were more strongly correlated with soil properties than those in the bulk soil; despite its lower complexity, the rhizosphere soil contains a higher abundance of bacteria which are beneficial for cabbage growth compared with the bulk soil.

Growth phase estimation for abundant bacterial populations sampled longitudinally from human stool metagenomes

Longitudinal sampling of the stool has yielded important insights into the ecological dynamics of the human gut microbiome. However, human stool samples are available approximately once per day, while commensal population doubling times are likely on the order of minutes-to-hours. Despite this mismatch in timescales, much of the prior work on human gut microbiome time series modeling has assumed that day-to-day fluctuations in taxon abundances are related to population growth or death rates, which is likely not the case. Here, we propose an alternative model of the human gut as a stationary system, where population dynamics occur internally and the bacterial population sizes measured in a bolus of stool represent a steady-state endpoint of these dynamics. We formalize this idea as stochastic logistic growth. We show how this model provides a path toward estimating the growth phases of gut bacterial populations in situ. We validate our model predictions using an in vitro Escherichia coli growth experiment. Finally, we show how this method can be applied to densely-sampled human stool metagenomic time series data. We discuss how these growth phase estimates may be used to better inform metabolic modeling in flow-through ecosystems, like animal guts or industrial bioreactors.

Inferring microbial co-occurrence networks from amplicon data: a systematic evaluation

Microbes commonly organize into communities consisting of hundreds of species involved in complex interactions with each other. 16S ribosomal RNA (16S rRNA) amplicon profiling provides snapshots that reveal the phylogenies and abundance profiles of these microbial communities. These snapshots, when collected from multiple samples, can reveal the co-occurrence of microbes, providing a glimpse into the network of associations in these communities. However, the inference of networks from 16S data involves numerous steps, each requiring specific tools and parameter choices. Moreover, the extent to which these steps affect the final network is still unclear. In this study, we perform a meticulous analysis of each step of a pipeline that can convert 16S sequencing data into a network of microbial associations. Through this process, we map how different choices of algorithms and parameters affect the co-occurrence network and identify the steps that contribute substantially to the variance. We further determine the tools and parameters that generate robust co-occurrence networks and develop consensus network algorithms based on benchmarks with mock and synthetic data sets. The Microbial Co-occurrence Network Explorer, or MiCoNE (available at https://github.com/segrelab/MiCoNE) follows these default tools and parameters and can help explore the outcome of these combinations of choices on the inferred networks. We envisage that this pipeline could be used for integrating multiple data sets and generating comparative analyses and consensus networks that can guide our understanding of microbial community assembly in different biomes. IMPORTANCE Mapping the interrelationships between different species in a microbial community is important for understanding and controlling their structure and function. The surge in the high-throughput sequencing of microbial communities has led to the creation of thousands of data sets containing information about microbial abundances. These abundances can be transformed into co-occurrence networks, providing a glimpse into the associations within microbiomes. However, processing these data sets to obtain co-occurrence information relies on several complex steps, each of which involves numerous choices of tools and corresponding parameters. These multiple options pose questions about the robustness and uniqueness of the inferred networks. In this study, we address this workflow and provide a systematic analysis of how these choices of tools affect the final network and guidelines on appropriate tool selection for a particular data set. We also develop a consensus network algorithm that helps generate more robust co-occurrence networks based on benchmark synthetic data sets.

Stability of human gut microbiome: Comparison of ecological modelling and observational approaches

The gut microbiome plays a pivotal role in the human body, and perturbations in its composition have been linked to various disorders. Stability is an essential property of a healthy human gut microbiome, which allows it to maintain its functional richness under the external influences. This property has been explored through two distinct methodologies - mathematical modelling based on ecological principles and statistical analysis drawn from observations in interventional studies. Here we conducted a meta-analysis aimed to compare the two approaches utilising the data from 9 interventional and time series studies encompassing 3512 gut microbiome profiles obtained via 16S rRNA gene sequencing. By employing the previously published compositional Lotka-Volterra method, we modelled the dynamics of the microbial community and evaluated ecological stability measures. These measures were compared to those based on observed microbiome changes. There was a substantial correlation between the outcomes of the two approaches. Particularly, local stability assessed within the ecological paradigm was positively correlated with observational stability measures accounting for the compositional nature of microbiome data. Additionally, we were able to reproduce the previously reported inverse relationship between the community's robustness to microorganism loss and local stability, attributed to the distinct impacts of coefficient characterising the network decomposition on these two stability assessments. Our findings demonstrate harmonisation between the ecological and observational approaches to microbiome analysis, advancing the understanding of healthy gut microbiome concept. This paves the way to develop efficient microbiome-targeting interventions for disease prevention and treatment.

Top-down identification of keystone taxa in the microbiome

Keystone taxa in ecological communities are native taxa that play an especially important role in the stability of their ecosystem. However, we still lack an effective framework for identifying these taxa from the available high-throughput sequencing without the notoriously difficult step of reconstructing the detailed network of inter-specific interactions. In addition, while most microbial interaction models assume pair-wise relationships, it is yet unclear whether pair-wise interactions dominate the system, or whether higher-order interactions are relevant. Here we propose a top-down identification framework, which detects keystones by their total influence on the rest of the taxa. Our method does not assume a priori knowledge of pairwise interactions or any specific underlying dynamics and is appropriate to both perturbation experiments and metagenomic cross-sectional surveys. When applied to real high-throughput sequencing of the human gastrointestinal microbiome, we detect a set of candidate keystones and find that they are often part of a keystone module - multiple candidate keystone species with correlated occurrence. The keystone analysis of single-time-point cross-sectional data is also later verified by the evaluation of two-time-points longitudinal sampling. Our framework represents a necessary advancement towards the reliable identification of these key players of complex, real-world microbial communities.

Arsenic-Containing Medicine Treatment Disturbed the Human Intestinal Microbial Flora

Human intestinal microbiome plays vital role in maintaining intestinal homeostasis and interacting with xenobiotics. Few investigations have been conducted to understand the effect of arsenic-containing medicine exposure on gut microbiome. Most animal experiments are onerous in terms of time and resources and not in line with the international effort to reduce animal experiments. We explored the overall microbial flora by 16S rRNA genes analysis in fecal samples from acute promyelocytic leukemia (APL) patients treated with arsenic trioxide (ATO) plus all-trans retinoic acid (ATRA). Gut microbiomes were found to be overwhelmingly dominated by Firmicutes and Bacteroidetes after taking medicines containing arsenic in APL patients. The fecal microbiota composition of APL patients after treatment showed lower diversity and uniformity shown by the alpha diversity indices of Chao, Shannon, and Simpson. Gut microbiome operational taxonomic unit (OTU) numbers were associated with arsenic in the feces. We evaluated Bifidobacterium adolescentis and Lactobacillus mucosae to be a keystone in APL patients after treatment. Bacteroides at phylum or genus taxonomic levels were consistently affected after treatment. In the most common gut bacteria Bacteroides fragilis, arsenic resistance genes were significantly induced by arsenic exposure in anaerobic pure culture experiments. Without an animal model, without taking arsenicals passively, the results evidence that arsenic exposure by drug treatment is not only associated with alterations in intestinal microbiome development at the abundance and diversity level, but also induced arsenic biotransformation genes (ABGs) at the function levels which may even extend to arsenic-related health outcomes in APL.

Construction and Modeling of a Coculture Microplate for Real-Time Measurement of Microbial Interactions

The dynamic structures of microbial communities emerge from the complex network of interactions between their constituent microorganisms. Quantitative measurements of these interactions are important for understanding and engineering ecosystem structure. Here, we present the development and application of the BioMe plate, a redesigned microplate device in which pairs of wells are separated by porous membranes. BioMe facilitates the measurement of dynamic microbial interactions and integrates easily with standard laboratory equipment. We first applied BioMe to recapitulate recently characterized, natural symbiotic interactions between bacteria isolated from the Drosophila melanogaster gut microbiome. Specifically, the BioMe plate allowed us to observe the benefit provided by two Lactobacillus strains to an Acetobacter strain. We next explored the use of BioMe to gain quantitative insight into the engineered obligate syntrophic interaction between a pair of Escherichia coli amino acid auxotrophs. We integrated experimental observations with a mechanistic computational model to quantify key parameters associated with this syntrophic interaction, including metabolite secretion and diffusion rates. This model also allowed us to explain the slow growth observed for auxotrophs growing in adjacent wells by demonstrating that, under the relevant range of parameters, local exchange between auxotrophs is essential for efficient growth. The BioMe plate provides a scalable and flexible approach for the study of dynamic microbial interactions. IMPORTANCE Microbial communities participate in many essential processes from biogeochemical cycles to the maintenance of human health. The structure and functions of these communities are dynamic properties that depend on poorly understood interactions among different species. Unraveling these interactions is therefore a crucial step toward understanding natural microbiota and engineering artificial ones. Microbial interactions have been difficult to measure directly, largely due to limitations of existing methods to disentangle the contribution of different organisms in mixed cocultures. To overcome these limitations, we developed the BioMe plate, a custom microplate-based device that enables direct measurement of microbial interactions, by detecting the abundance of segregated populations of microbes that can exchange small molecules through a membrane. We demonstrated the possible application of the BioMe plate for studying both natural and artificial consortia. BioMe is a scalable and accessible platform that can be used to broadly characterize microbial interactions mediated by diffusible molecules.

New perspectives into the vaginal microbiome with systems biology

The vaginal microbiome (VMB) is critical to female reproductive health; however, the mechanisms associated with optimal and non-optimal states remain poorly understood due to the complex community structure and dynamic nature. Quantitative systems biology techniques applied to the VMB have improved understanding of community composition and function using primarily statistical methods. In contrast, fewer mechanistic models that use a priori knowledge of VMB features to develop predictive models have been implemented despite their use for microbiomes at other sites, including the gastrointestinal tract. Here, we explore systems biology approaches that have been applied in the VMB, highlighting successful techniques and discussing new directions that hold promise for improving understanding of health and disease.

Identifying keystone species in microbial communities using deep learning

Previous studies suggested that microbial communities harbor keystone species whose removal can cause a dramatic shift in microbiome structure and functioning. Yet, an efficient method to systematically identify keystone species in microbial communities is still lacking. This is mainly due to our limited knowledge of microbial dynamics and the experimental and ethical difficulties of manipulating microbial communities. Here, we propose a Data-driven Keystone species Identification (DKI) framework based on deep learning to resolve this challenge. Our key idea is to implicitly learn the assembly rules of microbial communities from a particular habitat by training a deep learning model using microbiome samples collected from this habitat. The well-trained deep learning model enables us to quantify the community-specific keystoneness of each species in any microbiome sample from this habitat by conducting a thought experiment on species removal. We systematically validated this DKI framework using synthetic data generated from a classical population dynamics model in community ecology. We then applied DKI to analyze human gut, oral microbiome, soil, and coral microbiome data. We found that those taxa with high median keystoneness across different communities display strong community specificity, and many of them have been reported as keystone taxa in literature. The presented DKI framework demonstrates the power of machine learning in tackling a fundamental problem in community ecology, paving the way for the data-driven management of complex microbial communities.

A symbiotic physical niche in Drosophila melanogaster regulates stable association of a multi-species gut microbiota

The gut is continuously invaded by diverse bacteria from the diet and the environment, yet microbiome composition is relatively stable over time for host species ranging from mammals to insects, suggesting host-specific factors may selectively maintain key species of bacteria. To investigate host specificity, we used gnotobiotic Drosophila, microbial pulse-chase protocols, and microscopy to investigate the stability of different strains of bacteria in the fly gut. We show that a host-constructed physical niche in the foregut selectively binds bacteria with strain-level specificity, stabilizing their colonization. Primary colonizers saturate the niche and exclude secondary colonizers of the same strain, but initial colonization by Lactobacillus species physically remodels the niche through production of a glycan-rich secretion to favor secondary colonization by unrelated commensals in the Acetobacter genus. Our results provide a mechanistic framework for understanding the establishment and stability of a multi-species intestinal microbiome.

Predicting Neurodegenerative Disease Using Prepathology Gut Microbiota Composition: a Longitudinal Study in Mice Modeling Alzheimer's Disease Pathologies

The gut microbiota-brain axis is suspected to contribute to the development of Alzheimer's disease (AD), a neurodegenerative disease characterized by amyloid-β plaque deposition, neurofibrillary tangles, and neuroinflammation. To evaluate the role of the gut microbiota-brain axis in AD, we characterized the gut microbiota of female 3xTg-AD mice modeling amyloidosis and tauopathy and wild-type (WT) genetic controls. Fecal samples were collected fortnightly from 4 to 52 weeks, and the V4 region of the 16S rRNA gene was amplified and sequenced on an Illumina MiSeq. RNA was extracted from the colon and hippocampus, converted to cDNA, and used to measure immune gene expression using reverse transcriptase quantitative PCR (RT-qPCR). Diversity metrics were calculated using QIIME2, and a random forest classifier was applied to predict bacterial features that are important in predicting mouse genotype. Gene expression of glial fibrillary acidic protein (GFAP; indicating astrocytosis) was elevated in the colon at 24 weeks. Markers of Th1 inflammation (il6) and microgliosis (mrc1) were elevated in the hippocampus. Gut microbiota were compositionally distinct early in life between 3xTg-AD mice and WT mice (permutational multivariate analysis of variance [PERMANOVA], 8 weeks, P = 0.001, 24 weeks, P = 0.039, and 52 weeks, P = 0.058). Mouse genotypes were correctly predicted 90 to 100% of the time using fecal microbiome composition. Finally, we show that the relative abundance of Bacteroides species increased over time in 3xTg-AD mice. Taken together, we demonstrate that changes in bacterial gut microbiota composition at prepathology time points are predictive of the development of AD pathologies. IMPORTANCE Recent studies have demonstrated alterations in the gut microbiota composition in mice modeling Alzheimer's disease (AD) pathologies; however, these studies have only included up to 4 time points. Our study is the first of its kind to characterize the gut microbiota of a transgenic AD mouse model, fortnightly, from 4 weeks of age to 52 weeks of age, to quantify the temporal dynamics in the microbial composition that correlate with the development of disease pathologies and host immune gene expression. In this study, we observed temporal changes in the relative abundances of specific microbial taxa, including the genus Bacteroides, that may play a central role in disease progression and the severity of pathologies. The ability to use features of the microbiota to discriminate between mice modeling AD and wild-type mice at prepathology time points indicates a potential role of the gut microbiota as a risk or protective factor in AD.

Controlling the human microbiome

We coexist with a vast number of microbes that live in and on our bodies. Those microbes and their genes are collectively known as the human microbiome, which plays important roles in human physiology and diseases. We have acquired extensive knowledge of the organismal compositions and metabolic functions of the human microbiome. However, the ultimate proof of our understanding of the human microbiome is reflected in our ability to manipulate it for health benefits. To facilitate the rational design of microbiome-based therapies, there are many fundamental questions to be addressed at the systems level. Indeed, we need a deep understanding of the ecological dynamics associated with such a complex ecosystem before we rationally design control strategies. In light of this, this review discusses progress from various fields, e.g., community ecology, network science, and control theory, that are helping us make progress toward the ultimate goal of controlling the human microbiome.

Analysis of strain, sex, and diet-dependent modulation of gut microbiota reveals candidate keystone organisms driving microbial diversity in response to American and ketogenic diets

Background The gut microbiota is modulated by a combination of diet, host genetics, and sex effects. The magnitude of these effects and interactions among them is important to understanding inter-individual variability in gut microbiota. In a previous study, mouse strain-specific responses to American and ketogenic diets were observed along with several QTL for metabolic traits. In the current study, we searched for genetic variants underlying differences in the gut microbiota in response to American and ketogenic diets, which are high in fat and vary in carbohydrate composition, between C57BL/6J (B6) and FVB/NJ (FVB) mouse strains. Results Genetic mapping of microbial features revealed 18 loci under the QTL model (i.e., marginal effects that are not specific to diet or sex), 12 loci under the QTL by diet model, and 1 locus under the QTL by sex model. Multiple metabolic and microbial features map to the distal part of Chr 1 and Chr 16 along with eigenvectors extracted from principal coordinate analysis of measures of β-diversity. Bilophila , Ruminiclostridium 9 , and Rikenella (Chr 1) were identified as sex and diet independent QTL candidate keystone organisms and Rikenelleceae RC9 Gut Group (Chr 16) was identified as a diet-specific, candidate keystone organism in confirmatory factor analyses of traits mapping to these regions. For many microbial features, irrespective of which QTL model was used, diet or the interaction between diet and a genotype were the strongest predictors of the abundance of each microbial trait. Sex, while important to the analyses, was not as strong of a predictor for microbial abundances. Conclusions These results demonstrate that sex, diet, and genetic background have different magnitudes of effects on inter-individual differences in gut microbiota. Therefore, Precision Nutrition through the integration of genetic variation, microbiota, and sex affecting microbiota variation will be important to predict response to diets varying in carbohydrate composition.

Ecological dynamics imposes fundamental challenges in community-based microbial source tracking

Quantifying the contributions of possible environmental sources ("sources") to a specific microbial community ("sink") is a classical problem in microbiology known as microbial source tracking (MST). Solving the MST problem will not only help us understand how microbial communities were formed, but also have far-reaching applications in pollution control, public health, and forensics. MST methods generally fall into two categories: target-based methods (focusing on the detection of source-specific indicator species or chemicals); and community-based methods (using community structure to measure similarity between sink samples and potential source environments). As next-generation sequencing becomes a standard community-assessment method in microbiology, numerous community-based computational methods, referred to as MST solvers hereafter have been developed and applied to various real datasets to demonstrate their utility across different contexts. Yet, those MST solvers do not consider microbial interactions and priority effects in microbial communities. Here, we revisit the performance of several representative MST solvers. We show compelling evidence that solving the MST problem using existing MST solvers is impractical when ecological dynamics plays a role in community assembly. In particular, we clearly demonstrate that the presence of either microbial interactions or priority effects will render the MST problem mathematically unsolvable for MST solvers. We further analyze data from fecal microbiota transplantation studies, finding that the state-of-the-art MST solvers fail to identify donors for most of the recipients. Finally, we perform community coalescence experiments to demonstrate that the state-of-the-art MST solvers fail to identify the sources for most of the sinks. Our findings suggest that ecological dynamics imposes fundamental challenges in MST. Interpretation of results of existing MST solvers should be done cautiously.

Beyond allergic progression: From molecules to microbes as barrier modulators in the gut-lung axis functionality

The "epithelial barrier hypothesis" states that a barrier dysfunction can result in allergy development due to tolerance breakdown. This barrier alteration may come from the direct contact of epithelial and immune cells with the allergens, and indirectly, through deleterious effects caused by environmental changes triggered by industrialization, pollution, and changes in the lifestyle. Apart from their protective role, epithelial cells can respond to external factors secreting IL-25 IL-33, and TSLP, provoking the activation of ILC2 cells and a Th2-biased response. Several environmental agents that influence epithelial barrier function, such as allergenic proteases, food additives or certain xenobiotics are reviewed in this paper. In addition, dietary factors that influence the allergenic response in a positive or negative way will be also described here. Finally, we discuss how the gut microbiota, its composition, and microbe-derived metabolites, such as short-chain fatty acids, alter not only the gut but also the integrity of distant epithelial barriers, focusing this review on the gut-lung axis.

An evaluation of homeostatic plasticity for ecosystems using an analytical data science approach

The natural world is constantly changing, and planetary boundaries are issuing severe warnings about biodiversity and cycles of carbon, nitrogen, and phosphorus. In other views, social problems such as global warming and food shortages are spreading to various fields. These seemingly unrelated issues are closely related, but it can be said that understanding them in an integrated manner is still a step away. However, progress in analytical technologies has been recognized in various fields and, from a microscopic perspective, with the development of instruments including next-generation sequencers (NGS), nuclear magnetic resonance (NMR), gas chromatography-mass spectrometry (GC/MS), and liquid chromatography-mass spectrometry (LC/MS), various forms of molecular information such as genome data, microflora structure, metabolome, proteome, and lipidome can be obtained. The development of new technology has made it possible to obtain molecular information in a variety of forms. From a macroscopic perspective, the development of environmental analytical instruments and environmental measurement facilities such as satellites, drones, observation ships, and semiconductor censors has increased the data availability for various environmental factors. Based on these background, the role of computational science is to provide a mechanism for integrating and understanding these seemingly disparate data sets. This review describes machine learning and the need for structural equations and statistical causal inference of these data to solve these problems. In addition to introducing actual examples of how these technologies can be utilized, we will discuss how to use these technologies to implement environmentally friendly technologies in society.

Arsenic shapes the microbial community structures in tungsten mine waste rocks

Tungsten (W) is a critical material that is widely used in military applications, electronics, lighting technology, power engineering and the automotive and aerospace industries. In recent decades, overexploitation of W has generated large amounts of mine waste rocks, which generate elevated content of toxic elements and cause serious adverse effects on ecosystems and public health. Microorganisms are considered important players in toxic element migrations from waste rocks. However, the understanding of how the microbial community structure varies in W mine waste rocks and its key driving factors is still unknown. In this study, high-throughput sequencing methods were used to determine the microbial community profiles along a W content gradient in W mine waste rocks. We found that the microbial community structures showed clear differences across the different W levels in waste rocks. Notably, arsenic (As), instead of W and nutrients, was identified as the most important predictor influencing microbial diversity. Furthermore, our results also showed that As is the most important environmental factor that regulates the distribution patterns of ecological clusters and keystone ASVs. Importantly, we found that the dominant genera have been regulated by As and were widely involved in As biogeochemical cycling in waste rocks. Taken together, our results have provided useful information about the response of microbial communities to W mine waste rocks.

Inference of dynamic interaction networks: A comparison between Lotka-Volterra and multivariate autoregressive models

Networks are ubiquitous throughout biology, spanning the entire range from molecules to food webs and global environmental systems. Yet, despite substantial efforts by the scientific community, the inference of these networks from data still presents a problem that is unsolved in general. One frequent strategy of addressing the structure of networks is the assumption that the interactions among molecular or organismal populations are static and correlative. While often successful, these static methods are no panacea. They usually ignore the asymmetry of relationships between two species and inferences become more challenging if the network nodes represent dynamically changing quantities. Overcoming these challenges, two very different network inference approaches have been proposed in the literature: Lotka-Volterra (LV) models and Multivariate Autoregressive (MAR) models. These models are computational frameworks with different mathematical structures which, nevertheless, have both been proposed for the same purpose of inferring the interactions within coexisting population networks from observed time-series data. Here, we assess these dynamic network inference methods for the first time in a side-by-side comparison, using both synthetically generated and ecological datasets. Multivariate Autoregressive and Lotka-Volterra models are mathematically equivalent at the steady state, but the results of our comparison suggest that Lotka-Volterra models are generally superior in capturing the dynamics of networks with non-linear dynamics, whereas Multivariate Autoregressive models are better suited for analyses of networks of populations with process noise and close-to linear behavior. To the best of our knowledge, this is the first study comparing LV and MAR approaches. Both frameworks are valuable tools that address slightly different aspects of dynamic networks.

Bacterial Keystone Taxa Regulate Carbon Metabolism in the Earthworm Gut

As important ecosystem engineers in soils, earthworms strongly influence carbon cycling through their burrowing and feeding activities. Earthworms do not perform these roles in isolation, because their intestines create a special habitat favorable for complex bacterial communities. However, how the ecological functioning of these earthworm-microbe interactions regulates carbon cycling remains largely unknown. To fill this knowledge gap, we investigated the bacterial community structure and carbon metabolic activities in the intestinal contents of earthworms and compared them to those of the adjacent soils in a long-term fertilization experiment. We discovered that earthworms harbored distinct bacterial communities compared to the surrounding soil under different fertilization conditions. The bacterial diversity was significantly larger in the adjacent soils than that in the earthworm gut. Three statistically identified keystone taxa in the bacterial networks, namely, Solirubrobacterales, Ktedonobacteraceae, and Jatrophihabitans, were shared across the earthworm gut and adjacent soil. Environmental factors (pH and organic matter) and keystone taxa were important determinants of the bacterial community composition in the earthworm gut. Both PICRUSt2 (Phylogenetic Investigation of Communities by Reconstruction of Unobserved States) and FAPROTAX (Functional Annotation of Prokaryotic Taxa) predicted that carbon metabolism was significantly higher in adjacent soil than in the earthworm gut, which was consistent with the average well color development obtained by the Biolog assay. Structural equation modeling combined with correlation analysis suggested that pH, organic matter, and potential keystone taxa exhibited significant relationships with carbon metabolism. This study deepens our understanding of the mechanisms underlying keystone taxa regulating carbon cycling in the earthworm gut. IMPORTANCE The intestinal microbiome of earthworms is a crucial component of the soil microbial community and nutrient cycling processes. If we could elucidate the role of this microbiome in regulating soil carbon metabolism, we would make a crucial contribution to understanding the ecological role of these gut bacterial taxa and to promoting sustainable agricultural development. However, the ecological functioning of these earthworm-microbe interactions in regulating carbon cycling has so far not been fully investigated. In this study, we revealed, first, that the bacterial groups of Solirubrobacterales, Ktedonobacteraceae, and Jatrophihabitans were core keystone taxa across the earthworm gut and adjacent soil and, second, that the environmental factors (pH and organic carbon) and keystone taxa strongly affected the bacterial community composition and exhibited close correlations with microbial carbon metabolism. Our results provide new insights into the community assembly of the earthworm gut microbiome and the ecological importance of potential keystone taxa in regulating carbon cycling dynamics.

Reconstruction of a Soil Microbial Network Induced by Stress Temperature

The microbial community is viewed as a network of diverse microorganisms connected by various interspecific interactions. While the stress gradient hypothesis (SGH) predicts that positive interactions are favored in more stressful environments, the prediction has been less explored in complex microbial communities due to the challenges of identifying interactions. Here, by applying a nonlinear time series analysis to the amplicon-based diversity time series data of the soil microbiota cultured under less stressful (30°C) or more stressful (37°C) temperature conditions, we show how the microbial network responds to temperature stress. While the genera that persisted only under the less stressful condition showed fewer positive effects, the genera that appeared only under the more stressful condition received more positive effects, in agreement with SGH. However, temperature difference also induced reconstruction of the community network, leading to an increased proportion of negative interactions at the whole-community level. The anti-SGH pattern can be explained by the stronger competition caused by increased metabolic rate and population densities. IMPORTANCE By combining amplicon-based diversity survey with recently developed nonlinear analytical tools, we successfully determined the interaction networks of more than 150 natural soil microbial genera under less or more temperature stress and explored the applicability of the stress gradient hypothesis to soil microbiota, shedding new light on the well-known hypothesis.

Decomposing predictability to identify dominant causal drivers in complex ecosystems

Ecosystems are complex systems of various physical, biological, and chemical processes. Since ecosystem dynamics are composed of a mixture of different levels of stochasticity and nonlinearity, handling these data is a challenge for existing methods of time series-based causal inferences. Here, we show that, by harnessing contemporary machine learning approaches, the concept of Granger causality can be effectively extended to the analysis of complex ecosystem time series and bridge the gap between dynamical and statistical approaches. The central idea is to use an ensemble of fast and highly predictive artificial neural networks to select a minimal set of variables that maximizes the prediction of a given variable. It enables decomposition of the relationship among variables through quantifying the contribution of an individual variable to the overall predictive performance. We show how our approach, EcohNet, can improve interaction network inference for a mesocosm experiment and simulated ecosystems. The application of the method to a long-term lake monitoring dataset yielded interpretable results on the drivers causing cyanobacteria blooms, which is a serious threat to ecological integrity and ecosystem services. Since performance of EcohNet is enhanced by its predictive capabilities, it also provides an optimized forecasting of overall components in ecosystems. EcohNet could be used to analyze complex and hybrid multivariate time series in many scientific areas not limited to ecosystems.

Effect of milk protein hydrolysate supplementation on protein energy malnutrition-induced gut dysbiosis

Dairy proteins in the diet are beneficial for the growth of probiotics; however, what is unknown is the gut-mediated immune responses under protein energy malnutrition (PEM) and if dairy protein hydrolysates can be effective as dietary interventions. This study compares the composition of the gut microbiota of rats with moderate protein deficiency (M.PEM) and severe protein deficiency (S.PEM) induced by feeding 5% and 1% hypoprotein diets, followed by replenishment with buffalo and whey protein hydrolysates. Fecal samples were collected, and the composition of the gut bacteria was analyzed by whole genome sequencing using long-read sequencing. Gene expression studies of the immunomodulatory cytokines involved and quantification of sIgA were carried out. IL-6 and IFN-γ were downregulated by about 0.17 ± 0.06 and 0.12 ± 0.10 fold when supplemented with whey protein hydrolysate in SP-RWC rats and by about 0.02 ± 0.06 and 0.35 ± 0.12 fold when using buffalo milk hydrolysate. The percentage of Firmicutes decreased in M.PEM and S.PEM rats (33.57%, 28.83 versus 47.73% of control at 3 weeks) but increased upon protein replenishment for all three protein sources at the end of nine weeks. The percentage of Bacteroidetes increased to 31.03% in S.PEM-induced rats as against 28.17% in control rats. The relative abundance of Lactobacillus sp. decreased in M.PEM and S.PEM rats while it showed the opposite effect upon protein replenishment. Gut microbiota modulated the pathogenesis of PEM differentially based on protein intervention along with a significant increase in the relative abundance of the keystone Lactobacillus genus.

Species abundance correlations carry limited information about microbial network interactions

Unraveling the network of interactions in ecological communities is a daunting task. Common methods to infer interspecific interactions from cross-sectional data are based on co-occurrence measures. For instance, interactions in the human microbiome are often inferred from correlations between the abundances of bacterial phylogenetic groups across subjects. We tested whether such correlation-based methods are indeed reliable for inferring interaction networks. For this purpose, we simulated bacterial communities by means of the generalized Lotka-Volterra model, with variation in model parameters representing variability among hosts. Our results show that correlations can be indicative for presence of bacterial interactions, but only when measurement noise is low relative to the variation in interaction strengths between hosts. Indication of interaction was affected by type of interaction network, process noise and sampling under non-equilibrium conditions. The sign of a correlation mostly coincided with the nature of the strongest pairwise interaction, but this is not necessarily the case. For instance, under rare conditions of identical interaction strength, we found that competitive and exploitative interactions can result in positive as well as negative correlations. Thus, cross-sectional abundance data carry limited information on specific interaction types. Correlations in abundance may hint at interactions but require independent validation.

The bacteria in and on our body (the human microbiome) largely determine how our body functions, and whether we stay healthy or get sick. These bacteria do not live on their own, but interact among each other and with their human host. Finding out which bacteria interact with each other is cumbersome, but patterns of joint occurrence between species might provide a clue to their ecological dependencies. We investigated whether correlations in species abundance can be used for the purpose of ecological network reconstruction. We simulated different bacterial communities with known interactions according to a theoretical population model. After having collected virtual samples from our simulated data, we performed a correlation analysis and then compared the correlation network with our known interaction network. We found that correlations can be informative for underlying interactions, but ecological conclusions should be drawn carefully. An obvious limitation of correlation analysis is that direction of interaction cannot be recovered from co-occurrence data, making correlations insensitive for detection of asymmetric interactions. In addition, we found that competitive and exploitative interactions can induce positive as well as negative correlations. We recommend careful interpretation and validation when inferring networks from cross-sectional abundance data.

Data-driven causal analysis of observational biological time series

Complex systems are challenging to understand, especially when they defy manipulative experiments for practical or ethical reasons. Several fields have developed parallel approaches to infer causal relations from observational time series. Yet, these methods are easy to misunderstand and often controversial. Here, we provide an accessible and critical review of three statistical causal discovery approaches (pairwise correlation, Granger causality, and state space reconstruction), using examples inspired by ecological processes. For each approach, we ask what it tests for, what causal statement it might imply, and when it could lead us astray. We devise new ways of visualizing key concepts, describe some novel pathologies of existing methods, and point out how so-called ‘model-free’ causality tests are not assumption-free. We hope that our synthesis will facilitate thoughtful application of methods, promote communication across different fields, and encourage explicit statements of assumptions. A video walkthrough is available (Video 1 or https://youtu.be/AlV0ttQrjK8).

Computational approach to modeling microbiome landscapes associated with chronic human disease progression

A microbial community is a dynamic system undergoing constant change in response to internal and external stimuli. These changes can have significant implications for human health. However, due to the difficulty in obtaining longitudinal samples, the study of the dynamic relationship between the microbiome and human health remains a challenge. Here, we introduce a novel computational strategy that uses massive cross-sectional sample data to model microbiome landscapes associated with chronic disease development. The strategy is based on the rationale that each static sample provides a snapshot of the disease process, and if the number of samples is sufficiently large, the footprints of individual samples populate progression trajectories, which enables us to recover disease progression paths along a microbiome landscape by using computational approaches. To demonstrate the validity of the proposed strategy, we developed a bioinformatics pipeline and applied it to a gut microbiome dataset available from a Crohn’s disease study. Our analysis resulted in one of the first working models of microbial progression for Crohn’s disease. We performed a series of interrogations to validate the constructed model. Our analysis suggested that the model recapitulated the longitudinal progression of microbial dysbiosis during the known clinical trajectory of Crohn’s disease. By overcoming restrictions associated with complex longitudinal sampling, the proposed strategy can provide valuable insights into the role of the microbiome in the pathogenesis of chronic disease and facilitate the shift of the field from descriptive research to mechanistic studies.

The delineation of system dynamics of a microbial community can provide a wealth of insights into the roles of the microbiome in the pathogenesis of chronic disease. However, due to the difficulty in obtaining longitudinal samples, most existing microbiome studies have been cross-sectional and largely descriptive. Here, we present a novel computational strategy that leverages massive static sample data to model microbiome landscapes associated with chronic disease development. To demonstrate the validity of the proposed strategy, we applied it to a gut microbiome dataset available from a Crohn’s disease study and constructed one of the first microbial progression models of the disease. We performed a series of interrogations on the constructed model. Our analysis suggested that the constructed model recapitulated the longitudinal progression of microbial dysbiosis during the known clinical trajectory of Crohn’s disease. By overcoming the sampling restrictions inherent to slowly progressive diseases, our approach is potentially widely applicable in many different studies across body sites, diseases, and other conditions.

Predicting the humification degree of multiple organic solid waste during composting using a designated bacterial community

Microbes are the drivers for disposing of organic solid waste (OSW) during aerobic fermentation. Notwithstanding, the significance of microbes is underestimated in numerous studies on aerobic fermentation product assessments. Here, we investigated the humification degree (HD), and the humic acid content was assessed in terms of the bacterial community. The bacterial communities were useful indicators for making predictions and even correctly determined the categories of OSWs with 94% accuracy. The bacterial codes can also provide a better prediction of HD. Our results demonstrate that the bacteria code is a reliable biological method to assess HD effectively. Bacterial codes can be used as ecological and biological indicators to evaluate the quality of aerobic fermentation of different materials.