The pathobiology of human fungal infections.

Human fungal infections are a historically neglected area of disease research, yet they cause more than 1.5 million deaths every year. Our understanding of the pathophysiology of these infections has increased considerably over the past decade, through major insights into both the host and pathogen factors that contribute to the phenotype and severity of these diseases. Recent studies are revealing multiple mechanisms by which fungi modify and manipulate the host, escape immune surveillance and generate complex comorbidities. Although the emergence of fungal strains that are less susceptible to antifungal drugs or that rapidly evolve drug resistance is posing new threats, greater understanding of immune mechanisms and host susceptibility factors is beginning to offer novel immunotherapeutic options for the future. In this Review, we provide a broad and comprehensive overview of the pathobiology of human fungal infections, focusing specifically on pathogens that can cause invasive life-threatening infections, highlighting recent discoveries from the pathogen, host and clinical perspectives. We conclude by discussing key future challenges including antifungal drug resistance, the emergence of new pathogens and new developments in modern medicine that are promoting susceptibility to infection.

Revealing the genome of the microsporidian Vairimorpha bombi, a potential driver of bumble bee declines in North America.

Pollinators are vital for food security and the maintenance of terrestrial ecosystems. Bumblebees are important pollinators across northern temperate, arctic, and alpine ecosystems, yet are in decline across the globe. Vairimorpha bombi is a parasite belonging to the fungal class Microsporidia that has been implicated in the rapid decline of bumblebees in North America, where it may be an emerging infectious disease. To investigate the evolutionary basis of pathogenicity of V. bombi, we sequenced and assembled its genome using Oxford Nanopore and Illumina technologies and performed phylogenetic and genomic evolutionary analyses. The genome assembly for V. bombi is 4.73 Mb, from which we predicted 1,870 protein-coding genes and 179 tRNA genes. The genome assembly has low repetitive content and low GC content. V. bombi's genome assembly is the smallest of the Vairimorpha and closely related Nosema genera, but larger than those found in the Encephalitozoon and Ordospora sister clades. Orthology and phylogenetic analysis revealed 18 core conserved single-copy microsporidian genes including the histone acetyltransferase (HAT) GCN5. Surprisingly, V. bombi was unique to the microsporidia in not encoding the second predicted HAT ESA1. The V. bombi genome assembly annotation included 265 unique genes (i.e. not predicted in other microsporidia genome assemblies), 20% of which encode a secretion signal, which is a significant enrichment. Intriguingly, of the 36 microsporidian genomes we analyzed, 26 also had a significant enrichment of secreted signals encoded by unique genes, ranging from 6 to 71% of those predicted genes. These results suggest that microsporidia are under selection to generate and purge diverse and unique genes encoding secreted proteins, potentially contributing to or facilitating infection of their diverse hosts. Furthermore, V. bombi has 5/7 conserved spore wall proteins (SWPs) with its closest relative V. ceranae (that primarily infects honeybees), while also uniquely encoding four additional SWPs. This gene class is thought to be essential for infection, providing both environmental protection and recognition and uptake into the host cell. Together, our results show that SWPs and unique genes encoding a secretion signal are rapidly evolving in the microsporidia, suggesting that they underpin key pathobiological traits including host specificity and pathogenicity.

Strain and temperature dependent aggregation of Candida auris is attenuated by inhibition of surface amyloid proteins.

Candida auris is a multi-drug resistant human fungal pathogen that has become a global threat to human health due to its drug resistant phenotype, persistence in the hospital environment and propensity for patient to patient spread. Isolates display variable aggregation that may affect the relative virulence of strains. Therefore, dissection of this phenotype has gained substantial interest in recent years. We studied eight clinical isolates from four different clades (I-IV); four of which had a strongly aggregating phenotype and four of which did not. Genome analysis identified polymorphisms associated with loss of cell surface proteins were enriched in weakly-aggregating strains. Additionally, we identified down-regulation of chitin synthase genes involved in the synthesis of the chitinous septum. Characterisation of the cells revealed no ultrastructural defects in cytokinesis or cell separation in aggregating isolates. Strongly and weakly aggregating strains did not differ in net surface charge or in cell surface hydrophobicity. The capacity for aggregation and for adhesion to polystyrene microspheres were also not correlated. However, aggregation and extracellular matrix formation were all increased at higher growth temperatures, and treatment with the amyloid protein inhibitor Thioflavin-T markedly attenuated aggregation. Genome analysis further indicated strain specific differences in the genome content of GPI-anchored proteins including those encoding genes with the potential to form amyloid proteins. Collectively our data suggests that aggregation is a complex strain and temperature dependent phenomenon that may be linked in part to the ability to form extracellular matrix and cell surface amyloids.

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.

Two-speed genome evolution drives pathogenicity in fungal pathogens of animals.

The origins and evolution of virulence in amphibian-infecting chytrids Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are largely unknown. Here, we use deep nanopore sequencing of Bsal and comparative genomics against 21 high-quality genome assemblies that span the fungal Chytridiomycota. We discover that Bsal has the most repeat-rich genome of the Chytridiomycota, comprising 40.9% repetitive elements; this genome has expanded to more than 3× the length of its conspecific Bd, with autonomous and fully functional LTR/Gypsy elements contributing significantly to the expansion. The M36 metalloprotease virulence factors are highly expanded (n = 177) in Bsal, most of which (53%) are flanked by transposable elements, suggesting they have a repeat-associated expansion. We find enrichment upstream of M36 metalloprotease genes of three novel repeat families belonging to the repeat superfamily of LINEs that are implicated with gene copy number variations. Additionally, Bsal has a highly compartmentalized genome architecture, with virulence factors enriched in gene-sparse/repeat-rich compartments, while core conserved genes are enriched in gene-rich/repeat-poor compartments. Genes upregulated during infection are primarily found in the gene-sparse/repeat-rich compartment in both Bd and Bsal. Furthermore, genes with signatures of positive selection in Bd are enriched in repeat-rich regions, suggesting these regions are a cradle for the evolution of chytrid pathogenicity. These are the hallmarks of two-speed genome evolution, and this study provides evidence of two-speed genomes in an animal pathogen, shedding light on the evolution of fungal pathogens of vertebrates driving global declines and extinctions.

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.

HaplotypeTools: a toolkit for accurately identifying recombination and recombinant genotypes.

BACKGROUND: Identifying haplotypes is central to sequence analysis in diploid or polyploid genomes. Despite this, there remains a lack of research and tools designed for physical phasing and its downstream analysis. RESULTS: HaplotypeTools is a new toolset to phase variant sites using VCF and BAM files and to analyse phased VCFs. Phasing is achieved via the identification of reads overlapping ≥ 2 heterozygous positions and then extended by additional reads, a process that can be parallelized across a computer cluster. HaplotypeTools includes various utility scripts for downstream analysis including crossover detection and phylogenetic placement of haplotypes to other lineages or species. HaplotypeTools was assessed for accuracy against WhatsHap using simulated short and long reads, demonstrating higher accuracy, albeit with reduced haplotype length. HaplotypeTools was also tested on real Illumina data to determine the ancestry of hybrid fungal isolate Batrachochytrium dendrobatidis (Bd) SA-EC3, finding 80% of haplotypes across the genome phylogenetically cluster with parental lineages BdGPL (39%) and BdCAPE (41%), indicating those are the parental lineages. Finally, ~ 99% of phasing was conserved between overlapping phase groups between SA-EC3 and either parental lineage, indicating mitotic gene conversion/parasexuality as the mechanism of recombination for this hybrid isolate. HaplotypeTools is open source and freely available from https://github.com/rhysf/HaplotypeTools under the MIT License. CONCLUSIONS: HaplotypeTools is a powerful resource for analyzing hybrid or recombinant diploid or polyploid genomes and identifying parental ancestry for sub-genomic regions.

Drivers of methicillin-resistant Staphylococcus aureus (MRSA) lineage replacement in China.

BACKGROUND: Methicillin-resistant Staphylococcus aureus (MRSA) is a major nosocomial pathogen subdivided into lineages termed sequence types (STs). Since the 1950s, successive waves of STs have appeared and replaced previously dominant lineages. One such event has been occurring in China since 2013, with community-associated (CA-MRSA) strains including ST59 largely replacing the previously dominant healthcare-associated (HA-MRSA) ST239. We previously showed that ST59 isolates tend to have a competitive advantage in growth experiments against ST239. However, the underlying genomic and phenotypic drivers of this replacement event are unclear. METHODS: Here, we investigated the replacement of ST239 using whole-genome sequencing data from 204 ST239 and ST59 isolates collected in Chinese hospitals between 1994 and 2016. We reconstructed the evolutionary history of each ST and considered two non-mutually exclusive hypotheses for ST59 replacing ST239: antimicrobial resistance (AMR) profile and/or ability to colonise and persist in the environment through biofilm formation. We also investigated the differences in cytolytic activity, linked to higher virulence, between STs. We performed an association study using the presence and absence of accessory virulence genes. RESULTS: ST59 isolates carried fewer AMR genes than ST239 and showed no evidence of evolving towards higher AMR. Biofilm production was marginally higher in ST59 overall, though this effect was not consistent across sub-lineages so is unlikely to be a sole driver of replacement. Consistent with previous observations of higher virulence in CA-MRSA STs, we observed that ST59 isolates exhibit significantly higher cytolytic activity than ST239 isolates, despite carrying on average fewer putative virulence genes. Our association study identified the chemotaxis inhibitory protein (chp) as a strong candidate for involvement in the increased virulence potential of ST59. We experimentally validated the role of chp in increasing the virulence potential of ST59 by creating Δchp knockout mutants, confirming that ST59 can carry chp without a measurable impact on fitness. CONCLUSIONS: Our results suggest that the ongoing replacement of ST239 by ST59 in China is not associated to higher AMR carriage or biofilm production. However, the increase in ST59 prevalence is concerning since it is linked to a higher potential for virulence, aided by the carriage of the chp gene.

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.

Pathogen and host genetics underpinning cryptococcal disease.

Cryptococcosis is a severe fungal disease causing 220,000 cases of cryptococcal meningitis yearly. The etiological agents of cryptococcosis are taxonomically grouped into at least two species complexes belonging to the genus Cryptococcus. All of these yeasts are environmentally ubiquitous fungi (often found in soil, leaves and decaying wood, tree hollows, and associated with bird feces especially pigeon guano). Infection in a range of animals including humans begins following inhalation of spores or aerosolized yeasts. Recent advances provide fundamental insights into the factors from both the pathogen and its hosts which influence pathogenesis and disease. The complex interactions leading to disease in mammalian hosts have also updated from the availability of better genomic tools and datasets. In this review, we discuss recent genetic research on Cryptococcus, covering the epidemiology, ecology, and evolution of Cryptococcus pathogenic species. We also discuss the insights into the host immune response obtained from the latest genetic modified host models as well as insights from monogenic disorders in humans. Finally we highlight outstanding questions that can be answered in the near future using bioinformatics and genomic tools.

Nine Things Genomics Can Tell Us About Candida auris.

Candida auris is a recently emerged multidrug-resistant fungal pathogen causing severe illness in hospitalized patients. C. auris is most closely related to a few environmental or rarely observed but cosmopolitan Candida species. However, C. auris is unique in the concern it is generating among public health agencies for its rapid emergence, difficulty to treat, and the likelihood for further and more extensive outbreaks and spread. To date, five geographically distributed and genetically divergent lineages have been identified, none of which includes isolates that were collected prior to 1996. Indeed, C. auris' ecological niche(s) and emergence remain enigmatic, although a number of hypotheses have been proposed. Recent genomic and transcriptomic work has also identified a variety of gene and chromosomal features that may have conferred C. auris with several important clinical phenotypes including its drug-resistance and growth at high temperatures. In this review we discuss nine major lines of enquiry into C. auris that big-data technologies and analytical approaches are beginning to answer.

A New Lineage of Cryptococcus gattii (VGV) Discovered in the Central Zambezian Miombo Woodlands.

We discovered a new lineage of the globally important fungal pathogen Cryptococcus gattii on the basis of analysis of six isolates collected from three locations spanning the Central Miombo Woodlands of Zambia, Africa. All isolates were from environments (middens and tree holes) that are associated with a small mammal, the African hyrax. Phylogenetic and population genetic analyses confirmed that these isolates form a distinct, deeply divergent lineage, which we name VGV. VGV comprises two subclades (A and B) that are capable of causing mild lung infection with negligible neurotropism in mice. Comparing the VGV genome to previously identified lineages of C. gattii revealed a unique suite of genes together with gene loss and inversion events. However, standard URA5 restriction fragment length polymorphism (RFLP) analysis could not distinguish between VGV and VGIV isolates. We therefore developed a new URA5 RFLP method that can reliably identify the newly described lineage. Our work highlights how sampling understudied ecological regions alongside genomic and functional characterization can broaden our understanding of the evolution and ecology of major global pathogens.IMPORTANCECryptococcus gattii is an environmental pathogen that causes severe systemic infection in immunocompetent individuals more often than in immunocompromised humans. Over the past 2 decades, researchers have shown that C. gattii falls within four genetically distinct major lineages. By combining field work from an understudied ecological region (the Central Miombo Woodlands of Zambia, Africa), genome sequencing and assemblies, phylogenetic and population genetic analyses, and phenotypic characterization (morphology, histopathological, drug-sensitivity, survival experiments), we discovered a hitherto unknown lineage, which we name VGV (variety gattii five). The discovery of a new lineage from an understudied ecological region has far-reaching implications for the study and understanding of fungal pathogens and diseases they cause.

Evaluation of High-Resolution Melting for Rapid Differentiation of Phytophthora Hybrids and Their Parental Species.

Phytophthora species hybrids have been repeatedly reported as causing damaging diseases to cultivated and wild plants. Two known hybrids, P. andina and P. × pelgrandis, are pathogens of Solanaceae and ornamentals, respectively, although the extent of their host ranges are unknown. P. andina emerged from hybridization of P. infestans and an unidentified related species, whereas P. × pelgrandis emerged from P. nicotianae and P. cactorum. Considering that hybrids and parental species can coexist in the same regions and to distinguish them usually requires cloning or whole genome sequencing, we aimed to develop a rapid tool to distinguish them. Specifically, we used high-resolution melting (HRM) assays to differentiate genotypes based on their amplicon melting profiles. We designed primers for P. × pelgrandis and parental species based on available sequences of P. nicotianae and P. cactorum nuclear genes containing polymorphisms between species. For P. andina, heterozygous sites from Illumina short reads were used for the same purpose. We identified multiple amplicons exhibiting differences in melting curves between parental species and hybrids. We propose HRM as a rapid method for differentiation of P. andina and P. × pelgrandis hybrids from parental species that could be employed to advance research on these pathogens.

Rapid phenotypic evolution in multidrug-resistant Klebsiella pneumoniae hospital outbreak strains.

Carbapenem-resistant Klebsiella pneumoniae (CRKP) increasingly cause high-mortality outbreaks in hospital settings globally. Following a patient fatality at a hospital in Beijing due to a blaKPC-2-positive CRKP infection, close monitoring was put in place over the course of 14 months to characterize all blaKPC-2-positive CRKP in circulation in the hospital. Whole genome sequences were generated for 100 isolates from blaKPC-2-positive isolates from infected patients, carriers and the hospital environment. Phylogenetic analyses identified a closely related cluster of 82 sequence type 11 (ST11) isolates circulating in the hospital for at least a year prior to admission of the index patient. The majority of inferred transmissions for these isolates involved patients in intensive care units. Whilst the 82 ST11 isolates collected during the surveillance effort all had closely related chromosomes, we observed extensive diversity in their antimicrobial resistance (AMR) phenotypes. We were able to reconstruct the major genomic changes underpinning this variation in AMR profiles, including multiple gains and losses of entire plasmids and recombination events between plasmids, including transposition of blaKPC-2. We also identified specific cases where variation in plasmid copy number correlated with the level of phenotypic resistance to drugs, suggesting that the number of resistance elements carried by a strain may play a role in determining the level of AMR. Our findings highlight the epidemiological value of whole genome sequencing for investigating multi-drug-resistant hospital infections and illustrate that standard typing schemes cannot capture the extraordinarily fast genome evolution of CRKP isolates.

Genomic insights into multidrug-resistance, mating and virulence in Candida auris and related emerging species.

Candida auris is an emergent multidrug-resistant fungal pathogen causing increasing reports of outbreaks. While distantly related to C. albicans and C. glabrata, C. auris is closely related to rarely observed and often multidrug-resistant species from the C. haemulonii clade. Here, we analyze near complete genome assemblies for the four C. auris clades and three related species, and map intra- and inter-species rearrangements across the seven chromosomes. Using RNA-Seq-guided gene predictions, we find that most mating and meiosis genes are conserved and that clades contain either the MTLa or MTLα mating loci. Comparing the genomes of these emerging species to those of other Candida species identifies genes linked to drug resistance and virulence, including expanded families of transporters and lipases, as well as mutations and copy number variants in ERG11. Gene expression analysis identifies transporters and metabolic regulators specific to C. auris and those conserved with related species which may contribute to differences in drug response in this emerging fungal clade.

Transcriptional Heterogeneity of Cryptococcus gattii VGII Compared with Non-VGII Lineages Underpins Key Pathogenicity Pathways.

Cryptococcus gattii is a pathogenic yeast of humans and other animals which causes disease predominantly in immunocompetent hosts. Infection begins when aerosolized yeast or spores enter the body, triggering an immune response, including engulfment by macrophages. To understand the early transcriptional signals in both the yeast and its mammalian host, we performed a time-course dual-transcriptome sequencing (RNA-seq) experiment for four lineages of C. gattii (lineages VGI to IV) interacting with mouse macrophages at 1, 3, and 6 h postinfection. Comparisons of in vitro to ex vivo gene expression levels indicated that lineage VGII is transcriptionally divergent from non-VGII lineages, including differential expression of genes involved in capsule synthesis, capsule attachment, and ergosterol production. Several paralogous genes demonstrated subfunctionalization between lineages, including upregulation of capsule biosynthesis-related gene CAP2 and downregulation of CAP1 in VGIII. Isolates also compensate for lineage-specific gene losses by overexpression of genetically similar paralogs, including overexpression of capsule gene CAS3 in VGIV, which have lost the CAS31 gene. Differential expression of one in five C. gattii genes was detected following coincubation with mouse macrophages; all isolates showed high induction of oxidative-reduction functions and downregulation of capsule attachment genes. We also found that VGII switches expression of two laccase paralogs (from LAC1 to LAC2) during coincubation of macrophages. Finally, we found that mouse macrophages respond to all four lineages of C. gattii by upregulating FosB/Jun/Egr1 regulatory proteins at early time points. This report highlights the evolutionary breadth of expression profiles among the lineages of C. gattii and the diversity of transcriptional responses at this host-pathogen interface.IMPORTANCE The transcriptional profiles of related pathogens and their responses to host-induced stresses underpin their pathogenicity. Expression differences between related pathogens during host interaction can indicate when and how these genes contribute to virulence, ultimately informing new and improved treatment strategies for those diseases. In this paper, we compare the transcriptional profiles of five isolates representing four lineages of C. gattii in rich media. Our analyses identified key processes, including those involving cell capsule, ergosterol production, and melanin, that are differentially expressed between lineages, and we found that VGII has the most distinct profile in terms of numbers of differentially expressed genes. All lineages have also undergone subfunctionalization for several paralogs, including capsule biosynthesis and attachment genes. Most genes appeared downregulated during coincubation with macrophages, with the largest decrease observed for capsule attachment genes, which appeared to be coordinated with a stress response, as all lineages also upregulated oxidative stress response genes. Furthermore, VGII upregulated many genes that are linked to ergosterol biosynthesis and switched from expression of the laccase LAC1 to expression of LAC2 ex vivo Finally, we saw a pronounced increase in the FosB/Jun/Egr1 regulatory proteins at early time points in bone marrow-derived macrophages, marking a role in the host response to C. gattii This work highlights the dynamic roles of key C. gattii virulence genes in response to macrophages.

Global analysis of mutations driving microevolution of a heterozygous diploid fungal pathogen.

Candida albicans is a heterozygous diploid yeast that is a commensal of the human gastrointestinal tract and a prevalent opportunistic pathogen. Here, whole-genome sequencing was performed on multiple C. albicans isolates passaged both in vitro and in vivo to characterize the complete spectrum of mutations arising in laboratory culture and in the mammalian host. We establish that, independent of culture niche, microevolution is primarily driven by de novo base substitutions and frequent short-tract loss-of-heterozygosity events. An average base-substitution rate of ∼1.2 × 10-10 per base pair per generation was observed in vitro, with higher rates inferred during host infection. Large-scale chromosomal changes were relatively rare, although chromosome 7 trisomies frequently emerged during passaging in a gastrointestinal model and was associated with increased fitness for this niche. Multiple chromosomal features impacted mutational patterns, with mutation rates elevated in repetitive regions, subtelomeric regions, and in gene families encoding cell surface proteins involved in host adhesion. Strikingly, de novo mutation rates were more than 800-fold higher in regions immediately adjacent to emergent loss-of-heterozygosity tracts, indicative of recombination-induced mutagenesis. Furthermore, genomes showed biased patterns of mutations suggestive of extensive purifying selection during passaging. These results reveal how both cell-intrinsic and cell-extrinsic factors influence C. albicans microevolution, and provide a quantitative picture of genome dynamics in this heterozygous diploid species.