Genome-Wide Association Analysis for Triazole Resistance in Aspergillus fumigatus

Azole resistance mechanisms and population structure of the human pathogen Aspergillus fumigatus on retail plant products

Aspergillus fumigatus is a ubiquitous saprotroph and human-pathogenic fungus that is life-threatening to the immunocompromised. Triazole-resistant A. fumigatus was found in patients without prior treatment with azoles, leading researchers to conclude that resistance had developed in agricultural environments where azoles are used against plant pathogens. Previous studies have documented azole-resistant A. fumigatus across agricultural environments, but few have looked at retail plant products. Our objectives were to determine if azole-resistant A. fumigatus is prevalent in retail plant products produced in the United States (U.S.), as well as to identify the resistance mechanism(s) and population genetic structure of these isolates. Five hundred twenty-five isolates were collected from retail plant products and screened for azole resistance. Twenty-four isolates collected from compost, soil, flower bulbs, and raw peanuts were pan-azole resistant. These isolates had the TR34/L98H, TR46/Y121F/T289A, G448S, and H147Y cyp51A alleles, all known to underly pan-azole resistance, as well as WT alleles, suggesting that non-cyp51A mechanisms contribute to pan-azole resistance in these isolates. Minimum spanning networks showed two lineages containing isolates with TR alleles or the F46Y/M172V/E427K allele, and discriminant analysis of principle components identified three primary clusters. This is consistent with previous studies detecting three clades of A. fumigatus and identifying pan-azole-resistant isolates with TR alleles in a single clade. We found pan-azole resistance in U.S. retail plant products, particularly compost and flower bulbs, which indicates a risk of exposure to these products for susceptible populations and that highly resistant isolates are likely distributed worldwide on these products.IMPORTANCEAspergillus fumigatus has recently been designated as a critical fungal pathogen by the World Health Organization. It is most deadly to people with compromised immune systems, and with the emergence of antifungal resistance to multiple azole drugs, this disease carries a nearly 100% fatality rate without treatment or if isolates are resistant to the drugs used to treat the disease. It is important to determine the relatedness and origins of resistant A. fumigatus isolates in the environment, including plant-based retail products, so that factors promoting the development and propagation of resistant isolates can be identified.

Pan-azole- and multi-fungicide-resistant Aspergillus fumigatus is widespread in the United States

Aspergillus fumigatus is an important global fungal pathogen of humans. Azole drugs are among the most effective treatments for A. fumigatus infection. Azoles are also widely used in agriculture as fungicides against fungal pathogens of crops. Azole-resistant A. fumigatus has been increasing in Europe and Asia for two decades where clinical resistance is thought to be driven by agricultural use of azole fungicides. The most prevalent mechanisms of azole resistance in A. fumigatus are tandem repeats (TR) in the cyp51A promoter coupled with mutations in the coding region which result in resistance to multiple azole drugs (pan-azole resistance). Azole-resistant A. fumigatus has been isolated from patients in the United States (U.S.), but little is known about its environmental distribution. To better understand the distribution of azole-resistant A. fumigatus in the U.S., we collected isolates from agricultural sites in eight states and tested 202 isolates for sensitivity to azoles. We found azole-resistant A. fumigatus in agricultural environments in seven states showing that it is widespread in the U.S. We sequenced environmental isolates representing the range of U.S. sample sites and compared them with publicly available environmental worldwide isolates in phylogenetic, principal component, and ADMIXTURE analyses. We found worldwide isolates fell into three clades, and TR-based pan-azole resistance was largely in a single clade that was strongly associated with resistance to multiple agricultural fungicides. We also found high levels of gene flow indicating recombination between clades highlighting the potential for azole-resistance to continue spreading in the U.S.IMPORTANCEAspergillus fumigatus is a fungal pathogen of humans that causes over 250,000 invasive infections each year. It is found in soils, plant debris, and compost. Azoles are the first line of defense antifungal drugs against A. fumigatus. Azoles are also used as agricultural fungicides to combat other fungi that attack plants. Azole-resistant A. fumigatus has been a problem in Europe and Asia for 20 years and has recently been reported in patients in the United States (U.S.). Until this study, we did not know much about azole-resistant A. fumigatus in agricultural settings in the U.S. In this study, we isolated azole-resistant A. fumigatus from multiple states and compared it to isolates from around the world. We show that A. fumigatus which is resistant to azoles and to other strictly agricultural fungicides is widespread in the U.S.

Genomic diversity of the pathogenic fungus Aspergillus fumigatus in Japan reveals the complex genomic basis of azole resistance

Aspergillus fumigatus is a pathogenic fungus with a global distribution. The emergence of azole-resistant A. fumigatus (ARAf) other than the TR-mutants is a problem in Japan. Additionally, the genetic diversity of A. fumigatus strains in Japan remains relatively unknown. Here we show the diversity in the A. fumigatus strains isolated in Japan as well as the complexity in the global distribution of the pathogenic strains. First, we analyzed the genome sequences of 171 strains from Japan as well as the antifungal susceptibility of these strains. Next, we conducted a population analysis of 876 strains by combining the available genomic data for strains isolated worldwide, which were grouped in six clusters. Finally, a genome-wide association study identified the genomic loci associated with ARAf strains, but not the TR-mutants. These results highlight the complexity of the genomic mechanism underlying the emergence of ARAf strains other than the TR-mutants.

Genome-wide transcriptional analyses of Clarireedia jacksonii isolates associated with multi-drug resistance

Emerging multidrug resistance (MDR) in Clarireedia spp. is a huge challenge to the management of dollar spot (DS) disease on turfgrass. Insight into the molecular basis of resistance mechanisms may help identify key molecular targets for developing novel effective chemicals. Previously, a MDR isolate (LT586) of C. jacksonii with significantly reduced sensitivities to propiconazole, boscalid, and iprodione, and a fungicide-sensitive isolate (LT15) of the same species were isolated from creeping bentgrass (Agrostis stolonifera L.). The present study aimed to further explore the molecular mechanisms of resistance by using genome-wide transcriptional analyses of the two isolates. A total of 619 and 475 differentially expressed genes (DEGs) were significantly down and upregulated in the MDR isolate LT586, compared with the sensitive isolate LT15 without fungicide treatment. Three hundreds and six and 153 DEGs showed significantly lower and higher expression in the MDR isolate LT586 than those in the sensitive isolate LT15, which were commonly induced by the three fungicides. Most of the 153 upregulated DEGs were xenobiotic detoxification-related genes and genes with transcriptional functions. Fifty and 17 upregulated DEGs were also commonly observed in HRI11 (a MDR isolate of the C. jacksonii) compared with the HRS10 (a fungicide-sensitive isolate of same species) from a previous study without and with the treatment of propiconazole, respectively. The reliability of RNA-seq data was further verified by qRT-PCR method using a few select potentially MDR-related genes. Results of this study indicated that there were multiple uncharacterized genes, possibly responsible for MDR phenotypes in Clarireedia spp., which may have important implications in understanding the molecular mechanisms underlying MDR resistance.

The human fungal pathogen Aspergillus fumigatus can produce the highest known number of meiotic crossovers

Sexual reproduction involving meiosis is essential in most eukaryotes. This produces offspring with novel genotypes, both by segregation of parental chromosomes as well as crossovers between homologous chromosomes. A sexual cycle for the opportunistic human pathogenic fungus Aspergillus fumigatus is known, but the genetic consequences of meiosis have remained unknown. Among other Aspergilli, it is known that A. flavus has a moderately high recombination rate with an average of 4.2 crossovers per chromosome pair, whereas A. nidulans has in contrast a higher rate with 9.3 crossovers per chromosome pair. Here, we show in a cross between A. fumigatus strains that they produce an average of 29.9 crossovers per chromosome pair and large variation in total map length across additional strain crosses. This rate of crossovers per chromosome is more than twice that seen for any known organism, which we discuss in relation to other genetic model systems. We validate this high rate of crossovers through mapping of resistance to the laboratory antifungal acriflavine by using standing variation in an undescribed ABC efflux transporter. We then demonstrate that this rate of crossovers is sufficient to produce one of the common multidrug resistant haplotypes found in the cyp51A gene (TR34/L98H) in crosses among parents harboring either of 2 nearby genetic variants, possibly explaining the early spread of such haplotypes. Our results suggest that genomic studies in this species should reassess common assumptions about linkage between genetic regions. The finding of an unparalleled crossover rate in A. fumigatus provides opportunities to understand why these rates are not generally higher in other eukaryotes.

Identifying genetic variants associated with amphotericin B (AMB) resistance in Aspergillus fumigatus via k-mer-based GWAS

Aspergillus fumigatus is one of the most common pathogenic fungi, which results in high morbidity and mortality in immunocompromised patients. Amphotericin B (AMB) is used as the core drug for the treatment of triazole-resistant A. fumigatus. Following the usage of amphotericin B drugs, the number of amphotericin B-resistant A. fumigatus isolates showed an increasing trend over the years, but the mechanism and mutations associated with amphotericin B sensitivity are not fully understood. In this study, we performed a k-mer-based genome-wide association study (GWAS) in 98 A. fumigatus isolates from public databases. Associations identified with k-mers not only recapitulate those with SNPs but also discover new associations with insertion/deletion (indel). Compared to SNP sites, the indel showed a stronger association with amphotericin B resistance, and a significant correlated indel is present in the exon region of AFUA_7G05160, encoding a fumarylacetoacetate hydrolase (FAH) family protein. Enrichment analysis revealed sphingolipid synthesis and transmembrane transport may be related to the resistance of A. fumigatus to amphotericin B. The expansion of variant types detected by the k-mer method increases opportunities to identify and exploit complex genetic variants that drive amphotericin B resistance, and these candidate variants help accelerate the selection of prospective gene markers for amphotericin B resistance screening in A. fumigatus.

Pulmonary Aspergillosis in Humboldt Penguins—Susceptibility Patterns and Molecular Epidemiology of Clinical and Environmental Aspergillus fumigatus Isolates from a Belgian Zoo, 2017–2022

Aspergillus fumigatus is the main causative agent of avian aspergillosis and results in significant health problems in birds, especially those living in captivity. The fungal contamination by A. fumigatus in the environment of Humboldt penguins (Spheniscus humboldti), located in a Belgian zoo, was assessed through the analysis of air, water, sand and nest samples during four non-consecutive days in 2021–2022. From these samples, potential azole-resistant A. fumigatus (ARAF) isolates were detected using a selective culture medium. A total of 28 veterinary isolates obtained after necropsy of Humboldt penguins and other avian species from the zoo were also included. All veterinary and suspected ARAF isolates from the environment were characterized for their azole-resistance profile by broth microdilution. Isolates displaying phenotypic resistance against at least one medical azole were systematically screened for mutations in the cyp51A gene. A total of 14 (13.6%) ARAF isolates were identified from the environment (n = 8) and from Humboldt penguins (n = 6). The TR34/L98H mutation was observed in all resistant environmental strains, and in two resistant veterinary strains. To the best of our knowledge, this is the first description of this mutation in A. fumigatus isolates from Humboldt penguins. During the period 2017–2022, pulmonary aspergillosis was confirmed in 51 necropsied penguins, which reflects a death rate due to aspergillosis of 68.0%, mostly affecting adults. Microsatellite polymorphism analysis revealed a high level of diversity among environmental and veterinary A. fumigatus isolates. However, a cluster was observed between one veterinary isolate and six environmental strains, all resistant to medical azoles. In conclusion, the environment of the Humboldt penguins is a potential contamination source of ARAF, making their management even more complex.

Heritability and gene functions associated with sclerotia formation of Rhizoctonia solani AG-7 using whole genome sequencing and genome-wide association study

Sclerotia are specialized fungal structures formed by pigmented and aggregated hyphae, which can survive under unfavourable environmental conditions and serve as the primary inocula for several phytopathogenic fungi including Rhizoctonia solani. Among 154 R. solani anastomosis group 7 (AG-7) isolates collected in fields, the sclerotia-forming capability regarding sclerotia number and sclerotia size varied in the fungal population, but the genetic makeup of these phenotypes remained unclear. As limited studies have focused on the genomics of R. solani AG-7 and the population genetics of sclerotia formation, this study completed the whole genome sequencing and gene prediction of R. solani AG-7 using the Oxford NanoPore and Illumina RNA sequencing. Meanwhile, a high-throughput image-based method was established to quantify the sclerotia-forming capability, and the phenotypic correlation between sclerotia number and sclerotia size was low. A genome-wide association study identified three and five significant SNPs associated with sclerotia number and size in distinct genomic regions, respectively. Of these significant SNPs, two and four showed significant differences in the phenotypic mean separation for sclerotia number and sclerotia size, respectively. Gene ontology enrichment analysis focusing on the linkage disequilibrium blocks of significant SNPs identified more categories related to oxidative stress for sclerotia number, and more categories related to cell development, signalling and metabolism for sclerotia size. These results indicated that different genetic mechanisms may underlie these two phenotypes. Moreover, the heritability of sclerotia number and sclerotia size were estimated for the first time to be 0.92 and 0.31, respectively. This study provides new insights into the heritability and gene functions related to the development of sclerotia number and sclerotia size, which could provide additional knowledge to reduce fungal residues in fields and achieve sustainable disease management.

Genomic and Molecular Identification of Genes Contributing to the Caspofungin Paradoxical Effect in Aspergillus fumigatus

Aspergillus fumigatus is a deadly opportunistic fungal pathogen responsible for ~100,000 annual deaths. Azoles are the first line antifungal agent used against A. fumigatus, but azole resistance has rapidly evolved making treatment challenging. Caspofungin is an important second-line therapy against invasive pulmonary aspergillosis, a severe A. fumigatus infection. Caspofungin functions by inhibiting β-1,3-glucan synthesis, a primary and essential component of the fungal cell wall. A phenomenon termed the caspofungin paradoxical effect (CPE) has been observed in several fungal species where at higher concentrations of caspofungin, chitin replaces β-1,3-glucan, morphology returns to normal, and growth rate increases. CPE appears to occur in vivo, and it is therefore clinically important to better understand the genetic contributors to CPE. We applied genomewide association (GWA) analysis and molecular genetics to identify and validate candidate genes involved in CPE. We quantified CPE across 67 clinical isolates and conducted three independent GWA analyses to identify genetic variants associated with CPE. We identified 48 single nucleotide polymorphisms (SNPs) associated with CPE. We used a CRISPR/Cas9 approach to generate gene deletion mutants for seven genes harboring candidate SNPs. Two null mutants, ΔAfu3g13230 and ΔAfu4g07080 (dscP), resulted in reduced basal growth rate and a loss of CPE. We further characterized the dscP phosphatase-null mutant and observed a significant reduction in conidia production and extremely high sensitivity to caspofungin at both low and high concentrations. Collectively, our work reveals the contribution of Afu3g13230 and dscP in CPE and sheds new light on the complex genetic interactions governing this phenotype. IMPORTANCE This is one of the first studies to apply genomewide association (GWA) analysis to identify genes involved in an Aspergillus fumigatus phenotype. A. fumigatus is an opportunistic fungal pathogen that causes hundreds of thousands of infections and ~100,000 deaths each year, and antifungal resistance has rapidly evolved in this species. A phenomenon called the caspofungin paradoxical effect (CPE) occurs in some isolates, where high concentrations of the drug lead to increased growth rate. There is clinical relevance in understanding the genetic basis of this phenotype, since caspofungin concentrations could lead to unintended adverse clinical outcomes in certain cases. Using GWA analysis, we identified several interesting candidate polymorphisms and genes and then generated gene deletion mutants to determine whether these genes were important for CPE. Two of these mutant strains (ΔAfu3g13230 and ΔAfu4g07080/ΔdscP) displayed a loss of the CPE. This study sheds light on the genes involved in clinically important phenotype CPE.

Using genomics to understand the mechanisms of virulence and drug resistance in fungal pathogens

Fungal pathogens pose an increasingly worrying threat to human health, food security and ecosystem diversity. To tackle fungal infections and improve current diagnostic and therapeutic tools it is necessary to understand virulence and antifungal drug resistance mechanisms in diverse species. Recent advances in genomics approaches have provided a suitable framework to understand these phenotypes, which ultimately depend on genetically encoded determinants. In this work, we review how the study of genome sequences has been key to ascertain the bases of virulence and drug resistance traits. We focus on the contribution of comparative genomics, population genomics and directed evolution studies. In addition, we discuss how different types of genomic mutations (small or structural variants) contribute to intraspecific differences in virulence or drug resistance. Finally, we review current challenges in the field and anticipate future directions to solve them. In summary, this work provides a short overview of how genomics can be used to understand virulence and drug resistance in fungal pathogens.

Genetic Structure and Triazole Antifungal Susceptibilities of Alternaria alternata from Greenhouses in Kunming, China

Alternaria alternata is an opportunistic human fungal pathogen and a ubiquitous phytopathogen capable of causing diseases to >100 agricultural crops and ornamental plants. To control plant diseases caused by A. alternata, triazole fungicides have been widely used both in open crop and vegetable fields and in indoor growth facilities such as greenhouses. At present, the effect of fungicide use on triazole resistance development in A. alternata populations is not known. Here, we isolated 237 A. alternata strains from nine greenhouses around metropolitan Kunming in Yunnan, southwest China, determined their genotypes using 10 short tandem repeat markers, and quantified their susceptibility to four triazoles (difenoconazole, tebuconazole, itraconazole, and voriconazole). Abundant allelic and genotypic diversities were detected among these A. alternata strains. Significantly, over 17% of the strains were resistant to difenoconazole, and both known and new drug-resistance mutations were found in the triazole target gene cyp51. Our findings of high-level genetic variation of A. alternata in greenhouses coupled with high-frequency fungicide resistance call for greater attention to continued monitoring and to developing alternative plant fungal disease management strategies in greenhouses. IMPORTANCE Alternaria alternata is among the most common fungi in our environments, such as indoor facilities, the soil, and outdoor air. It can cause diseases in >100 crop and ornamental plants. Furthermore, it can cause human infections. However, our understanding of its genetic diversity and antifungal susceptibility is very limited. Indeed, the critical threshold values for resistance have not been defined for most antifungal drugs in this species. Greenhouses are known to have heavy applications of agricultural fungicides. In this study, we analyzed strains of A. alternata from nine greenhouses near metropolitan Kunming in southwestern China. Our study revealed very high genetic diversity and identified strains with high MIC values against two agricultural and two medical triazole antifungals within each of the nine greenhouses. Our study calls for greater attention to this emerging threat to food security and human health.

Genetic Diversity and Azole Resistance Among Natural Aspergillus fumigatus Populations in Yunnan, China

The emergence and spread of azole resistance alleles in clinical and environmental isolates of Aspergillus fumigatus is a global human health concern and endangers the "One Health" approach in our fight against antifungal resistance (AFR) in this pathogen. A major challenge to combat AFR in A. fumigatus is the massive aerial dispersal ability of its asexual spores. Our recent fine-scale survey of greenhouse populations of A. fumigatus near Kunming, Yunnan, China, suggested that the use of azole fungicides for plant protection was likely a major driver of the high-frequency azole-resistant A. fumigatus (ARAF) in greenhouses. Here, we investigated the potential spread of those ARAF and the structure of geographic populations of A. fumigatus by analyzing 452 isolates from 19 geographic locations across Yunnan. We found lower frequencies of ARAF in these outdoor populations than those in greenhouses near Kunming, but there were abundant new alleles and new genotypes, including those associated with azole resistance, consistent with multiple independent origins of ARAF across Yunnan. Interestingly, among the four ecological niches, the sediments of a large lake near Kunming were found to have the highest frequency of ARAF (~ 43%). While most genetic variations were observed within the 19 local populations, statistically significant genetic differentiations were found between many subpopulations within Yunnan. Furthermore, similar to greenhouse populations, these outdoor populations of A. fumigatus in Yunnan were significantly different from those in other parts of the world. Our results call for increased attention to local and regional studies of this fungal pathogen to help develop targeted control strategies against ARAF.

Use of Bulk Segregant Analysis for Determining the Genetic Basis of Azole Resistance in the Opportunistic Pathogen Aspergillus fumigatus

A sexual cycle was described in 2009 for the opportunistic fungal pathogen Aspergillus fumigatus, opening up for the first time the possibility of using techniques reliant on sexual crossing for genetic analysis. The present study was undertaken to evaluate whether the technique 'bulk segregant analysis' (BSA), which involves detection of differences between pools of progeny varying in a particular trait, could be applied in conjunction with next-generation sequencing to investigate the underlying basis of monogenic traits in A. fumigatus. Resistance to the azole antifungal itraconazole was chosen as a model, with a dedicated bioinformatic pipeline developed to allow identification of SNPs that differed between the resistant progeny pool and resistant parent compared to the sensitive progeny pool and parent. A clinical isolate exhibiting monogenic resistance to itraconazole of unknown basis was crossed to a sensitive parent and F1 progeny used in BSA. In addition, the use of backcrossing and increasing the number in progeny pools was evaluated as ways to enhance the efficiency of BSA. Use of F1 pools of 40 progeny led to the identification of 123 candidate genes with SNPs distributed over several contigs when aligned to an A1163 reference genome. Successive rounds of backcrossing enhanced the ability to identify specific genes and a genomic region, with BSA of progeny (using 40 per pool) from a third backcross identifying 46 genes with SNPs, and BSA of progeny from a sixth backcross identifying 20 genes with SNPs in a single 292 kb region of the genome. The use of an increased number of 80 progeny per pool also increased the resolution of BSA, with 29 genes demonstrating SNPs between the different sensitive and resistant groupings detected using progeny from just the second backcross with the majority of variants located on the same 292 kb region. Further bioinformatic analysis of the 292 kb region identified the presence of a cyp51A gene variant resulting in a methionine to lysine (M220K) change in the CYP51A protein, which was concluded to be the causal basis of the observed resistance to itraconazole. The future use of BSA in genetic analysis of A. fumigatus is discussed.

Epidemiological Studies of Pan-Azole Resistant Aspergillus fumigatus Populations Sampled during Tulip Cultivation Show Clonal Expansion with Acquisition of Multi-Fungicide Resistance as Potential Driver

Pan-azole resistant isolates are found in clinical and environmental Aspergillus fumigatus (Af) populations. Azole resistance can evolve in both settings, with Af directly targeted by antifungals in patients and, in the environment, Af unintendedly exposed to fungicides used for material preservation and plant disease control. Resistance to non-azole fungicides, including methyl benzimidazole carbamates (MBCs), quinone outside inhibitors (QoIs) and succinate dehydrogenase inhibitors (SDHIs), has recently been reported. These fungicide groups are not used in medicine but can play an important role in the further spread of pan-azole resistant genotypes. We investigated the multi-fungicide resistance status and the genetic diversity of Af populations sampled from tulip field soils, tulip peel waste and flower compost heaps using fungicide sensitivity testing and a range of genotyping tools, including STRAf typing and sequencing of fungicide resistant alleles. Two major clones were present in the tulip bulb population. Comparisons with clinical isolates and literature data revealed that several common clonal lineages of TR34/L98H and TR46/Y121F/T289A strains that have expanded successfully in the environment have also acquired resistance to MBC, QoI and/or SDHI fungicides. Strains carrying multiple fungicide resistant alleles have a competitive advantage in environments where residues of multiple fungicides belonging to different modes of action are present.

Genetic Analyses of Amphotericin B Susceptibility in Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous saprophytic mold that can cause a range of clinical syndromes, from allergic reactions to invasive infections. Amphotericin B (AMB) is a polyene antifungal drug that has been used to treat a broad range of systemic mycoses since 1958, including as a primary treatment option against invasive aspergillosis in regions with high rates (≥10%) of environmental triazole resistance. However, cases of AMB-resistant A. fumigatus strains have been increasingly documented over the years, and high resistance rates were recently reported in Brazil and Canada. The objective of this study is to identify candidate mutations associated with AMB susceptibility using a genome-wide association analysis of natural strains, and to further investigate a subset of the mutations in their putative associations with differences in AMB minimum inhibitory concentration (MIC) and in growths at different AMB concentrations through the analysis of progeny from a laboratory genetic cross. Together, our results identified a total of 34 candidate single-nucleotide polymorphisms (SNPs) associated with AMB MIC differences-comprising 18 intergenic variants, 14 missense variants, one synonymous variant, and one non-coding transcript variant. Importantly, progeny from the genetic cross allowed us to identify putative SNP-SNP interactions impacting progeny growth at different AMB concentrations.

Intimate genetic relationships and fungicide resistance in multiple strains of Aspergillus fumigatus isolated from a plant bulb

Fungal infections are increasingly dangerous because of environmentally dispersed resistance to antifungal drugs. Azoles are commonly used antifungal drugs, but they are also used as fungicides in agriculture, which may enable enrichment of azole-resistant strains of the human pathogen Aspergillus fumigatus in the environment. Understanding of environmental dissemination and enrichment of genetic variation associated with azole resistance in A. fumigatus is required to suppress resistant strains. Here, we focused on eight strains of azole-resistant A. fumigatus isolated from a single tulip bulb for sale in Japan. This set includes strains with TR34 /L98H/T289A/I364V/G448S and TR46 /Y121F/T289A/S363P/I364V/G448S mutations in the cyp51A gene, which showed higher tolerance to several azoles than strains harbouring TR46 /Y121F/T289A mutation. The strains were typed by microsatellite typing, single nucleotide polymorphism profiles, and mitochondrial and nuclear genome analyses. The strains grouped differently using each typing method, suggesting historical genetic recombination among the strains. Our data also revealed that some strains isolated from the tulip bulb showed tolerance to other classes of fungicides, such as QoI and carbendazim, followed by related amino acid alterations in the target proteins. Considering spatial-temporal factors, plant bulbs are an excellent environmental niche for fungal strains to encounter partners, and to obtain and spread resistance-associated mutations.