Yeast isolation from soil
The above yeast isolation protocol was implemented to culture yeasts from soil samples originating from 53 locations in nine countries9,12. In total, 1,473 yeast strains were isolated from 3,826 soil samples. Given the different climatic conditions of the nine originating countries, the best incubation temperature for each country was determined based on its mean annual temperature (Table 1). Given the slower yeast growth at 14 °C, soil samples from Iceland were incubated on a roller drum for an additional 48 h (96-120 h total). Following 2-3 days of incubation on solid medium, microbial growth was visible on plates for all samples (Figure 1A). A randomly selected colony for each yeast morphology present on a plate was streaked for single colonies on fresh plates (Figure 1B).
The rate of successful yeast isolation differed between countries (Table 1). For example, 97 yeast strains were cultured from 562 Saudi Arabian soil samples (17.3%), whereas 261 yeasts were cultured from 300 Canadian soil samples (87%). A rarefaction analysis should be performed to determine if sufficient soil sampling has been conducted to obtain accurate estimates of the true yeast diversity in sampled locations. An example rarefaction analysis for each country was preformed where the Shannon diversity index was used as a measure of soil yeast diversity (Figure 2). The resulting rarefaction curves approached the saturation asymptote indicating that additional sampling was not likely to have yielded more yeast diversity.
Sequencing of the ITS region revealed that the 1,473 yeast isolates can be categorized into 134 distinct species. This included 90 known species and 44 potentially novel species. We applied a mixed-effects model to identify predictors of global culturable soil yeast diversity, as quantified by the Shannon diversity index31. In this model, mean annual precipitation, mean annual temperature, elevation, and distance to equator of the sampling locations were fitted as fixed effects, while sampling country was set as a random effect9. This model identified mean annual precipitation to be significantly correlated with the Shannon diversity index (p = 0.012), while no significant correlations were found between other predictor variables and the Shannon diversity index9.
To compare this culture-based protocol to culture-independent methods, we compared these findings to a previous study by Tedersoo and colleagues that investigated global diversity of soil fungi using metagenomics2. Tedersoo and colleagues performed high-throughput sequencing of the ITS region on DNA directly extracted from soil samples of 39 countries, four of which, namely Cameroon, Canada, China, and New Zealand, were also sampled in this study. We performed BLAST searches to identify ITS sequences present in both datasets32, and found 26% of the ITS sequences had a significant match (>98.41% nucleotide identity) in the metagenomics dataset. However, <3% matched a fungal sequence from the same country, while the remaining 23% matched a fungal sequence found in a different country9. We were more successful than the metagenomics study in annotating the ITS sequences with yeast species identity. Tedersoo and colleagues reported a total of 50,589 fungal operational taxonomic units (OTUs), with the number of yeast OTUs not known. Of these fungal OTUs, 33% were merely annotated as environmental_sequence (724, 1.4%), uncultured_soil_fungus (2,405, 4.8%), uncultured_ectomycorrhizal_fungus (1,407, 2.8%), or uncultured_fungus (11,898, 23.5%).
Aspergillus fumigatus isolation from soil
After 3 days of soil incubation at 50 °C in 1 mL of SDB in a microcentrifuge tube, mycelia may be found growing within the SDB, on the SDB to air boundary, and/or on the inside walls. Aspergillus fumigatus mycelia are typically found growing on the SDB to air boundary but can also be harvested from just below the SDB surface. After this step, the mycelia are transferred to MEA and grown for 3 days at 37 °C. Mycelial growth on plates may only form the green suede morphology typical of A. fumigatus (Figure 3A). More commonly, other thermotolerant molds may be present within the same soil and grow with A. fumigatus on the same plate or by themselves (Figure3B-D). Aspergillus fumigatus isolation rates may vary between geographic locations, similar to what has been found for soil yeasts. For example, within Vancouver, Canada, 251 isolates were obtained through this method from 540 soil samples (46.5%) (unpublished results). By contrast, within Cameroon, 51 A. fumigatus isolates were harvested from 495 soil samples (10.3%)33.
A. fumigatus conidiophore structure can be viewed and identified using light microscopy. Conidiophores have a ball on stick morphology (Figure 4). Additionally, A. fumigatus conidiophores are uniseriate, where the phialides, attached to the conidia chains, are attached directly to the spherical vesicle. In other Aspergillus species, conidiophores are biseriate, where the phialides are connected to metulae attached to the vesicle.
Several software programs are available to perform fragment analysis of the raw data from capillary electrophoresis, converting the capillary electrophoresis spectrum to fragment sizes. Figure 5 is a chromatogram generated by the program Osiris. The three channels visualizing the fragment lengths of A. fumigatus dinucleotide (2A, 2B, and 2C), trinucleotide (3A, 3B, and 3C), and tetranucleotide (4A, 4B, and 4C) STR loci are shown. In addition, PCR artifacts that occur if the sample DNA concentration during PCR is too high are also shown. The highest peaks of each color represent the fragment sizes of the three STR loci in this plot.
The genetic variation present between A. fumigatus STR genotypes can be visualized using the R package poppr. The R script bruvo.msn creates a pairwise Bruvo's genetic distances matrix between each pair of genotypes. This matrix is then used to generate a minimum spanning network (MSN). A high-quality MSN can then be generated using the R script plot_poppr_msn or imsn (Figure 6). A discriminatory analysis of principal components (DAPC) is another method to visualize the genetic relationships among strains (Figure 7). The R script dapc in the R package adegenet is used to perform DAPC and can be used with known or unknown group priors. With unknown group priors, it uses K-means clustering to identify the likely number of groups of individuals.

Figure 1: Yeast isolation from soil samples. (A) Microbial growth is visible on solid agar following 2-3 days of incubation. Yeast colonies can be seen interspersed among other fungal/bacterial colonies. One representative yeast colony for each morphological type on a plate is selected to streak for single colonies. (B) Three streaks are performed to obtain single colonies of yeast isolates. One plate was used to obtain each soil isolate. Please click here to view a larger version of this figure.

Figure 2: Rarefaction analysis of soil sampling. In each of the nine sampled countries, soil yeast diversity, as quantified by the Shannon diversity index, approaches saturation as the number of soil sample increases. This figure was adapted from 9. Please click here to view a larger version of this figure.

Figure 3: Aspergillus fumigatus morphology on Sabouraud dextroseagar. (A) Green suede-like A. fumigatus growth morphology. (B-D) A. fumigatus growth with other thermophilic molds present within the soil sample. A. fumigatus conidiation is visibly reduced in B and D. Please click here to view a larger version of this figure.

Figure 4: Photo of an Aspergillus fumigatus conidiophore under a light microscope. Scale bar = 10 µm. Please click here to view a larger version of this figure.

Figure 5: Osiris output of nine microsatellite loci from an Aspergillus fumigatus strain. The outputs correspond to the (A) Dinucleotide (2A, 2B, and 2C), (B) trinucleotide (3A, 3B, and 3C), and (C) tetranucleotide (4A, 4B, and 4C) STR loci. Forward primer 5` fluorescent labels are: A = 6-FAM; B = HEX; C = ATTO550. Multiplex reaction products were diluted 30x prior to capillary electrophoresis. The LIZ600 dye standard was used during capillary electrophoresis. Raw data was analyzed using the peak analysis software Osiris identifying potential peaks between 60 and 400 bp. Several PCR artifacts are off-scale peaks that cause fluorescence bleeding, stutter peaks, and N-1 peaks (C). Abbreviations: 6-FAM = 6-carboxyfluorescein; HEX = hexachlorofluorescein; STR = short tandem repeats; RFU = relative fluorescence units; BPS = base pairs. Please click here to view a larger version of this figure.

Figure 6: Minimum-spanning network showing the genetic relationship between MLGs of A. fumigatus from Iceland and Northwest Territories in Canada. Each node represents one MLG, where node size corresponds to the number of strains for each MLG. Nodes that are more genetically similar have darker and thicker edges, whereas nodes genetically distant have lighter and thinner edges. This figure was adapted from 34. Abbreviations: ISL = Iceland; NWT = Northwest Territories in Canada; MLG = multilocus genotype. Please click here to view a larger version of this figure.

Figure 7: Genetic clustering using Discriminant Analysis of Principal Components (DAPC) of Iceland, NWT, Eurasian, North American, and Oceanian A. fumigatus populations. Isolates were genotyped at nine microsatellite loci and clone-corrected, totaling 1,703 unique multilocus genotypes. Genotypes were colored according to geographic origins. This figure was adapted from 34. Abbreviations: DAPC = discriminant analysis of principal components; NWT = Northwest Territories; ISL = Iceland; CMR = Cameroon; CAN = Hamilton, Ontario, Canada; BEL = Belgium; FRA = France; DEU = Germany; IND = India; NLD = Netherlands; NOR = Norway; NZL = New Zealand; ESP = Spain; CHE = Switzerland; USA = United States. Please click here to view a larger version of this figure.
| Country | Incubation temperature (°C) | Soil samples | Yeast isolates | Known species/ Novel species |
| Cameroon | 30 | 493 | 110 | 10/9 |
| Canada | 23 | 300 | 261 | 34/12 |
| China | 23 | 340 | 230 | 23/5 |
| Costa Rica | 30 | 388 | 95 | 20/2 |
| France | 23 | 327 | 175 | 12/2 |
| Iceland | 14 | 316 | 211 | 11/0 |
| New Zealand | 23 | 610 | 155 | 14/4 |
| Peru | 23 | 490 | 139 | 30/9 |
| Saudi Arabia | 30 | 562 | 97 | 8/1 |
| Total | | 3826 | 1473 | 90/44 |
Table 1: Soil yeast isolation from nine countries in six continents. Incubation temperature for soil samples from each country was determined based on its mean annual temperature. Results presented here are adapted from 9, 12.