Bacterial and Fungal Gut Community Dynamics Over the First Five Years of Life in Predominantly Rural Communities in Ghana

doi:10.21203/rs.3.rs-78160/v1

Download PDF

Research article

Bacterial and Fungal Gut Community Dynamics Over the First Five Years of Life in Predominantly Rural Communities in Ghana

https://doi.org/10.21203/rs.3.rs-78160/v1

This work is licensed under a CC BY 4.0 License

Version 1

posted

You are reading this latest preprint version

Background

Bacterial and fungal microbiotas are increasingly recognized as important in health and disease starting early in life. However, microbiota composition has not yet been investigated in most rural, low-resource settings, and in such settings, bacterial and fungal microbiotas have not been compared. Thus, we applied 16S and ITS2 amplicon sequencing, respectively, to investigate bacterial and fungal fecal microbiotas in rural Ghanaian children from birth to 5 years of age. Corresponding maternal fecal and breastmilk microbiotas were additionally investigated.

Results

While bacterial communities differed systematically across the age spectrum in composition and diversity, the same was not observed for the fungal microbiota. We also identified a novel and dramatic change in the maternal postpartum microbiota. This change included much higher abundance of E. coli and much lower abundance of Prevotella in the first versus fourth week postpartum. While infants shared more bacterial taxa with their mother’s stool and breastmilk than with those of unrelated mothers, there were far fewer shared fungal taxa.

Conclusion

Given the known ability of commensal fungi to influence host health, the distinct pattern of their acquisition likely has important health consequences. Similarly, the dynamics of mothers’ bacterial microbiotas around the time of birth may have important consequences for their children’s health. Both topics require further study.

General Microbiology

Breastmilk

child

fungi

bacteria

The dynamics of the bacterial gut microbiota have been well studied, but mostly in resource-rich settings. These surveys have shown microbial diversity to increase with age and reach adult-like profiles by 2-3 years^1-6. The breastmilk bacterial microbiota has also been studied using culture-based^7-9 and culture-independent^{10, 11} approaches, which revealed its potential to seed the infant intestinal microbiota¹². These findings together with disease-specific surveys have identified the bacterial microbiota as a key factor shaping health in early life¹³.

In contrast, the fungal gut microbiota has been far less studied, despite emerging evidence that fungi are also key modulators of health and disease^{14, 15}. Examples include recent mouse studies demonstrating that fungi effectively calibrate mucosal immune responses with an emphasis on Th17 immunity and neutrophil function¹⁴ and that colonization with Candida albicans increases survival of systemic viral challenge¹⁵. Sequence-based surveys of infant fungal microbiotas have also revealed that a higher abundance of specific fungal taxa is associated with the development of allergic disease later in childhood¹⁶. The human fungal microbiota in early life has been surveyed in very few studies, and in select populations: Europe (Norway, Florence and Luxembourg) and America (Puerto Rico and Ecuador) ^16-20. Based on diverse study designs, Candida, Aspergillus, and Saccharomyces feature prominently in newborns and older infants. There is, however, no clear age-dependent trajectory of fungal microbiota taxonomic composition^18-20 or diversity^{19, 20} in the first two years of life, aside from lower diversity in 10-day-old newborns versus 3-month-olds or their mothers¹⁸. There were also no clear trajectories identified during the first month of life in a more granular assessment of the newborn fungal microbiota¹⁷.

The breastmilk fungal microbiota has only recently been profiled using a culture-independent approach²¹, revealing dominance by Malassesia spp. also commonly found on skin and compositional differences across diverse populations, which were confirmed in separate studies^{22, 23}. Beyond this, the fungal composition of human breastmilk has not been investigated.

Despite the importance of both the bacterial and fungal gut microbiota in influencing the developmental trajectory towards health or disease, to our knowledge these two have never been jointly investigated in early life, which is likely the most important period of development. Equally surprising, despite the well-established fact that the maternal stool and breastmilk microbiotas shape the newborn’s microbiota, no study so far has captured this interaction between mother and infant across the bacterial and fungal domains. We thus conducted a cross-sectional bacterial and fungal microbiota survey encompassing fecal microbiotas of children in the first 5 years of life plus fecal and breastmilk microbiotas of mothers of children 6 months and younger.

Study Participants

Study participant recruitment, Study participants were recruited in the Kintampo North Municipality, located in the former Brong Ahafo Region (now Bono-East Region) in the middle-belt of Ghana, and serving as a part of the Kintampo Health Research Centre (KHRC) study area. Prospective participants were selected from the database of the Kintampo Health and Demographic Surveillance System (KHDSS). Expectant mothers, or mothers of older children, were approached recruited into the study were approached during pregnancy. Inclusion criteria for this study were for participants to be (a) living in the study area, (b) healthy with no known congenital defects, (c) under 63 months, and (d) to have provided written informed consent. An exclusion criterion was refusal to provide informed consent. Study participant recruitment was stratified by the age groups outlined: the newborn period, (day of life; DOL 0-5; DOL 13-17; DOL 26-35); 3 months (DOL 83-115), 6 months (DOL 165-200), 1 year; month of life (MOL) 11-13; 2 years, MOL 22-26; 3 years, MOL 33-39; and 5 years, MOL 57-63.

Microbiota sample collection and analysis

Stool and Breastmilk Sample collection. Field workers supplied mothers with sterile containers and scoops to collect stool samples from soiled diapers. Stool samples were collected from the diapers of children with sterile plastic spoons transported to the clinical laboratory at the KHRC. Additionally, mothers of participants who were 0 to 1 or 4 to 5 weeks old and were exclusively breastfed were asked to provide stool and breast milk samples in supplied sterile containers. All samples collected in the study were kept in cold boxes with ice-packs and sent to the laboratory within 2 hours. At the clinical laboratory, stool and breastmilk samples were split into multiple 1.0-ml aliquots and stored at -80^oC until further analysis. Biospecimens were transported to the University of British Columbia on dry ice with temperature monitoring via World Courier Inc.

Stool DNA extraction for amplicon sequencing. DNA was extracted from swabs using the MagAttract PowerSoil DNA KF kit (Qiagen Cat. No. 27000-4-KF) using the manufacturer’s protocol for the KingFisher Flex platform with the following modification using approximately 200 mg bulk stool loaded into each 96-well plate with sterile wooden picks.

Breastmilk DNA extraction for amplicon sequencing. Total DNA was extracted from 0.7 – 1.0 µL breastmilk using the Qiagen DNease PowerSoil DNA extraction kit (Qiagen Cat. 12830-50) with the following modifications: frozen breastmilk samples were thawed on ice and transferred to 2.0 mL screw-top tubes (VWR Cat. 211-0440) and spun for 10 min at 4^oC at 20,000 g to pellet all particles. The lipid portion of the milk remained on top of the aqueous phase after centrifugation. The lipid layer was retained, and the aqueous phase between the pellet and whey layer carefully removed with a 200 µL micropipette. The contents of one Qiagen 0.1 mm glass bead tube (Qiagen 13118-50) was then added to each sample with 500 µL Bead Solution and 200 µL phenol:chloroform:IAA pH 7-8 (Ambion Cat. AM 9730, pH adjusted with included Tris buffer). 60 µL of C1 was then added, and samples were homogenized in a Fasprep bead mill for 30s at 5.5 m/s for two cycles with a 5 min wait in between. The remainder of the protocol was carried out according to manufacturer’s instructions and samples eluted in 2 x 30 µL (total of 60 µL).

16S amplicon sequencing

Samples were submitted to Microbiota Insights Inc. (Vancouver, BC, Canada) for PCR amplification and DNA sequencing. Library preparation was done using a previously published protocol²⁴ with the details and product information outlined below. 10 µl of the final product was used to normalize using the SepalPrep Normalization Prep Plate Kit (ThermoFisher Cat. A1051001) to 1-2 ng/µl and 5 µl of each normalized sample was pooled into a single library per 96-well plate. Library pools were further concentrated using the DNA Clean and Concentrator kit (Zymo Cat. D4013). A dilution series was performed for each of the pooled libraries for subsequent quality control steps. Each pool was analyzed using the Agilent Bioanalyzer using the High Sensitivity DS DNA assay (Agilent Cat. 5047-4626) to determine approximate fragment size, and to verify library integrity. Library pools with unintended amplicons were purified using the Qiagen QIAquick Gel Extraction Kit (Qiagen Cat. 28706). Pooled library concentrations were determined using the KAPA Library Quantification Kit for Illumina (KAPA Cat. KK4824). The final libraries were loaded at 8 pM, with an additional PhiX spike-in of 20%. The amplicon library was sequenced on the MiSeq using the MiSeq 500 Cycle V2 Reagent Kit (Illumina Cat. MS-102-2003)

Paired-end reads were assembled using MOTHUR version 1.37.2, following the MiSeq SOP²⁵ (https://www.mothur.org/wiki/MiSeq_SOP; accessed Aug 2018). Operational Taxonomic Units (OTUs) were clustered at 97% identity and classified using the GreenGenes 13_8_99 database²⁶.

ITS2 amplicon sequencing

Processing for ITS2 amplicon sequencing were identical to those used for 16S amplicon sequencing, except with the following differences: Primers and PCR conditions were previously described²⁷. Paired-end reads were trimmed at the non-overlapping ends, and high-quality reads classified using UNITE (v. 7.1)²⁸.

16S and ITS2 sequencing Quality control (QC).

The potential for contamination was addressed by co-sequencing DNA amplified from specimens and from four each of template-free PCR controls and extraction kit reagents processed the same way as the specimens (two for breastmilk specimens). Two positive controls, consisting of cloned SUP05 DNA, were also included (number of copies = 2*10^6). Breastmilk samples were extracted and PCR amplified separately from the stool samples but sequenced on the same lane. Given the different distribution of contaminating OTUs from breastmilk blanks compared to stool, QC was done for the breastmilk samples separately. Contaminating OTUs can come from two main sources: from PCR and extraction reagents, and carry-over from neighboring samples during the extraction and amplification process. OTUs were considered reagent contaminants and removed from the dataset if they were present in over 50% of blank samples and their count geometric mean plus one standard deviation was greater than that in the samples. Using this approach, OTUs with high counts in many samples (mainly including taxa commonly known to be highly abundant in these sample types) were not flagged for removal, despite their presence in blanks. For OTUs whose mean counts were high in the samples but also found in over 50% of blanks, the mean count plus standard deviation from the blanks was subtracted from each of the samples. Samples with total counts under 1000 were also removed. Once filtered, data from breastmilk and stool were recombined for further analysis.

Statistical Analysis

Alpha diversity. We used the phyloseq R package²⁹ to compute alpha and beta diversity. We estimated microbial diversity with abundance-dependent (Shannon) and independent (Observed richness) metrics, after subsampling the OTU table to account for unequal sampling depth. To test for relationships between age, we used the Wilcoxon test comparing each age bin to the Mothers of DOL 26-35 infants, as this was considered the adult microbiota group. P-values were adjusted using the Bonferroni correction and adjusted p-values under 0.05 considered significant.

Community composition. Beta diversity was computed using the Bray-Curtis index. Sample clustering was visualized using non-metric-multidimensional scaling (NMDS). To test whether child age (all age bins except maternal samples) explained community composition, we used the Adonis test from R package vegan³⁰ using age in days calculated between the time interval between participant’s date of birth and visit date. To test if community composition between 5-year-olds (age bin MOL 57-63) and younger children differed, we computed the Bray-Curtis distance of individuals in each bin to their distance to all individuals in the MOL 57-63 bin. The Wilcoxon rank-sum test was then used to compare distances between each bin to the median distance among the 5-year-olds to others in the same bin. P-values were adjusted using the Bonferroni correction and adjusted p-values under 0.05 considered significant.

Differential Abundance. To identify OTUs differentially abundant between mothers of 0-5 day old vs. 26-35 day-old infants, we used the wald test implemented in the R package DESeq2³¹. Differential abundance was considered significant if the Benjamini-Hochberg adjusted p-value was less than 0.05. To further visualize select bacterial genera differentially abundant between maternal groups and their abundance across the age spectrum, we aggregated OTU counts by genus-level assignment and confirmed that the genera differed between groups using the Wilcoxon rank-sum test.

Shared OTUs between mother and infant. For shared OTU analysis, OTUs were considered present if their relative abundance was over 0.01%. For stool and breast milk separately, we calculated the number of shared OTUs between newborns and their mothers and newborns and unrelated mothers. Because the number of unrelated pairs was much greater than related pairs, we randomly sub-sampled the unrelated pairs, selecting the same number as related pairs in each analysis. The median number of shared OTUs between related and unrelated pairs and 0-5 and 25-36 days post-partum was then compared using the Wilcoxon rank-sum test.

Sample summary: amplicon sequencing

After quality filtering, we retained a total of 2.9 million reads, representing 4492 bacterial OTUs. These OTUs represented 435 unique taxa classified at the Genus level or above. Five samples yielded fewer than 1000 reads and were excluded from downstream analysis. The remainder yielded between 2.6 and 31.2 thousand reads per sample, with most fecal samples yielding greater than 10 thousand (Supl. Figure 1A).

After quality filtering, we retained a total of 2.5 million reads representing 1772 fungal OTUs, representing 375 unique taxa classified at the Genus level or above. Fifty-nine samples contained fewer than 1000 reads and were excluded from downstream analysis. The remainder contained between 1.2 and 96.5 thousand reads per sample (Supl. Figure 1B). Samples with no detectable fungi over our count threshold were found in both infant and adult stool samples alike (Supl. Figure 1C). The proportion of samples in each age stratification with no detectable fungi was highly variable, but usually below 50%.

Bacterial and fungal alpha diversity across the age spectrum. To measure alpha diversity, we used both observed richness, the number of distinct OTUs present, and the Shannon Diversity index, which reflects both richness and evenness. We compared diversity of fecal communities in each age bin to that of the one-month postpartum mothers, representing an adult microbiota. Alpha diversity of fecal bacterial communities increased with age. Observed bacterial richness was lower than that of adults until 3 years, and Shannon diversity, until 2 years (Figure 1A). Fungal fecal microbiota diversity of only a few age bins differed from adult diversity, otherwise there was no clear trend for fungal diversity versus age, except that it was lower during the first week of life (Figure 1B). Interestingly, the mothers of those newborns, also had lower fungal richness than one-month postpartum mothers. Most samples, especially from after 2 years of age contained fewer fungal OTUs (less than 60 OTUs in most samples), while between 100 and 300 bacterial OTUs were detected during this time.

Bacterial and fungal communities across the age spectrum. Beta diversity (community composition) of fecal bacterial communities differed across the age spectrum and, expectedly, was distinct from breastmilk bacterial communities (Figure 2A-B). There was no such trend for fungal fecal communities. While child age contributed to 18% of the variance of fecal bacterial communities (PERMANOVA p = 0.001), it accounted for only 2% of variance in fecal fungal communities (PERMANOVA p = 0.015). We did not identify one NMDS axis that drove an age gradient for either community. To further assess community composition across the age spectrum, we determined the similarity of communities in each age bin to those at 5 years. Fecal bacterial communities from all age bins, except 3 years, differed from those at 5 years (Figure 2C). Expectedly, fecal fungal communities did not follow this trend.

Relative abundances of the 25 most abundant Genera clearly revealed bacterial taxonomic shifts over the age spectrum (Figure 3A). In keeping with the literature on bacteria^{4, 32}, newborns were dominated by Escherichia coli in the first week of life, and it remained abundant up to 6 months. Then, Bifidobacteria were dominant from 1 month to 1 year. Subsequently, the diversity of genera increased and Prevotella became dominant, similar to the adult microbiota. Bacteroides was abundant during the first year, until the increase in Prevotella. Children were either dominated by Bacteroides or Prevotella at any given time (Supl. Figure 2). Faecalibacterium was abundant after the first year. The fecal fungal microbiota was dominated by Candida (primarily Candida albicans and Candida tropicalis) in children of all ages (Figure 3B). Malassezia featured prominently in the first 3 months of life, while Aspergillus was prominent after the first year.

Maternal bacterial microbiota. Our study revealed dramatic differences in bacterial stool microbiotas of mothers in the first postpartum week versus one month postpartum (Figure 3A). We further explored this difference by measuring the Bray-Curtis distance between mothers’ microbiotas at both time points to their respective children at 1 week, finding that mothers in the first versus fourth post-partum week harbored a bacterial microbiota significantly more similar to that of their child (Figure 4A). We identified 68 OTUs differentially abundant between the two mother groups, representing 34 taxa at the Genus level or above (Figure 4B). Abundances of some taxa obviously contributed to the similarity of microbiotas of mothers 1 week postpartum and their children at 1 week. Thus, in both mothers and infants from one to four weeks after birth, Prevotella had increasing abundance, while Escherichia had decreasing abundance (Figure 4C,D). Other taxa, such as Faecalibacterium and Blautia, were more abundant in mothers 1 week versus 4 weeks postpartum, yet were virtually absent in their infants at 1 week. We also found that, at 1 week, newborns share more bacterial but not fungal OTUs with their mother versus with unrelated mothers, but this was no longer the case at 1 month (Supl. Figure 3A,B). OTUs most commonly shared by mother-infant pairs include members of Escherichia, Enterococcus, and Streptococcus, all of which are of higher relative abundance in newborns than in their mothers (Supl. Figure 3C). Only two fungal OTUs were commonly shared between mother-infant pairs (Supl. Figure 3D).

Bacterial and fungal communities in mother’s breastmilk. Breastmilk is a well-recognized modulator of the infant gut microbiota, serving both as a source of colonizing bacteria and bacterial growth substrates. The relationship of breastmilk and the fungal microbiota has not been well described. Thus, we compared the bacterial and fungal community composition of both breastmilk and fecal microbiotas. Breastmilk microbiotas were typically dominated by skin taxa of both bacteria (Streptococcus,Staphylococcus, or Corynebacterium; Figure 5A) and fungi (Malassezia, including M. restricta, M. globosa, and M. furfur; Figure 5B). Alpha diversity did not differ between breastmilk samples collected in the first or fourth post-partum week for either bacteria or fungi. Only minor differences in community composition were found between the first and fourth post-partum week, including only eight differentially abundant bacterial OTUs (Figure 5C) and no differentially abundant fungal OTUs. For bacterial communities, mother-infant pairs shared significantly more OTUs than unrelated pairs at both postpartum time points (Figure 5D). While fungal data was available for fewer pairs, we found a significant difference between related and unrelated pair sharing 4 weeks postpartum (Figure 5E). The most commonly shared bacterial OTUs belonged to typical skin commensal genera, Streptococcus and Staphylococcus, which had higher relative abundance in the breastmilk than in the newborn stool. Shared OTUs between mother and infant pairs also included gut commensals Bifidobacterium and Bacteroides, both of which were of higher relative abundance in the newborn’s stool than in their mother’s milk (Supl. Figure 4A). Shared fungal OTUs were much more limited, and did not have a clear pattern of differing abundance between mother and child (Supl. Figure 4B).

Here, we report for the first time the bacterial and fungal microbiota dynamics during early childhood. In conducting a direct comparison between bacterial and fungal gut communities, we demonstrated how the two follow very distinct developmental trajectories, having contrasting associations with the mother’s microbiota. We also found a previously unreported and surprising compositional change in maternal stool microbiotas during the postpartum period.

We found undetectable levels of fungi in a large subset of stool samples from children and mothers in a rural African setting. This finding is consistent with a previous study in Norway, an urban high-resource setting¹⁸. We additionally confirmed that the samples with undetectable fungi, including those from the youngest age groups, contained readily amplifiable bacterial genes providing further evidence that the lack of detectable fungi in these samples reflected their absence or extremely low abundance, rather than sample processing artefacts. However, lysis methods that successfully lysed bacteria may have failed to lyse fungal cells. Interestingly, most breastmilk samples were positive for fungi, especially those from 4 weeks postpartum, while we did not detect fungi in half of the corresponding breastfeeding infants.

As in our study, increasing bacterial^{4, 6}, and unchanging fungal^{19, 20} richness and diversity during infancy were previously reported, but never together in the same study. However, the timing of increases in bacterial richness and diversity in our study differed from those in other studies in both high- and low-resource regions. In our Ghanaian cohort, richness and diversity remained stable over the first year of life and dramatically increased thereafter. In studies of American and European cohorts, alpha diversity steadily increased from birth to two years⁶ or 3 months to 5 years⁴, respectively. Numbers of observed OTUs increased steadily from birth to 3 years among Amerindian and Malawian infants¹, and from birth to two years in an Indian setting³³. The distinct timing in our study may reflect lifestyles, diet or an environment unique to this study cohort in Ghana.

The bacterial succession we observed, with sequential dominance by E. coli, Bifidobacterium, and Prevotella, as well as sequential high abundance of Staphylococcus,Streptococcus,Bacteroides, and Faecalibacterium, was consistent with previous studies of newborn and adult microbiotas^{4, 32}. The Bacteroidetes phylum in older children and adults is typically dominated by either Bacteroides or Prevotella^{1, 34, 35}. Prevotella was dominant in adults in rural Ghana, as demonstrated among the one month post-partum women in our study. We found that Bacteroides was abundant in the Ghanaian infants prior to 1 year, after which, Prevotella colonized the children and became dominant. The ability of Bacteroides species to utilize human milk oligosaccharides (HMOs)³⁶as well as diet-derived polysaccharides³⁷ as an energy source, is thought to contribute to their unique ability to persist in both infant and adult guts. However, Prevotella appears to outcompete Bacteroides in adults with a fiber-rich diet ³⁸. This transition occurs early in life in this Ghanaian population, probably associated changes to the child’s diet at this time.

We found dominance of the gut fungal microbiota of some children by Candida sp., as has been found in some^{17, 20} but not other^{16, 19} studies, where additional fungal taxa, such as Aspergillus and Penicillium, featured prominently. We also detected an early increase in the relative abundance of Malassezia that was not seen in other studies of newborns^{17, 18}. Given the high relative abundance of this taxon in breastmilk in our study and another²¹, and evidence of its’ transmission from mother to newborn³⁹, absence of this taxon in other studies is puzzling but may reflect regional differences in maternal breastmilk fungal microbiotas, which have not been well-studied.

We were surprised to find substantial differences between the gut bacterial microbiotas of mothers 1 week versus 4 weeks postpartum. These dramatic differences were not due to cross-contamination between maternal and infant samples, as mothers in both groups were sampled using the same protocol and time points, and samples from both groups were processed simultaneously. Given the cross-sectional study design, we cannot definitively conclude that the microbiota changed in these women over this time. However, both the 3-month recruitment window and the greater microbiota similarity of mothers and their infants 1 week postpartum imply that such a change is the most likely interpretation of the data. Studies of the stool microbiota during pregnancy have yielded mixed results, with some studies showing no change in community composition during pregnancy^{40, 41} and others showing substantial shifts^{42, 43}. The latter two studies showed an increase in Proteobacteria during pregnancy, which is consistent dominance of Proteobacteria that we observed 1 week postpartum. However, our finding of a low relative abundance of Prevotella 1 week postpartum is novel, and was one of the greatest differences we observed in the postpartum period. One week postpartum, we also observed a high relative abundance of Faecalibacterium, associated with decreased inflammation and protection from inflammatory bowel diseases⁴⁴, and Blautia, negatively associated with type 1 diabetes⁴⁵ and gestational diabetes melittus (GDM) without prescribed dietary intervention ⁴³. Together, our findings suggest that differences in maternal microbiota around the time of birth reflect both increases in taxa common to newborns (Escherichia) and decreases in taxa absent in newborns (Prevotella), as well as increases in abundance of taxa not abundant in infants but with potential health benefits to the mother (Faecaibacterim and Blautia).

We did not identify fungal taxa that were differentially abundant in the gut microbiota 1 week versus 4 weeks postpartum. We also did not identify many commonly shared fungal OTUs between pairs. In contrast to the bacterial gut microbiota, the fungal gut microbiota of mothers does not appear to be substantially affected by postpartum time and does not appear to strongly influence colonization of the infant’s gut.

We found that common skin taxa dominated breastmilk bacterial (Streptococcus, Staphylococcus) and fungal (Malasezzia) communities, as previously shown^{11, 21, 46, 47}. When dominant in a mother’s breastmilk, OTUs of these taxa tended to be present at lower relative abundance in her infant’s gut microbiota, consistent with the breastmilk being a source for these organisms colonizing the infant gut. However, it is important to note that detection of their DNA in breastmilk and feces does not demonstrate that these OTUs colonize, or even survive these two environments. Further investigation is required to determine if these skin taxa have any function in breastmilk or the gut. OTUs of common gut taxa (Bifidobacterium, Bacteroides, Lactobacillus and Rothia) were also frequently shared between a mother’s breastmilk and her infant’s gut microbiota, as previously reported⁴⁶. However, the gut commensals tended to have higher relative abundance in the gut than the breastmilk, suggesting that the infant’s (or mother’s) gut was a source for these organisms colonizing the breastmilk. Overall, our study identifies strong and novel associations between the microbiotas of breastmilk and the infant gut, particularly in the case of bacteria versus fungi.

With its cross-sectional design, this study cannot definitively prove that differences among age groups represent microbiota changes over time; although, such changes are the most compelling interpretation of the results. Sample size, a lack of data on infant and mother health, and on mode of delivery, limited the analysis to effects of age and not other host factors that also may have influenced the microbiotas.

The goal of this study was to assess, for the first time, the bacterial and fungal gut microbiota composition of children in Ghana to help identify microbial signatures relevant to the health of women and their children. In doing so, we not only identified bacterial and fungal community dynamics shared with other study cohorts, but also identified characteristics not previously reported. Unique to the Ghanaian cohort, alpha diversity was low and stable during the first year of life, and Prevotella became the dominant bacterial genus between the first and second year. Also, the mother’s bacterial microbiota changed drastically between 1 week and 4 weeks postpartum. We also contrasted fungal and bacterial community composition across early life in a consolidated study. While mothers are a major source of bacteria colonizing newborns, the sources of colonizing fungi are less apparent.

Future Directions

This small cross-sectional study broadly assessed ecology mother-newborn bacterial and fungal microbiotas, together with the microbiota in children over the first five years in rural Ghana. Crucial insights from this study inform future approaches to identify microbial signatures relevant to the health of women and their children. Based on our findings, longitudinal studies will address the following questions:

What demographic or environmental factors explain differences in maternal microbiotas after childbirth? And are these changes associated with maternal health or disease?
Are there host microbiota (bacterial and fungal) signatures associated with increased risk for, or protection from severe disease, especially in the newborn period?
Are commensal fungi identified in mother-infant pairs transient, or do they represent a stable community of commensal microbes co-inhabiting the gut with bacteria?

These studies will inform novel interventions targeting the microbiota, to improve the health of mothers and their children.

Ethics approval and consent to participate

All work in this study was in accordance to ethical guidelines at the University of British Columbia and done under the approved ethics protocol number H11-01423. Separate ethical approval was obtained from the Kintampo Health Research Centre Institutional Ethics Committee (KHRCIEC_2016-14).

Consent for Publication

Not Applicable

Availability of data and material

The datasets supporting the conclusions of this article are available in the National Centre for Biotechnology Information Sequence Read Archive (NCBI SRA) accessible via the following URL: https://www.ncbi.nlm.nih.gov/bioproject/PRJNA658595.

During the peer-review process, codes produced to analyze these data are available through Dropbox, https://www.dropbox.com/sh/iedyqnghu4ugdku/AAB2z69zvRxjM6pinyG7nju0a?dl=0.

Competing interests

The authors declare that they have no competing interests.

Funding

This work was supported the Canadian Institutes of Health Research (PJT-148781) and the Ghana Health Service infrastructure at Kintampo Health Research Centre.

Authors’ contributions

SOA, UG, NA, TRK, KPA, YE and PP contributed to the study design. UG, YE, SK, and DD facilitated participant recruitment and sample collection. NA and WWM generated the amplicon sequencing data. NA performed the statistical analyses and prepared the manuscript. WWM and SOA substantially edited the manuscript. All authors contributed to editing the manuscript.

Acknowledgements

We are grateful to the Management and Staff of Kintampo Health Research Centre, who supported the conduct of this study without with very little funding. Our sincere thanks to the Kintampo community, especially the families that allowed and consented for mothers and their children to participate in this study.

Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiota viewed across age and geography. Nature 2012; 486:222-7.
Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO. Development of the human infant intestinal microbiota. PLoS Biol 2007; 5:e177.
Subramanian S, Huq S, Yatsunenko T, Haque R, Mahfuz M, Alam MA, et al. Persistent gut microbiota immaturity in malnourished Bangladeshi children. Nature 2014; 510:417-21.
Stewart CJ, Ajami NJ, O'Brien JL, Hutchinson DS, Smith DP, Wong MC, et al. Temporal development of the gut microbiota in early childhood from the TEDDY study. Nature 2018; 562:583-8.
Backhed F, Roswall J, Peng Y, Feng Q, Jia H, Kovatcheva-Datchary P, et al. Dynamics and Stabilization of the Human Gut Microbiota during the First Year of Life. Cell Host Microbe 2015; 17:690-703.
Bokulich NA, Chung J, Battaglia T, Henderson N, Jay M, Li H, et al. Antibiotics, birth mode, and diet shape microbiota maturation during early life. Sci Transl Med 2016; 8:343ra82.
Jost T, Lacroix C, Braegger C, Chassard C. Assessment of bacterial diversity in breast milk using culture-dependent and culture-independent approaches. Br J Nutr 2013; 110:1253-62.
Jimenez E, Delgado S, Maldonado A, Arroyo R, Albujar M, Garcia N, et al. Staphylococcus epidermidis: a differential trait of the fecal microbiota of breast-fed infants. BMC Microbiol 2008; 8:143.
Solis G, de Los Reyes-Gavilan CG, Fernandez 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.
Cabrera-Rubio R, Collado MC, Laitinen K, Salminen S, Isolauri E, Mira A. The human milk microbiota changes over lactation and is shaped by maternal weight and mode of delivery. Am J Clin Nutr 2012; 96:544-51.
Moossavi S, Sepehri S, Robertson B, Bode L, Goruk S, Field CJ, et al. Composition and Variation of the Human Milk Microbiota Are Influenced by Maternal and Early-Life Factors. Cell Host Microbe 2019; 25:324-35.e4.
Martin V, Maldonado-Barragan A, Moles L, Rodriguez-Banos M, Campo RD, Fernandez L, et al. Sharing of bacterial strains between breast milk and infant feces. J Hum Lact 2012; 28:36-44.
Arrieta MC, Stiemsma LT, Amenyogbe N, Brown EM, Finlay B. The intestinal microbiota in early life: health and disease. Frontiers in immunology 2014; 5:427.
Shao TY, Ang WXG, Jiang TT, Huang FS, Andersen H, Kinder JM, et al. Commensal Candida albicans Positively Calibrates Systemic Th17 Immunological Responses. Cell Host Microbe 2019; 25:404-17.e6.
Jiang TT, Shao TY, Ang WXG, Kinder JM, Turner LH, Pham G, et al. Commensal Fungi Recapitulate the Protective Benefits of Intestinal Bacteria. Cell Host Microbe 2017; 22:809-16.e4.
Arrieta MC, Arevalo A, Stiemsma L, Dimitriu P, Chico ME, Loor S, et al. Associations between infant fungal and bacterial dysbiosis and childhood atopic wheeze in a nonindustrialized setting. J Allergy Clin Immunol 2018; 142:424-34.e10.
Ward TL, Dominguez-Bello MG, Heisel T, Al-Ghalith G, Knights D, Gale CA. Development of the Human Mycobiome over the First Month of Life and across Body Sites. mSystems 2018; 3.
Schei K, Avershina E, Oien T, Rudi K, Follestad T, Salamati S, et al. Early gut mycobiota and mother-offspring transfer. Microbiota 2017; 5:107.
Strati F, Di Paola M, Stefanini I, Albanese D, Rizzetto L, Lionetti P, et al. Age and Gender Affect the Composition of Fungal Population of the Human Gastrointestinal Tract. Front Microbiol 2016; 7:1227.
Wampach L, Heintz-Buschart A, Hogan A, Muller EEL, Narayanasamy S, Laczny CC, et al. Colonization and Succession within the Human Gut Microbiota by Archaea, Bacteria, and Microeukaryotes during the First Year of Life. Front Microbiol 2017; 8:738.
Boix-Amoros A, Puente-Sanchez F, du Toit E, Linderborg KM, Zhang Y, Yang B, et al. Mycobiome profiles in breast milk from healthy women depend on mode of delivery, geographic location and interaction with bacteria. Appl Environ Microbiol 2019.
Heisel T, Nyaribo L, Sadowsky MJ, Gale CA. Breastmilk and NICU surfaces are potential sources of fungi for infant mycobiomes. Fungal Genet Biol 2019; 128:29-35.
Moossavi S, Fehr K, Derakhshani H, Sbihi H, Robertson B, Bode L, et al. Human milk fungi: environmental determinants and inter-kingdom associations with milk bacteria in the CHILD Cohort Study. BMC Microbiol 2020; 20:146.
Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol 2013; 79:5112-20.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol 2009; 75:7537-41.
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. Applied and Environmental Microbiology 2006; 72:5069-72.
Gweon HS, Oliver A, Taylor J, Booth T, Gibbs M, Read DS, et al. PIPITS: an automated pipeline for analyses of fungal internal transcribed spacer sequences from the Illumina sequencing platform. Methods Ecol Evol 2015; 6:973-80.
Koljalg U, Larsson KH, Abarenkov K, Nilsson RH, Alexander IJ, Eberhardt U, et al. UNITE: a database providing web-based methods for the molecular identification of ectomycorrhizal fungi. New Phytol 2005; 166:1063-8.
Charlop-Powers Z, Brady SF. phylogeo: an R package for geographic analysis and visualization of microbiota data. Bioinformatics 2015; 31:2909-11.
Oksanen J, Blanchet G, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. 2017. https://CRAN.R-project.org/package=vegan,
Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol 2014; 15:550.
Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, et al. Succession of microbial consortia in the developing infant gut microbiota. Proc Natl Acad Sci U S A 2011; 108 Suppl 1:4578-85.
Chandel DS, Perez-Munoz ME, Yu F, Boissy R, Satpathy R, Misra PR, et al. Changes in the Gut Microbiota After Early Administration of Oral Synbiotics to Young Infants in India. J Pediatr Gastroenterol Nutr 2017; 65:218-24.
De Filippo C, Cavalieri D, Di Paola M, Ramazzotti M, Poullet JB, Massart S, et al. Impact of diet in shaping gut microbiota revealed by a comparative study in children from Europe and rural Africa. Proc Natl Acad Sci U S A 2010; 107:14691-6.
Lin A, Bik EM, Costello EK, Dethlefsen L, Haque R, Relman DA, et al. Distinct distal gut microbiota diversity and composition in healthy children from Bangladesh and the United States. PLoS One 2013; 8:e53838.
Marcobal A, Barboza M, Sonnenburg ED, Pudlo N, Martens EC, Desai P, et al. Bacteroides in the infant gut consume milk oligosaccharides via mucus-utilization pathways. Cell Host Microbe 2011; 10:507-14.
Flint HJ, Duncan SH, Scott KP, Louis P. Links between diet, gut microbiota composition and gut metabolism. Proc Nutr Soc 2015; 74:13-22.
Wu GD, Chen J, Hoffmann C, Bittinger K, Chen YY, Keilbaugh SA, et al. Linking long-term dietary patterns with gut microbial enterotypes. Science 2011; 334:105-8.
Nagata R, Nagano H, Ogishima D, Nakamura Y, Hiruma M, Sugita T. Transmission of the major skin microbiota, Malassezia, from mother to neonate. Pediatr Int 2012; 54:350-5.
Bisanz JE, Enos MK, PrayGod G, Seney S, Macklaim JM, Chilton S, et al. Microbiota at Multiple Body Sites during Pregnancy in a Rural Tanzanian Population and Effects of Moringa-Supplemented Probiotic Yogurt. Appl Environ Microbiol 2015; 81:4965-75.
DiGiulio DB, Callahan BJ, McMurdie PJ, Costello EK, Lyell DJ, Robaczewska A, et al. Temporal and spatial variation of the human microbiota during pregnancy. Proc Natl Acad Sci U S A 2015; 112:11060-5.
Koren O, Goodrich JK, Cullender TC, Spor A, Laitinen K, Backhed HK, et al. Host remodeling of the gut microbiota and metabolic changes during pregnancy. Cell 2012; 150:470-80.
Ferrocino I, Ponzo V, Gambino R, Zarovska A, Leone F, Monzeglio C, et al. Changes in the gut microbiota composition during pregnancy in patients with gestational diabetes mellitus (GDM). Sci Rep 2018; 8:12216.
Sokol H, Pigneur B, Watterlot L, Lakhdari O, Bermudez-Humaran LG, Gratadoux JJ, et al. Faecalibacterium prausnitzii is an anti-inflammatory commensal bacterium identified by gut microbiota analysis of Crohn disease patients. Proc Natl Acad Sci U S A 2008; 105:16731-6.
Murri M, Leiva I, Gomez-Zumaquero JM, Tinahones FJ, Cardona F, Soriguer F, et al. Gut microbiota in children with type 1 diabetes differs from that in healthy children: a case-control study. BMC Med 2013; 11:46.
Murphy K, Curley D, O'Callaghan TF, O'Shea CA, Dempsey EM, O'Toole PW, et al. The Composition of Human Milk and Infant Faecal Microbiota Over the First Three Months of Life: A Pilot Study. Sci Rep 2017; 7:40597.
Williams JE, Carrothers JM, Lackey KA, Beatty NF, York MA, Brooker SL, et al. Human Milk Microbial Community Structure Is Relatively Stable and Related to Variations in Macronutrient and Micronutrient Intakes in Healthy Lactating Women. J Nutr 2017; 147:1739-48.

msbacteriafungiBMCmicrobiolsupplement.pdf

Download PDF

Version 1

posted

You are reading this latest preprint version

Bacterial and Fungal Gut Community Dynamics Over the First Five Years of Life in Predominantly Rural Communities in Ghana

Status:

Version 1

Abstract

Figures

Introduction

Methods

Results

Discussion

Conclusion

Declarations

References

Supplementary Files

Status:

Version 1