Characterizing the role of RNA silencing components in Cryptococcus neoformans

Distinct evolutionary trajectories following loss of RNA interference in Cryptococcus neoformans

While increased mutation rates typically have negative consequences in multicellular organisms, hypermutation can be advantageous for microbes adapting to the environment. Previously, we identified two hypermutator Cryptococcus neoformans clinical isolates that rapidly develop drug resistance due to transposition of a retrotransposon, Cnl1. Cnl1-mediated hypermutation is caused by a nonsense mutation in the gene encoding a novel RNAi component, Znf3, combined with a tremendous transposon burden. To elucidate adaptative mechanisms following RNAi loss, two bioinformatic pipelines were developed to identify RNAi loss-of-function mutations in a collection of 387 sequenced C. neoformans isolates. Remarkably, several RNAi-loss isolates were identified that are not hypermutators and have not accumulated transposons. To test if these RNAi loss-of-function mutations can cause hypermutation, the mutations were introduced into a non-hypermutator strain with a high transposon burden, which resulted in a hypermutator phenotype. To further investigate if RNAi-loss isolates can become hypermutators, in vitro passaging was performed. Although no hypermutators were found in two C. neoformans RNAi-loss strains after short-term passage, hypermutation was observed in a passaged C. deneoformans strain with increased transposon burden. Consistent with a two-step evolution, when an RNAi-loss isolate was crossed with an isolate containing a high Cnl1 burden, F1 hypermutator progeny inheriting a high transposon burden were identified. In addition to Cnl1 transpositions, insertions of a novel gigantic DNA transposon KDZ1 (∼11 kb), contributed to hypermutation in the progeny. Our results suggest that RNAi loss is relatively common (7/387, ∼1.8%) and enables distinct evolutionary trajectories: hypermutation following transposon accumulation or survival without hypermutation.

Significance Statement

There is a dearth of antifungal drugs available to treat Cryptococcus neoformans , a human fungal pathogen of global impact. We previously identified natural hypermutators with a loss-of-function mutation in the RNAi machinery and transposon expansion. Here, we identified several novel natural isolates with RNAi defects, none of which are hypermutators or have undergone transposon expansion. Furthermore, we demonstrate that these isolates can lie on a pathway to hypermutation following introduction of a transposon burden. In addition, a novel DNA transposon class was discovered that contributes to antifungal drug resistance. These findings highlight the importance of transposons in driving rapid adaptation in the absence of RNAi and reveal distinct evolutionary trajectories following RNAi loss, a relatively common event in C. neoformans .

VmRDR2 of Valsa mali mediates the generation of VmR2-siR1 that suppresses apple resistance by RNA interference

Cross-kingdom RNA interference (RNAi) is a crucial mechanism in host-pathogen interactions, with RNA-dependent RNA polymerase (RdRP) playing a vital role in signal amplification during RNAi. However, the role of pathogenic fungal RdRP in siRNAs generation and the regulation of plant-pathogen interactions remains elusive. Using deep sequencing, molecular, genetic, and biochemical approaches, this study revealed that VmRDR2 of Valsa mali regulates VmR2-siR1 to suppress the disease resistance-related gene MdLRP14 in apple. Both VmRDR1 and VmRDR2 are essential for the pathogenicity of V. mali in apple, with VmRDR2 mediating the generation of endogenous siRNAs, including an infection-related siRNA, VmR2-siR1. This siRNA specifically degrades the apple intracellular LRR-RI protein gene MdLRP14 in a sequence-specific manner, and overexpression of MdLRP14 enhances apple resistance against V. mali, which can be suppressed by VmR2-siR1. Conversely, MdLRP14 knockdown reduces resistance. In summary, this study demonstrates that VmRDR2 contributes to the generation of VmR2-siR1, which silences the host's intracellular LRR protein gene, thereby inhibiting host resistance. These findings offer novel insights into the fungi-mediated pathogenicity mechanism through RNAi.

Illuminating the Cryptococcus neoformans species complex: unveiling intracellular structures with fluorescent-protein-based markers

Cryptococcus neoformans is a fungal pathogen of the top critical priority recognized by the World Health Organization. This clinically important fungus also serves as a eukaryotic model organism. A variety of resources have been generated to facilitate investigation of the C. neoformans species complex, including congenic pairs, well-annotated genomes, genetic editing tools, and gene deletion sets. Here, we generated a set of strains with all major organelles fluorescently marked. We tested short organelle-specific targeting sequences and successfully labeled the following organelles by fusing the targeting sequences with a fluorescence protein: the plasma membrane, the nucleus, the peroxisome, and the mitochondrion. We used native cryptococcal Golgi and late endosomal proteins fused with a fluorescent protein to label these two organelles. These fluorescence markers were verified via colocalization using organelle-specific dyes. All the constructs for the fluorescent protein tags were integrated in an intergenic safe haven region. These organelle-marked strains were examined for growth and various phenotypes. We demonstrated that these tagged strains could be employed to track cryptococcal interaction with the host in phagocytosis assays. These strains also allowed us to discover remarkable differences in the dynamics of proteins targeted to different organelles during sexual reproduction. Additionally, we revealed that "dormant" spores transcribed and synthesized their own proteins and trafficked the proteins to the appropriate subcellular compartments, demonstrating that spores are metabolically active. We anticipate that these newly generated fluorescent markers will greatly facilitate further investigation of cryptococcal biology and pathogenesis.

Comparative genomics of the closely related fungal genera Cryptococcus and Kwoniella reveals karyotype dynamics and suggests evolutionary mechanisms of pathogenesis

In exploring the evolutionary trajectories of both pathogenesis and karyotype dynamics in fungi, we conducted a large-scale comparative genomic analysis spanning the Cryptococcus genus, encompassing both global human fungal pathogens and nonpathogenic species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species, covering virtually all known diversity within these genera. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at preadaptive pathogenic potential, our analysis found evidence of gene gain (via horizontal gene transfer) and gene loss in pathogenic Cryptococcus species, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the 2 genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5, or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes showed reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Overall, our findings advance our understanding of genetic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Calcineurin contributes to RNAi-mediated transgene silencing and small interfering RNA production in the human fungal pathogen Cryptococcus neoformans

Adaptation to external environmental challenges at the cellular level requires rapid responses and involves relay of information to the nucleus to drive key gene expression changes through downstream transcription factors. Here, we describe an alternative route of adaptation through a direct role for cellular signaling components in governing gene expression via RNA interference-mediated small RNA production. Calcium-calcineurin signaling is a highly conserved signaling cascade that plays central roles in stress adaptation and virulence of eukaryotic pathogens, including the human fungal pathogen Cryptococcus neoformans. Upon activation in C. neoformans, calcineurin localizes to P-bodies, membraneless organelles that are also the site for RNA processing. Here, we studied the role of calcineurin and its substrates in RNAi-mediated transgene silencing. Our results reveal that calcineurin regulates both the onset and the reversion of transgene silencing. We found that some calcineurin substrates that localize to P-bodies also regulate transgene silencing but in opposing directions. Small RNA sequencing in mutants lacking calcineurin or its targets revealed a role for calcineurin in small RNA production. Interestingly, the impact of calcineurin and its substrates was found to be different in genome-wide analysis, suggesting that calcineurin may regulate small RNA production in C. neoformans through additional pathways. Overall, these findings define a mechanism by which signaling machinery induced by external stimuli can directly alter gene expression to accelerate adaptative responses and contribute to genome defense.

Calcineurin contributes to RNAi-mediated transgene silencing and small interfering RNA production in the human fungal pathogen Cryptococcus neoformans

Adaptation to external environmental challenges at the cellular level requires rapid responses and involves relay of information to the nucleus to drive key gene expression changes through downstream transcription factors. Here, we describe an alternative route of adaptation through a direct role for cellular signaling components in governing gene expression via RNA interference-mediated small RNA production. Calcium-calcineurin signaling is a highly conserved signaling cascade that plays central roles in stress adaptation and virulence of eukaryotic pathogens, including the human fungal pathogen Cryptococcus neoformans. Upon activation in C. neoformans, calcineurin localizes to P-bodies, membrane-less organelles that are also the site for RNA processing. Here, we studied the role of calcineurin and its substrates in RNAi-mediated transgene silencing. Our results reveal that calcineurin regulates both the onset and the reversion of transgene silencing. We found that some calcineurin substrates that localize to P-bodies also regulate transgene silencing but in opposing directions. Small RNA sequencing in mutants lacking calcineurin or its targets revealed a role for calcineurin in small RNA production. Interestingly, the impact of calcineurin and its substrates was found to be different in genome-wide analysis, suggesting that calcineurin may regulate small RNA production in C. neoformans through additional pathways. Overall, these findings define a mechanism by which signaling machinery induced by external stimuli can directly alter gene expression to accelerate adaptative responses and contribute to genome defense.

Comparative genomics of Cryptococcus and Kwoniella reveals pathogenesis evolution and contrasting karyotype dynamics via intercentromeric recombination or chromosome fusion

A large-scale comparative genomic analysis was conducted for the global human fungal pathogens within the Cryptococcus genus, compared to non-pathogenic Cryptococcus species, and related species from the sister genus Kwoniella. Chromosome-level genome assemblies were generated for multiple species of both genera, resulting in a dataset encompassing virtually all of their known diversity. Although Cryptococcus and Kwoniella have comparable genome sizes (about 19.2 and 22.9 Mb) and similar gene content, hinting at pre-adaptive pathogenic potential, our analysis found evidence in pathogenic Cryptococcus species of specific examples of gene gain (via horizontal gene transfer) and gene loss, which might represent evolutionary signatures of pathogenic development. Genome analysis also revealed a significant variation in chromosome number and structure between the two genera. By combining synteny analysis and experimental centromere validation, we found that most Cryptococcus species have 14 chromosomes, whereas most Kwoniella species have fewer (11, 8, 5 or even as few as 3). Reduced chromosome number in Kwoniella is associated with formation of giant chromosomes (up to 18 Mb) through repeated chromosome fusion events, each marked by a pericentric inversion and centromere loss. While similar chromosome inversion-fusion patterns were observed in all Kwoniella species with fewer than 14 chromosomes, no such pattern was detected in Cryptococcus. Instead, Cryptococcus species with less than 14 chromosomes, underwent chromosome reductions primarily through rearrangements associated with the loss of repeat-rich centromeres. Additionally, Cryptococcus genomes exhibited frequent interchromosomal translocations, including intercentromeric recombination facilitated by transposons shared between centromeres. Taken together, our findings advance our understanding of genomic changes possibly associated with pathogenicity in Cryptococcus and provide a foundation to elucidate mechanisms of centromere loss and chromosome fusion driving distinct karyotypes in closely related fungal species, including prominent global human pathogens.

Gene Silencing via RNA Interference in Cryptococcus

RNA interference (RNAi) is a molecular biology technique for silencing specific eukaryotic genes without altering the DNA sequence in the genome. The silencing effect occurs because of decreased levels of mRNA that then result in decreased protein levels for the gene. The specificity of the silencing is dependent upon the presence of sequence-specific double-stranded RNA (dsRNA) that activates the cellular RNAi machinery. This chapter describes the process of silencing a specific target gene in Cryptococcus using a dual promoter vector. The plasmid, pIBB103, was designed with two convergent GAL7 promoters flanking a ura5 fragment that acts as a reporter for efficient RNAi. The target gene fragment is inserted between the promoters to be transcribed from both directions leading to the production of dsRNA in cells that activate the RNAi pathway.

Simplified and effective RNA interference and CRISPR-Cas9 systems for Cryptococcus neoformans

The 3,4-dihydroxyphenylalanine (DOPA) melanin is one of the important virulence factors for Cryptococcus neoformans, which may trigger immune responses in the host. While the production of DOPA melanin is catalyzed by laccase that is predominantly encoded by LAC1 gene. Therefore, regulating the genetic expression of C. neoformans is conducive to exploring the impact of interested molecules on the host. In this work, we established two systems that were constructed quickly and easily for the knock-down/knock-out of LAC1 gene: RNA interference (RNAi) and clustered regularly interspaced short palindromic repeats CRISPR-Cas9. The RNAi system was constructed by pSilencer 4.1-CMV neo plasmid and short hairpin RNA to achieve effective transcriptional suppression. The CRISPR-Cas9 system was used the PNK003 vectors to obtain a stable albino mutant strain. The results of phenotype, quantitative real-time polymerase chain reaction, transmission electron microscope, and spectrophotometry were used to assess the ability of melanin production. As a result, the RNAi system displayed attenuation of transcriptional suppression when the transformants continuously passed on new plates. However, the transcriptional suppression of long loop in short hairpin RNA was more powerful and lasted longer. An albino strain produced by CRISPR-Cas9 was completely unable to synthesize melanin. In conclusion, strains with different capacities of melanin production were obtained by RNAi and CRISPR-Cas9 systems, which might be useful for exploring the linear relation between melanin and immunoreaction of the host. In addition, the two systems in this article might be convenient to quickly screen the possible trait-regulating genes of other serotypes of C. neoformans.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining (MAT) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By generating additional chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere-flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight toward the identification of a sexual cycle in Malassezia, with possible implications for pathogenicity.

Disruption of the Aspergillus fumigatus RNA interference machinery alters the conidial transcriptome

The RNA interference (RNAi) pathway has evolved numerous functionalities in eukaryotes, with many on display in Kingdom Fungi. RNAi can regulate gene expression, facilitate drug resistance, or even be altogether lost to improve growth potential in some fungal pathogens. In the WHO fungal priority pathogen, Aspergillus fumigatus, the RNAi system is known to be intact and functional. To extend our limited understanding of A. fumigatus RNAi, we first investigated the genetic variation in RNAi-associated genes in a collection of 217 environmental and 83 clinical genomes, where we found that RNAi components are conserved even in clinical strains. Using endogenously expressed inverted-repeat transgenes complementary to a conditionally essential gene (pabA) or a nonessential gene (pksP), we determined that a subset of the RNAi componentry is active in inverted-repeat transgene silencing in conidia and mycelium. Analysis of mRNA-seq data from RNAi double-knockout strains linked the A. fumigatus dicer-like enzymes (DclA/B) and RNA-dependent RNA polymerases (RrpA/B) to regulation of conidial ribosome biogenesis genes; however, surprisingly few endogenous small RNAs were identified in conidia that could explain this broad change. Although RNAi was not clearly linked to growth or stress response defects in the RNAi knockouts, serial passaging of RNAi knockout strains for six generations resulted in lineages with diminished spore production over time, indicating that loss of RNAi can exert a fitness cost on the fungus. Cumulatively, A. fumigatus RNAi appears to play an active role in defense against double-stranded RNA species alongside a previously unappreciated housekeeping function in regulation of conidial ribosomal biogenesis genes.

Frequent transitions in mating-type locus chromosomal organization in Malassezia and early steps in sexual reproduction

Fungi in the basidiomycete genus Malassezia are the most prevalent eukaryotic microbes resident on the skin of human and other warm-blooded animals and have been implicated in skin diseases and systemic disorders. Analysis of Malassezia genomes revealed that key adaptations to the skin microenvironment have a direct genomic basis, and the identification of mating/meiotic genes suggests a capacity to reproduce sexually, even though no sexual cycle has yet been observed. In contrast to other bipolar or tetrapolar basidiomycetes that have either two linked mating-type-determining ( MAT ) loci or two MAT loci on separate chromosomes, in Malassezia species studied thus far the two MAT loci are arranged in a pseudobipolar configuration (linked on the same chromosome but capable of recombining). By incorporating newly generated chromosome-level genome assemblies, and an improved Malassezia phylogeny, we infer that the pseudobipolar arrangement was the ancestral state of this group and revealed six independent transitions to tetrapolarity, seemingly driven by centromere fission or translocations in centromere- flanking regions. Additionally, in an approach to uncover a sexual cycle, Malassezia furfur strains were engineered to express different MAT alleles in the same cell. The resulting strains produce hyphae reminiscent of early steps in sexual development and display upregulation of genes associated with sexual development as well as others encoding lipases and a protease potentially relevant for pathogenesis of the fungus. Our study reveals a previously unseen genomic relocation of mating-type loci in fungi and provides insight towards the discovery of a sexual cycle in Malassezia , with possible implications for pathogenicity.

Significance Statement

Malassezia , the dominant fungal group of the mammalian skin microbiome, is associated with numerous skin disorders. Sexual development and yeast-to-hyphae transitions, governed by genes at two mating-type ( MAT ) loci, are thought to be important for fungal pathogenicity. However, Malassezia sexual reproduction has never been observed. Here, we used chromosome-level assemblies and comparative genomics to uncover unforeseen transitions in MAT loci organization within Malassezia , possibly related with fragility of centromeric-associated regions. Additionally, by expressing different MAT alleles in the same cell, we show that Malassezia can undergo hyphal development and this phenotype is associated with increased expression of key mating genes along with other genes known to be virulence factors, providing a possible connection between hyphal development, sexual reproduction, and pathogenicity.

Importance of Clinical Isolates in Cryptococcus neoformans Research

The human pathogenic fungus Cryptococcus neoformans is a global health concern. Previous research in the field has focused on studies using reference strains to identify virulence factors, generate mutant libraries, define genomic structures, and perform functional studies. In this review, we discuss the benefits and drawbacks of using reference strains to study C. neoformans, describe how the study of clinical isolates has expanded our understanding of pathogenesis, and highlight how studies using clinical isolates can further develop our understanding of the host-pathogen interaction during C. neoformans infection.

Exogenous double-stranded RNA inhibits the infection physiology of rust fungi to reduce symptoms in planta

Rust fungi (Pucciniales) are a diverse group of plant pathogens in natural and agricultural systems. They pose ongoing threats to the diversity of native flora and cause annual crop yield losses. Agricultural rusts are predominantly managed with fungicides and breeding for resistance, but new control strategies are needed on non-agricultural plants and in fragile ecosystems. RNA interference (RNAi) induced by exogenous double-stranded RNA (dsRNA) has promise as a sustainable approach for managing plant-pathogenic fungi, including rust fungi. We investigated the mechanisms and impact of exogenous dsRNA on rust fungi through in vitro and whole-plant assays using two species as models, Austropuccinia psidii (the cause of myrtle rust) and Coleosporium plumeriae (the cause of frangipani rust). In vitro, dsRNA either associates externally or is internalized by urediniospores during the early stages of germination. The impact of dsRNA on rust infection architecture was examined on artificial leaf surfaces. dsRNA targeting predicted essential genes significantly reduced germination and inhibited development of infection structures, namely appressoria and penetration pegs. Exogenous dsRNA sprayed onto 1-year-old trees significantly reduced myrtle rust symptoms. Furthermore, we used comparative genomics to assess the wide-scale amenability of dsRNA to control rust fungi. We sequenced genomes of six species of rust fungi, including three new families (Araucariomyceaceae, Phragmidiaceae, and Skierkaceae) and identified key genes of the RNAi pathway across 15 species in eight families of Pucciniales. Together, these findings indicate that dsRNA targeting essential genes has potential for broad-use management of rust fungi across natural and agricultural systems.

Genome-wide analysis of heat stress-stimulated transposon mobility in the human fungal pathogen Cryptococcus deneoformans

We recently reported transposon mutagenesis as a significant driver of spontaneous mutations in the human fungal pathogen Cryptococcus deneoformans during murine infection. Mutations caused by transposable element (TE) insertion into reporter genes were dramatically elevated at high temperatures (37° vs. 30°) in vitro, suggesting that heat stress stimulates TE mobility in the Cryptococcus genome. To explore the genome-wide impact of TE mobilization, we generated transposon accumulation lines by in vitro passage of C. deneoformans strain XL280α for multiple generations at both 30° and at the host-relevant temperature of 37°. Utilizing whole-genome sequencing, we identified native TE copies and mapped multiple de novo TE insertions in these lines. Movements of the T1 DNA transposon occurred at both temperatures with a strong bias for insertion between gene-coding regions. By contrast, the Tcn12 retrotransposon integrated primarily within genes and movement occurred exclusively at 37°. In addition, we observed a dramatic amplification in copy number of the Cnl1 (Cryptococcus neoformans LINE-1) retrotransposon in subtelomeric regions under heat-stress conditions. Comparing TE mutations to other sequence variations detected in passaged lines, the increase in genomic changes at elevated temperatures was primarily due to mobilization of the retroelements Tcn12 and Cnl1. Finally, we found multiple TE movements (T1, Tcn12, and Cnl1) in the genomes of single C. deneoformans isolates recovered from infected mice, providing evidence that mobile elements are likely to facilitate microevolution and rapid adaptation during infection.

The Cryptococcus gattii species complex: Unique pathogenic yeasts with understudied virulence mechanisms

Members of Cryptococcus gattii/neoformans species complex are the etiological agents of the potentially fatal human fungal infection cryptococcosis. C. gattii and its sister species cause disease in both immunocompetent and immunocompromised hosts, while the closely related species C. neoformans and C. deneoformans predominantly infect immunocompromised hosts. To date, most studies have focused on similarities in pathogenesis between these two groups, but over recent years, important differences have become apparent. In this review paper, we highlight some of the major phenotypic differences between the C. gattii and neoformans species complexes and justify the need to study the virulence and pathogenicity of the C. gattii species complex as a distinct cryptococcal group.

Uncontrolled transposition following RNAi loss causes hypermutation and antifungal drug resistance in clinical isolates of Cryptococcus neoformans

Cryptococcus neoformans infections cause approximately 15% of AIDS-related deaths owing to a combination of limited antifungal therapies and drug resistance. A collection of clinical and environmental C. neoformans isolates were assayed for increased mutation rates via fluctuation analysis, and we identified two hypermutator C. neoformans clinical isolates with increased mutation rates when exposed to the combination of rapamycin and FK506. Sequencing of drug target genes found that Cnl1 transposon insertions conferred the majority of resistance to rapamycin and FK506 and could also independently cause resistance to 5-fluoroorotic acid and the clinically relevant antifungal 5-flucytosine. Whole-genome sequencing revealed both hypermutator genomes harbour a nonsense mutation in the RNA-interference component ZNF3 and hundreds of Cnl1 elements organized into massive subtelomeric arrays on each of the fourteen chromosomes. Quantitative trait locus mapping in 28 progeny derived from a cross between a hypermutator and wild-type identified a locus associated with hypermutation that included znf3. CRISPR editing of the znf3 nonsense mutation abolished hypermutation and restored small-interfering-RNA production. We conclude that hypermutation and drug resistance in these clinical isolates result from RNA-interference loss and accumulation of Cnl1 elements.

Dynamic genome plasticity during unisexual reproduction in the human fungal pathogen Cryptococcus deneoformans

Genome copy number variation occurs during each mitotic and meiotic cycle and it is crucial for organisms to maintain their natural ploidy. Defects in ploidy transitions can lead to chromosome instability, which is a hallmark of cancer. Ploidy in the haploid human fungal pathogen Cryptococcus neoformans is exquisitely orchestrated and ranges from haploid to polyploid during sexual development and under various environmental and host conditions. However, the mechanisms controlling these ploidy transitions are largely unknown. During C. deneoformans (formerly C. neoformans var. neoformans, serotype D) unisexual reproduction, ploidy increases prior to the onset of meiosis, can be independent from cell-cell fusion and nuclear fusion, and likely occurs through an endoreplication pathway. To elucidate the molecular mechanisms underlying this ploidy transition, we identified twenty cell cycle-regulating genes encoding cyclins, cyclin-dependent kinases (CDK), and CDK regulators. We characterized four cyclin genes and two CDK regulator genes that were differentially expressed during unisexual reproduction and contributed to diploidization. To detect ploidy transition events, we generated a ploidy reporter, called NURAT, which can detect copy number increases via double selection for nourseothricin-resistant, uracil-prototrophic cells. Utilizing this ploidy reporter, we showed that ploidy transition from haploid to diploid can be detected during the early phases of unisexual reproduction. Interestingly, selection for the NURAT reporter revealed several instances of segmental aneuploidy of multiple chromosomes, which conferred azole resistance in some isolates. These findings provide further evidence of ploidy plasticity in fungi with significant biological and public health implications.

Ploidy is an intrinsic fundamental feature of all eukaryotic organisms, and ploidy variation and maintenance are critical to the organism survival and evolution. Fungi exhibit exquisite plasticity in ploidy variation in adaptation to various environmental stresses. For example, the haploid opportunistic human fungal pathogen C. deneoformans can generate diploid blastospores during unisexual reproduction and also forms polyploid titan cells during host infection; however, the mechanisms underlying these ploidy transitions are largely unknown. In this study, we elucidated the genetic regulatory circuitry governing ploidy duplication during C. deneoformans unisexual reproduction through the identification and characterization of cell cycle regulators that are differentially expressed during unisexual reproduction. We showed that four cyclin and two cyclin-dependent kinase regulator genes function in concert to orchestrate ploidy transition during unisexual reproduction. To trace and track ploidy transition events, we also generated a ploidy reporter and revealed the formation of segmental aneuploidy in addition to diploidization, illustrating the diverse mechanisms of genome plasticity in C. deneoformans.

Role of Dicer-Dependent RNA Interference in Regulating Mycoparasitic Interactions

Dicer-like proteins (DCLs) play a vital role in RNA interference (RNAi), by cleaving RNA filament into small RNAs. Although DCL-mediated RNAi can regulate interspecific communication between pathogenic/mutualistic organisms and their hosts, its role in mycoparasitic interactions is yet to be investigated. In this study, we deleted dcl genes in the mycoparasitic fungus Clonostachys rosea and characterize the functions of DCL-dependent RNAi in mycoparasitism. Deletion of dcl2 resulted in a mutant with reduced secondary metabolite production, antagonism toward the plant-pathogenic fungus Botrytis cinerea, and reduced ability to control Fusarium foot rot disease on wheat, caused by Fusarium graminearum. Transcriptome sequencing of the in vitro interaction between the C. rosea Δdcl2 strain and B. cinerea or F. graminearum identified the downregulation of genes coding for transcription factors, membrane transporters, hydrolytic enzymes, and secondary metabolites biosynthesis enzymes putatively involved in antagonistic interactions, in comparison with the C. rosea wild-type interaction. A total of 61 putative novel microRNA-like RNAs (milRNAs) were identified in C. rosea, and 11 were downregulated in the Δdcl2 mutant. In addition to putative endogenous gene targets, these milRNAs were predicted to target B. cinerea and F. graminearum virulence factor genes, which showed an increased expression during interaction with the Δdcl2 mutant incapable of producing the targeting milRNAs. In summary, this study constitutes the first step in elucidating the role of RNAi in mycoparasitic interactions, with important implications for biological control of plant diseases, and poses the base for future studies focusing on the role of cross-species RNAi regulating mycoparasitic interactions. IMPORTANCE Small RNAs mediated RNA interference (RNAi) known to regulate several biological processes. Dicer-like endoribonucleases (DCLs) play a vital role in the RNAi pathway by generating sRNAs. In this study, we investigated a role of DCL-mediated RNAi in interference interactions between mycoparasitic fungus Clonostachys rosea and the two fungal pathogens Botrytis cinerea and Fusarium graminearum (here called mycohosts). We found that the dcl mutants were not able to produce 11 sRNAs predicted to finetune the regulatory network of genes known to be involved in production of hydrolytic enzymes, antifungal compounds, and membrane transporters needed for antagonistic action of C. rosea. We also found C. rosea sRNAs putatively targeting known virulence factors in the mycohosts, indicating RNAi-mediated cross-species communication. Our study expanded the understanding of underlying mechanisms of cross-species communication during interference interactions and poses a base for future works studying the role of DCL-based cross-species RNAi in fungal interactions.

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.

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.

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.

Dicer promotes Atg8 expression through RNAi independent mechanism in Cryptococcus neoformans

ATG8 is one of the critical genes that participate in several essential autophagic steps. The expression of ATG8 must be exquisitely regulated to avoid physiological disorder and even cell death. However, the mechanisms of regulating ATG8 expression remain to be fully uncovered. In this investigation, we found that Dicer homologs in Cryptococcus neoformans could activate the expression of ATG8 independent of RNAi. Deletion of two Dicer homologs (DCR1 and DCR2) from C. neoformans, especially DCR2, led to significantly reduced Atg8 protein level, but deletion of other RNAi components did not result in the same phenotype. The autophagic flux, the numbers of autophagic bodies and the tolerance to glucose starvation of dcr2∆ were also significantly reduced. Further investigation showed that Dcr2 activates the expression of ATG8 through the promoter region, not the Open Reading Frame or 3' Untranslated Region. We also found that a similar phenomenon exists in mammalian cells, as DCR1 instead of AGO2 knockdown also reduced the expression of LC3, indicating that this mechanism may be conservative in eukaryotic cells. Therefore, a novel transcription activation mechanism was revealed in this paper.

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 Evolutionary Significance of RNAi in the Fungal Kingdom

RNA interference (RNAi) was discovered at the end of last millennium, changing the way scientists understood regulation of gene expression. Within the following two decades, a variety of different RNAi mechanisms were found in eukaryotes, reflecting the evolutive diversity that RNAi entails. The essential silencing mechanism consists of an RNase III enzyme called Dicer that cleaves double-stranded RNA (dsRNA) generating small interfering RNAs (siRNAs), a hallmark of RNAi. These siRNAs are loaded into the RNA-induced silencing complex (RISC) triggering the cleavage of complementary messenger RNAs by the Argonaute protein, the main component of the complex. Consequently, the expression of target genes is silenced. This mechanism has been thoroughly studied in fungi due to their proximity to the animal phylum and the conservation of the RNAi mechanism from lower to higher eukaryotes. However, the role and even the presence of RNAi differ across the fungal kingdom, as it has evolved adapting to the particularities and needs of each species. Fungi have exploited RNAi to regulate a variety of cell activities as different as defense against exogenous and potentially harmful DNA, genome integrity, development, drug tolerance, or virulence. This pathway has offered versatility to fungi through evolution, favoring the enormous diversity this kingdom comprises.

Advances in Molecular Tools and In Vivo Models for the Study of Human Fungal Pathogenesis

Pathogenic fungi represent an increasing infectious disease threat to humans, especially with an increasing challenge of antifungal drug resistance. Over the decades, numerous tools have been developed to expedite the study of pathogenicity, initiation of disease, drug resistance and host-pathogen interactions. In this review, we highlight advances that have been made in the use of molecular tools using CRISPR technologies, RNA interference and transposon targeted mutagenesis. We also discuss the use of animal models in modelling disease of human fungal pathogens, focusing on zebrafish, the silkworm, Galleria mellonella and the murine model.

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.

Transposon mobilization in the human fungal pathogen Cryptococcus is mutagenic during infection and promotes drug resistance in vitro

When transitioning from the environment, pathogenic microorganisms must adapt rapidly to survive in hostile host conditions. This is especially true for environmental fungi that cause opportunistic infections in immunocompromised patients since these microbes are not well adapted human pathogens. Cryptococcus species are yeastlike fungi that cause lethal infections, especially in HIV-infected patients. Using Cryptococcus deneoformans in a murine model of infection, we examined contributors to drug resistance and demonstrated that transposon mutagenesis drives the development of 5-fluoroorotic acid (5FOA) resistance. Inactivation of target genes URA3 or URA5 primarily reflected the insertion of two transposable elements (TEs): the T1 DNA transposon and the TCN12 retrotransposon. Consistent with in vivo results, increased rates of mutagenesis and resistance to 5FOA and the antifungal drugs rapamycin/FK506 (rap/FK506) and 5-fluorocytosine (5FC) were found when Cryptococcus was incubated at 37° compared to 30° in vitro, a condition that mimics the temperature shift that occurs during the environment-to-host transition. Inactivation of the RNA interference (RNAi) pathway, which suppresses TE movement in many organisms, was not sufficient to elevate TE movement at 30° to the level observed at 37°. We propose that temperature-dependent TE mobilization in Cryptococcus is an important mechanism that enhances microbial adaptation and promotes pathogenesis and drug resistance in the human host.

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.

Chromatin-Mediated Regulation of Genome Plasticity in Human Fungal Pathogens

Human fungal pathogens, such as Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans, are a public health problem, causing millions of infections and killing almost half a million people annually. The ability of these pathogens to colonise almost every organ in the human body and cause life-threating infections relies on their capacity to adapt and thrive in diverse hostile host-niche environments. Stress-induced genome instability is a key adaptive strategy used by human fungal pathogens as it increases genetic diversity, thereby allowing selection of genotype(s) better adapted to a new environment. Heterochromatin represses gene expression and deleterious recombination and could play a key role in modulating genome stability in response to environmental changes. However, very little is known about heterochromatin structure and function in human fungal pathogens. In this review, I use our knowledge of heterochromatin structure and function in fungal model systems as a road map to review the role of heterochromatin in regulating genome plasticity in the most common human fungal pathogens: Candida albicans, Aspergillus fumigatus and Cryptococcus neoformans.

Life Cycle of Cryptococcus neoformans

Cryptococcus neoformans is a ubiquitous environmental fungus and an opportunistic pathogen that causes fatal cryptococcal meningitis. Advances in genomics, genetics, and cellular and molecular biology of C. neoformans have dramatically improved our understanding of this important pathogen, rendering it a model organism to study eukaryotic biology and microbial pathogenesis. In light of recent progress, we describe in this review the life cycle of C. neoformans with a special emphasis on the regulation of the yeast-to-hypha transition and different modes of sexual reproduction, in addition to the impacts of the life cycle on cryptococcal populations and pathogenesis.

A Non-Dicer RNase III and Four Other Novel Factors Required for RNAi-Mediated Transposon Suppression in the Human Pathogenic Yeast Cryptococcus neoformans

The human pathogenic yeast Cryptococcus neoformans silences transposable elements using endo-siRNAs and an Argonaute, Ago1. Endo-siRNAs production requires the RNA-dependent RNA polymerase, Rdp1, and two partially redundant Dicer enzymes, Dcr1 and Dcr2, but is independent of histone H3 lysine 9 methylation. We describe here an insertional mutagenesis screen for factors required to suppress the mobilization of the C. neoformans HARBINGER family DNA transposon HAR1. Validation experiments uncovered five novel genes (RDE1-5) required for HAR1 suppression and global production of suppressive endo-siRNAs. The RDE genes do not impact transcript levels, suggesting the endo-siRNAs do not act by impacting target transcript synthesis or turnover. RDE3 encodes a non-Dicer RNase III related to S. cerevisiaeRnt1, RDE4 encodes a predicted terminal nucleotidyltransferase, while RDE5 has no strongly predicted encoded domains. Affinity purification-mass spectrometry studies suggest that Rde3 and Rde5 are physically associated. RDE1 encodes a G-patch protein homologous to the S. cerevisiaeSqs1/Pfa1, a nucleolar protein that directly activates the essential helicase Prp43 during rRNA biogenesis. Rde1 copurifies Rde2, another novel protein obtained in the screen, as well as Ago1, a homolog of Prp43, and numerous predicted nucleolar proteins. We also describe the isolation of conditional alleles of PRP43, which are defective in RNAi. This work reveals unanticipated requirements for a non-Dicer RNase III and presumptive nucleolar factors for endo-siRNA biogenesis and transposon mobilization suppression in C. neoformans.

Studying fungal pathogens of humans and fungal infections: fungal diversity and diversity of approaches

Seminal work by Louis Pasteur revealed the contribution of fungi - yeasts and microsporidia to agroindustry and disease in animals, respectively. More than 150 years later, the impact of fungi on human health and beyond is an ever-increasing issue, although often underestimated. Recent studies estimate that fungal infections, especially those caused by Candida, Cryptococcus and Aspergillus species, kill more than one million people annually. Indeed, these neglected infections are in general very difficult to cure and the associated mortality remains very high even when antifungal treatments exist. The development of new antifungals and diagnostic tools that are both necessary to fight fungal diseases efficiently, requires greater insights in the biology of the fungal pathogens of humans in the context of the infection, on their epidemiology, and on their role in the human mycobiota. We also need a better understanding of the host immune responses to fungal pathogens as well as the genetic basis for the increased sensitivity of some individuals to fungal infections. Here, we highlight some recent progress made in these different areas of research, in particular based on work conducted in our own laboratories. These progresses should lay the ground for better management of fungal infections, as they provide opportunities for better diagnostic, vaccination, the development of classical antifungals but also strategies for targeting virulence factors or the host.

Comprehensive analysis of full genome sequence and Bd-milRNA/target mRNAs to discover the mechanism of hypovirulence in Botryosphaeria dothidea strains on pear infection with BdCV1 and BdPV1

Pear ring rot disease, mainly caused by Botryosphaeria dothidea, is widespread in most pear and apple-growing regions. Mycoviruses are used for biocontrol, especially in fruit tree disease. BdCV1 (Botryosphaeria dothidea chrysovirus 1) and BdPV1 (Botryosphaeria dothidea partitivirus 1) influence the biological characteristics of B. dothidea strains. BdCV1 is a potential candidate for the control of fungal disease. Therefore, it is vital to explore interactions between B. dothidea and mycovirus to clarify the pathogenic mechanisms of B. dothidea and hypovirulence of B. dothidea in pear. A high-quality full-length genome sequence of the B. dothidea LW-Hubei isolate was obtained using Single Molecule Real-Time sequencing. It has high repeat sequence with 9.3% and DNA methylation existence in the genome. The 46.34 Mb genomes contained 14,091 predicted genes, which of 13,135 were annotated. B. dothidea was predicted to express 3833 secreted proteins. In bioinformatics analysis, 351 CAZy members, 552 transporters, 128 kinases, and 1096 proteins associated with plant-host interaction (PHI) were identified. RNA-silencing components including two endoribonuclease Dicer, four argonaute (Ago) and three RNA-dependent RNA polymerase (RdRp) molecules were identified and expressed in response to mycovirus infection. Horizontal transfer of the LW-C and LW-P strains indicated that BdCV1 induced host gene silencing in LW-C to suppress BdPV1 transmission. To investigate the role of RNA-silencing in B. dothidea defense, we constructed four small RNA libraries and sequenced B. dothidea micro-like RNAs (Bd-milRNAs) produced in response to BdCV1 and BdPV1 infection. Among these, 167 conserved and 68 candidate novel Bd-milRNAs were identified, of which 161 conserved and 20 novel Bd-milRNA were differentially expressed. WEGO analysis revealed involvement of the differentially expressed Bd-milRNA-targeted genes in metabolic process, catalytic activity, cell process and response to stress or stimulus. BdCV1 had a greater effect on the phenotype, virulence, conidiomata, vertical and horizontal transmission ability, and mycelia cellular structure biological characteristics of B. dothidea strains than BdPV1 and virus-free strains. The results obtained in this study indicate that mycovirus regulates biological processes in B. dothidea through the combined interaction of antiviral defense mediated by RNA-silencing and milRNA-mediated regulation of target gene mRNA expression.

Studying fungal pathogens of humans and fungal infections: fungal diversity and diversity of approaches

Seminal work by Louis Pasteur revealed the contribution of fungi-yeasts and microsporidia to agroindustry and disease in animals, respectively. More than 150 years later, the impact of fungi on human health and beyond is an ever-increasing issue, although often underestimated. Recent studies estimate that fungal infections, especially those caused by Candida, Cryptococcus and Aspergillus species, kill more than one million people annually. Indeed, these neglected infections are in general very difficult to cure and the associated mortality remains very high even when antifungal treatments exist. The development of new antifungals and diagnostic tools that are both necessary to fight fungal diseases efficiently, requires greater insights in the biology of the fungal pathogens of humans in the context of the infection, on their epidemiology, and on their role in the human mycobiota. We also need a better understanding of the host immune responses to fungal pathogens as well as the genetic basis for the increased sensitivity of some individuals to fungal infections. Here, we highlight some recent progress made in these different areas of research, in particular based on work conducted in our own laboratories. These progress should lay the ground for better management of fungal infections, as they provide opportunities for better diagnostic, vaccination, the development of classical antifungals but also strategies for targeting virulence factors or the host.

In silico characterization of microRNAs-like sequences in the genome of Paracoccidioides brasiliensis

Eukaryotic cells have different mechanisms of post-transcriptional regulation. Among these mechanisms, microRNAs promote regulation of targets by cleavage or degradation of the mRNA. Fungi of the Paracoccidioides complex are the etiological agents of the main systemic mycosis of Latin America. These fungi present a plasticity to adapt and survive in different conditions, and the presence of microRNAs-like molecules could be part of the mechanisms that provide such plasticity. MicroRNAs produced by the host influence the progression of this mycosis in the lungs besides regulating targets involved in apoptosis in macrophage, activation of T and B cells and the production of cytokines. Therefore, this work analyzed the presence of regions in the genome of this fungus with a potential to encode microRNAs-like molecules. Here we show by analysis of sequence similarity the presence of 18 regions, putatively coding for microRNAs-like molecules in the Paracoccidioides brasiliensis genome. We also described the conservation of dicer and argonaut proteins and the cognate transcripts induced in the yeast parasitic phase. This work represents a starting point for the analysis of the presence of those molecules in the morphological stages of the fungus and their role in fungal development.

Dicer-Like Proteins Regulate Sexual Development via the Biogenesis of Perithecium-Specific MicroRNAs in a Plant Pathogenic Fungus Fusarium graminearum

Ascospores act as the primary inoculum of Fusarium graminearum, which causes the destructive disease Fusarium head blight (FHB), or scab. MicroRNAs (miRNAs) have been reported in the F. graminearum vegetative stage, and Fgdcl2 is involved in microRNA-like RNA (milRNA) biogenesis but has no major impact on vegetative growth, abiotic stress or pathogenesis. In the present study, we found that ascospore discharge was decreased in the Fgdcl1 deletion mutant, and completely blocked in the double-deletion mutant of Fgdcl1 and Fgdcl2. Besides, more immature asci were observed in the double-deletion mutant. Interestingly, the up-regulated differentially expressed genes (DEGs) common to ΔFgdcl1 and ΔFgdcl1/2 were related to ion transmembrane transporter and membrane components. The combination of small RNA and transcriptome sequencing with bioinformatics analysis predicted 143 novel milRNAs in wild-type perithecia, and 138 of these milRNAs partly or absolutely depended on Fgdcl1, while only 5 novel milRNAs were still obtained in the Fgdcl1 and Fgdcl2 double-deletion mutant. Furthermore, 117 potential target genes were predicted. Overall, Fgdcl1 and Fgdcl2 genes were partly functionally redundant in ascospore discharge and perithecium-specific milRNA generation in F. graminearum, and these perithecium-specific milRNAs play potential roles in sexual development.

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

The “centromere paradox” refers to rapidly evolving and highly diverse centromere DNA sequences even in closely related eukaryotes. However, factors contributing to this rapid divergence are largely unknown. Here, we identified large regional, LTR retrotransposon-rich centromeres in a group of human fungal pathogens belonging to the Cryptococcus species complex. We provide evidence that loss-of-functional RNAi machinery and possibly cytosine DNA methylation trigger instability of the genome by activation of centromeric retrotransposons presumably suppressed by RNAi. We propose that RNAi, together with cytosine DNA methylation, serves as a critical determinant that maintains repetitive transposon-rich centromere structures. This study explores the direct link between RNAi and centromere structure evolution.

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.

The RNAi Universe in Fungi: A Varied Landscape of Small RNAs and Biological Functions

RNA interference (RNAi) is a conserved eukaryotic mechanism that uses small RNA molecules to suppress gene expression through sequence-specific messenger RNA degradation, translational repression, or transcriptional inhibition. In filamentous fungi, the protective function of RNAi in the maintenance of genome integrity is well known. However, knowledge of the regulatory role of RNAi in fungi has had to wait until the recent identification of different endogenous small RNA classes, which are generated by distinct RNAi pathways. In addition, RNAi research on new fungal models has uncovered the role of small RNAs and RNAi pathways in the regulation of diverse biological functions. In this review, we give an up-to-date overview of the different classes of small RNAs and RNAi pathways in fungi and their roles in the defense of genome integrity and regulation of fungal physiology and development, as well as in the interaction of fungi with biotic and abiotic environments.

Genomic CG dinucleotide deficiencies associated with transposable element hypermutation in Basidiomycetes, some lower fungi, a moss and a clubmoss

Many Basidiomycete genomes include substantial fractions that are deficient in CG dinucleotides, in extreme cases amounting to 70% of the genome. CG deficiency is variable and correlates with genome size and, more closely, with transposable element (TE) content. Many species have limited CG deficiency; it is therefore likely that there are other mechanisms that can control TE proliferation. Examination of TEs confirms that C-to-T transition mutations in CG dinucleotides may comprise a conspicuous proportion of differences between paired elements, however transition/transversion ratios are never as high as those due to RIP in some Ascomycetes, suggesting that repeat-associated CG mutation is not totally pervasive. This has allowed gene family expansion in Basidiomycetes, although CG transition differences are often prominent in paired gene family members, and are evidently responsible for destruction of some copies. A few lower fungal genomes exhibit similar evidence of repeat-associated CG mutation, as do the genomes of the two lower plants Physcomitrella patens and Selaginella moellendorffii, in both of which mutation parallels published methylation of CHG as well as CG nucleotides. In Basidiomycete DNA methylation has been reported to be largely confined to CG dinucleotides in repetitive DNA, but while methylation and mutation are evidently associated, it is not clear which is cause and which effect.

sRNAs as possible regulators of retrotransposon activity in Cryptococcus gattii VGII

Background

The absence of Argonaute genes in the fungal pathogen Cryptococcus gattii R265 and other VGII strains indicates that yeasts of this genotype cannot have a functional RNAi pathway, an evolutionarily conserved gene silencing mechanism performed by small RNAs. The success of the R265 strain as a pathogen that caused the Pacific Northwest and Vancouver Island outbreaks may imply that RNAi machinery loss could be beneficial under certain circumstances during evolution. As a result, a hypermutant phenotype would be created with high rates of genome retrotransposition, for instance. This study therefore aimed to evaluate in silicio the effect of retrotransposons and their control mechanisms by small RNAs on genomic stability and synteny loss of C. gattii R265 through retrotransposons sequence comparison and orthology analysis with other 16 C. gattii genomic sequences available.

Results

Retrotransposon mining identified a higher sequence count to VGI genotype compared to VGII, VGIII, and VGIV. However, despite the lower retrotransposon number, VGII exhibited increased synteny loss and genome rearrangement events. RNA-Seq analysis indicated highly expressed retrotransposons as well as sRNA production.

Conclusions

Genome rearrangement and synteny loss may suggest a greater retrotransposon mobilization caused by RNAi pathway absence, but the effective presence of sRNAs that matches retrotransposon sequences means that an alternative retrotransposon silencing mechanism could be active in genomic integrity maintenance of C. gattii VGII strains.

Electronic supplementary material

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

Small RNA-Based Antiviral Defense in the Phytopathogenic Fungus Colletotrichum higginsianum

Even though the fungal kingdom contains more than 3 million species, little is known about the biological roles of RNA silencing in fungi. The Colletotrichum genus comprises fungal species that are pathogenic for a wide range of crop species worldwide. To investigate the role of RNA silencing in the ascomycete fungus Colletotrichum higginsianum, knock-out mutants affecting genes for three RNA-dependent RNA polymerase (RDR), two Dicer-like (DCL), and two Argonaute (AGO) proteins were generated by targeted gene replacement. No effects were observed on vegetative growth for any mutant strain when grown on complex or minimal media. However, Δdcl1, Δdcl1Δdcl2 double mutant, and Δago1 strains showed severe defects in conidiation and conidia morphology. Total RNA transcripts and small RNA populations were analyzed in parental and mutant strains. The greatest effects on both RNA populations was observed in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains, in which a previously uncharacterized dsRNA mycovirus [termed Colletotrichum higginsianum non-segmented dsRNA virus 1 (ChNRV1)] was derepressed. Phylogenetic analyses clearly showed a close relationship between ChNRV1 and members of the segmented Partitiviridae family, despite the non-segmented nature of the genome. Immunoprecipitation of small RNAs associated with AGO1 showed abundant loading of 5’U-containing viral siRNA. C. higginsianum parental and Δdcl1 mutant strains cured of ChNRV1 revealed that the conidiation and spore morphology defects were primarily caused by ChNRV1. Based on these results, RNA silencing involving ChDCL1 and ChAGO1 in C. higginsianum is proposed to function as an antiviral mechanism.

Colletotrichum sp. comprises a diverse group of fungal pathogens that attack over 3000 plant species worldwide. Understanding the underlying mechanisms that govern fungal development and pathogenicity may enable more effective and sustainable approaches to crop disease management and control. In most organisms, RNA silencing is an important mechanism to control endogenous and exogenous RNA. RNA silencing utilizes small regulatory molecules (small RNAs) produced by proteins called Dicer (DCL), and exercise their function though effector proteins named Argonaute (AGO). Here, we investigated the role of RNA silencing machinery in the fungus Colletotrichum higginsianum, by generating deletions in genes encoding RNA silencing components. Severe defects were observed in both conidiation and conidia morphology in the Δdcl1, Δdcl1Δdcl2, and Δago1 strains. Analysis of transcripts and small RNAs revealed an uncharacterized dsRNA virus persistently infecting C. higginsianum. The virus was shown (1) to be de-repressed in the Δdcl1, Δdcl1Δdcl2 and Δago1 strains, and (2) to cause the conidiation and spore mutant phenotypes. Our results indicate that C. higginsianum employs RNA silencing as an antiviral mechanism to suppress viruses and their debilitating effects.

Transcriptional control of sexual development in Cryptococcus neoformans

Developmental processes are essential for the normal life cycles of many pathogenic fungi, and they can facilitate survival in challenging environments, including the human host. Sexual development of the human fungal pathogen Cryptococcus neoformans not only produces infectious particles (spores) but has also enabled the evolution of new disease-related traits such as drug resistance. Transcription factor networks are essential to the development and pathogenesis of C. neoformans, and a variety of sequence-specific DNA-binding proteins control both key developmental transitions and virulence by regulating the expression of their target genes. In this review we discuss the roles of known transcription factors that harbor important connections to both development and virulence. Recent studies of these transcription factors have identified a common theme in which metabolic, stress, and other responses that are required for sexual development appear to have been co-opted for survival in the human host, thus facilitating pathogenesis. Future work elucidating the connection between development and pathogenesis will provide vital insights into the evolution of complex traits in eukaryotes as well as mechanisms that may be used to combat fungal pathogens.

Gene Network Polymorphism Illuminates Loss and Retention of Novel RNAi Silencing Components in the Cryptococcus Pathogenic Species Complex

RNAi is a ubiquitous pathway that serves central functions throughout eukaryotes, including maintenance of genome stability and repression of transposon expression and movement. However, a number of organisms have lost their RNAi pathways, including the model yeast Saccharomyces cerevisiae, the maize pathogen Ustilago maydis, the human pathogen Cryptococcus deuterogattii, and some human parasite pathogens, suggesting there may be adaptive benefits associated with both retention and loss of RNAi. By comparing the RNAi-deficient genome of the Pacific Northwest Outbreak C. deuterogattii strain R265 with the RNAi-proficient genomes of the Cryptococcus pathogenic species complex, we identified a set of conserved genes that were lost in R265 and all other C. deuterogattii isolates examined. Genetic and molecular analyses reveal several of these lost genes play roles in RNAi pathways. Four novel components were examined further. Znf3 (a zinc finger protein) and Qip1 (a homolog of N. crassa Qip) were found to be essential for RNAi, while Cpr2 (a constitutive pheromone receptor) and Fzc28 (a transcription factor) are involved in sex-induced but not mitosis-induced silencing. Our results demonstrate that the mitotic and sex-induced RNAi pathways rely on the same core components, but sex-induced silencing may be a more specific, highly induced variant that involves additional specialized or regulatory components. Our studies further illustrate how gene network polymorphisms involving known components of key cellular pathways can inform identification of novel elements and suggest that RNAi loss may have been a core event in the speciation of C. deuterogattii and possibly contributed to its pathogenic trajectory.

Genome instability and mutations provoked by transposon movement are counteracted by novel defense mechanisms in organisms as diverse as fungi, plants, and mammals. In the human fungal pathogen Cryptococcus neoformans, an RNAi silencing pathway operates to defend the genome against mobile elements and transgene repeats. RNAi silencing pathways are conserved in the Cryptococcus pathogenic species complex and are mediated by canonical RNAi components. Surprisingly, several of these components are missing from all analyzed C. deuterogattii VGII strains, the molecular type responsible for the North American Pacific Northwest outbreak. To identify novel components of the RNAi pathways, we surveyed the reference genomes of C. deuterogattii, C. gattii, C. neoformans, and C. deneoformans. We identified 14 otherwise conserved genes missing in R265, including the RDP1, AGO1, and DCR1 canonical RNAi components, and focused on four potentially novel RNAi components: ZNF3, QIP1, CPR2, and FZC28. We found that Znf3 and Qip1 are both required for mitotic- and sex-induced silencing, while Cpr2 and Fzc28 contribute to sex-induced but not mitosis-induced silencing. Our studies reveal elements of RNAi pathways that operate to defend the genome during sexual development and vegetative growth and illustrate the power of network polymorphisms to illuminate novel components of biological pathways.

What makes Cryptococcus gattii a pathogen?

Cryptococcosis is an invasive fungal infection of humans and other animals, typically caused by the species Cryptococcus neoformans in patients with impaired immunity. However, there is growing recognition of the importance of the related species C. gattii in causing infections in apparently immunocompetent individuals. In particular, an ongoing outbreak of cryptococcal disease in the Pacific Northwest region, which started in 1999, has driven an intense research effort into this previously neglected pathogen. Here, we discuss some of the recent discoveries in this organism from the Pacific Northwest region and highlight areas for future investigation.

Genome Evolution and Innovation across the Four Major Lineages of Cryptococcus gattii

Cryptococcus gattii is a fungal pathogen of humans, causing pulmonary infections in otherwise healthy hosts. To characterize genomic variation among the four major lineages of C. gattii (VGI, -II, -III, and -IV), we generated, annotated, and compared 16 de novo genome assemblies, including the first for the rarely isolated lineages VGIII and VGIV. By identifying syntenic regions across assemblies, we found 15 structural rearrangements, which were almost exclusive to the VGI-III-IV lineages. Using synteny to inform orthology prediction, we identified a core set of 87% of C. gattii genes present as single copies in all four lineages. Remarkably, 737 genes are variably inherited across lineages and are overrepresented for response to oxidative stress, mitochondrial import, and metal binding and transport. Specifically, VGI has an expanded set of iron-binding genes thought to be important to the virulence of Cryptococcus, while VGII has expansions in the stress-related heat shock proteins relative to the other lineages. We also characterized genes uniquely absent in each lineage, including a copper transporter absent from VGIV, which influences Cryptococcus survival during pulmonary infection and the onset of meningoencephalitis. Through inclusion of population-level data for an additional 37 isolates, we identified a new transcontinental clonal group that we name VGIIx, mitochondrial recombination between VGII and VGIII, and positive selection of multidrug transporters and the iron-sulfur protein aconitase along multiple branches of the phylogenetic tree. Our results suggest that gene expansion or contraction and positive selection have introduced substantial variation with links to mechanisms of pathogenicity across this species complex.

The genetic differences between phenotypically different pathogens provide clues to the underlying mechanisms of those traits and can lead to new drug targets and improved treatments for those diseases. In this paper, we compare 16 genomes belonging to four highly differentiated lineages of Cryptococcus gattii, which cause pulmonary infections in otherwise healthy humans and other animals. Half of these lineages have not had their genomes previously assembled and annotated. We identified 15 ancestral rearrangements in the genome and over 700 genes that are unique to one or more lineages, many of which are associated with virulence. In addition, we found evidence for recent transcontinental spread, mitochondrial genetic exchange, and positive selection in multidrug transporters. Our results suggest that gene expansion/contraction and positive selection are diversifying the mechanisms of pathogenicity across this species complex.

Characterization of RNA silencing components in the plant pathogenic fungus Fusarium graminearum

The RNA interference (RNAi) plays a critical role in gene regulation in a variety of eukaryotic organisms. However, the role of RNAi remains largely unclear in plant pathogenic fungi. In this study, we explored the roles of core components of the RNAi pathway in Fusarium graminearum, the major causal agent of wheat head blight. Our results demonstrated that the hairpin RNA (hpRNA) can efficiently silence the expression level of target gene, and the argonaute protein FgAgo1 and dicer protein FgDicer2 are important in this silencing process. RNAi machinery was not involved in growth, abiotic stress and pathogenesis in F. graminearum under tested conditions. We firstly applied high-throughput sequencing technology to elucidate small RNA (17-40 nucleotides) (sRNA) transcriptome in F. graminearum, and found that a total of forty-nine micro-like-RNA (milRNA) candidates were identified in the wild-type and ∆FgDICER2, and twenty-four of them were FgDicer2-dependent. Fg-milRNA-4 negatively regulated expression of its target gene. Taken together, our results indicated that the hpRNA-induced gene silencing was a valuable genetic tool for exploring gene function in F. graminearum. FgAgo1 and FgDicer2 proteins played a critical role in the hpRNA mediated gene silencing process. In addition, FgDicer2 was involved in sRNA transcription and milRNA generation in this fungus.

Epigenetics as an emerging tool for improvement of fungal strains used in biotechnology

Filamentous fungi are today a major source of industrial biotechnology for the production of primary and secondary metabolites, as well as enzymes and recombinant proteins. All of them have undergone extensive improvement strain programs, initially by classical mutagenesis and later on by genetic manipulation. Thereby, strategies to overcome rate-limiting or yield-reducing reactions included manipulating the expression of individual genes, their regulatory genes, and also their function. Yet, research of the last decade clearly showed that cells can also undergo heritable changes in gene expression that do not involve changes in the underlying DNA sequences (=epigenetics). This involves three levels of regulation: (i) DNA methylation, (ii) chromatin remodeling by histone modification, and (iii) RNA interference. The demonstration of the occurrence of these processes in fungal model organisms such as Aspergillus nidulans and Neurospora crassa has stimulated its recent investigation as a tool for strain improvement in industrially used fungi. This review describes the progress that has thereby been obtained.

New technology and resources for cryptococcal research

Rapid advances in molecular biology and genome sequencing have enabled the generation of new technology and resources for cryptococcal research. RNAi-mediated specific gene knock down has become routine and more efficient by utilizing modified shRNA plasmids and convergent promoter RNAi constructs. This system was recently applied in a high-throughput screen to identify genes involved in host-pathogen interactions. Gene deletion efficiencies have also been improved by increasing rates of homologous recombination through a number of approaches, including a combination of double-joint PCR with split-marker transformation, the use of dominant selectable markers and the introduction of Cre-Loxp systems into Cryptococcus. Moreover, visualization of cryptococcal proteins has become more facile using fusions with codon-optimized fluorescent tags, such as green or red fluorescent proteins or, mCherry. Using recent genome-wide analytical tools, new transcriptional factors and regulatory proteins have been identified in novel virulence-related signaling pathways by employing microarray analysis, RNA-sequencing and proteomic analysis.