Multilocus resistance evolution to azole fungicides in fungal plant pathogen populations

Evolution of fungicide resistance is a major threat to food production in agricultural ecosystems. Fungal pathogens rapidly evolved resistance to all classes of fungicides applied to the field. Resistance to the commonly used azole fungicides is thought to be driven mainly by mutations in a gene (CYP51) encoding a protein of the ergosterol biosynthesis pathway. However, some fungi gained azole resistance independently of CYP51 mutations and the mechanisms leading to CYP51-independent resistance are poorly understood. We used whole-genome sequencing and genome-wide association studies (GWAS) to perform an unbiased screen of azole resistance loci in Rhynchosporium commune, the causal agent of the barley scald disease. We assayed cyproconazole resistance in 120 isolates collected from nine populations worldwide. We found that mutations in highly conserved genes encoding the vacuolar cation channel YVC1, a transcription activator, and a saccharopine dehydrogenase made significant contributions to fungicide resistance. These three genes were not previously known to confer resistance in plant pathogens. However, YVC1 is involved in a conserved stress response pathway known to respond to azoles in human pathogenic fungi. We also performed GWAS to identify genetic polymorphism linked to fungal growth rates. We found that loci conferring increased fungicide resistance were negatively impacting growth rates, suggesting that fungicide resistance evolution imposed costs. Analyses of population structure showed that resistance mutations were likely introduced into local populations through gene flow. Multilocus resistance evolution to fungicides shows how pathogen populations can evolve a complex genetic architecture for an important phenotypic trait within a short time span.