Rice bran supplementation modulates growth, microbiota and metabolome in weaning infants: a clinical trial in Nicaragua and Mali

Rice bran supplementation provides nutrients, prebiotics and phytochemicals that enhance gut immunity, reduce enteric pathogens and diarrhea, and warrants attention for improvement of environmental enteric dysfunction (EED) in children. EED is a subclinical condition associated with stunting due to impaired nutrient absorption. This study investigated the effects of rice bran supplementation on weight for age and length for age z-scores (WAZ, LAZ), EED stool biomarkers, as well as microbiota and metabolome signatures in weaning infants from 6 to 12 months old that reside in Nicaragua and Mali. Healthy infants were randomized to a control (no intervention) or a rice bran group that received daily supplementation with increasing doses at each month (1–5 g/day). Stool microbiota were characterized using 16S rDNA amplicon sequencing. Stool metabolomes were analyzed using ultra-high-performance liquid-chromatography tandem mass-spectrometry. Statistical comparisons were completed at 6, 8, and 12 months of age. Daily consumption of rice bran was safe and feasible to support changes in LAZ from 6–8 and 8–12 months of age in Nicaragua and Mali infants when compared to control. WAZ was significantly improved only for Mali infants at 8 and 12 months. Mali and Nicaraguan infants showed major differences in the overall gut microbiota and metabolome composition and structure at baseline, and thus each country cohort demonstrated distinct microbial and metabolite profile responses to rice bran supplementation when compared to control. Rice bran is a practical dietary intervention strategy that merits development in rice-growing regions that have a high prevalence of growth stunting due to malnutrition and diarrheal diseases. Rice is grown as a staple food, and the bran is used as animal feed or wasted in many low- and middle-income countries where EED and stunting is prevalent.


Study participant characteristics.
To study the effect of daily rice bran supplementation, monthly stool samples from 47 Nicaraguan and 48 Malian children were collected (average of 7 samples per child, total of 567 samples). The flow and number of infants from study recruitment to study completion is shown in Fig. 1. Baseline participant characteristics for Nicaragua and Mali are shown in Table 2. We collected information on demographic factors and infants' household characteristics. In Nicaragua, 54.2% of infants were born via caesarean section in the control group and 30.4% in the rice bran group. All participants from Mali were delivered vaginally. For breastfeeding status, 96% of the control group and 83% in the rice bran group were consuming breast milk at six months old in Nicaragua and all children in the Mali group were consuming breast milk at the beginning and throughout the study. Of the 95 infants enrolled, 52 received antibiotics with 87 total antibiotic courses in the six-month period. Most courses consisted of systemic antibiotics given orally, with some delivered by injection for respiratory, skin, ear or diarrheal infections.
Anthropometric growth measurements. Anthropometric data were collected using standardized procedures across study sites at 6, 8, and 12 months of age for each child and included length-for-age Z-score (LAZ), weight-for-age Z-score (WAZ) and weight-for-length Z-score (WLZ). Table 3 reports the z-scores analyzed for each country by repeated measures and adjusted by treatment (rice bran) and the rate of growth was significantly changed by age groups 6-8 months and 8-12 months. Significant differences were observed for anthropometric measures between treatment and ages in both cohorts. In Nicaraguan infants who consumed rice bran, LAZ was significantly changed over time (p-value = 0.0002) and for WLZ at 6-8 months (p-value < 0.0001). The Malian www.nature.com/scientificreports www.nature.com/scientificreports/ infants who consumed rice bran had significant growth results for WAZ (6-8 months, p-value = 0.0001, 8-12 months, p-value = 0.0175) and WLZ (6-8 months, p-value = 0.0141, 8-12 months, p-value = 0.0134). Figure 2 displays LAZ, WAZ and WLZ for each age group and by country. The significant changes observed in LAZ at 8 months and 12 months in Nicaraguan infants who consumed rice bran was compared to control group ( Fig. 2A, p < 0.01). No significant differences were observed for WAZ and WLZ with this control group comparison (Fig. 2B,C).

Diarrheal episodes and stool environmental enteric dysfunction (eeD) biomarkers.
For the Mali cohort, the diarrheal disease incidence was 33% in the rice bran intervention group compared to 79% in the control group over the 6-month study period. There was one rice bran intervention participant who experienced a single repeat episode of diarrhea, and 5 control group participants experienced between 2 and 4 repeat episodes of diarrhea during the 6-month trial. There were a higher number of diarrheal episodes in the Mali control group infants when compared to rice bran intervention group. No differences were detected for diarrheal episodes for the Nicaragua rice bran and control groups. We next evaluated four EED stool biomarkers at 6, 8 and 12 months of age using ELISA (see materials and methods). A significant decrease in AAT was observed at 12 months of age (p = 0.0368) in Nicaragua infants who consumed rice bran compared to control (  [43][44][45][46][47] . Sequences were preprocessed for quality assurance and classified into operational taxonomic units (OTUs) (see materials and methods), and the results were integrated to construct family and genus-level composition profiles for all 192 samples from both countries. No major differences were detected in alpha diversity indices (Observed, Shannon, InvSimpson and Richness) calculated for rice bran group and control group at 8 and 12 months (Table S1). A complete country-level separation in the overall gut microbial community composition was depicted for beta diversity using Nonmetric Multidimensional Scaling (NMDS) plot based on the Bray-Curtis distance measure (Fig. 3A). The magnitude of age and country-level differences in infant gut microbiota composition was expected, as differences were previously reported that compared these geographical regions 21 . The differences in microbiota composition at baseline provided compelling rationale for separately evaluating the control and rice bran responsive microbial communities in infants at each age and for each country. Figure 3B shows the NMDS plot separated by country. This figure highlights differences in the microbiota between two time periods, 8 and 12 months, which is more pronounced in the Malian samples. There is overlap of microbiota within each time period that indicates putative similarity between microbial community structures during growth. The overlap was observed to a greater level at 8 months of age as compared to 12 months that may illustrate microbial adaptation to new exposures 33 . Figure 4 illustrates the OTU with at least 2 log-fold change between the rice bran and control groups per country and per age group. These OTU were ordered based on statistical significance (FDR-adjusted p-value, from bottom (lowest p-value) to top). See additional details provided in Table S2 for Nicaragua and Table S3 for Mali. The effects of daily dietary rice bran supplementation on the infant microbiota were evaluated at 8 months of age for each country because this was the time point in which measurable growth differences (LAZ) were already observed for the rice bran group compared to control.
Gut microbiota changes that were reflective of a dietary response to rice bran in both the Mali and Nicaragua cohorts involved 12 OTUs at either 8 or 12 months of age when compared to control. The highest area of overlap occurred for both Lactobacillaceae_Lactobacillus_Otu0024 and the Lactobacillus_Otu0053 (Table S4). There was additional overlap in both countries with respect to microbiota changes following rice bran intake that included Bifidobacterium, Faecalibacterium, and Lachnospiriaceae. infant stool metabolomes. A total of 309 stool samples were analyzed for metabolomics from this 6-month prospective study to evaluate the effects of rice bran supplementation compared to control. Stool metabolite analysis at 8 months of age in Nicaraguan and Malian infants resulted in the detection of 1449 biochemicals, of which 1016 metabolites had confirmed names and 433 compounds were of unknown structural identity (see Table S5). ANOVA contrasts and Welch's two-sample t-test were used to identify biochemicals that significantly differed between experimental groups at 8-months of age. Table 5 lists the fold differences in the stool metabolite relative abundance between study diet groups at 8 months of age. There are 39 (Nicaragua) and 44 (Mali) stool metabolites with significant fold differences between children consuming rice bran for two months when compared to control. There were also 33 unknown metabolites for Nicaragua and 31 unknown metabolites for Mali that showed significant differences between groups (data shown in Table S5). Significant fold differences occurred for 15 amino acids, 2 peptides, 3 carbohydrates, 9 lipids, 1 cofactor and vitamin, and 9 xenobiotics (six of these metabolites were considered to be food components/plant-derived) in Nicaraguan children who consumed rice bran compared to control. At 8 months of age, the stool from Mali infants had 6 amino acids, 1 energy metabolite, 14 lipids, 6 cofactor and vitamins, 5 nucleotides and 12 xenobiotics (7 classified as food components/ plant-derived) that were significantly different between rice bran and control (Table 5).
There were 62 stool metabolites from the Nicaraguan cohort at 8 months that had lower relative abundances, and there were 10 stool metabolites with significantly increased fold differences in abundance between groups. The stool metabolites of food and nutritional relevance that were associated with increasing intake of rice bran were classified to the tryptophan metabolism pathway (e.g. indolepropionate), as well as monoacylglycerol (1-linolenoylglycerol) and diacylglycerol metabolic (linoleoyl-linolenoyl-glycerol) pathways.
In contrast to Nicaragua, there were 54 distinct stool metabolites from the Mali cohort that increased with significant fold differences at 8 months between rice bran and control infants. Other stool metabolites with gut health and nutritional importance that originated from rice bran food intake, included alpha-tocotrienol (vitamin E component), pyridoxine (vitamin B6), ferulic acid 4-sulfate, tyrosol, and N-acetyl sphingosine (Table 5). An estimated false discovery rate (q-value) was calculated to account for the multiple comparisons across metabolites that are typical of metabolomic-based studies.

Discussion
This study demonstrated that rice bran supplementation was feasible and safe for weaning infants with strong compliance in two distinct LMIC countries. Rice bran was well tolerated in a dose-escalation, supplementation trial regimen during the first 6 months of weaning, and was without side effects or adverse interactions. The distinct agro-ecological climates, dietary patterns, and birth methods between infants in West Africa and Latin America were considered highly relevant sources of variation that added scientific rigor to the generalizability of these clinical trial results with respect to safety, feasibility, and tolerability. However, when evaluating the infant gut microbiome and stool metabolite responses to increasing rice bran intake, this study showed distinct country and age specific variability in weaning infants. Rice bran supplementation in the diet supported the healthy growth of Nicaraguan and Malian infants with changes detected in LAZ between 6-8 months and 8-12 months of age. Nicaraguan infants had reduced gut permeability as measured by AAT at 12 months. Malian infants also showed changes over time for WAZ. The 95 healthy infants enrolled in this study had slightly higher WAZ, LAZ and WLZ scores than the documented country-wide averages that use the WHO scoring index 12 .
There is a growing body of scientific evidence for a strong relationship between EED and growth deficits in children 7,10,48 . EED biomarkers were selected in this study because high concentrations of stool AAT and MPO were associated with decreased growth in children 48 , Naylor et al. found that elevated AAT levels were associated with decreased oral rotavirus vaccine response 49 , and were associated with seroconversion to rotavirus vaccination in Nicaraguan infants 37 . Rice bran was shown to reduce AAT in neonatal pigs challenged with human Rotavirus infection 24 and was a finding of translational importance to this clinical trial because rice bran intake improved AAT levels in Nicaraguan infants ( Table 4). The lack of any statistically significant differences in EED stool biomarkers for Mali infants may be due to the higher number of overall diarrheal episodes that this cohort www.nature.com/scientificreports www.nature.com/scientificreports/ experienced, the level of variability within individuals, and breastfeeding status for all infants throughout the study period. Our findings did concur with EED biomarker concentrations reported in the MAL-ED cohort 36 . EED biomarkers merit continuous review for relevance with anthropometric measures due to extensive global variability in concentrations reported across studies 6,50-52 .
A study limitation and possible confounder of rice bran effects to modulate the infant gut microbiota was the percentage of exclusively breastfed infants and the number of vaginal deliveries being lower in the Nicaragua cohort. This was in contrast to the Mali cohort that had 100% of children who were breastfed and had vaginal delivery. These are key considerations to evaluating microbiota that has been previously characterized as less mature 53 , and with varied structure by geographical location 54 , diet deficiencies [55][56][57][58][59] , environmental exposures 54,60,61 and host factors 60,62 . Infants from each country showed favorable changes in both the stool microbiota and metabolome over time (with increasing age and growth), and between rice bran and control groups at 8 months of age.
Given that dietary rice bran intake has been previously shown to promote beneficial stool microbial communities, such as native gut probiotics in mice 25,63 , pigs 23,24,31 and adults 27,64,65 , differential abundance testing at the genus level assignment between rice bran and control fed infants was conducted at 8 months and 12 months of age. In Nicaragua, the Lactobacillus, Lachnospiraceae, Bifidobacterium, Ruminococcaceae and Veillonella were identified as responsive following rice bran consumption compared to an age and control matched group. These taxa have recognized saccharolytic mechanisms of action [66][67][68] , are known to produce and promote crossfeeding of short chain fatty acids 69,70 , as well as provide competitive inhibition of pathogen colonization 71,72 . The microbial enrichment of these communities has important implications for assessing how gut microbes metabolize rice bran that can support changes to infant LAZ growth outcomes. In Mali, the observed increase in the relative abundance of Lactobacillus in rice bran fed infants was consistent with prior studies of rice bran feeding to young animals. This finding was considered alongside evidence for human milk oligosaccharides that also promote Lactobacillus in breastfed infants 73 . Rice bran metabolism by Lactobacillus spp. was shown to produce a suite of small molecules that may be absorbed by the host or detected in stool. Stool metabolites from infants fed rice bran showed significant fold differences amongst several essential amino acids, cofactors and vitamins, lipids, phytochemicals, and in energy metabolism pathways when compared to the control. The stool metabolites that were highly likely to have originated from the rice bran consumed by both Nicaragua and Malian infants provided additional confirmation of compliance to the dietary intervention 74 . As predicted, we observed and reported exceptionally distinct profiles for both the stool microbiota and metabolome composition of infants between Mali and Nicaragua at all ages 41,54,75 and therefore separately discussed the response to rice bran supplementation by each region. For example, the nearly 5-fold increase of stool indolepropionate in rice bran infants compared to control from Nicaragua represents a tryptophan metabolite produced by the gut microbiota that may influence the developing immune system and intestinal homeostasis 76,77 . Increased levels of N-acetylmethionine (nutritionally and metabolically equal to L-methionine) and N-formylmethionine in Malian infants also represented rice bran derived amino acids required for normal growth and development 78 . There are several cofactors and vitamins from rice bran supplementation, such as alpha and gamma-tocotrienol and pyridoxine (vitamin B6) that merit attention for demonstrating multiple health benefits such as synthesis of amino acids and neurotransmitter precursors, as well as preventing anemia and skin problems 74,79,80 . Additional microbial digested food components in the stool metabolome that come from rice bran were ferulic acid 4-sulfate, indoleacetylaspartate, and Tyrosol. This trial had only supplemented   www.nature.com/scientificreports www.nature.com/scientificreports/ the rice bran for a 6 month window and according to WHO, growth assessment should be standardized and compared globally over the first two years of life 81 . We put forth that rice bran metabolism by host and gut microbes between 12-36 months of age should be captured in future assessments of early influences on growth velocity 81,82 . The focus on changes between 8 and 12 months in this study showed modifications in gut microbial communities and metabolites by rice bran intake, and suggests there will be long-term impact on the overall microbiota composition as it continues to develop and mature 83 .
This was the first randomized controlled trial of rice bran supplementation in LMIC infants and provides compelling rationale for continued follow-up investigation of rice bran supplementation for reducing risk of malnutrition, as well as for eliciting changes during child growth that protect against enteric pathogens and diarrheal www.nature.com/scientificreports www.nature.com/scientificreports/ diseases. The dose and feasibility outcomes from this study support development of rice bran based complementary weaning foods. Incorporating rice bran from local rice production and processing facilities should be a priority in subsequent trial designs. The development of sustainable, local and affordable food products for consumption by weaning infants and throughout childhood is particularly critical in regions where political instability continuously threaten the regions food and nutritional security.

Materials and Methods
Study design. A 6-month, prospective, randomized-controlled, phase one dose escalation dietary intervention was conducted in a cohort of weaning infants residing in León, Nicaragua and in the community of Dioro, Mali, West Africa. Nicaraguan infants were recruited from public health rosters provided by the local Health Ministry from Perla Maria and Sutiava health sectors, and Malian infants were recruited from the Dioro Community Health Center. To be eligible, male and female infants were screened between 4-5 months of age for health status, and then followed weekly for diarrhea episodes. Participants were excluded if they had experienced diarrhea or received antibiotic treatment within the previous month; had known allergies, or  www.nature.com/scientificreports www.nature.com/scientificreports/ immune-compromising conditions (e.g. parasitic or malarial infections); had previously been hospitalized; and/ or enrolled in a malnutrition treatment program. In Nicaragua, all eligible participants received 3 doses of the rotavirus vaccine per regular administration through the Immunization Program 82 . Rotavirus vaccination was not yet administrated to the Mali cohort. All Malian participants received vitamin A supplementation upon enrollment. Dietary intervention with rice bran started when infants were 6 months of age because WHO guidelines promote exclusive breastfeeding for the first six months of life 84,85 .
The required ethical board reviews and approvals were completed for Mali and Nicaragua as provided by the Internal Review Board (IRB) of the Colorado State University Research Integrity and the Compliance Review office (protocol ID# 14-5233 H Nicaragua, 15-5744 H Mali). In Mali, the Institut National de Recherche en Santé Publique (National Institute of Research in Public Health, FWA 00000892) approved the intervention, which occurred between October 2015 and May of 2016 and registered at clinicaltrial.gov as (NCT02557373) on 23 September 2015. Ethical review and approvals for the Nicaraguan intervention that occurred between March 2015 and October 2015 were provided by the IRBs of the Universidad Nacional Autónoma de Nicaragua -León, University of North Carolina at Chapel Hill, and Virginia Polytechnic Institute and State University and registered at clinicaltrial.gov on 26 November 2105 as (NCT02615886). Written informed consent was obtained from the infant's parent or responsible guardian prior to any data collection. Infant participants that met the eligibility criteria were randomized within each health sector, and sex (Nicaragua) and geographic location of household and sex (Mali) to either rice bran or control group (see Fig. S1 for enrollment details). Randomization was completed using sequential enrollment for each site independently. Participants were randomized by CSU, enrolled and assigned to groups by study coordinators in each country. The participants were not blinded. Complete study protocol is available online (http:// csu-cvmbs.colostate.edu/academics/erhs/Pages/elizabeth-ryan-lab-global-health.aspx).

Rice bran packaging for consumption. The United States Department of Agriculture-Agricultural
Research Service (USDA-ARS) Dale Bumpers National Rice Research Center provided rice bran that was polished from the U.S. variety, Calrose. Rice bran is prone to fat oxidation and heat-stabilization was performed to increase shelf-life by heating the bran at 100 degrees Celsius for five minutes to inactivate the lipase/lipoxygenase enzymes that cause rancidity 86 . The rice bran was sifted to remove any debris (rice husk, rice grain). Packaging of the rice bran was completed by Western Innovations, Inc. (Denver, CO) where 22 kg of rice bran was weighed into 1 g increments, separated into water-proof sachets, and heat-sealed to ensure the rice bran would be administered with accurate doses to infants.
Fourteen sachets (1 g/sachet) were filled into a 4″ × 3″ × 2″ box that was labeled for study participants and included nutrient information. These boxes were stored in a cool, dark, dry place until they were provided to study participants. nicaragua and Mali intervention. The study team (doctor, nurse and study coordinator) in Nicaragua and the community health workers (CHWs) in Mali provided a 2-week supply of rice bran at each routine home visit and instructed the participant's parent or guardian to add the daily amount of rice bran to the participant's food. At 6-7 months of age, participants in the rice bran group consumed 1 g of rice bran/day (1 sachet). Between the ages of 7-8 months, participants consumed 2 g of rice bran/day (2 sachets). At 8-10 months of age, participants consumed 3 g of rice bran/day (3 sachets). The amount increased to 4 g of rice bran/day (4 sachets) from 10-11 months, and 5 g of rice bran/day (5 sachets) from 11-12 months of age, respectively. The rice bran was added to appropriate weaning foods, such as rice cereal, yogurt, fruit and natural juices, vegetables, and soups. At the beginning of the intervention (six months of age), infant's parents or guardians were instructed and monitored daily for one week by study personnel to ensure that guardians knew how to administer and record the amount of rice bran consumed. Compliance to the rice bran intervention was calculated from records that had the dose/ amount of rice bran consumed circled in daily increments of none (0%), half (50%), or all (100%). The study team also collected any unused boxes or sachets during these visits. Participants in the control group did not receive any rice bran and there were no reports of brown rice intake during the 6-month study duration.
In Nicaragua, study personnel visited all infants weekly. In Mali, the CHWs visited each participant's household daily for the duration of the 6-month study to assess compliance and diarrhea episodes. If a participant had a diarrhea episode, the study team would collect a stool sample, and collect information that included the diarrhea onset date, how long the episode lasted, numbers of bowel movements within 24 hours, any associated signs and symptoms (e.g. nausea, vomiting, fever), if any other family members had diarrhea, and if any treatment was provided (e.g. antibiotics, rehydration).
The study team in Nicaragua collected data for control group participants at 6, 8, and 12 months old, and rice bran group every month. The anthropometric measures (weight and length) were collected via a portable stadiometer and weighing balance. Mali participants visited the Community Health Clinic every month. Length was measured in supine position using a reclining length-board. Length was collected to the nearest centimeter and weight to the nearest 0.1 kg. Anthropometric measures were calculated for LAZ, WAZ, and WLZ scores following the World Health Organizations (WHO) child growth standards using the WHO Anthro software (version 3.2.2) 87 .
Diapers were provided to all study participants. Stool was collected directly from soiled diapers. Freshly collected stool was diluted 20-fold and homogenized in a sterile pre-reduced anaerobic saline −0.1 M potassium phosphate buffer (pH 7.2) containing 20% glycerol (vol/vol). Four aliquot suspensions were prepared in 15 mL falcon tubes, transported on dry ice to the UNAN-León-Center of Infectious Diseases Laboratories (and liquid nitrogen in Bamako,Mali), immediately transferred to a −80 °C freezer, shipped in a liquid nitrogen chilled dry shipping dewar to Colorado State University, where they were relocated into a −80 °C freezer prior to analysis.
A study questionnaire was completed by the participant's caretaker (e.g. mother, father, or grandparent) to assess for duration of breastfeeding, types of and timing of introductions to complementary foods, as well as www.nature.com/scientificreports www.nature.com/scientificreports/ antibiotic use. The breastfeeding questions included whether or not the child was receiving breast milk, and/or had the child been receiving received formula. The complementary feeding history included a list of common Nicaraguan and Malian foods that are normally introduced to infants during weaning. Infant's parents or guardians recorded how often the infant consumed each of the eleven foods. The questionnaire also recorded if a participant had received treatment with antibiotics since the last visit, the reason for taking the antibiotic, the name of the antibiotic, as well as the length of time the participant had been taking the antibiotic. A household survey was also completed at the beginning of the trial to collect mother's education level, drinking water source, household flooring type, and animals present in the household. Analysis of breastfeeding and formula feeding patterns, complementary feeding practices, and associations with nutritional status at 6-months old (i.e. baseline) were previously reported for Nicaragua 88 . Monthly visits to the Mali community health clinic provided monitoring for malnutrition and severe adverse events; no adverse events were reported in the rice bran intervention group. There was one participant death reported in the control group (respiratory infection) and another withdrew to receive malnutrition treatment in the second month of the study. Diarrheal episodes were recorded, and a sample was collected in both countries using the same protocol.
Stool analysis for eeD markers. Stool biomarkers were selected to report gut inflammation and epithelial integrity as indicators of EED. These included neopterin (NEO), myeloperoxidase (MPO), calprotectin, (CAL) and alpha-1 antitrypsin (AAT) 89 . Suspended stool samples from 6, 8, and 12-month collections were centrifuged at 3,000 RPM to remove debris, following agitation, and the remaining supernatant was used for Enzyme-Linked-Immunosorbant-Assay (ELISA) determination of EED biomarker concentrations. Laboratory analysis protocols included in commercial kits were followed. Concentrations of CAL were determined at a 1:360 final dilution factor (Eagle Biosciences-Nashua, NH. Ref: CAL35-K01). Samples were diluted to 1:100 for determination of NEO concentrations (GenWay Biotech Inc-San Diego, CA, USA). MPO concentrations were determined at a 1:500 dilution factor (Immundiagnostik AG-Bensheim, Germany). Samples were diluted to 1:12,500 for determination of AAT concentrations (Immuchrom GMBH-Heppenheim, Germany), and dilution factors accounted for stool suspension ratios (20-fold). Final concentrations were determined from averages of three replicate assays and duplicate optical density readings, and interpolated using Graphpad 6.0 according to standards measured on each 96-well plate.
Stool microbiota analysis for nicaragua and Mali. The infant stool was collected at 6, 8 and 12 months of age from diapers and placed in a 1:19 ratio with Phosphate Buffered Saline + Glycerol solution. Diarrhea samples were collected using the same protocol. Suspended stool samples were vortexed before centrifuging at 3000 RPM to separate the stool debris. The remaining supernatant was used for Enzyme-Linked-Immunosorbant-Assay (ELISA) determination of EED biomarkers whereas; DNA was extracted for 16S microbial analysis from the stool pellet. DNA extraction was conducted using MoBio PowerSoil Kit (Reference number 12888, MoBio Laboratories Inc., Solana Beach, CA). PCR amplification of 390 bp amplicons was done in 50 µl reaction using Fischer Hot Start Master Mix and EMP standard protocols [43][44][45][46][47] . SPRI magnetic beads were used to purify DNA, and flourimetric quantification of Sybr Green tags was used to confirm adequate concentration of DNA. The pooled library was created with 50 ng DNA per sample and quantified using Kapa Kit (Kapa Biosystems). The pooled library was run on Illumina-MiSeq with 15% PhiX mock library to reduce discrepancies in read clustering, using the Illumina V2 500 cycle kit (2 × 250/250 paired-end reads).
Microbiota data processing and analysis. Sequence data were processed using mothur 90 version 1.39.5 and using a custom pipeline that provides an adjustment on the developers' standard operating procedure (SOP) for OTU calling and taxonomic classification of MiSeq data first presented in Kozich, et al. 91 . For alignment and classification within this SOP we used the SILVA database 92 version 128. Clustering, for OTU identification, was performed using VSEARCH using the distance based greedy clustering (DGC) option as implanted in mothur and utilizing 0.97 sequence similarity cutoff. We also used a cutoff of one read that was subtracted from all OTU read counts to guard against overestimation of sample richness. Rarefaction curves were generated using the package vegan 93 as implemented in R version 3.4.4 94 to assess diversity and suitability of depth of coverage per sample. The resulting OTU table was utilized in further data analyses as follows.
Exploring the data: Bar-graphs for relative abundance data based on the resulting OTU table were generated using the ggplot2 95 package in R. These plots were generated for the data at the genus and the family levels and meant to describe the microbial community structure per sampled infant and per time point under each of the treatment levels.
Data were normalized using cumulative sum scaling (CSS) 96 prior to beta diversity and log-fold change analyses. Nonmetric Multidimensional Scaling (NMDS) 97 was used on the OTU level to assess possible trends and clustering in the microbial community structure comparing the two countries, the treatment conditions and the two time points, using the vegan package and utilizing Bray-Curtis dissimilarity 97 . Data were separated per country and the metagenomeSeq. 98 package in R 94 was used to fit a zero inflated normal (ZIN) model to test for log-fold change differences between the rice bran treatment and control per age group. Benjamini and Hochberg's 99 false discovery rate (FDR) method was used to correct for multiple testing and compute the adjusted p-values used to determine significance of differences in the log-fold change of OTU abundance. www.nature.com/scientificreports www.nature.com/scientificreports/ to ten recovery standards were added prior to the first step in the extraction process for quality control purposes. Extraction was performed using 80% ice-cold methanol under vigorous shaking for 2 min (Glen Mills GenoGrinder 2000) followed by centrifugation to remove protein, dissociate small molecules bound to protein or trapped in the precipitated protein matrix. Each stool extract was divided into five fractions: two for analysis by two separate reverse phase UPLC-MS/MS methods with positive ion mode electrospray ionization, 1 for analysis by reverse phase UPLC-MS/MS methods with negative ion mode electrospray ionization, 1 for hydrophilic interaction liquid chromatography UPLC-MS/MS with negative ion mode electrospray ionization, and 1 sample for backup. All samples were placed briefly on Concentration Evaporator (TurboVap ® Zymark) to remove organic solvent. UPLC-MS/MS methods utilized a Waters ACQUITY ultra-performance liquid chromatography and a Thermo Scientific Q-Exactive high resolution/accurate mass spectrometer interfaced with a heated electrospray ionization (HESI-II) source and Orbitrap mass analyzer operated at 35,000 mass resolution. Raw data was extracted, peak-identified and processed for quality control using Metabolon's hardware and software. Statistical analysis. Statistical analyses for anthropometric measures (length, weight, LAZ, WAZ, and WLZ) and stool EED biomarkers were completed using SAS 9.4 (Cary, NC, USA). The sample size was calculated for achieving greater than 85% power and based on expected changes in selected stool metabolites following dietary rice bran consumption for one month 27 . Normality was evaluated by visual inspection. For anthropometric variables, two-sample t-tests were used to compare means for the 2 treatment groups (rice bran and control) separately at birth and 6 months (prior to start of treatment). A repeated measures analysis was performed for each response variable separately using SAS Proc Mixed. Specifically, treatment (rice bran or control) and age (6, 8 or 12 months), and treatment-age interaction were included in the model as fixed effects. The participant was included as a random effect to account for repeated measures. At each age, treatment groups were compared using contrasts of the model. A similar repeated measures analysis was used EED biomarkers, but log transformation was used to satisfy model assumptions. For stool metabolites, Welch's two-sample t-test was used to analyze statistical significance between groups' stool metabolites, after participating in the 6-month dietary trial. A p-value of ≤ 0.05 was used for statistical significance. An estimated false discovery rate (q-value) was calculated to account for the multiple comparisons across metabolites that are typical of metabolomic-based studies. one sentence summary. Dietary rice bran supplementation during infant weaning from 6-12 months of age improved growth outcomes and supported metabolism by the gut microbiota.

Data Availability
16S sequence data were submitted to the National Center for Biotechnology Information SRA under accession no. (SRP159269) and Bio-project (PRJNA488807).