Skip to main content

Fungal diversity in the gut microbiome of young South African children



The fungal microbiome, or mycobiome, is a poorly described component of the gut ecosystem and little is known about its structure and development in children. In South Africa, there have been no culture-independent evaluations of the child gut mycobiota. This study aimed to characterise the gut mycobiota and explore the relationships between fungi and bacteria in the gut microbiome of children from Cape Town communities.


Stool samples were collected from children enrolled in the TB-CHAMP clinical trial. Internal transcribed spacer 1 (ITS1) gene sequencing was performed on a total of 115 stool samples using the Illumina MiSeq platform. Differences in fungal diversity and composition in relation to demographic, clinical, and environmental factors were investigated, and correlations between fungi and previously described bacterial populations in the same samples were described.


Taxa from the genera Candida and Saccharomyces were detected in all participants. Differential abundance analysis showed that Candida spp. were significantly more abundant in children younger than 2 years compared to older children. The gut mycobiota was less diverse than the bacterial microbiota of the same participants, consistent with the findings of other human microbiome studies. The variation in richness and evenness of fungi was substantial, even between individuals of the same age. There was significant association between vitamin A supplementation and higher fungal alpha diversity (p = 0.047), and girls were shown to have lower fungal alpha diversity (p = 0.003). Co-occurrence between several bacterial taxa and Candida albicans was observed.


The dominant fungal taxa in our study population were similar to those reported in other paediatric studies; however, it remains difficult to identify the true core gut mycobiota due to the challenges set by the low abundance of gut fungi and the lack of true gut colonising species. The connection between the microbiota, vitamin A supplementation, and growth and immunity warrants exploration, especially in populations at risk for micronutrient deficiencies. While we were able to provide insight into the gut mycobiota of young South African children, further functional studies are necessary to explain the role of the mycobiota and the correlations between bacteria and fungi in human health.

Peer Review reports


Fungi are poorly studied organisms of the gut ecosystem, and it is estimated that they only represent 0.001 – 0.1% of the organisms in the gastrointestinal tract. Nonetheless, this still equates to about a billion organisms that have substantially larger genomes than that of prokaryotes in the gut [1]. Fungi therefore have the potential to contribute a substantial amount of genetic material and metabolic products to the microbiome. These fungi and their genetic products were first referred to as the mycobiome by Ghannoum and colleagues in 2010 [2].

The mycobiota is a known reservoir for potential pathogens, also known as pathobionts, and there is increasing evidence linking the gut mycobiota to intestinal inflammatory diseases [3,4,5,6]. The disruption of fungal communities by the overgrowth of certain pathobionts has also been linked to disease in early life [7], including type 1 diabetes [6], diarrhoea-associated intestinal candidiasis [8] and focal intestinal perforation in neonates [9]. However, despite the recent rise in interest, the gut mycobiota has not been as widely explored as the bacterial microbiota. The composition and development of healthy gut mycobiota during early life has not yet been fully established and the transient nature of the mycobiota makes this difficult to do.

The Human Microbiome Project found that the gut mycobiota of healthy adult volunteers was dominated by the phylum Ascomycota and, to a lesser extent, Basidiomycota [10]. Candida spp. (Ascomycota), Saccharomyces spp. (Ascomycota), and Malassezia spp. (Basidiomycota) were the most commonly detected fungal genera. When compared with the bacterial microbiota in the same individuals, the mycobiota were found to be less diverse, and more variable over time [10]. Only a limited number of studies have investigated the infant mycobiota, and thus far only in developed nations. Studies conducted in Europe and the United States of America have identified Candida and Debaryomyces (Ascomycota) as the most abundant genera [11, 12], with Candida being especially prevalent in preterm infants [13]. Saccharomyces cerevisiae was found to be more abundant after the first year of life [12, 14].

The mycobiota is affected by factors such as physiology and lifestyle, the immune system and to a lesser extent, host genotype. Interactions with the bacterial microbiome may also influence the mycobiota [15]. Some of the more notable examples of fungal-bacterial correlation, summarized by Cui and colleagues [15], include the occurrence of Mycobacterium superinfection and aspergillosis, and the suppression of fungal growth by Pseudomonas aeruginosa in cystic fibrosis patients. The Human Microbiome Project described significant strong negative correlations between Penicillium and Faecalibacterium, and Saccharomyces and Lachnospiraceae in the healthy human gut [10]. The interactions between bacteria and fungi in the gut could influence human health and it is therefore important to describe the mycobiota and bacteriota together, as well as any significant correlations found between bacteria and fungi in the population.

In South Africa, the mycobiota has mainly been studied in wine-making; however, a sequencing-based study was recently performed in healthy adult volunteers [16]. This study reported that genera belonging to the Ascomycota, Pichia, Candida, and Cladosporium, dominated the mycobiota and that diet and geographical location played the most important roles in differentiating the gut fungi between those from urban and rural communities. There is a distinct lack of mycobiota data from low resource settings, and especially in paediatric populations. We set out to characterise the diversity and community structure of the gut mycobiota of young South African children, investigate the relationships between clinical, environmental, and socio-economic factors and the mycobiota, and identify potential relationships between bacteria and fungi in this population.


Participant demographics

From November 2017 to July 2019, baseline stool samples were collected from 116 of 218 children enrolled in the Cape Town arm of the TB-CHAMP trial. DNA was extracted successfully from 115/116 samples; one sample and matched participant data were therefore excluded from further analysis. The children included in the study had a median age of 32 months (interquartile range [IQR]:15—43), 52.2% were boys and only one child was HIV-positive. The full set of participant demographics and exposures to different clinical and environmental factors that were investigated, as outlined in the methods section, have been published as part of a study describing the bacterial microbiota in the same samples, based on 16S rRNA gene sequencing [17].

Fungal microbiota composition

Fungi were detected in all participant stool samples, but during paired-end analysis it was noted that Saccharomyces was not detected in the mock community control as expected, nor in any of the samples. This was attributed to the unusually long (480 bp) internal transcribed spacer 1 (ITS1) region in Saccharomyces spp., which did not allow overlap and merge of the paired-end reads (2 × 200 bp). Forward read only analysis successfully identified Saccharomyces in the mock community and was able to detect more taxa than paired-end analysis (Supplementary data and Table S1); therefore, further analysis was performed using the forward reads only. For the forward reads, 11,351,965 sequencing reads were produced; after clean-up and chimaera removal 9,180,236 remained and were used to produce feature data.

The phylum Ascomycota was present in the highest abundance (94.8%), followed by Basidiomycota (4.8%) (Fig. 1). The remainder of the detected phyla were present in low abundance (< 1%) and/or detected in less than 5 samples.

Fig. 1
figure 1

Taxonomic profiles of participant samples at phylum level. Each vertical bar represents a single participant sample. Features unassigned at phylum level were filtered prior to visualization. The brackets on the x-axis show the different age groups (A to E) and samples within each group are ordered by increasing age

After filtering rare features that appeared in less than 5 samples, 67 populations could be identified, 53 to genus level, 6 to family level, 3 order level, one to class level and 4 to only phylum level. Despite the variety of fungi detected, the gut mycobiota was shown to be dominated by two genera, Candida and Saccharomyces, regardless of age group (Fig. 2).

Fig. 2
figure 2

The relative abundance of the top five taxa in different age groups. Taxa that could not be classified to phylum level and features appearing in less than 5 samples were filtered prior to analysis

At genus level, differential abundance analysis showed that Candida was significantly more abundant in children under 2 years, whereas Cyberlindnera was significantly more abundant in children over 2 years. While Saccharomyces was also shown to be more abundant in children older than 3 years, and Debaryomyces in children older than 4 years, the ANCOM results were not significant between any of the age groups. No other significant changes in abundance were noted for the tested factors.

Alpha and beta diversity

Nineteen samples were excluded from diversity analyses after rarefying to a sequence count depth of 5262. At this level of rarefaction, the Good’s coverage for all samples were > 0.97. For analysis of demographic, clinical and environmental factors, samples with no data were automatically excluded. As shown by Kruskal–Wallis testing, only four of the clinical, environmental, and socio-economic factors had significant associations with fungal alpha diversity based on Shannon’s H or observed features (OF) (Fig. 3, Supplementary Tables S2 and S3). Girls were found to have significantly lower alpha diversity (p = 0.003) and those older than 6 months who had received vitamin A supplementation in the preceding six months, had higher alpha diversity (p = 0.047). However, when considering only OF, these factors were not significant (girls vs boys, p = 0.06; vitamin A supplementation, p = 0.61). Children who had received antibiotic treatment in the two weeks before sample collection had lower richness as measured by OF (p = 0.03), but this was not significant with Shannon’s H (p = 0.069). Children who had been born with normal vaginal delivery had more observed features than those born by Caesarean section (p = 0.02).

Fig. 3
figure 3

Significant differences in fungal alpha diversity as measured by A Shannon’s H alpha diversity and B observed features, based on Kruskal–Wallis testing (p < 0.05). A: Sex (Shannon). B Vitamin A supplementation in children > 6 months (Shannon). C Antibiotic exposure < 2 weeks before sample collection (OF). D Mode of birth (OF)

The alpha diversity of fungi in the gut was consistently lower and had a larger range than that of their bacterial counterparts in all the age groups, based on Shannon’s H metric (Fig. 4). Bacterial diversity was significantly lower in the younger age groups, and stabilised in children from the age of 3 years, whereas the differences in fungal alpha diversity between age groups were not significant when considering either Shannon’s H or OF (Supplementary Table S4).

Fig. 4
figure 4

A comparison of bacterial and fungal diversity in children as measured by Shannon’s H alpha diversity. After rarefaction, 5 samples were excluded from bacterial diversity analysis (n = 110), and 19 samples from the fungal analysis (n = 96)

Of the tested factors, only age contributed to fungal community dissimilarity, using the Bray–Curtis metric, Jaccard distance and PERMANOVA for significance evaluation. For both metrics, the two youngest age groups (0 to < 1 and 1 to < 2) were significantly dissimilar (adjusted p value < 0.05) to the two oldest age groups (3 to < 4 and 4 to < 5). R2 values as reported by the adonis plug-in are shown in Supplementary Table S5.

Correlations between bacteria and fungi in the gut

Both positive and negative bacterial-fungal correlations were observed in this population, although none were strong (r >  ± 0.7) (Fig. 5). The Ruminococcus gnavus group and members of the Enterobacteriaceae family exhibited positive correlations with almost all the represented fungal taxa, the strongest of which was between the R. gnavus group and Candida albicans (r = 0.451). Features belonging to C. albicans exhibited the most correlations with bacterial taxa, which is not surprising as it was the most abundant fungal species in this population. The strongest negative bacterial-fungal correlation was between Agathobacter and C. albicans (r = -0.388). Fungal-fungal correlations were mainly represented by co-occurrence, of which the association between Candida parapsilosis and the unidentified Clavispora spp. was the strongest (r = 0.518). The strongest positive correlations observed in this population were between the bacterial genera Agathobacter and Roseburia (r = 0.645), and Streptococcus and Veillonella (r = 0.613).

Fig. 5
figure 5

Correlations between bacterial and fungal taxa in the gut. Red represents negative correlations and blue represents positive correlations. Only significant correlations (p < 0.05) have been coloured


This study investigated the mycobiota composition in young South African children enrolled in the TB-CHAMP clinical trial. The richness and evenness of mycobiota in the gut varied substantially, and age was the main factor contributing to differences in fungal diversity.

Even though fungi are present in very low abundance in the gut, fungal taxa were detected in all the stool samples tested in the present study, and Candida and Saccharomyces dominated the gut mycobiota. Culture-independent studies previously performed in children identified Candida and Debaryomyces as the most abundant genera [12, 14], and similarly to the present study, Saccharomyces has been found at a higher relative abundance after infancy [12, 14, 18]. Many other taxa were detected in the present study, but they were sparsely distributed between samples, and > 80% of the mycobiota abundance in each age group was comprised of only 5 taxa. This sparseness was also seen in a cohort of healthy adults, where individual fungal operational taxonomic units (OTUs) were only detected in 3% of participants on average [19]. The majority of the most commonly detected genera in our study were similar to those found in higher abundance in adults, including Candida, Saccharomyces, Aspergillus, Malassezia, Cyberlindnera, Penicillium, and Cladosporium [10, 19]. A recent review summarised evidence that suggests that mode of birth and feeding practises both influence mycobiota composition in early life [20]. In this study, it was noted that children born by Caesarean section had less observed features than those born by normal vaginal delivery; however, no differentially abundant features were identified, and the fungal gut communities were not significantly dissimilar between these groups. Another study has suggested that fungi may be acquired from various sources such as the mother’s breastmilk, parental skin or even the hospital or home environments, due to their ubiquitous nature [14].

Fungal alpha diversity was significantly lower in girls, in contrast to other reported findings. A study in neonates found that the diversity and structure of fungi in the gut were unaffected by sex, whereas another study that included both adults and children found that the gut mycobiota in females had increased richness and abundance [11, 13]. This disparity in diversity between girls and boys was not detected in the bacterial microbiota of our study population [17]. The lower degree of diversity in girls may be a consequence of different diets, play activities and socialisation, medication use or immunological responses. We were not able to evaluate this due to limitations in our clinical data, but this important and novel finding requires further evaluation. To our knowledge, the association of vitamin A supplementation in children older than six months with increased fungal diversity has not been reported elsewhere and it is unclear what the underlying cause may be. Evidence suggests that deficiencies in micronutrients such as vitamin A may disrupt microbiota development, but published works focus on the bacterial microbiota [21,22,23]. Unlike fungal diversity, bacterial richness in this study population was not associated with vitamin A supplementation, but the bacterial community compositions trended toward significant dissimilarity [17]. There were no substantial differences in fungal alpha diversity between age groups; however, the fungal community compositions were shown to differ significantly between the younger (< 2 years) and older (4 to < 5 years) age groups, mostly attributed to the shift in the dominant genera, Candida and Saccharomyces. This shift may be ascribed to dietary changes, as older children are exposed to a larger variety of food and more yeast-containing food, like bread.

Comparable to the observations made in adults by the Human Microbiome Project, fungal gut diversity was much lower than bacterial diversity in the same participants [10]. A common theory for the relatively lower fungal gut diversity is that the gut environment is not conducive to the growth of most fungal taxa. Resident fungi must be able to grow at 37 °C to colonize the gut; therefore, many of the detected fungi are considered transient and are introduced through diet and inhalation [1, 24]. In fact, it has been suggested that the predominant “true” gut fungi are composed mainly of species from the genus Candida, and that other commonly identified fungi such as Saccharomyces cerevisiae and Cladosporium may only be present due to their abundance in food and drink, and the environment, respectively [25]. David and colleagues found that Candida, Debaryomyces, Penicillium and Scopulariopsis in the human gut were foodborne fungi associated with an animal-based diet [26]. In early life, some transient species such as Malassezia may be introduced through breastfeeding, but are later overtaken by members of the order Saccharomycetales [3].

In the present study, significant bacterial-fungal correlations were not very strong, and both positive and negative correlations were observed, which is different to what has been reported for healthy adults where positive correlations were more common [19, 27]. The strongest of the observed correlations was the co-occurrence of Candida albicans and Ruminococcus gnavus group. Similar correlations between Candida and Ruminococcus have been described in the gut [28]. The negative correlations previously found between Penicillium and Faecalibacterium, and Saccharomyces and the family Lachnospiraceae in the healthy adult gut [8] were not replicated in this study. Several fungi also displayed positive correlations with the family Enterobacteriaceae in this study; this correlation has recently been attributed to cooperation between these bacteria and fungi to facilitate colonisation [29]. Overall, fungal-fungal and bacterial-fungal correlations were weaker and less common than bacterial-bacterial correlations in the present study, but there is not enough evidence to determine whether this is a common trend in the human gut.

Despite extensive data collection efforts, the present study was limited by incomplete data for several clinical and epidemiological factors, and the lack of comprehensive dietary history for the participants. Another limitation of our study was the use of only the ITS1 region, which has a large variation in read length. While most fungi have short ITS1 regions (< 260 bp), some of the genera detected commonly in humans, such as Saccharomyces, have long regions that require longer paired-end sequencing reads (2 × 300 bp or more). Shorter regions could also potentially be over-amplified during the sequencing library preparation. An advantage of using ITS1 over ITS2 is the fact that ITS2 may be biased toward amplifying the Ascomycota, which are already extremely dominant in the gut [30]. ITS-based studies are also hampered by copy-number variation, which limits the ability to accurately quantify the abundance of species. We were not able to identify the core gut mycobiota due to the uncertainty of the colonisation potential of the detected fungi, and the low abundance of fungi in the gut. Despite these shortcomings, it was possible to assign taxonomy to at least the genus level in most cases, which could be attributed to the highly variable nature and subsequent discriminatory power of the ITS1 region.


This study provides insight into the gut mycobiota of young African children and shows that the dominant taxa in this population are similar to those reported in other settings. Only age, sex and vitamin A supplementation were shown to be associated with the fungal diversity and composition. The connection between the microbiota, vitamin A supplementation, and growth and immunity warrants exploration, especially in populations at risk for micronutrient deficiencies. The low abundance of most of the detected fungal taxa and the lack of true gut colonising species among fungi remains an issue for mycobiota studies, due to the resultant difficulties in establishing the core mycobiota and determining which fungal taxa are truly affected by clinical and environmental factors. Functional studies are also necessary to explain the role of the mycobiota and the correlations between bacteria and fungi in human health.


Study design

This was a sub-study of the ongoing phase III, cluster randomised, double-blinded, placebo-controlled tuberculosis child multidrug-resistant preventive therapy (TB-CHAMP) trial ( TB-CHAMP aims to assess the efficacy and safety of levofloxacin as preventive therapy in children exposed to multidrug-resistant (MDR) tuberculosis (TB) within the household. This sub-study was a continuation of previous work that explored the bacterial gut microbiota of the children at their baseline visit, before randomisation into the treatment or placebo groups [17].

Population and participant eligibility

The study population has been described previously [17]. Briefly, children from the Cape Town metropolitan area under the age of 5 years who were residing in the same household as an enrolled adult index case (≥ 18 years of age) with bacteriologically confirmed pulmonary MDR-TB were eligible. The following were considered exclusion criteria: TB treatment in the previous 12 months; isoniazid or any fluoroquinolone treatment for ≥ 14 days at screening and known exposure to an isoniazid-susceptible TB index case.

Sample and participant data collection

One stool sample per child was collected without preservative in 25 mL faecal containers with spoons (Lasec, South Africa) at the baseline visit. Samples were transported to the laboratory, where they were homogenised and stored at -80ºC. A pilot study previously determined that the storage and transport conditions in this setting did not significantly affect microbiota composition, including fungi [31]. Demographic, clinical, and environmental data were collected during this visit using physical examinations, case report forms and lifestyle questionnaires. This included, but was not limited to, age, sex, HIV status, medical history (including mode of birth, breastfeeding habits, antibiotic use, deworming, and hospitalisation), anthropometric measurements, and exposure to environmental factors such as cigarette smoke, indoor cooking fires, pets, and day-care. Drinking water supply, ablution type and housing structure type were also recorded as a surrogate for socio-economic status. For age-related analyses, participants were grouped by year of life: A (0 to < 1 years) (n = 24), B (1 to < 2 years) (n = 25), C (2 to < 3 years) (n = 19), D (3 to < 4 years) (n = 24), and E (4 to < 5 years) (n = 23).

DNA extraction and sequencing

DNA was extracted as part of a previous study investigating the bacterial gut microbiota of the study participants [17]. Briefly, the QIAamp PowerFecal DNA Isolation Kit (Qiagen, Germany), which includes a bead-beating step, was used according to the manufacturer’s instructions and stored at -20 °C. Sequencing was performed at the Centre for Proteomic and Genomic Research (CPGR) in Cape Town, on the Illumina MiSeq platform. Targeted sequencing was performed using the Illumina Fungal Metagenomic Demonstrated Sequencing Protocol and internal transcribed spacer 1 (ITS1) primer pool [32, 33] (Table 1). The amplicon PCR was adapted by increasing the number of cycles to 35 and the forward and reverse primer pools were added at concentrations of 2 μM and 1 μM, respectively. The MiSeq Reagent v3 Kit (600 cycles) was used to generate sequencing libraries for 2 × 200 bp paired-end read sequencing. Sequencing libraries were spiked with 10% of a 5 pM PhiX sequencing control, and the ZymoBIOMICS Microbial Community DNA standard (Zymo Research, USA) and batched negative extraction controls were included in the run.

Table 1 ITS1 primer pool used for Illumina sequencing

Microbiome analysis and statistical testing

Sequence analysis and taxonomic classification

Sequence analysis was performed using the Quantitative Insights Into Microbial Ecology (QIIME2 2020.8) bioinformatics platform [34] in conjunction with other bioinformatics packages described below, on the Stellenbosch University high performance computing cluster (HPC) 2 ( ITS1 regions were extracted using ITSxpress [35] in order to remove read-through sequences from shorter ITS1 amplicons and improve read merging in paired-end sequences. The dada2 plug-in [36] was used to perform read error correction, quality filtering and chimaera removal, and feature tables were generated using amplicon sequence variants (ASVs). Taxonomy was assigned to the sequence features outputted from dada2, using a classifier trained with the feature-classifier plug-in [37, 38] on the UNITE fungal database version 8.2 ( Forward read analysis was performed as described for paired-end reads, except for the ITSxpress step, which was replaced with adapter- and read-through trimming using the cutadapt plug-in [39].

Microbial community diversity

Diversity analyses were performed at the feature level, and rarefaction, supported by Good’s coverage, was performed prior to analysis. The evenness and richness in participant samples (alpha diversity) was demonstrated using the Shannon’s H metric [40] and observed features (OF), and the significance between groups, based on clinical, demographic, and socio-economic factors, was determined using Kruskal–Wallis pairwise tests [41]. Spearman correlation analysis was performed for numerical data [42]. The Shannon alpha diversity of fungal communities was compared to that of bacterial communities in the same participant samples, investigated/described previously [17]. Microbial community dissimilarity (beta diversity) was determined with principal coordinate analysis (PCoA) based on the Bray–Curtis dissimilarity metric [43] and dissimilarity based on presence/absence was evaluated with Jaccard distance [44]. Statistical differences were determined by PERMANOVA [45]. For Bray–Curtis similarity, the adonis plug-in [46] was also implemented to report R2 values to show the strength of associations between tested factors. Benjamini–Hochberg False Discovery Rate (BH-FDR) multiple test correction was performed where appropriate and the threshold for significance was 0.05 for all p values and adjusted p values (q values) [47].

Differential abundance

Differentially abundant features for factors such as age, medication use, method of birth, breastfeeding, and exposure to pets and smoke were investigated at genus level using Analysis of composition of microbiomes (ANCOM) [48]. ANCOM performs pairwise testing and generates a test statistic (W), with a significance threshold, which shows how many times a hypothesis could be rejected for a specific genus. As recommended, only features present ≥ 20 times overall and in at least 25% of the samples were used for ANCOM analyses, to avoid noise from less abundant features.

Correlations between bacterial and fungal microbiota

FastSpar, an implementation of the Sparse Correlations for Compositional data (SparCC) method, was used to investigate the correlations between bacterial and fungal species in the gut environment and to test the significance of the correlations [49, 50]. Samples with less than 500 feature counts were excluded, and only features that were present in more than one sample and made up at least 0.5% of the overall abundance were included in correlation analysis. The correlations were plotted in R using corrplot version 0.89 [51]. FastSpar produces p values following a permutation-based approach and bootstrapping; these p values and raw correlation values can be found in Supplementary Tables S6 and S7.

Availability of data and materials

The sequencing data analysed during the current study are available through the NCBI Sequence Read Archive under BioProject PRJNA692358.


  1. Auchtung TA, Fofanova TY, Stewart CJ, Nash AK, Wong MC, Gesell JR, et al. Investigating colonization of the healthy adult gastrointestinal tract by fungi. MSphere. 2018;3(2):1–16.

    Article  Google Scholar 

  2. Ghannoum MA, Jurevic RJ, Mukherjee PK, Cui F, Sikaroodi M, Naqvi A, et al. Characterization of the oral fungal microbiome (Mycobiome) in healthy individuals. PLoS Pathog. 2010;6(1): e1000713.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Fujimura KE, Sitarik AR, Havstad S, Lin DL, Levan S, Fadrosh D, et al. Neonatal gut microbiota associates with childhood multisensitized atopy and T cell differentiation. Nat Med. 2016;22(10):1187–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Hoarau G, Mukherjee PK, Gower-Rousseau C, Hager C, Chandra J, Retuerto MA, et al. Bacteriome and mycobiome interactions underscore microbial dysbiosis in familial Crohn’s disease. MBio. 2016;7(5):1–11.

    Article  Google Scholar 

  5. Sokol H, Leducq V, Aschard H, Pham H-P, Jegou S, Landman C, et al. Fungal microbiota dysbiosis in IBD. Gut. 2017;66(6):1039–48.

    Article  CAS  PubMed  Google Scholar 

  6. Honkanen J, Vuorela A, Muthas D, Orivuori L, Luopajärvi K, Tejesvi MVG, et al. Fungal dysbiosis and intestinal inflammation in children with beta-cell autoimmunity. Front Immunol. 2020;11:1–14.

    Article  CAS  Google Scholar 

  7. Ward TL, Knights D, Gale CA. Infant fungal communities: current knowledge and research opportunities. BMC Med. 2017;15(1):30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Ezeonu IM, Ntun NW, Ugwu KO. Intestinal candidiasis and antibiotic usage in children: case study of Nsukka, South Eastern Nigeria. Afr Health Sci. 2018;17(4):1178.

    Article  Google Scholar 

  9. Coates EW. distinctive distribution of pathogens associated with peritonitis in neonates with focal intestinal perforation compared with necrotizing enterocolitis. Pediatrics. 2005;116(2):e241–6.

    Article  PubMed  Google Scholar 

  10. Nash AK, Auchtung TA, Wong MC, Smith DP, Gesell JR, Ross MC, et al. The gut mycobiome of the human microbiome project healthy cohort. Microbiome. 2017;5(1):153.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Strati F, Di Paola M, Stefanini I, Albanese D, Rizzetto L, Lionetti P, et al. Age and gender affect the composition of fungal population of the human gastrointestinal tract. Front Microbiol. 2016;7:1–16.

    Article  Google Scholar 

  12. James SA, Phillips S, Telatin A, Baker D, Ansorge R, Clarke P, et al. Preterm infants harbour a rapidly changing mycobiota that includes candida pathobionts. J Fungi. 2020;6(4):273.

    Article  Google Scholar 

  13. Willis KA, Purvis JH, Myers ED, Aziz MM, Karabayir I, Gomes CK, et al. Fungi form interkingdom microbial communities in the primordial human gut that develop with gestational age. FASEB J. 2019;33(11):12825–37.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Schei K, Avershina E, Øien T, Rudi K, Follestad T, Salamati S, et al. Early gut mycobiota and mother-offspring transfer. Microbiome. 2017;5(1):107.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Cui L, Morris A, Ghedin E. The human mycobiome in health and disease. Genome Med. 2013;5(7):63.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Kabwe MH, Vikram S, Mulaudzi K, Jansson JK, Makhalanyane TP. The gut mycobiota of rural and urban individuals is shaped by geography. BMC Microbiol. 2020;20(1):257.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. NelVanZyl K, Whitelaw AC, Hesseling AC, Seddon JA, Demers AM, Newton-Foot M. Association between clinical and environmental factors and the gut microbiota profiles in young South African children. Sci Rep. 2021;11(1):15895.

    Article  CAS  Google Scholar 

  18. Zou R, Wang Y, Duan M, Guo M, Zhang Q, Zheng H. Dysbiosis of gut fungal microbiota in children with autism spectrum disorders. J Autism Dev Disord. 2021;51(1):267–75.

    Article  PubMed  Google Scholar 

  19. Mahnic A, Rupnik M. Different host factors are associated with patterns in bacterial and fungal gut microbiota in Slovenian healthy cohort. PLoS ONE. 2018;13(12): e0209209.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Zhang L, Zhan H, Xu W, Yan S, Ng SC. The role of gut mycobiome in health and diseases. Therap Adv Gastroenterol. 2021;14:1–18.

    Article  CAS  Google Scholar 

  21. Hibberd MC, Wu M, Rodionov DA, Li X, Cheng J, Griffin NW, et al. The effects of micronutrient deficiencies on bacterial species from the human gut microbiota. Sci Transl Med. 2017;9(390):aal4069.

    Article  Google Scholar 

  22. Iyer N, Vaishnava S. Vitamin A at the interface of host–commensal–pathogen interactions. PLOS Pathog. 2019;15(6): e1007750.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Huda MN, Ahmad SM, Kalanetra KM, Taft DH, Alam MJ, Khanam A, et al. Neonatal vitamin A supplementation and vitamin A status are associated with gut microbiome composition in Bangladeshi infants in early infancy and at 2 years of age. J Nutr. 2019;149(6):1075–88.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Hallen-Adams HE, Suhr MJ. Fungi in the healthy human gastrointestinal tract. Virulence. 2017;8(3):352–8.

    Article  CAS  PubMed  Google Scholar 

  25. Suhr MJ, Hallen-Adams HE. The human gut mycobiome: pitfalls and potentials–a mycologists perspective. Mycologia. 2015;107(6):1057–73.

    Article  CAS  PubMed  Google Scholar 

  26. David LA, Maurice CF, Carmody RN, Gootenberg DB, Button JE, Wolfe BE, et al. Diet rapidly and reproducibly alters the human gut microbiome. Nature. 2014;505(7484):559–63.

    Article  CAS  PubMed  Google Scholar 

  27. Seelbinder B, Chen J, Brunke S, Vazquez-Uribe R, Santhaman R, Meyer A-C, et al. Antibiotics create a shift from mutualism to competition in human gut communities with a longer-lasting impact on fungi than bacteria. Microbiome. 2020;8(1):133.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Mirhakkak MH, Schäuble S, Klassert TE, Brunke S, Brandt P, Loos D, et al. Metabolic modeling predicts specific gut bacteria as key determinants for Candida albicans colonization levels. ISME J. 2021;15(5):1257–70.

    Article  CAS  PubMed  Google Scholar 

  29. Rao C, Coyte KZ, Bainter W, Geha RS, Martin CR, Rakoff-Nahoum S. Multi-kingdom ecological drivers of microbiota assembly in preterm infants. Nature. 2021;591(7851):633–8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Bellemain E, Carlsen T, Brochmann C, Coissac E, Taberlet P, Kauserud H. ITS as an environmental DNA barcode for fungi: an in silico approach reveals potential PCR biases. BMC Microbiol. 2010;10(1):189.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Nel Van Zyl K, Whitelaw AC, Newton-Foot M. The effect of storage conditions on microbial communities in stool. PLoS One. 2020;15(1):e0227486.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Illumina, Inc. Fungal sequencing and classification with the ITS Metagenomics Protocol [Internet]. San Diego: Illumina, Inc.; 2019 [Cited 2022 Aug 16]. Available from:

  33. Illumina, Inc. Fungal Metagenomic Sequencing Demonstrated Protocol [Internet]. Illumina, Inc.; 2019 [Cited 2022 Aug 16]. Available from:

  34. Bolyen E, Rideout JR, Dillon MR, Bokulich NA, Abnet CC, Al-Ghalith GA, et al. Reproducible, interactive, scalable and extensible microbiome data science using QIIME 2. Nat Biotechnol. 2019;37(8):852–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. Rivers AR, Weber KC, Gardner TG, Liu S, Armstrong SD. ITSxpress: Software to rapidly trim internally transcribed spacer sequences with quality scores for marker gene analysis. F1000Res. 2018;7:1418.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Callahan BJ, McMurdie PJ, Rosen MJ, Han AW, Johnson AJA, Holmes SP. DADA2: High-resolution sample inference from Illumina amplicon data. Nat Methods. 2016;13(7):581–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: machine learning in python. J Mach Learn Res. 2011;12(85):2825–30.

    Google Scholar 

  38. Bokulich NA, Kaehler BD, Rideout JR, Dillon M, Bolyen E, Knight R, et al. Optimizing taxonomic classification of marker-gene amplicon sequences with QIIME 2’s q2-feature-classifier plugin. Microbiome. 2018;6(1):90.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 2011;17(1):10.

    Article  Google Scholar 

  40. Shannon CE, Weaver W W. The mathematical theory of communication. Champaign, Illinois: The University of Illinois Press; 1949.

    Google Scholar 

  41. Kruskal WH, Wallis WA. Use of ranks in one-criterion variance analysis. J Am Stat Assoc. 1952;47(260):583–621.

    Article  Google Scholar 

  42. Spearman C. The proof and measurement of association between two things. Am J Psychol. 1904;15(1):72–101.

    Article  Google Scholar 

  43. Bray JR, Curtis JT. an ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27(4):325–49.

    Article  Google Scholar 

  44. Jaccard P. Nouvelles recherches sur la distribution florale. Bull la Société vaudoise des Sci Nat. 1908;44:223–70.

    Google Scholar 

  45. Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26(1):32–46.

    Article  Google Scholar 

  46. Jari Oksanen, F. Guillaume Blanchet, Michael Friendly, Roeland Kindt, Pierre Legendre, Dan McGlinn, Peter R. Minchin, R. B. O’Hara, Gavin L. Simpson, Peter Solymos, M. Henry H. Stevens, Eduard Szoecs and HW. Vegan: Community Ecology Package. 2018.

  47. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B. 1995;57(1):289–300.

    Article  Google Scholar 

  48. Mandal S, Van Treuren W, White RA, Eggesbø M, Knight R, Peddada SD. Analysis of composition of microbiomes: a novel method for studying microbial composition. Microb Ecol Heal Dis. 2015;26:27663.

    Article  Google Scholar 

  49. Watts SC, Ritchie SC, Inouye M, Holt KE. FastSpar: rapid and scalable correlation estimation for compositional data. Bioinformatics. 2019;35(6):1064–6.

    Article  CAS  PubMed  Google Scholar 

  50. Friedman J, Alm EJ. Inferring correlation networks from genomic survey data. PLoS Comput Biol. 2012;8(9): e1002687.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  51. Wei T, Simko V. Corrplot: visualization of a correlation matrix. 2016.

    Google Scholar 

Download references


The authors are sincerely grateful to the participants and staff of the TB-CHAMP clinical trial for collecting samples and data. We would like to thank Ms Anneen van Deventer, Dr Sue Purchase and Ms Elize Batist from the Desmond Tutu TB Centre for their valuable assistance. The Medical Research Council Clinical Trials Unit (MRC-CTU), UK, provided the data collected from participants. The authors would like to express their gratitude for the support given by the Division of Medical Microbiology Manuscript Writing Group.


This study was supported by grants funded through the NHLS Research Trust of South Africa (MNF received GRANT004_ 94632 and GRANT004_94679) and the Harry Crossley Foundation (KN received a HCF research grant). KN was supported by a postgraduate study bursary from the Harry Crossley Foundation. JAS is supported by a Clinician Scientist Fellowship jointly funded by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement (MR/R007942/1). AMD is supported by the National Institute of Allergy and Infectious Diseases (NIAID) of the National Institutes of Health (NIH) under Award Number UM1AI06716. The views expressed in this article are those of the authors and do not necessarily represent the official views of the National Institutes of Health. The TB-CHAMP trial is funded by the Joint Global Health Trials Scheme of the Department for International Development, UK (DFID), the Wellcome Trust and The Medical Research Council (MRC UK) (grant MR/M007340/1), and the South African Medical Research Council (SA MRC) Strategic Health Innovation Partnerships (SHIP) (grant S003805). This project was made possible thanks to Unitaid’s funding and support. Unitaid accelerates access to innovative health products and lays the foundations for their scale-up by countries and partners. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Author information

Authors and Affiliations



KN performed the analysis and wrote the main manuscript text. KN, AW, ACH, JAS, AMD and MNF contributed to the conceptualization of the study. All authors reviewed the manuscript. The author(s) read and approved the final manuscript.

Corresponding author

Correspondence to K Nel Van Zyl.

Ethics declarations

Ethics approval and consent to participate

The TB-CHAMP trial was approved by the South African Health Products Regulatory Authority (20160128), local health authorities and the Stellenbosch University Health Research Ethics Committee (SU-HREC, M16/02/009). This sub-study involved participants from the trial’s Cape Town sites only and ethical approval was obtained from the SU-HREC (S18/02/031). Written informed consent was obtained from the parents or legal guardians of the children. The research from both the parent trial and this sub-study were conducted as outlined by Stellenbosch University’s Policy for responsible research conduct in human participants and guidelines set out by the World Health Organization (WHO) and the Declaration of Helsinki.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Nel Van Zyl, K., Whitelaw, A.C., Hesseling, A.C. et al. Fungal diversity in the gut microbiome of young South African children. BMC Microbiol 22, 201 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: