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

Spray-induced gene silencing to control plant pathogenic fungi: A step-by-step guide

RNA interference (RNAi)-based control technologies are gaining popularity as potential alternatives to synthetic fungicides in the ongoing effort to manage plant pathogenic fungi. Among these methods, spray-induced gene silencing (SIGS) emerges as particularly promising due to its convenience and feasibility for development. This approach is a new technology for plant disease management, in which double-stranded RNAs (dsRNAs) targeting essential or virulence genes are applied to plants or plant products and subsequently absorbed by plant pathogens, triggering a gene silencing effect and the inhibition of the infection process. Spray-induced gene silencing has demonstrated efficacy in laboratory settings against various fungal pathogens. However, as research progressed from the laboratory to the greenhouse and field environments, novel challenges arose, such as ensuring the stability of dsRNAs and their effective delivery to fungal targets. Here, we provide a practical guide to SIGS for the control of plant pathogenic fungi. This guide outlines the essential steps and considerations needed for designing and assessing dsRNA molecules. It also addresses key challenges inherent to SIGS, including delivery and stability of dsRNA molecules, and how nanoencapsulation of dsRNAs can aid in overcoming these obstacles. Additionally, the guide underscores existing knowledge gaps that warrant further research and aims to provide assistance to researchers, especially those new to the field, encouraging the advancement of SIGS for the control of a broad range of fungal pathogens.

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 a gene encoding an RNA interference (RNAi) component, ZNF3, combined with a tremendous transposon burden. To elucidate adaptive mechanisms following RNAi loss, two bioinformatic pipelines were developed to identify RNAi loss-of-function (LOF) 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 whether these RNAi LOF mutations can cause hypermutation, the mutations were introduced into a nonhypermutator strain with a high transposon burden, which resulted in a hypermutator phenotype. To further investigate whether 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 Cryptococcus deneoformans strain with an 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 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.

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.

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.

Identification of an active RNAi pathway in Candida albicans

RNA interference (RNAi) is a fundamental regulatory pathway with a wide range of functions, including regulation of gene expression and maintenance of genome stability. Although RNAi is widespread in the fungal kingdom, well-known species, such as the model yeast Saccharomyces cerevisiae, have lost the RNAi pathway. Until now evidence has been lacking for a fully functional RNAi pathway in Candida albicans, a human fungal pathogen considered critically important by the World Health Organization. Here, we demonstrated that the widely used C. albicans reference strain (SC5314) contains an inactivating missense mutation in the gene encoding for the central RNAi component Argonaute. In contrast, most other C. albicans isolates contain a canonical Argonaute protein predicted to be functional and RNAi-active. Indeed, using high-throughput small and long RNA sequencing combined with seamless CRISPR/Cas9-based gene editing, we demonstrate that an active C. albicans RNAi machinery represses expression of subtelomeric gene families. Thus, an intact and functional RNAi pathway exists in C. albicans, highlighting the importance of using multiple reference strains when studying this dangerous pathogen.

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.

tRNA-derived fragments as new players in regulatory processes in yeast

For a very long time, RNA molecules were treated as transistory molecules, by which the genetic information flows from DNA to proteins; the model proposed in the 1960s accepted that proteins are both the products and the regulators of gene expression. Since then, thousands of reports proved that RNAs should be thought about as the factors that do control gene expression. The pervasive transcription has been reported in many eukaryotic organisms, illustrating a highly interwoven transcriptome organization that includes hundreds of previously unknown noncoding RNAs. The key roles of noncoding RNAs (microRNAs and small interfering RNAs) in gene expression regulation are no longer surprising, as are new classes of noncoding RNAs constantly being discovered. Transfer RNAs (tRNAs) are the second most abundant type of RNAs in the cell. Advances in high-throughput sequencing technologies exposed the existence of functional, regulatory tRNA-derived RNA fragments (tRFs), generated from precursor and mature tRNAs. These tRF molecules have been found to play central roles during stress and different pathological conditions. Herein, we present the critical assessment of the discoveries made in the field of tRNA-derived fragments in the past 15 years in various pathogenic and nonpathogenic yeast species.

ZNF3 regulates proliferation, migration and invasion through MMP1 and TWIST in colorectal cancer

Colorectal cancer (CRC) is a malignant tumor with a high incidence and mortality worldwide. Currently, the underlying molecular mechanisms of CRC are still unclear. Zinc finger protein 3 (ZNF3) is a zinc-finger transcription factor that has been reported as a candidate for breast cancer prognosis, suggesting its involvement in the regulation of tumorigenesis. However, the association between ZNF3 and CRC remains unknown. To investigate the role of ZNF3 in CRC, we first analyze the correlation between ZNF3 expression and CRC, and the results demonstrate that ZNF3 is highly expressed in CRC tissue and cells, which is associated with the age of CRC patients. In vitro studies show that ZNF3 overexpression promotes CRC cell migration. Compared to control cells, knockdown of ZNF3 markedly suppresses CRC cell proliferation, migration and invasion and promotes G0/G1 phase cell cycle arrest. The expressions of the EMT-related markers TWIST and MMP1 are significantly decreased when ZNF3 is silenced. Additionally, overexpression of MMP1 and TWIST exacerbates CRC cell proliferation, accelerates the S phase cell cycle in ZNF3-knockdown SW480 cells, and increases cell migration and invasion through Transwell chambers. These data suggest that ZNF3 is involved in cellular proliferation, migration and invasion by regulating MMP1 and TWIST in CRC cells.

Similar Characteristics of siRNAs of Plant Viruses Which Replicate in Plant and Fungal Hosts

Simple Summary

RNA silencing in fungi was shown to confer antiviral defense against plant viruses. In this study, using high-throughput sequencing and bioinformatic analyses, we showed that small interfering RNAs (siRNAs) of cucumber mosaic virus and tobacco mosaic virus (TMV) which replicated in phytopathogenic fungi Rhizoctonia solani and Fusarium graminearum had similarities with viral siRNAs produced in plant hosts in regard to the size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini. Additionally, our results also determined that both F. graminearum DCL1 and DCL2 were involved in the production of TMV siRNAs. Thus, the fungal RNA silencing machineries have adaptive capabilities to recognize and process the genome of invading plant viruses.

Abstract

RNA silencing is a host innate antiviral mechanism which acts via the synthesis of viral-derived small interfering RNAs (vsiRNAs). We have previously reported the infection of phytopathogenic fungi by plant viruses such as cucumber mosaic virus (CMV) and tobacco mosaic virus (TMV). Furthermore, fungal RNA silencing was shown to suppress plant virus accumulation, but the characteristics of plant vsiRNAs associated with the antiviral response in this nonconventional host remain unknown. Using high-throughput sequencing, we characterized vsiRNA profiles in two plant RNA virus–fungal host pathosystems: CMV infection in phytopathogenic fungus Rhizoctonia solani and TMV infection in phytopathogenic fungus Fusarium graminearum. The relative abundances of CMV and TMV siRNAs in the respective fungal hosts were much lower than those in the respective experimental plant hosts, Nicotiana benthamiana and Nicotiana tabacum. However, CMV and TMV siRNAs in fungi had similar characteristics to those in plants, particularly in their size distributions, proportion of plus and minus senses, and nucleotide preference for the 5′ termini of vsiRNAs. The abundance of TMV siRNAs largely decreased in F. graminearum mutants with a deletion in either dicer-like 1 (dcl1) or dcl2 genes which encode key proteins for the production of siRNAs and antiviral responses. However, deletion of both dcl1 and dcl2 restored TMV siRNA accumulation in F. graminearum, indicating the production of dcl-independent siRNAs with no antiviral function in the absence of the dcl1 and dcl2 genes. Our results suggest that fungal RNA silencing recognizes and processes the invading plant RNA virus genome in a similar way as in plants.

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.

RNA-based therapeutics to treat human fungal infections

In recent decades, RNA-based therapeutics have transitioned from a near impossibility to a compelling treatment alternative for genetic disorders and infectious diseases. The mRNA vaccines against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are truly groundbreaking, and new adaptations are already being proposed to fight other microbes. Unfortunately, the potential of RNA-based therapeutics to treat human fungal infections has remained mostly absent from the conversation, despite the fact that invasive fungal infections kill as many per year as tuberculosis and even more than malaria. Here, we argue that RNA-based therapeutics should be investigated for the treatment of human fungal infections and discuss several major roadblocks and potential circumventions that may allow for the realization of RNA-based therapies against human fungal pathogens.

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.

Cryptococcus neoformans: Sex, morphogenesis, and virulence

Cryptococcus neoformans is a dimorphic fungus that causes lethal meningoencephalitis mainly in immunocompromised individuals. Different morphotypes enable this environmental fungus and opportunistic pathogen to adapt to different natural niches and exhibit different levels of pathogenicity in various hosts. It is well-recognized that C. neoformans undergoes bisexual or unisexual reproduction in vitro to generate genotypic, morphotypic, and phenotypic diversity, which augments its ability for adaptation. However, if and how sexual reproduction and the meiotic machinery exert any direct impact on the infection process is unclear. This review summarizes recent discoveries on the regulation of cryptococcal life cycle and morphogenesis, and how they impact cryptococcal pathogenicity. The potential role of the meiotic machinery on ploidy regulation during cryptococcal infection is also discussed. This review aims to stimulate further investigation on links between fungal morphogenesis, sexual reproduction, and virulence.

Analysis of tRNA-derived RNA fragments (tRFs) in Cryptococcus spp.: RNAi-independent generation and possible compensatory effects in a RNAi-deficient genotype

Small RNAs (sRNAs) are key factors in the regulation of gene expression. Recently, a new class of regulatory sRNAs derived from tRNAs was described, the tRNA-derived RNA fragments (tRFs). Such RNAs range in length from 14 to 30 nucleotides and are produced from both mature and primary tRNA transcripts, with very specific cleavage sites along the tRNA sequence. Although several mechanisms have been proposed for how tRFs mediate regulation of gene expression, the exact mechanism of tRF biogenesis and its dependency upon the RNAi pathway remain unclear. Cryptococcus gattii and Cryptococcus neoformans are basidiomycetous yeasts and important human pathogens. While C. neoformans is RNAi proficient, C. gattii VGII has lost essential RNAi genes. Here, we sought to identify the tRF production profile in C. gattii VGII and C. neoformans in order to assess the RNAi-dependency of tRF production in these fungal species. We developed a RNA-sequencing-based tRF prediction workflow designed to improve the currently available prediction tools. Using this methodology, we were able to identify tRFs in both organisms. Despite the loss of the RNAi pathway, C. gattii VGII displayed a number of identified tRFs that did not differ significantly from those observed in C. neoformans. The analysis of predicted tRF targets revealed that a higher number of targets was found for C. gattii VGII tRFs compared to C. neoformans tRFs. These results support the idea that tRFs are at least partially independent of the canonical RNAi machinery, raising questions about possible compensatory roles of alternative regulatory RNAs in the absence of a functional RNAi pathway.

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.

Role of Non-coding RNAs in Fungal Pathogenesis and Antifungal Drug Responses

Purpose of Review

Non-coding RNAs (ncRNAs), including regulatory small RNAs (sRNAs) and long non-coding RNAs (lncRNAs), constitute a significant part of eukaryotic genomes; however, their roles in fungi are just starting to emerge. ncRNAs have been shown to regulate gene expression in response to varying environmental conditions (like stress) and response to chemicals, including antifungal drugs. In this review, I highlighted recent studies focusing on the functional roles of ncRNAs in pathogenic fungi.

Recent Findings

Emerging evidence suggests sRNAs (small RNAs) and lncRNAs (long non-coding RNAs) play an important role in fungal pathogenesis and antifungal drug response. Their roles include posttranscriptional gene silencing, histone modification, and chromatin remodeling. Fungal pathogens utilize RNA interference (RNAi) mechanisms to regulate pathogenesis-related genes and can also transfer sRNAs inside the host to suppress host immunity genes to increase virulence. Hosts can also transfer sRNAs to induce RNAi in fungal pathogens to reduce virulence. Additionally, sRNAs and lncRNAs also regulate gene expression in response to antifungal drugs increasing resistance (and possibly tolerance) to drugs.

Summary

Herein, I discuss what is known about ncRNAs in fungal pathogenesis and antifungal drug responses. Advancements in genomic technologies will help identify the ncRNA repertoire in fungal pathogens, and functional studies will elucidate their mechanisms. This will advance our understanding of host-fungal interactions and potentially help develop better treatment strategies.

Zinc Finger Proteins in the Human Fungal Pathogen Cryptococcus neoformans

Zinc is one of the essential trace elements in eukaryotes and it is a critical structural component of a large number of proteins. Zinc finger proteins (ZNFs) are zinc-finger domain-containing proteins stabilized by bound zinc ions and they form the most abundant proteins, serving extraordinarily diverse biological functions. In recent years, many ZNFs have been identified and characterized in the human fungal pathogen Cryptococcus neoformans, a fungal pathogen causing fatal meningitis mainly in immunocompromised individuals. It has been shown that ZNFs play important roles in the morphological development, differentiation, and virulence of C. neoformans. In this review, we, first, briefly introduce the ZNFs and their classification. Then, we explain the identification and classification of the ZNFs in C. neoformans. Next, we focus on the biological role of the ZNFs functionally characterized so far in the sexual reproduction, virulence factor production, ion homeostasis, pathogenesis, and stress resistance in C. neoformans. We also discuss the perspectives on future function studies of ZNFs in C. neoformans.

RNA interference in the oleaginous yeast Rhodosporidium toruloides

The red yeast Rhodosporidium toruloides is an excellent microbial host for production of carotenoids, neutral lipids and valuable enzymes. In recent years, genetic tools for gene expression and gene disruption have been developed for this red yeast. However, methods remain limited in terms of fine-tuning gene expression. In this study, we first demonstrated successful implementation of RNA interference (RNAi) in R. toruloides NP11, which was applied to down-regulate the expression of autophagy related gene 8 (ATG8), and fatty acid synthase genes (FAS1 and FAS2), respectively. Compared with the control strain, RNAi-engineered strains showed a silencing efficiency ranging from 11% to 92%. The RNAi approach described here ensures selective inhibition of the target gene expression, and should expand our capacity in the genetic manipulation of R. toruloides for both fundamental research and advanced cell factory development.

The Evolution of Sexual Reproduction and the Mating-Type Locus: Links to Pathogenesis of Cryptococcus Human Pathogenic Fungi

Cryptococcus species utilize a variety of sexual reproduction mechanisms, which generate genetic diversity, purge deleterious mutations, and contribute to their ability to occupy myriad environmental niches and exhibit a range of pathogenic potential. The bisexual and unisexual cycles of pathogenic Cryptococcus species are stimulated by properties associated with their environmental niches and proceed through well-characterized signaling pathways and corresponding morphological changes. Genes governing mating are encoded by the mating-type (MAT) loci and influence pathogenesis, population dynamics, and lineage divergence in Cryptococcus. MAT has undergone significant evolutionary changes within the Cryptococcus genus, including transition from the ancestral tetrapolar state in nonpathogenic species to a bipolar mating system in pathogenic species, as well as several internal reconfigurations. Owing to the variety of established sexual reproduction mechanisms and the robust characterization of the evolution of mating and MAT in this genus, Cryptococcus species provide key insights into the evolution of sexual reproduction.

Transcription factor Znf2 coordinates with the chromatin remodeling SWI/SNF complex to regulate cryptococcal cellular differentiation

Cellular differentiation is instructed by developmental regulators in coordination with chromatin remodeling complexes. Much information about their coordination comes from studies in the model ascomycetous yeasts. It is not clear, however, what kind of information that can be extrapolated to species of other phyla in Kingdom Fungi. In the basidiomycete Cryptococcus neoformans, the transcription factor Znf2 controls yeast-to-hypha differentiation. Through a forward genetic screen, we identified the basidiomycete-specific factor Brf1. We discovered Brf1 works together with Snf5 in the SWI/SNF chromatin remodeling complex in concert with existent Znf2 to execute cellular differentiation. We demonstrated that SWI/SNF assists Znf2 in opening the promoter regions of hyphal specific genes, including the ZNF2 gene itself. This complex also supports Znf2 to fully associate with its target regions. Importantly, our findings revealed key differences in composition and biological function of the SWI/SNF complex in the two major phyla of Kingdom Fungi.

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.

E Pluribus Unum: The Fungal Kingdom as a Rosetta Stone for Biology and Medicine

THE Genetics Society of America's (GSA's) Edward Novitski Prize recognizes a single experimental accomplishment or a body of work in which an exceptional level of creativity, and intellectual ingenuity, has been used to design and execute scientific experiments to solve a difficult problem in genetics. The 2019 recipient is Joseph Heitman, who is recognized for his work on fungal pathogens of humans and for ingenious experiments using yeast to identify the molecular targets of widely used immunosuppressive drugs. The latter work, part of Heitman's postdoctoral research, proved to be a seminal contribution to the discovery of the conserved Target of Rapamycin (TOR) pathway. In his own research group, a recurring theme has been the linking of fundamental insights in fungal biology to medically important problems. His studies have included defining fungal mating-type loci, including their evolution and links to virulence, and illustrating convergent transitions from outcrossing to inbreeding in fungal pathogens of plants and animals. He has led efforts to establish new genetic and genomic methods for studying pathogenesis in Cryptococcus species. Heitman's group also discovered unisexual reproduction, a novel mode of fungal reproduction with implications for pathogen evolution and the origins of sexual reproduction.

Genetic and genomic evolution of sexual reproduction: echoes from LECA to the fungal kingdom

Sexual reproduction is vastly diverse and yet highly conserved across the eukaryotic domain. This ubiquity suggests that the last eukaryotic common ancestor (LECA) was sexual. It is hypothesized that several critical processes in sexual reproduction, including cell fusion and meiosis, were acquired during the evolution from the first eukaryotic common ancestor (FECA) to the sexual LECA. However, it is challenging to delineate the exact origin and evolution of sexual reproduction given that both FECA and LECA are extinct. Studies of diverse eukaryotes have helped to shed light on this sexual evolutionary trajectory, revealing that a primordial sexual ploidy cycle likely involved endoreplication followed by concerted chromosome loss and that cell-cell fusion, meiosis, and sex determination later arose to shape modern sexual reproduction. Despite the general conservation of sexual reproduction processes throughout eukaryotes, modern sexual cycles are immensely diverse and complex. This diversity and complexity has become readily apparent in the fungal kingdom with the recent rapid expansion of whole-genome sequencing. This abundance of data, the variety of genetic tools available to manipulate and characterize fungi, and the thorough characterization of many fungal sexual cycles make the fungal kingdom an excellent forum, in which to study the conservation and diversification of sexual reproduction.

Epigenetic mechanisms of drug resistance in fungi

The emergence of drug-resistant fungi poses a continuously increasing threat to human health. Despite advances in preventive care and diagnostics, resistant fungi continue to cause significant mortality, especially in immunocompromised patients. Therapeutic resources are further limited by current usage of only four major classes of antifungal drugs. Resistance against these drugs has already been observed in pathogenic fungi requiring the development of much needed newer antifungal drugs. Epigenetic changes such as DNA or chromatin modifications alter gene expression levels in response to certain stimuli, including interaction with the host in the case of fungal pathogens. These changes can confer resistance to drugs by altering the expression of target genes or genes encoding drug efflux pumps. Multiple pathogens share many of these epigenetic pathways; thus, targeting epigenetic pathways might also identify drug target candidates for the development of broad-spectrum antifungal drugs. In this review, we discuss the importance of epigenetic pathways in mediating drug resistance in fungi as well as in the development of anti-fungal drugs.

Estimating the Intra-taxa Diversity, Population Genetic Structure, and Evolutionary Pathways of Cryptococcus neoformans and Cryptococcus gattii

Members of the Cryptococcus complex, includes Cryptococcus neoformans (most common fungal infection of the brain) and Cryptococcus gattii (high-impact emerging pathogen worldwide). Currently, the fungal multilocus sequence typing database (Fungal MLST Database) constitutes a valuable data repository of the genes used for molecular typing of these pathogens. We analyzed the data available in the Fungal MLST Database for seven housekeeping genes, with the aim to evaluate its contribution in the description of intra-taxa diversity, population genetic structure, and evolutionary patterns. Although the Fungal MLST Database has a greater number of reports for C. neoformans (n = 487) than for C. gattii (n = 344), similar results were obtained for both species in terms of allelic diversity. Phylogenetic reconstructions revealed grouping by molecular type in both species and allowed us to propose differences in evolutionary patterns (gradualism in the case of C. neoformans and punctuated evolution in the case of C. gattii). In addition, C. neoformans showed a population genetic structure consisting of 37 clonal complexes (CCs; CC1 being predominant), high crosslinking [without sequence type (ST) grouping by molecular type], marked divergence events in phylogenetic analysis, and few introgression events (mainly between VNI and VNIV). By contrast, C. gattii showed 50 CCs (with greater homogeneity in ST number by CC) and clustering by molecular type with marked crosslinking events in phylogenetic networks being less evident. Understanding relationships at the molecular level for species of the Cryptococcus complex, based on the sequences of the housekeeping genes, provides information for describing the evolutionary history of these emerging pathogens.

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.

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.

As humans, we often think of genetic mutations as being bad. Over the past several decades we have seen health warnings issued on a variety of environmental exposures, from cigarettes to tanning beds, and with good reason because they cause mutations. For multicellular organisms like humans, these mutations are strongly associated with cancer. But in bacteria, this is not true. In fact, the rate at which mutations occur sometimes increases to help bacteria cope with stressful environments.

Unlike bacteria, humans are eukaryotes – the name given to organisms whose cells contain different compartments separated by membranes, such as the nucleus of the cell. For years, we have assumed that eukaryotic microbes, like fungi and parasites, act more like humans than like bacteria because work in budding yeast (another eukaryote) has suggested this to be the case. However, recent work in disease-causing fungi has shown that, much like bacteria, elevated mutation rates may help them to respond to stress. This could also enable fungi to become resistant to drugs used to treat fungal infections.

Cryptococcus deuterogattii is a fungus that causes human diseases including meningoencephalitis and a lung infection called pulmonary cryptococcosis. An ongoing outbreak of the fungus began in the Pacific Northwest of Canada in the late 1990s and emerged in the United States in 2006/2007. Among isolates closely related to those fungi causing the outbreak, three were found that appear to have a specific mutation in their DNA mismatch repair pathway, meaning that they may also experience a higher mutation rate. These strains are also less able to cause disease than others.

Billmyre et al. now demonstrate experimentally that all three isolates have a specific DNA mismatch repair defect, and show that these fungi experience elevated mutation rates, resulting in what is known as a hypermutator state. Furthermore, whole genome sequencing and phylogenetic analysis showed that these hypermutator strains are derived from the outbreak-causing fungi, and that their reduced ability to cause disease is likely a result of accumulating mutations and the loss of the ability to grow at the higher temperatures found in the human body.

Fungal infections are difficult to treat, in part because there are a limited number of available drugs. Elevated mutation rates will likely increase how often and how rapidly fungi develop resistance to these drugs. Understanding how commonly fungi exhibit a hypermutator state that could impact the development of drug resistance will therefore be important for treating patients with fungal infections, which account for millions of infections and hundreds of thousands of deaths annually worldwide.

RNA Interference in Fungi: Retention and Loss

RNA interference (RNAi) is a mechanism conserved in eukaryotes, including fungi, that represses gene expression by means of small noncoding RNAs (sRNAs) of about 20 to 30 nucleotides. Its discovery is one of the most important scientific breakthroughs of the past 20 years, and it has revolutionized our perception of the functioning of the cell. Initially described and characterized in Neurospora crassa, the RNAi is widespread in fungi, suggesting that it plays important functions in the fungal kingdom. Several RNAi-related mechanisms for maintenance of genome integrity, particularly protection against exogenous nucleic acids such as mobile elements, have been described in several fungi, suggesting that this is the main function of RNAi in the fungal kingdom. However, an increasing number of fungal sRNAs with regulatory functions generated by specific RNAi pathways have been identified. Several mechanistic aspects of the biogenesis of these sRNAs are known, but their function in fungal development and physiology is scarce, except for remarkable examples such as Mucor circinelloides, in which specific sRNAs clearly regulate responses to environmental and endogenous signals. Despite the retention of RNAi in most species, some fungal groups and species lack an active RNAi mechanism, suggesting that its loss may provide some selective advantage. This article summarizes the current understanding of RNAi functions in the fungal kingdom.

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.

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.