The RNAi machinery controls distinct responses to environmental signals in the basal fungus Mucor circinelloides
© Nicolás et al.; licensee BioMed Central. 2015
Received: 7 November 2014
Accepted: 6 March 2015
Published: 25 March 2015
RNA interference (RNAi) is a conserved mechanism of genome defence that can also have a role in the regulation of endogenous functions through endogenous small RNAs (esRNAs). In fungi, knowledge of the functions regulated by esRNAs has been hampered by lack of clear phenotypes in most mutants affected in the RNAi machinery. Mutants of Mucor circinelloides affected in RNAi genes show defects in physiological and developmental processes, thus making Mucor an outstanding fungal model for studying endogenous functions regulated by RNAi. Some classes of Mucor esRNAs map to exons (ex-siRNAs) and regulate expression of the genes from which they derive. To have a broad picture of genes regulated by the silencing machinery during vegetative growth, we have sequenced and compared the mRNA profiles of mutants in the main RNAi genes by using RNA-seq. In addition, we have achieved a more complete phenotypic characterization of silencing mutants.
Deletion of any main RNAi gene provoked a deep impact in mRNA accumulation at exponential and stationary growth. Genes showing increased mRNA levels, as expected for direct ex-siRNAs targets, but also genes with decreased expression were detected, suggesting that, most probably, the initial ex-siRNA targets regulate the expression of other genes, which can be up- or down-regulated. Expression of 50% of the genes was dependent on more than one RNAi gene in agreement with the existence of several classes of ex-siRNAs produced by different combinations of RNAi proteins. These combinations of proteins have also been involved in the regulation of different cellular processes. Besides genes regulated by the canonical RNAi pathway, this analysis identified processes, such as growth at low pH and sexual interaction that are regulated by a dicer-independent non-canonical RNAi pathway.
This work shows that the RNAi pathways play a relevant role in the regulation of a significant number of endogenous genes in M. circinelloides during exponential and stationary growth phases and opens up an important avenue for in-depth study of genes involved in the regulation of physiological and developmental processes in this fungal model.
RNA silencing or RNA interference (RNAi) is a conserved defence mechanism against invasive nucleic acids, such as viruses, transposons or transgenes, in a wide spectrum of eukaryotic organisms, including fungi [1,2]. Double stranded RNA (dsRNA) molecules, derived from those exogenous sequences, are processed by the RNase III Dicer to produce short interfering RNAs (siRNAs), which are bound to an Argonaute protein within the RNA-induced silencing complex (RISC). These siRNAs serve as a guide to identify complementary target RNA molecules for silencing or destruction [3,4]. In many organisms, siRNA-induced silencing requires RNA-dependent RNA polymerases (RdRPs) to generate dsRNA from single-stranded RNA (ssRNA) or to amplify siRNA signals . In metazoans, RNAi also has a role in the regulation of endogenous functions through several classes of endogenous small RNA (esRNAs) molecules, which are generated from genome-encoded precursors. These endogenous pathways play many fundamental roles, including regulation of mRNA accumulation and translation, chromatin silencing, programmed DNA rearrangements and genome surveillance . Different classes of regulatory esRNAs have been also described in fungi, although information on their functional roles is very scarce .
Mucor circinelloides, a basal fungus belonging to the order Mucorales, has become a good model organism for the study of different molecular processes in the fungal kingdom, including light responses  and gene silencing , mainly due to the availability of a large number of molecular tools and to its evolutionary distance from other fungal model organisms, such as Neurospora crassa. More recently, M. circinelloides is attracting special attention as a causal agent of mucormycosis, an emerging fungal infection, not very common but often lethal, caused by various species of the order Mucorales . Although typically affects immunocompromised patients , various circumstances such as the use of certain antifungal drugs  and some natural disasters  have increased the number of cases in risk populations and have boosted the interest in further study about pathogenesis of this fungus [11-13]. Of particular interest is the recent discovery of a new epigenetic mechanism for developing transient resistance to an antifungal drug via an RNAi-mediated pathway, this epigenetic mechanism being particularly enhanced in pathogenic strains of M. circinelloides .
Exogenously-induced RNA silencing in M. circinelloides is associated with the accumulation of two size classes of siRNAs, 21 and 25 nt long, which are differentially accumulated during the vegetative growth . Only one of the two dicer genes that have been identified in M. circinelloides, dcl-2, is required for transgene-induced silencing and accumulation of the two size classes of siRNAs, both when silencing is induced by sense and inverted repeat transgenes [16,17]. Unlike N. crassa, RNA silencing in M. circinelloides is associated with an amplification step that generates secondary siRNAs corresponding to target sequences by the RNA-dependent RNA polymerase activity of the rdrp-2 gene product . A functionally distinct rdrp gene, rdrp-1, is essential for initiation of silencing by sense transgenes by producing antisense RNA transcripts derived from the transgene, but it is not involved in the amplification of the silencing signal . Finally, only one of the three ago genes identified in M. circinelloides, ago-1, is required for transgene-induced gene silencing, whatever the nature of the silencing trigger is, as shown by the negative silencing phenotype of ago-1 − mutants expressing sense- or inverted-repeat transgenes. Since neither primary nor secondary siRNAs are detected in those mutants, it has been suggested that Ago-1 is required for production/stability of siRNAs .
As in metazoans, the M. circinelloides RNAi pathway also has a role in the regulation of endogenous genes through several classes of esRNA molecules, which are generated from genome-encoded precursors . Deep sequencing of small RNAs endogenously accumulated in the wild type strain, dicer − and rdrp − mutants identified a number of esRNAs that map to exons and regulate the expression of many protein coding genes. These esRNAs, named exonic-siRNAs (ex-siRNAs), can be classified in different classes based on the silencing proteins required for their biogenesis. In addition to its role in silencing exogenous sequences, ago-1 is also required for the production of all of these ex-siRNAs . A large group of them (Class II), including 222 exons, is dcl-2-dependent and also requires the rdrp-1 gene product, whereas a small group of only nine exons (Class I), which is also dcl-2-dependent, does not require the rdrp-1 gene product but most of them requires RdRP-2 . These two dcl-2-dependent classes of ex-siRNAs are down-regulated in the ago-1 − mutant and are specifically bound to Ago-1, suggesting that Ago-1 is involved in the biogenesis/stability of these regulatory ex-siRNAs. Binding to Ago-1 indicates that they are functional siRNAs produced by a canonical RNAi pathway to suppress the expression of the corresponding target genes. In fact, validation experiments demonstrated that lack of detection of specific ex-siRNAs of these classes in the ago-1 − mutant was associated with an increase of mRNA accumulation of the corresponding protein coding genes. Thus, these ex-siRNAs regulate the expression of the protein coding genes from which they are derived . Classes III and IV of ex-siRNAs do not specifically bind Ago-1, although they are down-regulated in the ago-1 − mutant, suggesting that Ago-1 participates in the biogenesis of these ex-siRNAs . Class III ex-siRNAs (88 exons) can be produced either by Dcl-1 or Dcl-2 and requires RdRP-1 and RdRP-2 for their biogenesis, whereas class IV is a tiny group of ex-siRNAs that requires the dcl-1 gene product . The biogenesis requirements and structural characteristics of classes III and IV have suggested that these ex-siRNAs are not produced by a canonical RNAi pathway, a frequent situation in filamentous fungi .
Knowledge of the endogenous functions regulated by esRNAs in fungi has been hampered by the lack of clear phenotypes in the majority of mutants affected in the RNAi machinery. However, M. circinelloides silencing mutants display diverse phenotypes . The dcl-1 − mutants are affected in their development, showing a significant reduction in the vegetative growth rate and alterations in hyphal morphology, which are observed in the leading edge of the fungal colony . Also, a light-regulated process as spore production is significantly reduced in dcl-2 − , dcl-1 − /dcl-2 − , and ago-1 − mutants relative to the wild-type strain [17,19]. Besides a defect in sporulation, ago-1 − and other silencing mutants, such as dcl-2 − and rdrp-2 − , share a common phenotype related to the autolytic response to nutritional stress . The autolysis of filamentous fungi is not a simple cell necrosis phenomenon but is an active and well regulated process where many enzymatic activities are involved . The accelerated lysis phenotype shared by these silencing mutants suggests that this process could be controlled by gene/s regulated by ex-siRNAs. Vegetative growth, hyphal morphology, sporulation and autolysis are developmental responses to endogenous and environmental signals in which many genes are involved, thus revealing the importance of the RNAi machinery in the endogenous regulation of complex developmental processes. The involvement of the RNA silencing pathways in the response to environmental signals is supported by the ability of M. circinelloides to adapt to the environment through RNAi-dependent epimutations , pointing out the relevance of the RNAi mechanism in controlling phenotypic plasticity.
We have previously demonstrated that M. circinelloides ex-siRNAs regulate the expression of the protein coding genes from which they are derived [19,20]. To have a comprehensive picture of genes regulated by the silencing machinery during the exponential and stationary growth phases, we have addressed the characterization of the mRNA profiles shown by mutants in the RNAi genes. In addition, we have achieved a more complete phenotypic characterization of silencing mutants, mainly on aspects related to the response to different environmental stimuli. The results obtained allow getting a deeper insight into the role of the RNAi machinery in the regulation of endogenous processes in M. circinelloides.
Differential response to environmental signals of mutants in RNAi genes
Transcript profiles of the RNAi machinery mutants
Genes differentially expressed in silencing mutants
Growth time (h)
Average Log 2 fold change a
Average Log 2 fold change a
The role of ago-1 in genes regulated by dcl-1/dcl-2 was also observed in the up-regulated genes, since 32 out of 35 genes that increased their expression in the dcl-1Δ/dcl-2Δ mutant at 24 h showed the same behavior in the ago-1Δ strain. However, in this case, most of these genes (59.4%) were also up-regulated in the rdrp-2Δ strain (Figure 4, Additional file 1: Table S1) and only a small proportion of them (1%) also required rdrp-1, with a single gene showing increased expression exclusively in dcl-1Δ/dcl-2Δ, ago-1Δ and rdrp-1Δ (Figure 4, Additional file 1: Table S1). These results indicate that, besides Dicer and Ago-1, the RdRP-2 protein plays a prominent role in the up-regulation of genes during the exponential growth and suggest the participation of class I ex-siRNAs in this regulation, since those are the proteins required for the biogenesis of this ex-siRNA class.
The expression patterns observed at 48 h were very similar to those observed at 24 h, although the number of regulated genes was smaller (Figure 4, Table 1, Additional file 1: Table S1). However, there were some significant differences mainly concentrated in the repressed genes. Thus, the proportion of down-regulated genes in dcl-1Δ/dcl-2Δ that were also down-regulated in ago-1Δ was smaller (43.8%) than at 24 h and most of them (66.7%) were regulated by rdrp-2. This is in contrast to what happens at 24 h, where most of the repressed genes were regulated by rdrp-1. During stationary growth, only 38.1% of the genes repressed in dcl-1Δ/dcl-2Δ and ago-1Δ were also repressed in rdrp-1Δ, being all of them also altered in rdrp-2Δ (Figure 4, Additional file 1: Table S1). These results suggest that rdrp-2 plays a more relevant role than rdrp-1 in down regulation of genes during stationary phase (Figure 4, Additional file 1: Table S1). Moreover, changes in expression were not as large as those observed at exponential growth, with a small number of genes having Log2 fold change values lower than −3 (Table 1, Additional file 1: Table S1). On the other hand and similarly to what was observed at 24 h, most of genes (61.1%) induced in dcl-1Δ/dcl-2Δ and ago-1Δ mutants at 48 h were also induced in rdrp-2Δ, even at higher levels than in exponential phase (Table 1) and only a small proportion of them (11.1%) were also up-regulated in rdrp-1Δ (Figure 4, Additional file 1: Table S1). No genes showing increased expression exclusively in dcl-1Δ/dcl-2Δ, ago-1Δ and rdrp-1Δ were found (Figure 4, Additional file 1: Table S1). The major expression pattern among the up-regulated genes (that is, dcl-1/2, ago-1 and rdrp-2-dependent) is expected from a direct regulation by ex-siRNAs of class I, as it is the case for two of the identified genes (ID115001, ID82197) (Additional file 1: Table S1), which correspond to class I ex-siRNA-producing loci . Finally, an important proportion of genes (307 out of 537 at 24 h and 152 out of 225 at 48 h) were differentially expressed in one or more mutants but not in dcl-1Δ/dcl-2Δ, suggesting that they are regulated by a dicer-independent RNAi pathway (; our unpublished results).
Validation of RNA-seq data
Putative functions of genes regulated by the RNAi genes
Analysis of the presumed functions of genes showing the highest differential expression in the RNAi mutants relative to wild-type identified genes that may be involved in growth, stress responses, and autolysis in M. circinelloides (Additional file 1: Table S1). Focusing on particular genes, we identified several down-regulated in RNAi mutants during exponential growth that code for proteins presumably involved in cell wall biogenesis and modification, including two putative chitin synthases (ID145794, ID104542) and four mannosyltransferases (ID14092, ID112156, ID 86104, ID 150168) (Additional file 1: Table S1), which are cell wall modifying proteins that are crucial for cell wall integrity and viability [24,25]. Also 10 genes involved in cell division were identified among the down-regulated ones, including 4 putative septins (ID156101, ID35281, ID104636, ID152534) (Additional file 1: Table S1). Down-regulation of those genes in RNAi mutants could be responsible for the phenotypic defects on growth and vegetative development shown by those mutants. Furthermore, some of the genes that were up-regulated in RNAi mutants during the exponential growth code for proteins containing domains associated with regulation of gene expression or signal transduction (Additional file 1: Table S1), which may explain the high number of deregulated genes in silencing mutants. Finally, several genes differently expressed at 48 h code for proteins probably related to the activity of transposons (ID113076, which contains a retrovirus zinc finger-like and retrotransposon gag protein domains), sexual interaction (ID107961, which codes for a short-chain dehydrogenase), sporulation, (ID156744, which is annotated as a carbohydrate esterase family 4 protein that deacetylates substrates such as chitin to chitosan, a component of the fungal cell wall/spore cortex), cellular aging (ID114253, a hydroxy-acyl-CoA dehydrogenase probably involved in cellular aging in rats ) and stress response (142978, a small heat-shock protein of the Hsp26/Hsp42 family) (Additional file 1: Table S1). Confirmation of the putative role of these and other genes on the different phenotypes shown by silencing mutants would require specific gene expression analysis under the condition of interest, in which the different phenotypes were clearly observed (see Discussion).
M. circinelloides represents an outstanding model for studying the function of the RNAi pathway in the control of physiological processes and development, since mutants in genes of this pathway show altered phenotypes, contrarily to what usually happens in other fungal models . Phenotypes associated with mutations in genes involved in sRNA biogenesis have been described in few filamentous fungi, and in all cases the phenotypic changes have been moderate. Particularly, mutants in dicer genes of Magnaporthe oryzae and T. atroviride show a reduction in growth rate compared with the wild type [23,28]. In addition, the T. atroviride mutants are also affected in the production of vegetative spores . The diverse phenotypes observed in the M. circinelloides RNAi mutants are probably a consequence of the altered expression of genes involved in different processes, since the ex-siRNAs generated by the RNAi pathway regulate the expression of a large number of protein coding genes [19,20]. Besides those from which they derive, these ex-siRNAs could also regulate in trans other genes with sequences completely or partially complementary to a particular ex-siRNA molecule. Given the vast number of ex-siRNAs with different sequences accumulated in the wild-type strain and the lack of knowledge about the required sequence complementarity, identification of all putative target genes of the ex-siRNAs by in silico analysis is an unaffordable task. One approach to identify genes controlled by the RNAi pathway that could be involved in fungal physiology and development is to compare the mRNA profiles of mutants in the RNAi pathway with that of the wild type strain. This work has used RNA-seq to characterize the mRNA profiles, at both exponential and stationary phase, of dcl-1Δ/dcl-2Δ, ago-1Δ, rdrp-1Δ and rdrp-2Δ mutants, which are affected in the key components of the RNAi pathway in M. circinelloides [17-19]. These conditions make possible to connect ex-siRNA classes produced during vegetative growth by specific combination of silencing proteins with gene clusters that show similar expression patterns in specific groups of silencing mutants, allowing the identification of ex-siRNA targets. Moreover, as production of asexual spores occurs during both exponential and stationary phase, the differential expression of specific genes in mutants affected in this process could support the role of those genes in the control of asexual sporulation.
Deletion of any of the silencing genes provoked a deep impact in mRNA accumulation at both growth phases, although the number of genes with altered expression was higher at exponential phase. In total, nearly 700 genes significantly changed their expression at least in one RNAi mutant and, interestingly, not only by increasing their mRNAs levels, as expected for a direct regulation by ex-siRNAs, but also by decreasing their expression. This suggests that the initial targets of ex-siRNAs can regulate the expression of other genes, which can be up- or down-regulated. In fact, only a small proportion (8.75%) of the genes up-regulated in the mutants corresponded to direct targets of ex-siRNAs, i.e. the ex-siRNAs-producing loci. This seems to be a common situation in fungi, since only a low proportion of up-regulated mRNAs in T. atroviride dicer mutants were correlated with down-regulation of the corresponding esRNAs . However, the possibility that RNAi components could regulate gene expression by an unknown mechanism independent of ex-siRNAs cannot be excluded.
Expression of around 50% of the genes was altered in more than one mutant, indicating that several silencing proteins participate in the regulation of many genes. This is an expected result because ex-siRNAs are produced by different combination of silencing proteins, and therefore the mRNA accumulation patterns should replicate those of the ex-siRNA ones . Dicer enzymes are involved in the biogenesis of the four classes of ex-siRNAs , and the mRNA accumulation patterns also revealed a main role for Dicer enzymes. However, the expression of a large number of genes was Dicer independent, supporting previous observations about the existence of a dcl-1/dcl-2-independent mechanism of ex-siRNA biogenesis , which is currently under study. On the other hand, almost all (86%) genes regulated by dcl-1/dcl-2 are also regulated by ago-1, suggesting that the basic components of the canonical RNAi machinery are involved in this regulation (Figure 4).
A large proportion of the genes regulated by dcl-1/dcl-2 and ago-1 were also regulated by one or both rdrp genes at both exponential and stationary phase. The function of rdrp genes is particularly interesting, since their contribution to the regulation changes with the growth phase and the type of regulation. Most of the genes that are repressed in exponential phase in the dcl-1Δ/dcl-2Δ and ago-1Δ mutants were also repressed in the rdrp-1Δ mutant and may represent indirect targets of class II ex-siRNAs (dcl-2, ago-1, rdrp-1-dependent) (Figures 3 and 4). The class II ex-siRNAs derive mostly from genes coding proteins of unknown function which, in many cases, contain domains involved in signal transduction or in processing genetic information , and therefore, may act on a wide range of secondary response genes. A small proportion of the genes regulated by dcl-1/dcl-2, ago-1 and rdrp-1 were also regulated by rdrp-2 and would represent secondary targets of class III (dcl-1/dcl-2, ago-1, rdrp-1, rdrp-2-dependent) or IV (dcl-1, ago-1, rdrp-1, rdrp-2-dependent) ex-siRNAs [19,20]. Class IV is a tiny group of ex-siRNAs that derive from only five exons and it is expected to have slight role in regulation, whereas class III ex-siRNAs derive from a large group of genes that are expressed at a high level in the wild type strain , which anticipates a significant role in controlling gene expression. The prevalent role of rdrp-1 in down-regulation during exponential phase was not restricted to genes regulated by Dicer proteins but was also observed in Dicer-independent regulated genes (Figure 4).
The impact of rdrp-1 in gene regulation seems to decrease with the growth phase, because most of the genes that are repressed during stationary phase in the dcl-1Δ/dcl-2Δ and ago-1Δ mutants were repressed in the rdrp-2Δ mutant but not in rdrp-1Δ (Figure 4). Only a few genes were down-regulated in the rdrp-1Δ mutant and all of them were also repressed in rdrp-2Δ (Figure 4). Therefore, most of genes down-regulated in dcl-1Δ/dcl-2Δ and ago-1Δ mutants at stationary phase could be indirect targets of class I ex-siRNAs, which is defined by their dependency on dcl-2, ago-1 and, in most cases, rdrp-2. Besides that, some of those genes could be targets of classes III or IV, as it was confirmed by northern hybridization for three of them (ID138264, ID156744, ID113332) (Figure 5).
Contrary to what happens with the down-regulated genes, a similar expression pattern was observed for up-regulated genes at both exponential and stationary phase. Most of the genes with increased expression in the dcl-1Δ/dcl-2Δ and ago-1Δ mutants were also up-regulated in rdrp-2Δ both at exponential and stationary phase (Figure 4). Gene rdrp-1 plays a minor role in the regulation of those genes, since only few of them increased their expression in both rdrp-1Δ and rdrp-2Δ mutants. These patterns suggest that they are direct targets of class I ex-siRNAs and, to a lesser extent, classes III or IV. The regulation pattern corresponding to these ex-siRNA classes was confirmed by northern-blot experiments (ID114253, ID82197, ID142978) (Figure 5).
Previous results and those of this work demonstrated that RNAi genes are involved in the control of several physiological and developmental processes, which include vegetative growth , production of asexual spores, accelerated autolysis in response to nutritional stress [17,19], growth at different pH conditions (Figure 1C) and sexual interaction (Figure 2). These processes are regulated by specific set of genes. Thus, asexual spore production is controlled mainly by dcl-2, ago-1 and rdrp-2, with a possible role for rdrp-1 gene in some particular conditions, such as poor nutrient environments. The same set of genes are also involved in accelerated autolysis , suggesting that both processes could be controlled by genes regulated directly or indirectly by ex-siRNAs, particularly class I, although regulation of those processes by RNAi components independently of ex-siRNAs cannot be excluded. In fungi, both cellular processes are connected with nutrient sensing of the cells, particularly carbon source levels in the medium [22,23]. Therefore, class I ex-siRNAs may be involved in fine tuning of genes involved in nutrient sensing pathways. Although only 5 different genes were identified as producers of this class of ex-siRNAs , the number of genes identified in our RNA-seq analysis that fit the class I ex-siRNA-dependent pattern was much higher, suggesting an amplification of the regulation. Thus, 26 genes showed mRNA patterns at 24 h, 48 h or at both growth times (ID114253 and ID77287) that correspond to regulation by class I ex-siRNAs (Additional file 6: Table S4). Interestingly, the mRNA levels of two of the 5 class I ex-siRNAs producer loci (ID82197 and ID115001) were increased at 48 h, one of them (ID82197) being validated by northern-blot hybridization (Figure 5). The identification of genes with a class I ex-siRNA expression pattern gives a strong support to the RNA-seq data and suggests that some of the identified genes may be involved in asexual spore production. This group of genes included a putative methyltransferase (ID90667) whose down-regulation in mutants with defective asexual sporulation was confirmed by Northern blot analysis (Figure 5). Methyltransferases have been shown to control several aspects of fungal development , including asexual sporulation in several fungal species [30-34]. In addition to this validated gene, other genes with similar expression pattern code for protein with putative functions or domains (Additional file 6: Table S4) that suggest that they could regulate sporulation. Thus, two genes (ID77287 and ID105388) code for proteins containing domains usually present in regulatory proteins, suggesting that they could be transcription factors controlling the expression of downstream genes. In addition, one gene (ID109934) codes for a predicted Rho GTPase-activating protein. Rho GTPases have been implicated in cellular integrity and spore wall maturation in Schizosacharomyces pombe [35,36], and the M. circinelloides Rho protein Cdc42 was found associated with chitosomes, which are the major reservoir of the chitin synthases required for the synthesis of the fungal cell wall . Moreover, G proteins signaling pathways have been involved in sporulation in several fungi, such as Aspergillus nidulans and N. crassa . Thus, modulation of the expression of this protein in the dcl-2Δ, ago-1Δ and rdrp-2Δ mutants could be responsible of the defects in sporulation shown by these mutants.
Other genes with the class I ex-siRNAs expression pattern (ID114253 and ID184709) are connected with lipid metabolism, which may play an important role in glucose starvation . In particular, protein ID184709 is a putative lipase suggested to supply substrates for gluconeogenesis and energy production, a process necessary to support the morphological changes required for sporulation or to facilitate the energy storage in spores . Taken together, the putative functions of several genes regulated by dcl-2, ago-1 and rdrp-2 genes, presumably through the action of class I ex-siRNAs, may explain the defects in asexual spore production shown by the corresponding null mutants. Detailed and individual functional analyses of each of those genes by deletion or overexpression would allow identification of the specific genes responsible for this phenotype.
The two processes regulated by RNAi identified in this work, growth at low pH and sexual interaction, are only controlled by rdrp-1 and rdrp-2, suggesting that a different ex-siRNA class is involved. These ex-siRNAs would be produced by a Dicer-independent mechanism in which Ago-1 would not participate in their biogenesis or function, although other M. circinelloides Ago proteins (Ago-2 and/or Ago-3) could be involved. Alternatively, it can be also possible that RdRP-1 and RdRP-2 control these biological processes independently of the regulation by ex-siRNAs. The RNA-seq analysis revealed few genes regulated by both rdrp-1 and rdrp-2 during vegetative growth (Additional file 7: Table S5), although some of them could have a role in signal transduction (ID80421), gene regulation (ID164671), or post-translational modification (ID79209, an F-box domain protein). Further transcriptomic experiments performed under low pH conditions or during sexual interaction would probably shed light on rdrp regulation of these processes. Particularly interesting is the putative role of rdrp genes in the regulation of sexual interaction. However, efforts to analyze the transcriptomes of the wild type and mutant strains during sexual interaction have been hampered by the long time required for interaction to occur, which makes difficult to isolate the high quality RNA samples required for sequencing. In fact, RNA degradation starts several days before sexual interaction takes place. A similar problem has been found to identify differentially expressed genes with a possible role in the accelerated autolytic response shown by several silencing mutants, since autolysis induced by nutritional stress takes place after a long period of incubation . Works are in progress to solve those problems and identify genes differentially expressed under the conditions of interest. Alternatively, detection of ex-siRNAs generated by an RdRP-dependent Dicer-independent RNAi pathway could enable the identification of target genes involved in the response to low pH and sexual interaction in M. circinelloides.
Results obtained in this work reveal a role of canonical and non-canonical RNAi pathways in the regulation of a significant number of endogenous genes in M. circinelloides both at exponential and stationary growth phase. Cellular processes that respond to nutrient levels in the environment, such as autolysis and asexual sporulation, are regulated by the dicer-dependent RNAi pathway, whereas processes regulated by other signals, such as mating and low pH, are controlled by a non-canonical dicer-independent pathway. Several genes whose expression is modulated by each pathway code for proteins that could be involved in the regulation of specific physiological and developmental processes in response to environmental signals in this fungal model. Further analysis will identify the specific genes that are causative of the different phenotypes shown by the RNAi mutants.
Strains and growth conditions
The strains used in this study are listed in Additional file 8: Table S6. The leucine auxotroph R7B , derived from M. circinelloides f. lusitanicus CBS 277.49 , was used as the wild-type strain. All the mutant strains (Additional file 8: Table S6) derived from MU402, a uracil auxotroph derivative of R7B . Cultures were grown in complete medium YPG , minimal YNB medium  or semi-complex MMC medium . Media were supplemented with L-leucine (20 μg/ml) or uridine (200 μg/ml), when required. The pH was adjusted to 4.5 for mycelial growth. For RNA extraction, 2.5 · 105 spores were grown at 26°C in solid YPG during 24 or 48 h under light conditions. Illumination conditions were as previously described .
Quantification of sexual mating and zygospores formation
Sexual interactions were induced in complete yeast and dextrose agar (YPD) media, in which spores of the mutant strains ((−) mating type) were co-inoculated in the middle of the agar plates along with the wild type strain (NRRL3631, (+) mating type), placed approximately 2 cm apart. Plates were incubated at room temperature, under dark conditions during 20 days. Under these conditions, a dark line appears in the center of the plate as a consequence of the formation of zygospores in the contact zone. This line was sliced in portions of 1 cm2 and fixed in 10% formaldehyde during 10 hours. After fixation, samples were frozen and sliced using a cryotome to produce sections of 30 μm. Zygospores from twelve sections were counted by optical microscopy (bright field 10X) for each interaction.
RNA extraction, library preparation and sequencing
Total RNA was extracted using RNeasy®Plant Mini Kit (Qiagen, Hilden, Germany) following the manufacturer instructions. In order to maximize target coverage, equal amounts of total RNA from three biological replicates of each strain were pooled for RNA-seq library construction. Both library preparation and sequencing was performed at Baseclear (Leiden, The Netherlands). The cDNA library was sequenced using Illumina HiSEQ2500 sequencer with a one lane, 50 cycles, single-read sequencing strategy. The total number of reads for each sample was between 12 and 21 million, more than 95% of them mapping onto the M. circinelloides genome (Additional file 9: Table S7).
Validation of RNA-seq data by northern blot
Total RNA was isolated using Trizol reagent following the instruction of the supplier (Gibco-BRL). Standard recombinant DNA manipulations were performed as described in Sambrook and Russell  and Ausubel et al. . Northern blots were hybridized to radioactively labelled probes and washed as described previously . Probes of around 500 bp were generated by PCR amplification using specific primers for each gene (Additional file 3: Table S3). The PCR fragments were labelled with [α-32P]dCTP using Ready-to-Go DNA labelling beads (Amersham Pharmacia Biotech), following the instructions of the supplier. Membranes were exposed to Kodak Phosphor Screen SD230 for quantification. After exposure to the membrane, the screen was scanned on a Molecular Imager FX reader (BioRad) and signals were quantified using ImageJ software (National Institute of Health). The 25S rRNA gene was used as an endogenous control to normalize for differences in the amount of RNA between samples .
Gene based expression counts were calculated for annotated transcripts from the Joint Genome Institute (JGI) M. circinelloides genome using RSEM  and rounded to the nearest whole number. Expression values were then imported into R and differential expression analysis between wild-type and mutant samples was performed using DESeq . To avoid infinite values, a value of 1 was added to the normalized count value of each gene with zero value before log2 transformation. Hierarchical clustering was performed using R function hclust with Manhattan distance.
Availability of supporting data
The raw reads of M. circinelloides mRNA are deposited in the NCBI GenBank as Sequence Read Archive (SRA) under the following accession numbers: SRX731621 (wild type 24 h), SRX731623 (wild type 48 h), SRX731630 (dcl-1Δ/dcl-2Δ 24 h), SRX731637 (dcl-1Δ/dcl-2Δ 48 h), SRX731624 (ago-1Δ 24 h), SRX731629 (ago-1Δ 48 h), SRX731639 (rdrp-1Δ 24 h), SRX731640 (rdrp-1Δ 48 h), SRX731642 (rdrp-2Δ 24 h) and SRX731641 (rdrp-2Δ 48 h).
We thank to the Joint Genome Institute for early access to the M. circinelloides genome sequence. The work conducted by the U.S. Department of Energy Joint Genome Institute is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231. We also thank J. A. Madrid for technical assistance. This work was supported by the Spanish Ministerio de Economía y Competitividad (BFU2012-32246, co-financed by FEDER).
- Buchon N, Vaury C. RNAi: a defensive RNA-silencing against viruses and transposable elements. Heredity. 2006;96:195–202.View ArticlePubMedGoogle Scholar
- Nicolás FE, Ruiz-Vázquez RM. Functional diversity of RNAi-associated sRNAs in fungi. Int J Mol Sci. 2013;14:15348–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Carthew RW, Sontheimer EJ. Origins and mechanisms of miRNAs and siRNAs. Cell. 2009;136:642–55.View ArticlePubMed CentralPubMedGoogle Scholar
- Ghildiyal M, Zamore PD. Small silencing RNAs: an expanding universe. Nat Rev Genet. 2009;10:94–108.View ArticlePubMed CentralPubMedGoogle Scholar
- Garre V, Nicolás FE, Torres-Martínez S, Ruiz-Vázquez RM. The RNAi Machinery in Mucorales: The Emerging Role of Endogenous Small RNAs. In: Sesma A, von der Haar T, editors. Fungal RNA Biology. Switzerland: Springer International Publishing; 2014. p. 291–313.View ArticleGoogle Scholar
- Corrochano LM, Garre V. Photobiology in the Zygomycota: multiple photoreceptor genes for complex responses to light. Fungal Genet Biol. 2010;47:893–9.View ArticlePubMedGoogle Scholar
- Roden MM, Zaoutis TE, Buchanan WL, Knudsen TA, Sarkisova TA, Schaufele RL, et al. Epidemiology and outcome of zygomycosis: a review of 929 reported cases. Clin Infect Dis. 2005;41:634–53.View ArticlePubMedGoogle Scholar
- Chayakulkeeree M, Ghannoum M, Perfect J. Zygomycosis: the reemerging fungal infection. Eur J Clin Microbiol Infect Dis. 2006;25:215–29.View ArticlePubMedGoogle Scholar
- Bitar D, Van Cauteren D, Lanternier F, Dannaoui E, Che D, Dromer F, et al. Increasing Incidence of Zygomycosis (Mucormycosis), France, 1997–2006. Emerg Infect Dis. 2009;15:1395–401.View ArticlePubMed CentralPubMedGoogle Scholar
- Fanfair R, Benedict K, Bos J. Necrotizing cutaneous mucormycosis after a tornado in Joplin, Missouri, in 2011. N Engl J Med. 2012;367:2214–25.View ArticleGoogle Scholar
- Li CH, Cervantes M, Springer DJ, Boekhout T, Ruiz-Vazquez RM, Torres-Martinez SR, et al. Sporangiospore size dimorphism is linked to virulence of Mucor circinelloides. PLoS Pathog. 2011;7:e1002086.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee SC, Li A, Calo S, Heitman J. Calcineurin plays key roles in the dimorphic transition and virulence of the human pathogenic zygomycete Mucor circinelloides. PLoS Pathog. 2013;9:e1003625.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee SC, Billmyre RB, Li A, Carson S, Sykes SM, Huh EY, et al. Analysis of a food-borne fungal pathogen outbreak: virulence and genome of a Mucor circinelloides isolate from yogurt. MBio. 2014;5:e01390–14.PubMed CentralPubMedGoogle Scholar
- Calo S, Shertz-Wall C, Lee AC, Bastidas RJ, Nicolás FE, Granek JA, et al. Antifungal drug resistance evoked via RNAi-dependent epimutations. Nature. 2014;513:555–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Nicolás FE, Torres-Martínez S, Ruiz-Vázquez RM. Two classes of small antisense RNAs in fungal RNA silencing triggered by non-integrative transgenes. EMBO J. 2003;22:3983–91.View ArticlePubMed CentralPubMedGoogle Scholar
- Nicolás FE, de Haro JP, Torres-Martínez S, Ruiz-Vázquez RM. 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–16.View ArticlePubMedGoogle Scholar
- de Haro JP, Calo S, Cervantes M, Nicolás FE, Torres-Martínez S, Ruiz-Vázquez RM. 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–97.View ArticlePubMed CentralPubMedGoogle Scholar
- Calo S, Nicolás FE, Vila A, Torres-Martínez S, Ruiz-Vázquez RM. Two distinct RNA-dependent RNA polymerases are required for initiation and amplification of RNA silencing in the basal fungus Mucor circinelloides. Mol Microbiol. 2012;83:379–94.View ArticlePubMedGoogle Scholar
- Cervantes M, Vila A, Nicolás FE, Moxon S, de Haro JP, Dalmay T, et al. 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.View ArticlePubMed CentralPubMedGoogle Scholar
- Nicolas FE, Moxon S, de Haro JP, Calo S, Grigoriev IV, Torres-Martínez S, et al. 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:5535–41.View ArticlePubMed CentralPubMedGoogle Scholar
- Lee HC, Li L, Gu W, Xue Z, Crosthwaite SK, Pertsemlidis A, et al. Diverse pathways generate microRNA-like RNAs and Dicer-independent small interfering RNAs in fungi. Mol Cell. 2019;38:803–14.View ArticleGoogle Scholar
- Emri T, Molnár Z, Szilágyi M, Pócsi I. Regulation of autolysis in Aspergillus nidulans. Appl Biochem Biotechnol. 2008;151:211–20.View ArticlePubMedGoogle Scholar
- Carreras-Villaseñor N, Esquivel-Naranjo EU, Villalobos-Escobedo JM, 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.View ArticlePubMedGoogle Scholar
- Willer T, Brand M, Sipiczki M, Strahl S. Protein O-mannosylation is crucial for cell wall integrity, septation and viability in fission yeast. Mol Microbiol. 2005;57:156–70.View ArticlePubMedGoogle Scholar
- Munro CA, Bates S, Buurman ET, Hughes HB, MacCallum DM, Bertram G, et al. Mnt1p and Mnt2p of Candida albicans are partially redundant alpha-1,2-mannosyltransferases that participate in O-linked mannosylation and are required for adhesion and virulence. J Biol Chem. 2005;280:1051–60.View ArticlePubMed CentralPubMedGoogle Scholar
- Ivell R, Balvers M, Anand RJ, Paust HJ, McKinnell C, Sharpe R. Differentiation-dependent expression of 17β-hydroxysteroid dehydrogenase, type 10, in the rodent testis: effect of aging in Leydig cells. Endocrinology. 2003;144:3130–7.View ArticlePubMedGoogle Scholar
- Chang SS, Zhang Z, Liu Y. RNA interference pathways in fungi: mechanisms and functions. Annu Rev Microbiol. 2012;66:305–23.View ArticlePubMedGoogle Scholar
- Kadotani N, Nakayashiki H, Tosa Y, Mayama S. One of the two Dicer-like proteins in the filamentous fungi Magnaporthe oryzae genome is responsible for hairpin RNA-triggered RNA silencing and related siRNA accumulation. J Biol Chem. 2004;279:44467–74.View ArticlePubMedGoogle Scholar
- Sarikaya-Bayram Ö, Palmer JM, Keller N, Braus GH, Bayram Ö. One Juliet and four Romeos: VeA and its methyltransferases. Front Microbiol. 2015;6:1.View ArticlePubMed CentralPubMedGoogle Scholar
- Karimi Aghcheh R, Druzhinina IS, Kubicek CP. The putative protein methyltransferase LAE1 of Trichoderma atroviride is a key regulator of asexual development and mycoparasitism. PLoS One. 2013;8:e67144.View ArticlePubMedGoogle Scholar
- Bok JW, Balajee SA, Marr KA, Andes D, Nielsen KF, Frisvad JC, et al. LaeA, a regulator of morphogenetic fungal virulence factors. Eukaryot Cell. 2005;4:1574–82.View ArticlePubMed CentralPubMedGoogle Scholar
- Chang PK, Scharfenstein LL, Ehrlich KC, Wie Q, Bhatnagar D, Ingber BF. Effects of laeA deletion on Aspergillus flavus conidial development and hydrophobicity may contribute to loss of aflatoxin production. Fungal Biol. 2012;116:298–307.View ArticlePubMedGoogle Scholar
- Hoff B, Kamerewerd J, Sig IC, Mitterbauer R, Zadra I, Kürnsteiner H, et al. Two components of a velvet-like complex control hyphal morphogenesis, conidiophore development, and penicillin biosynthesis in Penicillium chrysogenum. Eukaryot Cell. 2010;9:1236–50.View ArticlePubMed CentralPubMedGoogle Scholar
- Karimi-Aghcheh R, Bok JW, Phatale PA, Smith KM, Baker SE, Lichius A, et al. Functional analyses of Trichoderma reesei LAE1 reveal conserved and contrasting roles of this regulator. G3. 2013;3:369–78.View ArticlePubMed CentralPubMedGoogle Scholar
- Pérez P, Cansado J. Cell integrity signaling and response to stress in fission yeast. Curr Protein Pept Sci. 2010;11:680–92.View ArticlePubMedGoogle Scholar
- Rincón SA, Santos B, Pérez P. Fission yeast Rho5p GTPase is a functional paralogue of Rho1p that plays a role in survival of spores and stationary-phase cells. Eukaryot Cell. 2006;5:435–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Moreno-Jimenez R, Garcia-Soto J, Martinez-Cadena G. Small GTP-binding proteins are associated with chitosomes and vesicles carrying glucose oxidase from Mucor circinelloides. Microbiology. 2008;154:842–51.View ArticlePubMedGoogle Scholar
- Park H, Yu J. Genetic control of asexual sporulation in filamentous fungi. Curr Opin Microbiol. 2012;15:669–77.View ArticlePubMedGoogle Scholar
- Roncero MIG. Enrichment method for the isolation of auxotrophic mutants of Mucor circinelloides using the polyene antibiotic N-glycosyl-polyfungin. Carlsberg Res Commun. 1984;49:685–90.View ArticleGoogle Scholar
- Schipper MAA. On Mucor circinelloides, Mucor racemosus and related species. Studies Mycol. 1976;12:1–40.Google Scholar
- Bartnicki-García S, Nickerson WJ. Nutrition, growth and morphogenesis of Mucor rouxii. J Bacteriol. 1962;84:841–58.PubMed CentralPubMedGoogle Scholar
- Lasker BA, Borgia PT. High-frequency heterokaryon formation by Mucor racemosus. J Bacteriol. 1980;141:565–9.PubMed CentralPubMedGoogle Scholar
- Quiles-Rosillo MD, Torres-Martínez S, Garre V. cigA, a light-inducible gene involved in vegetative growth in Mucor circinelloides is regulated by the carotenogenic repressor crgA. Fungal Genet Biol. 2003;38:122–32.View ArticlePubMedGoogle Scholar
- Sambrook J, Russell DW. Molecular Cloning: A Laboratory Manual. 3rd ed. New York: Cold Spring Harbor Laboratory Press; 2001.Google Scholar
- Ausubel FM, Brent R, Kingston RE, Moore DD, Seidmann JG, Smith JA, et al. Currents Protocols in Molecular Biology. New York: John Wiley and Sons; 1987.Google Scholar
- Li B, Dewey CN. RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011;12:323.View ArticlePubMed CentralPubMedGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.