Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes

We used whole-genome sequencing to analyze a collection of 35 fluconazole-resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two persistent clinical lineages were identified. We identified copy number variation (CNV) of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show that the annotated telomeric gene CDR1B is actually an artifactual in silico fusion of two highly similar neighboring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased the expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

A secondary mechanism of action for triazole antifungals in Aspergillus fumigatus mediated by hmg1

Triazole antifungals function as ergosterol biosynthesis inhibitors and are frontline therapy for invasive fungal infections, such as invasive aspergillosis. The primary mechanism of action of triazoles is through the specific inhibition of a cytochrome P450 14-α-sterol demethylase enzyme, Cyp51A/B, resulting in depletion of cellular ergosterol. Here, we uncover a clinically relevant secondary mechanism of action for triazoles within the ergosterol biosynthesis pathway. We provide evidence that triazole-mediated inhibition of Cyp51A/B activity generates sterol intermediate perturbations that are likely decoded by the sterol sensing functions of HMG-CoA reductase and Insulin-Induced Gene orthologs as increased pathway activity. This, in turn, results in negative feedback regulation of HMG-CoA reductase, the rate-limiting step of sterol biosynthesis. We also provide evidence that HMG-CoA reductase sterol sensing domain mutations previously identified as generating resistance in clinical isolates of Aspergillus fumigatus partially disrupt this triazole-induced feedback. Therefore, our data point to a secondary mechanism of action for the triazoles: induction of HMG-CoA reductase negative feedback for downregulation of ergosterol biosynthesis pathway activity. Abrogation of this feedback through acquired mutations in the HMG-CoA reductase sterol sensing domain diminishes triazole antifungal activity against fungal pathogens and underpins HMG-CoA reductase-mediated resistance.

Comparative genomics of Cryptococcus and Kwoniella reveals pathogenesis evolution and contrasting karyotype dynamics via intercentromeric recombination or chromosome fusion

A large-scale comparative genomic analysis was conducted for the global human fungal pathogens within the Cryptococcus genus, compared to non-pathogenic Cryptococcus species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species of both genera, resulting in a dataset encompassing virtually all of their known diversity. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at pre-adaptive pathogenic potential, our analysis found evidence in pathogenic Cryptococcus species of specific examples of gene gain (via horizontal gene transfer) and gene loss, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the two genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5 or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes, underwent chromosome reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Taken together, our findings advance our understanding of genomic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Analysis of clinical Candida parapsilosis isolates reveals copy number variation in key fluconazole resistance genes

We used whole-genome sequencing to analyse a collection of 35 fluconazole resistant and 7 susceptible Candida parapsilosis isolates together with coverage analysis and GWAS techniques to identify new mechanisms of fluconazole resistance. Phylogenetic analysis shows that although the collection is diverse, two probable outbreak groups were identified. We identified copy number variation of two genes, ERG11 and CDR1B, in resistant isolates. Two strains have a CNV at the ERG11 locus; the entire ORF is amplified in one, and only the promoter region is amplified in the other. We show the annotated telomeric gene CDR1B is actually an artefactual in silico fusion of two highly similar neighbouring CDR genes due to an assembly error in the C. parapsilosis CDC317 reference genome. We report highly variable copy numbers of the CDR1B region across the collection. Several strains have increased expansion of the two genes into a tandem array of new chimeric genes. Other strains have experienced a deletion between the two genes creating a single gene with a reciprocal chimerism. We find translocations, duplications, and gene conversion across the CDR gene family in the C. parapsilosis species complex, showing that it is a highly dynamic family.

Genetic diversity and microevolution in clinical Cryptococcus isolates from Cameroon

Cryptococcal meningitis is the second most common cause of death in people living with HIV/AIDS, yet we have a limited understanding of how cryptococcal isolates change over the course of infection. Cryptococcal infections are environmentally acquired, and the genetic diversity of these infecting isolates can also be geographically linked. Here, we employ whole genome sequences for 372 clinical Cryptococcus isolates from 341 patients with HIV-associated cryptococcal meningitis obtained via a large clinical trial, across both Malawi and Cameroon, to enable population genetic comparisons of isolates between countries. We see that isolates from Cameroon are highly clonal, when compared to those from Malawi, with differential rates of disruptive variants in genes with roles in DNA binding and energy use. For a subset of patients (22) from Cameroon, we leverage longitudinal sampling, with samples taken at days 7 and 14 post-enrollment, to interrogate the genetic changes that arise over the course of infection, and the genetic diversity of isolates within patients. We see disruptive variants arising over the course of infection in several genes, including the phagocytosis-regulating transcription factor GAT204. In addition, in 13% of patients sampled longitudinally, we see evidence for mixed infections. This approach identifies geographically linked genetic variation, signatures of microevolution, and evidence for mixed infections across a clinical cohort of patients affected by cryptococcal meningitis in Central Africa.

Impact of pathogen genetics on clinical phenotypes in a population of Talaromyces marneffei from Vietnam

Talaromycosis, a severe and invasive fungal infection caused by Talaromyces marneffei, is difficult to treat and impacts those living in endemic regions of Southeast Asia, India, and China. While 30% of infections result in mortality, our understanding of the genetic basis of pathogenesis for this fungus is limited. To address this, we apply population genomics and genome-wide association study approaches to a cohort of 336 T. marneffei isolates collected from patients who enrolled in the Itraconazole vs Amphotericin B for Talaromycosis trial in Vietnam. We find that isolates from northern and southern Vietnam form two distinct geographical clades, with isolates from southern Vietnam associated with increased disease severity. Leveraging longitudinal isolates, we identify multiple instances of disease relapse linked to unrelated strains, highlighting the potential for multistrain infections. In more frequent cases of persistent talaromycosis caused by the same strain, we identify variants arising over the course of patient infections that impact genes predicted to function in the regulation of gene expression and secondary metabolite production. By combining genetic variant data with patient metadata for all 336 isolates, we identify pathogen variants significantly associated with multiple clinical phenotypes. In addition, we identify genes and genomic regions under selection across both clades, highlighting loci undergoing rapid evolution, potentially in response to external pressures. With this combination of approaches, we identify links between pathogen genetics and patient outcomes and identify genomic regions that are altered during T. marneffei infection, providing an initial view of how pathogen genetics affects disease outcomes.

Genome Sequence of Lichtheimia ornata, an Emerging Opportunistic Mucorales Pathogen

Lichtheimia ornata is an emerging opportunistic Mucorales pathogen that is associated with fatal infections in immunocompromised individuals. While these environmentally acquired infections have rarely been reported to date, cases were noted in a recent analysis of coronavirus disease 2019 (COVID-19)-associated mucormycosis in India. Here, we report the annotated genome sequence of the environmental isolate CBS 291.66.

Impact of pathogen genetics on clinical phenotypes in a population of Talaromyces marneffei from Vietnam

Talaromycosis, a severe and invasive fungal infection caused by Talaromyces marneffei , is difficult to treat and impacts those living in endemic regions of southeast Asia, India, and China. While 30% of infections result in mortality, our understanding of the genetic basis of pathogenesis for this fungus is limited. To address this, we apply population genomics and genome wide association study approaches to a cohort of 336 T. marneffei isolates collected from patients who enrolled in the Itraconazole versus Amphotericin B for Talaromycosis (IVAP) trial in Vietnam. We find that isolates from northern and southern Vietnam form two distinct geographical clades, with isolates from southern Vietnam associated with increased disease severity. Leveraging longitudinal isolates, we identify multiple instances of disease relapse linked to unrelated strains, highlighting the potential for multi-strain infections. In more frequent cases of persistent talaromycosis caused by the same strain, we identify variants arising over the course of patient infections that impact genes predicted to function in the regulation of gene expression and secondary metabolite production. By combining genetic variant data with patient metadata for all 336 isolates, we identify pathogen variants significantly associated with multiple clinical phenotypes. In addition, we identify genes and genomic regions under selection across both clades, highlighting loci undergoing rapid evolution, potentially in response to external pressures. With this combination of approaches, we identify links between pathogen genetics and patient outcomes and identify genomic regions that are altered during T. marneffei infection, providing an initial view of how pathogen genetics affects disease outcomes.

Comparing genomic variant identification protocols for Candida auris

Genomic analyses are widely applied to epidemiological, population genetic and experimental studies of pathogenic fungi. A wide range of methods are employed to carry out these analyses, typically without including controls that gauge the accuracy of variant prediction. The importance of tracking outbreaks at a global scale has raised the urgency of establishing high-accuracy pipelines that generate consistent results between research groups. To evaluate currently employed methods for whole-genome variant detection and elaborate best practices for fungal pathogens, we compared how 14 independent variant calling pipelines performed across 35 Candida auris isolates from 4 distinct clades and evaluated the performance of variant calling, single-nucleotide polymorphism (SNP) counts and phylogenetic inference results. Although these pipelines used different variant callers and filtering criteria, we found high overall agreement of SNPs from each pipeline. This concordance correlated with site quality, as SNPs discovered by a few pipelines tended to show lower mapping quality scores and depth of coverage than those recovered by all pipelines. We observed that the major differences between pipelines were due to variation in read trimming strategies, SNP calling methods and parameters, and downstream filtration criteria. We calculated specificity and sensitivity for each pipeline by aligning three isolates with chromosomal level assemblies and found that the GATK-based pipelines were well balanced between these metrics. Selection of trimming methods had a greater impact on SAMtools-based pipelines than those using GATK. Phylogenetic trees inferred by each pipeline showed high consistency at the clade level, but there was more variability between isolates from a single outbreak, with pipelines that used more stringent cutoffs having lower resolution. This project generated two truth datasets useful for routine benchmarking of C. auris variant calling, a consensus VCF of genotypes discovered by 10 or more pipelines across these 35 diverse isolates and variants for 2 samples identified from whole-genome alignments. This study provides a foundation for evaluating SNP calling pipelines and developing best practices for future fungal genomic studies.

Comparison of long- and short-read metagenomic assembly for low-abundance species and resistance genes

Recent technological and computational advances have made metagenomic assembly a viable approach to achieving high-resolution views of complex microbial communities. In previous benchmarking, short-read (SR) metagenomic assemblers had the highest accuracy, long-read (LR) assemblers generated the most contiguous sequences and hybrid (HY) assemblers balanced length and accuracy. However, no assessments have specifically compared the performance of these assemblers on low-abundance species, which include clinically relevant organisms in the gut. We generated semi-synthetic LR and SR datasets by spiking small and increasing amounts of Escherichia coli isolate reads into fecal metagenomes and, using different assemblers, examined E. coli contigs and the presence of antibiotic resistance genes (ARGs). For ARG assembly, although SR assemblers recovered more ARGs with high accuracy, even at low coverages, LR assemblies allowed for the placement of ARGs within longer, E. coli-specific contigs, thus pinpointing their taxonomic origin. HY assemblies identified resistance genes with high accuracy and had lower contiguity than LR assemblies. Each assembler type's strengths were maintained even when our isolate was spiked in with a competing strain, which fragmented and reduced the accuracy of all assemblies. For strain characterization and determining gene context, LR assembly is optimal, while for base-accurate gene identification, SR assemblers outperform other options. HY assembly offers contiguity and base accuracy, but requires generating data on multiple platforms, and may suffer high misassembly rates when strain diversity exists. Our results highlight the trade-offs associated with each approach for recovering low-abundance taxa, and that the optimal approach is goal-dependent.

Multimodal analysis of the COVID-19-associated mucormycosis outbreak in Delhi, India indicates the convergence of clinical and environmental risk factors

BACKGROUND

The aetiology of the major outbreak of COVID-19-associated mucormycosis (CAM) in India in spring 2021 remains incompletely understood. Herein, we provide a multifaceted and multi-institutional analysis of clinical, pathogen-related, environmental and healthcare-related factors during CAM outbreak in the metropolitan New Delhi area.

METHODS

We reviewed medical records of all patients diagnosed with biopsy-proven CAM (n = 50) at 7 hospitals in the New Delhi, and NCR area in April-June 2021. Two multivariate logistic regression models were used to compare clinical characteristics of CAM cases with COVID-19-hospitalised contemporary patients as controls (n = 69). Additionally, meteorological parameters and mould spore concentrations in outdoor air were analysed. Selected hospital fomites were cultured. Mucorales isolates from CAM patients were analysed by ITS sequencing and whole-genome sequencing (WGS).

RESULTS

Independent risk factors for CAM identified by multivariate analysis were previously or newly diagnosed diabetes mellitus, active cancer and severe COVID-19 infection. Supplemental oxygen, remdesivir therapy and ICU admission for COVID-19 were associated with reduced CAM risk. The CAM incidence peak was preceded by an uptick in environmental spore concentrations in the preceding 3-4 weeks that correlated with increasing temperature, high evaporation and decreasing relative humidity. Rhizopus was the most common genus isolated, but we also identified two cases of the uncommon Mucorales, Lichtheimia ornata. WGS found no clonal population of patient isolates. No Mucorales were cultured from hospital fomites.

CONCLUSIONS

An intersection of host and environmental factors contributed to the emergence of CAM. Surrogates of access to advanced COVID-19 treatment were associated with lower CAM risk.

Genome-Wide Identification of Variants Associated with Antifungal Drug Resistance

Genomic approaches are widely applied to study the genetic basis of antifungal drug resistance in clinical isolates and experimental studies. Whole-genome sequencing of clinical isolates can comprehensively identify mutations associated with drug resistance and their frequency across fungal populations. In addition, genome comparison of serially collected isolates, such as from patient samples or in vitro drug selection experiments, will identify a small number of changes that can be evaluated for association with drug resistance. Here, we provide a detailed protocol for the computational analysis of genome sequences to identify variants associated with drug resistance.

Genomic Variation across a Clinical Cryptococcus Population Linked to Disease Outcome

Cryptococcus neoformans is the causative agent of cryptococcosis, a disease with poor patient outcomes that accounts for approximately 180,000 deaths each year. Patient outcomes may be impacted by the underlying genetics of the infecting isolate; however, our current understanding of how genetic diversity contributes to clinical outcomes is limited. Here, we leverage clinical, in vitro growth and genomic data for 284 C. neoformans isolates to identify clinically relevant pathogen variants within a population of clinical isolates from patients with human immunodeficiency virus (HIV)-associated cryptococcosis in Malawi. Through a genome-wide association study (GWAS) approach, we identify variants associated with the fungal burden and the growth rate. We also find both small and large-scale variation, including aneuploidy, associated with alternate growth phenotypes, which may impact the course of infection. Genes impacted by these variants are involved in transcriptional regulation, signal transduction, glycosylation, sugar transport, and glycolysis. We show that growth within the central nervous system (CNS) is reliant upon glycolysis in an animal model and likely impacts patient mortality, as the CNS yeast burden likely modulates patient outcome. Additionally, we find that genes with roles in sugar transport are enriched in regions under selection in specific lineages of this clinical population. Further, we demonstrate that genomic variants in two genes identified by GWAS impact virulence in animal models. Our approach identifies links between the genetic variation in C. neoformans and clinically relevant phenotypes and animal model pathogenesis, thereby shedding light on specific survival mechanisms within the CNS and identifying the pathways involved in yeast persistence. IMPORTANCE Infection outcomes for cryptococcosis, most commonly caused by C. neoformans, are influenced by host immune responses as well as by host and pathogen genetics. Infecting yeast isolates are genetically diverse; however, we lack a deep understanding of how this diversity impacts patient outcomes. To better understand both clinical isolate diversity and how diversity contributes to infection outcomes, we utilize a large collection of clinical C. neoformans samples that were isolated from patients enrolled in a clinical trial across 3 hospitals in Malawi. By combining whole-genome sequence data, clinical data, and in vitro growth data, we utilize genome-wide association approaches to examine the genetic basis of virulence. Genes with significant associations display virulence attributes in both murine and rabbit models, demonstrating that our approach can identify potential links between genetic variants and patho-biologically significant phenotypes.

The molecular and genetic basis of antifungal resistance in the emerging fungal pathogen Candida auris

Fungal infections are responsible for significant morbidity and mortality. Resistance to the limited number of agents in the antifungal armamentarium among pathogenic fungi represents a growing public health threat. Particularly concerning is the emerging fungal pathogen Candida auris that frequently exhibits resistance to the triazole class of antifungals and amphotericin B, and for which isolates resistant to all of the major antifungal classes have been reported. In this brief review, we provide an overview of what is currently known about the molecular and genetic basis for antifungal resistance in this fungal pathogen.

First Patient-to-Patient Intrahospital Transmission of Clade I Candida auris in France Revealed after a Two-Month Incubation Period

Candida auris is a recently described emerging pathogen in hospital settings. Five genetic clades have been delineated, with each clade being isolated from specific geographic regions. We here describe the first transmission between 2 patients (P0 and P1) of a clade I C. auris strain imported into our burn intensive care unit from the Middle East. The strains have been investigated with whole-genome sequencing, which validated the high similarity of the genomes between isolates from P0 and P1. We repeatedly screened the two patients and contact patients (i.e., other patients present in the same hospital ward at the time of the first positive sample from P0 or P1; n = 49; 268 tests) with fungal culture and a C. auris-specific quantitative PCR assay to assess transmission patterns. We observed that P1 developed C. auris colonization between 41 and 61 days after potential exposure to P0 contamination, despite three negative screening tests as recommended by our national authorities. This study illustrates that transmission of C. auris between patients can lead to long-term incubation times before the detection of colonization. The recommended screening strategy may not be optimal and should be improved in the light of our findings. IMPORTANCE While large outbreaks of C. auris in hospital settings have been described, few clear cases of direct transmission have been documented. We here investigated the transmission of C. auris clade I between two patients with a 41- to 61-day delay between exposure and the development of colonization. This may lead to changes in the recommendations concerning treatment of C. auris cases, as an incubation period of this length is one of the first to be reported.

S2.5d Exploring multidrug resistance, fitness compensation, and collateral sensitivity in Candida auris : Fight fire with fire?

S2.5 Rare yeasts, September 21, 2022, 3:00 PM - 4:30 PM

Candida auris (C. auris) is a recently emerged human fungal pathogen of growing concern due to its ability to acquire extensive multidrug resistance (MDR) to all four antifungal drug classes. The unprecedented extent of MDR in C. auris, suggests accelerated resistance evolution, novel mechanisms of resistance, and/or potential fitness compensation. Despite being the first fungus to be officially considered an urgent antimicrobial resistance threat by the CDC (US), insights into the resistance mechanisms and evolutionary dynamics of C. auris are still scarce.

 

By using high-throughput in vitro experimental evolution with various antifungal drugs, we have obtained a library of resistant strains from four different clades. Through both genome and targeted sequencing, we have discovered novel mutations, especially for polyene resistance, which indicate new mechanisms of resistance and fitness compensation. For the validation of mutations, we have optimized a recyclable CRISPR/Cas9 tool for C. auris based on the C. albicans HIS-FLP system.

 

By mapping drug susceptibility responses of evolved strains across a library of several antifungals and repurposed drugs, we have discovered trends of cross-resistance and collateral sensitivity. Both phenomena have been extensively studied in tumors and bacteria but remain unexplored in fungi. In the light of these observations, we explore novel treatment schemes that prevent antifungal drug resistance development in C. auris and other pathogenic fungi.

The importance of antimicrobial resistance in medical mycology

Prior to the SARS-CoV-2 pandemic, antibiotic resistance was listed as the major global health care priority. Some analyses, including the O'Neill report, have predicted that deaths due to drug-resistant bacterial infections may eclipse the total number of cancer deaths by 2050. Although fungal infections remain in the shadow of public awareness, total attributable annual deaths are similar to, or exceeds, global mortalities due to malaria, tuberculosis or HIV. The impact of fungal infections has been exacerbated by the steady rise of antifungal drug resistant strains and species which reflects the widespread use of antifungals for prophylaxis and therapy, and in the case of azole resistance in Aspergillus, has been linked to the widespread agricultural use of antifungals. This review, based on a workshop hosted by the Medical Research Council and the University of Exeter, illuminates the problem of antifungal resistance and suggests how this growing threat might be mitigated.

A novel rRNA hybridization-based approach to rapid, accurate Candida identification directly from blood culture

Invasive fungal infections are increasingly common and carry high morbidity and mortality, yet fungal diagnostics lag behind bacterial diagnostics in rapidly identifying the causal pathogen. We previously devised a fluorescent hybridization-based assay to identify bacteria within hours directly from blood culture bottles without subculture, called phylogeny-informed rRNA-based strain identification (Phirst-ID). Here, we adapt this approach to unambiguously identify 11 common pathogenic Candida species, including C. auris, with 100% accuracy from laboratory culture (33 of 33 strains in a reference panel, plus 33 of 33 additional isolates tested in a validation panel). In a pilot study on 62 consecutive positive clinical blood cultures from two hospitals that showed yeast on Gram stain, Candida Phirst-ID matched the clinical laboratory result for 58 of 59 specimens represented in the 11-species reference panel, without misclassifying the 3 off-panel species. It also detected mixed Candida species in 2 of these 62 specimens, including the one discordant classification, that were not identified by standard clinical microbiology workflows; in each case the presence of both species was validated by both clinical and experimental data. Finally, in three specimens that grew both bacteria and yeast, we paired our prior bacterial probeset with this new Candida probeset to detect both pathogen types using Phirst-ID. This simple, robust assay can provide accurate Candida identification within hours directly from blood culture bottles, and the conceptual approach holds promise for pan-microbial identification in a single workflow.

LAY SUMMARY

Candida bloodstream infections cause considerable morbidity and mortality, yet slow diagnostics delay recognition, worsening patient outcomes. We develop and validate a novel molecular approach to accurately identify Candida species directly from blood culture one day faster than standard workflows.

Population genetics and microevolution of clinical Candida glabrata reveals recombinant sequence types and hyper-variation within mitochondrial genomes, virulence genes, and drug targets

Candida glabrata is the second most common etiological cause of worldwide systemic candidiasis in adult patients. Genome analysis of 68 isolates from 8 hospitals across Scotland, together with 83 global isolates, revealed insights into the population genetics and evolution of C. glabrata. Clinical isolates of C. glabrata from across Scotland are highly genetically diverse, including at least 19 separate sequence types that have been recovered previously in globally diverse locations, and 1 newly discovered sequence type. Several sequence types had evidence for ancestral recombination, suggesting transmission between distinct geographical regions has coincided with genetic exchange arising in new clades. Three isolates were missing MATα1, potentially representing a second mating type. Signatures of positive selection were identified in every sequence type including enrichment for epithelial adhesins thought to facilitate fungal adhesin to human epithelial cells. In patent microevolution was identified from 7 sets of recurrent cases of candidiasis, revealing an enrichment for nonsynonymous and frameshift indels in cell surface proteins. Microevolution within patients also affected epithelial adhesins genes, and several genes involved in drug resistance including the ergosterol synthesis gene ERG4 and the echinocandin target FKS1/2, the latter coinciding with a marked drop in fluconazole minimum inhibitory concentration. In addition to nuclear genome diversity, the C. glabrata mitochondrial genome was particularly diverse, with reduced conserved sequence and conserved protein-encoding genes in all nonreference ST15 isolates. Together, this study highlights the genetic diversity within the C. glabrata population that may impact virulence and drug resistance, and 2 major mechanisms generating this diversity: microevolution and genetic exchange/recombination.

Inter-species geographic signatures for tracing horizontal gene transfer and long-term persistence of carbapenem resistance

BACKGROUND

Carbapenem-resistant Enterobacterales (CRE) are an urgent global health threat. Inferring the dynamics of local CRE dissemination is currently limited by our inability to confidently trace the spread of resistance determinants to unrelated bacterial hosts. Whole-genome sequence comparison is useful for identifying CRE clonal transmission and outbreaks, but high-frequency horizontal gene transfer (HGT) of carbapenem resistance genes and subsequent genome rearrangement complicate tracing the local persistence and mobilization of these genes across organisms.

METHODS

To overcome this limitation, we developed a new approach to identify recent HGT of large, near-identical plasmid segments across species boundaries, which also allowed us to overcome technical challenges with genome assembly. We applied this to complete and near-complete genome assemblies to examine the local spread of CRE in a systematic, prospective collection of all CRE, as well as time- and species-matched carbapenem-susceptible Enterobacterales, isolated from patients from four US hospitals over nearly 5 years.

RESULTS

Our CRE collection comprised a diverse range of species, lineages, and carbapenem resistance mechanisms, many of which were encoded on a variety of promiscuous plasmid types. We found and quantified rearrangement, persistence, and repeated transfer of plasmid segments, including those harboring carbapenemases, between organisms over multiple years. Some plasmid segments were found to be strongly associated with specific locales, thus representing geographic signatures that make it possible to trace recent and localized HGT events. Functional analysis of these signatures revealed genes commonly found in plasmids of nosocomial pathogens, such as functions required for plasmid retention and spread, as well survival against a variety of antibiotic and antiseptics common to the hospital environment.

CONCLUSIONS

Collectively, the framework we developed provides a clearer, high-resolution picture of the epidemiology of antibiotic resistance importation, spread, and persistence in patients and healthcare networks.

Obligate sexual reproduction of a homothallic fungus closely related to the Cryptococcus pathogenic species complex

<title>eLife digest</title>.

Fungi are enigmatic organisms that flourish in soil, on decaying plants, or during infection of animals or plants. Growing in myriad forms, from single-celled yeast to multicellular molds and mushrooms, fungi have also evolved a variety of strategies to reproduce. Normally, fungi reproduce in one of two ways: either they reproduce asexually, with one individual producing a new individual identical to itself, or they reproduce sexually, with two individuals of different 'mating types' contributing to produce a new individual. However, individuals of some species exhibit 'homothallism' or self-fertility: these individuals can produce reproductive cells that are universally compatible, and therefore can reproduce sexually with themselves or with any other cell in the population.

Homothallism has evolved multiple times throughout the fungal kingdom, suggesting it confers advantage when population numbers are low or mates are hard to find. Yet some homothallic fungi been overlooked compared to heterothallic species, whose mating types have been well characterised. Understanding the genetic basis of homothallism and how it evolved in different species can provide insights into pathogenic species that cause fungal disease.

With that in mind, Passer, Clancey et al. explored the genetic basis of homothallism in Cryptococcus depauperatus, a close relative of C. neoformans, a species that causes fungal infections in humans. A combination of genetic sequencing techniques and experiments were applied to analyse, compare, and manipulate C. depauperatus' genome to see how this species evolved self-fertility.

Passer, Clancey et al. showed that C. depauperatus evolved the ability to reproduce sexually by itself via a unique evolutionary pathway. The result is a form of homothallism never reported in fungi before. C. depauperatus lost some of the genes that control mating in other species of fungi, and acquired genes from the opposing mating types of a heterothallic ancestor to become self-fertile.

Passer, Clancey et al. also found that, unlike other Cryptococcus species that switch between asexual and sexual reproduction, C. depauperatus grows only as long, branching filaments called hyphae, a sexual form. The species reproduces sexually with itself throughout its life cycle and is unable to produce a yeast (asexual) form, in contrast to other closely related species.

This work offers new insights into how different modes of sexual reproduction have evolved in fungi. It also provides another interesting case of how genome plasticity and evolutionary pressures can produce similar outcomes, homothallism, via different evolutionary paths. Lastly, assembling the complete genome of C. depauperatus will foster comparative studies between pathogenic and non-pathogenic Cryptococcus species.

MycoSNP: A Portable Workflow for Performing Whole-Genome Sequencing Analysis of Candida auris

Candida auris is an urgent public health threat characterized by high drug-resistant rates and rapid spread in healthcare settings worldwide. As part of the C. auris response, molecular surveillance has helped public health officials track the global spread and investigate local outbreaks. Here, we describe whole-genome sequencing analysis methods used for routine C. auris molecular surveillance in the United States; methods include reference selection, reference preparation, quality assessment and control of sequencing reads, read alignment, and single-nucleotide polymorphism calling and filtration. We also describe the newly developed pipeline MycoSNP, a portable workflow for performing whole-genome sequencing analysis of fungal organisms including C. auris.

Generating Complete Genome Assemblies of Candida auris

Genomic studies of Candida auris are underpinned by the generation of high-quality genome assemblies. These reference genomes have been essential for investigations of the evolution and epidemiology of this emerging fungal pathogen. In addition to genomic epidemiology studies of local outbreaks and analysis of the global emergence of this species, comparisons of genomes of isolates from the five major clades have revealed differences in gene content and genomic structure. Here, we provide a detailed protocol for generating complete genome assemblies for C. auris.

Genome-wide analysis of experimentally evolved Candida auris reveals multiple novel mechanisms of multidrug-resistance

Candida auris is globally recognized as an opportunistic fungal pathogen of high concern, due to its extensive multidrug-resistance (MDR). Still, molecular mechanisms of MDR are largely unexplored. This is the first account of genome wide evolution of MDR in C. auris obtained through serial in vitro exposure to azoles, polyenes and echinocandins. We show the stepwise accumulation of multiple novel mutations in genes known and unknown in antifungal drug resistance, albeit almost all new for C. auris. Echinocandin resistance was accompanied by a codon deletion in FKS1hot spot 1 and a substitution in FKS1 ‘novel’ hot spot 3. Mutations in ERG3 and CIS2 further increased the echinocandin MIC. Decreased azole susceptibility was linked to a mutation in transcription factor TAC1b and overexpression of the drug efflux pump Cdr1; a segmental duplication of chromosome 1 containing ERG11; and a whole chromosome 5 duplication, which contains TAC1b. The latter was associated with increased expression of ERG11, TAC1band CDR2, but not CDR1. The simultaneous emergence of nonsense mutations in ERG3 and ERG11 was shown to decrease amphotericin B susceptibility, accompanied with fluconazole cross resistance. A mutation in MEC3, a gene mainly known for its role in DNA damage homeostasis, further increased the polyene MIC. Overall, this study shows the alarming potential and diversity for MDR development in C. auris, even in a clade until now not associated with MDR (clade II),hereby stressing its clinical importance and the urge for future research.

Phenotypic Characterization and Comparative Genomics of the Melanin-Producing Yeast Exophiala lecanii-corni Reveals a Distinct Stress Tolerance Profile and Reduced Ribosomal Genetic Content

The black yeast Exophiala lecanii-corni of the order Chaetothyriales is notable for its ability to produce abundant quantities of DHN-melanin. While many other Exophiala species are frequent causal agents of human infection, E. lecanii-corni CBS 102400 lacks the thermotolerance requirements that enable pathogenicity, making it appealing for use in targeted functional studies and biotechnological applications. Here, we report the stress tolerance characteristics of E. lecanii-corni, with an emphasis on the influence of melanin on its resistance to various forms of stress. We find that E. lecanii-corni has a distinct stress tolerance profile that includes variation in resistance to temperature, osmotic, and oxidative stress relative to the extremophilic and pathogenic black yeast Exophiala dermatitidis. Notably, the presence of melanin substantially impacts stress resistance in E. lecanii-corni, while this was not found to be the case in E. dermatitidis. The cellular context, therefore, influences the role of melanin in stress protection. In addition, we present a detailed analysis of the E. lecanii-corni genome, revealing key differences in functional genetic content relative to other ascomycetous species, including a significant decrease in abundance of genes encoding ribosomal proteins. In all, this study provides insight into how genetics and physiology may underlie stress tolerance and enhances understanding of the genetic diversity of black yeasts.

In vivo emergence of high-level resistance during treatment reveals the first identified mechanism of amphotericin B resistance in Candida auris

OBJECTIVE

Candida auris has emerged as a healthcare-associated and multidrug-resistant fungal pathogen of great clinical concern. While as many as 50% of C. auris clinical isolates are reported to be resistant to amphotericin B, to date, no mechanisms contributing to this resistance have been identified. Here we describe the clinical case in which high-level amphotericin B resistance was acquired in vivo during therapy and undertake molecular and genetic studies to identify and characterize the genetic determinant of resistance.

METHODS

Whole genome sequencing was performed on four C. auris isolates obtained from a single patient case. Cas9-mediated genetic manipulations were then used to generate mutant strains harboring mutations of interest, and these strains were subsequently subjected to amphotericin B susceptibility testing, and comprehensive sterol profiling.

RESULTS

A novel mutation in C. auris sterol-methyltransferase gene, ERG6, was found to be associated with amphotericin B resistance, and this mutation alone conferred a >32-fold increase in amphotericin B resistance. Comprehensive sterol profiling revealed an abrogation of ergosterol biosynthesis and a corresponding accumulation of cholesta-type sterols in isolates and strains harboring the clinically-derived ERG6 mutation.

CONCLUSIONS

Together these findings definitively demonstrate mutations in C. auris ERG6 as the first identified mechanism of clinical amphotericin B resistance in C. auris and represent a significant step forward in the understanding of antifungal resistance in this emerging public health threat.

Diversity, multifaceted evolution, and facultative saprotrophism in the European Batrachochytrium salamandrivorans epidemic

While emerging fungi threaten global biodiversity, the paucity of fungal genome assemblies impedes thoroughly characterizing epidemics and developing effective mitigation strategies. Here, we generate de novo genomic assemblies for six outbreaks of the emerging pathogen Batrachochytrium salamandrivorans (Bsal). We reveal the European epidemic currently damaging amphibian populations to comprise multiple, highly divergent lineages demonstrating isolate-specific adaptations and metabolic capacities. In particular, we show extensive gene family expansions and acquisitions, through a variety of evolutionary mechanisms, and an isolate-specific saprotrophic lifecycle. This finding both explains the chytrid's ability to divorce transmission from host density, producing Bsal's enigmatic host population declines, and is a key consideration in developing successful mitigation measures.

Blastomycosis in Africa and the Middle East: A Comprehensive Review of Reported Cases and Reanalysis of Historical Isolates Based on Molecular Data

BACKGROUND

Blastomycosis has been reported from countries in Africa and the Middle East, but a decades-long debate has persisted regarding whether this is the same disease known in North America and caused by Blastomyces dermatitidis and Blastomyces gilchristii.

METHODS

We reviewed published cases of human and veterinary blastomycosis from Africa and the Middle East. We abstracted epidemiological and clinical features of cases, including sites of disease, diagnosis, management, outcomes, and, where available, genetic and antigenic typing of case isolates. In addition, we sequenced nucleic acids from 9 clinical isolates from Africa deposited in global collections as B. dermatitidis; for 5, we sequenced the internal transcribed spacer regions, and for the other 4 we sequenced the whole genomes.

RESULTS

We identified 172 unique human patients with blastomycosis, including 159 patients from 25 African countries and 12 patients from 5 Middle Eastern countries, and also identified 7 reports of veterinary blastomycosis. In humans, cutaneous disease predominated (n = 100/137, 73%), followed by pulmonary (n = 73/129, 57%) and osteoarticular involvement (n = 61/128, 48%). Unusual direct microscopy/histopathological presentations included short hyphal fragments in tissues (n = 23/129, 18%). There were 34 genotyped case isolates that comprised 4 species: Blastomyces percursus (n = 22, 65%), from 8 countries throughout all regions; Blastomyces emzantsi (n = 9, 26%), from South Africa; B. dermatitidis (n = 1, 3%), from the Democratic Republic of Congo; and B. gilchristii (n = 2, 6%), from South Africa and Zimbabwe.

CONCLUSIONS

Blastomycosis occurs throughout Africa and the Middle East and is caused predominantly by B. percursus and, at least in South Africa, B. emzantsi, resulting in distinct clinical and pathological patterns of disease.

Diagnostic Allele-Specific PCR for the Identification of Candida auris Clades

Candida auris is an opportunistic pathogenic yeast that emerged worldwide during the past decade. This fungal pathogen poses a significant public health threat due to common multidrug resistance (MDR), alarming hospital outbreaks, and frequent misidentification. Genomic analyses have identified five distinct clades that are linked to five geographic areas of origin and characterized by differences in several phenotypic traits such as virulence and drug resistance. Typing of C. auris strains and the identification of clades can be a powerful tool in molecular epidemiology and might be of clinical importance by estimating outbreak and MDR potential. As C. auris has caused global outbreaks, including in low-income countries, typing C. auris strains quickly and inexpensively is highly valuable. We report five allele-specific polymerase chain reaction (AS-PCR) assays for the identification of C. auris and each of the five described clades of C. auris based on conserved mutations in the internal transcribed spacer (ITS) rDNA region and a clade-specific gene cluster. This PCR method provides a fast, cheap, sequencing-free diagnostic tool for the identification of C. auris, C. auris clades, and potentially, the discovery of new clades.

Whole-genome sequencing of multiple isolates of Puccinia triticina reveals asexual lineages evolving by recurrent mutations

The wheat leaf rust fungus, Puccinia triticina Erikss., is a worldwide pathogen of tetraploid durum and hexaploid wheat. Many races of P. triticina differ for virulence to specific leaf rust resistance genes and are found in most wheat-growing regions of the world. Wheat cultivars with effective leaf rust resistance exert selection pressure on P. triticina populations for virulent race types. The objectives of this study were to examine whole-genome sequence data of 121 P. triticina isolates and to gain insight into race evolution. The collection included isolates comprising of many different race phenotypes collected worldwide from common and durum wheat. One isolate from wild wheat relative Aegilops speltoides and two from Ae. cylindrica were also included for comparison. Based on 121,907 informative variants identified relative to the reference Race 1-1 genome, isolates were clustered into 11 major lineages with 100% bootstrap support. The isolates were also grouped based on variation in 1311 predicted secreted protein genes. In gene-coding regions, all groups had high ratios of nonsynonymous to synonymous mutations and nonsense to readthrough mutations. Grouping of isolates based on two main variation principle components for either genome-wide variation or variation just within the secreted protein genes, indicated similar groupings. Variants were distributed across the entire genome, not just within the secreted protein genes. Our results suggest that recurrent mutation and selection play a major role in differentiation within the clonal lineages.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary-based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results show that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5–10 million reads per RNA-seq replicate. We also showed that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

Phylogenomic Analysis of a 55.1-kb 19-Gene Dataset Resolves a Monophyletic Fusarium that Includes the Fusarium solani Species Complex

Scientific communication is facilitated by a data-driven, scientifically sound taxonomy that considers the end-user's needs and established successful practice. In 2013, the Fusarium community voiced near unanimous support for a concept of Fusarium that represented a clade comprising all agriculturally and clinically important Fusarium species, including the F. solani species complex (FSSC). Subsequently, this concept was challenged in 2015 by one research group who proposed dividing the genus Fusarium into seven genera, including the FSSC described as members of the genus Neocosmospora, with subsequent justification in 2018 based on claims that the 2013 concept of Fusarium is polyphyletic. Here, we test this claim and provide a phylogeny based on exonic nucleotide sequences of 19 orthologous protein-coding genes that strongly support the monophyly of Fusarium including the FSSC. We reassert the practical and scientific argument in support of a genus Fusarium that includes the FSSC and several other basal lineages, consistent with the longstanding use of this name among plant pathologists, medical mycologists, quarantine officials, regulatory agencies, students, and researchers with a stake in its taxonomy. In recognition of this monophyly, 40 species described as genus Neocosmospora were recombined in genus Fusarium, and nine others were renamed Fusarium. Here the global Fusarium community voices strong support for the inclusion of the FSSC in Fusarium, as it remains the best scientific, nomenclatural, and practical taxonomic option available.

Clade-specific chromosomal rearrangements and loss of subtelomeric adhesins in Candida auris

Candida auris is an emerging fungal pathogen of rising concern due to global spread, the ability to cause healthcare-associated outbreaks, and antifungal resistance. Genomic analyses revealed that early contemporaneously detected cases of C. auris were geographically stratified into four major clades. While Clades I, III, and IV are responsible for ongoing outbreaks of invasive and multidrug-resistant infections, Clade II, also termed the East Asian clade, consists primarily of cases of ear infection, is often susceptible to all antifungal drugs, and has not been associated with outbreaks. Here, we generate chromosome-level assemblies of twelve isolates representing the phylogenetic breadth of these four clades and the only isolate described to date from Clade V. This Clade V genome is highly syntenic with those of Clades I, III, and IV, although the sequence is highly divergent from the other clades. Clade II genomes appear highly rearranged, with translocations occurring near GC-poor regions, and large subtelomeric deletions in most chromosomes, resulting in a substantially different karyotype. Rearrangements and deletion lengths vary across Clade II isolates, including two from a single patient, supporting ongoing genome instability. Deleted subtelomeric regions are enriched in Hyr/Iff-like cell-surface proteins, novel candidate cell wall proteins, and an ALS-like adhesin. Cell wall proteins from these families and other drug-related genes show clade-specific signatures of selection in Clades I, III, and IV. Subtelomeric dynamics and the conservation of cell surface proteins in the clades responsible for global outbreaks causing invasive infections suggest an explanation for the different phenotypes observed between clades.

Best practices on the differential expression analysis of multi-species RNA-seq

Advances in transcriptome sequencing allow for simultaneous interrogation of differentially expressed genes from multiple species originating from a single RNA sample, termed dual or multi-species transcriptomics. Compared to single-species differential expression analysis, the design of multi-species differential expression experiments must account for the relative abundances of each organism of interest within the sample, often requiring enrichment methods and yielding differences in total read counts across samples. The analysis of multi-species transcriptomics datasets requires modifications to the alignment, quantification, and downstream analysis steps compared to the single-species analysis pipelines. We describe best practices for multi-species transcriptomics and differential gene expression.

Amoeba Predation of Cryptococcus neoformans Results in Pleiotropic Changes to Traits Associated with Virulence

Amoeboid predators, such as amoebae, are proposed to select for survival traits in soil microbes such as Cryptococcus neoformans; these traits can also function in animal virulence by defeating phagocytic immune cells, such as macrophages. Consistent with this notion, incubation of various fungal species with amoebae enhanced their virulence, but the mechanisms involved are unknown. In this study, we exposed three strains of C. neoformans (1 clinical and 2 environmental) to predation by Acanthamoeba castellanii for prolonged times and then analyzed surviving colonies phenotypically and genetically. Surviving colonies comprised cells that expressed either pseudohyphal or yeast phenotypes, which demonstrated variable expression of traits associated with virulence, such as capsule size, urease production, and melanization. Phenotypic changes were associated with aneuploidy and DNA sequence mutations in some amoeba-passaged isolates, but not in others. Mutations in the gene encoding the oligopeptide transporter (CNAG_03013; OPT1) were observed among amoeba-passaged isolates from each of the three strains. Isolates derived from environmental strains gained the capacity for enhanced macrophage toxicity after amoeba selection and carried mutations on the CNAG_00570 gene encoding Pkr1 (AMP-dependent protein kinase regulator) but manifested reduced virulence in mice because they elicited more effective fungal-clearing immune responses. Our results indicate that C. neoformans survival under constant amoeba predation involves the generation of strains expressing pleiotropic phenotypic and genetic changes. Given the myriad potential predators in soils, the diversity observed among amoeba-selected strains suggests a bet-hedging strategy whereby variant diversity increases the likelihood that some will survive predation.IMPORTANCECryptococcus neoformans is a ubiquitous environmental fungus that is also a leading cause of fatal fungal infection in humans, especially among immunocompromised patients. A major question in the field is how an environmental yeast such as C. neoformans becomes a human pathogen when it has no need for an animal host in its life cycle. Previous studies showed that C. neoformans increases its pathogenicity after interacting with its environmental predator amoebae. Amoebae, like macrophages, are phagocytic cells that are considered an environmental training ground for pathogens to resist macrophages, but the mechanism by which C. neoformans changes its virulence through interactions with protozoa is unknown. Our study indicates that fungal survival in the face of amoeba predation is associated with the emergence of pleiotropic phenotypic and genomic changes that increase the chance of fungal survival, with this diversity suggesting a bet-hedging strategy to ensure that some forms survive.

Genome-Wide Analysis of Experimentally Evolved Candida auris Reveals Multiple Novel Mechanisms of Multidrug Resistance

Candida auris is globally recognized as an opportunistic fungal pathogen of high concern, due to its extensive multidrug resistance (MDR). Still, molecular mechanisms of MDR are largely unexplored. This is the first account of genome-wide evolution of MDR in C. auris obtained through serial in vitro exposure to azoles, polyenes, and echinocandins. We show the stepwise accumulation of copy number variations and novel mutations in genes both known and unknown in antifungal drug resistance. Echinocandin resistance was accompanied by a codon deletion in FKS1 hot spot 1 and a substitution in FKS1 "novel" hot spot 3. Mutations in ERG3 and CIS2 further increased the echinocandin MIC. Decreased azole susceptibility was linked to a mutation in transcription factor TAC1b and overexpression of the drug efflux pump Cdr1, a segmental duplication of chromosome 1 containing ERG11, and a whole chromosome 5 duplication, which contains TAC1b The latter was associated with increased expression of ERG11, TAC1b, and CDR2 but not CDR1 The simultaneous emergence of nonsense mutations in ERG3 and ERG11 was shown to decrease amphotericin B susceptibility, accompanied with fluconazole cross-resistance. A mutation in MEC3, a gene mainly known for its role in DNA damage homeostasis, further increased the polyene MIC. Overall, this study shows the alarming potential for and diversity of MDR development in C. auris, even in a clade until now not associated with MDR (clade II), stressing its clinical importance and the urge for future research.IMPORTANCECandida auris is a recently discovered human fungal pathogen and has shown an alarming potential for developing multi- and pan-resistance toward all classes of antifungals most commonly used in the clinic. Currently, C. auris has been globally recognized as a nosocomial pathogen of high concern due to this evolutionary potential. So far, this is the first study in which the stepwise progression of multidrug resistance (MDR) in C. auris is monitored in vitro Multiple novel mutations in known resistance genes and genes previously not or vaguely associated with drug resistance reveal rapid MDR evolution in a C. auris clade II isolate. Additionally, this study shows that in vitro experimental evolution can be a powerful tool to discover new drug resistance mechanisms, although it has its limitations.

Genomic insights into the host specific adaptation of the Pneumocystis genus

Pneumocystis jirovecii, the fungal agent of human Pneumocystis pneumonia, is closely related to macaque Pneumocystis. Little is known about other Pneumocystis species in distantly related mammals, none of which are capable of establishing infection in humans. The molecular basis of host specificity in Pneumocystis remains unknown as experiments are limited due to an inability to culture any species in vitro. To explore Pneumocystis evolutionary adaptations, we have sequenced the genomes of species infecting macaques, rabbits, dogs and rats and compared them to available genomes of species infecting humans, mice and rats. Complete whole genome sequence data enables analysis and robust phylogeny, identification of important genetic features of the host adaptation, and estimation of speciation timing relative to the rise of their mammalian hosts. Our data reveals insights into the evolution of P. jirovecii, the sole member of the genus able to infect humans.

Cissé, Ma et al. utilize genomic data from Pneumocystis species infecting macaques, rabbit, dogs and rats to investigate the molecular basis of host specificity in Pneumocystis. Their analyses provide insight to the specific adaptations enabling the infection of humans by P. jirovecii.

Genomic epidemiology of a Cryptococcus neoformans case cluster in Glasgow, Scotland, 2018

In 2018, a cluster of two cases of cryptococcosis occurred at the Queen Elizabeth University Hospital (QEUH) in Glasgow, Scotland (UK). It was postulated that these cases may have been linked to pigeon droppings found on the hospital site, given there have been previous reports of Cryptococcus neoformans associated with pigeon guano. Although some samples of pigeon guano taken from the site yielded culturable yeast from genera related to Cryptococcus, they have since been classified as Naganishia or Papiliotrema spp., and no isolates of C. neoformans were recovered from either the guano or subsequent widespread air sampling. In an attempt to further elucidate any possible shared source of the clinical isolates, we used whole-genome sequencing and phylogenetic analysis to examine the relationship of the two Cryptococcus isolates from the QEUH cases, along with two isolates from sporadic cases treated at a different Glasgow hospital earlier in 2018. Our work demonstrated that these four clinical isolates were not clonally related; while all isolates were from the VNI global lineage and of the same mating type (MATα), the genotypes of the two QEUH isolates were separated by 1885 base changes and belonged to different sub-lineages, recently described as the intercontinental sub-clades VNIa-93 and VNIa-5. In contrast, one of the two sporadic 2018 clinical isolates was determined to belong to the VNIb sub-lineage and the other classified as a VNIV/VNI hybrid. Our work demonstrated that the two 2018 QEUH isolates and the two prior C. neoformans clinical isolates were all genetically distinct. It was not possible to determine whether the QEUH genotypes stemmed from independent sources or from the same source, i.e. pigeons carrying different genotypes, but it should be noted that whilst members of allied genera within the Tremellomycetes were isolated from the hospital environment, there were no environmental isolations of C. neoformans.

Comparative genomics of white and opaque cell states supports an epigenetic mechanism of phenotypic switching in Candida albicans

Several Candida species can undergo a heritable and reversible transition from a 'white' state to a mating proficient 'opaque' state. This ability relies on highly interconnected transcriptional networks that control cell-type-specific gene expression programs over multiple generations. Candida albicans, the most prominent pathogenic Candida species, provides a well-studied paradigm for the white-opaque transition. In this species, a network of at least eight transcriptional regulators controls the balance between white and opaque states that have distinct morphologies, transcriptional profiles, and physiological properties. Given the reversible nature and the high frequency of white-opaque transitions, it is widely assumed that this switch is governed by epigenetic mechanisms that occur independently of any changes in DNA sequence. However, a direct genomic comparison between white and opaque cells has yet to be performed. Here, we present a whole-genome comparative analysis of C. albicans white and opaque cells. This analysis revealed rare genetic changes between cell states, none of which are linked to white-opaque switching. This result is consistent with epigenetic mechanisms controlling cell state differentiation in C. albicans and provides direct evidence against a role for genetic variation in mediating the switch.

Landscape of gene expression variation of natural isolates of Cryptococcus neoformans in response to biologically relevant stresses

Cryptococcus neoformans is an opportunistic fungal pathogen that at its peak epidemic levels caused an estimated million cases of cryptococcal meningitis per year worldwide. This species can grow in diverse environmental (trees, soil and bird excreta) and host niches (intracellular microenvironments of phagocytes and free-living in host tissues). The genetic basic for adaptation to these different conditions is not well characterized, as most experimental work has relied on a single reference strain of C. neoformans. To identify genes important for yeast infection and disease progression, we profiled the gene expression of seven C. neoformans isolates grown in five representative in vitro environmental and in vivo conditions. We characterized gene expression differences using RNA-Seq (RNA sequencing), comparing clinical and environmental isolates from two of the major lineages of this species, VNI and VNBI. These comparisons highlighted genes showing lineage-specific expression that are enriched in subtelomeric regions and in lineage-specific gene clusters. By contrast, we find few expression differences between clinical and environmental isolates from the same lineage. Gene expression specific to in vivo stages reflects available nutrients and stresses, with an increase in fungal metabolism within macrophages, and an induction of ribosomal and heat-shock gene expression within the subarachnoid space. This study provides the widest view to date of the transcriptome variation of C. neoformans across natural isolates, and provides insights into genes important for in vitro and in vivo growth stages.

Genomic diversity of the human pathogen Paracoccidioides across the South American continent

Paracoccidioidomycosis (PCM) is a life-threatening systemic mycosis widely reported in the Gran Chaco ecosystem. The disease is caused by different species from the genus Paracoccidioides, which are all endemic to South and Central America. Here, we sequenced and analyzed 31 isolates of Paracoccidioides across South America, with particular focus on isolates from Argentina and Paraguay. The de novo sequenced isolates were compared with publicly available genomes. Phylogenetics and population genomics revealed that PCM in Argentina and Paraguay is caused by three distinct Paracoccidioides genotypes, P. brasiliensis (S1a and S1b) and P. restrepiensis (PS3). P. brasiliensis S1a isolates from Argentina are frequently associated with chronic forms of the disease. Our results suggest the existence of extensive molecular polymorphism among Paracoccidioides species, and provide a framework to begin to dissect the connection between genotypic differences in the pathogen and the clinical outcomes of the disease.

Understanding the Emergence of Multidrug-Resistant Candida: Using Whole-Genome Sequencing to Describe the Population Structure of Candida haemulonii Species Complex

The recent emergence of a multidrug-resistant yeast, Candida auris, has drawn attention to the closely related species from the Candida haemulonii complex that include C. haemulonii, Candida duobushaemulonii, Candida pseudohaemulonii, and the recently identified Candida vulturna. Here, we used antifungal susceptibility testing and whole-genome sequencing (WGS) to investigate drug resistance and genetic diversity among isolates of C. haemulonii complex from different geographic areas in order to assess population structure and the extent of clonality among strains. Although most isolates of all four species were genetically distinct, we detected evidence of the in-hospital transmission of C. haemulonii and C. duobushaemulonii in one hospital in Panama, indicating that these species are also capable of causing outbreaks in healthcare settings. We also detected evidence of the rising azole resistance among isolates of C. haemulonii and C. duobushaemulonii in Colombia, Panama, and Venezuela linked to substitutions in ERG11 gene as well as amplification of this gene in C. haemulonii in isolates in Colombia suggesting the presence of evolutionary pressure for developing azole resistance in this region. Our results demonstrate that these species need to be monitored as possible causes of outbreaks of invasive infection.

Threats Posed by the Fungal Kingdom to Humans, Wildlife, and Agriculture

The fungal kingdom includes at least 6 million eukaryotic species and is remarkable with respect to its profound impact on global health, biodiversity, ecology, agriculture, manufacturing, and biomedical research. Approximately 625 fungal species have been reported to infect vertebrates, 200 of which can be human associated, either as commensals and members of our microbiome or as pathogens that cause infectious diseases. These organisms pose a growing threat to human health with the global increase in the incidence of invasive fungal infections, prevalence of fungal allergy, and the evolution of fungal pathogens resistant to some or all current classes of antifungals. More broadly, there has been an unprecedented and worldwide emergence of fungal pathogens affecting animal and plant biodiversity. Approximately 8,000 species of fungi and Oomycetes are associated with plant disease. Indeed, across agriculture, such fungal diseases of plants include new devastating epidemics of trees and jeopardize food security worldwide by causing epidemics in staple and commodity crops that feed billions. Further, ingestion of mycotoxins contributes to ill health and causes cancer. Coordinated international research efforts, enhanced technology translation, and greater policy outreach by scientists are needed to more fully understand the biology and drivers that underlie the emergence of fungal diseases and to mitigate against their impacts. Here, we focus on poignant examples of emerging fungal threats in each of three areas: human health, wildlife biodiversity, and food security.

Tracing the Evolutionary History and Global Expansion of Candida auris Using Population Genomic Analyses

Candida auris has emerged globally as a multidrug-resistant yeast that can spread via nosocomial transmission. An initial phylogenetic study of isolates from Japan, India, Pakistan, South Africa, and Venezuela revealed four populations (clades I, II, III, and IV) corresponding to these geographic regions. Since this description, C. auris has been reported in more than 30 additional countries. To trace this global emergence, we compared the genomes of 304 C. auris isolates from 19 countries on six continents. We found that four predominant clades persist across wide geographic locations. We observed phylogeographic mixing in most clades; clade IV, with isolates mainly from South America, demonstrated the strongest phylogeographic substructure. C. auris isolates from two clades with opposite mating types were detected contemporaneously in a single health care facility in Kenya. We estimated a Bayesian molecular clock phylogeny and dated the origin of each clade within the last 360 years; outbreak-causing clusters from clades I, III, and IV originated 36 to 38 years ago. We observed high rates of antifungal resistance in clade I, including four isolates resistant to all three major classes of antifungals. Mutations that contribute to resistance varied between the clades, with Y132F in ERG11 as the most widespread mutation associated with azole resistance and S639P in FKS1 for echinocandin resistance. Copy number variants in ERG11 predominantly appeared in clade III and were associated with fluconazole resistance. These results provide a global context for the phylogeography, population structure, and mechanisms associated with antifungal resistance in C. aurisIMPORTANCE In less than a decade, C. auris has emerged in health care settings worldwide; this species is capable of colonizing skin and causing outbreaks of invasive candidiasis. In contrast to other Candida species, C. auris is unique in its ability to spread via nosocomial transmission and its high rates of drug resistance. As part of the public health response, whole-genome sequencing has played a major role in characterizing transmission dynamics and detecting new C. auris introductions. Through a global collaboration, we assessed genome evolution of isolates of C. auris from 19 countries. Here, we described estimated timing of the expansion of each C. auris clade and of fluconazole resistance, characterized discrete phylogeographic population structure of each clade, and compared genome data to sensitivity measurements to describe how antifungal resistance mechanisms vary across the population. These efforts are critical for a sustained, robust public health response that effectively utilizes molecular epidemiology.

Tracking a Global Threat: a New Genotyping Method for Candida auris

Over the past decade, Candida auris has emerged as an urgent threat to public health. Initially reported from cases of ear infections in Japan and Korea, C. auris has since been detected around the world. While whole-genome sequencing has been extensively used to trace the genetic relationships of the global emergence and local outbreaks, a recent report in mBio describes a targeted genotyping method as a rapid and inexpensive method for classifying C. auris isolates (T. de Groot, Y. Puts, I. Berrio, A. Chowdhary, and J. F. Meis, mBio 11:e02971-19, https://doi.org/10.1128/mBio.02971-19, 2020).

Diversity and Complexity of the Large Surface Protein Family in the Compacted Genomes of Multiple Pneumocystis Species

Pneumocystis, a major opportunistic pathogen in patients with a broad range of immunodeficiencies, contains abundant surface proteins encoded by a multicopy gene family, termed the major surface glycoprotein (Msg) gene superfamily. This superfamily has been identified in all Pneumocystis species characterized to date, highlighting its important role in Pneumocystis biology. In this report, through a comprehensive and in-depth characterization of 459 msg genes from 7 Pneumocystis species, we demonstrate, for the first time, the phylogeny and evolution of conserved domains in Msg proteins and provide a detailed description of the classification, unique characteristics, and phylogenetic relatedness of five Msg families. We further describe, for the first time, the relative expression levels of individual msg families in two rodent Pneumocystis species, the substantial variability of the msg repertoires in P. carinii from laboratory and wild rats, and the distinct features of the expression site for the classic msg genes in Pneumocystis from 8 mammalian host species. Our analysis suggests multiple functions for this superfamily rather than just conferring antigenic variation to allow immune evasion as previously believed. This study provides a rich source of information that lays the foundation for the continued experimental exploration of the functions of the Msg superfamily in Pneumocystis biology.IMPORTANCEPneumocystis continues to be a major cause of disease in humans with immunodeficiency, especially those with HIV/AIDS and organ transplants, and is being seen with increasing frequency worldwide in patients treated with immunodepleting monoclonal antibodies. Annual health care associated with Pneumocystis pneumonia costs ∼$475 million dollars in the United States alone. In addition to causing overt disease in immunodeficient individuals, Pneumocystis can cause subclinical infection or colonization in healthy individuals, which may play an important role in species preservation and disease transmission. Our work sheds new light on the diversity and complexity of the msg superfamily and strongly suggests that the versatility of this superfamily reflects multiple functions, including antigenic variation to allow immune evasion and optimal adaptation to host environmental conditions to promote efficient infection and transmission. These findings are essential to consider in developing new diagnostic and therapeutic strategies.