Search

ABSTRACT 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.

ABSTRACT 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.

ABSTRACT Histoplasma capsulatum, a dimorphic fungal pathogen, is the most common cause of fungal respiratory infections in immunocompetent hosts. Histoplasma is endemic in the Ohio and Mississippi River Valleys in the United States and is also distributed worldwide. Previous studies have revealed at least eight clades, each specific to a geographic location: North American classes 1 and 2 (NAm 1 and NAm 2), Latin American groups A and B (LAm A and LAm B), Eurasian, Netherlands, Australian and African, and an additional distinct lineage (H81) comprised of Panamanian isolates. Previously assembled Histoplasma genomes are highly fragmented, with the highly repetitive G217B (NAm 2) strain, which has been used for most whole-genome-scale transcriptome studies, assembled into over 250 contigs. In this study, we set out to fully assemble the repeat regions and characterize the large-scale genome architecture of Histoplasma species. We resequenced five Histoplasma strains (WU24 [NAm 1], G217B [NAm 2], H88 [African], G186AR [Panama], and G184AR [Panama]) using Oxford Nanopore Technologies long-read sequencing technology. Here, we report chromosomal-level assemblies for all five strains, which exhibit extensive synteny among the geographically distant Histoplasma isolates. The new assemblies revealed that RYP2, a major regulator of morphology and virulence, is duplicated in G186AR. In addition, we mapped previously generated transcriptome data sets onto the newly assembled chromosomes. Our analyses revealed that the expression of transposons and transposon-embedded genes are upregulated in yeast phase compared to mycelial phase in the G217B and H88 strains. This study provides an important resource for fungal researchers and further highlights the importance of chromosomal-level assemblies in analyzing high-throughput data sets. IMPORTANCE Histoplasma species are dimorphic fungi causing significant morbidity and mortality worldwide. These fungi grow as mold in the soil and as budding yeast within the human host. Histoplasma can be isolated from soil in diverse regions, including North America, South America, Africa, and Europe. Phylogenetically distinct species of Histoplasma have been isolated and sequenced. However, for the commonly used strains, genome assemblies have been fragmented, leading to underutilization of genome-scale data. This study provides chromosome-level assemblies of the commonly used Histoplasma strains using long-read sequencing technology. Comparative analysis of these genomes shows largely conserved gene order within the chromosomes. Mapping existing transcriptome data on these new assemblies reveals clustering of transcriptionally coregulated genes. The results of this study highlight the importance of obtaining chromosome-level assemblies in understanding the biology of human fungal pathogens.

ABSTRACT Cryptococcus neoformans is a major human central nervous system (CNS) fungal pathogen causing considerable morbidity and mortality. In this study, we provide the widest view to date of the yeast transcriptome directly from the human subarachnoid space and within cerebrospinal fluid (CSF). We captured yeast transcriptomes from C. neoformans of various genotypes in 31 patients with cryptococcal meningoencephalitis as well as several Cryptococcus gattii infections. Using transcriptome sequencing (RNA-seq) analyses, we compared the in vivo yeast transcriptomes to those from other environmental conditions, including in vitro growth on nutritious media or artificial CSF as well as samples collected from rabbit CSF at two time points. We ranked gene expressions and identified genetic patterns and networks across these diverse isolates that reveal an emphasis on carbon metabolism, fatty acid synthesis, transport, cell wall structure, and stress-related gene functions during growth in CSF. The most highly expressed yeast genes in human CSF included those known to be associated with survival or virulence and highlighted several genes encoding hypothetical proteins. From that group, a gene encoding the CMP1 putative glycoprotein (CNAG_06000) was selected for functional studies. This gene was found to impact the virulence of Cryptococcus in both mice and the CNS rabbit model, in agreement with a recent study also showing a role in virulence. This transcriptional analysis strategy provides a view of regulated yeast genes across genetic backgrounds important for human CNS infection and a relevant resource for the study of cryptococcal genes, pathways, and networks linked to human disease. IMPORTANCE Cryptococcus is the most common fungus causing high-morbidity and -mortality human meningitis. This encapsulated yeast has a unique propensity to travel to the central nervous system to produce disease. In this study, we captured transcriptomes of yeasts directly out of the human cerebrospinal fluid, the most concerning site of infection. By comparing the RNA transcript levels with other conditions, we gained insights into how the basic machinery involved in metabolism and environmental responses enable this fungus to cause disease at this body site. This approach was applied to clinical isolates with diverse genotypes to begin to establish a genotype-agnostic understanding of how the yeast responds to stress. Based on these results, future studies can focus on how these genes and their pathways and networks can be targeted with new therapeutics and possibly classify yeasts with bad infection outcomes.

Candida aurisC. auris

Cryptococcus neoformansC. neoformans

ABSTRACT Candida auris has emerged as a multidrug-resistant pathogen of great clinical concern. Approximately 90% of clinical C. auris isolates are resistant to fluconazole, the most commonly prescribed antifungal agent, and yet it remains unknown what mechanisms underpin this fluconazole resistance. To identify novel mechanisms contributing to fluconazole resistance in C. auris, fluconazole-susceptible C. auris clinical isolate AR0387 was passaged in media supplemented with fluconazole to generate derivative strains which had acquired increased fluconazole resistance in vitro. Comparative analyses of comprehensive sterol profiles, [3H]fluconazole uptake, sequencing of C. auris genes homologous to genes known to contribute to fluconazole resistance in other species of Candida, and relative expression levels of C. auris ERG11, CDR1, and MDR1 were performed. All fluconazole-evolved derivative strains were found to have acquired mutations in the zinc-cluster transcription factor-encoding gene TAC1B and to show a corresponding increase in CDR1 expression relative to the parental clinical isolate, AR0387. Mutations in TAC1B were also identified in a set of 304 globally distributed C. auris clinical isolates representing each of the four major clades. Introduction of the most common mutation found among fluconazole-resistant clinical isolates of C. auris into fluconazole-susceptible isolate AR0387 was confirmed to increase fluconazole resistance by 8-fold, and the correction of the same mutation in a fluconazole-resistant isolate, AR0390, decreased fluconazole MIC by 16-fold. Taken together, these data demonstrate that C. auris can rapidly acquire resistance to fluconazole in vitro and that mutations in TAC1B significantly contribute to clinical fluconazole resistance. IMPORTANCE Candida auris is an emerging multidrug-resistant pathogen of global concern, known to be responsible for outbreaks on six continents and to be commonly resistant to antifungals. While the vast majority of clinical C. auris isolates are highly resistant to fluconazole, an essential part of the available antifungal arsenal, very little is known about the mechanisms contributing to resistance. In this work, we show that mutations in the transcription factor TAC1B significantly contribute to clinical fluconazole resistance. These studies demonstrated that mutations in TAC1B can arise rapidly in vitro upon exposure to fluconazole and that a multitude of resistance-associated TAC1B mutations are present among the majority of fluconazole-resistant C. auris isolates from a global collection and appear specific to a subset of lineages or clades. Thus, identification of this novel genetic determinant of resistance significantly adds to the understanding of clinical antifungal resistance in C. auris.

ABSTRACT 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.

ABSTRACT 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. IMPORTANCE Pneumocystis 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.

ABSTRACT 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. auris. IMPORTANCE 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.

ABSTRACT 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).

ABSTRACT 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. IMPORTANCE Cryptococcus 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.

Cryptococcus species are life-threatening human fungal pathogens that cause cryptococcal meningoencephalitis in both immunocompromised and healthy hosts. The natural environmental niches of Cryptococcus include pigeon (Columba livia) guano, soil, and a variety of tree species such as Eucalyptus camaldulensis, Ceratonia siliqua, Platanus orientalis, and Pinus spp. Genetic and genomic studies of extensive sample collections have provided insights into the population distribution and composition of different Cryptococcus species in geographic regions around the world. However, few such studies examined Cryptococcus in Turkey. We sampled 388 Olea europaea (olive) and 132 E. camaldulensis trees from seven locations in coastal and inland areas of the Aegean region of Anatolian Turkey in September 2016 to investigate the distribution and genetic diversity present in the natural Cryptococcus population. We isolated 84 Cryptococcus neoformans strains (83 MATα and 1 MATa) and 3 Cryptococcus deneoformans strains (all MATα) from 87 (22.4% of surveyed) O. europaea trees; a total of 32 C. neoformans strains were isolated from 32 (24.2%) of the E. camaldulensis trees, all of which were MATα. A statistically significant difference was observed in the frequency of C. neoformans isolation between coastal and inland areas (P < 0.05). Interestingly, the MATaC. neoformans isolate was fertile in laboratory crosses with VNI and VNB MATα tester strains and produced robust hyphae, basidia, and basidiospores, thus suggesting potential sexual reproduction in the natural population. Sequencing analyses of the URA5 gene identified at least five different genotypes among the isolates. Population genetics and genomic analyses revealed that most of the isolates in Turkey belong to the VNBII lineage of C. neoformans, which is predominantly found in southern Africa; these isolates are part of a distinct minor clade within VNBII that includes several isolates from Zambia and Brazil. Our study provides insights into the geographic distribution of different C. neoformans lineages in the Mediterranean region and highlights the need for wider geographic sampling to gain a better understanding of the natural habitats, migration, epidemiology, and evolution of this important human fungal pathogen.

ABSTRACT Speciation is a central mechanism of biological diversification. While speciation is well studied in plants and animals, in comparison, relatively little is known about speciation in fungi. One fungal model is the Cryptococcus genus, which is best known for the pathogenic Cryptococcus neoformans/Cryptococcus gattii species complex that causes >200,000 new human infections annually. Elucidation of how these species evolved into important human-pathogenic species remains challenging and can be advanced by studying the most closely related nonpathogenic species, Cryptococcus amylolentus and Tsuchiyaea wingfieldii. However, these species have only four known isolates, and available data were insufficient to determine species boundaries within this group. By analyzing full-length chromosome assemblies, we reappraised the phylogenetic relationships of the four available strains, confirmed the genetic separation of C. amylolentus and T. wingfieldii (now Cryptococcus wingfieldii), and revealed an additional cryptic species, for which the name Cryptococcus floricola is proposed. The genomes of the three species are ∼6% divergent and exhibit significant chromosomal rearrangements, including inversions and a reciprocal translocation that involved intercentromeric ectopic recombination, which together likely impose significant barriers to genetic exchange. Using genetic crosses, we show that while C. wingfieldii cannot interbreed with any of the other strains, C. floricola can still undergo sexual reproduction with C. amylolentus. However, most of the resulting spores were inviable or sterile or showed reduced recombination during meiosis, indicating that intrinsic postzygotic barriers had been established. Our study and genomic data will foster additional studies addressing fungal speciation and transitions between nonpathogenic and pathogenic Cryptococcus lineages. IMPORTANCE The evolutionary drivers of speciation are critical to our understanding of how new pathogens arise from nonpathogenic lineages and adapt to new environments. Here we focus on the Cryptococcus amylolentus species complex, a nonpathogenic fungal lineage closely related to the human-pathogenic Cryptococcus neoformans/Cryptococcus gattii complex. Using genetic and genomic analyses, we reexamined the species boundaries of four available isolates within the C. amylolentus complex and revealed three genetically isolated species. Their genomes are ∼6% divergent and exhibit chromosome rearrangements, including translocations and small-scale inversions. Although two of the species (C. amylolentus and newly described C. floricola) were still able to interbreed, the resulting hybrid progeny were usually inviable or sterile, indicating that barriers to reproduction had already been established. These results advance our understanding of speciation in fungi and highlight the power of genomics in assisting our ability to correctly identify and discriminate fungal species.

ABSTRACT Candida auris is an emerging fungal pathogen and a serious global health threat as the majority of clinical isolates display elevated resistance to currently available antifungal drugs. Despite the increased prevalence of C. auris infections, the mechanisms governing drug resistance remain largely elusive. In diverse fungi, the evolution of drug resistance is enabled by the essential molecular chaperone Hsp90, which stabilizes key regulators of cellular responses to drug-induced stress. Hsp90 also orchestrates temperature-dependent morphogenesis in Candida albicans, a key virulence trait. However, the role of Hsp90 in the pathobiology of C. auris remains unknown. In order to study regulatory functions of Hsp90 in C. auris, we placed HSP90 under the control of a doxycycline-repressible promoter to enable transcriptional repression. We found that Hsp90 is essential for growth in C. auris and that it enables tolerance of clinical isolates with respect to the azoles, which inhibit biosynthesis of the membrane sterol ergosterol. High-level azole resistance was independent of Hsp90 but dependent on the ABC transporter CDR1, deletion of which resulted in abrogated resistance. Strikingly, we discovered that C. auris undergoes a morphogenetic transition from yeast to filamentous growth in response to HSP90 depletion or cell cycle arrest but not in response to other cues that induce C. albicans filamentation. Finally, we observed that this developmental transition is associated with global transcriptional changes, including the induction of cell wall-related genes. Overall, this report provides a novel insight into mechanisms of drug tolerance and resistance in C. auris and describes a developmental transition in response to perturbation of a core regulator of protein homeostasis. IMPORTANCE Fungal pathogens pose a serious threat to public health. Candida auris is an emerging fungal pathogen that is often resistant to commonly used antifungal drugs. However, the mechanisms governing drug resistance and virulence in this organism remain largely unexplored. In this study, we adapted a conditional expression system to modulate the transcription of an essential gene, HSP90, which regulates antifungal resistance and virulence in diverse fungal pathogens. We showed that Hsp90 is essential for growth in C. auris and is important for tolerance of the clinically important azole antifungals, which block ergosterol biosynthesis. Further, we established that the Cdr1 efflux transporter regulates azole resistance. Finally, we discovered that C. auris transitions from yeast to filamentous growth in response to Hsp90 inhibition, accompanied by global transcriptional remodeling. Overall, this work provides a novel insight into mechanisms regulating azole resistance in C. auris and uncovers a distinct developmental program regulated by Hsp90.

ABSTRACT Pathogenic species of Cryptococcus cause hundreds of thousands of deaths annually. Considerable phenotypic variation is exhibited during infection, including increased capsule size, capsule shedding, giant cells (≥15 μm), and micro cells (≤1 μm). We examined 70 clinical isolates of Cryptococcus neoformans and Cryptococcus tetragattii from HIV/AIDS patients in Botswana to determine whether the capacity to produce morphological variants was associated with clinical parameters. Isolates were cultured under conditions designed to simulate in vivo stresses. Substantial variation was seen across morphological and clinical data. Giant cells were more common in C. tetragattii, while micro cells and shed capsule occurred in C. neoformans only. Phenotypic variables fell into two groups associated with differing symptoms. The production of “large” phenotypes (greater cell and capsule size and giant cells) was associated with higher CD4 count and was negatively correlated with intracranial pressure indicators, suggesting that these are induced in early stage infection. “Small” phenotypes (micro cells and shed capsule) were associated with lower CD4 counts, negatively correlated with meningeal inflammation indicators, and positively correlated with intracranial pressure indicators, suggesting that they are produced later during infection and may contribute to immune suppression and promote proliferation and dissemination. These trends persisted at the species level, indicating that they were not driven by association with particular Cryptococcus species. Isolates possessing giant cells, micro cells, and shed capsule were rare, but strikingly, they were associated with patient death (P = 0.0165). Our data indicate that pleomorphism is an important driver in Cryptococcus infection. IMPORTANCE Cryptococcosis results in hundreds of thousands of deaths annually, predominantly in sub-Saharan Africa. Cryptococcus is an encapsulated yeast, and during infection, cells have the capacity for substantial morphological changes, including capsule enlargement and shedding and variations in cell shape and size. In this study, we examined 70 Cryptococcus isolates causing meningitis in HIV/AIDS patients in Botswana in order to look for associations between phenotypic variation and clinical symptoms. Four variant phenotypes were seen across strains: giant cells of ≥15 µm, micro cells of ≤1 µm, shed extracellular capsule, and irregularly shaped cells. We found that “large” and “small” phenotypes were associated with differing disease symptoms, indicating that their production may be important during the disease process. Overall, our study indicates that Cryptococcus strains that can switch on cell types under different situations may be more able to sustain infection and resist the host response.

ABSTRACT 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.

ABSTRACT Rhizopus delemar is an invasive fungal pathogen responsible for the frequently fatal disease mucormycosis. Germination, a crucial mechanism by which infectious spores of Rhizopus delemar cause disease, is a key developmental process that transforms the dormant spore state into a vegetative one. The molecular mechanisms that underpin this transformation may be key to controlling mucormycosis; however, the regulation of germination remains poorly understood. This study describes the phenotypic and transcriptional changes that take place over the course of germination. This process is characterized by four distinct stages: dormancy, isotropic swelling, germ tube emergence, and hyphal growth. Dormant spores are shown to be transcriptionally unique, expressing a subset of transcripts absent in later developmental stages. A large shift in the expression profile is prompted by the initiation of germination, with genes involved in respiration, chitin, cytoskeleton, and actin regulation appearing to be important for this transition. A period of transcriptional consistency can be seen throughout isotropic swelling, before the transcriptional landscape shifts again at the onset of hyphal growth. This study provides a greater understanding of the regulation of germination and highlights processes involved in transforming Rhizopus delemar from a single-cellular to multicellular organism. IMPORTANCE Germination is key to the growth of many organisms, including fungal spores. Mucormycete spores exist abundantly within the environment and germinate to form hyphae. These spores are capable of infecting immunocompromised individuals, causing the disease mucormycosis. Germination from spore to hyphae within patients leads to angioinvasion, tissue necrosis, and often fatal infections. This study advances our understanding of how spore germination occurs in the mucormycetes, identifying processes we may be able to inhibit to help prevent or treat mucormycosis.

ABSTRACT Pneumocystis species are opportunistic mammalian pathogens that cause severe pneumonia in immunocompromised individuals. These fungi are highly host specific and uncultivable in vitro. Human Pneumocystis infections present major challenges because of a limited therapeutic arsenal and the rise of drug resistance. To investigate the diversity and demographic history of natural populations of Pneumocystis infecting humans, rats, and mice, we performed whole-genome and large-scale multilocus sequencing of infected tissues collected in various geographic locations. Here, we detected reduced levels of recombination and variations in historical demography, which shape the global population structures. We report estimates of evolutionary rates, levels of genetic diversity, and population sizes. Molecular clock estimates indicate that Pneumocystis species diverged before their hosts, while the asynchronous timing of population declines suggests host shifts. Our results have uncovered complex patterns of genetic variation influenced by multiple factors that shaped the adaptation of Pneumocystis populations during their spread across mammals. IMPORTANCE Understanding how natural pathogen populations evolve and identifying the determinants of genetic variation are central issues in evolutionary biology. Pneumocystis, a fungal pathogen which infects mammals exclusively, provides opportunities to explore these issues. In humans, Pneumocystis can cause a life-threatening pneumonia in immunosuppressed individuals. In analysis of different Pneumocystis species infecting humans, rats, and mice, we found that there are high infection rates and that natural populations maintain a high level of genetic variation despite low levels of recombination. We found no evidence of population structuring by geography. Our comparisons of the times of divergence of these species to their respective hosts suggest that Pneumocystis may have undergone recent host shifts. The results demonstrate that Pneumocystis strains are widely disseminated geographically and provide a new understanding of the evolution of these pathogens.

Cryptococcus species are the causative agents of cryptococcal meningitis, a significant source of mortality in immunocompromised individuals. Initial work on the molecular epidemiology of this fungal pathogen utilized genotyping approaches to describe the genetic diversity and biogeography of two species, Cryptococcus neoformans and Cryptococcus gattii. Whole genome sequencing of representatives of both species resulted in reference assemblies enabling a wide array of downstream studies and genomic resources. With the increasing availability of whole genome sequencing, both species have now had hundreds of individual isolates sequenced, providing fine-scale insight into the evolution and diversification of Cryptococcus and allowing for the first genome-wide association studies to identify genetic variants associated with human virulence. Sequencing has also begun to examine the microevolution of isolates during prolonged infection and to identify variants specific to outbreak lineages, highlighting the potential role of hyper-mutation in evolving within short time scales. We can anticipate that further advances in sequencing technology and sequencing microbial genomes at scale, including metagenomics approaches, will continue to refine our view of how the evolution of Cryptococcus drives its success as a pathogen.

ABSTRACT The pathogenic species of Cryptococcus are a major cause of mortality owing to severe infections in immunocompromised as well as immunocompetent individuals. Although antifungal treatment is usually effective, many patients relapse after treatment, and in such cases, comparative analyses of the genomes of incident and relapse isolates may reveal evidence of determinative, microevolutionary changes within the host. Here, we analyzed serial isolates cultured from cerebrospinal fluid specimens of 18 South African patients with recurrent cryptococcal meningitis. The time between collection of the incident isolates and collection of the relapse isolates ranged from 124 days to 290 days, and the analyses revealed that, during this period within the patients, the isolates underwent several genetic and phenotypic changes. Considering the vast genetic diversity of cryptococcal isolates in sub-Saharan Africa, it was not surprising to find that the relapse isolates had acquired different genetic and correlative phenotypic changes. They exhibited various mechanisms for enhancing virulence, such as growth at 39°C, adaptation to stress, and capsule production; a remarkable amplification of ERG11 at the native and unlinked locus may provide stable resistance to fluconazole. Our data provide a deeper understanding of the microevolution of Cryptococcus species under pressure from antifungal chemotherapy and host immune responses. This investigation clearly suggests a promising strategy to identify novel targets for improved diagnosis, therapy, and prognosis. IMPORTANCE Opportunistic infections caused by species of the pathogenic yeast Cryptococcus lead to chronic meningoencephalitis and continue to ravage thousands of patients with HIV/AIDS. Despite receiving antifungal treatment, over 10% of patients develop recurrent disease. In this study, we collected isolates of Cryptococcus from cerebrospinal fluid specimens of 18 patients at the time of their diagnosis and when they relapsed several months later. We then sequenced and compared the genomic DNAs of each pair of initial and relapse isolates. We also tested the isolates for several key properties related to cryptococcal virulence as well as for their susceptibility to the antifungal drug fluconazole. These analyses revealed that the relapsing isolates manifested multiple genetic and chromosomal changes that affected a variety of genes implicated in the pathogenicity of Cryptococcus or resistance to fluconazole. This application of comparative genomics to serial clinical isolates provides a blueprint for identifying the mechanisms whereby pathogenic microbes adapt within patients to prolong disease.

ABSTRACT Cryptococcosis is a potentially lethal disease of humans/animals caused by Cryptococcus neoformans and Cryptococcus gattii. Distinction between the two species is based on phenotypic and genotypic characteristics. Recently, it was proposed that C. neoformans be divided into two species and C. gattii into five species based on a phylogenetic analysis of 115 isolates. While this proposal adds to the knowledge about the genetic diversity and population structure of cryptococcosis agents, the published genotypes of 2,606 strains have already revealed more genetic diversity than is encompassed by seven species. Naming every clade as a separate species at this juncture will lead to continuing nomenclatural instability. In the absence of biological differences between clades and no consensus about how DNA sequence alone can delineate a species, we recommend using “Cryptococcus neoformans species complex” and “C. gattii species complex” as a practical intermediate step, rather than creating more species. This strategy recognizes genetic diversity without creating confusion.

ABSTRACT The Paracoccidioides genus includes two species of thermally dimorphic fungi that cause paracoccidioidomycosis, a neglected health-threatening human systemic mycosis endemic to Latin America. To examine the genome evolution and the diversity of Paracoccidioides spp., we conducted whole-genome sequencing of 31 isolates representing the phylogenetic, geographic, and ecological breadth of the genus. These samples included clinical, environmental and laboratory reference strains of the S1, PS2, PS3, and PS4 lineages of P. brasiliensis and also isolates of Paracoccidioides lutzii species. We completed the first annotated genome assemblies for the PS3 and PS4 lineages and found that gene order was highly conserved across the major lineages, with only a few chromosomal rearrangements. Comparing whole-genome assemblies of the major lineages with single-nucleotide polymorphisms (SNPs) predicted from the remaining 26 isolates, we identified a deep split of the S1 lineage into two clades we named S1a and S1b. We found evidence for greater genetic exchange between the S1b lineage and all other lineages; this may reflect the broad geographic range of S1b, which is often sympatric with the remaining, largely geographically isolated lineages. In addition, we found evidence of positive selection for the GP43 and PGA1 antigen genes and genes coding for other secreted proteins and proteases and lineage-specific loss-of-function mutations in cell wall and protease genes; these together may contribute to virulence and host immune response variation among natural isolates of Paracoccidioides spp. These insights into the recent evolutionary events highlight important differences between the lineages that could impact the distribution, pathogenicity, and ecology of Paracoccidioides. IMPORTANCE Characterization of genetic differences between lineages of the dimorphic human-pathogenic fungus Paracoccidioides can identify changes linked to important phenotypes and guide the development of new diagnostics and treatments. In this article, we compared genomes of 31 diverse isolates representing the major lineages of Paracoccidioides spp. and completed the first annotated genome sequences for the PS3 and PS4 lineages. We analyzed the population structure and characterized the genetic diversity among the lineages of Paracoccidioides, including a deep split of S1 into two lineages (S1a and S1b), and differentiated S1b, associated with most clinical cases, as the more highly recombining and diverse lineage. In addition, we found patterns of positive selection in surface proteins and secreted enzymes among the lineages, suggesting diversifying mechanisms of pathogenicity and adaptation across this species complex. These genetic differences suggest associations with the geographic range, pathogenicity, and ecological niches of Paracoccidioides lineages.

ABSTRACT Cryptococcus gattii is a fungal pathogen of humans, causing pulmonary infections in otherwise healthy hosts. To characterize genomic variation among the four major lineages of C. gattii (VGI, -II, -III, and -IV), we generated, annotated, and compared 16 de novo genome assemblies, including the first for the rarely isolated lineages VGIII and VGIV. By identifying syntenic regions across assemblies, we found 15 structural rearrangements, which were almost exclusive to the VGI-III-IV lineages. Using synteny to inform orthology prediction, we identified a core set of 87% of C. gattii genes present as single copies in all four lineages. Remarkably, 737 genes are variably inherited across lineages and are overrepresented for response to oxidative stress, mitochondrial import, and metal binding and transport. Specifically, VGI has an expanded set of iron-binding genes thought to be important to the virulence of Cryptococcus, while VGII has expansions in the stress-related heat shock proteins relative to the other lineages. We also characterized genes uniquely absent in each lineage, including a copper transporter absent from VGIV, which influences Cryptococcus survival during pulmonary infection and the onset of meningoencephalitis. Through inclusion of population-level data for an additional 37 isolates, we identified a new transcontinental clonal group that we name VGIIx, mitochondrial recombination between VGII and VGIII, and positive selection of multidrug transporters and the iron-sulfur protein aconitase along multiple branches of the phylogenetic tree. Our results suggest that gene expansion or contraction and positive selection have introduced substantial variation with links to mechanisms of pathogenicity across this species complex. IMPORTANCE The genetic differences between phenotypically different pathogens provide clues to the underlying mechanisms of those traits and can lead to new drug targets and improved treatments for those diseases. In this paper, we compare 16 genomes belonging to four highly differentiated lineages of Cryptococcus gattii, which cause pulmonary infections in otherwise healthy humans and other animals. Half of these lineages have not had their genomes previously assembled and annotated. We identified 15 ancestral rearrangements in the genome and over 700 genes that are unique to one or more lineages, many of which are associated with virulence. In addition, we found evidence for recent transcontinental spread, mitochondrial genetic exchange, and positive selection in multidrug transporters. Our results suggest that gene expansion/contraction and positive selection are diversifying the mechanisms of pathogenicity across this species complex.

ABSTRACT Food-borne pathogens are ongoing problems, and new pathogens are emerging. The impact of fungi, however, is largely underestimated. Recently, commercial yogurts contaminated with Mucor circinelloides were sold, and >200 consumers became ill with nausea, vomiting, and diarrhea. Mucoralean fungi cause the fatal fungal infection mucormycosis, whose incidence has been continuously increasing. In this study, we isolated an M. circinelloides strain from a yogurt container, and multilocus sequence typing identified the strain as Mucor circinelloides f. circinelloides. M. circinelloides f. circinelloides is the most virulent M. circinelloides subspecies and is commonly associated with human infections, whereas M. circinelloides f. lusitanicus and M. circinelloides f. griseocyanus are less common causes of infection. Whole-genome analysis of the yogurt isolate confirmed it as being close to the M. circinelloides f. circinelloides subgroup, with a higher percentage of divergence with the M. circinelloides f. lusitanicus subgroup. In mating assays, the yogurt isolate formed sexual zygospores with the (−) M. circinelloides f. circinelloides tester strain, which is congruent with its sex locus encoding SexP, the (+) mating type sex determinant. The yogurt isolate was virulent in murine and wax moth larva host systems. In a murine gastromucormycosis model, Mucor was recovered from fecal samples of infected mice for up to 10 days, indicating that Mucor can survive transit through the GI tract. In interactions with human immune cells, M. circinelloides f. lusitanicus induced proinflammatory cytokines but M. circinelloides f. circinelloides did not, which may explain the different levels of virulence in mammalian hosts. This study demonstrates that M. circinelloides can spoil food products and cause gastrointestinal illness in consumers and may pose a particular risk to immunocompromised patients. IMPORTANCE The U.S. FDA reported that yogurt products were contaminated with M. circinelloides, a mucoralean fungal pathogen, and >200 consumers complained of symptoms, including vomiting, nausea, and diarrhea. The manufacturer voluntarily withdrew the affected yogurt products from the market. Compared to other food-borne pathogens, including bacteria, viruses, and parasites, less focus has been placed on the risk of fungal pathogens. This study evaluates the potential risk from the food-borne fungal pathogen M. circinelloides that was isolated from the contaminated commercial yogurt. We successfully cultured an M. circinelloides isolate and found that the isolate belongs to the species M. circinelloides f. circinelloides, which is often associated with human infections. In murine and insect host models, the isolate was virulent. While information disseminated in the popular press would suggest this fungal contaminant poses little or no risk to consumers, our results show instead that it is capable of causing significant infections in animals.

ABSTRACT An outbreak of the fungal pathogen Cryptococcus gattii began in the Pacific Northwest (PNW) in the late 1990s. This outbreak consists of three clonal subpopulations: VGIIa/major, VGIIb/minor, and VGIIc/novel. Both VGIIa and VGIIc are unique to the PNW and exhibit increased virulence. In this study, we sequenced the genomes of isolates from these three groups, as well as global isolates, and analyzed a total of 53 isolates. We found that VGIIa/b/c populations show evidence of clonal expansion in the PNW. Whole-genome sequencing provided evidence that VGIIb originated in Australia, while VGIIa may have originated in South America, and these were likely independently introduced. Additionally, the VGIIa outbreak lineage may have arisen from a less virulent clade that contained a mutation in the MSH2 ortholog, but this appears to have reverted in the VGIIa outbreak strains, suggesting that a transient mutator phenotype may have contributed to adaptation and evolution of virulence in the PNW outbreak. PNW outbreak isolates share genomic islands, both between the clonal lineages and with global isolates, indicative of sexual recombination. This suggests that VGII C. gattii has undergone sexual reproduction, either bisexual or unisexual, in multiple locales contributing to the production of novel, virulent subtypes. We also found that the genomes of two basal VGII isolates from HIV+ patients contain an introgression tract spanning three genes. Introgression substantially contributed to intra-VGII polymorphism and likely occurred through sexual reproduction with VGI. More broadly, these findings illustrate how both microevolution and sexual reproduction play central roles in the development of infectious outbreaks from avirulent or less virulent progenitors. IMPORTANCE Cryptococcus gattii is the causative agent responsible for ongoing infections in the Pacific Northwest of the United States and western Canada. The incidence of these infections increased dramatically in the 1990s and remains elevated. These infections are attributable to three clonal lineages of C. gattii, VGIIa, VGIIb, and VGIIc, with only VGIIa identified once previously in the Pacific Northwest prior to the start of the outbreak, albeit in a less virulent form. This study addresses the origin and emergence of this outbreak, using whole-genome sequencing and comparison of both outbreak and global isolates. We show that VGIIa arose mitotically from a less virulent clonal group, possibly via the action of a mutator phenotype, while VGIIb was likely introduced from Australia, and VGIIc appears to have emerged in the United States or in an undersampled locale via sexual reproduction. This work shows that multiple processes can contribute to the emergence of an outbreak.

ABSTRACT The major cause of athlete’s foot is Trichophyton rubrum, a dermatophyte or fungal pathogen of human skin. To facilitate molecular analyses of the dermatophytes, we sequenced T. rubrum and four related species, Trichophyton tonsurans, Trichophyton equinum, Microsporum canis, and Microsporum gypseum. These species differ in host range, mating, and disease progression. The dermatophyte genomes are highly colinear yet contain gene family expansions not found in other human-associated fungi. Dermatophyte genomes are enriched for gene families containing the LysM domain, which binds chitin and potentially related carbohydrates. These LysM domains differ in sequence from those in other species in regions of the peptide that could affect substrate binding. The dermatophytes also encode novel sets of fungus-specific kinases with unknown specificity, including nonfunctional pseudokinases, which may inhibit phosphorylation by competing for kinase sites within substrates, acting as allosteric effectors, or acting as scaffolds for signaling. The dermatophytes are also enriched for a large number of enzymes that synthesize secondary metabolites, including dermatophyte-specific genes that could synthesize novel compounds. Finally, dermatophytes are enriched in several classes of proteases that are necessary for fungal growth and nutrient acquisition on keratinized tissues. Despite differences in mating ability, genes involved in mating and meiosis are conserved across species, suggesting the possibility of cryptic mating in species where it has not been previously detected. These genome analyses identify gene families that are important to our understanding of how dermatophytes cause chronic infections, how they interact with epithelial cells, and how they respond to the host immune response. IMPORTANCE Athlete’s foot, jock itch, ringworm, and nail infections are common fungal infections, all caused by fungi known as dermatophytes (fungi that infect skin). This report presents the genome sequences of Trichophyton rubrum, the most frequent cause of athlete’s foot, as well as four other common dermatophytes. Dermatophyte genomes are enriched for four gene classes that may contribute to the ability of these fungi to cause disease. These include (i) proteases secreted to degrade skin; (ii) kinases, including pseudokinases, that are involved in signaling necessary for adapting to skin; (iii) secondary metabolites, compounds that act as toxins or signals in the interactions between fungus and host; and (iv) a class of proteins (LysM) that appear to bind and mask cell wall components and carbohydrates, thus avoiding the host’s immune response to the fungi. These genome sequences provide a strong foundation for future work in understanding how dermatophytes cause disease.