Physiological stressors and invasive plant infections alter the small RNA transcriptome of the rice blast fungus, Magnaporthe oryzae
- Vidhyavathi Raman†1,
- Stacey A Simon1, 2,
- Amanda Romag1, 3,
- Feray Demirci1, 2,
- Sandra M Mathioni1, 4,
- Jixian Zhai1, 2,
- Blake C Meyers1, 2 and
- Nicole M Donofrio1Email author
© Raman et al.; licensee BioMed Central Ltd. 2013
Received: 6 August 2012
Accepted: 2 May 2013
Published: 12 May 2013
The rice blast fungus, Magnaporthe oryzae is a destructive pathogen of rice and other related crops, causing significant yield losses worldwide. Endogenous small RNAs (sRNAs), including small interfering RNAs (siRNAs) and microRNAs (miRNAs) are critical components of gene regulation in many eukaryotic organisms. Recently several new species of sRNAs have been identified in fungi. This fact along with the availability of genome sequence makes M. oryzae a compelling target for sRNA profiling. We have examined sRNA species and their biosynthetic genes in M. oryzae, and the degree to which these elements regulate fungal stress responses. To this end, we have characterized sRNAs under different physiological stress conditions, which had not yet been examined in this fungus.
The resulting libraries are composed of more than 37 million total genome matched reads mapping to intergenic regions, coding sequences, retrotransposons, inverted, tandem, and other repeated regions of the genome with more than half of the small RNAs arising from intergenic regions. The 24 nucleotide (nt) size class of sRNAs was predominant. A comparison to transcriptional data of M. oryzae undergoing the same physiological stresses indicates that sRNAs play a role in transcriptional regulation for a small subset of genes. Support for this idea comes from generation and characterization of mutants putatively involved in sRNAs biogenesis; our results indicate that the deletion of Dicer-like genes and an RNA-Dependent RNA Polymerase gene increases the transcriptional regulation of this subset of genes, including one involved in virulence.
Various physiological stressors and in planta conditions alter the small RNA profile of the rice blast fungus. Characterization of sRNA biosynthetic mutants helps to clarify the role of sRNAs in transcriptional control.
RNA silencing, the RNA-mediated suppression of a target gene in a sequence specific manner, has been described in many organisms as post-transcriptional gene silencing, RNA interference, and quelling, in plants, animals, and fungi, respectively [1, 2]. These silencing pathways are involved in transposon silencing, viral defense, DNA elimination, heterochromatin formation, and post-transcriptional repression of cellular genes [3, 4]. The critical steps of these pathways include production of small RNAs (sRNAs) of 19–27 nucleotide (nt) from structured or double-stranded RNA by RNase III-like endonucleases called Dicers, followed by loading of these sRNAs into Argonaute-containing complexes to form RNA-induced silencing complexes (RISC) that guide the cleavage of target transcripts [3, 5].
Magnaporthe oryzae (hereafter M. oryzae), a filamentous fungal pathogen that causes rice blast disease, leads to severe yield losses in cultivated rice worldwide. The Magnaporthe-rice pathosystem is an excellent model to study the molecular basis of pathogenesis, in part because of the economic importance and molecular and genetic tractability of both host and pathogen, availability of genome sequences and transcriptome data, and the similarity of the rice blast fungus with other cereal diseases in terms of appressorium formation and intra-cellular tissue invasion –. The blast fungus infects rice plants at all stages and all tissues – leaves, stems, nodes, panicles and roots [8, 9]. Leaf infection starts with reversible attachment of conidia to the host surface. The release of spore tip mucilage from conidia forms an irreversible attachment to the surface . Conidia then germinate and form infection structures called appressoria at the end of the germ tubes. Appressoria generate enormous amounts of turgor pressure to penetrate the outer surface of the plant, which is followed by tissue invasion and colonization resulting in disease lesions .
The carefully controlled and timed life cycle of M. oryzae warrants examination into both transcriptional and post-transcriptional gene regulation. To date, endogenous small RNA species have been examined in several fungal species including the budding yeast, Saccharomyces cerevisiae, the model ascomycete fungus and close relative of M. oryzae, Neurospora crassa and M. oryzae itself (; reviewed in –). The M. oryzae genome encodes putative RNAi pathway components consisting of two Dicer-like proteins (DCL), three Argonautes (AGO), three RNA-dependent RNA polymerases (RdRP) and one exportin-5 [11, 13, 14]. In 2004, Kadotani et al. showed that Magnaporthe Dicer-like-2 (MoDcl2) is necessary for generating siRNAs in hairpin dsRNA-mediated RNA silencing, unlike their counterparts in Neurospora crassa in which both of the Dicers are redundantly involved in the silencing process . While the Kadotani et al. study (2004) did not find a clear role for the Magnaporthe Dicer-like-1 (MoDcl1) in the generation of siRNAs, the function of this gene is yet to be discovered. Exportin-5, which plays a role in nuclear export of pre-miRNAs in animals, was identified and deleted in M. oryzae and shown to have reduced pathogenicity on roots, and to complement its yeast homolog, the nucleocytoplasmic transporter Msn5 –.
The complexity of RNAi components in this fungus suggests the possible existence of multiple sRNA pathways. Though miRNAs have not been identified in M. oryzae, the published data have demonstrated numerous components of the M. oryzae sRNA repertoire, including small interfering (siRNA) -induced silencing of transgenes , the presence of viral derived siRNAs , methylguanosine-capped and polyadenylated small RNAs (CPA-sRNAs)  and sRNAs matching to repetitive elements, tRNA loci, rRNAs, protein coding genes, snRNAs and intergenic regions (; this study). Moreover, several new classes of small RNAs have been described recently in fungi, including QDE-2-interacting siRNAs (qiRNAs), microRNA-like small RNAs (milRNAs) and Dicer-independent small interfering RNAs (disiRNAs) of N. crassa. The existence of RNAi has been demonstrated in many fungal species including pathogenic species –. In Mucor circinelloides, production of two classes of small antisense RNAs requires a functional Dcl-2 gene .
In the M. oryzae life cycle, transition from the biotrophic to the necrotrophic stage is likely induced by environmental cues that the fungus encounters during the pathogenic process . Conditions encountered by the fungus during its growth in rice leaves, particularly during disease symptom expression, include nitrogen starvation, and this is regarded as one of the inductive cues for disease symptom expression during rice plant infections [26, 27]. Starvation is also considered a key factor for inducing the germ tube tip to differentiate into an appressorium . These and similar studies prompted us to examine whether post-transcriptional changes occur in M. oryzae during exposure to different physiological stressors, as well as during several invasive stages of pathogenic growth. Recent studies on genome-wide transcriptional changes demonstrate that in vitro conditions such as nitrogen and carbon deprivation are comparable to genome-wide changes that occur during invasive plant growth . The role of post-transcriptional control during these stresses, however, has not yet been studied. We examined small RNA profiles of M. oryzae during stress conditions matching those of Mathioni et al. , namely, nitrogen starvation, carbon starvation, minimal media, and paraquat-induced oxidative stress. We also examined small RNA fungal profiles in planta at 72 and 96 hours post inoculation (hpi), as well as mock inoculated plants. Small RNAs from certain conditions showed association with specific genomic loci; for example we observed an abundance of sRNAs associated with retrotransposons during nitrogen starvation. We also identified a small subset of genes whose transcriptional down-regulation during several of the stress conditions showed association with high numbers of sRNAs. To further examine whether this represented a true link between transcriptional and post-transcriptional control, we generated four small RNA biosynthetic fungal mutants, and performed quantitative real-time RT-PCR on this down-regulated subset of genes. These results revealed a previously unknown role for the Magnaporthe Dicer-Like1 (MoDcl1) gene.
Small RNA profiling of M. oryzaestress libraries
Description, total sequences, genome-matched reads and distinct reads for each sRNA library
Distinct genome- matched readsa
Mycelium, re-inoculated into complete medium for 16 hours
Mycelium, re-inoculated into carbon starved medium for 16 hours
Mycelium, re-inoculated into minimal medium for 16 hours
Mycelium, nitrogen-starved for 16 hours
Mycelium, paraquat-treated for 24 hours
Leaf, 0 hr control for Mg-treated libraries, Nipponbare
Leaf, 72 hpi, Magnaporthe infected, Nipponbare
Leaf, 96 hpi, Magnaporthe infected, Nipponbare
Cluster analysis of mycelial and in plantalibraries
Cluster analysis of all small RNA libraries compared between two conditions
No. of genes associated
CM vs CS
CM vs MM
CS vs MM
NS vs CM
CM vs PQ
Infected Rice leaf libraries
LMg0 vs 72
LMg0 vs 96
LMg72 vs 96
sRNA associate with different genomic regions, depending upon the fungal stress
Total, and distinct genome-matched sRNAs among all libraries for different classes of genomic elements
Antisense to protein-coding genes
• Low complexity
• Simple repeat
Noncoding RNAs (tRNAs or rRNAs)
Antisense to noncoding RNAs
To determine whether our libraries identified repeat-associated siRNAs (known as rasiRNAs), we analyzed six different classes of repeats, and found that the majority of the sRNAs emanated from rRNA loci (89%) followed by retrotransposons (LTR/Gypsy and LINE/Tad1; 11%). Small RNAs that matched to transposons (DNA/TcMar-Pogo), simple repeats, and low complexity repeats occupied only a small fraction of the total sRNAs. Though rRNA-associated sRNAs are more abundant compared to retrotransposons-associated sRNAs, the distinct reads associated with the former were less (43%) compared to the latter (52%). The majority of sRNAs matching to protein-coding regions were in a sense orientation (77%; Table 3), and a similar phenomenon was observed for non-coding RNAs (tRNAs or rRNAs; 87%).
To further characterize the sRNAs associated with retrotransposons, representative loci from major repeat classes in M. oryzae such as LTR/Gypsy and LINE/Tad1 retrotransposons were analyzed for their size distribution, changes under different physiological conditions, strand specificity and 5’ nt identity. Strand orientation of sRNAs was also included in our analysis because abundant antisense siRNAs are known to be associated with repeats  and are affected by the loss of RNAi machinery . Also, studying strand specificity is important because sRNAs originating from overlapping gene pairs or natural cis-antisense gene pairs have been shown to play an important role in salt tolerance in Arabidopsis  and gene regulation in mouse oocytes . LTR/Gypsy retrotransposons showed major peaks at 20–23 nt and abundant sRNAs in the antisense orientation (Additional file 2: Figure S2A and B). These sRNAs differentially accumulated under different physiological stress conditions; for example, one class of LTR/GYPSY retrotransposon called GYMAG2 showed a spike in sRNA accumulation during nitrogen and oxidative stress conditions that was roughly 28% and 50% higher than any other condition (Additional file 2: Figure S2D). Analysis of the next most abundant major repeat class of retrotransposons, LINE/Tad1, revealed that its major peak for sRNAs based on abundance, was at 20 nt followed by 23/24 nt (Additional file 3: Figure S3A). Similar to the LTR/Gypsy retrotransposon class, sRNAs matching to LINE/Tad1 repeats also showed differential accumulation in various physiological stress conditions, but unlike the LTR class, LINE/Tad1 sRNAs did not show a spike under any one particular condition (Additional file 3: Figure S3C). Unlike LTR/Gypsy, sRNAs matching to LINE/Tad1 did not show abundant sRNAs in the antisense orientation (Additional file 3: Figure S3B).
We observed that many of the annotated tRNAs and 5S rRNAs represent a significant portion of intergenic region -associated sRNAs. In order to study these genomic features, we selected random loci of each feature and analyzed their size distribution. The results showed a bimodal distribution for sRNAs originating from a majority of the rRNA loci, with size peaks at 20 and 23/24 nt (Additional file 4: Figure S4). Analysis of random tRNA loci revealed the presence of varied size distribution and increased sRNAs under the PQ condition (Additional file 4: Figure S4). We also observed strand bias for sRNAs originating from 5S rRNA and tRNA. All analyzed 5S rRNA and tRNA loci showed sRNAs on the sense strand (data not shown). sRNAs originating from mitochondrial sequences occupied only a small proportion, comprising 0.6% and 1.2% of total and distinct genome-matched reads, respectively (Table 3).
Stress conditions altered the proportion of individual classes of sRNAs associated with different genomic features (Additional file 5: Figure S5A). The PQ library revealed a notable reduction in both total and distinct abundance sRNAs matching to protein coding regions, compared to all the other libraries. Distinct abundance for the PQ library also showed an increase in sRNAs associated with intergenic regions, and with repeats. Another prominent difference was observed in the CS library; there was an increase in sRNAs associated with LTR Gypsy transposons, and a decrease in sRNAs associated with rRNA repeats (Additional file 5: Figure S5C). An increase in LTR retrotransposons was also observed during infection of rice plants with M. oryzae (Additional file 5: Figure S5D).
5’ nucleotide preference
Correlation between transcriptional and post-transcriptional regulation under stress conditions
Features of sRNAs associated with transcriptionally down-regulated genes under different physiological stress conditions
MGG_01439.7 inorganic phosphate transporter PHO84
MGG_04470.7 nucleolar complex protein 14
MGG_01439.7 inorganic phosphate transporter PHO84
MGG_16711.7 C Hypothetical protein
Features of sRNAs associated with fungal genes down-regulated during infection
MGG_02390.7 conserved hypothetical protein
MGG_06609.7 acetyl-CoA hydrolase
MGG_09965.7 conserved hypothetical protein
MGG_08843.7 magnesium transporter ALR2
MGG_01991.7 betaine aldehyde dehydrogenase
MGG_06868.7 acetolactate synthase catalytic subunit
MGG_04428.7 zinc finger transcription factor ace1
MGG_01596.7 DNA damage response protein kinase DUN1
has Tandem repeats on it
Only a small fraction of RNA species (2-3%) are translated into proteins in higher organisms, while the remainder includes non-coding yet functional species such as small RNAs [39, 40]. Small RNAs in eukaryotes have some well-established roles in modulating gene expression at the transcriptional and post-transcriptional levels, as well as in shaping the organization and modification of chromatin . There is a robust body of knowledge on sRNAs in the model fungus Neurospora crassa, and in the budding yeast (; reviewed in [10, 42]). Recently, several detailed studies on sRNAs have emerged in fungal pathogens, including M. oryzae amongst others [11, 12, 20, 21, 37]. The current study presented here, adds to this growing information by analyzing sRNAs under various physiological stress conditions and during plant infection. Understanding the epigenetic response of plant pathogens to different environmental conditions may elucidate how these pathogens so successfully adapt to new climates and/or new hosts.
Our results demonstrate that the majority of the sRNAs originate from intergenic regions followed by repeats, protein coding regions and mitochondrial regions, when sRNAs from all libraries were taken together. Our results are supported by those of Nunes et al. , who reported abundant sRNAs matching to intergenic regions, transposable elements, rRNA, tRNA, and protein-coding genes, from mycelial and appressorial libraries. Transposons account for about 7% of the M. oryzae genome (including 5.4% contribution by retrotransposon; ); our data suggests that only 4.5% of the total sRNA population matches to transposons (11% of the repeats population), which indicates RNAi may not act as the only mechanism to suppress transposons in M. oryzae. This observation is further supported by the hypothesis of Murata et al. , which supports the operation of a transposon- or repeat-silencing mechanism distinct from both transcriptional gene silencing and siRNA-mediated gene silencing.
After spore germination and appressorium formation, M. oryzae generates a penetration peg that enters into the plant’s epidermal cell. Primary infection hyphae colonize this cell followed by formation of bulbous secondary hyphae and spreading to adjacent epidermal cells (the biotrophic phase; ). At about 72 hpi, secondary hyphae spread into mesophyll cells and at this time approximately 10% of total biomass is fungal , which could account for the reduced amount of total and distinct genome matched reads in in planta libraries, compared to the mycelial libraries. Under favorable disease conditions, visible lesions are formed by 96 hpi . The drastic changes undergone by the pathogen during invasive, later stage growth are reflected in the increased number of genes showing differential abundance for sRNAs during 72 and 96 hpi, as determined by the cluster analysis (Table 2). During these later invasive stages, the fungus is believed to be coping with nutrient-limited environments [27, 29]. In the in planta libraries, we observed 562 distinct reads associated with LMg72 compared to LMg96 after the removal of LMg0, which represent rice-only reads (Additional file 9: Figure S7). Furthermore, out of the ~8000-9000 sRNA reads found in both the LMg72 and LMg96 libraries, only 2132 overlapped, indicating that many unique small RNAs are generated during these invasive, later growth stages. We took advantage of previous transcriptional data generated from our lab in 2011 to determine whether clusters with sRNAs have a negative correlation with significantly differentially expressed genes in the microarray study. Ninety-eight such genes from the microarray study associated with sRNAs and 37.7% of these showed a negative correlation indicating that these particular genes may be regulated epigenetically. The remaining ~60% of the microarray genes showed a positive correlation with the sRNA data. While these are relatively low numbers from a genome-wide perspective, the microarray data was generated from a very early version of the M. oryzae genome annotation, which has changed considerably in the last four years; furthermore, the 98 genes from the microarray study were based upon highly stringent criteria used both in this study, as well as the study from 2011, limiting the number of such genes we can identify.
Small RNAs differentially associate with retrotransposons, in a stress-dependent manner
Living organisms adapt to the environmental stress conditions by undergoing several phenotypic, physiologic and molecular changes, which requires precise regulation at the transcriptional, post-transcriptional and translational levels. These changes include but are not limited to activation and repression of several genes , signaling cascades , changes in molecular chaperones , activation of transposable elements [49, 50] and chromatin changes and epigenetic modifications . Several nutrient stress responsive sRNAs, have been identified in plants [52, 53], bacteria  and Caenorhabditis elegans.
In pathogenic fungal species, efficient pathogen nutrition is essential for successful colonization and fungal fitness . Transcriptome studies on M. oryzae and many other fungi have demonstrated that nutrient acquisition is a tightly regulated process, involving genes that function in amino acid metabolism and uptake, global nitrogen regulation and pathogenicity [27, 57]–. In the present study, we found that repeat-associated sRNAs generally and LTR-retrotransposons specifically were induced by oxidative stress and nitrogen starvation, respectively (Additional file 5: Figure S5). This observation may be similar to Arabidopsis roots undergoing phosphate starvation, which induces a high level of DICER-LIKE1-dependent sRNAs derived from the long terminal repeat of a retrotransposon and a group of 19 nt sRNAs that correspond to the 5’ end of tRNA . In rice, some repeat-derived sRNAs are down-regulated by abscisic acid and abiotic stresses accompanied by an increased accumulation of their precursors .
In M. oryzae, MAGGY retrotransposons are activated by heat shock and oxidative stress along with up-regulation of MAGGY RNA  and an EST matching to a non-LTR-retrotransposon was identified in an M. oryzae-infected rice leaf library . Similarly, our results revealed a spike in sRNAs associated with an LTR class during oxidative stress. Several ESTs matching to LTR-retrotransposons have also been identified in the in planta stage of the wheat fungal pathogen, Mycosphaerella graminicola, and are less abundant compared to their mycelial stage . Similarly, our data revealed that distinct sRNAs reads matching to repeats in in planta growth were ~10% higher than in the mycelial libraries (Additional file 5: Figure S5A and B). Reactivation of transposable elements during stress conditions have been found in many organisms from fission yeast to plants and animals (reviewed in ). Recently, a widely-studied tobacco transposon called TNT1 was found to be reactivated during infection of tobacco . While we did see an increase in sRNAs associated with LTR/Gypsy retrotransposons during nitrogen starvation compared to the non-stressed condition, overall, there was not a striking difference in these types of associations in the case of the LINE/Tad1 class of retrotransposon, where sRNA abundance was not markedly different in the stress versus non-stressed conditions. Together, our data and the wealth of information on transposon-associated small RNAs in other organisms suggest that retrotransposons might be suppressed via epigenetic mechanisms. Additional testing will be required to determine what underlies the increase in sRNAs associated with LTR elements during nitrogen starvation.
Interestingly, we noted an abundance of sRNAs on the antisense strands of LTR/Gypsy elements, while there was a roughly equal distribution of them on either strand of the LINE/Tad1 elements. This was different from the recent study by Nunes et al.  on sRNA species in M. oryzae, which showed a more equal distribution of siRNAs on both LTR-associated strands, and may be attributable to the differing conditions used in each study and differing sequencing depth. Association of small RNAs with one strand or the other has been well-studied in Drosophila, for example, antisense siRNAs are abundant for repeat regions . In general, transcription of transposons is activated by stress conditions; for example, carbon or nitrogen stresses activate the transcription of Drifter, a novel low copy transposon in Fusarium oxysporum. Some of the transposons are transcribed in both directions  and environmental stress-regulated antisense transcription of Ty1 LTR-retrotransposons is observed in S. cerevisiae.
Derivation of the sRNAs in M. oryzae: antisense and 5’ nucleotide position
Around 25% of the tags that matched to genes are in antisense orientation (Table 3). Interestingly, one-fifth of the annotated genes of this fungus encode both sense and antisense transcripts . The sources of the dsRNAs for Dicer processing include inverted repeats, bidirectional transcription and antisense transcripts from various loci . There is a possibility that overlapping transcripts in an antisense orientation produced from a pair of neighboring genes on opposite DNA strands (convergent overlapping gene pairs or natural cis-antisense gene pairs) form dsRNAs, which may be further processed into sRNAs [33, 34]. Recently, antisense transcripts were also found in other fungal species like the ectomycorrizal ascomycete Tuber melanosporum and basidiomycete mushroom Schizophyllus commune. It is also proven that siRNAs matching to mRNAs in S. castellii and in the RNAi reconstituted strain of S. cerevisiae were indeed produced from dsRNA formed by antisense and overlapping transcription [5, 72]. Even in M. circinelloides, the majority of the exonic-siRNAS (ex-siRNAs) are antisense to mRNA . Together, these results suggest that at least part of the siRNAs that matched to genes are produced from dsRNA formed by overlapping sense and antisense transcripts in M. oryzae. We observed that both 5S rRNA and tRNA loci were producing more of sense strand specific sRNAs whereas pseudo-tRNA loci have both sense and antisense sRNAs (data not shown). It has been shown that pseudo-tRNA of Bacillus cereus is detected to have a function outside protein synthesis . Moreover, a subset of pseudogenes generates endo-siRNAs from double-stranded RNAs or directly in mouse oocytes . The same scenario may exist in M. oryzae, because around 40% of pseudo-tRNAs are present either within a gene or opposite strand to it.
In the RNAi pathway, the Dicer-cleaved sRNAs are bound by Argonautes, and based on their 3’ complementarity Argonautes exert homology dependent transcriptional cleavage or post-transcriptional silencing . The binding of sRNAs to AGO is mainly dependent on structural features of precursor duplexes in animals [76, 77] whereas in plants, it is determined by the 5’ nt of the sRNA. In Arabidopsis AGO2 and AGO4 preferentially recruit sRNAs with 5’ A, likewise AGO1 with U and AGO5 with C . Abundant U at the 5’ end observed in all the mycelial libraries in our study may be biologically significant (Figure 3), since this kind of preference was observed in other organisms including the single cell green alga Chlamydomonas. Importantly, a preference for U (T) and suppression of C was observed in M. oryzae in a recent study, which our data supports . In the fungal kingdom, 5’ preference for U is observed in qde-2 (Argonaute) associated sRNAs and milRNAs, and disiRNAs of N. crassa; ago-1 associated sRNAs of S. pombe; sRNAs of S. castelleii and K. polysporus; and Ex-siRNAs of Mucor circinelloides.
A role for Dicer1 (MoDcl1) in M. oryzae?
Though Dicer-like (DCL) proteins are evolutionarily conserved in eukaryotes, their numbers vary among organisms. In the multiple DCL protein-containing organisms such as Drosophila and A. thaliana, they function in distinct RNA silencing pathways with overlapping specificity . Both M. oryzae and its close relative N. crassa encode two Dicers. However, MoDcl2 in M. oryzae functions in siRNA production under normal conditions, whereas N. crassa Dicers, orthologs of the M. oryzae Dicers, act redundantly. While it is likely that MoDcl1 does not have the ability to render or exert a silencing role under normal conditions either triggered by hairpin RNA-expressing transgenes or endogenous repetitive sequences, it is capable of producing siRNAs when overexpressed in an MoDcl2 deletion background. Moreover, its mRNA levels were up-regulated by 10- to 15- fold during the sexual stage. It has been proven that both MoDcl1 and MoDcl2 define functionally diversified RNAi pathways which may have arisen from transcriptional control and protein specialization .
In M. oryzae, deletion of MoDcl2 leads to increased abundance of MAGGY mRNA when it is introduced into the native genome . With these reports in mind, it was expected that deletion of MoDcl2 in our experiments would increase the transcript abundance of genes which are known to possess more sRNAs in one or more libraries, in comparison with mycelia grown in complete media (control condition). Conversely, our results showed there is an overall reduction in transcript levels of these genes in the Δmodcl2 mutant and an increased level in the Δmodcl1 mutant. Interestingly, the known virulence gene ACE1 was among the few genes that fit this profile. ACE1 is tightly regulated at the transcriptional level; transcripts are only found in the penetration peg of the fungus, a specialized hypha that allows ingress through the surface of the leaf and into the plant cell . We postulate that strong control of this gene might be governed at the post-transcriptional level. Further studies will reveal whether these and the other genes tested show the same profile when the biosynthetic mutants are subjected to the environmental stress conditions, or during pathogenicity. Several studies have shown that MoDcl2 is exclusively involved in siRNA production for transgenes and active transposons in M. oryzae, while less active transposons are not silenced by this same mechanism [43, 80, 82]. Moreover, MoDcl1 is also involved in silencing when it is over-expressed in a Δmodcl2 mutant background. Together, these findings along with results presented here suggest the hypothesis that MoDcl1, at least in M. oryzae, might define the existence of a MoDcl2-independent RNAi mechanism. The recent discovery of Argonaute associated Dicer independent siRNAs (disiRNAs) in N. crassa, the first discovered Dicer-independent siRNAs in the fungal kingdom, supports this hypothesis. Surprisingly, in M. circinelloides, Dicer2 plays a critical role in transgene-induced silencing , but redundancy for the two Dicer genes was observed for production of endogenous sRNAs with DCL2 being the primary protein responsible for producing the majority of endogenous sRNAs . Taken together, these studies suggest that M. oryzae presumably has similar redundant Dicer-based mechanisms for the production of sRNAs.
In conclusion, various stress-inducing and in planta conditions alter the sRNA profile of M. oryzae, with the main effect exerted on loci associated with protein coding regions. A subset of these genes showed an association with the transcriptome data, as well as a differential accumulation of their transcripts between the two Dicer-deficient mutant backgrounds. Stress conditions repressed the sRNAs associated with protein coding regions, and marginally increased the intergenic- and repeat- associated sRNAs. Considering our results and previous reports, questions of Dicer-independent mechanisms, and genetic redundancy of multiple copies of sRNA biosynthetic genes in the genome, warrant a detailed dissection of sRNA generating pathway(s) through the use of various sRNA biosynthetic mutants. Further studies on in planta interaction using sRNA biosynthetic mutants would provide better understanding of molecular mechanisms behind host-pathogen interaction which in turn useful for developing novel disease control strategies.
Fungal strain and growth conditions
M. oryzae wild type strain 70–15 was used in this study because of the availability of genome sequence and transcriptome data. Filter paper stocks stored at -20°C were inoculated onto Oatmeal agar plates (OM: oatmeal 50 g/L and agar 15 g/L) for germination and to establish starter cultures. Fungus grown on OM plates were transferred to complete medium (CM: 10 g sucrose, 6 g yeast extract, 6 g casamino acid and 1 mL trace elements/L) plates and subsequently to complete liquid medium. Liquid grown fungus was harvested under sterile conditions and washed with sterile distilled water. The resulting fungal mat was carefully divided into equal pieces for further exposure to physiological conditions in liquid such as complete medium (CM), minimal medium (MM), carbon starved medium (CS: MM lacking carbon source), nitrogen starved medium (NS: MM lacking nitrogen source) and a reactive oxygen species-rich environment generated by amendment of Paraquat (PQ) to the CM, as described in detail by Mathioni et al. . For the in planta libraries, twenty-day-old rice (Nipponbare) seedlings were mock inoculated with 0.2% gelatin and leaves were harvested immediately (LMg0). Conidia from M. oryzae 70–15 grown on OM plates for 10 days (d) were harvested using 0.2% gelatin. Concentration of the conidial suspension was adjusted to 1 × 105 spores/mL and 3 mL of the suspension was spray inoculated onto rice plants; leaves were harvested after 72 hours post-inoculation (hpi) and 96 hpi respectively for LMg72 and LMg96 libraries.
Small RNA library construction and data analysis
Total RNA was isolated using Trizol reagent (Invitrogen, New York, USA) according to manufacturer’s instructions. The quality of the RNA was checked by agarose gel electrophoresis and NanoDrop ND-1000. High molecular weight RNA was removed by PEG precipitation, and low molecular weight RNA was size-fractionated on a polyacrylamide gel. Gel portions corresponding to 20–30 nt were recovered and small RNA libraries were constructed by following the protocol by Lu et al.  with modifications using standard Illumina sRNA adapters. The quality of the resulting libraries was checked by subjecting an aliquot of each library to TOPO® cloning (Invitrogen, USA), and approximately 96 clones from each library were sequenced by the Sanger method. The resulting sequences were checked for the presence of adapters, size of cDNAs, and their genome match. The libraries were then sequenced using Illumina’s Genome Analyzer II (GAII) at the Delaware Biotechnology Institute. Adapters were removed from the resulting sequences and matched to the M. oryzae genome (version 6 assembly) from the Broad Institute, and reads with sequencing errors were filtered out using PERL scripts. tRNA, rRNA, snRNA or snoRNAs -matched reads were removed from the total reads, and genome-matched reads were obtained using the program Bowtie, accepting only perfect matches. From the total genome matched reads, distinct genome matched reads were calculated. The abundance of each read was normalized to 5 million units of transcripts per 5 million (TP5M) to account for the difference in the sequencing depth and to facilitate accurate comparisons across libraries. To compare the different libraries, the entire genome was divided arbitrarily into 500 bp windows/clusters based on the proximity-based algorithm (; described in the Results section). Clusters showing ≥ 10-fold difference between 2 or more libraries were selected for further cluster analysis. Repeat analysis was performed using programs namely, RepeatMasker , etandem and einverted .
Targeted disruption of genes involved in sRNA biogenesis and mutant screening
Gene replacement cassettes were constructed using adaptamer-mediated PCR . Typically, 1.3 kb of upstream and downstream sequence of each target gene was amplified with primers that contained adaptamer sequences. A 1.4 kb fragment containing the hygromycin B phosphotransferese gene with the trpC promoter from A. nidulans was amplified from plasmid pCB1003 using the adaptamer sequence attached to the forward HPHF and reverse HPHR primer set (Additional file 11: Table S3). Using nested primers, three individually amplified fragments of upstream and downstream sequence and hygromycin resistance gene were combined and amplified together to construct a hygromycin cassette for gene replacement of approximately 3.3 kb in length. The hygromycin cassette was transformed into 70–15 protoplasts as previously described . Gene disruption mutants were identified by PCR screening using primers outside the flanking regions and gene specific primers (Additional file 11: Table S3) and further confirmed by Southern blot analysis (Additional file 10: Figure S8). For deletion of MoDcl1 and MoDcl2, gene replacement cassettes were constructed as described previously with the modification of BAR gene, amplified from pCB1530, being used as the selectable marker  in place of HYG. Gene replacement cassettes for MoDcl1 and MoDcl2 were transformed into protoplasts of MoDcl2 and MoDcl1 knockouts, respectively.
Quantitative reverse transcription-polymerase chain reaction (qRT-PCR)
Genomic DNA was removed from total RNA by TURBO® DNase treatment (Ambion, New York, USA). Reverse transcription reactions were performed with 500 ng of RNA in 20 μl reaction using High Capacity RNA-to-cDNA Kit (Applied Biosystems, California, USA) by following the manufacturer’s instructions. Real time PCR was done using 5 PRIME SYBR mastermix for Real-time PCR (5 PRIME, Maryland, USA) in an Eppendorf Mastercycler using 2 μl of diluted reverse transcribed templates. All samples were repeated three times. Relative expression values were calculated using the 2-ΔΔCT method with GADPH (MGG_01084.6) as a housekeeping control. Primers used are listed in Additional file 12: Table S4.
GEO accession numbers
The data obtained from small RNA deep sequencing studies were deposited in the Gene Expression Omnibus database at NCBI (http://www.ncbi.nlm.nih.gov/geo/). The accession number for GEO is GSE43277). The data can also be found at http://mpss.udel.edu/mg_sRNA.
The authors thank James Sweigard for providing the plasmids pCB1003 and pCB1530 for generating the single and double knock out mutants. Further, they wish to thank members of the Donofrio lab for guidance with the wet lab work, and members of the Meyers group for help with additional bioinformatics. NMD would like to thank Ralph Dean for valuable discussions on data interpretation and generation.
This project is supported by a United States Department of Agriculture – Agriculture and Food Research Initiative grant to ND and BCM (# 2008-35319-04465).
- Romano N, Macino G: Quelling: transient inactivation of gene expression in Neurospora crassa by transformation with homologous sequences. Mol Microbiol. 1992, 6 (22): 3343-3353. 10.1111/j.1365-2958.1992.tb02202.x.View ArticlePubMedGoogle Scholar
- Baulcombe D: RNA silencing. Curr Biol. 2002, 12 (3): R82-4. 10.1016/S0960-9822(02)00665-6.View ArticlePubMedGoogle Scholar
- Tomari Y, Zamore PD: Perspective: machines for RNAi. Genes Dev. 2005, 19 (5): 517-529. 10.1101/gad.1284105.View ArticlePubMedGoogle Scholar
- Baulcombe D: RNA silencing in plants. Nature. 2004, 431 (7006): 356-363. 10.1038/nature02874.View ArticlePubMedGoogle Scholar
- Drinnenberg IA, Weinberg DE, Xie KT, Mower JP, Wolfe KH, Fink GR, Bartel DP: RNAi in budding yeast. Science. 2009, 326 (5952): 544-550. 10.1126/science.1176945.PubMed CentralView ArticlePubMedGoogle Scholar
- Dean RA, Talbot NJ, Ebbole DJ, Farman ML, Mitchell TK, Orbach MJ, Thon M, Kulkarni R, Xu JR, Pan H, Read ND, Lee YH, Carbone I, Brown D, Oh YY, Donofrio N, Jeong JS, Soanes DM, Djonovic S, Kolomiets E, Rehmeyer C, Li W, Harding M, Kim S, Lebrun MH, Bohnert H, Coughlan S, Butler J, Calvo S, Ma LJ, Nicol R, Purcell S, Nusbaum C, Galagan JE, Birren BW: The genome sequence of the rice blast fungus Magnaporthe grisea. Nature. 2005, 434 (7036): 980-986. 10.1038/nature03449.View ArticlePubMedGoogle Scholar
- Ebbole DJ: Magnaporthe as a model for understanding host-pathogen interactions. Annu Rev Phytopathol. 2007, 45: 437-456. 10.1146/annurev.phyto.45.062806.094346.View ArticlePubMedGoogle Scholar
- Wilson RA, Talbot NJ: Under pressure: investigating the biology of plant infection by Magnaporthe oryzae. Nat Rev Microbiol. 2009, 7 (3): 185-195. 10.1038/nrmicro2032.View ArticlePubMedGoogle Scholar
- Sesma A, Osbourn AE: The rice leaf blast pathogen undergoes developmental processes typical of root-infecting fungi. Nature. 2004, 431 (7008): 582-586. 10.1038/nature02880.View ArticlePubMedGoogle Scholar
- Fulci V, Macino G: Quelling: post-transcriptional gene silencing guided by small RNAs in Neurospora crassa. Curr Opin Microbiol. 2007, 10 (2): 199-203. 10.1016/j.mib.2007.03.016.View ArticlePubMedGoogle Scholar
- Kadotani N, Nakayashiki H, Tosa Y, Mayama S: RNA silencing in the phytopathogenic fungus Magnaporthe oryzae. Mol Plant Microbe Interact. 2003, 16 (9): 769-776. 10.1094/MPMI.2003.16.9.769.View ArticlePubMedGoogle Scholar
- Nunes CC, Gowda M, Sailsbery J, Xue M, Chen F, Brown DE, Oh Y, Mitchell TK, Dean RA: Diverse and tissue-enriched small RNAs in the plant pathogenic fungus. Magnaporthe oryzae. BMC Genomics. 2011, 12: 288-10.1186/1471-2164-12-288.View ArticlePubMedGoogle Scholar
- Nakayashiki H: RNA silencing in fungi: mechanisms and applications. FEBS Lett. 2005, 579 (26): 5950-5957. 10.1016/j.febslet.2005.08.016.View ArticlePubMedGoogle Scholar
- Murphy D, Dancis B, Brown JR: The evolution of core proteins involved in microRNA biogenesis. BMC Evol Biol. 2008, 8: 92-10.1186/1471-2148-8-92.PubMed CentralView ArticlePubMedGoogle Scholar
- Catalanotto C, Pallotta M, ReFalo P, Sachs MS, Vayssie L, Macino G, Cogoni C: Redundancy of the two dicer genes in transgene-induced posttranscriptional gene silencing in Neurospora crassa. Mol Cell Biol. 2004, 24 (6): 2536-2545. 10.1128/MCB.24.6.2536-2545.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Kaffman A, Rank NM, O'Neill EM, Huang LS, O'Shea EK: The receptor Msn5 exports the phosphorylated transcription factor Pho4 out of the nucleus. Nature. 1998, 396 (6710): 482-486. 10.1038/24898.View ArticlePubMedGoogle Scholar
- Bohnsack MT, Czaplinski K, Gorlich D: Exportin 5 is a RanGTP-dependent dsRNA-binding protein that mediates nuclear export of pre-miRNAs. RNA. 2004, 10 (2): 185-191. 10.1261/rna.5167604.PubMed CentralView ArticlePubMedGoogle Scholar
- Tucker SL, Besi MI, Galhano R, Franceschetti M, Goetz S, Lenhert S, Osbourn A, Sesma A: Common genetic pathways regulate organ-specific infection-related development in the rice blast fungus. Plant Cell. 2010, 22 (3): 953-972. 10.1105/tpc.109.066340.PubMed CentralView ArticlePubMedGoogle Scholar
- Himeno M, Maejima K, Komatsu K, Ozeki J, Hashimoto M, Kagiwada S, Yamaji Y, Namba S: Significantly low level of small RNA accumulation derived from an encapsidated mycovirus with dsRNA genome. Virology. 2010, 396 (1): 69-75. 10.1016/j.virol.2009.10.008.View ArticlePubMedGoogle Scholar
- Gowda M, Nunes CC, Sailsbery J, Xue M, Chen F, Nelson CA, Brown DE, Oh Y, Meng S, Mitchell T, Hagedorn CH, Dean RA: Genome-wide characterization of methylguanosine-capped and polyadenylated small RNAs in the rice blast fungus Magnaporthe oryzae. Nucleic Acids Res. 2010, 38 (21): 7558-7569. 10.1093/nar/gkq583.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee HC, Li L, Gu W, Xue Z, Crosthwaite SK, Pertsemlidis A, Lewis ZA, Freitag M, Selker EU, Mello CC, Liu Y: Diverse pathways generate microRNA-like RNAs and Dicer-independent small interfering RNAs in fungi. Mol Cell. 2010, 38 (6): 803-814. 10.1016/j.molcel.2010.04.005.PubMed CentralView ArticlePubMedGoogle Scholar
- Matityahu A, Hadar Y, Dosoretz CG, Belinky PA: Gene silencing by RNA Interference in the white rot fungus Phanerochaete chrysosporium. Appl Environ Microbiol. 2008, 74 (17): 5359-5365. 10.1128/AEM.02433-07.PubMed CentralView ArticlePubMedGoogle Scholar
- Mouyna I, Henry C, Doering TL, Latge JP: Gene silencing with RNA interference in the human pathogenic fungus Aspergillus fumigatus. FEMS Microbiol Lett. 2004, 237 (2): 317-324.PubMedGoogle Scholar
- Janus D, Hoff B, Kuck U: Evidence for Dicer-dependent RNA interference in the industrial penicillin producer Penicillium chrysogenum. Microbiology. 2009, 155 (Pt 12): 3946-3956.View ArticlePubMedGoogle Scholar
- de Haro JP, Calo S, Cervantes M, Nicolas FE, Torres-Martinez S, Ruiz-Vazquez 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 (10): 1486-1497. 10.1128/EC.00191-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Talbot NJ, McCafferty HRK, Ma M, Moore K, Hamer JE: Nitrogen starvation of the rice blast fungus Magnaporthe grisea may act as an environmental cue for disease symptom expression. Physiol Mol Plant Pathol. 1997, 50 (3): 179-195. 10.1006/pmpp.1997.0081.View ArticleGoogle Scholar
- Donofrio NM, Oh Y, Lundy R, Pan H, Brown DE, Jeong JS, Coughlan S, Mitchell TK, Dean RA: Global gene expression during nitrogen starvation in the rice blast fungus. Magnaporthe grisea. Fungal Genet Biol. 2006, 43 (9): 605-617. 10.1016/j.fgb.2006.03.005.View ArticlePubMedGoogle Scholar
- Liu XH, Lu JP, Lin FC: Autophagy during conidiation, conidial germination and turgor generation in Magnaporthe grisea. Autophagy. 2007, 3 (5): 472-473.View ArticlePubMedGoogle Scholar
- Mathioni SM, Belo A, Rizzo CJ, Dean RA, Donofrio NM: Transcriptome profiling of the rice blast fungus during invasive plant infection and in vitro stresses. BMC Genomics. 2011, 12: 49-10.1186/1471-2164-12-49.PubMed CentralView ArticlePubMedGoogle Scholar
- Lu C, Tej SS, Luo S, Haudenschild CD, Meyers BC, Green PJ: Elucidation of the small RNA component of the transcriptome. Science. 2005, 309 (5740): 1567-1569. 10.1126/science.1114112.View ArticlePubMedGoogle Scholar
- Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. 1996-2010, http://www.repeatmasker.org,Google Scholar
- Vagin VV, Sigova A, Li C, Seitz H, Gvozdev V, Zamore PD: A distinct small RNA pathway silences selfish genetic elements in the germline. Science. 2006, 313 (5785): 320-324. 10.1126/science.1129333.View ArticlePubMedGoogle Scholar
- Borsani O, Zhu J, Verslues PE, Sunkar R, Zhu JK: Endogenous siRNAs derived from a pair of natural cis-antisense transcripts regulate salt tolerance in Arabidopsis. Cell. 2005, 123 (7): 1279-1291. 10.1016/j.cell.2005.11.035.PubMed CentralView ArticlePubMedGoogle Scholar
- Watanabe T, Totoki Y, Toyoda A, Kaneda M, Kuramochi-Miyagawa S, Obata Y, Chiba H, Kohara Y, Kono T, Nakano T, Surani MA, Sakaki Y, Sasaki H: Endogenous siRNAs from naturally formed dsRNAs regulate transcripts in mouse oocytes. Nature. 2008, 453 (7194): 539-543. 10.1038/nature06908.View ArticlePubMedGoogle Scholar
- Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, Wu L, Li S, Zhou H, Long C, Chen S, Hannon GJ, Qi Y: Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5' terminal nucleotide. Cell. 2008, 133 (1): 116-127. 10.1016/j.cell.2008.02.034.PubMed CentralView ArticlePubMedGoogle Scholar
- Buhler M, Spies N, Bartel DP, Moazed D: TRAMP-mediated RNA surveillance prevents spurious entry of RNAs into the Schizosaccharomyces pombe siRNA pathway. Nat Struct Mol Biol. 2008, 15 (10): 1015-1023. 10.1038/nsmb.1481.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicolas FE, Moxon S, de Haro JP, Calo S, Grigoriev IV, Torres-Martinez S, Moulton V, Ruiz-Vazquez RM, 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 (16): 5535-5541. 10.1093/nar/gkq301.PubMed CentralView ArticlePubMedGoogle Scholar
- Fudal I, Collemare J, Bohnert HU, Melayah D, Lebrun MH: Expression of Magnaporthe grisea avirulence gene ACE1 is connected to the initiation of appressorium-mediated penetration. Eukaryot Cell. 2007, 6 (3): 546-554. 10.1128/EC.00330-05.PubMed CentralView ArticlePubMedGoogle Scholar
- Perkins DO, Jeffries C, Sullivan P: Expanding the 'central dogma': the regulatory role of nonprotein coding genes and implications for the genetic liability to schizophrenia. Mol Psychiatry. 2005, 10 (1): 69-78. 10.1038/sj.mp.4001577.View ArticlePubMedGoogle Scholar
- Eddy SR: Non-coding RNA genes and the modern RNA world. Nat Rev Genet. 2001, 2 (12): 919-929. 10.1038/35103511.View ArticlePubMedGoogle Scholar
- van Wolfswinkel JC, Ketting RF: The role of small non-coding RNAs in genome stability and chromatin organization. J Cell Sci. 2010, 123 (Pt 11): 1825-1839.View ArticlePubMedGoogle Scholar
- Dang Y, Yang Q, Xue Z, Liu Y: RNA interference in fungi: pathways, functions, and applications. Eukaryot Cell. 2011, 10 (9): 1148-1155. 10.1128/EC.05109-11.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (18): 5987-5994. 10.1093/nar/gkm646.PubMed CentralView ArticlePubMedGoogle Scholar
- Kankanala P, Czymmek K, Valent B: Roles for rice membrane dynamics and plasmodesmata during biotrophic invasion by the blast fungus. Plant Cell. 2007, 19 (2): 706-724. 10.1105/tpc.106.046300.PubMed CentralView ArticlePubMedGoogle Scholar
- Parker D, Beckmann M, Enot DP, Overy DP, Rios ZC, Gilbert M, Talbot N, Draper J: Rice blast infection of Brachypodium distachyon as a model system to study dynamic host/pathogen interactions. Nat Protoc. 2008, 3 (3): 435-445. 10.1038/nprot.2007.499.View ArticlePubMedGoogle Scholar
- Causton HC, Ren B, Koh SS, Harbison CT, Kanin E, Jennings EG, Lee TI, True HL, Lander ES, Young RA: Remodeling of yeast genome expression in response to environmental changes. Mol Biol Cell. 2001, 12 (2): 323-337.PubMed CentralView ArticlePubMedGoogle Scholar
- Shinozaki K, Yamaguchi-Shinozaki K: Molecular responses to dehydration and low temperature: differences and cross-talk between two stress signaling pathways. Curr Opin Plant Biol. 2000, 3 (3): 217-223.View ArticlePubMedGoogle Scholar
- Feder ME, Hofmann GE: Heat-shock proteins, molecular chaperones, and the stress response: evolutionary and ecological physiology. Annu Rev Physiol. 1999, 61: 243-282. 10.1146/annurev.physiol.61.1.243.View ArticlePubMedGoogle Scholar
- Capy P, Gasperi G, Biemont C, Bazin C: Stress and transposable elements: co-evolution or useful parasites?. Heredity (Edinb). 2000, 85 (Pt 2): 101-106.View ArticleGoogle Scholar
- Grandbastien M: Activation of plant retrotransposons under stress conditions. Trends in plant science. 1998, 3: 181-187. 10.1016/S1360-1385(98)01232-1.View ArticleGoogle Scholar
- Pecinka A, Scheid OM: Stress-Induced Chromatin Changes: a Critical View on their Heritability. Plant Cell Physiol. 2012, 10.1093/pcp/pcs044.Google Scholar
- Chiou TJ: The role of microRNAs in sensing nutrient stress. Plant Cell Environ. 2007, 30 (3): 323-332. 10.1111/j.1365-3040.2007.01643.x.View ArticlePubMedGoogle Scholar
- Sunkar R, Chinnusamy V, Zhu J, Zhu JK: Small RNAs as big players in plant abiotic stress responses and nutrient deprivation. Trends Plant Sci. 2007, 12 (7): 301-309. 10.1016/j.tplants.2007.05.001.View ArticlePubMedGoogle Scholar
- Moller T, Franch T, Udesen C, Gerdes K, Valentin-Hansen P: Spot 42 RNA mediates discoordinate expression of the E. coli galactose operon. Genes Dev. 2002, 16 (13): 1696-1706. 10.1101/gad.231702.PubMed CentralView ArticlePubMedGoogle Scholar
- Holtz J, Pasquinelli AE: Uncoupling of lin-14 mRNA and protein repression by nutrient deprivation in Caenorhabditis elegans. RNA. 2009, 15 (3): 400-405. 10.1261/rna.1258309.PubMed CentralView ArticlePubMedGoogle Scholar
- Divon HH, Fluhr R: Nutrition acquisition strategies during fungal infection of plants. FEMS Microbiol Lett. 2007, 266 (1): 65-74. 10.1111/j.1574-6968.2006.00504.x.View ArticlePubMedGoogle Scholar
- Marzluf GA: Genetic regulation of nitrogen metabolism in the fungi. Microbiol Mol Biol Rev. 1997, 61 (1): 17-32.PubMed CentralPubMedGoogle Scholar
- Froeliger EH, Carpenter BE: NUT1, a major nitrogen regulatory gene in Magnaporthe grisea, is dispensable for pathogenicity. Mol Gen Genet. 1996, 251 (6): 647-656.PubMedGoogle Scholar
- Lau G, Hamer JE: Regulatory Genes Controlling MPG1 Expression and Pathogenicity in the Rice Blast Fungus Magnaporthe grisea. Plant Cell. 1996, 8 (5): 771-781.PubMed CentralView ArticlePubMedGoogle Scholar
- Hsieh LC, Lin SI, Shih AC, Chen JW, Lin WY, Tseng CY, Li WH, Chiou TJ: Uncovering small RNA-mediated responses to phosphate deficiency in Arabidopsis by deep sequencing. Plant Physiol. 2009, 151 (4): 2120-2132. 10.1104/pp.109.147280.PubMed CentralView ArticlePubMedGoogle Scholar
- Yan Y, Zhang Y, Yang K, Sun Z, Fu Y, Chen X, Fang R: Small RNAs from MITE-derived stem-loop precursors regulate abscisic acid signaling and abiotic stress responses in rice. Plant J. 2011, 65 (5): 820-828. 10.1111/j.1365-313X.2010.04467.x.View ArticlePubMedGoogle Scholar
- Ikeda K, Nakayashiki H, Takagi M, Tosa Y, Mayama S: Heat shock, copper sulfate and oxidative stress activate the retrotransposon MAGGY resident in the plant pathogenic fungus Magnaporthe grisea. Mol Genet Genomics. 2001, 266 (2): 318-325. 10.1007/s004380100560.View ArticlePubMedGoogle Scholar
- Kim S, Ahn IP, Lee YH: Analysis of genes expressed during rice-Magnaporthe grisea interactions. Mol Plant Microbe Interact. 2001, 14 (11): 1340-1346. 10.1094/MPMI.2001.14.11.1340.View ArticlePubMedGoogle Scholar
- Keon J, Antoniw J, Rudd J, Skinner W, Hargreaves J, Hammond-Kosack K: Analysis of expressed sequence tags from the wheat leaf blotch pathogen Mycosphaerella graminicola(anamorph Septoria tritici). Fungal Genet Biol. 2005, 42 (5): 376-389. 10.1016/j.fgb.2004.12.005.View ArticlePubMedGoogle Scholar
- Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8 (4): 272-285.View ArticlePubMedGoogle Scholar
- Grandbastien MA, Audeon C, Bonnivard E, Casacuberta JM, Chalhoub B, Costa AP, Le QH, Melayah D, Petit M, Poncet C, Tam SM, Van Sluys MA, Mhiri C: Stress activation and genomic impact of Tnt1 retrotransposons in Solanaceae. Cytogenet Genome Res. 2005, 110 (1–4): 229-241.View ArticlePubMedGoogle Scholar
- Kimura M, Yamaguchi I: Convergent transcription units and their promoters at both ends of pot2, an inverted repeat transposon from the rice blast fungus. J Biochem. 1998, 124 (2): 268-273. 10.1093/oxfordjournals.jbchem.a022106.View ArticlePubMedGoogle Scholar
- Servant G, Pinson B, Tchalikian-Cosson A, Coulpier F, Lemoine S, Pennetier C, Bridier-Nahmias A, Todeschini AL, Fayol H, Daignan-Fornier B, Lesage P: Tye7 regulates yeast Ty1 retrotransposon sense and antisense transcription in response to adenylic nucleotides stress. Nucleic Acids Res. 2012, 40 (12): 5271-5282. 10.1093/nar/gks166.PubMed CentralView ArticlePubMedGoogle Scholar
- Gowda M, Venu RC, Li H, Jantasuriyarat C, Chen S, Bellizzi M, Pampanwar V, Kim H, Dean RA, Stahlberg E, Wing R, Soderlund C, Wang GL: Magnaporthe grisea infection triggers RNA variation and antisense transcript expression in rice. Plant Physiol. 2007, 144 (1): 524-533. 10.1104/pp.107.095653.PubMed CentralView ArticlePubMedGoogle Scholar
- Tisserant E, Da Silva C, Kohler A, Morin E, Wincker P, Martin F: Deep RNA sequencing improved the structural annotation of the Tuber melanosporum transcriptome. New Phytol. 2011, 189 (3): 883-891. 10.1111/j.1469-8137.2010.03597.x.View ArticlePubMedGoogle Scholar
- Ohm RA, de Jong JF, Lugones LG, Aerts A, Kothe E, Stajich JE, de Vries RP, Record E, Levasseur A, Baker SE, Bartholomew KA, Coutinho PM, Erdmann S, Fowler TJ, Gathman AC, Lombard V, Henrissat B, Knabe N, Kues U, Lilly WW, Lindquist E, Lucas S, Magnuson JK, Piumi F, Raudaskoski M, Salamov A, Schmutz J, Schwarze FW, vanKuyk PA, Horton JS, Grigoriev IV, Wosten HA: Genome sequence of the model mushroom Schizophyllum commune. Nat Biotechnol. 2010, 28 (9): 957-963. 10.1038/nbt.1643.View ArticlePubMedGoogle Scholar
- Drinnenberg IA, Fink GR, Bartel DP: Compatibility with killer explains the rise of RNAi-deficient fungi. Science. 2011, 333 (6049): 1592-10.1126/science.1209575.PubMed CentralView ArticlePubMedGoogle Scholar
- Ataide SF, Rogers TE, Ibba M: The CCA anticodon specifies separate functions inside and outside translation in Bacillus cereus. RNA Biol. 2009, 6 (4): 479-487. 10.4161/rna.6.4.9332.PubMed CentralView ArticlePubMedGoogle Scholar
- Tam OH, Aravin AA, Stein P, Girard A, Murchison EP, Cheloufi S, Hodges E, Anger M, Sachidanandam R, Schultz RM, Hannon GJ: Pseudogene-derived small interfering RNAs regulate gene expression in mouse oocytes. Nature. 2008, 453 (7194): 534-538. 10.1038/nature06904.PubMed CentralView ArticlePubMedGoogle Scholar
- Irvine DV, Zaratiegui M, Tolia NH, Goto DB, Chitwood DH, Vaughn MW, Joshua-Tor L, Martienssen RA: Argonaute slicing is required for heterochromatic silencing and spreading. Science. 2006, 313 (5790): 1134-1137. 10.1126/science.1128813.View ArticlePubMedGoogle Scholar
- Steiner FA, Hoogstrate SW, Okihara KL, Thijssen KL, Ketting RF, Plasterk RH, Sijen T: Structural features of small RNA precursors determine Argonaute loading in Caenorhabditis elegans. Nat Struct Mol Biol. 2007, 14 (10): 927-933. 10.1038/nsmb1308.View ArticlePubMedGoogle Scholar
- Tomari Y, Du T, Zamore PD: Sorting of Drosophila small silencing RNAs. Cell. 2007, 130 (2): 299-308. 10.1016/j.cell.2007.05.057.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhao T, Li G, Mi S, Li S, Hannon GJ, Wang XJ, Qi Y: A complex system of small RNAs in the unicellular green alga Chlamydomonas reinhardtii. Genes Dev. 2007, 21 (10): 1190-1203. 10.1101/gad.1543507.PubMed CentralView ArticlePubMedGoogle Scholar
- Jaskiewicz L, Filipowicz W: Role of Dicer in posttranscriptional RNA silencing. Curr Top Microbiol Immunol. 2008, 320: 77-97. 10.1007/978-3-540-75157-1_4.PubMedGoogle Scholar
- Kadotani N, Murata T, Quoc NB, Adachi Y, Nakayashiki H: Transcriptional control and protein specialization have roles in the functional diversification of two dicer-like proteins inMagnaporthe oryzae. Genetics. 2008, 180 (2): 1245-1249. 10.1534/genetics.108.093922.PubMed CentralView ArticlePubMedGoogle Scholar
- Bohnert HU, Fudal I, Dioh W, Tharreau D, Notteghem JL, Lebrun MH: A putative polyketide synthase/peptide synthetase from Magnaporthe grisea signals pathogen attack to resistant rice. Plant Cell. 2004, 16 (9): 2499-2513. 10.1105/tpc.104.022715.PubMed CentralView 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 small interfering RNA accumulation. J Biol Chem. 2004, 279 (43): 44467-44474. 10.1074/jbc.M408259200.View ArticlePubMedGoogle Scholar
- Lu C, Meyers BC, Green PJ: Construction of small RNA cDNA libraries for deep sequencing. Methods. 2007, 43 (2): 110-117. 10.1016/j.ymeth.2007.05.002.View ArticlePubMedGoogle Scholar
- Rice P, Longden I, Bleasby A: EMBOSS: the European Molecular Biology Open Software Suite. Trends Genet. 2000, 16 (6): 276-277. 10.1016/S0168-9525(00)02024-2.View ArticlePubMedGoogle Scholar
- Reid RJ, Sunjevaric I, Keddache M, Rothstein R: Efficient PCR-based gene disruption in Saccharomyces strains using intergenic primers. Yeast. 2002, 19 (4): 319-328. 10.1002/yea.817.View ArticlePubMedGoogle Scholar
- Sweigard JA, Chumley FG, Valent B: Disruption of a Magnaporthe grisea cutinase gene. Mol Gen Genet. 1992, 232 (2): 183-190.PubMedGoogle Scholar
- Sweigard JA, Chumley F, Carroll AM, Farrall L, Valent B: A series of vectors for fungal transformation. Fungal Genetics Newsletter. 1997, 44: 52-53.Google Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.