Predicting microbial interactions with approaches based on flux balance analysis: an evaluation.

BACKGROUND: Given a genome-scale metabolic model (GEM) of a microorganism and criteria for optimization, flux balance analysis (FBA) predicts the optimal growth rate and its corresponding flux distribution for a specific medium. FBA has been extended to microbial consortia and thus can be used to predict interactions by comparing in-silico growth rates for co- and monocultures. Although FBA-based methods for microbial interaction prediction are becoming popular, a systematic evaluation of their accuracy has not yet been performed. RESULTS: Here, we evaluate the accuracy of FBA-based predictions of human and mouse gut bacterial interactions using growth data from the literature. For this, we collected 26 GEMs from the semi-curated AGORA database as well as four previously published curated GEMs. We tested the accuracy of three tools (COMETS, Microbiome Modeling Toolbox and MICOM) by comparing growth rates predicted in mono- and co-culture to growth rates extracted from the literature and also investigated the impact of different tool settings and media. We found that except for curated GEMs, predicted growth rates and their ratios (i.e. interaction strengths) do not correlate with growth rates and interaction strengths obtained from in vitro data. CONCLUSIONS: Prediction of growth rates with FBA using semi-curated GEMs is currently not sufficiently accurate to predict interaction strengths reliably.

Starvation responses impact interaction dynamics of human gut bacteria Bacteroides thetaiotaomicron and Roseburia intestinalis.

Bacterial growth often alters the environment, which in turn can impact interspecies interactions among bacteria. Here, we used an in vitro batch system containing mucin beads to emulate the dynamic host environment and to study its impact on the interactions between two abundant and prevalent human gut bacteria, the primary fermenter Bacteroides thetaiotaomicron and the butyrate producer Roseburia intestinalis. By combining machine learning and flow cytometry, we found that the number of viable B. thetaiotaomicron cells decreases with glucose consumption due to acid production, while R. intestinalis survives post-glucose depletion by entering a slow growth mode. Both species attach to mucin beads, but only viable cell counts of B. thetaiotaomicron increase significantly. The number of viable co-culture cells varies significantly over time compared to those of monocultures. A combination of targeted metabolomics and RNA-seq showed that the slow growth mode of R. intestinalis represents a diauxic shift towards acetate and lactate consumption, whereas B. thetaiotaomicron survives glucose depletion and low pH by foraging on mucin sugars. In addition, most of the mucin monosaccharides we tested inhibited the growth of R. intestinalis but not B. thetaiotaomicron. We encoded these causal relationships in a kinetic model, which reproduced the observed dynamics. In summary, we explored how R. intestinalis and B. thetaiotaomicron respond to nutrient scarcity and how this affects their dynamics. We highlight the importance of understanding bacterial metabolic strategies to effectively modulate microbial dynamics in changing conditions.

Deciphering Interactions Within a 4-Strain Riverine Bacterial Community.

The dynamics of a community of four planktonic bacterial strains isolated from river water was followed in R2 broth for 72 h in batch experiments. These strains were identified as Janthinobacterium sp., Brevundimonas sp., Flavobacterium sp. and Variovorax sp. 16S rRNA gene sequencing and flow cytometry analyses were combined to monitor the change in abundance of each individual strain in bi-cultures and quadri-culture. Two interaction networks were constructed that summarize the impact of the strains on each other's growth rate in exponential phase and carrying capacity in stationary phase. The networks agree on the absence of positive interactions but also show differences, implying that ecological interactions can be specific to particular growth phases. Janthinobacterium sp. was the fastest growing strain and dominated the co-cultures. However, its growth rate was negatively affected by the presence of other strains 10 to 100 times less abundant than Janthinobacterium sp. In general, we saw a positive correlation between growth rate and carrying capacity in this system. In addition, growth rate in monoculture was predictive of carrying capacity in co-culture. Taken together, our results highlight the necessity to take growth phases into account when measuring interactions within a microbial community. In addition, evidence that a minor strain can greatly influence the dynamics of a dominant one underlines the necessity to choose population models that do not assume a linear dependency of interaction strength to abundance of other species for accurate parameterization from such empirical data.

Thermocaproicibacter melissae gen. nov., sp. nov., a thermophilic chain-elongating bacterium, producing n-caproate from polymeric carbohydrates.

Strain MDTJ8T is a chain-elongating thermophilic bacterium isolated from a thermophilic acidogenic anaerobic digestor treating human waste while producing the high commodity chemical n-caproate. The strain grows and produces formate, acetate, n-butyrate, n-caproate and lactate from mono-, di- and polymeric saccharides at 37-60 °C (optimum, 50-55 °C) and at pH 5.0-7.0 (optimum, pH 6.5). The organism is an obligate anaerobe, is motile and its cells form rods (0.3-0.5×1.0-3.0 µm) that stain Gram-positive and occur primarily as chains. Phylogenetic analysis of both the 16S rRNA gene and full genome sequence shows that strain MDTJ8T belongs to a group that consists of mesophylic chain-elongating bacteria within the family Oscillospiraceae, being nearest to Caproicibacter fermentans EA1T (94.8 %) and Caproiciproducens galactitolivorans BS-1T (93.7 %). Its genome (1.96 Mbp) with a G+C content of 49.6 mol% is remarkably smaller than those of other chain-elongating bacteria of the family Oscillospiraceae. Pairwise average nucleotide identity and DNA-DNA hybridization values between strain MDJT8T and its mesophilic family members are less than 70 and 35 %, respectively, while pairwise average amino acid identity values are less than 68 %. In addition, strain MDJT8T uses far less carbohydrate and non-carbohydrate substrates compared to its nearest family members. The predominant cellular fatty acids of strain MDTJ8T are C14 : 0, C14 : 0 DMA (dimethyl acetal) and C16 : 0, while its polar lipid profile shows three unidentified glycophospholipids, 11 glycolipids, 13 phospholipids and six unidentified lipids. No respiratory quinones and polyamines are detected. Based on its phylogenetic, genotypic, morphological, physiological, biochemical and chemotaxonomic characteristics, strain MDTJ8T represents a novel species and novel genus of the family Oscillospiraceae and Thermocaproicibacter melissae gen. nov., sp. nov. is proposed as its name. The type strain is MDTJ8T (=DSM 114174T=LMG 32615T=NCCB 100883T).

Metabolic models of human gut microbiota: Advances and challenges.

The human gut is a complex ecosystem consisting of hundreds of microbial species interacting with each other and with the human host. Mathematical models of the gut microbiome integrate our knowledge of this system and help to formulate hypotheses to explain observations. The generalized Lotka-Volterra model has been widely used for this purpose, but it does not describe interaction mechanisms and thus does not account for metabolic flexibility. Recently, models that explicitly describe gut microbial metabolite production and consumption have become popular. These models have been used to investigate the factors that shape gut microbial composition and to link specific gut microorganisms to changes in metabolite concentrations found in diseases. Here, we review how such models are built and what we have learned so far from their application to human gut microbiome data. In addition, we discuss current challenges of these models and how these can be addressed in the future.

Technical versus biological variability in a synthetic human gut community.

Synthetic communities grown in well-controlled conditions are an important tool to decipher the mechanisms driving community dynamics. However, replicate time series of synthetic human gut communities in chemostats are rare, and it is thus still an open question to what extent stochasticity impacts gut community dynamics. Here, we address this question with a synthetic human gut bacterial community using an automated fermentation system that allows for a larger number of biological replicates. We collected six biological replicates for a community initially consisting of five common gut bacterial species that fill different metabolic niches. After an initial 12 hours in batch mode, we switched to chemostat mode and observed the community to stabilize after 2-3 days. Community profiling with 16S rRNA resulted in high variability across replicate vessels and high technical variability, while the variability across replicates was significantly lower for flow cytometric data. Both techniques agree on the decrease in the abundance of Bacteroides thetaiotaomicron, accompanied by an initial increase in Blautia hydrogenotrophica. These changes occurred together with reproducible metabolic shifts, namely a fast depletion of glucose and trehalose concentration in batch followed by a decrease in formic acid and pyruvic acid concentrations within the first 12 hours after the switch to chemostat mode. In conclusion, the observed variability in the synthetic bacterial human gut community, as assessed with 16S rRNA gene sequencing, is largely due to technical variability. The low variability seen in HPLC and flow cytometry data suggests a highly deterministic system.

Isolation and characterization of a thermophilic chain elongating bacterium that produces the high commodity chemical n-caproate from polymeric carbohydrates.

A thermophilic chain elongating bacterium, strain MDTJ8, was isolated from a thermophilic acidogenic anaerobic digestor producing n-caproate from human waste, growing optimally at 50-55 °C and pH 6.5. 16S rRNA gene analysis suggests that MDTJ8 represents a new species/genus within a group currently composed of mesophilic chain elongators of the Oscillospiraceae family. Genome analysis showed that strain MDTJ8 contains homologues of genes encoding for chain elongation and energy conservation but also indicated n-caproate production from carbohydrates including polymeric substances. This was confirmed by culturing experiments in which MDTJ8 converted, at pH 6.5 and 55 °C, mono-, di- and polymeric carbohydrates (starch and hemicellulose) to n-caproate reaching concentrations up to 283 mg/L and accounting for up to 10 % of the measured fermentation products. MDTJ8 is the first axenic organism that thermophilically performs chain elongation, opening doors to understand and intensify thermophilic bioprocesses targeting anaerobic digestion towards the production of the value-added chemical n-caproate.

Tick microbial associations at the crossroad of horizontal and vertical transmission pathways.

BACKGROUND: Microbial communities can affect disease risk by interfering with the transmission or maintenance of pathogens in blood-feeding arthropods. Here, we investigated whether bacterial communities vary between Ixodes ricinus nymphs which were or were not infected with horizontally transmitted human pathogens. METHODS: Ticks from eight forest sites were tested for the presence of Borrelia burgdorferi sensu lato, Babesia spp., Anaplasma phagocytophilum, and Neoehrlichia mikurensis by quantitative polymerase chain reaction (qPCR), and their microbiomes were determined by 16S rRNA amplicon sequencing. Tick bacterial communities clustered poorly by pathogen infection status but better by geography. As a second approach, we analysed variation in tick microorganism community structure (in terms of species co-infection) across space using hierarchical modelling of species communities. For that, we analysed almost 14,000 nymphs, which were tested for the presence of horizontally transmitted pathogens B. burgdorferi s.l., A. phagocytophilum, and N. mikurensis, and the vertically transmitted tick symbionts Rickettsia helvetica, Rickettsiella spp., Spiroplasma ixodetis, and Candidatus Midichloria mitochondrii. RESULTS: With the exception of Rickettsiella spp., all microorganisms had either significant negative (R. helvetica and A. phagocytophilum) or positive (S. ixodetis, N. mikurensis, and B. burgdorferi s.l.) associations with M. mitochondrii. Two tick symbionts, R. helvetica and S. ixodetis, were negatively associated with each other. As expected, both B. burgdorferi s.l. and N. mikurensis had a significant positive association with each other and a negative association with A. phagocytophilum. Although these few specific associations do not appear to have a large effect on the entire microbiome composition, they can still be relevant for tick-borne pathogen dynamics. CONCLUSIONS: Based on our results, we propose that M. mitochondrii alters the propensity of ticks to acquire or maintain horizontally acquired pathogens. The underlying mechanisms for some of these remarkable interactions are discussed herein and merit further investigation. Positive and negative associations between and within horizontally and vertically transmitted symbionts.

Quantitative microbial population study reveals geographical differences in bacterial symbionts of Ixodes ricinus.

BACKGROUND: Ixodes ricinus ticks vector pathogens that cause serious health concerns. Like in other arthropods, the microbiome may affect the tick's biology, with consequences for pathogen transmission. Here, we explored the bacterial communities of I. ricinus across its developmental stages and six geographic locations by the 16S rRNA amplicon sequencing, combined with quantification of the bacterial load. RESULTS: A wide range of bacterial loads was found. Accurate quantification of low microbial biomass samples permitted comparisons to high biomass samples, despite the presence of contaminating DNA. The bacterial communities of ticks were associated with geographical location rather than life stage, and differences in Rickettsia abundance determined this association. Subsequently, we explored the geographical distribution of four vertically transmitted symbionts identified in the microbiome analysis. For that, we screened 16,555 nymphs from 19 forest sites for R. helvetica, Rickettsiella spp., Midichloria mitochondrii, and Spiroplasma ixodetis. Also, the infection rates and distributions of these symbionts were compared to the horizontally transmitted pathogens Borrelia burgdorferi sensu lato, Anaplasma phagocytophilum, and Neoehrlichia mikurensis. The infection rates of all vertically transmitted symbionts differed between the study sites, and none of the symbionts was present in all tested ticks suggesting a facultative association with I. ricinus. The proportions in which symbionts occurred in populations of I. ricinus were highly variable, but geographically close study sites expressed similar proportions. These patterns were in contrast to what we observed for horizontally transmitted pathogens. Lastly, nearly 12% of tested nymphs were free of any targeted microorganisms, which is in line with the microbiome analyses. CONCLUSIONS: Our results show that the microbiome of I. ricinus is highly variable, but changes gradually and ticks originating from geographically close forest sites express similar bacterial communities. This suggests that geography-related factors affect the infection rates of vertically transmitted symbionts in I. ricinus. Since some symbionts, such as R. helvetica can cause disease in humans, we propose that public health investigations consider geographical differences in its infection rates.

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.

The virota and its transkingdom interactions in the healthy infant gut.

SignificanceMicrobes colonizing the infant gut during the first year(s) of life play an important role in immune system development. We show that after birth the (nearly) sterile gut is rapidly colonized by bacteria and their viruses (phages), which often show a strong cooccurrence. Most viruses infecting the infant do not cause clinical signs and their numbers strongly increase after day-care entrance. The infant diet is clearly reflected by identification of plant-infecting viruses, whereas fungi and parasites are not part of a stable gut microbiota. These temporal high-resolution baseline data about the gut colonization process will be valuable for further investigations of pathogenic viruses, dynamics between phages and their bacterial host, as well as studies investigating infants with a disturbed microbiota.

Fast quantification of gut bacterial species in cocultures using flow cytometry and supervised classification.

A bottleneck for microbial community experiments with many samples and/or replicates is the fast quantification of individual taxon abundances, which is commonly achieved through sequencing marker genes such as the 16S rRNA gene. Here, we propose a new approach for high-throughput and high-quality enumeration of human gut bacteria in a defined community, combining flow cytometry and supervised classification to identify and quantify species mixed in silico and in defined communities in vitro. We identified species in a 5-species in silico community with an F1 score of 71%. In addition, we demonstrate in vitro that our method performs equally well or better than 16S rRNA gene sequencing in two-species cocultures and agrees with 16S rRNA gene sequencing data on the most abundant species in a four-species community. We found that shape and size differences alone are insufficient to distinguish species, and that it is thus necessary to exploit the multivariate nature of flow cytometry data. Finally, we observed that variability of flow cytometry data across replicates differs between gut bacterial species. In conclusion, the performance of supervised classification of gut species in flow cytometry data is species-dependent, but is for some combinations accurate enough to serve as a faster alternative to 16S rRNA gene sequencing.

Regime transition Shapes the Composition, Assembly Processes, and Co-occurrence Pattern of Bacterioplankton Community in a Large Eutrophic Freshwater Lake.

At certain nutrient concentrations, shallow freshwater lakes are generally characterized by two contrasting ecological regimes with disparate patterns of biodiversity and biogeochemical cycles: a macrophyte-dominated regime (MDR) and a phytoplankton-dominated regime (PDR). To reveal ecological mechanisms that affect bacterioplankton along the regime shift, Illumina MiSeq sequencing of the 16S rRNA gene combined with a novel network clustering tool (Manta) were used to identify patterns of bacterioplankton community composition across the regime shift in Taihu Lake, China. Marked divergence in the composition and ecological assembly processes of bacterioplankton community was observed under the regime shift. The alpha diversity of the bacterioplankton community consistently and continuously decreased with the regime shift from MDR to PDR, while the beta diversity presents differently. Moreover, as the regime shifted from MDR to PDR, the contribution of deterministic processes (such as environmental selection) to the assembly of bacterioplankton community initially decreased and then increased again as regime shift from MDR to PDR, most likely as a consequence of differences in nutrient concentration. The topological properties, including modularity, transitivity and network diameter, of the bacterioplankton co-occurrence networks changed along the regime shift, and the co-occurrences among species changed in structure and were significantly shaped by the environmental variables along the regime transition from MDR to PDR. The divergent environmental state of the regimes with diverse nutritional status may be the most important factor that contributes to the dissimilarity of bacterioplankton community composition along the regime shift.

Fast and flexible analysis of linked microbiome data with mako.

Mako is a software tool that converts microbiome data and networks into a graph database and visualizes query results, thus allowing users without programming knowledge to carry out network-based queries. Mako is accompanied by a database compiled from 60 microbiome studies that is easily extended with the user's own data. We illustrate mako's strengths by enumerating association partners linked to propionate production and comparing frequencies of different network motifs across habitat types.

Successional Stages in Infant Gut Microbiota Maturation.

Disturbances in the primary colonization of the infant gut can result in lifelong consequences and have been associated with a range of host conditions. Although early-life factors have been shown to affect infant gut microbiota development, our current understanding of human gut colonization in early life remains limited. To gain more insights into the unique dynamics of this rapidly evolving ecosystem, we investigated the microbiota over the first year of life in eight densely sampled infants (n = 303 total samples). To evaluate the gut microbiota maturation transition toward an adult configuration, we compared the microbiome composition of the infants to that of the Flemish Gut Flora Project (FGFP) population (n = 1,106). We observed the infant gut microbiota to mature through three distinct, conserved stages of ecosystem development. Across these successional gut microbiota maturation stages, the genus predominance was observed to shift from Escherichia over Bifidobacterium to Bacteroides. Both disease and antibiotic treatment were observed to be associated occasionally with gut microbiota maturation stage regression, a transient setback in microbiota maturation dynamics. Although the studied microbiota trajectories evolved to more adult-like constellations, microbiome community typing against the background of the FGFP cohort clustered all infant samples within the (in adults) potentially dysbiotic Bacteroides 2 (Bact2) enterotype. We confirmed the similarities between infant gut microbial colonization and adult dysbiosis. Profound knowledge about the primary gut colonization process in infants might provide crucial insights into how the secondary colonization of a dysbiotic adult gut can be redirected. IMPORTANCE After birth, microbial colonization of the infant intestinal tract is important for health later in life. However, this initial process is highly dynamic and influenced by many factors. Studying this process in detail requires a dense longitudinal sampling effort. In the current study, the bacterial microbiota of >300 stool samples was analyzed from 8 healthy infants, suggesting that the infant gut microbial population matures along a path involving distinct microbial constellations and that the timing of these transitions is infant specific and can temporarily retrace upon external events. We also showed that the infant microbial populations show similarities to suboptimal bacterial populations in the guts of adults. These insights are crucial for a better understanding of the dynamics and characteristics of a "healthy gut microbial population" in both infants and adults and might allow the identification of intervention targets in cases of microbial disturbances or disease.

Gut Microbiome Profiling Uncovers a Lower Abundance of Butyricicoccus in Advanced Stages of Chronic Kidney Disease.

Chronic kidney disease (CKD) is characterized by the accumulation of uremic toxins which exert deleterious effects on various organ systems. Several of these uremic toxins originate from the bacterial metabolization of aromatic amino acids in the colon. This study assessed whether the gut microbial composition varies among patients in different stages of CKD. Uremic metabolites were quantified by UPLC/fluorescence detection and microbial profiling by 16S rRNA amplicon sequencing. Gut microbial profiles of CKD patients were compared among stages 1-2, stage 3 and stages 4-5. Although a substantial inter-individual difference in abundance of the top 15 genera was observed, no significant difference was observed between groups. Bristol stool scale (BSS) correlated negatively with p-cresyl sulfate and hippuric acid levels, irrespective of the intake of laxatives. Butyricicoccus, a genus with butyrate-generating properties, was decreased in abundance in advanced stages of CKD compared to the earlier stages (p = 0.043). In conclusion, in this cross-sectional study no gradual differences in the gut microbial profile over the different stages of CKD were observed. However, the decrease in the abundance of Butyricicoccus genus with loss of kidney function stresses the need for more in-depth functional exploration of the gut microbiome in CKD patients not on dialysis.

Temporal variability in quantitative human gut microbiome profiles and implications for clinical research.

While clinical gut microbiota research is ever-expanding, extending reference knowledge of healthy between- and within-subject gut microbiota variation and its drivers remains essential; in particular, temporal variability is under-explored, and a comparison with cross-sectional variation is missing. Here, we perform daily quantitative microbiome profiling on 713 fecal samples from 20 Belgian women over six weeks, combined with extensive anthropometric measurements, blood panels, dietary data, and stool characteristics. We show substantial temporal variation for most major gut genera; we find that for 78% of microbial genera, day-to-day absolute abundance variation is substantially larger within than between individuals, with up to 100-fold shifts over the study period. Diversity, and especially evenness indicators also fluctuate substantially. Relative abundance profiles show similar but less pronounced temporal variation. Stool moisture, and to a lesser extent diet, are the only significant host covariates of temporal microbiota variation, while menstrual cycle parameters did not show significant effects. We find that the dysbiotic Bact2 enterotype shows increased between- and within-subject compositional variability. Our results suggest that to increase diagnostic as well as target discovery power, studies could adopt a repeated measurement design and/or focus analysis on community-wide microbiome descriptors and indices.

Disentangling environmental effects in microbial association networks.

BACKGROUND: Ecological interactions among microorganisms are fundamental for ecosystem function, yet they are mostly unknown or poorly understood. High-throughput-omics can indicate microbial interactions through associations across time and space, which can be represented as association networks. Associations could result from either ecological interactions between microorganisms, or from environmental selection, where the association is environmentally driven. Therefore, before downstream analysis and interpretation, we need to distinguish the nature of the association, particularly if it is due to environmental selection or not. RESULTS: We present EnDED (environmentally driven edge detection), an implementation of four approaches as well as their combination to predict which links between microorganisms in an association network are environmentally driven. The four approaches are sign pattern, overlap, interaction information, and data processing inequality. We tested EnDED on networks from simulated data of 50 microorganisms. The networks contained on average 50 nodes and 1087 edges, of which 60 were true interactions but 1026 false associations (i.e., environmentally driven or due to chance). Applying each method individually, we detected a moderate to high number of environmentally driven edges-87% sign pattern and overlap, 67% interaction information, and 44% data processing inequality. Combining these methods in an intersection approach resulted in retaining more interactions, both true and false (32% of environmentally driven associations). After validation with the simulated datasets, we applied EnDED on a marine microbial network inferred from 10 years of monthly observations of microbial-plankton abundance. The intersection combination predicted that 8.3% of the associations were environmentally driven, while individual methods predicted 24.8% (data processing inequality), 25.7% (interaction information), and up to 84.6% (sign pattern as well as overlap). The fraction of environmentally driven edges among negative microbial associations in the real network increased rapidly with the number of environmental factors. CONCLUSIONS: To reach accurate hypotheses about ecological interactions, it is important to determine, quantify, and remove environmentally driven associations in marine microbial association networks. For that, EnDED offers up to four individual methods as well as their combination. However, especially for the intersection combination, we suggest using EnDED with other strategies to reduce the number of false associations and consequently the number of potential interaction hypotheses. Video abstract.