Nonrandom Distribution of Azole Resistance across the Global Population of Aspergillus fumigatus

The emergence of azole resistance in the pathogenic fungus Aspergillus fumigatus has continued to increase, with the dominant resistance mechanisms, consisting of a 34-nucleotide tandem repeat (TR(34))/L98H and TR(46)/Y121F/T289A, now showing a structured global distribution. Using hierarchical clustering and multivariate analysis of 4,049 A. fumigatus isolates collected worldwide and genotyped at nine microsatellite loci using analysis of short tandem repeats of A. fumigatus (STRAf), we show that A. fumigatus can be subdivided into two broad clades and that cyp51A alleles TR(34)/L98H and TR(46)/Y121F/T289A are unevenly distributed across these two populations. Diversity indices show that azole-resistant isolates are genetically depauperate compared to their wild-type counterparts, compatible with selective sweeps accompanying the selection of beneficial mutations. Strikingly, we found that azole-resistant clones with identical microsatellite profiles were globally distributed and sourced from both clinical and environmental locations, confirming that azole resistance is an international public health concern. Our work provides a framework for the analysis of A. fumigatus isolates based on their microsatellite profile, which we have incorporated into a freely available, user-friendly R Shiny application (AfumID) that provides clinicians and researchers with a method for the fast, automated characterization of A. fumigatus genetic relatedness. Our study highlights the effect that azole drug resistance is having on the genetic diversity of A. fumigatus and emphasizes its global importance upon this medically important pathogenic fungus.