Next Article in Journal
Cellular Basis of Embryonic Hematopoiesis and Its Implications in Prenatal Erythropoiesis
Next Article in Special Issue
The RNAi Mechanism Regulates a New Exonuclease Gene Involved in the Virulence of Mucorales
Previous Article in Journal
Anti-Stem Cell Property of Pterostilbene in Gastrointestinal Cancer Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

The Evolutionary Significance of RNAi in the Fungal Kingdom

by
Carlos Lax
1,
Ghizlane Tahiri
1,
José Alberto Patiño-Medina
2,
José T. Cánovas-Márquez
1,
José A. Pérez-Ruiz
1,
Macario Osorio-Concepción
1,
Eusebio Navarro
1 and
Silvia Calo
3,*
1
Department of Genetics and Microbiology, Faculty of Biology, University of Murcia, 30100 Murcia, Spain
2
Instituto de Investigaciones Químico Biológicas, Universidad Michoacana de San Nicolás de Hidalgo, Ciudad Universitaria, Morelia, Michoacán CP 58030, Mexico
3
School of Natural and Exact Sciences, Pontificia Universidad Católica Madre y Maestra, 51033 Santiago de los Caballeros, Dominican Republic
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(24), 9348; https://doi.org/10.3390/ijms21249348
Submission received: 31 October 2020 / Revised: 23 November 2020 / Accepted: 24 November 2020 / Published: 8 December 2020
(This article belongs to the Special Issue Functional and Mechanistic Diversity of RNA Silencing in Fungi)

Abstract

:
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.

Graphical Abstract

1. Introduction

RNA interference (RNAi) or RNA silencing has been deeply studied in the last two decades, as its discovery entailed a revolution in the understanding of the regulation of gene expression. This RNAi pathway, broadly conserved in eukaryotes (Table 1), uses small interfering RNAs (siRNAs) to suppress gene expression of homologous sequences. These siRNAs, of 20–30 nucleotides (nt) long, are produced from double-stranded RNA (dsRNA) by an RNase III called Dicer (Dcr) and loaded into an RNA-induced silencing complex (RISC), which contains an Argonaute protein (Ago) that drives the selective degradation of homologous messenger RNAs (mRNA), as well as translational or transcriptional repression of target sequences. Moreover, in fungi and other organisms, an RNA-dependent RNA polymerase (Rdp or RdRP) generates dsRNA from certain single-stranded RNA (ssRNA) or from the target messenger RNA, activating or amplifying the silencing response, respectively [1,2].
Fungi have proven to be excellent model organisms for the study of the RNAi pathway, since many of the discoveries accomplished in these organisms were later extended to higher eukaryotes. In fact, one of the first RNA silencing phenomena reported was found in Neurospora crassa, which is an essential model organism to study modern genetics. N. crassa has developed different RNAi mechanisms, but the two that were originally found are the best described. The first is called quelling, a post-transcriptional gene silencing (PTGS) guided by siRNAs [15], which suppresses transposons and virus infections [16]. Quelling is triggered by the introduction of transgenes homologous to an endogenous gene. After the transgene is transcribed, it follows a canonical pathway to activate silencing, carried out by an Rdp protein (QDE-1), two Dicer-like proteins (Dcl1 or Dcl2), and an Argonaute (QDE-2). The second mechanism is called meiotic silencing of unpaired DNA (MSUD) and is involved in silencing genes that are not paired with their partner on the homologous chromosome during meiosis [17]. This mechanism is present not only in N. crassa but also in other ascomycetes, such as Gibberella zeae, and operates during prophase I [18,19]. Some of the elements involved in the canonical RNAi pathway, such as Dcl1, are necessary for this other mechanism, as well as a MSUD-specific Rdp (SAD-1), a second Argonaute (SMS-2), and the helicase SAD-3. These proteins form a multiprotein complex located at the perinuclear region that acts generating MSUD-associated siRNAs (masiRNAs) [20]. After these first discoveries in N. crassa, the RNAi mechanism was found in several other fungi, such as Schizosaccharomyces pombe [21], Cryptococcus [22], or Mucor [23].
When RNAi was discovered, it was thought to be a defense system against exogenous and potentially harmful DNA, including transposons, virus, and transgenes. However, very soon, its involvement in other cellular functions, such as genome integrity or gene regulation, was found. Recent studies in fungal pathogens stated that the RNAi pathway is also implicated in development, drug tolerance, and virulence. Thus, fungi have exploited RNAi to tune their cellular processes, reaching unsuspected limits. The evolutionary relevance of the retention or loss of this regulatory pathway, as well as its diversity of functions in fungi (Table 2), is evaluated in this review.

2. Defense against Viruses

The first predicted function of RNAi as a defense against viral infections would explain the conservation of the pathway through the evolution of eukaryotes, since viral infections affect every phylum on the tree of life. RNAi drives virus inactivation through the production of virus-derived small interfering RNAs (vsRNAs) originated from viral dsRNA [70]. This RNAi-based mechanism has been widely explored in the ascomycete filamentous fungus Cryphonectria parasitica, which is considered a model for the analysis of virus–fungus interactions because it can be infected by five different virus families. Furthermore, several studies have shown that RNAi regulates virus infection in this fungus [24,25]. C. parasitica contains two dicer-like and four Argonaute-like genes. Dcl2 and Ago2 are the primary components of the antiviral defense response, since mutants defective in dcl2 or ago2 genes are highly susceptible to mycovirus infection and are affected in the generation of vsRNAs [24,25]. Moreover, there is an induction of dcl2 and ago2 expression regulated at the transcriptional level by viral infection [24]. In other fungi, such as N. crassa, dsRNA also induces the expression of RNAi genes, supporting the protective role of the RNAi mechanism against virus [71]. On the other hand, viruses have evolved to express RNAi suppressors [72]. The mycovirus Cryphonectria hypovirus 1 (CHV1), which can also infect C. parasitica, expresses a Papain-like protease. This protein, known as p29, inhibits the expression of dcl2 and ago2 in the fungus, thereby avoiding immunity [24].
Furthermore, the mycoviruses that infect the filamentous fungi of the genus Aspergillus have been extensively studied, contributing to our understanding of viral diversity and leading to numerous discoveries about virus infection in fungi. In addition to its impact in agriculture and food industry, Aspergillus has great relevance in medicine as the causative agent of the fungal infection Aspergillosis. In numerous important Aspergillus species, such as A. nidulans, A. fumigatus, A. niger, and A. flavus, RNAi seems to be involved in virus–fungus interactions [73]. Specifically, in A. nidulans, siRNAs derived from Aspergillus virus 341 were detected, indicating that RNAi acts as a defense mechanism against the virus [26]. A. nidulans has genes encoding one Argonaute and one Dicer. Although two rdrp genes are found in the fungus genome, neither is required for small RNA (sRNA) production, suggesting the absence of silencing amplification [1]. Similarly to viruses that infect C. parasitica, infection of A. nidulans with Aspergillus virus 1816 suppresses RNAi, suggesting the existence of an RNA silencing suppressor encoded by the viral genome [26]. Thus, evolution of the virus–host interaction includes tuning of the host RNAi mechanism for both the virus and the host’s benefit.

3. Control of Transposable Elements

Transposable elements (TEs) are described as DNA sequences that have the ability to change their position within a genome. For details on TE classification, we encourage readers to refer to the review by Bourque et al. [74]. Although TEs were initially considered junk DNA, they have been associated with several important activities since their discovery, including centromere function, genome reorganization, and gene expression regulation [75,76,77]. TEs are also considered “selfish” DNA because they try to be perpetuated whilst the host tries to curtail their spread and, thus, their consequences on genome integrity. As a result, many organisms have developed mechanisms to ensure the control of TE activity [78,79]. Some of those mechanisms are well described in the literature and include DNA methylation [80], histone methylation [81], and heterochromatin-inducing protein [82]. Another mechanism that controls TE spread, probably the most ancient of all mentioned, is RNAi. Small RNAs associated with proteins can act at the transcriptional or post-transcriptional level against TE activity. The role of RNAi in TE repression has been well characterized in the plant kingdom [83] and other organisms, such as Drosophila melanogaster and Caenorhabditis elegans [84,85]. Noteworthy, in C. elegans and Chlamydomonas reinhardtii, some genes have proven to be essential for both RNAi and TE control mechanisms [86,87,88], while, in mammals, the existence of active siRNA molecules against LINE-1 transposons (Class 1, non-LTR) has been reported [89]. Studies on the RNAi role in TE control in fungi have served to further understand these mechanisms. In N. crassa, the repression of LINE1-like retrotransposons is exclusively achieved by post-transcriptional gene silencing and, unlike other eukaryotic species, DNA methylation appears not to be involved [27]. This TE control mechanism shares common elements with quelling (QDE-2 and Dicer) but is independent of QDE-1 or QDE-3, suggesting a diversification of pathways depending on the target. In addition to TE control during mitotic growth, MSUD acts by suppressing repetitive elements during meiosis [19,28].
RNAi mechanisms for TE regulation have also been described in several other fungi. Yamanaka et al. [29] reported the existence of siRNAs targeting the Tf2 retrotransposon in S. pombe, where a correlation between siRNAs and methylation of histone H3 in lysine 9 (H3K9me) was also detected. In the phytopathogenic fungus Magnaporthe oryzae, approximately 10% of all the small RNAs detected in the mycelial stage mapped to repetitive DNA, especially LTR-retrotransposons [30], and RNAi components were proven to be in control of the LTR-retrotransposon MAGGY [31,32]. Previous studies in M. oryzae already confirmed that DNA methylation was not implicated in MAGGY repression and suggested the influence of gene silencing [33]. In other fungi with active RNAi machinery, the possible role of RNAi in TE control has also been suggested, i.e., Aspergillus [90], or discarded, as occurred in the case of Trichoderma atroviride due to the very low number of TEs present on its genome and the absence of small RNAs mapping to these regions [44].
In Mucor lusitanicus, formerly known as Mucor circinelloides f. lusitanicus [91], sRNAs pairing with transposons and derived from the canonical pathway have been detected (see below for a more detailed description of RNAi pathways of this fungus) [43,92]. A recent study discovered the presence of TE sequences surrounding the centromeric area of M. lusitanicus, which seem to be present in all Mucoromycotina species and, thus, were termed as Grem-LINE1 (genomic retrotransposable element of Mucoromycota) [34]. The canonical RNAi pathway negatively regulates theses sequences, as mutant strains lacking RNAi typical components (Ago1, Dcl1, and Dcl2) show a decrease in sRNA levels targeting Grem-LINE and, consequently, increased Grem-LINE mRNA expression. Surprisingly, mutant strains lacking proteins (RdRP1 and R3B2) of a different pathway known as noncanonical RNAi pathway (NCRIP) show enhanced levels of Grem-LINE1 sRNAs compared to the wild-type strain [64]. Therefore, the NCRIP regulates the silencing of these Grem-LINE1 through the modulation of the canonical mechanism [64].
The pathogenic fungus Cryptococcus neoformans has also developed a mechanism to maintain genome integrity by controlling TE activity via RNAi, especially during sexual reproduction. Its sequenced genome revealed the presence of genes encoding the core components of RNAi: one Rdp (Rdp1), two Dicers (Dcr1 and Dcr2), and one Argonaute (Ago1) [93]. Furthermore, a significant accumulation of sRNAs mapping to TE sequences was found in strains undergoing sexual reproduction. However, in a Δrdp1 mutant strain, these sRNAs were lacking and a higher expression of retrotransposons was reported [35,36]. Thus, this RNAi mechanism was named sex-induced silencing (SIS) and is significantly more active than silencing during mitotic growth [36]. Moreover, by using the FKBP12-encoding gene FRR1 as a transposon trap, the Δrdp1 mutant strain showed an increased transposition rate during both opposite sex (a-α) and unisexual reproduction (α-α), which evidences the strategic role of SIS as a TE control mechanism [37]. Moreover, the process underlying target sequence recognition in this basidiomycete has been studied in detail. Transcripts targeted by RNAi encode suboptimal splicing signals that result in an abnormally high occupancy of spliceosomes. Then, the mRNA stalled in the spliceosomes is used for siRNAs production by a spliceosome-coupled and nuclear RNAi complex (SCANR), as well as the lariat debranching enzyme (Dbr1) [38]. These siRNAs are then used to silence TE. Therefore, the presence of intron-containing transposon-derived mRNA in those spliceosomes evidences TE silencing.
Recently, a set of five novel genes (rde1–5) required for the repression of the DNA transposons HAR1 was also identified in C. neoformans [39]. Rde1 and Rde2 proteins appear to be associated with and located in the nucleolus. Rde1 also associates with Ago1 and Prp43, a helicase involved in disassembling stalled and post-catalytic spliceosomes. Rde4 presents a predicted terminal nucleotidyl transferase domain that resembles polyA polymerases and terminal-uridylyl transferases (TUTases). Since Rde4 is implicated in siRNAs biogenesis, it may be possible for TUTases to have a role of TE control in this fungus, similarly to LINE-1 repression in humans [94]. Interestingly, Rde3, a non-Dicer RNAse III related to Saccharomyces cerevisiae Rnt1 protein, could be interacting with Rde5 (function unknown) and with RNA quality control enzymes, such as Rrp6 (component of the nuclear exosome) and Rnh1 (RNase H that cleaves RNA-DNA duplexes) [39]. Moreover, in addition to siRNAs, microRNAs (miRNAs) targeting transposons and pseudogenes have also been detected in Cryptococcus, suggesting a possible participation of these molecules in TE control [95].
Lastly, analogously to M. lusitanicus, centromere sequences in C. neoformans, C. deneoformans, and C. deuterogattii are swamped by transposons [96]. Therefore, by controlling the expression of TEs at centromeres, RNAi could play an important regulatory role in the structural evolution of centromeres. Likewise, centromeres of Puccinia graminis f. sp. tritici (Pgt) are highly repetitive [40]. Pgt produces a strong wave of sRNAs derived from centromeric regions in the later stages of infection, probably to maintain genome stability during cell division. These centromeric sRNAs induce silencing of TEs both inside and outside the centromeres, as well as affect neighboring genes [40]. A similar phenomenon occurs in other phytopathogen species, where the surrounding of effector genes (term used when referring to genes implicated in the infection process) by repetitive sequences has been reported. RNAi suppresses not only the transposon sequences but also the effector genes, allowing fungi to avoid detection during infection, as plants deploy intracellular immune receptors involved in the recognition of those effector proteins [97,98]. On the other hand, the presence of transposons flanking effector genes may also confer an evolutive advantage. Due to the mobile nature of TE, the effector genes could be translocated to highly dynamic regions of chromosomes, modifying their regulation and allowing rapid adaptation to changing conditions [99].

4. Regulation of Endogenous Genes

As described above, the development of RNAi mechanisms represents an evolutive advantage regarding the defense against exogenous nucleic acids. However, specialization of those RNAi mechanisms has led to the establishment of novel post-transcriptional regulation networks of endogenous genes according to the use of Rdp, Dicer, and Argonaute proteins. The human pathogenic fungus M. lusitanicus, for instance, shows an intricate RNAi mechanism as a function the interplay of the silencing proteins in three different pathways, named the canonical, epimutational, and noncanonical RNAi pathways. The crosstalk of the RNAi pathways creates a complex network that regulates both basic cellular activities, such as metabolism or vegetative growth, and elaborated mechanisms, including sexual reproduction and pathogenesis [23,100,101].
M. lusitanicus encodes two Dicer-like (Dcl1-2), three Argonaute (Ago1-3), and three Rdp (RdRP1-3) proteins. RdRP1 is essential for the generation of antisense RNA from transgene transcripts. Triggering dsRNA molecules produced by RdRP1 are then processed by Dcl2 into two different classes of siRNAs, 21 and 25 nt long. Afterward, Ago1 binds to the siRNA molecules and provokes degradation of target transcripts. Lastly, RdRP2 amplifies the silencing signal. This basic RNAi machinery acts as a defense mechanism against exogenous nucleic acid sequences, such as transgenes, viruses, and transposons (see above), but also regulates the expression of endogenous genes through exonic siRNAs (ex-siRNAs), which are produced from exons. This new type of siRNA has a post-transcriptional effect on mRNA [92] and regulates the expression of the genes they are produced from and other genes in trans, with complementary sequences or regulated by the initial targets [41]. Four different classes have been described for the 324 exonic siRNAs found so far, depending on the proteins involved in their biogenesis: class I (Dcl2 and RdRP2), class II (Dcl2 and RdRP1), class III (Dcl1, Dcl2, and RdRP1 or RdRP2), and class IV (Dcl1 and RdRP1 or RdRP2) [92]. The protein Ago1 is implicated in all four classes, for their biogenesis and/or stability [102], and the ribonuclease III, R3B2, participates in the generation of most ex-siRNAs of class II, III, and IV, and some of class I [43]. Nearly 700 genes are regulated by these ex-siRNAs, which are involved in several functions, such as carbohydrate, lipid, or amino-acid metabolism, energy production and conversion, and even defense mechanisms. Thus, ex-siRNAs have an important role in M. lusitanicus development and physiological processes [41]. In fact, mutant strains affected in components of the silencing machinery required for biogenesis of ex-siRNAs show reduced asexual sporulation [41,102,103], growth, and altered hyphal morphology [104], and accelerated autolysis as a response to nutrient starvation [102]. Likewise, in T. atroviride, sRNAs similar to M. lusitanicus ex-siRNAs are essential for growth and development [44]. Fusarium graminearum also produces ex-siRNAs, which are involved specifically in the ascospore formation [45].
The most recent RNAi mechanism discovered in M. lusitanicus is known as NCRIP because none of the known Dicer and Argonaute proteins are involved. A new group of sRNAs, named rdrp-dependent dicer-independent sRNAs (rdRNAs), was found during the analysis of the double Dicer mutant (Δdcl1/dcl2) of M. lusitanicus [43]. The rdRNAs are degradation products of the mRNA without a defined size but with a strong bias for uracil in the penultimate position. Furthermore, they are produced by R3B2, an atypical RNase III-like protein with one RNase III domain and two double-strand RNA-binding domains, together with RdRP1, RdRP3, and the helicase RnhA. Thus, these new sRNAs are not random degradation products of the fungal transcriptome, and NCRIP drives its genetic control through the degradation of the mRNA from its target genes [42,43]. This mechanism was initially described as a gene regulation pathway implicated in oxidative stress and sexual interaction [23,43]. Additionally, the latest studies established that NCRIP regulates more than 25% of the genome, 3187 genes mostly linked to cellular metabolism, germination, and development, indicating a relevant role in saprophytic growth conditions [64]. Furthermore, the transcriptomic profile of M. lusitanicus during phagocytosis is similar to that of mutant strains lacking NCRIP [64,101], suggesting that this mechanism is repressing the genetic program of defense against phagocytosis. Thus, the NCRIP also has a relevant role in the pathogenesis of this filamentous fungus (see below for a more detailed description). On the other hand, the epimutational RNAi pathway, which is also reviewed more thoroughly later on, allows in M. lusitanicus the transitory silencing of specific mRNA during a stress challenge [42,55,56]. The fact that the inactivation of the NCRIP increments the production of epimutants shows the competition of both RNAi pathways over the messengers with high expression rates [42,43]. This, together with the phagocytosis results, exposes again the relevant role of the NCRIP repressing the defense genetic programs during optimal growth of M. lusitanicus. Lastly, the R3B2 protein, which produces rdRNA, is also required by the canonical RNAi pathway of M. lusitanicus. As described before, the atypical RNase III is involved in the production of most of the ex-siRNAs, together with RdRP, Dicer, and Argonaute proteins [41,43]. The relationship of the NCRIP with the canonical RNAi and the epimutation pathways reveals the intricate interaction of proteins developed through evolution that allows M. lusitanicus refined control over mRNA levels to confront environmental challenges [64,100].
On another note, in N. crassa, Dicer-independent small interfering RNAs (disiRNAs) emerged as a new type of silencing pathway, without the implication of the classic RNAi proteins for its generation [48]. High levels of disiRNAs result in high levels of DNA methylation and promoter-specific RNA transcripts, suggesting a regulation of endogenous gene expression by DNA methylation, which was, thus, named disiRNA loci DNA methylation (DLDM) [50]. This mechanism has a high dynamism and an on/off pattern, where antisense transcription provokes the stalling of RNA polymerase II and the exonuclease ERI-1 has a main role in both DLDM and disiRNA biogenesis [49]. On the other hand, microRNAs (miRNAs) are small RNAs which are present in many organisms, regulating the expression of endogenous genes by complementarity with their target mRNAs [105]. N. crassa produces microRNA-like RNAs (milRNAs), which have significant similarities with miRNAs related to their precursors, their production by Dicer, and their function regulating cellular activities [48]. Different combinations of Dicer proteins, the exonuclease QIP (quelling inducer protein), the RNase III domain-containing protein MRPL3, and the Argonaute QDE-2, are involved in the diverse mechanisms (at least four) that produce these milRNAs [48]. In addition, it has been demonstrated that RNA polymerase II and III participate in their transcription. This remarkable role of RNA polymerase III in the production of milRNAs seems to be one of the differences between milRNAs and miRNAs [106]. milRNAs have also been studied in other organisms. Coprinopsis cinerea, a model mushroom, produces milRNAs implicated in development and fruiting body formation [47]; moreover, in F. graminearum, milRNAs were initially found not only in asexual phases [107] but also in sexual stages later on [46], where they seem to have an important role in the regulation of sexual development. As described above, F. graminearum also produces ex-siRNAs [45]. The participation of Argonaute and Dicer proteins is essential for the biogenesis of both kinds of sRNAs in this fungus, FgDcl1 and FgAgo2, to generate ex-siRNAs, and FgDcl1/2 and FgAgo1 to generate milRNAs [45,46].
In Magnaporthe oryzae, two different sRNAs have been observed: transfer RNA fragments (tRFs) and circular RNAs (circRNAs). tRFs have been detected in bacteria, fungi, plants, and animals [108]. In M. oryzae, tRFs were identified in spores and mycelia, but their specific function has not yet been identified [30]. These tRNA fragments are involved in many different cell activities in human and other organisms, including the inhibition of translation [109], repression of endogenous genes involved in cell proliferation and DNA damage response [110], and even the suppression of antiviral responses by a respiratory syncytial virus [111]. tRFs have also been observed in the plant pathogen Phytophthora infestans, produced with the implication of PiAgo1, but not PiDcl1. High levels of tRFs in this fungus are linked to plant infection [52]. On the other hand, circRNAs, which repress miRNAs in animals and plants, were recently discovered in M. oryzae. circRNAs have been studied in mycelium and conidium of this fungus, and they are produced from genes involved in metabolism and normal growth, as well as in the biogenesis of lipid, glycogen, and amino-acid storages, which may be used for plant infection. However, circRNAs are still under study, because their general functions are not entirely known [51].
Thus, the huge variety of RNAi-based pathways that regulate endogenous gene expression in fungi supports the idea that RNAi appeared early in evolution. RNA silencing may even have played an important role in the differentiation of eukaryotes, as it allows a fine-tuning of gene expression, which would change drastically with few alterations in the implicated proteins.

5. Heterochromatin Formation

Heterochromatin constitutes a highly condensed state of DNA. It is considered to have no transcriptional activity due to the limited access of the regulatory proteins to the promoter regions [112]. Generally, heterochromatin is concentrated in the telomeric, centromeric, ribosomal, and mating type regions of the eukaryotic chromosome [113]. Heterochromatin assembly is strictly regulated for accurate chromosome segregation, maintenance of telomere integrity, transcriptional silencing, and transposon control [112,114]. S. pombe is an excellent study model to discern and understand the epigenetic mechanism regulating heterochromatin assembly [21]. Furthermore, S. pombe’s genome presents a single copy of the genes encoding RNAi proteins: Argonaute (Ago1), Dicer (Dcr1), and RNA-dependent RNA polymerase (Rdp1) [53,115]. Heterochromatin formation in S. pombe is triggered by the production of siRNAs derived from centromeric regions with numerous repeats [53,115]. These siRNAs are loaded into the RNA-induced transcriptional silencing complex (RITS) through the Argonaute protein, leading the complex to complementary nascent RNA transcripts from the initial centromeric repeats [116,117]. Chromatin-associated RITS promotes the recruitment of the histone methyltransferase Clr4, a subunit of the cryptic loci regulator complex (CLRC), which is involved in nucleosome methylation and propagation of heterochromatin. The Clr4 enzyme, therefore, catalyzes the methylation of the histone H3 (H3K9me) at the centromeric repeats, inducing heterochromatin assembly [115,117,118,119]. Absence of Clr4 results in the interruption of both events, histone methylation and heterochromatin silencing [117,120]. Additionally, the H3K9me marks at the centromeric regions function as binding sites for chromodomain proteins, such as Chp1 (a component of RITS), Swi6 (a heterochromatin protein 1 homolog), and Clr4 itself. Consequently, H3K9me promotes the recruitment of chromatin-modifying proteins and RNAi components, which stimulate siRNA amplification and the spread of H3K9me domains to adjacent centromeric regions [53,118,121]. Thus, the coupling between RITS/RNAi and CLRC functions is essential for the formation and maintenance of heterochromatin, as well as centromeric silencing in S. pombe [121,122]. Nevertheless, heterochromatin formation to maintain epigenetic silencing is a highly complex process and could integrate alternate pathways, as cells lacking the RNAi components maintain basal levels of H3K9me and drive the formation of heterochromatin [123].
Remarkably, heterochromatin formation and the RNAi pathway can also regulate the epigenetic inheritance of gene silencing in S. pombe [54]. siRNA amplification derived from hairpin RNA promotes the establishment of the H3K9 methylation at a target euchromatic locus through the recruitment of chromatin-modifying enzymes, leading to heterochromatin formation and gene silencing. Once established, the H3K9me-mediated epigenetic silencing can be stably inherited through multiple generations, even in the absence of the primary RNA [54,119,124]. The H3K9me contained in the parental histones is recognized by chromatin regulatory proteins, which spread the methylation to newly incorporated histones at the target locus [54,117]. Hence, it has been proposed that H3K9me can act as a molecular mark and confer a cellular memory of gene silencing events in past generations [124]. However, this phenomenon only occurs under the same conditions in which the original silencing was acquired. Thus, the epigenetic memory may promote survival of this species in hostile environments [124]. In conclusion, the development of a crosstalk between the RNAi and the H3K9me epigenetic mechanism has been crucial in the evolution of gene silencing in S. pombe [54,119].

6. Adaptation to Stressful Conditions

The epimutational pathway in M. lusitanicus was discovered after the emergence of isolates resistant to the antifungal drug FK506 with no apparent mutations in the target genes [55]. FK506 interacts with the peptidylprolyl isomerase FKBP12, forming a complex that inhibits Calcineurin, a phosphatase activated by calcium-calmodulin necessary for dimorphism and virulence of M. lusitanicus. Thus, inhibition of Calcineurin after treatment with this drug enforces Mucor to grow as yeast [125]. However, exposure to FK506 produced the emergence of several Mucor isolates with mycelial growth. Most of these FK506-resistant isolates harbored mutations in the fkbA gene, which encodes FKBP12, or in the cnaA or cnbR genes, which encode the subunits of Calcineurin, and presented stable drug resistance. However, several FK506-resistant strains did not have mutations in any of those genes. They also showed resistance to rapamycin, which has the same FKBP12 target, but were not resistant to any other tested drug. These FK506-resistant isolates had a complete loss of fkbA expression, but fkbA levels were restored to wild-type levels, and the isolates became drug-sensitive after several generations of vegetative growth on FK506-free media. Thus, this process was reversible, and it appeared solely epigenetic. The isolates, called epimutants, produced siRNAs from the mature mRNA of fkbA gene, with an average of 21–24 nt in length and a prevalence of uracil at the 5′ end. Both are typical features of Mucor siRNAs that interact with Argonaute proteins [55,56]. Epimutants developing resistance to other antifungal agents, such as 5-fluoroorotic acid (5-FOA), have also been isolated. In this case, siRNAs were produced from the exonic regions of the pyrG or pyrF genes, which encode the proteins that transform 5-FOA into a toxic compound for the fungus. Again, the mRNA levels from the target gene were significantly lower in the epimutants than in the wild-type strain [56]. Therefore, the epimutation process does not appear to occur at a specific gene locus, suggesting it might constitute a general mechanism that generates phenotypic plasticity in Mucor by silencing key genes and allowing rapid and reversible adaptation to environmental stresses [42,126].
This mechanism of epimutation requires the action of the canonical RNAi pathway, since the deletion of some components of this pathway, such as the Dicer (Dcl1 or Dcl2), Argonaute (Ago1), or RNA-dependent RNA polymerase (RdRP2) proteins, as well as QIP or a Sad-3-like Helicase (RnhA), abolished the ability to produce drug resistance through epimutations [42,55]. On the other hand, deletions of rdrp1, rdrp3, or r3b2 genes increased the production of drug-resistant epimutants. Thus, the NCRIP inhibits the mechanism of generation of epimutants [42]. The implication of both RNAi pathways (canonical and NCRIP) in the epimutant process indicates the existence of a balanced regulation of both mechanisms that produces the activation of the RNAi-dependent epimutant pathway under stress conditions or its repression when the mRNA degradation mechanism operates under nonstress conditions [42]. In addition, Mucor circinelloides strains (formerly known as Mucor circinelloides f. circinelloides and a close relative of M. lusitanicus) isolated from human patients displayed an overly activated epimutation pathway, suggesting that the selective pressure of the host–pathogen interaction might promote the phenotypic plasticity provided by this mechanism [42,55]. As NCRIP inhibits the generation of epimutation, it may be interesting to analyze if this mechanism is repressed in these strains.
Moreover, after the infection of immunocompetent mice with an epimutant FK506-resistant Mucor strain, epimutation was found to be maintained in vivo during infection. However, reversion of epimutation to a wild-type drug-sensitive phenotype was markedly higher in isolates recovered from the brain than from other evaluated organs [127]. Furthermore, in vivo passage of a wild-type strain of Mucor led to increased epimutant-driven induction of drug resistance, especially in isolates collected from the brain. This induction of epimutation did not occur during the infection itself but aroused several days after the recovery of the strains from the different organs and under exposure to drug selection. Therefore, the brain environment does not appear to be the direct inducer of epimutation, but rather prepares the fungus to respond quickly to stress, including antifungal drug exposure [127]. Mucor is a ubiquitous fungus that can be found mostly in soil, but is also able to infect and grow in both invertebrate and vertebrate hosts. The adaptation to this variety of environments may be favored by this quick change at the epigenetic level [127].

7. Pathogenesis

RNAi has also been found to play an important role in pathogenesis, more thoroughly studied in plant pathogens. Many crops with worldwide importance are susceptible to being infected by pathogenic fungi, which translates into economic losses. Thus, alternative methods of infection control have been investigated, allowing a deeper understanding of fungal pathogenesis and the involvement of RNAi. Some of those pathogenic fungi, such as Colletotrichum gloeosporioides or M. oryzae, have active RNAi pathways implicated in growth and conidiation, which influences their pathogenicity. Mutants lacking one or two RNAi genes in those species usually lose their ability to penetrate the leaves [57,58]. Similarly, Sclerotinia sclerotiorum, a devastating plant fungal pathogen that causes up to 100% yield losses in crop production and affects a wide array of crops, encodes two dicer (dcr1-2) and two argonaute (ago2,4) genes that have been functionally characterized. RNA-silencing-deficient mutants, Δago2 mutant [128] and Δdcr1/dcr2 double mutant [59], displayed significantly debilitated growth and reduced virulence. Other species, such as Valsa mali, which severely damages apple production, and Fusarium oxysporum f. sp. lycopersici, which causes wilt disease in tomato, show significantly decreased virulence when they lack the Argonaute gene ago2 [60,61]. In Penicillium italicum, a citrus fruit pathogen that generate citrus blue mold, two Dicer genes (dcr1-2) have been characterized. Silencing of the dcr2 gene, but not dcr1, incapacitates the fungus to infect [62]. In the same way, the plant fungal pathogen F. graminearum, the major causal agent of wheat head blight, encodes two Dicer, two Argonaute, and five Rdp proteins [107]. FgDcr1 and FgAgo2, which were essential proteins for the generation of ex-siRNAs, are also positive regulators for F. graminearum virulence, since mutation of the encoding genes resulted in reduced Fusarium head blight development. Moreover, F. graminearum produces several mycotoxins that help the fungus to colonize the plant tissues. When the concentration of some of the mycotoxins was analyzed in wheat spikes infected with RNAi mutants, they showed a lower accumulation than in spikes infected with the wild-type, suggesting that the fungal RNAi mechanism may promote F. graminearum mycotoxin production [63].
RNAi is also involved in virulence of animal pathogens. In M. lusitanicus, knockout mutants in rdrp1 and r3b2, two NCRIP key proteins, showed an augmented oxidative stress tolerance in vitro and a reduced pathogenic potential in vivo, when tested in an immunosuppressed mice model. Analysis of the transcriptomic profiles of NCRIP mutants identified a substantial number of differentially expressed genes in both Δrdrp1 and Δr3b2 mutants, compared to the wild-type strain. Remarkably, most of the fungal genes regulated by phagocytosis (involved in functions such as cell motility, cytoskeleton, inorganic ion transport, and metabolism) are under the control of NCRIP, suggesting that this RNAi-based mechanism is a master regulator of the response of the pathogen to phagocytosis and, thus, essential for virulence [64].
An important feature of fungal pathogenesis is the mechanism called cross-kingdom RNAi, which has evolved to regulate the host–pathogen interaction. The existence of sRNA trafficking between the host and the pathogen and silencing target genes of the counterparty in trans was first discovered in plants, but afterward extended to mammal systems. Some fungal pathogens have been found to produce sRNAs that function as RNA effectors to suppress host immunity. Phytophthora sojae, an oomycete soybean pathogen, produces two effectors, PSR1 and PSR2, that target the host RNAi processes to attenuate the plant’s immune response and to promote the infection. PSR1 inhibits the host biogenesis of both miRNAs and siRNAs, whereas PSR2 specifically impairs siRNAs accumulation [65]. Likewise, Botrytis cinerea, a necrotrophic fungal pathogen, produces small RNAs (Bc-sRNAs) during infection that hijack the host plant RNAi machinery to induce the silencing of host immunity genes in tomato and the Arabidopsis thaliana B. cinerea dcl1/dcl2 double mutant that could no longer produce Bc-sRNAs exhibited reduced virulence, whereas the A. thaliana ago1 mutant, which lost RNAi function, developed resistance to B. cinerea [66]. Thus, both the pathogen and the host’s RNAi mechanisms are crucial for this cross-kingdom phenomenon to work. Similarly, Puccinia striiformis f. sp. tritici (Pst), one of the most destructive pathogens of wheat, produces a microRNA-like RNA 1 (Pst-milR1), which suppresses wheat defenses during wheat–Pst interaction. Silencing of the Pst-milR1 precursor resulted in increased wheat resistance to Pst infection [67]. Moreover, in the cereal powdery mildew pathogens Blumeria graminis f. sp. hordei (Bgh) and Blumeria graminis f. sp. tritici (Bgt), sRNA-seq data derived from infected wheat revealed six putative sRNAs from Bgt and 15 putative sRNAs from Bgh. Those sRNAs had predicted targets exclusively present in plants, with associated functions that would alter primary metabolism, suggesting potential cross-kingdom RNA transference [68]. Lastly, in the oomycete biotrophic pathogen Hyaloperonospora arabidopsidis, 34 sRNAs with the capacity to translocate into plant cells and suppress host target genes were found. However, Arabidopsis siRNA biogenesis mutant strains displayed increased H. arabidopsidis growth, demonstrating the important role of siRNAs in plant immunity [69]. Therefore, the host–pathogen interaction is constantly evolving, and RNAi has become an advantage for both the host and the pathogen. A fragile balance will decide which succeeds, infection or immunity.

8. Loss of RNAi

In essence, RNAi has crucial regulatory and defense roles in eukaryotes, suggesting that this key mechanism has been positively selected through evolution in plants, nematodes, animals, and fungi. Yet, some members of the fungal kingdom have lost key components of the RNAi pathway [129], resulting in its inactivation (Table 3). A hypothesis to explain this contradiction could be that those species may have other defensive mechanisms more advantageous than RNAi. Alternatively, perhaps, this RNA-based mechanism constitutes a disadvantage for them, forcing the survival of the RNAi-deficient species. Loss of RNAi affects species of the subphyla Ustilaginomycotina, Saccharomycotina, Wallemiomycetes, and some members of the phylum Microsporidia [130]. Those organisms lack key components of the RNAi machinery: Dicer, Argonaute, and Rdp proteins. Ustilago maydis, for instance, has lost genes encoding the main components of this RNA degradation pathway, even when related species such as Ustilago hordei have a functional RNAi mechanism [131]. Furthermore, genes involved in heterochromatin formation, DNA methylation, and repeat-induced point mutation (RIP), a mechanism that inhibits TE activity, are also found in U. hordei but have no homologs in U. maydis [131]. Therefore, this fungus lacks the usual mechanisms to maintain genome integrity. The absence of typical defense mechanisms might be partially compensated by its very active recombination system, which could facilitate genome stability [131]. Moreover, a correlation between loss of RNAi and length of centromeres has been observed among Ustilago and Cryptococcus species [96]. While the centromeric regions of the RNA-proficient species C. neoformans and C. deneoformans are composed of full-length transposable elements [96], C. deuterogattii, deficient in RNAi, lacks full-length TE sequences and, thus, features a shorter centromeric region [96,132]. Similarly, U. maydis harbors shorter centromeres and fewer transposons than U. hordei [96]. Those findings establish a relationship between RNAi and centromere structure evolution. On the other hand, the loss of RNAi in C. deuterogattii is likely associated with the evolution of an increased virulence since free movement of retrotransposable elements could provide rapid adaptative responses to the stressful environment host cells entail through a major genome plasticity [132].
Lastly, some budding yeast, including the human pathogen Candida albicans and Naumovozyma castellii (previously known as Saccharomyces castellii and close relative of S. cerevisiae), have lost some components of the RNAi machinery, but have replaced them with noncanonical RNAi proteins to produce dsRNA and short hairpin RNA (shRNA) as regulatory and defensive elements against transposable elements [133]. However, S. cerevisiae lacks the key RNAi components and, thus, the ability of TE inhibition by RNAi. To understand why RNAi has been negatively selected in these organisms, we should take into account their evolutionary context. For example, S. cerevisiae can be infected by mycoviruses known as “Killer viruses”, which are dsRNA viruses that can be cytoplasmically inherited. This endemic virus encodes for a toxin that kills the neighboring cells lacking the “Killer system” (not infected by the mycovirus), while conferring immunity to those making the toxin. Restoring the RNAi mechanism in S. cerevisiae, by introducing the Dicer and Ago proteins from the related species N. castellii, contributes to genome stability through TE silencing [133], but negatively affects the Killer system. The acquirement and retention of Killer virus, thus, favored the evolution of the RNAi-deficient species [134]. Therefore, the role and even the presence of RNAi differ across all fungal kingdom, as it has evolved adapting to the particularities and needs of each species. Fungi have exploited RNAi to conquer different niches, adapt to diverse stresses, acquire drug tolerance, and avoid or even modulate the host immune response during infection. According to these RNAi-mediated roles, it would not be wrong to say that RNAi has played a very significant role in the evolutionary success of fungi as a kingdom, as well as individual organisms.

Author Contributions

C.L. wrote the Control of Transposable Elements section; G.T. wrote Defense against virus and Loss of RNAi; J.A.P.-M. wrote the Pathogenesis section; J.T.C.-M. and J.A.P.-R. wrote the Regulation of Endogenous Genes section; M.O.-C. wrote Heterochromatin formation; E.N. wrote Adaptation to Stressful Conditions; S.C. wrote the introduction, made the figures and tables and gave coherence to the manuscript. All the authors revised and approved the final draft. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

We would like to thank Victoriano Garre and Vikas Yadav for their critical reading and valuable comments. We also thank Juan García for lending us his graphic design skills.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

RNAiRNA interference
dsRNAdouble-stranded RNA
siRNAssmall interfering RNAs
sRNAssmall RNAs
RISCRNA-induced silencing complex
RITSRNA-induced transcriptional silencing complex
ssRNAsingle-stranded RNA
RdRP/RdpRNA-dependent RNA polymerase
PTGSPost-transcriptional gene silencing
MSUDMeiotic silencing of unpaired DNA
masiRNAsMSUD-associated siRNAs
vsRNAsvirus-derived small interfering RNAs
TEtransposable elements
SCANRspliceosome-coupled and nuclear RNAi
TUTasesterminal-uridylyl transferases
miRNAsmicro RNAs
milRNAsmicroRNA-like RNAs
NCRIPnoncanonical RNAi pathway
ex-siRNAsexonic siRNAs
rdRNAsRdRP-dependent dicer-independent sRNAs
disiRNAsDicer-independent small interfering RNAs
DLDMdisiRNA loci DNA methylation
tRFstRNA fragments
circRNAscircular RNAs
CLRCcryptic loci regulator complex
H3K9memethylation of the histone 3 lysin 9
5-FOA5-fluoroorotic acid

References

  1. Chang, S.S.; Zhang, Z.; Liu, Y. RNA interference pathways in fungi: Mechanisms and functions. Annu. Rev. Microbiol. 2012, 66, 305–323. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ghildiyal, M.; Zamore, P.D. Small silencing RNAs: An expanding universe. Nat. Rev. Genet. 2009, 10, 94–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Cerutti, H.; Casas-Mollano, J.A. On the origin and functions of RNA-mediated silencing: From protists to man. Curr. Genet. 2006, 50, 91–99. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Matveyev, A.V.; Alves, J.M.P.; Serrano, M.G.; Lee, V.; Lara, A.M.; Barton, W.A.; Costa-Martins, A.G.; Beverley, S.M.; Camargo, E.P.; Teixeira, M.M.G.; et al. The Evolutionary Loss of RNAi Key Determinants in Kinetoplastids as a Multiple Sporadic Phenomenon. J. Mol. Evol. 2017, 84, 104–115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Ketting, R.F. The Many Faces of RNAi. Dev. Cell 2011, 20, 148–161. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Kataoka, K.; Mochizuki, K. Programmed DNA elimination in tetrahymena: A small RNA-mediated genome surveillance mechanism. Adv. Exp. Med. Biol. 2011, 722, 156–173. [Google Scholar] [CrossRef] [Green Version]
  7. Nicolás, F.E.; Torres-Martínez, S.; Ruiz-Vázquez, R.M. Loss and Retention of RNA Interference in Fungi and Parasites. PLoS Pathog. 2013, 9, e1003089. [Google Scholar] [CrossRef] [Green Version]
  8. Sarkies, P.; Selkirk, M.E.; Jones, J.T.; Blok, V.; Boothby, T.; Goldstein, B.; Hanelt, B.; Ardila-Garcia, A.; Fast, N.M.; Schiffer, P.M.; et al. Ancient and Novel Small RNA Pathways Compensate for the Loss of piRNAs in Multiple Independent Nematode Lineages. PLoS Biol. 2015, 13, e1002061. [Google Scholar] [CrossRef] [Green Version]
  9. Tomoyasu, Y.; Miller, S.C.; Tomita, S.; Schoppmeier, M.; Grossmann, D.; Bucher, G. Exploring systemic RNA interference in insects: A genome-wide survey for RNAi genes in Tribolium. Genome Biol. 2008, 9, R10. [Google Scholar] [CrossRef] [Green Version]
  10. Dowling, D.; Pauli, T.; Donath, A.; Meusemann, K.; Podsiadlowski, L.; Petersen, M.; Peters, R.S.; Mayer, C.; Liu, S.; Zhou, X.; et al. Phylogenetic origin and diversification of RNAi pathway genes in insects. Genome Biol. Evol. 2016, 8, 3784–3793. [Google Scholar] [CrossRef] [Green Version]
  11. Hedil, M.; Kormelink, R. Viral RNA silencing suppression: The enigma of bunyavirus NSs proteins. Viruses 2016, 8, 208. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Muhammad, T.; Zhang, F.; Zhang, Y.; Liang, Y. RNA Interference: A Natural Immune System of Plants to Counteract Biotic Stressors. Cells 2019, 8, 38. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Wei, W.; Ba, Z.; Gao, M.; Wu, Y.; Ma, Y.; Amiard, S.; White, C.I.; Danielsen, J.M.R.; Yang, Y.G.; Qi, Y. A role for small RNAs in DNA double-strand break repair. Cell 2012, 149, 101–112. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Svoboda, P. Renaissance of mammalian endogenous RNAi. FEBS Lett. 2014, 588, 2550–2556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Romano, N.; Macino, G. Quelling: Transient inactivation of gene expression in Neurospora crassa by transformation with homologous sequences. Mol. Microbiol. 1992, 6, 3343–3353. [Google Scholar] [CrossRef]
  16. Cogoni, C.; Macino, G. Gene silencing in Neurospora crassa requires a protein homologous to RNA-dependent RNA polymerase. Nature 1999, 399, 166–169. [Google Scholar] [CrossRef]
  17. Nakayashiki, H. RNA silencing in fungi: Mechanisms and applications. FEBS Lett. 2005, 579, 5950–5957. [Google Scholar] [CrossRef] [Green Version]
  18. Son, H.; Min, K.; Lee, J.; Raju, N.B.; Lee, Y.W. Meiotic silencing in the homothallic fungus Gibberella zeae. Fungal Biol. 2011, 12, 1290–1302. [Google Scholar] [CrossRef]
  19. Wang, Y.; Smith, K.M.; Taylor, J.W.; Freitag, M.; Stajich, J.E. Endogenous small RNA mediates meiotic silencing of a novel DNA transposon. G3 Genes Genomes Genet. 2015, 5, 1949–1960. [Google Scholar] [CrossRef] [Green Version]
  20. Decker, L.M.; Boone, E.C.; Xiao, H.; Shanker, B.S.; Boone, S.F.; Kingston, S.L.; Lee, S.A.; Hammond, T.M.; Shiu, P.K.T. Complex formation of RNA silencing proteins in the perinuclear region of Neurospora crassa. Genetics 2015, 199, 1017–1021. [Google Scholar] [CrossRef] [Green Version]
  21. Hirano, Y.; Asakawa, H.; Sakuno, T.; Haraguchi, T.; Hiraoka, Y. Nuclear Envelope Proteins Modulating the Heterochromatin Formation and Functions in Fission Yeast. Cells 2020, 9, 1908. [Google Scholar] [CrossRef] [PubMed]
  22. Billmyre, R.B.; Calo, S.; Feretzaki, M.; Wang, X.; Heitman, J. RNAi function, diversity, and loss in the fungal kingdom. Chromosom. Res. 2013, 21, 561–572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Torres-Martínez, S.; Ruiz-Vázquez, R.M. The RNAi Universe in Fungi: A Varied Landscape of Small RNAs and Biological Functions. Annu. Rev. Microbiol. 2017, 71, 371–391. [Google Scholar] [CrossRef] [PubMed]
  24. Sun, Q.; Choi, G.H.; Nuss, D.L. A single Argonaute gene is required for induction of RNA silencing antiviral defense and promotes viral RNA recombination. Proc. Natl. Acad. Sci. USA 2009, 106, 17927–17932. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Zhang, X.; Nuss, D.L. A host dicer is required for defective viral RNA production and recombinant virus vector RNA instability for a positive sense RNA virus. Proc. Natl. Acad. Sci. USA 2008, 105, 16749–16754. [Google Scholar] [CrossRef] [Green Version]
  26. Hammond, T.M.; Andrewski, M.D.; Roossinck, M.J.; Keller, N.P. Aspergillus mycoviruses are targets and suppressors of RNA silencing. Eukaryot. Cell 2008, 7, 350–357. [Google Scholar] [CrossRef] [Green Version]
  27. Nolan, T.; Braccini, L.; Azzalin, G.; De Toni, A.; Macino, G.; Cogoni, C. The post-transcriptional gene silencing machinery functions independently of DNA methylation to repress a LINE1-like retrotransposon in Neurospora crassa. Nucleic Acids Res. 2005, 33, 1564–1573. [Google Scholar] [CrossRef]
  28. Shiu, P.K.T.; Raju, N.B.; Zickler, D.; Metzenberg, R.L. Meiotic silencing by unpaired DNA. Cell 2001, 107, 905–916. [Google Scholar] [CrossRef] [Green Version]
  29. Yamanaka, S.; Mehta, S.; Reyes-Turcu, F.E.; Zhuang, F.; Fuchs, R.T.; Rong, Y.; Robb, G.B.; Grewal, S.I.S. RNAi triggered by specialized machinery silences developmental genes and retrotransposons. Nature 2013, 493, 557–560. [Google Scholar] [CrossRef] [Green Version]
  30. Nunes, C.C.; Gowda, M.; Sailsbery, J.; Xue, M.; Chen, F.; Brown, D.E.; Oh, Y.Y.; Mitchell, T.K.; Dean, R.A. Diverse and tissue-enriched small RNAs in the plant pathogenic fungus, Magnaporthe oryzae. BMC Genom. 2011, 12, 288. [Google Scholar] [CrossRef] [Green Version]
  31. Kadotani, N.; Nakayashiki, H.; Tosa, Y.; Mayama, S. RNA silencing in the phytopathogenic fungus Magnaporthe oryzae. Mol. Plant-Microbe Interact. 2003, 16, 769–776. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Murata, T.; Kadotani, N.; Yamaguchi, M.; Tosa, Y.; Mayama, S.; Nakayashiki, H. siRNA-dependent and -independent post-transcriptional cosuppression of the LTR-retrotransposon MAGGY in the phytopathogenic fungus Magnaporthe oryzae. Nucleic Acids Res. 2007, 35, 5987–5994. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Nakayashiki, H.; Ikeda, K.; Hashimoto, Y.; Tosa, Y.; Mayama, S. Methylation is not the main force repressing the retrotransposon MAGGY in Magnaporthe grisea. Nucleic Acids Res. 2001, 29, 1278–1284. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Navarro-Mendoza, M.I.; Pérez-Arques, C.; Panchal, S.; Nicolás, F.E.; Mondo, S.J.; Ganguly, P.; Pangilinan, J.; Grigoriev, I.V.; Heitman, J.; Sanyal, K.; et al. Early Diverging Fungus Mucor circinelloides Lacks Centromeric Histone CENP-A and Displays a Mosaic of Point and Regional Centromeres. Curr. Biol. 2019, 29, 3791–3802.e6. [Google Scholar] [CrossRef] [Green Version]
  35. Janbon, G.; Maeng, S.; Yang, D.H.; Ko, Y.J.; Jung, K.W.; Moyrand, F.; Floyd, A.; Heitman, J.; Bahn, Y.S. Characterizing the role of RNA silencing components in Cryptococcus neoformans. Fungal Genet. Biol. 2010, 47, 1070. [Google Scholar] [CrossRef] [Green Version]
  36. Wang, X.; Hsueh, Y.P.; Li, W.; Floyd, A.; Skalsky, R.; Heitman, J. Sex-induced silencing defends the genome of Cryptococcus neoformans via RNAi. Genes Dev. 2010, 24, 2566–2582. [Google Scholar] [CrossRef] [Green Version]
  37. Wang, X.; Darwiche, S.; Heitman, J. Sex-induced silencing operates during opposite-sex and unisexual reproduction in Cryptococcus neoformans. Genetics 2013, 193, 1163–1174. [Google Scholar] [CrossRef] [Green Version]
  38. Dumesic, P.A.; Natarajan, P.; Chen, C.; Drinnenberg, I.A.; Schiller, B.J.; Thompson, J.; Moresco, J.J.; Yates, J.R.; Bartel, D.P.; Madhani, H.D. Stalled spliceosomes are a signal for RNAi-mediated genome defense. Cell 2013, 152, 957–968. [Google Scholar] [CrossRef] [Green Version]
  39. Burke, J.E.; Longhurst, A.D.; Natarajan, P.; Rao, B.; Liu, J.; Sales-Lee, J.; Mortensen, Y.; Moresco, J.J.; Diedrich, J.K.; Yates, J.R.; et al. A non-dicer RNase III and four other novel factors required for RNAi-Mediated transposon suppression in the human pathogenic yeast cryptococcus neoformans. G3 Genes Genomes Genet. 2019, 9, 2235–2244. [Google Scholar] [CrossRef] [Green Version]
  40. Sperschneider, J.; Jones, A.W.; Nasim, J.; Xu, B.; Jacques, S.; Upadhyaya, N.M.; Mago, R.; Figueroa, M.; Singh, K.B.; Stone, E.A.; et al. The stem rust fungus Puccinia graminis f. sp. tritici induces centromeric small RNAs during late infection that direct genome-wide DNA methylation. bioRxiv 2020, 469338. [Google Scholar] [CrossRef] [Green Version]
  41. Nicolás, F.E.; Vila, A.; Moxon, S.; Cascales, M.D.; Torres-Martínez, S.; Ruiz-Vázquez, R.M.; Garre, V. The RNAi machinery controls distinct responses to environmental signals in the basal fungus Mucor circinelloides. BMC Genom. 2015, 16, 237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Calo, S.; Nicolás, F.E.; Lee, S.C.; Vila, A.; Cervantes, M.; Torres-Martinez, S.; Ruiz-Vazquez, R.M.; Cardenas, M.E.; Heitman, J. A non-canonical RNA degradation pathway suppresses RNAi-dependent epimutations in the human fungal pathogen Mucor circinelloides. PLoS Genet. 2017, 13, e1006686. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Trieu, T.A.; Calo, S.; Nicolás, F.E.; Vila, A.; Moxon, S.; Dalmay, T.; Torres-Martínez, S.; Garre, V.; Ruiz-Vázquez, R.M. A Non-canonical RNA Silencing Pathway Promotes mRNA Degradation in Basal Fungi. PLoS Genet. 2015, 11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Carreras-Villaseñor, N.; Esquivel-Naranjo, E.U.; Villalobos-Escobedo, J.M.; Abreu-Goodger, C.; Herrera-Estrella, A. The RNAi machinery regulates growth and development in the filamentous fungus Trichoderma atroviride. Mol. Microbiol. 2013, 89, 96–112. [Google Scholar] [CrossRef] [PubMed]
  45. Son, H.; Park, A.R.; Lim, J.Y.; Shin, C.; Lee, Y.W. Genome-wide exonic small interference RNA-mediated gene silencing regulates sexual reproduction in the homothallic fungus Fusarium graminearum. PLoS Genet. 2017, 13, e1006595. [Google Scholar] [CrossRef]
  46. Zeng, W.; Wang, J.; Wang, Y.; Lin, J.; Fu, Y.; Xie, J.; Jiang, D.; Chen, T.; Liu, H.; Cheng, J. Dicer-like proteins regulate sexual development via the biogenesis of perithecium-specific microRNAs in a plant pathogenic fungus Fusarium graminearum. Front. Microbiol. 2018, 9, 818. [Google Scholar] [CrossRef] [Green Version]
  47. Lau, A.Y.T.; Cheng, X.; Cheng, C.K.; Nong, W.; Cheung, M.K.; Chan, R.H.F.; Hui, J.H.L.; Kwan, H.S. Discovery of microRNA-like RNAs during early fruiting body development in the model mushroom Coprinopsis cinerea. PLoS ONE 2018, 13, e0198234. [Google Scholar] [CrossRef] [Green Version]
  48. Lee, H.C.; Li, L.; Gu, W.; Xue, Z.; Crosthwaite, S.K.; Pertsemlidis, A.; Lewis, Z.A.; Freitag, M.; Selker, E.U.; Mello, C.C.; et al. Diverse pathways generate microRNA-like RNAs and Dicer-independent small interfering RNAs in fungi. Mol. Cell 2010, 38, 803–814. [Google Scholar] [CrossRef] [Green Version]
  49. Dang, Y.; Cheng, J.; Sun, X.; Zhou, Z.; Liu, Y. Antisense transcription licenses nascent transcripts to mediate transcriptional gene silencing. Genes Dev. 2016, 30, 2417–2432. [Google Scholar] [CrossRef] [Green Version]
  50. Dang, Y.; Li, L.; Guo, W.; Xue, Z.; Liu, Y. Convergent Transcription Induces Dynamic DNA Methylation at disiRNA Loci. PLoS Genet. 2013, 9, e1003761. [Google Scholar] [CrossRef] [Green Version]
  51. Yuan, J.; Wang, Z.; Xing, J.; Yang, Q.; Chen, X.L. Genome-wide Identification and characterization of circular RNAs in the rice blast fungus Magnaporthe oryzae. Sci. Rep. 2018, 8, 6757. [Google Scholar] [CrossRef] [PubMed]
  52. Åsman, A.K.M.; Vetukuri, R.R.; Jahan, S.N.; Fogelqvist, J.; Corcoran, P.; Avrova, A.O.; Whisson, S.C.; Dixelius, C. Fragmentation of tRNA in Phytophthora infestans asexual life cycle stages and during host plant infection. BMC Microbiol. 2014, 14, 308. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Martienssen, R.; Moazed, D. RNAi and heterochromatin assembly. Cold Spring Harb. Perspect. Biol. 2015, 7, a019323. [Google Scholar] [CrossRef] [PubMed]
  54. Yu, R.; Wang, X.; Moazed, D. Epigenetic inheritance mediated by coupling of RNAi and histone H3K9 methylation. Nature 2018, 558, 615–619. [Google Scholar] [CrossRef]
  55. Calo, S.; Shertz-Wall, C.; Lee, S.C.; Bastidas, R.J.; Nicolás, F.E.; Granek, J.A.; Mieczkowski, P.; Torres-Martínez, S.; Ruiz-Vázquez, R.M.; Cardenas, M.E.; et al. Antifungal drug resistance evoked via RNAi-dependent epimutations. Nature 2014, 513. [Google Scholar] [CrossRef] [Green Version]
  56. Chang, Z.; Billmyre, R.B.; Lee, S.C.; Heitman, J. Broad antifungal resistance mediated by RNAi-dependent epimutation in the basal human fungal pathogen Mucor circinelloides. PLoS Genet. 2019, 15, e1007957. [Google Scholar] [CrossRef] [Green Version]
  57. Wang, Q.; An, B.; Hou, X.; Guo, Y.; Luo, H.; He, C. Dicer-like proteins regulate the growth, conidiation, and pathogenicity of Colletotrichum gloeosporioides from Hevea brasiliensis. Front. Microbiol. 2018, 8, 2621. [Google Scholar] [CrossRef] [Green Version]
  58. Raman, V.; Simon, S.A.; Demirci, F.; Nakano, M.; Meyers, B.C.; Donofrio, N.M. Small RNA functions are required for growth and development of magnaporthe oryzae. Mol. Plant-Microbe Interact. 2017, 30, 517–530. [Google Scholar] [CrossRef] [Green Version]
  59. Mochama, P.; Jadhav, P.; Neupane, A.; Marzano, S.Y.L. Mycoviruses as triggers and targets of RNA silencing in white mold fungus Sclerotinia sclerotiorum. Viruses 2018, 10, 214. [Google Scholar] [CrossRef] [Green Version]
  60. Feng, H.; Xu, M.; Liu, Y.; Gao, X.; Yin, Z.; Voegele, R.T.; Huang, L. The distinct roles of Argonaute protein 2 in the growth, stress responses and pathogenicity of the apple tree canker pathogen. For. Pathol. 2017, 47, e12354. [Google Scholar] [CrossRef]
  61. Jo, S.M.; Ayukawa, Y.; Yun, S.H.; Komatsu, K.; Arie, T. A putative RNA silencing component protein FoQde-2 is involved in virulence of the tomato wilt fungus Fusarium oxysporum f. sp. lycopersici. J. Gen. Plant Pathol. 2018, 84, 395–398. [Google Scholar] [CrossRef]
  62. Yin, C.; Zhu, H.; Jiang, Y.; Shan, Y.; Gong, L. Silencing Dicer-Like Genes Reduces Virulence and sRNA Generation in Penicillium italicum, the Cause of Citrus Blue Mold. Cells 2020, 9, 363. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  63. Gaffar, F.Y.; Imani, J.; Karlovsky, P.; Koch, A.; Kogel, K.-H. Different Components of the RNA Interference Machinery Are Required for Conidiation, Ascosporogenesis, Virulence, Deoxynivalenol Production, and Fungal Inhibition by Exogenous Double-Stranded RNA in the Head Blight Pathogen Fusarium graminearum. Front. Microbiol. 2019, 10, 1662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Pérez-Arques, C.; Navarro-Mendoza, M.I.; Murcia, L.; Navarro, E.; Garre, V.; Nicolás, F.E. A non-canonical RNAi pathway controls virulence and genome stability in Mucorales. PLoS Genet. 2020, 16, e1008611. [Google Scholar] [CrossRef]
  65. Qiao, Y.; Liu, L.; Xiong, Q.; Flores, C.; Wong, J.; Shi, J.; Wang, X.; Liu, X.; Xiang, Q.; Jiang, S.; et al. Oomycete pathogens encode RNA silencing suppressors. Nat. Genet. 2013, 45, 330–333. [Google Scholar] [CrossRef] [Green Version]
  66. Weiberg, A.; Wang, M.; Lin, F.M.; Zhao, H.; Zhang, Z.; Kaloshian, I.; Huang, H.D.; Jin, H. Fungal small RNAs suppress plant immunity by hijacking host RNA interference pathways. Science 2013, 342, 118–123. [Google Scholar] [CrossRef] [Green Version]
  67. Wang, B.; Sun, Y.; Song, N.; Zhao, M.; Liu, R.; Feng, H.; Wang, X.; Kang, Z. Puccinia striiformis f. sp. tritici microRNA-like RNA 1 (Pst-milR1), an important pathogenicity factor of Pst, impairs wheat resistance to Pst by suppressing the wheat pathogenesis-related 2 gene. New Phytol. 2017, 215, 338–350. [Google Scholar] [CrossRef] [Green Version]
  68. Kusch, S.; Frantzeskakis, L.; Thieron, H.; Panstruga, R. Small RNAs from cereal powdery mildew pathogens may target host plant genes. Fungal Biol. 2018, 122, 1050–1063. [Google Scholar] [CrossRef]
  69. Dunker, F.; Trutzenberg, A.; Rothenpieler, J.S.; Kuhn, S.; Pröls, R.; Schreiber, T.; Tissier, A.; Kemen, A.; Kemen, E.; Hückelhoven, R.; et al. Oomycete small RNAs bind to the plant RNA-induced silencing complex for virulence. Elife 2020, 9, e56096. [Google Scholar] [CrossRef]
  70. Segers, G.C.; Zhang, X.; Deng, F.; Sun, Q.; Nuss, D.L. Evidence that RNA silencing functions as an antiviral defense mechanism in fungi. Proc. Natl. Acad. Sci. USA 2007, 104, 12902–12906. [Google Scholar] [CrossRef] [Green Version]
  71. Choudhary, S.; Lee, H.-C.; Maiti, M.; He, Q.; Cheng, P.; Liu, Q.; Liu, Y. A Double-Stranded-RNA Response Program Important for RNA Interference Efficiency. Mol. Cell. Biol. 2007, 27, 3995–4005. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Segers, G.C.; Van Wezel, R.; Zhang, X.; Hong, Y.; Nuss, D.L. Hypovirus papain-like protease p29 suppresses RNA silencing in the natural fungal host and in a heterologous plant system. Eukaryot. Cell 2006, 5, 896–904. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Kotta-Loizou, I.; Coutts, R.H.A. Mycoviruses in Aspergilli: A comprehensive review. Front. Microbiol. 2017, 8, 1699. [Google Scholar] [CrossRef] [Green Version]
  74. Bourque, G.; Burns, K.H.; Gehring, M.; Gorbunova, V.; Seluanov, A.; Hammell, M.; Imbeault, M.; Izsvák, Z.; Levin, H.L.; Macfarlan, T.S.; et al. Ten things you should know about transposable elements. Genome Biol. 2018, 19, 199. [Google Scholar] [CrossRef] [PubMed]
  75. Klein, S.J.; O’Neill, R.J. Transposable elements: Genome innovation, chromosome diversity, and centromere conflict. Chromosom. Res. 2018, 26, 5–23. [Google Scholar] [CrossRef] [Green Version]
  76. Hirsch, C.D.; Springer, N.M. Transposable element influences on gene expression in plants. Biochim. Biophys. Acta Gene Regul. Mech. 2017, 1860, 157–165. [Google Scholar] [CrossRef] [Green Version]
  77. Diehl, A.G.; Ouyang, N.; Boyle, A.P. Transposable elements contribute to cell and species-specific chromatin looping and gene regulation in mammalian genomes. Nat. Commun. 2020, 11, 1796. [Google Scholar] [CrossRef] [Green Version]
  78. Friedli, M.; Trono, D. The Developmental Control of Transposable Elements and the Evolution of Higher Species. Annu. Rev. Cell Dev. Biol. 2015, 31, 429–451. [Google Scholar] [CrossRef]
  79. Castanera, R.; López-Varas, L.; Borgognone, A.; LaButti, K.; Lapidus, A.; Schmutz, J.; Grimwood, J.; Pérez, G.; Pisabarro, A.G.; Grigoriev, I.V.; et al. Transposable Elements versus the Fungal Genome: Impact on Whole-Genome Architecture and Transcriptional Profiles. PLoS Genet. 2016, 12, e1006108. [Google Scholar] [CrossRef] [Green Version]
  80. Bourc’his, D.; Bestor, T.H. Meiotic catastrophe and retrotransposon reactivation in male germ cells lacking Dnmt3L. Nature 2004, 431, 96–99. [Google Scholar] [CrossRef]
  81. Bulut-Karslioglu, A.; DeLaRosa-Velázquez, I.A.; Ramirez, F.; Barenboim, M.; Onishi-Seebacher, M.; Arand, J.; Galán, C.; Winter, G.E.; Engist, B.; Gerle, B.; et al. Suv39h-Dependent H3K9me3 Marks Intact Retrotransposons and Silences LINE Elements in Mouse Embryonic Stem Cells. Mol. Cell 2014, 55, 277–290. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  82. Quenneville, S.; Turelli, P.; Bojkowska, K.; Raclot, C.; Offner, S.; Kapopoulou, A.; Trono, D. The KRAB-ZFP/KAP1 System Contributes to the Early Embryonic Establishment of Site-Specific DNA Methylation Patterns Maintained during Development. Cell Rep. 2012, 2, 766–773. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  83. Bourc’his, D.; Voinnet, O. A small-RNA perspective on gametogenesis, fertilization, and early zygotic development. Science 2010, 330, 617–622. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Aravin, A.A.; Naumova, N.M.; Tulin, A.V.; Vagin, V.V.; Rozovsky, Y.M.; Gvozdev, V.A. Double-stranded RNA-mediated silencing of genomic tandem repeats and transposable elements in the D. melanogaster germline. Curr. Biol. 2001, 11, 1017–1027. [Google Scholar] [CrossRef] [Green Version]
  85. Sijen, T.; Plasterk, R.H.A. Transposon silencing in the Caenorhabditis elegans germ line by natural RNAi. Nature 2003, 426, 310–314. [Google Scholar] [CrossRef] [PubMed]
  86. Ketting, R.F.; Haverkamp, T.H.A.; Van Luenen, H.G.A.M.; Plasterk, R.H.A. mut-7 of C. elegans, required for transposon silencing and RNA interference, is a homolog of werner syndrome helicase and RNaseD. Cell 1999, 99, 133–141. [Google Scholar] [CrossRef] [Green Version]
  87. Tabara, H.; Sarkissian, M.; Kelly, W.G.; Fleenor, J.; Grishok, A.; Timmons, L.; Fire, A.; Mello, C.C. The rde-1 gene, RNA interference, and transposon silencing in C. elegans. Cell 1999, 99, 123–132. [Google Scholar] [CrossRef] [Green Version]
  88. Wu-Scharf, D.; Jeong, B.R.; Zhang, C.; Cerutti, H. Transgene and transposon silencing in Chlamydomonas reinhardtii by a DEAH-Box RNA helicase. Science 2000, 290, 1159–1162. [Google Scholar] [CrossRef] [Green Version]
  89. Soifer, H.S.; Zaragoza, A.; Peyvan, M.; Behlke, M.A.; Rossi, J.J. A potential role for RNA interference in controlling the activity of the human LINE-1 retrotransposon. Nucleic Acids Res. 2005, 33, 846–856. [Google Scholar] [CrossRef] [Green Version]
  90. Hammond, T.M.; Bok, J.W.; Andrewski, M.D.; Reyes-Domínguez, Y.; Scazzocchio, C.; Keller, N.P. RNA silencing gene truncation in the filamentous fungus Aspergillus nidulans. Eukaryot. Cell 2008, 7, 339–349. [Google Scholar] [CrossRef] [Green Version]
  91. Wagner, L.; De Hoog, S.; Alastruey-Izquierdo, A.; Voigt, K.; Kurzai, O.; Walther, G. A revised species concept for opportunistic Mucor species reveals species-specific antifungal susceptibility profiles. Antimicrob. Agents Chemother. 2019, 63, e00653-19. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Nicolas, F.E.; Moxon, S.; de Haro, J.P.; Calo, S.; Grigoriev, I.V.; Torres-MartÍnez, S.; Moulton, V.; Ruiz-Vázquez, R.M.; Dalmay, T. Endogenous short RNAs generated by Dicer 2 and RNA-dependent RNA polymerase 1 regulate mRNAs in the basal fungus Mucor circinelloides. Nucleic Acids Res. 2010, 38. [Google Scholar] [CrossRef] [PubMed]
  93. Loftus, B.J.; Fung, E.; Roncaglia, P.; Rowley, D.; Amedeo, P.; Bruno, D.; Vamathevan, J.; Miranda, M.; Anderson, I.J.; Fraser, J.A.; et al. The genome of the basidiomycetous yeast and human pathogen Cryptococcus neoformans. Science 2005, 307, 1321–1324. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Warkocki, Z.; Krawczyk, P.S.; Adamska, D.; Bijata, K.; Garcia-Perez, J.L.; Dziembowski, A. Uridylation by TUT4/7 Restricts Retrotransposition of Human LINE-1s. Cell 2018, 174, 1537–1548.e29. [Google Scholar] [CrossRef] [Green Version]
  95. Jiang, N.; Yang, Y.; Janbon, G.; Pan, J.; Zhu, X. Identification and Functional Demonstration of miRNAs in the Fungus Cryptococcus neoformans. PLoS ONE 2012, 7, e52734. [Google Scholar] [CrossRef] [Green Version]
  96. Yadav, V.; Sun, S.; Billmyre, R.B.; Thimmappa, B.C.; Shea, T.; Lintner, R.; Bakkeren, G.; Cuomo, C.A.; Heitman, J.; Sanyal, K. RNAi is a critical determinant of centromere evolution in closely related fungi. Proc. Natl. Acad. Sci. USA 2018, 115, 3108–3113. [Google Scholar] [CrossRef] [Green Version]
  97. Fahlgren, N.; Bollmann, S.R.; Kasschau, K.D.; Cuperus, J.T.; Press, C.M.; Sullivan, C.M.; Chapman, E.J.; Hoyer, J.S.; Gilbert, K.B.; Grünwald, N.J.; et al. Phytophthora Have Distinct Endogenous Small RNA Populations That Include Short Interfering and microRNAs. PLoS ONE 2013, 8, e77181. [Google Scholar] [CrossRef]
  98. Vetukuri, R.R.; Åsman, A.K.M.; Tellgren-Roth, C.; Jahan, S.N.; Reimegård, J.; Fogelqvist, J.; Savenkov, E.; Söderbom, F.; Avrova, A.O.; Whisson, S.C.; et al. Evidence for Small RNAs Homologous to Effector-Encoding Genes and Transposable Elements in the Oomycete Phytophthora infestans. PLoS ONE 2012, 7, e51399. [Google Scholar] [CrossRef] [Green Version]
  99. Grandaubert, J.; Lowe, R.G.T.; Soyer, J.L.; Schoch, C.L.; Van De Wouw, A.P.; Fudal, I.; Robbertse, B.; Lapalu, N.; Links, M.G.; Ollivier, B.; et al. Transposable element-assisted evolution and adaptation to host plant within the Leptosphaeria maculans-Leptosphaeria biglobosa species complex of fungal pathogens. BMC Genomics 2014, 15, 891. [Google Scholar] [CrossRef] [Green Version]
  100. Lax, C.; Pérez-arques, C.; Navarro-mendoza, M.I.; Cánovas-márquez, J.T.; Tahiri, G.; Pérez-ruiz, J.A.; Osorio-concepción, M.; Murcia-flores, L.; Navarro, E.; Garre, V.; et al. Genes, pathways, and mechanisms involved in the virulence of mucorales. Genes 2020, 11, 317. [Google Scholar] [CrossRef] [Green Version]
  101. Nicolás, F.E.; Murcia, L.; Navarro, E.; Navarro-Mendoza, M.I.; Pérez-Arques, C.; Garre, V. Mucorales species and macrophages. J. Fungi 2020, 6, 94. [Google Scholar] [CrossRef] [PubMed]
  102. Cervantes, M.; Vila, A.; Nicolas, F.E.; Moxon, S.; de Haro, J.P.; Dalmay, T.; Torres-Martinez, S.; Ruiz-Vazquez, R.M. A single argonaute gene participates in exogenous and endogenous RNAi and controls cellular functions in the basal fungus Mucor circinelloides. PLoS ONE 2013, 8, e69283. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  103. De Haro, J.P.; Calo, S.; Cervantes, M.; Nicolás, F.E.; Torres-Martinez, S.; Ruiz-Vázquez, R.M. A single dicer gene is required for efficient gene silencing associated with two classes of small antisense RNAs in mucor circinelloides. Eukaryot. Cell 2009, 8, 1486–1497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  104. Nicolas, F.E.; de Haro, J.P.; Torres-Martinez, S.; Ruiz-Vazquez, R.M. Mutants defective in a Mucor circinelloides dicer-like gene are not compromised in siRNA silencing but display developmental defects. Fungal Genet. Biol. 2007, 44, 504–516. [Google Scholar] [CrossRef]
  105. Ambros, V. The functions of animal microRNAs. Nature 2004, 431, 350–355. [Google Scholar] [CrossRef]
  106. Yang, Q.; Li, L.; Xue, Z.; Ye, Q.; Zhang, L.; Li, S.; Liu, Y. Transcription of the Major Neurospora crassa microRNA-Like Small RNAs Relies on RNA Polymerase III. PLoS Genet. 2013, 9, e1003227. [Google Scholar] [CrossRef] [Green Version]
  107. Chen, Y.; Gao, Q.; Huang, M.; Liu, Y.; Liu, Z.; Liu, X.; Ma, Z. Characterization of RNA silencing components in the plant pathogenic fungus Fusarium graminearum. Sci. Rep. 2015, 5, 12500. [Google Scholar] [CrossRef] [Green Version]
  108. Lee, Y.S.; Shibata, Y.; Malhotra, A.; Dutta, A. A novel class of small RNAs: tRNA-derived RNA fragments (tRFs). Genes Dev. 2009, 23, 2639–2649. [Google Scholar] [CrossRef] [Green Version]
  109. Sobala, A.; Hutvagner, G. Small RNAs derived from the 5′ end of tRNAs can inhibit protein translation in human cells. RNA Biol. 2013, 10, 553–563. [Google Scholar] [CrossRef] [Green Version]
  110. Maute, R.L.; Schneider, C.; Sumazin, P.; Holmes, A.; Califano, A.; Basso, K.; Dalla-Favera, R. TRNA-derived microRNA modulates proliferation and the DNA damage response and is down-regulated in B cell lymphoma. Proc. Natl. Acad. Sci. USA 2013, 110, 1404–1409. [Google Scholar] [CrossRef] [Green Version]
  111. Deng, J.; Ptashkin, R.N.; Chen, Y.; Cheng, Z.; Liu, G.; Phan, T.; Deng, X.; Zhou, J.; Lee, I.; Lee, Y.S.; et al. Respiratory Syncytial Virus Utilizes a tRNA Fragment to Suppress Antiviral Responses Through a Novel Targeting Mechanism. Mol. Ther. 2015, 23, 1622–1629. [Google Scholar] [CrossRef] [Green Version]
  112. Liu, J.; Ali, M.; Zhou, Q. Establishment and evolution of heterochromatin. Ann. N. Y. Acad. Sci. 2020, 9. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  113. Allshire, R.C.; Madhani, H.D. Ten principles of heterochromatin formation and function. Nat. Rev. Mol. Cell Biol. 2018, 19, 229–244. [Google Scholar] [CrossRef]
  114. Taglini, F.; Chapman, E.; van Nues, R.; Theron, E.; Bayne, E.H. Mkt1 is required for RNAi-mediated silencing and establishment of heterochromatin in fission yeast. Nucleic Acids Res. 2020, 48, 1239–1253. [Google Scholar] [CrossRef] [Green Version]
  115. Volpe, T.A.; Kidner, C.; Hall, I.M.; Teng, G.; Grewal, S.I.; Martienssen, R.A. Regulation of heterochromatic silencing and histone H3 lysine-9 methylation by RNAi. Science 2002, 297, 1833–1837. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  116. Verdel, A.; Jia, S.; Gerber, S.; Sugiyama, T.; Gygi, S.; Grewal, S.I.; Moazed, D. RNAi-mediated targeting of heterochromatin by the RITS complex. Science 2004, 303, 672–676. [Google Scholar] [CrossRef] [Green Version]
  117. Bhattacharjee, S.; Roche, B.; Martienssen, R.A. RNA-induced initiation of transcriptional silencing (RITS) complex structure and function. RNA Biol. 2019, 16, 1133–1146. [Google Scholar] [CrossRef] [PubMed]
  118. Jih, G.; Iglesias, N.; Currie, M.A.; Bhanu, N.V.; Paulo, J.A.; Gygi, S.P.; Garcia, B.A.; Moazed, D. Unique roles for histone H3K9me states in RNAi and heritable silencing of transcription. Nature 2017, 547, 463–467. [Google Scholar] [CrossRef] [Green Version]
  119. Duempelmann, L.; Skribbe, M.; Bühler, M. Small RNAs in the Transgenerational Inheritance of Epigenetic Information. Trends Genet. 2020, 36, 203–214. [Google Scholar] [CrossRef] [Green Version]
  120. Thon, G.; Hansen, K.R.; Altes, S.P.; Sidhu, D.; Singh, G.; Verhein-Hansen, J.; Bonaduce, M.J.; Klar, A.J.S. The Clr7 and Clr8 directionality factors and the Pcu4 cullin mediate heterochromatin formation in the fission yeast Schizosaccharomyces pombe. Genetics 2005, 171, 1583–1595. [Google Scholar] [CrossRef] [Green Version]
  121. Hall, I.M.; Shankaranarayana, G.D.; Noma, K.; Ayoub, N.; Cohen, A.; Grewal, S.I. Establishment and maintenance of a heterochromatin domain. Science 2002, 297, 2232–2237. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  122. Goto, D.B.; Nakayama, J.I. RNA and epigenetic silencing: Insight from fission yeast. Dev. Growth Differ. 2012, 54, 129–141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  123. Sadaie, M.; Kawaguchi, R.; Ohtani, Y.; Arisaka, F.; Tanaka, K.; Shirahige, K.; Nakayama, J. Balance between Distinct HP1 Family Proteins Controls Heterochromatin Assembly in Fission Yeast. Mol. Cell. Biol. 2008, 28, 6973–6988. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  124. Duempelmann, L.; Mohn, F.; Shimada, Y.; Oberti, D.; Andriollo, A.; Lochs, S.; Bühler, M. Inheritance of a Phenotypically Neutral Epimutation Evokes Gene Silencing in Later Generations. Mol. Cell 2019, 74, 534–541. [Google Scholar] [CrossRef] [PubMed]
  125. Lee, S.C.; Li, A.; Calo, S.; Inoue, M.; Tonthat, N.K.; Bain, J.M.; Louw, J.; Shinohara, M.L.; Erwig, L.P.; Schumacher, M.A.; et al. Calcineurin orchestrates dimorphic transitions, antifungal drug responses and host-pathogen interactions of the pathogenic mucoralean fungus Mucor circinelloides. Mol. Microbiol. 2015, 97, 844–865. [Google Scholar] [CrossRef] [Green Version]
  126. Chang, Z.; Yadav, V.; Lee, S.C.; Heitman, J. Epigenetic mechanisms of drug resistance in fungi. Fungal Genet. Biol. 2019, 132, 103253. [Google Scholar] [CrossRef]
  127. Chang, Z.; Heitman, J. Drug-resistant epimutants exhibit organ-specific stability and induction during murine infections caused by the human fungal pathogen Mucor circinelloides. MBio 2019, 10, e02579-19. [Google Scholar] [CrossRef] [Green Version]
  128. Neupane, A.; Feng, C.; Mochama, P.K.; Saleem, H.; Lee Marzano, S.Y. Roles of Argonautes and Dicers on Sclerotinia sclerotiorum Antiviral RNA Silencing. Front. Plant Sci. 2019, 10, 976. [Google Scholar] [CrossRef] [Green Version]
  129. Nicolás, F.E.; Garre, V. RNA Interference in Fungi: Retention and Loss. In The Fungal Kingdom; Heitman, J., Howlett, B.J., Crous, P.W., Stukenbrock, E.H., James, T.Y., Gow, N.A.R., Eds.; American Society for Microbiology: Washington, DC, USA, 2017; pp. 657–672. [Google Scholar]
  130. Choi, J.; Kim, K.T.; Jeon, J.; Wu, J.; Song, H.; Asiegbu, F.O.; Lee, Y.H. FunRNA: A fungi-centered genomics platform for genes encoding key components of RNAi. BMC Genom. 2014, 15, S14. [Google Scholar] [CrossRef] [Green Version]
  131. Laurie, J.D.; Ali, S.; Linning, R.; Mannhaupt, G.; Wong, P.; Güldener, U.; Münsterkötter, M.; Moore, R.; Kahmann, R.; Bakkeren, G.; et al. Genome comparison of barley and maize smut fungi reveals targeted loss of RNA silencing components and species-specific presence of transposable elements. Plant Cell 2012, 24, 1733–1745. [Google Scholar] [CrossRef] [Green Version]
  132. Feretzaki, M.; Billmyre, R.B.; Clancey, S.A.; Wang, X.; Heitman, J. Gene Network Polymorphism Illuminates Loss and Retention of Novel RNAi Silencing Components in the Cryptococcus Pathogenic Species Complex. PLoS Genet. 2016, 12, e1005868. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  133. Drinnenberg, I.A.; Weinberg, D.E.; Xie, K.T.; Mower, J.P.; Wolfe, K.H.; Fink, G.R.; Bartel, D.P. RNAi in budding yeast. Science 2009, 326, 544–550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  134. Drinnenberg, I.A.; Fink, G.R.; Bartel, D.P. Compatibility with killer explains the rise of RNAi-deficient fungi. Science 2011, 333, 1592. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Table 1. RNA interference (RNAi) mechanisms in eukaryotes.
Table 1. RNA interference (RNAi) mechanisms in eukaryotes.
Eukaryotic GroupEvolutionary AdvantageLoss of RNAiEffector ProteinssiRNAsReferences
FungiDefense against viruses
Control of TE
Regulation of endogenous genes
Heterochromatin formation
Adaptation to stressful conditions
Pathogenesis
Saccharomyces cerevisiae
Ustilago maydis
Cryptococcus deuterogatti
Ago
Dicer
Rdp/RdRP
siRNA
masiRNA
milRNA
ex-siRNA
endo-siRNA
rdRNA
circRNA
tRF
disiRNA
Table 2 and Table 3
ProtistsDefense against viruses
Control of TE
Heterochromatin formation
DNA elimination
Trypanosoma cruzi
Leishmania major
Leishmania donovani
Leishmania tarentolae
Plasmodium falciparum
Piwi
Dicer
RdRP
siRNA
scnRNA
[3,4,5,6,7]
NematodesDefense against viruses
Control of TE
Heterochromatin formation
Regulation of endogenous genes
Ago/Piwi
Dicer
Rrf/Ego
miRNA
piRNA
endo-siRNA
exo-siRNA
[3,5,8]
InsectsDefense against viruses
Control of TE
Heterochromatin formation
Regulation of endogenous genes
Ago/Piwi
Dicer/Loquacious/
Drosha/Pasha/R2D2
siRNA
miRNA
piRNA
endo-siRNA
[3,5,8,9,10]
PlantsDefense against pathogens
Control of TE
Heterochromatin formation
Regulation of endogenous genes
Repair DNA double-strand breaks (DSB)
Ago
Dcl
Rdr
siRNA
miRNA
hp-siRNA
hc-siRNA or hetsiRNA
tasiRNA
t-siRNA
lsiRNA
nat-siRNA
phasiRNA
easiRNA
diRNA
[3,8,11,12,13]
MammalsDefense against viruses
Control of TE
Heterochromatin formation
Regulation of endogenous genes
Repair DNA double-strand breaks
Ago/Piwi
Dicer/Drosha
siRNA
miRNA
piRNA
endo-siRNA
diRNA
[3,5,8,13,14]
TE: transposable elements; Ago: Argonaute; Rdp/RdRP/Rdr: RNA-dependent RNA polymerase; Dcl: Dicer-like protein; siRNA: small interfering RNA; masiRNA: MSUD-associated siRNAs; MSUD: meiotic silencing of unpaired DNA; milRNA: microRNA-like RNAs; ex-siRNA: exonic siRNAs; endo-siRNA: endogenous siRNA; rdRNA: rdrp-dependent dicer-independent sRNAs; circRNA: circular RNAs; tRF: transfer RNA fragments; disiRNA: Dicer-independent siRNAs; scnRNA: small scan RNA; miRNA: micro RNA; piRNA: Piwi-interacting RNA; exo-siRNA: exogenously-derived siRNA; hp-siRNA: hairpin-derived siRNA; hc-siRNA or hetsiRNA: heterochromatic siRNA; tasiRNA: trans-acting siRNA; t-siRNA: transgene-derived siRNA; lsiRNA: long siRNA; nat-siRNA: natural antisense siRNA; phasiRNA: phased siRNA; easiRNA: epigenetically-activated siRNA; diRNA: DSB-induced siRNA.
Table 2. Classification of RNA interference (RNAi) mechanisms found in fungi.
Table 2. Classification of RNA interference (RNAi) mechanisms found in fungi.
FunctionRegulation ofMain RNAi Proteins siRNAs RNAi MechanismFungi SpeciesReferences
Defense against virusesCryphonectria hypovirus 1Dcl2, Ago2vsRNAs Cryphonectria parasitica[24,25]
Aspergillus virus 1816DclB, RsdAvsRNAs Aspergillus nidulans[26]
Control of Transposable ElementLINE1-like retrotransposonsQDE-2, Dcl1, Dcl2 QuellingNeurospora crassa[27]
DNA transposon
Sly1-1
Dcl1, SAD-1-5, SMS-2masiRNAsMSUDNeurospora crassa[19,28]
Tf2 retrotransposonAgo1, Rdp1, Dcr1, Clr4 RNA-induced transcriptional silencing (RITS)Schizosaccharomyces pombe[29]
LTR-retrotransposon MAGGYMoDcl2, MoAgo1, MoAgo3LTR-siRNAs Magnaporthe oryzae[30,31,32,33]
Grem-LINE1Ago2, Dcl1, Dcl2 Canonical Mucor lusitanicus[34]
RetrotransposonsRdp1, Dcr1, Dcr2, Ago1, Qip1, Gwc1, Srr1endo-siRNAsSex-induced silencing (SIS)
Spliceosome-Coupled And Nuclear RNAi complex (SCANR)
Cryptococcus neoformans[35,36,37,38]
DNA transposon
HAR1
Rde1-5endo-siRNAs Cryptococcus neoformans[39]
Young TEAgo1, Dicer3, RdRP1/3Late wave sRNAs Puccinia graminis f. sp. triticci[40]
Regulation of endogenous genes RdRP1, RdRP2, Dcl1, Dcl2, Ago1, R3B2ex-siRNAsCanonicalMucor lusitanicus[23,41]
Canonical RNAi, phagocytosis responseRdRP1, RdRP3, R3B2, RnhArdRNAsNon-canonical (NCRIP)Mucor lusitanicus[42,43]
Growth and developmentDcr2, Rdr3ex-siRNAs Trichoderma atroviride[44]
Ascospore formationFgDcl1, FgAgo2ex-siRNAs Fusarium graminearum[45]
Sexual developmentFgDcl1/2, FgAgo1milRNAs Fusarium graminearum[46]
Development and fruiting body formation milRNAs Coprinopsis cinerea[47]
Dicer, QIP, MRPL3, QDE-2milRNAs Neurospora crassa[48]
ERI-1, QDE-2disiRNAsdisiRNA loci DNA methylation (DLDM)Neurospora crassa[49,50]
miRNA circRNAs Magnaporthe oryzae[51]
tRFs Magnaporthe oryzae[30]
Plant infectionPiAgo1tRFs Phytophthora infestans[52]
Heterochromatin formationCentromeric regionsDcr1, Ago1, Rdp1, Tas3, Chp1Centromeric siRNAs Schizosaccharomyces pombe[53,54]
Adaptation to stressful conditionsFKBP12
PyrG, PyrF
Dcr1-2, Ago1, RdRP2, QIP, RnhA EpimutationMucor circinelloides
Mucor lusitanicus
[55,56]
PathogenesisAbility to infect leavesDcr1, Dcr2 Colletotrichum gloeosporioides[57]
Ability to infect leavesAgo3, Rdp1 Magnaportheoryzae[58]
VirulenceDcr1-2, Ago2 Sclerotinia sclerotiorum[59]
Resistance to H2O2 and virulenceAgo2 Valsa mali[60]
VirulenceAgo2 Fusarium oxysporum f. sp. lycopersici[61]
InfectionDcr2 Penicillium italicum[62]
InfectionFgDcl1, FgAgo2 Fusarium graminearum[63]
Resistance to H2O2 and virulenceRdRP1, R3B2rdRNAsNCRIPMucor lusitanicus[64]
Host miRNA and siRNA PSR1 and PSR2Cross-kingdom RNAiPhytophthora sojae[65]
Host immunity genes Bc-sRNAsCross-kingdom RNAiBotrytis cinerea[66]
Host defenses Pst-milR1Cross-kingdom RNAiPuccinia striiformis f. sp. tritici[67]
Host metabolism genes Bgh-sRNAs
Bgt-sRNAs
Cross-kingdom RNAiBlumeria graminis f. sp. hordei
Blumeria graminis f. sp. tritici
[68]
Host genes Cross-kingdom RNAiHyaloperonospora arabidopsidis[69]
Grem-LINE1: genomic retrotransposable element of Mucoromycota; TE: Transposable Elements; miRNA: micro RNA; siRNA: small interfering RNA; Dcr: Dicer; Dcl: Dicer-like protein; Ago: Argonaute; Rdp/RdRP/Rdr: RNA-dependent RNA polymerase; vsRNAs: virus-derived siRNAs; masiRNA: MSUD-associated siRNAs; MSUD: meiotic silencing of unpaired DNA; endo-siRNA: endogenous siRNA; ex-siRNA: exonic siRNAs; rdRNA: rdrp-dependent dicer-independent sRNAs; milRNA: microRNA-like RNAs; disiRNA: Dicer-independent siRNAs; circRNA: circular RNAs; tRF: transfer RNA fragments; sRNAs: small RNAs; NCRIP: noncanonical RNAi pathway.
Table 3. Loss of the RNAi mechanism among fungi species.
Table 3. Loss of the RNAi mechanism among fungi species.
SpeciesPutative AdvantageLost ProteinsSubstitutesReferences
Ustilago maydisRetention of Killer virusAgo1, RdRP1-3, Dcr1Partially by recombination system[131]
Crytpococcus deuterogattiIncreased virulence by TE activityAgo1, Ago2, Dcr1, Rdp1 [96]
Candida albicans DicerNoncanonical CaDcr1[133]
Naumovozyma castellii DicerNoncanonical NcaDcr1[133]
Saccharomyces cerevisiaeRetention of Killer virusAgo1, Dcr1 [134]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lax, C.; Tahiri, G.; Patiño-Medina, J.A.; Cánovas-Márquez, J.T.; Pérez-Ruiz, J.A.; Osorio-Concepción, M.; Navarro, E.; Calo, S. The Evolutionary Significance of RNAi in the Fungal Kingdom. Int. J. Mol. Sci. 2020, 21, 9348. https://doi.org/10.3390/ijms21249348

AMA Style

Lax C, Tahiri G, Patiño-Medina JA, Cánovas-Márquez JT, Pérez-Ruiz JA, Osorio-Concepción M, Navarro E, Calo S. The Evolutionary Significance of RNAi in the Fungal Kingdom. International Journal of Molecular Sciences. 2020; 21(24):9348. https://doi.org/10.3390/ijms21249348

Chicago/Turabian Style

Lax, Carlos, Ghizlane Tahiri, José Alberto Patiño-Medina, José T. Cánovas-Márquez, José A. Pérez-Ruiz, Macario Osorio-Concepción, Eusebio Navarro, and Silvia Calo. 2020. "The Evolutionary Significance of RNAi in the Fungal Kingdom" International Journal of Molecular Sciences 21, no. 24: 9348. https://doi.org/10.3390/ijms21249348

APA Style

Lax, C., Tahiri, G., Patiño-Medina, J. A., Cánovas-Márquez, J. T., Pérez-Ruiz, J. A., Osorio-Concepción, M., Navarro, E., & Calo, S. (2020). The Evolutionary Significance of RNAi in the Fungal Kingdom. International Journal of Molecular Sciences, 21(24), 9348. https://doi.org/10.3390/ijms21249348

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop