ABSTRACT
With the advent of next-generation sequencing and microbial community characterization, we are beginning to understand the key factors that shape early-life microbial colonization and associated health outcomes. Studies characterizing infant microbial colonization have focused mostly on bacteria in the microbiome and have largely neglected fungi (the mycobiome), despite their relevance to mucosal infections in healthy infants. In this pilot study, we characterized the skin, oral, and anal mycobiomes of infants over the first month of life (n = 17) and the anal and vaginal mycobiomes of mothers (n = 16) by internal transcribed spacer 2 (ITS2) amplicon sequencing. We found that infant mycobiomes differed by body site, with the infant mycobiomes at the anal sites being different from those at the skin and oral sites. The relative abundances of body site-specific taxa differed by birth mode, with significantly more Candida albicans fungi present on the skin of vaginally born infants on day 30 and significantly more Candida orthopsilosis fungi present in the oral cavity of caesarean section-born infants throughout the first month of life. We found the mycobiomes within individual infants to be variable over the first month of life, and vaginal birth did not result in infant mycobiomes that were more similar to the mother’s vaginal mycobiome. Therefore, although vertical transmission of specific fungal isolates from mother to infant has been reported, it is likely that other sources (environment, other caregivers) also contribute to early-life mycobiome establishment. Thus, future longitudinal studies of mycobiome and bacterial microbiome codevelopment, with dense sampling from birth to beyond the first month of life, are warranted.
IMPORTANCE Humans are colonized by diverse fungi (mycobiome), which have received much less study to date than colonizing bacteria. We know very little about the succession of fungal colonization in early life and whether it may relate to long-term health. To better understand fungal colonization and its sources, we studied the skin, oral, and anal mycobiomes of healthy term infants and the vaginal and anal mycobiomes of their mothers. Generally, infants were colonized by few fungal taxa, and fungal alpha diversity did not increase over the first month of life. There was no clear community maturation over the first month of life, regardless of body site. Key body-site-specific taxa, but not overall fungal community structures, were impacted by birth mode. Thus, additional studies to characterize mycobiome acquisition and succession throughout early life are needed to form a foundation for research into the relationship between mycobiome development and human disease.
INTRODUCTION
Early-life microbial colonization plays an important role in the healthy development of an infant. Extensive epidemiological surveys of the human microbiome have improved our understanding of how factors such as birth mode, diet, and exposure to antibiotics shape the bacterial microbiome in early life (1). Recent studies have also begun to elucidate key members of the bacterial community that may play specific roles in the health outcomes of infants. For example, the genera Lachnospira, Veillonella, Faecalibacterium, and Rothia have been shown to be protective against asthma development if present in the fecal bacterial microbiome early in life (3 months of age) (2). In comparison to the infant bacterial microbiome, the structure and potential connections to health outcomes of the infant fungal microbiome (mycobiome) have received much less attention.
Recent characterization of the infant mycobiome through amplicon sequencing of the internal transcribed spacer (ITS) regions of the rRNA locus (ITS region 1 [ITS1] and/or ITS2) has shown infant gastrointestinal tracts to be colonized predominantly with Saccharomyces, Candida, and Malassezia genera (3–7). Similarly to colonizing bacteria, colonizing fungi have been linked to health outcomes, including early-life (1 to 11 months of age) fecal colonization, with Candida and Rhodotorula being linked to atopy and asthma (6). Additionally, targeted approaches have demonstrated direct transmission of specific fungal species from mother to infant (8, 9). However, how the very early stages of fungal colonization occur and whether there is birth mode-specific fungal community transmission from mother to child remain unknown. To gain a better understanding of early fungal community establishment, we characterized the skin, oral, and anal mycobiomes of a cohort of infants (n = 17) over the first 30 days of life through amplicon sequencing of ITS2. For comparison purposes, maternal vaginal and anal mycobiomes were also characterized. We measured differences in mycobiomes with respect to body site and maturational dynamics and for impact of birth mode on community composition.
RESULTS
Infant mycobiomes vary by body site.To understand how early infant mycobiomes develop according to body site, we analyzed fungal DNA isolated from swabs of the skin, oral cavity, and anus (a mix of skin and feces) of infants at regular intervals over the first 30 days of life. Using principal-coordinate analysis of weighted UniFrac distances, infant anal mycobiomes clustered separately from oral and skin mycobiomes (Fig. 1a) (P = 0.01, R2 = 0.039). Clustering of anal samples was driven, in part, by the increased relative abundances of Candida albicans in anal mycobiomes versus skin mycobiomes (P = 0.003) and oral mycobiomes (P = 0.015), as well as increased abundances of Candida parapsilosis in anal samples versus skin samples (P = 0.020) (false-discovery rate [FDR] adjusted) (Fig. 1b to e). Oral mycobiome alpha diversity (observed species) was significantly lower than that of skin and anal mycobiomes using all samples across the first 30 days of life (P = 0.002 and 0.001, respectively) (Fig. 2a). However, comparing body sites at a specific time point (see Fig. S1a in the supplemental material), alpha diversity was not consistently lower in the oral mycobiome.
FIG S1
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Infant mycobiomes vary by body site. (a) Principal-coordinate analysis of weighted UniFrac distances for infant skin, oral, and anal mycobiomes over the first 30 days of life. Box plots shown along each axis represent the median and interquartile range and indicate the distribution of samples along the given axis. Each point represents a single sample and is colored by body site as follows: skin, yellow (n = 58 samples); oral, teal (n = 56 samples); anal, pink (n = 60 samples). PERMANOVA values, R2 values, and P values are shown. (b and c) The same principal-coordinate analysis colored by the relative abundances of (b) Candida albicans and (c) C. parapsilosis, with anal samples denoted with a solid border. (d and e) Relative abundances of (d) C. albicans and (e) C. parapsilosis within the skin, oral, and anal mycobiomes of infants, as assessed by Wilcoxon rank sum tests with false-discovery-rate correction.
Specific taxa within the infant skin, oral, and anal mycobiomes. (a) Numbers of observed taxa within the skin, oral, and anal mycobiomes of infants, as assessed by Student’s t tests. (b) Relative abundances of fungal taxa within each body site. Each bar represents an individual sample, and samples are ordered by infant. The union of the 10 most abundant and 10 most common taxa is shown. “Other” represents taxa whose relative abundance is <10%. A full taxon legend is located in Fig. S1c.
Overall, infant mycobiomes tended to be dominated by relatively few taxa (Fig. 2b). The taxa of highest prevalence and relative abundance across all infant body sites sampled were C. parapsilosis, C. tropicalis, Saccharomyces cerevisiae, C. albicans, C. orthopsilosis, Cryptococcus pseudolongus, Cladosporium velox, Debaryomyces renaii, D. hansenii, Hanseniaspora uvarum, and Issatchenkia orientalis (also known as C. krusei [10, 11]) (Fig. 2b; see also Tables S1 and S2 in the supplemental material). Among the body sites, the most abundant and/or prevalent taxa were C. tropicalis, C. parapsilosis, S. cerevisiae, C. albicans, and C. orthopsilosis in the infant skin mycobiome; C. parapsilosis, C. tropicalis, S. cerevisiae, C. orthopsilosis, C. albicans, and Cladosporium velox in the infant oral mycobiome; and C. parapsilosis, C. tropicalis, C. albicans, S. cerevisiae, C. orthopsilosis, and Cryptococcus pseudolongus in the infant anal mycobiome (Tables S1 and S2). Similarly, the maternal mycobiomes (vaginal and anal) were also often dominated by a single taxon (Fig. S1b), with C. albicans being the most abundant (>20% in relative abundance) in the vaginal mycobiome and C. albicans, S. cerevisiae, and C. parapsilosis being the most abundant in the maternal anal mycobiome (Table S1). Of note, all of the major fungal taxa presented here have been previously reported to be human-associated taxa, and many are opportunistic pathogens in the setting of a compromised immune system (12).
TABLE S1
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
TABLE S2
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Infant mycobiomes are individualized and variable.Skin mycobiomes showed within-infant similarity over time, demonstrated by lower within-infant weighted UniFrac distances than between-infant distances (P = 0.041), as well as significant clustering of weighted UniFrac distances according to infant (permutational multivariate analysis of variance [PERMANOVA]; P = 0.054, R2 = 0.372) (Fig. 3a and b). The oral and anal mycobiomes of infants, however, exhibited high intraindividual variability for beta diversities over time with a lack of clustering by infant over time and similar average within-infant and between-infant weighted UniFrac distances (P > 0.05) (Fig. 3c to f). No individual infant’s mycobiome showed a clear trajectory toward a mature or distinct state during the first 30 days of life, regardless of body site, as demonstrated by the lack of a decrease in UniFrac distances from the same infant's day 30 sample over time (P > 0.05) (Fig. S2a to c).
FIG S2
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Infant mycobiome dynamics during the first 30 days of life. Data represent results of principal-coordinate analysis (a, c, and e) of weighted UniFrac (W-Unifrac) distances and median within-infant and between-infant weighted UniFrac distances (b, d, and f) for (a and b) skin, (c and d) oral, and (e and f) anal mycobiomes over time. Principal-coordinate analysis data (individual infant data noted by distinct colored shape; see legend) were tested with PERMANOVA; R2 values and P values are shown. Distances were compared using a Wilcoxon rank sum test (n = 17 infants).
Infant and adult (maternal) mycobiomes were similar to each other with respect to alpha and beta diversity measures. The total numbers of observed taxa were similar for infant and adult samples for all infant sample collection time points, regardless of body site (Fig. S1a). Over the first 30 days of life, alpha diversities (Shannon index values) of infant mycobiomes did not significantly change (Fig. 4). For skin mycobiomes, there was a trend toward increasing diversity over time, although the data did not reach statistical significance (P = 0.070) (Fig. 4a). Infant anal mycobiome beta diversity was similar to that of adults (P = 0.369, R2 = 0.014) (Fig. S3a), with infant mycobiomes maintaining a similar distance from adult mycobiomes over the first 30 days of life (P = 0.121) (Fig. S3b).
FIG S3
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Infant alpha diversity over time. Shannon index values corresponding to fungal OTUs from (a) skin, (b) oral, and (c) anal mycobiomes are presented (n = 17 infants). P values represent group-wise averages of permuted Spearman correlation test statistics controlled by subject. The dotted lines represent a summary linear model that controls for subject, with R2 reported.
Birth mode affects the relative abundances of specific taxa but not overall fungal community structures.To examine whether early-life mycobiomes are impacted by birth mode, we first analyzed skin fungal communities for differences in alpha and beta diversity and for differentiated taxa in infants born vaginally compared to those born by caesarean section. Note that the effect of caesarean section represents a composite of the effects of lack of labor, perinatal antibiotics, and bypassing of the vaginal canal (Table 1). As shown by weighted UniFrac distances and principal-coordinate analysis, infant skin mycobiomes did not cluster according to birth mode (P = 0.613, R2 = 0.009) (Fig. 5a), and an infant’s skin mycobiome was not more similar to the mother’s vaginal or anal mycobiome than to a randomly chosen mother’s mycobiome (Fig. S4a). Similar findings were also obtained using Bray Curtis and unweighted UniFrac distances (data not shown).
FIG S4
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
Participant dataa
Infant skin and oral mycobiomes by birth mode. (a) Principal-coordinate analysis of weighted UniFrac distances of infant skin mycobiomes colored by birth mode (vaginal = blue, n = 30 samples; caesarean section = green, n = 28 samples). PERMANOVA values, R2 values, and P values are shown. (b and d) Weighted UniFrac distances of infant (b) skin and (d) oral mycobiomes from the mother’s vaginal mycobiome by birth mode (paired analysis), as assessed by a Wilcoxon rank sum test (n = 7 vaginal birth families and 9 caesarean section birth families across all time points). (c and e) Relative abundances of (c) Candida albicans in infant skin on day 30 according to birth mode and (e) C. orthopsilosis in the oral cavity across all time points, as assessed with a Wilcoxon rank sum test (false discovery rate corrected; n = 5 vaginally born infants, n = 6 caesarean section-born infants).
To test whether this lack of differentiation by birth mode was due to high vaginal mycobiome variability across mothers, we tested for differences in the average distance of an infant’s skin sample data from the mother’s vaginal sample data. Infants born vaginally did not have a skin mycobiome that was more similar to the mother’s vaginal mycobiome than those born by caesarean section, although a nonsignificant trend was observed (P = 0.22) (Fig. 5b). This held true regardless of the infant sampling day or the beta diversity metric used. Skin mycobiome alpha diversities were also not significantly different with respect to birth mode for any time point within the first 30 days of life but did trend toward increasing over time in vaginally born infants but not in caesarean section-born infants (P = 0.071, and P = 0.331, respectively) (Fig. S5a). Interestingly, C. albicans, a common fungus inhabiting the vagina (13), had significantly higher relative abundance in skin mycobiomes of 30-day-old infants born vaginally than in skin mycobiomes of those born by caesarean section (FDR-adjusted P = 0.004) (Fig. 5c).
FIG S5
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
For oral mycobiomes, birth mode did not significantly impact the alpha diversity trajectory over time (caesarean section P = 0.238; vaginal P = 0.873) (Fig. S5a) or beta diversity clustering (P = 0.261, R2 = 0.022) (Fig. S5b). The oral mycobiome of an infant was also not more similar to that of the mother than to that of a randomly chosen mother (Fig. S4b) or to the mother’s vaginal mycobiome among the infants born vaginally than among those born by caesarean section (P = 0.917) (Fig. 5d). Caesarean section-born infants, however, did have a significantly higher relative abundance of Candida orthopsilosis than infants born vaginally (FDR-adjusted P = 0.001) (Fig. 5e).
For anal mycobiomes, alpha diversity was significantly higher among the infants born by caesarean section than among those born vaginally, but only for day 21 of life (P < 0.001), and did not significantly increase or decrease over time regardless of birth mode (P = 0.099 for birth by caesarean section; P = 0.973 for vaginal birth) (Fig. S5a). Birth mode did not significantly affect the beta diversity of the infant anal mycobiome (P = 0.97, R2 = 0.002) (Fig. S5c), and the anal mycobiome of an infant was not more similar to the mother’s mycobiome than to that of a randomly chosen mother (Fig. S4c) or to the mother’s vaginal mycobiome among the infants born vaginally than among those born by caesarean section (P = 0.513) (Fig. S5d).
DISCUSSION
Infant mycobiomes and body site differentiation.Overall, the early infant mycobiomes consisted of a few (generally <16 per sample) unique operational taxonomic units (OTUs) and taxa, with individual mycobiomes often being dominated by one specific taxon (Fig. 2b). The finding of low fungal richness reported here is supported by a recent study of the infant fecal mycobiome which detected low fungal biomass and, on average, fewer than 12 fungal OTUs per sample over the first 2 years of life (7, 14). The most prevalent and abundant taxa across all infant body sites were primarily Candida species, with C. parapsilosis and C. tropicalis having the highest combined relative abundances (37%, 39%, and 56% of the skin, oral, and anal fungi, respectively; see Table S1 in the supplemental material). C. albicans and S. cerevisiae were also prominent (up to 15% relative abundance), depending on the body site. The detection of these species was expected, as they are well-known human-associated fungi. Previous analyses of skin-associated fungi showed adults to be predominantly colonized with the genus Malassezia, as well as being colonized with other genera, such as Candida and Saccharomyces (15, 16). Here, Malassezia was also detected in the infant skin mycobiome but accounted for only 2% of the relative abundance. This result is supported by the finding of lower relative abundance of Malassezia in the skin of children than in that of adults as shown by next-generation sequencing (17), PCR-based (18, 19), and culture-based (20, 21) approaches. The oral mycobiome of adults has previously been reported to be predominantly composed of species of the Candida (22) and Malassezia (23) genera. Previous exploration of the oral mycobiome of infants, however, has been restricted to the use of Candida-specific culturing, demonstrating the presence of numerous Candida species such as C. albicans, C. parapsilosis, C. krusei, C. guilliermondii, C. geocandidum, and C. tropicalis (24, 25). Similarly, we found the Candida species C. parapsilosis, C. tropicalis, and C. orthopsilosis, as well as S. cerevisiae, to be the most relatively abundant and prevalent fungal species in infant oral mycobiomes during the first month of life (Tables S1 and S2). Among the anal samples, C. parapsilosis, C. tropicalis, C. albicans, and S. cerevisiae were common and had the highest relative abundances (Table S1), consistent with previous studies of the infant fecal mycobiome, which reported Saccharomyces and Candida to be the predominant genera as well as the presence of other minor members (3, 5–7). Because the data presented here are based on relative abundance, which could potentially be impacted by ITS2 copy number, future studies highlighting the relative and total abundances of fungal taxa in infants are warranted.
The issue of whether or not immediate (at birth) and early-life body site differentiation of the bacterial microbiome exists is controversial. Some reports state that there is a detectable difference immediately after birth, when the infant meconium is significantly different in bacterial composition from skin and oral sites (26). Other reports, however, state that infant skin, oral, and rectal bacterial microbiomes cluster primarily according to delivery mode (27). Differentiation of the bacterial microbiome continues over time for skin, oral, and anal sites, with some studies showing detectable differences among all three sites as early as 4 to 6 weeks after birth (26, 28). Here, early-life anal mycobiomes of infants were found to be different from oral and skin mycobiomes and this was driven, in part, by the significantly higher relative abundances of C. albicans and C. parapsilosis in anal samples than in oral and skin samples (Fig. 1). Higher relative abundances of these species in the anal mycobiome may have been due, in part, to transition of the gut from an aerobic to an anaerobic state during the first weeks of life, favoring the initial expansion of opportunistic Candida species (7). For example, C. albicans is the most common fungal colonizer and pathogen of infants, likely due to its abilities to rapidly adapt to environmental changes through yeast-to-hypha morphogenesis, to adhere to human epithelial and endothelial cells, and to form biofilms (29, 30). Additionally, expansion of abundances of C. albicans and C. parapsilosis in the skin and oral mycobiomes may be limited by the presence of specific resident bacteria in these areas, as demonstrated by the ability of abundant oral taxa Streptococcus and Lactobacillus (31) to prevent the growth of (32) and to out-compete (33–35) C. albicans. Due to the small sample size and lack of bacterial sequencing in this pilot study, the mechanism determining fungal mycobiome differentiation by body site should be further explored with larger studies that include both bacterial and fungal sequencing.
Birth mode comparisons.The bacterial microbiome is highly variable with respect to birth mode; beta diversity differences are detectable immediately following birth (for skin and oral sites), and the initially higher fecal alpha diversity seen with caesarean section-born infants decreases over the first month of life to below that seen with vaginally born infants (26, 36). Many studies have also reported higher levels of specific taxa, such as Bacteroides, in the fecal bacterial microbiome of vaginally born infants than in the bacterial microbiome of those born by caesarean section (27, 36, 37). For the mycobiome, we did not detect differences in beta diversity with respect to birth mode over the first month of life regardless of sampling site (Fig. 5; see also Fig. S5 in the supplemental material). Mycobiome alpha diversities also did not significantly increase or decrease over the first month of life, regardless of birth mode and body site (Fig. S5a), which may have been due to the use of a relatively consistent food source of either breast milk or formula. Interestingly, the skin mycobiome of vaginally born infants did contain a higher relative abundance of C. albicans than that of caesarean section-born infants at 30 days of life (Fig. 5c), suggesting that either birth mode or perinatal antibiotics affect colonization by specific fungal species. A delayed increase in vagina-associated taxa in vaginally born infants was also previously seen for bacteria, albeit at much later times, where the number of maternal vagina- and infant stool-shared operational taxonomic units (OTUs) for vaginally born infants peaked between 18 and 24 months of life (36). Additionally, vertical transmission of specific C. albicans strains from mother to infant has been shown previously using culture-based and sequence-based (restriction fragment length polymorphisms) approaches (8). Due to the small sample size used here, future longitudinal studies of infants with larger cohorts that extend beyond 30 days after birth are needed to definitively determine if abundances of vagina-associated fungal taxa, including C. albicans, increase in the mycobiomes of infants born vaginally. Of note, a recent report showed low rates of C. albicans transmission from mother to infant, suggesting environmental as well as maternal sources for infant mycobiome establishment (38). Our results support the need for analyses for other mycobiome sources, as the average distance of an infant’s mycobiome from that of its mother was not significantly lower for infants born vaginally than for those born by caesarean section (Fig. 5b and d; see also Fig. S5d). Because maternal samples in this study were taken immediately following birth (after passage of amniotic fluid, the baby, and the placenta in cases of vaginal delivery), similarities between mother and infant mycobiomes could have been impacted by these factors. If mothers had been sampled just prior to delivery, perhaps the similarities between the mother and infant mycobiomes would have been higher.
Vertical transmission of skin-associated fungal species of the Malassezia genus from mothers to their caesarean section-born infants as determined using a targeted PCR-based approach was previously shown (9). We did not detect a significant increase in the relative abundances of skin-associated Malassezia species in the skin mycobiomes of infants born by caesarean section compared to vaginally born infants (data not shown). Although the maternal skin mycobiome was not analyzed here, based on the extensive evidence indicating the presence of Malassezia in the adult skin mycobiome (15, 16, 39, 40), the lack of increased skin-associated Malassezia abundance in caesarean section-born infants is not likely due to low maternal Malassezia levels. It is possible that Malassezia remains in low relative abundance within the first month of life, suggesting that, similarly to C. albicans colonization, Malassezia colonization may represent a mixture of maternal and environmental colonization sources. Additionally, because Malassezia species are basidiomycetous fungi, their low relative abundance in our results may have been due to our use of the ITS2 region, which may be biased toward the amplification of ascomycetous fungi (41, 42), for taxonomic classification.
The oral mycobiome of infants born by caesarean section contained a higher relative abundance of C. orthopsilosis than that of vaginally born infants over the first month of life (Fig. 5e). C. orthopsilosis was previously characterized as C. parapsilosis due to morphological similarities of the two species and has been reported to be a cause of invasive infections (43–45). The higher relative abundance of C. orthopsilosis in caesarean section-born infants shown in our results may have been due to their exposure to maternal antibiotics (Table 1). Because the growth of fungi is likely influenced by interkingdom interactions with bacterial community members (6), it is possible that birth mode-specific and differentially abundant bacteria may impact fungal community characteristics. For example, specific bacterial taxa, including Streptococcaceae and Lactobacillales, have been reported to have higher relative abundances in the oral cavity of infants with no maternal antibiotic exposure at the time of delivery (46), and species within these families have been shown to suppress biofilm formation by specific Candida species (47–50). Therefore, antibiotic administration may have allowed expansion of C. orthopsilosis in the oral cavity, or, because Candida species are often found in the hospital environment, antibiotic administration may have permitted higher rates of C. orthopsilosis acquisition from the environment (51–54). Due to the small number of subjects in this study, an analysis of bacterial (31) and fungal microbiome interactions was not completed and would likely have been underpowered with respect to detection of significant correlations. Longitudinal studies with large infant cohorts are needed in the future to understand the role of interkingdom relationships in microbiome and mycobiome development during infancy.
Fungal variability within young infants.In this study of early-life mycobiomes, we observed that the fungal communities of the skin, oral, and anal cavities were variable and that there was no clear progression toward a different or mature infant mycobiome composition (Fig. S2) during the first month of life. Additionally, for anal mycobiomes, early-life fungal compositions did not differ from those of adult females (Fig. S3). The anal and oral mycobiomes were diverse across infants and time, where within-infant beta diversity was not significantly lower than between-infant beta diversity (Fig. 3). In addition, mycobiome alpha diversity did not increase over time, regardless of body site, demonstrating a lack of new species acquisition and limited richness increases within the first month of life (Fig. 4). Unfortunately, there have been few longitudinal studies of very early (first month of life) bacterial microbiomes in humans that allow a comparison with our mycobiome results. One study, focused on bacterial communities of feces, observed that alpha diversity increases over the first month of life (55), in contrast to our findings determined with fungal communities. Longer-term studies of fecal microbiome development have shown that bacterial communities mature toward an adult-like state during the first 2 to 3 years of life (36, 56). Expanded longitudinal studies to similarly characterize mycobiome development are warranted. In particular, dense sampling in the earliest months of life, along with fungal and bacterial sequencing, would allow detection of key fungus-bacterium interactions that might prevent or enable colonization with specific taxa.
In summary, this report presents characteristics of mycobiome development during the first month of life. Initial body site differences, including the anal mycobiome having significantly more C. albicans and C. parapsilosis fungi than the skin and oral mycobiomes (Fig. 1a to c) and the oral mycobiome being significantly less diverse than the skin and anal mycobiomes (Fig. 2a), were detected. The skin community composition (beta diversity) was more stable than that of the oral and anal mycobiomes over time, where the skin mycobiome of each of the infants was more similar to that of the same infant over time than to those of other infants (Fig. 3a and b). Alpha diversity did not significantly increase over time (Fig. 4a), and we did not observe community composition maturation of the mycobiome over the first 30 days of life (Fig. S2a). Birth mode impacted the relative abundances of body site-specific taxa, including an increase in levels of skin-associated C. albicans in vaginally born infants and an increase in levels of orally associated C. orthopsilosis in infants born by caesarean section (Fig. 5c and e). These findings highlight the need for additional studies to better characterize the sources of and variability in fungal colonization early in life.
MATERIALS AND METHODS
Participant data.Mothers (n = 17) who were pronounced to be healthy by their physicians were recruited from the University Hospital Puerto Rico Medical Center under approval from the Institutional Review Boards of the University of Puerto Rico Medical Sciences campus (approval A9710112) and from the Rio Piedras campus (approval 1011-107). Informed consent was obtained from all participants. All participants received standard-of-care treatment, such as prophylaxis antibiotics for mothers delivering by caesarean section and/or by vaginal birth if positive for group B Streptococcus (GBS) (Table 1).
Sample collection and sequencing.Samples from infants (n = 17) were collected with sterile swabs immediately following birth and again at days 3, 7, 14, 21, 27 and 30. Infant sample sites included the oral mucosa, forehead, and anal cavity. Samples from vaginal and anal sites of mothers (n = 16) were also collected immediately following the birth of the infant (Table 1). Of note, both infant and maternal anal samples were composed of both skin and fecal matter. Samples were placed on ice and frozen at −80°C within 2 h of collection until further processing was performed. DNA was extracted from thawed samples using a Mo Bio Powersoil kit (Mo Bio Laboratories, Inc., Carlsbad, CA). The ITS2 region of the fungal genome was amplified by PCR as previously described (5) with the following modifications. ITS2-specific oligonucleotide primers were modified to include a barcode of 6 bp for both the forward and the reverse primers to improve the amount of multiplexing performed on a single sequencing run. In addition, degenerate base pairs were included on the 5′ end of each primer to improve the quality of sequencing through increased length diversity and reductions in base pair homogeneity during photo acquisition. Finally, Kapa HiFi polymerase (Kapa Biosystems, Woburn, MA) was used to improve sequencing quality and reduce rates of PCR-related errors (57). ITS2 amplicons were generated under the following conditions: 95°C for 5 min; 30 cycles of 98°C for 20 s, 65°C for 15 s, and 72°C for 15 s; and a final extension of 72°C for 5 min. Amplicons were purified using a QIAQuick PCR purification kit (Qiagen, Germantown, MD) and eluted into 30 μl of PCR-grade water. Purified amplicons were quantitated with a Qubit 2.0 fluorometer (Invitrogen, Eugene, OR) using a Qubit double-stranded DNA (dsDNA) high-sensitivity (HS) assay kit (Invitrogen) and pooled for Illumina library construction using a TruSeq Nano kit (Illumina, San Diego, CA). Sequencing was performed at the University of Minnesota Genomics Center on an Illumina MiSeq system using a 2 × 250-bp paired-end version 2 MiSeq reagent kit (Illumina) as previously described.
Mycobiome characterization.Raw sequencing reads were preprocessed with SHI7 using the TruSeq3-2 adaptor trimming option (58). Primers were removed by trimming the first 25 bases of each read and truncating the sequences to a maximum read length of 150 bases. Reads shorter than 150 bases were dropped (see Table S3 in the supplemental material). Using a validated protocol, NINJA-OPS was used to align preprocessed reads against the UNITE v7 singleton-exclusive dynamic fungal ITS database release (31 January 2016) for NINJA-OPS using default options (59, 60). The resulting operational taxonomic unit (OTU) table was filtered to keep only those samples with at least 50 aligned reads and OTUs occurring in two or more samples (Table S3). OTU counts were converted to relative abundance data for all downstream analyses. Taxa correlating in relative abundance at a level above 95% across samples were collapsed into one taxon, decreasing the number of taxa from 123 to 96 across all samples. Alpha and beta diversity were calculated using QIIME v.1.8.0 (61), with the use of the dynamic 100 ghost-tree (0116_s) for phylogenetic metrics (62). Only one sample per individual, per time point, and per body site was retained for subsequent analysis, as some participants had multiple samples collected per day. For oral samples, only those taken prior to feeding were used for analysis, except in cases where only a sampling collected following feeding was available. Samples from feces and those listed as inoculated were dropped from the analysis due to inadequate sample size. Subsequent analysis of the data was performed in R, with the code available for download at the following URL: https://github.com/TonyaWard/fungal_infant1.
TABLE S3
Copyright © 2018 Ward et al.This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.
To ensure that the levels of OTUs reported here were valid, a secondary analysis was performed using de novo OTU picking with the default setting of UCLUST (63) through QIIME v.1.8.0. The numbers of OTUs found with this method are reported in Table S2. Of note, a mean of 94% of de novo OTUs per sample failed to align to a known reference sequence.
Statistical analyses.Statistical analyses were performed in R. PERMANOVA was applied to test for differences in beta diversity centroids across sample types (Fig. 1a, 3, and 5a; see also Fig. S3a and S4b and c), and Wilcoxon rank sum tests were used to test for differences in the mean beta diversity distance from one sample type to another, including intra- and interinfant distances (Fig. 3 and 5b and d; see also Fig. S2). Testing for changes in alpha diversity (Fig. 4) and distance to self over time (Fig. S2) was performed using a permutation test of within-subject Spearman correlations. Student’s t tests were used to test for differences in alpha diversity across birth modes (within a time point and body site; Fig. S1) and across body sites (Fig. 2a). Testing for differentiated taxa across sample types was performed using a Kruskal-Wallis test for comparisons of data from more than two groups (across infant body sites, for example) or Wilcoxon rank sum tests for comparisons of data from two groups (delivery mode; Fig. 5c and e). Body site comparisons between specific sites were preformed pairwise with Wilcoxon rank sum tests with false-discovery-rate correction (Fig. 1d and e). Analyses comparing infant samples to maternal samples were not paired unless noted otherwise in the figure legend (Fig. 5; see also Fig. S4 and S5).
Accession number(s).Sequencing data are available under BioProject accession number PRJNA393442 within NCBI.
ACKNOWLEDGMENTS
M.G.D.-B. and C.A.G. developed the study, and M.G.D.-B. provided the samples. T.H. prepared fungal amplicons for sequencing. T.L.W., T.H., G.A.-G., and D.K. performed the data analysis. All authors assisted with the analysis design and data interpretation. T.L.W. and C.A.G. prepared the manuscript with edits from all authors. All authors read and approved the final manuscript.
This work was supported in part by National Center for Advancing Translational Sciences award UL1TR000114 (T.L.W.), National Institutes of Health awards R01 DK090989 (M.G.D.-B.) and R01 AI121383 (D.K.), the C&D Research Fund (M.G.D.-B.), and Minnesota Partnership for Biotechnology and Medical Genomics award 13.28 (D.K.). The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health.
FOOTNOTES
- Received October 17, 2017.
- Accepted February 13, 2018.
- Copyright © 2018 Ward et al.
This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.
