Tomato Spotted Wilt Virus Benefits Its Thrips Vector by Modulating Metabolic and Plant Defense Pathways in Tomato

Abstract

Several plant viruses modulate vector fitness and behavior in ways that may enhance virus transmission. Previous studies have documented indirect, plant-mediated effects of tomato spotted wilt virus (TSWV) infection on the fecundity, growth and survival of its principal thrips vector, Frankliniella occidentalis , the western flower thrips. We conducted thrips performance and preference experiments combined with plant gene expression, phytohormone and total free amino acid analyses to determine if systemically-infected tomato plants modulate primary metabolic and defense-related pathways to culminate into a more favorable environment for the vector. In a greenhouse setting, we documented a significant increase in the number of offspring produced by F. occidentalis on TSWV-infected tomato plants compared to mock-inoculated plants, and in choice test assays, females exhibited enhanced settling on TSWV-infected leaves. Microarray analysis combined with phytohormone signaling pathway analysis revealed reciprocal modulation of key phytohormone pathways under dual attack, possibly indicating a coordinated and dampening defense against the vector on infected plants. TSWV infection, alone or in combination with thrips, suppressed genes associated with photosynthesis and chloroplast function thereby significantly impacting primary metabolism of the host plant, and hierarchical cluster and network analyses revealed that many of these genes were co-regulated with phytohormone defense signaling genes. TSWV infection increased expression of genes related to protein synthesis and degradation which was reflected in the increased total free amino acid content in virus-infected plants that harbored higher thrips populations. These results suggest coordinated gene networks that regulate plant primary metabolism and defense responses rendering virus-infected plants more conducive for vector colonization, an outcome that is potentially beneficial to the vector and the virus when considered within the context of the complex transmission biology of TSWV. To our knowledge this is the first study to identify global transcriptional networks that underlie the TSWV-thrips interaction as compared to a single mechanistic approach. Findings of this study increase our fundamental knowledge of host plant-virus-vector interactions and identifies underlying mechanisms of induced host susceptibility to the insect vector.

Introduction

Most plant-pathogenic viruses depend exclusively on insects for transmission to plants (Hogenhout et al., 2008). It would therefore seem to be evolutionarily advantageous for a plant virus to modify vector behavior and performance in ways that enhance its likelihood of acquisition and dissemination in the landscape [reviewed inEigenbrode et al. (2018)]. Indeed, there are cases that demonstrate the direct effect of some plant viruses to alter vector performance. For example, acquisition of barley yellow dwarf virus (BYDV) altered aphid vector preference for non-infected compared to BYDV-infected plants (Ingwell et al., 2012) and acquisition of tomato mottle virus increased oviposition of its whitefly vector (McKenzie, 2002). Indirect or plant-mediated effects due to virus infection also play a role in vector-plant interactions. For example, plants infected with potato leafroll virus and BYDV were not only more attractive to aphid vectors, but also increased nutrient quality that enhanced vector survival and fecundity (Eigenbrode et al., 2002;Jimenez-Martinez et al., 2004;Ngumbi et al., 2007). These studies highlight the widespread occurrence of vector manipulation by plant viruses across several pathosystems for aphids (Fereres et al., 1989;Blua et al., 1994;Eigenbrode et al., 2002;Mauck et al., 2010), whiteflies (Guo et al., 2010), and thrips (Carter, 1939;Maris et al., 2004;Belliure et al., 2005;Abe et al., 2012).

Among thrips-transmitted tospoviruses, Tomato spotted wilt orthotospovirus (order Bunyavirales, family Tospoviridae, genus Orthotospovirus) is ranked among the top 10 most economically important plant viruses worldwide (Rybicki, 2015). The virus is transmitted in a circulative-persistent manner exclusively by thrips, the most efficient of which is the western flower thrips, Frankliniella occidentalis (Pergande) (Whitfield et al., 2005,2015;Rotenberg and Whitfield, 2018). The thrips-tospovirus relationship is unique in that adult thrips are only able to transmit TSWV if acquisition occurs in the first instar and early second larval thrips stages (Van De Wetering et al., 1996). The virus then replicates in the midgut and adjacent tissues, and eventually reaches the primary salivary glands (Ullman et al., 1992;Montero-Astúa et al., 2016;Rotenberg and Whitfield, 2018). After infection of the salivary glands, adults during feeding release the virus into viable plant cells along with the saliva. Hence, TSWV replicates in both the vector and the host plant, offering the opportunity for both direct and indirect (i.e., plant-mediated) effects on the vector. Studies of direct effects are rare (Belliure et al., 2005;Shrestha et al., 2012) but there are numerous reports of indirect effects of orthotospoviruses on thrips vectors. For example, TSWV infection has been shown to increase vector fecundity, development, population growth, and survival (Carter, 1939;Maris et al., 2004;Belliure et al., 2005;Shrestha et al., 2012;Ogada et al., 2013) and longevity (Ogada et al., 2013) on infected plants relative to uninfected plants. Although most documented virus effects are positive, there are reports of negative (DeAngelis et al., 1993;Stumpf and Kennedy, 2005) and neutral (Wijkamp et al., 1996) effects of TSWV on thrips vectors as well. These reports suggest that TSWV infection can alter plant physiology to benefit vector fitness.

The literature on the molecular mechanisms underlying the tospovirus-thrips vector-plant interaction is emerging. Most studies have focused on defense-related signaling pathways as a potential mechanism in the model plants, Arabidopsis thaliana and Nicotiana benthamiana.Abe et al. (2012)demonstrated the role of antagonistic crosstalk between salicylic acid (SA) and jasmonic acid (JA)-responses in the TSWV-thrips interaction using Arabidopsis mutants. The authors suggested that increased performance of thrips on TSWV-infected plants was caused by a reduction of the JA-regulated plant defense, which was suppressed by an increase in SA-regulated plant defense responses in virus-infected plants. More recently,Wu et al. (2019)showed that the TSWV non-structural protein, NSs, suppressed biosynthesis of monoterpenes, which are known to repel F. occidentalis by directly interacting with MYCs, key regulators of the JA signaling pathway. To date, there have been four studies on transcriptional responses (microarray and RNA-Seq) to TSWV in host plants which report changes in plant immune defenses and metabolism in virus-infected plants (Catoni et al., 2009;Nachappa et al., 2013;Padmanabhan et al., 2019;Xu et al., 2020). However, these studies did not analyze global transcriptional changes in response to the thrips vector and the combination treatment of virus and vector.

The roles of the three major phytohormones, SA, JA, and ethylene (ET) in defense responses to pathogens and insects is well-established (Walling, 2000;Stout et al., 2006;Howe and Jander, 2008;Bari and Jones, 2009;Pieterse et al., 2009;Verhage et al., 2010). However, research in the past decade demonstrates that plant defense responses is more than just SA and JA/ET pathways, with more coordination and integration of a range of hormones including, abscisic acid (ABA), auxins, brassinosteroids, cytokinins and gibberellins (De Bruyne et al., 2014;Shigenaga and Argueso, 2016;Yue et al., 2016;Berens et al., 2017;Yang et al., 2019). For instance, ABA can act as both a positive and negative regulator of disease resistance. ABA can suppress SA-mediated defenses, and plant susceptibility to pathogens can increase following exogenous applications of ABA (Asselbergh et al., 2008). In certain instances, exogenous application of ABA can have the opposite effect on SA-mediated defenses resulting in increased resistance to pathogens (Ton and Mauch-Mani, 2004;Wiese et al., 2004;Melotto et al., 2006). The antagonism of SA-mediated defenses by ABA may be explained in part by the positive effect of ABA on JA biosynthesis (Adie et al., 2007). These phytohormones not only mediate immunity, but also growth and development (Pieterse et al., 2009;De Bruyne et al., 2014;Berens et al., 2017;Züst and Agrawal, 2017). For example, the classic growth hormone-auxin acts in an antagonistic manner with SA during plant defense, whereas auxin and JA signaling act synergistically. Moreover, some pathogens either produce auxin themselves or increase plant auxin biosynthesis to manipulate plant defense responses and development [reviewed inKazan and Manners (2009)]. The evolutionary conservation of the intricate network of phytohormone signaling pathways likely enables plants to effectively balance trade-offs between defense and growth (Züst and Agrawal, 2017).

In the current study, we sought to elucidate the plant-mediated mechanisms underlying the interaction between TSWV and its insect vector, F. occidentalis, on the plant host, tomato, Solanum lycopersicum L. First, we performed replicated greenhouse and laboratory experiments to confirm that TSWV altered vector performance and behavior on TSWV-infected plants. To characterize plant molecular mechanisms, microarray analyses using the Affymetrix Tomato GeneChip^® were performed to analyze gene expression in tomato plants receiving one of four treatments: (i) mock-inoculated, or (ii) systemically-infected with TSWV, and subsequently infested (iii) with or (iv) without non-viruliferous thrips. The microarray represents myriad genes associated with primary metabolism, including chloroplast function, cell wall modification, protein synthesis, as well as defense- and stress-related genes associated with phytohormone signaling pathways (JA, SA, ET, ABA, and auxin). In addition, we analyzed phytohormone levels and total free amino acid content in plants exposed to single and dual challengers. This is the first comprehensive analysis of plant-mediated mechanisms (genes and metabolites) that potentially improve the quality of TSWV-infected plants for its thrips vectors.

Materials and Methods

Greenhouse Experimental Design and Treatment Structure

The greenhouse experiment consisted of four treatments: (1) TSWV infection alone; (2) thrips infestation alone; (3) TSWV infection and thrips infestation, and (4) mock-inoculated, healthy controls. Each treatment was replicated three or four times (i.e., three or four plants) in a randomized complete block design with one plant per treatment per block. The experiment was conducted three times (i.e., three biological replications).

TSWV and Thrips Sources

Tomato spotted wilt virus (isolate TSWV-MT2) was maintained by mechanical inoculations on caged tomato plants (S. lycopersicum cv Moneymaker) under greenhouse conditions as per (Rotenberg et al., 2009). A colony of thrips was maintained on green bean pods (Phaseolus vulgaris) as described byBautista et al. (1995)at ambient room temperatures of 24 ± 1°C and a photoperiod of 16:8 h (L:D) light: dark cycle.

Plant Growth Conditions, Virus Inoculations, and Thrips Release

Tomato plants (cv Moneymaker) were grown in 6-inch pots filled with Metro mix^® potting soil and housed in one thrips-proof-screened (No Thrips Insect^® Screen, BioQuip Products, Rancho Dominguez, CA) cage in a greenhouse room. Plants were fertilized once a week with Miracle Gro^®-Water Soluble All Purpose Plant Food (24-8-16) NPK. The temperatures in the greenhouse ranged from 23 to 25°C and the photoperiod was 16:8 h (L:D).

To generate TSWV-infected plants, 3-week-old plants were mechanically-inoculated with TSWV from infected plant tissue or mock-inoculated with buffer and healthy plant tissue. Virus inoculum was prepared by grinding two to three young symptomatic tomato leaves in ice-cold 5–10 ml of inoculation buffer (10 mM sodium sulfite and 5% wt/vol celite) using a pre-chilled mortar and pestle. Inoculum was applied and dragged lightly with a cotton swab over the surface of all fully-expanded leaves on the plant. Approximately 2 weeks after inoculation, plants were visually inspected for TSWV symptoms, i.e., stunting and deformation, chlorotic ring spots, mosaic patterns, and leaf-bronzing, and the most uniform group of symptomatic plants were chosen for each experiment.

Individual 5-week-old symptomatic and mock-inoculated control plants were moved to single-plant cages constructed from 19-liter cardboard ice cream buckets (38 cm tall × 26 cm diameter) with four 14 cm × 27 cm apertures cut into the side walls. The four apertures were covered with No-Thrips Insect Screen^TM (Green-Tek, Inc., Edgerton, WI, United States) and sealed with silicone. The top of the container was covered by thrips-proof-screen secured by rubber bands. These single-plant cages prevented cross-contamination between treatments. Cohorts of 40 adult females from the laboratory colony, 7-day post eclosion were transferred into 1.5 ml microcentrifuge tubes and then allowed to escape in a synchronous manner from one tube placed at the base of each plant.

Thrips Feeding Damage, Performance and Settling Assays

Thrips feeding leaf damage index (LDI) was determined using a modified visual rating scale according toMcKinney (1923). All leaves from each plant were evaluated visually and grouped into six classes with values of 0, 1, 3, 5, 7, and 9, which corresponded to percentages of total leaf area damaged by thrips feeding of 0, <10, 10–25, 25–50, 50–75, and >75%, respectively. In addition, the number of thrips feeding lesions per plant were counted.

Thrips performance on virus-infected plants or mock-inoculated plants was defined as the number of offspring that emerged from leaf tissue at 7-day post release of adult thrips. For each biological replicate, 2-sample rank tests (i.e., Mann–Whitney) using Minitab v.14 (Minitab, Inc., State College, PA, United States) was performed to determine if the median count of offspring obtained from TSWV-infected plants was similar to that obtained from healthy, mock-inoculated plants. Data from all three experiments were pooled together and a two-way analysis of variance in Minitab using a GLM with treatment and biological replicate as fixed effects.

To determine adult thrips settling preference, assays using detached-leaflets were conducted as per (Nachappa et al., 2013). Same-age leaflets were taken from mock-inoculated plants and plants with TSWV infection alone at the termination of each greenhouse experiment (i.e., 6-week-old plants). Leaflet pairings consisted of choice and non-choice tests using pairs of leaflets in 15-cm diameter Petri-dishes with the lids of each dish fitted with thrips-proof screen to allow ventilation. Each Petri dish contained 15 ml of 1.5% water agar onto which were placed the adaxial surfaces of each paired leaflet. The agar prevented desiccation of leaflets while allowing thrips to move across the surface and choose leaflets. Ten adult female thrips were placed with a small paint brush in the center of the Petri dish, equidistant from each leaflet, and lids were sealed with Parafilm. The assay was conducted under laboratory conditions with a 16:8 h photoperiod and ambient temperatures of 24–25°C. Adult thrips preference (choice of virus infected vs. mock controls) was determined by counting the number of thrips present on each of the paired leaflet every hour for the first 6 h, and then at 12, 24, 36, 48, and 72 h after the release. For each biological replicate, we performed 1-Sample Wilcoxon sign rank tests using Minitab to determine if the median paired difference in accumulation on virus-infected vs. non-infected healthy tissue was significantly greater than zero at each time-point.

Leaf Tissue Sampling for Gene Expression, Phytohormone and Amino Acid Analyses

Preliminary observations revealed that thrips preferred to feed on older (basal) leaves, so leaflets were harvested from 10th–11th youngest leaf, counting down from the top to the base of the plant. From each of the three experiments (biological replications), two same-age, paired leaflets (on each side of leaf rachis) were harvested. We harvested leaflets exhibiting similar leaf damage ratings (thrips alone average LDI = 2; TSWV + thrips average LDI = 2.66) between the two thrips treatments (thrips alone and TSWV + thrips) to reduce the confounding effect of variation due to amount of feeding (Supplementary Table 1). In the current study, thrips were not localized or caged to a specific leaf rather allowed to colonize the whole plant. Given their thigmotactic behavior, i.e., preferring to hide in small crevices on plant surfaces, and high mobility of adult thrips, it is difficult to correlate feeding damage with number of thrips. Hence, we used LDI and number of lesions as a more reliable estimate of the effect of thrips feeding on plants rather than sighting the insects. Moreover, there was a significant positive correlation between LDI and lesions per plant (r = 0.739, P < 0.0001). One leaflet was processed for gene expression analysis (i.e., microarray hybridizations and RT-qPCR analysis) and the other was processed to determine phytohormone contents. Leaflet samples were flash-frozen in liquid nitrogen and stored at −80°C. To obtain enough leaf tissue (approximately 5 g) required for determination of total free amino acid content, a third leaf immediately basal to the leaflets chosen for microarray and phytohormone analyses was simultaneously harvested and freeze-dried.

RNA Isolation, Amplification, Labeling and Generation of GeneChip Data

Total RNA isolation and cDNA synthesis was performed as perNachappa et al. (2013). Briefly, RNA was isolated from frozen leaflet samples (100 mg of tissue) using the Qiagen RNeasy Plant Mini Kit (Qiagen, Valencia, CA, United States) following manufacturer’s protocol. RNA quantity was determined by NanoDrop spectrophotometer (NanoDrop Technologies, Wilmington, DE, United States) and quality was assessed with the 2100 Bioanalyzer using Nanochip technology (Agilent Technologies, Inc., Palo Alto, CA, United States). Total RNA was pooled for each treatment (3–4 experimental replicates, i.e., plants) within a biological replication. Pooled RNA samples were subjected to cRNA synthesis, labeling, and hybridization to Affymetrix Tomato Genome Arrays (GeneChip)^® at Kansas State University, Integrated Genomics Center. Each GeneChip^® contains more than 10,000 probe sets for over 9,200 genes with each gene being represented by at least one probe set containing 25-mer oligonucleotides. Hybridization intensities of scanned microarrays for each of the three biological replicates were generated with Gene Chip^® Operating System, GCOS (Affymetrix Inc.). Global scaling was applied for each GeneChip to adjust the Target Intensity (TGT) Value to an arbitrary target of 500 so that hybridization intensity of all chips was equivalent. In addition, expressed genes were identified by GCOS, using a detection algorithm and assigned a present, marginal, or absent call for genes represented by each probe set on the array (GeneChip Expression Analysis Technical Manual). Microarray data files (.CEL) were analyzed using GeneSpring 10.1 (Silicon Genetics, Redwood, CA, United States) and normalized using RMA (Robust Multichip Average) algorithm. Differentially-expressed genes were selected using 2 criteria: (i) an expression ratio of at least ±2-fold change and (ii) P ≤ 0.05 in ANOVA tests comparing log2 (normalized hybridization intensity) of treatment to the mock control. Differentially-expressed genes were assigned functional annotations with Blast2GO software (Conesa and Götz, 2008) and classified by GO-biological process, cellular component and molecular function using default parameters and an E-value cut-off of 10^–5. The microarray experiment design details and raw microarray data is available at ArrayExpress under the accession numberE-MTAB-9294.

Hierarchical Cluster Analysis

We used the heatmap () base function in R to create a heatmap of a sub-set of 369 differentially expressed genes selected for their membership in candidate pathways (defense, phytohormones, photosynthesis, cell wall organization, and protein metabolism). The colors range from red to green, showing the expression of genes in each treatment group ranging from highly similar (green) to less similar (red). The heatmap was plotted from the adjacency matrix of the network, using the R heatmap () function. The Pearson correlation coefficients depicts correlation between genes and treatments based on the average fold change values from 3 replicates for each treatment.

Phytohormone Signaling Pathway Analysis

We analyzed changes in tomato phytohormone pathways for SA, JA, ET, ABA, and AUX-responses associated with each treatment using the pathway analysis method developed byStudham and MacIntosh (2012). This method provides a broad indication (pathway score) of relative upregulation or downregulation of individual phytohormone pathways based on the magnitude (fold change), direction of change (up or down), and significance of differentially-expressed genes (compared to mock treatment) associated with each pathway for each treatment. Pathway scores for each phytohormone pathway can be compared to provide a comprehensive picture of treatment-associated, host plant responses and to identify dominant pathways affected by treatment. To perform the pathway analysis, annotations for probe sets in the Tomato GeneChip^® were obtained manually based on literature and from Blast2GO biological processes¹ (Conesa and Götz, 2008). Additional annotations for probe sets that were not assigned a function were obtained by identifying the EC numbers using the annot8r program². The EC numbers were used to identify pathway membership by querying UniProt³ and KEGG (Kanehisa et al., 2011). The genes associated with different functional categories have different roles such as biosynthesis, response, regulation, and others. Each role was assigned different weights indicating its relative importance in determining the induction or suppression of a hormone/protein (Studham and MacIntosh, 2012). A cumulative pathway score was then calculated for the subset of genes that exhibited microarray fold changes ≥ 1.0 and p-values ≤ 0.05 for each phytohormone pathway.

Co-expression Network Analysis

Gene co-expression network analysis is used to describe the correlation patterns among genes across microarray and other multidimensional expression data sets (e.g., RNAseq). The weighted gene correlation network analysis method (WGCNA,Langfelder and Horvath, 2008) can be used to identify clusters (modules) of highly correlated genes, and each module provides information about the pairwise relationships (correlations) between genes. In addition, community analysis of gene co-expression networks can identify communities (clusters and modules) of nodes. The nodes in a particular community have a higher likelihood of connecting to each other than to nodes from other communities (Barabási, 2016), and these relationships can provide additional information about groups of genes that are likely to be expressed together under particular experimental conditions. The WGCNA R package⁴ (Langfelder and Horvath, 2008) was used to analyze a sub-set of differentially expressed genes in candidate functional categories of interest [i.e., photosynthesis, plant hormone-related (ABA, AUX, ET, JA, and SA), protein metabolism and turnover, cell wall and defense functions] across four treatments (mock-infected, TSWV, thrips, and TSWV + thrips). The sample code in the WGCNA tutorial⁵ (Langfelder and Horvath, 2008) was adapted to perform WGCNA on our expression data set. The input data to WGCNA consisted of averaged three replicates for each of the four treatments. Code from the first tutorial on WGCNA network analysis⁶ was adapted and used to construct the correlation network. The nodes and edges of the WGCNA network were exported to files using the export network to Cytoscape (Shannon et al., 2003) function in the WGCNA package. The network figure was created using the R igraph package⁷. Different community detection methods were applied to the network data to determine the optimal method and number of community modules, and the leading eigen method was selected (Newman, 2006). As our goal was to visualize the network nodes clearly, we tried most of the igraph network layouts and compared the resulting network images. The ‘layout_with_fr’ igraph layout was chosen for the figure because it produced the clearest visualization of the network. This layout uses the force-directed layout algorithm of Fruchterman and Reingold to place the network vertices on a plane (Fruchterman and Reingold, 1991). The resulting network included five communities of nodes, each community representing a functional category (see list above).

Confirmation of Microarray Data Using Reverse Transcription Quantitative –PCR (RT-qPCR)

We targeted six genes associated with the SA, JA, and antiviral small-RNA-mediated gene silencing pathways. These genes were BGL2 and NPR1 (SA pathway) and OPR3, AOS and CI (JA pathway), and RNA-directed RNA polymerase 1 (RDR1). We chose elongation factor 1-alpha (leEF1) as the internal reference gene for normalization because expression of this gene was found to be invariant with TSWV infection or herbivore challenge in the microarray experiments and had been previously shown to be stably-expressed in Moneymaker tomato systemically-infected with tobacco rattle virus (Rotenberg et al., 2006). Target and leEF1 primer pair sequences, their corresponding melting temperatures, and real-time PCR efficiencies are indicated inTable 1. The normalized abundance of TSWV nucleocapsid (N) RNA compared to leEF1 was also determined to estimate virus titer in leaf tissue using primers tested previously (Rotenberg et al., 2009). We selected one plant per treatment (mock included) per biological replicate (i.e., 24 RNA samples in total) of the greenhouse experiment that represented the average fecundity and/or TSWV symptom severity for a given treatment. Subsamples of total RNA isolated from leaflet tissue used in the microarray hybridization experiment were treated with DNase using the rigorous DNA removal procedure of the Turbo DNA-free kit (Applied Biosystems Inc., Carlsbad, CA, United States) and cDNA was synthesized from 1 μg DNA-free RNA using the iSCRIPT cDNA synthesis kit (Bio-Rad, Hercules, CA, United States). Real-time PCR master mixes were prepared using iQ sybr Green Mix (Bio-Rad) according to manufacturer’s specifications and final reaction (20 μl) concentrations of 200 nM of each primer. Reactions were performed in duplicate using the iCycler iQ Thermal Cycler with a 96 ml × 0.2 ml reaction module and iCycler iQ software (Bio-Rad).

TABLE 1

Real-time quantitative reverse transcriptase – PCR primer pair sequences and corresponding PCR efficiencies.

Primer name	Gene name/accession number	Primer sequence (5′–3′, forward/reverse)	aPCR efficiency
AOS	Allene oxide synthase/AJ271093	ATCGTCTTATCGTGTTAGTATTC/	1.98
GATGATGATGGTGATTGTGAT
BGL2	Beta-1,3-glucanase/M80604	CTTGTTGGGCTTCTAATCC/	1.91
CTTGATCCGATGGTAAATTATTG
CI	Cathepsin D inhibitor protein/X73986	GCGTTAGGTGGTGATGTA/	1.97
GAATTGTAGGTCCATTAGTTGAT
leEF-1	Elongation factor -1 alpha/X14449	GATTGGTGGTATTGGAACTGTC/	1.97
AGCTTCGTGGTGCATCTC
NPR1	Non-pathogenesis related protein	GATAAGTCCTTGCCTCAT/	2.00
1/NM_001247629	AATGCTCTATGTATCCTCTT
OPR3	12-oxophytodienoate reductase/AJ24255	GGTGGTTACGATAGAGAAGA/	1.91
GGATAATCAGTATAGCCAACAAT
RDR1	RNA-directed RNA polymerase 1/Y10403	GCGACCTTCACAAGAGAT/	1.80
TCATAATGCCACCACTAAGT
TSWV-Nb	TSWV nucleocapsid gene/AF306490	GCTTCCCACCCTTTGATTC/	1.90
ATAGCCAAGACAACACTGATC
TSWV-NSs	TSWV non-structural protein (silencing suppressor)/NC_002051.1:89-1483	ACTCTGTTCTGGCACTATCTG/	2.03
GCTGGAATCGGTCTGTAATAT

^aPCR efficiencies were calculated as 10^–1/slope fromPfaffl (2001). ^bPrimer sequences obtained fromRotenberg et al. (2009).

The relative abundance of target RNA was determined for each virus and herbivore treatment compared to the mock control. The relative expression ratio (RER) equation (Pfaffl, 2001) was calculated as follows: RER = E _target ^{ΔCttarget[control–treatment]}/E _ref ^{ΔCtref[control–treatment]} where E refers to the PCR primer efficiency for target or internal reference (ref; leEF1) genes and ΔCt is the difference in Ct-values (i.e., threshold cycle values automatically calculated by the IQ software) between a treatment and mock control. The average Ct value (n = 3) obtained for the mock control for each target and reference gene was determined and used in the RER calculations. To estimate virus titer, the normalized abundance of TSWV N or NSs RNA (genomic and transcript RNA) was calculated using the Pfaffl inverse equation (Pfaffl, 2001): E _ref ^Ctref/E _N ^CtN as described previously (Rotenberg et al., 2009).

Phytohormone Analysis

Frozen tomato leaflet samples from each biological replicate were sent to The Donald Danforth Proteomics & Mass Spectrometry Facility, St. Louis, MO, United States for chemical extraction and quantification of SA, JA, ABA, jasmonyl-isoleucine (JA-Ile), and 12-oxo-phytodienoic acid (OPDA) using liquid chromatography–electrospray tandem mass spectrometry as per (Pan et al., 2008). Analysis of variance was performed on log₁₀-transformed phytohormone contents (ng analyte g fresh weight^–1) in Minitab using a GLM that included treatment and biological replicate as two fixed factors and their interaction term as the random factor. The analyses revealed no apparent main effect or treatment interaction for any of the phytohormones measured due to biological replicate; therefore, data was combined for the three biological replicates and one-way ANOVA was performed to determine the main effect of treatment on phytohormone concentration. Pairwise treatment comparisons were performed using Tukey’s family error rate (P ≤ 0.05).

Total Free Amino Acid Analyses

Frozen tomato leaf samples obtained from each experiment were sent to the Ruminant Nutrition Laboratory at Kansas State University for extraction and quantification of total free amino acids. Total free amino acid content was analyzed using a modified protocol described (El Fahaam et al., 1990). Briefly, whole leaves (excluding rachis and petioles) were harvested from tomato plants and freeze-dried in an oven/desiccator at 80°C for 2–3 days or until no change in weight was recorded. Approximately 0.1 g of dry tissue was extracted with 10 ml of 70% hot ethanol and centrifuged at 2,500 g for 5 min. The dry residue was dissolved in 2.5 ml of 0.1 N HCl and kept at −20°C until assayed. Colorimetric procedures were adapted to Technicon Autoanalyzer II for simultaneous determination of total free amino acid content based on an internal standard (leucine) in plant tissue samples <modified fromBroderick and Kang (1980)>. Total free amino acid content data were analyzed using a similar statistical model as the phytohormone analysis; however, biological replicate had a significant main or treatment interaction effect. Therefore, data from biological replicates were interpreted separately.

Results

Performance and Settling Behavior of Thrips

The effect of TSWV infection of tomato plants on thrips performance (number of offspring) was evaluated 7-days after adult females were released onto individual plants. Our data revealed significant differences in the number of offspring (first and second instar larvae) produced on virus-infected plants compared to mock-inoculated plants (F = 11.78, df = 1, P = 0.003) (Figure 1A). There was no impact of time or biological replication (F = 1.73, df = 2, P = 0.21) and the interaction between number of offspring and biological replication (F = 0.03, df = 2, P = 0.97) on thrips population. On average, there were twofold more thrips offspring on virus-infected plants compared to mock-inoculated plants (Mean ± SE: 25.75 ± 2.02 and 13.50 ± 2.67, respectively). Leaf damage caused by thrips was quantified and no apparent differences were detected between thrips on healthy or virus-infected plants (leaf damage index: P = 0.34; number of lesions: P = 0.48, Kruskal–Wallis test;Supplementary Table 1). In Petri dish assays, thrips adults were given a choice between TSWV-infected and mock-inoculated leaflets harvested from the greenhouse experiments. By 3 to 4-h post-release, there were significantly more thrips adults (Z = 10.0, P = 0.03 and Z = 3.0, P = 0.02, respectively) associated with TSWV-infected leaflets compared to mock-inoculated leaflets (Figure 1B). This trend persisted over the course of the 72-h experiment. There were no apparent differences between the numbers of thrips observed between leaflets in the non-choice situation (P > 0.2 for all time-points, data not shown).

[Image Omitted: See PDF]

Single and Combined Effects of TSWV and Thrips on Tomato Global Gene Expression Profiles

Tomato microarray hybridizations were performed to describe and quantify single and combined effects of virus and vector feeding on transcription-level expression. Collectively, of the 10,209 probe sets (i.e., 9,200 unique coding sequences) represented on the Tomato GeneChip, 1,722 sequences were differentially-expressed in plants challenged by the various treatments compared to mock-inoculated, healthy plants (P < 0.05, regardless of magnitude of fold change). Of these sequences, 307, 171, and 424 genes were significantly expressed by at least twofold and P ≤ 0.05 in TSWV, thrips, and TSWV + thrips challenged plants compared to the mock-inoculated controls (Supplementary Tables 2–4, respectively). Venn diagrams depicting the number of unique and shared genes that were differentially-expressed among treatments revealed several patterns (Figures 2A,BandSupplementary Table 5). First, only a small proportion of genes were unique to individual challengers. Second, systemic infection of tomato plants by TSWV contributed the most to global gene expression changes in both the positive (up-regulation) and negative (down-regulation) direction as compared to changes induced by thrips feeding. Third, the majority of genes differentially-expressed in response to thrips alone were down-regulated (64%). Fourth, the combined effect of TSWV and thrips resulted in a greater proportion of down-regulated genes that were unique to dual challenger (i.e., 32% down vs. 16% up) (Figures 2A,B).

[Image Omitted: See PDF]

The differentially-expressed genes were functionally-classified by Gene Ontology (GO) terms into biological process, cellular component and molecular function with relevance to plant health and responses. Overall, a greater proportion of genes responded to TSWV infection alone and combined treatment compared to thrips feeding alone in all three GO categories (Figures 3A–C). The GO-biological process most represented by the differential expression was in response to stimulus, including biotic and abiotic stimulus (Figure 3A). Within this category, a large percentage of genes were up-regulated in response to virus infection alone (74%) and the combination treatment (64%), whereas thrips feeding down-regulated a larger percentage of stress-related genes (63%). The second largest category was photosynthesis, and a majority of photosynthesis-associated genes (90%) were down-regulated in all three treatments. This trend is also reflected in GO cellular component category where chloroplast-associated genes were over-represented, and a large percentage (95%) were down-regulated across all three treatments (Figure 3B). Another category of potential interest was cell wall organization, which was overrepresented in GO-cellular component and GO- biological process as well (Figures 3A,B). Virus infection down-regulated 25% of cell wall related genes, but some genes were also up-regulated (14%). With regards to GO-molecular function, ATP binding genes were over-expressed, which agrees with the GO-CC category where cell membrane-associated genes were overrepresented (Figure 3C). There was no trend in ATP binding genes in response to virus infection, in contrast, thrips feeding resulted in down-regulation of ATP binding genes (100%). The protein binding category was overrepresented, with virus infection in the single and combined treatment resulting in induction of protein binding genes (69 and 87%, respectively), whereas thrips feeding resulted in suppression of the genes (88%).

[Image Omitted: See PDF]

Single and Combined Effects of TSWV and Thrips on Defense-Related and Primary Metabolic Processes Inferred From Gene Expression

Defense Signaling Pathways

We examined the microarray hybridization data for possible interaction between SA and JA pathways. The magnitude and direction of differential expression of all known phytohormone genes in the Tomato GeneChip are shown inSupplementary Table 6. TSWV infection alone significantly up-regulated the majority (59%) of signature SA-responsive genes compared to mock-inoculated plants (Figure 4A). These include genes involved both upstream and downstream of SA synthesis such as chorismate mutase and pathogenesis-related proteins, respectively. Genes encoding the proteins NPR1 and NPR3 that have been shown to interact with SA were also up-regulated in virus-infected plants. We observed suppression of JA-genes in response to virus treatment; 29% suppression and 18% up-regulation in virus-infected plants (Figure 4B). This suggests that the strength of the SA-induced suppression of JA-genes was not as widespread. Moreover, most of the down-regulated JA genes belonged to the protease inhibitor category, specifically wound-induced proteinase inhibitors. A large percentage (60%) of ET-associated genes including those involved in ET synthesis such as 1-aminocyclopropane-1-carboxylate oxidase and ET signaling such as ethylene-responsive transcription factor 5 were up-regulated in TSWV-infected plants (Supplementary Table 6). Feeding by thrips alone resulted in up-regulation of JA-related genes (70%), but no major impact was observed on ET genes (Figure 4BandSupplementary Table 6). Interestingly, genes in all three signaling pathways were significantly up-regulated in TSWV and thrips dual treatment (Supplementary Table 6). In addition to signaling pathways, genes involved in general stress responses such as heat shock proteins, GSTs, and reactive oxygen species (ROS) were differentially-regulated in TSWV-infected tomato, and in most cases up-regulation occurred regardless of the presence of thrips (Supplementary Table 6). We also found transcription factors such as WRKYs, Myb family, bZIP family and Mitogen-activated kinase 6 that initiate defense responses to be up-regulated in response to virus infection.

[Image Omitted: See PDF]

We analyzed other phytohormones that are known to interact with SA and JA/ET, namely ABA and auxin and that were differentially-expressed in our microarray analysis. Virus infection alone did not significantly impact ABA-related genes (66% showed similar expression to mock) and 33% ABA-related genes were down-regulated (Supplementary Table 6). Thrips alone and the dual treatment did not differ significantly in the expression from mock (83 and 60% similarity, respectively) (Supplementary Table 6). With regards to auxin, TSWV infection alone and in combination with thrips down-regulated 40 and 60% of auxin genes, whereas thrips activity had similar percentage of up-regulated and down-regulated genes (40% for both) (Supplementary Table 6).

Photosynthesis-Related Processes

Genes involved in photosynthesis were largely down-regulated in all three treatments (Figure 5AandSupplementary Table 7). Virus infection, both alone and in combination with thrips feeding, had a greater impact on photosynthesis-related genes with 54 and 70% of the genes being suppressed compared to mock-inoculated plants, respectively (Supplementary Table 7). These include ribulose bisphosphate carboxylase, ATP binding protein, chaperone, and chlorophyll a/b-binding proteins. Plants challenged with thrips also repressed 37% of photosynthesis-related genes compared to the control.

[Image Omitted: See PDF]

Cell Wall Organization

Virus infection altered gene expression in both upward and downward direction. Among the up-regulated genes (29%), cell wall degradation genes such as expansin, some cellulases and beta-1, 4-glucanases that are also involved in the SA-pathway (Figure 5BandSupplementary Table 7). In contrast, cell wall-related genes that were significantly down-regulated (25%) were cell wall modification enzymes as xyloglucan endotransglycosylase/hydrolases and pectinesterase. In contrast, thrips feeding alone and combined treatment resulted in up-regulation of a large percentage of cell wall genes (37% for both;Supplementary Table 7). This perturbation of cell wall gene expression shows that thrips feeding has a major impact on cell wall genes even in the absence of virus infection in host tissues.

Protein Synthesis and Degradation

A majority of genes involved in protein synthesis such as constituents of 40S, 50S, and 60S ribosomal subunits and genes associated with protein degradation such as those that encode 26S proteasome and involved in ubiquitination were induced in response to TSWV and the combination treatment (64 and 70%, respectively;Figure 5CandSupplementary Table 8). The presence of thrips alone did not have an impact on expression of protein metabolism genes compared to mock-inoculated plants (Figure 5C).

Hierarchical Cluster Analysis

Focusing on the 369 DE genes associated with pathways of interest, the cluster analysis revealed more similar global gene expression patterns between TSWV infection alone and TSWV + thrips, and likewise, mock-inoculated and thrips alone treatments produced similar patterns (Figure 6andSupplementary Table 8). Genes that were consistently up-regulated by virus alone or the dual treatment were protein metabolism genes such as E3 ubiquitin-ligase, 40S, 50S, and 60S ribosomal subunits (clusters 1, 5, and 6), but also involved genes with broad function in defense and phytohormone pathways, such as NPR1, PR-5 and EDS1 (Supplementary Table 8). In contrast, virus alone and the dual treatment consistently down-regulated photosynthesis and cell wall functional categories such as chlorophyll a-b binding proteins, rhodanese-like domain-containing proteins and xyloglucan endotransglycosylase/hydrolases, respectively (clusters 2, 3, and 4). Interestingly, TSWV alone and the mock treatment showed similar trends in expression (down-regulation) of the genes in cluster 7, which was predominantly comprised of JA and ET pathway genes (Figure 6). Moreover, in this cluster of genes, thrips activity alone resulted in an expression pattern (up-regulation) that was more similar to that of the dual treatment, which included defense and phytohormone related genes. Thrips activity largely down-regulated genes related to protein metabolism, photosynthesis and cell wall functional categories (clusters 1, 4, and 5). These results provide insights into correlation between biochemical pathways, including photosynthesis, protein metabolism, cell wall biogenesis and defense during single and dual attack by TSWV and the thrips vector.

[Image Omitted: See PDF]

Pathway Analysis

Pathway analysis was performed on the microarray data to determine relative activities of five different phytohormone signaling pathways in modulating the TSWV-thrips interaction as perStudham and MacIntosh (2012). The pathway analysis uses q-values for each gene, magnitude and direction of fold change and assigns a role score based on Blast2GO and KEGG functional categories, then sums up their contributions to get the cumulative pathway score. Hence, it differs from an ANOVA that compares log2 (normalized hybridization intensity) of individual genes in the treatment relative to the control (Supplementary Table 6). Four of the five pathways (SA, JA, ABA, and AUX) appeared to be modulated by TSWV-infection alone or in combination with thrips (Table 2), however, plants infected with TSWV alone exhibited an overall negative effect (pathway score = −58.28) on the JA pathway. Thrips infestation on both non-infected and TSWV-infected plants activated the JA pathway, indication that the negative effect of TSWV alone on JA gene expression (Figure 4B) was neutralized by the large positive effect of thrips feeding and/or oviposition on the JA pathway (pathway score = 40.02 and 56.99). Tomato plants singly-challenged by pathogen or pest exhibited an apparent negative co-regulation or cross-talk between the SA and JA pathways, however, under dual challenge, both SA and JA pathways were generally up-regulated, suggesting that other factors, possibly the large ABA effect in the dual treatment (pathway score = 76.36), modulated the SA-JA crosstalk. It was also apparent that virus infection, regardless of thrips infestation, stimulated the ABA pathway and suppressed AUX pathway-associated genes, and had a moderate effect on the ET pathway even after infestation with thrips, indication that the observed thrips-only effect on the ET pathway (Table 2, pathway score = 40.01) was neutralized by TSWV infection (pathway score = 0.64) in the dual treatment. In total, our analysis revealed reciprocal modulation of key phytohormone pathways under dual attack. The phytohormone related genes used for this analysis are listed inSupplementary Table 9.

TABLE 2

Relative activity of phytohormones in the TSWV- Frankliniella occidentalis interaction.

Treatment	Salicylic acid	Jasmonic acid	Ethylene	Abscisic acid	Auxin
TSWV	62.25	–58.28	24.07	27.05	–77.83
Thrips	–38.07	40.02	40.01	–16.8	6.57
TSWV + thrips	37.62	56.99	0.64	76.36	–48.63

Phytohormone pathway scores indicate induction (positive), suppression (negative), or insignificant effect (less than 1) of signaling pathways in response to virus infection, thrips infestation and the combination treatment.

Gene Co-expression Network Analysis

Network analysis was performed on WGCNA-generated modules of co-expressed genes (Supplementary Table 10) to visualize connections among functional categories of interest (color coded) and differentially-expressed, annotated genes (Supplementary Figure 1), with the knowledge that the modular structure of complex networks plays a critical role in their functionality. The result of our network analysis indicated that the gene co-expression network contains five tightly connected groups of gene sequences, whose expression is highly correlated. While photosynthesis is the predominant functional category in all of the modules, each of the groups contain genes from other functional categories, such as plant hormone-related (ABA, AUX, ET, JA, and SA), protein metabolism, cell wall organization and defense functions, thus illustrating the interconnectedness among primary biological processes significantly perturbed by the virus and vector, and co-regulation of genes involved in primary metabolism (photosynthesis, cell wall organization, and protein synthesis) and defense responses.

Reverse Transcription-Quantitative PCR (RT-qPCR) Validation of Microarray Hybridization Data

Six genes associated with the SA, JA, and antiviral small-RNA-mediated gene silencing pathways that were determined to be differentially-expressed by the microarray hybridization experiment were further validated by reverse transcription-quantitative PCR (RT-qPCR) (Table 3). We included RNA-directed RNA polymerase 1 (RDR1), a key siRNA-pathway gene involved in the amplification of virus derived siRNAs that target viral dsRNAs for degradation, to examine antiviral defense. Overall, the average relative expression ratios of SA and JA marker genes, and RDR1 in response to virus infection alone or in combination with thrips mirrored the direction (positive or negative) of expression for the microarray analyses (Table 3). Furthermore, pairwise comparisons of treatment averages obtained for the microarray and real-time RT-qPCR analyses revealed similar patterns among treatments (Table 3). The normalized abundance (NA) of TSWV nucleocapsid (N) and silencing suppressor (NSs) RNA relative to leEF1 was also determined to estimate virus titer in leaf tissue. For virus infected plants (± thrips), there was a significant correlation (Pearson) (r = 0.996, P < 0.0001) between normalized abundance of N RNA and NSs RNA. Analysis of variance revealed significantly lower virus titers in leaves from TSWV-infected plants infested with F. occidentalis compared to plants infected with TSWV alone (Figure 7) (NSs RNA: P = 0.006; N RNA: P = 0.067), indicating that thrips infestation on infected leaves for 1 week significantly influenced virus accumulation. There were significant correlations (n = 6) between TSWV titer and Log₂RER for NPR1 (NSs: r = 0.941, P = 0.005; N: r = 0.935, P = 0.006), and Cathepsin D (CI) inhibitor protein (NSs: r = −0.877, P = 0.022; N: r = −0.897, P = 0.015), an indication that expression of SA- and JA-associated transcripts may be quantitatively associated with the extent of TSWV infection (seeSupplementary Table 11for log2RER values for plant RNA expression (OPR3, AOS, CI, NPR1, BGL, and RDR1) and normalized abundance values for viral RNAs (N, NSs) and the corresponding correlation matrix).

TABLE 3

Reverse transcription quantitative-PCR (RT-qPCR) validation of differential genes.

Log2(fold change intensity)			Log2(relative expression ratio)
Gene title	Probe set ID	TSWV	Thrips	TSWV + thrips	TSWV	Thrips	TSWV + thrips
SA pathway
Beta 1,3 glucanase	Les.3673.1.S1_at	3.6a	−0.9a	3.4a	9.2a	3.3b	7.3a
Non-expressor of pathogenesis-related 1	Les.5940.1.S1_at	0.7a	−0.4b	0.6a	2.3a	−0.2b	1.0a
JA pathway
Allene oxide synthase	Les.13.1.S1_at	−1.4a	0.6a	−0.04a	1.0a	0.3a	1.5a
Cathepsin D inhibitor protein	Les.3740.1.S1_at	−2.4b	0.9a	2.4a	−2.8b	4.2ab	5.1a
12-oxophytodienoate reductase	Les.22.1.S1_at	0.8a	0.2a	0.8a	2.4a	−0.2b	1.4ab
RNAi pathway
RNA-directed RNA polymerase 1	Les.61.1.S1_at	1.1a	−0.04b	1.1a	3.8a	1.2b	3.5a

Relative transcript abundance in tomato plants systemically-infected with TSWV and/or infested with thrips. Values represent mean of three biological replicates. Different letters indicate statistical difference between treatments at P ≤ 0.05.[Image Omitted: See PDF]

TSWV Infection Altered Phytohormone Levels

To determine single and combined effects of virus infection and thrips activity on the levels of the signal molecules (i.e., phytohormones), we quantified the levels of SA, JA, JA-Ile, and OPDA in leaf tissues 7-days post thrips release. Same-age leaflets that were immediately basal (older) to the leaves harvested for microarray analysis were used for phytohormone analysis. TSWV infection alone or in combination with thrips enhanced SA content in tomato leaves (Table 4). Thrips feeding had no apparent effect on any of the phytohormones at the time of sampling; however, the combined effect of TSWV and thrips resulted in significantly higher leaf contents of JA. There were no apparent differences among treatments with regards to JA-Ile or OPDA at the time of sampling (Table 4). TABLE 4

Phytohormone content expressed as Log₁₀(ng analyte g fresh weight^–1) in leaflet samples from tomato plants systemically-infected with TSWV and/or infested with thrips 1 week-post thrips release.

Phytohormone	Mock	TSWV	Thrips	TSWV + thrips
Salicylic acid	1.72 ± 0.12b	3.15 ± 0.18a	1.85 ± 0.13b	2.90 ± 0.22a
Jasmonic acid	0.16 ± 0.06b	0.29 ± 0.14b	0.25 ± 0.11b	0.67 ± 0.14a
JA-isoleucine	0.05 ± 0.48a	0.26 ± 0.36a	0.24 ± 0.50a	0.92 ± 0.36a
OPDA (12-oxo-phytodienoic acid)	0.60 ± 0.39a	1.15 ± 0.34a	0.89 ± 0.34a	1.33 ± 0.31a

Values represent mean ± standard deviation (n = 4 plants). Different letters indicate statistical difference between treatments at P ≤ 0.05.

TSWV Infection Increased Total Free Amino Acid Content

As a measure of plant quality to the insect vector, we measured total free amino acid content to determine if TSWV infection altered the nutritional status of the host. There was a significant effect of time (biological replicate) (F = 10.87, df = 2, P = 0.0002) and treatment (F = 4.25, df = 3, P = 0.01) but not the interaction term (F = 1.11, df = 6, P = 0.37) on the total free amino acid content. Overall, TSWV infection alone and in combination with thrips feeding harbored greater total free amino acid content (Mean ± SE: 288.64 ± 110.87 and 290.12 ± 64.85, respectively) compared to mock-inoculated and thrips –fed plants (Mean ± SE: 164.03 ± 20.56 and 181.55 ± 25.23, respectively) (Figure 8).[Image Omitted: See PDF]

Discussion

Plants encounter and respond to multiple and often co-occurring biotic stresses such as pathogens and insect herbivores in a myriad of ways. Moreover, plant pathogens and insects that share the same host are likely to interact and these interactions may have positive, negative or neutral consequences (Belliure et al., 2005,2010;Thaler et al., 2010;Pan et al., 2013;Eigenbrode et al., 2018). For instance, previous research by us and others showed that a plant virus, TSWV, enhanced survival and oviposition of a non-vector herbivore, spider mites (Belliure et al., 2010;Nachappa et al., 2013). In this study, we explored plant responses that underlie the impact of TSWV on its thrips vector, F. occidentalis. While previous research identified plant defense response as a key molecular mechanism through which TSWV affects its vector (Abe et al., 2012;Wu et al., 2019), other pathways could be involved (Nachappa et al., 2013). Here, we demonstrate that virus infection alters the expression of coordinated networks of regulatory genes controlling primary metabolic pathways and defense responses thereby rendering virus-infected plants more suitable hosts for its insect vectors. To our knowledge, this is the first study to identify global transcriptional networks that modulate TSWV-thrips interaction in tomato and provides information on key pathway players.

Orthotospoviruses depend solely on the thrips vector for transmission (Rotenberg and Whitfield, 2018). Several studies documenting changes in thrips settling behavior/host preference (Maris et al., 2004), feeding behavior (Stafford et al., 2011) or feeding ability (Shrestha et al., 2012) and/or performance (Maris et al., 2004;Belliure et al., 2005;Stumpf and Kennedy, 2005;Inoue and Sakurai, 2006;Stumpf and Kennedy, 2007;Shrestha et al., 2012;Ogada et al., 2013) due to direct or indirect effects of the virus illustrate the possibility of positive and negative outcomes (reviewed inEigenbrode et al. (2018)). In the current study, we found that tomato plants infected with TSWV enhanced performance of thrips compared to healthy or mock-inoculated plants. Previous studies have shown that virus infection increases free amino acid content in infected plants (Markkula and Laurema, 1964;Shrestha et al., 2012), which is known to impact vitellogenesis in thrips and also serve as nutrients for the developing eggs (Klowden, 2013). In agreement with these results, we found that virus infection alone or in combination with thrips resulted in 2.5 and 1.9 times greater total free amino acid content compared to mock-inoculated plants. This suggests that the increased population of thrips is influenced by the increased concentrations of total free amino acids in TSWV-infected plants. In choice assays, we found that TSWV infection increased attractiveness of the host plant for thrips vectors. Recently,Wu et al. (2019)showed that TSWV infection enhances plant attractiveness to the thrips vector by suppressing synthesis of volatile plant terpenes, which is known to repel herbivores including thrips. There were two oligos related to terpene synthesis in the tomato GeneChip, monoterpene synthase 1 and sesquiterpene synthase 1, but they were not differentially-expressed (P < 0.05, fold change >2) between the treatments (data not shown). Taken together, we hypothesize that increased aggregation and population growth of thrips on virus-infected plants can potentially increase the number of viruliferous vectors, which in turn could be expected to increase virus secondary spread by thrips dispersal to neighboring plants. While it is well documented that TSWV spread in various vegetable crops grown in Australia (Coutts et al., 2004) and the south-eastern United States (Camann et al., 1995;Gitaitis et al., 1998) is primarily monocyclic in nature, i.e., primary flight of migratory viruliferous adults settling on plants and slow progression of within-field spread during the season, a 6-year epidemiological study of TSWV and F. occidentalis on processing tomatoes grown in the Central Valley of California points to the likelihood and impact of secondary spread of TSWV by viruliferous adults arising from larvae produced on early-season crops to neighboring late-season tomato crops (Batuman et al., 2020). Consistent with this scenario, attraction of female thrips to TSWV-infected plants would promote oviposition on these plants, and emerging larvae, the requisite stage for TSWV acquisition and subsequent plant inoculation as adults, would thrive on these plants, and at or around crop harvest, migrate as adults to inoculate later season crops.

Because we were primarily interested in separating indirect from direct effects of virus on adult thrips settling preference, we used non-viruliferous females in the current study. However, as stated above, one expected scenario in a field or greenhouse setting is the migration of viruliferous thrips settling on naïve or healthy plants. In this case, one might expect that plants would respond to thrips probing and feeding prior to establishment of a localized, and subsequently systemic viral infection, and this may have implications regarding host response to dual attack. Indeed, there is ample evidence that the type of the attacker (feeding guilds or host specialization) (Inbar et al., 1999;Kaplan and Denno, 2007;Poelman et al., 2008) and order of attack (Inbar et al., 1999;Viswanathan et al., 2007;Poelman et al., 2008;Erb et al., 2011) by two organisms can influence their plant-mediated interactions. A recent meta-analysis of published pest and/or pathogen, plant-mediated interaction studies revealed that in general, attack by a pathogen prior to herbivore attack had no significant effect on herbivore performance or the predicted outcome of JA-SA crosstalk (Moreira et al., 2018), an outcome that differs considerably from our findings. In addition, because very few studies have tested the reciprocal order (herbivore first, pathogen second), this meta-analysis study could not resolve the outcome. As such, the timing and relative occurrence of non-viruliferous and viruliferous thrips relative to orthotospovirus delivery, localized infection and systemic spread, and how these spatiotemporal events may coordinate plant host responses warrants future research to disentangle the complexity and dynamic nature of vector-transmitted plant diseases.

Plants infected with viruses are often found to be more suitable hosts for insect vectors than uninfected plants, specifically those transmitted in a persistent mode [reviewed inEigenbrode et al. (2018)]. Until now, most studies have focused on virus-induced suppression of anti-herbivore defenses (Thaler et al., 2010;Abe et al., 2012;Casteel et al., 2012;Mauck et al., 2019;Wu et al., 2019). A recent study found that plants mount several layers of defense to resist attack from TSWV infection [reviewed inZhu et al. (2019)]. Consistent with these findings, our study demonstrates that virus infection up-regulated a suite of genes related to plant innate immunity and defense response. In the TSWV-thrips interaction, virus infection has been shown to increase the anti-pathogen response (SA-related defenses) which suppresses the anti-herbivore response (JA-related defenses) by exploiting the antagonistic crosstalk between SA-JA plant defenses. This attracts and benefits the vector thrips compared to uninfected plants (Abe et al., 2012). Similarly, we found that virus infection induced a majority of SA-regulated genes and SA content in the plants, but only 29% of the JA-related responses in virus-infected plants were repressed in these plants, which suggests that the inhibition by SA of JA responsive genes is transient. SA-mediated suppression of JA- responsive gene expression is thought to mainly occur downstream of the JA biosynthesis pathway (Leon-Reyes et al., 2010). This may explain the down-regulation of JA-inducible genes such as wound-induced proteinase inhibitors, and the lack of down-regulation in JA biosynthesis genes such as OPR3, LOX AND AOS in TSWV-infected plants. In the current study, we used mechanical leaf-rub inoculation to infect plants with TSWV and thrips were released 2-weeks post virus-inoculation, by which time we expected that wound-related responses to be attenuated. Wounding and insect feeding activate JA-regulated wound response genes such as proteinase inhibitors (Green and Ryan, 1972;Walling, 2000;Wasternack et al., 2006) which could prime the plant to respond more strongly against thrips. However, previous studies found that early wound-response gene RNA levels are up-regulated 0.5 to 2 h after injury and late wound-response gene RNA levels increase from 4 to 8 h (Ryan, 2000;Chen et al., 2008;Scranton et al., 2013). Nevertheless, future experiments could include an undisturbed or healthy control to differentiate between host response to mock-inoculation and thrips feeding. It is likely that timing and magnitude of responses play a major role in orchestrating SA-JA antagonism. For example,Abe et al. (2012)found that the magnitude of suppression of JA-related genes was reduced at 7-day compared to 14-day post-TSWV infection. In the current study, gene expression was measured 3 weeks post-virus-infection, which may be one reason for attenuation of plant responses. Moreover, the host plants were different in the two studies, tomato (current study) versus Arabidopsis (Abe et al., 2012), which may in part explain the inconsistencies between the two studies. Further experimentation utilizing tomato mutants of the SA- and JA-signaling pathways would be useful to explore this result.Koornneef et al. (2008)showed that the antagonistic effect of SA on JA signaling was evident when SA was applied simultaneously with MeJA; however, when SA was applied more than 30 h prior to the onset of the JA response, SA-mediated suppression of JA was not observed. These results suggest that SA-JA crosstalk is transient and depends on the timing, magnitude and order of elicitation. Given the dynamic nature of TSWV infection and symptom development (chlorosis, stunting, and wilt), it is crucial to analyze gene expression profiles during early versus late stages of infection or disease development as it relates to vector performance. While the roles of SA and JA/ET in plant defenses are well-established, these hormones also affect a myriad of development processes such as growth repression, flower development and fertility (JA), flowering time (SA), and germination, senescence and fruit ripening (ET) that may influence outcomes of virus infection. This suggests an intimate interplay exists between phytohormone regulation and primary metabolism [reviewed inRojas et al. (2014)andBerens et al. (2017)].

Other phytohormones such as ET, ABA, auxins, brassinosteroids, cytokinins and gibberellins also play in role in plant-pathogen interactions (De Bruyne et al., 2014;Shigenaga and Argueso, 2016;Yue et al., 2016;Berens et al., 2017;Yang et al., 2019). In the current study, we found that virus infection alone and in combination with thrips up-regulated ABA-related genes and down-regulated auxin-related genes. In general, ABA acts antagonistically with SA, hence it may be beneficial for the virus to induce ABA genes, resulting in suppression of anti-pathogen or SA-mediated defenses (Mauch-Mani and Mauch, 2005;Asselbergh et al., 2008). The role of ABA in response to insect feeding is not clear; however, studies have shown that aphids induce ABA as a decoy strategy to suppress SA- and JA-related defenses (Studham and MacIntosh, 2013;Nachappa et al., 2016). Thrips feeding down-regulated ABA-related genes. Auxin regulates plant growth and development, and is also involved in stress responses [reviewed inKazan and Manners (2009)]. Auxin can affect disease outcomes directly and indirectly. Direct interaction may involve interaction with other phytohormones including SA, JA, and ET. For example, auxin suppresses immunity against the bacterial pathogen Pseudomonas syringae via SA suppression (Wang et al., 2007). Indirect effects of auxin may involve changes in plant growth and development and thereby outcomes of disease resistance. Virus-infected plants often show developmental abnormalities such as stunting and leaf curling, which resemble auxin mutants leading to the conclusion that virus infection alters host auxin signaling. We found that TSWV infection alone and the dual treatment suppressed auxin genes potentially resulting in TSWV-related developmental abnormalities; however, the impact of such developmental alterations on the insect vector is not known.

Although at different time points post-infection, similar trends documented in the transcriptome profiles of tomato and Arabidopsis during TSWV were consistent with our findings (Catoni et al., 2009;Padmanabhan et al., 2019;Xu et al., 2020). We collected plant samples after thrips infestation (3 weeks after virus inoculation) but the other published work suggests that at the time of first thrips deposition on plants, 2 weeks after TSWV-inoculation, these important defense responses and plant barriers were compromised and promoted thrips performance on plants infected with TSWV. Notably, we found that virus infection alone and the combined treatment of virus and thrips down-regulated the majority of the photosynthesis-related genes such as ribulose bisphosphate carboxylase, ATP binding protein, chaperone, and chlorophyll a/b-binding proteins. Although considerable progress has been made in understanding the plant defense responses, very little is known about the role of primary metabolic pathways associated with plant growth and development in regulating defense responses. It is hypothesized that the energy saved by down regulation of primary metabolism specifically photosynthesis and chlorophyll biosynthesis is diverted and used for defense responses (Less et al., 2011;Rojas et al., 2014). Moreover, the chloroplast houses several important steps in the synthesis of phytohormones involved in defense, such as SA, JA, and ABA [reviewed inKangasjärvi et al. (2012)]. Hence, chloroplast-related proteins potentially crosstalk with defense-related proteins. The chloroplast is also the engine of plant growth and plays a crucial role in symptom development. For example, the development of chlorotic or yellowing symptoms of virus-infected plants is likely due to suppression of chloroplast pigment genes. Previous research found that TSWV-infected lettuce plants were more attractive for thrips vectors compared to healthy plants because of the yellow color of the infected plants (Yudin et al., 1987).Catoni et al. (2009)also found that genes related to photosynthesis were suppressed in tomato leaves 14 dpi. More recently,Xu et al. (2020)studied gene expression in response to TSWV infection in Arabidopsis across different development stages (9, 12, and 15 dpi). They too found that genes such as rubisco and Chaperonin and others involved in the photosynthesis pathway were largely suppressed in response to TSWV infection.

In this study, we found that virus infection altered genes related to cell wall organization or biogenesis. Most cell-wall-modification enzymes such as xyloglucan endochitinase, endotransglycosylase/hydrolases and pectinesterase expression levels were decreased in virus-infected plants. In contrast, cell wall degradation genes such as expansin, some cellulases and beta-1, 3-glucanases were upregulated. In agreement with our findings,Xu et al. (2020)found that cellulose synthesis genes were down-regulated whereas cellulases and beta-1,3-glucanases were up-regulated to promote cell wall degradation. Interestingly, beta-1, 3-glucanase is a SA-inducible PR gene that results in SA accumulation during pathogen attack. In a study on transcriptome analysis of a TSWV resistant tomato line (Sw-7) and susceptible line challenged with TSWV infection, the authors found cell wall-related genes to be down-regulated in the resistant line, Sw-7 and only three out of the five pectinesterases were up-regulated (Padmanabhan et al., 2019). The authors also reported up-regulation of pathogenesis-related proteins, PR-1 and PR-5 (osmotin) that modulate callose and lignin deposition in the cell wall leading to restricted virus movement in the resistant line. Indeed, our study also found increased expression of PR-1 and PR-5 in virus-infected plants. These findings provide evidence supporting the dual regulatory role of cell wall organization genes in plant structure and defense. It is possible that suppression of cell wall genes renders the plant more susceptible to the penetration of mouthparts of thrips vectors.

Another functional category of interest was protein synthesis and degradation. Plants are known to accelerate protein metabolism to resist virus infection and spread (Verchot, 2016). In the current study, genes related to protein synthesis and degradation were induced in response to TSWV infection alone and in combination with thrips feeding. Specifically, genes associated with ribosome 60S or 40S subunits were induced as were genes associated with protein degradation such as those that encode 26S proteasome and in the ubiquitination pathway. These results are in agreement with a previous report of increased up-regulation of protein degradation genes during TSWV infection (Xu et al., 2020). Analysis of total free amino acids revealed increased AA content in TSWV-infected plants compared to mock and thrips infested plant. This increase may in part be due to increased protein metabolism in virus-infected plants. It is likely that induction or suppression of specific amino acids is involved in plant defense against pathogens and insects. However, since we did not analyze individual amino acids, we cannot make conclusions about the role of amino acids in defense. There are reports of the involvement of particular amino acids in defense. For example, the lht1 (lysine histidine transporter 1) mutant of Arabidopsis has significantly reduced contents of glutamine, alanine, and proline in comparison with wild-type plants and showed enhanced resistance to diverse bacterial, fungal, and oomycete pathogens (Liu et al., 2010).

The positive effect of TSWV on F. occidentalis described here is consistent with the work of other researchers, but we cannot overlook the broad acting effects of TSWV on non-vector and other trophic levels. Indeed, previous research showed that TSWV infection had positive effects on non-vectors as well (Belliure et al., 2010;Nachappa et al., 2013) by modulating the same pathways namely defense response, photosynthesis, cell wall metabolism and protein metabolism (Nachappa et al., 2013). Hence, changes in these key pathways may also be exploited by non-vector herbivores. Interestingly, one study found that TSWV-infected female thrips were more predaceous on spider mite eggs compared to uninfected thrips (Stafford-Banks et al., 2014). While this behavior is unlikely to increase virus transmission directly, TSWV infection appears to indirectly enhance fitness of thrips vector by improving oviposition of the spider mite prey. In contrast,Pan et al. (2013)found that TSWV infection decreased fitness of whiteflies on peppers. Further research is required to determine whether these mechanisms are different from those identified in our study. Virus infection can also alter plant volatiles that can have consequences for both vector and non-vector herbivores. For example, significantly more aphids settled onto barley yellow dwarf virus-infected than non-infected plants (Jimenez-Martinez et al., 2004). In contrast, fungus gnat adults preferred non-infected plants compared to white clover mosaic virus-infected plants based on their volatile blends (van Molken et al., 2012). These results suggest that the consequences of virus infection on vector and non-vectors and the third trophic level depends on the specific ecological setting, i.e., the species involved, the order in which they attack, the time span of the interaction, the frequency and spatial scales at which pathogens and herbivores co-occur (Roossinck, 2011;van Molken et al., 2012).

In the current study, TSWV appeared to be the prominent driver of plant responses, with some modulation by the thrips vector. With a few specific exceptions, the combination of TSWV and thrips mirrored global gene expression patterns of TSWV infection alone, which suggests that TSWV infection has a significant effect on plant physiology compared to thrips. However, this was likely due, in part, to the length of time allowed for virus accumulation and symptom development (2 weeks) prior to release of female thrips on these plants. Interestingly, we found that thrips – on infected plants for only 1 week prior to sampling – had a negative effect on virus titer as measured by real time qRT-PCR using two viral genes, TSWV-N and NSs. The expression of canonical virus-activated pathways was similar for TSWV alone and the dual treatment leading us to hypothesize that (i) the dual attack may further compromise plant health and thus creates a less suitable host, or (ii) the presence of the insect vector alters the virus composition in the plant host through a yet unidentified pathway or mechanism. Further exploration of this phenomenon is warranted and it will be interesting to determine if thrips perception by the plant alters the virus in a way to promote acquisition similar to transmission morphs described for non-circulative viruses (Gutiérrez et al., 2013;Martinière et al., 2013;Berthelot et al., 2019).

Data Availability Statement

The datasets generated for this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author Contributions

PN, DM, JN, AW, and DR conceived and designed the experiments. PN and DR performed the experiments. PN, DR, and JC analyzed the data. PN, DM, JN, AW, and DR contributed to reagents, materials, and analysis tools. PN and DR wrote the manuscript. All authors contributed to the article and approved the submitted version.

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.

We thank Alexandria Leach-Kieffaber and Vamsi J. Nalam for their assistance with the RT-qPCR experiments.

Supplementary Material

The Supplementary Material for this article can be found online at:https://www.frontiersin.org/articles/10.3389/fpls.2020.575564/full#supplementary-material

Supplementary Figure 1

Weighted gene co-expression network (WGCNA) plotted using the ‘layout_with_fr’ layout in igraph. Each node represents a probe set id (sequence) from the gene chip. The five clusters of nodes following leading eigen community detection are illustrated by the different colored nodes. The membership of nodes in each community is provided in Supplementary Table 10.

Supplementary Table 1

Quantification of thrips leaf feeding damage using a visual rating scale according toMcKinney (1923)and the number of thrips feeding lesions per plant.

Supplementary Table 2

List of differentially-expressed genes (fold change ≥2.0 andP< 0.05) in tomato plants systemically-infected with TSWV. NA means non-annotated, no BLAST hits.

Supplementary Table 3

List of differentially -expressed genes (fold change ≥2.0 andP< 0.05) in tomato plants infested withF. occidentalis. NA means non-annotated, no BLAST hits.

Supplementary Table 4

List of differentially -expressed genes (fold change ≥2.0 andP< 0.05) in tomato plants systemically-infected with TSWV and infested withF. occidentalis. NA means non-annotated, no BLAST hits.

Supplementary Table 5

List of differentially-expressed genes associated with Venn diagram. NA means non-annotated, no BLAST hits.

Supplementary Table 6

List of differentially-expressed genes associated with plant defenses and phytohormone pathways. NA means non-annotated, no BLAST hits. Different letters indicate statistical difference between treatments atP≤ 0.05.

Supplementary Table 7

List of differentially-expressed genes associated with photosynthesis, cell wall organization and protein metabolism. NA means non-annotated, no BLAST hits. Different letters indicate statistical difference between treatments atP≤ 0.05.

Supplementary Table 8

List of differentially-expressed genes associated with specific clusters in the heatmap or hierarchical cluster analysis. NA means non-annotated, no BLAST hits.

Supplementary Table 10

List of differentially-expressed genes used to generate phytohormone pathway scores. NA means non-annotated, no BLAST hits.

Supplementary Table 10

List of differentially-expressed genes associated with each cluster of nodes or community in the WGCNA analysis and gene co-expression network.

Supplementary Table 11

(A)TSWV titer and defense gene expression in leaf tissue from tomato plants infected with TSWV alone or in combination withFrankliniella occidentalis(thrips) from three biological replicate microarray experiments(B)and associated correlation matrix.

References

AbeH.TomitakaY.ShimodaT.SeoS.SakuraiT.KugimiyaS.(2012).Antagonistic plant defense system regulated by phytohormones assists interactions among vector insect, thrips and a Tospovirus.Plant Cell Physiol.53204–212.10.1093/pcp/pcr17322180600

AdieB. A. T.Pérez-PérezJ.Pérez-PérezM. M.GodoyM.Sánchez-SerranoJ.-J.SchmelzE. A.(2007).ABA is an essential signal for plant resistance to pathogens affecting ja biosynthesis and the activation of defenses in Arabidopsis.Plant Cell191665–1681.10.1105/tpc.106.04804117513501

AsselberghB.De VleesschauwerD.HöfteM.(2008).Global switches and fine-tuning-ABA modulates plant pathogen defense.Mol. Plant Microbe Interact.21709–719.10.1094/mpmi-21-6-070918624635

BarabásiA.-L.(2016).Network Science.Cambridge:Cambridge university press.

BariR.JonesJ. D.(2009).Role of plant hormones in plant defence responses.Plant Mol. Biol.69473–488.10.1007/s11103-008-9435-019083153

BatumanO.TuriniT. A.LestrangeM.StoddardS.MiyaoG.AegerterB. J.(2020).Development of an IPM strategy for thrips and tomato spotted wilt virus in processing tomatoes in the Central Valley of California.Pathogens9:636.10.3390/pathogens908063632764311

BautistaR. ÑMauR. F. L.ChoJ. J.CusterD. M.(1995).Potential of tomato spotted wilt tospovirus plant hosts in Hawaii as virus reservoirs for transmission byFranklinieila occidentalis(Thysanoptera: Thripidae).Phytopathology15:32.

BelliureB.JanssenA.MarisP. C.PetersD.SabelisM. W.(2005).Herbivore arthropods benefit from vectoring plant viruses.Ecol. Lett.870–79.10.1111/j.1461-0248.2004.00699.x

BelliureB.SabelisM. W.JanssenA.(2010).Vector and virus induce plant responses that benefit a non-vector herbivore.Basic Appl. Ecol.11162–169.10.1016/j.baae.2009.09.004

BerensM. L.BerryH. M.MineA.ArguesoC. T.TsudaK.(2017).Evolution of hormone signaling networks in plant defense.Annu. Rev. Phytopathol.55401–425.10.1146/annurev-phyto-080516-03554428645231

BerthelotE.DucoussoM.MaciaJ.-L.BogaertF.BaeckerV.ThébaudG.(2019).Turnip mosaic virus is a second example of a virus using transmission activation for plant-to-plant propagation by aphids.J. Virol.93:e01822-18.10.1128/JVI.01822-1830760573

BluaM. J.PerringT. M.MadoreM. A.(1994).Plant virus-induced changes in aphid population development and temporal fluctuations in plant nutrients.J. Chem. Ecol.20691–707.10.1007/bf0205960724242121

BroderickG. A.KangJ. H.(1980).Automated simultaneous determination of ammonia and total amino acids in ruminal fluid and in vitro media.J. Dairy Sci.6364–75.10.3168/jds.s0022-0302(80)82888-8

CamannM.CulbreathA.PickeringJ.ToddJ.DemskiJ.(1995).Spatial and temporal patterns of spotted wilt epidemics in peanut.Phytopathology85879–885.10.1094/phyto-85-879

CarterW.(1939).Populations ofThrips tabaci, with special reference to virus transmission.J. Anim. Ecol.8261–276.10.2307/1234

CasteelC. L.HansenA. K.WallingL. L.PaineT. D.(2012).Manipulation of plant defense responses by the tomato psyllid (Bactericerca cockerelli) and its associated endosymbiont candidatus liberibacter psyllaurous.PLoS One7:e35191.10.1371/journal.pone.003519122539959

CatoniM.MiozziL.FiorilliV.LanfrancoL.AccottoG. P.(2009).Comparative analysis of expression profiles in shoots and roots of tomato systemically infected by Tomato spotted wilt virus reveals organ-specific transcriptional responses.Mol. Plant Microbe Interact.221504–1513.10.1094/mpmi-22-12-150419888816

ChenH.Gonzales-VigilE.HoweG. A.(2008). “Action of plant defensive enzymes in the insect midgut,” inInduced Plant Resistance to Herbivory,ed.SchallerA.(Houten:Springer Verlag),271–284.10.1007/978-1-4020-8182-8_13

ConesaA.GötzS.(2008).Blast2GO: a comprehensive suite for functional analysis in plant genomics.Int. J. Plant Genomics2008:619832.

CouttsB. A.Thomas-CarrollM. L.JonesR. A. C.(2004).Patterns of spread of Tomato spotted wilt virus in field crops of lettuce and pepper: spatial dynamics and validation of control measures.Ann. Appl. Biol.145231–245.10.1111/j.1744-7348.2004.tb00379.x

De BruyneL.HöfteM.De VleesschauwerD.(2014).Connecting growth and defense: the emerging roles of brassinosteroids and gibberellins in plant innate immunity.Mol. Plant7943–959.10.1093/mp/ssu05024777987

DeAngelisJ.SetherD.RossignolP.(1993).Survival, development, and reproduction in western flower thrips (Thysanoptera: Thripidae) exposed to impatiens necrotic spot virus.Environ. Entomol.221308–1312.10.1093/ee/22.6.1308

EigenbrodeS. D.Bosque-PérezN. A.DavisT. S.(2018).Insect-borne plant pathogens and their vectors: ecology, evolution, and complex interactions.Annu. Rev. Entomol.63169–191.10.1146/annurev-ento-020117-04311928968147

EigenbrodeS. D.DingH.ShielP.BergerP. H.(2002).Volatiles from potato plants infected with potato leafroll virus attracts and arrest the virus vector, Myzus persicae (Homoptera:Aphididae).Proc. R. Soc. Lond. B.269455–460.10.1098/rspb.2001.190911886636

El FahaamY. M.FeglaG. I.WagihE. E.El-KaryoniH. A.(1990).Biochemical changes in lettuce plants infected with lettuce mosaic virus.Agric. Sci.2271–277.

ErbM.RobertC. A.HibbardB. E.TurlingsT. C.(2011).Sequence of arrival determines plant−mediated interactions between herbivores.J. Ecol.997–15.10.1111/j.1365-2745.2010.01757.x

FereresA.ListerR. M.ArayaJ. E.FosterJ. E.(1989).Development and reproduction of the English grain aphid (Homoptera:Aphididae) on wheat cultivars infected with Barley yellow dwarf virus.Environ. Entomol.18388–393.10.1093/ee/18.3.388

FruchtermanT. M.ReingoldE. M.(1991).Graph drawing by force−directed placement.Softw. Pract. Exp.211129–1164.10.1002/spe.4380211102

GitaitisR.DowlerC.ChalfantR.(1998).Epidemiology of tomato spotted wilt in pepper and tomato in southern Georgia.Plant Dis.82752–756.10.1094/pdis.1998.82.7.75230856944

GreenT. R.RyanC. A.(1972).Wound-induced proteinase inhibitor in plant leaves: a possible defense mechanism against insects.Science175776–777.10.1126/science.175.4023.77617836138

GuoJ.-Y.YeG.-Y.DongS.-Z.LiuS.-S.(2010).An invasive whitefly feeding on a virus-infected plant increased its egg production and realized fecundity.PLoS One5:e11713.10.1371/journal.pone.001171320676356

GutiérrezS.MichalakisY.Van MunsterM.BlancS.(2013).Plant feeding by insect vectors can affect life cycle, population genetics and evolution of plant viruses.Funct. Ecol.27610–622.10.1111/1365-2435.12070

HogenhoutS. A.AmmarE.-D.WhitfieldA. E.RedinbaughM. G.(2008).Insect vector interactions with persistently transmitted viruses.Annu. Rev. Phytopathol.46327–359.10.1146/annurev.phyto.022508.09213518680428

HoweG. A.JanderG.(2008).Plant immunity to insect herbivores.Annu. Rev. Plant Biol.5941–46.10.1146/annurev.arplant.59.032607.09282518031220

InbarM.DoostdarH.LeibeeG. L.MayerR. T.(1999).The role of plant rapidly induced responses in asymmetric interspecific interactions among insect herbivores.J. Chem. Ecol.251961–1979.

IngwellL. L.EigenbrodeS. D.Bosque-PérezN. A.(2012).Plant viruses alter insect behavior to enhance their spread.Sci. Rep.2:578.

InoueT.SakuraiT.(2006).Infection of Tomato spotted wilt virus (TSWV) shortens the life span of thelytokous Thrips tabaci (Thysanoptera: Thripidae).Appl. Entomol. Zool.41239–246.10.1303/aez.2006.239

Jimenez-MartinezE. S.Bosque-PerezN. A.BergerP. H.ZemetraR. S.(2004).Life history of the bird cherry-oat aphid, Rhopalosiphum padi (Homoptera:Aphididae), on transgenic and untrnasformed wheat challenged with Barley yellow dwarf virus.J. Ecol. Entomol.97203–212.10.1093/jee/97.2.20315154437

KanehisaM.GotoS.SatoY.FurumichiM.TanabeM.(2011).KEGG for integration and interpretation of large-scale molecular data sets.Nucleic Acids Res.40D109–D114.

KangasjärviS.NeukermansJ.LiS.AroE.-M.NoctorG.(2012).Photosynthesis, photorespiration, and light signalling in defence responses.J. Exp. Bot.631619–1636.10.1093/jxb/err40222282535

KaplanI.DennoR. F.(2007).Interspecific interactions in phytophagous insects revisited: a quantitative assessment of competition theory.Ecol. Lett.10977–994.10.1111/j.1461-0248.2007.01093.x17855811

KazanK.MannersJ. M.(2009).Linking development to defense: auxin in plant–pathogen interactions.Trends Plant Sci.14373–382.10.1016/j.tplants.2009.04.00519559643

KlowdenM. J.(2013).Physiological Systems in Insects.Cambridge, MA:Academic press.

KoornneefA.Leon-ReyesA.RitsemaT.VerhageA.Den OtterF. C.Van LoonL.(2008).Kinetics of salicylate-mediated suppression of jasmonate signaling reveal a role for redox modulation.Plant Physiol.1471358–1368.10.1104/pp.108.12139218539774

LangfelderP.HorvathS.(2008).WGCNA: an R package for weighted correlation network analysis.BMC Bioinform.9:559.10.1186/1471-2105-9-55919114008

Leon-ReyesA.Van der DoesD.De LangeE.DelkerC.WasternackC.Van WeesS.(2010).Salicylate-mediated suppression of jasmonate-responsive gene expression in Arabidopsis is targeted downstream of the jasmonate biosynthesis pathway.Planta2321423–1432.10.1007/s00425-010-1265-z20839007

LessH.AngeloviciR.TzinV.GaliliG.(2011).Coordinated gene networks regulating Arabidopsis plant metabolism in response to various stresses and nutritional cues.Plant Cell231264–1271.10.1105/tpc.110.08286721487096

LiuG.JiY.BhuiyanN.PilotG.SelvarajG.ZouJ.(2010).Amino acid homeostasis nodulates salicylic acid–associated redox status and defense responses in arabidopsis.Plant Cell223845–3863.10.1105/tpc.110.07939221097712

MarisP. C.JoostenN. N.GoldbachR. W.PetersD.(2004).Tomato spotted wilt virus infection improves host suitability for its vector Frankliniella occidentalis.Phytopathology94706–711.10.1094/PHYTO.2004.94.7.70618943902

MarkkulaM.LauremaS.(1964).Changes in the concentration of free amino acids in plants induced by virus disease and the reproduction of aphids.Ann. Agric. Fenn.3265–271.

MartinièreA.BakA.MaciaJ.-L.LautredouN.GarganiD.DoumayrouJ.(2013).A virus responds instantly to the presence of the vector on the host and forms transmission morphs.eLife2:e00183.

Mauch-ManiB.MauchF.(2005).The role of abscisic acid in plant–pathogen interactions.Curr. Opin. Plant Biol.8409–414.10.1016/j.pbi.2005.05.01515939661

MauckK. E.De MoraesC. M.MescherM. C.(2010).Deceptive chemical signals induced by a plant virus attract insect vectors to inferior hosts.Proc. Natl. Acad. Sci. U.S.A.1073600–3605.10.1073/pnas.090719110720133719

MauckK. E.KenneyJ.ChesnaisQ.(2019).Progress and challenges in identifying molecular mechanisms underlying host and vector manipulation by plant viruses.Curr. Opin. insect Sci.337–18.10.1016/j.cois.2019.01.00131358199

McKenzieC.(2002).Effect of tomato mottle virus (ToMoV) on Bemisia tabaci biotype B (Homoptera: Aleyrodidae) oviposition and adult survivorship on healthy tomato.Florida Entomol.85367–369.10.1653/0015-4040(2002)085[0367:eotmvt]2.0.co;2

McKinneyH.(1923).A new system of grading plant diseases.J. Agric. Res26195–218.

MelottoM.UnderwoodW.KoczanJ.NomuraK.HeS. Y.(2006).Plant stomata function in innate immunity against bacterial invasion.Cell126969–980.

Montero-AstúaM.UllmanD. E.WhitfieldA. E.(2016).Salivary gland morphology, tissue tropism and the progression of tospovirus infection in Frankliniella occidentalis.Virology49339–51.10.1016/j.virol.2016.03.00326999025

MoreiraX.Abdala−RobertsL.CastagneyrolB.(2018).Interactions between plant defence signalling pathways: Evidence from bioassays with insect herbivores and plant pathogens.J. Ecol.1062353–2364.10.1111/1365-2745.12987

NachappaP.CulkinC. T.SayaP. M.HanJ.NalamV. J.(2016).Water stress modulates soybean aphid performance, feeding behavior and virus transmission in soybean.Front. Plant Sci.7:552.10.3389/fpls.2016.0055227200027

NachappaP.MargoliesD. C.NecholsJ. R.WhitfieldA. E.RotenbergD.(2013).Tomato spotted wilt virus benefits a non-vector arthropod, Tetranychus urticae, by modulating different plant responses in tomato.PLoS One8:e75909.10.1371/journal.pone.007590924058708

NewmanM. E.(2006).Finding community structure in networks using the eigenvectors of matrices.Phys. Rev. E74:036104.

NgumbiE.EigenbrodeS. D.Bosque-PerezN. A.DingH.RodriguezA.(2007).Myzus persicae is arrested more by blends than by individual compounds elevated in headspace of PLRV-infected potato.J. Chem. Ecol.331733–1747.10.1007/s10886-007-9340-z17680312

OgadaP. A.MaissE.PoehlingH. M.(2013).Influence of tomato spotted wilt virus on performance and behaviour of western flower thrips (Frankliniella occidentalis).J. Appl. Entomol.137488–498.10.1111/jen.12023

PadmanabhanC.MaQ.ShekastebandR.StewartK. S.HuttonS. F.ScottJ. W.(2019).Comprehensive transcriptome analysis and functional characterization of PR-5 for its involvement in tomato Sw-7 resistance to tomato spotted wilt tospovirus.Sci. Rep.9:7673.10.1038/s41598-019-44100-x31114006

PanH.ChenG.LiF.WuQ.WangS.XieW.(2013).Tomato spotted wilt virus infection reduces the fitness of a nonvector herbivore on pepper.J. Econ. Entomol.106924–928.10.1603/ec1236523786083

PanX.WeltiR.WangX.(2008).Simultaneous quantifictaion of major phytohormones and related compounds in crude plant extracts by liquid chromatography-electrospray.Phytochemistry691773–1781.10.1016/j.phytochem.2008.02.00818367217

PfafflM. W.(2001).A new mathematical model for relative quantification in real-time RT-PCR.Nucleic Acids Res.29:e45.

PieterseC.ReyesA.-L.Van Der EntS.Van WeesS.(2009).Networking by small-molecule hormones in plant immunity.Nat. Chem. Ecol.5308–316.10.1038/nchembio.116423314914

PoelmanE. H.GaliartR.RaaijmakersC. E.Van LoonJ. J. A.Van DamN. M.(2008).Performance of specialist and generalist herbivores feeding on cabbage cultivars is not explained by glucosinolate profiles.Entomol. Exp. Appl.127218–228.10.1111/j.1570-7458.2008.00700.x

RojasC.Senthil-KumarM.TzinV.MysoreK.(2014).Regulation of primary plant metabolism during plant-pathogen interactions and its contribution to plant defense.Front. Plant Sci.5:17.10.3389/fpls.2014.0001724575102

RoossinckM. J.(2011).The good viruses: viral mutualistic symbioses.Nat. Rev. Microbiol.999–108.10.1038/nrmicro249121200397

RotenbergD.Krishna KumarN. K.UllmanD. E.Montero-AstúaM.WillisD. K.GermanT. L.(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-040419271982

RotenbergD.ThompsonT. S.GermanT. L.WillisD. K.(2006).Methods for effective real-time RT-PCR analysis of virus-induced gene silencing.J. Virol. Methods13849–59.10.1016/j.jviromet.2006.07.01716959330

RotenbergD.WhitfieldA. E.(2018).Molecular interactions between tospoviruses and thrips vectors.Curr. Opin. Virol.33191–197.10.1016/j.coviro.2018.11.00730500681

RyanC. A.(2000).The systemin signaling pathway: differential activation of plant defensive genes.Biochim. Biophys. Acta1477112–121.10.1016/s0167-4838(99)00269-1

RybickiE. P.(2015).A Top Ten list for economically important plant viruses.Arch. Virol.16017–20.10.1007/s00705-014-2295-925430908

ScrantonM. A.FowlerJ. H.GirkeT.WallingL. L.(2013).Microarray analysis of tomato’s early and late wound response reveals new regulatory targets for leucine aminopeptidase A.PLoS One8:e77889.10.1371/journal.pone.007788924205013

ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.(2003).Cytoscape: a software environment for integrated models of biomolecular interaction networks.Genome Res.132498–2504.10.1101/gr.123930314597658

ShigenagaA. M.ArguesoC. T.(2016).No hormone to rule them all: interactions of plant hormones during the responses of plants to pathogens.Semin. Cell Dev. Biol.56174–189.10.1016/j.semcdb.2016.06.00527312082

ShresthaA.SrinivasanR.RileyD. G.CulbreathA. K.(2012).Direct and indirect effects of a thrips-transmitted Tospovirus on the preference and fitness of its vector,Frankliniella fusca.Entomol. Exp. Appl.145260–271.10.1111/eea.12011

StaffordC. A.WalkerG. P.UllmanD. E.(2011).Infection with a plant virus modifies vector feeding behavior.Proc. Natl. Acad. Sci U.S.A1089350–9355.10.1073/pnas.110077310821606372

Stafford-BanksC. A.YangL. H.McmunnM. S.UllmanD. E.(2014).Virus infection alters the predatory behavior of an omnivorous vector.Oikos1231384–1390.10.1111/oik.01148

StoutM. J.ThalerJ. S.ThommaB. P. H. J.(2006).Plant-mediated interactions between pathogenic microorganisms and herbivorous arthropods.Annu. Rev. Entomol.51663–689.10.1146/annurev.ento.51.110104.15111716332227

StudhamM. E.MacIntoshG. C.(2012).Phytohormone signaling pathway analysis method for comparing hormone responses in plant-pest interactions.BMC Res. Notes5:392.10.1186/1756-0500-5-39222846705

StudhamM. E.MacIntoshG. C.(2013).Multiple phytohormone signals control the transcriptional response to soybean aphid infestation in susceptible and resistant soybean plants.Mol. Plant Microbe Interact.26116–129.10.1094/mpmi-05-12-0124-fi22992001

StumpfC. F.KennedyG. G.(2005).Effects of tomato spotted wilt virus (TSWV) isolates, host plants, and temperature on survival, size, and development time ofFrankliniella fusca.Entomol. Exp. Appl.114215–225.10.1111/j.1570-7458.2005.00251.x

StumpfC. F.KennedyG. G.(2007).Effects of tomato spotted wilt virus isolates, host plants, and temperature on survival, size, and development time ofFrankliniella occidentalis.Entomol. Exp. Appl.123139–147.10.1111/j.1570-7458.2007.00541.x

ThalerJ. S.AgrawalA. A.HalitschkeR.(2010).Salicylate-mediated interactions between pathogens and herbivores.Ecology911075–1082.10.1890/08-2347.1

TonJ.Mauch-ManiB.(2004).β-amino-butyric acid-induced resistance against necrotrophic pathogens is based on ABA-dependent priming for callose.Plant J.38119–130.10.1111/j.1365-313x.2004.02028.x15053765

UllmanD. E.ChoJ. J.MauR. F.WestcotD. M.CusterD. M.(1992).A midgut barrier to tomato spotted wilt virus acquisition by adult western flower thrips.Phytopathology821333–1333.10.1094/phyto-82-1333

Van De WeteringF.GoldbachR.PetersD.(1996).Tomato spotted wilt tospovirus ingestion by first instar larvae ofFrankliniella occidentalisis a prerequisite for transmission.Phytopathology86900–905.10.1094/phyto-86-900

van MolkenT.De CaluweH.HordijkC. A.Leon-ReyesA.SnoerenT. A.Van DamN. M.(2012).Virus infection decreases the attractiveness of white clover plants for a non-vectoring herbivore.Oecologia170433–444.10.1007/s00442-012-2322-z22526939

VerchotJ.(2016).Plant virus infection and the ubiquitin proteasome machinery: arms race along the endoplasmic reticulum.Viruses8:314.10.3390/v811031427869775

VerhageA.Van WeesS.PieterseC.(2010).Plant immunity:Its the hormone talking, but what do they say?Plant Physiol.154536–540.10.1104/pp.110.16157020921180

ViswanathanD.LifchitsO.ThalerJ.(2007).Consequences of sequential attack for resistance to herbivores when plants have specific induced responses.Oikos1161389–1399.10.1111/j.0030-1299.2007.15882.x

WallingL.(2000).The myraid of plant responses to herbivores.J. Plant Growth Regul.19195–216.10.1007/s00344000002611038228

WangD.Pajerowska-MukhtarK.CullerA. H.DongX.(2007).Salicylic acid inhibits pathogen growth in plants through repression of the auxin signaling pathway.Curr. Biol.171784–1790.10.1016/j.cub.2007.09.02517919906

WasternackC.StenzelI.HauseB.HauseG.KutterC.MaucherH.(2006).The wound response in tomato – Role of jasmonic acid.J. Plant Physiol.163297–306.10.1016/j.jplph.2005.10.01416368162

WhitfieldA. E.FalkB. W.RotenbergD.(2015).Insect vector-mediated transmission of plant viruses.Virology479-480278–289.10.1016/j.virol.2015.03.02625824478

WhitfieldA. E.UllmanD. E.GermanT. L.(2005).Tospovirus–thrips interactions.Annu. Rev. Phytopathol.43459–489.10.1146/annurev.phyto.43.040204.14001716078892

WieseJ.KranzT.SchubertS.(2004).Induction of pathogen resistance in barley by abiotic stress.Plant Biol.6529–536.10.1055/s-2004-82117615375723

WijkampI.GoldbachR.PetersD.(1996).Propagation of tomato spotted wilt virus inFrankliniella occidentalisdoes neither result in pathological effects nor in transovarial passage of the virus.Entomol. Exp. Appl.81285–292.10.1046/j.1570-7458.1996.00098.x

WuX.XuS.ZhaoP.ZhangX.YaoX.SunY.(2019).The Orthotospovirus nonstructural protein NSs suppresses plant MYC-regulated jasmonate signaling leading to enhanced vector attraction and performance.PLoS Pathog.15:e1007897.10.1371/journal.ppat.100789731206553

XuM.ChenJ.HuangY.ShenD.SunP.XuY.(2020).Dynamic transcriptional profiles ofArabidopsis thalianaInfected by Tomato spotted wilt virus.Phytopathology110153–163.10.1094/phyto-06-19-0199-fi31544594

YangJ.DuanG.LiC.LiuL.HanG.ZhangY.(2019).The crosstalks between jasmonic acid and other plant hormone signaling highlight the involvement of jasmonic acid as a core component in plant response to biotic and abiotic stresses.Front. Plant Sci.10:1349.10.3389/fpls.2019.0134931681397

YudinL. S.MitchellW. C.ChoJ. J.(1987).Color preference of thrips (Thysanoptera:Thripidae) with reference to aphids (Homoptera: Aphididae) and leafminers in Hawaiian lettuce farms.J. Econ. Entomol.8051–55.10.1093/jee/80.1.51

YueX.LiX. G.GaoX.-Q.ZhaoX. Y.DongY. X.ZhouC.(2016).The Arabidopsis phytohormone crosstalk network involves a consecutive metabolic route and circular control units of transcription factors that regulate enzyme-encoding genes.BMC Syst. Biol.10:87.10.1186/s12918-016-0333-927590055

ZhuM.Van GrinsvenI. L.KormelinkR.TaoX.(2019).Paving the way to tospovirus infection: multilined interplays with plant innate immunity.Annu. Rev. Phytopathol.5741–62.10.1146/annurev-phyto-082718-10030930893008

ZüstT.AgrawalA. A.(2017).Trade-offs between plant growth and defense against insect herbivory: an emerging mechanistic synthesis.Annu. Rev. Plant Biol.68513–534.10.1146/annurev-arplant-042916-04085628142282