Clinical Predictors of Recurrence After Primary Clostridioides difficile Infection: A Prospective Cohort Study.

BACKGROUND: Recurrent Clostridioides difficile infection (CDI) is a major public health threat. While clinical prediction tools exist, they do not incorporate the newest Infectious Diseases Society of America guidelines. METHODS: This was a prospective longitudinal study of patients experiencing their first episode of uncomplicated CDI. Patients were followed from diagnosis through 8 weeks post-completion of their anti-CDI therapy to assess recurrence. Stool was collected at diagnosis and weekly for 8 weeks following treatment. Recurrence was defined as diarrhea as well as a positive stool test by toxin EIA (EIA) for C. difficile. Fisher's exact test for binary variables and Student's t test for continuous variables were performed. Cox regression was performed to assess for predictors of CDI recurrence. RESULTS: Seventy-five patients were enrolled between August 1, 2015, and September 1, 2018. Mean age 58.1 years ± 15.5, 69.3% female, 74.7% were white, 11.3% had baseline irritable bowel syndrome, and 54.7% were actively using PPIs. Over the 8-week follow-up period, 22 patients developed a confirmed CDI recurrence. Univariate predictors of recurrence included treatment with metronidazole (40.9% vs 15.1%, p = 0.03), initially diagnosis by EIA (77.3% vs 43.4%, p = 0.007) and platelet count (206 ± 72.1 vs 270.9 ± 114.8, p = 0.03). A Cox regression model revealed primary diagnosis by EIA (HR 3.39, 95% CI 1.23, 9.31, p = 0.018) and treatment with metronidazole (HR 3.27 95% CI 1.31-8.19, p = 0.01) remain predictors for CDI recurrence. CONCLUSION: In a large prospective longitudinal cohort of uncomplicated CDI patients, treatment with metronidazole and diagnosis via EIA were the most robust predictors of CDI recurrence.

Improving microbial fitness in the mammalian gut by in vivo temporal functional metagenomics.

Elucidating functions of commensal microbial genes in the mammalian gut is challenging because many commensals are recalcitrant to laboratory cultivation and genetic manipulation. We present Temporal FUnctional Metagenomics sequencing (TFUMseq), a platform to functionally mine bacterial genomes for genes that contribute to fitness of commensal bacteria in vivo. Our approach uses metagenomic DNA to construct large-scale heterologous expression libraries that are tracked over time in vivo by deep sequencing and computational methods. To demonstrate our approach, we built a TFUMseq plasmid library using the gut commensal Bacteroides thetaiotaomicron (Bt) and introduced Escherichia coli carrying this library into germfree mice. Population dynamics of library clones revealed Bt genes conferring significant fitness advantages in E. coli over time, including carbohydrate utilization genes, with a Bt galactokinase central to early colonization, and subsequent dominance by a Bt glycoside hydrolase enabling sucrose metabolism coupled with co-evolution of the plasmid library and E. coli genome driving increased galactose utilization. Our findings highlight the utility of functional metagenomics for engineering commensal bacteria with improved properties, including expanded colonization capabilities in vivo.

A microbiota signature associated with experimental food allergy promotes allergic sensitization and anaphylaxis.

BACKGROUND: Commensal microbiota play a critical role in maintaining oral tolerance. The effect of food allergy on the gut microbial ecology remains unknown. OBJECTIVE: We sought to establish the composition of the gut microbiota in experimental food allergy and its role in disease pathogenesis. METHODS: Food allergy-prone mice with a gain-of-function mutation in the IL-4 receptor α chain (Il4raF709) and wild-type (WT) control animals were subjected to oral sensitization with chicken egg ovalbumin (OVA). Enforced tolerance was achieved by using allergen-specific regulatory T (Treg) cells. Community structure analysis of gut microbiota was performed by using a high-density 16S rDNA oligonucleotide microarrays (PhyloChip) and massively parallel pyrosequencing of 16S rDNA amplicons. RESULTS: OVA-sensitized Il4raF709 mice exhibited a specific microbiota signature characterized by coordinate changes in the abundance of taxa of several bacterial families, including the Lachnospiraceae, Lactobacillaceae, Rikenellaceae, and Porphyromonadaceae. This signature was not shared by similarly sensitized WT mice, which did not exhibit an OVA-induced allergic response. Treatment of OVA-sensitized Il4raF709 mice with OVA-specific Treg cells led to a distinct tolerance-associated signature coincident with the suppression of the allergic response. The microbiota of allergen-sensitized Il4raF709 mice differentially promoted OVA-specific IgE responses and anaphylaxis when reconstituted in WT germ-free mice. CONCLUSION: Mice with food allergy exhibit a specific gut microbiota signature capable of transmitting disease susceptibility and subject to reprogramming by enforced tolerance. Disease-associated microbiota may thus play a pathogenic role in food allergy.

A Spatial Multi-Modal Dissection of Host-Microbiome Interactions within the Colitis Tissue Microenvironment.

The intricate and dynamic interactions between the host immune system and its microbiome constituents undergo dynamic shifts in response to perturbations to the intestinal tissue environment. Our ability to study these events on the systems level is significantly limited by in situ approaches capable of generating simultaneous insights from both host and microbial communities. Here, we introduce Microbiome Cartography (MicroCart), a framework for simultaneous in situ probing of host features and its microbiome across multiple spatial modalities. We demonstrate MicroCart by comprehensively investigating the alterations in both gut host and microbiome components in a murine model of colitis by coupling MicroCart with spatial proteomics, transcriptomics, and glycomics platforms. Our findings reveal a global but systematic transformation in tissue immune responses, encompassing tissue-level remodeling in response to host immune and epithelial cell state perturbations, and bacterial population shifts, localized inflammatory responses, and metabolic process alterations during colitis. MicroCart enables a deep investigation of the intricate interplay between the host tissue and its microbiome with spatial multiomics.

Inferring dynamic signatures of microbes in complex host ecosystems.

The human gut microbiota comprise a complex and dynamic ecosystem that profoundly affects host development and physiology. Standard approaches for analyzing time-series data of the microbiota involve computation of measures of ecological community diversity at each time-point, or measures of dissimilarity between pairs of time-points. Although these approaches, which treat data as static snapshots of microbial communities, can identify shifts in overall community structure, they fail to capture the dynamic properties of individual members of the microbiota and their contributions to the underlying time-varying behavior of host ecosystems. To address the limitations of current methods, we present a computational framework that uses continuous-time dynamical models coupled with Bayesian dimensionality adaptation methods to identify time-dependent signatures of individual microbial taxa within a host as well as across multiple hosts. We apply our framework to a publicly available dataset of 16S rRNA gene sequences from stool samples collected over ten months from multiple human subjects, each of whom received repeated courses of oral antibiotics. Using new diversity measures enabled by our framework, we discover groups of both phylogenetically close and distant bacterial taxa that exhibit consensus responses to antibiotic exposure across multiple human subjects. These consensus responses reveal a timeline for equilibration of sub-communities of micro-organisms with distinct physiologies, yielding insights into the successive changes that occur in microbial populations in the human gut after antibiotic treatments. Additionally, our framework leverages microbial signatures shared among human subjects to automatically design optimal experiments to interrogate dynamic properties of the microbiota in new studies. Overall, our approach provides a powerful, general-purpose framework for understanding the dynamic behaviors of complex microbial ecosystems, which we believe will prove instrumental for future studies in this field.

MDITRE: Scalable and Interpretable Machine Learning for Predicting Host Status from Temporal Microbiome Dynamics.

Longitudinal microbiome data sets are being generated with increasing regularity, and there is broad recognition that these studies are critical for unlocking the mechanisms through which the microbiome impacts human health and disease. However, there is a dearth of computational tools for analyzing microbiome time-series data. To address this gap, we developed an open-source software package, Microbiome Differentiable Interpretable Temporal Rule Engine (MDITRE), which implements a new highly efficient method leveraging deep-learning technologies to derive human-interpretable rules that predict host status from longitudinal microbiome data. Using semi-synthetic and a large compendium of publicly available 16S rRNA amplicon and metagenomics sequencing data sets, we demonstrate that in almost all cases, MDITRE performs on par with or better than popular uninterpretable machine learning methods, and orders-of-magnitude faster than the prior interpretable technique. MDITRE also provides a graphical user interface, which we show through case studies can be used to derive biologically meaningful interpretations linking patterns of microbiome changes over time with host phenotypes. IMPORTANCE The human microbiome, or collection of microbes living on and within us, changes over time. Linking these changes to the status of the human host is crucial to understanding how the microbiome influences a variety of human diseases. Due to the large scale and complexity of microbiome data, computational methods are essential. Existing computational methods for linking changes in the microbiome to the status of the human host are either unable to scale to large and complex microbiome data sets or cannot produce human-interpretable outputs. We present a new computational method and software package that overcomes the limitations of previous methods, allowing researchers to analyze larger and more complex data sets while producing easily interpretable outputs. Our method has the potential to enable new insights into how changes in the microbiome over time maintain health or lead to disease in humans and facilitate the development of diagnostic tests based on the microbiome.

Gut metabolites predict Clostridioides difficile recurrence.

BACKGROUND: Clostridioides difficile infection (CDI) is the most common hospital acquired infection in the USA, with recurrence rates > 15%. Although primary CDI has been extensively linked to gut microbial dysbiosis, less is known about the factors that promote or mitigate recurrence. Moreover, previous studies have not shown that microbial abundances in the gut measured by 16S rRNA amplicon sequencing alone can accurately predict CDI recurrence. RESULTS: We conducted a prospective, longitudinal study of 53 non-immunocompromised participants with primary CDI. Stool sample collection began pre-CDI antibiotic treatment at the time of diagnosis, and continued up to 8 weeks post-antibiotic treatment, with weekly or twice weekly collections. Samples were analyzed using (1) 16S rRNA amplicon sequencing, (2) liquid chromatography/mass-spectrometry metabolomics measuring 1387 annotated metabolites, and (3) short-chain fatty acid profiling. The amplicon sequencing data showed significantly delayed recovery of microbial diversity in recurrent participants, and depletion of key anaerobic taxa at multiple time-points, including Clostridium cluster XIVa and IV taxa. The metabolomic data also showed delayed recovery in recurrent participants, and moreover mapped to pathways suggesting distinct functional abnormalities in the microbiome or host, such as decreased microbial deconjugation activity, lowered levels of endocannabinoids, and elevated markers of host cell damage. Further, using predictive statistical/machine learning models, we demonstrated that the metabolomic data, but not the other data sources, can accurately predict future recurrence at 1 week (AUC 0.77 [0.71, 0.86; 95% interval]) and 2 weeks (AUC 0.77 [0.69, 0.85; 95% interval]) post-treatment for primary CDI. CONCLUSIONS: The prospective, longitudinal, and multi-omic nature of our CDI recurrence study allowed us to uncover previously unrecognized dynamics in the microbiome and host presaging recurrence, and, in particular, to elucidate changes in the understudied gut metabolome. Moreover, we demonstrated that a small set of metabolites can accurately predict future recurrence. Our findings have implications for development of diagnostic tests and treatments that could ultimately short-circuit the cycle of CDI recurrence, by providing candidate metabolic biomarkers for diagnostics development, as well as offering insights into the complex microbial and metabolic alterations that are protective or permissive for recurrence. Video Abstract.

Discovering homotypic binding events at high spatial resolution.

MOTIVATION: Clusters of protein-DNA interaction events involving the same transcription factor are known to act as key components of invertebrate and mammalian promoters and enhancers. However, detecting closely spaced homotypic events from ChIP-Seq data is challenging because random variation in the ChIP fragmentation process obscures event locations. RESULTS: The Genome Positioning System (GPS) can predict protein-DNA interaction events at high spatial resolution from ChIP-Seq data, while retaining the ability to resolve closely spaced events that appear as a single cluster of reads. GPS models observed reads using a complexity penalized mixture model and efficiently predicts event locations with a segmented EM algorithm. An optional mode permits GPS to align common events across distinct experiments. GPS detects more joint events in synthetic and actual ChIP-Seq data and has superior spatial resolution when compared with other methods. In addition, the specificity and sensitivity of GPS are superior to or comparable with other methods. AVAILABILITY: http://cgs.csail.mit.edu/gps.

In vivo commensal control of Clostridioides difficile virulence.

Leveraging systems biology approaches, we illustrate how metabolically distinct species of Clostridia protect against or worsen Clostridioides difficile infection in mice by modulating the pathogen's colonization, growth, and virulence to impact host survival. Gnotobiotic mice colonized with the amino acid fermenter Paraclostridium bifermentans survive infection with reduced disease severity, while mice colonized with the butyrate-producer, Clostridium sardiniense, succumb more rapidly. Systematic in vivo analyses revealed how each commensal alters the gut-nutrient environment to modulate the pathogen's metabolism, gene regulatory networks, and toxin production. Oral administration of P. bifermentans rescues conventional, clindamycin-treated mice from lethal C. difficile infection in a manner similar to that of monocolonized animals, thereby supporting the therapeutic potential of this commensal species. Our findings lay the foundation for mechanistically informed therapies to counter C. difficile disease using systems biology approaches to define host-commensal-pathogen interactions in vivo.

High-resolution temporal profiling of the human gut microbiome reveals consistent and cascading alterations in response to dietary glycans.

BACKGROUND: Dietary glycans, widely used as food ingredients and not directly digested by humans, are of intense interest for their beneficial roles in human health through shaping the microbiome. Characterizing the consistency and temporal responses of the gut microbiome to glycans is critical for rationally developing and deploying these compounds as therapeutics. METHODS: We investigated the effect of two chemically distinct glycans (fructooligosaccharides and polydextrose) through three clinical studies conducted with 80 healthy volunteers. Stool samples, collected at dense temporal resolution (~ 4 times per week over 10 weeks) and analyzed using shotgun metagenomic sequencing, enabled detailed characterization of participants' microbiomes. For analyzing the microbiome time-series data, we developed MC-TIMME2 (Microbial Counts Trajectories Infinite Mixture Model Engine 2.0), a purpose-built computational tool based on nonparametric Bayesian methods that infer temporal patterns induced by perturbations and groups of microbes sharing these patterns. RESULTS: Overall microbiome structure as well as individual taxa showed rapid, consistent, and durable alterations across participants, regardless of compound dose or the order in which glycans were consumed. Significant changes also occurred in the abundances of microbial carbohydrate utilization genes in response to polydextrose, but not in response to fructooligosaccharides. Using MC-TIMME2, we produced detailed, high-resolution temporal maps of the microbiota in response to glycans within and across microbiomes. CONCLUSIONS: Our findings indicate that dietary glycans cause reproducible, dynamic, and differential alterations to the community structure of the human microbiome.

High-resolution computational models of genome binding events.

Direct physical information that describes where transcription factors, nucleosomes, modified histones, RNA polymerase II and other key proteins interact with the genome provides an invaluable mechanistic foundation for understanding complex programs of gene regulation. We present a method, joint binding deconvolution (JBD), which uses additional easily obtainable experimental data about chromatin immunoprecipitation (ChIP) to improve the spatial resolution of the transcription factor binding locations inferred from ChIP followed by DNA microarray hybridization (ChIP-Chip) data. Based on this probabilistic model of binding data, we further pursue improved spatial resolution by using sequence information. We produce positional priors that link ChIP-Chip data to sequence data by guiding motif discovery to inferred protein-DNA binding sites. We present results on the yeast transcription factors Gcn4 and Mig2 to demonstrate JBD's spatial resolution capabilities and show that positional priors allow computational discovery of the Mig2 motif when a standard approach fails.

MITRE: inferring features from microbiota time-series data linked to host status.

Longitudinal studies are crucial for discovering causal relationships between the microbiome and human disease. We present MITRE, the Microbiome Interpretable Temporal Rule Engine, a supervised machine learning method for microbiome time-series analysis that infers human-interpretable rules linking changes in abundance of clades of microbes over time windows to binary descriptions of host status, such as the presence/absence of disease. We validate MITRE's performance on semi-synthetic data and five real datasets. MITRE performs on par or outperforms conventional difficult-to-interpret machine learning approaches, providing a powerful new tool enabling the discovery of biologically interpretable relationships between microbiome and human host ( https://github.com/gerberlab/mitre/ ).

Dynamic Modulation of the Gut Microbiota and Metabolome by Bacteriophages in a Mouse Model.

The human gut microbiome is comprised of densely colonizing microorganisms including bacteriophages, which are in dynamic interaction with each other and the mammalian host. To address how bacteriophages impact bacterial communities in the gut, we investigated the dynamic effects of phages on a model microbiome. Gnotobiotic mice were colonized with defined human gut commensal bacteria and subjected to predation by cognate lytic phages. We found that phage predation not only directly impacts susceptible bacteria but also leads to cascading effects on other bacterial species via interbacterial interactions. Metabolomic profiling revealed that shifts in the microbiome caused by phage predation have a direct consequence on the gut metabolome. Our work provides insight into the ecological importance of phages as modulators of bacterial colonization, and it additionally suggests the potential impact of gut phages on the mammalian host with implications for their therapeutic use to precisely modulate the microbiome.

Impaired glucose metabolism and altered gut microbiome despite calorie restriction of ob/ob mice.

BACKGROUND: Growing evidence supports the role of gut microbiota in obesity and its related disorders including type 2 diabetes. Ob/ob mice, which are hyperphagic due to leptin deficiency, are commonly used models of obesity and were instrumental in suggesting links between gut microbiota and obesity. Specific changes in their gut microbiota such as decreased microbial diversity and increased Firmicutes to Bacteroidetes ratio have been suggested to contribute to obesity via increased microbiota capacity to harvest energy. However, the differential development of ob/ob mouse gut microbiota compared to wild type microbiota and the role of hyperphagia in their metabolic impairment have not been investigated thoroughly. RESULTS: We performed a 10-week long study in ob/ob (n = 12) and wild type control (n = 12) mice fed ad libitum. To differentiate effects of leptin deficiency from hyperphagia, we pair-fed an additional group of ob/ob mice (n = 11) based on the food consumption of control mice. Compared to control mice, ob/ob mice fed ad libitum exhibited compromised glucose metabolism and increased body fat percentage. Pair-fed ob/ob mice exhibited even more compromised glucose metabolism and maintained strikingly similar high body fat percentage at the cost of lean body mass. Acclimatization of the microbiota to our facility took up to 5 weeks. Leptin deficiency impacted gut microbial composition, explaining 18.3% of the variance. Pair-feeding also altered several taxa, although the overall community composition at the end of the study was not significantly different. We found 24 microbial taxa associations with leptin deficiency, notably enrichment of members of Lactobacillus and depletion of Akkermansia muciniphila. Microbial metabolic functions related to energy harvest, including glycan degradation, phosphotransferase systems and ABC transporters, were enriched in the ob/ob mice. Taxa previously reported as relevant for obesity were associated with body weight, including Oscillibacter and Alistipes (both negatively correlated) and Prevotella (positively correlated). CONCLUSIONS: Leptin deficiency caused major changes in the mouse gut microbiota composition. Several microbial taxa were associated with body composition. Pair-fed mice maintained a pre-set high proportion of body fat despite reduced calorie intake, and exhibited more compromised glucose metabolism, with major implications for treatment options for genetically obese individuals.

Engineered Interspecies Amino Acid Cross-Feeding Increases Population Evenness in a Synthetic Bacterial Consortium.

In nature, microbes interact antagonistically, neutrally, or beneficially. To shed light on the effects of positive interactions in microbial consortia, we introduced metabolic dependencies and metabolite overproduction into four bacterial species. While antagonistic interactions govern the wild-type consortium behavior, the genetic modifications alleviated antagonistic interactions and resulted in beneficial interactions. Engineered cross-feeding increased population evenness, a component of ecological diversity, in different environments, including in a more complex gnotobiotic mouse gut environment. Our findings suggest that metabolite cross-feeding could be used as a tool for intentionally shaping microbial consortia in complex environments.IMPORTANCE Microbial communities are ubiquitous in nature. Bacterial consortia live in and on our body and in our environment, and more recently, biotechnology is applying microbial consortia for bioproduction. As part of our body, bacterial consortia influence us in health and disease. Microbial consortium function is determined by its composition, which in turn is driven by the interactions between species. Further understanding of microbial interactions will help us in deciphering how consortia function in complex environments and may enable us to modify microbial consortia for health and environmental benefits.

Author Correction: Microbiota therapy acts via a regulatory T cell MyD88/RORγt pathway to suppress food allergy.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Microbiota therapy acts via a regulatory T cell MyD88/RORγt pathway to suppress food allergy.

The role of dysbiosis in food allergy (FA) remains unclear. We found that dysbiotic fecal microbiota in FA infants evolved compositionally over time and failed to protect against FA in mice. Infants and mice with FA had decreased IgA and increased IgE binding to fecal bacteria, indicative of a broader breakdown of oral tolerance than hitherto appreciated. Therapy with Clostridiales species impacted by dysbiosis, either as a consortium or as monotherapy with Subdoligranulum variabile, suppressed FA in mice as did a separate immunomodulatory Bacteroidales consortium. Bacteriotherapy induced expression by regulatory T (Treg) cells of the transcription factor ROR-γt in a MyD88-dependent manner, which was deficient in FA infants and mice and ineffectively induced by their microbiota. Deletion of Myd88 or Rorc in Treg cells abrogated protection by bacteriotherapy. Thus, commensals activate a MyD88/ROR-γt pathway in nascent Treg cells to protect against FA, while dysbiosis impairs this regulatory response to promote disease.

Automated discovery of functional generality of human gene expression programs.

An important research problem in computational biology is the identification of expression programs, sets of co-expressed genes orchestrating normal or pathological processes, and the characterization of the functional breadth of these programs. The use of human expression data compendia for discovery of such programs presents several challenges including cellular inhomogeneity within samples, genetic and environmental variation across samples, uncertainty in the numbers of programs and sample populations, and temporal behavior. We developed GeneProgram, a new unsupervised computational framework based on Hierarchical Dirichlet Processes that addresses each of the above challenges. GeneProgram uses expression data to simultaneously organize tissues into groups and genes into overlapping programs with consistent temporal behavior, to produce maps of expression programs, which are sorted by generality scores that exploit the automatically learned groupings. Using synthetic and real gene expression data, we showed that GeneProgram outperformed several popular expression analysis methods. We applied GeneProgram to a compendium of 62 short time-series gene expression datasets exploring the responses of human cells to infectious agents and immune-modulating molecules. GeneProgram produced a map of 104 expression programs, a substantial number of which were significantly enriched for genes involved in key signaling pathways and/or bound by NF-kappaB transcription factors in genome-wide experiments. Further, GeneProgram discovered expression programs that appear to implicate surprising signaling pathways or receptor types in the response to infection, including Wnt signaling and neurotransmitter receptors. We believe the discovered map of expression programs involved in the response to infection will be useful for guiding future biological experiments; genes from programs with low generality scores might serve as new drug targets that exhibit minimal "cross-talk," and genes from high generality programs may maintain common physiological responses that go awry in disease states. Further, our method is multipurpose, and can be applied readily to novel compendia of biological data.

Quorum Sensing Can Be Repurposed To Promote Information Transfer between Bacteria in the Mammalian Gut.

The gut microbiome is intricately involved with establishing and maintaining the health of the host. Engineering of gut microbes aims to add new functions and expand the scope of control over the gut microbiome. To create systems that can perform increasingly complex tasks in the gut, it is necessary to harness the ability of the bacteria to communicate in the gut environment. Interestingly, acyl-homoserine lactone (acyl-HSL)-mediated Gram-negative bacterial quorum sensing, a widely used mode of intercellular signaling system in nature, has not been identified in normal healthy mammalian gut. It remains unknown whether the gut bacteria that do not natively use quorum sensing can be engineered to successfully signal to other bacteria using acyl-HSLs in the gut environment. Here, we repurposed quorum sensing to create an information transfer system between native gut Escherichia coli and attenuated Salmonella enterica serovar Typhimurium. Specifically, we functionalized one species with inducible signal production and the other with signal detection and recording using genomically integrated circuits. The information transfer system demonstrated successful intra- and interspecies signaling in the murine gut. This study provides a basis for further understanding of interbacterial interactions in an otherwise hard-to-study environment as well as a basis for further investigation of the potential of acyl-HSLs as intercellular signaling molecules of engineered gut consortia.

Computer-guided design of optimal microbial consortia for immune system modulation.

Manipulation of the gut microbiota holds great promise for the treatment of diseases. However, a major challenge is the identification of therapeutically potent microbial consortia that colonize the host effectively while maximizing immunologic outcome. Here, we propose a novel workflow to select optimal immune-inducing consortia from microbiome compositicon and immune effectors measurements. Using published and newly generated microbial and regulatory T-cell (Treg) data from germ-free mice, we estimate the contributions of twelve Clostridia strains with known immune-modulating effect to Treg induction. Combining this with a longitudinal data-constrained ecological model, we predict the ability of every attainable and ecologically stable subconsortium in promoting Treg activation and rank them by the Treg Induction Score (TrIS). Experimental validation of selected consortia indicates a strong and statistically significant correlation between predicted TrIS and measured Treg. We argue that computational indexes, such as the TrIS, are valuable tools for the systematic selection of immune-modulating bacteriotherapeutics.