Skip to main content
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems
  • Log in
  • My alerts
  • My Cart

Main menu

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
  • ASM Journals
    • Antimicrobial Agents and Chemotherapy
    • Applied and Environmental Microbiology
    • Clinical Microbiology Reviews
    • Clinical and Vaccine Immunology
    • EcoSal Plus
    • Infection and Immunity
    • Journal of Bacteriology
    • Journal of Clinical Microbiology
    • Journal of Microbiology & Biology Education
    • Journal of Virology
    • mBio
    • Microbiology and Molecular Biology Reviews
    • Microbiology Resource Announcements
    • Microbiology Spectrum
    • Molecular and Cellular Biology
    • mSphere
    • mSystems

User menu

  • Log in
  • My alerts
  • My Cart

Search

  • Advanced search
mSystems
publisher-logosite-logo

Advanced Search

  • Home
  • Articles
    • Latest Articles
    • Special Issues
    • COVID-19 Special Collection
    • Editor's Picks
    • Special Series: Sponsored Minireviews and Video Abstracts
    • Archive
  • Topics
    • Applied and Environmental Science
    • Ecological and Evolutionary Science
    • Host-Microbe Biology
    • Molecular Biology and Physiology
    • Novel Systems Biology Techniques
    • Early-Career Systems Microbiology Perspectives
  • For Authors
    • Getting Started
    • Submit a Manuscript
    • Scope
    • Editorial Policy
    • Submission, Review, & Publication Processes
    • Organization and Format
    • Errata, Author Corrections, Retractions
    • Illustrations and Tables
    • Nomenclature
    • Abbreviations and Conventions
    • Publication Fees
    • Ethics
  • About the Journal
    • About mSystems
    • Editor in Chief
    • Board of Editors
    • For Reviewers
    • For the Media
    • For Librarians
    • For Advertisers
    • Alerts
    • RSS
    • FAQ
Research Article | Clinical Science and Epidemiology

Development of the Human Mycobiome over the First Month of Life and across Body Sites

Tonya L. Ward, Maria Gloria Dominguez-Bello, Tim Heisel, Gabriel Al-Ghalith, Dan Knights, Cheryl A. Gale
Catherine Lozupone, Editor
Tonya L. Ward
aBioTechnology Institute, University of Minnesota, St. Paul, Minnesota, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Maria Gloria Dominguez-Bello
bDepartments of Biochemistry and Microbiology and Anthropology, Rutgers University, New Brunswick, New Jersey, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Tim Heisel
cDepartment of Pediatrics, University of Minnesota, Minneapolis, Minnesota, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Gabriel Al-Ghalith
dBioinformatics and Computational Biology, University of Minnesota, Minneapolis, Minnesota, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Gabriel Al-Ghalith
Dan Knights
aBioTechnology Institute, University of Minnesota, St. Paul, Minnesota, USA
eDepartment of Computer Science and Engineering, University of Minnesota, Minneapolis, Minnesota, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Cheryl A. Gale
cDepartment of Pediatrics, University of Minnesota, Minneapolis, Minnesota, USA
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for Cheryl A. Gale
Catherine Lozupone
University of Colorado Denver
Roles: Editor
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
DOI: 10.1128/mSystems.00140-17
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

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 

Alpha diversity and taxa of maternal and infant mycobiomes over the first 30 days of life. (a) Total numbers of observed species in each body site per day (for the infants, total n = 58 skin samples, total n = 56 oral samples, and total n = 60 anal samples; for the mothers, n = 16 vaginal and 15 anal samples) were compared within each sampling day by the use of pairwise Wilcox tests and false-discovery-rate correction. No tests yielded a significant difference. (b) Relative abundances of fungal taxa within each maternal body site. Each bar represents an individual sample, and samples are ordered by mother. (c) Full taxon legend for Fig. 2b. “Other” represents taxa whose relative abundance is <10%. Download FIG S1, PDF file, 0.3 MB.
Copyright © 2018 Ward et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 1 
  • Open in new tab
  • Download powerpoint
FIG 1 

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.

FIG 2 
  • Open in new tab
  • Download powerpoint
FIG 2 

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 

Top five most abundant infant and maternal fungal taxa by body site. Download TABLE S1, DOCX file, 0.1 MB.
Copyright © 2018 Ward et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

TABLE S2 

Top five most prevalent infant and maternal fungal taxa by body site. Download TABLE S2, DOCX file, 0.1 MB.
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 

Infant mycobiome maturation over the first 30 days of life. Data represent the distance (weighted UniFrac) of each infant mycobiome from the corresponding day 30 mycobiome, according to body site. (a) Skin mycobiomes. (b) Oral mycobiomes. (c) Anal mycobiomes. n = 11 infants. P value data 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. Download FIG S2, TIF file, 1 MB.
Copyright © 2018 Ward et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 3 
  • Open in new tab
  • Download powerpoint
FIG 3 

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 

Diversity in infant and maternal anal mycobiomes. (a) Principal-coordinate analysis of weighted UniFrac distances of maternal (gray) and infant (pink) anal mycobiomes over the first 30 days of life. PERMANOVA values, R2 values, and P values are shown. (b) Distance (weighted UniFrac) of each infant anal mycobiome from the corresponding adult anal mycobiome, over the first 30 days of life (n = 17 infants and 16 mothers). The dotted lines represent a summary linear model that controls for subject, with R2 reported. The P value data represent group-wise averages of permuted Spearman correlation test statistics controlled by subject. Download FIG S3, TIF file, 0.9 MB.
Copyright © 2018 Ward et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

FIG 4 
  • Open in new tab
  • Download powerpoint
FIG 4 

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 

Infant mycobiome similarity to maternal mycobiome. Data represent the distance (weighted UniFrac) of each infant’s sample from the mother’s anal or vaginal sample (paired; orange) or from a permuted randomly chosen mother's sample (random; blue) for the infant’s skin (a), oral (b), and anal (c) mycobiomes, tested with Wilcoxon rank sum tests (n = 16). Download FIG S4, TIF file, 1.1 MB.
Copyright © 2018 Ward et al.

This content is distributed under the terms of the Creative Commons Attribution 4.0 International license.

View this table:
  • View inline
  • View popup
TABLE 1 

Participant dataa

FIG 5 
  • Open in new tab
  • Download powerpoint
FIG 5 

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 

Infant mycobiome characteristics according to birth mode. (a) Data corresponding to the Shannon index values determined for the skin, oral, and anal mycobiomes for caesarean section (green) and vaginally born (blue) infants over time are shown (n = 17 infants). The dotted lines represent a summary linear model that controls for subject, with R2 reported. P value data represent group-wise averages of permuted Spearman correlation test statistics controlled by subject. t tests were also performed on each sampling day to test for differences across birth modes. The only significant difference within day, after false-discovery-rate correction, is shown for the anal samples on day 21 (dotted box). (b) Principal-coordinate analysis of oral weighted UniFrac distances of vaginally born (blue, n = 26) and caesarean-section-born (green, n = 30) infants over the first 30 days of life. (c) Principal-coordinate analysis of anal weighted UniFrac distances of vaginally born (blue, n = 28) and caesarean-section-born (green, n = 32) infants over the first 30 days of life. (d) Weighted UniFrac distance of infant anal mycobiomes from the mother’s vaginal mycobiome by birth mode (paired analysis), as assessed by a Wilcox test (n = 7 vaginal birth families and 9 caesarean section birth families, with samples from all time points). PERMANOVA values, R2 values, and P values are shown for principal-coordinate analysis. Download FIG S5, PDF file, 0.5 MB.
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 

Sequencing processing outcomes. Download TABLE S3, DOCX file, 0.1 MB.
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.

REFERENCES

  1. 1.↵
    1. Mueller NT,
    2. Bakacs E,
    3. Combellick J,
    4. Grigoryan Z,
    5. Dominguez-Bello MG
    . 2015. The infant microbiome development: mom matters. Trends Mol Med 21:109–117. doi:10.1016/j.molmed.2014.12.002.
    OpenUrlCrossRefPubMed
  2. 2.↵
    1. Arrieta MC,
    2. Stiemsma LT,
    3. Dimitriu PA,
    4. Thorson L,
    5. Russell S,
    6. Yurist-Doutsch S,
    7. Kuzeljevic B,
    8. Gold MJ,
    9. Britton HM,
    10. Lefebvre DL,
    11. Subbarao P,
    12. Mandhane P,
    13. Becker A,
    14. McNagny KM,
    15. Sears MR,
    16. Kollmann T; CHILD Study,
    17. Mohn WW,
    18. Turvey SE,
    19. Finlay BB
    . 2015. Early infancy microbial and metabolic alterations affect risk of childhood asthma. Sci Transl Med 7:307ra152. doi:10.1126/scitranslmed.aab2271.
    OpenUrlAbstract/FREE Full Text
  3. 3.↵
    1. LaTuga MS,
    2. Ellis JC,
    3. Cotton CM,
    4. Goldberg RN,
    5. Wynn JL,
    6. Jackson RB,
    7. Seed PC
    . 2011. Beyond bacteria: a study of the enteric microbial consortium in extremely low birth weight infants. PLoS One 6:e27858. doi:10.1371/journal.pone.0027858.
    OpenUrlCrossRefPubMed
  4. 4.↵
    1. Strati F,
    2. Di Paola M,
    3. Stefanini I,
    4. Albanese D,
    5. Rizzetto L,
    6. Lionetti P,
    7. Calabrò A,
    8. Jousson O,
    9. Donati C,
    10. Cavalieri D,
    11. De Filippo C
    . 2016. Age and gender affect the composition of fungal population of the human gastrointestinal tract. Front Microbiol 7:1227. doi:10.3389/fmicb.2016.01227.
    OpenUrlCrossRef
  5. 5.↵
    1. Heisel T,
    2. Podgorski H,
    3. Staley CM,
    4. Knights D,
    5. Sadowsky MJ,
    6. Gale CA
    . 2015. Complementary amplicon-based genomic approaches for the study of fungal communities in humans. PLoS One 10:e0116705. doi:10.1371/journal.pone.0116705.
    OpenUrlCrossRefPubMed
  6. 6.↵
    1. Fujimura KE,
    2. Sitarik AR,
    3. Havstad S,
    4. Lin DL,
    5. Levan S,
    6. Fadrosh D,
    7. Panzer AR,
    8. LaMere B,
    9. Rackaityte E,
    10. Lukacs NW,
    11. Wegienka G,
    12. Boushey HA,
    13. Ownby DR,
    14. Zoratti EM,
    15. Levin AM,
    16. Johnson CC,
    17. Lynch SV
    . 2016. Neonatal gut microbiota associates with childhood multisensitized atopy and T cell differentiation. Nat Med 22:1187–1191. doi:10.1038/nm.4176.
    OpenUrlCrossRefPubMed
  7. 7.↵
    1. Wampach L,
    2. Heintz-Buschart A,
    3. Hogan A,
    4. Muller EEL,
    5. Narayanasamy S,
    6. Laczny CC,
    7. Hugerth LW,
    8. Bindl L,
    9. Bottu J,
    10. Andersson AF,
    11. de Beaufort C,
    12. Wilmes P
    . 2017. Colonization and succession within the human gut microbiome by archaea, bacteria, and microeukaryotes during the first year of life. Front Microbiol 8:738. doi:10.3389/fmicb.2017.00738.
    OpenUrlCrossRef
  8. 8.↵
    1. Bliss JM,
    2. Basavegowda KP,
    3. Watson WJ,
    4. Sheikh AU,
    5. Ryan RM
    . 2008. Vertical and horizontal transmission of Candida albicans in very low birth weight infants using DNA fingerprinting techniques. Pediatr Infect Dis J 27:231–235. doi:10.1097/INF.0b013e31815bb69d.
    OpenUrlCrossRefPubMed
  9. 9.↵
    1. Nagata R,
    2. Nagano H,
    3. Ogishima D,
    4. Nakamura Y,
    5. Hiruma M,
    6. Sugita T
    . 2012. Transmission of the major skin microbiota, Malassezia, from mother to neonate. Pediatr Int Off J Jpn Pediatr Soc 54:350–355. doi:10.1111/j.1442-200X.2012.03563.x.
    OpenUrlCrossRef
  10. 10.↵
    1. Kurtzman CP,
    2. Robnett CJ,
    3. Basehoar-Powers E
    . 2008. Phylogenetic relationships among species of Pichia, Issatchenkia and Williopsis determined from multigene sequence analysis, and the proposal of Barnettozyma gen. nov., Lindnera gen. nov. and Wickerhamomyces gen. nov. FEMS Yeast Res 8:939–954. doi:10.1111/j.1567-1364.2008.00419.x.
    OpenUrlCrossRefPubMedWeb of Science
  11. 11.↵
    1. Borman AM,
    2. Linton CJ,
    3. Oliver D,
    4. Palmer MD,
    5. Szekely A,
    6. Johnson EM
    . 2010. Rapid molecular identification of pathogenic yeasts by pyrosequencing analysis of 35 nucleotides of internal transcribed spacer 2. J Clin Microbiol 48:3648–3653. doi:10.1128/JCM.01071-10.
    OpenUrlAbstract/FREE Full Text
  12. 12.↵
    1. Wheeler ML,
    2. Limon JJ,
    3. Underhill DM
    . 2017. Immunity to commensal fungi: detente and disease. Annu Rev Pathol 12:359–385. doi:10.1146/annurev-pathol-052016-100342.
    OpenUrlCrossRef
  13. 13.↵
    1. Drell T,
    2. Lillsaar T,
    3. Tummeleht L,
    4. Simm J,
    5. Aaspõllu A,
    6. Väin E,
    7. Saarma I,
    8. Salumets A,
    9. Donders GGG,
    10. Metsis M
    . 2013. Characterization of the vaginal micro- and mycobiome in asymptomatic reproductive-age Estonian women. PLoS One 8:e54379. doi:10.1371/journal.pone.0054379.
    OpenUrlCrossRefPubMed
  14. 14.↵
    1. Schei K,
    2. Avershina E,
    3. Øien T,
    4. Rudi K,
    5. Follestad T,
    6. Salamati S,
    7. Ødegård RA
    . 2017. Early gut mycobiota and mother-offspring transfer. Microbiome 5:107. doi:10.1186/s40168-017-0319-x.
    OpenUrlCrossRef
  15. 15.↵
    1. Findley K,
    2. Oh J,
    3. Yang J,
    4. Conlan S,
    5. Deming C,
    6. Meyer JA,
    7. Schoenfeld D,
    8. Nomicos E,
    9. Park M; NIH Intramural Sequencing Center Comparative Sequencing Program,
    10. Kong HH,
    11. Segre JA
    . 2013. Topographic diversity of fungal and bacterial communities in human skin. Nature 498:367–370. doi:10.1038/nature12171.
    OpenUrlCrossRefPubMedWeb of Science
  16. 16.↵
    1. Zhang E,
    2. Tanaka T,
    3. Tajima M,
    4. Tsuboi R,
    5. Nishikawa A,
    6. Sugita T
    . 2011. Characterization of the skin fungal microbiota in patients with atopic dermatitis and in healthy subjects. Microbiol Immunol 55:625–632. doi:10.1111/j.1348-0421.2011.00364.x.
    OpenUrlCrossRefPubMed
  17. 17.↵
    1. Jo JH,
    2. Deming C,
    3. Kennedy EA,
    4. Conlan S,
    5. Polley EC,
    6. Ng WI; NISC Comparative Sequencing Program,
    7. Segre JA,
    8. Kong HH
    . 2016. Diverse human skin fungal communities in children converge in adulthood. J Invest Dermatol 136:2356–2363. doi:10.1016/j.jid.2016.05.130.
    OpenUrlCrossRef
  18. 18.↵
    1. Jang SJ,
    2. Lim SH,
    3. Ko JH,
    4. Oh BH,
    5. Kim SM,
    6. Song YC,
    7. Yim SM,
    8. Lee YW,
    9. Choe YB,
    10. Ahn KJ
    . 2009. The investigation on the distribution of Malassezia yeasts on the normal Korean skin by 26S rDNA PCR-RFLP. Ann Dermatol 21:18–26. doi:10.5021/ad.2009.21.1.18.
    OpenUrlCrossRefPubMed
  19. 19.↵
    1. Sugita T,
    2. Suzuki M,
    3. Goto S,
    4. Nishikawa A,
    5. Hiruma M,
    6. Yamazaki T,
    7. Makimura K
    . 2010. Quantitative analysis of the cutaneous Malassezia microbiota in 770 healthy Japanese by age and gender using a real-time PCR assay. Med Mycol 48:229–233.
    OpenUrlPubMed
  20. 20.↵
    1. Faergemann J,
    2. Fredriksson T
    . 1980. Age incidence of Pityrosporum orbiculare on human skin. Acta Derm Venereol 60:531–533.
    OpenUrlPubMedWeb of Science
  21. 21.↵
    1. Gupta AK,
    2. Kohli Y
    . 2004. Prevalence of Malassezia species on various body sites in clinically healthy subjects representing different age groups. Med Mycol 42:35–42. doi:10.1080/13693780310001610056.
    OpenUrlCrossRefPubMed
  22. 22.↵
    1. Ghannoum MA,
    2. Jurevic RJ,
    3. Mukherjee PK,
    4. Cui F,
    5. Sikaroodi M,
    6. Naqvi A,
    7. Gillevet PM
    . 2010. Characterization of the oral fungal microbiome (mycobiome) in healthy individuals. PLoS Pathog 6:e1000713. doi:10.1371/journal.ppat.1000713.
    OpenUrlCrossRefPubMed
  23. 23.↵
    1. Dupuy AK,
    2. David MS,
    3. Li L,
    4. Heider TN,
    5. Peterson JD,
    6. Montano EA,
    7. Dongari-Bagtzoglou A,
    8. Diaz PI,
    9. Strausbaugh LD
    . 2014. Redefining the human oral mycobiome with improved practices in amplicon-based taxonomy: discovery of Malassezia as a prominent commensal. PLoS One 9:e90899. doi:10.1371/journal.pone.0090899.
    OpenUrlCrossRefPubMed
  24. 24.↵
    1. Kleinegger CL,
    2. Lockhart SR,
    3. Vargas K,
    4. Soll DR
    . 1996. Frequency, intensity, species, and strains of oral Candida vary as a function of host age. J Clin Microbiol 34:2246–2254.
    OpenUrlAbstract/FREE Full Text
  25. 25.↵
    1. Russell C,
    2. Lay KM
    . 1973. Natural history of Candida species and yeasts in the oral cavities of infants. Arch Oral Biol 18:957–962. doi:10.1016/0003-9969(73)90176-3.
    OpenUrlCrossRefPubMedWeb of Science
  26. 26.↵
    1. Chu DM,
    2. Ma J,
    3. Prince AL,
    4. Antony KM,
    5. Seferovic MD,
    6. Aagaard KM
    . 2017. Maturation of the infant microbiome community structure and function across multiple body sites and in relation to mode of delivery. Nat Med 23:314–326. doi:10.1038/nm.4272.
    OpenUrlCrossRef
  27. 27.↵
    1. Dominguez-Bello MG,
    2. Costello EK,
    3. Contreras M,
    4. Magris M,
    5. Hidalgo G,
    6. Fierer N,
    7. Knight R
    . 2010. Delivery mode shapes the acquisition and structure of the initial microbiota across multiple body habitats in newborns. Proc Natl Acad Sci U S A 107:11971–11975. doi:10.1073/pnas.1002601107.
    OpenUrlAbstract/FREE Full Text
  28. 28.↵
    1. Costello EK,
    2. Carlisle EM,
    3. Bik EM,
    4. Morowitz MJ,
    5. Relman DA
    . 2013. Microbiome assembly across multiple body sites in low-birthweight infants. mBio 4:e00782-13. doi:10.1128/mBio.00782-13.
    OpenUrlAbstract/FREE Full Text
  29. 29.↵
    1. Chin VK,
    2. Lee TY,
    3. Rusliza B,
    4. Chong PP
    . 2016. Dissecting Candida albicans infection from the perspective of C. albicans virulence and omics approaches on host-pathogen interaction: a review. Int J Mol Sci 17:1643. doi:10.3390/ijms17101643.
    OpenUrlCrossRefPubMed
  30. 30.↵
    1. Mayer FL,
    2. Wilson D,
    3. Hube B
    . 2013. Candida albicans pathogenicity mechanisms. Virulence 4:119–128. doi:10.4161/viru.22913.
    OpenUrlCrossRefPubMedWeb of Science
  31. 31.↵
    1. Dominguez-Bello MG,
    2. De Jesus-Laboy KM,
    3. Shen N,
    4. Cox LM,
    5. Amir A,
    6. Gonzalez A,
    7. Bokulich NA,
    8. Song SJ,
    9. Hoashi M,
    10. Rivera-Vinas JI,
    11. Mendez K,
    12. Knight R,
    13. Clemente JC
    . 2016. Partial restoration of the microbiota of cesarean-born infants via vaginal microbial transfer. Nat Med 22:250–253. doi:10.1038/nm.4039.
    OpenUrlCrossRefPubMed
  32. 32.↵
    1. Barbosa JO,
    2. Rossoni RD,
    3. Vilela SFG,
    4. de Alvarenga JA,
    5. Velloso Mdos S,
    6. Prata MC,
    7. Jorge AO,
    8. Junqueira JC
    . 2016. Streptococcus mutans can modulate biofilm formation and attenuate the virulence of Candida albicans. PLoS One 11:e0150457. doi:10.1371/journal.pone.0150457.
    OpenUrlCrossRef
  33. 33.↵
    1. Parolin C,
    2. Marangoni A,
    3. Laghi L,
    4. Foschi C,
    5. Ñahui Palomino RA,
    6. Calonghi N,
    7. Cevenini R,
    8. Vitali B
    . 2015. Isolation of vaginal Lactobacilli and characterization of anti-Candida activity. PLoS One 10:e0131220. doi:10.1371/journal.pone.0131220.
    OpenUrlCrossRef
  34. 34.↵
    1. Sharma A,
    2. Srivastava S
    . 2014. Anti-Candida activity of spent culture filtrate of Lactobacillus plantarum strain LR/14. J Mycol Med 24:e25–e34. doi:10.1016/j.mycmed.2013.11.001.
    OpenUrlCrossRef
  35. 35.↵
    1. Köhler GA,
    2. Assefa S,
    3. Reid G
    . 2012. Probiotic interference of Lactobacillus rhamnosus GR-1 and Lactobacillus reuteri RC-14 with the opportunistic fungal pathogen Candida albicans. Infect Dis Obstet Gynecol 2012:e636474. doi:10.1155/2012/636474.
    OpenUrlCrossRef
  36. 36.↵
    1. Bokulich NA,
    2. Chung J,
    3. Battaglia T,
    4. Henderson N,
    5. Jay M,
    6. Li H,
    7. D Lieber AD,
    8. Wu F,
    9. Perez-Perez GI,
    10. Chen Y,
    11. Schweizer W,
    12. Zheng X,
    13. Contreras M,
    14. Dominguez-Bello MG,
    15. Blaser MJ
    . 2016. Antibiotics, birth mode, and diet shape microbiome maturation during early life. Sci Transl Med 8:343ra82. doi:10.1126/scitranslmed.aad7121.
    OpenUrlAbstract/FREE Full Text
  37. 37.↵
    1. Yassour M,
    2. Vatanen T,
    3. Siljander H,
    4. Hämäläinen AM,
    5. Härkönen T,
    6. Ryhänen SJ,
    7. Franzosa EA,
    8. Vlamakis H,
    9. Huttenhower C,
    10. Gevers D,
    11. Lander ES,
    12. Knip M; DIABIMMUNE Study Group,
    13. Xavier RJ
    . 2016. Natural history of the infant gut microbiome and impact of antibiotic treatment on bacterial strain diversity and stability. Sci Transl Med 8:343ra81. doi:10.1126/scitranslmed.aad0917.
    OpenUrlAbstract/FREE Full Text
  38. 38.↵
    1. Payne MS,
    2. Cullinane M,
    3. Garland SM,
    4. Tabrizi SN,
    5. Donath SM,
    6. Bennett CM,
    7. Amir LH
    . 2016. Detection of Candida spp. in the vagina of a cohort of nulliparous pregnant women by culture and molecular methods: is there an association between maternal vaginal and infant oral colonisation? Aust N Z J Obstet Gynaecol 56:179–184. doi:10.1111/ajo.12409.
    OpenUrlCrossRef
  39. 39.↵
    1. Leung MHY,
    2. Chan KCK,
    3. Lee PKH
    . 2016. Skin fungal community and its correlation with bacterial community of urban Chinese individuals. Microbiome 4:46. doi:10.1186/s40168-016-0192-z.
    OpenUrlCrossRef
  40. 40.↵
    1. Park HK,
    2. Ha MH,
    3. Park SG,
    4. Kim MN,
    5. Kim BJ,
    6. Kim W
    . 2012. Characterization of the fungal microbiota (mycobiome) in healthy and dandruff-afflicted human scalps. PLoS One 7:e32847. doi:10.1371/journal.pone.0032847.
    OpenUrlCrossRefPubMed
  41. 41.↵
    1. Bellemain E,
    2. Carlsen T,
    3. Brochmann C,
    4. Coissac E,
    5. Taberlet P,
    6. Kauserud H
    . 2010. ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiol 10:189. doi:10.1186/1471-2180-10-189.
    OpenUrlCrossRefPubMed
  42. 42.↵
    1. Bazzicalupo AL,
    2. Bálint M,
    3. Schmitt I
    . 2013. Comparison of ITS1 and ITS2 rDNA in 454 sequencing of hyperdiverse fungal communities. Fungal Ecol 6:102–109. doi:10.1016/j.funeco.2012.09.003.
    OpenUrlCrossRef
  43. 43.↵
    1. Tavanti A,
    2. Davidson AD,
    3. Gow NAR,
    4. Maiden MCJ,
    5. Odds FC
    . 2005. Candida orthopsilosis and Candida metapsilosis spp. nov. to replace Candida parapsilosis groups II and III. J Clin Microbiol 43:284–292. doi:10.1128/JCM.43.1.284-292.2005.
    OpenUrlAbstract/FREE Full Text
  44. 44.↵
    1. Bonfietti LX,
    2. Martins Mdos A,
    3. Szeszs MW,
    4. Pukiskas SB,
    5. Purisco SU,
    6. Pimentel FC,
    7. Pereira GH,
    8. Silva DC,
    9. Oliveira L,
    10. Carvalho Melhem MDS
    . 2012. Prevalence, distribution and antifungal susceptibility profiles of Candida parapsilosis, Candida orthopsilosis and Candida metapsilosis bloodstream isolates. J Med Microbiol 61:1003–1008. doi:10.1099/jmm.0.037812-0.
    OpenUrlCrossRefPubMedWeb of Science
  45. 45.↵
    1. Lovero G,
    2. Borghi E,
    3. Balbino S,
    4. Cirasola D,
    5. De Giglio O,
    6. Perdoni F,
    7. Caggiano G,
    8. Morace G,
    9. Montagna MT
    . 2016. Molecular identification and echinocandin susceptibility of Candida parapsilosis complex bloodstream isolates in Italy, 2007–2014. PLoS One 11:e0150218. doi:10.1371/journal.pone.0150218.
    OpenUrlCrossRef
  46. 46.↵
    1. Gomez-Arango LF,
    2. Barrett HL,
    3. McIntyre HD,
    4. Callaway LK,
    5. Morrison M,
    6. Dekker Nitert M
    . 2017. Antibiotic treatment at delivery shapes the initial oral microbiome in neonates. Sci Rep 7:43481. doi:10.1038/srep43481.
    OpenUrlCrossRef
  47. 47.↵
    1. Chandra J,
    2. Mukherjee PK
    . 2015. Candida biofilms: development, architecture, and resistance. Microbiol Spectr 3:10. doi:10.1128/microbiolspec.MB-0020-2015.
    OpenUrlCrossRef
  48. 48.↵
    1. Reid G,
    2. Tieszer C,
    3. Lam D
    . 1995. Influence of Lactobacilli on the adhesion of Staphylococcus aureus and Candida albicans to fibers and epithelial cells. J Ind Microbiol 15:248–253. doi:10.1007/BF01569832.
    OpenUrlCrossRefPubMedWeb of Science
  49. 49.↵
    1. Webb BC,
    2. Willcox MD,
    3. Thomas CJ,
    4. Harty DW,
    5. Knox KW
    . 1995. The effect of sodium hypochlorite on potential pathogenic traits of Candida albicans and other Candida species. Oral Microbiol Immunol 10:334–341. doi:10.1111/j.1399-302X.1995.tb00163.x.
    OpenUrlCrossRefPubMedWeb of Science
  50. 50.↵
    1. Holmes AR,
    2. Cannon RD,
    3. Jenkinson HF
    . 1995. Interactions of Candida albicans with bacteria and salivary molecules in oral biofilms. J Ind Microbiol 15:208–213. doi:10.1007/BF01569827.
    OpenUrlCrossRefPubMedWeb of Science
  51. 51.↵
    1. Sabino R,
    2. Sampaio P,
    3. Carneiro C,
    4. Rosado L,
    5. Pais C
    . 2011. Isolates from hospital environments are the most virulent of the Candida parapsilosis complex. BMC Microbiol 11:180. doi:10.1186/1471-2180-11-180.
    OpenUrlCrossRefPubMed
  52. 52.↵
    1. Sautour M,
    2. Dalle F,
    3. Olivieri C,
    4. L’ollivier C,
    5. Enderlin E,
    6. Salome E,
    7. Chovelon I,
    8. Vagner O,
    9. Sixt N,
    10. Fricker-Pap V,
    11. Aho S,
    12. Fontaneau O,
    13. Cachia C,
    14. Bonnin A
    . 2009. A prospective survey of air and surface fungal contamination in a medical mycology laboratory at a tertiary care university hospital. Am J Infect Control 37:189–194. doi:10.1016/j.ajic.2008.06.009.
    OpenUrlCrossRefPubMed
  53. 53.↵
    1. van Asbeck EC,
    2. Huang YC,
    3. Markham AN,
    4. Clemons KV,
    5. Stevens DA
    . 2007. Candida parapsilosis fungemia in neonates: genotyping results suggest healthcare workers hands as source, and review of published studies. Mycopathologia 164:287–293. doi:10.1007/s11046-007-9054-3.
    OpenUrlCrossRefPubMedWeb of Science
  54. 54.↵
    1. Welsh RM,
    2. Bentz ML,
    3. Shams A,
    4. Houston H,
    5. Lyons A,
    6. Rose LJ,
    7. Litvintseva AP
    . 2017. Survival, persistence, and isolation of the emerging multidrug-resistant pathogenic yeast Candida auris on a plastic health care surface. J Clin Microbiol 55:2996–3005. doi:10.1128/JCM.00921-17.
    OpenUrlAbstract/FREE Full Text
  55. 55.↵
    1. Hill CJ,
    2. Lynch DB,
    3. Murphy K,
    4. Ulaszewska M,
    5. Jeffery IB,
    6. O’Shea CA,
    7. Watkins C,
    8. Dempsey E,
    9. Mattivi F,
    10. Tuohy K,
    11. Ross RP,
    12. Ryan CA,
    13. O’Toole PW,
    14. Stanton C
    . 2017. Evolution of gut microbiota composition from birth to 24 weeks in the INFANTMET Cohort. Microbiome 5:21. doi:10.1186/s40168-017-0240-3.
    OpenUrlCrossRef
  56. 56.↵
    1. Yatsunenko T,
    2. Rey FE,
    3. Manary MJ,
    4. Trehan I,
    5. Dominguez-Bello MG,
    6. Contreras M,
    7. Magris M,
    8. Hidalgo G,
    9. Baldassano RN,
    10. Anokhin AP,
    11. Heath AC,
    12. Warner B,
    13. Reeder J,
    14. Kuczynski J,
    15. Caporaso JG,
    16. Lozupone CA,
    17. Lauber C,
    18. Clemente JC,
    19. Knights D,
    20. Knight R,
    21. Gordon JI
    . 2012. Human gut microbiome viewed across age and geography. Nature 486:222–227. doi:10.1038/nature11053.
    OpenUrlCrossRefPubMedWeb of Science
  57. 57.↵
    1. Gohl DM,
    2. Vangay P,
    3. Garbe J,
    4. MacLean A,
    5. Hauge A,
    6. Becker A,
    7. Gould TJ,
    8. Clayton JB,
    9. Johnson TJ,
    10. Hunter R,
    11. Knights D,
    12. Beckman KB
    . 2016. Systematic improvement of amplicon marker gene methods for increased accuracy in microbiome studies. Nat Biotechnol 34:942–949. doi:10.1038/nbt.3601.
    OpenUrlCrossRef
  58. 58.↵
    1. Al-Ghalith GA,
    2. Ang K,
    3. Hillmann B,
    4. Knights D
    . 2017. SHI7: a streamlined short-read iterative trimming pipeline. https://zenodo.org. doi:10.5281/zenodo.808832.
    OpenUrlCrossRef
  59. 59.↵
    1. Al-Ghalith GA,
    2. Montassier E,
    3. Ward HN,
    4. Knights D
    . 2016. NINJA-OPS: fast accurate marker gene alignment using concatenated ribosomes. PLoS Comput Biol 12:e1004658. doi:10.1371/journal.pcbi.1004658.
    OpenUrlCrossRefPubMed
  60. 60.↵
    1. Cline LC,
    2. Song Z,
    3. Al-Ghalith GA,
    4. Knights D,
    5. Kennedy PG
    . 2017. Moving beyond de novo clustering in fungal community ecology. New Phytol 216:629–634. doi:10.1111/nph.14752.
    OpenUrlCrossRef
  61. 61.↵
    1. Caporaso JG,
    2. Kuczynski J,
    3. Stombaugh J,
    4. Bittinger K,
    5. Bushman FD,
    6. Costello EK,
    7. Fierer N,
    8. Peña AG,
    9. Goodrich JK,
    10. Gordon JI,
    11. Huttley GA,
    12. Kelley ST,
    13. Knights D,
    14. Koenig JE,
    15. Ley RE,
    16. Lozupone CA,
    17. McDonald D,
    18. Muegge BD,
    19. Pirrung M,
    20. Reeder J,
    21. Sevinsky JR,
    22. Turnbaugh PJ,
    23. Walters WA,
    24. Widmann J,
    25. Yatsunenko T,
    26. Zaneveld J,
    27. Knight R
    . 2010. QIIME allows analysis of high-throughput community sequencing data. Nat Methods 7:335–336. doi:10.1038/nmeth.f.303.
    OpenUrlCrossRefPubMedWeb of Science
  62. 62.↵
    1. Fouquier J,
    2. Rideout JR,
    3. Bolyen E,
    4. Chase J,
    5. Shiffer A,
    6. McDonald D,
    7. Knight R,
    8. Caporaso JG,
    9. Kelley ST
    . 2016. Ghost-tree: creating hybrid-gene phylogenetic trees for diversity analyses. Microbiome 4:11. doi:10.1186/s40168-016-0153-6.
    OpenUrlCrossRef
  63. 63.↵
    1. Edgar RC
    . 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460–2461. doi:10.1093/bioinformatics/btq461.
    OpenUrlCrossRefPubMedWeb of Science
PreviousNext
Back to top
Download PDF
Citation Tools
Development of the Human Mycobiome over the First Month of Life and across Body Sites
Tonya L. Ward, Maria Gloria Dominguez-Bello, Tim Heisel, Gabriel Al-Ghalith, Dan Knights, Cheryl A. Gale
mSystems Mar 2018, 3 (3) e00140-17; DOI: 10.1128/mSystems.00140-17

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Print
Alerts
Sign In to Email Alerts with your Email Address
Email

Thank you for sharing this mSystems article.

NOTE: We request your email address only to inform the recipient that it was you who recommended this article, and that it is not junk mail. We do not retain these email addresses.

Enter multiple addresses on separate lines or separate them with commas.
Development of the Human Mycobiome over the First Month of Life and across Body Sites
(Your Name) has forwarded a page to you from mSystems
(Your Name) thought you would be interested in this article in mSystems.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Share
Development of the Human Mycobiome over the First Month of Life and across Body Sites
Tonya L. Ward, Maria Gloria Dominguez-Bello, Tim Heisel, Gabriel Al-Ghalith, Dan Knights, Cheryl A. Gale
mSystems Mar 2018, 3 (3) e00140-17; DOI: 10.1128/mSystems.00140-17
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Top
  • Article
    • ABSTRACT
    • INTRODUCTION
    • RESULTS
    • DISCUSSION
    • MATERIALS AND METHODS
    • ACKNOWLEDGMENTS
    • FOOTNOTES
    • REFERENCES
  • Figures & Data
  • Info & Metrics
  • PDF

KEYWORDS

ITS2
fungi
infant
microbiome
mycobiome

Related Articles

Cited By...

About

  • About mSystems
  • Author Videos
  • Board of Editors
  • Policies
  • Overleaf Pilot
  • For Reviewers
  • For the Media
  • For Librarians
  • For Advertisers
  • Alerts
  • RSS
  • FAQ
  • Permissions
  • Journal Announcements

Authors

  • ASM Author Center
  • Submit a Manuscript
  • Author Warranty
  • Types of Articles
  • Getting Started
  • Ethics
  • Contact Us

Follow #mSystemsJ

@ASMicrobiology

       

 

ASM Journals

ASM journals are the most prominent publications in the field, delivering up-to-date and authoritative coverage of both basic and clinical microbiology.

About ASM | Contact Us | Press Room

 

ASM is a member of

Scientific Society Publisher Alliance

Copyright © 2021 American Society for Microbiology | Privacy Policy | Website feedback

Online ISSN: 2379-5077