In-depth sequencing of the siRNAs associated with peach latent mosaic viroid infection
© Bolduc et al; licensee BioMed Central Ltd. 2010
Received: 23 October 2009
Accepted: 16 February 2010
Published: 16 February 2010
It has been observed that following viroid infection, there is an accumulation of viroid-derived siRNAs in infected plants. Some experimental results suggest that these small RNAs may be produced by the plant defense system to protect it from infection, indicating that viroids can elicit the RNA-silencing pathways. The objective of this study is to identify in the peach latent mosaic viroid (PLMVd), a model RNA genome, the regions that are most susceptible to RNA interference machinery.
The RNA isolated from an infected tree have been used to sequence in parallel viroid species and small non-coding RNA species. Specifically, PLMVd RNAs were amplified, cloned and sequenced according to a conventional approach, while small non-coding RNAs were determined by high-throughput sequencing. The first led to the typing of 18 novel PLMVd variants. The second provided a library of small RNAs including 880 000 sequences corresponding to PLMVd-derived siRNAs, which makes up 11.2% of the sequences of the infected library. These siRNAs contain mainly 21-22 nucleotide RNAs and are equivalently distributed between the plus and the minus polarities of the viroid. They cover the complete viroid genome, although the amount varies depending on the regions. These regions do not necessarily correlate with the double-stranded requirement to be a substrate for Dicer-like enzymes. We noted that some sequences encompass the hammerhead self-cleavage site, indicating that the circular conformers could be processed by the RNA-silencing machinery. Finally, a bias in the relative abundance of the nature of the 5' nucleotides was observed (A, U >> G, C).
The approach used provided us a quantitative representation of the PLMVd-derived siRNAs retrieved from infected peach trees. These siRNAs account for a relatively large proportion of the small non-coding RNAs. Surprisingly, the siRNAs from some regions of the PLMVd genome appear over-represented, although these regions are not necessarily forming sufficiently long double-stranded structures to satisfy Dicer-like criteria for substrate specificity. Importantly, this large library of siRNAs gave several hints as to the components of the involved silencing machinery.
Viroids are small (240-460 nucleotides, nt), single-stranded circular RNA pathogens that self-replicate in the plants causing important economical losses in agricultural industries . Their genomes do not encode any protein, so they rely exclusively on host proteins. They are divided into two families: the Pospiviroidae and the Avsunviroidae. The members of the first family, for which the typical representative is the potato spindle tuber viroid (PSTVd), share a conserved central region, a nuclear localization and an asymmetric rolling circle replication that appears to involve host enzymes. Members of the Avsunviroidae, including the peach latent mosaic viroid (PLMVd), which is the viroid model in the present study, have been shown to accumulate in the chloroplasts . The latter viroids possess hammerhead self-cleaving RNA motifs and replicate according to a symmetric rolling circle mechanism.
Plants accumulate a multitude of small non-coding RNAs (i.e. 18-25 nucleotides in length). The actions of most of these small RNAs result in gene repression through what is called post-transcriptional gene silencing (PTGS) to protect the plants from biotic stresses, like viroid or virus infection. Besides their role in response to pathogens, these small RNAs are also implicated in development and genome integrity . Based on the detection of small RNAs, it has been established that viroids induce PTGS [4–8]. In addition, it has been demonstrated that a small hairpin derived from PSTVd was sufficient to cause symptoms in its natural host, providing a potential explanation for the pathogenicity mechanism . Moreover the compact viroid structure makes them resistant to RISC degradation  at least in the case of PSTVd. However, contradictory reports have shown that, under some circumstances, viroid-derived small RNAs could protect the plants from viroid infection [11, 12], making the association between RNA silencing and pathogenesis questionable.
In the case of PLMVd, some aspects of the RNA-silencing mechanism, including how it might trigger Dicer-like enzyme activity, were studied using a cell-free wheat germ assay . In addition, we recently reported the results of a sequencing effort directed towards the siRNAs associated with PLMVd infection in peach leaves . Specifically, a fraction of small, non-coding RNAs (ncRNAs) was isolated, cloned and sequenced according to the classical procedure for ncRNAs, and any PLMVd-derived siRNAs that were found were then analyzed. This work yielded several interesting observations, including the fact that PLMVd strands of both polarities may enter the peach RNA silencing pathways. These observations should be considered with caution as they are based on the analysis of only 60 PLMVd-derived siRNAs, which out of 284 ncRNAs detected in infected plant leaves, corresponds to only 3.9× coverages of the viroid genome. Moreover, since RNA silencing in plant includes a phenomenon of amplification by an RNA-dependent RNA polymerase , it was surprising to note that some regions of the genome were not covered by siRNAs.
High-throughput sequencing (HTS) is a method that provides a relatively unbiased and direct digital readout of cDNA sequences generated from an RNA sample [15, 16]. The unparalleled ability of HTS to yield quantitative information on transcript abundance should help to give a more representative picture of the PLMVd-derived siRNA population than what is available right now. With the goal of obtaining a clearer view of the PLMVd-derived siRNAs, the sequencing of the small RNAs associated with this infection was repeated using HTS. Concurrently, the viroid species present in RNA samples used were also sequenced. The parallel sequencing of the viroid and derived siRNAs retrieved within a unique RNA sample should remove potential bias.
Results and Discussion
Analysis of the HTS data for the RNA obtained from both the healthy and the PLMVd-infected peach leaf RNA samples.
Unfiltered total reads
38 035 624
10 502 809
12 476 493
7 862 905
Reads homologous to rRNAs
Reads homologous to tRNAs
Reads homologous to PLMVd
PLMVd (+) polarity
PLMVd (-) polarity
Subsequently, both libraries were analyzed for the presence of PLMVd-derived sequences via a homology search using the novel variant viroid sequences. Only 2 reads, corresponding to only 16 × 10-5 % of the total number (see Table 1), from the library derived from the healthy sample were found to be homologous with PLMVd sequences. This low number simply represents the background noise. Conversely, 884 234 sequences from the library derived from the infected sample were found to be homologous with the eighteen PLMVd sequences. This abundance of PLMVd-associated small RNAs corresponds to 11.2% of the reads of the infected library, and is significantly greater than what is observed for two Pospiviroidae family viroids (i.e. <5% for both the PSTVd  and the citrus exocortis viroid (CEVd)  that are known to accumulate in the nucleus). Among the sequences complementary to PLMVd, no apparent polarity bias was detected (i.e. (+) 423 951 and (-) 460 283). The relative abundances of the PLMVd replication intermediates of both polarities present in peach cells was previously estimated to be almost identical in a Siberian C cultivar (i.e. 2.2 ng and 2.0 ng for the (+) and (-) polarities respectively, per mg of tissues . This is the same order of magnitude, allowing for experimental error, as was previously suggested for other PLMVd variants in a distinct peach cultivar . This result confirms that if the PLMVd genome or its replicative intermediates, enter the plant siRNA pathways, both strands are susceptible to be implicated.
Relatively few species larger than 23 nt in size (Figure 2B) were found in the PLMVd-associated sequence subset. In fact, most of the reads were either 21 or 22 nt in size (765 766 sequences, representing 86.6% of all of the PLMVd derived sequences). Since Dicer-like enzymes cleave double-stranded RNA structures, it is tempting to assume that the viroid RNA genome (or replication intermediates) can be substrates of those enzymes in vivo as shown for PSTVd . Moreover, in vitro studies have shown that incubation of PLMVd RNA generates siRNAs in a cell-free wheat-germ extract . According to previous data using a genetic approach in Arabidopsis thaliana , more than one Dicer-like enzyme can produce siRNAs of 21-22 nt, rendering difficult and premature the identification of which one is involved with PLMVd silencing.
Clearly, the accumulated siRNAs differ significantly for the two PLMVd polarities. The secondary structures adopted by the PLMVd strands of both polarities [18, 25] differ, and may have an important impact on this phenomenon. It is likely that the accessibility of the double-stranded structure to the Dicer-like enzyme is an important criteria to consider. Surprisingly, both the P6 and P7 stems of the (+) polarity, which interact through the P8 pseudoknot that should reduce their accessibility to the enzyme, appear to constitute a preferred region as evidenced by the high abundance of siRNAs derived from this region. Also the P11 stem is a long double stranded RNA structure that could satisfy the Dicer-like enzymes for cleavage; yet, our results show that it is not necessarily the case. For example, only one strand of the P11 stem of both polarities was observed to preferentially accumulate siRNAs, while the other did not. One possible explanation might be that host proteins interact with the latter strand, thereby preventing Dicer-like binding. Interestingly, it is not the same strand of the P11 stem of both polarities that shows preferential accumulation. The reason for this phenomenon remains obscure. Moreover, it has been reported that disease symptoms induced by the hop stunt viroid (HSVd) in Nicotina benthamiana were independent of viroid accumulation levels, and instead relied on the involvement of the RNA-dependent RNA polymerase 6 (RDR6) . Based on this, it is possible that Dicer plays a relatively important role, but not as important as RDR6 potentially could, and that the over-represented regions in the pool of PLMVd derived small RNAs may actually reflect the action of RDR6 and other as yet unknown factors more than the activity of Dicer-like proteins. If this is indeed the case, we can speculate that the viroids may potentially exert their pathogenicity through the tasiRNA pathway involving RDR6 and DCL-4, as it has been previously hypothesized . Another aspect of the problem might be the preferential degradation of some siRNAs by the cellular machinery
Recent reports have shown that the sorting of small RNAs into Arabidopsis argonaute (AGO) complexes is influenced by the 5' terminal nucleotide [28–30]. In order to attribute the potential implication of specific AGO complexes, the nature of the 5' nucleotide of PLMVd-derived siRNAs of 21-22 nt was analyzed. PLMVd-derived RNA of 21 nt in size contain 48.7% of the time a U at the 5' end and 24.9% a C (i.e. C + U = 73.6%). Similarly, the PLMVd-derived siRNA of 22 nt harbor a U in the 5' end 39.4% of the time and a C in 29.2% (C + U = 68.6%). When the fraction of PLMVd-derived siRNA was removed from the library and the same analysis was repeated, the same bias for U and C was observed (C + U = 70.3% (19.0% + 51.3%) and 63.7% (23.1% + 40.6%) for the 21 and 22 nt ncRNA in size, respectively; data not shown). Conversely, the analysis of the library of small non-coding RNA obtained from the healthy leaves revealed a different bias for the nucleotide identity retrieved at the 5' end. The small RNAs of 21 as well as 22 nt in size showed a significant preference for a A and G (C + U = 43.0% and 36.8% for the 21 and 22 nt ncRNA in size, respectively, data not shown). In the latter case, the A was the most over-represented nucleotide in 5' (~ 40%). AGO1 and AGO5 were shown to preferentially recruit small RNAs beginning by a U and a C, respectively [29, 30]. Taking into account this demonstration, the bias towards a U or a C in the 5' terminal nucleotide suggests that PLMVd-derived small RNAs, if they are implicated in the plant RNA silencing pathways, may be recruited by AGO1 and AGO5, if homologous protein complexes exist in the Prunus genus. More generally, it is tempting to suggest that the active AGO proteins in the infected plant cell differ from those retrieved in the healthy ones. However, this does not exclude that other AGO proteins may be involved in the binding of those small RNAs since plants may encode up to ten different AGO proteins with redundant functions .
Lastly, in the previous analysis, the use of a minimal set of siRNAs , led to the proposal that a potential bias exists for the transportation of the circular conformers into the cytoplasm where they trigger further RNA silencing. This hypothesis is based on the observation of a slightly larger proportion of siRNAs that overlapped the self-cleavage site of PLMVd, a species that is by definition derived only from circular conformers. According to the alignment illustrated in Figure 3, this is no longer the case. The regions that include the hammerhead self-cleavage sites are not part of the hotspots on the structure. However, PLMVd-derived siRNAs of both polarities are retrieved from the infected plant. Also, no RNA silencing components within the chloroplast have been described so far. These components are localized mainly in the cytoplasm. So the viroid needs to be transported into that compartment to be processed by the RNA silencing machinery. These results suggest that, at least, the circular PLMVd conformer of both polarities is transported into the cytoplasm from the chloroplast. This might also be the case for the linear monomers that are the most abundant replication intermediates . Importantly, this study clearly shows that any conclusion based on a relatively small number of sequences should be taken with considerable caution.
For the first time, all of the RNA isolated from an infected tree have been used in parallel to sequence viroid species and small non coding RNA species. On the one hand, this work led to the determination of 18 novel PLMVd variants. On the other hand, the use of high throughput sequencing provided the largest library of siRNA reported yet for a viroid infection. It was shown that the entire PLMVd genome is retrieved in the siRNA population. Moreover, PLMVd strands of both polarities appeared susceptible for targeting, although specific regions for each one are over-represented in terms of siRNAs. Those regions, which are not the same for each polarity, do not necessarily correlate with double-stranded regions that could be substrates for Dicer-like enzymes. The nature of the first 5' nucleotide in the 21-22 nt population revealed a bias toward a C or a U in viroid-derived siRNAs. In fact this bias is also observed in the non-PLMVd-siRNAs. Unlike the infected library, the majority of the healthy sequences has A or G at their 5'. Importantly, the large libraries of siRNAs created in this study led to several observations concerning the silencing mechanism associated with PLMVd infection.
At the same time this manuscript was being submitted, another study describing deep sequencing of small RNAs associated to PLMVd infection was released . The authors have showed that PLMVd of both polarities can generate predominantly small 21-22 nt long RNAs in vivo. These small RNAs map the PLMVd genome adopting a hotspot profile which is different from one polarity to the other (plus vs minus), although their libraries of PLMVd-associated siRNAs were smaller. However, the Di Serio et al. study has the advantage of reporting data isolated from two peach plants. The present study offers two important differences: Firstly, the fact that the sequencing of the small ncRNAs was also performed on a healthy plant ensures that there are no species showing homologies with PLMVd. Secondly, the PLMVd sequences used to identify viroid-associated siRNAs were obtained from the same sample as the small RNAs instead of using sequences reported several years ago . It has been shown that PLMVd is a quasi-species and infection using a plasmid was sufficient for the observation of several sequences in a single plant . There is another recent study that includes high throughput sequencing of small ncRNAs from grapevine infected by hop stunt viroid and grapevine yellow speckle viroid 1, two viroids from the other family (i.e. Pospiviroidae family) . Clearly, high-throughput sequencing has become the method of choice for the identification of viroid-associated siRNAs.
RNA extraction and PLMVd analysis
Total RNA was extracted from the leaves of peach (Siberian C cultivar) and nectarine trees grown on the West Coast of Canada (Sidney, British Columbia) using Trizol (Invitrogen) according to the manufacturer's instructions. PLMVd can infect nectarine as well as peach trees . The sample from the nectarine plant should be considered as a negative control in which there is no direct comparison of the ncRNA library with the one from the PLMVd-infected peach plant. These RNA samples were quantified by UV spectroscopy, and their qualities verified by 1% agarose gel electrophoresis.
For PLMVd cloning, 5 μg of total RNA were incubated in the presence of 3 U of RQ1 DNAse (Promega) at 37°C for 1 h in a total volume of 20 μL containing 40 mM Tris-HCl pH 8.0, 10 mM MgSO4 and 1 mM CaCl2. Following two phenol-chloroform extractions and an ethanol precipitation, the RNA was reverse transcribed for 1 h at 45°C in the presence of the reverse primer 5'-AGTGCTCCGAATAGGGCACCCCAAGGTGG-3' using SuperScript III (Invitrogen) in 20 μL reactions as recommended by the manufacturer. Upon completion of the reverse transcription reactions, 10 μg of RNase A were added and the mixtures incubated for 15 min at room temperature. Ten microliters of the reverse transcriptase reactions were then used in PCR amplification reactions performed using purified Pfu DNA polymerase. The reactions were carried out in a final volume of 100 μL containing 2 μM of each DNA primer, 2 mM MgSO4, 200 μM dNTP, 10 mM Tris-HCl (pH 8.85), 25 mM KCl and 5 mM (NH4)2SO4. The reactions were incubated at 94°C for 1 min prior to the addition of the enzyme, and were then subjected to 25 amplification cycles (30 sec at 94°C, 60 sec at 50°C and 30 sec at 72°C). The reverse primer was the same as that used previously in the reverse transcription reactions, while the forward primer was 5'-GCAGTTCCCGATAGAAAGGCTAAGCACC-3'. The amplifications were analyzed by loading 10 μL aliquots on a 1% agarose gel. The gel band corresponding to a fragment of 338 nt in length was cut out and the DNA extracted using Spin-X columns (Corning Science Products). The resulting DNA was then incubated in the presence of Taq DNA polymerase (New England Biolabs) for 20 min at 72°C and ligated into pGEM-T vector (Promega). After transformation into E.coli, 23 clones were sequenced at the Platefome de séquençage et de génotypage des genomes (CRCHUL) (Québec, Canada). All sequences were deposited to NCBI, accession numbers GQ499305 to GQ499322.
Isolation of small RNAs and high-throughput sequencing
For each library (healthy and PLMVd-infected) 10 μg of total RNA was size fractionated on a 15% denaturing polyacrylamide gel and the 18-30 nt fraction excised from the gel. Following elution of the RNA, 5' and 3' adapters were added to the RNA species using the Digital Gene Expression for Small RNA Sample Prep Kit (Illumina). After reverse transcription, the samples were PCR amplified, purified on 10% polyacrylamide gels and then eluted according to the manufacturer's protocol. High-throughput Illumina Solexa sequencing was performed at The Center for Applied Genomics (Toronto, Canada).
Unfiltered sequence reads from both sets were trimmed of the 3' adapter. Reads that did not contain the 3' linker, and those shorter than 16 nt, were discarded. The resulting reads were aligned on the PLMVd sequences, as well as on other control RNA sequences, using the Exonerate software . No mismatches were permitted in the pairwise sequence comparison. All small RNA reads are available at the GEO database with accession number GSE18764.
This work was supported by a grant from the Natural Sciences and Engineering Research Council of Canada (NSERC Canada, grant number 155219-07) to J-PP. The RNA group is supported by grants from the Université de Sherbrooke and the Canadian Institutes of Health Research (CIHR, grant number PRG-80169). J-PP holds the Canada Research Chair in Genomics and Catalytic RNA. J-PP is a member of the Centre de Recherche Clinique Étienne-Lebel.
- Tsagris EM, Martínez de Alba AE, Gozmanova M, Kalantidis K: Viroids. Cell Microbiol. 2008, 10: 2168-2179. 10.1111/j.1462-5822.2008.01231.xView ArticlePubMedGoogle Scholar
- Bussière F, Lehoux J, Thompson DA, Skrzeczkowski J, Perreault JP: Subcellular localization and rolling circle replication of peach latent mosaic viroid: hallmarks of group A viroids. J Virol. 1999, 73: 6353-6360.PubMedPubMed CentralGoogle Scholar
- Bonnet E, Peer Van de Y, Rouzé P: The small RNA world of plants. New Phytol. 2006, 171: 451-468.View ArticlePubMedGoogle Scholar
- Papaefthimiou I, Hamilton A, Denti M, Baulcombe D, Tsagris M, Tabler M: Replicating potato spindle tuber viroid RNA is accompanied by short RNA fragments that are characteristic of post-transcriptional gene silencing. Nucleic Acids Res. 2001, 29: 2395-2400. 10.1093/nar/29.11.2395View ArticlePubMedPubMed CentralGoogle Scholar
- Itaya A, Folimonov A, Matsuda Y, Nelson RS, Ding B: Potato spindle tuber viroid as inducer of RNA silencing in infected tomato. Mol Plant Microbe Interac. 2001, 14: 1332-1334. 10.1094/MPMI.2001.14.11.1332..View ArticleGoogle Scholar
- Martinez de Alba AE, Flores R, Hernandez C: Two chloroplastic viroids induce the accumulation of small RNAs associated with posttranscriptional gene silencing. J Virol. 2002, 76: 13094-13096. 10.1128/JVI.76.24.13094-13096.2002View ArticlePubMedPubMed CentralGoogle Scholar
- Denti MA, Boutla A, Tsagris M, Tabler M: Short interfering RNAs specific for potato spindle tuber viroid are found in the cytoplasm but not in the nucleus. Plant J. 2004, 37: 762-769. 10.1111/j.1365-313X.2004.02001.xView ArticlePubMedGoogle Scholar
- Landry P, Thompson D, Perreault JP: The role of viroids in gene silencing: the model case of peach latent mosaic viroid. Can J Plant Pathol. 2004, 26: 31-38.View ArticleGoogle Scholar
- Wang MB, Bian XY, Wu LM, Liu LX, Smith NA, Isenegger D, Wu RM, Masuta C, Vance VB, Watson JM, Rezaian A, Dennis ES, Waterhouse PM: On the role of RNA silencing in the pathogenecity and evolution of viroids and viral satellites. Proc Natl Acad Sci USA. 2004, 101: 3275-3280. 10.1073/pnas.0400104101View ArticlePubMedPubMed CentralGoogle Scholar
- Itaya A, Zhong X, Bundschuh R, Qi Y, Wang Y, Takeda R, Harris AR, Molina C, Nelson RS, Ding B: A structured viroid RNA serves as a substrate for dicer-like cleavage to produce biologically active small RNAs but is resistant to RNA-induced silencing complex-mediated degradation. J Virol. 2007, 81: 2980-2994. 10.1128/JVI.02339-06View ArticlePubMedPubMed CentralGoogle Scholar
- Carbonell A, Martínez de Alba AE, Flores R, Gago S: Double-stranded RNA interferes in a sequence-specific manner with the infection of representative members of the two viroid families. Virology. 2008, 371: 44-53. 10.1016/j.virol.2007.09.031View ArticlePubMedGoogle Scholar
- Schwind N, Zwiebel M, Itaya A, Ding B, Wang MB, Krczal G, Wassenegger M: RNAi-mediated resistance to potato spindle tuber viroid in transgenic tomato expressing a viroid hairpin RNA construct. Mol Plant Pathol. 2009, 10: 459-469. 10.1111/j.1364-3703.2009.00546.xView ArticlePubMedGoogle Scholar
- Landry P, Perreault JP: Identification of a peach latent mosaic viroid hairpin able to act as a Dicer-like substrate. J Virol. 2005, 79: 6540-6543. 10.1128/JVI.79.10.6540-6543.2005View ArticlePubMedPubMed CentralGoogle Scholar
- St-Pierre P, Hassen IF, Thompson D, Perreault JP: Characterization of the siRNAs associated with peach latent mosaic viroid infection. Virology. 2009, 383: 178-182. 10.1016/j.virol.2008.11.008View ArticlePubMedGoogle Scholar
- Hall N: Advanced sequencing technologies and their wider impact in microbiology. J Exp Biol. 2007, 210 (Pt 9): 1518-1525. 10.1242/jeb.001370View ArticlePubMedGoogle Scholar
- Morin RD, O'Connor MD, Griffith M, Kuchenbauer F, Delaney A, Prabhu AL, Zhao Y, McDonald H, Zeng T, Hirst M, Eaves CJ, Marra MA: Application of massively parallel sequencing to microRNA profiling and discovery in human embryonic stem cells. Genome Res. 2008, 18: 610-621. 10.1101/gr.7179508View ArticlePubMedPubMed CentralGoogle Scholar
- Rocheleau L, Pelchat M: The Subviral RNA Database: a toolbox for viroids, the hepatitis delta virus and satellite RNAs research. BMC Microbiol. 2006, 6: 24- 10.1186/1471-2180-6-24View ArticlePubMedPubMed CentralGoogle Scholar
- Bussière F, Ouellet J, Côté F, Lévesque D, Perreault JP: Mapping in solution shows the peach latent mosaic viroid to possess a new pseudoknot in a complex, branched secondary structure. J Virol. 2000, 74: 2647-2654. 10.1128/JVI.74.6.2647-2654.2000View ArticlePubMedPubMed CentralGoogle Scholar
- Slater GS, Birney E: Automated generation of heuristics for biological sequence comparison. BMC Bioinformatics. 2005, 6: 31- 10.1186/1471-2105-6-31View ArticlePubMedPubMed CentralGoogle Scholar
- Kawaji H, Nakamura M, Takahashi Y, Sandelin A, Katayama S, Fukuda S, Daub CO, Kai C, Kawai J, Yasuda J, Carninci P, Hayashizaki Y: Hidden layers of human small RNAs. BMC Genomics. 2008, 9: 157- 10.1186/1471-2164-9-157View ArticlePubMedPubMed CentralGoogle Scholar
- Xie Z, Johansen LK, Gustafson AM, Kasschau KD, Lellis AD, Zilberman D, Jacobsen SE, Carrington JC: Genetic and functional diversification of small RNA pathways in plants. PLoS Biol. 2004, 2: E104- 10.1371/journal.pbio.0020104View ArticlePubMedPubMed CentralGoogle Scholar
- Xie Z, Qi X: Diverse small RNA-directed silencing pathways in plants. Biochim Biophys Acta. 2008, 1779: 720-724.View ArticlePubMedGoogle Scholar
- Nobuta K, Lu C, Shrivastava R, Pillay M, De Paoli E, Accerbi M, Arteaga-Vazquez M, Sidorenko L, Jeong DH, Yen Y, Green PJ, Chandler VL, Meyers BC: Distinct size distribution of endogeneous siRNAs in maize: Evidence from deep sequencing in the mop1-1 mutant. Proc Natl Acad Sci USA. 2008, 105: 14958-14963. 10.1073/pnas.0808066105View ArticlePubMedPubMed CentralGoogle 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: 135-146. 10.1016/j.virol.2007.05.011View ArticlePubMedGoogle Scholar
- Dubé A, Baumstark T, Bisaillon M, Perreault JP: The RNA strands of the plus and minus polarities of Peach Latent Mosaic Viroid fold into different structures. RNA. 2010, 16: 463-473. 10.1261/rna.1826710View ArticlePubMedPubMed CentralGoogle Scholar
- Gómez G, Martínez G, Pallás V: Viroid-induced symptoms in Nicotiana benthamiana plants are dependent on RDR6 activity. Plant physiol. 2008, 148: 414-423. 10.1104/pp.108.120808View ArticlePubMedPubMed CentralGoogle Scholar
- Gómez G, Martínez G, Pallás V: Interplay between viroid-induced pathogenesis and RNA silencing pathways. Trends Plant Sci. 2009, 14: 264-269. 10.1016/j.tplants.2009.03.002View ArticlePubMedGoogle Scholar
- Montgomery TA, Howell MD, Cuperus JT, Li D, Hansen JE, Alexander AL, Chapman EJ, Fahlgren N, Allen E, Carrington JC: Specificity of ARGONAUTE7-miR390 interaction and dual functionality in TAS3 trans-acting siRNA formation. Cell. 2008, 133: 128-141. 10.1016/j.cell.2008.02.033View ArticlePubMedGoogle Scholar
- Mi S, Cai T, Hu Y, Chen Y, Hodges E, Ni F, Wu L, Li S, Zhou H, Long C, Chen S, Hannon GJ, Qi Y: Sorting of small RNAs into Arabidopsis argonaute complexes is directed by the 5' terminal nucleotide. Cell. 2008, 133: 116-127. 10.1016/j.cell.2008.02.034View ArticlePubMedPubMed CentralGoogle 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: 493-500. 10.1093/pcp/pcn043View ArticlePubMedGoogle Scholar
- Di Serio F, Gisel A, Navarro B, Delgado S, Martinez de Alba A-E, Donvito G, Flores R: Deep sequencing of the small RNAs derived from two symptomatic variants of a chlroplastic viroid: implications for their genesis and for pathogenesis. PLoS One. 2009, 4: e7539- 10.1371/journal.pone.0007539View ArticlePubMedPubMed CentralGoogle Scholar
- Ambros S, Hernandez C, Desvignes JC, Flores R: Genomic structure of three phenotypically different isolates of peach latent mosaic viroid: Implications of the existence of constraints limiting the heterogeneity of viroid quasispecies. J Virol. 1998, 72: 7397-7406.PubMedPubMed CentralGoogle Scholar
- Navarro B, Pantaleo V, Gisel A, Moxon S, Dalmay T, Bisztray G, Di Serio F, Burgyán J: 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: e7686- 10.1371/journal.pone.0007686View ArticlePubMedPubMed CentralGoogle Scholar
- Skrzeczkowski LJ, Howell WE, Mink GI: Occurrence of peach latent mosaic viroid in commercial peach and nectarine cultivars in the U.S. Plant Dis. 1996, 80: 823-View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.