- Research article
- Open Access
Physiological stressors and invasive plant infections alter the small RNA transcriptome of the rice blast fungus, Magnaporthe oryzae
BMC Genomics volume 14, Article number: 326 (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
We analyzed small RNA (sRNA) size profiles of libraries constructed from M. oryzae mycelial tissue, and from susceptible rice leaves infected with M. oryzae at two time-points. For the former, M. oryzae was subjected to the in vitro stresses of nitrogen starvation (NS), carbon starvation (CS), minimal media (MM), and oxidative-inducing conditions via paraquat (PQ). Mycelia were also grown on a rich, complete media (CM) for comparative purposes. Growth conditions for these in vitro stresses mimicked those of Mathioni et al. , and hereafter will be referred to as the “mycelial” libraries. In order to compare the in vitro stress sRNA profiles to small RNAs generated during in planta growth, we also isolated fungal-induced rice lesions at 72 and 96 hpi. These will be hereafter referred to as the “in planta” libraries. Small RNA libraries were generated and quality checked before performing deep sequencing using Illumina’s Genome Analyzer II (GAII) system. This resulted in 58,511,280 and 28,941,539 sequences from the mycelial and in planta libraries, respectively (Table 1). The sequences were aligned to the M. oryzae genome assembly version 6, downloaded from the Broad Institute (http://www.broadinstitute.org/). Total genome-matched reads (mapping to one or more sites in the M. oryzae genome) and distinct genome-matched reads were calculated. There were 34,652,258 and 2,414,112 total genome-matched reads, and 927,501 and 54,394 distinct genome-matched reads for mycelial and in planta libraries, respectively, (Table 1). For the in planta libraries, M. oryzae sequences were parsed out from rice sequences by alignment to the fungal genome. An informatics workflow describing this process is presented in Additional file 1: Figure S1, along with an example of the output from the pipeline as applied to a mycelial library.
The size distribution of total genome-matched reads for mycelial libraries were bimodal with peaks at 23/24 nucleotides (nt) and 26/27 nt, while distinct genome-matched reads showed a normal distribution curve with the peak at 23/24 nt (Figure 1). The abundance of 24 nt reads represented 16.77% of the total population. In planta libraries also had a similar pattern of abundance, but the major peak was at 26 nt followed by 24 nt for the total genome-matched reads and 22/23 nt for the distinct genome-matched reads (Figure 2). Almost 12% of the sRNAs are represented by the 24 nt class. Collectively, all mycelial libraries showed the majority of sRNA peaks at 24 followed by 26 nt, whereas in planta libraries showed the majority of sRNAs peaking at 26 followed by 24 nt.
Cluster analysis of mycelial and in plantalibraries
We performed a cluster analysis on each M. oryzae supercontig, based upon the proximity-based algorithm of Lu et al.  for a comparative analysis of the mycelial and in planta small RNA libraries. This clustering algorithm allowed us to determine the abundance of sRNAs and their distribution across the M. oryzae genome by dividing the genome into a series of 500 base pair (bp) windows (clusters) and calculating the abundance of sRNAs for each cluster. The clusters are defined by fixed coordinates, which allows comparative analyses of multiple libraries and the identification of differentially regulated clusters (DRCs) in the genome. In our study, DRCs represent regions of the genome with a sRNA abundance that differs by 10-fold or more in pairwise comparisons between the libraries. All mycelial libraries were compared to CM, which resulted in a total of 477 DRCs. We examined the total number of sRNA clusters for each library and identified clusters that map to genic (which includes UTRs, introns and exons) and intergenic regions, various forms of repeat regions and transposable element regions. Out of the 477 total sRNA clusters, 251 were from genic regions, 213 were from intergenic regions and 20 were from repeat regions and transposable element regions (Table 2). The CM vs. CS comparison revealed four DRCs, indicating carbon starvation had a minimal impact on the accumulation of sRNAs. NS vs. CM and PQ vs. CM had the greatest number of intergenic and genic-associated clusters. The genic and intergenic-associated clusters were 52 and 96 for NS vs. CM, and 179 and 69 for PQ vs. CM respectively. The NS vs. CM comparison resulted in the highest number (5) of sRNAs mapping to retrotransposons and transposons, versus other comparisons. Results from the in planta libraries showed that the majority of clusters associated with genes and intergenic regions, as opposed to other genomic regions, and this was similar to what was seen in the CM vs. PQ comparison. Overall, the mycelial libraries showed more small RNAs mapping to transposons and retrotransposons, compared to the in planta libraries; for the latter libraries, with the exception of two clusters mapping to retrotransposons, almost everything else associated with genic regions (Table 2).
sRNA associate with different genomic regions, depending upon the fungal stress
In total, taking data from all the libraries together, the majority of sRNAs were found to be associated with intergenic regions (54%), followed by repeat-containing regions (41%; Table 3). On the other hand, comparisons across libraries resulted in more DRCs associated with genes (Table 3). Although more total sRNAs mapped to intergenic regions compared to protein-coding regions, protein-coding regions contained a greater number of distinct sRNAs (Table 2). This difference may be explained by the presence of repeats which have not yet been annotated, and which were not identified by the RepeatMasker repeat identification program (; http://www.repeatmasker.org). Though approximately 50% of the genome codes for proteins and 7% is occupied by repeats , our sRNA data suggests only 4% of the total sRNAs sequenced matched to protein coding regions.
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
Nucleotide preference in the 5’ position is the determining factor in plants for binding of sRNAs to specific AGO proteins . While fungi do not appear to share this exact characteristic, preference for particular 5’ nucleotides have been observed [12, 21, 36, 37]. In M. oryzae, analysis of different libraries for their 5’ nt preference revealed that sRNAs starting with U were abundant in both mycelial and in planta libraries, and there was a suppression for C in mycelial libraries and for G in in planta libraries (Figure 3). The 5’ U preference was maintained in individual mycelial libraries – CM, MM, NS and PQ, but the CS library contained more sRNAs starting with G, closely followed by U. The suppression of C at the 5’ position was also maintained in most of the mycelial libraries except for MM and PQ (data not shown). The 5’ end nucleotide preference for the LTR class was U, and for LINE/Tad1 was both A and U (data not shown). The general pattern of abundant U and less G for sRNAs from in planta libraries was reflected in the LMg0 and LMg96, while LMg72 library had abundant A and less G (data not shown).
Correlation between transcriptional and post-transcriptional regulation under stress conditions
We hypothesized that transcriptional down-regulation of genes may correlate with a plethora of sRNAs accumulating either around that gene in the intergenic spaces, or within the gene itself. In order to gain a genome-wide perspective on how transcriptional and post-transcriptional gene regulation might be related, we compared our sRNA data to transcriptional genome data generated from M. oryzae under the same “stress” conditions as Mathioni et al.  in their microarray study. From this dataset, we identified 98 genes associated with clusters that showed differential accumulation of sRNAs (≥ 10-fold compared to CM control), 28 and 9 of which showed a negative correlation with microarray data under PQ and NS conditions, respectively (Additional file 6: Figure S6; Additional file 7: Table S1; Additional file 8: Table S2; Figure 4). Further analysis of four of these genes showed that they have sRNAs originating from exons or introns in either the sense or antisense orientations (Table 4). Of these, only MGG_16711.7 (earlier annotated as MGG_05869.6) coding for a hypothetical protein containing an importin-beta N-terminal domain, had many sRNAs originating from its intron.
For the in planta libraries, we were most interested in sRNAs that were exclusively present at 72 hpi and/or 96 hpi (Additional file 9: Figure S7), indicating that these fungal sRNAs might be involved in invasive growth. In order to determine whether any of these sRNAs were negatively correlated with genes from the microarray study as described above, we first selected genes associated with sRNA clusters from the comparisons of LMg0 vs. LMg72, LMg0 vs. LMg96 and LMg72 vs. LMg96 that showed a ≥ 10-fold difference. Next, their corresponding loci from the CM vs. LMg72 data were compared to the microarray study . This resulted in 12 genes that were both down-regulated in the microarray, associated with large quantities of sRNAs (Table 5). For each gene, abundant sRNAs present in LMg72 and/or LMg96 were analyzed for their position, and checked for their absence in non-inoculated rice sRNA libraries (http://mpss.udel.edu/rice_sRNA), indicating that they were originating from the fungal, and not the rice, genome. The resultant sRNAs were largely found to be associated with exons, 5’ UTR or 3’ UTRs in a sense orientation (Table 5). Interestingly, one sRNA matching to MGG_08843 coding for the magnesium transporter, ALR2, was abundant at 72 hpi and present in an antisense orientation to the 5’UTR; Mathioni et al.  demonstrated down-regulation of this gene during carbon starvation, as well as in planta growth in rice and barley. Another gene includes the well-characterized transcription factor, ACE1 (MGG_04428.7); this gene is of interest, as it is known to be tightly regulated and involved in appressorial formation .
In order to confirm whether these genes might truly be under post-transcriptional regulation, we generated four sRNA biosynthetic mutants, and utilized real-time qRT-PCR to determine the expression of seven genes identified from the above analysis, namely MGG_08843.6, MGG_01596.6, MGG_04428.6, MGG_06609.6, MGG_05869.6, MGG_01439.6 and MGG_04470.6. Mutants tested included gene replacement deletions in the following loci: Dicer 1 and Dicer 2 (MGG_01541 and MGG_12357, previously generated and characterized by Kadotani et al. 2002, and hereafter referred to as MoDcl-1 and MoDcl-2), one RNA-dependent RNA Polymerase (RdRP MGG_13453; hereafter referred to as RdRP) and a double-knockout mutant in both Dicers (see Additional file 10: Figure S8 for deletion construct details, confirmation of deletions and Southern hybridization, revealing one disruption cassette per mutant line). We hypothesized that suppression of gene expression would be lifted if the genes for generating sRNAs in the fungus were no longer present, and sRNAs were no longer made. Six of the seven genes analyzed showed 1.5 to 1.9-fold increased expression over the wild type in the MoDcl1 mutant. The MGG_01596 gene, encoding a DNA damage response protein DUN1, increased in expression in the ∆modcl1, modcl2/modcl1 and in the ∆mordrp-1 and was reduced in the ∆modcl2 mutant backgrounds. The MGG_04428 gene encoding the transcription factor ACE1 was increased ~1.7-fold in the ∆modcl1 background (Figure 5).
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.
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.
Baulcombe D: RNA silencing. Curr Biol. 2002, 12 (3): R82-4. 10.1016/S0960-9822(02)00665-6.
Tomari Y, Zamore PD: Perspective: machines for RNAi. Genes Dev. 2005, 19 (5): 517-529. 10.1101/gad.1284105.
Baulcombe D: RNA silencing in plants. Nature. 2004, 431 (7006): 356-363. 10.1038/nature02874.
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.
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.
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.
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.
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.
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.
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.
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.
Nakayashiki H: RNA silencing in fungi: mechanisms and applications. FEBS Lett. 2005, 579 (26): 5950-5957. 10.1016/j.febslet.2005.08.016.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Liu XH, Lu JP, Lin FC: Autophagy during conidiation, conidial germination and turgor generation in Magnaporthe grisea. Autophagy. 2007, 3 (5): 472-473.
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.
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.
Smit AFA, Hubley R, Green P: RepeatMasker Open-3.0. 1996-2010, http://www.repeatmasker.org,
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.
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.
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.
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.
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.
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.
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.
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.
Eddy SR: Non-coding RNA genes and the modern RNA world. Nat Rev Genet. 2001, 2 (12): 919-929. 10.1038/35103511.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Grandbastien M: Activation of plant retrotransposons under stress conditions. Trends in plant science. 1998, 3: 181-187. 10.1016/S1360-1385(98)01232-1.
Pecinka A, Scheid OM: Stress-Induced Chromatin Changes: a Critical View on their Heritability. Plant Cell Physiol. 2012, 10.1093/pcp/pcs044.
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.
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.
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.
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.
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.
Marzluf GA: Genetic regulation of nitrogen metabolism in the fungi. Microbiol Mol Biol Rev. 1997, 61 (1): 17-32.
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.
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.
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.
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.
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.
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.
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.
Slotkin RK, Martienssen R: Transposable elements and the epigenetic regulation of the genome. Nat Rev Genet. 2007, 8 (4): 272-285.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
Sweigard JA, Chumley FG, Valent B: Disruption of a Magnaporthe grisea cutinase gene. Mol Gen Genet. 1992, 232 (2): 183-190.
Sweigard JA, Chumley F, Carroll AM, Farrall L, Valent B: A series of vectors for fungal transformation. Fungal Genetics Newsletter. 1997, 44: 52-53.
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).
The authors declare that they have no competing interests.
BCM and ND designed the study and NMD and VR wrote the manuscript. VR and SAS generated the libraries. FD and JZ performed bioinformatic analyses. MR generated one mutant and SMM provided guidance on microarray comparisons. Mutant generation and molecular analyses were performed by VR SAS, VR, NMD, BCM, SMM and MR edited the manuscript. All authors read and approved the final manuscript.
Vidhyavathi Raman contributed equally to this work.
Electronic supplementary material
Additional file 1: Figure S1: Informatics pipeline used for the analysis sRNA libraries of M. oryzae different libraries. We used a pipeline with different filtering steps to analyze the raw reads obtained from Illumina. The filtering steps (left column) are summarized here and fully described in the main text. Results of each filtering step to the M. oryzae mycelial complete media (CM) library are shown in the right column, as an example. The colors of the numbers indicate either total (black) or distinct (red) sequences. (PPTX 69 KB)
Additional file 2: Figure S2: sRNAs associate with the retrotransposon class, LTR/Gypsy. Size distribution of representative loci from the major retrotransposon classes LTR/Gypsy in mycelial libraries (A); One representative locus of GYMAG1_I-int LTR/Gypsy on chromosome 6 under wild type (complete media) conditions (B); abundance of 16 different LTR/Gypsy retrotransposons on sense (open box) and antisense (filled box) strands of all mycelial libraries (C); abundance of sRNAs from LTR/Gypsy classes under different stresses in the mycelial libraries (D). (CM = complete media; CS = carbon starved; MM = minimal media; NS = nitrogen starved; PQ = paraquat). (PPTX 357 KB)
Additional file 3: Figure S3: sRNAs associate with the retrotransposon class, LINE/Tad1. Size distribution of representative LINE/Tad1 from different chromosomes in mycelial libraries (A); abundance of sRNAs from six LINE/Tad1 loci on sense (open box) and antisense (closed box) strands of all mycelial libraries (B); and abundance of sRNAs from LINE/Tad1 classes under different stresses from the mycelial libraries (C). (CM = complete media; CS = carbon starved; MM = minimal media; NS = nitrogen starved; PQ = paraquat). (PPTX 78 KB)
Additional file 4: Figure S4: Characterization of sRNAs associated with 5S rRNA and tRNA. sRNA size distribution (A and B) and abundance under different stress conditions (C and D) for 5S rRNA (left panels) and tRNA (right panels). (PPTX 75 KB)
Additional file 5: Figure S5: sRNAs associated with different genomic regions. Fraction of small RNA reads that associate with different genomic loci in mycelial libraries (A) and in in planta libraries (B). Proportion of individual repeat classes to total repeats in mycelial libraries (C) and in in planta libraries (D). (CM = complete media; CS = carbon starved; MM = minimal media; NS = nitrogen starved; PQ = paraquat; LMg0 = mock inoculated rice; LMg72 = 72 hpi; LMg96 = 96 hpi). (PPTX 74 KB)
Additional file 6: Figure S6: Screenshot representing sRNAs associated with MGG_01439, the putative inorganic phosphate transporter, and its neighboring intergenic regions under CS (carbon starved) and NS (nitrogen starved) conditions. (PPTX 196 KB)
Additional file 7: Table S1: List of genes showing negative correlation between sRNA and microarray data under PQ condition. (XLSX 11 KB)
Additional file 8: Table S2: List of genes showing negative correlation between sRNA and microarray data under NS condition. (XLSX 9 KB)
Additional file 9: Figure S7: Venn diagram depicting the clustering of sRNAs from LMg0, LMg72 and LMg96 libraries. The center gray area likely represents small RNAs from rice that were also found in the inoculated samples. (LMg0 = mock inoculated rice; LMg72 = 72 hpi; LMg96 = 96 hpi). (PPTX 89 KB)
Additional file 10: Figure S8: Generation of knock-out mutants in sRNA pathways. Schematic representation of the chromosomal locus of wild type (WT; top panel) and the deletion mutant MoDcl1 (middle panel), as well as MoDcl2 in the MoDcl1 mutant background (lower panel) (A); Confirmation of the deletion of MoDcl1 and the double mutant, using flanking, gene-specific and hygromycin/BAR-specific primers. KO = knockout; ECT = ectopic (B); Southern analysis of ∆modcl1, ∆modcl2 and ∆modcl2/modcl1 mutants. 1kb ladder (1), 70–15 (2), ∆modcl2 KO (3), and Ectopic (5) mutants genomic DNA digested with HindIII, and 70–15 (7), ∆modcl1 KO (8) and Ectopic (11) mutants genomic DNA digested with EcoRI and HindIII (C); ∆modcl2/modcl1 KO1 (2), KO2 (3), Ectopic (4) mutants genomic DNA digested with EcoRI (D). All knock-out mutants were performed in the same way, with the same confirmations. (PPTX 226 KB)
About this article
Cite this article
Raman, V., Simon, S.A., Romag, A. et al. Physiological stressors and invasive plant infections alter the small RNA transcriptome of the rice blast fungus, Magnaporthe oryzae. BMC Genomics 14, 326 (2013). https://doi.org/10.1186/1471-2164-14-326
- Small RNA
- Magnaporthe oryzae
- Rice blast fungus