CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair.

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

CryptoCEN: A Co-Expression Network for Cryptococcus neoformans reveals novel proteins involved in DNA damage repair.

Elucidating gene function is a major goal in biology, especially among non-model organisms. However, doing so is complicated by the fact that molecular conservation does not always mirror functional conservation, and that complex relationships among genes are responsible for encoding pathways and higher-order biological processes. Co-expression, a promising approach for predicting gene function, relies on the general principal that genes with similar expression patterns across multiple conditions will likely be involved in the same biological process. For Cryptococcus neoformans, a prevalent human fungal pathogen greatly diverged from model yeasts, approximately 60% of the predicted genes in the genome lack functional annotations. Here, we leveraged a large amount of publicly available transcriptomic data to generate a C. neoformans Co-Expression Network (CryptoCEN), successfully recapitulating known protein networks, predicting gene function, and enabling insights into the principles influencing co-expression. With 100% predictive accuracy, we used CryptoCEN to identify 13 new DNA damage response genes, underscoring the utility of guilt-by-association for determining gene function. Overall, co-expression is a powerful tool for uncovering gene function, and decreases the experimental tests needed to identify functions for currently under-annotated genes.

The wtf meiotic driver gene family has unexpectedly persisted for over 100 million years.

Meiotic drivers are selfish elements that bias their own transmission into more than half of the viable progeny produced by a driver+/driver- heterozygote. Meiotic drivers are thought to exist for relatively short evolutionary timespans because a driver gene or gene family is often found in a single species or in a group of very closely related species. Additionally, drivers are generally considered doomed to extinction when they spread to fixation or when suppressors arise. In this study, we examine the evolutionary history of the wtf meiotic drivers first discovered in the fission yeast Schizosaccharomyces pombe. We identify homologous genes in three other fission yeast species, S. octosporus, S. osmophilus, and S. cryophilus, which are estimated to have diverged over 100 million years ago from the S. pombe lineage. Synteny evidence supports that wtf genes were present in the common ancestor of these four species. Moreover, the ancestral genes were likely drivers as wtf genes in S. octosporus cause meiotic drive. Our findings indicate that meiotic drive systems can be maintained for long evolutionary timespans.

Genome-wide quantification of contributions to sexual fitness identifies genes required for spore viability and health in fission yeast.

Numerous genes required for sexual reproduction remain to be identified even in simple model species like Schizosaccharomyces pombe. To address this, we developed an assay in S. pombe that couples transposon mutagenesis with high-throughput sequencing (TN-seq) to quantitatively measure the fitness contribution of nonessential genes across the genome to sexual reproduction. This approach identified 532 genes that contribute to sex, including more than 200 that were not previously annotated to be involved in the process, of which more than 150 have orthologs in vertebrates. Among our verified hits was an uncharacterized gene, ifs1 (important for sex), that is required for spore viability. In two other hits, plb1 and alg9, we observed a novel mutant phenotype of poor spore health wherein viable spores are produced, but the spores exhibit low fitness and are rapidly outcompeted by wild type. Finally, we fortuitously discovered that a gene previously thought to be essential, sdg1 (social distancing gene), is instead required for growth at low cell densities and can be rescued by conditioned medium. Our assay will be valuable in further studies of sexual reproduction in S. pombe and identifies multiple candidate genes that could contribute to sexual reproduction in other eukaryotes, including humans.

Drug Resistance and Evolvability in an Emerging Human Fungal Pathogen.

Over the past decade, Candida auris has emerged as a highly transmissible human fungal pathogen. Because of its ability to transmit between patients in hospitals and its ability to rapidly develop drug resistance, C. auris presents unique challenges. However, at a genetic and genomic level we still understand relatively little about how drug resistance develops in this pathogen. Burrack et al. use experimental evolution and whole-genome sequencing to identify mutations correlated with fluconazole resistance in C. auris. They identify interesting genomic features, including highly plastic subtelomeric regions and whole chromosomal and segmental aneuploidies. Excitingly, they also identify the first example of a hypermutator strain in C. auris. In comparison with the model human fungal pathogen Candida albicans, C. auris is more likely to undergo mutation and less likely to undergo copy number variation in response to drug selection, which may be linked to differences in base ploidy level.

Diverse mating phenotypes impact the spread of wtf meiotic drivers in Schizosaccharomyces pombe.

Meiotic drivers are genetic elements that break Mendel's law of segregation to be transmitted into more than half of the offspring produced by a heterozygote. The success of a driver relies on outcrossing (mating between individuals from distinct lineages) because drivers gain their advantage in heterozygotes. It is, therefore, curious that Schizosaccharomyces pombe, a species reported to rarely outcross, harbors many meiotic drivers. To address this paradox, we measured mating phenotypes in S. pombe natural isolates. We found that the propensity for cells from distinct clonal lineages to mate varies between natural isolates and can be affected both by cell density and by the available sexual partners. Additionally, we found that the observed levels of preferential mating between cells from the same clonal lineage can slow, but not prevent, the spread of a wtf meiotic driver in the absence of additional fitness costs linked to the driver. These analyses reveal parameters critical to understanding the evolution of S. pombe and help explain the success of meiotic drivers in this species.

The fission yeast, Schizosaccharomyces pombe, is a haploid organism, meaning it has a single copy of each of its genes. S. pombe cells generally carry one copy of each chromosome and can reproduce clonally by duplicating these chromosomes and then dividing into two cells. However, when the yeast are starving, they can reproduce sexually. This involves two cells mating by fusing together to create a 'diploid zygote', which contains two copies of each gene. The zygote then undergoes 'meiosis', a special type of cell division in which the zygote first duplicates its genome and then divides twice. This results in four haploid spores which are analogous to sperm and eggs in humans that each contain one copy of the genome. The spores will grow and divide normally when conditions improve. The genes carried by each of the haploid spores depend on the cells that formed the zygote. If the two 'parent' yeast had the same version or 'allele' of a gene, all four spores will have it in their genome. However, if the two parents have different alleles, only 50% of the offspring will carry each version. Although this is usually the case, there are certain alleles, called meiotic drivers, that are transmitted to all offspring even in situations where it is only carried by one parent. Meiotic drivers can be found in many organisms, including mammals, but their behavior is easiest to study in yeast. Meiotic drivers known as killers achieve this by disposing of any 'sister' spores that do not inherit the same allele of this gene. This 'killing' can only happen when only one of the 'parents' carries the driver. This scenario is thought to rarely occur in species that inbreed, as inbreeding leads to both gene copies being the same. However, this does not appear to be the case for S. pombe, which contain a whole family of killer meiotic drivers, the wtf genes, despite also being reported to mainly inbreed. To investigate this contradiction, López Hernández et al. isolated several genetically distinct populations of S.pombe. These isolates were grown together to determine how often the each one would outcross (mate with an individual from a different population) or inbreed. The results found that levels of inbreeding varied between isolates. Next, López Hernández et al. used mathematical modelling and experimental evolution analyses to study how wtf drivers spread amongst these populations. This revealed that wtf genes spread faster in populations with more outcrossing. In some instances, the wtf driver was linked to a gene that could harm the population. In these cases, López Hernández et al. found than inbreeding could purge these drivers and stop them from spreading the dangerous alleles through the population. López Hernández et al. establish a simple experimental system to model driver evolution and experimentally demonstrate how key parameters, such as outcrossing rates, affect the spread of these genes. Understanding how meiotic drivers spread is important, as these systems could potentially be used to modify populations important to humans, such as crops or disease vectors.

A Novel Resistance Pathway for Calcineurin Inhibitors in the Human-Pathogenic Mucorales Mucor circinelloides.

Mucormycosis is an emerging lethal fungal infection in immunocompromised patients. Mucor circinelloides is a causal agent of mucormycosis and serves as a model system to understand genetics in Mucorales. Calcineurin is a conserved virulence factor in many pathogenic fungi, and calcineurin inhibition or deletion of the calcineurin regulatory subunit (CnbR) in Mucor results in a shift from hyphal to yeast growth. We analyzed 36 calcineurin inhibitor-resistant or bypass mutants that exhibited hyphal growth in the presence of calcineurin inhibitors or in the yeast-locked cnbRΔ mutant background without carrying any mutations in known calcineurin components. We found that a majority of the mutants had altered sequence in a gene, named here bycA (bypass of calcineurin). bycA encodes an amino acid permease. We verified that both the bycAΔ single mutant and the bycAΔ cnbRΔ double mutant are resistant to calcineurin inhibitor FK506, thereby demonstrating a novel mechanism of resistance against calcineurin inhibitors. We also found that the level of expression of bycA was significantly higher in the wild-type strain treated with FK506 and in the cnbRΔ mutants but was significantly lower in the wild-type strain without FK506 treatment. These findings suggest that bycA is a negative regulator of hyphal growth and/or a positive regulator of yeast growth in Mucor and that calcineurin suppresses expression of the bycA gene at the mRNA level to promote hyphal growth. BycA is involved in the Mucor hypha-yeast transition as our data demonstrate positive correlations among bycA expression, protein kinase A activity, and Mucor yeast growth. Also, calcineurin, independently of its role in morphogenesis, contributes to virulence traits, including phagosome maturation blockade, host cell damages, and proangiogenic growth factor induction during interactions with hosts.IMPORTANCEMucor is intrinsically resistant to most known antifungals, which makes mucormycosis treatment challenging. Calcineurin is a serine/threonine phosphatase that is widely conserved across eukaryotes. When calcineurin function is inhibited in Mucor, growth shifts to a less virulent yeast growth form, which makes calcineurin an attractive target for development of new antifungal drugs. Previously, we identified two distinct mechanisms through which Mucor can become resistant to calcineurin inhibitors involving Mendelian mutations in the gene for FKBP12, including mechanisms corresponding to calcineurin A or B subunits and epimutations silencing the FKBP12 gene. Here, we identified a third novel mechanism where loss-of-function mutations in the amino acid permease corresponding to the bycA gene contribute to resistance against calcineurin inhibitors. When calcineurin activity is absent, BycA can activate protein kinase A (PKA) to promote yeast growth via a cAMP-independent pathway. Our data also show that calcineurin activity contributes to host-pathogen interactions primarily in the pathogenesis of Mucor.

5-fluorocytosine resistance is associated with hypermutation and alterations in capsule biosynthesis in Cryptococcus.

Patients infected with the fungal pathogen Cryptococcus are most effectively treated with a combination of 5-fluorocytosine (5FC) and amphotericin B. 5FC acts as a prodrug, which is converted into toxic 5-fluorouracil (5FU) upon uptake into fungal cells. However, the pathogen frequently develops resistance through unclear mechanisms. Here we show that resistance to 5FC in Cryptococcus deuterogattii is acquired more frequently in isolates with defects in DNA mismatch repair that confer an elevated mutation rate. We use whole genome sequencing of 16 independent isolates to identify mutations associated with 5FC resistance in vitro. We find mutations in known resistance genes (FUR1 and FCY2) and in a gene UXS1, previously shown to encode an enzyme that converts UDP-glucuronic acid to UDP-xylose for capsule biosynthesis, but not known to play a role in 5FC metabolism. Mutations in UXS1 lead to accumulation of UDP-glucuronic acid and alterations in nucleotide metabolism, which appear to suppress toxicity of both 5FC and its toxic derivative 5FU.

Roles for Stress Response and Cell Wall Biosynthesis Pathways in Caspofungin Tolerance in Cryptococcus neoformans.

Limited antifungal diversity and availability are growing problems for the treatment of fungal infections in the face of increasing drug resistance. The echinocandins, one of the newest classes of antifungal drugs, inhibit production of a crucial cell wall component. However, these compounds do not effectively inhibit the growth of the opportunistic fungal pathogen Cryptococcus neoformans, despite potent inhibition of the target enzyme in vitro Therefore, we performed a forward genetic screen to identify cellular processes that mediate the relative tolerance of this organism to the echinocandin drug caspofungin. Through these studies, we identified 14 genetic mutants that enhance caspofungin antifungal activity. Rather than directly affecting caspofungin antifungal activity, these mutations seem to prevent the activation of various stress-induced compensatory cellular processes. For example, the pfa4Δ mutant has defects in the palmitoylation and localization of many of its target proteins, including the Ras1 GTPase and the Chs3 chitin synthase, which are both required for caspofungin tolerance. Similarly, we have confirmed the link between caspofungin treatment and calcineurin signaling in this organism, but we suggest a deeper mechanism in which caspofungin tolerance is mediated by multiple pathways downstream of calcineurin function. In summary, we describe here several pathways in C. neoformans that contribute to the complex caspofungin tolerance phenotype in this organism.

Gastrointestinal microbiota alteration induced by Mucor circinelloides in a murine model.

Mucor circinelloides is a pathogenic fungus and etiologic agent of mucormycosis. In 2013, cases of gastrointestinal illness after yogurt consumption were reported to the US FDA, and the producer found that its products were contaminated with Mucor. A previous study found that the Mucor strain isolated from an open contaminated yogurt exhibited virulence in a murine systemic infection model and showed that this strain is capable of surviving passage through the gastrointestinal tract of mice. In this study, we isolated another Mucor strain from an unopened yogurt that is closely related but distinct from the first Mucor strain and subsequently examined if Mucor alters the gut microbiota in a murine host model. DNA extracted from a ten-day course of stool samples was used to analyze the microbiota in the gastrointestinal tracts of mice exposed via ingestion of Mucor spores. The bacterial 16S rRNA gene and fungal ITS1 sequences obtained were used to identify taxa of each kingdom. Linear regressions revealed that there are changes in bacterial and fungal abundance in the gastrointestinal tracts of mice which ingested Mucor. Furthermore, we found an increased abundance of the bacterial genus Bacteroides and a decreased abundance of the bacteria Akkermansia muciniphila in the gastrointestinal tracts of exposed mice. Measurements of abundances show shifts in relative levels of multiple bacterial and fungal taxa between mouse groups. These findings suggest that exposure of the gastrointestinal tract to Mucor can alter the microbiota and, more importantly, illustrate an interaction between the intestinal mycobiota and bacteriota. In addition, Mucor was able to induce increased permeability in epithelial cell monolayers in vitro, which might be indicative of unstable intestinal barriers. Understanding how the gut microbiota is shaped is important to understand the basis of potential methods of treatment for gastrointestinal illness. How the gut microbiota changes in response to exposure, even by pathogens not considered to be causative agents of food-borne illness, may be important to how commercial food producers prevent and respond to contamination of products aimed at the public. This study provides evidence that the fungal microbiota, though understudied, may play an important role in diseases of the human gut.

Broad antifungal resistance mediated by RNAi-dependent epimutation in the basal human fungal pathogen Mucor circinelloides.

Mucormycosis-an emergent, deadly fungal infection-is difficult to treat, in part because the causative species demonstrate broad clinical antifungal resistance. However, the mechanisms underlying drug resistance in these infections remain poorly understood. Our previous work demonstrated that one major agent of mucormycosis, Mucor circinelloides, can develop resistance to the antifungal agents FK506 and rapamycin through a novel, transient RNA interference-dependent mechanism known as epimutation. Epimutations silence the drug target gene and are selected by drug exposure; the target gene is re-expressed and sensitivity is restored following passage without drug. This silencing process involves generation of small RNA (sRNA) against the target gene via core RNAi pathway proteins. To further elucidate the role of epimutation in the broad antifungal resistance of Mucor, epimutants were isolated that confer resistance to another antifungal agent, 5-fluoroorotic acid (5-FOA). We identified epimutant strains that exhibit resistance to 5-FOA without mutations in PyrF or PyrG, enzymes which convert 5-FOA into the active toxic form. Using sRNA hybridization as well as sRNA library analysis, we demonstrate that these epimutants harbor sRNA against either pyrF or pyrG, and further show that this sRNA is lost after reversion to drug sensitivity. We conclude that epimutation is a mechanism capable of targeting multiple genes, enabling Mucor to develop resistance to a variety of antifungal agents. Elucidation of the role of RNAi in epimutation affords a fuller understanding of mucormycosis. Furthermore, it improves our understanding of fungal pathogenesis and adaptation to stresses, including the evolution of drug resistance.

A High-Resolution Map of Meiotic Recombination in Cryptococcus deneoformans Demonstrates Decreased Recombination in Unisexual Reproduction.

Multiple species within the basidiomycete genus Cryptococcus cause cryptococcal disease. These species are estimated to affect nearly a quarter of a million people leading to ∼180,000 mortalities, annually. Sexual reproduction, which can occur between haploid yeasts of the same or opposite mating type, is a potentially important contributor to pathogenesis as recombination can generate novel genotypes and transgressive phenotypes. However, our quantitative understanding of recombination in this clinically important yeast is limited. Here, we describe genome-wide estimates of recombination rates in Cryptococcus deneoformans and compare recombination between progeny from α-α unisexual and a-α bisexual crosses. We find that offspring from bisexual crosses have modestly higher average rates of recombination than those derived from unisexual crosses. Recombination hot and cold spots across the C. deneoformans genome are also identified and are associated with increased GC content. Finally, we observed regions genome-wide with allele frequencies deviating from the expected parental ratio. These findings and observations advance our quantitative understanding of the genetic events that occur during sexual reproduction in C. deneoformans, and the impact that different forms of sexual reproduction are likely to have on genetic diversity in this important fungal pathogen.

RNAi is a critical determinant of centromere evolution in closely related fungi.

The centromere DNA locus on a eukaryotic chromosome facilitates faithful chromosome segregation. Despite performing such a conserved function, centromere DNA sequence as well as the organization of sequence elements is rapidly evolving in all forms of eukaryotes. The driving force that facilitates centromere evolution remains an enigma. Here, we studied the evolution of centromeres in closely related species in the fungal phylum of Basidiomycota. Using ChIP-seq analysis of conserved inner kinetochore proteins, we identified centromeres in three closely related Cryptococcus species: two of which are RNAi-proficient, while the other lost functional RNAi. We find that the centromeres in the RNAi-deficient species are significantly shorter than those of the two RNAi-proficient species. While centromeres are LTR retrotransposon-rich in all cases, the RNAi-deficient species lost all full-length retroelements from its centromeres. In addition, centromeres in RNAi-proficient species are associated with a significantly higher level of cytosine DNA modifications compared with those of RNAi-deficient species. Furthermore, when an RNAi-proficient Cryptococcus species and its RNAi-deficient mutants were passaged under similar conditions, the centromere length was found to be occasionally shortened in RNAi mutants. In silico analysis of predicted centromeres in a group of closely related Ustilago species, also belonging to the Basidiomycota, were found to have undergone a similar transition in the centromere length in an RNAi-dependent fashion. Based on the correlation found in two independent basidiomycetous species complexes, we present evidence suggesting that the loss of RNAi and cytosine DNA methylation triggered transposon attrition, which resulted in shortening of centromere length during evolution.

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

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

Natural mismatch repair mutations mediate phenotypic diversity and drug resistance in Cryptococcus deuterogattii.

Pathogenic microbes confront an evolutionary conflict between the pressure to maintain genome stability and the need to adapt to mounting external stresses. Bacteria often respond with elevated mutation rates, but little evidence exists of stable eukaryotic hypermutators in nature. Whole genome resequencing of the human fungal pathogen Cryptococcus deuterogattii identified an outbreak lineage characterized by a nonsense mutation in the mismatch repair component MSH2. This defect results in a moderate mutation rate increase in typical genes, and a larger increase in genes containing homopolymer runs. This allows facile inactivation of genes with coding homopolymer runs including FRR1, which encodes the target of the immunosuppresive antifungal drugs FK506 and rapamycin. Our study identifies a eukaryotic hypermutator lineage spread over two continents and suggests that pathogenic eukaryotic microbes may experience similar selection pressures on mutation rate as bacterial pathogens, particularly during long periods of clonal growth or while expanding into new environments.

Importance of Resolving Fungal Nomenclature: the Case of Multiple Pathogenic Species in the Cryptococcus Genus.

Cryptococcosis is a major fungal disease caused by members of the Cryptococcus gattii and Cryptococcus neoformans species complexes. After more than 15 years of molecular genetic and phenotypic studies and much debate, a proposal for a taxonomic revision was made. The two varieties within C. neoformans were raised to species level, and the same was done for five genotypes within C. gattii. In a recent perspective (K. J. Kwon-Chung et al., mSphere 2:e00357-16, 2017, https://doi.org/10.1128/mSphere.00357-16), it was argued that this taxonomic proposal was premature and without consensus in the community. Although the authors of the perspective recognized the existence of genetic diversity, they preferred the use of the informal nomenclature "C. neoformans species complex" and "C. gattii species complex." Here we highlight the advantage of recognizing these seven species, as ignoring these species will impede deciphering further biologically and clinically relevant differences between them, which may in turn delay future clinical advances.

Fungal genome and mating system transitions facilitated by chromosomal translocations involving intercentromeric recombination.

Species within the human pathogenic Cryptococcus species complex are major threats to public health, causing approximately 1 million annual infections globally. Cryptococcus amylolentus is the most closely known related species of the pathogenic Cryptococcus species complex, and it is non-pathogenic. Additionally, while pathogenic Cryptococcus species have bipolar mating systems with a single large mating type (MAT) locus that represents a derived state in Basidiomycetes, C. amylolentus has a tetrapolar mating system with 2 MAT loci (P/R and HD) located on different chromosomes. Thus, studying C. amylolentus will shed light on the transition from tetrapolar to bipolar mating systems in the pathogenic Cryptococcus species, as well as its possible link with the origin and evolution of pathogenesis. In this study, we sequenced, assembled, and annotated the genomes of 2 C. amylolentus isolates, CBS6039 and CBS6273, which are sexual and interfertile. Genome comparison between the 2 C. amylolentus isolates identified the boundaries and the complete gene contents of the P/R and HD MAT loci. Bioinformatic and chromatin immunoprecipitation sequencing (ChIP-seq) analyses revealed that, similar to those of the pathogenic Cryptococcus species, C. amylolentus has regional centromeres (CENs) that are enriched with species-specific transposable and repetitive DNA elements. Additionally, we found that while neither the P/R nor the HD locus is physically closely linked to its centromere in C. amylolentus, and the regions between the MAT loci and their respective centromeres show overall synteny between the 2 genomes, both MAT loci exhibit genetic linkage to their respective centromere during meiosis, suggesting the presence of recombinational suppressors and/or epistatic gene interactions in the MAT-CEN intervening regions. Furthermore, genomic comparisons between C. amylolentus and related pathogenic Cryptococcus species provide evidence that multiple chromosomal rearrangements mediated by intercentromeric recombination have occurred during descent of the 2 lineages from their common ancestor. Taken together, our findings support a model in which the evolution of the bipolar mating system was initiated by an ectopic recombination event mediated by similar repetitive centromeric DNA elements shared between chromosomes. This translocation brought the P/R and HD loci onto the same chromosome, and further chromosomal rearrangements then resulted in the 2 MAT loci becoming physically linked and eventually fusing to form the single contiguous MAT locus that is now extant in the pathogenic Cryptococcus species.

Elucidation of the calcineurin-Crz1 stress response transcriptional network in the human fungal pathogen Cryptococcus neoformans.

Calcineurin is a highly conserved Ca2+/calmodulin-dependent serine/threonine-specific protein phosphatase that orchestrates cellular Ca2+ signaling responses. In Cryptococcus neoformans, calcineurin is activated by multiple stresses including high temperature, and is essential for stress adaptation and virulence. The transcription factor Crz1 is a major calcineurin effector in Saccharomyces cerevisiae and other fungi. Calcineurin dephosphorylates Crz1, thereby enabling Crz1 nuclear translocation and transcription of target genes. Here we show that loss of Crz1 confers phenotypes intermediate between wild-type and calcineurin mutants, and demonstrate that deletion of the calcineurin docking domain results in the inability of Crz1 to translocate into the nucleus under thermal stress. RNA-sequencing revealed 102 genes that are regulated in a calcineurin-Crz1-dependent manner at 37°C. The majority of genes were down-regulated in cna1Δ and crz1Δ mutants, indicating these genes are normally activated by the calcineurin-Crz1 pathway at high temperature. About 58% of calcineurin-Crz1 target genes have unknown functions, while genes with known or predicted functions are involved in cell wall remodeling, calcium transport, and pheromone production. We identified three calcineurin-dependent response element motifs within the promoter regions of calcineurin-Crz1 target genes, and show that Crz1 binding to target gene promoters is increased upon thermal stress in a calcineurin-dependent fashion. Additionally, we found a large set of genes independently regulated by calcineurin, and Crz1 regulates 59 genes independently of calcineurin. Given the intermediate crz1Δ mutant phenotype, and our recent evidence for a calcineurin regulatory network impacting mRNA in P-bodies and stress granules independently of Crz1, calcineurin likely acts on factors beyond Crz1 that govern mRNA expression/stability to operate a branched transcriptional/post-transcriptional stress response network necessary for fungal virulence. Taken together, our findings reveal the core calcineurin-Crz1 stress response cascade is maintained from ascomycetes to a pathogenic basidiomycete fungus, but its output in C. neoformans appears to be adapted to promote fungal virulence.

Identification of small RNAs in extracellular vesicles from the commensal yeast Malassezia sympodialis.

Malassezia is the dominant fungus in the human skin mycobiome and is associated with common skin disorders including atopic eczema (AE)/dermatitis. Recently, it was found that Malassezia sympodialis secretes nanosized exosome-like vesicles, designated MalaEx, that carry allergens and can induce inflammatory cytokine responses. Extracellular vesicles from different cell-types including fungi have been found to deliver functional RNAs to recipient cells. In this study we assessed the presence of small RNAs in MalaEx and addressed if the levels of these RNAs differ when M. sympodialis is cultured at normal human skin pH versus the elevated pH present on the skin of patients with AE. The total number and the protein concentration of the released MalaEx harvested after 48 h culture did not differ significantly between the two pH conditions nor did the size of the vesicles. From small RNA sequence data, we identified a set of reads with well-defined start and stop positions, in a length range of 16 to 22 nucleotides consistently present in the MalaEx. The levels of small RNAs were not significantly differentially expressed between the two different pH conditions indicating that they are not influenced by the elevated pH level observed on the AE skin.