Next Article in Journal
Genome-Wide Association Mapping of Crown and Brown Rust Resistance in Perennial Ryegrass
Next Article in Special Issue
Life Cycle and Phylogeography of True Truffles
Previous Article in Journal
Haplotype-Based Single-Step GWAS for Yearling Temperament in American Angus Cattle
Previous Article in Special Issue
Diversity, Dispersal and Mode of Reproduction of Amanita exitialis in Southern China
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Genetic Diversity and Dispersal of Aspergillus fumigatus in Arctic Soils

Department of Biology, McMaster University, 1280 Main Street West, Hamilton, ON L8S 4K1, Canada
*
Author to whom correspondence should be addressed.
Submission received: 2 December 2021 / Revised: 20 December 2021 / Accepted: 21 December 2021 / Published: 22 December 2021
(This article belongs to the Special Issue Population Genetics of Fungi)

Abstract

:
Aspergillus fumigatus is a saprophytic mold and an opportunistic pathogen with a broad geographic and ecological distribution. A. fumigatus is the most common etiological agent of aspergillosis, affecting over 8,000,000 individuals worldwide. Due to the rising number of infections and increasing reports of resistance to antifungal therapy, there is an urgent need to understand A. fumigatus populations from local to global levels. However, many geographic locations and ecological niches remain understudied, including soil environments from arctic regions. In this study, we isolated 32 and 52 A. fumigatus strains from soils in Iceland and the Northwest Territories of Canada (NWT), respectively. These isolates were genotyped at nine microsatellite loci and the genotypes were compared with each other and with those in other parts of the world. Though significantly differentiated from each other, our analyses revealed that A. fumigatus populations from Iceland and NWT contained evidence for both clonal and sexual reproductions, and shared many alleles with each other and with those collected from across Europe, Asia, and the Americas. Interestingly, we found one triazole-resistant strain containing the TR34 /L98H mutation in the cyp51A gene from NWT. This strain is closely related to a triazole-resistant genotype broadly distributed in India. Together, our results suggest that the northern soil populations of A. fumigatus are significantly influenced by those from other geographic regions.

1. Introduction

Aspergillus fumigatus is a thermotolerant ascomycete mould with a ubiquitous presence around the world. Its primary ecological niche is within decaying plant matter and soil. However, it is also a common opportunistic fungal pathogen capable of infecting both immunocompetent and immunocompromised humans. Over the past three decades, there have been a rising number of reports of antifungal-resistant A. fumigatus strains, which have significantly impacted the treatment of patients with A. fumigatus infections [1,2,3,4,5]. A. fumigatus infections, collectively known as aspergillosis, encompass a range of illnesses from asthma to invasive aspergillosis. Invasive aspergillosis, the most severe form of aspergillosis, is estimated to occur at over 200,000 cases annually worldwide [6]. Depending on the underlying conditions of patients and the effectiveness of antifungal management, the mortality of invasive aspergillosis ranges from 40% to 90%. For both treatment and prophylaxis, triazole antifungals such as itraconazole, voriconazole, posaconazole, and isavuconazole are often used for frontline therapy. Triazole antifungals target the enzyme lanosterol 14α-demethylase encoded by the gene cyp51A. This enzyme is required for the biosynthesis of ergosterol, an essential sterol in the cytoplasmic membrane of fungal cells. Resistance to triazoles is commonly conferred by mutations within the cyp51A gene, inhibiting triazole binding and/or causing overexpression of the enzyme [7]. However, non-cyp51A mediated triazole-resistance mechanisms also exist, such as the overexpression of abc and cdr efflux pumps [8]. Aside from the triazoles, two other antifungal classes, the polyene amphotericin B and the echinocandins, are used for salvage therapy of aspergillosis, especially in the case of triazole-resistance [9].
Over the past 30 years, there have been rising incidences of triazole resistant A. fumigatus infections worldwide, including the identifications of triazole-resistant strains in triazole naïve patients [10,11,12,13]. These results suggest the importance of environmental populations of A. fumigatus to patients and to the clinical populations of this species. Consequently, it is extremely important to understand the environmental populations of A. fumigatus. Indeed, an increasing number of environmental populations from different geographic regions have been surveyed to aid in monitoring drug resistance rates and identifying/tracking resistant A. fumigatus genotypes. The results so far suggest that agricultural use of triazole fungicides can contribute to the development of triazole resistant strains, which subsequently infect patients [14,15,16,17,18]. In these and other studies, a panel of nine highly polymorphic short tandem repeat (STR) loci are often used to analyze the relationships among strains and populations. Based on genotype information at these nine loci on thousands of A. fumigatus isolates, divergent lineages as well as high levels of gene flow between geographic populations have been identified [18,19,20]. For example, Ashu et al. [19] revealed eight genetic populations in their global sample with all eight genetic populations being broadly distributed across multiple countries and continents. However, the samples from Cameroon in central western Africa was genetically unique, different from those analyzed so far from Eurasia, Oceania, and the Americas [18].
Among the currently analyzed geographic populations of A. fumigatus, almost all have come from tropical, subtropical, and temperate regions. Little is known about geographic populations of A. fumigatus from cold climates. In this study, we isolated and analyzed strains of A. fumigatus from two high latitude regions: Iceland and the Northwest Territories (NWT) in Canada. Due to their long-distance separation, we hypothesized that populations of A. fumigatus from these two geographic regions should be genetically different from each other. In addition, due to their unique climates, we expect that these two A. fumigatus geographic populations would be genetically different from those reported from other geographic regions in temperate, subtropical, and tropical climates. Finally, since both geographic regions have very limited agriculture and no history of agricultural fungicide usage, we expected that all isolates from the soil in both regions should be susceptible to the common clinical triazole drugs itraconazole and voriconazole. Our results showed that some of the expectations were met. However, unexpected results were also identified.

2. Materials and Methods

2.1. Soil Collection, A. fumigatus Isolation, and STR Genotyping

Soil samples were obtained from Iceland and the Northwest Territories in Canada. All soil samples were obtained from at least 1 cm below the soil surface. In Iceland, a total of 314 soil samples were collected from the following six sites: Dimmuborgir (60), Thingvellir (51), Skaftafell (52), Myvatn Lake (42), Landbrotalaug (60), and Reykjavik University (49). In Northwest Territories (NWT) in Canada, 220 soil samples were obtained from the following four sites within and near Yellowknife: Yellowknife Downtown (80), 30 km North of Yellowknife (50), 30 km Northwest of Yellowknife (60), and in Yellowknife Old Town (30). The numbers in parenthesis represent the number of soil samples taken from each site.
To isolate A. fumigatus, approximately 0.2 g of soil from each soil sample was put in 1 mL of Sabouraud dextrose (SD) broth containing 50 mg/L chloramphenicol within a 1.5 mL microcentrifuge tube. The soil suspension was mixed through vertexing and incubated at 50 °C for 3 days to select for A. fumigatus growth. For the soil samples that failed the initial attempt, two additional attempts were made, each with a different incubation temperature, at 42 °C and 30 °C, respectively. After incubation, mycelia grown on the surface of the broth from each incubation temperature was then transferred to SD agar containing 50 mg/L chloramphenicol and incubated at 37 °C for 3 days. Colonies resembling A. fumigatus were sub-cultured to a new agar medium and their conidia were harvested with sterile 30% glycerol and stored at −80 °C. To confirm their identity, these putative A. fumigatus isolates were genotyped at the mating type locus using A. fumigatus species-specific and mating-type specific primer pairs [21]. Strains confirmed as A. fumigatus were then analyzed for their genotypes at nine highly polymorphic short tandem repeat (STR) loci specific for A. fumigatus, following protocols previously described by De Valk et al. [22]. The only modification was the fluorescent tag for primers STRAf2C, STRAf 3C, and STRAf 4C that was changed from the TET probe to ATTO550 (IDT) to provide better fragment size detection and scoring. STR fragments at the nine loci were separated using capillary electrophoresis at McMaster University’s MoBix Lab and analyzed using the STR analysis software Osiris [23].

2.2. Allelic and Genotypic Diversities

Allelic diversity was calculated for each of the nine STR loci within both the Iceland and NWT populations using the Excel add-in GenAlEx 6.5 [24,25]. GenAlEx 6.5 was also used to compute unbiased haploid genetic diversity (uh). The uh diversity measures the probability that two individuals will be different within the population using the sum of allele frequencies within the populations and adjusts for differences in sample size. Additionally, the number of effective alleles (Ne) was calculated. To determine whether the Iceland and NWT populations differ significantly in their mean uh and mean Ne, two-tailed Student’s T-tests were conducted.

2.3. Clonality and Recombination

To assess the potential evidence for clonality and/or recombination in these two regional populations, we conducted two tests. The first test was the overall linkage disequilibrium within each of the two geographic populations using the standardized index of association ( r d ¯ ) and the second was the proportion of loci that are phylogenetically compatible. The null hypothesis of the index of association test is random recombination. The function poppr from the R package poppr was used to calculate r d ¯ [26]. Statistical significance was determined through 999 permutations. To account for the influence of clonal genotypes on r d ¯ , both populations were clone-corrected using the poppr function clonecorrect and the r d ¯ was recalculated. The proportion of phylogenetically compatible pairs of loci was calculated using the software MULTILOCUS v1.3b [27]. Two loci are phylogenetically compatible if all observed genotypes can be accounted for by mutation alone without inferring homoplasy or recombination. This test is also called the four-gamete test. Statistical significance was determined through 1000 randomizations.

2.4. Genetic Relationships among Strains from the Two Arctic Populations

To assess the proportion of genetic variation present within and between the two arctic populations, analysis of molecular variance (AMOVA) was conducted using GenAlEx 6.5 [24,25]. In a panmictic population, most of the variance observed will arise within the samples. In contrast, a high-level among-sample variance would suggest the presence of genetically differentiated populations. The null hypothesis for AMOVA is that there is no genetic difference between the populations tested. For genetic distance calculation among strains, the parameter haploid-SSR was used to incorporate the stepwise mutation model. The smaller the difference between allele repeat numbers at each locus, the smaller the genetic distance between the alleles and consequently, their multilocus genotypes (MLGs) would be more similar to each other. Genetic differentiation between the two arctic populations was calculated using ϕpt, an analogue of Fst. Significance was determined through 999 permutations.
To visualize genetic distance between MLGs from Iceland and NWT, a minimum spanning network (MSN) was generated by calculating Bruvo’s genetic distance between strains using the R package poppr [26]. Bruvo’s genetic distance is specific for STR genotypes and incorporates the stepwise mutation model. Genetic distance was calculated for each pair of MLGs, and the resulting matrix is visualized in the MSN. Minimum edge genetic distance was set to 0.05. In addition, a multivariate analysis, Discriminant Analysis of Principal Components (DAPC), implemented by adegenet package in R, was used to cluster MLGs genotypes in relation to their geographic origins [28]. In this analysis, DAPC first transforms the data using principal component analysis (PCA) to reduce the number of variables and allowing retention of the variables that have the greatest contribution to the variation within the dataset. This is followed by discriminate analysis (DA) to cluster the MLGs by optimizing the between-group variation and reducing within-group variation.

2.5. Genetic Relationships between Arctic Samples and Those from Other Regions

To compare our arctic STR dataset with those from other geographic and climatic regions, STR data from 12 populations within Eurasia, North America, and Oceania, were included, totaling 2339 MLGs [19,29]. Specifically, the number of MLGs within these 12 populations were: Belgium (108), Hamilton, Canada (195), Cameroon (51), France (66), Germany (85), India (94), Netherlands (1082), Norway (203), New Zealand (104), Spain (180), Switzerland (70), and USA (101).
In our comparisons, the number of private alleles that were present only in Iceland and NWT populations but were absent in these 12 populations was determined. To analyze the level of differentiation between these two arctic populations and those from other geographic regions, we used the same parameters and AMOVA function as described above. In addition to the overall AMOVA, pairwise population ϕpt was estimated among the 14 geographic populations within this dataset. Statistical significance was determined through 999 permutations. An MSN was generated and DAPC was conducted as described in the above section.

2.6. Triazole Susceptibility Testing

Antifungal susceptibility testing was conducted following the reference protocols described in the Clinical and Laboratory Standards Institute (CLSI) document “M38 Reference Method for Broth Dilution Antifungal Susceptibility Testing of Filamentous Fungi” 3rd edition [30]. Here, two triazoles, itraconazole and voriconazole, were used. Two Candida strains, C. parapsilosis ATCC® 22019 and C. krusei ATCC® 6258, were used as references. Briefly, spore suspensions were adjusted to an optical density between 0.09 and 0.13 at 530 nm. The adjusted spore suspensions were then diluted 50× into RPMI1640 at pH 7 ± 0.1. The two triazole antifungals were dissolved to 200 mg/L in DMSO. A two-fold serial dilution series was conducted for each antifungal to create a 3 to 200 mg/L concentration range in addition to a DMSO only control. Each concentration was then diluted 50× in RPMI 1640. Lastly, 100 μL from both the spore suspension and triazole solutions were added to 96-well cell culture plates to test the susceptibility of each isolate from NWT and Iceland at the final drug concentrations 0.03 to 2 mg/L. The 96-well cell culture plates were then incubated at 35 °C without agitation. The minimum inhibitory concentrations (MIC) for both antifungals were visually determined after 48 h. Strains that grew at 2 mg/L were further tested at the antifungal concentrations between 0.25 to 16 mg/L.

2.7. Cyp51A Sequencing

The cyp51A gene was sequenced for a strain that had an itraconazole and/or voriconazole MIC of 2 or higher to identify putative triazole resistance mutation(s) in this gene. Following the PCR protocol described by Mellado et al. [31], five pairs of primers were used to amplify the cyp51A region and its 5′UTR. Amplified PCR products were sent for Sanger sequencing at McMaster University’s MoBix Lab. To identify mutations within the cyp51A gene, the obtained sequences were aligned to the cyp51A gene from the A. fumigatus reference strain af293 using the program MEGA-X [32].

3. Results

3.1. Isolation Rates of A. fumigatus from Arctic Soils

In total, 52 (23.6%) and 32 (10.2%) A. fumigatus isolates were obtained from 220 NWT and 314 Iceland soils samples, respectively. In both regions, there were considerable variations among sites in their isolation rates of A. fumigatus (Table 1). Following isolation, all isolates were genotyped at nine microsatellite loci.

3.2. Local Genetic Diversity within Iceland and NWT

The STR genotypes were obtained for isolates from both Iceland and NWT A. fumigatus populations. Because of the small sample sizes and unequal sample size distributions among local populations from within both Iceland and NWT (Table 1), we elected to not compare the local populations within each of the two regions but instead focused on comparing the two arctic populations. We found that the two regional populations had similar uh between Iceland and NWT where no significant difference was observed (Table 2; p = 0.36). For the effective number of alleles, the NWT population had a higher value than that from Iceland, but the difference was statistically insignificant (p = 0.13). Similarly, the number of observed alleles was higher in the NWT population but after being adjusted for sample size effects, the difference was statistically nonsignificant (p = 0.08). Together, the results indicated similar levels of allelic and genotypic diversities between the Iceland and NWT soil populations of A. fumigatus.

3.3. Clonality and Recombination within Iceland and Northwest Territories

To investigate evidence for clonality and recombination within the two arctic A. fumigatus populations, the r d ¯ and the proportion of phylogenetically compatible loci were calculated for each population. The r d ¯ values of Iceland and NWT were 0.36 (p = 0.001) and 0.16 (p = 0.001), respectively, rejecting the null hypothesis of random mating. The hypothesis of random mating was rejected even after clonal corrections where only 20 and 40 unique MLGs within Iceland and NWT populations, respectively, were retained for analyses. Specifically, after clone-correction, the r d ¯ values of the Iceland and NWT populations were 0.19 (p = 0.001) and 0.13 (p = 0.001), respectively. However, phylogenetic compatibility analyses revealed that 55.6% (p = 0.001) and 5.6% (p = 0.001) pairs of loci within Iceland and NWT were phylogenetically compatible. The high proportions of phylogenetically incompatible pairs of loci, at 44.4% and 94.4%, respectively, for the Iceland and NWT populations, are consistent with recombination. Overall, these results indicate that both the Iceland and NWT soil populations of A. fumigatus contain signatures of both clonal and recombining modes of reproduction.

3.4. Relationships between Iceland and NWT Samples

We analyzed the relationships among A. fumigatus isolates from Iceland and NWT. Bruvo’s genetic distance was calculated between isolates of both populations and visualized through an MSN (Figure 1). Overall, while we observed some geographic clustering among MLGs, isolates from the two arctic populations were intermixed. For example, two different MLGs were shared between Iceland and NWT populations. One of these two MLGs was shared among six strains with three isolates from each of the two regional populations. The other was shared among four isolates (one from Iceland and three from NWT). A DAPC of the Iceland and NWT samples yielded comparable results, which both showed some geographic clustering as well as overlaps between these two regions (Figure 2). Results from AMOVA indicated that while the majority (86.8%) of the genetic variations were found within the two regional populations, 13.2% of the total observed genetic variance was found between these two populations (ϕpt = 0.132, p = 0.001). Together, these results are consistent with an overall statistically significant genetic differentiation between these two arctic populations.

3.5. Relationship between the Arctic Populations to Additional Global A. fumigatus Populations

To determine how the arctic populations fit in the global context, the allelic and genotypic diversities of 12 other geographic populations totaling 2339 A. fumigatus MLGs were compared to the Iceland and NWT A. fumigatus populations. Within the 2423 MLGs, the number of private alleles present only within Iceland and NWT was 2 and 1 respectively. In Iceland, these two private alleles were #116 and #117 at the 3A locus. In NWT, the private allele was #57 at the 3B locus. Together, the allelic distribution results suggest limited novel genetic diversity in the arctic soil samples.
To determine the genetic relationships between the 2 arctic populations and the 12 other geographic populations, ϕpt was calculated between all pairwise population combinations (Table 3). For this analysis, the dataset was also clone-corrected to remove the influence of multiple clonal MLGs in determining genetic differentiation between populations. Overall, our analyses revealed that the two arctic populations of A. fumigatus were significantly different from most other geographic populations (Table 3). Specifically, before clonal correction, only the Belgian and French populations showed insignificant difference to the Iceland population. After clonal correction, only the Indian and French populations were not significantly differentiated from the Iceland population while the Norwegian and New Zealand populations were not significantly different from the NWT population.
We also analyzed the genotypic similarities among individual MLGs from the 14 geographic populations. An MSN of all 2416 strains was generated to visualize the genetic distance among the MLGs (Figure 3). To highlight the distribution of the arctic isolates, the 12 previously surveyed A. fumigatus populations were coloured white. Overall, while the MLGs from the arctic regions showed some geographic clustering, these arctic genotypes were broadly embedded in the global genotype network. This overall pattern was similarly supported by DAPC analyses (Figure 4) where MLGs from Iceland and NWT clustered with those from Eurasia and the United States but were different from the Hamilton, Ontario, Canada population as well as from the Cameroonian population.

3.6. Susceptibility Testing

Antifungal susceptibility testing using voriconazole and itraconazole was conducted for all isolates from the two arctic populations (Table 4). One isolate resistant to both triazoles was identified within NWT but no resistant strain was identified from Iceland (Figure 5). Our sequence analysis revealed that this triazole-resistant strain NCY6_13_2 had the TR34/L98H mutation within the cyp51A gene. Interestingly, this strain had a near identical genotype at the nine STR loci to a clonal population of triazole-resistant strains discovered from India in 2012, differing by only one repeat unit at the 3A locus and with identical alleles at the other eight STR loci (Table 5) [16]. Furthermore, all the A. fumigatus isolates in this Indian clone contained the identical TR34/L98H mutation within the cyp51A gene as the strain from NWT in our sample.

4. Discussion

In this study, we obtained and genotyped A. fumigatus isolates from two northern climatic regions, Iceland and Northwest Territories, Canada. This is the first study to investigate A. fumigatus population structures at high latitudes and cold climates. A total of 32 (20 unique MLGs) and 52 isolates (40 unique MLGs) were obtained from Iceland and NWT, respectively. Both populations contained abundant allelic and genotypic diversities at the nine microsatellite loci, with evidence of recombination within both populations. However, the observed recombination was not due to random mating. Our analyses demonstrated that while Iceland and NWT shared several identical/similar multilocus genotypes, the two regional populations were significantly differentiated from each other. Similarly, while many alleles were shared, these two arctic soil populations were significantly differentiated from geographic samples from other parts of the world. Interestingly, one triazole resistant strain was identified from NWT and this strain had a near identical MLG as the clonal triazole-resistant population previously identified in India. Below we discuss the implications and relevance of this work to previous studies.
The high level of allelic and genotypic diversity observed within both the Iceland and NWT populations resemble that seen within most regional A. fumigatus populations [18,29,33,34,35]. The high allelic and genotypic diversities within both populations (Table 2) suggested that either these two populations were ancient or that there had been multiple migrations from other regions to these two arctic regions. Indeed, the high allelic and genotypic diversities in both regions suggested no evidence of recent selective sweep in either Iceland or NWT. In comparison, a large and broadly distributed MLG within India was found to be caused by a selective sweep of a highly drug-resistant genotype [16]. However, several MLG were over-represented by two or more isolates in both arctic samples, suggesting that local selection could happen in both populations. At present, due to the relatively small and uneven population sizes among local populations within each of these two regions, it is difficult to determine what selective force(s) might be responsible for their over-representation. The over-representation of certain genotypes is consistent with clonal reproduction of A. fumigatus in nature. However, evidence for recombination was also found in both regional populations. Our results are similar to those reported previously for other regional populations of A. fumigatus [18,29,33,34,35] as well as other human pathogenic and non-pathogenic fungi [36].
Within individual fungal species, both latitudinal differences and geographic distances have been found to correlate with population genetic differences [37,38]. For example, a latitudinal gradient-based differentiation was observed among populations of the chili pepper anthracnose pathogen Colletotrichum acutatum [37]. In the rust fungus Melampsora laricipopulina, a gradual dispersal of genetic elements across the European populations was observed, which led to higher genetic variations between populations at greater geographic distances [38]. In contrast, the significant genetic differentiation of the Canadian and Icelandic populations of M. laricipopulina was likely caused by founder effects [38]. Regarding A. fumigatus, geographic A. fumigatus populations have been shown to be shaped by both historical differentiation and high contemporary gene flow due to anthropogenic influences [19,29]. Our results from the Icelandic and NWT samples of A. fumigatus are consistent with previous findings. Specifically, we observed both allele and genotype-sharing between these two regions but also found an overall significant genetic differentiation both between these two arctic samples as well as between the arctic samples and the geographic populations from other regions of the world.
At present, while the close relationship between the NWT and Norwegian samples might be due to similar selection pressure, the reasons for the low/insignificant level of differentiation between Icelandic A. fumigatus population and those from Belgium, France, and Indian populations as well as between the NWT A. fumigatus population and the New Zealand population are unknown. As previously discussed by Korfanty et al. [29], two types of factors may contribute to gene flow between distant countries. The first type is natural factors such as wind. However, given the geographic distance between most of these areas, while possible, direct wind-aided dispersal of A. fumigatus conidia is unlikely, especially across hemispheres such as between New Zealand and NWT in Canada (~12,694 km). The second type is anthropogenic factors such as human travel and commercial trade. For both Iceland and NWT, the European Union is a major exporter and importer of goods to and from both regions, contributing to 30% of exports from NWT and 53% of Iceland’s total trade in goods in 2019 [39,40]. Similarly, India was the top importer of NWT goods in 2019, where Asian specific exports from NWT began to increase in 2011. Within Iceland in 2019, 682,700 tourists (34.4%) originated from the European Union [41]. Similarly, Yellowknife in NWT has been a popular destination for observing northern lights, which have attracted tourists from around the world, including from many Asian countries such as China and India. However, the number of tourists from all destinations declined substantially since the start of the SARS-Covid-19 pandemic. Given the resilience of A. fumigatus spores to a variety of environmental stressors and their abundance within the atmosphere [42], spores could have followed these anthropogenic routes and migrated to and from Iceland and NWT.
One triazole resistant A. fumigatus strain was identified in the NWT population. This triazole-resistant strain contained the TR34/L98H mutation in the cyp51A gene. This mutation has been reported from environmental and clinical samples from many countries [12,43,44,45], including the selective sweep observed in a widespread Indian genotype of A. fumigatus in 2012 [16]. The resistant strain identified in our study had a near identical MLG to the Indian population. However, how the strain migrated into the NWT region is unclear. A recent long distance dispersal event may have occurred, which introduced this strain to NWT from India, potentially by an anthropogenic event. Unfortunately, the lack of clinical research on A. fumigatus and aspergillosis within both Iceland and Northern regions in Canada makes it difficult to determine the clinical significance of triazole-resistant aspergillosis in both northern regions and the potential relationship between the triazole-resistant strain isolated here from the soil and the clinical strains of A. fumigatus in NWT. Systematic analyses of clinical aspergillosis and susceptibility testing of the strains responsible for such infections in both regions are needed in order to assess the frequency of triazole-resistant aspergillosis and the relationship between environmental and clinical populations of A. fumigatus in both Iceland and Northern Canada.

5. Conclusions

This study investigated the population structure of the opportunistic fungal pathogen A. fumigatus obtained from northern climates above 60° latitude. Prior to this study, limited isolation and genotypic data has been described from soil environments in these northern regions [46,47,48]. Here we characterized A. fumigatus diversity present within Iceland and Northwest Territories in Canada. Our study revealed significant genetic diversity with evidence of both clonal and sexual reproduction within both regional populations. While evidence for gene flow was found between these two regional arctic populations as well as between them and other geographic populations, unique alleles and genotypes were also found in both arctic populations. The finding of a triazole-resistant strain from the arctic with a nearly identical genotype to strains in India demonstrates the complexity of managing drug resistance in A. fumigatus. Our study highlights the importance of both genetically characterizing A. fumigatus populations from a diversity of geographic and ecological niches and identifying the level of triazole susceptibility within these populations to help track the spread of drug-resistant strains. Indeed, the discovery of high genetic diversity in cold environments coupled with the known high-temperature tolerance of A. fumigatus calls for greater effort for understanding the thermo-adaptability of this organism, with implications even for space travel and contamination of Mars from the Mars missions [46,49].

Author Contributions

Conceptualization, J.X.; methodology, G.A.K.; software, G.A.K.; formal analysis, G.A.K.; investigation, G.A.K. and M.D.; resources, H.J., H.Y. and J.X.; writing—original draft preparation, G.A.K.; writing—review and editing, G.A.K., M.D., H.J., H.Y. and J.X.; visualization, G.A.K.; supervision, J.X.; project administration, J.X.; funding acquisition, J.X. and G.A.K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by grants from the Natural Sciences and Engineering Research Council of Canada (Grant No. CRDPJ 474638-14), and by the Institute of Infectious Diseases Research (IIDR) Antibiotic Resistance Initiative and the Faculty of Science’s Global Science Initiative of McMaster University. The APC was funded by IIDR. G.A.K. is supported by MITACS and NSERC PGS-D Scholarships.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Verweij, P.E.; Lucas, J.A.; Arendrup, M.C.; Bowyer, P.; Brinkmann, A.J.F.; Denning, D.W.; Dyer, P.S.; Fisher, M.C.; Geenen, P.L.; Gisi, U.; et al. The one health problem of azole resistance in Aspergillus fumigatus: Current insights and future research agenda. Fungal Biol. Rev. 2020, 34, 202–214. [Google Scholar] [CrossRef]
  2. Zhang, M.; Feng, C.-L.; Chen, F.; He, Q.; Su, X.; Shi, Y. Triazole resistance in Aspergillus fumigatus clinical isolates obtained in Nanjing, China. Chin. Med. J. 2017, 130, 665–668. [Google Scholar] [CrossRef] [PubMed]
  3. Resendiz Sharpe, A.; Lagrou, K.; Meis, J.F.; Chowdhary, A.; Lockhart, S.R.; Verweij, P.E. Triazole resistance surveillance in Aspergillus fumigatus. Med. Mycol. 2018, 56, S83–S92. [Google Scholar] [CrossRef] [Green Version]
  4. Bueid, A.; Howard, S.J.; Moore, C.B.; Richardson, M.D.; Harrison, E.; Bowyer, P.; Denning, D.W. Azole antifungal resistance in Aspergillus fumigatus: 2008 and 2009. J. Antimicrob. Chemother. 2010, 65, 2116–2118. [Google Scholar] [CrossRef] [PubMed]
  5. Abdolrasouli, A.; Petrou, M.A.; Park, H.; Rhodes, J.L.; Rawson, T.M.; Moore, L.S.P.; Donaldson, H.; Holmes, A.H.; Fisher, M.C.; Armstrong-James, D. Surveillance for azole-resistant Aspergillus fumigatus in a centralized diagnostic mycology service, London, United Kingdom, 1998–2017. Front. Microbiol. 2018, 9, 2234. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Bongomin, F.; Gago, S.; Oladele, R.; Denning, D. Global and multi-national prevalence of fungal diseases—Estimate precision. J. Fungi 2017, 3, 57. [Google Scholar] [CrossRef] [PubMed]
  7. Liu, M.; Zheng, N.; Li, D.; Zheng, H.; Zhang, L.; Ge, H.; Liu, W. cyp51A -based mechanism of azole resistance in Aspergillus fumigatus: Illustration by a new 3D structural model of Aspergillus fumigatus CYP51A protein. Med. Mycol. 2016, 54, 400–408. [Google Scholar] [CrossRef] [Green Version]
  8. Moye-Rowley, W.S. Multiple mechanisms contribute to the development of clinically significant azole resistance in Aspergillus fumigatus. Front. Microbiol. 2015, 6, 70. [Google Scholar] [CrossRef]
  9. Dockrell, D.H. Salvage therapy for invasive aspergillosis. J. Antimicrob. Chemother. 2008, 61, i41–i44. [Google Scholar] [CrossRef]
  10. Dauchy, C.; Bautin, N.; Nseir, S.; Reboux, G.; Wintjens, R.; Le Rouzic, O.; Sendid, B.; Viscogliosi, E.; Le Pape, P.; Arendrup, M.C.; et al. Emergence of Aspergillus fumigatus azole resistance in azole-naïve patients with chronic obstructive pulmonary disease and their homes. Indoor Air 2018, 28, 298–306. [Google Scholar] [CrossRef] [Green Version]
  11. Astvad, K.M.T.; Jensen, R.H.; Hassan, T.M.; Mathiasen, E.G.; Thomsen, G.M.; Pedersen, U.G.; Christensen, M.; Hilberg, O.; Arendrup, M.C. First detection of TR46/Y121F/T289A and TR34/L98H alterations in Aspergillus fumigatus isolates from azole-naive patients in Denmark despite negative findings in the environment. Antimicrob. Agents Chemother. 2014, 58, 5096–5101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Berger, S.; El Chazli, Y.; Babu, A.F.; Coste, A.T. Azole Resistance in Aspergillus fumigatus: A consequence of antifungal use in agriculture? Front. Microbiol. 2017, 8, 1024. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Pelaez, T.; Monteiro, M.C.; Garcia-Rubio, R.; Bouza, E.; Gomez-Lopez, A.; Mellado, E. First detection of Aspergillus fumigatus azole-resistant strain due to Cyp51A TR46/Y121F/T289A in an azole-naive patient in Spain. New Microbes New Infect. 2015, 6, 33–34. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Stensvold, C.R.; Jørgensen, L.N.; Arendrup, M.C. Azole-Resistant invasive aspergillosis: Relationship to agriculture. Curr. Fungal Infect. Rep. 2012, 6, 178–191. [Google Scholar] [CrossRef]
  15. Ren, J.; Jin, X.; Zhang, Q.; Zheng, Y.; Lin, D.; Yu, Y. Fungicides induced triazole-resistance in Aspergillus fumigatus associated with mutations of TR46/Y121F/T289A and its appearance in agricultural fields. J. Hazard. Mater. 2017, 326, 54–60. [Google Scholar] [CrossRef]
  16. Chowdhary, A.; Kathuria, S.; Xu, J.; Sharma, C.; Sundar, G.; Singh, P.K.; Gaur, S.N.; Hagen, F.; Klaassen, C.H.; Meis, J.F. Clonal expansion and emergence of environmental multiple-triazole-resistant Aspergillus fumigatus strains carrying the TR34/L98H mutations in the cyp51A gene in India. PLoS ONE 2012, 7, e52871. [Google Scholar] [CrossRef]
  17. Chowdhary, A.; Kathuria, S.; Xu, J.; Meis, J.F. Emergence of Azole-Resistant Aspergillus fumigatus Strains due to agricultural azole use creates an increasing threat to human health. PLoS Pathog. 2013, 9, e1003633. [Google Scholar] [CrossRef]
  18. Ashu, E.E.; Korfanty, G.A.; Xu, J. Evidence of unique genetic diversity in Aspergillus fumigatus isolates from Cameroon. Mycoses 2017, 60, 739–748. [Google Scholar] [CrossRef]
  19. Ashu, E.E.; Hagen, F.; Chowdhary, A.; Meis, J.F.; Xu, J. Global population genetic analysis of Aspergillus fumigatus. mSphere 2017, 2, e00019-17. [Google Scholar] [CrossRef] [Green Version]
  20. Sewell, T.R.; Zhu, J.; Rhodes, J.; Hagen, F.; Meis, J.F.; Fisher, M.C.; Jombart, T. Nonrandom distribution of azole resistance across the global population of Aspergillus fumigatus. MBio 2019, 10, e00392-19. [Google Scholar] [CrossRef] [Green Version]
  21. Paoletti, M.; Rydholm, C.; Schwier, E.U.; Anderson, M.J.; Szakacs, G.; Lutzoni, F.; Debeaupuis, J.-P.; Latgé, J.-P.; Denning, D.W.; Dyer, P.S. Evidence for sexuality in the opportunistic fungal pathogen Aspergillus fumigatus. Curr. Biol. 2005, 15, 1242–1248. [Google Scholar] [CrossRef] [PubMed]
  22. de Valk, H.A.; Meis, J.F.G.M.; Curfs, I.M.; Muehlethaler, K.; Mouton, J.W.; Klaassen, C.H.W. Use of a novel panel of nine short tandem repeats for exact and high-resolution fingerprinting of Aspergillus fumigatus isolates. J. Clin. Microbiol. 2005, 43, 4112–4120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Goor, R.M.; Neall, L.F.; Hoffman, D.; Sherry, S.T. A Mathematical approach to the analysis of multiplex DNA profiles. Bull. Math. Biol. 2011, 73, 1909. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Peakall, R.; Smouse, P.E. genalex 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes 2006, 6, 288–295. [Google Scholar] [CrossRef]
  25. Peakall, R.; Smouse, P.E. GenAlEx 6.5: Genetic analysis in Excel. Population genetic software for teaching and research—An update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef] [Green Version]
  26. Kamvar, Z.N.; Tabima, J.F.; Grünwald, N.J. Poppr: An R package for genetic analysis of populations with clonal, partially clonal, and/or sexual reproduction. PeerJ 2014, 2, e281. [Google Scholar] [CrossRef] [Green Version]
  27. Agapow, P.-M.; Burt, A. Indices of multilocus linkage disequilibrium. Mol. Ecol. Notes 2001, 1, 101–102. [Google Scholar] [CrossRef]
  28. Jombart, T.; Ahmed, I. adegenet 1.3-1: New tools for the analysis of genome-wide SNP data. Bioinformatics 2011, 27, 3070–3071. [Google Scholar] [CrossRef] [Green Version]
  29. Korfanty, G.A.; Teng, L.; Pum, N.; Xu, J. Contemporary gene flow is a major force shaping the Aspergillus fumigatus population in Auckland, New Zealand. Mycopathologia 2019, 184, 479–492. [Google Scholar] [CrossRef]
  30. Alexander, B.D.; Procop, G.W.; Dufresne, P.; Espinel-Ingroff, A.; Fuller, J.; Ghannoum, M.A.; Hanson, K.E.; Holliday, D.; Holliday, N.M.; Ostrosky-Zeichner, L.; et al. Reference Method for Broth Dilution Antifungal Susceptibility Testing of Filamentous Fungi, 3rd ed.; Clinical and Laboratory Standards Institute: Wayne, PA, USA, 2017; Volume 37. [Google Scholar]
  31. Mellado, E.; Diaz-Guerra, T.M.; Cuenca-Estrella, M.; Rodriguez-Tudela, J.L. Identification of two different 14-α sterol demethylase- related genes (cyp51A and cyp51B) in Aspergillus fumigatus and other Aspergillus species. J. Clin. Microbiol. 2001, 39, 2431–2438. [Google Scholar] [CrossRef] [Green Version]
  32. Kumar, S.; Stecher, G.; Li, M.; Knyaz, C.; Tamura, K. MEGA X: Molecular evolutionary genetics analysis across computing platforms. Mol. Biol. Evol. 2018, 35, 1547. [Google Scholar] [CrossRef] [PubMed]
  33. Ashu, E.E.; Kim, G.Y.; Roy-Gayos, P.; Dong, K.; Forsythe, A.; Giglio, V.; Korfanty, G.; Yamamura, D.; Xu, J. Limited evidence of fungicide-driven triazole-resistant Aspergillus fumigatus in Hamilton, Canada. Can. J. Microbiol. 2018, 64, 119–130. [Google Scholar] [CrossRef] [Green Version]
  34. Zhou, D.; Wang, R.; Li, X.; Peng, B.; Yang, G.; Zhang, K.Q.; Zhang, Y.; Xu, J. Genetic diversity and azole resistance among natural Aspergillus fumigatus populations in Yunnan, China. Microb. Ecol. 2021, 7, 1–17. [Google Scholar] [CrossRef] [PubMed]
  35. Zhou, D.; Korfanty, G.A.; Mo, M.; Wang, R.; Li, X.; Li, H.; Li, S.; Wu, J.-Y.; Zhang, K.-Q.; Zhang, Y.; et al. Extensive genetic diversity and widespread azole resistance in greenhouse populations of Aspergillus fumigatus in Yunnan, China. mSphere 2021, 6, e00066-21. [Google Scholar] [CrossRef] [PubMed]
  36. Xu, J. Is Natural population of Candida tropicalis sexual, parasexual, and/or asexual? Front. Cell. Infect. Microbiol. 2021, 11, 751676. [Google Scholar] [CrossRef]
  37. Diao, Y.; Zhang, C.; Xu, J.; Lin, D.; Liu, L.; Mtung’e, O.G.; Liu, X. Genetic differentiation and recombination among geographic populations of the fungal pathogen Colletotrichum truncatum from chili peppers in China. Evol. Appl. 2015, 8, 108–118. [Google Scholar] [CrossRef]
  38. Barrès, B.; Halkett, F.; Dutech, C.; Andrieux, A.; Pinon, J.; Frey, P. Genetic structure of the poplar rust fungus Melampsora larici-populina: Evidence for isolation by distance in Europe and recent founder effects overseas. Infect. Genet. Evol. 2008, 8, 577–587. [Google Scholar] [CrossRef]
  39. van den Berg, R. Library of Parliament. Available online: https://lop.parl.ca/sites/PublicWebsite/default/en_CA/ResearchPublications/TradeAndInvestment/2020517E (accessed on 21 November 2021).
  40. European Commission. Available online: https://ec.europa.eu/trade/policy/countries-and-regions/countries/iceland/ (accessed on 21 November 2021).
  41. Ferdamalastofa, Iceland Tourist Board. Available online: https://www.ferdamalastofa.is/en/recearch-and-statistics/numbers-of-foreign-visitors (accessed on 21 November 2021).
  42. Kwon-Chung, K.J.; Sugui, J.A. Aspergillus fumigatus—What makes the species a ubiquitous human fungal pathogen? PLoS Pathog. 2013, 9, e1003743. [Google Scholar] [CrossRef]
  43. Wu, C.-J.; Wang, H.-C.; Lee, J.-C.; Lo, H.-J.; Dai, C.-T.; Chou, P.-H.; Ko, W.-C.; Chen, Y.-C. Azole-resistant Aspergillus fumigatus isolates carrying TR 34 /L98H mutations in Taiwan. Mycoses 2015, 58, 544–549. [Google Scholar] [CrossRef]
  44. Etienne, K.A.; Berkow, E.L.; Gade, L.; Nunnally, N.; Lockhart, S.R.; Beer, K.; King Jordan, I.; Rishishwar, L.; Litvintseva, A.P. Genomic diversity of azole-resistant Aspergillus fumigatus in the United States. MBio 2021, 12, e01803-21. [Google Scholar] [CrossRef]
  45. Fan, Y.; Wang, Y.; Xu, J. Comparative genome sequence analyses of geographic samples of Aspergillus fumigatus—relevance for Amphotericin B resistance. Microorganisms 2020, 8, 1673. [Google Scholar] [CrossRef]
  46. Schuerger, A.C.; Lee, P. Microbial ecology of a crewed rover traverse in the arctic: Low microbial dispersal and implications for planetary protection on human mars missions. Astrobiology 2015, 15, 478–491. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Land, C.J.; Sostaric, B.; Fuchs, R.; Lundstrom, H.; Hult, K. Intratracheal exposure of rats to Aspergillus fumigatus spores isolated from sawmills in Sweden. Appl. Environ. Microbiol. 1989, 55, 2856–2860. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Viegas, C.; Almeida, B.; Aranha Caetano, L.; Afanou, A.; Straumfors, A.; Veríssimo, C.; Gonçalves, P.; Sabino, R. Algorithm to assess the presence of Aspergillus fumigatus resistant strains: The case of Norwegian sawmills. Int. J. Environ. Health Res. 2020, 1–9. [Google Scholar] [CrossRef] [PubMed]
  49. Wong Sak Hoi, J.; Beau, R.; Latgé, J.P. A novel dehydrin-like protein from Aspergillus fumigatus regulates freezing tolerance. Fungal Genet. Biol. 2012, 49, 210–216. [Google Scholar] [CrossRef]
Figure 1. Minimum-spanning network showing the genetic relationship between MLGs of A. fumigatus from Iceland and Northwest Territory in Canada. The genetic distance between MLGs was calculated using Bruvo’s genetic distance from the nine microsatellite loci that incorporates the stepwise mutation model. Each node represents one or more identical MLGs, 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.
Figure 1. Minimum-spanning network showing the genetic relationship between MLGs of A. fumigatus from Iceland and Northwest Territory in Canada. The genetic distance between MLGs was calculated using Bruvo’s genetic distance from the nine microsatellite loci that incorporates the stepwise mutation model. Each node represents one or more identical MLGs, 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.
Genes 13 00019 g001
Figure 2. Discriminant analysis of principal components (DAPC) of Iceland and NWT A. fumigatus populations representing the first discriminant function in an individual density plot. Isolates were genotyped at nine microsatellite loci and clone corrected, totaling 60 multilocus genotypes. NWT—Northwest Territories, ISL—Iceland.
Figure 2. Discriminant analysis of principal components (DAPC) of Iceland and NWT A. fumigatus populations representing the first discriminant function in an individual density plot. Isolates were genotyped at nine microsatellite loci and clone corrected, totaling 60 multilocus genotypes. NWT—Northwest Territories, ISL—Iceland.
Genes 13 00019 g002
Figure 3. Minimum-spanning network showing the genetic relationship between MLGs from Iceland and Northwest Territory to 12 other A. fumigatus geographic populations. The genetic distance between MLGs was calculated using Bruvo’s genetic distance based on the nine microsatellite loci that incorporates the stepwise mutation model. Each node represents one or more identical MLGs, where the node size corresponds to the number of identical MLGs. Nodes that are more genetically similar have darker and thicker edges, whereas nodes genetically distant have lighter and thinner edges. 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, and USA—United States.
Figure 3. Minimum-spanning network showing the genetic relationship between MLGs from Iceland and Northwest Territory to 12 other A. fumigatus geographic populations. The genetic distance between MLGs was calculated using Bruvo’s genetic distance based on the nine microsatellite loci that incorporates the stepwise mutation model. Each node represents one or more identical MLGs, where the node size corresponds to the number of identical MLGs. Nodes that are more genetically similar have darker and thicker edges, whereas nodes genetically distant have lighter and thinner edges. 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, and USA—United States.
Genes 13 00019 g003
Figure 4. 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 1703 unique multilocus genotypes. Populations were clustered according to population of origin. 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, and USA—United States.
Figure 4. 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 1703 unique multilocus genotypes. Populations were clustered according to population of origin. 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, and USA—United States.
Genes 13 00019 g004
Figure 5. Radial growth and voriconazole susceptibility of triazole resistant A. fumigatus strain NCY6_13_2. (A) Growth morphology of A. fumigatus strain NCY6_13_2 on Malt extract agar after 5-day incubation at 35°C. (B) Voriconazole susceptibility of NCY6_13_2. Two control strains, Candida parasilosis ATCC® 22019 and C. krusei ATCC® 6258, were used to verify the susceptibility results. A negative control that was not inoculated with any microbial strain was included in the bottom row. Susceptibility endpoints are circled in red for each strain.
Figure 5. Radial growth and voriconazole susceptibility of triazole resistant A. fumigatus strain NCY6_13_2. (A) Growth morphology of A. fumigatus strain NCY6_13_2 on Malt extract agar after 5-day incubation at 35°C. (B) Voriconazole susceptibility of NCY6_13_2. Two control strains, Candida parasilosis ATCC® 22019 and C. krusei ATCC® 6258, were used to verify the susceptibility results. A negative control that was not inoculated with any microbial strain was included in the bottom row. Susceptibility endpoints are circled in red for each strain.
Genes 13 00019 g005
Table 1. Soil samples and A. fumigatus isolates obtained from Iceland and Northwest Territories (NWT).
Table 1. Soil samples and A. fumigatus isolates obtained from Iceland and Northwest Territories (NWT).
Country/TerritoryLocal SiteNumber of Soil SamplesNumber of Isolates (MAT1-1:MAT1-2) *
IcelandDimmuborgir601 (0:1)
Thingvellir518 (3:5)
Skaftafell 526 (4:2)
Myvatn Lake4210 (7:3)
Landbrotalaug601 *
Reykjavik University 496 (4:2)
Total31232 (18:13) *
Canada/NWTYellowknife Downtown8043 (24:17) *
30 km North of Yellowknife500
30 km Northwest of Yellowknife603 (0:3)
Yellowknife Old Town 306 (1:5)
Total22052 (25:25) *
* 1 isolate from Iceland and 2 from NWT were unable to be genotyped at the MAT locus.
Table 2. Number of alleles and allelic diversity for the nine microsatellite loci of 32 and 52 Aspergillus fumigatus isolates from Iceland and Northwest Territories, respectively.
Table 2. Number of alleles and allelic diversity for the nine microsatellite loci of 32 and 52 Aspergillus fumigatus isolates from Iceland and Northwest Territories, respectively.
PopulationLocusN 1Na 2Ne 3Uh 4
Iceland2A3274.920.82
2B3254.700.81
2C3286.100.86
3A32136.020.86
3B3264.660.81
3C32106.240.87
4A2952.870.68
4B3262.940.68
4C30105.630.85
Mean31.447.784.900.81
Standard Error0.380.910.430.03
Northwest Territories2A51144.430.79
2B52103.990.76
2C52116.930.87
3A501913.740.95
3B50114.630.80
3C501912.890.94
4A5293.250.71
4B5283.220.70
4C52145.880.85
Mean51.2212.786.550.82
Standard Error0.321.351.340.03
1 N = sample size; 2 Na = number of unique alleles; 3 Ne = number of effective alleles = 1/(Σpi2) where pi is the frequency of the ith allele for the population; 4 uh = unbiased diversity = (N/(N−1)) × (1−Σpi2).
Table 3. Pairwise differentiations between Iceland and Northwest Territories A. fumigatus populations to each other and 12 A. fumigatus populations from Eurasia, Oceania, and North America. NWT/ISL row represents the Iceland or Northwest Territories population being compared to the other. 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, and USA—United States.
Table 3. Pairwise differentiations between Iceland and Northwest Territories A. fumigatus populations to each other and 12 A. fumigatus populations from Eurasia, Oceania, and North America. NWT/ISL row represents the Iceland or Northwest Territories population being compared to the other. 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, and USA—United States.
Pairwise ϕpt
IcelandNorthwest Territories
CountryOriginal DatasetClone CorrectedOriginal DatasetClone Corrected
NWT/ISL0.132 **0.083 *0.132 **0.083 *
CMR0.580 ***0.629 ***0.598 ***0.588 ***
CAN0.412 ***0.439 ***0.337 ***0.327 ***
BEL0.0200.048 *0.178 ***0.041 **
FRA0.0280.0120.070 **0.050 *
DEU0.170 ***0.192 ***0.285 ***0.206 ***
IND0.130 ***0.0520.213 ***0.056 *
NLD0.054 **0.069 *0.088 ***0.062 **
NOR0.152 ***0.050 *0.037 **0.014
NZL0.067 **0.123 **0.051 **0.024
ESP0.082 **0.074 *0.122 ***0.069 ***
CHE0.089 **0.100 **0.063 ***0.036 *
USA0.162 ***0.146 **0.129 ***0.051 **
* p-value = 0.05, ** p-value = 0.01, *** p-value = 0.001.
Table 4. Triazole antifungal susceptibility of the 32 and 52 A. fumigatus isolates from Iceland and NWT, respectively, to itraconazole and voriconazole. Numbers in the table refer to the number of strains with the respective minimum inhibitory concentrations for each of the two drugs in each geographic population.
Table 4. Triazole antifungal susceptibility of the 32 and 52 A. fumigatus isolates from Iceland and NWT, respectively, to itraconazole and voriconazole. Numbers in the table refer to the number of strains with the respective minimum inhibitory concentrations for each of the two drugs in each geographic population.
Country/TerritoryTriazoleMinimum Inhibitory Concentration (mg/L)
0.1250.250.51216
IcelandItraconazole0024800
Voriconazole2291000
Canada/NWTItraconazole00351601
Voriconazole02724001
Table 5. Short tandem repeat genotype of the A. fumigatus strain with the TR34/L98H mutation in the triazole-target gene cyp51A obtained from Northwest Territories in Canada. Alleles at eight of the nine STR loci in this strain were identical to those of the large Indian triazole-resistant clonal population as reported by Chowdhary et al. [16].
Table 5. Short tandem repeat genotype of the A. fumigatus strain with the TR34/L98H mutation in the triazole-target gene cyp51A obtained from Northwest Territories in Canada. Alleles at eight of the nine STR loci in this strain were identical to those of the large Indian triazole-resistant clonal population as reported by Chowdhary et al. [16].
StrainRegion2A2B2C3A3B3C4A4B4C
NCY6-13_2NWT142093181081028
1042/09India142093191081028
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Korfanty, G.A.; Dixon, M.; Jia, H.; Yoell, H.; Xu, J. Genetic Diversity and Dispersal of Aspergillus fumigatus in Arctic Soils. Genes 2022, 13, 19. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13010019

AMA Style

Korfanty GA, Dixon M, Jia H, Yoell H, Xu J. Genetic Diversity and Dispersal of Aspergillus fumigatus in Arctic Soils. Genes. 2022; 13(1):19. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13010019

Chicago/Turabian Style

Korfanty, Gregory A., Mykaelah Dixon, Haoran Jia, Heather Yoell, and Jianping Xu. 2022. "Genetic Diversity and Dispersal of Aspergillus fumigatus in Arctic Soils" Genes 13, no. 1: 19. https://0-doi-org.brum.beds.ac.uk/10.3390/genes13010019

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop