| Computational Biology | Research Article Food-breastmilk combinations alter the colonic microbiome of weaning infants: an in silico study Vitor G. da Silva,1,2 Nick W. Smith,1 Jane A. Mullaney,1,2,3 Clare Wall,2,4 Nicole C. Roy,1,2,5 Warren C. McNabb1,2 AUTHOR AFFILIATIONS See affiliation list on p. 16. ABSTRACT The introduction of solid foods to infants, also known as weaning, is a critical point for the development of the complex microbial community inhabiting the human colon, impacting host physiology in infancy and later in life. This research investigated in silico the impact of food-breastmilk combinations on growth and metabolite production by colonic microbes of New Zealand weaning infants using the metagenome-scale metabolic model named Microbial Community. Eighty-nine foods were individually combined with breastmilk, and the 12 combinations with the strongest influence on the microbial production of short-chain fatty acids (SCFAs) and branched- chain fatty acids (BCFAs) were identified. Fiber-rich and polyphenol-rich foods, like pumpkin and blackcurrant, resulted in the greatest increase in predicted fluxes of total SCFAs and individual fluxes of propionate and acetate when combined, respectively, with breastmilk. Identified foods were further combined with other foods and breastmilk, resulting in 66 multiple food-breastmilk combinations. These combinations altered in silico the impact of individual foods on the microbial production of SCFAs and BCFAs, suggesting that the interaction between the dietary compounds composing a meal is the key factor influencing colonic microbes. Blackcurrant combined with other foods and breastmilk promoted the greatest increase in the production of acetate and total SCFAs, while pork combined with other foods and breastmilk decreased the production of total BCFAs. IMPORTANCE Little is known about the influence of complementary foods on the colonic microbiome of weaning infants. Traditional in vitro and in vivo microbiome methods are limited by their resource-consuming concerns. Modeling approaches represent a promising complementary tool to provide insights into the behavior of microbial communities. This study evaluated how foods combined with other foods and human milk affect the production of short-chain fatty acids and branched-chain fatty acids by colonic microbes of weaning infants using a rapid and inexpensive in silico approach. Foods and food combinations identified here are candidates for future experimental investigations, helping to fill a crucial knowledge gap in infant nutrition. KEYWORDS gut microbiome, food, infant, in silico T he human gastrointestinal tract is colonized by a complex microbial community, with the greatest concentration and diversity being found in the large intestine, or colon (1). Diet is a well-known key factor shaping the composition and function of colonic microbes throughout human life (2–4). In turn, colonic microbes impact host physiology and are associated with health and disease biomarkers (5). One mechanism by which the colonic microbiota affects host health is the production of metabolites that are later absorbed by the host (6). Organic acids, such as short-chain fatty acids (SCFAs) and branched-chain fatty acids (BCFAs), are among the most studied microbial metabolites September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 1 Editor Daniel Garrido, Pontificia Universidad Catolica de Chile, Santiago, Santiago, Chile Address correspondence to Warren C. McNabb, w.mcnabb@massey.ac.nz. Vitor G. da Silva, Nick W. Smith, Jane A. Mullaney, Clare Wall, Nicole C. Roy, and Warren C. McNabb contributed equally to this article. Author order was determined in order of increasing seniority. The authors declare no conflict of interest. See the funding table on p. 16. Received 22 April 2024 Accepted 22 July 2024 Published 27 August 2024 Copyright © 2024 Geniselli da Silva et al. This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license. D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://crossmark.crossref.org/dialog/?doi=10.1128/msystems.00577-24&domain=pdf&date_stamp=2024-08-27 https://doi.org/10.1128/msystems.00577-24 https://creativecommons.org/licenses/by/4.0/ produced in the colon, and decreases in their production have been associated with disease, particularly gastrointestinal and metabolic disorders (7, 8). One of the challenges of investigating the colonic microbiome is that the composi­ tion of this microbial community is unique to each individual, making it difficult to define an ideal microbiota in terms of composition (9, 10). On the other hand, the functional capacity of the microbiota is similar among individuals with the same health status (11), and colonic microbes are functionally redundant, meaning that different taxa can perform the same metabolic functions (12). Thus, focusing on what the colonic microbiota produces, rather than which microbes compose it, may be more relevant to host health implications, also reducing the complexity of microbiome investigations. Animal and human studies have demonstrated that imbalances in the colonic microbiota, or dysbiosis, during infancy are associated with an increased risk for several diseases later in life (13–16). Therefore, adequate nutrition from an early postnatal age is crucial to prevent or limit dysbiosis. The introduction of solid foods represents a window of opportunity for establishing long-term beneficial host-microbiota interactions that may influence host physiology later in life (17, 18). However, among the studies that examined the impact of solid foods on the colonic microbiota of infants (19–24), only a few assessed the microbial production of organic acids (25, 26), limiting our understand­ ing of how complementary foods affect colonic microbes in this crucial life stage. Studies conducted in vitro or using animal models have investigated the effect of foods on the colonic microbiome of infants in the absence of human milk (27–30), not accurately reflecting how complementary foods are introduced to weaning infants. Furthermore, there is no evidence in vitro, in animals (such as mice and piglets), or in silico of how food-breastmilk combinations affect the colonic microbiota in early life. Mathematical models are a rapid and inexpensive complementary strategy to study microbial communities, being able to generate hypotheses that can be further tested by experimental approaches. Genome-scale metabolic models stand out for using genome information to infer the metabolism of microorganisms, having the major advantage of predicting microbial growth rates and fluxes of produced metabolites without prior definition of kinetic parameters, whose experimental determination is technically challenging (31). Recently, a metagenome-scale community metabolic model named Microbial Community (MICOM) was proposed to predict how diets shape the composition and function of the human colonic microbiota (29). So far, in silico investigations of the colonic microbiome using metagenome-scale metabolic models have mainly focused on adults (32–35). This study used this modelling approach to identify foods with the strongest impact on the production of organic acids by colonic microbes of weaning infants when combined with breastmilk. SCFAs and BCFAs were chosen due to their association with the host health. Insights generated by this research help the design of future experiments, advancing the understanding of the relationship between comple­ mentary foods and colonic microbes in a decisive stage of human life. RESULTS Predicted fluxes of SCFAs and BCFAs This study aimed to verify in silico the effect of complementary foods on the colonic microbiota of New Zealand weaning infants. For that, a literature survey was conducted to identify 89 foods from various food groups that are commonly introduced to weaning infants, particularly in New Zealand (36–40). These foods were combined in silico with breastmilk to create combinations representative of those consumed by 6-month-old infants (15% food and 85% breastmilk, by caloric intake). We further evaluated in silico the effect of these food-breastmilk combinations on the 10 most abundant genera composing the average fecal microbiota of New Zealand infants at weaning age (5–12 months of age), using previously published data (28). Candidate foods with the strongest impact on the microbial production of organic acids (SCFAs and BCFAs) were identified. These candidate foods were chosen to provide a diversity of model outcomes, including Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 2 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 foods that increased, decreased, or had little effect in total SCFAs (acetate, propionate, and butyrate combined) or total BCFAs (isobutyrate and isovalerate combined) fluxes compared to breastmilk alone. Importantly, the employed modeling approach results in repeatable outcomes when simulations are performed with the same initial conditions. Consequently, statistical analyses could not be performed to assess any potential variation in microbial growth or metabolite production resulting from the use of different food-breastmilk combinations. Another limitation of this study was using a 1% relative abundance threshold to select the microbial genera included in the simulations, resulting in a small-scale community that lacked low-abundant taxa and may not fully represent the metabolic potential of the colonic microbiota of weaning infants. Thus, the findings presented are preliminary and require further experimental validation. A graphical representation of the modeling workflow is depicted in Fig. 1. Combining foods with breastmilk altered the fluxes of major SCFAs and BCFAs compared to breastmilk alone (Table S1). Twenty-nine food-breastmilk combinations increased the flux of total SCFAs (32% of the evaluated combinations), and 64 food- breastmilk combinations reduced the flux of total BCFAs (72% of the combinations). Twelve food-breastmilk combinations that, compared to breastmilk alone, resulted in the greatest increase, decrease, or similar values of total SCFAs and BCFAs fluxes were identified (Fig. 2). The combination of pumpkin with breastmilk promoted the greatest increase in the flux of total SCFAs (11.7%) and the second greatest decrease in the flux of total BCFAs (40.2%), compared to breastmilk alone. Raspberries-breastmilk and blackcurrant-breast­ milk increased the production of total SCFAs (6.4 and 6.2%, respectively) with no change in the flux of total BCFAs. Similarly, soybean-breastmilk and sweet potato-breastmilk elevated fluxes of total SCFAs (4.0 and 3.0%, respectively) without changing total BCFA production. The increase in the flux of total SCFAs obtained with the pumpkin-breastmilk combination was driven by heightened fluxes of propionate and butyrate. On the other hand, the other four food-breastmilk combinations increasing total SCFA flux increased the flux of acetate. Indeed, combining pumpkin with breastmilk resulted in the largest increase in propionate (65.0%), while breastmilk combined with blackcurrant in the largest increase in acetate (13.3%). The couscous-breastmilk combination resulted in the highest increase in butyrate (104.8%) but decreased acetate and propionate individually. In addition, it caused the largest reduction in total BCFA flux (61.7%). Split peas-breastmilk reduced the flux of FIG 1 Workflow to predict the influence of food-breastmilk combinations on growth rates of colonic microbes of weaning infants and produced fluxes of organic acids using the model MICOM. Image created using Biorender.com. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 3 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 total BCFAs by 26.9% while increasing the production of acetate (12.8%) and butyrate (14.8%). On the other hand, the combination of strawberry with breastmilk promoted the greatest reduction in the total SCFA flux (60.9%) also reducing total BCFA flux (33.9%), while shrimp-breastmilk had the greatest increase in the production of BCFAs (22.9%). Conversely, black beans, pork, and chickpeas promoted little to no change in the fluxes of SCFAs and BCFAs when individually combined with breastmilk. Predicted microbial growth rates The modeling approach makes predictions about microbial growth assuming that the system is in a steady state, meaning that there is no accumulation of substrates. In addition, MICOM employs a linearization strategy that correlates the predicted growth rates of individual taxa to their relative abundance. As a result, the profile of faster- growing microbes was not expected to change much under different food-breastmilk combinations. High-abundance genera were expected to have the greatest growth rates, while low-abundance genera were expected to have slower growth rates. As expected, the high-abundance genera Bifidobacterium, Bacteroides, and Bacil­ lus had higher growth rates for most of the food-breastmilk combinations, while FIG 2 Heatmap of food-breastmilk combinations with the greatest influence on predicted fluxes of SCFAs and BCFAs. Fluxes of organic acids are expressed in relative variation in comparison to breastmilk alone. Cells are colored according to intensity, with the highest values in green and the lowest values in red. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 4 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 low-abundance genera, such as Lacticaseibacillus and Streptococcus, had slower growth rates (Table S2). On the other hand, low-abundance genera also had negligible growth (growth rates lower than 10−6). This result suggests a failure to respect the trade-off between community growth and individual growth, resulting in a numerical problem that was not optimally solved (numerical instabilities, see Materials and Methods section). Food-breastmilk combinations with the greatest influence on the fluxes of SCFAs and BCFAs impacted the microbial growth rates of high-abundance genera in different ways (Fig. 3). Compared to breastmilk alone, little to no changes in the growth rates of the genus Bacteroides were observed for these breastmilk-food combinations. On the other hand, most of the food-breastmilk combinations reduced the growth of the Bifidobacte­ rium and Lactobacillus genera. The greatest reduction in these genera was observed for couscous-breastmilk (17% reduction), followed by shrimp-breastmilk (16%). In addition, a few food-breastmilk combinations altered the growth of Prevotella, Collinsella, and Bacillus genera. The pork-breastmilk combination had the strongest influence on the growth rates of these genera, increasing it by 48% for the Prevotella genus and 35% for the Bacillus genus, while decreasing growth rates by 53% for the Collinsella genus. Impact of multiple food-breastmilk combinations on colonic microbes of weaning infants As complementary foods are normally introduced to infants in combination with other foods rather than consumed individually, the identified foods with the strongest impact on microbial SCFAs and BCFAs fluxes were combined with other identified foods and breastmilk for additional simulations. Sixty-six combinations were generated (multiple food-breastmilk combinations), composed of 7.5% of a first food item, 7.5% of a second food item, and 85% of breastmilk by caloric intake. In comparison to breastmilk alone, 24 multiple food-breastmilk combinations increased the flux of total SCFAs (36% of the FIG 3 Heatmap of microbial growth rates for food-breastmilk combinations with the greatest influence on predicted fluxes of SCFAs and BCFAs. Microbial growth rates are expressed in relative variation in comparison to breastmilk alone. Cells are colored according to intensity, with the highest values in green and the lowest values in red. Genera with predicted negligible growth are not presented. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 5 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 combinations), while 47 reduced the flux of total BCFAs (70% of the combinations, Table S3). The combination of blackcurrant with breastmilk and soybean, strawberries, or sweet potato promoted the highest increases in acetate production, leading to the greatest increases in total SCFA flux by 11.7%, 6.9%, and 6.6%, respectively, as compared to breastmilk (Fig. 4). Pumpkin combined with breastmilk and couscous, or split peas resulted in the greatest increase in propionate (9.4% and 8.4%, respectively) but did not alter total SCFA flux. On the other hand, the combination of black beans-pumpkin- breastmilk had the greatest increase in butyrate flux (410.9%). However, the same combination had the greatest reduction in isobutyrate, isovaler­ ate, and resulting total BCFA flux (84.5%). The combination chickpea-split peas-breast­ milk followed behind, with a total BCFA flux reduction of 75.6%. Pork combined with breastmilk and other foods also decreased total flux of BCFAs; reduction was observed for the pork-couscous-breastmilk (26.4%), pork-chickpea-breastmilk (36.5%), and pork-blackcurrant-breastmilk combinations (47.7%). Combining foods with other foods and breastmilk shifted the way foods influenced the production of total SCFAs and total BCFAs when combined with breastmilk only. Among the foods individually combined with breastmilk, pumpkin promoted the greatest increase in the flux of total SCFAs, while combining pumpkin with other foods and breastmilk, such as pumpkin-split peas-breastmilk and couscous-pumpkin-breast­ milk, promoted little to no changes in the flux of total SCFAs (Fig. 5). In turn, among the multiple food-breastmilk combinations, it was blackcurrant that stood out for promoting the greatest increases in total SCFA flux when combined with breastmilk and soybean, strawberries, or sweet potato. Most concerning, the combination of strawberries-breast­ milk promoted the smallest total SCFA flux, but when strawberries were combined with blackcurrant and breastmilk, it resulted in the second highest increase in total SCFA flux (Fig. 5). Nevertheless, the influence of pumpkin and blackcurrant on the produc­ tion of propionate and acetate, respectively, was maintained when combining those foods with other foods and breastmilk. Although the combinations of FIG 4 Heatmap of multiple food-breastmilk combinations with the greatest influence on predicted fluxes of SCFAs and BCFAs. Fluxes of organic acids are expressed in relative variation in comparison to breastmilk alone. Cells are colored according to intensity, with the highest values in green and the lowest values in red. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 6 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 couscous-pumpkin-breastmilk and pumpkin-split peas-breastmilk did not change total SCFA flux, they resulted in the greatest increase in propionate. Similarly, the combi­ nations of blackcurrant-soybean-breastmilk, blackcurrant-strawberries-breastmilk, and blackcurrant-sweet potato-breastmilk had the highest increases in acetate flux. As expected, high-abundance genera, like Bifidobacterium, Bacteroides, and Bacillus, grew faster, while low-abundance genera, such as Lacticaseibacillus and Streptococcus, grew slower under multiple food-breastmilk combinations (Table S4). Negligible growth (growth rates near zero) for the low-abundance genera was also observed. For instance, 55 out of 89 multiple-food breastmilk combinations did not promote the growth of at least one genus included in the simulations, indicating that the trade-off between maximal community growth and individual growth was not respected, probably due to numerical issues. Overall, multiple food-breastmilk combinations tended toward decreasing the growth of the Bifidobacterium and Lactobacillus genera. Among the combinations with the greatest influence on the production of total SCFAs and BCFAs, blackcurrant-soybean- breastmilk resulted in the greatest increase in the growth rates of the Bacteroides genus (13.8%, Fig. 6). Blackcurrant-sweet potato-breastmilk promoted the greatest relative increases in the growth rates of the genera Prevotella (36.4%), Bacillus (28.5%), and Lactobacillus (17.0%), while black beans-blackcurrant-breastmilk the greatest relative decreases for these same genera (78.9, 76.5, and 76.7%, respectively). On the other hand, the combinations of chickpea-split peas-breastmilk, black beans-blackcurrant-breastmilk, and black beans-pumpkin-breastmilk decreased the growth of all high-abundance genera. DISCUSSION In this study, a metagenome-scale metabolic modeling approach was employed to identify foods and food combinations with the greatest impact on total SCFA and BCFA production by the colonic microbiota of New Zealand weaning infants. This is the first in silico screening of the influence of complementary foods on the colonic microbes of weaning infants. Another originality of this study was the in silico evaluation of foods combined with breastmilk and other foods, allowing a better representation of how complementary foods are introduced to weaning infants. FIG 5 Comparison of SCFA and BCFA fluxes when foods are combined individually with breastmilk or with other foods and breastmilk. Food-breastmilk combinations are ordered according to their predicted total SCFA flux, in which combinations resulting in higher fluxes are shown on the left. Fluxes for breastmilk alone (control) are shown on the right. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 7 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 A total of 155 food-breastmilk combinations (89 single food-breastmilk and 66 multiple-food breastmilk combinations) were investigated, and the 12 individual foods and 12 food combinations with the strongest influence on predicted fluxes of total SCFAs and BCFAs, when combined with breastmilk, were identified. Consistent with nutritional recommendations for weaning infants (41), foods included in the simulations were from diverse food groups, including legumes, fruits, vegetables, and animal proteins, while sweets and fast-food products were excluded. For instance, legumes and meats are recommended for infants as sources of protein and iron, while fresh fruits and vegetables are sources of certain vitamins, other minerals, and dietary fiber (39, 41, 42). MICOM assumes a steady state and employs a two-step linearization strategy that predicts growth rates based on the relative abundance of individual taxa (32). Further­ more, the food-breastmilk combinations consisted of 85% breastmilk by caloric intake, resulting in a similar nutritional profile. As a result, the predicted profile of faster-growing microbes was expected to remain consistent across different food-breastmilk combina­ tions, with high-abundance genera having higher growth rates than low-abundance genera. As expected, the highest growth rates were observed for highly abundant genera, while less abundant genera had slower growth. However, predicted microbial growth was not homogeneous under different dietary conditions, and growth rates near zero were also predicted for the less abundant genera Limosilactobacillus, Lacticaseibacil­ lus, Streptococcus, and Veillonella. This indicates that the solver failed to optimally solve the linear problem, most likely affecting outcomes of predicted growth rates and fluxes of organic acids. FIG 6 Heatmap of microbial growth rates for multiple food-breastmilk combinations with the greatest influence on predicted fluxes of SCFAs and BCFAs. Microbial growth rates are expressed in relative variation in comparison to breastmilk alone. Cells are colored according to intensity, with the highest values in green and the lowest values in red. Genera with predicted negligible growth are not presented. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 8 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 Compared to breastmilk (control), most of the food-breastmilk combinations decreased the growth rate of the Bifidobacterium and Collinsella genera, agreeing with findings from longitudinal studies showing the positive association between breastmilk consumption and the relative abundance of these genera (2, 17, 20). On the other hand, only small increases (up to 3.9%) in the growth rate of the genus Bacteroides, which are associated with complex carbohydrates and protein degradation, were predicted in silico. This contrasts with the observed increased fecal relative abundance of the Bacteroides genus when solid foods are introduced to infants (3, 43). Nevertheless, our ability to compare in silico predictions of microbial growth rates with in vitro or in vivo observed microbial abundance is limited. It is a technical challenge to accurately measure the growth rates of individual colonic microbes experimentally (44) and available data are rare, particularly for under-investigated groups like wean­ ing infants. Furthermore, changes in microbiota composition are not as informative as changes in microbiota function, considering that healthy individuals tend to have distinct microbiota composition but with similar functionality (9, 11). Thus, the focus of the study was on how food-breastmilk combinations affect the microbial production of health-relevant organic acids, particularly SCFAs and BCFAs, in the colonic microbiota of infants, rather than its microbial composition. As compared to breastmilk alone, pumpkin, raspberries, blackcurrant, soybeans, and sweet potato individually combined with breastmilk resulted in the greatest increases in the flux of total SCFAs. Among these foods, pumpkin and blackcurrant also stood out for promoting the greatest increase in the production of propionate and acetate, respectively. These findings suggest beneficial alterations to the colonic microbiome of infants at weaning, as SCFAs are associated with health benefits (45–48). SCFAs are primarily produced in the colon through the microbial fermentation of non-absorbed carbohydrates like dietary fiber and resistant starch. Acetate is the most abundant SCFA, especially during the first months of life, having anti-inflamma- tory properties and contributing to epithelial barrier integrity (49, 50). As the colonic microbiota matures during the weaning period, the proportion of propionate and butyrate in the colon increases (49). These organic acids, in addition to supporting colonic barrier integrity, provide various other health benefits, such as promoting metabolic and neurological health, regulating appetite, and preventing colon cancer (51–56). Butyrate also serves as an energy source for colon cells (57). The in silico predicted increase in SCFA production can be attributed to the high content of dietary fiber and phytochemicals in these foods, as suggested by other microbiome investigations in vitro, in animals, and clinical trials. For instance, fecal microbes of adults highly fermented pumpkin skin in vitro, resulting in increased production of propionate and total SCFAs (58), while adding pumpkin polyphenols into the high-fat diet of type 2 diabetic rats increased the colonic concentration of butyrate and total SCFAs in digesta (59). Soybean, which contains insoluble dietary fiber, increased the fecal content of acetate, propionate, and butyrate in mice on a high-fat diet (60), while sweet potato, which contains non-starch polysaccharides, increased the fecal content of acetate, propionate, and isobutyrate in diarrheic mice (61). Fecal fermentations of sweet potatoes using inoculum from healthy adults increased the production of acetate and total SCFAs due to the presence of fibers (62, 63), and also increased acetate production, likely due to anthocyanins found in the purple variety (64). Berries contain bioactive polyphenols and flavonoids, which influence the com­ position and function of colonic microbes. For example, blackcurrant anthocyanins decreased the ratio of Bacillota/Bacteroidota phyla in murine models (65) and increased caecal and serum concentrations of propionate and butyrate (66, 67). The fecal fermentation of different raspberry compounds, including phenolic extract and total dietary fiber, showed that SCFA production was mainly driven by polyphenol content and to a lesser extent by fiber content (68). Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-24 9 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 When combined with breastmilk, couscous resulted in the greatest increase in butyrate production, although decreased total SCFA flux. This observation may be justified by the presence of dietary fibers in durum wheat, in particular arabinoxylan (69). Interventions on adults consuming wheat arabinoxylan reported increased fecal butyrate concentration (70, 71). Similarly, the fecal fermentation of wheat cereal using inoculum from six healthy weaning infants increased the production of butyrate, although the extent of the increase varied among individuals (29). Couscous, pumpkin, strawberries, and split peas also reduced total BCFA flux when combined with breastmilk. BCFAs, like isobutyrate and isovalerate, are produced in the colon through microbial fermentation of non-absorbed amino acids. Less is known about their influence on host health but evidence suggests a link with metabolic functions (72). BCFAs are biomarkers for protein fermentation in the distal colon (73), which can generate potential deleterious metabolites, such as ammonia and phenols (74). Colonic microbial production of BCFAs was negatively correlated with dietary insoluble fiber intake (75), which could partially explain the decreased BCFAs production observed in silico when adding the above fiber-rich foods to breastmilk. For instance, an intervention with yellow pea fiber decreased the fecal concentration of the BCFA isovalerate in overweight adults (76). Nevertheless, an in vitro fecal fermen­ tation of 22 plant sources of fiber, including whole cereals, seeds, and pulses, using inoculum from healthy adults, found no changes in the production of BCFAs (77). In this in silico study, other fiber-rich foods like sweet potato, black beans, and chickpeas also did not change the flux of BCFAs when combined with breastmilk, suggesting that other dietary compounds affect BCFA production. Indeed, a recent murine model study demonstrated that the protein source is a key factor affecting the fecal BCFA content (78). On the other hand, when strawberries and shrimp were individually combined with breastmilk, they, respectively, showed the greatest decrease in flux of total SCFAs and the greatest increase in total BCFA flux in silico. These observations suggest a potential deleterious alteration in the function of colonic microbes. Evidence demonstrated that reduced SCFA fecal concentration in infancy is associated with an increased risk of diseases and allergies (16, 79, 80) while increased content of BCFAs in infants’ feces was linked to weight gain (81). These imbalances in organic acid production by the colonic microbiota during infancy may affect later life. In adults, decreased fecal concentration of SCFAs has been linked with colonic dysbiosis and diseases, such as encephalitis, Parkinson, and type 2 diabetes (46, 82, 83). Currently, there is a lack of studies evaluating how the consumption of shrimp affects colonic microbes, while evidence found for strawberries contrasted with the in silico observations. Strawberry interventions did not change the fecal SCFA content in healthy adults (84) but increased the caecal production of acetate, propionate, and butyrate in mice with colitis (85). Importantly, combining foods with other foods and with breastmilk shifted their influence on colonic microbes. The combination of pumpkin-breastmilk promoted the highest flux of total SCFAs but, when pumpkin was combined with other foods and breastmilk, resulted in little or no change in total SCFA fluxes. Among the multi­ ple food-breastmilk combinations, blackcurrant stood out for promoting the greatest increases in SCFA production. This observation suggests that the interaction between dietary compounds composing a meal has a stronger influence on the metabolism of colonic microbes than individual foods. As foods are rarely consumed individually, dietary patterns rather than individual foods are more likely to promote alterations in colonic microbes that may affect host health (86, 87). A recent in vitro study evaluated how 32 foods are fermented by the fecal micro­ biota of New Zealand infants at weaning age. After 24 h of fermentation, blackcurrant, pumpkin, and sweet potato increased total SCFA production, which is consistent with the in silico observations reported here (28). Similarly, pumpkin and blackcurrant increased the production of acetate after 24 h of fermentation using fecal inoculum from New Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2410 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 Zealand weaning infants (27). However, none of these in vitro studies combined foods with human milk or with other foods and human milk, providing limited information about how foods may influence the colonic microbiota of infants when added to their pre-existing dietary patterns. In vivo approaches can better evaluate the relationship between dietary patterns and the colonic microbiota, measuring eventual host health outcomes. For example, a series of studies investigated the colonic microbes of malnourished infants using animal models to identify food combinations that promoted the growth of bacterial taxa associated with healthy colonic microbial development during weaning (30, 88, 89). Observational trials evaluating dietary patterns in weaning infants associated increased dietary diversity with increased microbial diversity and consequently the stabilization of the colonic microbiome (90), and with increased SCFA production (3), highlighting the importance of introducing infants to meals composed of diverse foods. Nevertheless, clinical trials evaluating the relationship between diet and colonic microbiome are resource- and time-consuming, also facing other limitations, such as ethical concerns, low patient recruitment and adherence to the intervention, and poor accuracy of food questionnaires. In this scenario, modeling is a useful complementary tool to evaluate conditions that cannot be investigated with traditional in vitro and in vivo methods due to technical or logistical limitations. However, this in silico approach has limitations. To simulate the impact of foods on the colonic microbiota of weaning infants, average fecal microbiota data were used rather than data at the individual level. Repeating the simulations with the same dietary conditions produced similar outcomes. Consequently, it was impossible to conduct statistical analyses to assess potential differences in the effect of various food-breastmilk combinations on microbial growth and metabolite production. To reduce computational demands, a threshold of 1% of relative abundance was used to select the microbial genera included in the simulations. This resulted in a small-scale microbial community that excluded genera usually found in low abundance in the colon of weaning infants but that have key metabolic functions, such as the butyrate producers Clostridium and Faecalibacterium. The design of food-breastmilk combinations was limited by the accuracy and availability of data in food composition databases, which predominantly focus on Western-type foods, lacking information on the diversity of food varieties and cook­ ing/preparation methods (for a review on the limitations of food composition databases, see reference (91)). This study prioritized foods produced or available in New Zealand. However, the Virtual Metabolic Human database (92) is based on the USDA National Nutrient Database for Standard Reference Release 28 (93) and does not cover all foods consumed by weaning infants in New Zealand. For example, traditional indigenous foods like kūmara (sweet potato variety) are not covered in the database and had to be replaced by the most similar food available in the database. Furthermore, data available in the VMH database do not account for individual variability in breastmilk composition. Our study did not consider key factors influencing the composition of the colonic microbiota in infants during the first months of life, such as the mode of delivery, breastfeeding versus infant formula feeding, and the mother’s diet (94, 95). Since our simulations focused on the weaning period (typically between 5 and 12 months of age), we assumed that these early-life factors would have a lesser impact on the infant microbiome compared to dietary changes introduced during weaning (consumption of solids and less intake of breastmilk or infant formula). Nonetheless, metagenome-scale metabolic models create personalized simulations based on input data. This methodol­ ogy is adaptable and can be applied to other infant populations, potentially accounting for varying early-life influences on the microbiome in future investigations. Another drawback is that the accuracy of the simulations presented here is depend­ ent on the quality of the microbial metabolic reconstructions (a mathematical repre­ sentation of the biochemical reactions that a microorganism can perform), which may contain missing information (31) and greatly affect the prediction capability of Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2411 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 genome-scale metabolic models (96). Microbial relative abundance data were obtained from 16S rRNA sequencing of fecal samples. This is the most common method to assess colonic microbes of weaning infants but 16S rRNA sequencing provides resolution suitable only to the genus rank (97), while the use of feces does not accurately represent the ratios of microbial communities found in the colonic mucosa and in the proximal colon (98). Relative abundances are not independent by nature, meaning that the relative abundance of one microbe is affected by the abundance of any other microbe in the community. As a result, this may lead to a false representation of the actual microbial community. Although the metabolic reconstructions of colonic microbes used in the simulations were validated using independent data sets (99), our simulations using MICOM were not validated experimentally (due to resource and ethical constraints) or using an external data set (due to a lack of available data). A recent study using MICOM reported agree­ ment between fluxes of propionate and butyrate predicted in silico and those estimated in vitro (100). The package stood out among seven other static metabolic modeling tools in a qualitative and quantitative assessment (101) but did not accurately correlate growth rates predicted in silico with those observed in vitro (96). Furthermore, numerical instabilities are common in MICOM when solving quadratic programming problems in large community models (32). Near zero growth rates were observed in our simulations for the less abundant genera, indicating instabilities that are likely to alter predicted fluxes of SCFAs and BCFAs. The lack of research evaluating the effect of food-breastmilk combinations on the colonic microbiome of weaning infants and the differences between predicted in silico outcomes (microbial growth rates and fluxes of metabolites) and experimental outcomes (microbial relative abundance and concentration of metabolites) limits even a qualitative comparison between in silico observations reported here and published in vitro or in vivo results. To calculate the fluxes of microbial metabolites through flux balance analysis, a mass steady state was assumed, implying that there is no accumulation of substrates in the intracellular space (see review (102)). However, this assumption strongly differs from in vitro static conditions, where resources are depleted, and microbial products accumulate over time. Due to this assumption, parameters used in the simulation, such as substrate fluxes and microbial growth rates, are fixed and thus may not represent the rapid changes of colonic microbes in response to diet (4, 103). Dynamic flux balance analysis may be a promising alternative to better represent in vivo conditions (104). However, this strategy is computationally intensive when dealing with complex microbial communities (105), hence limiting its use to screen the effect of a broad range of food combinations on colonic microbes. Finally, in silico predictions may diverge from the behavior of colonic microbes observed in vitro or in vivo. Both unabsorbed carbohydrates and amino acids contribute to colonic microbial production of SCFAs, but colonic microbes preferentially ferment carbohydrates (106, 107). Acetate is produced in larger quantities by most colonic commensals, while propionate and butyrate are produced in lesser amounts by only a few genera normally through cross-feeding interactions (108, 109). Reduced dietary carbohydrate intake, even if replaced with protein, ultimately decreases the production of butyrate (110). Thus, one could expect that fiber-rich plant foods would increase the in silico production of SCFAs when combined with breastmilk. However, that was not the case for all food combinations here evaluated. This result may be justified by the highly individualized and varied response of colonic microbes to dietary compounds (111), particularly dietary fibers like resistant starch (112, 113), but also impacted by the limitations cited above. Therefore, further experimental investigation is essential to validate the foods and food combinations identified in this study. Nevertheless, data generated by this research provide a direction for future food-microbiome investigations Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2412 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 and show the potential of modeling approaches to complement in vitro and in vivo techniques. Currently, there is a lack of knowledge about how solid foods affect the colonic microbiome of weaning infants. This study evaluated for the first time how food-breast­ milk combinations affect the colonic microbial production of SCFAs and BCFAs using a metagenome-scale metabolic modeling approach. By quickly and inexpensively generating insights in silico, this study helps the design of future research in vitro and in vivo, contributing to fill a crucial knowledge gap in infant nutrition. Furthermore, our in silico observations suggest that the interaction of foods composing a meal has a key influence on the colonic microbial production of SCFAs and BCFAs. This encourages future microbiome investigations to focus on the combined effect of foods on colonic microbes instead of focusing on the effect of individual food items. MATERIALS AND METHODS Software Simulations were performed in Python (version 3.9) using the package MICOM(32) (version 0.32.5) and the integrated development environment Spyder (version 5). The solver CPLEX Optimisation Studio (IBM ILOG, version 22.1) was employed under an academic license. Modeling workflow Simulations used the relative abundance of the genera composing the fecal microbiota of weaning infants, metabolic reconstructions for each taxon (a list of the biochemi­ cal reactions performed by a microorganism), and the fluxes of dietary compounds of food-breastmilk combinations (a list of nutrients composing a diet, in which their respective quantities are expressed in units of time) as inputs. Outputs were the predicted individual microbial growth rates and the fluxes of metabolites produced by the microbiota, particularly SCFAs and BCFAs. The average relative abundance of the fecal microbiota of 14 New Zealand infants aged between 5 and 12 months, obtained from a previous in vitro study (28), was used to represent the colonic microbiota of weaning infants (while minimizing variations in microbiota composition between individuals). Raw 16S rRNA sequencing data (Illumina MiSeq, 2 × 250 bp paired-end reads) for the fermentation control (water, at t = 0 h) were imported from NCBI archives (BioProject PRJNA669972) using the plugin q2-fondue (114). A total of 241,611 paired-end demultiplexed reads (~80,000 sequences/sample) were imported as “.qza” artefact in QIIME2 (115) (version 2023.2), resulting in a quality score of approximately 37. DADA2 (116) was used to denoise and filter the reads, which were trimmed at 250 bp. Taxonomy assignment was carried out using the q2-feature- classifier plugin and the Greengenes2 database (117), in which amplicon sequence variants were collapsed at the genus level (see Table S5). To reduce the numerical instability (a failure in optimally solving a numerical problem) and the processing time of the simulations, taxa were filtered to include only genera with at least 1% relative abundance. As proposed by MICOM’s authors (32), low-abundance taxa can be discarded as they are unlikely to affect the overall production of microbial metabolites but increase the computational cost of the simulations. Eleven genera remained, accounting for 95% of the relative abundance of the microbial community. The Assembly of Gut Organisms through Reconstruction and Analysis version 2 (AGORA2) (99) metabolic reconstructions were employed in the simulations. AGORA2 reconstructions were available for 10 of these genera, representing 94% of the relative abundance used as input in the simulations. The genus rank was chosen because it suited both the resolution of 16S rRNA sequencing (97) and pan models of AGORA2 metabolic reconstructions (which are available on the MICOM’s GitHub page). Pan models of the AGORA2 metabolic reconstructions at the genus rank were built by Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2413 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 pooling individual metabolic reconstructions for strains of colonic microbes into higher taxonomic ranks. Fluxes of dietary compounds were generated according to a pre-defined workflow, which is available on the MICOM GitHub page. Foods included in the simulations were identified through a literature survey. Those corresponded to food items that are commonly introduced to infants at weaning (6 to 12 months old) in New Zealand (38, 40). To account for cultural variations in dietary habits, complementary foods introduced to infants in other geographical locations such as North America (37) and Europe (36) were also included. In alignment with current nutrition recommendations (39, 41), included food items were from various food groups, such as vegetables, fruits, nuts, legumes, meats, cereals, and dairy products, while food groups not recommended to infants, such as sweets and fast-food products, were excluded. To represent how solid foods are consumed by weaning infants, identified food items were combined, individually or in pairs, with mature breastmilk using the “Design a diet” function of the Virtual Metabolic Human database (92) (VMH) (see Table S6 for the list of foods). Food varieties are prioritized for New Zealand production or local availability. When New Zealand varieties were not included in the VMH database (which is based on the USDA National Nutrient Database), they were replaced with similar available varieties. Previous evidence estimated that at 6 months of age, the average weight of infants is 7.6 kg (118) and their energy requirement is 85 kcal/kg/day (119), of which 85% is supplied by breastmilk (119). Therefore, it was assumed that diets for 6-month-old infants would consist of 85% of breastmilk and 15% of complementary foods, totaling 608 kcal/day. In total, dietary fluxes of 155 food-breastmilk combinations and two controls (infant formula only and breastmilk only) were imported from the VMH database and processed in their original unit, mmol/day, as this reduced the numerical instability of the simulations. Host-secreted mucin cores and the bile acids glycocholic acid and taurocholic acid were then added to the dietary fluxes (values of 1 mmol/h) to better represent the substrates available for microbes in the human colon. The next step of the MICOM workflow used the human metabolic network Recon3D (120) resource to identify the dietary compounds that are absorbed in the small intestine. The flux of these com­ pounds was then multiplied by a factor of 0.2 to account for their intestinal absorption. Finally, AGORA2 (99) metabolic reconstructions were used to identify and add the missing nutrients necessary to allow growth rates of at least 0.01 /h for all microbial species of colonic microbes comprised by these metabolic reconstructions. This step is necessary because the dietary fluxes provided by food databases often lack essential cofactors for microbial growth, such as vitamins and minerals. After supplementing the fluxes with the minimal additional substrates (which are added in mmol/h), the final fluxes were in mmol/h. Modelling assumptions This modeling approach was based on flux balance analysis under a mass steady-state assumption, meaning that there is no accumulation of substances inside the microbial cells (102). A constrained linear programming problem (a mathematical problem) for the fluxes (v) was created using a stoichiometric matrix (S), which is a matrix representing all biochemical reactions (as columns) and involved metabolites (as rows) performed by the microbial community. The objective of the flux balance analysis was to maximize biomass reaction (vbm) such that there is no accumulation of substrates in the system (S*v = 0). The steady-state assumption represents the exponential phase of bacterial growth, in which growth rates can be assumed constant, and the linear problem describing fluxes of metabolites can be solved by a programming solver. The community growth rate (µc) is determined by the sum of growth rates of individual microbes (µi) weighted by their relative abundances (ai), as described by the following equation: Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2414 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 μc = i aiμi MICOM also distributes growth across all members of the microbial community by limiting the maximal community growth with a cooperative trade-off (a value between 0 and 1). The optimal trade-off between maximal microbial community growth and maximal individual microbial growth was determined for each food-breastmilk combination. This allows most microbes present in the community to grow, rather than only the dominant members, thus better representing in vivo conditions in which multiple microbes are found in different abundances in the human colon (11). The cooperative trade-off (α) was calculated in two steps, initially determining the maximal community growth rate (µc max) and then calculating the optimal individual growth rates (µi opt), which are positively correlated with the relative abundance (a) of individual taxa (32): μiopt = α μcmaxaT a  a MICOM assumes that the relative abundance of microbes present in a sample corresponds to their relative biomass, in which the biomass reaction is normalized to produce 1 g, predicting the production of microbial metabolites in mmoles/g.h (of microbial biomass in dry weight) (for further description see reference (32)). The total flux of a given microbial metabolite (vtot m) is calculated as the sum of the fluxes of this metabolite produced by individual microbes (vi m) weighted by their relative abundances (ai), as described by the equation below: vtotm = i ai vim Criteria for selecting food-breastmilk combinations To identify food and food combinations having the strongest influence on the colonic microbiome of New Zealand weaning infants, this study focused on the production of microbial metabolites (microbial function) rather than changes in microbial composition. This is due to the potential of functional analyses to be more informative about how diets shape colonic microbes. Indeed, evidence linked variations in the genetic content of colonic commensals between individuals with their ability to metabolize dietary compounds and produce bioactive metabolites (121, 122). Predicted microbial fluxes of the major SCFAs (acetate, propionate, and butyrate) and BCFAs (isovalerate and isobutyrate) were chosen due to the importance of these organic acids on host physiology and because their concentration is known to be influenced by host dietary patterns during weaning (49, 123). As criteria of choice, selected foods and food combinations were representative of different model outputs, corresponding to candidates that, when combined with breastmilk, resulted in the greatest increase, decrease, or did not change the fluxes of total SCFAs and BCFAs, in comparison to breastmilk alone. Foods commonly consumed by weaning infants in New Zealand and produced in the country were prioritized. ACKNOWLEDGMENTS V.G.S. is supported by a Riddet Institute PhD Scholarship. PhD stipend and the project are funded by the High-Value Nutrition National Science Challenge. This study was funded by the New Zealand Ministry for Business, Innovation and Employment (MBIE, grant no. 3710040) through the High-Value Nutrition National Science Challenge. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2415 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 AUTHOR AFFILIATIONS 1Riddet Institute, Massey University, Palmerston North, New Zealand 2High-Value Nutrition National Science Challenge, Auckland, New Zealand 3AgResearch, Palmerston North, New Zealand 4Department of Nutrition and Dietetics, The University of Auckland, Auckland, New Zealand 5Department of Human Nutrition, University of Otago, Dunedin, New Zealand AUTHOR ORCIDs Vitor G. da Silva http://orcid.org/0000-0001-9218-4473 Warren C. McNabb http://orcid.org/0000-0003-2514-6551 FUNDING Funder Grant(s) Author(s) High-Value Nutrition National Science Challenge Vitor G. da Silva Jane A. Mullaney Clare Wall Nicole C. Roy Warren C. McNabb ADDITIONAL FILES The following material is available online. Supplemental Material Supplemental tables (mSystems00577-24-s0001.pdf). Tables S1 to S6. REFERENCES 1. Sender R, Fuchs S, Milo R. 2016. Revised estimates for the number of human and bacteria cells in the body. PLoS Biol 14:e1002533. https:// doi.org/10.1371/journal.pbio.1002533 2. Laursen MF, Andersen LBB, Michaelsen KF, Mølgaard C, Trolle E, Bahl MI, Licht TR. 2016. Infant gut microbiota development is driven by transition to family foods independent of maternal obesity. mSphere 1:e00069-15. https://doi.org/10.1128/mSphere.00069-15 3. Gómez-Martín M, Saturio S, Arboleya S, Herrero-Morín D, Calzón M, López T, González S, Gueimonde M. 2022. Association between diet and fecal microbiota along the first year of life. Food Res Int 162:111994. https://doi.org/10.1016/j.foodres.2022.111994 4. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, Ling AV, Devlin AS, Varma Y, Fischbach MA, Biddinger SB, Dutton RJ, Turnbaugh PJ. 2014. Diet rapidly and reproducibly alters the human gut microbiome. Nature New Biol 505:559–563. https://doi.org/10.1038/ nature12820 5. Manor O, Dai CL, Kornilov SA, Smith B, Price ND, Lovejoy JC, Gibbons SM, Magis AT. 2020. Health and disease markers correlate with gut microbiome composition across thousands of people. Nat Commun 11:5206. https://doi.org/10.1038/s41467-020-18871-1 6. Ghosh S, Whitley CS, Haribabu B, Jala VR. 2021. Regulation of intestinal barrier function by microbial metabolites. Cell Mol Gastroenterol Hepatol 11:1463–1482. https://doi.org/10.1016/j.jcmgh.2021.02.007 7. Frank DN, St. Amand AL, Feldman RA, Boedeker EC, Harpaz N, Pace NR. 2007. Molecular-phylogenetic characterization of microbial community imbalances in human inflammatory bowel diseases. Proc Natl Acad Sci USA 104:13780–13785. https://doi.org/10.1073/pnas.0706625104 8. Zhao L, Lou H, Peng Y, Chen S, Fan L, Li X. 2020. Elevated levels of circulating short-chain fatty acids and bile acids in type 2 diabetes are linked to gut barrier disruption and disordered gut microbiota. Diabetes Res Clin Pract 169:108418. https://doi.org/10.1016/j.diabres. 2020.108418 9. Olsson LM, Boulund F, Nilsson S, Khan MT, Gummesson A, Fagerberg L, Engstrand L, Perkins R, Uhlén M, Bergström G, Tremaroli V, Bäckhed F. 2022. Dynamics of the normal gut microbiota: a longitudinal one-year population study in Sweden. Cell Host Microbe 30:726–739. https://doi. org/10.1016/j.chom.2022.03.002 10. Rinninella E, Raoul P, Cintoni M, Franceschi F, Miggiano GAD, Gasbarrini A, Mele MC. 2019. What is the healthy gut microbiota composition? A changing ecosystem across age, environment, diet, and diseases. Microorganisms 7:14. https://doi.org/10.3390/microorganisms7010014 11. Human Microbiome Project Consortium. 2012. Structure, function and diversity of the healthy human microbiome. Nat New Biol 486:207–214. https://doi.org/10.1038/nature11234 12. Tian L, Wang X-W, Wu A-K, Fan Y, Friedman J, Dahlin A, Waldor MK, Weinstock GM, Weiss ST, Liu Y-Y. 2020. Deciphering functional redundancy in the human microbiome. Nat Commun 11:6217. https:// doi.org/10.1038/s41467-020-19940-1 13. Cox LM, Yamanishi S, Sohn J, Alekseyenko AV, Leung JM, Cho I, Kim SG, Li H, Gao Z, Mahana D, Zárate Rodriguez JG, Rogers AB, Robine N, Loke P, Blaser MJ. 2014. Altering the intestinal microbiota during a critical developmental window has lasting metabolic consequences. Cell 158:705–721. https://doi.org/10.1016/j.cell.2014.05.052 14. Cho I, Yamanishi S, Cox L, Methé BA, Zavadil J, Li K, Gao Z, Mahana D, Raju K, Teitler I, Li H, Alekseyenko AV, Blaser MJ. 2012. Antibiotics in early life alter the murine colonic microbiome and adiposity. Nature New Biol 488:621–626. https://doi.org/10.1038/nature11400 15. Miller SA, Wu RKS, Oremus M. 2018. The association between antibiotic use in infancy and childhood overweight or obesity: a systematic review and meta‐analysis. Obes Rev 19:1463–1475. https://doi.org/10. 1111/obr.12717 Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2416 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1128/msystems.00577-24 https://doi.org/10.1371/journal.pbio.1002533 https://doi.org/10.1128/mSphere.00069-15 https://doi.org/10.1016/j.foodres.2022.111994 https://doi.org/10.1038/nature12820 https://doi.org/10.1038/s41467-020-18871-1 https://doi.org/10.1016/j.jcmgh.2021.02.007 https://doi.org/10.1073/pnas.0706625104 https://doi.org/10.1016/j.diabres.2020.108418 https://doi.org/10.1016/j.chom.2022.03.002 https://doi.org/10.3390/microorganisms7010014 https://doi.org/10.1038/nature11234 https://doi.org/10.1038/s41467-020-19940-1 https://doi.org/10.1016/j.cell.2014.05.052 https://doi.org/10.1038/nature11400 https://doi.org/10.1111/obr.12717 https://doi.org/10.1128/msystems.00577-24 16. Arrieta M-C, Stiemsma LT, Dimitriu PA, Thorson L, Russell S, Yurist- Doutsch S, Kuzeljevic B, Gold MJ, Britton HM, Lefebvre DL, Subbarao P, Mandhane P, Becker A, McNagny KM, Sears MR, Kollmann T, Investigators CS, Mohn WW, Turvey SE, Finlay BB. 2015. Early infancy microbial and metabolic alterations affect risk of childhood asthma. Sci Transl Med 7:307ra152. https://doi.org/10.1126/scitranslmed.aab2271 17. Stewart CJ, Ajami NJ, O’Brien JL, Hutchinson DS, Smith DP, Wong MC, Ross MC, Lloyd RE, Doddapaneni H, Metcalf GA, et al. 2018. Temporal development of the gut microbiome in early childhood from the TEDDY study. Nat New Biol 562:583–588. https://doi.org/10.1038/ s41586-018-0617-x 18. Thompson AL, Monteagudo-Mera A, Cadenas MB, Lampl ML, Azcarate- Peril MA. 2015. Milk- and solid-feeding practices and daycare attendance are associated with differences in bacterial diversity, predominant communities, and metabolic and immune function of the infant gut microbiome. Front Cell Infect Microbiol 5:3. https://doi.org/ 10.3389/fcimb.2015.00003 19. Qasem W, Azad MB, Hossain Z, Azad E, Jorgensen S, Castillo San Juan S, Cai C, Khafipour E, Beta T, Roberts LJ, Friel J. 2017. Assessment of complementary feeding of Canadian infants: effects on microbiome & oxidative stress, a randomized controlled trial. BMC Pediatr 17:54. https: //doi.org/10.1186/s12887-017-0805-0 20. Oyedemi OT, Shaw S, Martin JC, Ayeni FA, Scott KP. 2022. Changes in the gut microbiota of Nigerian infants within the first year of life. PLoS One 17:e0265123. https://doi.org/10.1371/journal.pone.0265123 21. Ordiz MI, Janssen S, Humphrey G, Ackermann G, Stephenson K, Agapova S, Divala O, Kaimila Y, Maleta K, Zhong C, Knight R, Trehan I, Tarr PI, Rusconi B, Manary MJ. 2020. The effect of legume supplementa­ tion on the gut microbiota in rural Malawian infants aged 6 to 12 months. Am J Clin Nutr 111:884–892. https://doi.org/10.1093/ajcn/ nqaa011 22. Bierut T, Duckworth L, Grabowsky M, Ordiz MI, Laury ML, Callaghan- Gillespie M, Maleta K, Manary MJ. 2021. The effect of bovine colostrum/egg supplementation compared with corn/soy flour in young Malawian children: a randomized, controlled clinical trial. Am J Clin Nutr 113:420–427. https://doi.org/10.1093/ajcn/nqaa325 23. Plaza-Diaz J, Bernal MJ, Schutte S, Chenoll E, Genovés S, Codoñer FM, Gil A, Sanchez-Siles LM. 2021. Effects of whole-grain and sugar content in infant cereals on gut microbiota at weaning: a randomized trial. Nutrients 13:1496. https://doi.org/10.3390/nu13051496 24. Krebs NF, Sherlock LG, Westcott J, Culbertson D, Hambidge KM, Feazel LM, Robertson CE, Frank DN. 2013. Effects of different complementary feeding regimens on iron status and enteric microbiota in breastfed infants. J Pediatr 163:416–423. https://doi.org/10.1016/j.jpeds.2013.01. 024 25. Zambrana LE, McKeen S, Ibrahim H, Zarei I, Borresen EC, Doumbia L, Boré A, Cissoko A, Douyon S, Koné K, Perez J, Perez C, Hess A, Abdo Z, Sangaré L, Maiga A, Becker-Dreps S, Yuan L, Koita O, Vilchez S, Ryan EP. 2019. Rice bran supplementation modulates growth, microbiota and metabolome in weaning infants: a clinical trial in Nicaragua and Mali. Sci Rep 9:13919. https://doi.org/10.1038/s41598-019-50344-4 26. Tang M, Ma C, Weinheimer-Haus EM, Robertson CE, Kofonow JM, Berman LM, Waljee A, Zhu J, Frank DN, Krebs NF. 2022. Different gut microbiota in U.S. formula-fed infants consuming a meat vs. dairy- based complementary foods: a randomized controlled trial. Front Nutr 9:1063518. https://doi.org/10.3389/fnut.2022.1063518 27. Parkar SG, Frost JKT, Rosendale D, Stoklosinski HM, Jobsis CMH, Hedderley DI, Gopal P. 2021. The sugar composition of the fibre in selected plant foods modulates weaning infants’ gut microbiome composition and fermentation metabolites in vitro. Sci Rep 11:9292. https://doi.org/10.1038/s41598-021-88445-8 28. Parkar SG, Rosendale DI, Stoklosinski HM, Jobsis CMH, Hedderley DI, Gopal P. 2021. Complementary food ingredients alter infant gut microbiome composition and metabolism in vitro. Microorganisms 9:2089. https://doi.org/10.3390/microorganisms9102089 29. Gamage H, Tetu SG, Chong RWW, Ashton J, Packer NH, Paulsen IT. 2017. Cereal products derived from wheat, sorghum, rice and oats alter the infant gut microbiota in vitro. Sci Rep 7:14312. https://doi.org/10.1038/ s41598-017-14707-z 30. Gehrig JL, Venkatesh S, Chang H-W, Hibberd MC, Kung VL, Cheng J, Chen RY, Subramanian S, Cowardin CA, Meier MF, et al. 2019. Effects of microbiota-directed foods in gnotobiotic animals and undernourished children. Science 365:eaau4732. https://doi.org/10.1126/science. aau4732 31. Thiele I, Palsson BØ. 2010. A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat Protoc 5:93–121. https:// doi.org/10.1038/nprot.2009.203 32. Diener C, Gibbons SM, Resendis-Antonio O. 2020. MICOM: metage­ nome-scale modeling to infer metabolic interactions in the gut microbiota. mSystems 5:e00606-19. https://doi.org/10.1128/mSystems. 00606-19 33. Henson MA. 2021. Interrogation of the perturbed gut microbiota in gouty arthritis patients through in silico metabolic modeling. Eng Life Sci 21:489–501. https://doi.org/10.1002/elsc.202100003 34. Hertel J, Harms AC, Heinken A, Baldini F, Thinnes CC, Glaab E, Vasco DA, Pietzner M, Stewart ID, Wareham NJ, Langenberg C, Trenkwalder C, Krüger R, Hankemeier T, Fleming RMT, Mollenhauer B, Thiele I. 2019. Integrated analyses of microbiome and longitudinal metabolome data reveal microbial-host interactions on sulfur metabolism in Parkinson’s disease. Cell Rep 29:1767–1777. https://doi.org/10.1016/j.celrep.2019. 10.035 35. Bauer E, Thiele I. 2018. From metagenomic data to personalized in silico microbiotas: predicting dietary supplements for Crohn's disease. NPJ Syst Biol Appl 4:27. https://doi.org/10.1038/s41540-018-0063-2 36. Betoko A, Charles M-A, Hankard R, Forhan A, Bonet M, Saurel-Cubizolles M-J, Heude B, de Lauzon-Guillain B, EDEN mother-child cohort study group. 2013. Infant feeding patterns over the first year of life: influence of family characteristics. Eur J Clin Nutr 67:631–637. https://doi.org/10. 1038/ejcn.2012.200 37. Bailey RL, Stang JS, Davis TA, Naimi TS, Schneeman BO, Dewey KG, Donovan SM, Novotny R, Kleinman RE, Taveras EM, Bazzano L, Snetselaar LG, de Jesus J, Casavale KO, Stoody EE, Goldman JD, Moshfegh AJ, Rhodes DG, Herrick KA, Koegel K, Perrine CG, Pannucci T. 2022. Dietary and complementary feeding practices of us infants, 6 to 12 months: a narrative review of the federal nutrition monitoring data. J Acad Nutr Diet 122:2337–2345. https://doi.org/10.1016/j.jand.2021.10. 017 38. Gontijo de Castro T, Gerritsen S, Santos LP, Marchioni DML, Morton SMB, Wall C. 2022. Child feeding indexes measuring adherence to New Zealand nutrition guidelines: development and assessment. Matern Child Nutr 18:e13402. https://doi.org/10.1111/mcn.13402 39. New Zealand Ministry of Health. 2021. Healthy eating guidelines for New Zealand babies and toddlers (0–2 years old). Wellington: Ministry of Health 40. New Zealand Ministry of Social Development. 2018. Infant feeding in New Zealand: adherence to food and nutrition guidelines among the growing up in New Zealand cohort. Ministry of Social Development 41. Alvisi P, Brusa S, Alboresi S, Amarri S, Bottau P, Cavagni G, Corradini B, Landi L, Loroni L, Marani M, Osti IM, Povesi-Dascola C, Caffarelli C, Valeriani L, Agostoni C. 2015. Recommendations on complementary feeding for healthy, full-term infants. Ital J Pediatr 41:36. https://doi.org/ 10.1186/s13052-015-0143-5 42. New Zealand Ministry of Health. 2008. Food and nutrition guidelines for healthy infants and toddlers (aged 0-2): background paper. Ministry of Health, Wellington. 43. Aguilar-Lopez M, Wetzel C, MacDonald A, Ho TTB, Donovan SM. 2021. Human milk-based or bovine milk-based fortifiers differentially impact the development of the gut microbiota of preterm infants. Front Pediatr 9:719096. https://doi.org/10.3389/fped.2021.719096 44. Koch BJ, McHugh TA, Hayer M, Schwartz E, Blazewicz SJ, Dijkstra P, van Gestel N, Marks JC, Mau RL, Morrissey EM, Pett‐Ridge J, Hungate BA. 2018. Estimating taxon‐specific population dynamics in diverse microbial communities. Ecosphere 9:e02090. https://doi.org/10.1002/ ecs2.2090 45. Friščić J, Dürholz K, Chen X, Engdahl C, Möller L, Schett G, Zaiss MM, Hoffmann MH. 2021. Dietary derived propionate regulates pathogenic fibroblast function and ameliorates experimental arthritis and inflammatory tissue priming. Nutrients 13:1643. https://doi.org/10. 3390/nu13051643 46. Sanna S, van Zuydam NR, Mahajan A, Kurilshikov A, Vich Vila A, Võsa U, Mujagic Z, Masclee AAM, Jonkers DMAE, Oosting M, Joosten LAB, Netea MG, Franke L, Zhernakova A, Fu J, Wijmenga C, McCarthy MI. 2019. Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2417 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1126/scitranslmed.aab2271 https://doi.org/10.1038/s41586-018-0617-x https://doi.org/10.3389/fcimb.2015.00003 https://doi.org/10.1186/s12887-017-0805-0 https://doi.org/10.1371/journal.pone.0265123 https://doi.org/10.1093/ajcn/nqaa011 https://doi.org/10.1093/ajcn/nqaa325 https://doi.org/10.3390/nu13051496 https://doi.org/10.1016/j.jpeds.2013.01.024 https://doi.org/10.1038/s41598-019-50344-4 https://doi.org/10.3389/fnut.2022.1063518 https://doi.org/10.1038/s41598-021-88445-8 https://doi.org/10.3390/microorganisms9102089 https://doi.org/10.1038/s41598-017-14707-z https://doi.org/10.1126/science.aau4732 https://doi.org/10.1038/nprot.2009.203 https://doi.org/10.1128/mSystems.00606-19 https://doi.org/10.1002/elsc.202100003 https://doi.org/10.1016/j.celrep.2019.10.035 https://doi.org/10.1038/s41540-018-0063-2 https://doi.org/10.1038/ejcn.2012.200 https://doi.org/10.1016/j.jand.2021.10.017 https://doi.org/10.1111/mcn.13402 https://doi.org/10.1186/s13052-015-0143-5 https://doi.org/10.3389/fped.2021.719096 https://doi.org/10.1002/ecs2.2090 https://doi.org/10.3390/nu13051643 https://doi.org/10.1128/msystems.00577-24 Causal relationships among the gut microbiome, short-chain fatty acids and metabolic diseases. Nat Genet 51:600–605. https://doi.org/10. 1038/s41588-019-0350-x 47. Kaye DM, Shihata WA, Jama HA, Tsyganov K, Ziemann M, Kiriazis H, Horlock D, Vijay A, Giam B, Vinh A, Johnson C, Fiedler A, Donner D, Snelson M, Coughlan MT, Phillips S, Du X-J, El-Osta A, Drummond G, Lambert GW, Spector TD, Valdes AM, Mackay CR, Marques FZ. 2020. Deficiency of prebiotic fiber and insufficient signaling through gut metabolite-sensing receptors leads to cardiovascular disease. Circulation 141:1393–1403. https://doi.org/10.1161/CIRCULATIONAHA. 119.043081 48. Li Q, Chen H, Zhang M, Wu T, Liu R. 2019. Altered short chain fatty acid profiles induced by dietary fiber intervention regulate AMPK levels and intestinal homeostasis. Food Funct 10:7174–7187. https://doi.org/10. 1039/C9FO01465A 49. Łoniewska B, Fraszczyk-Tousty M, Tousty P, Skonieczna-Żydecka K, Maciejewska-Markiewicz D, Łoniewski I. 2023. Analysis of fecal short- chain fatty acids (SCFAs) in healthy children during the first two years of life: an observational prospective cohort study. Nutrients 15:367. https:/ /doi.org/10.3390/nu15020367 50. Deleu S, Arnauts K, Deprez L, Machiels K, Ferrante M, Huys GRB, Thevelein JM, Raes J, Vermeire S. 2023. High acetate concentration protects intestinal barrier and exerts anti-inflammatory effects in organoid-derived epithelial monolayer cultures from patients with ulcerative colitis. Int J Mol Sci 24:768. https://doi.org/10.3390/ ijms24010768 51. Byrne CS, Chambers ES, Alhabeeb H, Chhina N, Morrison DJ, Preston T, Tedford C, Fitzpatrick J, Irani C, Busza A, Garcia-Perez I, Fountana S, Holmes E, Goldstone AP, Frost GS. 2016. Increased colonic propionate reduces anticipatory reward responses in the human striatum to high- energy foods. Am J Clin Nutr 104:5–14. https://doi.org/10.3945/ajcn. 115.126706 52. Todesco T, Rao AV, Bosello O, Jenkins DJ. 1991. Propionate lowers blood glucose and alters lipid metabolism in healthy subjects. Am J Clin Nutr 54:860–865. https://doi.org/10.1093/ajcn/54.5.860 53. Grüter T, Mohamad N, Rilke N, Blusch A, Sgodzai M, Demir S, Pedreiturria X, Lemhoefer K, Gisevius B, Haghikia A, Fisse AL, Motte J, Gold R, Pitarokoili K. 2023. Propionate exerts neuroprotective and neuroregenerative effects in the peripheral nervous system. Proc Natl Acad Sci U S A 120:e2216941120. https://doi.org/10.1073/pnas. 2216941120 54. Ryu TY, Kim K, Han T-S, Lee M-O, Lee J, Choi J, Jung KB, Jeong E-J, An DM, Jung C-R, Lim JH, Jung J, Park K, Lee M-S, Kim M-Y, Oh SJ, Hur K, Hamamoto R, Park D-S, Kim D-S, Son M-Y, Cho H-S. 2022. Human gut- microbiome-derived propionate coordinates proteasomal degradation via HECTD2 upregulation to target EHMT2 in colorectal cancer. ISME J 16:1205–1221. https://doi.org/10.1038/s41396-021-01119-1 55. Oncel S, Safratowich BD, Lindlauf JE, Liu Z, Palmer DG, Briske-Anderson M, Zeng H. 2024. Efficacy of butyrate to inhibit colonic cancer cell growth is cell type-specific and apoptosis-dependent. Nutrients 16:529. https://doi.org/10.3390/nu16040529 56. Zhang Y, Xu S, Qian Y, Mo C, Ai P, Yang X, Xiao Q. 2023. Sodium butyrate ameliorates gut dysfunction and motor deficits in a mouse model of Parkinson’s disease by regulating gut microbiota. Front Aging Neurosci 15:1099018. https://doi.org/10.3389/fnagi.2023.1099018 57. Donohoe DR, Garge N, Zhang X, Sun W, O’Connell TM, Bunger MK, Bultman SJ. 2011. The microbiome and butyrate regulate energy metabolism and autophagy in the mammalian colon. Cell Metab 13:517–526. https://doi.org/10.1016/j.cmet.2011.02.018 58. Gual-Grau A, Guirro M, Crescenti A, Boqué N, Arola L. 2022. In vitro fermentability of a broad range of natural ingredients by fecal microbiota from lean and obese individuals: potential health benefits . Int J Food Sci Nutr 73:195–209. https://doi.org/10.1080/09637486.2021. 1954144 59. Liu G, Liang L, Yu G, Li Q. 2018. Pumpkin polysaccharide modifies the gut microbiota during alleviation of type 2 diabetes in rats. Int J Biol Macromol 115:711–717. https://doi.org/10.1016/j.ijbiomac.2018.04.127 60. Wang B, Yu H, He Y, Wen L, Gu J, Wang X, Miao X, Qiu G, Wang H. 2021. Effect of soybean insoluble dietary fiber on prevention of obesity in high-fat diet fed mice via regulation of the gut microbiota. Food Funct 12:7923–7937. https://doi.org/10.1039/d1fo00078k 61. Bie N, Duan S, Meng M, Guo M, Wang C. 2021. Regulatory effect of non- starch polysaccharides from purple sweet potato on intestinal microbiota of mice with antibiotic-associated diarrhea. Food Funct 12:5563–5575. https://doi.org/10.1039/d0fo03465g 62. de Albuquerque TMR, Magnani M, Lima MDS, Castellano LRC, de Souza EL. 2021. Effects of digested flours from four different sweet potato (Ipomoea batatas L.) root varieties on the composition and metabolic activity of human colonic microbiota in vitro. J Food Sci 86:3707–3719. https://doi.org/10.1111/1750-3841.15852 63. Cao Y, Tian B, Zhang Z, Yang K, Cai M, Hu W, Guo Y, Xia Q, Wu W. 2022. Positive effects of dietary fiber from sweet potato [Ipomoea batatas (L.) Lam.] peels by different extraction methods on human fecal microbiota in vitro fermentation. Front Nutr 9:986667. https://doi.org/10.3389/fnut. 2022.986667 64. Zhang X, Yang Y, Wu Z, Weng P. 2016. The modulatory effect of anthocyanins from purple sweet potato on human intestinal microbiota in vitro. J Agric Food Chem 64:2582–2590. https://doi.org/ 10.1021/acs.jafc.6b00586 65. Cao L, Lee SG, Melough MM, Sakaki JR, Maas KR, Koo SI, Chun OK. 2020. Long-term blackcurrant supplementation modified gut microbiome profiles in mice in an age-dependent manner: an exploratory study. Nutrients 12:290. https://doi.org/10.3390/nu12020290 66. Yang H-J, Zhang T, Wu X-G, Kim M-J, Kim Y-H, Yang E-S, Yoon Y-S, Park S. 2021. Aqueous blackcurrant extract improves insulin sensitivity and secretion and modulates the gut microbiome in non-obese type 2 diabetic rats. Antioxidants (Basel) 10:756. https://doi.org/10.3390/ antiox10050756 67. Jakobsdottir G, Blanco N, Xu J, Ahrné S, Molin G, Sterner O, Nyman M. 2013. Formation of short-chain Fatty acids, excretion of anthocyanins, and microbial diversity in rats fed blackcurrants, blackberries, and raspberries. J Nutr Metab 2013:202534. https://doi.org/10.1155/2013/ 202534 68. Núñez-Gómez V, Periago MJ, Navarro-González I, Campos-Cava MP, Baenas N, González-Barrio R. 2021. Influence of raspberry and its dietary fractions on the in vitro activity of the colonic microbiota from normal and overweight subjects. Plant Foods Hum Nutr 76:494–500. https://doi.org/10.1007/s11130-021-00923-6 69. Pritchard JR, Lawrence GJ, Larroque O, Li Z, Laidlaw HK, Morell MK, Rahman S. 2011. A survey of β‐glucan and arabinoxylan content in wheat. J Sci Food Agric 91:1298–1303. https://doi.org/10.1002/jsfa.4316 70. Walton GE, Lu C, Trogh I, Arnaut F, Gibson GR. 2012. A randomised, double-blind, placebo controlled cross-over study to determine the gastrointestinal effects of consumption of arabinoxylan-oligosacchar­ ides enriched bread in healthy volunteers. Nutr J 11:36. https://doi.org/ 10.1186/1475-2891-11-36 71. Hald S, Schioldan AG, Moore ME, Dige A, Lærke HN, Agnholt J, Bach Knudsen KE, Hermansen K, Marco ML, Gregersen S, Dahlerup JF. 2016. Effects of arabinoxylan and resistant starch on intestinal microbiota and short-chain fatty acids in subjects with metabolic syndrome: a randomised crossover study. PLoS One 11:e0159223. https://doi.org/10. 1371/journal.pone.0159223 72. Heimann E, Nyman M, Pålbrink A-K, Lindkvist-Petersson K, Degerman E. 2016. Branched short-chain fatty acids modulate glucose and lipid metabolism in primary adipocytes. Adipocyte 5:359–368. https://doi. org/10.1080/21623945.2016.1252011 73. Macfarlane GT, Gibson GR, Beatty E, Cummings JH. 1992. Estimation of short-chain fatty acid production from protein by human intestinal bacteria based on branched-chain fatty acid measurements. FEMS Microbiol Ecol 10:81–88. https://doi.org/10.1111/j.1574-6941.1992. tb00002.x 74. Richardson AJ, McKain N, Wallace RJ. 2013. Ammonia production by human faecal bacteria, and the enumeration, isolation and characteri­ zation of bacteria capable of growth on peptides and amino acids. BMC Microbiol 13:6. https://doi.org/10.1186/1471-2180-13-6 75. Rios-Covian D, González S, Nogacka AM, Arboleya S, Salazar N, Gueimonde M, de Los Reyes-Gavilán CG. 2020. An overview on fecal branched short-chain fatty acids along human life and as related with body mass index: associated dietary and anthropometric factors. Front Microbiol 11:973. https://doi.org/10.3389/fmicb.2020.00973 76. Mayengbam S, Lambert JE, Parnell JA, Tunnicliffe JM, Nicolucci AC, Han J, Sturzenegger T, Shearer J, Mickiewicz B, Vogel HJ, Madsen KL, Reimer Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2418 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1038/s41588-019-0350-x https://doi.org/10.1161/CIRCULATIONAHA.119.043081 https://doi.org/10.1039/C9FO01465A https://doi.org/10.3390/nu15020367 https://doi.org/10.3390/ijms24010768 https://doi.org/10.3945/ajcn.115.126706 https://doi.org/10.1093/ajcn/54.5.860 https://doi.org/10.1073/pnas.2216941120 https://doi.org/10.1038/s41396-021-01119-1 https://doi.org/10.3390/nu16040529 https://doi.org/10.3389/fnagi.2023.1099018 https://doi.org/10.1016/j.cmet.2011.02.018 https://doi.org/10.1080/09637486.2021.1954144 https://doi.org/10.1016/j.ijbiomac.2018.04.127 https://doi.org/10.1039/d1fo00078k https://doi.org/10.1039/d0fo03465g https://doi.org/10.1111/1750-3841.15852 https://doi.org/10.3389/fnut.2022.986667 https://doi.org/10.1021/acs.jafc.6b00586 https://doi.org/10.3390/nu12020290 https://doi.org/10.3390/antiox10050756 https://doi.org/10.1155/2013/202534 https://doi.org/10.1007/s11130-021-00923-6 https://doi.org/10.1002/jsfa.4316 https://doi.org/10.1186/1475-2891-11-36 https://doi.org/10.1371/journal.pone.0159223 https://doi.org/10.1080/21623945.2016.1252011 https://doi.org/10.1111/j.1574-6941.1992.tb00002.x https://doi.org/10.1186/1471-2180-13-6 https://doi.org/10.3389/fmicb.2020.00973 https://doi.org/10.1128/msystems.00577-24 RA. 2019. Impact of dietary fiber supplementation on modulating microbiota-host-metabolic axes in obesity. J Nutr Biochem 64:228–236. https://doi.org/10.1016/j.jnutbio.2018.11.003 77. Calatayud M, Van den Abbeele P, Ghyselinck J, Marzorati M, Rohs E, Birkett A. 2021. Comparative effect of 22 dietary sources of fiber on gut microbiota of healthy humans in vitro; Front Nutr 8:700571. https://doi. org/10.3389/fnut.2021.700571 78. Choi B-Y, Daniel N, Houde VP, Ouellette A, Marcotte B, Varin TV, Vors C, Feutry P, Ilkayeva O, Ståhlman M, St-Pierre P, Bäckhed F, Tremblay A, White PJ, Marette A. 2021. Feeding diversified protein sources exacerbates hepatic insulin resistance via increased gut microbial branched-chain fatty acids and mTORC1 signaling in obese mice. Nat Commun 12:3377. https://doi.org/10.1038/s41467-021-23782-w 79. Cheng HY, Chan JCY, Yap GC, Huang C-H, Kioh DYQ, Tham EH, Loo EXL, Shek LPC, Karnani N, Goh A, Van Bever HPS, Teoh OH, Chan YH, Lay C, Knol J, Yap F, Tan KH, Chong Y-S, Godfrey KM, Chan ECY, Lee BW, Ta LDH. 2022. Evaluation of stool short chain fatty acids profiles in the first year of life with childhood atopy-related outcomes. Front Allergy 3:873168. https://doi.org/10.3389/falgy.2022.873168 80. Sasaki M, Suaini NHA, Afghani J, Heye KN, O’Mahony L, Venter C, Lauener R, Frei R, Roduit C. 2024. Systematic review of the association between short‐chain fatty acids and allergic diseases. Allergy 79:1789–1811. https://doi.org/10.1111/all.16065 81. Pekmez CT, Larsson MW, Lind MV, Vazquez Manjarrez N, Yonemitsu C, Larnkjaer A, Bode L, Mølgaard C, Michaelsen KF, Dragsted LO. 2020. Breastmilk lipids and oligosaccharides influence branched short-chain fatty acid concentrations in infants with excessive weight gain. Mol Nutr Food Res 64:e1900977. https://doi.org/10.1002/mnfr.201900977 82. Xu R, Tan C, He Y, Wu Q, Wang H, Yin J. 2020. Dysbiosis of gut microbiota and short-chain fatty acids in encephalitis: a Chinese pilot study. Front Immunol 11:1994. https:​/​/​doi.org/​10.3389/​fimmu.2020. 01994 83. Chen S-J, Chen C-C, Liao H-Y, Lin Y-T, Wu Y-W, Liou J-M, Wu M-S, Kuo C- H, Lin C-H. 2022. Association of fecal and plasma levels of short-chain fatty acids with gut microbiota and clinical severity in patients with Parkinson disease. Neurology (ECronicon) 98:e848–e858. https://doi. org/10.1212/WNL.0000000000013225 84. Ezzat-Zadeh Z, Henning SM, Yang J, Woo SL, Lee R-P, Huang J, Thames G, Gilbuena I, Tseng C-H, Heber D, Li Z. 2021. California strawberry consumption increased the abundance of gut microorganisms related to lean body weight, health and longevity in healthy subjects. Nutr Res 85:60–70. https://doi.org/10.1016/j.nutres.2020.12.006 85. Han Y, Song M, Gu M, Ren D, Zhu X, Cao X, Li F, Wang W, Cai X, Yuan B, Goulette T, Zhang G, Xiao H. 2019. Dietary intake of whole strawberry inhibited colonic inflammation in dextran-sulfate-sodium-treated mice via restoring immune homeostasis and alleviating gut microbiota dysbiosis. J Agric Food Chem 67:9168–9177. https://doi.org/10.1021/ acs.jafc.8b05581 86. De Angelis M, Ferrocino I, Calabrese FM, De Filippis F, Cavallo N, Siragusa S, Rampelli S, Di Cagno R, Rantsiou K, Vannini L, Pellegrini N, Lazzi C, Turroni S, Lorusso N, Ventura M, Chieppa M, Neviani E, Brigidi P, O’Toole PW, Ercolini D, Gobbetti M, Cocolin L. 2020. Diet influences the functions of the human intestinal microbiome. Sci Rep 10:4247. https:// doi.org/10.1038/s41598-020-61192-y 87. Meslier V, Laiola M, Roager HM, De Filippis F, Roume H, Quinquis B, Giacco R, Mennella I, Ferracane R, Pons N, Pasolli E, Rivellese A, Dragsted LO, Vitaglione P, Ehrlich SD, Ercolini D. 2020. Mediterranean diet intervention in overweight and obese subjects lowers plasma cholesterol and causes changes in the gut microbiome and metabo­ lome independently of energy intake. Gut 69:1258–1268. https://doi. org/10.1136/gutjnl-2019-320438 88. Raman AS, Gehrig JL, Venkatesh S, Chang H-W, Hibberd MC, Subramanian S, Kang G, Bessong PO, Lima AAM, Kosek MN, Petri WA, Rodionov DA, Arzamasov AA, Leyn SA, Osterman AL, Huq S, Mostafa I, Islam M, Mahfuz M, Haque R, Ahmed T, Barratt MJ, Gordon JI. 2019. A sparse covarying unit that describes healthy and impaired human gut microbiota development. Science 365:eaau4735. https://doi.org/10. 1126/science.aau4735 89. Chen RY, Mostafa I, Hibberd MC, Das S, Mahfuz M, Naila NN, Islam MM, Huq S, Alam MA, Zaman MU, Raman AS, Webber D, Zhou C, Sundaresan V, Ahsan K, Meier MF, Barratt MJ, Ahmed T, Gordon JI. 2021. A microbiota-directed food intervention for undernourished children. N Engl J Med 384:1517–1528. https://doi.org/10.1056/NEJMoa2023294 90. Homann C-M, Rossel CAJ, Dizzell S, Bervoets L, Simioni J, Li J, Gunn E, Surette MG, de Souza RJ, Mommers M, Hutton EK, Morrison KM, Penders J, van Best N, Stearns JC. 2021. Infants’ first solid foods: impact on gut microbiota development in two intercontinental cohorts. Nutrients 13:2639. https://doi.org/10.3390/nu13082639 91. Delgado A, Issaoui M, Vieira MC, Saraiva de Carvalho I, Fardet A. 2021. Food composition databases: does it matter to human health? Nutrients 13:2816. https://doi.org/10.3390/nu13082816 92. Noronha A, Modamio J, Jarosz Y, Guerard E, Sompairac N, Preciat G, Daníelsdóttir AD, Krecke M, Merten D, Haraldsdóttir HS, et al. 2019. The Virtual Metabolic Human database: integrating human and gut microbiome metabolism with nutrition and disease. Nucleic Acids Res 47:D614–D624. https://doi.org/10.1093/nar/gky992 93. Haytowitz DB, Ahuja JKC, Showell B, Somanchi M, Nickle M, Nguyen QA, Williams JR, Roseland JM, Khan M, Patterson KY, Exler J, Wasswa-Kintu S, Thomas R, Pehrsson PR. 2015. Composition of foods raw, processed, prepared USDA national nutrient database for standard reference, release 28. Nutrient Data Laboratory, Beltsville Human Nutrition Research Center, ARS, USDA. 94. Rutayisire E, Huang K, Liu Y, Tao F. 2016. The mode of delivery affects the diversity and colonization pattern of the gut microbiota during the first year of infants' life: a systematic review. BMC Gastroenterol 16:86. https://doi.org/10.1186/s12876-016-0498-0 95. Taylor R, Keane D, Borrego P, Arcaro K. 2023. Effect of maternal diet on maternal milk and breastfed infant gut microbiomes: a scoping review. Nutrients 15:1420. https://doi.org/10.3390/nu15061420 96. Joseph C, Zafeiropoulos H, Bernaerts K, Faust K. 2024. Predicting microbial interactions with approaches based on flux balance analysis: an evaluation. BMC Bioinformatics 25:36. https://doi.org/10.1186/ s12859-024-05651-7 97. Johnson JS, Spakowicz DJ, Hong B-Y, Petersen LM, Demkowicz P, Chen L, Leopold SR, Hanson BM, Agresta HO, Gerstein M, Sodergren E, Weinstock GM. 2019. Evaluation of 16S rRNA gene sequencing for species and strain-level microbiome analysis. Nat Commun 10:5029. https://doi.org/10.1038/s41467-019-13036-1 98. Flynn KJ, Ruffin MT IV, Turgeon DK, Schloss PD. 2018. Spatial variation of the native colon microbiota in healthy adults. Cancer Prev Res (Phila) 11:393–402. https://doi.org/10.1158/1940-6207.CAPR-17-0370 99. Heinken A, Hertel J, Acharya G, Ravcheev DA, Nyga M, Okpala OE, Hogan M, Magnúsdóttir S, Martinelli F, Nap B, Preciat G, Edirisinghe JN, Henry CS, Fleming RMT, Thiele I. 2023. Genome-scale metabolic reconstruction of 7,302 human microorganisms for personalized medicine. Nat Biotechnol 41:1320–1331. https://doi.org/10.1038/ s41587-022-01628-0 100. Quinn-Bohmann N, Wilmanski T, Sarmiento KR, Levy L, Lampe JW, Gurry T, Rappaport N, Ostrem EM, Venturelli OS, Diener C, Gibbons SM. 2023. Microbial community-scale metabolic modeling predicts personalized short chain fatty acid production profiles in the human gut. bioRxiv:2023.02.28.530516. https://doi.org/10.1101/2023.02.28. 530516 101. Scott WT, Benito-Vaquerizo S, Zimmermann J, Bajić D, Heinken A, Suarez-Diez M, Schaap PJ. 2023. A structured evaluation of genome- scale constraint-based modeling tools for microbial consortia. PLoS Comput Biol 19:e1011363. https://doi.org/10.1371/journal.pcbi. 1011363 102. Orth JD, Thiele I, Palsson BØ. 2010. What is flux balance analysis? Nat Biotechnol 28:245–248. https://doi.org/10.1038/nbt.1614 103. Dürholz K, Hofmann J, Iljazovic A, Häger J, Lucas S, Sarter K, Strowig T, Bang H, Rech J, Schett G, Zaiss MM. 2020. Dietary short-term fiber interventions in arthritis patients increase systemic SCFA levels and regulate inflammation. Nutrients 12:3207. https://doi.org/10.3390/ nu12103207 104. Henson MA, Hanly TJ. 2014. Dynamic flux balance analysis for synthetic microbial communities. IET Syst Biol 8:214–229. https://doi.org/10. 1049/iet-syb.2013.0021 105. Popp D, Centler F. 2020. μBialSim: constraint-based dynamic simulation of complex microbiomes. Front Bioeng Biotechnol 8:574. https://doi. org/10.3389/fbioe.2020.00574 106. Yang J, Martínez I, Walter J, Keshavarzian A, Rose DJ. 2013. In vitro characterization of the impact of selected dietary fibers on fecal Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2419 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1016/j.jnutbio.2018.11.003 https://doi.org/10.3389/fnut.2021.700571 https://doi.org/10.1038/s41467-021-23782-w https://doi.org/10.3389/falgy.2022.873168 https://doi.org/10.1111/all.16065 https://doi.org/10.1002/mnfr.201900977 https://doi.org/10.3389/fimmu.2020.01994 https://doi.org/10.1212/WNL.0000000000013225 https://doi.org/10.1016/j.nutres.2020.12.006 https://doi.org/10.1021/acs.jafc.8b05581 https://doi.org/10.1038/s41598-020-61192-y https://doi.org/10.1136/gutjnl-2019-320438 https://doi.org/10.1126/science.aau4735 https://doi.org/10.1056/NEJMoa2023294 https://doi.org/10.3390/nu13082639 https://doi.org/10.3390/nu13082816 https://doi.org/10.1093/nar/gky992 https://doi.org/10.1186/s12876-016-0498-0 https://doi.org/10.3390/nu15061420 https://doi.org/10.1186/s12859-024-05651-7 https://doi.org/10.1038/s41467-019-13036-1 https://doi.org/10.1158/1940-6207.CAPR-17-0370 https://doi.org/10.1038/s41587-022-01628-0 https://doi.org/10.1101/2023.02.28.530516 https://doi.org/10.1371/journal.pcbi.1011363 https://doi.org/10.1038/nbt.1614 https://doi.org/10.3390/nu12103207 https://doi.org/10.1049/iet-syb.2013.0021 https://doi.org/10.3389/fbioe.2020.00574 https://doi.org/10.1128/msystems.00577-24 microbiota composition and short chain fatty acid production. Anaerobe 23:74–81. https://doi.org/10.1016/j.anaerobe.2013.06.012 107. Smith EA, Macfarlane GT. 1998. Enumeration of amino acid fermenting bacteria in the human large intestine: effects of pH and starch on peptide metabolism and dissimilation of amino acids. FEMS Microbiol Ecol 25:355–368. https://doi.org/10.1111/j.1574-6941.1998.tb00487.x 108. Duncan SH, Barcenilla A, Stewart CS, Pryde SE, Flint HJ. 2002. Acetate utilization and butyryl coenzyme A (CoA):acetate-CoA transferase in butyrate-producing bacteria from the human large intestine. Appl Environ Microbiol 68:5186–5190. https://doi.org/10.1128/AEM.68.10. 5186-5190.2002 109. Reichardt N, Duncan SH, Young P, Belenguer A, McWilliam Leitch C, Scott KP, Flint HJ, Louis P. 2014. Phylogenetic distribution of three pathways for propionate production within the human gut microbiota. ISME J 8:1323–1335. https://doi.org/10.1038/ismej.2014.14 110. Duncan SH, Belenguer A, Holtrop G, Johnstone AM, Flint HJ, Lobley GE. 2007. Reduced dietary intake of carbohydrates by obese subjects results in decreased concentrations of butyrate and butyrate- producing bacteria in feces. Appl Environ Microbiol 73:1073–1078. https://doi.org/10.1128/AEM.02340-06 111. Johnson AJ, Vangay P, Al-Ghalith GA, Hillmann BM, Ward TL, Shields- Cutler RR, Kim AD, Shmagel AK, Syed AN, Walter J, Menon R, Koecher K, Knights D. 2019. Daily sampling reveals personalized diet-microbiome associations in humans. Cell Host Microbe 25:789–802. https://doi.org/ 10.1016/j.chom.2019.05.005 112. Venkataraman A, Sieber JR, Schmidt AW, Waldron C, Theis KR, Schmidt TM. 2016. Variable responses of human microbiomes to dietary supplementation with resistant starch. Microbiome 4:33. https://doi. org/10.1186/s40168-016-0178-x 113. Yang J, Rose DJ. 2014. Long-term dietary pattern of fecal donor correlates with butyrate production and markers of protein fermenta­ tion during in vitro fecal fermentation. Nutr Res 34:749–759. https://doi. org/10.1016/j.nutres.2014.08.006 114. Ziemski M, Adamov A, Kim L, Flörl L, Bokulich NA. 2022. Reproducible acquisition, management and meta-analysis of nucleotide sequence (meta)data using q2-fondue. Bioinformatics 38:5081–5091. https://doi. org/10.1093/bioinformatics/btac639 115. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, Alexander H, Alm EJ, Arumugam M, Asnicar F, et al. 2019. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol 37:852–857. https://doi.org/10.1038/s41587- 019-0209-9 116. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. 2016. DADA2: high-resolution sample inference from Illumina amplicon data. Nat Methods 13:581–583. https://doi.org/10.1038/nmeth.3869 117. McDonald D, Jiang Y, Balaban M, Cantrell K, Zhu Q, Gonzalez A, Morton JT, Nicolaou G, Parks DH, Karst S, Albertsen M, Hugenholtz P, DeSantis T, Song SJ, Bartko A, Havulinna AS, Jousilahti P, Cheng S, Inouye M, Niiranen T, Jain M, Salomaa V, Lahti L, Mirarab S, Knight R. 2022. Greengenes2 enables a shared data universe for microbiome studies. Bioinformatics. https://doi.org/10.1101/2022.12.19.520774 118. Growth charts - WHO child growth standards. 2022. Available from: https://www.cdc.gov/growthcharts/who_charts.htm. Retrieved 31 May 2023. 119. Heinig MJ, Nommsen LA, Peerson JM, Lonnerdal B, Dewey KG. 1993. Energy and protein intakes of breast-fed and formula-fed infants during the first year of life and their association with growth velocity: the DARLING study. Am J Clin Nutr 58:152–161. https://doi.org/10.1093/ ajcn/58.2.152 120. Brunk E, Sahoo S, Zielinski DC, Altunkaya A, Dräger A, Mih N, Gatto F, Nilsson A, Preciat Gonzalez GA, Aurich MK, Prlić A, Sastry A, Danielsdot­ tir AD, Heinken A, Noronha A, Rose PW, Burley SK, Fleming RMT, Nielsen J, Thiele I, Palsson BO. 2018. Recon3D enables a three-dimensional view of gene variation in human metabolism. Nat Biotechnol 36:272–281. https://doi.org/10.1038/nbt.4072 121. Kerimi A, Kraut NU, da Encarnacao JA, Williamson G. 2020. The gut microbiome drives inter- and intra-individual differences in metabolism of bioactive small molecules. Sci Rep 10:19590. https://doi.org/10.1038/ s41598-020-76558-5 122. Zhu A, Sunagawa S, Mende DR, Bork P. 2015. Inter-individual differences in the gene content of human gut bacterial species. Genome Biol 16:82. https://doi.org/10.1186/s13059-015-0646-9 123. Łoniewski I, Skonieczna-Żydecka K, Stachowska L, Fraszczyk-Tousty M, Tousty P, Łoniewska B. 2022. Breastfeeding affects concentration of faecal short chain fatty acids during the first year of life: results of the systematic review and meta-analysis. Front Nutr 9:939194. https://doi. org/10.3389/fnut.2022.939194 Research Article mSystems September 2024 Volume 9 Issue 9 10.1128/msystems.00577-2420 D ow nl oa de d fr om h ttp s: //j ou rn al s. as m .o rg /jo ur na l/m sy st em s on 0 8 O ct ob er 2 02 4 by 1 25 .2 37 .8 2. 15 5. https://doi.org/10.1016/j.anaerobe.2013.06.012 https://doi.org/10.1111/j.1574-6941.1998.tb00487.x https://doi.org/10.1128/AEM.68.10.5186-5190.2002 https://doi.org/10.1038/ismej.2014.14 https://doi.org/10.1128/AEM.02340-06 https://doi.org/10.1016/j.chom.2019.05.005 https://doi.org/10.1186/s40168-016-0178-x https://doi.org/10.1016/j.nutres.2014.08.006 https://doi.org/10.1093/bioinformatics/btac639 https://doi.org/10.1038/s41587-019-0209-9 https://doi.org/10.1038/nmeth.3869 https://doi.org/10.1101/2022.12.19.520774 https://www.cdc.gov/growthcharts/who_charts.htm https://doi.org/10.1093/ajcn/58.2.152 https://doi.org/10.1038/nbt.4072 https://doi.org/10.1038/s41598-020-76558-5 https://doi.org/10.1186/s13059-015-0646-9 https://doi.org/10.3389/fnut.2022.939194 https://doi.org/10.1128/msystems.00577-24 Food-breastmilk combinations alter the colonic microbiome of weaning infants: an in silico study RESULTS Predicted fluxes of SCFAs and BCFAs Predict