- Research article
- Open Access
Comparative analysis among the small RNA populations of source, sink and conductive tissues in two different plant-virus pathosystems
© Herranz et al.; licensee BioMed Central. 2015
- Received: 6 August 2014
- Accepted: 6 February 2015
- Published: 22 February 2015
In plants, RNA silencing plays a fundamental role as defence mechanism against viruses. During last years deep-sequencing technology has allowed to analyze the sRNA profile of a large variety of virus-infected tissues. Nevertheless, the majority of these studies have been restricted to a unique tissue and no comparative analysis between phloem and source/sink tissues has been conducted. In the present work, we compared the sRNA populations of source, sink and conductive (phloem) tissues in two different plant virus pathosystems. We chose two cucurbit species infected with two viruses very different in genome organization and replication strategy; Melon necrotic spot virus (MNSV) and Prunus necrotic ringspot virus (PNRSV).
Our findings showed, in both systems, an increase of the 21-nt total sRNAs together with a decrease of those with a size of 24-nt in all the infected tissues, except for the phloem where the ratio of 21/24-nt sRNA species remained constant. Comparing the vsRNAs, both PNRSV- and MNSV-infected plants share the same vsRNA size distribution in all the analyzed tissues. Similar accumulation levels of sense and antisense vsRNAs were observed in both systems except for roots that showed a prevalence of (+) vsRNAs in both pathosystems. Additionally, the presence of overrepresented discrete sites along the viral genome, hot spots, were identified and validated by stem-loop RT-PCR. Despite that in PNRSV-infected plants the presence of vsRNAs was scarce both viruses modulated the host sRNA profile.
We compare for the first time the sRNA profile of four different tissues, including source, sink and conductive (phloem) tissues, in two plant-virus pathosystems. Our results indicate that antiviral silencing machinery in melon and cucumber acts mainly through DCL4. Upon infection, the total sRNA pattern in phloem remains unchanged in contrast to the rest of the analyzed tissues indicating a certain tissue-tropism to this polulation. Independently of the accumulation level of the vsRNAs both viruses were able to modulate the host sRNA pattern.
- Small RNAs
In most eukaryotes RNA silencing is an essential mechanism which controls gene expression in a wide range of biological processes. In addition, RNA silencing plays an important role as defence against viruses in plants, invertebrates and possibly also in mammals [1-4]. Double-stranded RNA (dsRNA) is the factor triggering the process and upon infection the host silencing machinery is induced by dsRNAs generated during viral replication or secondary structures derived from the pathogen genome [5-7]. These viral dsRNAs are processed by RNase III enzymes belonging to the Dicer family (Dicer-like or DCL in plants). In Arabidopsis thaliana, four classes of DCL enzymes (DCL1 to DCL4) have been reported, DCL1 mainly contributes to the production of miRNAs  and DCL4, DCL2 and DCL3 are involved in the processing of viral genomes yielding viral-sRNAs (vsRNAs) of 21-, 22- and 24-nt, respectively [9-11]. In infected plants, hierarchical roles of DCL4 and DCL2 have been established; 21-nt vsRNAs are the most abundant class followed by 22-nt vsRNAs [12-16]. Subsequent loading of vsRNAs to Argonaute (AGO)-containing complexes or RISCs directs the degradation of both genomic and subgenomic viral RNAs [17,18]. In plants, this antiviral silencing response can be amplified by cellular RNA-dependent RNA polymerases (RDRs) which use as templates cleaved viral RNAs to generate new dsRNA substrates for Dicer enzymes .
In plants, a series of grafting experiments demonstrated that RNA silencing has a systemic nature and that the long-distance transmission of the putative sequence-specific signaling molecule must be transported through the phloem .
Little is known about the nature of the systemic signal that mediates the process and how it is transmitted from the site of initiation to the rest of the plant. However, the detection of sRNAs in the phloem of different plant species [21-24] as well as the identification of phloem RNA binding proteins with capacity to bind and translocate sRNAs [21,25] have led to the assumption that the systemic signal must be a ribonucleoprotein complex. In addition, Tournier et al.  demonstrated that the silencing signal is transported from source to sink tissues following the direction of phloem flow. Whereas some lines of evidence have indicated a role for sRNAs with different sizes in long-distance transmission of RNA silencing [27-29], other studies have suggested that the signal could be either a sRNA precursor or produced from dsRNA by a DCL-independent mechanism .
During last years deep-sequencing analyses in different plant/virus systems have allowed to study the vsRNA profile of a high variety of infected tissues [7,31-42]. Nevertheless, with very few exceptions [32,35] this type of studies have been restricted to a unique tissue and no comparisons among phloem and source/sink tissues have been made. Therefore, the specific aims of this work were: 1) to compare the populations of viral-derived sRNAs among source, sink and conductive (phloem) tissues and 2) to compare these populations in two different pathosystems sharing cucurbit species as a host. To do that, we selected two plant hosts infected with viruses very different in genome organization and replication strategy. Thus, inoculated cotyledon as a source, root and symptomatic primary leaf as recipient tissues and phloem sap from melon and cucumber plants infected with Melon necrotic spot virus (MNSV) and Prunus necrotic ringspot virus (PNRSV), respectively, were analyzed.
Melon necrotic ringspot virus (MNSV) is a plant virus belonging to the genus Carmovirus within the family Tombusviridae [43,44] which is present in cucurbit crops worldwide. The MNSV genome is a single-stranded RNA molecule of 4.3 kb with positive polarity encoding at least five different proteins [44,45]. The open reading frame (ORF) at the 5′ end terminates in an amber codon yielding two proteins involved in replication, p29 and p89. Cell-to-cell viral movement is supported by two proteins, p7A and p7B, encoded by two small centrally-located ORFs [46-49]. The ORF at the 3′end encodes the coat protein (CP) p42 which is also involved in systemic transport of the virus and is a symptom determinant . By contrast, Prunus necrotic ringspot virus (PNRSV) is a positive single-stranded RNA virus, member of the genus Ilarvirus in the family Bromoviridae and with a tripartite genome. RNAs 1 and 2 encode the replicase subunits P1 and P2. RNA 3 is bicistronic and contains the coding sequences of the putative movement protein (MP) gene and the CP gene. CP synthesis occurs via a subgenomic monocistronic mRNA (RNA 4) [51,52].
Our findings showed, in both systems, an increase of the 21-nt total sRNAs together with a decrease of those with a size of 24-nt in all the infected tissues, except for the phloem where the ratio of 21/24-nt sRNAs remained constant. Comparing the vsRNAs, both PNRSV and MNSV infected plants share the same vsRNA size distribution in all the analyzed tissues. Although in PNRSV-infected plants the percentage of vsRNAs was one hundred times less than that observed for the MNSV- infected plants, both viruses were able to modulate the host sRNA profile. In addition, similar accumulation levels of sense and antisense vsRNAs were observed in both systems except for roots that showed a prevalence of (+) vsRNAs in both pathosystems. The biological significance of these results is discussed.
Viruses, host plants and tissues used for construction of sRNA libraries
Number of reads
To confirm that our virus-infected samples were indeed infected and to compare virus accumulation in melon and cucumber, respectively, we performed Northern blot analysis of total RNA from cotyledon, leaf, root and phloem samples of MNSV and PNRSV-infected plants using a virus-specific probe (Additional file 2: Figure S1). Analysis of serial dilutions of each infected sample by dot-blot (Additional file 3: Figure S2) and quantification of the signal using the Java image processing programe ImageJ indicated similar virus concentrations in leaves and the highest difference in cotyledons with a ratio of 7 considering the higher dilution of the virus (Additional file 4: Table S2).
The total sRNA size profile differs in healthy and MNSV infected melon plants
The sRNA profile of phloem sap was considerably different from the profile of the other tissues. As it is shown in Figure 1 the majority of sRNAs in healthy and infected phloem have a length of 24-nt (35%). The band corresponding to this size is clearly visible in the PAGE-SN (Figure 2C, lanes 13 and 14). Moreover, in both the infected and healthy plants, a higher amount of small RNAs with a 16-nt size was observed (Figure 1 and Figure 2C, lanes 13 and 14). In contrast to the other tissues, the 21-nt sRNA population did not significantly increase upon infection in phloem sap. In addition, in this tissue, the 21/24 ratio is maintained regardless of whether the plant is infected or not (Additional file 5: Figure S3A).
The sRNA size profile differs in healthy and PNRSV-infected cucumber plants
Analysis of the sRNA profile in healthy and PNRSV-infected cucumber plants revealed that the total sRNAs size distribution pattern was very similar to that observed for melon plants. Compared with healthy tissues, infected cotyledon and leaf samples exhibited increased levels of 21-nt sRNAs (from around 27% to 45% in both tissues) together with a lower amount of 24-nt (from 25% to 12% in cotyledon and 35% to 14% in leaf) species (Figure 1C and 1D). It is noticeable that in infected root- although there is also an increase-the amount of sRNAs with a size of 21-nt was smaller (Figure 1D) and species with a size between 16 and 20-nt were much more represented in both healthy and infected root. This is also clearly visible in the PAGE-SN (Figure 2B, lanes 9 and 12) where a band pattern appears below the 21-nt size. As seen in cotyledons and leaves, virus infection produced a three-fold reduction of the sRNAs with a size of 24-nt (approximately 13% in mock-inoculated vs 4% in PNRSV-infected roots).
The sRNA profile of cucumber phloem was different from the profile of the other tissues. While phloem of healthy cucumber plants contained mostly sRNAs with a size of 24 (25%), 23 (21%) and 21 (15%) nucleotides (Figure 1C), the ratio 21/24-nt sRNAs remained constant in infected phloem tissue (Figure 1D and Additional file 5: Figure S3B). Although the signal was very weak, we could confirm the presence of the 24-nt sRNAs by polyacrylamide gel analysis (Figure 2C, lanes 15 and 16), the rest of the sizes were hardly visible. As for melon phloem tissue, a more intense double band at the bottom of the gel appeared as well.
The percentage of vsRNAs in PNRSV infected cucumber plants is below 1%
vsRNAs are mainly 21 nucleotides in length
In cucumber the results were similar, approximately 50% of the vsRNAs from PNRSV infected leaves, cotyledons, roots and phloem had a length of 21 nucleotides, however, they were hardly detected in infected cotyledon, leaf and root tissues (Figure 2B, lower panel). This low vsRNA signal confirms the low percentage of these species found in PNRSV-infected cucumber plants by deep-sequencing analysis.
Based on these two observations, we decided to test whether or not the concentration of vsRNAs in infected cucumber plants was under the detection limit analyzing serial dilutions of a 21-nt transcript with known concentration by hybridization with the isotopic riboprobe. With this approach we were able to detect up to 15 ng of the specific sRNA (Figure 2D). Considering the percentages of vsRNAs in all infected cucumber tissues, calculations were made with cotyledon which has the highest concentration representing 0.8% of the total RNA. 4 μg of sample containing thus 32 ng of vsRNAs were loaded into the gel and therefore, the amount of vsRNAs in the cucumber samples was close to the detection limit and clearly under this limit for the phloem samples. This result explains the weak or null signals obtained in the Northern blot analysis and validates the low percentage of vsRNAs observed in infected cucumber plants in comparison with the infected melon plants. Although it is reasonable to assume that accumulation levels of vsRNAs could correlate with the accumulation levels of viral RNAs this seems not to be the case since, as stated above, both pathosystems accumulated similar viral RNAs levels (Additional file 2: Figure S1).
The antiviral silencing machinery acts mainly through DCL4 in this plant genus
Knowing that vsRNAs in melon and cucumber plants were mainly 21-nt in size, we decided to analyze expression levels of the different DCL transcripts by real-time qRT-PCR. A genome-wide analysis using protein BLAST search of the corresponding databases (http://www.phytozome.net and http://www.melonomics.net) identified the four homologues of A. thaliana DCL1 to 4 in both C. sativus and C. melo genome (melon accession numbers and cucumber locus names are provided in Methods section). Additional file 6: Figure S4 and Additional file 7: Figure S5 show the differential expression of DCL1 to 4 transcripts in different tissues of both melon and cucumber, respectively. As the expression of these enzymes has not been described in phloem to date, this tissue was excluded from the analysis. According to the results, DCL4 seemed to be the most important antiviral silencing enzyme as its expression level clearly increased in all the infected tissues and in both pathosystems. Regarding the other DCLs, the patterns that compared both hosts differed. Upon MNSV infection, we observed a significant increase in the expression level of the DCL3 transcript in infected cotyledon (Additional file 6: Figure S4). However in the PNRSV-infected cucumber plants, infection modulated the expression of DCL3 in leaves and roots, and also of DCL2 in cotyledon which led to higher levels of transcripts (Additional file 7: Figure S5).
Mapping of vsRNAs against the viral genome suggests tissue-specific differences in the origin of the viral sRNAs
In A. thaliana, the 5′ terminal nucleotide of the sRNAs influences their selective loading into specific AGOs. Thus, investigating the 5′-terminal nucleotide of vsRNAs can give insight into vsRNA function. Figures 5E and 6E shows the percentage of vsRNAs of 21 to 24 nt with 5′ terminal Uracil, Cytosine, Guanine, or Adenine in the different tissues. We found a similar percentage for each nucleotide independently of the size, tissue and plant-virus system (Figures 5E and 6E). This result suggests the involvement of different Ago complexes in all analyzed tissues, including phloem.
Presence of selected vsRNAs “hot-spots”was validated by stem-loop RT-PCR
As described above, both MNSV and PNRSV vsRNAs mapped along the whole corresponding virus genome. However, PNRSV vsRNA distribution was more heterogeneous with a higher accumulation along the RNA 3 (Figure 6). Additionally and in both systems, vsRNAs accumulated specifically at discrete sites along the virus genome. Three hot-spots of vsRNA accumulation located at the N-terminal half of the p42 MNSV ORF were especially prevalent. Two of them were of sense polarity and mainly composed of a unique 21-nt length vsRNA. The third MNSV hot-spot was of antisense polarity and composed of vsRNAs of 21- and 22-nt. According to their 5′ nucleotide position on the genome, these hot-spots were named as vsRNA 3078, 3099 and 3160, respectively. The relative abundance of vsRNA 3078 and 3160 was comparable in all the tissues (0.1-0.12% for vsRNA 3078 and 0.27-0.35% for vsRNA 3160) except in phloem sap were both, vsRNA 3078 and 3160, were selectively enriched becoming the most abundant vsRNAs (0.6% and 0.59% respectively). For vsRNA 3099, the highest relative abundance was observed in cotyledon (0.37%) and root (0.32%) decreasing in phloem (0.25%) and leaf (0.1%) (Figure 5). For PNRSV, two noteworthy hot-spots were observed, one antisense vsRNA of 21-nt at the 3′ UTR region of the RNA3 (vsRNA 103) which was mainly abundant in cotyledon (2.3%) and an interesting sense vsRNA of 25-nt located at the C-terminus of the coat protein (vsRNA 1691) with a high number of reads in phloem (7.3%) (Figure 6).
RNA silencing controls countless biological processes including, in the case of plants, defence against biotic and abiotic stresses. When a plant is infected by a virus, the silencing machinery has to combine efforts to guarantee the correct development of the plant on the one hand and an appropriate defence against the pathogen on the other. This involves, among other aspects, a silencing signal with nucleotide-sequence specificity that spreads over long distances through the phloem tissue. Next generation sequencing has allowed a considerable progress in identifying and characterizing the hallmark molecules in this process, the vsRNAs. However, in most of these studies, the characterization of the vsRNAs has been usually restricted to a particular tissue. In addition, when phloem tissue was analyzed, its vsRNA population has not been compared, with very few exceptions, with sink and/or source tissues. In this work we deep-sequenced the sRNA population of four tissues; cotyledon as source, symptomatic-primary leaf and root as sink tissues and phloem as conductive tissue. In addition, we compared two different plant-virus systems, MNSV and PNRSV, infecting melon and cucumber plants, respectively.
Sequence analysis of the data revealed a remarkable difference between the total number of vsRNA reads in both systems. Consistent with our data, Donaire et al.  found high variability in vsRNA reads when comparing a variety of plant-virus interactions. The authors attributed this variability to differences in virus accumulation and replication, efficiency of the RNA machinery and to the mode of action of the viral silencing suppressor (VSR). Similar to that reported here, the same percentage of vsRNAs (around 57%) was obtained in MNSV infected cotyledons from melon plants. Curiously, the profile of vsRNAs of the Cucumber mosaic virus (CMV), a member belonging to the same family than PNRSV, showed also a lower percentage of virus-derived sRNAs (around 14%) in Arabidopsis thaliana.
Our study did not reveal a significant correlation between virus accumulation and number of vsRNAs. Thus, other factors like differences in virus counter defense strategies in the respective host plants likely accounts for the differences in vsRNA accumulation in the plant-virus systems tested. Little is known about the MNSV and PNRSV VSRs. For MNSV, the two viral proteins, p42 and p7B have been described to have VSR activity . They were classified as weak because they delayed but did not prevent PTGS in transient expression experiments on GFP-transgenic plants. The high number of MNSV-derived sRNAs indicated that MNSV RNAs serve as highly effective substrates for DCL enzymes. The fact that large amounts of vsRNAs coexist with high titers of MNSV genome accumulation indicates that MNSV VSRs could act downstream of sRNA production in the PTGS pathway. One possibility could be that the MNSV VSRs avoid RISC assembly through vsRNA sequestering. Alternatively, MNSV replication and spreading could occur at rates that outcompete the capacity of the plant silencing machinery for its degradation. Further experiments to determine how these two proteins work at the molecular level and whether or not they interfere with the silencing signal spread are necessary. For PNRSV, no viral protein acting as a VSR has been described so far. Our finding that only very low levels of PNRSV sRNAs accumulate upon infection indicates that this virus encodes a VSR acting upstream of vsRNA production or uses another non-conventional mechanism allowing the virus to escape silencing. Examples of such mechanisms have been previously described; Cauliflower mosaic virus (CMV), for instance, uses decoy RNAs to evade silencing , Red clover mosaic virus (RCNMV) suppresses RNA silencing recruiting DCL1 proteins or their homologs to the replication complexes  and highly structured viral RNAs of human adenovirus sequester Dicer and suppress antiviral defence . Like PNRSV, Magnaporthe oryzae virus 2 (MoV2), a dsRNA mycovirus, accumulates very low amounts of vsRNAs (0.5%). None of the viral proteins exhibit RNA silencing suppressor activity but MoV2 seems to evade silencing by replicating within viral particles . Similarly, the PNRSV genome might be protected from degradation through interaction with proteins or early encapsidation. Interestingly, PNRSV, like other ilarviruses and Alfalfa mosaic virus (AMV), require a specific and efficient interaction of their CPs with the viral genome for replication and translation . It would be interesting to know if the low accumulation level of vsRNAs observed from PNRSV is a general rule for the rest of ilarviruses and AMV.
Importantly, despite the fact that PNRSV infection hardly produced any vsRNAs, total sRNA profiles MNSV and PNRSV infected plants showed a clear increase of the 21-nt species in all the tissues being less evident in phloem. A sharp decline of 24-nt sRNA relative abundance was also observed in all the tissues except for the phloem sap where the amount remains similar and even higher than the 21-nt species. This indicates that virus infection modulates dicer activities independent of whether the viral RNA is targeted by silencing or not.
Similar to our findings, a 2.5-3 fold reduction of 24-nt sRNAs relative abundance compared to healthy plants was detected in Citrus tristeza virus (CTV)-infected Mexican lime and sweet orange but not in sour orange . Interestingly, in the two former citrus species, CTV invades adjacent phloem-associated cells at distal regions of the plant producing small infection clusters. In less-susceptible sour orange species, the virus is able to exit sieve elements but cannot spread to adjacent cells resulting in isolated single infected cells and less than 1/10 of Mexican lime virus concentration . A similar effect on 24-nt sRNAs concentration occurs in the synergistic sweet potato virus disease (SPVD) . In SPVD, a phloem-limited crinivirus (Sweet potato chlorotic stunt virus, SPCSV) increases 600-fold the titer of an unrelated potyvirus (Sweet potato feathery mottle virus, SPFMV) in non-phloem tissue resulting in highly symptomatic plants. It was hypothesized that SPCSV interferes with the systemic phloem dependent signaling required for SPFMV resistance . Only mobile 24-nt sRNAs have been associated with RNA-dependent DNA methylation (RdDM) in recipient cells that may trigger transcriptional gene silencing (TGS) . Whether local reduction of the 24-nt abundance allows high titers of viral accumulation by affecting fundamental mechanisms related with gene expression of viral resistance in distal parts needs to be clarified.
From our results and previous studies, it appears that all size classes of sRNA are mobile. However, when we compared libraries from healthy and infected samples, the most remarkable new observation was that whereas for the rest of the tissues the ratio between the 21/24-nt sRNAs increases, it is maintained in phloem (Additional file 5: Figure S3). In infected samples, both types of sRNAs have to be loaded into sieve tubes against a concentration gradient suggesting the existence of a mechanism that actively selects the sRNAs to be transported. Analysis of the sRNA composition in phloem has been previously carried out in different systems [21,23,32] but only Yoo et al.  identified the small RNA species in phloem of Cucumber yellows closterovirus (CuYV) infected pumpkin plants. In contrast to what we report here, the authors observed a distribution pattern clearly different from the healthy one. However, this apparent contradiction needs to be carefully considered since CuYV vsRNAs, mostly of 20 and 21-nt, represented 57% of total sRNAs. In our study, due to the low percentage of vsRNAs in infected phloem (10% in melon and close to 0% in cucumber), the bulk of the sRNA population in phloem belongs to the endogenous sRNA group. Consequently the size distribution analysis of the host genome sRNAs showed the same pattern than for the total sRNAs (Additional file 8 Figure S6 and Additional file 9: Figure S7). To better understand host responses to infection and possibly to gain insight into transported host sRNAs it will be interesting to see which host sRNAs exhibit modulated abundance in the different tissues upon infection. Additionally, the fact that despite an overall increase of 21-nt sRNAs and decrease of 24-nt sRNAs upon infection in leaves, cotyledons and roots the ratio between 21-nt to 24-nt sRNAs remains constant in phloem tissue could suggest that 24-nt sRNAs are the main transported sRNA species in these plants. In addition, as the percentage of 24-nt species in phloem is similar upon infection (Additional file 8: Figure S6 and Additional file 9: Figure S7), a qualitative change of these small molecules could alert the distal parts of the plant. Some lines of evidence have indicated a role for sRNAs with different sizes in long-distance transmission of RNA silencing [27-29]. Alternatively, viruses may have found a way to selectively inhibit transport of 21-nt sRNA species.
The antiviral silencing machinery in Arabidopsis acts mainly through DCL4, which produces 21-nt small RNAs . As mentioned above, and consistently with this, vsRNAs in melon and cucumber plants were mainly 21-nt in size, which correlated with a clear increase in the expression level of the DCL4 transcript in all the tissues, and indicates that dicer functions and dicing patterns are conserved between species. Accumulation of the 21-nt long vsRNAs in both virus/plant systems indicates that neither of the two unrelated viruses used in this study impaired the DCL4 function. Thus DCL4 seemed to be the most important antiviral enzyme in this plant genus. In all tissues, our results also indicated that whereas MNSV vsRNAs derived equally from whole genomic RNA, PNRSV vsRNAs mapped at a higher percentage to RNA3. It is noteworthy that all the tissues in both systems showed equivalent amounts of sense and antisense vsRNAs, except for roots, where the number of vsRNAs with a positive polarity increased. The bias to the genomic sense strand in roots could support a model by which the secondary structures within viral single strand RNA contribute to vsRNA generation. But why does this happen only in roots? A previous work by our group revealed high levels of MNSV accumulation in infected melon roots . This differential tropism was also observed recently for citrus species infected with Citrus tristeza virus (CTV) . In MNSV-infected plants, the concentration of dsRNAs, as a measure of replication, was 5-fold lower in roots than in cotyledons, despite the amount of virus being comparable in both tissues. This observation could consequently imply a greater accumulation of viral genomic RNA, which could serve as a substrate for DCL cleavage, and would support the observed increase in sense vsRNAs in our study. According to Gosalvez-Bernal et al. , one would imagine roots to be like a prelude to the virus being loaded into the internal phloem; in other words, a bottleneck where the replication of the virus would decrease, and would accumulate while waiting to enter the vascular system in order to spread. Alternatively, Andika et al.  provided evidence that the RNA silencing mechanism is less effective in roots than in leaves, which could explain the substantial viral RNA accumulation in this tissue. However, and as stated above, the accumulation levels of vsRNAs in MNSV-infected roots were similar to the other tissues, while other mechanisms, which differed from RNA silencing, might account for this differential tissue tropism. Further experiments that investigate this increase in sense vsRNAs in MNSV- and PNRSV-infected roots, as well as the molecular mechanism that underlies reduced viral replication in this tissue, will provide new insights into this phenomenon. According to previous work, sRNA loading into specific AGO complexes is conditioned by the 5′ terminal nucleotide [67-69]. In the present work, none of the vsRNA 5′ends showed any prevalence. High variability in this aspect [36,70] further evidences the great complexity of this mechanism and the wide range of AGO complexes, presumably with different roles, which might be involved in the process.
Finally, our comparative sRNA profiling analysis has revealed that MNSV and PNRSV infections both lead to the production of mainly 21-nt vsRNAs, which is consistent with DCL4 being the main antiviral silencing component in plants. The fact that PNRSV infection hardly led to the production of vsRNAs indicates that this virus has evolved mechanisms to avoid being targeted by the host silencing machinery. However, both viruses strongly altered the total sRNA profile, which reveals that virus infection modulates the host gene expression. In future analyses, it will be interesting to see whether these changes are host defense responses or whether the virus actively uses silencing to down-regulate defense. The different pattern of sRNAs found in phloem tissue is probably one of the most interesting observations to have emerged from our comprehensive sRNA analysis. Given that viral or host sRNAs are transported over long distances, it will be interesting to analyze the identity of the sRNA population in phloem compared to other tissues, and to also analyze the target gene expression in sink tissues.
We have compared the sRNA profile of four different tissues, including source, sink and conductive (phloem) tissues, in two different pathosystems for the first time. Our results indicate that the antiviral silencing machinery in melon and cucumber acts mainly through DCL4. One of the most interesting results to have emerged from our sRNA analysis is that the total sRNA pattern in phloem remains unchanged upon infection, unlike the other analyzed tissues. In addition, and independently of the accumulation level of vsRNAs, both viruses were able to modulate the host sRNA pattern.
Plant material and RNA extraction
Cotyledons of both Cucumis melo cv. Galia and Cucumis sativus cv. Supermarketer plants were inoculated with purified virions and crude virus-containing extracts from infected plants respectively, 9 days after germination. As a control, mock-inoculated plants were grown in parallel and under the same greenhouse conditions (16 hours at 24°C and 8 hours at 18°C in daylight conditions). Inoculated-cotyledons, root, phloem and primary leaf from mock-inoculated and symptomatic infected plants were harvested 15 days after inoculation (15 dpi).
To collect the exudates, the first petiole and the stem internode between the first and second leaf were cut with a sterile razor blade and with the help of a pipet the phloem was collected in an Eppendorf-tube containing 1 ml of TRIzol.
Total RNA extraction was identical for all the samples and carried out by using conventional TRIzol (Sigma-Aldrich) protocol according to the manufacturer’s instructions. The quality of the samples was verified by Northern-blot analysis using specific digoxigenin-labelled riboprobes for each virus and evaluating both A260/230 and A260/280 absorption ratios (Nanodrop- Thermo Scientific).
Small RNA PEG fractionation
To enrich the samples in sRNAs, 3 volumes of polyethylene glycol (MW 8000) and NaCl to final concentrations of 5% and 500 mM respectively were added. After incubation on ice for 30 minutes, the samples were centrifuged at 15000 rpm for 15 min at 4°C. The upper phase was transferred to a new tube and 4 volumes of phenol were added. Then the tubes were vortexed thoroughly and centrifuged during 2 minutes at maximum speed. Afterwards the upper phase was precipitated and the pellet containing sRNAs was resuspended in 50 μl of water.
1 μg of each fractionated sample was separated in a 17% polyacrylamide/7 M urea gel for the subsequent silver staining. For this purpose, the gel was soaked for 45 minutes in ethanol/H2O/acetic acid (ratio 50:40:1) followed by a second 45 minutes incubation in the same solution with a ratio of 10:89:1. Thereafter the gel was dipped in 150 ml of H2O with 0.3 g AgNO3 for 45 minutes. Then, the gel was rinsed 3 times with sterile water and developed in a formaldehyde solution (0.5 M KOH, 1.2 ml 37% formaldehyde, in 150 ml H2O). The reaction was stopped by washing with sterile water.
Isotopic Northern-blot analysis
In order to detect vsRNAs, melon (2 μg) and cucumber (4 μg) fractionated RNA samples were separated in a 17% polyacrylamide/ 7 M urea gel in 0.5X TBE, transferred onto a positively-charged nylon membrane using the Bio-Rad Trans-Blot Cell and covalently UV cross-linked to the membrane (700 × 100 μJ/cm2). Detection of the vsRNAs was carried out using specific riboprobes labeled with (α-32P) ATP. After overnight hybridization, the membranes were washed twice with 2X SSC plus 0.1% SDS for 10 min at room temperature, and once with 0.1X SSC plus 0.1% SDS at 55°C for 15 min, and examined with a bioimage analyzer (Fujifilm FLA-5100).
sRNA sequence processing
Production and sequencing of the libraries were carried out by the biotechnology company GenoScreen (http://www.genoscreen.com). Small RNAs were fractionated from total RNA by acrylamide gel purification. Single strand ligation of 3′ and 5′ adaptors was done before a second acrylamide gel purification. Reverse transcription and PCR amplification were performed to generate the cDNA colonies template library. To verify quality and ensure that there was no contamination, libraries were then titrated by a 1x50bp run on HiSeq2000.
A total of 16 sRNA libraries were sequenced on an Illumina genomic DNA analyzer in 1 HiSeq 2000® channel. Adapter trimming and cleaning of the reads was carried out by the same company. Most of the inserts found had a length ranging from 18 to 26 bases.
Each sRNA library was aligned with the corresponding plant (Melonomics: melon_genome_pseudomolecules_V3.5.fasta http://www.melonomics.net Cucumber genome database: ACHR00000000.1 http://www.cucumber.genomics.org.cn) and known virus genomes (MNSV Al: DQ339157 or PNRSV NcM1.NctSp.mur1: AJ306818) by the Bioinformatics Service at the IBMCP (http://www.ibmcp.upv.es). Count of matching sequence reads were normalized to the total number of reads after last filtering step and given in parts per million
Hot-spot validation by Stem-loop RT-PCR
We used the protocol described by Varkonyi-Gasic, et al.  with some modifications. For the pulse reverse transcription reaction we mixed 0.5 μM of appropriate stem-loop RT primer (Additional file 10: Table S3), 1 mM dNTPs and 100 ng of the corresponding RNA in a volume of 14.9 μl and incubated at 65°C during 5 minutes and then on ice. Then, 200 U of RevertAid Reverse Transcriptase (Thermo Scientific; http://www.thermoscientific.com/en/home.html) and 40 U RNase Inhibitor were added in a final volume of 20 μl and pulsed RT was performed as follow: 30 min at 16°C, pulsed RT of 60 cycles at 30°C for 30 s, 42°C for 30 s and 50°C for 1 s. Then samples were incubated at 85°C for 5 min. PCR was carried out with the following components: 0.2 μM of each primer (Additional file 10: Table S3), 0.2% PVP, 0.2 mM dNTP’s, 0.5 mM MgCl2 and 2 U GoTaq (Promega; http://www.promega.es). The PCR reaction was performed as follows: 2 min at 94°C, 25 cycles 15 s at 94°C and 1 min at 60°C.
Identification of AtDCL1-4 genes homologs in both Cucumis sativus and Cucumis melo genomes
Arabidopsis thaliana DCL proteins (DCL1; At1g01040, DCL2; At3g03300, DCL3; At3g43920 and DCL4; At5g 20320) were blasted against the corresponding database (http://www.phytozome.net and http://www.melonomics.net) to identify the four homologues in C. sativus and C. melo. Accession numbers: Cucsa.260100, Cucsa.356600, Cucsa.055950, Cucsa.350750 (C. sativus) and MELO3C005929P1/30P1, MELO3C010042P1, MELO3C011495P1, MELO3C010254P1 (C. melo).
Real-time qRT-PCR analysis
To analyze the differential expression of DCL transcripts, qRT-PCR was performed with 100 ng of total RNA using one step SYBR PrimeScript RT-PCR Kit II (Takara) following the manufacturer’s recommendations in an Applied Biosystems 7500 Real-Time PCR System. Each biological replicate was assayed in triplicate. Gene-specific oligonucleotide primers were designed using Primer Express® version 3.0 software (Applied Biosystems). Primer information is shown in Additional file 11: Table S4. Expression levels for target genes were normalized to Elongation Factor 1-alpha (EF1 α) and phosphatase 2A regulatory subunit (PP2A) and fold expression changes compared to the healthy controls calculated using the ΔΔthreshold cycle (Ct) method.
We thank Dr A. Niehl for critical reading and helpful comments on the manuscript. This work was funded by a supporting program for the research from the Universidad Politécnica de Valencia (PAID-05-10), a grant BIO2011-25018 from the Spanish granting agency Direccion General de Investigacion Cientifica and the PROMETEO program 2011/003 from the Generalitat Valenciana. MCH is the recipient of a contract from JAE-DOC program of the CSIC, JAN is the recipient of a postdoctoral contract from the Ministerio de Educación y Ciencia of Spain.
- Pumplin N, Voinnet O. RNA silencing suppression by plant pathogens: defence, counter-defence and counter-counter-defence. Nat Rev Microbiol. 2013;11(11):745–60.View ArticlePubMedGoogle Scholar
- Brodersen P, Voinnet O. The diversity of RNA silencing pathways in plants. Trends Genet. 2006;22(5):268–80.View ArticlePubMedGoogle Scholar
- Ghildiyal M, Zamore PD. Small silencing RNAs: an expanding universe. Nat Rev Genet. 2009;10(2):94–108.View ArticlePubMed CentralPubMedGoogle Scholar
- Ciaudo C, Jay F, Okamoto I, Chen CJ, Sarazin A, Servant N, et al. RNAi-dependent and independent control of LINE1 accumulation and mobility in mouse embryonic stem cells. PLoS Genet. 2013;9(11):e1003791.View ArticlePubMed CentralPubMedGoogle Scholar
- Ding SW, Voinnet O. Antiviral immunity directed by small RNAs. Cell. 2007;130(3):413–26.View ArticlePubMed CentralPubMedGoogle Scholar
- Szittya G, Moxon S, Pantaleo V, Toth G, Rusholme Pilcher RL, Moulton V, et al. Structural and functional analysis of viral siRNAs. PLoS Pathog. 2010;6(4):e1000838.View ArticlePubMed CentralPubMedGoogle Scholar
- Donaire L, Wang Y, Gonzalez-Ibeas D, Mayer KF, Aranda MA, Llave C. Deep-sequencing of plant viral small RNAs reveals effective and widespread targeting of viral genomes. Virology. 2009;392(2):203–14.View ArticlePubMedGoogle Scholar
- Voinnet O. Origin, biogenesis, and activity of plant microRNAs. Cell. 2009;136(4):669–87.View ArticlePubMedGoogle Scholar
- Liu Q, Feng Y, Zhu Z. Dicer-like (DCL) proteins in plants. Funct Integr Genomics. 2009;9(3):277–86.View ArticlePubMedGoogle Scholar
- Henderson IR, Zhang X, Lu C, Johnson L, Meyers BC, Green PJ, et al. Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning. Nat Genet. 2006;38(6):721–5.View ArticlePubMedGoogle Scholar
- Margis R, Fusaro AF, Smith NA, Curtin SJ, Watson JM, Finnegan EJ, et al. The evolution and diversification of Dicers in plants. FEBS Lett. 2006;580(10):2442–50.View ArticlePubMedGoogle Scholar
- Deleris A, Gallego-Bartolome J, Bao J, Kasschau KD, Carrington JC, Voinnet O. Hierarchical action and inhibition of plant Dicer-like proteins in antiviral defense. Science. 2006;313(5783):68–71.View ArticlePubMedGoogle Scholar
- Blevins T, Rajeswaran R, Shivaprasad PV, Beknazariants D, Si-Ammour A, Park HS, et al. Four plant Dicers mediate viral small RNA biogenesis and DNA virus induced silencing. Nucleic Acids Res. 2006;34(21):6233–46.View ArticlePubMed CentralPubMedGoogle Scholar
- Bouche N, Lauressergues D, Gasciolli V, Vaucheret H. An antagonistic function for Arabidopsis DCL2 in development and a new function for DCL4 in generating viral siRNAs. EMBO J. 2006;25(14):3347–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Moissiard G, Voinnet O. RNA silencing of host transcripts by cauliflower mosaic virus requires coordinated action of the four Arabidopsis Dicer-like proteins. Proc Natl Acad Sci U S A. 2006;103(51):19593–8.View ArticlePubMed CentralPubMedGoogle Scholar
- Qu F, Ye X, Morris TJ. Arabidopsis DRB4, AGO1, AGO7, and RDR6 participate in a DCL4-initiated antiviral RNA silencing pathway negatively regulated by DCL1. Proc Natl Acad Sci U S A. 2008;105(38):14732–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Vaucheret H. Plant ARGONAUTES. Trends Plant Sci. 2008;13(7):350–8.View ArticlePubMedGoogle Scholar
- Hutvagner G, Simard MJ. Argonaute proteins: key players in RNA silencing. Nat Rev Mol Cell Biol. 2008;9(1):22–32.View ArticlePubMedGoogle Scholar
- Voinnet O. Use, tolerance and avoidance of amplified RNA silencing by plants. Trends Plant Sci. 2008;13(7):317–28.View ArticlePubMedGoogle Scholar
- Palauqui JC, Elmayan T, Pollien JM, Vaucheret H. Systemic acquired silencing: transgene-specific post-transcriptional silencing is transmitted by grafting from silenced stocks to non-silenced scions. EMBO J. 1997;16(15):4738–45.View ArticlePubMed CentralPubMedGoogle Scholar
- Yoo BC, Kragler F, Varkonyi-Gasic E, Haywood V, Archer-Evans S, Lee YM, et al. A systemic small RNA signaling system in plants. Plant Cell. 2004;16(8):1979–2000.View ArticlePubMed CentralPubMedGoogle Scholar
- Buhtz A, Pieritz J, Springer F, Kehr J. Phloem small RNAs, nutrient stress responses, and systemic mobility. BMC Plant Biol. 2010;10:64.View ArticlePubMed CentralPubMedGoogle Scholar
- Buhtz A, Springer F, Chappell L, Baulcombe DC, Kehr J. Identification and characterization of small RNAs from the phloem of Brassica napus. Plant J. 2008;53(5):739–49.View ArticlePubMedGoogle Scholar
- Rodriguez-Medina C, Atkins CA, Mann AJ, Jordan ME, Smith PM. Macromolecular composition of phloem exudate from white lupin (Lupinus albus L.). BMC Plant Biol. 2011;11:36.View ArticlePubMed CentralPubMedGoogle Scholar
- Pallas V, Gomez G. Phloem RNA-binding proteins as potential components of the long-distance RNA transport system. Frontiers in Plant Science. 2013;4:130.PubMed CentralPubMedGoogle Scholar
- Tournier B, Tabler M, Kalantidis K. Phloem flow strongly influences the systemic spread of silencing in GFP Nicotiana benthamiana plants. Plant J. 2006;47(3):383–94.View ArticlePubMedGoogle Scholar
- Hamilton A, Voinnet O, Chappell L, Baulcombe D. Two classes of short interfering RNA in RNA silencing. EMBO J. 2002;21(17):4671–9.View ArticlePubMed CentralPubMedGoogle Scholar
- Voinnet O. MicroRNA and autophagy--C. elegans joins the crew. EMBO Rep. 2013;14(6):485–7.View ArticlePubMed CentralPubMedGoogle Scholar
- Dunoyer P, Schott G, Himber C, Meyer D, Takeda A, Carrington JC, et al. Small RNA duplexes function as mobile silencing signals between plant cells. Science. 2010;328(5980):912–6.View ArticlePubMedGoogle Scholar
- Brosnan CA, Mitter N, Christie M, Smith NA, Waterhouse PM, Carroll BJ. Nuclear gene silencing directs reception of long-distance mRNA silencing in Arabidopsis. Proc Natl Acad Sci U S A. 2007;104(37):14741–6.View ArticlePubMed CentralPubMedGoogle Scholar
- Silva TF, Romanel EA, Andrade RR, Farinelli L, Osteras M, Deluen C, et al. Profile of small interfering RNAs from cotton plants infected with the polerovirus Cotton leafroll dwarf virus. BMC Mol Biol. 2011;12:40.View ArticlePubMed CentralPubMedGoogle Scholar
- Martinez G, Donaire L, Llave C, Pallas V, Gomez G. High-throughput sequencing of Hop stunt viroid-derived small RNAs from cucumber leaves and phloem. Mol Plant Pathol. 2010;11(3):347–59.View ArticlePubMedGoogle Scholar
- Donaire L, Barajas D, Martinez-Garcia B, Martinez-Priego L, Pagan I, Llave C. Structural and genetic requirements for the biogenesis of tobacco rattle virus-derived small interfering RNAs. J Virol. 2008;82(11):5167–77.View ArticlePubMed CentralPubMedGoogle Scholar
- Qi X, Bao FS, Xie Z. Small RNA deep sequencing reveals role for Arabidopsis thaliana RNA-dependent RNA polymerases in viral siRNA biogenesis. PLoS One. 2009;4(3):e4971.View ArticlePubMed CentralPubMedGoogle Scholar
- Pantaleo V, Saldarelli P, Miozzi L, Giampetruzzi A, Gisel A, Moxon S, et al. Deep sequencing analysis of viral short RNAs from an infected Pinot Noir grapevine. Virology. 2010;408(1):49–56.View ArticlePubMedGoogle Scholar
- Lin KY, Cheng CP, Chang BC, Wang WC, Huang YW, Lee YS, et al. Global analyses of small interfering RNAs derived from Bamboo mosaic virus and its associated satellite RNAs in different plants. PLoS One. 2010;5(8):e11928.View ArticlePubMed CentralPubMedGoogle Scholar
- Navarro B, Pantaleo V, Gisel A, Moxon S, Dalmay T, Bisztray G, et al. Deep sequencing of viroid-derived small RNAs from grapevine provides new insights on the role of RNA silencing in plant-viroid interaction. PLoS One. 2009;4(11):e7686.View ArticlePubMed CentralPubMedGoogle Scholar
- Martin R, Arenas C, Daros JA, Covarrubias A, Reyes JL, Chua NH. Characterization of small RNAs derived from Citrus exocortis viroid (CEVd) in infected tomato plants. Virology. 2007;367(1):135–46.View ArticlePubMedGoogle Scholar
- St-Pierre P, Hassen IF, Thompson D, Perreault JP. Characterization of the siRNAs associated with peach latent mosaic viroid infection. Virology. 2009;383(2):178–82.View ArticlePubMedGoogle Scholar
- Di Serio F, Gisel A, Navarro B, Delgado S, de Alba AE M, Donvito G, et al. Deep sequencing of the small RNAs derived from two symptomatic variants of a chloroplastic viroid: implications for their genesis and for pathogenesis. PLoS One. 2009;4(10):e7539.View ArticlePubMed CentralPubMedGoogle Scholar
- Li R, Gao S, Hernandez AG, Wechter WP, Fei Z, Ling KS. Deep sequencing of small RNAs in tomato for virus and viroid identification and strain differentiation. PLoS One. 2012;7(5):e37127.View ArticlePubMed CentralPubMedGoogle Scholar
- Hu Q, Hollunder J, Niehl A, Korner CJ, Gereige D, Windels D, et al. Specific impact of tobamovirus infection on the Arabidopsis small RNA profile. PLoS One. 2011;6(5):e19549.View ArticlePubMed CentralPubMedGoogle Scholar
- Hibi T, Furuki I. Melon Necrotic Spot Virus. In: CMI: AAB Descriptions of Plants Viruses N° 302. Kew, UK: Commonwealth Mycological Institute; 1985.Google Scholar
- Riviere CJ, Rochon DM. Nucleotide sequence and genomic organization of melon necrotic spot virus. J Gen Virol. 1990;71(Pt 9):1887–96.View ArticlePubMedGoogle Scholar
- Diaz JA, Nieto C, Moriones E, Truniger V, Aranda MA. Molecular characterization of a Melon necrotic spot virus strain that overcomes the resistance in melon and nonhost plants. Mol Plant Microbe Interact. 2004;17(6):668–75.View ArticlePubMedGoogle Scholar
- Navarro JA, Genoves A, Climent J, Sauri A, Martinez-Gil L, Mingarro I, et al. RNA-binding properties and membrane insertion of Melon necrotic spot virus (MNSV) double gene block movement proteins. Virology. 2006;356(1–2):57–67.View ArticlePubMedGoogle Scholar
- Genoves A, Navarro JA, Pallas V. A self-interacting carmovirus movement protein plays a role in binding of viral RNA during the cell-to-cell movement and shows an actin cytoskeleton dependent location in cell periphery. Virology. 2009;395(1):133–42.View ArticlePubMedGoogle Scholar
- Genoves A, Navarro JA, Pallas V. The Intra- and intercellular movement of Melon necrotic spot virus (MNSV) depends on an active secretory pathway. Mol Plant Microbe Interact. 2010;23(3):263–72.View ArticlePubMedGoogle Scholar
- Serra-Soriano M, Pallas V, Navarro JA. A model for transport of a viral membrane protein through the early secretory pathway: minimal sequence and endoplasmic reticulum lateral mobility requirements. Plant J. 2014;77(6):863–79.View ArticlePubMedGoogle Scholar
- Genoves A, Navarro JA, Pallas V. Functional analysis of the five melon necrotic spot virus genome-encoded proteins. J Gen Virol. 2006;87(Pt 8):2371–80.View ArticlePubMedGoogle Scholar
- Pallas V, Aparicio F, Herranz MC, Amari K, Sanchez-Pina MA, Myrta A, et al. Ilarviruses of Prunus spp.: a continued concern for fruit trees. Phytopathology. 2012;102(12):1108–20.View ArticlePubMedGoogle Scholar
- Pallas V, Aparicio F, Herranz MC, Sanchez-Navarro JA, Scott SW. The molecular biology of ilarviruses. Adv Virus Res. 2013;87:139–81.View ArticlePubMedGoogle Scholar
- Varkonyi-Gasic E, Wu R, Wood M, Walton EF, Hellens RP. Protocol: a highly sensitive RT-PCR method for detection and quantification of microRNAs. Plant Methods. 2007;3:12.View ArticlePubMed CentralPubMedGoogle Scholar
- Blevins T, Rajeswaran R, Aregger M, Borah BK, Schepetilnikov M, Baerlocher L, et al. Massive production of small RNAs from a non-coding region of Cauliflower mosaic virus in plant defense and viral counter-defense. Nucleic Acids Res. 2011;39(12):5003–14.View ArticlePubMed CentralPubMedGoogle Scholar
- Takeda A, Tsukuda M, Mizumoto H, Okamoto K, Kaido M, Mise K, et al. A plant RNA virus suppresses RNA silencing through viral RNA replication. EMBO J. 2005;24(17):3147–57.View ArticlePubMed CentralPubMedGoogle Scholar
- Andersson MG, Haasnoot PC, Xu N, Berenjian S, Berkhout B, Akusjarvi G. Suppression of RNA interference by adenovirus virus-associated RNA. J Virol. 2005;79(15):9556–65.View ArticlePubMed CentralPubMedGoogle Scholar
- Himeno M, Maejima K, Komatsu K, Ozeki J, Hashimoto M, Kagiwada S, et al. Significantly low level of small RNA accumulation derived from an encapsidated mycovirus with dsRNA genome. Virology. 2010;396(1):69–75.View ArticlePubMedGoogle Scholar
- Aparicio F, Vilar M, Perez-Paya E, Pallas V. The coat protein of prunus necrotic ringspot virus specifically binds to and regulates the conformation of its genomic RNA. Virology. 2003;313(1):213–23.View ArticlePubMedGoogle Scholar
- Ruiz-Ruiz S, Navarro B, Gisel A, Pena L, Navarro L, Moreno P, et al. Citrus tristeza virus infection induces the accumulation of viral small RNAs (21-24-nt) mapping preferentially at the 3′-terminal region of the genomic RNA and affects the host small RNA profile. Plant Mol Biol. 2011;75(6):607–19.View ArticlePubMedGoogle Scholar
- Folimonova SY, Folimonov AS, Tatineni S, Dawson WO. Citrus tristeza virus: survival at the edge of the movement continuum. J Virol. 2008;82(13):6546–56.View ArticlePubMed CentralPubMedGoogle Scholar
- Kreuze JF, Perez A, Untiveros M, Quispe D, Fuentes S, Barker I, et al. Complete viral genome sequence and discovery of novel viruses by deep sequencing of small RNAs: a generic method for diagnosis, discovery and sequencing of viruses. Virology. 2009;388(1):1–7.View ArticlePubMedGoogle Scholar
- Karyeija RF, Kreuze JF, Gibson RW, Valkonen JP. Synergistic interactions of a potyvirus and a phloem-limited crinivirus in sweet potato plants. Virology. 2000;269(1):26–36.View ArticlePubMedGoogle Scholar
- Melnyk CW, Molnar A, Bassett A, Baulcombe DC. Mobile 24 nt small RNAs direct transcriptional gene silencing in the root meristems of Arabidopsis thaliana. Curr Biol. 2011;21(19):1678–83.View ArticlePubMedGoogle Scholar
- Gosalvez-Bernal B, Genoves A, Navarro JA, Pallas V, Sanchez-Pina MA. Distribution and pathway for phloem-dependent movement of Melon necrotic spot virus in melon plants. Mol Plant Pathol. 2008;9(4):447–61.View ArticlePubMedGoogle Scholar
- Harper SJ, Cowell SJ, Robertson CJ, Dawson WO. Differential tropism in roots and shoots infected by Citrus tristeza virus. Virology. 2014;460–461:91–9.View ArticlePubMedGoogle Scholar
- Andika IB, Kondo H, Tamada T. Evidence that RNA silencing-mediated resistance to beet necrotic yellow vein virus is less effective in roots than in leaves. Mol Plant Microbe Interact. 2005;18(3):194–204.View ArticlePubMedGoogle Scholar
- Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, et al. Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5' terminal nucleotide. Cell. 2008;133(1):116–27.View ArticlePubMed CentralPubMedGoogle Scholar
- Takeda A, Iwasaki S, Watanabe T, Utsumi M, Watanabe Y. The mechanism selecting the guide strand from small RNA duplexes is different among argonaute proteins. Plant Cell Physiol. 2008;49(4):493–500.View ArticlePubMedGoogle Scholar
- Wu L, Zhang Q, Zhou H, Ni F, Wu X, Qi Y. Rice MicroRNA effector complexes and targets. Plant Cell. 2009;21(11):3421–35.View ArticlePubMed CentralPubMedGoogle Scholar
- Xu Y, Huang L, Fu S, Wu J, Zhou X. Population diversity of rice stripe virus-derived siRNAs in three different hosts and RNAi-based antiviral immunity in Laodelphgax striatellus. PLoS One. 2012;7(9):e46238.View ArticlePubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.