Prevotella copri and microbiota members mediate the beneficial effects of a therapeutic food for malnutrition.

Microbiota-directed complementary food (MDCF) formulations have been designed to repair the gut communities of malnourished children. A randomized controlled trial demonstrated that one formulation, MDCF-2, improved weight gain in malnourished Bangladeshi children compared to a more calorically dense standard nutritional intervention. Metagenome-assembled genomes from study participants revealed a correlation between ponderal growth and expression of MDCF-2 glycan utilization pathways by Prevotella copri strains. To test this correlation, here we use gnotobiotic mice colonized with defined consortia of age- and ponderal growth-associated gut bacterial strains, with or without P. copri isolates closely matching the metagenome-assembled genomes. Combining gut metagenomics and metatranscriptomics with host single-nucleus RNA sequencing and gut metabolomic analyses, we identify a key role of P. copri in metabolizing MDCF-2 glycans and uncover its interactions with other microbes including Bifidobacterium infantis. P. copri-containing consortia mediated weight gain and modulated energy metabolism within intestinal epithelial cells. Our results reveal structure-function relationships between MDCF-2 and members of the gut microbiota of malnourished children with potential implications for future therapies.

Bioactive glycans in a microbiome-directed food for children with malnutrition.

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Here we analyse biospecimens from a randomized, controlled trial of a microbiome-directed complementary food (MDCF-2) that produced superior rates of weight gain compared with a calorically more dense conventional ready-to-use supplementary food in 12-18-month-old Bangladeshi children with moderate acute malnutrition4. We reconstructed 1,000 bacterial genomes (metagenome-assembled genomes (MAGs)) from the faecal microbiomes of trial participants, identified 75 MAGs of which the abundances were positively associated with ponderal growth (change in weight-for-length Z score (WLZ)), characterized changes in MAG gene expression as a function of treatment type and WLZ response, and quantified carbohydrate structures in MDCF-2 and faeces. The results reveal that two Prevotella copri MAGs that are positively associated with WLZ are the principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilizing the component glycans of MDCF-2. The predicted specificities of carbohydrate-active enzymes expressed by their polysaccharide-utilization loci are correlated with (1) the in vitro growth of Bangladeshi P. copri strains, possessing varying degrees of polysaccharide-utilization loci and genomic conservation with these MAGs, in defined medium containing different purified glycans representative of those in MDCF-2, and (2) the levels of faecal carbohydrate structures in the trial participants. These associations suggest that identifying bioactive glycan structures in MDCFs metabolized by growth-associated bacterial taxa will help to guide recommendations about their use in children with acute malnutrition and enable the development of additional formulations.

Prevotella copri-related effects of a therapeutic food for malnutrition.

Preclinical and clinical studies are providing evidence that the healthy growth of infants and children reflects, in part, healthy development of their gut microbiomes1-5. This process of microbial community assembly and functional maturation is perturbed in children with acute malnutrition. Gnotobiotic animals, colonized with microbial communities from children with severe and moderate acute malnutrition, have been used to develop microbiome-directed complementary food (MDCF) formulations for repairing the microbiomes of these children during the weaning period5. Bangladeshi children with moderate acute malnutrition (MAM) participating in a previously reported 3-month-long randomized controlled clinical study of one such formulation, MDCF-2, exhibited significantly improved weight gain compared to a commonly used nutritional intervention despite the lower caloric density of the MDCF6. Characterizing the 'metagenome assembled genomes' (MAGs) of bacterial strains present in the microbiomes of study participants revealed a significant correlation between accelerated ponderal growth and the expression by two Prevotella copri MAGs of metabolic pathways involved in processing of MDCF-2 glycans1. To provide a direct test of these relationships, we have now performed 'reverse translation' experiments using a gnotobiotic mouse model of mother-to-offspring microbiome transmission. Mice were colonized with defined consortia of age- and ponderal growth-associated gut bacterial strains cultured from Bangladeshi infants/children in the study population, with or without P. copri isolates resembling the MAGs. By combining analyses of microbial community assembly, gene expression and processing of glycan constituents of MDCF-2 with single nucleus RNA-Seq and mass spectrometric analyses of the intestine, we establish a principal role for P. copri in mediating metabolism of MDCF-2 glycans, characterize its interactions with other consortium members including Bifidobacterium longum subsp. infantis, and demonstrate the effects of P. copri-containing consortia in mediating weight gain and modulating the activities of metabolic pathways involved in lipid, amino acid, carbohydrate plus other facets of energy metabolism within epithelial cells positioned at different locations in intestinal crypts and villi. Together, the results provide insights into structure/function relationships between MDCF-2 and members of the gut communities of malnourished children; they also have implications for developing future prebiotic, probiotic and/or synbiotic therapeutics for microbiome restoration in children with already manifest malnutrition, or who are at risk for this pervasive health challenge.

Bioactive glycans in a microbiome-directed food for malnourished children.

Evidence is accumulating that perturbed postnatal development of the gut microbiome contributes to childhood malnutrition1-4. Designing effective microbiome-directed therapeutic foods to repair these perturbations requires knowledge about how food components interact with the microbiome to alter its expressed functions. Here we use biospecimens from a randomized, controlled trial of a microbiome-directed complementary food prototype (MDCF-2) that produced superior rates of weight gain compared to a conventional ready-to-use supplementary food (RUSF) in 12-18-month-old Bangladeshi children with moderate acute malnutrition (MAM)4. We reconstructed 1000 bacterial genomes (metagenome-assembled genomes, MAGs) present in their fecal microbiomes, identified 75 whose abundances were positively associated with weight gain (change in weight-for-length Z score, WLZ), characterized gene expression changes in these MAGs as a function of treatment type and WLZ response, and used mass spectrometry to quantify carbohydrate structures in MDCF-2 and feces. The results reveal treatment-induced changes in expression of carbohydrate metabolic pathways in WLZ-associated MAGs. Comparing participants consuming MDCF-2 versus RUSF, and MDCF-2-treated children in the upper versus lower quartiles of WLZ responses revealed that two Prevotella copri MAGs positively associated with WLZ were principal contributors to MDCF-2-induced expression of metabolic pathways involved in utilization of its component glycans. Moreover, the predicted specificities of carbohydrate active enzymes expressed by polysaccharide utilization loci (PULs) in these two MAGs correlate with the (i) in vitro growth of Bangladeshi P. copri strains, possessing differing degrees of PUL and overall genomic content similarity to these MAGs, cultured in defined medium containing different purified glycans representative of those in MDCF-2, and (ii) levels of carbohydrate structures identified in feces from clinical trial participants. In the accompanying paper5, we use a gnotobiotic mouse model colonized with age- and WLZ-associated bacterial taxa cultured from this study population, and fed diets resembling those consumed by study participants, to directly test the relationship between P. copri, MDCF-2 glycan metabolism, host ponderal growth responses, and intestinal gene expression and metabolism. The ability to identify bioactive glycan structures in MDCFs that are metabolized by growth-associated bacterial taxa will help guide recommendations about use of this MDCF for children with acute malnutrition representing different geographic locales and ages, as well as enable development of bioequivalent, or more efficacious, formulations composed of culturally acceptable and affordable ingredients.

Human Milk Oligosaccharide Utilization in Intestinal Bifidobacteria Is Governed by Global Transcriptional Regulator NagR.

Bifidobacterium longum subsp. infantis is a prevalent beneficial bacterium that colonizes the human neonatal gut and is uniquely adapted to efficiently use human milk oligosaccharides (HMOs) as a carbon and energy source. Multiple studies have focused on characterizing the elements of HMO utilization machinery in B. longum subsp. infantis; however, the regulatory mechanisms governing the expression of these catabolic pathways remain poorly understood. A bioinformatic regulon reconstruction approach used in this study implicated NagR, a transcription factor from the ROK family, as a negative global regulator of gene clusters encoding lacto-N-biose/galacto-N-biose (LNB/GNB), lacto-N-tetraose (LNT), and lacto-N-neotetraose (LNnT) utilization pathways in B. longum subsp. infantis. This conjecture was corroborated by transcriptome profiling upon nagR genetic inactivation and experimental assessment of binding of recombinant NagR to predicted DNA operators. The latter approach also implicated N-acetylglucosamine (GlcNAc), a universal intermediate of LNT and LNnT catabolism, and its phosphorylated derivatives as plausible NagR transcriptional effectors. Reconstruction of NagR regulons in various Bifidobacterium lineages revealed multiple potential regulon expansion events, suggesting evolution from a local regulator of GlcNAc catabolism in ancestral bifidobacteria to a global regulator controlling the utilization of mixtures of GlcNAc-containing host glycans in B. longum subsp. infantis and Bifidobacterium bifidum. IMPORTANCE The predominance of bifidobacteria in the gut of breastfed infants is attributed to the ability of these bacteria to metabolize human milk oligosaccharides (HMOs). Thus, individual HMOs such as lacto-N-tetraose (LNT) and lacto-N-neotetraose (LNnT) are considered promising prebiotics that would stimulate the growth of bifidobacteria and confer multiple health benefits to preterm and malnourished children suffering from impaired (stunted) gut microbiota development. However, the rational selection of HMO-based prebiotics is hampered by the incomplete knowledge of regulatory mechanisms governing HMO utilization in target bifidobacteria. This study describes NagR-mediated transcriptional regulation of LNT and LNnT utilization in Bifidobacterium longum subsp. infantis. The elucidated regulatory network appears optimally adapted to simultaneous utilization of multiple HMOs, providing a rationale to add HMO mixtures (rather than individual components) to infant formulas. The study also provides insights into the evolutionary trajectories of complex regulatory networks controlling carbohydrate metabolism in bifidobacteria.

Priority effects shape the structure of infant-type Bifidobacterium communities on human milk oligosaccharides.

Bifidobacteria are among the first colonizers of the infant gut, and human milk oligosaccharides (HMOs) in breastmilk are instrumental for the formation of a bifidobacteria-rich microbiota. However, little is known about the assembly of bifidobacterial communities. Here, by applying assembly theory to a community of four representative infant-gut associated Bifidobacterium species that employ varied strategies for HMO consumption, we show that arrival order and sugar consumption phenotypes significantly affected community formation. Bifidobacterium bifidum and Bifidobacterium longum subsp. infantis, two avid HMO consumers, dominate through inhibitory priority effects. On the other hand, Bifidobacterium breve, a species with limited HMO-utilization ability, can benefit from facilitative priority effects and dominates by utilizing fucose, an HMO degradant not utilized by the other bifidobacterial species. Analysis of publicly available breastfed infant faecal metagenome data showed that the observed trends for B. breve were consistent with our in vitro data, suggesting that priority effects may have contributed to its dominance. Our study highlights the importance and history dependency of initial community assembly and its implications for the maturation trajectory of the infant gut microbiota.

Microbial liberation of N-methylserotonin from orange fiber in gnotobiotic mice and humans.

Plant fibers in byproduct streams produced by non-harsh food processing methods represent biorepositories of diverse, naturally occurring, and physiologically active biomolecules. To demonstrate one approach for their characterization, mass spectrometry of intestinal contents from gnotobiotic mice, plus in vitro studies, revealed liberation of N-methylserotonin from orange fibers by human gut microbiota members including Bacteroides ovatus. Functional genomic analyses of B. ovatus strains grown under permissive and non-permissive N-methylserotonin "mining" conditions revealed polysaccharide utilization loci that target pectins whose expression correlate with strain-specific liberation of this compound. N-methylserotonin, orally administered to germ-free mice, reduced adiposity, altered liver glycogenesis, shortened gut transit time, and changed expression of genes that regulate circadian rhythm in the liver and colon. In human studies, dose-dependent, orange-fiber-specific fecal accumulation of N-methylserotonin positively correlated with levels of microbiome genes encoding enzymes that digest pectic glycans. Identifying this type of microbial mining activity has potential therapeutic implications.

Bifidobacterium infantis treatment promotes weight gain in Bangladeshi infants with severe acute malnutrition.

Disrupted development of the gut microbiota is a contributing cause of childhood malnutrition. Bifidobacterium longum subspecies infantis is a prominent early colonizer of the infant gut that consumes human milk oligosaccharides (HMOs). We found that the absolute abundance of Bifidobacterium infantis is lower in 3- to 24-month-old Bangladeshi infants with severe acute malnutrition (SAM) compared to their healthy age-matched counterparts. A single-blind, placebo-controlled trial (SYNERGIE) was conducted in 2- to 6-month-old Bangladeshi infants with SAM. A commercial U.S. donor-derived B. infantis strain (EVC001) was administered daily with or without the HMO lacto-N-neotetraose for 28 days. This intervention increased fecal B. infantis abundance in infants with SAM, although to levels still 10- to 100-fold lower than in untreated healthy controls. EVC001 treatment promoted weight gain that was associated with reduced intestinal inflammation markers in infants with SAM. We cultured fecal B. infantis strains from Bangladeshi infants and colonized gnotobiotic mice with these cultured strains. The gnotobiotic mice were fed a diet representative of that consumed by 6-month-old Bangladeshi infants, with or without HMO supplementation. One B. infantis strain, Bg_2D9, expressing two gene clusters involved in uptake and utilization of N-glycans and plant-derived polysaccharides, exhibited superior fitness over EVC001. The fitness advantage of Bg_2D9 was confirmed in a gnotobiotic mouse model of mother-to-infant gut microbiota transmission where dams received a pretreatment fecal community from a SAM infant in the SYNERGIE trial. Whether Bg_2D9 is superior to EVC001 for treating malnourished infants who consume a diet with limited breastmilk requires further clinical testing.

Milk glycan metabolism by intestinal bifidobacteria: insights from comparative genomics.

Bifidobacteria are early colonizers of the human neonatal gut and provide multiple health benefits to the infant, including inhibiting the growth of enteropathogens and modulating the immune system. Certain Bifidobacterium species prevail in the gut of breastfed infants due to the ability of these microorganisms to selectively forage glycans present in human milk, specifically human milk oligosaccharides (HMOs) and N-linked glycans. Therefore, these carbohydrates serve as promising prebiotic dietary supplements to stimulate the growth of bifidobacteria in the guts of children suffering from impaired gut microbiota development. However, the rational formulation of milk glycan-based prebiotics requires a detailed understanding of how bifidobacteria metabolize these carbohydrates. Accumulating biochemical and genomic data suggest that HMO and N-glycan assimilation abilities vary remarkably within the Bifidobacterium genus, both at the species and strain levels. This review focuses on the delineation and genome-based comparative analysis of differences in respective biochemical pathways, transport systems, and associated transcriptional regulatory networks, providing a foundation for genomics-based projection of milk glycan utilization capabilities across a rapidly growing number of sequenced bifidobacterial genomes and metagenomic datasets. This analysis also highlights remaining knowledge gaps and suggests directions for future studies to optimize the formulation of milk-glycan-based prebiotics that target bifidobacteria.

Transcriptional Regulation of Plant Biomass Degradation and Carbohydrate Utilization Genes in the Extreme Thermophile Caldicellulosiruptor bescii.

Extremely thermophilic bacteria from the genus Caldicellulosiruptor can degrade polysaccharide components of plant cell walls and subsequently utilize the constituting mono- and oligosaccharides. Through metabolic engineering, ethanol and other industrially important end products can be produced. Previous experimental studies identified a variety of carbohydrate-active enzymes in model species Caldicellulosiruptor saccharolyticus and Caldicellulosiruptor bescii, while prior transcriptomic experiments identified their putative carbohydrate uptake transporters. We investigated the mechanisms of transcriptional regulation of carbohydrate utilization genes using a comparative genomics approach applied to 14 Caldicellulosiruptor species. The reconstruction of carbohydrate utilization regulatory network includes the predicted binding sites for 34 mostly local regulators and point to the regulatory mechanisms controlling expression of genes involved in degradation of plant biomass. The Rex and CggR regulons control the central glycolytic and primary redox reactions. The identified transcription factor binding sites and regulons were validated with transcriptomic and transcription start site experimental data for C. bescii grown on cellulose, cellobiose, glucose, xylan, and xylose. The XylR and XynR regulons control xylan-induced transcriptional response of genes involved in degradation of xylan and xylose utilization. The reconstructed regulons informed the carbohydrate utilization reconstruction analysis and improved functional annotations of 51 transporters and 11 catabolic enzymes. Using gene deletion, we confirmed that the shared ATPase component MsmK is essential for growth on oligo- and polysaccharides but not for the utilization of monosaccharides. By elucidating the carbohydrate utilization framework in C. bescii, strategies for metabolic engineering can be pursued to optimize yields of bio-based fuels and chemicals from lignocellulose. IMPORTANCE To develop functional metabolic engineering platforms for nonmodel microorganisms, a comprehensive understanding of the physiological and metabolic characteristics is critical. Caldicellulosiruptor bescii and other species in this genus have untapped potential for conversion of unpretreated plant biomass into industrial fuels and chemicals. The highly interactive and complex machinery used by C. bescii to acquire and process complex carbohydrates contained in lignocellulose was elucidated here to complement related efforts to develop a metabolic engineering platform with this bacterium. Guided by the findings here, a clearer picture of how C. bescii natively drives carbohydrate utilization is provided and strategies to engineer this bacterium for optimal conversion of lignocellulose to commercial products emerge.

Transcriptional control of central carbon metabolic flux in Bifidobacteria by two functionally similar, yet distinct LacI-type regulators.

Bifidobacteria resident in the gastrointestinal tract (GIT) are subject to constantly changing environmental conditions, which require rapid adjustments in gene expression. Here, we show that two predicted LacI-type transcription factors (TFs), designated AraQ and MalR1, are involved in regulating the central, carbohydrate-associated metabolic pathway (the so-called phosphoketolase pathway or bifid shunt) of the gut commensal Bifidobacterium breve UCC2003. These TFs appear to not only control transcription of genes involved in the bifid shunt and each other, but also seem to commonly and directly affect transcription of other TF-encoding genes, as well as genes related to uptake and metabolism of various carbohydrates. This complex and interactive network of AraQ/MalR1-mediated gene regulation provides previously unknown insights into the governance of carbon metabolism in bifidobacteria.

A sparse covarying unit that describes healthy and impaired human gut microbiota development.

Characterizing the organization of the human gut microbiota is a formidable challenge given the number of possible interactions between its components. Using a statistical approach initially applied to financial markets, we measured temporally conserved covariance among bacterial taxa in the microbiota of healthy members of a Bangladeshi birth cohort sampled from 1 to 60 months of age. The results revealed an "ecogroup" of 15 covarying bacterial taxa that provide a concise description of microbiota development in healthy children from this and other low-income countries, and a means for monitoring community repair in undernourished children treated with therapeutic foods. Features of ecogroup population dynamics were recapitulated in gnotobiotic piglets as they transitioned from exclusive milk feeding to a fully weaned state consuming a representative Bangladeshi diet.

Effects of microbiota-directed foods in gnotobiotic animals and undernourished children.

To examine the contributions of impaired gut microbial community development to childhood undernutrition, we combined metabolomic and proteomic analyses of plasma samples with metagenomic analyses of fecal samples to characterize the biological state of Bangladeshi children with severe acute malnutrition (SAM) as they transitioned, after standard treatment, to moderate acute malnutrition (MAM) with persistent microbiota immaturity. Host and microbial effects of microbiota-directed complementary food (MDCF) prototypes targeting weaning-phase bacterial taxa underrepresented in SAM and MAM microbiota were characterized in gnotobiotic mice and gnotobiotic piglets colonized with age- and growth-discriminatory bacteria. A randomized, double-blind controlled feeding study identified a lead MDCF that changes the abundances of targeted bacteria and increases plasma biomarkers and mediators of growth, bone formation, neurodevelopment, and immune function in children with MAM.

Micronutrient Requirements and Sharing Capabilities of the Human Gut Microbiome.

The human gut microbiome harbors a diverse array of metabolic pathways contributing to its development and homeostasis via a complex web of diet-dependent metabolic interactions within the microbial community and host. Genomics-based reconstruction and predictive modeling of these interactions would provide a framework for diagnostics and treatment of dysbiosis-related syndromes via rational selection of therapeutic prebiotics and dietary nutrients. Of particular interest are micronutrients, such as B-group vitamins, precursors of indispensable metabolic cofactors, that are produced de novo by some gut bacteria (prototrophs) but must be provided exogenously in the diet for many other bacterial species (auxotrophs) as well as for the mammalian host. Cross-feeding of B vitamins between prototrophic and auxotrophic species is expected to strongly contribute to the homeostasis of microbial communities in the distal gut given the efficient absorption of dietary vitamins in the upper gastrointestinal tract. To confidently estimate the balance of microbiome micronutrient biosynthetic capabilities and requirements using available genomic data, we have performed a subsystems-based reconstruction of biogenesis, salvage and uptake for eight B vitamins (B1, B2, B3, B5, B6, B7, B9, and B12) and queuosine (essential factor in tRNA modification) over a reference set of 2,228 bacterial genomes representing 690 cultured species of the human gastrointestinal microbiota. This allowed us to classify the studied organisms with respect to their pathway variants and infer their prototrophic vs. auxotrophic phenotypes. In addition to canonical vitamin pathways, several conserved partial pathways were identified pointing to alternative routes of syntrophic metabolism and expanding a microbial vitamin "menu" by several pathway intermediates (vitamers) such as thiazole, quinolinate, dethiobiotin, pantoate. A cross-species comparison was applied to assess the extent of conservation of vitamin phenotypes at distinct taxonomic levels (from strains to families). The obtained reference collection combined with 16S rRNA gene-based phylogenetic profiles was used to deduce phenotype profiles of the human gut microbiota across in two large cohorts. This analysis provided the first estimate of B-vitamin requirements, production and sharing capabilities in the human gut microbiome establishing predictive phenotype profiling as a new approach to classification of microbiome samples. Future expansion of our reference genomic collection of metabolic phenotypes will allow further improvement in coverage and accuracy of predictive phenotype profiling of the human microbiome.

Comparative Genomics Reveals the Regulatory Complexity of Bifidobacterial Arabinose and Arabino-Oligosaccharide Utilization.

Members of the genus Bifidobacterium are common inhabitants of the human gastrointestinal tract. Previously it was shown that arabino-oligosaccharides (AOS) might act as prebiotics and stimulate the bifidobacterial growth in the gut. However, despite the rapid accumulation of genomic data, the precise mechanisms by which these sugars are utilized and associated transcription control still remain unclear. In the current study, we used a comparative genomic approach to reconstruct arabinose and AOS utilization pathways in over 40 bacterial species belonging to the Bifidobacteriaceae family. The results indicate that the gene repertoire involved in the catabolism of these sugars is highly diverse, and even phylogenetically close species may differ in their utilization capabilities. Using bioinformatics analysis we identified potential DNA-binding motifs and reconstructed putative regulons for the arabinose and AOS utilization genes in the Bifidobacteriaceae genomes. Six LacI-family transcriptional factors (named AbfR, AauR, AauU1, AauU2, BauR1 and BauR2) and a TetR-family regulator (XsaR) presumably act as local repressors for AOS utilization genes encoding various α- or ß-L-arabinofuranosidases and predicted AOS transporters. The ROK-family regulator AraU and the LacI-family regulator AraQ control adjacent operons encoding putative arabinose transporters and catabolic enzymes, respectively. However, the AraQ regulator is universally present in all Bifidobacterium species including those lacking the arabinose catabolic genes araBDA, suggesting its control of other genes. Comparative genomic analyses of prospective AraQ-binding sites allowed the reconstruction of AraQ regulons and a proposed binary repression/activation mechanism. The conserved core of reconstructed AraQ regulons in bifidobacteria includes araBDA, as well as genes from the central glycolytic and fermentation pathways (pyk, eno, gap, tkt, tal, galM, ldh). The current study expands the range of genes involved in bifidobacterial arabinose/AOS utilization and demonstrates considerable variations in associated metabolic pathways and regulons. Detailed comparative and phylogenetic analyses allowed us to hypothesize how the identified reconstructed regulons evolved in bifidobacteria. Our findings may help to improve carbohydrate catabolic phenotype prediction and metabolic modeling, while it may also facilitate rational development of novel prebiotics.