Parallel evolution under constraint shapes echinocandin resistance in Candida auris

Drug resistance emerges repeatedly in outbreaks of Candida fungal pathogens, but little is known about its origins or persistence. Here, we investigated the evolutionary processes shaping echinocandin resistance in Candida auris, a globally emerging and predominantly clonal fungal pathogen. Genome-wide association across over 600 isolates identified mutations in the {beta}-1,3-glucan synthase gene FKS1 as the most significant driver of resistance to an echinocandin drug. Ancestral reconstruction of this population traced shared resistance mutations among small groups typically consisting of 2-3 closely related isolates, but clusters could include up to 16 isolates. Nearly all resistant clusters consisted of isolates collected in the same year and region, consistent with local transmission. To further examine population-level selection, we measured adaptive signatures in FKS1 and the highly diverged paralog FKS2 across 22,000 genomes. This revealed excess nonsynonymous polymorphisms in FKS1, primarily due to independent, recurrent mutations at resistance hotspots, consistent with parallel evolution and incomplete fixation of adaptive alleles. In FKS2, there is no evidence of hotspots and little support for diversifying selection. Together, these results indicate that resistance mutations emerge under strong genetic constraint, with adaptation restricted to only one FKS homolog and predominantly at mutational hotspots.