Evolutionary Persistence of DNA Methylation for Millions of Years after Ancient Loss of a De Novo Methyltransferase

Cytosine methylation of DNA is a widespread modification of DNA that plays numerous critical roles. In the yeast Cryptococcus neoformans, CG methylation occurs in transposon-rich repeats and requires the DNA methyltransferase Dnmt5. We show that Dnmt5 displays exquisite maintenance-type specificity in vitro and in vivo and utilizes similar in vivo cofactors as the metazoan maintenance methylase Dnmt1. Remarkably, phylogenetic and functional analysis revealed that the ancestral species lost the gene for a de novo methylase, DnmtX, between 50-150 mya. We examined how methylation has persisted since the ancient loss of DnmtX. Experimental and comparative studies reveal efficient replication of methylation patterns in C. neoformans, rare stochastic methylation loss and gain events, and the action of natural selection. We propose that an epigenome has been propagated for >50 million years through a process analogous to Darwinian evolution of the genome.

Complete Genome Sequences for Two Talaromyces marneffei Clinical Isolates from Northern and Southern Vietnam

Talaromyces marneffei is a thermally dimorphic fungus endemic in China and Southeast Asia that causes fatal infections in immunocompromised individuals, particularly in patients with advanced HIV disease. Here, we report the complete genome sequences of two clinical isolates from northern and southern Vietnam.

Cryptococcus neoformans Recovered From Olive Trees (Olea europaea) in Turkey Reveal Allopatry With African and South American Lineages

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 MAT a) 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 MAT a C. 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.

Genetic and Genomic Analyses Reveal Boundaries between Species Closely Related to Cryptococcus Pathogens

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.

Coordinated host-pathogen transcriptional dynamics revealed using sorted subpopulations and single macrophages infected with Candida albicans

The outcome of fungal infections depends on interactions with innate immune cells. Within a population of macrophages encountering Candida albicans, there are distinct host-pathogen trajectories; however, little is known about the molecular heterogeneity that governs these fates. Here we developed an experimental system to separate interaction stages and single macrophage cells infected with C. albicans from uninfected cells and assessed transcriptional variability in the host and fungus. Macrophages displayed an initial up-regulation of pathways involved in phagocytosis and proinflammatory response after C. albicans exposure that declined during later time points. Phagocytosed C. albicans shifted expression programs to survive the nutrient poor phagosome and remodeled the cell wall. The transcriptomes of single infected macrophages and phagocytosed C. albicans displayed a tightly coordinated shift in gene expression co-stages and revealed expression bimodality and differential splicing that may drive infection outcome. This work establishes an approach for studying host-pathogen trajectories to resolve heterogeneity in dynamic populations.

The outcomes of the interactions between individual host cells and pathogens are heterogeneous. Here, the authors assess transcriptional variability in both host and pathogen during infection of macrophages with the fungus Candida albicans, using sorted subpopulations and single macrophages.

Pathways of Pathogenicity: Transcriptional Stages of Germination in the Fatal Fungal Pathogen Rhizopus delemar

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.

Genomic characterization of recurrent mold infections in thoracic transplant recipients

Invasive mold disease in thoracic organ transplant recipients is a well-recognized complication, but the long-term persistence of molds within the human body and evasion of host defenses has not been well-described. We present 2 cases of invasive mold disease (Verruconis gallopava and Aspergillus fumigatus) in thoracic transplant recipients who had the same mold cultured years prior to the invasive disease presentation. The paired isolates from the index and recurrent infections in both patients were compared using whole-genome sequencing to determine if the same strain of mold caused both the index and recurrent infections. In Case 1, the isolates were found to be of the same strain indicating that the initial colonizing isolate identified pre-transplant eventually caused invasive mold disease post-transplant while in Case 2, the 2 isolates were not of the same strain. These results demonstrate the distinct possibility of molds both persisting within the human body for years prior to invasive mold disease or the long-term risk of recurrent, persistent infection with more than one strain. Further studies of long-term molecular epidemiology of IMD and risk factors for mold persistence in transplant recipients are encouraged.

Titan cells formation in Cryptococcus neoformans is finely tuned by environmental conditions and modulated by positive and negative genetic regulators

The pathogenic fungus Cryptococcus neoformans exhibits morphological changes in cell size during lung infection, producing both typical size 5 to 7 μm cells and large titan cells (> 10 μm and up to 100 μm). We found and optimized in vitro conditions that produce titan cells in order to identify the ancestry of titan cells, the environmental determinants, and the key gene regulators of titan cell formation. Titan cells generated in vitro harbor the main characteristics of titan cells produced in vivo including their large cell size (>10 μm), polyploidy with a single nucleus, large vacuole, dense capsule, and thick cell wall. Here we show titan cells derived from the enlargement of progenitor cells in the population independent of yeast growth rate. Change in the incubation medium, hypoxia, nutrient starvation and low pH were the main factors that trigger titan cell formation, while quorum sensing factors like the initial inoculum concentration, pantothenic acid, and the quorum sensing peptide Qsp1p also impacted titan cell formation. Inhibition of ergosterol, protein and nucleic acid biosynthesis altered titan cell formation, as did serum, phospholipids and anti-capsular antibodies in our settings. We explored genetic factors important for titan cell formation using three approaches. Using H99-derivative strains with natural genetic differences, we showed that titan cell formation was dependent on LMP1 and SGF29 genes. By screening a gene deletion collection, we also confirmed that GPR4/5-RIM101, and CAC1 genes were required to generate titan cells and that the PKR1, TSP2, USV101 genes negatively regulated titan cell formation. Furthermore, analysis of spontaneous Pkr1 loss-of-function clinical isolates confirmed the important role of the Pkr1 protein as a negative regulator of titan cell formation. Through development of a standardized and robust in vitro assay, our results provide new insights into titan cell biogenesis with the identification of multiple important factors/pathways.

Whole-Genome Analysis Illustrates Global Clonal Population Structure of the Ubiquitous Dermatophyte Pathogen Trichophyton rubrum

Dermatophytes include fungal species that infect humans, as well as those that also infect other animals or only grow in the environment. The dermatophyte species Trichophyton rubrum is a frequent cause of skin infection in immunocompetent individuals. While members of the T. rubrum species complex have been further categorized based on various morphologies, their population structure and ability to undergo sexual reproduction are not well understood. In this study, we analyze a large set of T. rubrum and T. interdigitale isolates to examine mating types, evidence of mating, and genetic variation. We find that nearly all isolates of T. rubrum are of a single mating type, and that incubation with T. rubrum "morphotype" megninii isolates of the other mating type failed to induce sexual development. While the region around the mating type locus is characterized by a higher frequency of SNPs compared to other genomic regions, we find that the population is remarkably clonal, with highly conserved gene content, low levels of variation, and little evidence of recombination. These results support a model of recent transition to asexual growth when this species specialized to growth on human hosts.

Population genomics and the evolution of virulence in the fungal pathogen Cryptococcus neoformans

Cryptococcus neoformans is an opportunistic fungal pathogen that causes approximately 625,000 deaths per year from nervous system infections. Here, we leveraged a unique, genetically diverse population of C. neoformans from sub-Saharan Africa, commonly isolated from mopane trees, to determine how selective pressures in the environment coincidentally adapted C. neoformans for human virulence. Genome sequencing and phylogenetic analysis of 387 isolates, representing the global VNI and African VNB lineages, highlighted a deep, nonrecombining split in VNB (herein, VNBI and VNBII). VNBII was enriched for clinical samples relative to VNBI, while phenotypic profiling of 183 isolates demonstrated that VNBI isolates were significantly more resistant to oxidative stress and more heavily melanized than VNBII isolates. Lack of melanization in both lineages was associated with loss-of-function mutations in the BZP4 transcription factor. A genome-wide association study across all VNB isolates revealed sequence differences between clinical and environmental isolates in virulence factors and stress response genes. Inositol transporters and catabolism genes, which process sugars present in plants and the human nervous system, were identified as targets of selection in all three lineages. Further phylogenetic and population genomic analyses revealed extensive loss of genetic diversity in VNBI, suggestive of a history of population bottlenecks, along with unique evolutionary trajectories for mating type loci. These data highlight the complex evolutionary interplay between adaptation to natural environments and opportunistic infections, and that selection on specific pathways may predispose isolates to human virulence.

Advances in Cryptococcus genomics: insights into the evolution of pathogenesis

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.

Novel taxa of thermally dimorphic systemic pathogens in the Ajellomycetaceae (Onygenales)

Recent discoveries of novel systemic fungal pathogens with thermally dimorphic yeast-like phases have challenged the current taxonomy of the Ajellomycetaceae, a family currently comprising the genera Blastomyces, Emmonsia, Emmonsiellopsis, Helicocarpus, Histoplasma, Lacazia and Paracoccidioides. Our morphological, phylogenetic and phylogenomic analyses demonstrated species relationships and their specific phenotypes, clarified generic boundaries and provided the first annotated genome assemblies to support the description of two new species. A new genus, Emergomyces, accommodates Emmonsia pasteuriana as type species, and the new species Emergomyces africanus, the aetiological agent of case series of disseminated infections in South Africa. Both species produce small yeast cells that bud at a narrow base at 37°C and lack adiaspores, classically associated with the genus Emmonsia. Another novel dimorphic pathogen, producing broad-based budding cells at 37°C and occurring outside North America, proved to belong to the genus Blastomyces, and is described as Blastomyces percursus.

Genomic innovations linked to infection strategies across emerging pathogenic chytrid fungi

To understand the evolutionary pathways that lead to emerging infections of vertebrates, here we explore the genomic innovations that allow free-living chytrid fungi to adapt to and colonize amphibian hosts. Sequencing and comparing the genomes of two pathogenic species of Batrachochytrium to those of close saprophytic relatives reveals that pathogenicity is associated with remarkable expansions of protease and cell wall gene families, while divergent infection strategies are linked to radiations of lineage-specific gene families. By comparing the host–pathogen response to infection for both pathogens, we illuminate the traits that underpin a strikingly different immune response within a shared host species. Our results show that, despite commonalities that promote infection, specific gene-family radiations contribute to distinct infection strategies. The breadth and evolutionary novelty of candidate virulence factors that we discover underscores the urgent need to halt the advance of pathogenic chytrids and prevent incipient loss of biodiversity.

Batrachochytrium dendrobatidis and B. salamandrivorans are both important pathogens of amphibians, but they differ in their host ranges, infection strategies, and host immune responses. Here, Farrer and colleagues compare their genomes and transcriptomes to identify the genetic basis of these differences.

Comparative Genomics Integrated with Association Analysis Identifies Candidate Effector Genes Corresponding to Lr20 in Phenotype-Paired Puccinia triticina Isolates from Australia

Leaf rust is one of the most common and damaging diseases of wheat, and is caused by an obligate biotrophic basidiomycete, Puccinia triticina (Pt). In the present study, 20 Pt isolates from Australia, comprising 10 phenotype-matched pairs with contrasting pathogenicity for Lr20, were analyzed using whole genome sequencing. Compared to the reference genome of the American Pt isolate 1-1 BBBD Race 1, an average of 404,690 single nucleotide polymorphisms (SNPs) per isolate was found and the proportion of heterozygous SNPs was above 87% in the majority of the isolates, demonstrating a high level of polymorphism and a high rate of heterozygosity. From the genome-wide SNPs, a phylogenetic tree was inferred, which consisted of a large clade of 15 isolates representing diverse presumed clonal lineages including 14 closely related isolates and the more diverged isolate 670028, and a small clade of five isolates characterized by lower heterozygosity level. Principle component analysis detected three distinct clusters, corresponding exactly to the two major subsets of the small clade and the large clade comprising all 15 isolates without further separation of isolate 670028. While genome-wide association analysis identified 302 genes harboring at least one SNP associated with Lr20 virulence (p < 0.05), a Wilcoxon rank sum test revealed that 36 and 68 genes had significant (p < 0.05) and marginally significant (p < 0.1) differences in the counts of non-synonymous mutations between Lr20 avirulent and virulent groups, respectively. Twenty of these genes were predicted to have a signal peptide without a transmembrane segment, and hence identified as candidate effector genes corresponding to Lr20. SNP analysis also implicated the potential involvement of epigenetics and small RNA in Pt pathogenicity. Future studies are thus warranted to investigate the biological functions of the candidate effectors as well as the gene regulation mechanisms at epigenetic and post-transcription levels. Our study is the first to integrate phenotype-genotype association with effector prediction in Pt genomes, an approach that may circumvent some of the technical difficulties in working with obligate rust fungi and accelerate avirulence gene identification.

Transcriptional analysis of mating and pre-infection stages of the anther smut, Microbotryum lychnidis-dioicae

Microbotryum lychnidis-dioicae is an obligate biotrophic parasite of the wildflower species Silene latifolia. This dikaryotic fungus, commonly known as an anther smut, requires that haploid, yeast-like sporidia of opposite mating types fuse and differentiate into dikaryotic hyphae that penetrate host tissue as part of the fungal life cycle. Mating occurs under conditions of cool temperatures and limited nutrients. Further development requires host cues or chemical mimics, including a variety of lipids, e.g. phytols. To identify global changes in transcription associated with developmental shifts, RNA-Seq was conducted at several in vitro stages of fungal propagation, i.e. haploid cells grown independently on rich and nutrient-limited media, mated cells on nutrient-limited media as well as a time course of such mated cells exposed to phytol. Comparison of haploid cells grown under rich and nutrient-limited conditions identified classes of genes probably associated with general nutrient availability, including components of the RNAi machinery. Some gene enrichment patterns comparing the nutrient-limited and mated transcriptomes suggested gene expression changes associated with the mating programme (e.g. homeodomain binding proteins, secreted proteins, proteins unique to M. lychnidis-dioicae¸ multicopper oxidases and RhoGEFs). Analysis for phytol treatment compared with mated cells alone allowed identification of genes likely to be involved in the dikaryotic switch (e.g. oligopeptide transporters). Gene categories of particular note in all three conditions included those in the major facilitator superfamily, proteins containing PFAM domains of the secretory lipase family as well as proteins predicted to be secreted, many of which have the hallmarks of fungal effectors with potential roles in pathogenicity.

Paracoccidioides spp. catalases and their role in antioxidant defense against host defense responses

Dimorphic human pathogenic fungi interact with host effector cells resisting their microbicidal mechanisms. Yeast cells are able of surviving within the tough environment of the phagolysosome by expressing an antioxidant defense system that provides protection against host-derived reactive oxygen species (ROS). This includes the production of catalases (CATs). Here we identified and analyzed the role of CAT isoforms in Paracoccidioides, the etiological agent of paracoccidioidomycosis. Firstly, we found that one of these isoforms was absent in the closely related dimorphic pathogen Coccidioides and dermatophytes, but all of them were conserved in Paracoccidioides, Histoplasma and Blastomyces species. We probed the contribution of CATs in Paracoccidioides by determining the gene expression levels of each isoform through quantitative RT-qPCR, in both the yeast and mycelia phases, and during the morphological switch (transition and germination), as well as in response to oxidative agents and during interaction with neutrophils. PbCATP was preferentially expressed in the pathogenic yeast phase, and was associated to the response against exogenous H₂O₂. Therefore, we created and analyzed the virulence defects of a knockdown strain for this isoform, and found that CATP protects yeast cells from H₂O₂ generated in vitro and is relevant during lung infection. On the other hand, CATA and CATB seem to contribute to ROS homeostasis in Paracoccidioides cells, during endogenous oxidative stress. CAT isoforms in Paracoccidioides might be coordinately regulated during development and dimorphism, and differentially expressed in response to different stresses to control ROS homeostasis during the infectious process, contributing to the virulence of Paracoccidioides.

Microevolutionary traits and comparative population genomics of the emerging pathogenic fungus Cryptococcus gattii

Emerging fungal pathogens cause an expanding burden of disease across the animal kingdom, including a rise in morbidity and mortality in humans. Yet, we currently have only a limited repertoire of available therapeutic interventions. A greater understanding of the mechanisms of fungal virulence and of the emergence of hypervirulence within species is therefore needed for new treatments and mitigation efforts. For example, over the past decade, an unusual lineage of Cryptococcus gattii, which was first detected on Vancouver Island, has spread to the Canadian mainland and the Pacific Northwest infecting otherwise healthy individuals. The molecular changes that led to the development of this hypervirulent cryptococcal lineage remain unclear. To explore this, we traced the history of similar microevolutionary events that can lead to changes in host range and pathogenicity. Here, we detail fine-resolution mapping of genetic differences between two highly related Cryptococcus gattii VGIIc isolates that differ in their virulence traits (phagocytosis, vomocytosis, macrophage death, mitochondrial tubularization and intracellular proliferation). We identified a small number of single site variants within coding regions that potentially contribute to variations in virulence. We then extended our methods across multiple lineages of C. gattii to study how selection is acting on key virulence genes within different lineages.This article is part of the themed issue 'Tackling emerging fungal threats to animal health, food security and ecosystem resilience'.

Simultaneous Emergence of Multidrug-Resistant Candida auris on 3 Continents Confirmed by Whole-Genome Sequencing and Epidemiological Analyses

BACKGROUND

Candida auris, a multidrug-resistant yeast that causes invasive infections, was first described in 2009 in Japan and has since been reported from several countries.

METHODS

To understand the global emergence and epidemiology of C. auris, we obtained isolates from 54 patients with C. auris infection from Pakistan, India, South Africa, and Venezuela during 2012-2015 and the type specimen from Japan. Patient information was available for 41 of the isolates. We conducted antifungal susceptibility testing and whole-genome sequencing (WGS).

RESULTS

Available clinical information revealed that 41% of patients had diabetes mellitus, 51% had undergone recent surgery, 73% had a central venous catheter, and 41% were receiving systemic antifungal therapy when C. auris was isolated. The median time from admission to infection was 19 days (interquartile range, 9-36 days), 61% of patients had bloodstream infection, and 59% died. Using stringent break points, 93% of isolates were resistant to fluconazole, 35% to amphotericin B, and 7% to echinocandins; 41% were resistant to 2 antifungal classes and 4% were resistant to 3 classes. WGS demonstrated that isolates were grouped into unique clades by geographic region. Clades were separated by thousands of single-nucleotide polymorphisms, but within each clade isolates were clonal. Different mutations in ERG11 were associated with azole resistance in each geographic clade.

CONCLUSIONS

C. auris is an emerging healthcare-associated pathogen associated with high mortality. Treatment options are limited, due to antifungal resistance. WGS analysis suggests nearly simultaneous, and recent, independent emergence of different clonal populations on 3 continents. Risk factors and transmission mechanisms need to be elucidated to guide control measures.

A Multistate Toggle Switch Defines Fungal Cell Fates and Is Regulated by Synergistic Genetic Cues

Heritable epigenetic changes underlie the ability of cells to differentiate into distinct cell types. Here, we demonstrate that the fungal pathogen Candida tropicalis exhibits multipotency, undergoing stochastic and reversible switching between three cellular states. The three cell states exhibit unique cellular morphologies, growth rates, and global gene expression profiles. Genetic analysis identified six transcription factors that play key roles in regulating cell differentiation. In particular, we show that forced expression of Wor1 or Efg1 transcription factors can be used to manipulate transitions between all three cell states. A model for tristability is proposed in which Wor1 and Efg1 are self-activating but mutually antagonistic transcription factors, thereby forming a symmetrical self-activating toggle switch. We explicitly test this model and show that ectopic expression of WOR1 can induce white-to-hybrid-to-opaque switching, whereas ectopic expression of EFG1 drives switching in the opposite direction, from opaque-to-hybrid-to-white cell states. We also address the stability of induced cell states and demonstrate that stable differentiation events require ectopic gene expression in combination with chromatin-based cues. These studies therefore experimentally test a model of multistate stability and demonstrate that transcriptional circuits act synergistically with chromatin-based changes to drive cell state transitions. We also establish close mechanistic parallels between phenotypic switching in unicellular fungi and cell fate decisions during stem cell reprogramming.

Genome Diversity, Recombination, and Virulence across the Major Lineages of Paracoccidioides

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.

Lower prevalence but similar fitness in a parasitic fungus at higher radiation levels near Chernobyl

Nuclear disasters at Chernobyl and Fukushima provide examples of effects of acute ionizing radiation on mutations that can affect the fitness and distribution of species. Here, we investigated the prevalence of Microbotryum lychnidis-dioicae, a pollinator-transmitted fungal pathogen of plants causing anther-smut disease in Chernobyl, its viability, fertility and karyotype variation, and the accumulation of nonsynonymous mutations in its genome. We collected diseased flowers of Silene latifolia from locations ranging by more than two orders of magnitude in background radiation, from 0.05 to 21.03 μGy/h. Disease prevalence decreased significantly with increasing radiation level, possibly due to lower pollinator abundance and altered pollinator behaviour. Viability and fertility, measured as the budding rate of haploid sporidia following meiosis from the diploid teliospores, did not vary with increasing radiation levels and neither did karyotype overall structure and level of chromosomal size heterozygosity. We sequenced the genomes of twelve samples from Chernobyl and of four samples collected from uncontaminated areas and analysed alignments of 6068 predicted genes, corresponding to 1.04 × 10(7)  base pairs. We found no dose-dependent differences in substitution rates (neither dN, dS, nor dN/dS). Thus, we found no significant evidence of increased deleterious mutation rates at higher levels of background radiation in this plant pathogen. We even found lower levels of nonsynonymous substitution rates in contaminated areas compared to control regions, suggesting that purifying selection was stronger in contaminated than uncontaminated areas. We briefly discuss the possibilities for a mechanistic basis of radio resistance in this nonmelanized fungus.

Hsf1 and Hsp90 orchestrate temperature-dependent global transcriptional remodelling and chromatin architecture in Candida albicans

Fever is a universal response to infection, and opportunistic pathogens such as Candida albicans have evolved complex circuitry to sense and respond to heat. Here we harness RNA-seq and ChIP-seq to discover that the heat shock transcription factor, Hsf1, binds distinct motifs in nucleosome-depleted promoter regions to regulate heat shock genes and genes involved in virulence in C. albicans. Consequently, heat shock increases C. albicans host cell adhesion, damage and virulence. Hsf1 activation depends upon the molecular chaperone Hsp90 under basal and heat shock conditions, but the effects are opposite and in part controlled at the level of Hsf1 expression and DNA binding. Finally, we demonstrate that Hsp90 regulates global transcription programs by modulating nucleosome levels at promoters of stress-responsive genes. Thus, we describe a mechanism by which C. albicans responds to temperature via Hsf1 and Hsp90 to orchestrate gene expression and chromatin architecture, thereby enabling thermal adaptation and virulence.

The Dynamic Genome and Transcriptome of the Human Fungal Pathogen Blastomyces and Close Relative Emmonsia

Three closely related thermally dimorphic pathogens are causal agents of major fungal diseases affecting humans in the Americas: blastomycosis, histoplasmosis and paracoccidioidomycosis. Here we report the genome sequence and analysis of four strains of the etiological agent of blastomycosis, Blastomyces, and two species of the related genus Emmonsia, typically pathogens of small mammals. Compared to related species, Blastomyces genomes are highly expanded, with long, often sharply demarcated tracts of low GC-content sequence. These GC-poor isochore-like regions are enriched for gypsy elements, are variable in total size between isolates, and are least expanded in the avirulent B. dermatitidis strain ER-3 as compared with the virulent B. gilchristii strain SLH14081. The lack of similar regions in related species suggests these isochore-like regions originated recently in the ancestor of the Blastomyces lineage. While gene content is highly conserved between Blastomyces and related fungi, we identified changes in copy number of genes potentially involved in host interaction, including proteases and characterized antigens. In addition, we studied gene expression changes of B. dermatitidis during the interaction of the infectious yeast form with macrophages and in a mouse model. Both experiments highlight a strong antioxidant defense response in Blastomyces, and upregulation of dioxygenases in vivo suggests that dioxide produced by antioxidants may be further utilized for amino acid metabolism. We identify a number of functional categories upregulated exclusively in vivo, such as secreted proteins, zinc acquisition proteins, and cysteine and tryptophan metabolism, which may include critical virulence factors missed before in in vitro studies. Across the dimorphic fungi, loss of certain zinc acquisition genes and differences in amino acid metabolism suggest unique adaptations of Blastomyces to its host environment. These results reveal the dynamics of genome evolution and of factors contributing to virulence in Blastomyces.

Dimorphic fungal pathogens including Blastomyces are the cause of major fungal diseases in North and South America. The genus Emmonsia includes species infecting small mammals as well as a newly emerging pathogenic species recently reported in HIV-positive patients in South Africa. Here, we synthesize both genome sequencing of four isolates of Blastomyces and two species of Emmonsia as well as deep sequencing of Blastomyces RNA to draw major new insights into the evolution of this group and the pathogen response to infection. We investigate the trajectory of genome evolution of this group, characterizing the phylogenetic relationships of these species, a remarkable genome expansion that formed large isochore-like regions of low GC content in Blastomyces, and variation of gene content, related to host interaction, among the dimorphic fungal pathogens. Using RNA-Seq, we profile the response of Blastomyces to macrophage and mouse pulmonary infection, identifying key pathways and novel virulence factors. The identification of key fungal genes involved in adaptation to the host suggests targets for further study and therapeutic intervention in Blastomyces and related dimorphic fungal pathogens.

Sex and parasites: genomic and transcriptomic analysis of Microbotryum lychnidis-dioicae, the biotrophic and plant-castrating anther smut fungus

Background

The genus Microbotryum includes plant pathogenic fungi afflicting a wide variety of hosts with anther smut disease. Microbotryum lychnidis-dioicae infects Silene latifolia and replaces host pollen with fungal spores, exhibiting biotrophy and necrosis associated with altering plant development.

Results

We determined the haploid genome sequence for M. lychnidis-dioicae and analyzed whole transcriptome data from plant infections and other stages of the fungal lifecycle, revealing the inventory and expression level of genes that facilitate pathogenic growth. Compared to related fungi, an expanded number of major facilitator superfamily transporters and secretory lipases were detected; lipase gene expression was found to be altered by exposure to lipid compounds, which signaled a switch to dikaryotic, pathogenic growth. In addition, while enzymes to digest cellulose, xylan, xyloglucan, and highly substituted forms of pectin were absent, along with depletion of peroxidases and superoxide dismutases that protect the fungus from oxidative stress, the repertoire of glycosyltransferases and of enzymes that could manipulate host development has expanded. A total of 14 % of the genome was categorized as repetitive sequences. Transposable elements have accumulated in mating-type chromosomal regions and were also associated across the genome with gene clusters of small secreted proteins, which may mediate host interactions.

Conclusions

The unique absence of enzyme classes for plant cell wall degradation and maintenance of enzymes that break down components of pollen tubes and flowers provides a striking example of biotrophic host adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1660-8) contains supplementary material, which is available to authorized users.