Analysis of the genome and transcriptome of Cryptococcus neoformans var. grubii reveals complex RNA expression and microevolution leading to virulence attenuation

Hsp90 interaction networks in fungi-tools and techniques

Heat-shock protein 90 (Hsp90) is a central regulator of cellular proteostasis. It stabilizes numerous proteins that are involved in fundamental processes of life, including cell growth, cell-cycle progression and the environmental response. In addition to stabilizing proteins, Hsp90 governs gene expression and controls the release of cryptic genetic variation. Given its central role in evolution and development, it is important to identify proteins and genes that interact with Hsp90. This requires sophisticated genetic and biochemical tools, including extensive mutant collections, suitable epitope tags, proteomics approaches and Hsp90-specific pharmacological inhibitors for chemogenomic screens. These usually only exist in model organisms, such as the yeast Saccharomyces cerevisiae. Yet, the importance of other fungal species, such as Candida albicans and Cryptococcus neoformans, as serious human pathogens accelerated the development of genetic tools to study their virulence and stress response pathways. These tools can also be exploited to map Hsp90 interaction networks. Here, we review tools and techniques for Hsp90 network mapping available in different fungi and provide a summary of existing mapping efforts. Mapping Hsp90 networks in fungal species spanning >500 million years of evolution provides a unique vantage point, allowing tracking of the evolutionary history of eukaryotic Hsp90 networks.

Intron distribution and emerging role of alternative splicing in fungi

Spliceosomal introns are noncoding sequences that are spliced from pre-mRNA. They are ubiquitous in eukaryotic genomes, although the average number of introns per gene varies considerably between different eukaryotic species. Fungi are diverse in terms of intron numbers ranging from 4% to 99% genes with introns. Alternative splicing is one of the most common modes of posttranscriptional regulation in eukaryotes, giving rise to multiple transcripts from a single pre-mRNA and is widespread in metazoans and drives extensive proteome diversity. Earlier, alternative splicing was considered to be rare in fungi, but recently, increasing numbers of studies have revealed that alternative splicing is also widespread in fungi and has been implicated in the regulation of fungal growth and development, protein localization and the improvement of survivability, likely underlying their unique capacity to adapt to changing environmental conditions. However, the role of alternative splicing in pathogenicity and development of drug resistance is only recently gaining attention. In this review, we describe the intronic landscape in fungi. We also present in detail the newly discovered functions of alternative splicing in various cellular processes and outline areas particularly in pathogenesis and clinical drug resistance for future studies that could lead to the development of much needed new therapeutics.

Comparative analysis of RNA enrichment methods for preparation of Cryptococcus neoformans RNA sequencing libraries

RNA sequencing (RNA-Seq) experiments focused on gene expression involve removal of ribosomal RNA (rRNA) because it is the major RNA constituent of cells. This process, called RNA enrichment, is done primarily to reduce cost: without rRNA removal, deeper sequencing must be performed to compensate for the sequencing reads wasted on rRNA. The ideal RNA enrichment method removes all rRNA without affecting other RNA in the sample. We tested the performance of three RNA enrichment methods on RNA isolated from Cryptococcus neoformans, a fungal pathogen of humans. We find that the RNase H depletion method is more efficient in depleting rRNA and more specific in recapitulating non-rRNA levels present in unenriched controls than the commonly-used Poly(A) isolation method. The RNase H depletion method is also more effective than the Ribo-Zero depletion method as measured by rRNA depletion efficiency and recapitulation of protein-coding RNA levels present in unenriched controls, while the Ribo-Zero depletion method more closely recapitulates annotated non-coding RNA (ncRNA) levels. Finally, we leverage these data to accurately map the C. neoformans mitochondrial rRNA genes, and also demonstrate that RNA-Seq data generated with the RNase H and Ribo-Zero depletion methods can be used to explore novel C. neoformans long non-coding RNA genes.

Mechanics of microtubule organizing center clustering and spindle positioning in budding yeast Cryptococcus neoformans

The dynamic process of mitotic spindle assembly depends on multitudes of inter-dependent interactions involving kinetochores (KTs), microtubules (MTs), spindle pole bodies (SPBs), and molecular motors. Before forming the mitotic spindle, multiple visible microtubule organizing centers (MTOCs) coalesce into a single focus to serve as an SPB in the pathogenic budding yeast, Cryptococcus neoformans. To explain this unusual phenomenon in the fungal kingdom, we propose a "search and capture" model, in which cytoplasmic MTs (cMTs) nucleated by MTOCs grow and capture each other to promote MTOC clustering. Our quantitative modeling identifies multiple redundant mechanisms mediated by a combination of cMT-cell cortex interactions and inter-cMT coupling to facilitate MTOC clustering within the physiological time limit as determined by time-lapse live-cell microscopy. Besides, we screen various possible mechanisms by computational modeling and propose optimal conditions that favor proper spindle positioning-a critical determinant for timely chromosome segregation. These analyses also reveal that a combined effect of MT buckling, dynein pull, and cortical push maintains spatiotemporal spindle localization.

Docusate-Based Ionic Liquids of Anthelmintic Benzimidazoles Show Improved Pharmaceutical Processability, Lipid Solubility, and in Vitro Activity against Cryptococcus neoformans

As the existing therapeutic modalities for the treatment of cryptococcal meningitis (CM) have suboptimal efficacy, repurposing existing drugs for the treatment of CM is of great interest. The FDA-approved anthelmintic benzimidazoles, albendazole, mebendazole, and flubendazole, have demonstrated potent but variable in vitro activity against Cryptococcus neoformans, the predominant fungal species responsible for CM. We performed molecular docking studies to ascertain the interaction of albendazole, mebendazole, and flubendazole with a C. neoformans β-tubulin structure, which revealed differential binding interactions and explained the different in vitro efficacies reported previously and observed in this investigation. Despite their promising in vitro efficacy, the repurposing of anthelmintic benzimidazoles for oral CM therapy is significantly hampered due to their high crystallinity, poor pharmaceutical processability, low and pH-dependent solubility, and drug precipitation upon entering the intestine, all of which result in low and variable oral bioavailability. Here, we demonstrate that the anthelmintic benzimidazoles can be transformed into partially amorphous low-melting ionic liquids (ILs) with a simple metathesis reaction using amphiphilic sodium docusate as a counterion. In vitro efficacy studies on a laboratory reference and a clinical isolate of C. neoformans showed 2- to 4-fold lower IC90 values for docusate-based ILs compared to the pure anthelmintic benzimidazoles. Furthermore, using a C. neoformans strain with green fluorescent protein (GFP)-tagged β-tubulin and albendazole and its docusate IL as model candidates, we showed that the benzimidazoles and their ILs reduce the viability of C. neoformans by interfering with its microtubule assembly. Unlike pure anthelmintic benzimidazoles, the docusate-based ILs showed excellent solubility in organic solvents and >30-fold higher solubility in bioavailability-enhancing lipid vehicles. Finally, the docusate ILs were successfully incorporated into SoluPlus, a self-assembling biodegradable polymer, which upon dilution with water formed polymeric micelles with a size of <100 nm. Thus, the development of docusate-based ILs represents an effective approach to improve the physicochemical properties and potency of anthelmintic benzimidazoles to facilitate their repurposing and preclinical development for CM therapy.

Epigenetic dynamics of centromeres and neocentromeres in Cryptococcus deuterogattii

Deletion of native centromeres in the human fungal pathogen Cryptococcus deuterogattii leads to neocentromere formation. Native centromeres span truncated transposable elements, while neocentromeres do not and instead span actively expressed genes. To explore the epigenetic organization of neocentromeres, we analyzed the distribution of the heterochromatic histone modification H3K9me2, 5mC DNA methylation and the euchromatin mark H3K4me2. Native centromeres are enriched for both H3K9me2 and 5mC DNA methylation marks and are devoid of H3K4me2, while neocentromeres do not exhibit any of these features. Neocentromeres in cen10Δ mutants are unstable and chromosome-chromosome fusions occur. After chromosome fusion, the neocentromere is inactivated and the native centromere of the chromosome fusion partner remains as the sole, active centromere. In the present study, the active centromere of a fused chromosome was deleted to investigate if epigenetic memory promoted the re-activation of the inactive neocentromere. Our results show that the inactive neocentromere is not re-activated and instead a novel neocentromere forms directly adjacent to the deleted centromere of the fused chromosome. To study the impact of transcription on centromere stability, the actively expressed URA5 gene was introduced into the CENP-A bound regions of a native centromere. The introduction of the URA5 gene led to a loss of CENP-A from the native centromere, and a neocentromere formed adjacent to the native centromere location. Remarkably, the inactive, native centromere remained enriched for heterochromatin, yet the integrated gene was expressed and devoid of H3K9me2. A cumulative analysis of multiple CENP-A distribution profiles revealed centromere drift in C. deuterogattii, a previously unreported phenomenon in fungi. The CENP-A-binding shifted within the ORF-free regions and showed a possible association with a truncated transposable element. Taken together, our findings reveal that neocentromeres in C. deuterogattii are highly unstable and are not marked with an epigenetic memory, distinguishing them from native centromeres.

Linear eukaryotic chromosomes require a specific chromosomal region, the centromere, where the macromolecular kinetochore protein complex assembles. In most organisms, centromeres are located in gene-free, repeat-rich chromosomal regions and may or may not be associated with heterochromatic epigenetic marks. We report that the native centromeres of the human fungal pathogen Cryptococcus deuterogattii are enriched with heterochromatin marks. Deleting a centromere in C. deuterogattii results in formation of neocentromeres that span genes. In some cases, neocentromeres are unstable leading to chromosome-chromosome fusions and neocentromere inactivation. These neocentromeres, unlike native centromeres, lack any of the tested heterochromatic marks or any epigenetic memory. We also found that neocentromere formation can be triggered not only by deletion of the native centromere but also by disrupting its function via insertion of a gene. These results show that neocentromere dynamics in this fungal pathogen are unique among organisms studied so far. Our results also revealed key differences between epigenetics of native centromeres between C. deuterogattii and its sister species, C. neoformans. These finding provide an opportunity to test and study the evolution of centromeres, as well as neocentromeres, in this species complex and how it might contribute to their genome evolution.

Structural features of Cryptococcus neoformans bifunctional GAR/AIR synthetase may present novel antifungal drug targets

Cryptococcus neoformans is a fungus that causes life-threatening systemic mycoses. During infection of the human host, this pathogen experiences a major change in the availability of purines; the fungus can scavenge the abundant purines in its environmental niche of pigeon excrement, but must employ de novo biosynthesis in the purine-poor human CNS. Eleven sequential enzymatic steps are required to form the first purine base, IMP, an intermediate in the formation of ATP and GTP. Over the course of evolution, several gene fusion events led to the formation of multifunctional purine biosynthetic enzymes in most organisms, particularly the higher eukaryotes. In C. neoformans, phosphoribosyl-glycinamide synthetase (GARs) and phosphoribosyl-aminoimidazole synthetase (AIRs) are fused into a bifunctional enzyme, while the human ortholog is a trifunctional enzyme that also includes GAR transformylase. Here we functionally, biochemically, and structurally characterized C. neoformans GARs and AIRs to identify drug targetable features. GARs/AIRs are essential for de novo purine production and virulence in a murine inhalation infection model. Characterization of GARs enzymatic functional parameters showed that C. neoformans GARs/AIRs have lower affinity for substrates glycine and PRA compared with the trifunctional metazoan enzyme. The crystal structure of C. neoformans GARs revealed differences in the glycine- and ATP-binding sites compared with the Homo sapiens enzyme, while the crystal structure of AIRs shows high structural similarity compared with its H. sapiens ortholog as a monomer but differences as a dimer. The alterations in functional and structural characteristics between fungal and human enzymes could potentially be exploited for antifungal development.

Epitope-Based Immunoinformatic Approach on Heat Shock 70 kDa Protein Complex of Cryptococcus neoformans var. grubii

INTRODUCTION

Cryptococcosis is a ubiquitous opportunistic fungal disease caused by Cryptococcus neoformans var. grubii. It has high global morbidity and mortality among HIV patients and non-HIV carriers with 99% and 95%, respectively. Furthermore, the increasing prevalence of undesired toxicity profile of antifungal, multidrug-resistant organisms and the scarcity of FDA-authorized vaccines were the hallmark in the present days. This study was undertaken to design a reliable epitope-based peptide vaccine through targeting highly conserved immunodominant heat shock 70 kDa protein of Cryptococcus neoformans var. grubii that covers a considerable digit of the world population through implementing a computational vaccinology approach.

MATERIALS AND METHODS

A total of 38 sequences of Cryptococcus neoformans var. grubii's heat shock 70 kDa protein were retrieved from the NCBI protein database. Different prediction tools were used to analyze the aforementioned protein at the Immune Epitope Database (IEDB) to discriminate the most promising T-cell and B-cell epitopes. The proposed T-cell epitopes were subjected to the population coverage analysis tool to compute the global population's coverage. Finally, the T-cell projected epitopes were ranked based on their binding scores and modes using AutoDock Vina software. Results and Discussion. The epitopes (ANYVQASEK, QSEKPKNVNPVI, SEKPKNVNPVI, and EKPKNVNPVI) had shown very strong binding affinity and immunogenic properties to B-cell. (FTQLVAAYL, YVYDTRGKL) and (FFGGKVLNF, FINAQLVDV, and FDYALVQHF) exhibited a very strong binding affinity to MHC-I and MHC-II, respectively, with high population coverage for each, while FYRQGAFEL has shown promising results in terms of its binding profile to MHC-II and MHC-I alleles and good strength of binding when docked with HLA-C∗12:03. In addition, there is massive global population coverage in the three coverage modes. Accordingly, our in silico vaccine is expected to be the future epitope-based peptide vaccine against Cryptococcus neoformans var. grubii that covers a significant figure of the entire world citizens.

Virulence Traits and Population Genomics of the Black Yeast Aureobasidium melanogenum

The black yeast-like fungus Aureobasidium melanogenum is an opportunistic human pathogen frequently found indoors. Its traits, potentially linked to pathogenesis, have never been systematically studied. Here, we examine 49 A. melanogenum strains for growth at 37 °C, siderophore production, hemolytic activity, and assimilation of hydrocarbons and human neurotransmitters and report within-species variability. All but one strain grew at 37 °C. All strains produced siderophores and showed some hemolytic activity. The largest differences between strains were observed in the assimilation of hydrocarbons and human neurotransmitters. We show for the first time that fungi from the order Dothideales can assimilate aromatic hydrocarbons. To explain the background, we sequenced the genomes of all 49 strains and identified genes putatively involved in siderophore production and hemolysis. Genomic analysis revealed a fairly structured population of A.melanogenum, raising the possibility that some phylogenetic lineages have higher virulence potential than others. Population genomics indicated that the species is strictly clonal, although more than half of the genomes were diploid. The existence of relatively heterozygous diploids in an otherwise clonal species is described for only the second time in fungi. The genomic and phenotypic data from this study should help to resolve the non-trivial taxonomy of the genus Aureobasidium and reduce the medical hazards of exploiting the biotechnological potential of other, non-pathogenic species of this genus.

Cryptococcus extracellular vesicles properties and their use as vaccine platforms

Whereas extracellular vesicle (EV) research has become commonplace in different biomedical fields, this field of research is still in its infancy in mycology. Here we provide a robust set of data regarding the structural and compositional aspects of EVs isolated from the fungal pathogenic species Cryptococcus neoformans, C. deneoformans and C. deuterogattii. Using cutting-edge methodological approaches including cryogenic electron microscopy and cryogenic electron tomography, proteomics, and flow cytometry, we revisited cryptococcal EV features and suggest a new EV structural model, in which the vesicular lipid bilayer is covered by mannoprotein-based fibrillar decoration, bearing the capsule polysaccharide as its outer layer. About 10% of the EV population is devoid of fibrillar decoration, adding another aspect to EV diversity. By analysing EV protein cargo from the three species, we characterized the typical Cryptococcus EV proteome. It contains several membrane-bound protein families, including some Tsh proteins bearing a SUR7/PalI motif. The presence of known protective antigens on the surface of Cryptococcus EVs, resembling the morphology of encapsulated virus structures, suggested their potential as a vaccine. Indeed, mice immunized with EVs obtained from an acapsular C. neoformans mutant strain rendered a strong antibody response in mice and significantly prolonged their survival upon C. neoformans infection.

Genome-Wide Analysis of Nutrient Signaling Pathways Conserved in Arbuscular Mycorrhizal Fungi

Arbuscular mycorrhizal (AM) fungi form a mutualistic symbiosis with a majority of terrestrial vascular plants. To achieve an efficient nutrient trade with their hosts, AM fungi sense external and internal nutrients, and integrate different hierarchic regulations to optimize nutrient acquisition and homeostasis during mycorrhization. However, the underlying molecular networks in AM fungi orchestrating the nutrient sensing and signaling remain elusive. Based on homology search, we here found that at least 72 gene components involved in four nutrient sensing and signaling pathways, including cAMP-dependent protein kinase A (cAMP-PKA), sucrose non-fermenting 1 (SNF1) protein kinase, target of rapamycin kinase (TOR) and phosphate (PHO) signaling cascades, are well conserved in AM fungi. Based on the knowledge known in model yeast and filamentous fungi, we outlined the possible gene networks functioning in AM fungi. These pathways may regulate the expression of downstream genes involved in nutrient transport, lipid metabolism, trehalase activity, stress resistance and autophagy. The RNA-seq analysis and qRT-PCR results of some core genes further indicate that these pathways may play important roles in spore germination, appressorium formation, arbuscule longevity and sporulation of AM fungi. We hope to inspire further studies on the roles of these candidate genes involved in these nutrient sensing and signaling pathways in AM fungi and AM symbiosis.

A mutation in C. neoformans mitochondrial NADH dehydrogenase results in increased virulence in mice

Cryptococcus neoformans: (H99W) was serially passaged in the invertebrate wax moth Galleria mellonella fifteen times to study how fungal virulence evolves under selection and whether those adaptations affect virulence. The G. mellonella passaged strain (P15) and the pre-passage H99W strains were used to infect three different host models of C. neoformans: C. elegans, G. mellonella, and Balb/c mice. While there was no difference in survival in the invertebrate models, P15 killed mice faster than H99W through both intratracheal and intravenous routes of infection and mice infected intravenously with P15 showed higher fungal burden in the brain. Characterization of the major virulence factors of C. neoformans found that P15 had increased capsule size, GXM release, and melanization. Whole genome sequencing of P15 and H99W revealed two mutations in P15, an insertion in the promoter region of NADH dehydrogenase (CNAG_09000) and an insertion in the LMP1 gene (CNAG_06765). Both ATP production and metabolic rate were higher in P15 compared to H99W. Quantitative RT-PCR suggested that the increased ATP was due to increased RNA levels of NADH dehydrogenase. Thus, adaptation to growth in hemocytes resulted in increased production of ATP, increased metabolic rate, and increased virulence in mice. This was likely due to differential expression of virulence factors, which skewed the host immune response to a less efficient Th2 response, with higher levels of IL-4, IL-10, and TNF-α in the brain. Overall, serial passage experiments have increased our understanding of how this yeast evolves under innate immune selection pressure.

CNS-infecting pathogens Escherichia coli and Cryptococcus neoformans exploit the host Pdlim2 for intracellular traversal and exocytosis in the blood-brain barrier

Microbial penetration of the blood-brain barrier, a prerequisite for development of the central nervous system (CNS) infection, involves microbial invasion, intracellular traversal and exocytosis. Microbial invasion of the blood-brain barrier has been investigated, but the molecular basis for microbial traversal and exit from the blood-brain barrier remains unknown. We performed transcriptome analysis of the human brain microvascular endothelial cell (HBMEC) infected with Escherichia coli and Cryptococcus neoformans, representative bacterial and fungal pathogens common in CNS infection. Among the upregulated targets in response to E. coli and C. neoformans infection, PDLIM2 was knocked down by shRNA in HBMEC for further investigation. We demonstrated that Pdlim2 specifically regulated the microbial traversal and exit from HBMEC by assessing microbial invasion, transcytosis, intracellular multiplication and egression. Additionally, the defective exocytosis of internalized E. coli from the PDLIM2 shRNA knockdown cell was restored by treatment with a calcium ionophore (ionomycin). Moreover, we performed the proximity-dependent biotin labeling with the biotin ligase BioID2 and identified 210 potential Pdlim2-interactors. Among the nine enriched Pdlim2-interactors in response to both E. coli and C. neoformans infection, we selected MPRIP and showed that HBMEC with knockdown of MPRIP mimicked the phenotype of PDLIM2 knockdown cell. These results suggest that the CNS-infecting microbes hijack Pdlim2 and Mprip for intracellular traversal and exocytosis in the blood-brain barrier.

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

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

The evolving species concepts used for yeasts: from phenotypes and genomes to speciation networks

Here we review how evolving species concepts have been applied to understand yeast diversity. Initially, a phenotypic species concept was utilized taking into consideration morphological aspects of colonies and cells, and growth profiles. Later the biological species concept was added, which applied data from mating experiments. Biophysical measurements of DNA similarity between isolates were an early measure that became more broadly applied with the advent of sequencing technology, leading to a sequence-based species concept using comparisons of parts of the ribosomal DNA. At present phylogenetic species concepts that employ sequence data of rDNA and other genes are universally applied in fungal taxonomy, including yeasts, because various studies revealed a relatively good correlation between the biological species concept and sequence divergence. The application of genome information is becoming increasingly common, and we strongly recommend the use of complete, rather than draft genomes to improve our understanding of species and their genome and genetic dynamics. Complete genomes allow in-depth comparisons on the evolvability of genomes and, consequently, of the species to which they belong. Hybridization seems a relatively common phenomenon and has been observed in all major fungal lineages that contain yeasts. Note that hybrids may greatly differ in their post-hybridization development. Future in-depth studies, initially using some model species or complexes may shift the traditional species concept as isolated clusters of genetically compatible isolates to a cohesive speciation network in which such clusters are interconnected by genetic processes, such as hybridization.

Proteomic analysis of serial isolates of Trichosporon asahii identifies host-specific adaptations using the TMT/MRM approach

The opportunistic fungal pathogen Trichosporon asahii (T. asahii) is an important causal agent of mortality in immunocompromised patients and associated with frequent relapses, even with sufficient antifungal treatment. Investigating the proteomes of initial and recurrent isolates may help to identify within-host adaptive changes. In this study, using tandem mass tag (TMT)-labeling combined with liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS) technology, we analyzed the proteomes of two T. asahii strains that were isolated 15 years apart from the same patient who suffered initial and recurrent episodes of systemic disseminated trichosporonosis. A total of 597 differentially expressed proteins were identified. Functional analysis showed that the increased proteins were primarily concentrated on peptide/protein/energy/drug metabolism and translation. Most of the results were determined to be consistent with the findings of phenotypic assays, such as tests for drug susceptibility, temperature growth, biofilm formation, melanization and paromomycin assays. Moreover, we performed multiple reaction monitoring (MRM) mass spectrometry to verify 27 candidate proteins, and the results of this experiment were also highly consistent with the results of the TMT analysis. Therefore, to the best of our knowledge, these data provide the first molecular evidence of how the T. asahii proteome changes related to host-specific adaptation during human infection. SIGNIFICANCE: Systemic infection with Trichosporon asahii (T. asahii) has recently been recognized as an important causal agent of mortality in immunocompromised patients. Although triazole treatment usually works efficiently in the early phase of infection, many patients relapse. Hence, comparative analyses of the proteomics of initial and recurrent isolates may reveal evidence of adaptive changes within the host. Our study demonstrates that the recurrent strain has undergone proteomic changes using tandem mass tag (TMT)-labeling combined with liquid chromatography-mass spectrometry/mass spectrometry (LC-MS/MS). Moreover, the results of phenotypic assays, including drug susceptibility, temperature growth, biofilm formation, melanization and paromomycin assays, were highly consistent with the proteomic changes, and multiple reaction monitoring (MRM) verification also showed similar trends to the TMT results. In summary, our study is the first to investigate the adaptation of T. asahii under pressure from antifungal chemotherapy and host immune responses.

Population genomic analysis of Cryptococcus Brazilian isolates reveals an African type subclade distribution

The genomes of a large number of Cryptococcus neoformans isolates have been sequenced and analyzed in recent years. These genomes have been used to understand the global population structure of this opportunistic pathogen. However, only a small number of South American isolates have been considered in these studies, and the population structure of C. neoformans in this part of the world remains elusive. Here, we analyzed the genomic sequences of 53 Brazilian Cryptococcus isolates and deciphered the C. neoformans population structure in this country. Our data reveal an African-like structure that suggested repeated intercontinental transports from Africa to South America. We also identified a mutator phenotype in one VNBII Brazilian isolate, exemplifying how fast-evolving isolates can shape the Cryptococcus population structure. Finally, phenotypic analyses revealed wide diversity but not lineage specificity in the expression of classical virulence traits within the set of isolates.

Moderate levels of 5-fluorocytosine cause the emergence of high frequency resistance in cryptococci

The antifungal agent 5-fluorocytosine (5-FC) is used for the treatment of several mycoses, but is unsuitable for monotherapy due to the rapid development of resistance. Here, we show that cryptococci develop resistance to 5-FC at a high frequency when exposed to concentrations several fold above the minimal inhibitory concentration. The genomes of resistant clones contain alterations in genes relevant as well as irrelevant for 5-FC resistance, suggesting that 5-FC may be mutagenic at moderate concentrations. Mutations in FCY2 (encoding a known permease for 5-FC uptake), FCY1, FUR1, UXS1 (encoding an enzyme that converts UDP-glucuronic acid to UDP-xylose) and URA6 contribute to 5-FC resistance. The uxs1 mutants accumulate UDP-glucuronic acid, which appears to down-regulate expression of permease FCY2 and reduce cellular uptake of the drug. Additional mutations in genes known to be required for UDP-glucuronic acid synthesis (UGD1) or a transcriptional factor NRG1 suppress UDP-glucuronic acid accumulation and 5-FC resistance in the uxs1 mutants.

Pathogenic fungi rapidly develop resistance to the antifungal agent 5-fluorocytosine (5-FC). Here, Chang et al. explore the mechanisms by which Cryptococcus develops 5-FC resistance at a high frequency, including mutations in several genes and altered levels of key metabolites.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Molecular Markers Reveal Epidemiological Patterns and Evolutionary Histories of the Human Pathogenic Cryptococcus

The human pathogenic Cryptococcus species are the main agents of fungal meningitis in humans and the causes of other diseases collectively called cryptococcosis. There are at least eight evolutionary divergent lineages among these agents, with different lineages showing different geographic and/or ecological distributions. In this review, we describe the main strain typing methods that have been used to analyze the human pathogenic Cryptococcus and discuss how molecular markers derived from the various strain typing methods have impacted our understanding of not only cryptococcal epidemiology but also its evolutionary histories. These methods include serotyping, multilocus enzyme electrophoresis, electrophoretic karyotyping, random amplified polymorphic DNA, restriction fragment length polymorphism, PCR-fingerprinting, amplified fragment length polymorphism, multilocus microsatellite typing, single locus and multilocus sequence typing, matrix-assisted laser desorption/ionization time of flight mass spectrometry, and whole genome sequencing. The major findings and the advantages and disadvantages of each method are discussed. Together, while controversies remain, these strain typing methods have helped reveal (i) the broad phylogenetic pattern among these agents, (ii) the centers of origins for several lineages and their dispersal patterns, (iii) the distributions of genetic variation among geographic regions and ecological niches, (iv) recent hybridization among several lineages, and (v) specific mutations during infections within individual patients. However, significant challenges remain. Multilocus sequence typing and whole genome sequencing are emerging as the gold standards for continued strain typing and epidemiological investigations of cryptococcosis.

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

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

Transcriptomic response of Cryptococcus neoformans to ecologically relevant nitrogen concentrations

Nitrogen availability is vital for the growth and survival of Cryptococcus neoformans in the natural environment. Two major ecological reservoirs were previously described for C. neoformans, namely, pigeon guano and the woody debris of various tree species. In contrast to the abundance of available nitrogen in guano, C. neoformans must adapt to severely limited nitrogen conditions within arboreal ecological niches. Previously, we demonstrated the role of nitrogen limitation in the production of cryptococcal virulence factors and drug tolerance. The genetic response underlying this adaptation to nitrogen deficiency, however, remains to be determined. Therefore, in the present study we investigated the transcriptomic response of C. neoformans to ecologically relevant nitrogen concentrations using RNA-sequencing. Our data revealed that low nitrogen conditions modulate the expression of numerous virulence genes in C. neoformans. Among these were, CTR4 and CGP1, which showed highly significant modulation under low nitrogen conditions. Furthermore, data analysis revealed the upregulation of antifungal tolerance-related genes in low nitrogen conditions, including genes involved in ergosterol biosynthetic processes and cell wall integrity. Overall, our findings provide insight into the survival of C. neoformans in nitrogen-poor ecological niches and suggest that pre-adaptation to these conditions may influence the pathobiology of this yeast.

Expression Patterns in Reductive Iron Assimilation and Functional Consequences during Phagocytosis of Lichtheimia corymbifera, an Emerging Cause of Mucormycosis

Iron is an essential micronutrient for most organisms and fungi are no exception. Iron uptake by fungi is facilitated by receptor-mediated internalization of siderophores, heme and reductive iron assimilation (RIA). The RIA employs three protein groups: (i) the ferric reductases (Fre5 proteins), (ii) the multicopper ferroxidases (Fet3) and (iii) the high-affinity iron permeases (Ftr1). Phenotyping under different iron concentrations revealed detrimental effects on spore swelling and hyphal formation under iron depletion, but yeast-like morphology under iron excess. Since access to iron is limited during pathogenesis, pathogens are placed under stress due to nutrient limitations. To combat this, gene duplication and differential gene expression of key iron uptake genes are utilized to acquire iron against the deleterious effects of iron depletion. In the genome of the human pathogenic fungus L. corymbifera, three, four and three copies were identified for FRE5, FTR1 and FET3 genes, respectively. As in other fungi, FET3 and FTR1 are syntenic and co-expressed in L. corymbifera. Expression of FRE5, FTR1 and FET3 genes is highly up-regulated during iron limitation (Fe-), but lower during iron excess (Fe+). Fe- dependent upregulation of gene expression takes place in LcFRE5 II and III, LcFTR1 I and II, as well as LcFET3 I and II suggesting a functional role in pathogenesis. The syntenic LcFTR1 I–LcFET3 I gene pair is co-expressed during germination, whereas LcFTR1 II- LcFET3 II is co-expressed during hyphal proliferation. LcFTR1 I, II and IV were overexpressed in Saccharomyces cerevisiae to represent high and moderate expression of intracellular transport of Fe3+, respectively. Challenge of macrophages with the yeast mutants revealed no obvious role for LcFTR1 I, but possible functions of LcFTR1 II and IVs in recognition by macrophages. RIA expression pattern was used for a new model of interaction between L. corymbifera and macrophages.

The interplay of phenotype and genotype in Cryptococcus neoformans disease

Cryptococcus neoformans is an opportunistic fungal pathogen that causes life-threatening meningitis primarily in immunocompromised individuals. In order to survive and proliferate during infection, C. neoformans must adapt to a variety of stresses it encounters within the host. Patient outcome depends on the interaction between the pathogen and the host. Understanding the mechanisms that C. neoformans uses to facilitate adaptation to the host and promote pathogenesis is necessary to better predict disease severity and establish proper treatment. Several virulence phenotypes have been characterized in C. neoformans, but the field still lacks a complete understanding of how genotype and phenotype contribute to clinical outcome. Furthermore, while it is known that C. neoformans genotype impacts patient outcome, the mechanisms remain unknown. This lack of understanding may be due to the genetic heterogeneity of C. neoformans and the extensive phenotypic variation observed between and within isolates during infection. In this review, we summarize the current understanding of how the various genotypes and phenotypes observed in C. neoformans correlate with human disease progression in the context of patient outcome and recurrence. We also postulate the mechanisms underlying the genetic and phenotypic changes that occur in vivo to promote rapid adaptation in the host.

The Transcription Factor Pdr802 Regulates Titan Cell Formation and Pathogenicity of Cryptococcus neoformans

The pathogenic yeast Cryptococcus neoformans presents a worldwide threat to human health, especially in the context of immunocompromise, and current antifungal therapy is hindered by cost, limited availability, and inadequate efficacy. After the infectious particle is inhaled, C. neoformans initiates a complex transcriptional program that integrates cellular responses and enables adaptation to the host lung environment.

Cryptococcus neoformans is a ubiquitous, opportunistic fungal pathogen that kills almost 200,000 people worldwide each year. It is acquired when mammalian hosts inhale the infectious propagules; these are deposited in the lung and, in the context of immunocompromise, may disseminate to the brain and cause lethal meningoencephalitis. Once inside the host, C. neoformans undergoes a variety of adaptive processes, including secretion of virulence factors, expansion of a polysaccharide capsule that impedes phagocytosis, and the production of giant (Titan) cells. The transcription factor Pdr802 is one regulator of these responses to the host environment. Expression of the corresponding gene is highly induced under host-like conditions in vitro and is critical for C. neoformans dissemination and virulence in a mouse model of infection. Direct targets of Pdr802 include the quorum sensing proteins Pqp1, Opt1, and Liv3; the transcription factors Stb4, Zfc3, and Bzp4, which regulate cryptococcal brain infectivity and capsule thickness; the calcineurin targets Had1 and Crz1, important for cell wall remodeling and C. neoformans virulence; and additional genes related to resistance to host temperature and oxidative stress, and to urease activity. Notably, cryptococci engineered to lack Pdr802 showed a dramatic increase in Titan cells, which are not phagocytosed and have diminished ability to directly cross biological barriers. This explains the limited dissemination of pdr802 mutant cells to the central nervous system and the consequently reduced virulence of this strain. The role of Pdr802 as a negative regulator of Titan cell formation is thus critical for cryptococcal pathogenicity.

Genetic contribution to high temperature tolerance in Cryptococcus neoformans

The human fungal pathogen Cryptococcus neoformans relies on a complex signaling network for the adaptation and survival at the host temperature. Protein phosphatase calcineurin is central to proliferation at 37°C but its exact contributions remain ill-defined. To better define genetic contributions to the C. neoformans temperature tolerance, 4031 gene knockouts were screened for genes essential at 37°C and under conditions that keep calcineurin inactive. Identified 83 candidate strains, potentially sensitive to 37°C, were subsequently subject to technologically simple yet robust assay, in which cells are exposed to a temperature gradient. This has resulted in identification of 46 genes contributing to the maximum temperature at which C. neoformans can proliferate (Tmax). The 46 mutants, characterized by a range of Tmax on drug-free media, were further assessed for Tmax under conditions that inhibit calcineurin, which led to identification of several previously uncharacterized knockouts exhibiting synthetic interaction with the inhibition of calcineurin. A mutant that lacked septin Cdc11 was among those with the lowest Tmax and failed to proliferate in the absence of calcineurin activity. To further define connections with calcineurin and the role for septins in high temperature growth, the 46 mutants were tested for cell morphology at 37°C and growth in the presence of agents disrupting cell wall and cell membrane. Mutants sensitive to calcineurin inhibition were tested for synthetic lethal interaction with deletion of the septin-encoding CDC12 and the localization of the septin Cdc3-mCherry. The analysis described here pointed to previously uncharacterized genes that were missed in standard growth assays indicating that the temperature gradient assay is a valuable complementary tool for elucidating the genetic basis of temperature range at which microorganisms proliferate.

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

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

A conserved Zn2Cys6 transcription factor, identified in a spontaneous mutant from in vitro passaging, is involved in pathogenicity of the blackleg fungus Leptosphaeria maculans

Continuous passaging in vitro can lead to the accumulation of changes in DNA sequence that potentially affect the properties of microbes, making them different from the original isolates. The identification of such genetic alterations is rare in fungi. A set of insertional mutants in the plant pathogenic fungus Leptosphaeria maculans, all derived from the same transformation experiment, had independent Agrobacterium T-DNA insertions and reduced pathogenicity on canola (Brassica napus). None of the insertions co-segregated in progeny from crosses with the reduction in pathogenicity. Genome sequences of three strains were analysed, and a mutation identified in a gene (ptf1, for pathogenicity-associated transcription factor 1) encoding a putative Zn₂(II)Cys₆ transcription factor. Homologs are found in other ascomycetes, and are required for pathogenicity by Fusarium graminearum, Fusarium oxysporum and Magnaporthe oryzae. The mutation in the L. maculans ptf1 gene co-segregates in progeny from crosses with the reduction in pathogenicity, a strain with an independent mutant allele isolated using CRISPR-Cas9 editing has reduced pathogenicity, and addition of wild type copies of the gene restores pathogenicity. Thus, this work defines a base pair substitution that occurred during in vitro passaging of a fungus that contributed to an attenuation of pathogenicity.

Transcription factor-driven alternative localization of Cryptococcus neoformans superoxide dismutase

Cryptococcus neoformans is an opportunistic fungal pathogen whose pathogenic lifestyle is linked to its ability to cope with fluctuating levels of copper (Cu), an essential metal involved in multiple virulence mechanisms, within distinct host niches. During lethal cryptococcal meningitis in the brain, C. neoformans senses a Cu-deficient environment and is highly dependent on its ability to scavenge trace levels of Cu from its host and adapt to Cu scarcity to successfully colonize this niche. In this study, we demonstrate for this critical adaptation, the Cu-sensing transcription factor Cuf1 differentially regulates the expression of the SOD1 and SOD2 superoxide dismutases in novel ways. Genetic and transcriptional analysis reveals Cuf1 specifies 5’-truncations of the SOD1 and SOD2 mRNAs through specific binding to Cu responsive elements within their respective promoter regions. This results in Cuf1-dependent repression of the highly abundant SOD1 and simultaneously induces expression of two isoforms of SOD2, the canonical mitochondrial targeted isoform and a novel alternative cytosolic isoform, from a single alternative transcript produced specifically under Cu limitation. The generation of cytosolic Sod2 during Cu limitation is required to maintain cellular antioxidant defense against superoxide stress both in vitro and in vivo. Further, decoupling Cuf1 regulation of Sod2 localization compromises the ability of C. neoformans to colonize organs in murine models of cryptococcosis. Our results provide a link between transcription factor–mediated alteration of protein localization and cell proliferation under stress, which could impact tissue colonization by a fungal pathogen.

Factors enforcing the species boundary between the human pathogens Cryptococcus neoformans and Cryptococcus deneoformans

Hybridization has resulted in the origin and variation in extant species, and hybrids continue to arise despite pre- and post-zygotic barriers that limit their formation and evolutionary success. One important system that maintains species boundaries in prokaryotes and eukaryotes is the mismatch repair pathway, which blocks recombination between divergent DNA sequences. Previous studies illuminated the role of the mismatch repair component Msh2 in blocking genetic recombination between divergent DNA during meiosis. Loss of Msh2 results in increased interspecific genetic recombination in bacterial and yeast models, and increased viability of progeny derived from yeast hybrid crosses. Hybrid isolates of two pathogenic fungal Cryptococcus species, Cryptococcus neoformans and Cryptococcus deneoformans, are isolated regularly from both clinical and environmental sources. In the present study, we sought to determine if loss of Msh2 would relax the species boundary between C. neoformans and C. deneoformans. We found that crosses between these two species in which both parents lack Msh2 produced hybrid progeny with increased viability and high levels of aneuploidy. Whole-genome sequencing revealed few instances of recombination among hybrid progeny and did not identify increased levels of recombination in progeny derived from parents lacking Msh2. Several hybrid progeny produced structures associated with sexual reproduction when incubated alone on nutrient-rich medium in light, a novel phenotype in Cryptococcus. These findings represent a unique, unexpected case where rendering the mismatch repair system defective did not result in increased meiotic recombination across a species boundary. This suggests that alternative pathways or other mismatch repair components limit meiotic recombination between homeologous DNA and enforce species boundaries in the basidiomycete Cryptococcus species.

Several mechanisms enforce species boundaries by either preventing the formation of hybrid zygotes, known as pre-zygotic barriers, or preventing the viability and fecundity of hybrids, known as post-zygotic barriers. Despite these barriers, interspecific hybrids form at an appreciable frequency, such as hybrid isolates of the human fungal pathogenic species, Cryptococcus neoformans and Cryptococcus deneoformans, which are regularly isolated from both clinical and environmental sources. C. neoformans x C. deneoformans hybrids are typically highly aneuploid, sterile, and display phenotypes intermediate to those of either parent, although self-fertile isolates and transgressive phenotypes have been observed. One important mechanism known to enforce species boundaries or lead to incipient speciation is the DNA mismatch repair system, which blocks recombination between divergent DNA sequences during meiosis. The aim of this study was to determine if genetically deleting the DNA mismatch repair component Msh2 would relax the species boundary between C. neoformans and C. deneoformans. Progeny derived from C. neoformans x C. deneoformans crosses in which both parental strains lacked Msh2 had higher viability, and unlike previous studies in Saccharomyces, these Cryptococcus hybrid progeny had higher levels of aneuploidy and no observable increase in meiotic recombination at the whole-genome level.

Treatment strategies for cryptococcal infection: challenges, advances and future outlook

Cryptococcus spp., in particular Cryptococcus neoformans and Cryptococcus gattii, have an enormous impact on human health worldwide. The global burden of cryptococcal meningitis is almost a quarter of a million cases and 181,000 deaths annually, with mortality rates of 100% if infections remain untreated. Despite these alarming statistics, treatment options for cryptococcosis remain limited, with only three major classes of drugs approved for clinical use. Exacerbating the public health burden is the fact that the only new class of antifungal drugs developed in decades, the echinocandins, displays negligible antifungal activity against Cryptococcus spp., and the efficacy of the remaining therapeutics is hampered by host toxicity and pathogen resistance. Here, we describe the current arsenal of antifungal agents and the treatment strategies employed to manage cryptococcal disease. We further elaborate on the recent advances in our understanding of the intrinsic and adaptive resistance mechanisms that are utilized by Cryptococcus spp. to evade therapeutic treatments. Finally, we review potential therapeutic strategies, including combination therapy, the targeting of virulence traits, impairing stress response pathways and modulating host immunity, to effectively treat infections caused by Cryptococcus spp. Overall, understanding of the mechanisms that regulate anti-cryptococcal drug resistance, coupled with advances in genomics technologies and high-throughput screening methodologies, will catalyse innovation and accelerate antifungal drug discovery.

Pleiotropy and epistasis within and between signaling pathways defines the genetic architecture of fungal virulence

Cryptococcal disease is estimated to affect nearly a quarter of a million people annually. Environmental isolates of Cryptococcus deneoformans, which make up 15 to 30% of clinical infections in temperate climates such as Europe, vary in their pathogenicity, ranging from benign to hyper-virulent. Key traits that contribute to virulence, such as the production of the pigment melanin, an extracellular polysaccharide capsule, and the ability to grow at human body temperature have been identified, yet little is known about the genetic basis of variation in such traits. Here we investigate the genetic basis of melanization, capsule size, thermal tolerance, oxidative stress resistance, and antifungal drug sensitivity using quantitative trait locus (QTL) mapping in progeny derived from a cross between two divergent C. deneoformans strains. Using a "function-valued" QTL analysis framework that exploits both time-series information and growth differences across multiple environments, we identified QTL for each of these virulence traits and drug susceptibility. For three QTL we identified the underlying genes and nucleotide differences that govern variation in virulence traits. One of these genes, RIC8, which encodes a regulator of cAMP-PKA signaling, contributes to variation in four virulence traits: melanization, capsule size, thermal tolerance, and resistance to oxidative stress. Two major effect QTL for amphotericin B resistance map to the genes SSK1 and SSK2, which encode key components of the HOG pathway, a fungal-specific signal transduction network that orchestrates cellular responses to osmotic and other stresses. We also discovered complex epistatic interactions within and between genes in the HOG and cAMP-PKA pathways that regulate antifungal drug resistance and resistance to oxidative stress. Our findings advance the understanding of virulence traits among diverse lineages of Cryptococcus, and highlight the role of genetic variation in key stress-responsive signaling pathways as a major contributor to phenotypic variation.

Comparative Analysis of Carbohydrate Active Enzymes in the Flammulina velutipes var. lupinicola Genome

The purpose of this study was to determine the genome sequence of Flammulina velutipes var. lupinicola based on next-generation sequencing (NGS) and to identify the genes encoding carbohydrate-active enzymes (CAZymes) in the genome. The optimal assembly (71 kmer) based on ABySS de novo assembly revealed a total length of 33,223,357 bp (49.53% GC content). A total of 15,337 gene structures were identified in the F. velutipes var. lupinicola genome using ab initio gene prediction method with Funannotate pipeline. Analysis of the orthologs revealed that 11,966 (96.6%) out of the 15,337 predicted genes belonged to the orthogroups and 170 genes were specific for F. velutipes var. lupinicola. CAZymes are divided into six classes: auxiliary activities (AAs), glycosyltransferases (GTs), carbohydrate esterases (CEs), polysaccharide lyases (PLs), glycoside hydrolases (GHs), and carbohydrate-binding modules (CBMs). A total of 551 genes encoding CAZymes were identified in the F. velutipes var. lupinicola genome by analyzing the dbCAN meta server database (HMMER, Hotpep, and DIAMOND searches), which consisted of 54-95 AAs, 145-188 GHs, 55-73 GTs, 6-19 PLs, 13-59 CEs, and 7-67 CBMs. CAZymes can be widely used to produce bio-based products (food, paper, textiles, animal feed, and biofuels). Therefore, information about the CAZyme repertoire of the F. velutipes var. lupinicola genome will help in understanding the lignocellulosic machinery and in-depth studies will provide opportunities for using this fungus for biotechnological and industrial applications.

Registered report protocol: Quantitative analysis of septin Cdc10-associated proteome in Cryptococcus neoformans

Cryptococcus neoformans is a pathogenic basidiomycetous yeast that primarily infects immunocompromised individuals. C. neoformans can thrive during infections due to its three main virulence-related characteristics: the ability to grow at host temperature (37°C), formation of carbohydrate capsule, and its ability to produce melanin. C. neoformans strains lacking septin proteins Cdc3 or Cdc12 are viable at 25°C; however, they fail to proliferate at 37°C and are avirulent in the murine model of infection. The basis of septin contribution to growth at host temperature remains unknown. Septins are a family of conserved filament-forming GTPases with roles in cytokinesis and morphogenesis. In the model organism Saccharomyces cerevisiae septins are essential. S. cerevisiae septins form a higher order complex at the mother-bud neck to scaffold over 80 proteins, including those involved in cell wall organization, cell polarity, and cell cycle control. In C. neoformans, septins also form a complex at the mother-bud neck but the septin interacting proteome in this species remains largely unknown. Moreover, it remains possible that septins play other roles important for high temperature stress that are independent of their established role in cytokinesis. Therefore, we propose to perform a global analysis of septin Cdc10 binding partners in C. neoformans, including those that are specific to high temperature stress. This analysis will shed light on the underlying mechanism of survival of this pathogenic yeast during infection and can potentially lead to the discovery of novel drug targets.

An intergenic "safe haven" region in Cryptococcus neoformans serotype D genomes

Cryptococcus neoformans is an opportunistic human fungal pathogen and serves as a model organism for studies of eukaryotic microbiology and microbial pathogenesis. C. neoformans species complex is classified into serotype A, serotype D, and AD hybrids, which are currently considered different subspecies. Different serotype strains display varied phenotypes, virulence, and gene regulation. Genetic investigation of important pathways is often performed in both serotype A and D reference strains in order to identify diversification or conservation of the interrogated signaling network. Many genetic tools have been developed for C. neoformans serotype A reference strain H99, including the gene free "safe haven" (SH) regions for DNA integration identified based on genomic features. However, no such a genomic safe haven region has been identified in serotype D strains. Here, capitalizing on the available genomic, transcriptomic, and chromatin data, we identified an intergenic region named as SH3 for the serotype D reference strains JEC21 and XL280. We also designed a sgRNA and a vector facilitating any alien gene integration into SH3 through a CRISPR-Cas9 system. We found that gene inserted in this region complemented the corresponding gene deletion mutant. Fluorescent reporter gene inserted in SH3 can also be expressed efficiently. Insertion in SH3 itself did not alter the expression of adjacent genes and did not affect the growth or mating of C. neoformans. Thus, SH3 provides a resource for genetic manipulations in serotype D strains and will facilitate comparative analyses of gene functions in this species complex. In addition, the incorporation of the multi-omic data in our selection of the safe haven region could help similar studies in other organisms.

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

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

Microbe Profile: Cryptococcus neoformans species complex

Cryptococcus neoformans is a lethal fungus disguised in a polysaccharide coat. It can remain dormant in the host for decades prior to reactivation, causing systemic cryptococcosis in humans and other mammals. Cryptococcus deploys a multitude of traits to adapt to and survive within the host, including immunosuppression, an ability to replicate intra- and extra-cellularly in phagocytes, changes in morphology and ploidy, a predilection to infect the CNS, and the capacity to utilize neurotransmitters and unique carbon sources available in the brain. These pathogenic strategies displayed by this fungus might have evolved through its interactions with microbial predators in the environment.

Serial Passage of Cryptococcus neoformans in Galleria mellonella Results in Increased Capsule and Intracellular Replication in Hemocytes, but Not Increased Resistance to Hydrogen Peroxide

To gain insight into how pathogens adapt to new hosts, Cryptococcus neoformans (H99W) was serially passaged in Galleria mellonella. The phenotypic characteristics of the passaged strain (P15) and H99W were evaluated. P15 grew faster in hemolymph than H99W, in vitro and in vivo, suggesting that adaptation had occurred. However, P15 was more susceptible to hydrogen peroxide in vitro, killed fewer mouse macrophages, and had less fungal burden in human ex vivo macrophages than H99W. Analysis of gene expression changes during Galleria infection showed only a few different genes involved in the reactive oxygen species response. As P15 sheds more GXM than H99W, P15 may have adapted by downregulating hemocyte hydrogen peroxide production, possibly through increased capsular glucuronoxylomannan (GXM) shedding. Hemocytes infected with P15 produced less hydrogen peroxide, and hydrogen peroxide production in response to GXM-shedding mutants was correlated with shed GXM. Histopathological examination of infected larvae showed increased numbers and sizes of immune nodules for P15 compared to H99W, suggesting an enhanced, but functionally defective, response to P15. These results could explain why this infection model does not always correlate with murine models. Overall, C. neoformans’ serial passage in G. mellonella resulted in a better understanding of how this yeast evolves under selection.

Intracellular Cryptococcus neoformans disrupts the transcriptome profile of M1- and M2-polarized host macrophages

Macrophages serve as a first line of defense against infection with the facultative intracellular pathogen, Cryptococcus neoformans (Cn). However, the ability of these innate phagocytic cells to destroy ingested Cn is strongly influenced by polarization state with classically (M1) activated macrophages better able to control cryptococcal infections than alternatively (M2) activated cells. While earlier studies have demonstrated that intracellular Cn minimally affects the expression of M1 and M2 markers, the impact on the broader transcriptome associated with these states remains unclear. To investigate this, an in vitro cell culture model of intracellular infection together with RNA sequencing-based transcriptome profiling was used to measure the impact of Cn infection on gene expression in both polarization states. The gene expression profile of both M1 and M2 cells was extensively altered to become more like naive (M0) macrophages. Gene ontology analysis suggested that this involved changes in the activity of the Janus kinase-signal transducers and activators of transcription (JAK-STAT), p53, and nuclear factor-κB (NF-κB) pathways. Analyses of the principle polarization markers at the protein-level also revealed discrepancies between the RNA- and protein-level responses. In contrast to earlier studies, intracellular Cn was found to increase protein levels of the M1 marker iNos. In addition, common gene expression changes were identified that occurred post-Cn infection, independent of polarization state. This included upregulation of the transcriptional co-regulator Cited1, which was also apparent at the protein level in M1-polarized macrophages. These changes constitute a transcriptional signature of macrophage Cn infection and provide new insights into how Cn impacts gene expression and the phenotype of host phagocytes.

The Role of Testosterone and Gibberellic Acid in the Melanization of Cryptococcus neoformans

Cryptococcus neoformans, a spore-producing pathogenic yeast, affects immunocompromised individuals causing meningoencephalitis. Once C. neoformans is introduced via the respiratory tract, it is engulfed by macrophages and other phagocytes. One of C. neoformans's primary virulence factors is the pigment melanin, which is formed in the cell wall and protects the yeast against UV radiation and oxidizing agents produced by macrophages during phagocytosis. To better understand the observed sex bias (3:1; male:female) in C. neoformans infections, the phenotype of various virulence factors was determined in the presence of exogenous sex hormones. C. neoformans melanized faster in the presence of testosterone than it did in the presence of estradiol. Using a combination of RNA sequencing analysis and ELISA results, we identified a growth hormone, gibberellic acid (GA), produced in C. neoformans that was highly upregulated in the presence of testosterone. A variety of knockout strains of genes involved in the GA biosynthesis pathway showed significantly reduced melanization in the presence of testosterone. Additionally, inhibitors of GA also reduced melanization in the presence of testosterone. Thus, these data suggest that the gibberellic biosynthesis pathway is involved in melanization in C. neoformans, and the melanization difference observed in the presence of testosterone may be due to increased production of GA, which may partly explain the sex bias observed in C. neoformans infections.

Identification of Genomewide Alternative Splicing Events in Sequential, Isogenic Clinical Isolates of Candida albicans Reveals a Novel Mechanism of Drug Resistance and Tolerance to Cellular Stresses

The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

Alternative splicing (AS)—a process by which a single gene gives rise to different protein isoforms in eukaryotes—has been implicated in many basic cellular processes, but little is known about its role in drug resistance and fungal pathogenesis. The most common human fungal pathogen, Candida albicans, has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Here, we report AS regulating drug resistance in C. albicans. Comparative RNA-sequencing of two different sets of sequential, isogenic azole-sensitive and -resistant isolates of C. albicans revealed differential expression of splice isoforms of 14 genes. One of these was the superoxide dismutase gene SOD3, which contains a single intron. The sod3Δ/Δ mutant was susceptible to the antifungals amphotericin B (AMB) and menadione (MND). While AMB susceptibility was rescued by overexpression of both the spliced and unspliced SOD3 isoforms, only the spliced isoform could overcome MND susceptibility, demonstrating the functional relevance of this splicing in developing drug resistance. Furthermore, unlike AMB, MND inhibits SOD3 splicing and acts as a splicing inhibitor. Consistent with these observations, MND exposure resulted in increased levels of unspliced SOD3 isoform that are unable to scavenge reactive oxygen species (ROS), resulting in increased drug susceptibility. Collectively, these observations suggest that AS is a novel mechanism for stress adaptation and overcoming drug susceptibility in C. albicans.

IMPORTANCE The emergence of resistance in Candida albicans, an opportunistic pathogen, against the commonly used antifungals is becoming a major obstacle in its treatment. The necessity to identify new drug targets demands fundamental insights into the mechanisms used by this organism to develop drug resistance. C. albicans has introns in 4 to 6% of its genes, the functions of which remain largely unknown. Using the RNA-sequencing data from isogenic pairs of azole-sensitive and -resistant isolates of C. albicans, here, we show how C. albicans uses modulations in mRNA splicing to overcome antifungal drug stress.

A Transcriptional Regulatory Map of Iron Homeostasis Reveals a New Control Circuit for Capsule Formation in Cryptococcus neoformans

Iron is essential for the growth of the human fungal pathogen Cryptococcus neoformans within the vertebrate host, and iron sensing contributes to the elaboration of key virulence factors, including the formation of the polysaccharide capsule. C. neoformans employs sophisticated iron acquisition and utilization systems governed by the transcription factors Cir1 and HapX. However, the details of the transcriptional regulatory networks that are governed by these transcription factors and connections to virulence remain to be defined. Here, we used chromatin immunoprecipitation followed by next-generation sequencing (ChIP-seq) and transcriptome analysis (RNA-seq) to identify genes directly regulated by Cir1 and/or HapX in response to iron availability. Overall, 40 and 100 genes were directly regulated by Cir1, and 171 and 12 genes were directly regulated by HapX, under iron-limited and replete conditions, respectively. More specifically, we found that Cir1 directly controls the expression of genes required for iron acquisition and metabolism, and indirectly governs capsule formation by regulating specific protein kinases, a regulatory connection not previously revealed. HapX regulates the genes responsible for iron-dependent pathways, particularly under iron-depleted conditions. By analyzing target genes directly bound by Cir1 and HapX, we predicted the binding motifs for the transcription factors and verified that the purified proteins bind these motifs in vitro Furthermore, several direct target genes were coordinately and reciprocally regulated by Cir1 and HapX, suggesting that these transcription factors play conserved roles in the response to iron availability. In addition, biochemical analyses revealed that Cir1 and HapX are iron-containing proteins, implying that the regulatory networks of Cir1 and HapX may be influenced by the incorporation of iron into these proteins. Taken together, our identification of the genome-wide transcriptional networks provides a detailed understanding of the iron-related regulatory landscape, establishes a new connection between Cir1 and kinases that regulate capsule, and underpins genetic and biochemical analyses that reveal iron-sensing mechanisms for Cir1 and HapX in C. neoformans.

The Novel J-Domain Protein Mrj1 Is Required for Mitochondrial Respiration and Virulence in Cryptococcus neoformans

Cryptococcus neoformans is the causative agent of cryptococcal meningitis, a disease responsible for ∼15% of all HIV-related deaths. Unfortunately, development of antifungal drugs is challenging because potential targets are conserved between humans and C. neoformans. In this context, we characterized a unique J-domain protein, Mrj1, which lacks orthologs in humans. We showed that Mrj1 was required for normal mitochondrial respiration and that mutants lacking Mrj1 were deficient in growth, capsule elaboration, and virulence. Furthermore, we were able to phenocopy the defects in growth and capsule elaboration by inhibiting respiration. This result suggests that the role of Mrj1 in mitochondrial function was responsible for the observed virulence defects and reinforces the importance of mitochondria to fungal pathogenesis. Mitochondria are difficult to target, as their function is also key to human cells; however, Mrj1 presents an opportunity to target a unique fungal protein required for mitochondrial function and virulence in C. neoformans.

The opportunistic fungal pathogen Cryptococcus neoformans must adapt to the mammalian environment to establish an infection. Proteins facilitating adaptation to novel environments, such as chaperones, may be required for virulence. In this study, we identified a novel mitochondrial co-chaperone, Mrj1 (mitochondrial respiration J-domain protein 1), necessary for virulence in C. neoformans. The mrj1Δ and J-domain-inactivated mutants had general growth defects at both routine laboratory and human body temperatures and were deficient in the major virulence factor of capsule elaboration. The latter phenotype was associated with cell wall changes and increased capsular polysaccharide shedding. Accordingly, the mrj1Δ mutant was avirulent in a murine model of cryptococcosis. Mrj1 has a mitochondrial localization and co-immunoprecipitated with Qcr2, a core component of complex III of the electron transport chain. The mrj1 mutants were deficient in mitochondrial functions, including growth on alternative carbon sources, growth without iron, and mitochondrial polarization. They were also insensitive to complex III inhibitors and hypersensitive to an alternative oxidase (AOX) inhibitor, suggesting that Mrj1 functions in respiration. In support of this conclusion, mrj1 mutants also had elevated basal oxygen consumption rates which were completely abolished by the addition of the AOX inhibitor, confirming that Mrj1 is required for mitochondrial respiration through complexes III and IV. Furthermore, inhibition of complex III phenocopied the capsule and cell wall defects of the mrj1 mutants. Taken together, these results indicate that Mrj1 is required for normal mitochondrial respiration, a key aspect of adaptation to the host environment and virulence.

Epidemiological investigation for grouped cases of Trichosporon asahii using whole genome and IGS1 sequencing

BACKGROUND

Trichosporonosis is a rare invasive infection in humans mainly due to Trichosporon asahii, and especially recovered from patients having haematological malignancy. Since 2012, IGS1 region sequencing is used as a genotyping method to distinguish isolates, with high frequency of one haplotype worldwide and a geographic specificity for some haplotypes.

OBJECTIVES

We compared the IGS1 genotyping method and whole genome sequencing (WGS) to study the relationship between clinical isolates involved in two grouped cases in France.

METHODS

IGS1 sequencing and antifungal susceptibility testing were performed for 54 clinical isolates. Clinical data for 28 isolates included in surveillance programs were analysed. Whole genome was sequenced for 32 clinical isolates and the type strain.

RESULTS

All isolates were intrinsically resistant to flucytosine, while voriconazole had the most potent in vitro activity. The majority of the isolates was recovered from patients with haematological malignancies (42.86%), with a high proportion of children (<15 yrs-old, 32.14%) and a high mortality rate at three months (46.15%). Based on the WGS analysis, isolates exhibiting IGS1 haplotype 1, 3 and 7 belonged to different clades. Five isolates recovered during the first grouped cases had the same IGS1 haplotype and shared 99% of SNPs similarity. For the second grouped cases, four isolates had 98.7% of SNPs similarity while the isolate recovered 4 years earlier was totally unlinked.

CONCLUSIONS

We confirmed the usefulness of IGS1 sequencing for grouped cases infection of T. asahii. We underlined its limitation for the study of population structure and the utility of WGS analysis for the study of epidemiologically unrelated isolates.

Polymorphic centromere locations in the pathogenic yeast Candida parapsilosis

Centromeres pose an evolutionary paradox: strongly conserved in function but rapidly changing in sequence and structure. However, in the absence of damage, centromere locations are usually conserved within a species. We report here that isolates of the pathogenic yeast species Candida parapsilosis show within-species polymorphism for the location of centromeres on two of its eight chromosomes. Its old centromeres have an inverted-repeat (IR) structure, whereas its new centromeres have no obvious structural features but are located within 30 kb of the old site. Centromeres can therefore move naturally from one chromosomal site to another, apparently spontaneously and in the absence of any significant changes in DNA sequence. Our observations are consistent with a model in which all centromeres are genetically determined, such as by the presence of short or long IRs or by the ability to form cruciforms. We also find that centromeres have been hotspots for genomic rearrangements in the C. parapsilosis clade.

Synergistic and antagonistic drug interactions in the treatment of systemic fungal infections

Invasive fungal infections cause 1.6 million deaths annually, primarily in immunocompromised individuals. Mortality rates are as high as 90% due to limited treatments. The azole class antifungal, fluconazole, is widely available and has multi-species activity but only inhibits growth instead of killing fungal cells, necessitating long treatments. To improve treatment, we used our novel high-throughput method, the overlap² method (O2M) to identify drugs that interact with fluconazole, either increasing or decreasing efficacy. We identified 40 molecules that act synergistically (amplify activity) and 19 molecules that act antagonistically (decrease efficacy) when combined with fluconazole. We found that critical frontline beta-lactam antibiotics antagonize fluconazole activity. A promising fluconazole-synergizing anticholinergic drug, dicyclomine, increases fungal cell permeability and inhibits nutrient intake when combined with fluconazole. In vivo, this combination doubled the time-to-endpoint of mice with Cryptococcus neoformans meningitis. Thus, our ability to rapidly identify synergistic and antagonistic drug interactions can potentially alter the patient outcomes.

Individuals with weakened immune systems – such as recipients of organ transplants – can fall prey to illnesses caused by fungi that are harmless to most people. These infections are difficult to manage because few treatments exist to fight fungi, and many have severe side effects. Antifungal drugs usually slow the growth of fungi cells rather than kill them, which means that patients must remain under treatment for a long time, or even for life.

One way to boost efficiency and combat resistant infections is to combine antifungal treatments with drugs that work in complementary ways: the drugs strengthen each other’s actions, and together they can potentially kill the fungus rather than slow its progression. However, not all drug combinations are helpful. In fact, certain drugs may interact in ways that make treatment less effective. This is particularly concerning because people with weakened immune systems often take many types of medications.

Here, Wambaugh et al. harnessed a new high-throughput system to screen how 2,000 drugs (many of which already approved to treat other conditions) affected the efficiency of a common antifungal called fluconazole. This highlighted 19 drugs that made fluconazole less effective, some being antibiotics routinely used to treat patients with weakened immune systems.

On the other hand, 40 drugs boosted the efficiency of fluconazole, including dicyclomine, a compound currently used to treat inflammatory bowel syndrome. In fact, pairing dicyclomine and fluconazole more than doubled the survival rate of mice with severe fungal infections. The combined treatment could target many species of harmful fungi, even those that had become resistant to fluconazole alone.

The results by Wambaugh et al. point towards better treatments for individuals with serious fungal infections. Drugs already in circulation for other conditions could be used to boost the efficiency of fluconazole, while antibiotics that do not decrease the efficiency of this medication should be selected to treat at-risk patients.

Diverse transposable element landscapes in pathogenic and nonpathogenic yeast models: the value of a comparative perspective

Genomics and other large-scale analyses have drawn increasing attention to the potential impacts of transposable elements (TEs) on their host genomes. However, it remains challenging to transition from identifying potential roles to clearly demonstrating the level of impact TEs have on genome evolution and possible functions that they contribute to their host organisms. I summarize TE content and distribution in four well-characterized yeast model systems in this review: the pathogens Candida albicans and Cryptococcus neoformans, and the nonpathogenic species Saccharomyces cerevisiae and Schizosaccharomyces pombe. I compare and contrast their TE landscapes to their lifecycles, genomic features, as well as the presence and nature of RNA interference pathways in each species to highlight the valuable diversity represented by these models for functional studies of TEs. I then review the regulation and impacts of the Ty1 and Ty3 retrotransposons from Saccharomyces cerevisiae and Tf1 and Tf2 retrotransposons from Schizosaccharomyces pombe to emphasize parallels and distinctions between these well-studied elements. I propose that further characterization of TEs in the pathogenic yeasts would enable this set of four yeast species to become an excellent set of models for comparative functional studies to address outstanding questions about TE-host relationships.

Centromere deletion in Cryptococcus deuterogattii leads to neocentromere formation and chromosome fusions

The human fungal pathogen Cryptococcus deuterogattii is RNAi-deficient and lacks active transposons in its genome. C. deuterogattii has regional centromeres that contain only transposon relics. To investigate the impact of centromere loss on the C. deuterogattii genome, either centromere 9 or 10 was deleted. Deletion of either centromere resulted in neocentromere formation and interestingly, the genes covered by these neocentromeres maintained wild-type expression levels. In contrast to cen9∆ mutants, cen10∆ mutant strains exhibited growth defects and were aneuploid for chromosome 10. At an elevated growth temperature (37°C), the cen10∆ chromosome was found to have undergone fusion with another native chromosome in some isolates and this fusion restored wild-type growth. Following chromosomal fusion, the neocentromere was inactivated, and the native centromere of the fused chromosome served as the active centromere. The neocentromere formation and chromosomal fusion events observed in this study in C. deuterogattii may be similar to events that triggered genomic changes within the Cryptococcus/Kwoniella species complex and may contribute to speciation throughout the eukaryotic domain.

Centromere scission drives chromosome shuffling and reproductive isolation

A fundamental characteristic of eukaryotic organisms is the generation of genetic variation via sexual reproduction. Conversely, significant large-scale genome structure variations could hamper sexual reproduction, causing reproductive isolation and promoting speciation. The underlying processes behind large-scale genome rearrangements are not well understood and include chromosome translocations involving centromeres. Recent genomic studies in the Cryptococcus species complex revealed that chromosome translocations generated via centromere recombination have reshaped the genomes of different species. In this study, multiple DNA double-strand breaks (DSBs) were generated via the CRISPR/Cas9 system at centromere-specific retrotransposons in the human fungal pathogen Cryptococcus neoformans The resulting DSBs were repaired in a complex manner, leading to the formation of multiple interchromosomal rearrangements and new telomeres, similar to chromothripsis-like events. The newly generated strains harboring chromosome translocations exhibited normal vegetative growth but failed to undergo successful sexual reproduction with the parental wild-type strain. One of these strains failed to produce any spores, while another produced ∼3% viable progeny. The germinated progeny exhibited aneuploidy for multiple chromosomes and showed improved fertility with both parents. All chromosome translocation events were accompanied without any detectable change in gene sequences and thus suggest that chromosomal translocations alone may play an underappreciated role in the onset of reproductive isolation and speciation.

Pathogen and host genetics underpinning cryptococcal disease

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