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

Hao-Wei Chang, Evan Lee, Yi Wang, Cyrus Zhou, Kali Pruss, Suzanne Henrissat, Robert Chen, Clara Kao, Matthew Hibberd, Hannah Lynn, et al.

► To cite this version:

Hao-Wei Chang, Evan Lee, Yi Wang, Cyrus Zhou, Kali Pruss, et al.. Prevotella copri and microbiota members mediate the beneficial effects of a therapeutic food for malnutrition. Nature Microbiology, 2024, 9 (4), pp.922-937. 10.1038/s41564-024-01628-7. hal-04579760

HAL Id: hal-04579760

https://amu.hal.science/hal-04579760v1

Submitted on 6 Sep 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

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

Received: 20 November 2023

Accepted: 31 January 2024

Published online: 19 March 2024

<figure><img src="image_3.png" /></figure>

Check for updates

Hao-Wei Chang Ⓓ ¹,²,³,¹⁴, Evan M. Lee¹,²,¹⁴, Yi Wang¹,²,³,¹⁴, Cyrus Zhou Ⓓ ¹,²,
Kali M. Pruss Ⓓ ¹,²,³, Suzanne Henrissat¹,²,⁴, Robert Y. Chen¹,², Clara Kao Ⓓ ¹,²,

Matthew C. Hibberd Ⓓ ¹,²,³, Hannah M. Lynn Ⓓ ¹,², Daniel M. Webber Ⓓ ¹,²,³,

Marie Crane Ⓓ ¹,², Jiye Cheng Ⓓ ¹,², Dmitry A. Rodionov⁵,

Aleksandr A. Arzamasov Ⓓ⁵, Juan J. Castillo Ⓓ⁶, Garret Couture⁶, Ye Chen Ⓓ ¹,²,⁶,

Nikita P. Balcazo Jr. Ⓓ⁶, Carlito B. Lebrilla Ⓓ⁶, Nicolas Terrapon Ⓓ⁴,

Bernard Henrissat Ⓓ ⁷,⁸, Olga Ilkayeva Ⓓ ⁹,¹⁰,¹¹, Michael J. Muehlbauer⁹,¹⁰,

Christopher B. Newgard Ⓓ ⁹,¹⁰,¹¹,¹², Ishita Mostafa¹³, Subhasish Das¹³,

Mustafa Mahfuz¹³, Andrei L. Osterman Ⓓ⁵, Michael J. Barratt¹,²,³,

Tahmeed Ahmed¹³ & Jeffrey I. Gordon Ⓓ ¹,²,³ ☑

<figure><img src="image_4.png" /></figure>

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.

Microbial colonization of the infant gut begins at birth in a process that is influenced by mode of delivery, exposure to maternal microbes from various sources (vaginal, skin, faecal and breast milk), environmental microbes and first foods¹,². By the third postnatal year, the process of

microbial community assembly ('maturation') is largely complete, with communities from healthy children at this age showing configurations that resemble those of adult members of the same family²⁻⁴. By contrast, undernourished children exhibit delayed development of their

<figure><img src="image_5.png" /></figure>
microbiota; transplantation of these 'immature' communities into germ-free mice produces impairments in growth and metabolism in recipient animals⁵,⁶. In preclinical models, perturbations in the small intestinal microbiota and protein-deficient diets have both been shown to produce enteropathies that share histologic and pathologic features of environmental enteric dysfunction, a chronic condition affecting undernourished children that is characterized by gut/systemic inflammation and impaired nutrient absorption⁷⁻⁹. Taken together, these findings highlight the role of the intestinal microbiota and its collection of genes (microbiome) in fostering healthy postnatal growth.

We previously performed a genome-resolved metagenomic analysis of faecal samples serially collected from 12- to 18-month-old Bangladeshi children with moderate acute malnutrition enrolled in a randomized controlled clinical trial. The trial compared the effects of administering a microbiota-directed complementary food prototype (MDCF-2) versus a ready-to-use supplementary food (RUSF) on host physiology and the composition and expressed functions of the childrens' microbiomes¹⁰,¹¹. During the 3 month intervention, MDCF-2 produced significantly greater increases in ponderal growth (defined by rate of change in weight for length (height), expressed as WLZ scores), even though MDCF-2 had a 15% lower caloric density. MDCF-2 also promoted significantly greater increases in levels of plasma proteins positively associated with WLZ, including biomarkers and mediators of musculoskeletal and central nervous system development¹⁰. We also identified metagenome-assembled genomes (MAGs) whose abundances were associated with WLZ scores. Cel-lulose, galactan, arabinan, xylan and mannan represent the principal non-starch polysaccharides in MDCF-2 (ref. 11). Among the 75 MAGs found to be positively correlated with WLZ, two *Prevotella copri* MAGs dominated the expression of carbohydrate utilization pathways that target MDCF-2 polysaccharides. The expression of these pathways and the faecal concentration of glycan metabolism products were significantly positively correlated with the magnitude of the childrens' improvement in WLZ¹¹. The two WLZ-associated *P. copri* MAGs (Bg0018 and Bg0019) share ten functionally conserved polysaccharide utilization loci (PULs), including seven that are completely conserved. The degree of representation of these seven PULs among the 11 *P. copri* MAGs identified in study participants was predictive of each MAG's strength of association with WLZ.

*Bifidobacterium longum* subsp. *infantis* (hereafter referred to as *B. infantis*) is a prominent early colonizer of the infant gut and a principal consumer of human milk oligosaccharides¹². Studies of healthy versus malnourished Bangladeshi children showed that *B. infantis* is depleted or absent in the microbiota of infants/children with severe acute malnutrition (SAM). A randomized controlled clinical trial showed that administration of a commercial *B. infantis* strain to infants with SAM improved ponderal growth and reduced levels of biomarkers of gut inflammation¹³. Follow-up preclinical studies revealed that the combination of a Bangladeshi *B. infantis* strain (Bg2D9) and a commercial *B. infantis* strain from a US donor promoted weight gain when introduced into gnotobiotic mice colonized with a pretreatment uncultured faecal microbiota from a Bangladeshi infant who had been a participant in the clinical trial¹³.

In this Article, we use gnotobiotic mice to further examine the role of *P. copri* and other community members, including *B. infantis*, in mediating MDCF-2 metabolism, as well as the effects of MDCF-2 on host physiology, notably weight gain and intestinal function. Defined collections of genome-sequenced bacterial strains, cultured from Bangladeshi children, were sequentially introduced into gnotobiotic female mice (dams), with subsequent transmission of these strains to their pups during the suckling and weaning periods. The strains represented microorganisms whose prominence changes at different stages of postnatal gut microbial community assembly in healthy Bangladeshi children ('age-discriminatory' taxa)⁴⁻⁶,¹⁴. They also included WLZ-correlated bacterial strains selected based on their

shared features with MAGs identified in the MDCF-2 clinical study¹⁰,¹¹. Pups were subjected to the same dietary sequence of exclusive milk feeding (from the dam) followed by weaning onto an MDCF-2 supplemented diet. As well as characterizing changes in microbial community composition and function in response to MDCF-2, we also looked at facets of host physiology. Substantial metabolic investments are needed to support the normal daily replacement of large numbers of gut epithelial cells¹⁵ as well as the functions they normally express in their differentiating and differentiated states. Therefore, we tested the hypothesis that epithelial cell gene expression and metabolism would be sensitive reporters of functional differences associated with defined consortia with or without *P. copri* MAG-representing strains. We report a central role for *P. copri* strains closely resembling the WLZ-associated MAGs in metabolizing glycans present in MDCF-2, plus their capacity, in the context of other microbiota members such as *B. longum* subsp. *infantis* and MDCF-2, to promote weight gain and influence expression of metabolic functions in enterocytes.

## Results

### Dam-pup transmission of age- and WLZ-associated bacteria

We first designed a defined human gut microbial community that reflected the developing gut microbiota of children enrolled in the clinical study¹¹. We selected 20 bacterial strains, 16 of which were cultured from the faecal microbiota of 6- to 24-month-old Bangladeshi children living in Mirpur, the urban slum where the previously reported randomized controlled MDCF-2 clinical trial had been performed (Supplementary Table 1a). They included strains initially identified by the close correspondence of their 16S ribosomal RNA gene sequences to (1) a group of taxa that reflect gut microbiota development in healthy Bangladeshi children⁴,¹⁴ and (2) taxa whose abundances had statistically significant associations (positive or negative) with the rate of weight gain (β-WLZ)⁶,¹⁰. The relatedness of these strains to the 1,000 MAGs assembled from faecal samples obtained from all participants in the clinical study¹¹ was determined by average nucleotide sequence identity (ANI) scores, alignment coverage parameters¹⁶,¹⁷ and their encoded metabolic pathways (Supplementary Table 1b–d). Encoded metabolic pathways for carbohydrate utilization, amino acid and vitamin/cofactor biosynthesis, and fermentation in MAGs and cultured strains were defined by in silico reconstructions; the results are described in the form of ‘binary phenotype’ scores denoting pathway presence or absence¹¹.

To test how *P. copri* colonization with a strain resembling the two WLZ-associated MAGs, Bg0018 and Bg0019, affected microbial community composition and expressed functions, dietary glycan degradation and host metabolism, we selected the Bangladeshi *P. copri* strain PS131.S11 (abbreviated *P. copri* Bg131). This strain was chosen because of its phylogenetic similarity to Bg0018 and Bg0019 (Extended Data Fig. 1a), the concordance of its metabolic pathway representation with these MAGs (Supplementary Table 1c,d) and its representation of five of the ten functionally conserved PULs shared by Bg0018 and Bg0019. These five PULs are predicted to be involved in degradation of starch, β-glucan, pectin, pectic galactan and xylan (Supplementary Table 1e,f). An additional arabinogalactan-targeting PUL was found adjacent to a conserved PUL-targeting starch, although it did not meet criteria for conservation with the corresponding PUL in MAGs Bg0018 and Bg0019 (Supplementary Table 1f).

To assess the specificity of responses of *P. copri* to MDCF-2, we additionally included an isolate from another *Prevotella* species, *Prevotella stercorea*. No *P. stercorea* MAGs were associated with WLZ in the clinical study, and the cultured isolate did not share any of the PULs present in MAGs Bg0018 and Bg0019 or *P. copri* Bg131. Instead, the PULs in the *P. stercorea* isolate contain glycoside hydrolases that mainly target animal-derived glycans (Supplementary Table 1g). Therefore, we hypothesized this isolate would show lower fitness on the MDCF-2 plant glycan-based diet.
Initial attempts to mono-colonize mice with *P. copri* revealed that it was a poor colonizer on its own (Extended Data Fig. 1b,c). To ensure that *B. infantis* was well represented at the earliest stages of assembly of the defined community so that later colonizers such as *P. copri* could establish themselves, the collection of cultured isolates included two strains of *B. infantis* recovered from Bangladeshi children–*B. infantis* Bg463 and *B. infantis* Bg2D9. The Bg463 strain had been used in our earlier preclinical studies that led to the development of MDCF-2<sup>6,14</sup>. *B. infantis* Bg2D9 had shown greater fitness (absolute abundance) than Bg463 when administered to just-weaned germ-free mice consuming a diet representative of that consumed by 6-month-old children living in Mirpur<sup>13</sup>; this superior fitness was attributed to additional carbohydrate utilization pathways that the strain possesses<sup>13</sup>.

We used the 20-strain collection to perform a three-arm, fixed diet study that involved 'successive' waves of maternal colonization with four different bacterial consortia (Fig. 1a–c). The sequence of introduction of taxa into dams was designed to emulate temporal features of the normal postnatal development of the human gut community; for example, consortia 1 and 2 consisted of strains that are prominent colonizers of healthy infants/children in the first postnatal year (including the *B. infantis* isolates), while those in consortium 3 are prominent during weaning in the second postnatal year<sup>4–6,14</sup>. This dam-to-pup colonization strategy also helped overcome the technical challenge of reliable delivery of bacterial consortia to newborn pups via oral gavage.

Dually housed germ-free dams were switched from a standard breeder chow to a 'weaning-diet' supplemented with MDCF-2 on postpartum day 2, 2 days before initiation of the colonization sequence. This weaning diet was formulated to emulate the diets consumed by children in the clinical trial during MDCF-2 treatment (Methods; Fig. 1a and Supplementary Table 2). Pups in all three arms were subjected to a diet sequence that began with exclusive milk feeding (from the nursing dam) followed by a weaning period where pups had access to the weaning phase diet supplemented with MDCF-2. Pups were weaned at postnatal day 24 (P24), after which time they received MDCF-2 alone ad libitum until P53 when they were euthanized.

To identify factors that promote stable *P. copri* colonization, we compared the effects of the two different *B. infantis* strains, Bg2D9 and Bg463, on *P. copri* and the other cultured strains representing age- and WLZ-associated MAGs. We included *B. infantis* Bg2D9 in consortium 1 of arm 1 of the three-arm experiment and *B. infantis* Bg463 in consortium 1 of arm 2 (Fig. 1b,c). This consortium consisted of five 'early' infant gut community colonizers and was administered to dams on postpartum day 4. Dams in these two arms subsequently received the following gavages: (1) on postpartum day 7, *P. copri* and *P. stercorea*; (2) on postpartum day 10, additional age-discriminatory and WLZ-associated taxa and (3) on postpartum day 21, three strains–*P. copri*, *P. stercorea* and *Faecalibacterium prausnitzii* (Fig. 1b). At this last time point, the three strains were given by oral gavage to both the dams and their offspring to help promote successful colonization. To test the effects of *P. copri* colonization, we included arm 3, which was a replicate of arm 2 but without *Prevotella* in consortia 2 and 4. Our rationale for the timing of the first three gavages was based on the diet sequence (gavage 1 of early colonizers at a time (P4) when mice were exclusively consuming

the dam's milk, gavage 2 as the pups were just beginning to consume the human weaning (complementary food) diet, gavage 3 later during this period of 'complementary feeding' and the fourth gavage to help to ensure a consistent level of *P. copri* colonization at the end of weaning (and subsequently through the post-weaning period)).

## Effects of *B. infantis* on host weight and microbiota

Gnotobiotic mice colonized with *B. infantis* Bg2D9 exhibited a significantly greater increase in weight gain between P23 (the first time point measured, 2 days after the final gavage) and P53 compared to mice in the two other experimental arms (Fig. 1d and Supplementary Table 3). By contrast, there was no significant difference in weight gain between animals colonized with *B. infantis* Bg463 with or without *Prevotella* species in arms 2 and 3.

To explore whether these differences in weight gain were associated with differences in microbial community composition, we used shotgun sequencing of community DNA to quantify the relative and absolute abundances of all administered strains in faecal samples collected from dams on postpartum days 21, 24 and 35, faecal samples collected from their offspring on P21, P24, P35 and P53, and caecal and ileal contents collected from offspring at P53 (n = 2 dams and 5–8 pups analysed per arm; Supplementary Table 4a). Colonization of each consortium member was highly consistent among all animals in each treatment group (Supplementary Table 4b–d). *B. infantis* Bg2D9 successfully colonized pups at P21 in arm 1. By contrast, *B. infantis* Bg463 colonized at 5–8 orders-of-magnitude lower absolute abundance levels in arms 2 and 3 (Supplementary Table 4c). These differences between the groups were sustained through P53. Consistent with the role of *B. infantis* as an early pre-weaning colonizer, both strains decreased in abundance between P21 and P53 (Fig. 1e). Exposure to *B. infantis* Bg2D9 in Arm 1 was associated with levels of *P. copri* colonization that were 3 orders of magnitude greater in the pre-weaning period (P21) than in arm 2 mice exposed to *B. infantis* Bg463 (Fig. 1f and Supplementary Table 4c). Administering the fourth gavage on P21 elevated the absolute abundance of faecal *P. copri* in arm 2 to a level comparable to arm 1; this level was sustained throughout the post-weaning period (P24 to P53) (Fig. 1f and Supplementary Table 4c). This effect of the fourth gavage was also evident in the ileal and caecal microbiota at the time of euthanasia (Supplementary Table 4d).

Based on these results, we directly tested the colonization dependency of *P. copri* on *B. infantis* in two independent experiments whose designs are outlined in Extended Data Fig. 1b. Dually housed germ-free dams were switched from standard breeder chow to the weaning Bangladeshi diet supplemented with MDCF-2 on postpartum day 2. On postpartum day 4, one group of dams was colonized with *B. infantis* Bg2D9, while the other group received a sham gavage. On postpartum days 7 and 10, both groups of gnotobiotic mice were gavaged with a consortium containing five *P. copri* strains. These five *P. copri* strains (Bg131, 1A8, 2C6, 2D7 and G8) were all isolated from faecal samples obtained from Bangladeshi children (Supplementary Table 5). Pups were separated from their dams at the completion of weaning, and their diet was switched to MDCF-2 until euthanasia on P42 (n = 9 mice per treatment group; 2 independent control experiments). At this time

**Fig. 1 |** Identifying factors that affect the efficiency of colonization of gnotobiotic dam–pup dyads with *P. copri* in the presence of other cultured age-discriminatory and WLZ-associated bacterial strains and the effects of colonization on pup weight gain. a, Energy contribution from different modules of the 'weaning diet supplemented with MDCF-2'. b,c, Study design (n = 2 dams and 8, 5 and 7 offspring for arms 1, 2 and 3, respectively). b, The timing of bacterial colonization of dams and diet switches. c, The gavages administered to members of each treatment arm. d, Body weights of the offspring of dams, normalized to P23. e, Absolute abundance of *B. infantis* Bg2D9 (arm 1) and *B. infantis* Bg463 (arm 2) in faecal samples obtained from pups. f, Absolute abundance of *P. copri* in faecal samples collected from pups in the indicated

treatment arms at the indicated postnatal time points. Inset: the absolute abundance of *P. copri* in faecal samples collected from pups at P21. g, Principal components analysis of absolute abundances of other community members in faecal samples obtained from pups at P21 and P53. Mean values ± s.d. are shown in d–f. Each dot in d–f represents an individual animal. P values were calculated using a linear mixed-effects model (d, Methods), a two-sided Mann–Whitney U test (f, inset) or PERMANOVA (g). Centroids are denoted by a coloured 'X'. Shaded ellipses represent the 95% confidence interval of the sample distribution. Each dot represents an individual animal. Data generated from all of the offspring were used in the analyses shown in d–g.
a

<figure><img src="image_6.png" /></figure>

b

<figure><img src="image_7.png" /></figure>

c

<table><thead><tr><th rowspan="2">Consortium</th><th>Arm 1</th><th>Arm 2</th><th>Arm 3</th></tr><tr><th>With <i>B. infantis</i> Bg2D9</th><th>With <i>B. infantis</i> Bg463</th><th>With <i>B. infantis</i> Bg463</th></tr></thead><tbody><tr><td>1</td><td>+</td><td>+</td><td>-</td></tr><tr><td>2</td><td>+</td><td>+</td><td>+</td></tr><tr><td>3</td><td>+</td><td>+</td><td>+</td></tr><tr><td>4</td><td>+</td><td>+</td><td>+(but without <i>Prevotella</i> sp.)</td></tr></tbody></table>

d

<figure><img src="image_8.png" /></figure>

e

<figure><img src="image_9.png" /></figure>

f

<figure><img src="image_10.png" /></figure>

g

<figure><img src="image_11.png" /></figure>
point, the absolute abundance of *P. copri* in faeces collected from mice that had received *B. infantis* Bg2D9 was 3 orders of magnitude higher than in animals never exposed to *B. infantis* (Extended Data Fig. 1c and Supplementary Table 6). There was no statistically significant difference in weight gain from P23 to P42 between the mono- and bi-colonization groups, although interpreting this finding is at least partially confounded by the increased caecal size and fluid content observed in animals mono-colonized with *P. copri*.

The effects of *B. infantis* on *P. copri* did not generalize to *P. stercorea*. Unlike *P. copri*, the absolute abundance of *P. stercorea* in faeces sampled on P21 and P24 was not significantly different in mice belonging to arms 1 and 2 (Supplementary Table 4c). Before weaning at P24, the absolute abundance of *P. stercorea* was 5 orders of magnitude lower than that of *P. copri*. Throughout the post-weaning period, the absolute abundance of *P. stercorea* remained similar in members of both treatment arms but 2 orders of magnitude below that of *P. copri*.

Over the course of the experiment, *B. infantis* Bg2D9 colonization resulted in a faecal community composition that was distinct from that of animals colonized with *B. infantis* Bg463, regardless of their *Prevotella* colonization (Fig. 1g). In animals colonized with *Prevotella* and either of the two *B. infantis* strains (arm 1 versus arm 2), these differences were observed as early as P21 and became more pronounced by the end of the experiment at P53 (Fig. 1g). In animals harbouring *Prevotella*-containing communities, colonization with *B. infantis* Bg2D9 compared to Bg463 significantly increased the fitness of three organisms in the P53 faecal community (*B. luti*, *D. longicatena*, *M. multacida*) and five organisms in the P53 caecal community (*B. breve*, *B. catenulatum*, *B. obeum*, *D. longicatena*, *M. multacida*), while the absolute abundances of the other community members were not significantly affected (Extended Data Fig. 2a–c and Supplementary Table 4c,d).

In animals colonized with *B. infantis* Bg463, the addition of *Prevotella* to the community did not result in significant differences in community composition at either P21 or P53 and only significantly increased the fitness of one organism (*Streptococcus gallolyticus*; Fig. 1g and Extended Data Fig. 2a). However, when comparing arms 1 and 3, the combination of *B. infantis* Bg2D9 and *Prevotella* colonization increased the fitness of a larger set of seven and six organisms in P53 faecal and caecal communities, respectively, including *B. catenulatum*, *Blautia obeum* and *Mitsuokella multacida*—three of the four organisms predicted to be capable of utilizing arabinose (Extended Data Fig. 2a,b); this suggests a potential synergistic interaction between the two organisms in mediating effects on community structure. Based on these results, we concluded that in the context of this preclinical model, (1) *B. infantis* Bg2D9 colonization was an important determinant of microbial community structure, including the fitness of *P. copri*; (2) communities containing *B. infantis* Bg2D9 were associated with augmented weight gain and (3) the temporal profile of community member fitness produced when *B. infantis* Bg2D9 was included more closely resembled that of children in the clinical study who, during the weaning period when MDCF-2 treatment was initiated, all had substantial levels of *P. copri*¹¹.

## Microbial metabolism of MDCF-2 glycans

Given these observed differences in microbial community structure, we used ultra-high performance liquid chromatography-triple quadrupole mass spectrometric (UHPLC-QqQ-MS)-based measurements of mono-saccharide and linkage content of glycans to analyse the metabolism of MDCF-2 across treatment groups. We sampled the caecum because we wanted to compare microbial gene expression with polysaccharide degradative capacity in a large gut habitat specialized for microbial fermentation¹⁸.

While *B. infantis* Bg2D9 colonization was an important determinant of microbial community structure, *Prevotella* colonization drove the degradation of MDCF-2 glycans, regardless of *B. infantis*

colonization (Extended Data Fig. 3). *Prevotella* colonization significantly reduced the levels of arabinose in caecal glycans (Extended Data Fig. 3a) and levels of the arabinose-containing linkages t-Arap, t-Araf, 2-Araf, 2,3-Araf and 3,4-Xylp/3,5-Araf (Extended Data Fig. 3b). These differences are supported by the fact that *P. copri* Bg131, unlike *P. stercorea*, contains PULs involved in processing arabinose-containing MDCF-2 glycans: that is, PUL27b specifies carbohydrate active enzymes (CAZymes) known or predicted to digest arabinogalactan, while PUL2 possesses a fucosidase that could target the terminal residues found in arabinogalactan II (Supplementary Table 1f). By contrast, there were no significant differences between animals colonized with *B. infantis* Bg2D9 versus *B. infantis* Bg463 in arms 1 and 2 for any of the monosaccharides or linkages measured (Supplementary Table 7). Together, these results indicate that *Prevotella*-containing communities show more complete degradation of branched arabinans and a greater degree of liberation of arabinose from MDCF-2 glycans.

To investigate how these changes in glycan utilization are associated with the expressed metabolic functions of community members, we performed microbial RNA sequencing (RNA-seq) on caecal contents (Supplementary Table 8). We first analysed the expression levels of PULs by *P. copri* and *P. stercorea*. In both *Prevotella*-containing arms (arms 1 and 2), *P. copri* PULs with predicted targets of starch and arabinoxalactan (PUL 27a and 27b, respectively) were the most significantly enriched for higher levels of expression by gene set enrichment analysis (GSEA) (Methods) (Extended Data Fig. 4a and Supplementary Table 8b).

In contrast to *P. copri*, only two *P. stercorea* PULs with predicted targets of α-mannan and N-linked glycans (PULs 5 and 7, respectively) were significantly enriched for higher expression (Extended Data Fig. 4b and Supplementary Table 8b). These two PULs had lower levels of expression than the *P. copri* PULs. Furthermore, unlike the observed reductions in the arabinose content of caecal contents collected from mice colonized with the *Prevotella*-containing consortia, there were no significant differences in mannose, N-acetylglucosamine or N-acetylgalactosamine, the primary components of glycans targeted by these *P. stercorea* PULs (Extended Data Fig. 4c). These results indicate that *P. copri*, not *P. stercorea*, is responsible for increased liberation of arabinan from MDCF-2 and contributes to degradation of polysaccharides represented in MDCF-2.

We then turned to the rest of the community to determine whether these changes in glycan metabolism had effects on the metabolic responses of other organisms. The *B. infantis* Bg2D9 transcriptome was distinct from *B. infantis* Bg463 (Extended Data Fig. 3c). Consistent with the UHPLC-QqQ-MS-based analysis of caecal glycans, arabinose utilization was among the most upregulated pathways in *B. catenulatum*, *M. multacida* and *B. obeum*—three of the four organisms predicted to be capable of utilizing arabinose (Extended Data Fig. 3d and Supplementary Table 8c,d). These three organisms were also significantly more abundant with *B. infantis* Bg2D9 colonization (Extended Data Fig. 2b). *B. catenulatum* upregulated arabinose utilization genes directly in response to *P. copri* colonization (arm 2 versus arm 3; Extended Data Fig. 3d). By contrast, *M. multacida* and *B. obeum* showed upregulation of arabinose utilization in response to *B. infantis* Bg2D9 colonization (arm 1 versus arm 2; Extended Data Fig. 3d). *M. multacida* and *B. obeum* also showed significant upregulation of almost all their genes involved in biosynthesis of the branched chain amino acids as well as glutamate and glutamine with *B. infantis* Bg2D9 colonization (Extended Data Fig. 3d and Supplementary Table 8d). Together, these findings suggest that while *B. infantis* Bg2D9 colonization affects the abundances of these organisms in the community, glycosidic activities (for example, arabinan degradation) associated with *P. copri* colonization are a primary determinant of their metabolic responses.

## *B. infantis* Bg2D9, *P. copri* and intestinal gene expression

We next tested whether the effects of *B. infantis* Bg2D9 and *P. copri* on microbial community structure and expressed metabolic functions
<figure><img src="image_12.png" /><figcaption>**Fig. 2** | snRNA-seq analysis and targeted mass spectrometric analysis of intestinal tissue and plasma collected from mice containing bacterial communities with or without *P. copri* and two different strains of *B. infantis*.</figcaption></figure>

Jejunal tissue samples collected from arm 1 (with *P. copri* and with *B. infantis* Bg2D9) and arm 3 (without *P. copri* and with *B. infantis* Bg463) at the end of the experiment (P53) described in Fig. 1a were analysed (n = 4 samples/treatment arm for a d and f). a, The number of Recon2 reactions with statistically significant differences in their predicted flux between mice in Arm 1 and Arm 3. TA, transit amplifying. b, The number of Recon2 reactions in each Recon2 subsystem that are predicted to have statistically significant differences in their activities between the two treatment groups. Colours denote values normalized to the sum of all statistically significantly different Recon2 reactions found in all selected cell clusters for a given Recon2 subsystem in each treatment group. c, Proportional representation of cell clusters identified by snRNA-seq. Asterisks denote 'statistically credible differences' as defined by scCODA (Supplementary Table 10c and Methods). Mean values ± s.d. are shown. d, Selected Recon2 reactions in enterocyte clusters distributed along the villus involved in the urea cycle and glutamine metabolism. e, Targeted mass spectrometric quantifications

of citrulline levels along the length of the gut and in plasma. Mean values ± s.d. and *P* values from the two-sided Mann–Whitney *U* test are shown. Each dot represents an individual animal (n = 8 and 7 for arms 1 and 3, respectively). **f**, Effect of colonization with bacterial consortia containing or lacking *P. copri* on extracellular transporters for monosaccharides, amino acids and dipeptides. Ala, alanine; Arg, arginine; Asp, aspartate; Cys, cysteine; Gal, galactose; Glc, glucose; Gln, glutamine; Glu, glutamate; Gly, glycine; His, histidine; Ile, isoleucine; Leu, leucine; Lys, lysine; Met, methionine; Orn, ornithine; Phe, phenylalanine; Pro, proline; Sar, sarcosine; Ser, serine; Thr, threonine; Trp, tryptophan; Tyr, tyrosine; Val, valine. These transporters were selected, and the spatial information of their expressed region along the length of the villus was assigned based on published experimental evidence²⁹. Arrows in **d** and **f** indicate the 'forward' direction of each Recon2 reaction. The Wilcoxon rank-sum test was used to evaluate the statistical significance of the net reaction scores (**a**, **b**, **d** and **e**) between the two treatment groups. *P* values were calculated from Wilcoxon rank-sum tests and adjusted for multiple comparisons (Benjamini–Hochberg method); *q* < 0.05 was used as the cut-off for statistical significance.

were associated with metabolic changes in epithelial cells in portions of the small intestine that are dedicated to nutrient absorption. For this analysis, we used small intestinal samples from arms 1 and 3 for further analysis. This is because (1) the combination of *B. infantis* Bg2D9 and *P. copri* mediated a set of effects greater than those mediated by either organism alone and (2) successful colonization with both

*B. infantis* and *P. copri* before and through the weaning transition at P21 in arm 1 better represented the microbial communities of children in the clinical study¹¹. Because of the different *B. infantis* strains used and the greater fitness and greater expression of PULs targeting MDCF-2 glycans shown by *P. copri* but not *P. stercorea*, we refer to arm 1 as 'with *B. infantis* Bg2D9' and with *P. copri*' and arm 3 as 'with *B. infantis* Bg463'
and without *P. copri* for comparisons of the effect of the combination of these organisms.

There were no significant differences in jejunum villus height and crypt depth in mice 'with *B. infantis* Bg2D9 and with *P. copri* and 'with *B. infantis* Bg463 and without *P. copri* (n = 8 and 7 animals, respectively; Supplementary Table 9). Single-nucleus RNA sequencing (snRNA-seq) was used to investigate whether these two colonization states produced differences in expressed functions in jejunal tissue collected from P53 animals (n = 4 per treatment arm; Fig. 2a-c, Extended Data Fig. 5 and Supplementary Table 10). Cell clusters were assigned to the four principal intestinal epithelial cell lineages (enterocytic, goblet, enteroendocrine and Paneth cell) as well as to vascular endothelial cells, lymphatic endothelial cells, smooth muscle cells and enteric neurons (Extended Data Fig. 5a,b). Marker gene analysis allowed us to further subdivide the enterocytic lineage into three clusters: 'villus base', 'mid-villus' and 'villus tip'. The majority of all statistically significant differentially expressed genes (3,651 of 5,765; 63.3%) were assigned to these three enterocyte clusters (Extended Data Fig. 5c and Supplementary Table 10b).

We used NicheNet¹⁹ to identify potential ligand–receptor interactions between receiver and sender cells from our snRNA-seq dataset. We designated six of the epithelial cell clusters (crypt stem cells, proliferating transit amplifying/stem cells, villus base, mid-villus and villus tip enterocytes and goblet cells) as ‘receiver cells’, and all other clusters (both epithelial and mesenchymal) were designated ‘sender cells’. Extended Data Fig. 6 shows ‘bona fide ligand–receptor interactions’¹⁹ that are altered between the two colonization conditions for each receiver cell cluster. Ligands identified include those known to affect cell proliferation (igf-1), cell adhesion (cadm1, cadm3, cdh3, lama2, npnt), zonation of epithelial cell function/differentiation along the length of the villus (bmp4, bmp5) and immune responses (cadm1, il15, tgfbi, tnc) (Extended Data Fig. 6). Among all receiver cell clusters, crypt stem cells showed the highest number of altered bona fide ligand–receptor interactions. For example, igf-1 is known to enhance intestinal epithelial regeneration²⁰. We found that colonization with the *P. copri*-containing consortium was associated with markedly elevated expression of *igf-1* in goblet and lymphatic endothelial sender cells that signal to crypt stem cell receivers.

We subsequently used Compass²¹ and the Recon2²² database of metabolic reactions to generate in silico predictions of metabolic flux in different cell clusters: (1) stem cell and proliferating transit amplifying cell clusters positioned in crypts of Lieberkühn, (2) the three villus-associated enterocyte clusters and (3) the goblet cell cluster. Figure 2 shows the predicted metabolic flux differences (Methods) for enterocytes distributed along the length of the villus and in goblet cells. In clusters belonging to the enterocyte lineage, the number of statistically significant differences is greatest in villus base enterocytes and decreases towards the villus tip (Fig. 2a). Compared to mice in arm 3, those in arm 1 had greater predicted increases in the activities of subsystems related to energy metabolism, the metabolism of carbohydrates, amino acids and fatty acids, and various transporters, in villus base and mid-villus enterocytes (Fig. 2b and Extended Data Fig. 7).

While enterocytes prioritize glutamine as their primary energy source, they are also able to utilize fatty acids and glucose. We observed an increase in reactions related to fatty acid oxidation that occur in the villus enterocytes of mice in arm 1 compared to those in arm 3 extended to their crypts of Lieberkühn (Fig. 2b). Fatty acid oxidation has been linked to intestinal stem cell maintenance and regeneration²³. Mice colonized with *P. copri* and *B. infantis* Bg2D9 exhibited increases in the proportional representation of crypt stem cells and proliferating transit amplifying/stem cells but not in their villus-associated enterocytic clusters (Fig. 2c and Supplementary Table 10c). Compared to mice colonized with *B. infantis* Bg463 and lacking *P. copri*, those colonized with *B. infantis* Bg2D9 and *P. copri* also had predicted increases in energy metabolism in their goblet cells, as judged by the activities of

subsystems involved in glutamate metabolism, the urea cycle, fatty acid oxidation and glycolysis (Fig. 2b).

Citrulline is generally poorly represented in human diets; it is predominantly synthesized via metabolism of glutamine in small intestinal enterocytes and transported into the circulation²⁴ where it can serve as a quantitative biomarker of metabolically active enterocyte mass. Plasma levels are indicative of the absorptive capacity of the small intestine; they are lower in undernourished children and were increased in Bangladeshi children with moderate acute malnutrition after treatment with MDCF-2²⁵⁻²⁷. Compass predicted that mice harbouring *B. infantis* Bg2D9 and *P. copri* exhibit statistically significant increases in reactions involved in citrulline synthesis in villus base and mid-villus enterocyte clusters (*q* < 0.05 (adjusted *P* value); Wilcoxon ranked sum test; Fig. 2d). Follow-up targeted mass spectrometric analysis confirmed that citrulline was significantly increased in jejunal, ileal and colonic tissue segments, as well as in the plasma of mice in arm 1 compared to arm 3 (Fig. 2e and Supplementary Table 10d).

The presence of *P. copri* and *B. infantis* Bg2D9 was also associated with significantly greater predicted activities in the transport of nine amino acids (including the essential amino acids leucine, isoleucine, valine and phenylalanine), dipeptides and monosaccharides (glucose and galactose) in villus base and mid-villus enterocytes (Fig. 2f). These predictions suggest a greater absorptive capacity for these important growth-promoting nutrients, which are known to be transported within the jejunum at the base and middle regions of villi²⁴.

## *P. copri* effects on host metabolism and weight gain

We repeated the experiment described above with a larger number of animals (4 dually housed germ-free dams yielding 18–19 viable pups per arm). To examine whether the weight gain phenotype and metabolic alterations observed in the experiment described above could be attributed to the presence or absence of *P. copri* in the microbial community, we administered *B. infantis* Bg2D9 in both arms of this repeat experiment. Outside of this change, the same cultured strains, the same sequence of their introduction and the same sequence of diet switches were applied (Extended Data Fig. 8a). Reproducible colonization of consortium members within each arm was confirmed by quantifying their absolute abundances in caecal samples collected at the time of euthanasia (P53; Supplementary Table 11). As in the previous experiment, animals in the arm containing *P. copri* exhibited significantly greater weight gain between P23 and P53 than those in the no *P. copri* arm (Extended Data Fig. 8b and Supplementary Table 12).

We used targeted mass spectrometry to quantify levels of 20 amino acids, 19 biogenic amines and 66 acylcarnitines in the jejunum, colon, gastrocnemius, quadriceps, heart muscle and liver of the two groups of mice. In addition, we quantified the 66 acylcarnitines in their plasma. The results are described in Supplementary Table 13 and Extended Data Fig. 8c–e. Consistent with the previous experiment, citrulline, a biomarker for metabolically active enterocyte biomass, was significantly elevated in the jejunums of mice belonging to the with-P. copri group (Extended Data Fig. 8c and Supplementary Table 13a). We observed significant elevations of acylcarnitines derived from palmitic acid (C16:0), stearic acid (C18:0), oleic acid (C18:1), linoleic acid (C18:2) and linolenic acid (C18:3) in the jejunums of *P. copri*-colonized animals (Extended Data Fig. 8d); these are the major fatty acids found in soybean oil triglycerides²⁸, which is the principal source of lipids in MDCF-2. These acylcarnitine chain lengths were found at higher abundances than all other medium or long-chain acylcarnitine species in our jejunal samples, indicating their role as primary dietary lipid energy sources (Supplementary Table 13b). Elevation of these species suggests increased transport and β-oxidation of long-chain dietary lipids in the jejunums of the *P. copri*-colonized animals.

Analysis of colonic tissue showed significant elevation of C16:0, C18:1 and C18:2 acylcarnitines in *P. copri*-colonized animals, suggesting that β-oxidation is also elevated in tissue compartments not directly
involved in lipid absorption (*P* < 0.01; Mann–Whitney *U* test) (Extended Data Fig. 8e and Supplementary Table 13b). This finding was matched by a significant elevation in plasma levels of non-esterified fatty acids in *P. copri*-colonized animals, which would support fatty acid β-oxidation in peripheral tissues (Extended Data Fig. 8f and Supplementary Table 13c). In addition, colonic (and jejunal) levels of C3 and C4 acyl-carnitines known to be derived from branched-chain amino acid catabolism were significantly elevated in the *P. copri*-colonized animals (Extended Data Fig. 8d,e and Supplementary Table 13b).

## P. copri colonization and weight gain with another diet

We next looked at the effects of *P. copri* colonization in the context of a 'control' diet representative of that typically consumed by 18-month-old Bangladeshi children living in Mirpur ('Mirpur-18 diet')⁴. The design was similar to that used for the experiments described in Fig. 1 and Extended Data Fig. 8 with two exceptions: (1) *B. infantis* Bd2D9 was used in both groups and (2) on P24, pups from different litters were mixed and randomly assigned to two diet treatment groups, MDCF-2 and Mirpur-18 (*n* = 2 dams and 12 pups per group). The absolute abundances of community members were quantified in caecal contents collected at the time of euthanasia on P53. While the absolute abundance of *P. copri* Bd131 was not significantly different between the two diet groups, there were statistically significant differences between the abundances of 11 of the 19 community members (Extended Data Fig. 9). Nonetheless, we proceeded to test whether the increased weight gain phenotype associated with the presence of *P. copri* in the community was evident in the Mirpur-18 diet context. To do so, we repeated the dam-to-pup microbial transmission experiment where all animals were weaned onto the Mirpur-18 diet but where one group had received *P. copri* Bd131 (*n* = 8 animals) and the other had not (*n* = 9). All animals were euthanized on P53. *P. copri* successfully colonized mice and was maintained throughout the experiment at levels comparable to previous experiments (10.4 ± 0.1 log₁₀ genome equivalents per gram of caecal contents at P53). Importantly, there was no statistically significant difference in weight gain between the two groups (*P* = 0.297; linear mixed-effects model (Methods)). These findings provide evidence that the effect of *P. copri* on weight gain in this preclinical model is diet dependent.

## Tests of *P. copri* isolates resembling MAGs Bd0018 and Bd0019

We previously characterized five additional faecal *P. copri* strains that we cultured from Bangladeshi children living in Mirpur¹¹. Two of these strains (BgdS_2 and BgF_2) had greater genomic similarity to MAGs Bd0018 and Bd0019 than the other isolates, including *P. copri* Bd131, as quantified by phylogenetic distance, PUL content and the representation of metabolic pathways (Fig. 3a and Supplementary Table 14a–d); for example, nine of the ten functionally conserved PULs in Bd0018/Bg0019 were present in *P. copri* BdS_2 and BgF_2 as were 53 of 55 carbohydrate utilization pathways (Fig. 3a and Supplementary Table 14a–d). In addition, results of in vitro growth assays conducted in defined medium supplemented with different glycans represented in MDCF-2 disclosed that strain BgF5_2 showed stronger preference than Bd131 for glycans enriched in and/or unique to MDCF-2 compared to ready-to-use supplementary food (that is, arabinan, arabinoxylan, galactan and galactomannan)¹¹.

To directly determine whether pre-weaning colonization with *P. copri* strains resembling MAGs Bd0018 and Bd0019 is sufficient to promote growth and produce the metabolic effects described above, we performed an experiment whose design (Fig. 3b) was similar to our previous experiments (Fig. 1b and Extended Data Fig. 8a) but with several modifications. First, because of their greater genomic similarity to WLZ-associated MAGs Bd0018 and Bd0019, we used *P. copri* BdS_2 and BgF5_2 in place of *P. copri* Bd131. Second, to control for differences in *B. infantis* strain used in the initial experiment (Fig. 1), both arms received *B. infantis* Bd2D9 in this 'third' experiment (as was the case

in the second experiment described in Extended Data Fig. 8). Third, because *P. stercorea* had colonized at a lower abundance than *P. copri* and did not express CAZymes related to MDCF-2 glycans, it was not included in the second gavage mixture, which now only contained *P. copri*. Fourth, given that *B. infantis* Bd2D9 promoted pre-weaning colonization of *P. copri* in the initial experiment, we omitted the fourth gavage, previously administered at the end of the weaning period, that had included *P. copri* and *P. stercorea*. As before, the control group of animals were those that did not receive *P. copri* (*n* = 2 dams and 13 pups per treatment group).

Shotgun sequencing of DNA isolated from caecal contents collected at the time of euthanasia (P53) confirmed that animals in the experimental group had been colonized with both *P. copri* isolates as well as all other members of the defined consortia (Supplementary Table 15). Even though strain BdS_2 grew much more poorly than BgF5_2 when cultured in defined medium¹¹, in animals colonized with both isolates, the BdS_2 strain was present at higher absolute abundance than the BgF5_2 strain (Fig. 3c) (their relative abundances were 37.8 ± 4.4% and 15.5 ± 1.0%, respectively, versus 31 ± 6.6% and 24 ± 8.0% for *P. copri* Bd131 in the experiments described in Fig. 1b and Extended Data Fig. 8a). Comparing the two groups disclosed that colonization with BdS_2 and BgF5_2 augmented community biomass without displacing other bacteria (Fig. 3c and Supplementary Table 15b).

We observed a significantly greater increase in body weight between P23 and P53 in mice colonized with *P. copri* BdS_2 and BgF5_2 compared to those without *P. copri* (Fig. 3d and Supplementary Table 16). The difference in the mean percentage increase in postweaning weight between the experimental and control groups (24%) was comparable to that documented in the two previous experiments shown in Fig. 1b and Extended Data Fig. 8a (25% and 13%, respectively). As in these previous experiments, the weight difference was not attributable to differences in caecal size.

Mass spectrometry confirmed that preweaning colonization with *P. copri* affected intestinal lipid metabolism and was a major determinant of MDCF-2 glycan degradation. Targeted LC-MS of ileal and colonic tissue revealed a significant elevation of long-chain acylcarnitines corresponding to soybean oil lipids (Extended Data Fig. 10 and Supplementary Table 17), consistent with changes observed in the experiment described in Extended Data Fig. 8. Microbial RNA-seq of caecal contents revealed that among all PUL genes, those present in the three conserved PULs with predicted arabinan/starch, pectin and pectic galactan substrates were significantly enriched for higher levels of expression (Methods; Fig. 3e and Supplementary Table 18). UHPLC-QqQ-MS of monosaccharides in glycans present in caecal contents indicated that the presence of *P. copri* BdS_2 and BgF5_2 resulted in significantly lower levels of arabinose, consistent with our previous observations using *P. copri* Bd131, as well as galactose (a finding specific to this experiment) (Fig. 3f and Supplementary Table 19). Colonization with *P. copri* BdS_2 and BgF5_2 also significantly lowered levels of all arabinose-containing glycosidic linkages measured, as well as three galactose-containing linkages (Fig. 3g,h and Supplementary Table 19). Together, these data indicate that the PUL content of these two isolates is associated with enhanced degradation of MDCF-2 glycans compared to the *P. copri* Bd131-containing microbial community. Targeted UPHLC-QqQ-MS measurements of all 20 amino acids and 7 B vitamins also revealed that compared to the control group, colonization with *P. copri* BdS_2 and BgF5_2 was associated with significantly higher caecal levels of two essential amino acids (tryptophan, lysine), seven non-essential amino acids (glutamate, glutamine, aspartate, asparagine, arginine, proline, glycine) and pantothenic acid (vitamin B5) (Supplementary Table 20).

Based on all of these experiments, we concluded that (1) pre-weaning colonization with *P. copri* augments weight gain in the context of the MDCF-2 diet, (2) the presence of specific strains of this species is a major determinant/effector of MDCF-2 glycan degradation
<figure><img src="image_13.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>g</figcaption></figure>

<figure><img src="image_19.png" /><figcaption>h</figcaption></figure>

<figure><img src="image_20.png" /></figure>

**Fig. 3** | Testing the effects of pre-weaning colonization with two *P. copri* strains closely related to MAGs Bg0018 and Bg0019 on host weight gain and MDCF-2 glycan degradation. **a,** Comparison of PULs highly conserved in the two *P. copri* MAGs with their representation in the three cultured *P. copri* strains. **b,** Study design (n = 2 dams and 13 offspring per treatment arm). **c,** Absolute abundance of *P. copri* strains and total bacterial load in caecal contents collected at the end of the experiment (P53). Exact *P* values for comparisons of BgD5_2 and BgF5_2 and total bacterial load are 2 × 10⁻⁵ and 2 × 10⁻⁵, respectively. **d,** Body weights of the offspring of dams, normalized to P23. The *P* value for the group difference is *P* = 4 × 10⁻⁵ (linear mixed effects model (Methods)). **e,** GSEA of expression of PULs shared by *P. copri* BgD5_2 and BgF5_2 in the caecal contents of animals. Benjamini–Hochberg adjusted *P* values were calculated using GSEA ranking genes by their mean log₂ TPM across the *P. copri* colonized samples, with each PUL comprising a gene set against the background of all predicted PUL genes. Violin plots show the log₂ TPM of all genes assigned to any of the 22

predicted PULs in each isolate (n = 201 genes) in each of the samples, split to show homologues of consensus PUL 17 (arabinan, starch; n = 22 genes), PUL 4 (pectin; n = 13 genes) and PUL 16 (pectic galactan; n = 15 genes) in colour compared to the remainder of all PUL genes in grey. Internal box plots show the median (circle) and quartiles (box boundaries) for all genes assigned to PULs. P = 1 × 10⁻⁴ for PUL 17. **f,** UHPLC-QqQ-MS analysis of total arabinose and galactose in glycans present in caecal contents collected at P53. The *P* value for both arabinose and galactose is 2 × 10⁻⁵. **g,h,** UHPLC-QqQ-MS of glycosidic linkages containing arabinose (g) and galactose (h) in caecal contents. The exact *P* values for t-Araf, 2-Araf, 2,3-Araf, 3,4-Xylp/3,5-Araf and 5-Araf (g) are 2 × 10⁻⁵, 8 × 10⁻⁵, 2 × 10⁻⁵, 2 × 10⁻⁵ and 2 × 10⁻⁵, respectively. The exact *P* values for 2,4,6-galactose, 3,4,6-galactose and 4-galactose are 3 × 10⁻⁵, 2 × 10⁻⁵ and 2 × 10⁻⁵, respectively. Mean values ± s.d. are shown. *P* values were calculated using a two-sided Mann–Whitney *U* test (c,f–h). Each dot in **b–h** represents an individual animal.
and (3) incorporating these strains into the gut community changes intestinal cellular metabolism.

Discussion

Accessing tissue from different regions of the human intestine and extra-intestinal sites represents a major challenge when trying to characterize the mechanisms by which microbiome-targeted nutritional interventions impact the microbiota and human physiology at a molecular, cellular and systems level. In this study, we illustrate a 'reverse translation' strategy that can be used to address this challenge. Gnotobiotic mice were colonized with defined consortia of age- and WLZ-associated bacterial strains cultured from faecal samples collected from children living in a Bangladeshi community where the prevalence of malnutrition is high. *P. copri* was represented by cultured isolates whose genomic features, including PULs and metabolic pathways involved in carbohydrate utilization, are highly similar to MAGs associated with improved weight gain in the clinical trial. Dam-to-pup transmission of these communities occurred in the context of a sequence of diets that re-enacted those consumed by children enrolled in a clinical study of a MDCF-2. Microbial RNA-seq and targeted mass spectrometry of glycosidic linkages present in intestinal contents provided evidence that *P. copri* plays a key role in the metabolism of polysaccharides contained in MDCF-2. Consistent with this, *P. copri*-associated weight gain in the preclinical model was dependent on consumption of MDCF-2; this phenotype was not observed when a diet commonly consumed by Bangladeshi children was administered, despite comparable levels of *P. copri* colonization in the two diet contexts. snRNA-seq and targeted mass spectrometry of the intestine indicated that colonization with the consortium that contains the combination of *B. infantis* Bg2D9 and *P. copri* increases the uptake and metabolism of lipids (including those fatty acids that are most prominently represented in the soybean oil that comprises the principal lipid component of MDCF-2). Additional effects on uptake and metabolism of amino acids (including essential amino acids) and monosaccharides were predicted and in select cases validated by mass spectrometric assays. These effects on nutrient processing and energy metabolism involve proliferating epithelial progenitors in the crypts as well as their descendant lineages distributed along the villus. snRNA-seq revealed discrete spatial features of these effects, with populations of enterocytes positioned at the base-, mid- and tip regions of villi manifesting distinct patterns of differential expression of a number of metabolic functions.

Inspired by the results of the clinical trial, the goal of our 'reverse translation' experiments was to ascertain the impact of the presence or absence of *P. copri* in a model that emulated postnatal gut microbial community assembly and exposure to MDCF-2. The current study raises several questions that have both mechanistic and therapeutic implications. We were not able to successfully mono-colonize mice with our cultured Bangladeshi *P. copri* strains. Therefore, we could not directly test the effects of these strains in vivo on MDCF-2 glycan metabolism, weight gain and/or gut epithelial biology in the absence of other potential microbial interactions. Moreover, findings from the current study together with our findings from direct analysis of the faecal microbiomes of participants in the clinical trial¹¹ indicate that degradation of MDCF-2 glycans is necessary for promoting weight gain, albeit involving the actions of downstream metabolite(s)/signalling events that remain to be fully characterized. Additional work, involving systematic manipulation of the composition of the bacterial consortia introduced into dams (and subsequently transmitted to their offspring) will be required to ascertain the extent to which *P. copri* has (1) direct effects on intestinal epithelial gene expression, host metabolism and weight gain versus (2) effects of other community members that are dependent upon its presence or absence. If *P. copri* has direct effects on the host, it remains to be determined whether the mediators of these effects are the direct products of MDCF-2 glycan metabolism or the products of other metabolic pathways in *P. copri* whose activities

are regulated by biotransformation of these glycans, or other MDCF-2 components. Future studies are also needed to disambiguate the extent to which the observed effects of *P. copri* on gut epithelial carbohydrate, lipid and amino acid metabolism contribute to weight gain. The spatial features of metabolic pathway expression documented by snRNA-seq must be characterized further. This effort will be technically challenging; for example, it could require (1) documenting the distribution of *P. copri* and other community members along the crypt-villus axis, (2) advancing methods for spatial transcriptomics²⁹,³⁰ so that they can be (simultaneously) applied to both epithelial cell lineages and microbial community members and (3) using in situ mass spectrometry to directly characterize the metabolic profiles of discrete gut cell populations.

The relationship between prominent initial colonization by *B. infantis* Bg2D9 and the capacity of *P. copri* to subsequently colonize also needs further investigation. *B. infantis* Bg2D9 contains several genomic loci not represented in most other cultured *B. infantis* strains, which could enhance its ability to utilize a variety of dietary carbohydrates¹³. In principle, these loci could support increased fitness of *B. infantis* Bg2D9 in malnourished children whose consumption of breast milk is low. Given that Bangladeshi infants and young children with SAM have markedly lower levels or even completely lack *B. infantis* compared to their healthy counterparts¹³, the *B. infantis*–*P. copri* interaction documented in this preclinical study provides a rationale for testing the effects of first administering *B. infantis* Bg2D9 to individuals with SAM and subsequently MDCF-2 to restore age-appropriate microbiome configurations and promote healthy growth.

In summary, this and our companion study¹¹ illustrate one approach for identifying members of a gut microbial community that function as principal metabolizers of dietary components as well as key effectors of host biological responses. The results can provide a starting point for developing microbiome-based diagnostics for stratification of populations of undernourished children who are candidates for treatment with MDCF and for monitoring their treatment responses including in adaptive clinical trial designs. Another potential return on investment for this approach is a knowledge base needed for (1) creating ‘next generation’ MDCFs composed of (already identified) bioactive glycans but from alternative food sources that may be more readily available, affordable and culturally acceptable for populations living in different geographic locales; (2) making more informed decisions about dosing of an MDCF for undernourished children as a function of their stage of development (age) and disease severity and (3) evolving policies about complementary feeding practices based on insights about how food components impact the fitness and expressed beneficial functions of growth-promoting elements of a child’s microbiome.

Methods

Ethics statement

The studies reported complied with all applicable ethical regulations. Bacterial strains were cultured from faecal samples collected with informed consent, under protocols approved by the International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b) Ethical Review Committee. Material transfer agreements between icddr,b and Washington University in St. Louis were established for the use of these samples. Gnotobiotic mouse experiments were performed following Institutional Animal Care and Use Committee and Institutional Biological and Chemical Safety Committee protocols approved by the Washington University Animal Studies and Environmental Health and Safety Committee.

Bacterial genome sequencing and annotation

Monocultures of each isolate were grown overnight at 37 °C in Wilkins-Chalgren Anaerobe Broth (Oxoid, catalogue number CM0643) in a Coy Chamber under anaerobic conditions (atmosphere; 75% N₂, 20% CO₂ and 5% H₂) without shaking. Cells were recovered by centrifugation
(5,000 × g for 10 min at 4 °C). High-molecular-weight genomic DNA was purified (MagAttract HMW DNA kit, Qiagen) following the manufacturer's protocol, and the amount was quantified (Qubit fluorometer). The sample was passed up and down through a 29-gauge needle 6–8 times, and the fragment size distribution was determined (~30 kbp; TapeStation, Agilent).

Fragmented genomic DNA (400–1,000 ng) was prepared for long-read sequencing using a SMRTbell Express Template Prep Kit 2.0 (Pacific Biosciences) adapted to a deep 96-well plate (Fisher Scientific) format. All DNA handling and transfer steps were performed with wide-bore, genomic DNA pipette tips (ART). Barcoded adapters were ligated to A-tailed fragments (overnight incubation at 20 °C), and damaged or partial SMRTbell templates were subsequently removed (SMRTbell Enzyme Cleanup Kit). High-molecular-weight templates were purified (volume of added undiluted AMPure beads = 0.45 times the volume of the DNA solution). Libraries prepared from different strains were pooled (3–6 libraries per pool). A second round of size selection was then performed; AMPure beads were diluted to a final concentration of 40% (v/v) with SMRTbell elution buffer with the resulting mixture added at 2.2 times the volume of the pooled libraries. DNA was eluted from the AMPure beads with 12 µl of SMRTbell elution buffer. Pooled libraries were quantified (Qubit), their size distribution was assessed (TapeStation) and sequenced (Sequel System, Sequel Binding Kit 3.0 and Sequencing Primer v4 (Pacific Biosystems)). The resulting reads were demultiplexed, and Q20 circular consensus sequencing reads were generated (Cromwell workflow configured in SMRT Link software). Genomes were assembled using Flye³¹ (v2.8.1) with HiFi-error set to 0.003, min-overlap set at 2,000 and other options set to default. Genome quality was evaluated using CheckM³² (v1.1.3) (Supplementary Table 1a, 5 and 14a).

We applied Prokka³³ (v1.14) to identify potential open reading frames (ORFs) in each assembled genome. Additional functional annotation of these ORFs using a 'subsystems' approach adapted from the SEED genome annotation platform³⁴ was performed as described in our companion study¹¹. We assigned functions to 9,820 ORFs in 20 isolate genomes using a collection of mcSEED metabolic subsystems that capture the core metabolism of 98 nutrients/metabolites in four major categories (amino acids, vitamins, carbohydrates and fermentation products) projected over 2,856 annotated human gut bacterial genomes³⁵⁻³⁷. In silico reconstructions of selected mcSEED metabolic pathways were based on functional gene annotation and prediction using homology-based methods and genome context analysis. Reconstructions were represented as a binary phenotype matrix (BPM) where assignment of a '1' for amino acids and B vitamins denotes a predicted prototroph and a '0' denotes an auxotroph. For carbohydrates, '1' and '0' refer to a strain's predicted ability or inability, respectively, to utilize the indicated mono-, di- or oligosaccharide, while for fermentation end products, a '1' and '0' indicate a strain's predicted ability or inability to produce the indicated compound, respectively (Supplementary Tables 1c and 14d).

To calculate phylogenetic relationships between five *P. copri* isolates and MAGs Bg0018 and Bg0019, we first used CheckM³² (v1.1.3) to extract and align the amino acid sequences of 43 single-copy marker genes in each isolate or each of the two MAGs, plus an isolate genome sequence of *Bacteroides thetaiotaomicron* VPI-5482 (PATRIC accession number 226186.12). Concatenated marker gene sequences were analysed using fasttree³⁸ (v2.1.10) to construct a phylogenetic tree using the Jones–Taylor–Thornton model and CAT evolution rate approximation, followed by tree rescaling using the Gamma20 optimization. The tree was subsequently processed in R using ape³⁹ (v5.6-2) to root the tree with the *B. thetaiotaomicron* genome and extract phylogenetic distances between genomes, followed by ggtree⁴⁰ (v3.2.1) for tree plotting.

The similarity between the genomes of these strains and MAGs was quantified by calculating the ANI score with pyani⁴¹ (ANIm (ANI calculated with the MUMmer algorithm) implementation of ANI, v0.2.10).

We first calculated ANIm scores for all possible combinations between MAGs and the genomes of cultured bacterial strains and subsequently removed any MAG-strain genome combination with <10% alignment coverage¹⁷. For the remaining MAGs, a 'highly similar' genome in the collection of cultured bacterial strains was defined as having >94% ANIm score (Supplementary Table 1b). We then determined the degree of binary phenotype concordance between each genome in the collection of cultured bacterial strains and its 'highly similar' MAG. A binary phenotype concordance score was calculated by dividing the number of binary phenotypes¹¹ shared between a cultured strain's genome and a MAG by the total number of binary phenotypes annotated in the strain and MAG. A 'representative MAG' for each genome was defined as having a binary phenotype concordance score >90% (Supplementary Table 1c).

PULs were predicted based on the method described in ref. 42 and displayed with the Polysaccharide Utilization Loci Database (PULDB) interface⁴³. PULs were placed into three categories: (1) 'functionally conserved' (PULs containing shared ORFs encoding the same CAZymes and SusC/SusD proteins in the same organization in their respective genomes with ≥90% amino acid identity between proteins); (2) 'structurally distinct' (PULs present in respective genomes but where one or more CAZymes or one or both SusC/SusD proteins are missing or fragmented in a way likely to impact function, or where extra PUL elements are present) and (3) 'not conserved' (PULs present in respective genomes but with mutations likely to completely compromise function, or no PUL identified).

## Colonization and husbandry

Germ-free C57BL/6J mice were maintained in plastic flexible film isolators (Class Biologically Clean) at 23 °C under a strict 12 h light cycle (lights on at 0600 h). Autoclaved paper 'shepherd shacks' were kept in each cage to facilitate natural nesting behaviours and provide environmental enrichment.

A weaning diet containing MDCF-2 was formulated as described in the main text. Ingredients represented in the different diet modules were combined, and the mixture was dried, pelleted and sterilized by gamma irradiation (30–50 kGy). Sterility was confirmed by culturing the pellets in brain-heart infusion medium supplemented with 0.5% yeast extract (LYBHI medium⁴⁴) and in Wilkins-Chalgren Anaerobe Broth under aerobic and anaerobic conditions for 7 days at 37 °C followed by plating on LYBHI- and blood-agar plates. Nutritional analysis of each irradiated diet was performed by Nestlé Purina Analytical Laboratories (Supplementary Table 2c).

Pregnant C57BL/6J mice originating from trio matings were given ad libitum access to an autoclaved breeder chow (Purina Mills; Lab Diet 5021) throughout their pregnancy and until postpartum day 2. Key points about the experimental design of the gnotobiotic mouse experiments described in Figs. 1b and 3b and Extended Data Figs. 8a and 9a are as follows: (1) all bacterial strains were cultured in Wilkins-Chalgren Anaerobe Broth (except for *F. prausnitzii* which was cultured in LYBHI medium) and were collected after overnight growth at 37 °C (Supplementary Table 1a), (2) all gavage mixtures contained equivalent amounts (by OD₆₀₀) of their constituent bacterial strains except for *F. prausnitzii* which was concentrated 100-fold before preparing the gavage mixture, (3) each bacterial consortium was administered to the postpartum dams in a volume of 200 µl using an oral gavage needle (Cadence Science; catalogue number 7901), (4) the number of dams and pups per treatment group (2 dams and 7–8 pups per treatment group for the experiment described in Fig. 1b; 4 dams and 18–19 pups per treatment group for the experiment outlined in Extended Data Fig. 8a; 2 dams and 13 pups per treatment group for the experiment illustrated in Fig. 3b; 2 dams and 12 pups per treatment group for the experiment shown in Extended Data Fig. 9a), (5) half of the bedding was replaced with fresh bedding in each cage each day from postpartum day 1 to day 14, after which time bedding was changed every 7 days, (6) diets
were provided to mothers as well as to their weaning and post-weaning pups ad libitum and (7) biospecimens collected from mice when they were euthanized (without previous fasting) were snap frozen in liquid nitrogen and stored at -80 °C before use.

Pups were weighed on P23, P35 and P53 and normalized to the weight on P23. A linear mixed-effects model was used to evaluate the effect of different microbial communities on normalized mouse weight gain:

$$ \text{Normalized weight} \sim \beta_1 (\text{arm}) + \beta_2 (\text{postnatal day}) + (1|\text{mouse}) $$

## Defining the absolute abundances of bacterial strains in ileal, caecal and faecal communities

The absolute abundances of bacterial strains were determined using previously described methods with minor modifications$^{45,46}$. In brief, 3.3 × 10⁶ cells of *Alicyclobacillus acidiphilus* DSM 14558 and 1.49 × 10⁷ cells of *Agrobacterium radiobacter* DSM 30147 (ref. 45) were added to each weighed frozen sample before DNA isolation and preparation of barcoded libraries for shotgun sequencing. Sequencing was performed using an Illumina NextSeq instrument. Bacterial abundances were determined by assigning reads to each bacterial genome, followed by a normalization for genome uniqueness in the context of a given community$^{46}$. The resulting count table was imported into R$^{47}$ (v4.0.4). We calculated the absolute abundance of a given strain *i* in sample *j* in reference to the spike-in *A. acidiphilus* (*Aa*) and *A. radiobacter* (*Ar*) genomes using the following equation:

$$ \text{strain}_{i,j} = \left( \frac{\text{counts}_{i,j} \times \text{Aa cells added}_j}{\text{Aa counts}_{i,j} \times \text{sample weight}_j} + \frac{\text{counts}_{i,j} \times \text{Ar cells added}_j}{\text{Ar counts}_{i,j} \times \text{sample weight}_j} \right) \times 0.5 $$

The statistical significance of observed differences in the abundance of a given strain between experimental groups in ileal, caecal and faecal samples was determined by using the Kruskall-Wallis test followed by Dunn's test for each pairwise comparison among the three arms in the experiment described in Fig. 1 or the Mann-Whitney U test between the two arms in experiments described in Fig. 3 and Extended Data Figs. 1b and 9. The statistical significance of differences in the composition of communities in the three arms described in Fig. 1g was determined using permutational multivariate analysis of variance (PERMANOVA)$^{48}$ on sample projections onto principal components calculated from the log₁₀ absolute abundance profiles of the 16 organisms that were not *B. infantis* or *Prevotella* species. The statistical significance of observed differences in the abundance of a given strain across different treatment groups and time was tested using a linear mixed effects model within the R packages lme4$^{49}$ (v1.1-27) and lmerTest$^{50}$ (v3.1-3). The change in *P. copri* absolute abundance in faecal samples during the course of the experiment was determined by a linear mixed-effects model:

$$ \text{P. copri} \text{ absolute abundance} \sim \beta_1 (\text{arm}) + \beta_2 (\text{postnatal day}) + (1|\text{mouse}) $$

For the experiment described in Fig. 1, 96 faecal samples were sequenced (2.2 × 10⁶ ± 1.2 × 10⁵ unidirectional 75 nt reads per sample (mean ± s.d.)), along with 20 caecal samples (1.5 × 10⁶ ± 4.2 × 10⁵ unidirectional 75 nt reads per sample) and 20 ileal samples (1.5 × 10⁶ ± 9.1 × 10⁴ unidirectional 75 nt reads per sample) (Supplementary Table 4a). For the experiment described in Extended Data Fig. 1b, 37 faecal samples were sequenced (5.8 × 10⁶ ± 1.6 × 10⁶ unidirectional 75 nt reads per sample) (Supplementary Table 6a)), while for the experiment described in Extended Data Fig. 8a, 37 caecal samples were sequenced (1.3 × 10⁶ ± 1.3 × 10⁵ unidirectional 75 nt reads per sample (Supplementary Table 11a)). For the experiment described in Fig. 3b, 26 caecal samples were sequenced (2.9 × 10⁶ ± 5.6 × 10⁶ 150 nt paired-end reads per sample (Supplementary Table 15a)).

## Microbial RNA-seq

RNA was isolated⁶ from caecal contents collected at the end of the experiment. Complementary DNA libraries were generated from isolated RNA samples using the 'Total RNA Prep with Ribo-Zero Plus' kit (Illumina). Barcoded libraries were sequenced (Illumina NovaSeq instrument). For the experiment described in Fig. 1b, cDNA recovered from 20 different samples of caecal contents were each sequenced to a depth of 7.8 × 10⁷ ± 9.6 × 10⁶ bidirectional 150 nt reads (mean ± s.d.) (Supplementary Table 8a). For the experiment shown in Fig. 3b, DNA from 26 different samples of caecal contents were each sequenced to a depth of 6.5 × 10⁷ ± 2.1 × 10⁷ bidirectional 150 nt reads (Supplementary Table 18a). Raw reads were trimmed by using TrimGalore⁵¹ (v0.6.4). Trimmed reads longer than 100 bp were mapped to reference genomes with kallisto⁵² (v0.43.0).

To analyse expression of genes from the *P. copri* and *P. stercorea* PULs described in Fig. 3e and Extended Data Fig. 4a,b, transcripts per million (TPM) values were obtained by mapping reads, using kallisto, to their genomes (Figs. 1b and 3b). TPM normalized expression was used to control for differences in library depth and gene length. For Fig. 3e and Extended Data Fig. 4a,b, log₂ TPM with a pseudocount of 1 were visualized for all predicted PUL genes using the seaborn⁵³ (v0.12.1) Python package, splitting the set to show individual PULs in colour on the right side of each violin plot against the remainder of the PUL genes shown in grey on the left side of each violin plot. Benjamini-Hochberg adjusted GSEA P values were calculated with fgsea⁵⁴ (v1.20.0), ranking genes by their mean log₂ TPM across the *Prevotella*-colonized samples in a given arm. Each PUL comprised a gene set against the background of all PUL genes from a given isolate, with minimum and maximum gene set sizes of 5 and 50 genes, respectively (Supplementary Table 8b and Supplementary Table 18b).

For the analysis described in Fig. 3e, modifications were made to account for the high fraction of genes with 100% nucleotide identity between *P. copri* BgD5_2 and BgF5_2. The kallisto index was generated to represent the set of unique genes between the two isolates by including the entire *P. copri* BgF5_2 gene set and only the subset of *P. copri* BgD5_2 genes that did not share identical sequences in *P. copri* BgF5_2 (61 out of 3,066 genes). TPMs were filtered down to the set of unique genes in the *P. copri* BgD5_2 and BgF5_2 PULs (n = 201 unique genes, consisting of 195 genes with identical nucleotide sequences between the two isolates and 3 unique genes from each isolate (Supplementary Table 14c)).

For abundance-normalized differential expression analysis of microbial transcripts described in Extended Data Fig. 3d, paired metagenomic and meta-transcriptomic kallisto pseudocounts were generated by mapping reads from the caecal DNA and cDNA libraries described above to the specific set of bacteria administered in the arm from which the sample was derived. To prevent skewing of library size normalization, the resulting counts matrix then underwent filtering of counts for rRNA loci predicted by Prokka³³. MTXmodel³⁵ is a generalized linear model-based approach for testing for differential expression of genes in microbial communities that controls for false positives due to differences in underlying metagenomic abundances. MTXmodel was run using the paired metagenomic and meta-transcriptomic counts tables for each organism to achieve 'within-taxon-sum-scaling'⁵⁶. Each pairwise comparison of arms was run using the respective sets of samples with 'arm' as the single fixed effect in the generalized linear model design. Transcripts with statistically significant differences in their expression (that is, 'arm' coefficients that were significantly different from 0) after normalizing for absolute abundance were identified (q value (Benjamini-Hochberg adjusted P value) < 0.1) after multiple hypothesis correction was applied to the entire set of transcripts from a given organism. GSEA was performed for each organism with fgsea⁵⁴, ranking all genes surviving zero-filtering from MTXmodel differential expression testing by their estimated log₂ fold difference. Each mcSEED metabolic pathway in each organism was used as a gene set against the background of all
genes tested for differential expression, with minimum and maximum gene set sizes of 5 and 50 genes, respectively.

## Histomorphometric analysis of villus height and crypt depth

Jejunal and ileal segments were fixed in formalin and embedded vertically in paraffin; 5 µm-thick sections were prepared, and the sections were stained with haematoxylin and eosin. Slides were scanned (Nano-Zoomer instrument, Hamamatsu). For each animal, ten well-oriented crypt-villus units were selected from each intestinal segment for measurement of villus height and crypt depth using QuPath³⁷ (v0.3.2). Measurements were performed with the investigator blinded with respect to colonization group. A two-tailed Mann-Whitney U test was applied to the resulting datasets.

## snRNA-seq

Jejunal segments (1.5 cm in length) were collected from mice and snap frozen in liquid nitrogen (n = 4 animals per treatment group (2 male mice and 2 female mice); 2 treatment groups in total). The method for extracting nuclei was adapted from a previously described protocol for the pancreas⁵⁸. Briefly, tissues were thawed and minced in lysis buffer (25 mM citric acid, 0.25 M sucrose, 0.1% NP-40, and 1× protease inhibitor (Roche)). Nuclei were released from cells using a pestle duncer (Wheaton), washed 3 times with buffer (25 mM citric acid, 0.25 M sucrose and 1× protease inhibitor) and filtered successively through 100 µm, 70 µm, 40 µm, 20 µm and finally 5 µm diameter strainers (pluriSelect) to obtain single nuclei in resuspension buffer (25 mM KCl, 3 mM MgCl₂, 50 mM Tris, 1 mM DTT, 0.4 U µl⁻¹ rNase inhibitor (Sigma) and 0.4 U µl⁻¹ Superase inhibitor (ThermoFisher)). Approximately 10,000 nuclei per sample were subjected to gel bead-in-emulsion generation, reverse transcription and construction of libraries for sequencing according to the protocol provided in the 3' gene expression v3.1 kit manual (10× Genomics). Libraries were balanced, pooled and sequenced (Illumina NovaSeq S4; 3.23 × 10⁸ ± 1.39 × 10⁷ paired-end 150 nt reads per nucleus (mean ± s.d.) from jejunal samples). Read alignment, feature-barcode matrices and quality controls were processed by using the 10× Genomics CellRanger⁵⁹ 5.0 pipeline with the flag '-include-introns' to ensure that reads would be allowed to map to intronic regions of the mouse reference genome (GRCm38/mm10). Nuclei with over 2.5% reads from mitochondria-encoded genes or ribosomal protein genes were filtered out.

**Analysis of snRNA-seq datasets.** Sample integration, count normalization, cell clustering and marker gene identification was performed using Seurat 4.0⁶⁰. Briefly, filtered feature-barcode matrices outputted from CellRanger were imported as a Seurat object using CreateSeuratObject (min cells = 5, min features = 200). Each sample was normalized using SCTransform⁶¹,⁶² and integrated using SelectIntegrationFeatures. PrepSCTIntegration, FindIntegrationAnchors and IntegrateData from the Seurat software package. The integrated dataset, incorporating nuclei from all samples, was subjected to unsupervised clustering using FindNeighbors (dimensions = 1:30) and FindClusters (resolution = 1) from the Seurat package, which executes a shared nearest-neighbour graph clustering algorithm to identify putative cell clusters. Cell type assignment was performed manually based on expression of reported markers²⁹,⁶³,⁶⁴.

Cross-condition differential gene expression analysis was performed based on a 'pseudobulk' strategy: for each cell cluster, gene counts were aggregated to obtain sample-level counts; each pseudo-bulked sample served as an input for edgeR-based differential gene expression analysis⁶⁵,⁶⁶.

For NicheNet-based analysis¹⁹ (v1.1.0), all clusters in our snRNA-seq dataset were used as senders for crypt stem cells, proliferating transit amplifying/stem cells, villus base enterocytes, mid-villus enterocytes and villus tip enterocytes, plus goblet cells. We used the nichenet_seuratobj_aggregate (assay_oi = 'RNA') function with its default settings

to incorporate differential gene expression information from Seurat into our NicheNet analysis and to select bona fide ligand-receptor interactions.

Compass-based in silico metabolic flux analysis²¹ (v0.9.10.2) was performed using transcripts from each of six epithelial cell clusters (crypt stem cells, proliferating transit amplifying cells, villus-base, mid-villus and villus tip enterocytes, and goblet cells). The reaction scores calculated by Compass were filtered based on (1) the confidence levels of the Recon2 reactions and (2) the completeness of information for Recon2 reaction annotations. Only Recon2 reactions that are supported by biochemical evidence (defined by Recon2 as having a confidence level of 4; ref. 22) and that have complete enzymatic information for the reaction were advanced to the follow-on analysis (yield: 2,075 pass filter reactions in 83 Recon2 subsystems).

We subsequently calculated a 'metabolic flux difference' to determine whether the presence or absence of *P. copri* affected Compass-based predictions of metabolic activities at the Recon2 reaction level in the six cell clusters. The 'net reaction score' was calculated as follows:

$$c = c_f - c_r$$

where $c_f$ denotes the Compass score for a given reaction in the 'forward' direction, and, if the biochemical reaction is reversible, $c_r$ denotes the score for the 'reverse' reaction.

A Wilcoxon rank-sum test was used to test the significance of the net reaction score between the two treatment groups. *P* values from the Wilcoxon rank-sum tests were adjusted for multiple comparisons with the Benjamini-Hochberg method.

Cohen's *d* can be used to show the effect size of $c_f$ or $c_r$ for each reaction between two groups²¹ (in our case mice harbouring communities with and without *P. copri*). Briefly, Cohen's *d* of two groups, *j* and *k*, was calculated based on the following two equations⁶⁷ where *n*, *s* and *a* represent the number, the variance and the mean of the observations (in our case, the net reaction scores).

$$s_{\text{pool}} = \sqrt{\frac{(n_j - 1)s_j^2 - (n_k - 1)s_k^2}{n_j + n_k - 2}} \\ d = \frac{a_j - a_k}{s_{\text{pool}}}$$

If both $a_j$ and $a_k$ are non-negative numbers, a positive Cohen's *d* indicates that the mean of group *j* is greater than that of group *k*, whereas a negative Cohen's *d* means the mean of group *j* is smaller in that comparison. The magnitude of Cohen's *d* represents the effect size and is correlated with the difference between the means of the two groups. Because the mean of the net subsystem scores as well as the net reaction scores could be negative, we made the following adjustments to Cohen's *d* to preserve the concordance of sign and the order of group means. The adjusted Cohen's *d* represents the metabolic flux difference *m* and is defined as follows:

$$ \left\{ \begin{array}{l} a_j > 0; a_k < 0; |a_j| < |a_k| : m = -d \\ a_j < 0; a_k > 0; |a_j| > |a_k| : m = -d \\ a_j < 0; a_k < 0 : m = \frac{|a_j|-|a_k|}{s_{\text{pool}}} \\ \text{other} : m = d \end{array} \right. $$

scCODA⁶⁸ (v0.1.8) is a Bayesian probabilistic model for detecting 'statistically credible differences' in the proportional representation of cell clusters, identified from snRNA-seq datasets, between different treatment conditions. This method accounts for two main challenges when analysing snRNA-seq data: (1) low sample number and (2) the compositionality of the dataset (an increase in the proportional representation of a specific cell cluster will inevitably lead to decreases in
the proportional representation of all other cell clusters). Therefore, applying univariate statistical tests, such as a t-test, without accounting for this inherent negative correlation bias will result in reported false positives.

scCODA uses a Bayesian generalized linear multivariate regression model to describe the 'effect' of treatment groups on the proportional representation of each cell cluster; Hamiltonian Monte Carlo sampling is used to calculate the posterior inclusion probability of including the effect of treatment in the model. The type I error (false discovery) is derived from the posterior inclusion probability for each effect. The set of 'statistically credible effects' is the largest set of effects that can be chosen without exceeding a user-defined false discovery threshold α (α = 0.05 by default). We applied scCODA using default parameters, including choice of prior probability in the Bayesian model and the setting for Hamiltonian Monte Carlo sampling. The enteroendocrine cell cluster was used as the reference cluster in accordance with recommendations from the creators of scCODA to choose a cell cluster that has consistent proportional representation across samples.

## Mass spectrometry

**UHPLC-QqQ-MS of caecal glycosidic linkages and GC-MS of short-chain fatty acids.** UHPLC-QqQ-MS quantification of glycosidic linkages and monosaccharides present in caecal glycans was performed using methods described in the accompanying study¹¹. Levels of short-chain fatty acid levels in caecal contents were measured by GC-MS using a procedure outlined in ref. 6.

**LC-MS of acylcarnitines, amino acids and biogenic amines in host tissues.** Acylcarnitines were measured in jejunum, colon, liver, gastrocnemius, quadriceps and heart muscle, plus plasma according to ref. 69, while 20 amino acids plus 19 biogenic amines were quantified in jejunum, liver and muscle according to methods detailed in ref. 70. Plasma levels of non-esterified fatty acids were measured using a UniCel DxC600 clinical analyser (Beckman Coulter).

**Targeted mass spectrometry of caecal amino acids and B vitamins.** Methods for targeted LC-QqQ-MS of amino acids and B vitamins were adapted from a previous publication⁷¹. Caecal samples were extracted with ice-cold methanol, and a 200 µl aliquot was dried (vacuum centrifugation; LabConco CentriVap) and reconstituted with 200 µl of a solution containing 80% methanol in water. A 2 µl aliquot of extracted metabolites was then injected into an Agilent 1290 Infinity II UHPLC system coupled with an Agilent 6470 QqQ-MS operated in positive ion dynamic multiple reaction monitor mode. The native metabolites were separated on HILIC column (ACQUITY BEH Amide, 2.1 × 150 mm, 1.7 µm particle size, Waters) using a 20 min binary gradient with constant flow rate of 0.4 ml min⁻¹. The mobile phases were composed of 10 mM ammonium formate buffer in water with 0.125% formic acid (phase A) and 10 mM ammonium formate in 95% acetonitrile/H₂O (v/v) with 0.125% formic acid (phase B). The binary gradient was as follows: 0–8 min, 91–90% B; 8–14 min, 90–70% B; 15–15.1 min, 70–91% B; 15.1–20 min, 91% B. A pool of 20 amino acids and 7 B vitamins standards with known concentrations (amino acid pool: 0.1 ng ml⁻¹ to 100 µg ml⁻¹; B vitamin pool: 0.01 ng ml⁻¹ to 10 µg ml⁻¹) was injected along with the samples as an external calibration curve for absolute quantification.

## Statistics and reproducibility

Dams were randomly assigned to different treatment groups in all experiments. The sample size of treatment groups in all of the experiments was determined by the size of the litter. No statistical methods were used to pre-determine sample sizes, but our sample sizes are similar to those reported in previous publications⁶,¹³. Both sexes of the offspring were used. No data generated from the experiment were excluded. The statistical tests used in this study do not require data to be normally distributed and do not assume equal variance.

While data collection and analysis were not performed blinded to the experimental group, the initial experiment and validation experiments were performed independently by different members of our team. The snRNA-seq analysis and the subsequent targeted metabolomics quantification were also conducted by different scientists.

## Biological materials

Human specimens and bacterial strains cultured from faecal samples collected from Bangladeshi children are the property of icddr.b. Material transfer agreements exist between icddr.b and Washington University for the use of these samples. Requests for materials should be made to J.I.G.

## Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

## Data availability

Microbial community and long-read bacterial strain genome sequencing datasets, bacterial genome assemblies and microbial RNA-seq and snRNA-seq datasets have been deposited in the National Center for Biotechnology Information Sequence Read Archive (SRA) under study accession number PRJNA1067830. UHPLC-QqQ-MS datasets are available in Glycopost under study accession number GPST000392.

## Code availability

All software used were from publicly available sources.

## References

1. Bäckhed, F. et al. Dynamics and stabilization of the human gut microbiome during the first year of life. *Cell Host Microbe* **17**, 690–703 (2015).

2. Stewart, C. J. et al. Temporal development of the gut microbiome in early childhood from the TEDDY study. *Nature* **562**, 583–588 (2018).

3. Yatsunenko, T. et al. Human gut microbiome viewed across age and geography. *Nature* **486**, 222–227 (2012).

4. Subramanian, S. et al. Persistent gut microbiota immaturity in malnourished Bangladeshi children. *Nature* **510**, 417–421 (2014).

5. Blanton, L. V. et al. Gut bacteria that prevent growth impairments transmitted by microbiota from malnourished children. *Science* **351**, aad3311 (2016).

6. Gehrig, J. L. et al. Effects of microbiota-directed foods in gnotobiotic animals and undernourished children. *Science* **365**, eaau4732 (2019).

7. Brown, E. M. et al. Diet and specific microbial exposure trigger features of environmental enteropathy in a novel murine model. *Nat. Commun.* **6**, 7806 (2015).

8. Chen, R. Y. et al. Duodenal microbiota in stunted undernourished children with enteropathy. *N. Engl. J. Med.* **383**, 321–333 (2020).

9. Malique, A. et al. NAD⁺ precursors and bile acid sequestration treat preclinical refractory environmental enteric dysfunction. *Sci. Transl. Med.* **16**, eabq4145 (2024).

10. Chen, R. Y. et al. A microbiota-directed food intervention for undernourished children. *N. Engl. J. Med.* **384**, 1517–1528 (2021).

11. Hibberd, M. C. et al. Bioactive glycans in a microbiome-directed food for malnourished children. *Nature* **625**, 157–165 (2024).

12. Underwood, M. A., German, J. B., Lebrilla, C. B. & Mills, D. A. *Bifidobacterium longum* subspecies *infantis*: champion colonizer of the infant gut. *Pediatr. Res.* **77**, 229–235 (2015).

13. Barratt, M. J. et al. *Bifidobacterium longum* subsp. *infantis* strains for treating severe acute malnutrition in Bangladeshi infants. *Sci. Trans. Med.* **14**, eabk1107 (2022).

14. Raman, A. S. et al. A sparse covarying unit that describes healthy and impaired human gut microbiota development. *Science* **365**, eaau4735 (2019).
15. Sender, R. & Milo, R. The distribution of cellular turnover in the human body. *Nat. Med.* **27**, 45–48 (2021).

16. Richter, M. & Rosselló-Móra, R. Shifting the genomic gold standard for the prokaryotic species definition. *Proc. Natl Acad. Sci. USA* **106**, 19126–19131 (2009).

17. Olm, M. R., Brown, C. T., Brooks, B. & Banfield, J. F. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. *ISME J.* **11**, 2864–2868 (2017).

18. Karasov, W. H. & Douglas, A. E. Comparative digestive physiology. *Compr. Physiol.* **3**, 741 (2013).

19. Browaeys, R., Saelens, W. & Saeys, Y. NicheNet: modeling intercellular communication by linking ligands to target genes. *Nat. Methods* **17**, 159–162 (2020).

20. Zheng, Y. et al. Intestinal epithelial cell-specific IGF1 promotes the expansion of intestinal stem cells during epithelial regeneration and functions on the intestinal immune homeostasis. *Am. J. Physiol. Endocrinol. Metab.* **315**, E638–E649 (2018).

21. Wagner, A. et al. Metabolic modeling of single Th17 cells reveals regulators of autoimmunity. *Cell* **184**, 4168–4185 (2021).

22. Thiele, I. et al. A community-driven global reconstruction of human metabolism. *Nat. Biotechnol.* **31**, 419–425 (2013).

23. Mihaylova, M. M. Fasting activates fatty acid oxidation to enhance intestinal stem cell function during homeostasis and aging. *Cell Stem Cell* **22**, 769–778 (2018).

24. Crenn, P., Messing, B. & Cynober, L. Citrulline as a biomarker of intestinal failure due to enterocyte mass reduction. *Clin. Nutr.* **27**, 328–339 (2008).

25. Lanyero, B. Correlates of gut function in children hospitalized for severe acute malnutrition, a cross-sectional study in Uganda. *J. Pediatr. Gastroenterol. Nutr.* **69**, 292–298 (2019).

26. Guerrant, R. L. et al. Biomarkers of environmental enteropathy, inflammation, stunting, and impaired growth in children in northeast Brazil. *PLoS ONE* **11**, e0158772 (2016).

27. Mostafa, I. et al. Effect of gut microbiota-directed complementary food supplementation on fecal and plasma biomarkers of gut health and environmental enteric dysfunction in slum-dwelling children with moderate acute malnutrition. *Children* **11**, 69 (2024).

28. Clemente, T. E. & Cahoon, E. B. Soybean oil: genetic approaches for modification of functionality and total content. *Plant Physiol.* **151**, 1030–1040 (2009).

29. Moor, A. E. et al. Spatial reconstruction of single enterocytes uncovers broad zonation along the intestinal villus axis. *Cell* **175**, 1156–1167 (2018).

30. Fawkner-Corbett, D. et al. Spatiotemporal analysis of human intestinal development at single-cell resolution. *Cell* **184**, 810–826 (2021).

31. Kolmgorov, M., Yuan, J. & Pevzner, P. A. Assembly of long, error-prone reads using repeat graphs. *Nat. Biotechnol.* **37**, 540–546 (2019).

32. Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P. & Tyson, G. W. CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. *Genome Res.* **25**, 1043–1055 (2015).

33. Seemann, T. Prokka: rapid prokaryotic genome annotation. *Bioinformatics* **30**, 2068–2069 (2014).

34. Aziz, R. K. et al. SEED servers: high-performance access to the SEED genomes, annotations, and metabolic models. *PLoS ONE* **7**, e48053 (2012).

35. Rodionov, D. A. et al. Micronutrient requirements and sharing capabilities of the human gut microbiome. *Front. Microbiol.* **10**, 1316 (2019).

36. Frolova, M. S., Suvorova, I. A., Iablokov, S. N., Petrov, S. N. & Rodionov, D. A. Genomic reconstruction of short-chain fatty acid production by the human gut microbiota. *Front. Mol. Biosci.* **9**, 949563 (2022).

37. Ashniev, G. A., Petrov, S. N., Iablokov, S. N. & Rodionov, D. A. Genomics-based reconstruction and predictive profiling of amino acid biosynthesis in the human gut microbiome. *Microorganisms* **10**, 740 (2022).

38. Price, M. N. et al. FastTree 2—approximately maximum-likelihood trees for large alignments. *PLoS ONE* **5**, e9490 (2010).

39. Paradis, E. & Schliep, K. ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R. *Bioinformatics* **35**, 526–528 (2018).

40. Yu, G. Using ggtree to visualize data on tree-like structures. *Curr. Protoc. Bioinformatics* **69**, e96 (2020).

41. Pritchard, L., Glover, R. H., Humphris, S., Elphinstone, J. G. & Toth, I. K. Genomics and taxonomy in diagnostics for food security: soft-rotting enterobacterial plant pathogens. *Anal. Methods* **8**, 12–24 (2016).

42. Terrapon, N., Lombard, V., Gilbert, H. J. & Henrissat, B. Automatic prediction of polysaccharide utilization loci in Bacteroidetes species. *Bioinformatics* **31**, 647–655 (2015).

43. Terrapon, N. et al. PULDB: the expanded database of polysaccharide utilization loci. *Nucleic Acids Res.* **46**, D677–D683 (2018).

44. Sokol, H., Pigneur, B., Watterlot, L. & Langella, P. *Faecalibacterium prausnitzii* is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. *Proc. Natl Acad. Sci. USA* **105**, 16731–16736 (2008).

45. Stammler, F. et al. Adjusting microbiome profiles for differences in microbial load by spike-in bacteria. *Microbiome* **21**, 28 (2016).

46. McNulty, N. P. et al. Effects of diet on resource utilization by a model human gut microbiota containing *Bacteroides cellulosilyticus* WH2, a symbiont with an extensive glycobiome. *PLoS Biol.* **11**, e1001637 (2013).

47. R Core Team. R: a language and environment for statistical computing (R Foundation for Statistical Computing, 2020); https://www.R-project.org/

48. Anderson, M. J. Permutational Multivariate Analysis of Variance (PERMANOVA) (Wiley, 2017); https://doi.org/10.1002/9781118445112.stat07841

49. Bates, D., Mächler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 1–48 (2015).

50. Kuznetsova, A., Brockhoff, P. B. & Bojesen-Christensen, R. H. lmerTest: tests in linear mixed effects models. *J. Stat. Softw.* **82**, 1–12 (2017).

51. Krueger F. FelixKrueger/TrimGalore: v0.6.7. Zenodo https://doi.org/10.5281/zenodo.5127899 (2021).

52. Bray, N. L., Pimentel, H., Melsted, P. & Pachter, L. Near-optimal probabilistic RNA-seq quantification. *Nat. Biotechnol.* **34**, 525–527 (2016).

53. Waskom, M. L. seaborn: statistical data visualization. *J. Open Source Softw.* **6**, 3021 (2021).

54. Korotkevich, G. et al. Fast gene set enrichment analysis. Preprint at bioRxiv https://doi.org/10.1101/060012 (2021).

55. Zhang, Y. et al. Statistical approaches for differential expression analysis in metatranscriptomics. *Bioinformatics* **37**, i34–i41 (2021).

56. Klingenberg, H. & Meinicke, P. How to normalize metatranscriptomic count data for differential expression analysis. *Peer J* **5**, e3859 (2017).

57. Bankhead, P. et al. QuPath: open source software for digital pathology image analysis. *Sci. Rep.* **7**, 16878 (2017).

58. Tosti, L. et al. Single-nucleus and in situ RNA-sequencing reveal cell topographies in the human pancreas. *Gastroenterology* **160**, 1330–1344 (2021).

59. Zheng, G. et al. Massively parallel digital transcriptional profiling of single cells. *Nat. Commun.* **8**, 14049 (2017).

60. Hao, Y. et al. Integrated analysis of multimodal single-cell data. *Cell* **184**, 3573–3587 (2021).

61. Hafemeister, C. & Satija, R. Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. *Genome Biol.* **20**, 296 (2019).
62. Choudhary, S. & Satija, R. Comparison and evaluation of statistical error models for scRNA-seq. *Genome Biol.* **23**, 20 (2022).

63. Capdevila, C. et al. Single-cell transcriptional profiling of the intestinal epithelium. *Methods Mol. Biol.* **2171**, 129–153 (2020).

64. Haber, A. L. et al. A single-cell survey of the small intestinal epithelium. *Nature* **551**, 333–339 (2017).

65. Robinson, M. D., McCarthy, D. J. & Smyth, G. K. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics* **26**, 139–140 (2010).

66. Squair, J. W. et al. Confronting false discoveries in single-cell differential expression. *Nat. Commun.* **12**, 5692 (2021).

67. Cohen, J. Statistical Power Analysis for the Behavioral Sciences 2nd edn (Lawrence Erlbaum Associates, 1988).

68. Büttner, M., Ostner, J., Müller, C. L., Theis, F. J. & Schubert, B. scCODA is a Bayesian model for compositional single-cell data analysis. *Nat. Commun.* **12**, 6876 (2021).

69. An, J. et al. Hepatic expression of malonyl-CoA decarboxylase reverses muscle, liver and whole-animal insulin resistance. *Nat. Med.* **10**, 268–274 (2004).

70. Gray, N. et al. High-speed quantitative UPLC-MS analysis of multiple amines in human plasma and serum via precolumn derivatization with 6-aminoquinolyl-N-hydroxysuccinimidyl carbamate: application to acetaminophen-induced liver failure. *Anal. Chem.* **89**, 2478–2487 (2017).

71. Konieczna, L. et al. Bioanalysis of underivatized amino acids in non-invasive exhaled breath condensate samples using liquid chromatography coupled with tandem mass spectrometry. *J. Chromatogr. A* **1542**, 72–81 (2018).

## Acknowledgements

We thank D. O'Donnell and M. Karlsson for their invaluable assistance with mouse husbandry, M. Meier for his essential role in generating shotgun sequencing and microbial RNA-seq datasets, J. Guruge for her help with culturing and maintaining bacterial strains and S. Deng and J. Serugo for oversight of biospecimen archives. We also appreciate Washington University Advanced Imaging and Tissue Analysis Core of the Digestive Disease Research Core Center (NIDDK P30 DK052574) for the support in preparation of the histology slides. This work was supported by grants from the National Institutes of Health (NIH) (DK30292), the Washington University Personalized Medicine Initiative and the Bill & Melinda Gates Foundation (INV-016367). E.M.L. is supported by a NIH training grant (T32 HG000045). Y.W. is the recipient of a career development award from the NIH (K01 DK134840). C.Z. is the recipient of a NIH F30 predoctoral MD/ PhD fellowship (DK131866). K.M.P. is supported by a postdoctoral fellowship from the Helen Hay Whitney Foundation. C.K. is the recipient of a NIH F30 predoctoral MD/PhD fellowship (DK123838).

## Author contributions

H.-W.C., E.M.L., Y.W., C.Z. and J.I.G. designed the gnotobiotic mouse experiments, which were performed by H.-W.C., E.M.L. and C.Z. Mouse diets were formulated by H.-W.C. with assistance from M.J.B. and were based on the results of human diet surveys provided by I.M, S.D., M.M. and T.A. Bacterial genomes were sequenced and assembled

by D.M.W., H.-W.C. and M.C. with S.H., D.A.R., A.A.A., N.T., B.H. and A.L.O. contributing to genome annotation. Shotgun sequencing and microbial RNA-seq datasets were generated by H.-W.C. and E.M.L. and analysed together with M.C.H., H.M.L. and D.M.W. Mass spectrometry-based metabolomic studies were conducted by J.C., O.I., M.J.M. and C.B.N. Glycan analyses were performed by J.J.C., G.C., Y.C., N.P.B. and C.B.L. Morphometric analysis was conducted by Y.W. snRNA-seq datasets were generated and analysed by H.-W.C., Y.W., K.M.P., R.Y.C., C.K. and M.C. H.-W.C., E.M.L., Y.W. and J.I.G. wrote this manuscript with invaluable assistance from co-authors.

## Competing interests

A.L.O. and D.A.R. are co-founders of Phenobiome, a company pursuing development of computational tools for predictive phenotype profiling of microbial communities. C.B.L is a co-founder of Infinant Health, interVenn Bio and one.bio, companies involved in the characterization of glycans and developing carbohydrate applications for human health. The remaining authors of this paper declare no competing interests.

## Additional information

**Extended data** is available for this paper at
https://doi.org/10.1038/s41564-024-01628-7.

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41564-024-01628-7.

**Correspondence and requests for materials** should be addressed to Jeffrey I. Gordon.

**Peer review information** *Nature Microbiology* thanks the anonymous reviewer(s) for their contribution to the peer review of this work.

**Reprints and permissions information** is available at
www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2024

<footnote n="1">¹Edison Family Center for Genome Sciences and Systems Biology, Washington University School of Medicine, St. Louis, MO, USA. ²Center for Gut Microbiome and Nutrition Research, Washington University School of Medicine, St. Louis, MO, USA. ³Department of Pathology and Immunology, Washington University School of Medicine, St. Louis, MO, USA. ⁴Architecture et Fonction des Macromolécules Biologiques, CNRS, Aix-Marseille University, Marseille, France. ⁵Infectious and Inflammatory Disease Center, Sanford Burnham Prebys Medical Discovery Institute, La Jolla, CA, USA. ⁶Department of Chemistry, University of California, Davis, CA, USA. ⁷Department of Biotechnology and Biomedicine (DTU Bioengineering), Technical University of Denmark, Lyngby, Denmark. ⁸Department of Biological Sciences, King Abdulaziz University, Jeddah, Saudi Arabia. ⁹Sarah W. Stedman Nutrition and Metabolism Center, Duke University Medical Center, Durham, NC, USA. ¹⁰Duke Molecular Physiology Institute, Duke University Medical Center, Durham, NC, USA. ¹¹Department of Medicine, Duke University Medical Center, Durham, NC, USA. ¹²Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, NC, USA. ¹³International Centre for Diarrhoeal Disease Research, Bangladesh (icddr,b), Dhaka, Bangladesh. ¹⁴These authors contributed equally: Hao-Wei Chang, Evan M. Lee, Yi Wang. ✉e-mail: jgordon@wustl.edu</footnote>
<figure><img src="image_21.png" /></figure>

<figure><img src="image_22.png" /></figure>

<figure><img src="image_23.png" /><figcaption>**Extended Data Fig. 1** | Determining the relationship between *P. copri* colonization efficiency and pre-colonization with *B. longum* subsp. *infantis*.</figcaption></figure>

(a) Phylogenetic tree of cultured *P. copri* isolates used in the mouse studies described and the two MAGs positively associated with WLZ in the randomized controlled human study. The phylogenetic distance between each pair of comparisons is shown in the matrix. (b) Experimental design (n = 9 offspring per

group in experiment 1 and n = 10 offspring per group in experiment 2). Mice were weaned at P28 and P25 for experiments 1 and 2, respectively. (c) Total absolute abundances of *P. copri* strains in fecal samples collected from mice at P42. Mean values ± SD are shown. Each dot represents a separate animal. *P*-values were calculated using a two-sided Mann-Whitney U test and are 4x10⁻⁴ and 3x10⁻⁴ for experiments 1 and 2, respectively.
<figure><img src="image_24.png" /><figcaption>Extended Data Fig. 2 | Absolute abundances of other bacterial species in the defined community across different time points and locations.</figcaption></figure>

(a) Absolute abundances of organisms which were significantly higher with either *B. infantis* Bg2D9 (Arm 1 vs. Arm 2), or with the combination of *B. infantis* Bg2D9 and *Prevotella* (Arm 1 vs. Arm 3) in fecal samples collected at P53 (n = 8, 5, and 7 offspring for arms 1, 2, and 3, respectively). The adjusted *P*-values

for *B. obeum* (Arm 1 vs. Arm 3) and *D. longicatena* (Arm 1 vs. Arm 2) are 2x10⁻⁴ and 7x10⁻⁴, respectively. (b, c) Absolute abundances of the same organisms in cecal contents collected at P53 (panel b) and in fecal samples collected before weaning at P21 (panel c). Mean values ± SD are shown. Each dot represents an individual animal. *P*-values were calculated by the Kruskal-Wallis test followed by post-hoc Dunn's test with Bonferroni correction. N.S., not significant (*P* > 0.05).
<figure><img src="image_25.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_26.png" /></figure>

<figure><img src="image_27.png" /></figure>

<figure><img src="image_28.png" /></figure>

<figure><img src="image_29.png" /></figure>

<figure><img src="image_30.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_32.png" /></figure>

<figure><img src="image_33.png" /></figure>

<figure><img src="image_34.png" /></figure>

B. infantis Bg463
- Prevotella

**Extended Data Fig. 3** | Targeted mass spectrometric and microbial RNA-Seq analyses of consortia of cultured age-discriminatory and WLZ-associated bacteria strains that colonized gnotobiotic mice. Cecal contents collected at the end of the experiment described in Fig. 1a were analyzed. (a,b) UHPLC-QqQ-MS-based quantitation of levels of total arabinose (panel a) and arabinose-containing glycosidic linkages (panel b) in cecal glycans collected at P53. Abbreviations: Araf, arabinofuranose; Arap, arabinopyranose; Xylp, xylopyranose. Each dot represents an individual animal (n = 8, 5, and 7 offspring for arms 1, 2, and 3, respectively). Mean values ± s.d. are shown. *P* values were calculated by the Kruskal-Wallis test followed by post-hoc Dunn's test with Bonferroni correction for panels a and b. (c) PCA of profiles of normalized

meta-transcriptomic counts (see Methods). Centroids are denoted by a colored 'X' for each group. *P*-values were calculated by PERMANOVA. (d) MTXmodel abundance-normalized differential expression analysis of genes involved in specific carbohydrate utilization and amino acid biosynthetic pathways in the four arabinose-utilizing bacteria. Violin plots show the distribution of log₂ fold-differences for all expressed genes with metabolic pathway annotations in the indicated organism. Dots in panel d represent differential expression test results for individual genes involved in the corresponding pathway and are coloured if their Benjamini-Hochberg adjusted *P*-value is less than 0.1 (see Methods). Abbreviations: BCAA, branched-chain amino acid; Glu, glutamate; Gln, glutamine.
<figure><img src="image_35.png" /></figure>

**Extended Data Fig. 4** | Expression of *P. copri* and *P. stercorea* PULs and targeted mass spectrometric analysis of their predicted targets. (a,b) GSEA of expression of PULs shared by *P. copri* Bg131 (panel a) and *P. stercorea* (panel b) in the two *Prevotella*-containing arms of the experiment described in Fig. 1 (n = 8 and 5 offspring for arms 1 and 2 respectively). Benjamini-Hochberg adjusted *P*-values were calculated using GSEA ranking genes by their mean log~2~ TPM across the *Prevotella*-colonized samples in Arms 1 and 2, with each PUL comprising a gene set against the background of all predicted PUL genes. Violin plots show the

log~2~ TPM of all genes assigned to any of the PULs in each isolate; plots are split to show the indicated PUL. The exact adjusted *P*-value for PUL27a in Arm 2 (panel a) is 3x10~−4~, (c) UHPLC-QqQ-MS-based quantitation of levels of total mannose, N-acetylglucosamine, and N-acetylgalactosamine in cecal glycans collected at P53 (n = 8, 5, and 7 offspring for arms 1, 2, and 3 respectively). Mean values ± SD are shown. Each dot represents an individual animal. *P*-values were calculated by the Kruskal-Wallis test followed by post-hoc Dunn's test with Bonferroni correction.
<figure><img src="image_36.png" /></figure>

<figure><img src="image_37.png" /><figcaption>**Extended Data Fig. 5** | snRNA-Seq analysis of differential intestinal gene expression in mice colonized with bacterial consortia with or without *P. copri*. Jejunal tissue samples collected from 'w/ *P. copri* & w/ *B. infantis* Bg2D9' and 'w/o *P. copri* & w/ *B. infantis* Bg463' groups at the end of the experiment described in Fig. 1 were analyzed. (a) Dot plot of marker gene expression across</figcaption></figure>

<figure><img src="image_38.png" /></figure>

epithelial cell types. The average expression level and percentage of nuclei that express a given gene within a cell type are indicated by dot color and size, respectively. (b) Integrated UMAP plot for all jejunal nuclei isolated from animals in both arms (n = 4 mice/arm). (c) The number and directionality of statistically significant differentially expressed genes in each cell cluster.
<figure><img src="image_39.png" /><figcaption>**Extended Data Fig. 6** | *NicheNet-based analysis of the effects of P. copri colonization on cell-cell signaling activities.* Each row represents different sender cell clusters. Each column represents ligands expressed by these sender cells. Cells are colored based on the log₂-fold difference in expression of ligands</figcaption></figure>

in the sender cell clusters between mice in 'w/ P. copri & w/ B. infantis Bg2D9' and 'w/o P. copri & w/ B. infantis Bg463' groups from the experiment described in Fig. 1. Ligands (columns) are grouped based on receiver cell clusters and the indicated functions of downstream signaling pathways in these receiver cells.
<figure><img src="image_40.png" /></figure>

**Extended Data Fig. 7** | Normalized number of Recon2 reactions in Recon2 subsystems predicted to have statistically significant differences in their activities between mice in the 'w/ P. copri & w/ B. infantis Bg2D9' and 'w/o P. copri & w/ B. infantis Bg463' groups. See legend to Fig. 2b, which shows other affected subsystems, for details.
a

<figure><img src="image_41.png" /></figure>

b

<figure><img src="image_42.png" /></figure>

<figure><img src="image_43.png" /></figure>

<figure><img src="image_44.png" /></figure>

e

<figure><img src="image_45.png" /></figure>

d

**Extended Data Fig. 8** | Validating the effects of *P. copri* colonization on postnatal weight gain and host metabolism in gnotobiotic dam-pup dyads. (a) Study design (n = 4 dams and 18 and 19 offspring for arms 1 and 2, respectively). (b) Body weights of the offspring of dams, normalized to postnatal day 23 [linear mixed effects model (see Methods)]. (c-e) Targeted mass spectrometric analysis of jejunal citrulline (panel c) and acylcarnitine levels (panel d), plus colonic acylcarnitine levels (panel e). Exact *P*-values for jejunal C3

and jejunal C18:1 (panel d) are $1\times10^{-5}$ and $7\times10^{-4}$, respectively. Exact *P*-values for colonic C4/Ci4, C5, C16, and C18:1 (panel e) are $3\times10^{-6}$, $2\times10^{-6}$, $3\times10^{-4}$, and $7\times10^{-4}$, respectively. (f) Plasma levels of non-esterified fatty acids. Each dot represents a single animal. Mean values ± SD are shown for panels b-f. *P*-values were calculated from the linear mixed effect model (panel b) or two-sided Mann-Whitney U test (panels c-f). N.S., *P*-value > 0.05.
<figure><img src="image_46.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_47.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_48.png" /></figure>

**Extended Data Fig. 9** | Evaluating the effect of diet on the defined community

**in gnotobiotic dam-pup dyads.** (a) Experimental design (n = 2 dams and 12 offspring/diet treatment). (b) Principal component analysis showing the significant differences in community structure in the cecums of mice euthanized on P53 (P = 1x10⁻⁵; PERMANOVA). Ellipses represent 95% confidence intervals. (c) Absolute abundances of the defined community members in

cecal contents at P53. Exact *P*-values for *B. breve*, *B. catenulatum*, *B. infantis* Bg2D9, *D. formicigenerans*, *E. coli*, *F. prausnitzii*, *L. garavieae*, *L. ruminis*, and *M. multacida* are 4×10⁻⁵, 4×10⁻⁵, 4×10⁻⁵, 6×10⁻⁴, 4×10⁻⁵, 2×10⁻⁴, 4×10⁻⁵, 1×10⁻⁴, and 4×10⁻⁵, respectively. *P*-values were calculated by PERMANOVA (panel b) and a two-sided Mann-Whitney U test (panel c). N.S., *P*-value > 0.05. Mean values ± SD are shown. Each dot represents an individual animal.
<figure><img src="image_49.png" /></figure>

**Extended Data Fig. 10** | LC-MS of ileal and colonic acylcarnitines in gnotobiotic mice colonized with *P. copri* BgD5_2 and BgF5_2. (a) LC-MS of ileal acylcarnitines corresponding to soybean oil lipids. (b) LC-MS of colonic

<figure><img src="image_50.png" /></figure>

acylcarnitines corresponding to soybean oil lipids. Each dot represents an individual animal (n = 13 animals in each group). Mean values ± SD are shown. *P*-values were calculated using a two-sided Mann-Whitney U test for panel a and b.
# Reporting Summary

Nature Portfolio wishes to improve the reproducibility of the work that we publish. This form provides structure for consistency and transparency in reporting. For further information on Nature Portfolio policies, see our [Editorial Policies](#) and the [Editorial Policy Checklist](#).

## Statistics

For all statistical analyses, confirm that the following items are present in the figure legend, table legend, main text, or Methods section.

<figure><img src="image_51.png" /></figure>

Our web collection on [statistics for biologists](#) contains articles on many of the points above.

## Software and code

Policy information about [availability of computer code](#)

**Data collection** No software was used for data collection

**Data analysis** All software packages used are described in Methods and are publicly available; they include ape (v5.6-2), bcl2fastq (v1.8.4), Trim Galore (v0.6.4), Kallisto (v0.43.0), MTXmodel (v1.2.3), edgeR (v3.32.1), seaborn (v0.12.1), fgsea (v1.20.0), QuPath (v0.3.2), CellRanger (v5.0), Seurat (v4.0), NicheNet (v1.1.0), Compass (v0.9.10.2), scCODA (v0.1.8), Flye (v2.8.1), CheckM (v1.1.3), Prokka (v1.14), and pyani (v0.2.10). Code detailing the steps in the MAG assembly workflow and analyses of microbial RNA-Seq and glycan datasets are available from GitLab (https://gitlab.com/hibberdm/hibberd_webber_et_al_mdcf_poc_mags) and have been accessioned at Zenodo (DOI:10.5281/ zenodo.8000098). Code for annotation of bacterial genes and prediction of metabolic phenotypes is available from GitHub (https:// github.com/rodionovdima/PhenotypePredictor) and has been accessioned at Zenodo (DOI: 10.5281/zenodo.10049439).

For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers. We strongly encourage code deposition in a community repository (e.g. GitHub). See the Nature Portfolio [guidelines for submitting code & software](#) for further information.

## Data

Policy information about [availability of data](#)

All manuscripts must include a [data availability statement](#). This statement should provide the following information, where applicable:

- Accession codes, unique identifiers, or web links for publicly available datasets

- A description of any restrictions on data availability

- For clinical datasets or third party data, please ensure that the statement adheres to our [policy](#)

Microbial community and long-read bacterial strain genome sequencing datasets, bacterial genome assemblies, and microbial RNA-Seq and snRNA-Seq datasets
# Field-specific reporting

Please select the one below that is the best fit for your research. If you are not sure, read the appropriate sections before making your selection.

<figure><img src="image_52.png" /></figure>

Life sciences

<figure><img src="image_53.png" /></figure>

Behavioural & social sciences

<figure><img src="image_54.png" /></figure>

Ecological, evolutionary & environmental sciences

For a reference copy of the document with all sections, see [nature.com/documents/nr-reporting-summary-flat.pdf](https://nature.com/documents/nr-reporting-summary-flat.pdf)

# Life sciences study design

All studies must disclose on these points even when the disclosure is negative.

## Sample size

The number of animals studied per treatment group was based on the number of pups born to each litter. Sample sizes were based off past preclinical studies on effects of MDCF prototypes on microbial communities and the effects of B. infantis on weight gain (ref. 6 and 13). The exact number of animals used in each experiment is noted in the main text, figure legends and methods.

## Data exclusions

No data were excluded

## Replication

Two similar validation experiments were performed (Extended Data Fig. 8 and Fig. 3) which reproduced the key findings of weight gain with P. copri colonization and changes in host and microbial community metabolism. The experiment described in Extended Data Fig. 1 was repeated independently twice confirming the colonization dependency of P. copri on B. infantis.

## Randomization

Before initiation of the gnotobiotic experiments, dams were randomly assigned to each treatment group. Processing of samples from gnotobiotic mice alternated experimental groups to reduce potential technical biases.

## Blinding

All data were generated from all samples without knowledge of treatment group

# Reporting for specific materials, systems and methods

We require information from authors about some types of materials, experimental systems and methods used in many studies. Here, indicate whether each material, system or method listed is relevant to your study. If you are not sure if a list item applies to your research, read the appropriate section before selecting a response.

## Materials & experimental systems

<table><thead><tr><th>n/a</th><th>Involved in the study</th></tr></thead><tbody><tr><td>☑</td><td>Antibodies</td></tr><tr><td>☑</td><td>Eukaryotic cell lines</td></tr><tr><td>☑</td><td>Palaeontology and archaeology</td></tr><tr><td>☐</td><td>Animals and other organisms</td></tr><tr><td>☐</td><td>Human research participants</td></tr><tr><td>☑</td><td>Clinical data</td></tr><tr><td>☑</td><td>Dual use research of concern</td></tr></tbody></table>

## Methods

<table><thead><tr><th>n/a</th><th>Involved in the study</th></tr></thead><tbody><tr><td>☑</td><td>ChIP-seq</td></tr><tr><td>☑</td><td>Flow cytometry</td></tr><tr><td>☑</td><td>MRI-based neuroimaging</td></tr></tbody></table>

## Animals and other organisms

Policy information about studies involving animals; ARRIVE guidelines recommended for reporting animal research

### Laboratory animals

Germ-free C57BL/6J mice and their pups were used. Animals were euthanized at P42 for experiments testing colonization dependency of P. copri and B. infantis and at P53 for all other experiments.

Germ-free C57BL/6J mice were maintained in plastic flexible film isolators (Class Biologically Clean Ltd) at 23 °C under a strict 12-hour light cycle (lights on at 0600h). Autoclaved paper 'shepherd shacks' were kept in each cage to facilitate natural nesting behaviors and provide environmental enrichment. For all experiments, pregnant C57BL/6J mice originating from trio matings were given ad libitum access to an autoclaved breeder chow (Purina Mills; Lab Diet 5021) throughout their pregnancy and until postpartum day 2. Mice then underwent the sequence of diet changes described in each experiment's design with ad libitum food access.

### Wild animals

No wild animals were used

### Field-collected samples

No field-collected samples were used

### Ethics oversight

All gnotobiotic mouse experiments were performed following IACUC and IBC protocols that were approved by the Washington University Animal Studies and Environmental Health and Safety Committee.
## Human research participants

Policy information about studies involving human research participants

**Population characteristics**
The human biospecimens used for the gnotobiotic mouse experiments described in this manuscript are from an already reported clinical study whose population characteristics are detailed in ref. 6 (N Engl J Med 2021 (DOI: 10.1056/NEJMoa2023294)).

**Recruitment**
The human biospecimens used for the gnotobiotic mouse experiments described in this manuscript are from an already reported clinical study whose protocols for participant recruitment are described in ref. 6 (N Engl J Med 2021 (DOI: 10.1056/NEJMoa2023294)).

**Ethics oversight**
The human biospecimens used for the gnotobiotic mouse experiments described in this manuscript are from an already reported clinical study that was approved by the ethical review committee of the International Centre for Diarrhoeal Disease Research, Bangladesh as described in ref. 6 (N Engl J Med 2021 (DOI: 10.1056/NEJMoa2023294)).

Note that full information on the approval of the study protocol must also be provided in the manuscript.