- Open Access
Characterization of virus-derived small interfering RNAs in Apple stem grooving virus-infected in vitro-cultured Pyrus pyrifolia shoot tips in response to high temperature treatment
Virology Journalvolume 13, Article number: 166 (2016)
Heat treatment (known as thermotherapy) together with in vitro culture of shoot meristem tips is a commonly used technology to obtain virus-free germplasm for the effective control of virus diseases in fruit trees. RNA silencing as an antiviral defense mechanism has been implicated in this process. To understand if high temperature-mediated acceleration of the host antiviral gene silencing system in the meristem tip facilitates virus-derived small interfering RNAs (vsiRNA) accumulation to reduce the viral RNA titer in the fruit tree meristem tip cells, we used the Apple stem grooving virus (ASGV)–Pyrus pyrifolia pathosystem to explore the possible roles of vsiRNA in thermotherapy.
At first we determined the full-length genome sequence of the ASGV-Js2 isolate and then profiled vsiRNAs in the meristem tip of in vitro-grown pear (cv. ‘Jinshui no. 2’) shoots infected by ASGV-Js2 and cultured at 24 and 37 °C. A total of 7,495 and 7,949 small RNA reads were obtained from the tips of pear shoots cultured at 24 and 37 °C, respectively. Mapping of the vsiRNAs to the ASGV-Js2 genome revealed that they were unevenly distributed along the ASGV-Js2 genome, and that 21- and 22-nt vsiRNAs preferentially accumulated at both temperatures. The 5′-terminal nucleotides of ASGV-specific siRNAs in the tips cultured under different temperatures had a similar distribution pattern, and the nucleotide U was the most frequent. RT-qPCR analyses suggested that viral genome accumulation was drastically compromised at 37 °C compared to 24 °C, which was accompanied with the elevated levels of vsiRNAs at 37 °C. As plant Dicer-like proteins (DCLs), Argonaute proteins (AGOs), and RNA-dependent RNA polymerases (RDRs) are implicated in vsiRNA biogenesis, we also cloned the partial sequences of PpDCL2,4, PpAGO1,2,4 and PpRDR1 genes, and found their expression levels were up-regulated in the ASGV-infected pear shoots at 37 °C.
Collectively, these results showed that upon high temperature treatment, the ASGV-infected meristem shoot tips up-regulated the expression of key genes in the RNA silencing pathway, induced the biogenesis of vsiRNAs and inhibited viral RNA accumulation. This study represents the first report on the characterization of the vsiRNA population in pear plants infected by ASGV-Js2, in response to high temperature treatment.
Virus-induced gene silencing is an antiviral defense mechanism in plants, in which the vital elements involved are virus-derived small interfering RNAs (vsiRNAs) which are mainly 20–24 nucleotides (nt) in length [1, 2]. The vsiRNAs are produced by Dicer-like proteins (DCLs) from viral double-stranded RNA (dsRNA) replication intermediates during the viral replication process, or from highly structured single-stranded RNA molecules present in virus-infected host plants [1–4]. These vsiRNAs are templates generated in host cells that are loaded into Argonaute proteins (AGOs)-containing RNA-induced silencing complexes (RISCs), which then guide the target genomic and subgenomic viral mRNA in a sequence-specific manner to interfere with virus replication, translation and movement, and, in some cases, eliminate the viral infection [3, 5–11].
It is well known that DCLs, AGOs, and the RNA-dependent RNA polymerases (RDRs) participate in the antiviral silencing pathways as key silencing factors, and their RNA silencing activities are varied under different temperatures [4, 12–15]. In the model plant Arabidopsis, among the four relatively well-studied DCLs, DCL4 and DCL2 function in RNA silencing against RNA viruses by producing 21- and 22-nt vsiRNAs, respectively. It has been demonstrated that the activity of DCL2 in producing specific 22-nt vsiRNAs derived from Turnip crinkle virus (TCV) in Arabidopsis is enhanced by higher temperatures [16, 17]. DCLs-generated vsiRNAs are associated with specific AGO complexes, a process partially dependent on the 5'-terminal nucleotides. For instances, vsiRNAs with a 5'-terminal uridine or adenosine are recruited preferentially by AGO1 and AGO2 . Recent studies have also shown that AGO2 plays an antiviral role in the temperature-dependent survival of TCV- and Potato virus X (PVX)-infected Arabidopsis plants [17, 19, 20]. In addition, RDRs like RDR1, RDR2, or RDR6 are involved in the biogenesis of secondary vsiRNAs to further enhance the antiviral RNA silencing efficiency [21–23]. In N. benthamiana, silencing of RDR6 increases viral RNA accumulation and facilitates viruses to invade the meristem tissue . High temperatures intensify the RDR6 activity in the antiviral RNA-silencing defense response . It has been documented that RDR6 plays a tissue-specific role in the inhibition of Chinese wheat mosaic virus (CWMV) accumulation and vsiRNA biogenesis at higher temperatures [7, 15, 24, 25].
Apple stem grooving virus (ASGV), a member of the genus Capillovirus in the family Betaflexiviridae . The ASGV genome is a positive-sense ssRNA with 6.5 kb in length that contains two overlapping open reading frames (ORFs). The larger ORF1 encodes a polyprotein of 240 kDa, in which the N-terminal region contains replicase domains including methyltransferase (Met), papain-like protease (P-pro), NTP-binding helicase (Hel), and the RNA dependent RNA polymerase (RdRp), and the C-terminal region is the coat protein (CP) of 27 kDa . ORF2 is embedded within the ORF1 and encodes a movement protein (MP) of approximately 36 kDa . MP and CP may be produced through the 3′-coterminal subgenomic RNAs (sgRNAs), and CP expression from sgRNA is essential for ASGV systemic infection in the host [28, 29]. Phylogenetic analysis of 16 ASGV full-length genomic sequences clusters them into two groups with no correlations to host and geographical origins [30, 31]. ASGV infection is symptomless on most commercial cultivars of apple and pear, but does induce the typical symptoms of stem pitting and grooving on some cultivars of citrus, lily, kiwifruit, and pear [32–35]. In asymptomatic apple plantlets, ASGV infection induces global gene expression changes, suggesting that extensive host genome-wide gene expression changes do not necessarily lead to disease symptoms . In pear, ASGV infection often deteriorates fruit quality . In the past several years, an increasing incidence of ASGV infection was observed in the pear-growing areas of China, leading to substantial economic losses [37, 38]. High temperature in combination with shoot meristem tip culture is an effective way to obtain virus-free germplasm to control viral diseases of fruit trees [39, 40]. The absence of viruses in the shoot meristem tip tissues is of practical importance because virus-free clones can be generated from infected shoots by culturing excised meristem tips. The effect of temperature on the RNA-silencing activities in plants has been investigated. Accumulated evidence suggests that low temperature inhibits RNA silencing-mediated defense by limiting the generation of small interfering RNA (siRNA) molecules, and high temperature promotes this innate immunity via increasing siRNA accumulation levels [15, 41–43]. In virus-infected plants, viral symptoms disappear in new leaves at high temperatures, resulting from the higher temperature-mediated acceleration of the host antiviral gene silencing system in the meristem tip [9, 39, 44]. In a recent study, we found that thermotherapy elimination of ASGV from Asian pear is associated with the high temperature-induced mixed action of a number of miRNA-mediated target genes related to disease defense and hormone signal transduction pathways in the apical meristem of pear shoots . These data suggest that elevated temperatures may enhance vsiRNA-mediated antiviral gene silencing activity, which in turn reduces the accumulation of viral RNA in the infected meristem tip cells of in vitro-cultured pear shoots.
To explore the possible roles of vsiRNA in the interaction of ASGV and pear plants at higher temperatures, we determined the complete genomic sequence of the ASGV-Js2 isolate and profiled the vsiRNAs in the meristem tip of in vitro-cultured pear shoots at 24 and 37 °C. We analyzed the distribution of vsiRNAs along the viral genome. The corresponding accumulation levels in relation to viral titer and the expression profiles of vsiRNA were also characterized in the shoot tip tissues by RT-qPCR. Furthermore, we determined the PpDCL2,4, PpAGOs, and PpRDR1 mRNA sequences and assessed their relative accumulation levels in the ASGV-infected pear shoot tips at 24 and 37 °C. Collectively, these results assist in a better understanding of the roles of vsiRNAs in ASGV infection in vitro-grown pear meristem tips in response to high temperature.
In vitro-grown P. pyrifolia cv. ‘Jinshui no. 2’ shoots infected with ASGV or ASGV-free, were confirmed by RT-PCR (Additional file 1: Figure S1). ASGV-infected and ASGV-free ‘Jinshui no. 2’ shoots ~1 cm in length were cut and transferred to MS medium, and treated at a 37 °C thermotherapy chamber (16 h light of 1500 lux and 8 h dark) or cultured at 24 °C as controls.
Bioinformatics analysis of vsiRNA sequences
Clean small RNA reads were obtained from previously constructed sRNA libraries derived from T24 and T37, representing the mixed samples containing equal amounts of total RNAs isolated from the pear shoots treated for 1 and 5 days at 24 and 37 °C, respectively . The small RNA reads were then aligned to the ASGV-Js2 genome (serving as a reference sequence) obtained in this study (see below). The vsiRNA profiles along the viral sense and antisense genome were determined using the Perl scripts and Bowtie software with a stringency allowing no more than two mismatches per read [46, 47].
Amplification and cloning of the full-length ASGV genome from P. pyrifolia cv. ‘Jinshui no. 2’ shoots
The full-length genome sequence of the ASGV-Js2 isolate sequence was determined by standard PCR with specific primers (Additional file 2: Table S1) designed based on the full-length genomic sequences of 17 previously reported ASGV isolates, followed by 5'-RACE and 3'-RACE (Takara Biotechnology Company, Dalian, China) to obtain the 5'- and 3'-terminal regions of the viral genomic RNA. The PCR products were gel-purified and cloned into the pMD18-T vector (Takara). At least three independent overlapping clones in both orientations were sequenced at Jinsirui Biotechnology and Service Co. Ltd (Nanjing, Jiangsu province, China). The full-length sequence of the virus isolate was assembled using clones derived from overlapping RT-PCR fragments by the program Vector NTI 10.0 (Invitrogen, USA).
Identification and cloning of PpDCL2,4, PpAGO1,2,4, and PpRDR1 sequences from P. pyrifolia
Primers used for amplification were designed based on the pear (Pyrus bretschneideri Rehd.) genome database  using Oligo7  and are listed (Additional file 3: Table S2). PpDCL2,4, PpAGO1,2,4, and PpRDR1 sequences were amplified by RT-PCR using specific primers from total RNA as template. The PCR products were purified and cloned. The PpDCL2,4, PpAGO1,2,4, and PpRDR1 DNA sequences were deposited into GenBank and their accession numbers are given in Additional file 3: Table S2.
Phylogenetic tress Sequence analysis
Sequence similarity searches were performed with the NCBI BLASTN and BLASTX online programs. Pairwise alignments of nucleotide and amino acid sequence were made to determine sequence identity and similarity. Phylogenetic analyses based on multiple sequence alignments of the nucleotide sequences and the predicted protein sequences, respectively, were performed using the Neighbour-joining (NJ) method incorporated in the programs Clustal X 1.83 and MEGA6 software [50, 51].
Quantitative real-time PCR analyses
To detect the vsiRNA expression levels obtained from the high-throughput sRNA sequencing, real-time quantitative PCR (RT-qPCR) was performed by One Step microRNA Prime-Script sRNA cDNA Synthesis Kit with Poly (A) tails and a special oligo-dT adaptor according to supplier’s instructions (Takara) as described by Liu et al. . Total RNAs from 5 mm ASGV-infected or virus-free meristem tips cultured at 37 and 24 °C for 1 and 5 days were extracted using the CTAB method [45, 52]. All specific vsiRNA forward primers were designed based on the mature vsiRNA sequence (Additional file 4: Table S3). In addition, the expression levels of ASGV mp, PpDCL2,4, PpAGO1,2,4, and PpRDR1 were analyzed by real-time PCR. Two or three sets of primers were designed and estimated for each gene. The best primer set for each gene is provided in Additional file 4: Table S3. Each sample was a mixture of equal amounts of total RNA from the ASGV-infected or virus-free shoots treated at 24 and 37 °C for 1 and 5 days. RNA samples were digested with DNAase I, reverse-transcribed into cDNA, and then used as templates for RT-qPCR as described previously . The Actin gene was used as an internal reference gene for normalization of ASGV mp, PpDCL2,4, PpAGOs, and PpRDR1 gene expression levels. The RT-qPCR for each gene was performed for three replicates. The experimental data were analyzed essentially as described previously [45, 53].
Cloning and sequencing of the ASGV-Js2 genome from in vitro-grown shoots of P. pyrifolia cv. ‘Jinshui no. 2’
To profile vsiRNAs derived from the ASGV-Js2 isolate in pear, we needed to determine the full-length genome sequence of the ASGV-Js2 isolate. The full-length ASGV-Js2 genome was cloned and sequenced as described in Methods (Additional file 2: Table S1). The genomic sequence was deposited into GenBank (accession no. KU198289). The ASGV-Js2 viral genome is 6,497 bp in length and contains two overlapping ORFs (ORF1 and ORF2) (Fig. 1a), similar to the data reported previously [27, 28]. A phylogenetic analysis of ASGV-Js2 and 17 additional full-length ASGV genomic sequences available in GenBank showed that ASGV-Js2 defines a unique cluster, and the pairwise nucleotide identities range from 80.0 to 83 % within this phylogenetic ASGV group (Fig. 1b and Additional file 5: Table S4).
Characteristics of the ASGV-derived vsiRNA population in P. pyrifolia shoot tips in response to high temperature
Small RNAs obtained from previously constructed libraries of in vitro-cultured P. pyrifolia shoots infected by ASGV-Js2 at 24 and 37 °C were subjected to deep-sequencing on an Illumina platform by Liu et al. . To determine whether siRNAs identified in pear derived from ASGV-Js2 are responsive to high temperature treatment, the generated siRNA sequences were aligned against the ASGV genomic and the negative-sense genomic RNA sequences. A total of 7,495 and 7,949 reads ranging from 18 to 26 nt in length were identified to be ASGV-derived siRNAs from the 24 and 37 °C libraries, respectively. These vsiRNAs were mapped to the ASGV-Js2 genome sequence. An analysis of the length distribution of the vsiRNAs showed that the 21-nt vsiRNAs were the most abundant, representing 61 % (4,592 reads) and 62 % (4,928 reads) of total vsiRNAs from the 24 °C- and 37 °C-treated ASGV-infected pear shoot libraries, respectively (Fig. 2a). The next most abundant class was the 22 nt molecules, accounting for 31.5 % (2,362/7,495 reads) and 30 % (2,388/7,949 reads) of total vsiRNAs from the 24 and 37 °C libraries, respectively (Fig. 2a). Relatively more vsiRNAs (4,203 from the 24 °C library and 4,447 from the 37 °C libraries) were found to be derived from the positive strand (Fig. 2b). Fewer vsiRNAs (3,292 and 3,502) from both libraries were derived from the negative strand (Fig. 2b). The ratio of sense/antisense 21-nt v siRNAs was 1.21 and 1.25 for the 24 and 37 °C libraries, respectively (Fig. 2b). The relatively high proportion of 21- and 22-nt vsiRNAs confirmed the presence of DCL4- and DCL2-like Dicer ribonucleases in pear shoot tips that target viral RNAs to mediate vsiRNA biogenesis.
Previous data have shown that the 5'-terminal nucleotides partially influence the loading of siRNAs into specific AGO complexes in Arabidopsis, and may have roles in targeting vsiRNAs to AGOs . In this study, we analyzed the relative frequencies of the 5'-terminal nucleotides of vsiRNAs from the 24 and 37 °C libraries in an effort to understand potential interactions between vsiRNAs and AGOs in pear. Among all the vsiRNAs, U was found to be the most common nucleotide at the 5' end, with 55.32 % and 56.92 % derived from the ASGV positive and negative strands from in vitro-grown pear shoots at 24 °C, respectively (Fig. 3a). U was also the most abundant 5' terminal nucleotide (58.88 and 57.35 % corresponding to the ASGV positive and negative strands, respectively) in the 37 °C-treated shoots (Fig. 3b). In addition, detailed analyses of the relative frequencies of the 5'-terminal nucleotides for the 20 to 24-nt vsiRNAs revealed that the vsiRNAs ranging from 20 to 23-nt showed a preference for U, with 43.75 to 59.09 % and 40.77 to 60.77 % of those derived from the genomic ASGV RNA in the 24 and 37 °C libraries, respectively. However, the 24-nt vsiRNAs showed a preference for A at the 5' end, with 45.10 and 54.55 % in the 24 and 37 °C libraries (Fig. 4). Also, the 20 to 24-nt vsiRNA had a preference for U with 39.02–58.89 and 37.50–58.63 % derived from the negative ASGV strand in the 24 and 37 °C libraries, respectively (Fig. 4). These results support the idea that the 21- and 22-nt vsiRNAs with U as the 5'-terminal nucleotide might be preferentially loaded into AGO1 , consistent with the role described for AGO1 in defending against RNA virus infection in Arabidopsis .
Distribution of ASGV-derived siRNA abundance in P. pyrifolia shoot tips in response to high temperature
To determine the distribution pattern of the vsiRNAs along the ASGV-Js2 genome, vsiRNA sequences were mapped to the viral genome. The results showed that the vsiRNAs were distributed unevenly along both strands of the ASGV-Js2 genome, implying there are variations in the relative abundance of siRNAs targeting different regions of the ASGV-Js2 genome in the 24- and 37 °C-treated pear shoots (Fig. 5a and b). No apparent vsiRNA-production hotspots in the ASGV-Js2 genome were found from our analyses.
To validate the existence of the predicted vsiRNAs in the pear shoot meristem tip using small RNA sequence analysis, some vsiRNAs derived from different ASGV-Js2 genomic regions in the two libraries were analyzed by RT-qPCR. These vsiRNAs matched the ASGV-Js2 genome positive strand at positions 85–105, 2931–2951 and 4379–4399, designated vsiRNA85(+), vsiRNA2931(+), and vsiRNA4379(+), respectively. Similarly, some reads were mapped to the negative strand at positions 5839–5859 and 4625–4645, designated vsiRNA5839(−) and vsiRNA4625(−), respectively. RT-qPCR results discovered that the expression levels of vsiRNA2931(+) and vsiRNA4625(−) in the shoot meristem tip tissue at 37 °C were nearly equal to those at 24 °C (Fig. 6). In contrast, the expression levels of vsiRNA85(+), vsiRNA4379(+) and vsiRNA5839(−) were elevated at 37 °C, when compared to those from the 24 °C-treated shoot tips (Fig. 6).
The effects of high temperature on ASGV infection in P. pyrifolia shoot tips
ASGV infection may be evaluated by quantification of the level of ASGV genomic RNA and subgenomic RNA. As mentioned earlier, MP may be produced via subgenomic RNA. To investigate the effect of high temperature on ASGV infection in the meristem tip of in vitro-grown P. pyrifolia shoots, we determined the ASGV genomic RNA and its MP subgenomic RNA accumulation at 24 and 37 °C by RT-qPCR. The level of viral genomic RNA and MP subgenomic RNA was found to decrease by 50 % in response to the 37 °C treatment compared to 24 °C (Fig. 7). This indicated that high temperature drastically inhibited ASGV infection in the pear meristem tip.
Expression profiles of PpDCL2, PpDCL4, PpAGO1, PpAGO2, PpAGO4, and PpRDR1 mRNAs in ASGV-infected P. pyrifolia shoot tips at high temperature
In plants, the generation of small RNAs mainly depends on the activities of DCL, AGO, and RDR. Dicer and AGOs are the key factors involved in recognition of dsRNA and degrading target viral RNAs, whereas RDRs mediate the secondary accumulation of siRNAs through de novo synthesis [2, 4]. To study the expression of these genes in the ASGV-infected pear shoot tips in response to high temperature treatment, we cloned partial cDNA sequences of PpDCL2, PpDCL4, PpAGO1, PpAGO2, PpAGO4, and PpRDR1 mRNA and the resulting sequences were deposited into the GenBank database under accession numbers provided in Additional file 3: Table S2. Phylogenetic tree analysis showed that the deduced amino acid sequences of the partial PpDCL2, PpDCL4, PpAGO1, PpAGO2, PpAGO4, and PpRDR1 cDNA sequences share high levels of sequence conservation with their counterparts from Pyrus bretschneideri, Malus, Arabidopsis, Nicotiana, and Oryza sativa (Fig. 8a, b and c). These data support the notion that PpDCL2, PpDCL4, PpAGO1, PpAGO2, PpAGO4, and PpRDR1 genes in P. pyrifolia are homologs of the corresponding genes from Nicotiana, Arabidopsis, and Oryza sativa, and thus may function like DCLs, AGOs, and RDRs as demonstrated in those model plant species [2, 16, 17, 20].
Based on the cloned cDNA sequences of PpDCL2, PpDCL4, PpAGO1, PpAGO2, PpAGO4, and PpRDR1, primers were designed for RT-qPCR to determine their expression levels (Additional file 4: Table S3). We found that PpDCL2 and PpDCL4 mRNA increased dramatically in the ASGV-infected pear shoots at the higher temperature, as compared to the controls (Fig. 9a). Similarly, PpRDR1 expression was also up-regulated in ASGV-infected pear shots at 37 °C (Fig. 9c). Interestingly, it is ASGV infection rather than high temperature treatment significantly up-regulated PpAGO1 and PpAGO4 expression (Fig. 9b). In contrast, high temperature did induce higher levels of PpAGO2-specific mRNA accumulation (Fig. 9b).
Viral infection triggers the host gene silencing response, leading to vsiRNA accumulation [54–60]. Temperature plays a key role in the RNA silencing-mediated antiviral defense in plants due to its effect on the control of siRNA generation [15, 17, 42, 43]. Although vsiRNAs associated with virulent virus infections under various environmental conditions have been relatively well studied in model plant species, their involvement in latent virus infection in woody plants, especially fruit trees remains poorly characterized [61, 62]. Our study is the first report on the characterization of a population of vsiRNAs in pear plants infected by a latent virus, ASGV-Js2, in response to high temperature treatment.
In this study, small RNA sequencing data revealed the presence of vsiRNAs in the tip of the ASGV-infected P. pyrifolia shoots (Figs. 2, 3 and 4). We found that approximately 0.05 % (7,495/13,741,468 reads) and 0.06 % (7,949/14,071,933 reads) reads from the ASGV-infected samples matched the ASGV-Js2 genome from the 24 and 37 °C libraries, respectively. The results revealed that vsiRNAs accounted for a relatively small proportion of small RNAs in the ASGV-infected pear shoots when compared to the levels of vsiRNAs in other virus-host pathosystems in which virus-infected leaves are often used as materials for siRNA profiling studies [11, 57]. Therefore, the low levels of ASGV-derived vsiRNAs may be attributed to the meristem tips of in vitro-cultured pear shoots used in this study. It is well known that there are endogenous restrictions preventing viral genomes from moving into plant meristems. In this study, we also found a continuous and uneven distribution of plus- and minus-sense vsiRNAs throughout the ASGV-Js2 genome in pear at either 24 or 37 °C (Figs. 2b and 5). This is different from findings from several previous reports that vsiRNAs are increasingly distributed toward the 3' end of the viral genome [61–63]. It is not clear if the distribution pattern of vsiRNAs found in this study is also related to the particular tissues, e.g., pear shoots used in this study.
Real-time PCR analysis using primers specific for the MP coding region (Fig. 7) or RdRp sequence (data not shown) demonstrated a reduction in ASGV accumulation in the pear meristem tissues in response to high temperature treatment, similar to what was found in our previous study . Also the accumulation of vsiRNA4379 (+) and vsiRNA5839 (−) derived from the RdRp and CP regions of the ASGV genome increased accordingly (Figs. 5 and 6). These data suggest that vsiRNA abundance is negatively correlated with the levels of the ASGV viral RNA in response to high temperature (37 °C) treatment. This is consistent with the results obtained using other virus-host pathosystems . An exemption is the case of a DNA virus, Cucurbit leaf crumble virus (CLCV). The relative abundance of CLCV-derived siRNA is apparently positively correlated with viral titers in pumpkin . A possible explanation for this discrepancy is that each virus-host combination might reflect unique characteristics; specifically, a dynamic equilibrium established during viral infection may affect sRNA levels in different virus-host systems [54, 55, 57, 64]. It has been demonstrated that DCLs, RDRs, AGOs and other factors involved in RNA silencing also participate in antiviral defense in model plant species such as Arabidopsis, tobacco, and rice [2, 13, 15, 20, 23, 25, 65–67]. To explore the association of the corresponding homolog proteins with ASGV infection in pear in response to high temperature treatment, we cloned and obtained the partial sequences of PpDCL2, PpDCL4, PpAGO1, PpAGO2, PpAGO4, and PpRDR1 and determined their relative expression levels (Figs. 8 and 9). Overall these genes were up-regulated in the ASGV-infected pear shoots at 37 °C, which was accompanied with the reduced level of viral RNA (Fig. 7) and the elevated levels of vsiRNAs (Fig. 6). These data support that high temperature treatment may enhance the RNA silencing capacity in the pear meristem tissue via up-regulating the expression of key components of the antiviral pathway to cope with ASGV infection. Future study is directed to elucidate how temperature regulates gene silencing in the ASGV-infected pear shoots and if this is tissue- or ASGV-specific. Such work would help better understand and improve thermotherapy for the effective control of virus diseases in fruit trees.
This study represents the first report on the characterization of vsiRNA in the in vitro-grown ASGV-infected pear shoots in response to high temperature treatment. The profiles of vsiRNAs showed an uneven distribution along the ASGV-Js2 genome, and that 21- and 22-nt vsiRNAs preferentially accumulated when cultured at higher temperature. ASGV-specific siRNAs from all libraries had a similar distribution of 5'-terminal nucleotides. U was the most frequent among the 5’ terminal nucleotides, and its frequency was slightly higher at 37 °C. The expression levels of the viral mp gene and vsiRNAs were characterized by RT-qPCR. We also cloned PpDCL2,4, PpAGO1,2,4 and PpRDR1 partial sequences and examined their expression patterns, and found their expression levels were up-regulated in the ASGV-infected pear shoots at 37 °C. This up-regulation was accompanied with the reduced level of viral RNA and the elevated levels of vsiRNAs. Taken together these data suggest that high temperature may induce and enhance the RNA silencing capacity in the pear meristem tissue.
Hamilton AJ, Baulcombe DC. A species of small antisense RNA in posttranscriptional gene silencing in plants. Science. 1999;286:950–2.
Ding SW, Voinnet O. Antiviral immunity directed by small RNAs. Cell. 2007;130:413–26.
Miozzi L, Gambino G, Burgyan J, Pantaleo V. Genome-wide identification of viral and host transcripts targeted by viral siRNAs in Vitis vinifera. Mol Plant Pathol. 2013;14:30–43.
Schuck J, Gursinsky T, Pantaleo V, Burgyán J, Behrens SE. AGO/RISC mediated antiviral RNA silencing in a plant in vitro system. Nucleic Acids Res. 2013;41:5090–103.
Baumberger N, Baulcombe DC. Arabidopsis ARGONAUTE1 is an RNA slicer that selectively recruits microRNAs and short interfering RNAs. Proc Natl Acad Sci U S A. 2005;102:11928–33.
Llave C. Virus-derived small interfering RNAs at the core of plant-virus interactions. Trends Plant Sci. 2010;15:701–7.
Qu F, Ye XH, Hou GC, Sato S, Clemente TE, Morris TJ. RDR6 has a broad-spectrum but temperature-dependent antiviral defense role in Nicotiana benthamiana. J Virol. 2005;24:15209–17.
Shimura H, Pantaleo V, Ishihara T, Myojo N, Inaba J, Sueda K, et al. Viral satellite RNA induces yellow symptoms on tobacco by targeting a gene involved in chlorophyll biosynthesis using the RNA silencing machinery. PLoS Pathog. 2011;7:e1002021.
Smith NA, Eamens AL, Wang MB. Viral small interfering RNAs target host genes to mediate disease symptoms in plants. PLoS Pathog. 2011;7:e1002022.
Xia ZH, Peng J, Li YQ, Chen L, Li S, Zhou T, et al. Characterization of small interfering RNA derived from Sugarcane mosaic virus in infected maize plants by deep sequencing. PLoS One. 2014;9:e7013.
Guo C, Li L, Wang XF, Liang C. Alterations in siRNA and miRNA expression profiles detected by deep sequencing of transgenic rice with siRNA-mediated viral resistance. PLoS One. 2015;10:e0116175.
Zhang YJ, Wang W, Chen J, Liu JB, Xi MX, Shen FF. Identification of miRNAs and their targets in cotton inoculated with Verticillium dahlia by high-throughput sequencing and degradome analysis. Int J Mol Sci. 2015;16:14749–68.
Vaucheret H. Plant ARGONAUTES. Trends Plant Sci. 2008;13:350–8.
Del Toro FJ, Aguilar E, Hernández-Walias FJ, Tenllado F, Chung BN, Canto T. High temperature, high ambient CO2 affect the interactions between three positive-sense RNA viruses and a compatible host differentially, but not their silencing suppression efficiencies. PLoS One. 2015;10:e0136062.
Andika IB, Sun LY, Xiang R, Li JM, Chen JP. Root-Specific Role for Nicotiana benthamiana RDR6 in the Inhibition of Chinese wheat mosaic virus accumulation at higher temperatures. Mol Plant Microbe Interact. 2013;26:1165–75.
Garcia-Ruiz H, Takeda A, Chapman EJ, Sullivan CM, Fahlgren N, Brempelis KJ, et al. Arabidopsis RNA-dependent RNA polymerases and Dicer-like proteins in antiviral defense and small interfering RNA biogenesis during Turnip mosaic virus infection. Plant Cell. 2010;22:481–96.
Zhang X, Singh J, Li D, Qu F. Temperature-dependent survival of Turnip crinkle virus-infected Arabidopsis plants relies on an RNA silencing-based defense that requires DCL2, AGO2, and HEN1. J Virol. 2012;86:6847–68.
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:116–27.
Jaubert MJ, Bhattacharjee S, Mell AF, Perry KL, Moffett P. ARGONAUTE2 mediates RNA-silencing antiviral defenses against Potato virus X in Arabidopsis. Plant Physiol. 2011;156:1556–64.
Scholthof HB, Alvarado VY, Vega-Arreguin JC, Ciomperlik J, Odokonyero D, Brosseau C, et al. Identification of an ARGONAUTE for antiviral RNA silencing in Nicotiana benthamiana. Plant Physiol. 2011;156:1548–55.
Schwach F, Vaistij FE, Jones L, Baulcombe DC. An RNA-dependent RNA polymerase prevents meristem invasion by Potato virus X and is required for the activity but not the production of a systemic silencing signal. Plant Physiol. 2005;138:1842–52.
Hong W, Qian D, Sun RH, Jiang L, Wang Y, Wei CH, et al. OsRDR6 plays role in host defense against double-stranded RNA virus. Rice dwarf phytoreovirus Scientific Rep. 2015;5:11324.
Qi XP, Bao FS, Xie ZX. Small RNA deep sequencing reveals roles for Arabidophsis thaliana RNA-dependent RNA polymerases in viral siRNA biogenesis. PLoS One. 2009;4:e4971.
Qu F, Ye X, Morri 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:14732–7.
Jiang L, Qian D, Zheng H, Meng LY, Chen J, Le WJ, et al. RNA-dependent RNA polymerase 6 of rice (Oryza sativa) plays role in host defense against negative-strand RNA virus, Rice stripe virus. Virus Res. 2012;163:512–9.
Adams MJ, Candresse T, Hammond J, Kreuze JF, Martelli GP, Namba S, et al. Family Betaflexiviridae. In: King AMQ, Lefkowitz E, Adams MJ, Carstens EB, editors. Virus Taxonomy: Ninth Report of the International Committee on Taxonomy of Viruses. San Diego, USA: Elsevier; 2012. p. 920–41.
Yoshikawa N, Sasaki E, Kat M, Takahashi T. The nucleotide sequence of Apple stems grooving Capillovirus genome. Virology. 1992;191:98–105.
Tatineni S, Afunian MR, Gowda S, Hilf ME, Bar-Joseph M, Dawon WO. Characterization of the 5′- and 3′-terminal subgenomic RNAs produced by Capillovirus: Evidence for a CP subgenomic RNA. Virology. 2009;385:521–8.
Komatsu K, Hirata H, Fukagawa T, Yamaji Y, Okano Y, Ishikawa K, et al. Infection of capilloviruses requires subgenomic RNAs whose transcription is controlled by promoter-like sequences conserved among flexiviruses. Virus Res. 2012;167:8–15.
Liebenberg A, Moury B, Sabath N, Hell R, Kappis A, Jarausch W, et al. Molecular evolution of the genomic RNA of apple stem grooving capillovirus. J Mol Evol. 2012;75:92–101.
Chen H, Chen S, Li Y, Ye T, Hao L, Fan Z, et al. Phylogentic analysis and recombination events in full genome sequences of apple stem grooving virus. Acta Virological. 2014;58:309–16.
Inouye N, Maeda T, Mitsuhata K. Citrus tatter leaf virus isolated from lily. Ann Phytopathol Soc Jap. 1979;45:712–20.
Nemeth M.Virus, myvoplasma and rickettsia diseases of fruit trees. Akademiai Kiado. Budapest. 1986. p 841
Clover GRG, Pearson MN, Elliot DR, Tang Z, Smales TE, Alexander BJR. Characterization of a strain of apple stem grooving viruses in Actinidia chinesis from China. Plant Pathol. 2003;52:371–8.
Shim H, Min Y, Hong S, Kwon M, Kim H, Choi Y, et al. Nucleotide sequences of a Korean isolate of Apple stem grooving virus associated with black necrotic leaf spot disease on pear (Pyrus pyrofolia). Moleculars cells. 2004;18:192–9.
Chen SY, Ye T, Hao L, Chen H, Wang SJ, Fan ZF, et al. Infection of apple by Apple stem grooving virus leads to extensive alterations in gene expression patterns but no disease symptoms. PLoS One. 2014;9:e95239.
Wang GP, Hong N, Zhang ZP. Identification of virus species in pears cultivated in northern China. China Fruits. 1994;4:1–4.
Zheng YY, Wang GP, Hong N. The biological characteristics and molecular identification of some Apple stem grooving virus isolates. Acta Phytopathol Sinica. 2005;32:266–70.
Wang QC, Cuellar WJ, Rajamaki ML, Hirata Y, Valkone JT. Combined thermotherapy and cryotherapy for efficient virus eradication: Relation of virus distribution, subcellular changes, cell survival and viral RNA degradation in shoot tips. Mol Plant Pathol. 2008;9:237–50.
Hu GJ, Hong N, Wang LP, Hu HJ, Wang GP. Efficacy of virus elimination from in vitro-cultured sand pear (Pyrus pyrifolia) by chemotherapy combined with thermotherapy. Crop Prot. 2012;37:20–5.
Szittya G, Silhavy D, Molnar A, Havelda Z, Lovas AL, Banfalvi Z, et al. Low temperature inhibits RNA silencing-mediated defense by the control of siRNA generation. EMBO J. 2003;22:633–40.
Chellappan P, Vanitharani R, Ogbe F, Fauque CM. Effect of temperature on geminivirus-induced RNA-silencing in plants. Plant Physiol. 2005;138:1828–41.
Velázquez K, Renovell A, Comellas M, Serra P, García ML, Pina JA, et al. Effect of temperature on RNA silencing of negative-strand RNA plant virus: Citrus psorosis virus. Plant Pathol. 2010;59:982–90.
Mochizuk T, Ohki ST. Shoot meristem tissues of tobacco inoculated with Cucumber mosaic virus is infected with the virus and subsequently recovers from infection by RNA silencing. Gen Plant Pathol. 2004;70:363–6.
Juan L, Zhang XJ, Zhang FP, Hong N, Wang GP, Wang A, et al. Identification and characterization of microRNAs from in vitro-grown pear shoots infected with Apple stem grooving virus in response to high temperature using small RNA sequencing. BMC Genomics. 2015;16:945.
Zerbino DR, Birney E. Algorithms for de novo short read assembly using de Bruijn graphs. Genome Res. 2008;18:821–9.
Langmead B, Trapnel C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25.
Wu J, Wang ZW, Shi ZB, Zhang S, Ming R, Zhu SL, et al. The genome of the pear (Pyrus bretschneideri Rehd.). Genome Res. 2013;23:396–408.
Rychlik W. OLIGO 7 primer analysis software. Methods Mol Biol. 2007;402:35–60.
Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, et al. Multiple sequence alignment with the clustal series of programs. Nucleic Acids Res. 2003;31:3497–500.
Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. MEGA6: Molecular Evolutionary Genetics Analysis version 6.0. Mol Biol Evol. 2013;30:2725–9.
Angelini E, Clair D, Forgo M, Bertaccini A, Boudon-Padieu E. Flavescence dorée in France and Italy-Occurrence of closely related phytoplasma isolates and their near relationships to palatinate grapevine yellows and an alder phytoplasma. Vitis. 2001;40:79–86.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real–time quantitative PCR and the 2-ΔΔCT method. Methods. 2001;25:402–8.
Chellappan P, Vanitharani R, Pita J, Fauquet CM. Short interfering RNA accumulation correlates with host recovery in DNA virus-infected hosts, and gene silencing targets specific viral sequences. J Virol. 2004;78:7465–77.
Hagen C, Rojas MR, Kon T, Gilbertson RL. Recovery from Cucurbit leaf crumple virus (Family Geminiviridae, genus Begomovirus) infection is an adaptive antiviral response associated with changes in viral small RNAs. Phytopathology. 2008;98:1029–37.
Wang XB, Wu Q, Ito T, Cillo F, Li WX, Chen XM, et al. RNAi-mediated viral immunity requires amplification of virus-derived siRNAs in Arabidopsis thaliana. Proc Natl Acad Sci U S A. 2010;107:484–9.
Yang XL, Wang Y, Guo W, Xie Y, Xie Q, Fan LJ, et al. Characterization of Small Interfering RNAs Derived from the Geminivirus/Betasatellite Complex Using Deep Sequencing. PLoS One. 2011;6:e1692.
Xu Y, Huang LZ, Fu S, Wu JX, Zhou XP. Population diversity of Rice Stripe Virus-derived siRNAs in three different hosts and RNAi-Based antiviral immunity in Laodelphgax striatellus. PLoS One. 2012;9:e46238.
Czosnek H, Eybishtz A, Sade D, Gorovits R, Sobol I, Bejarano E, et al. Discovering Host Genes Involved in the Infection by the Tomato Yellow Leaf Curl Virus Complex and in the Establishment of Resistance to the Virus Using Tobacco Rattle Virus-based Post Transcriptional Gene Silencing. Viruses. 2013;5:998–1022.
Dey KK, Borth WB, Melzer MJ, Wang ML, Hu JS. Analysis of Pineapple mealybug wilt associated virus-1 and −2 for potential RNA silencing suppressors and pathogenicity factors. Viruses. 2015;7:969–95.
Alabi OJ, Zheng Y, Jagadeeswaran G, Sunkar R, Nadiu RA. High-throughput sequence analysis of small RNAs in grapevine (Vitis vinifera L.) affected by grapevine leafroll disese. Mol Plant Pathol. 2012;13:1060–76.
Visser M, Maree HJ, Rees DJ, Burger JT. High-throughput sequencing reveals small RNAs involved in ASGV infection. BMC Genomics. 2014;15:568.
Donaire L, Barajas D, Martínez-García B, Martínez-Priego L, Pagán I, Llave C. Structural and genetic requirement for the biogenesis of Tobacco rattle virus-derived small interfering RNAs. J Virol. 2008;82:5167–77.
Ribeiro SG, Lohuis H, Goldbach R, Prins M. Tomato chlorotic mottle virus is a target of RNA silencing but the presence of specific short interfering RNAs does not guarantee resistance in transgenic plants. J Virol. 2007;81:1563–73.
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:19593–8.
Takeda A, Iwasaki S, Watanabe T, Utumi 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.
Qu F. Antiviral role of plant-encoded RNA-dependent RNA polymerases revisited with deep sequencing of small interfering RNAs of virus origin. Mol Plant Microbe Interact. 2010;23:1248–52.
This work was supported by the National Natural Science Foundation of China (No.31201488) and (No.31471862), and the Earmarked Fund for Pear Modern Agro-industry Technology Research System (CARS-29-10).
LJ, ZXJ and YYK designed, carried out the experiments and analyzed the data. WLP wrote the draft manuscript. WGP, HN and WAM revised the manuscript. All authors read and approved the submitted version manuscript.
The authors declare that they have no competing interests.
1.2% agarose gel electrophoresis of RT-PCR products of ASGV-Js2 isolate. M: Marker II (TIANGEN Biotech, Beijing Co., Ltd.), CK+: P. pyrifolia cv. ‘HuangHua’ shoots as positive control, CK-: P .pyrifolia cv. ‘YuanHuang’ shoots as negative control. Lanes 1-6 show amplification products from virus-free plants (A) and ASGV-infected plants (B). (DOC 532 kb)
Oligonucleotide primers used for amplification of the ASGV-Js2 genome (DOC 45 kb)
Oligonucleotide primer sequences and amplicon characteristics for the PpDCL2,4, PpAGO1,2,4, and PpRDR1 genes. (DOC 35 kb)
Oligonucleotide primers used for qRT-PCR expression analysis of vsiRNAs, P. pyrifolia mRNAs, and the ASGV-Js2 mp gene from in vitro-grown shoots of P. pyrifolia. (DOC 42 kb)
Pairwise sequence similarities comparing the ASGV-Js2 full-length genome with 17 ASGV isolates sequences from GenBank. (DOC 66 kb)
About this article
- Pyrus pyrifolia
- Gene silencing
- Apple stem grooving virus
- High temperature
- Virus-derived small interfering RNA (vsiRNA)
- Argonaute (AGO)
- RNA dependent RNA polymerase (RDRs)
- Dicer-like (DCL)