Genomic Assembly of Clinical Candida glabrata (Nakaseomyces glabrata) Isolates Reveals within-Species Structural Plasticity and Association with In Vitro Antifungal Susceptibility

Abstract

The human pathogen Candida glabrata has become a global threat to human health, with infections characterized by high mortality and multidrug resistance. We have obtained nine fully assembled genomes from clinical isolates through a combination of short- and long-read sequencing approaches. The opportunistic human pathogen Candida glabrata has become an increasingly important threat to human health, with infections globally characterized by high mortality rates and multidrug resistance. To face this threat, more efficient diagnostic and therapeutic approaches are required, underpinning research to help define the intraspecies epidemiology, genetic variability, and therefore, diagnostic and therapeutic target stability. Previous comparative genetics studies conducted on limited numbers of strains only revealed partial resolution of chromosomal settings. In this study, by combining short- and long-read genome sequencing, phenotypic characterization, and comparative genomics over a large set of strains, we detected strict relationships between large chromosomal rearrangements and phylogenetic clades, genes subjected to different selective pressures, and new sets of genes associated with resistance to antifungals. Overall, these results not only provide a fundamental contribution to our knowledge of C. glabrata evolution and epidemiology but may also lay the foundations for the future development of tailored therapeutic approaches. IMPORTANCE The human pathogen Candida glabrata has become a global threat to human health, with infections characterized by high mortality and multidrug resistance. We have obtained nine fully assembled genomes from clinical isolates through a combination of short- and long-read sequencing approaches. The quality and completeness of such genomes and their subsequent comparison to the broadest set of genomes so far allowed us to pinpoint chromosomal rearrangements in several genomes and detect phylogenetic clades that were not associated with geographic location or isolation source. We identified a new set of genes associated with resistance to antifungals coding for adhesin or adhesin-like proteins, suggesting C. glabrata resists antifungals by forming aggregates or adhering to the host tissue. These results, which provide a fundamental contribution to our knowledge of C. glabrata evolution and epidemiology, may initiate the development of precision medicine interventions for patients with suspected or proven invasive fungal infections.