Abstract
Thrips-transmitted tospoviruses are economically important viruses affecting a wide range of field and horticultural crops worldwide. Tomato spotted wilt virus (TSWV) is the type member of the Tospovirus genus with a broad host range of more than 900 plant species. Interactions between these viruses and their plant hosts and insect vectors via RNAi pathways are likely a key determinant of pathogenicity. The current investigation, for the first time, compares biogenesis of small RNAs between the plant host and insect vector in the presence or absence of TSWV. Unique viral small interfering RNA (vsiRNA) profiles are evident for Arachis hypogaea (peanut) and Frankliniella fusca (thrips vector) following infection with TSWV. Differences between vsiRNA profiles for these plant and insect species, such as the relative abundance of 21 and 22 nt vsiRNAs and locations of alignment hotspots, reflect the diverse siRNA biosynthesis pathways of their respective kingdoms. The presence of unique vsiRNAs in F. fusca samples indicates that vsiRNA generation takes place within the thrips, and not solely through uptake via feeding on vsiRNAs produced in infected A. hypogaea. The study also shows key vsiRNA profile differences for TSWV among plant families, which are evident in the case of A. hypogaea, a legume, and members of Solanaceae (S. lycopersicum and Nicotiana benthamiana). Distinctively, overall small RNA (sRNA) biogenesis in A. hypogaea is markedly affected with an absence of the 24 nt sRNAs in TSWV-infected plants, possibly leading to wide-spread molecular and phenotypic perturbations specific to this species. These findings add significant information on the host–virus–vector interaction in terms of RNAi pathways and may lead to better crop and vector specific control strategies.
Introduction
Thrips and tospoviruses continue to be one of the biggest production constraints to numerous crops in both developed and developing parts of the world, with an economic cost in excess of US$1 billion annually (Pappu et al., 2009). Tospoviruses have been reported to infect more than 1090 species in 85 families of monocots and dicots (Mandal, 2012). Several informative reviews on different aspects of tospoviruses, the only plant-infecting viruses in the family Bunyaviridae, and their thrips vectors have been published in recent years (Jones, 2005; Whitfield et al., 2005; Persley et al., 2006; Lebas and Ochoa-Corona, 2007; Tsompana and Moyer, 2008; Pappu et al., 2009; Riley et al., 2011). Tomato spotted wilt virus (TSWV) is a member of the Tospovirus genus and is ranked second among the top ten economically important viruses in the world (Scholthof et al., 2011). It is one of the most important pathogens affecting peanut production in the south-eastern United States (Culbreath and Srinivasan, 2011). At least 10 species of thrips (Thysanoptera:Thripidae) have been implicated in TSWV transmission in a persistent and propagative manner; that is, the virus also replicates in the insect vector (Ullman et al., 1993; Whitfield et al., 2005; Pappu et al., 2009). Frankliniella fusca Hinds (tobacco thrips) is considered to be an important vector of TSWV on peanut (Riley et al., 2011). Use of resistant cultivars along with the integration of cultural and chemical tactics including control of the insect vector through insecticides are the primary methods used to manage TSWV. Understanding the interactions between the host, vector and virus at a molecular level is a crucial component in designing new management strategies and remains a highly active area of research.
RNA interference (RNAi), a key conserved mechanism in eukaryotes, plays a major role in innate defense against viruses and transposons (Baulcombe, 2004; Pelaez and Sanchez, 2013). The trigger for RNAi is double-stranded RNA (dsRNA), which may be generated during viral replication in the host, be it plants or insects. DsRNA is processed into 21–24 nt small interfering RNAs (siRNAs) in a complex series of interactions involving Dicer-like proteins (DCLs), Argonauts (AGOs), and RNA-dependent RNA polymerases (RdRps) (Ding and Voinnet, 2007). Another form of RNAi produces 21–24 nt microRNAs (miRNAs) from imperfect hairpin structures known as pre-miRNAs. MiRNAs often target mRNA transcripts and thus have a role in gene regulation (Mallory and Vaucheret, 2004). Together, miRNAs and siRNAs molecules are classified as small RNAs (sRNAs).
Virus specific siRNAs (vsiRNAs) are generated via RNAi pathways in virus-infected plants, and have been identified in diverse eukaryotic species (Ding, 2010; Wu et al., 2010). The presence of abundant vsiRNA hotspots across the viral genome reveal the likely initiation and phased loci of the host’s RNAi targeting mechanisms. Previous reports have demonstrated that vsiRNA hotspots are not uniformly distributed across viral genome segments and gene-coding regions, suggesting that the host has a preference for targeting specific regions, with consequences for both viral pathogenicity and host specificity (Donaire et al., 2008; Wang et al., 2010; Mitter et al., 2013).
Due to replication of TSWV in the thrips vector, infection-induced changes in sRNA abundance along with the location and magnitude of vsiRNA hotspots relative to the plant host are also of interest. Reports of TSWV vsiRNA hotspot position and magnitude in insect vector-derived samples are largely absent from the literature, though studies of vsiRNA hotspots related to other viruses and insect hosts have been carried out (Myles et al., 2008; Mueller et al., 2010; Hess et al., 2011). Similar to the plant response to viral infection, inhibition of viral replication via the insect’s dsRNA-triggered RNAi machinery is well documented in Drosophila melanogaster (Meigen; Kemp et al., 2013), Anopheles gambiae (Giles; Keene et al., 2004) and Aedes aegypti (L.; Campbell et al., 2008). A comparative study of Rice stripe virus-infected plant hosts Oryza sativa (L.) and Nicotiana benthamiana (Domin) and insect vector Laodelphax striatellus (Fallen) demonstrated the similarities and differences in both overall sRNA abundance and vsiRNA-profiles in these plant and insect species (Xu et al., 2012). Differences in the biogenesis of sRNAs were readily evident, with 21 and 24 nt sRNAs most prevalent in the plant species, and 22 and 26–27 nt sRNAs in the insect vector. However, vsiRNAs aligning to the viral genome were dominated by the 21 nt sRNAs in the plant samples, and 22 nt sRNAs in the insect derived sample. Significant variation in the location and magnitude of vsiRNA hotspots between the hosts was readily evident, which may be more broadly indicative of virus–vector–host relationships.
The use of next-generation sequencing has allowed for comprehensive analysis of sRNA populations in virus-infected samples, both generally and the specific sRNA subset directed against the virus (Szittya et al., 2010; Xu et al., 2012; Miozzi et al., 2013; Mitter et al., 2013). In this study, deep sequencing of the small RNA component of samples derived from TSWV-infected and non-infected Arachis hypogaea (L.) host plants and the F. fusca (Hinds) vector was performed. In silico analysis of these data demonstrates both common and unique features among the plant host and insect vector, as well as between A. hypogaea, a legume, and members of the Solanaceae infected with the same viral isolate.
Materials and Methods
Maintenance of Non-infected and TSWV-Infected A. hypogaea Plants
Arachis hypogaea seedlings (cv. Georgia Green) were germinated and maintained in a growth room at 25–30°C, 40–50% relative humidity (rh) and L14:D10 photoperiod. Subsequently, sprouted seeds were transplanted and maintained in thrips-proof cages in a greenhouse at 25–30°C, 80–90% rh, and L14:D10 photoperiod. TSWV-infected volunteer A. hypogaea plants were collected from Belflower Farm (Coastal Plain Experimental Station, Tifton, GA, USA) and maintained in a greenhouse under the same environmental conditions. Leaflets with TSWV symptoms were collected, placed in a Munger cage, into which non-viruliferous F. fusca female adults were released (Munger, 1942). After the next generation adults emerged, potentially viruliferous adults obtained from the cages were used to inoculate prepared A. hypogaea plants. Ten female adults (up to 2 days old) were released on 1-week-old plants that had been dusted with 0.05 g of pine (Pinus taeda L.) pollen. Each plant was enclosed in a Mylar film (Grafix, Cleveland, PA, USA) cage with a copper mesh top. Plants were maintained in thrip-proof cages in the greenhouse as described above. Subsequently 3 weeks post thrips inoculation, infection status of the plants was determined by real-time quantitative reverse transcription polymerase chain reaction (qRT-PCR) using TSWV N-gene-specific forward and reverse primers, 5′-GCTTCCCACCC TTTGATTC-3′ and 5′-ATAGCCAAGACAACACTGATC-3′, respectively (Rotenberg et al., 2009). Plasmids with TSWV N-gene inserts were used as external copy number standards.
Maintenance of Non-viruliferous and Viruliferous F. fusca
Adult F. fusca were initially collected from A. hypogaea blooms from Belflower farm. Subsequently, non-viruliferous and potentially viruliferous F. fusca colonies were separately reared on non-infected and TSWV-infected A. hypogaea leaflets, respectively, in Munger cages under the conditions described earlier. F. fusca reared for an entire generation (egg to adult) on TSWV-infected A. hypogaea leaflets alone were considered as potentially viruliferous. The TSWV infection status of subsamples from potentially viruliferous and non-viruliferous thrips used for this study was determined by qRT-PCR using the primers and conditions as described above. In both non-viruliferous and potentially viruliferous colonies, a majority of thrips (>95%) were females. Rotenberg et al. (2009) revealed that TSWV nucleocapsid RNA copies were two to three times greater in female thrips [Frankliniella occidentalis (Pergande)] when compared with male thrips. Although a different thrips species was the focus of this study [F. fusca (Hinds)], only females were used.
Total RNA Extraction and Deep Sequencing
Approximately 160 female adults (up to 2 days old) were separately pooled from viruliferous and non-viruliferous thrips colonies, and total RNA was extracted. Similarly, leaf tissue from three TSWV-infected or non-infected plants (about 8 weeks old) was used for plant total RNA extraction. A total of two independent biological replicates were included for each F. fusca treatment and three independent biological replicates for each A. hypogaea treatment. Total RNA from thrips and plant samples were extracted using Direct-zol RNA Miniprep kit (Zymo Research, Irvine, CA, USA) using the manufacturer’s protocol. Library construction and sRNA sequencing (50 bp, single end) were carried out by the Beijing Genomic Institute using the Illumina HiSeq sequencing system. Samples were multiplexed and sequenced over two lanes. All sequence files are publicly accessible via the NCBI Sequence Read Archive (uninfected A. hypogaea: SRR3727343, SRR3727344, SRR3727345; TSWV-infected A. hypogaea: SRR3727346, SRR3727347, SRR3727348; uninfected F. fusca: SRR3727364, SRR3727365; TSWV-infected F. fusca: SRR3727366, SRR3727368).
Bioinformatic Analysis of sRNA Sequence Files
Alignments of sRNA reads were carried out against the following TSWV RNA segment accessions: L RNA - NC_002052.1; M RNA - AY744483.1; S RNA - AY744475.1. For A. hypogaea sample alignments, reference sequences were obtained for A. ipaensis, a sequenced progenitor of A. hypogaea. These references included an annotated transcriptome (Araip.K30076.a1.M1) and mobile element data set (BB051914; Dash et al., 2016). Previously published raw sRNA sequence files from TSWV-infected S. lycopersicum (SRR947065) and N. benthamiana (SRR947064) were re-processed and used in this study for comparative purposes (Mitter et al., 2013). In this previous study, S. lycopersicum and N. benthamiana plants were mechanically inoculated with TSWV 3–4 weeks post-germination and leaf samples from TSWV positive plants were collected 17 days post inoculation.
Briefly, adapter sequences were removed from 50 bp single-end reads, which were then collapsed to identical sequences with maintained counts using the FASTX-Toolkit package. Only reads with lengths between 18 and 32 nt were retained for alignment normalization. Alignment of collapsed reads to reference sequences was performed using the SCRAM software package which allowed for exact matches to the reference sequence only (Fletcher, 2016). Normalization of aligned read count at each position was calculated as reads aligned per million reads between length 18 and 32 nt in the collapsed read file, allowing for robust quantitative comparison between alignments. Boxplots and sRNA length by count profiles were generated in R Studio using custom R scripts. Graphical representations of alignment profiles were produced either directly from the Python alignment software or via the Circos software package using tabulated output (Krzywinski et al., 2009). Scatter plots and Venn diagrams were generated in R Studio, and Euler diagrams in eulerAPE.
Results And Discussion
TSWV infection induced a wide range of symptoms in A. hypogaea plants in the greenhouse. Typical symptoms included chlorotic mottling on the adaxial leaflet surface (Figure 1A), concentric ringspots with green centers (Figure 1B), drooping of leaf petioles, and drying of terminals (Figure 1C). In contrast, viruliferous F. fusca displayed no obvious symptoms, with their external appearance indistinguishable from that of non-viruliferous F. fusca (Figure 1D).
FIGURE 1
sRNA Profiles in A. hypogaea and F. fusca as a Result of TSWV Infection
The size distributions of sRNAs in TSWV-infected and uninfected samples are shown in Figure 2A. The uninfected samples represent only the host sRNAs whereas the infected samples include total sRNAs (host and vsiRNAs). We have also aligned only vsiRNAs in the infected samples. The categorisation of 20–24 nt total sRNAs into host sRNAs (non-vsiRNAs) and vsiRNAs is presented in Supplementary Table S1. The bimodal size distributions of sRNA reads in healthy A. hypogaea and F. fusca are typical of their origin, with 21 and 24 nt sRNAs dominating in the plant host and 22 nt and 27–28 nt sRNAs most prevalent in the insect vector (Figure 2A). The relative abundance of 21 and 24 nt peaks apparent in the healthy A. hypogaea profile closely resembled that of N. benthamiana, a Solanaceae member, in a previous report (Mitter et al., 2013). Similarly, the distribution of sRNA reads in the infected and healthy F. fusca samples matched those in profiles of insect species as diverse as the silkworm Bombyx mori (L.; Jagadeeswaran et al., 2010) and migratory locust Locusta migratoria (Wei et al., 2009). L. striatellus infected with Rice stripe virus also displayed a bipartite profile, though a peak was evident at 26–27 nt rather than 27–28 nt, and contrastingly, the abundance of the sRNAs making up the larger sized peak was lower than that of the 22 nt sRNAs (Xu et al., 2012). The distinct profiles of the plant host and insect vector likely reflect the differences in sRNA biogenesis and function among the organisms. The use of female thrips may have influenced the prevalence of vsiRNAs observed, as they could be more abundant in female thrips than in male thrips (Rotenberg et al., 2009).
FIGURE 2
Remarkably, the impact of TSWV infection on the size distribution of total sRNAs in A. hypogaea is profound, with a significant reduction in total 24 nt sRNAs after infection (Figures 2A,B). Compared to uninfected plants, there is also an increase in 22 nt sRNAs, and to a lesser extent, 20 nt sRNAs (Figure 2B). Such a change in relative abundance of 22and 24 nt sRNAs due to TSWV infection is in contrast to an increase in 24 nt sRNAs and decrease in 22 nt sRNAs in other reported species, such as Bamboo mosaic virus-inoculated N. benthamiana plants (Lin et al., 2010). A slight decrease in 24 nt sRNAs in cotton plants infected with Cotton leaf roll dwarf virus indicates a similar pattern, but at a much reduced magnitude compared to TSWV-infected A. hypogaea (Silva et al., 2011). These data suggest that infection by TSWV significantly alters the productive balance of the plant’s sRNA biogenesis machinery. In contrast, infection of the F. fusca vector by TSWV results in minimal change in the relative abundance of the various sRNA discrete sizes (Figure 2A).
Alterations in the Endogenous sRNA Landscape of A. hypogaea Due to TSWV Infection
Due to their role in maintenance of genome stability via repression of transposable elements, the reduction in the number of 24 nt sRNAs in TSWV-infected A. hypogaea could have significant physiological consequences. Additionally, an increase in 22 nt sRNAs may lead to other further changes in endogenous sRNA abundance, as these sRNAs trigger biogenesis of phased sRNAs spanning complementary RNA sequences, along with silencing of the affected genes (McHale et al., 2013). Widespread perturbation of gene expression in infected plants, either directly or indirectly, may be a fundamental factor in the infected plant phenotype.
Analysis of sRNA alignments to repetitive and transposable elements within the A. hypogaea genome demonstrated a clear reduction in sRNAs derived from these elements in TSWV-infected plants (Figure 3). This shift likely has wide-ranging impacts, with transposition events potentially happening at a greater frequency, and the expression of nearby endogenous genes altered. In support of this prediction, previous studies have shown that the disruption of biogenesis of 24 nt sRNAs in plants had an extensive impact on the degree of DNA methylation along with the activity of targeted transposable elements, and the expression of neighboring genes (Jia et al., 2009; Groszmann et al., 2011).
FIGURE 3
Alignments to protein-coding sequences also demonstrated clear changes in the abundance of 21, 22, and 24 nt sRNAs targeting these sequences (Figure 4). Interestingly, an increase in the variance of coverage of coding sequences by 21, 22, and 24 nt sRNAs is particularly evident, pointing toward possible broad up- and down-regulation of gene expression in infected leaf tissue. A previous study has suggested that the TSWV-derived Non-structural protein (NSs) binds directly to intermediates of sRNA processing and potentially suppresses loading of sRNA into AGO complexes required for degradation or translational inhibition of complementary RNAs (Schnettler et al., 2010). Therefore, the disruption of sRNA levels observed in this study and the potential for NSs action on plant-generated sRNAs may be a means for perturbation of target RNA transcript levels and consequential physiological and developmental changes.
FIGURE 4
This report adds further evidence for the species-specific impact virus inoculation has on the dynamics of sRNA processing. The exact mechanism driving this shift in sRNA abundance and the reason it differs between species remains to be fully elucidated. DCL2, DCL3, and DCL4 homologs of plants all cleave dsRNA substrates, either endogenous or exogenous/viral derived dsRNAs, into 22, 24, and 21 nt sRNAs, respectively (Henderson et al., 2006; Bologna and Voinnet, 2014). Studies on Arabidopsis dcl mutants show DCL4 and DCL2 are partly redundant and act hierarchically in processing dsRNA substrates associated with post-transcriptional gene silencing (Parent et al., 2015). In contrast, DCL3 is accepted as the primary enzyme for generation of 24 nt siRNAs from dsRNAs produced by the Pol IV-RDR2 complex and associated with transcriptional silencing of repetitive elements and transposons (Gasciolli et al., 2005; Pontes et al., 2006; Wierzbicki et al., 2009; Johnson et al., 2014). Our data demonstrated that TSWV infection in A. hypogaea, but not the two Solanaceae species (Mitter et al., 2013), triggered a genome-wide depletion in 24 nt siRNAs associated with repetitive elements and transposons (Figure 3). This suggests that in A. hypogaea there was a TSWV-dependent suppression in expression levels, subcellular localization and/or the enzyme function of DCL3, or another upstream component of the transcriptional gene silencing pathway, such as Pol IV or RDR2 (Wierzbicki et al., 2009). Thus, inhibition of DCL3 or another component of the transcriptional gene silencing pathway is the most likely explanation for the shift in production of 24 nt siRNAs to 22 nt siRNAs in TSWV-infected A. hypogaea. Such a reaction to TSWV infection on the accumulation of 24 nt siRNAs clearly varies between plant species. Further expression analysis on components of the transcriptional gene silencing pathway may shed light on differential effects on TSWV infection in different host species.
Divergence in Relative Abundance of A. hypogaea and F. fusca vsiRNAs Aligned to Each TSWV Genome Segment
The portion of total sRNAs that align to the infecting viral genome has been shown to vary greatly depending on the virus and the plant host (Dehaan et al., 1991; Bolduc et al., 2010; Mitter et al., 2013; Hong et al., 2015). Approximately 16% of all reads in each infected A. hypogaea replicate aligned to the TSWV genome, while just 0.4% aligned to the virus in the infected F. fusca vector (Table 1). The abundance of vsiRNAs aligning to the TSWV genome in A. hypogaea is high compared to S. lycopersicum and N. benthamiana infected with the same viral isolate (6% and 0.02%, respectively), though in the range of other host–virus combinations such as peach infected with Peach latent mosaic viroid (11.2 %; Bolduc et al., 2010) and rice infected with Rice dwarf phytoreovirus (∼22%; Hong et al., 2015).
Table 1
| Sample | Genome Size (Mb/1C) | Total Reads (average if multiple replicates) | vsiRNA reads per million reads (21–22 nt) | vsiRNA percentage of total (%) |
|---|---|---|---|---|
| A. hypogaea | ∼2800 | 9179622 | 165298 | 16.53% |
| F. fusca | ∼400 | 21915376 | 4197 | 0.42% |
| S. lycopersicum | ∼950 | 19908563 | 62952 | 6.30% |
| N. benthamiana | ∼3000 | 18780917 | 145 | 0.01% |
Effect of the host organism on Tomato spotted wilt virus (TSWV)-specific vsiRNA (viral small interfering RNA) levels.
The mean percentage of total sRNA reads corresponding to the TSWV genome for infected Arachis hypogaea and Frankliniella fusca samples is displayed. S. lycopersicum and Nicotiana benthamiana alignment percentages are shown for comparison. All four species were infected with the same isolate of TSWV. The vsiRNA load in the A. hypogaea host is ∼30 times that of the F. fusca vector.
Tospoviruses are characterized by a tripartite genome structure consisting of large (L), medium (M) and small (S) single-stranded RNA segments with sense and ambi-sense coding frames. The L segment consists of a large open reading frame (ORF) that encodes RdRp (RNA-dependent RNA polymerase), which is required for viral replication (Dehaan et al., 1991). ORFs for Gn/Gc (precursor for the N-terminal G2 and C-terminal G1proteins) and NSm (non-structural protein) occur in the M segment. NSm is required for cell-to-cell viral movement (Soellick et al., 2000), and Gn/Gc for synthesis of glycoproteins present in the viruses’ lipid envelope (Snippe et al., 2007). The S segment encodes the N (nucleocapsid) gene, which encapsidates the RNA segments (Dehaan et al., 1990), and NSs (non-structural/suppressor of silencing protein), which is of particular importance in the defense against the plant’s RNAi defense strategy. Indeed, the degree of NSs protein accumulation in plant tissues has been associated with the severity of TSWV infections (Takeda et al., 2002). NSs can bind to both dsRNA and vsiRNA targets, allowing it to block antiviral RNA silencing before and after DCL-mediated dsRNA cleavage (Schnettler et al., 2010). Importantly, NSs has been shown to be effective in not only plant hosts, but also in insect vectors (Margaria et al., 2014), possibly through its affinity for longer dsRNA forms. Alignment of reads from infected plant samples to individual genome segments of TSWV has been previously reported, also with dissimilar results depending on the plant species and viral isolate (Mitter et al., 2013; Margaria et al., 2015).
In this study, TSWV-aligned vsiRNAs from A. hypogaea demonstrated a distribution distinct from vsiRNAs from the Solanaceae hosts, with a generally high portion aligned to the L segment and variable percentages aligning to the M and S segments, depending on the replicate. In contrast, F. fusca-derived vsiRNAs more closely matched the aligned vsiRNA distribution from S. lycopersicum and N. benthamiana (Figure 5; Table 2).
FIGURE 5
Table 2
| Genome segment | Expected percentage alignment based on size | Mean percentage of 21 – 22 nt vsiRNAs aligned – A. hypogaea | Mean percentage of 21 – 22 nt vsiRNAs aligned – F. fusca |
|---|---|---|---|
| L RNA | 54% | 39% | 24% |
| M RNA | 29% | 25% | 36% |
| S RNA | 18% | 36% | 40% |
Mean percentage of aligned reads for each TSWV segment along with the expected alignment percentage based on the size of the segment.
The abundance of aligned vsiRNAs in both A. hypogaea and F. fusca samples diverge from these expected percentages, demonstrating a segment-specific targeting preference.
The use of multiple fully independent biological replicates for A. hypogaea sRNA sequencing provides a detailed picture of the strand polarity of alignments to individual TSWV gene sequences (Figure 6). Significant variation amongst the A. hypogaea replicates indicates that there is little preference for a particular alignment orientation. This is in contrast to earlier reports for TSWV, where a greater abundance of vsiRNAs generally aligned to the sense strand of the virus in Solanaceae hosts (Mitter et al., 2013). The possibility exists that greater replication of sequencing for other TSWV-infected plant species may demonstrate similar variation. F. fusca, on the other hand, shows a sense preference for the Gn/Gc, RdRp and Nc genes, and an antisense preference for the NSm and NSs genes.
FIGURE 6
Due to vsiRNA being processed from dsRNA intermediates, it could be expected that vsiRNAs are equally distributed toward sense and antisense strands. The observed strand bias of vsiRNAs in F. fusca may be due to vsiRNAs binding to highly abundant complimentary viral transcripts, as has been suggested in previous studies (Smith et al., 2010; Harris et al., 2015). Such binding could potentially sequester and remove vsiRNA from the pool for sRNAs submitted for sequencing. If such an explanation were correct, it is predicted that the difference between strand biases observed for A. hypogea compared to F. fusca may be due to host-specific differences in the level of viral RNA transcript expression and/or availability.
vsiRNA Preference and Hotspot Locations Evident in Alignment Profiles
The location and magnitude of vsiRNA hotspots spanning the TSWV genome in A. hypogaea and F. fusca samples are partially similar, though key differences exist (Figure 7A). In general, 21 and 22 nt vsiRNA hotspots within each species are co-located, with 21 nt peaks of greater magnitude in A. hypogaea, and 22 nt peaks larger in F. fusca. Visualization of vsiRNA coverage via circos plots demonstrates that relative to F. fusca, hotspots are more evenly dispersed across each genome segment in A. hypogaea. The F. fusca plots show fewer easily discernible peaks in the large RNA segment, and a greater abundance in the sRNA segment. Hotspots often occur in coding regions of the genome for both species, with a distinct gap apparent between adjacent genes located on the sRNA segment. Whilst hotspots do sometimes occur in the same locus in both the plant host and insect vector, many are unique to each organism; reinforcing that RNAi targeting mechanisms do have species-specific preferences. Importantly, these differences also indicate that TSWV-specific sRNAs in the insect are not simply derived from the plant itself via feeding, and are unique to the insect vector.
FIGURE 7
A comparison of A. hypogaea vsiRNA hotspots with those previously reported for S. lycopersicum and N. benthamiana infected with the same TSWV isolate (Mitter et al., 2013) shows distinct differences between plant families (Figure 7B). Similar to F. fusca, the vsiRNA hotspots in the L RNA segments corresponding to the RdRP gene are lacking in S. lycopersicum and N. benthamiana, indicating a very specific focus for A. hypogaea. Accordingly, though the external symptoms of TSWV infection are comparable in severity between A. hypogaea and S. lycopersicum, this A. hypogaea-specific direction of vsiRNAs toward the TSWV’s RdRP sequence may have significant downstream impacts upon pathogenicity relative to the Solanaceae.
The biogenesis and specific effects of vsiRNAs that form hotspots in alignment profiles are an area of ongoing study. Highly expressed vsiRNAs are thought to be a result of advantageous structural conditions in the confined regions of dsRNAs that vsiRNAs are generated from (Zhang et al., 2015). Differences in vsiRNA hotspot prevalence among plant families infected with the same viral isolate indicate that the localized environment the virus inhabits varies enough to impact RNA secondary structure, and thus the relative makeup of vsiRNAs generated.
The presence and absence of discrete vsiRNAs that are unique or shared between host and vector, or between different plant species infected with TSWV, reinforce the interpretations of the hotspot analyses. Unique vsiRNAs are readily apparent in F. fusca samples, further supporting the indication that vsiRNA generation takes place within the thrips, and not solely through uptake via feeding on vsiRNAs produced in infected A. hypogaea leaf tissue (Figure 8A). Additionally, the differences between A. hypogaea and the Solanaceae remain pronounced (Figure 8B). Due to the low TSWV vsiRNA load in N. benthamiana, almost all vsiRNAs present are common to S. lycopersicum. However, a similar number of discrete vsiRNAs only show partial overlap between A. hypogaea and N. benthamiana/S. lycopersicum, reinforcing the idea of a distinctive RNAi-mediated response to TSWV infection for different host plant families.
FIGURE 8
Conclusion
The unique sRNA and TSWV-induced vsiRNA profiles of the host plant A. hypogaea and thrips vector F. fusca are of specific interest in understanding the mechanisms of viral infection and repression. The relative abundance of sRNAs in the non-infected plant and insect samples, along with the 21 and 22 nt vsiRNAs in TSWV-infected samples, reflect the different means of sRNA biosynthesis in these kingdoms. Phenotypic changes resulting from TSWV infection in the thrips vector are not readily apparent, and indeed the abundance of generated vsiRNAs is relatively low. In contrast, infection with the same viral isolate results in a large-scale changes in the peanut host, from overall reductions of 24 nt sRNAs and increases in 22 nt sRNAs, to up- and down-regulation of sRNAs directed against repetitive elements and endogenous genes. The mechanism behind these perturbations, along with their full impact on gene expression and genome stability, are yet to be elucidated, but indicate that A. hypogaea is a species of particular significance in virus/host interaction studies.
Funding
This work was supported by the National Peanut Board (grant: NPB-2521RF330498) at the University of Georgia, the College of Agricultural, Human and Natural Resource Sciences, Agricultural Research Center at Washington State University (project WNPO 0545), and by Horticulture Innovation Australia Limited at the University of Queensland.
Statements
Author contributions
SF, AS, BC, RS, HP, and NM conceived and designed experiments. AS and SF performed experiments. SF, AS, and JP analyzed the data. All authors contributed to the writing, discussion, and approval of the final manuscript. SF and AS are to be considered as equal joint first authors.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fpls.2016.01349
References
1
BaulcombeD. (2004). RNA silencing in plants.Nature431356–363. 10.1038/nature02874
2
BolducF.HoareauC.St-PierreP.PerreaultJ. P. (2010). In-depth sequencing of the siRNAs associated with peach latent mosaic viroid infection.BMC Mol. Biol.11:16. 10.1186/1471-2199-11-16
3
BolognaN. G.VoinnetO. (2014). The diversity, biogenesis, and activities of endogenous silencing small RNAs in Arabidopsis.Annu. Rev. Plant Biol.65473–503. 10.1146/annurev-arplant-050213-035728
4
CampbellC. L.KeeneK. M.BrackneyD. E.OlsonK. E.BlairC. D.WiluszJ.et al (2008). Aedes aegypti uses RNA interference in defense against Sindbis virus infection.BMC Microbiol.8:47. 10.1186/1471-2180-8-47
5
CulbreathA. K.SrinivasanR. (2011). Epidemiology of spotted wilt disease of peanut caused by Tomato spotted wilt virus in the southeastern U.S.Virus Res.159101–109. 10.1016/j.virusres.2011.04.014
6
DashS.CannonE. K. S.KalbererS. R.FarmerA. D.CannonS. B. (2016). “Chapter 8 – peanutbase and other bioinformatic resources for peanut A2” in Peanuts, edsWilsonR. F.StalkerH. T. (Urbana, IL: AOCS Press), 241–252.
7
DehaanP.KormelinkR.ResendeR. D.VanpoelwijkF.PetersD.GoldbachR. (1991). Tomato spotted wilt virus-L RNA encodes a putative RNA-polymerase.J. Gen. Virol.722207–2216. 10.1099/0022-1317-72-9-2207
8
DehaanP.WagemakersL.PetersD.GoldbachR. (1990). The S-RNA segment of Tomato spotted wilt virus has an ambisense character.J. Gen. Virol.711001–1007. 10.1099/0022-1317-71-5-1001
9
DingS.-W.VoinnetO. (2007). Antiviral immunity directed by small RNAs.Cell130413–426. 10.1016/j.cell.2007.07.039
10
DingS. W. (2010). RNA-based antiviral immunity.Nat. Rev. Immunol.10632–644. 10.1038/nri2824
11
DonaireL.BarajasD.Martinez-GarciaB.Martinez-PriegoL.PaganI.LlaveC. (2008). Structural and genetic requirements for the biogenesis of tobacco rattle virus-derived small interfering RNAs.J. Virol.825167–5177. 10.1128/JVI.00272-08
12
FletcherS. J. (2016). Small Complementary RnA Mapper : A Quick and Simple Small RNA Read Analysis Package [Online]. Available: http://carroll-lab.github.io/scram/
13
GasciolliV.MalloryA. C.BartelD. P.VaucheretH. (2005). Partially redundant functions of Arabidopsis DICER-like enzymes and a role for DCL4 in producing trans-acting siRNAs.Curr. Biol.151494–1500. 10.1016/j.cub.2005.07.024
14
GroszmannM.GreavesI. K.AlbertynZ. I.ScofieldG. N.PeacockW. J.DennisE. S. (2011). Changes in 24-nt siRNA levels in Arabidopsis hybrids suggest an epigenetic contribution to hybrid vigor.Proc. Natl. Acad. Sci. U.S.A.1082617–2622. 10.1073/pnas.1019217108
15
HarrisC. J.MolnarA.MullerS. Y.BaulcombeD. C. (2015). FDF-PAGE: a powerful technique revealing previously undetected small RNAs sequestered by complementary transcripts.Nucleic Acids Res.437590–7599. 10.1093/nar/gkv604
16
HendersonI. R.ZhangX.LuC.JohnsonL.MeyersB. C.GreenP. J.et al (2006). Dissecting Arabidopsis thaliana DICER function in small RNA processing, gene silencing and DNA methylation patterning.Nat. Genet.38721–725. 10.1038/ng1804
17
HessA. M.PrasadA. N.PtitsynA.EbelG. D.OlsonK. E.BarbacioruC.et al (2011). Small RNA profiling of Dengue virus-mosquito interactions implicates the PIWI RNA pathway in anti-viral defense.BMC Microbiol.11:45. 10.1186/1471-2180-11-45
18
HongW.QianD.SunR. H.JiangL.WangY.WeiC. H.et al (2015). OsRDR6 plays role in host defense against double-stranded RNA virus, Rice Dwarf Phytoreovirus.Sci. Rep.511324. 10.1038/srep11324
19
JagadeeswaranG.ZhengY.SumathipalaN.JiangH.ArreseE. L.SoulagesJ. L.et al (2010). Deep sequencing of small RNA libraries reveals dynamic regulation of conserved and novel microRNAs and microRNA-stars during silkworm development.BMC Genomics11:52. 10.1186/1471-2164-11-52
20
JiaY.LischD. R.OhtsuK.ScanlonM. J.NettletonD.SchnableP. S. (2009). Loss of RNA-dependent RNA polymerase 2 (RDR2) function causes widespread and unexpected changes in the expression of transposons, genes, and 24-nt small RNAs.PLoS Genet.5:e1000737. 10.1371/journal.pgen.1000737
21
JohnsonL. M.DuJ.HaleC. J.BischofS.FengS.ChodavarapuR. K.et al (2014). SRA- and SET-domain-containing proteins link RNA polymerase V occupancy to DNA methylation.Nature507124–128. 10.1038/nature12931
22
JonesD. R. (2005). Plant viruses transmitted by thrips.Eur. J. Plant Pathol.113119–157. 10.1007/s10658-005-2334-1
23
KeeneK. M.FoyB. D.Sanchez-VargasI.BeatyB. J.BlairC. D.OlsonK. E. (2004). RNA interference acts as a natural antiviral response to O’nyong-nyong virus (Alphavirus; Togaviridae) infection of Anopheles gambiae. Proc. Natl. Acad. Sci. U.S.A.10117240–17245. 10.1073/pnas.0406983101
24
KempC.MuellerS.GotoA.BarbierV.ParoS.BonnayF.et al (2013). Broad RNA interference-mediated antiviral immunity and virus-specific inducible responses in Drosophila.J. Immunol.190650–658. 10.4049/jimmunol.1102486
25
KrzywinskiM.ScheinJ.BirolI.ConnorsJ.GascoyneR.HorsmanD.et al (2009). Circos: an information aesthetic for comparative genomics.Genome Res.191639–1645. 10.1101/gr.092759.109
26
LebasB. S. M.Ochoa-CoronaF. M. (2007). “Impatiens necrotic spot virus,” in Characterization, Diagnosis and Management of Plant Viruses, Grain Crops and OrnamentalsVol. 4edsRaoG. P. B.BragardC.LebasB. S. M. (Houston, TX: Stadium Press LLC), 221–243.
27
LinK. Y.ChengC. P.ChangB. C.WangW. C.HuangY. W.LeeY. S.et al (2010). Global analyses of small interfering RNAs derived from Bamboo mosaic virus and its associated satellite RNAs in different plants.PLoS ONE5:e11928. 10.1371/journal.pone.0011928
28
MalloryA. C.VaucheretH. (2004). MicroRNAs: something important between the genes.Curr. Opin. Plant Biol.7120–125. 10.1016/j.pbi.2004.01.006
29
MandalS. K. (2012). Field evaluation of alternate use of insecticides against chilli thrips, Scirtothrips dorsalis (Hood).Ann. Plant Protect. Sci.2059–62.
30
MargariaP.BoscoL.VallinoM.CiuffoM.MautinoG. C.TavellaL.et al (2014). The NSs protein of tomato spotted wilt virus is required for persistent infection and transmission by Frankliniella occidentalis.J. Virol.885788–5802. 10.1128/JVI.00079-14
31
MargariaP.MiozziL.RosaC.AxtellM. J.PappuH. R.TurinaM. (2015). Small RNA profiles of wild-type and silencing suppressor-deficient tomato spotted wilt virus infected Nicotiana benthamiana.Virus Res.20830–38. 10.1016/j.virusres.2015.05.021
32
McHaleM.EamensA. L.FinneganE. J.WaterhouseP. M. (2013). A 22-nt artificial microRNA mediates widespread RNA silencing in Arabidopsis.Plant J.76519–529. 10.1111/tpj.12306
33
MiozziL.GambinoG.BurgyanJ.PantaleoV. (2013). Genome-wide identification of viral and host transcripts targeted by viral siRNAs in Vitis vinifera.Mol. Plant Pathol.1430–43. 10.1111/j.1364-3703.2012.00828.x
34
MitterN.KoundalV.WilliamsS.PappuH. (2013). Differential expression of Tomato spotted wilt virus-derived viral small RNAs in infected commercial and experimental host plants.PLoS ONE8:e76276. 10.1371/journal.pone.0076276
35
MuellerS.GaussonV.VodovarN.DeddoucheS.TroxlerL.PerotJ.et al (2010). RNAi-mediated immunity provides strong protection against the negative-strand RNA vesicular stomatitis virus in Drosophila.Proc. Nat. Acad. Sci. U.S.A.10719390–19395. 10.1073/pnas.1014378107
36
MungerF. (1942). A method for rearing citrus thrips in the laboratory.J. Econ. Entomol.35373–375. 10.1093/jee/35.3.373
37
MylesK. M.WileyM. R.MorazzaniE. M.AdelmanZ. N. (2008). Alphavirus-derived small RNAs modulate pathogenesis in disease vector mosquitoes.Proc. Nat. Acad. Sci. U.S.A.10519938–19943. 10.1073/pnas.0803408105
38
PappuH. R.JonesR. A.JainR. K. (2009). Global status of tospovirus epidemics in diverse cropping systems: successes achieved and challenges ahead.Virus Res.141219–236. 10.1016/j.virusres.2009.01.009
39
ParentJ. S.BouteillerN.ElmayanT.VaucheretH. (2015). Respective contributions of Arabidopsis DCL2 and DCL4 to RNA silencing.Plant J.81223–232. 10.1111/tpj.12720
40
PelaezP.SanchezF. (2013). Small RNAs in plant defense responses during viral and bacterial interactions: similarities and differences.Front. Plant Sci.4:343. 10.3389/fpls.2013.00343
41
PersleyD. M.ThomasJ. E.SharmanM. (2006). Tospoviruses: an Australian perspective.Austral. Plant Pathol.35161–180. 10.1071/AP06015
42
PontesO.LiC. F.Costa NunesP.HaagJ.ReamT.VitinsA.et al (2006). The Arabidopsis chromatin-modifying nuclear siRNA pathway involves a nucleolar RNA processing center.Cell12679–92. 10.1016/j.cell.2006.05.031
43
RileyD. G.JosephS. V.SrinivasanR.DiffieS. (2011). Thrips vectors of tospoviruses.J. Int. Pest Manag.11–10.
44
RotenbergD.Krishna KumarN. K.UllmanD. E.Montero-AstuaM.WillisD. K.GermanT. L.et al (2009). Variation in tomato spotted wilt virus titer in Frankliniella occidentalis and its association with frequency of transmission.Phytopathology99404–410. 10.1094/PHYTO-99-4-0404
45
SchnettlerE.HemmesH.HuismannR.GoldbachR.PrinsM.KormelinkR. (2010). Diverging affinity of tospovirus RNA silencing suppressor proteins, NSs, for various RNA duplex molecules.J. Virol.8411542–11554. 10.1128/JVI.00595-10
46
ScholthofK. B.AdkinsS.CzosnekH.PalukaitisP.JacquotE.HohnT.et al (2011). Top 10 plant viruses in molecular plant pathology.Mol. Plant Pathol.12938–954. 10.1111/j.1364-3703.2011.00752.x
47
SilvaT. F.RomanelE. A. C.AndradeR. R. S.FarinelliL.OsterasM.DeluenC.et al (2011). Profile of small interfering RNAs from cotton plants infected with the polerovirus Cotton leafroll dwarf virus.BMC Mol. Biol.12:40. 10.1186/1471-2199-12-40
48
SmithN. A.EamensA. L.WangM. B. (2010). The presence of high-molecular-weight viral RNAs interferes with the detection of viral small RNAs.RNA161062–1067. 10.1261/rna.2049510
49
SnippeM.SmeenkL.GoldbachR.KormelinkR. (2007). The cytoplasmic domain of tomato spotted wilt virus Gn glycoprotein is required for Golgi localisation and interaction with Gc.Virology363272–279. 10.1016/j.virol.2006.12.038
50
SoellickT. R.UhrigJ. F.BucherG. L.KellmannJ. W.SchreierP. H. (2000). The movement protein NSm of tomato spotted wilt tospovirus (TSWV): RNA binding, interaction with the TSWV N protein, and identification of interacting plant proteins.Proc. Nat. Acad. Sci. U.S.A.972373–2378. 10.1073/pnas.030548397
51
SzittyaG.MoxonS.PantaleoV.TothG.PilcherR. L. R.MoultonV.et al (2010). Structural and functional analysis of viral siRNAs.PLoS Pathog.6:e1000838. 10.1371/journal.ppat.1000838
52
TakedaA.SugiyamaK.NaganoH.MoriM.KaidoM.MiseK.et al (2002). Identification of a novel RNA silencing suppressor, NSs protein of Tomato spotted wilt virus.FEBS Lett.53275–79.
53
TsompanaM.MoyerJ. W. (2008). “Tospoviruses,” in Encyclopedia of Virology, 3rd Edn, edsMahyB. W. J.Van RegenmortelM. H. V. (Oxford: Elsevier Ltd), 157–162.
54
UllmanD. E.GermanT. L.SherwoodJ. L.WestcotD. M.CantoneF. A. (1993). Tospovirus replication in insect vector cells – immunocytochemical evidence that the nonstructural protein encoded by the S-RNA of Tomato spotted wilt tospovirus is present in thrips vector cells.Phytopathology83456–463. 10.1094/Phyto-83-456
55
WangX. B.WuQ. F.ItoT.CilloF.LiW. X.ChenX. M.et al (2010). RNAi-mediated viral immunity requires amplification of virus-derived siRNAs in Arabidopsis thaliana.Proc. Nat. Acad. Sci. U.S.A.107484–489. 10.1073/pnas.0904086107
56
WeiY.ChenS.YangP.MaZ.KangL. (2009). Characterization and comparative profiling of the small RNA transcriptomes in two phases of locust.Genome Biol.10 R6. 10.1186/gb-2009-10-1-r6
57
WhitfieldA. E.UllmanD. E.GermanT. L. (2005). Tospovirus-thrips interactions.Annu. Rev. Phytopathol.43459–489. 10.1146/annurev.phyto.43.040204.140017
58
WierzbickiA. T.ReamT. S.HaagJ. R.PikaardC. S. (2009). RNA polymerase V transcription guides ARGONAUTE4 to chromatin.Nat. Genet.41630–634. 10.1038/ng.365
59
WuQ. F.LuoY. J.LuR.LauN.LaiE. C.LiW. X.et al (2010). Virus discovery by deep sequencing and assembly of virus-derived small silencing RNAs.Proc. Nat. Acad. Sci. U.S.A.1071606–1611. 10.1073/pnas.0911353107
60
XuY.HuangL.FuS.WuJ.ZhouX. (2012). Population diversity of rice stripe virus-derived siRNAs in three different hosts and RNAi-based antiviral immunity in Laodelphax striatellus.PLoS ONE7:e46238. 10.1371/journal.pone.0046238
61
ZhangC.WuZ.LiY.WuJ. (2015). Biogenesis, function, and applications of virus-derived Small RNAs in plants.Front. Microbiol.6:1237. 10.3389/fmicb.2015.01237
Summary
Keywords
tomato spotted wilt virus, Arachis hypogea, Frankliniella fusca, plant defense, plant–insect–virus interactions, vsiRNA, small RNA
Citation
Fletcher SJ, Shrestha A, Peters JR, Carroll BJ, Srinivasan R, Pappu HR and Mitter N (2016) The Tomato Spotted Wilt Virus Genome Is Processed Differentially in its Plant Host Arachis hypogaea and its Thrips Vector Frankliniella fusca. Front. Plant Sci. 7:1349. doi: 10.3389/fpls.2016.01349
Received
03 June 2016
Accepted
22 August 2016
Published
07 September 2016
Volume
7 - 2016
Edited by
Agnieszka Ludwików, Adam Mickiewicz University in Poznań, Poland
Reviewed by
Agnieszka Kiełbowicz-Matuk, Institute of Plant Genetics of the Polish Academy of Science, Poland; Neeti Sanan-Mishra, International Center for Genetic Engineering and Biotechnology, India
Updates
Copyright
© 2016 Fletcher, Shrestha, Peters, Carroll, Srinivasan, Pappu and Mitter.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Neena Mitter, n.mitter@uq.edu.au
†These authors have contributed equally to this work.
This article was submitted to Plant Biotechnology, a section of the journal Frontiers in Plant Science
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.