Ethnic and diet-related differences in the healthy infant microbiome
Genome Medicinevolume 9, Article number: 32 (2017)
The infant gut is rapidly colonized by microorganisms soon after birth, and the composition of the microbiota is dynamic in the first year of life. Although a stable microbiome may not be established until 1 to 3 years after birth, the infant gut microbiota appears to be an important predictor of health outcomes in later life.
We obtained stool at one year of age from 173 white Caucasian and 182 South Asian infants from two Canadian birth cohorts to gain insight into how maternal and early infancy exposures influence the development of the gut microbiota. We investigated whether the infant gut microbiota differed by ethnicity (referring to groups of people who have certain racial, cultural, religious, or other traits in common) and by breastfeeding status, while accounting for variations in maternal and infant exposures (such as maternal antibiotic use, gestational diabetes, vegetarianism, infant milk diet, time of introduction of solid food, infant birth weight, and weight gain in the first year).
We demonstrate that ethnicity and infant feeding practices independently influence the infant gut microbiome at 1 year, and that ethnic differences can be mapped to alpha diversity as well as a higher abundance of lactic acid bacteria in South Asians and a higher abundance of genera within the order Clostridiales in white Caucasians.
The infant gut microbiome is influenced by ethnicity and breastfeeding in the first year of life. Ethnic differences in the gut microbiome may reflect maternal/infant dietary differences and whether these differences are associated with future cardiometabolic outcomes can only be determined after prospective follow-up.
The developing gastrointestinal microbiota in the first years of life is important for immune function, nutrient metabolism and protection from pathogens [1–3]. Microbial colonization of the infant gut proceeds through infancy and establishment of an adult-like microbiome is estimated to occur within the first 3 years . Identifying factors that shape the gut microbiome is currently an active area of research and early evidence suggests that host genetics  and early life exposures, including delivery method, antibiotics [6, 7], and diet, influence the infant gut microbiome [8, 9]. In addition to these established roles, the gut microbiota is emerging as a potentially important contributor to the development of non-communicable diseases (NCDs), having been associated with conditions such as obesity [10, 11], type 2 diabetes [12, 13], allergy and atopy , inflammatory bowel disease , and the development of colon cancer . The influence of the infant microbiome on the development of these conditions is of great clinical and economic interest as rates of NCDs in adults are increasing globally and by 2030 are predicted to account for 89% of all deaths in high income countries .
South Asians are people whose ancestors originate from the Indian subcontinent and they have among the highest rates of type 2 diabetes and premature cardiovascular disease (CVD) in the world. CVD risk factors, including adiposity, type 2 diabetes, and dyslipidemia, are higher among South Asians compared to white Caucasians of the same BMI . There is preliminary evidence that gut microbial composition in adults and children varies by age [4, 19], dietary intake [20, 21], ethnicity, geography [4, 22], and adoption of western lifestyles [19, 23]. Bacterial richness has been shown to increase with age and to be lower in residents of the United States compared with other populations . Core bacterial metabolic genes varied between these populations as well; however, the underlying reasons for ethnic and geographic differences in the microbiome have not been characterized. In this paper we investigate the associations of ethnicity and early life exposures with the gut microbiome among 1-year-old infants born and living in Canada while accounting for a diverse set of covariates that represent dietary differences as well as other exposures throughout infancy. This study explores the effect of ethnicity separately from region and provides a preliminary look at effects of ethnicity on the gut microbiota in early life.
Participants from two prospective Canadian birth cohorts were included in this gut microbiome substudy. The Canadian Healthy Infant Longitudinal Development study (CHILD) enrolled 3624 mainly white Caucasian mother–child pairs and most fathers from four Canadian centers (Vancouver, BC; Edmonton, AB; Winnipeg/Winkler-Morden, MB; and Toronto, ON) to investigate the root causes of allergy and asthma, including genetic and environmental triggers, and the ways in which they interact [24–26]. In this analysis, ethnicity refers to groups of people who have certain racial, cultural, religious, or other traits in common, whereas race refers to a person’s physical characteristics, such as bone structure, or skin, hair, or eye color [27, 28]. In the CHILD cohort, white Caucasian ancestry was confirmed by participants’ response to the question “To which ethnic or cultural group did your parents belong?” The South Asian Birth Cohort (START-Canada) enrolled 1012 South Asian mother–child pairs from the Brampton and Peel Region of Ontario to investigate the influence of diverse environmental exposures and genetics on early life adiposity, growth trajectory, and cardiometabolic factors . South Asian ethnicity was verified by the mother’s self-report of her and the father’s, and their parents’, ancestral origin being from India, Pakistan, Sri Lanka, or Bangladesh.
Harmonization of clinical data across cohorts was done by extracting them with the same definitions, where possible. When questions were not identical, we worked to extract the data from each cohort in such a way as to satisfy the same definition. Gestational diabetes mellitus was defined as having diabetes on the birth chart but no diabetes prior to pregnancy. A child was considered to have had formula in the first year if formula use was recorded at any time in the first year (from several questionnaires). In both cohorts, timing of infant weighing at 1 year was typically performed on the same day as 1-year stool collection (r2 > 0.93; median = 0 days; 95% confidence interval 0 to 2 days).
In this gut microbiome substudy, 1-year fecal samples from 173 white Caucasian infants in CHILD and 182 South Asian infants in START were used for the main analysis. An additional 77 samples from the CHILD cohort, from infants who are not white Caucasian, were used to explore trends found in the main analysis. For both cohorts, the collection of 1-year fecal samples was scheduled with the mother in advance. Stool collection was taken from a regular diaper in START and a specially lined diaper in CHILD . Mothers were instructed to record the time and date of the stool sample and place it in a sterile bag in the refrigerator for their scheduled appointment with the research nurse. Upon arrival, the nurse used depyrogenized stainless steel spatulas to divide the sample between four pre-labeled cryovials. The cryovials were then transported to the lab in a cooler, weighed, and stored at −80 °C or liquid nitrogen. START samples were stored at 4 °C for 2–4 h prior to freezing whereas CHILD samples were stored at 4 °C for an average of 14 ± 12 h.
DNA extraction, 16S rRNA gene sequencing, and analysis
DNA was extracted with a custom DNA extraction protocol described in . Briefly, 100–200 mg of stool was added to 2.8 mm and 0.1 mm glass beads (MoBio Laboratories Inc., Carlsbad, CA, USA) along with 800 μl of 200 mM sodium phosphate monobasic (pH 8) and 100 μl guanidinium thiocyanate EDTA N-lauroylsarkosine buffer (50.8 mM guanidine thiocyanate, 100 mM ethylenediaminetetraacetic acid, and 34 mM N-lauroylsarcosine). These were then homogenized in a PowerLyzer 24 Bench Top Homogenizer (MoBio Laboratories Inc.) for 3 min at 3000 RPM. Next, two enzymatic lysis steps were performed. First, the sample was incubated with 50 μl of 100 mg/ml lysozyme, 500 U mutanolysin, and 10 μl of 10 mg/ml RNase for 1 h at 37 °C. Next, the sample was incubated with 25 μl 25% sodium dodecyl sulphate, 25 μl of 20 mg/ml Proteinase K, and 62.5 μl of 5 M NaCl at 65 °C for 1 h. Next, debris was pelleted in a tabletop centrifuge at maximum speed for 5 min and the supernatant added to 900 μl of phenol:chloroform:isoamyl alcohol (25:24:1). The sample was then vortexed and centrifuged at maximum speed in a tabletop centrifuge for 10 min. The aqueous phase was removed and the sample run through the Clean and Concentrator-25 column (Zymo Research, Irvine, CA, USA) according to kit directions except for elution, which was done with 50 μl of ultrapure water and allowed to sit for 5 min before elution. The DNA was quantified using a Nanodrop 2000c Spectrophotometer . Amplification of the bacterial 16S rRNA gene v3 region (150 bp) tags was performed as previously described  with the following changes: 5 pmol of primer, 200 μM of each dNTP, 1.5 mM MgCl2, 2 μl of 10 mg/ml bovine serum albumin, and 1.25 U Taq polymerase (Life Technologies, Carlsbad, CA, USA) were used in a 50 μl reaction volume. The PCR program used was as follows: 94 °C for 2 min followed by 30 cycles of 94 °C for 30 s, 50 °C for 30 s, and 72 °C for 30 s, then a final extension step at 72 °C for 10 min. DNA extraction and PCR amplification of 16S rRNA gene v3 libraries were found to be reproducible using a set of five samples from each cohort (total of ten samples) that were extracted in triplicate (29 extractions since one extraction failed) and a subset of three extractions from each cohort amplified in triplicate for a total of 41 datasets (Additional file 1: Figure S1).
Illumina libraries were sequenced in the McMaster Genomics Facility with 250-bp sequencing in the forward and reverse directions on the Illumina MiSeq instrument. Custom, in-house Perl scripts were used to process Illumina sequences as previously described . Briefly, after sequence trimming and alignment, operational taxonomic units (OTU) were clustered using AbundantOTU+  with a threshold of 97%. Chimera checking was not done since we have shown that amplification of the short V3 region of the 16S rRNA gene leads to very few genuine chimeric sequences . Taxonomy for the representative sequence of each OTU was assigned using the Ribosomal Database Project classifier  with a minimum confidence cutoff of 0.8 against the Greengenes (2013 release) reference database . All OTUs classified as “Root:Other” (comprising 0.03% of the total reads sequenced) were then excluded as was one sample with <500 sequenced reads; however, singleton OTUs were not excluded. This resulted in a total of 41.4 million reads with a minimum of 2.0 × 103, maximum of 4.3 × 105, and a median of 9.0 × 104 reads per sample.
Bacterial community richness and diversity (alpha diversity) were calculated using the estimated species richness and Shannon diversity functions with the vegan package in R , using OTU abundances. Differences between bacterial communities in each sample (beta diversity) were quantified using the Bray–Curtis dissimilarity measure on relative abundance values of all bacterial genera and principal coordinate analysis was also done using the vegan package or the phyloseq package  in R.
Simple linear regression was used to determine the effect of ethnicity and breastfeeding on alpha diversity estimates. Permutational multivariate analysis of variance on Bray–Curtis dissimilarities of genus level relative abundances, done with the adonis function from the vegan package in R , was used to examine bacterial community differences associated with ethnicity after adjustment for potential covariates of ethnicity–microbiome associations.
Candidate covariates in the multivariable model were informed by the existing literature and assessed formally in univariable models against microbiome diversity (i.e., years mother lived in Canada, breastfeeding at time of collection, time since weaning, formula and cow’s milk use in the first year, time of introduction of solid foods, infant weight gain in the first year, birth weight, infant age at stool collection, and mode of delivery, gestational diabetes, mother’s antibiotic use during pregnancy and labor, and mother’s vegetarian status). Next, the candidate variables chosen above were used to separately predict dissimilarities with the same method as above. Those with p < 0.10 were subjected to a forward stepwise procedure. We then added the most significant covariates into the model in order of the proportion of variance explained, and stopped when the next most significant covariate was above the 0.05 threshold.
The association between genus level abundances and ethnicity and/or breastfeeding was determined through a multivariate algorithm adjusting for significant covariates performed with the Maaslin package in R [39, 40]. Briefly, covariates found to be significant (p < 0.05) predictors of the microbiome (described above) were included into a multivariate boosted, additive general linear model between covariate data and bacterial genus level abundances. P values were adjusted for multiple testing with the false discovery rate, reported as q values, and q < 0.05 was considered significant. Genera with a coefficient of variation >0.001 were included in Additional file 2: Table S1.
Table 1 shows the baseline demographic and anthropometric characteristics of the mothers and infants selected from CHILD (white Caucasians only) and START. Briefly, South Asian mothers lived in Canada for an average of 8 years versus a lifetime for white Caucasian mothers. Furthermore, South Asian mothers were younger, more likely to be vegetarian (34% versus 2%, p < 0.001), and to be diagnosed with gestational diabetes during pregnancy (14% versus 4%, p < 0.001) compared to white Caucasian mothers. There were no significant differences in the rates of Caesarian section between ethnic groups (18% in South Asian versus 15% in white Caucasian); however, white Caucasian mothers were more likely to receive antibiotics during pregnancy (8% versus 0.5%, p < 0.001) and South Asian mothers were more likely to receive antibiotics during labor (43% versus 34%, p < 0.05). South Asian infants were born earlier (39.1 weeks versus 39.5 weeks, p < 0.05), had lower birth weight (3.3 kg versus 3.5 kg, p < 0.001), and gained more weight in the first year of life (7.1 kg gained versus 6.4 kg gained, p < 0.001) than did white Caucasian infants. While both white Caucasian and South Asian mothers reported that they breastfed their infants at some point during the first year (97.1% versus 94.4%), a greater proportion of South Asian infants were still breastfeeding at the time of 1-year stool sample collection (43% versus 32%, p < 0.05). Additionally, there was more formula use during the first year (77% versus 65%, p < 0.001) and earlier introduction of solid food among South Asians (88% versus 50% from 3 to 6 months, 9.4% versus 40% from 6 to 9 months, p < 0.001). We suspect that more South Asian infant diets were vegetarian on account of the greater proportion of their mothers who identified as vegetarian (34% versus 2%, p < 0.001). Furthermore, there was no difference in age at time of stool collection (p = 0.39) between white Caucasians (12.3 ± 1.71 months) and South Asians (12.4 ± 1.69 months; Table 1).
Abundance of microorganisms within all samples
The v3 region of 16S rRNA genes was profiled from 355 participant stool samples collected at 1 year of age, 173 white Caucasians from the CHILD cohort and 182 from the START cohort. The range of alpha diversity estimates for each ethnicity separated by current breastfeeding at the time of sampling is illustrated in Fig. 1. Using simple linear regression, species richness estimates were found to be significantly affected by ethnicity after taking into account breastfeeding at time of collection (p < 0.05). Shannon diversity was significantly affected by ethnicity, taking into account breastfeeding at time of collection (p < 0.001), and likewise breastfeeding at time of collection within each ethnicity significantly affected Shannon diversity (p < 0.05). Further, when START samples, all collected within the Brampton/Peel region of Ontario Canada, were compared with each study center within the CHILD cohort (Vancouver, Edmonton, Toronto, and Winnipeg/Winkler-Morden) only Winnipeg/Winkler-Morden, MB had significantly lower species richness estimates (p < 0.05; Additional file 1: Figure S2). Although there was variability in Shannon diversity estimates across sample sites for CHILD, all sites were found to have significantly lower diversity than the START samples (p < 0.05; Additional file 1: Figure S2), while accounting for current breastfeeding. By including sample sites into the regression model the effect of current breastfeeding on Shannon diversity was no longer significant (p = 0.054).
Differences in the relative abundance of the dominant bacterial genera are presented in Additional file 1: Figure S3, broken down by ethnic group and breastfeeding status. Heterogeneity of samples can be seen in Additional file 1: Figure S3 as well as differences in genus level microbobial profiles between ethnic groups and breastfeeding status, differences that are explored in detail below.
Principal coordinate analysis of Bray–Curtis dissimilarities illustrates between-community differences in samples from white Caucasians and South Asian infants. Variation in the gut microbiome across geography has been observed in studies involving adults ; however, in our study the effect of ethnicity was larger than the effect of geographic location (Fig. 2) shown as the separation of the centroid for samples from South Asians from the centroids of samples from white Caucasians from all study centers. Also evident from the principal coordinate analysis, breastfeeding at time of collection affected the gut microbial profiles, although when stratified by currently breastfed and not currently breastfed infants, the strong effect of ethnicity persisted (Fig. 2). Several studies have found the infant gut microbiome to vary between infants born by Caesarean section and those born vaginally with the effect diminishing with age. Here, delivery method was not found to be a significant predictor of the structure of the gut microbiome in 1-year-old infants (Additional file 1: Figure S4). This may be because differences were no longer strong enough to be detected or because members of the phylum Bacteroidetes, often missing from the gut microbiome in Caesarean section delivered infants, were not abundant in our vaginally born infants (Additional file 1: Figure S2).
Association between ethnicity, milk diet, and solid food diet
In addition to ethnicity, 13 potential covariates were also associated with the microbiome in univariable regression analysis. These included mother’s years living in Canada, infant age, breastfeeding status at time of collection, time since weaning, vegetarian status, timing of introduction of solid foods, birth weight, infant weight gain in the first year, antibiotics during pregnancy, antibiotics during labor, formula use in the first year, formula use at collection, and cow’s milk in the first year (all p < 0.10; Table 2). We entered this set of covariates into a forward stepwise regression model to determine which factors remained significant and independently influenced the gut microbiome. Only ethnicity (p < 0.001), breastfeeding status (p < 0.001), infant age at stool collection (p < 0.01), and weight gain in the first year (p < 0.01) remained independently associated with the gut microbiome as a whole.
There was no statistically significant multiplicative interaction between ethnicity and breastfeeding (p = 0.23). Nevertheless, we acknowledge that such tests may be underpowered, and thus the results were also stratified by ethnicity and breastfeeding status in order to examine trends. Forward stepwise regression was conducted within white Caucasians and separately within South Asians (Table 3). This revealed that breastfeeding (p < 0.01) and infant age (p < 0.05) were independently associated with differences in the microbiome within each ethnic group, while antibiotic use during labor (p < 0.05) and weight gain in the first year (p < 0.05) remained independently associated with differences in the microbiome only in white Caucasians. Forward stepwise regression was also conducted separately within infants breastfed and not breastfed at the time of collection (Table 4), which indicated that ethnicity (p < 0.01) and the infant age (p < 0.05) remained independently associated with differences in the gut microbiome in both groups.
Differentially abundant genera within each group
Difference in the relative abundance of the dominant bacterial genera is presented as a taxa bar chart in Additional file 1: Figure S2, broken down by ethnic group and breastfeeding status. The relative abundance of individual bacterial genera was assessed for association with ethnicity and breastfeeding while accounting for infant age and weight gain in the first year. These covariates, which had survived the stepwise regression on the entire community, were included in the multivariate algorithm in order to strike a balance between overfitting the model and identifying the most comprehensive list of predictors. Taxa significantly associated with ethnicity, breastfeeding at time of collection, infant weight gain in the first year, and infant age (q value <0.05) are listed in Additional file 2: Table S1 and their abundance is illustrated in Fig. 3.
South Asians had higher abundances of several genera within the Actinobacteria (Bifidobacterium, Collinsella, Actinomyces, Atopobium) and of three unclassified genera compared to white Caucasians. Genera within the phylum Firmicutes within two distinct taxonomic groups were associated with ethnicity. Genera such as Streptococcus, Enterococcus, and Lactobacillus (class Bacilli, order Lactobacillales) were more abundant within South Asians whereas genera such as Blautia, Pseudobutyrivibrio, Ruminococcus, and Oscillospira (order Clostridiales) were more abundant in white Caucasians. The most differentially abundant genus were unclassified members of the Lachnospiraceae which were higher in white Caucasians. In order to investigate whether these differences were specific to each cohort or were indicative of true ethnic differences, five genera significantly associated with either white Caucasians or South Asians were plotted among the small number of South Asians recruited within the CHILD cohort (n = 6 that were not used for the previous microbiome analysis). Despite the small number available, the same trends were seen for the five genera plotted (Additional file 1: Figure S5).
Not surprisingly, breastfeeding status at the time of sample collection was strongly associated with the abundance of the genera Bifdobacterium (phylum Actinobacteria; Fig. 3). Several genera within the phylum Firmicutes were associated with breastfeeding at the time of collection; some were more abundant (Veillonella, Megasphaera, and Dialister) and others were less abundant (Blautia, unclassified Lachospiraceae, Clostridium, Ruminococcus, Coprobacillus, Lactococcus, as well as several unclassified genera within the Clostridiales and Erysipelotrichales).
Our results demonstrate that the gut microbiome of infants is influenced by ethnicity, infant age, weight gain, and breastfeeding. The gut microbiome has been proposed to influence the progression of chronic diseases and has been associated with adverse health outcomes . Development of the microbiome within the first years of life may influence long-term health, and can be affected by perinatal, genetic, and dietary factors, including solid foods and milk diet.
The distribution of a number of maternal and infant parameters differed between white Caucasian and South Asians (i.e., vegetarian status, gestational diabetes mellitus prevalence, timing of introduction of solid foods, antibiotic use during pregnancy, mode of delivery, etc.) and thus seemed likely candidates to explain the microbiome differences by ethnicity. However, when these variables were added as independent predictors of the gut microbiome composition in the multivariable model, none except breastfeeding status at the time of sampling, infant age, and weight gain in the first year improved the fit of the model (ethnicity R2 = 0.084 versus R2 = 0.082 with all additional variables; breastfeeding status R2 = 0.040 versus R2 = 0.032 with all additional variables). This suggests that these variables were largely captured by the higher order variables of interest, i.e., ethnicity and breastfeeding. Next, after taking into account these significant predictors (breastfeeding status, infant age at 1-year stool, and weight gain in the first year of life) we found that groups of bacterial genera which are phylogenetically distinct (i.e., within the order Lactobacillales versus Clostridales) were present at different abundances within each ethnic group. This suggests that different metabolic strategies are at work within the gut microbiome of South Asian and white Caucasian infants. Additionally, these bacterial taxa are good candidates to predict diet-related influences on the microbiome, microbial influences on host metabolism, and bacterial stimulation of the host immune system .
Several members of the lactic acid bacteria (LAB), specifically Bifidobacterium, Lactococcus, Streptococcus, and Enterococcus, were more abundant within South Asians after taking into account breastfeeding status at the time of collection, infant age, and weight gain in the first year. LAB break down mainly carbohydrates that are not absorbed by the host to produce acetate and lactate, both of which are used as energy sources by other microbial groups [43, 44]. Also the abundance of members of the Atopobium cluster of Actinobacteria (i.e., genera such as Collinsella and Atopobium) was higher in South Asians. This group of bacteria are saccharolytic (i.e., they break down small sugars)  and have been seen to decrease in abundance in the microbiome of individuals with a diet rich in whole grains . These genera have also been associated with higher levels of low-density lipoprotein in humans  and, along with other members of the Actinobacteria, have been associated with high hepatic levels of triglycerides and low hepatic levels of glycogen and glucose in mice . It is of interest to note that these observations are based on v3 16S rRNA gene data. Several studies of the infant gut microbiome, which employ amplification and sequencing of other variable regions of the same gene often report very low levels of Actinobacteria [6, 9, 49]. Members of this phylum, such as the Bifidobacteria, have been shown to dominate the infant gut microbiome [4, 50, 51], suggesting a possible primer bias against this group.
In contrast, white Caucasians showed higher abundances of members of the Firmicutes from the order Clostridiales, which have been shown to be increased in response to diets rich in animal protein  and high in fat . Products of bacterial fermentation of acetate and lactate, mentioned above, as well as non-digestible fiber and oligosaccharides by members of the Clostridiales seen here (Ruminococcus, Lachnospiraceae, and Oscillospira) include short chain fatty acids like butyrate, which is used by host cells as an energy source and can signal increased barrier function . Though also proposed to be chemoprotective, the relationship between luminal butyrate exposure and colorectal cancer in humans has been examined only indirectly in case-control studies . Nevertheless, these findings suggest different metabolic processes and immune stimuli at work within the South Asian and white Caucasian infant gastrointestinal tract, some of which may be explained by their heterogeneous diets.
When switching from a milk-based diet to a solid food diet, prior studies have shown a decrease in the abundance of Bifidobacterium along with an increase in members of the Firmicutes (such as Clostridium sp.) and Bacteroidetes [12, 20]. One study suggests that it is the cessation of breastfeeding that is required for maturation of the gut microbiota to occur with a decrease in Bifidobacterium and an increase in members of the Clostridiales only occurring after weaning . As expected, after adjustment for ethnicity, infant age, and weight gain in the first year, Bifidobacterium and Lactobacillus were significantly associated with breastfeeding. Additionally, an increase in the abundance of several genera within the phylum Firmicutes were associated with not being breastfed at the time of sampling.
Bifidobacteria, along with the LAB, are known to be abundant members of the microbiome of breastfeeding infants , whereas genera within the order Clostridiales are known to be more abundant within the gut of adults . Here bacterial profiles indicative of a breast milk diet were common among South Asians even those that were not breastfeeding at the time of collection, suggesting that these infants retain more of a breastfeeding microbiome than do white Caucasians of the same age. The reasons for this are unclear; however, dietary differences may be contributing. Our data show that equal proportions of infants in both groups were breastfed in the first year but does not capture breastfeeding frequency. It also shows that there was a much higher rate of formula use and an earlier introduction to solid food within South Asian than within white Caucasian infants. Because self-reported vegetarianism was more frequent in South Asians, it is possible that meat consumption hastens, or non-meat diets delay, changes induced within the infant gut microbiome during the switch to a solid food diet. It is important to note, however, that to our knowledge an analysis of the adult South Asian microbiome has not been reported, nor has a description of the maturation of the South Asian infant microbiome toward an adult-like composition; thus, our data must be interpreted within the context of the study, i.e., of South Asian infants born in Canada who consume a South Asian diet.
The underlying construct of “Ethnicity” brings together several biological and cultural factors, and it can be characterized using a number of different parameters (e.g., dietary habits, ancestral country of origin, etc.). In our multivariate model, ethnicity and breast feeding status remained independent and significant predictors of differences in the overall microbial communities (beta diversity), whereas vegetarian diet did not, which implies that the impact of ethnicity which incorporates some unique dietary patterns is not wholly explained by these dietary differences, as it also reflects other differences between the groups. After ensuring that these additional factors were potentially accounted for (i.e., years living in Canada, antibiotic use, timing of solid food introduction, etc.) we observed that breastfeeding, infant age, and weight gain in the first year significantly influenced the infant gut microbiome.
Strengths of our study include its relatively large size of nearly 200 infants from each of two different ethnic groups who have diverse dietary intakes; the availability of stool samples collected at similar times using similar methods; the high quality deep sequencing of the 16S rRNA gene for bacterial identification; a reliability analysis to demonstrate reproducibility of our methods; and detailed measurement of maternal and infant covariates. Limitations include incomplete data on maternal weight gain during and prior to pregnancy, which limits our ability to assess the influence of this important covariate on the infant gut microbiome; ethnicity in this study refers to the group a person self-identifies with and reflects a mix of cultural factors, including language, diet, religion, and ancestry—thus, ethnicity is a multidimensional construct which includes some within-group heterogeneity, and differences attributable to ethnicity may reflect a broad range of factors which are not purely biological; and the lack of a direct measure of infant dietary intake beyond feeding type at the time of stool collection.
The infant gut microbiome is influenced by ethnicity and breastfeeding in the first year of life. Ethnic differences in the gut microbiome may reflect maternal/infant dietary differences and whether these differences are associated with future cardiometabolic outcomes can only be determined after prospective follow-up.
Canadian Healthy Infant Longitudinal Development study
Lactic acid bacteria
Operational taxonomic unit
South Asian Birth Cohort
Falk PG, Hooper LV, Midtvedt T, Gordon JI. Creating and maintaining the gastrointestinal ecosystem: what we know and need to know from gnotobiology. Microbiol Mol Biol Rev. 1998;62:1157–70.
Guarner F, Malagelada J-R. Gut flora in health and disease. Lancet. 2003;361:512–9.
Newburg DS, Walker WA. Protection of the neonate by the innate immune system of developing gut and of human milk. Pediatr Res. 2007;61:2–8.
Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486:222–7.
Li Y, Oosting M, Deelen P, Ricaño-Ponce I, Smeekens S, Jaeger M, et al. Inter-individual variability and genetic influences on cytokine responses to bacteria and fungi. Nat Med. 2016;22:952–60.
Bokulich NA, Chung J, Battaglia T, Henderson N, Jay M, Li H, et al. Antibiotics, birth mode, and diet shape microbiome maturation during early life. Sci Transl Med. 2016;8:343ra82.
Azad MB, Konya T, Maughan H, Guttman DS, Field CJ, Chari RS, et al. Gut microbiota of healthy Canadian infants: profiles by mode of delivery and infant diet at 4 months. CMAJ. 2013;185:385–94.
Goodrich JK, Waters JL, Poole AC, Sutter JL, Koren O, Blekhman R, et al. Human genetics shape the gut microbiome. Cell. 2014;159:789–99.
Bäckhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, et al. Dynamics and stabilization of the human gut microbiome during the first year of life. Cell Host Microbe. 2015;17:852.
Ley RE, Bäckhed F, Turnbaugh P, Lozupone CA, Knight RD, Gordon JI. Obesity alters gut microbial ecology. Proc Natl Acad Sci U S A. 2005;102:11070–5.
Finucane MM, Sharpton TJ, Laurent TJ, Pollard KS. A taxonomic signature of obesity in the microbiome? Getting to the guts of the matter. PLoS One. 2014;9:e84689.
Laursen MF, Andersen LBB, Michaelsen KF, Mølgaard C, Trolle E, Bahl MI, et al. Infant gut microbiota development is driven by transition to family foods independent of maternal obesity. mSphere. 2016;1:e00069–15.
Qin J, Li Y, Cai Z, Li S, Zhu J, Zhang F, et al. A metagenome-wide association study of gut microbiota in type 2 diabetes. Nature. 2012;490:55–60.
Penders J, Gerhold K, Stobberingh EE, Thijs C, Zimmermann K, Lau S, et al. Establishment of the intestinal microbiota and its role for atopic dermatitis in early childhood. J Allergy Clin Immunol. 2013;132:601–607.e8.
Knights D, Lassen KG, Xavier RJ. Advances in inflammatory bowel disease pathogenesis: linking host genetics and the microbiome. Gut. 2013;62:1505–10.
Irrazábal T, Belcheva A, Girardin SE, Martin A, Philpott DJ. The multifaceted role of the intestinal microbiota in colon cancer. Mol Cell. 2014;54:309–20.
Nikolic IA, Stanciole AE, Zaydman M. Chronic emergency: why NCDs matter. Health, Nutrition, and Population Discussion Paper. 2011. http://www.ghd-net.org/sites/default/files/ChronicEmergencyWhyNCDsMatter.pdf. Accessed June 2016.
Rana A, de Souza RJ, Kandasamy S, Lear SA, Anand SS. Cardiovascular risk among South Asians living in Canada: a systematic review and meta-analysis. CMAJ Open. 2014;2:E183–91.
Clemente JC, Pehrsson EC, Blaser MJ, Sandhu K, Gao Z, Wang B, et al. The microbiome of uncontacted Amerindians. Sci Adv. 2015;1. doi:10.1126/sciadv.1500183
Fallani M, Amarri S, Uusijarvi A, Adam R, Khanna S, Aguilera M, et al. Determinants of the human infant intestinal microbiota after the introduction of first complementary foods in infant samples from five European centres. Microbiology. 2011;157:1385–92.
Zhang J, Guo Z, Xue Z, Sun Z, Zhang M, Wang L, et al. A phylo-functional core of gut microbiota in healthy young Chinese cohorts across lifestyles, geography and ethnicities. ISME J. 2015;9:1979–90.
Martínez I, Stegen JC, Maldonado-Gómez MX, Eren AM, Siba PM, Greenhill AR, et al. The gut microbiota of rural papua new guineans: composition, diversity patterns, and ecological processes. Cell Rep. 2015;11:527–38.
Schnorr SL. The diverse microbiome of the hunter-gatherer. Nature. 2015;518:S14–5.
Moraes TJ, Lefebvre DL, Chooniedass R, Becker AB, Brook JR, Denburg J, et al. The Canadian healthy infant longitudinal development birth cohort study: biological samples and biobanking. Paediatr Perinat Epidemiol. 2015;29:84–92.
Takaro TK, Scott JA, Allen RW, Anand SS, Becker AB, Befus AD, et al. The Canadian Healthy Infant Longitudinal Development (CHILD) birth cohort study: assessment of environmental exposures. J Expo Sci Environ Epidemiol. 2015;25:580–92.
Subbarao P, Anand SS, Becker AB, Befus AD, Brauer M, Brook JR, et al. The Canadian Healthy Infant Longitudinal Development (CHILD) Study: examining developmental origins of allergy and asthma. Thorax. 2015;70:998–1000.
Anand SS. Using ethnicity as a classification variable in health research: perpetuating the myth of biological determinism, serving socio-political agendas, or making valuable contributions to medical sciences? Ethn Health. 1999;4:241–4.
de Souza RJ, Anand SS. Cardiovascular disease in Asian Americans: unmasking heterogeneity. J Am Coll Cardiol. 2014;64:2495–7.
Anand SS, Vasudevan A, Gupta M, Morrison K, Kurpad A, Teo KK, et al. Rationale and design of South Asian birth cohort (START): a Canada-India collaborative study. BMC Public Health. 2013;13:79.
Stearns JC, Davidson CJ, McKeon S, Whelan FJ, Fontes ME, Schryvers AB, et al. Culture and molecular-based profiles show shifts in bacterial communities of the upper respiratory tract that occur with age. ISMEJ. 2015;9:1246–59.
Bartram AK, Lynch MDJ, Stearns JC, Moreno-Hagelsieb G, Neufeld JD. Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end illumina reads. Appl Environ Microbiol. 2011;77:3846–52.
Whelan FJ, Verschoor CP, Stearns JC, Rossi L, Johnstone J, Surette MG, et al. The loss of topography in the microbial communities of the upper respiratory tract in the elderly. Ann Am Thorac Soc. 2014;11:513–21.
Ye Y. Identification and quantification of abundant species from pyrosequences of 16S rRNA by consensus alignment. 2010 IEEE International Conference on Bioinformatics and Biomedicine (BIBM). New York: Institute of Electrical and Electronics Engineers (IEEE); 2010. pp. 153–7. doi:10.1109/BIBM.2010.5706555.
Stearns JC, Lynch MDL, Senadheera DB, Tenenbaum HC, Goldberg MB, Cvitkovitch DG, et al. Bacterial biogeography of the human digestive tract. Sci Rep. 2011;1:1–9.
Wang Q, George MG, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.
DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.
Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, et al. vegan: Community Ecology Package. 2015. http://CRAN.R-project.org/package = vegan.
McMurdie PJ, Holmes S. Phyloseq: An R Package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.
Morgan XC, Tickle TL, Sokol H, Gevers D, Devaney KL, Ward DV, et al. Dysfunction of the intestinal microbiome in inflammatory bowel disease and treatment. Genome Biol. 2012;13:R79.
Tickle T, Huttenhower C. Maaslin. Multivariate statistical framework that finds associations between clinical metadata and microbial community abundance or function. 2014. https://bitbucket.org/biobakery/maaslin
Dugas LR, Fuller M, Gilbert J, Layden BT. The obese gut microbiome across the epidemiologic transition. Emerg Themes Epidemiol. 2016;13:2.
Aron-Wisnewsky J, Clément K. The gut microbiome, diet, and links to cardiometabolic and chronic disorders. Nat Rev Nephrol. 2016;12:169–81.
Flint HJ, Duncan SH, Scott KP, Louis P. Links between diet, gut microbiota composition and gut metabolism. Proc Nutr Soc. 2015;74:13–22.
Duncan SH, Louis P, Flint HJ. Lactate-utilizing bacteria, isolated from human feces, that produce butyrate as a major fermentation product. Appl Environ Microbiol. 2004;70:5810–7.
Thorasin T, Hoyles L, McCartney AL. Dynamics and diversity of the “Atopobium cluster” in the human faecal microbiota, and phenotypic characterization of “Atopobium cluster” isolates. Microbiology. 2015;161:565–79.
Martínez I, Lattimer JM, Hubach KL, Case JA, Yang J, Weber CG, et al. Gut microbiome composition is linked to whole grain-induced immunological improvements. ISME J. 2013;7:269–80.
Lahti L, Salonen A, Kekkonen RA, Salojärvi J, Jalanka-Tuovinen J, Palva A, et al. Associations between the human intestinal microbiota, Lactobacillus rhamnosus GG and serum lipids indicated by integrated analysis of high-throughput profiling data. PeerJ. 2013;1:e32.
Claus SP, Ellero SL, Berger B, Krause L, Bruttin A, Molina J, et al. Colonization-induced host-gut microbial metabolic interaction. MBio. 2011;2:e00271–10.
Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO. Development of the human infant intestinal microbiota. PLoS Biol. 2007;5:e177.
Penders J, Thijs C, Vink C, Stelma FF, Snijders B, Kummeling I, et al. Factors influencing the composition of the intestinal microbiota in early infancy. Pediatrics. 2006;118:511–21.
Turroni F, Peano C, Pass DA, Foroni E, Severgnini M, Claesson MJ, et al. Diversity of bifidobacteria within the infant gut microbiota. PLoS One. 2012;7:e36957.
Wu GD, Chen J, Hoffmann C, Bittinger K, Chen Y-Y, Keilbaugh SA, et al. Linking long-term dietary patterns with gut microbial enterotypes. Science. 2011;334:105–8.
de La Serre CB, Ellis CL, Lee J, Hartman AL, Rutledge JC, Raybould HE. Propensity to high-fat diet-induced obesity in rats is associated with changes in the gut microbiota and gut inflammation. Am J Physiol Gastrointest Liver Physiol. 2010;299:G440–8.
Sengupta S, Muir JG, Gibson PR. Does butyrate protect from colorectal cancer? J Gastroenterol Hepatol. 2006;21:209–18.
Solís G, de Los Reyes-Gavilan CG, Fernández N, Margolles A, Gueimonde M. Establishment and development of lactic acid bacteria and bifidobacteria microbiota in breast-milk and the infant gut. Anaerobe. 2010;16:307–10.
Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature. 2012;489:220–30.
We are grateful to all the families who took part in this study, and the whole START and CHILD teams, which includes interviewers, nurses, computer and laboratory technicians, clerical workers, research scientists, volunteers, managers, and receptionists. In particular, we would like to thank Dipika Desai, Nora Abdalla, and Diana LeFebrevre for help with coordination of the studies as wells as Laura Rossi for her contributions to the technical aspects of this work.
This work was funded by a grant from the Canadian Institutes of Health Research (grant number FH6 129924). The CHILD Study was primarily funded by CIHR and the Allergy, Genes and Environment (AllerGen) Network of Centres of Excellence. The START study was part of a bilateral ICMR/CIHR funded program (grant number INC-109205) and from the Heart and Stroke Foundation (grant number NA7283). JCS holds the Endowed Farncombe Family Chair in Microbial Ecology and Bioinformatics at McMaster University. SSA holds a Canada Research Chair in Ethnicity and Cardiovascular Disease and the Michael G. DeGroote Heart and Stroke Foundation of Canada Chair in Population Health. MRS holds the AstraZeneca Endowed Chair in Respiratory Epidemiology. MGS holds a Canada Research Chair in Interdisciplinary Microbiome Research. MAZ holds a CIHR RCT Fellowship grant (MTP201410, MAZ).
Availability of data and materials
The datasets generated and/or analyzed during the current study are not publicly available, since the CHILD and START studies are bound by consent and cannot provide identifiable information to an outside group, but are available from the corresponding author on reasonable request.
JCS analyzed and interpreted the microbiome data. NCC, MAZ, RDS, SSA, and JB harmonized and assisted in the analysis of the subject data across cohorts and MF processed stool samples and contributed to technical aspects of the microbiome profiles within the lab of MGS. MG and SSA coordinated the collection of samples and data for START. MRS, ABB, PJS, PM, and SET coordinated the collection of samples and data for CHILD. JCS, MAZ, RDS, MS, and SSA were major contributors in writing the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
The Hamilton Integrated Research Ethics Board (HIREB) approved the research protocols for studies on human samples (CHILD, Malcolm Sears REB Project #07-2929; START, Sonia Anand HIREB Project # 10-640) and each participating parent gave signed informed consent. Our study conforms to the Declaration of Helsinki.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.