Genome variation of Sporothrix schenckii and Sporothrix brasiliensis

Sporotrichosis, a subcutaneous mycosis caused by dimorphic fungi of the Sporothrix genus, has become a major zoonotic epidemic in South America, primarily driven by Sporothrix brasiliensis. To elucidate the genomic basis of its emergence and antifungal adaptation, we analyzed whole-genome sequences from 95 Sporothrix isolates, integrating single-nucleotide polymorphism (SNP), copy number variation (CNV), and genome-wide association (GWAS) analyses. Comparative genomics revealed 610,242 SNPs within S. brasiliensis and 1,474,627 within S. schenckii, confirming a marked disparity in intraspecific diversity. Phylogenomic tree inference resolved six well-supported S. brasiliensis clades with limited internal divergence, reflecting recent population expansion, while S. schenckii displayed deep phylogeographic structure separating North and South American lineages. CNV profiling identified 158 affected genes in S. brasiliensis (60 gains, 98 losses) and 88 in S. schenckii (54 gains, 34 losses), concentrated near sub-telomeric regions. In S. brasiliensis, gains were enriched for kinases and intracellular trafficking functions, whereas losses involved genes related to translation and primary metabolism, suggesting regulatory reinforcement coupled with metabolic streamlining. A GWAS of itraconazole resistance identified 81 SNPs distributed across multiple scaffolds, with many located within genes related with transport, signaling, and redox balance, supporting a polygenic basis for azole response. The implication of these alleles with itraconazole resistance warrants experimental followup. Together, these results highlight distinct evolutionary strategies of closely related Sporothrix species and pinpoint genomic changes potentially associated with the emergence and drug tolerance of S. brasiliensis.