Transcriptome Proﬁling Unravels the Involvement of Phytohormones in Tomato Resistance to the Tomato Yellow Leaf Curl Virus (TYLCV)

: Tomato yellow leaf curl virus (TYLCV) is a serious pathogen transmitted by the whiteﬂy ( Bemisia tabaci ). Due to the quick spread of the virus, which is assisted by its vector, tomato yield and quality have suffered a crushing blow. Resistance to TYLCV has been intensively investigated in transmission, yet the mechanism of anti-TYLCV remains elusive. Herein, we conducted transcriptome proﬁling with a TYLCV-resistant cultivar (CLN2777A) and a susceptible line (Moneymaker) to identify the potential mechanism of resistance to TYLCV. Compared to the susceptible line, CLN2777A maintained a lower level of lipid peroxidation (LPO) after TYLCV infection. Through RNA-seq, over 1000 differentially expressed genes related to the metabolic process, cellular process, response to stimulus, biological regulation, and signaling were identiﬁed, indicating that the defense response was activated after the virus attack. Further analysis showed that TYLCV infection could induce the expression of the genes involved in salicylic and jasmonic acid biosynthesis and the signal transduction of phytohormones, which illustrated that phytohormones were essential for tomatoes to defend against TYLCV. These ﬁndings provide greater insight into the effective source of resistance for TYLCV control, indicating a potential molecular tool for the design of TYLCV-resistant tomatoes. results indicate that tomato plants of T2 probably regulate the biosynthesis of SA to defend against pathogen attack, making sure that plants produce more ZT and IAA for growth.


Introduction
Tomato yellow leaf curl virus (TYLCV), belonging to the genus Begomovirus (family Geminiviridae), is a devastating pathogen and can lead to great crop loss, especially in tropical and subtropical regions with high heat and humidity [1]. The migration of infected whitefly (Bemisia tabaci) is considered the main route of transmission. However, planting or grafting seedlings with TYLCV could also spread the virus [2][3][4][5]. In addition to vegetables such as tomatoes, peppers, and cucumbers, TYLCV can infect various crops, including tobacco, cassava, and cotton, and can cause great damage to their production [6]. As the main host of TYLCV, the typical symptom in tomatoes induced by the virus is arrested development with upward curling of the leaves and chlorosis, which results in a reduced yield and poor quality [7,8]. It is hard to completely eliminate TYLCV or control an outbreak of it as its vector, Bemisia tabaci, has a strong reproduction capability and migration rate. Hence, characterization of the TYLCV resistance gene to select and breed cultivars against the virus may be an ideal strategy to defend against the pathogen [2,3,9,10].
Plant hormones are a key part of plant-pathogen interactions. Exogenous application of salicylic acid (SA) can inhibit virus accumulation and improve the resistance of susceptible cultivars. Furthermore, the blocking of SA signaling results in more detectable TYLCV in tomatoes [32,33]. Jasmonic acid (JA) can regulate the biosynthesis of phytoalexin to avoid the virus transmitting from the vector to the plant [34]. The virulence factor βC1 of TYLCV directly interacts with MYC2, the fundamental transcription factor in JA signaling, and compromises its transcription activation, resulting in the suppression of terpene biosynthesis and subverting plant resistance [35]. To date, researchers have mainly focused on gene expression and function with samples after a long period of infection in TYLCV and tomato interaction [36][37][38]. To investigate the differential expression genes induced by TYLCV in a short time, we carried out RNA-seq with Moneymaker (MK; a susceptible line) and CLN2777A (T2; a resistant line) infecting them with virus after 24 and 48 h. Meanwhile, the activity of reactive oxygen species (ROS)-related enzymes and plant hormones was detected to uncover the mechanism of TYLCV resistance in tomatoes.

Plant Materials and TYLCV Inoculation
Moneymaker (MK; susceptible to TYLCV infection) and CLN2777A (T2; harboring the Ty-2 gene, resistant to TYLCV infection) were used for TYLCV inoculation. The germplasms of MK and T2 were obtained from Jiangsu Academy of Agricultural Sciences. The seeds were germinated at 25 • C, and then put in plastic boxes containing nutritional soil (soil/perlite/vermiculite = 2:1:1, v/v/v). Seedlings of the two-leaf stage were transferred into plastic pots. The temperature and relative humidity for all of the experiments were set at 25 • C and 70%-80% with a 16/8 h day/night photoperiod. For TYLCV inoculation, we used an Agrobacterium tumefaciens vector constructed by Zhang et al. [39], which contained the effective infection fragment of TYLCV. The method of inoculation of TYLCV was performed as described by Zhang et al. [39] with slight modifications. The first leaf of four-leaf stage plants was injected with 1 mL of A. tumefaciens (OD 600 = 0.6) harboring TYLCV or not (control (CK)). After inoculating with TYLCV for 24 and 48 h, the leaves of plants were harvested for ROS-related enzyme activity determination, RNA sequencing, phytohormones assay, and gene expression detection. Thirty seedlings were planted for each treatment of each tomato line, and the leaves of six random plants were collected at designated time points; then, all samples were frozen in liquid nitrogen and immediately stored at −80 • C.

Growth Parameter Measurement
After 26 d of TYLCV inoculation, plant height and fresh and dry weight (shoots and roots) were measured. The fresh weight was measured on a scale (±0.01 g) with an ambient temperature of 25 • C. Then, the samples were dried at 80 • C in an oven until the weight no longer changed, and the dry weight was measured on a scale (±0.01 g). For each measurement, one biological replicate containing six seedlings and three biological replicates was performed.

Enzyme Activity Assays
For detection of the enzyme activities, tomato leaves (0.2 g) were ground to homogenate in 1 mL of ice-cold phosphate buffer (pH = 7.8, 50 mM), and then centrifuged at 12,000× g for 20 min at 4 • C. SOD and POD activity was measured according to Macadam [40]. The activity of ascorbate oxidase (AAO) was detected by the method from Nakano [41]. Three biological replicates were performed for each measurement.

Library Construction and Illumina Sequencing
For library preparation, RNA of leaf samples collected at 0, 24, and 48 h post-TYLCV inoculation (hpi) were isolated, respectively, and the preparation strategy can be found in previous study [42]. First, total RNA was extracted using a Trizol reagent kit (Invitrogen, Carlsbad, CA, USA), and the extracted RNA concentration and quality were assessed by Agilent 2100 Bioanalyzer (Agilent Technologies, Palo Alto, CA, USA). Next, the eukaryotic mRNA was enriched by Oligo (dT)-attached magnetic beads, while the ribosomal RNA was removed by Ribo-zeroTM rRNA removal kit (Epicentre, Madison, WI, USA). Then the enriched mRNA was fragmented into short fragments using fragmentation buffer and reverse-transcribed into double-strand cDNA with random primers. The libraries were constructed using TruSeq RNA sample Preparation kit (Illumina, San Diego, CA, USA) following the manufacturer's recommendations, and then sequenced on Illumina HiSeq 2500 platform with 2 × 150 bp PE mode by Gene Denovo Biotechnology Co. (Guangzhou, China).

Read Mapping, Statistical Analysis, and Functional Annotation
Useless sequences such as sequencing adapters, rRNA reads, and short fragment reads were filtered out. Then, the clean reads were mapped to the Solanum lycopersicum V4.0 reference genome (https://solgenomics.net/organism/Solanum_lycopersicoides/ genome) on 22 April 2021 by HISAT2 software [43] (version 2.2.4, Johns Hopkins University, Baltimore, MD, USA, http://daehwankimlab.github.io/hisat2/, 22 April 2021) with "-rnastrandness RF" and other parameters set as a default. After this, the FPKM (fragments per kilobase of transcript per million mapped reads) values for known gene models were identified by StringTie software (version 1.3.1, Johns Hopkins University, Baltimore, MD, USA, http://ccb.jhu.edu/software/stringtie/index.shtml on 22 April 2021) according to Pertea et al. [44]. The threshold of the p-value was determined using the false discovery rate (FDR). In this study, differentially expressed genes (DEGs) were identified using edgeR [45], with the rules FDR ≤ 0.05 and fold change ≥2. The GO annotation of DEGs was analyzed by Blast2go software (http://www.blast2go.org) on 26 April 2021. KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway enrichment analysis was performed by a KEGG enrichment tool of the platform OmicShare (http://www.omicshare.com/ tools/Home/Soft/pathwaygsea) on 28 April 2021.

Phytohormone Content Assays
The phytohormone contents were detected as previously described [46], and the detection was performed at Genepioneer Biotechnologies Co. Ltd, Nanjing., China. Briefly, 0.2 g of the leaf sample was ground and extracted with a mixture buffer (2-propanol/H 2 O/ concentrated HCl, 2:1:0.002, v/v/v). The sample solution was loaded into the reverse-phase C18 Gemini HPLC column for HPLC-ESI-MS/MS analysis.

Validation of DEGs with Quantitative RT-PCR
Ten genes involved in the TYLCV resistance of tomatoes were selected to conduct quantitative RT-PCR validation. Primers were designed and their specificity for quantitative RT-PCR were evaluated using NCBI primer-blast (Table S1). First-strand cDNA synthesis, PCR amplifications, and calculations of the relative gene expression level were performed as previously described [47]. The gene α-Tubulin (Solyc04g077020.2) was used as the internal reference. Three independent biological replicates and two technical repeats per biological replicate were performed for each chosen gene.

TYLCV Leads to Tomato Development Arrest
To further investigate the influence of the virus on tomato growth, the TYLCVsusceptible line MK and the resistant line T2 were infected with the virus. The biomass, including plant weight and fresh and dry weight, was measured after 26 d of infection. The plant height of MK was reduced by 20%, while no significant difference was noticed in T2 invaded by TYLCV ( Figure 1A). Similarly, both the fresh and dry weights of MK infected by TYLCV were significantly decreased ( Figure 1B,C). In contrast, the decrease in the biomass of T2 showed no significant difference ( Figure 1B,C). These data suggest that TYLCV can repress tomato growth and reduce biomass.
The phytohormone contents were detected as previously described [46], and the detection was performed at Genepioneer Biotechnologies Co. Ltd, Nanjing., China. Briefly, 0.2 g of the leaf sample was ground and extracted with a mixture buffer (2-propanol/H2O/concentrated HCl, 2:1:0.002, v/v/v). The sample solution was loaded into the reverse-phase C18 Gemini HPLC column for HPLC-ESI-MS/MS analysis.

Validation of DEGs with quantitative RT-PCR
Ten genes involved in the TYLCV resistance of tomatoes were selected to conduct quantitative RT-PCR validation. Primers were designed and their specificity for quantitative RT-PCR were evaluated using NCBI primer-blast (Table S1). First-strand cDNA synthesis, PCR amplifications, and calculations of the relative gene expression level were performed as previously described [47]. The gene α-Tubulin (Solyc04g077020.2) was used as the internal reference. Three independent biological replicates and two technical repeats per biological replicate were performed for each chosen gene.

TYLCV leads to tomato development arrest
To further investigate the influence of the virus on tomato growth, the TYLCV-susceptible line MK and the resistant line T2 were infected with the virus. The biomass, including plant weight and fresh and dry weight, was measured after 26 d of infection. The plant height of MK was reduced by 20%, while no significant difference was noticed in T2 invaded by TYLCV (Figure 1a). Similarly, both the fresh and dry weights of MK infected by TYLCV were significantly decreased (Figure 1b, c). In contrast, the decrease in the biomass of T2 showed no significant difference (Figure 1b, c). These data suggest that TYLCV can repress tomato growth and reduce biomass. infected by TYLCV for 26 d. The data represent the mean ± SD (n = 3). Data were analyzed by an ANOVA in the SPSS software. * (p < 0.05) and ** (p < 0.01) indicate statistically significant differences.

TYLCV leads to lipid peroxidation stress at the initial time of infection
The accumulation of ROS is a stress reaction of a plant to disease. The amount of ROS could be elevated in tomatoes, breaking the state of its production and degradation balance under virus treatment for 2-14 d [33]. However, few explorations have focused on the initial time of symptom appearance in plants attacked by TYLCV. We checked the content of MDA along with the ROS-scavenging enzymes in response to TYLCV infection at early times (within 48 h). The results showed that more MDA accumulated from 24 h in the TYLCV-sensitive cultivar MK (S), while the process was delayed in T2 (R) after virus infection ( Figure 2a). Furthermore, the MDA content in TYLCV-attacked T2 maintained half of that in MK, which suggested that the virus-resistant lines could alleviate TYLCVinduced lipid peroxidation. All of the ROS-scavenging enzymes had higher activity in T2 than those in MK with the virus treatment at the same time point. Remarkably, the activity infected by TYLCV for 26 d. The data represent the mean ± SD (n = 3). Data were analyzed by an ANOVA in the SPSS software. * (p < 0.05) and ** (p < 0.01) indicate statistically significant differences.

TYLCV Leads to Lipid Peroxidation Stress at the Initial Time of Infection
The accumulation of ROS is a stress reaction of a plant to disease. The amount of ROS could be elevated in tomatoes, breaking the state of its production and degradation balance under virus treatment for 2-14 d [33]. However, few explorations have focused on the initial time of symptom appearance in plants attacked by TYLCV. We checked the content of MDA along with the ROS-scavenging enzymes in response to TYLCV infection at early times (within 48 h). The results showed that more MDA accumulated from 24 h in the TYLCV-sensitive cultivar MK (S), while the process was delayed in T2 (R) after virus infection ( Figure 2A). Furthermore, the MDA content in TYLCV-attacked T2 maintained half of that in MK, which suggested that the virus-resistant lines could alleviate TYLCVinduced lipid peroxidation. All of the ROS-scavenging enzymes had higher activity in T2 than those in MK with the virus treatment at the same time point. Remarkably, the activity of POD and SOD of T2 at 48 h increased by 60% and 90%, respectively, which was 30% and 150% higher than that in MK ( Figure 2B-D). of POD and SOD of T2 at 48 h increased by 60% and 90%, respectively, which was 30% and 150% higher than that in MK (Figure 2b-d). respectively. The data represent the mean ± SD (n = 3). Data were analyzed by an ANOVA in the SPSS software. * (p < 0.05) and ** (p < 0.01) indicate statistically significant differences.

Comparative transcriptome of tomatoes in response to TYLCV infection
To investigate the gene expression pattern at the initial time of virus inoculation, the leaves of MK and T2, before (0 h) and after (24 and 48 h) TYLCV inoculation, were extracted for RNA sequencing. The raw data are summarized in Table S2 and the gene expression level (estimated by FPKM values) is shown in Figure S1. To assess the data quality, Pearson's correlation coefficient between biological replicates was observed (R 2 > 0.91 for all samples ( Figure S2), which indicated that the biological replicates were reliable. Additionally, PCA was performed ( Figure 3a) and the results revealed that the 18 samples could be assigned to different groups, with PC1 and PC2 capturing most of the variance in the data. In addition, the biological replicates were projected closely in the space, which indicated a good correlation between replicates (Figure 3a).
The genes with the parameter of false discovery rate (FDR) below 0.05 and absolute fold change ≥ 2 were considered DEGs. To determine the gene expression changes resulting from TYLCV infection, the DEGs between the control samples (CK; 0 h) and the infected samples (24 and 48 h) were identified according to the criterion. A total of 938, 1010, 65, and 477 DEGs were found after 24/48 h of TYLCV infection in comparison to CK in MK and T2, respectively. Of these DEGs, 688, 695, 37, and 269 genes were up-regulated and 250, 315, 28, and 208 genes were down-regulated, respectively (Figure 3b and Table  S3). Interestingly, the up-regulated DEGs occupied 70%, which suggested that a batch of TYLCV-responsive genes were induced in the initial stage. Moreover, there were 20 DEGs common to two infection stages in MK and T2, compared to the uninfected at 0 h (Figure

Comparative Transcriptome of Tomatoes in Response to TYLCV Infection
To investigate the gene expression pattern at the initial time of virus inoculation, the leaves of MK and T2, before (0 h) and after (24 and 48 h) TYLCV inoculation, were extracted for RNA sequencing. The raw data are summarized in Table S2 and the gene expression level (estimated by FPKM values) is shown in Figure S1. To assess the data quality, Pearson's correlation coefficient between biological replicates was observed (R 2 > 0.91 for all samples ( Figure S2), which indicated that the biological replicates were reliable. Additionally, PCA was performed ( Figure 3A) and the results revealed that the 18 samples could be assigned to different groups, with PC1 and PC2 capturing most of the variance in the data. In addition, the biological replicates were projected closely in the space, which indicated a good correlation between replicates ( Figure 3A). The genes with the parameter of false discovery rate (FDR) below 0.05 and absolute fold change ≥ 2 were considered DEGs. To determine the gene expression changes resulting from TYLCV infection, the DEGs between the control samples (CK; 0 h) and the infected samples (24 and 48 h) were identified according to the criterion. A total of 938, 1010, 65, and 477 DEGs were found after 24/48 h of TYLCV infection in comparison to CK in MK and T2, respectively. Of these DEGs, 688, 695, 37, and 269 genes were up-regulated and 250, 315, 28, and 208 genes were down-regulated, respectively ( Figure 3B and Table S3). Interestingly, the up-regulated DEGs occupied 70%, which suggested that a batch of TYLCV-responsive genes were induced in the initial stage. Moreover, there were 20 DEGs common to two infection stages in MK and T2, compared to the uninfected at 0 h ( Figure 3C), which indicated that these 20 DEGs might be the hub genes in response to TYLCV infection.

Functional Annotation of DEGs
To gain an insight into the functions of the DEGs, GO and KEGG pathway enrichment analysis was applied. The graphical results revealed that a large number of DEGs were enriched in GO terms, such as the metabolic process, cellular process, response to stimulus, biological regulation, and signaling ( Figure 4). Notably, the number of DEGs involved in all the categories was much greater in the S line MK than that in the R line T2, which indicated that TYLCV might severely damage the biological processes in the S line. In addition, a few DEGs were found to participate in the immune system process, cell killing, and the rhythmic process. As for the molecular function class, several DEGs were classified into metabolic pathways, biosynthesis of secondary metabolites, plant-pathogen interaction, and MAPK signaling pathway. On top of this, a set of DEGs were involved in the phytohormone signaling pathway, e.g., zeatin biosynthesis (ZT), sesquiterpenoid and triterpenoid biosynthesis (SA), linolenic acid metabolism (JA), carotenoid biosynthesis (ABA), tryptophan metabolism (IAA), and plant hormone signal transduction. Hence, we further explored the DEGs involved in the metabolism of phytohormones and their signal pathways systematically.  Table S4.

Response of Plant Hormones to TYLCV in Tomatoes
Among the genes involved in the SA biosynthesis and signal pathway, SA carboxyl methyltransferase (SAMT, Solyc09g091550.3) and a set of SA-binding protein 2 homologs (Solyc02g065280. 3 Table S5). In the IAA-mediated signaling pathway, there were nine downregulated genes and 10 up-regulated genes after TYLCV inoculation ( Figure 5 and Table S5). In the case of the JA signaling pathway, six and two genes were down-regulated and up-regulated, respectively. The genes involved in the ABA signaling pathway are detailed in Table S5. Taken together, a KEGG map of all of the plant hormone signal transductions is displayed, in which the DEGs induced by TYLCV infection are highlighted ( Figure S3). These results indicated that TYLCV universally reprogramed the expression of the genes related to the biosynthesis and signaling transduction of plant hormones.  Table S5). In the IAA-mediated signaling pathway, there were nine down-regulated genes and 10 up-regulated genes after TYLCV inoculation ( Figure 5 and Table S5). In the case of the JA signaling pathway, six and two genes were down-regulated and up-regulated, respectively. The genes involved in the ABA signaling pathway are detailed in Table S5. Taken together, a KEGG map of all of the plant hormone signal transductions is displayed, in which the DEGs induced by TY-LCV infection are highlighted ( Figure S3). These results indicated that TYLCV universally reprogramed the expression of the genes related to the biosynthesis and signaling transduction of plant hormones. To further explore the roles of plant hormones in the crosstalk between TYLCV and tomatoes, the contents of several plant hormones were measured after 24/48 h of TYLCV infection ( Figure 6). The contents of zeatin and indole-3-acetic acid (IAA) were sharply increased after 48 h of TYLCV inoculation ( Figure 6). SA is essential to plants for acquiring resistance to many pathogens, whose content was up-regulated in the R line T2 after 48 h of TYLCV infection. However, SA was reduced in the S line MK ( Figure 6). Interestingly, in the first 48 h of infection of TYLCV, the contents of JA and ABA were also obviously affected. These results indicate that tomatoes probably regulate the biosynthesis of SA to To further explore the roles of plant hormones in the crosstalk between TYLCV and tomatoes, the contents of several plant hormones were measured after 24/48 h of TYLCV infection ( Figure 6). The contents of zeatin and indole-3-acetic acid (IAA) in TYLCVresistant cultivar (T2) were sharply increased after 48 h of TYLCV inoculation, while it is opposite in susceptible line (MK) (Figure 6). SA is essential to plants for acquiring resistance to many pathogens, whose content was up-regulated in the R line T2 after 48 h of TYLCV infection. However, SA was reduced in the S line MK ( Figure 6). Interestingly, in the first 48 h of infection of TYLCV, the contents of JA and ABA were also obviously affected. These results indicate that tomato plants of T2 probably regulate the biosynthesis of SA to defend against pathogen attack, making sure that plants produce more ZT and IAA for growth. R PEER REVIEW 9 of 14 defend against pathogen attack, making sure that plants produce more ZT and IAA for growth.

Quantitative RT-PCR validation of the RNA-seq expression
To validate the reliability of the RNA-seq data, we randomly selected nine DEGs to perform qRT-PCR. These genes included jasmonate ZIM-domain protein 1 (Solyc12g009220.2), LRR receptor-like serine/threonine-protein kinase GSO1 isoform X1 (Solyc08g016270.3), and transcription factor bHLH120-like (Solyc10g008270.3) (see Table  S1 for more details). The results showed a good correlation between the FPKM values and the qRT-PCR results (Figure 7), which suggested that the gene expression analysis by RNA-seq was reliable.

Quantitative RT-PCR Validation of the RNA-Seq Expression
To validate the reliability of the RNA-seq data, we randomly selected nine DEGs to perform qRT-PCR. These genes included jasmonate ZIM-domain protein 1 (Solyc12g009220.2), LRR receptor-like serine/threonine-protein kinase GSO1 isoform X1 (Solyc08g016270.3), and transcription factor bHLH120-like (Solyc10g008270.3) (see Table S1 for more details). The results showed a good correlation between the FPKM values and the qRT-PCR results (Figure 7), which suggested that the gene expression analysis by RNA-seq was reliable. Figure 7. Verification of nine randomly selected DEGs by qRT-PCR. Comparison of RNA-seq data (blue bar) with qRT-PCR data (red line). The normalized expression level (FPKM) of RNA-seq is indicated on the Y-axis to the left, and the relative qRT-PCR expression level is shown on the Y-axis to the right. α-Tubulin (Solyc04g077020.2) was used as the internal control. Both methods show similar gene expression trends. Three biological replicates were performed. Data were analyzed by an ANOVA in the SPSS software, and a, b, c, and d indicate statistically significant differences (p < 0.05) for the designated time point.

Discussion
Over the past decade, a great deal has been learned about viral vector-host interactions. The dramatic rise in the populations of Bemisia tabaci worldwide and the subsequent virus epidemics have led to increased research aimed at controlling the economic losses resulting from TYLCV. To block transmissible whitefly-borne geminiviruses, the host selection behavior of the vector has been studied [5,48]. Moreover, from an ecological perspective, the application of chemical insecticides and the improvement in cultural methods could protect tomatoes from attack by whitefly carrying the TYLCV [7][8][9]. To disturb the virus transferring to plants, the features of crops should be another key point to explore [5,8]. Generally, visual cues, such as plant color, and olfactory cues, such as plant volatiles, are the first barriers to prevent crops from whitefly [5]. Furthermore, the inbred lines or cultivars harboring the TYLCV-resistant locus or genes have been utilized for breeding and production. Among the germplasm, six genes have been identified, including Ty-1 to Ty-6 [11][12][13][14][15][16][17][18][19]. Ty-1 and its allelic gene Ty-3, which have been widely used in tomato breeding for a long time and have been studied several times concerning their siRNA [20], and now some new variants of TYLCV have emerged to break through resistance to this type of gene [23]. The understanding of Ty-4, Ty-5, and Ty-6 is still limited. Moreover, Ty-2, which was identified as TYNBS1, plays an increasingly important role in the breeding of resistant varieties, but its specific mechanism of anti-TYLCV remains unclear.
Data from plant physiology and molecular analyses demonstrated that TYLCV can repress tomato growth from the initial time of TYLCV inoculation through modulating

Discussion
Over the past decade, a great deal has been learned about viral vector-host interactions. The dramatic rise in the populations of Bemisia tabaci worldwide and the subsequent virus epidemics have led to increased research aimed at controlling the economic losses resulting from TYLCV. To block transmissible whitefly-borne geminiviruses, the host selection behavior of the vector has been studied [5,48]. Moreover, from an ecological perspective, the application of chemical insecticides and the improvement in cultural methods could protect tomatoes from attack by whitefly carrying the TYLCV [7][8][9]. To disturb the virus transferring to plants, the features of crops should be another key point to explore [5,8]. Generally, visual cues, such as plant color, and olfactory cues, such as plant volatiles, are the first barriers to prevent crops from whitefly [5]. Furthermore, the inbred lines or cultivars harboring the TYLCV-resistant locus or genes have been utilized for breeding and production. Among the germplasm, six genes have been identified, including Ty-1 to Ty-6 [11][12][13][14][15][16][17][18][19]. Ty-1 and its allelic gene Ty-3, which have been widely used in tomato breeding for a long time and have been studied several times concerning their siRNA [20], and now some new variants of TYLCV have emerged to break through resistance to this type of gene [23]. The understanding of Ty-4, ty-5, and Ty-6 is still limited. Moreover, Ty-2, which was identified as TYNBS1, plays an increasingly important role in the breeding of resistant varieties, but its specific mechanism of anti-TYLCV remains unclear.
Data from plant physiology and molecular analyses demonstrated that TYLCV can repress tomato growth from the initial time of TYLCV inoculation through modulating the ROS balance and plant hormone biosynthesis and signaling. As reported, TYLCV causes great damage to susceptible cultivars [1], but the process could be delayed when tomatoes possess the Ty-2 gene (Figure 1). Thus, Ty-2 has been widely introduced to tomato cultivars to combat TYLCV. Through mapping-based cloning, Ty-2 was identified as the NB-LRR family gene TYNBS1, providing the genetic and molecular basis for further studies [24]. Under biotic and abiotic stresses, ROS, superoxide, hydrogen peroxide, and nitric oxide were over accumulated in plants [49]. This process was considered to be harmful by damaging the protein structure, enhancing membrane lipid peroxidation, and producing MDA, leading to cell cytotoxicity [50]. To overcome ROS toxicity, antioxidant enzymes, such as SOD, POD, CAT, APX, and PPO, can clear ROS [51]. As for plant-pathogen defense, ROS may also orchestrate a plant's hypersensitive disease resistance response [51]. Combined with proteomics and metabolomics analysis with TYLCV R/S lines, the expression level of reactive oxygen species (ROS) compounds and antioxidative proteins remains high in tomatoes invaded by the virus [31]. In our study, the MDA content combined with ROS-eliminating enzymes suggested that Ty-2 could maintain the normal physiological condition for more time. Given that Ty-1 suppresses virus infection in tomatoes via the induction of RNA silencing [21], recent studies have shown that siRNAs are also involved in ROS signals to regulate plant immunity [52]. Hence, Ty-2 may have crosstalk with Ty-1 in the regulation of the host defense pathway. Coordinating the two genes could be a potential strategy for tomato TYLCV-resistant breeding.
Plant hormones play a crucial role in crop growth and virus defense. SA involves the production of ROS and the accumulation of pathogenesis-related protein, leading to characteristic hypersensitive responses (HRs), which restrict viruses locally [53]. In SAsignaling-active mutants such as cpr1 and cim10 [54] or the application of exogenous SA [33], the accumulation of TYLCV can be reduced, which can enhance the virus resistance. Our results showed that both SA biosynthesis and signal genes were dramatically induced by TYLCV at the time of infection ( Figure 5A), which illustrates that SA acts as the principal defender. JA is another key phytohormone that tackles viruses and vectors. A previous study showed that the C2 protein from TYLCV decreases the sensitivity toward JA and suppresses JA-mediated defenses [55]. Interestingly, our results argue that JA-related genes can still be up-regulated after TYLCV attack, presenting the drastic arms race between plants and the pathogen. Since plant resistance is usually promoted at the cost of development, growth hormones are also involved in plant immunity. When attacked by geminiviruses, brassinosteroid can induce an ROS burst to resist the virus in an SA-independent manner [56]. Similarly, cytokinin plays a role in ROS induction and synergistic association with SA to participate in the plant immunity [57]. Apart from these essential antivirus hormones, others-including zeatin and IAA-might have also been involved in TYLCV resistance in our results ( Figure 5).
RNA-seq is a powerful tool for identifying the key genes and pathways involved in particular biological processes [58]. Data from RNA-seq with tomatoes infected with TYLCV within 48 h showed that the virus has a big impact on the catalytic activity and metabolic process (Figure 4). This is consistent with the results that both biosynthetic genes and the contents of phytohormones, including SA, IAA, and zeatin, were up-regulated ( Figures 5 and 6). Moreover, the increase in the activity of ROS-scavenging enzymes revealed that they should be downstream of phytohormones to maintain the ROS balance participating in TYLCV resistance as per a previous report [53,57]. Meanwhile, JA, the major defense phytohormone, which could enhance the production of defense compounds against insect herbivores and pathogens [34], was also induced in the initial stage of TYLCV inoculation (24 h) ( Figure 6). This indicates that the contents of phytoalexins such as terpene regulated by JA [35] may be affected by the virus just after infection. For further studies, the metabolome-especially phytoalexin-should be detected to explain TYLCV resistance.

Conclusions
TYLCV is a common pathogen causing damage to tomato production. We used the susceptible line MK and the resistant line T2 to perform TYLCV inoculation. The infection of the virus promoted the accumulation of MDA, resulting in lipid peroxidation stress at the initial time of treatment. The activity of the ROS-scavenging enzymes was higher in T2 (R), suggesting that the ROS balance was important in the virus resistance. RNA-seq was carried out to detect the antivirus-related genes or pathways. Combined with qRT-PCR and the measurement of the phytohormones, both the biosynthesis and signal of SA, IAA, and zeatin were involved in the tomato defense against TYLCV. These insights will likely facilitate the genetic modification of tomatoes for desirable resistance traits and possibly many other species of plants infected by TYLCV.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/horticulturae8020143/s1, Figure S1: Distribution of gene expression (FPKM) level. Figure S2: Pearson correlation analysis of the FPKM values of leaf samples in tomato under TYLCV infection. Figure S3: The Kyoto encyclopedia of genes and genomes (KEGG) map of the signal transduction of different plant hormones. Table S1: Genes and primers for qRT-PCR. Table S2: Summary of aligned and mapped reads. Table S3: Differentially expressed genes in Moneymaker (MK) and CLN2777A (T2). Table S4: GO and KEGG pathway annotation of DEGs. Table S5: Expressed genes for phytohormones biosynthesis and signaling pathway.

Data Availability Statement:
The transcriptome data used in this study have been submitted to the NCBI database under the accession number PRJNA787385. These data have not been released due to unfinished research content but could be available upon request.