Introduction
Tomato yellow leaf curl virus (TYLCV) is a single-stranded circular DNA virus of the genus Begomovirus (family Geminiviridae). Since it was first discovered in Israel, TYLCV has spread to many other places, including the Middle East, Mediterranean coast, Africa, and Asia (Lu et al., 2023). TYLCV is one of the most detrimental viruses to tomato cultivation and it causes severe crop yield losses in tropical and subtropical regions. Because TYLCV is transmitted through whiteflies (Bemisia tabaci) (Gafni, 2003), chemical pesticides are often used to control this virus disease (Palumbo et al., 2001). However, the long-term use of pesticides is harmful to the environment and easily induces resistance in B. tabaci (Cahill et al., 1996). Given its advantages such as low cost and environmental friendliness, the use of beneficial microbes as biocontrol strains has attracted increasing interest in the control of TYLCV disease.
Members of Bacillus (Kumar et al., 2016), Pseudomonas (Shen et al., 2014), and Trichoderma (Elsharkawy et al., 2013; Vitti et al., 2015) have been used to control plant diseases by improving plant systemic resistance. There are two main pathways of systemic resistance in plants (Wang et al., 2009). Systemic acquired resistance (SAR) is dependent on salicylic acid and the regulatory protein NPR1 (Mou et al., 2003), whereas induced systemic resistance (ISR) is modulated by hormones such as jasmonic acid and ethylene. It has been shown that plant beneficial microbes can activate disease resistance pathways by regulating the expression of resistance-related genes (Elsharkawy et al., 2012) and the biosynthesis of hormones (Niu et al., 2011), ultimately enhancing plant resistance to the virus.
In addition to systemic resistance, RNA interference (RNAi) is recognized as an important defense mechanism in plants against viral infection (Kanakala and Ghanim, 2016). RNAi can be divided into small interfering RNA (siRNA)-mediated post-transcriptional gene silencing (PTGS) and transcriptional gene silencing (TGS), depending on the site and mechanism of action (Arif et al., 2013). PTGS occurs in the cytoplasm to precisely target and degrade mRNA transcripts of specific genes, which has been widely demonstrated to play an essential role in plant resistance against viral infection (El-Sappah et al., 2021). This process is mainly mediated by 21–24 nt virus-derived siRNA cut by Dicer-like (DCL) nucleases (Zheng et al., 2019). RNA-directed DNA methylation (RdDM)-mediated TGS is another important mechanism of plant defense against geminiviral infection, and it inhibits virus replication by increasing the methylation level of cytosine on the viral genome (Raja et al., 2008). Among them, methylation of CpG sites is crucial for plant defense against virus invasion, and viruses even evade plant defense by reducing the frequency of CpG sites in the genome (de Pradena et al., 2020; Loew et al., 2020).
The occurrence of PTGS and TGS involves the ARGONAUTE (AGO) family proteins, which are effectors of RNAi in plants (Li Z. C. et al., 2019). A growing number of AGO proteins have been identified to participate in the antiviral defense processes in plants. In particular, AGO4 is a major component of TGS that has been widely studied. In the classical RNA polymerase IV-RdDM pathway, the siRNA binds to AGO4 to form an AGO4-siRNA complex. This complex then binds to polymerase IV to take part in methylation of virus genomic sequences that are homologous to the siRNA under the action of domains-rearranged methyltransferase 2 (Matzke et al., 2015). Additionally, AGO3 has been reported to play a vital role in abscisic acid (ABA)-mediated antiviral defense (Alazem et al., 2014). Evidence suggests that both AGO3 and AGO7 play antiviral roles during TCV infection (Zheng et al., 2019). However, it is still unclear whether inoculation with biocontrol strains could modulate AGOs to enhance plant antiviral defense.
A previous study has demonstrated the robust effects of a bacterial strain, Bacillus amyloliquefaciens Ba13, for plant growth promotion and TYLCV disease control in tomato (Guo et al., 2019). How the presence of Ba. amyloliquefaciens Ba13 modulates the RNAi-based antiviral system in the host plant remains unclear. Here we verified the biocontrol effects of Ba. amyloliquefaciens Ba13 against TYLCV and explored the underlying mechanisms from the RNAi perspective. We analyzed the changes in RNAi-related gene expression in TYLCV-infected tomato leaves, as well as the methylation rate of TYLCV genome and the abundance of antiviral siRNAs in plant cells, after bacterial inoculation. The aim of this study was to decipher the possible role of Ba. amyloliquefaciens Ba13 in the antiviral RNAi system against TYLCV. The results could provide new evidence for microbial control of viral plant diseases.
Materials and methods
Experimental materials
The biocontrol strain, Ba. amyloliquefaciens Ba13 (SUB3558699; GenBank: MG846076), was isolated from a stratified old manured loessial soil under a corn-wheat rotation system in Yangling, Shaanxi Province, China. The strain was preserved by the Resource and Environmental Biology Laboratory in the College of Resources & Environment, Northwest A&F University (Yangling, Shaanxi Province, China). Seeds of tomato (Solanum lycopersicum L.) cultivar “Jinpeng 1” were purchased from Xi’an Jinpeng Seedling Co., Ltd. (Xi’an, Shaanxi Province, China). To mimic natural viral infection, the TYLCV pathogen was sourced from viruliferous whiteflies (B. tabaci) collected in the field in Yangling.
Seedling cultivation, bacterial inoculation, and viral infection
The pot experiment was carried out in a greenhouse on campus of Northwest A&F University. Plump tomato seeds were selected and surface sterilized with 75% ethanol for 30 s, followed by at least five washes with sterile distilled water. The sterilized seeds were placed evenly onto a filter paper in sterile Petri dishes (∼50 seeds per dish). To each dish, 2 mL of sterile water was added and the seeds were incubated in the dark at 26°C until germination. The germinated seeds were evenly sown into potted trays (50 holes per tray) containing seedling substrate and cultivated for 21 days. When true leaves formed, tomato seedlings with uniform growth were selected and transplanted into 60 pots with a diameter of 18.3 cm and a height of 12 cm (one seedling per pot). Each pot contained fertile field soil collected from a wheat-maize rotation system in Yangling. The soil was classified as Eum-Orthic Anthrosols (Staff, 2010). The soil was passed through a 2-mm sieve before use. The experimental used a randomized complete block design with three blocks, and each block contained 20 pots. The pots in each block were randomly and equally divided into two groups. One group was inoculated with Ba. amyloliquefaciens Ba13 (treatment group), while the other group was not inoculated (control group).
Ba. amyloliquefaciens Ba13 was cultured in beef extract–peptone broth at 37°C for 48 h with a shaking frequency of 180 rpm. The bacterial culture was diluted to 1 × 108 colony-forming units (CFU) mL–1 with sterile phosphate-buffered saline. For the treatment group, 100 mL of diluted cell culture per pot was applied to the seedlings 7 days post-transplantation. The control group received 100 mL of sterile phosphate-buffered saline per pot. Natural viral infection was performed 14 days after bacterial inoculation. To do this, both non-inoculated controls and Ba. amyloliquefaciens Ba13-inoculated plants were placed in a greenhouse with whiteflies. Each block was covered with a 150-mesh high-quality nylon net of 1.2 m × 1.2 m × 1.5 m (length × width × height), and 1000 whiteflies were released under the net for infection over 4 weeks. The average temperature of the greenhouse was 28°C and the average light period per day was 12 h.
Disease severity evaluation, growth measurement, and leaf sampling
The infection rate and disease severity of tomato plants in each group were determined at 14 and 28 days post-viral infection (dpvi). The disease severity was scored as follows (Sun et al., 2016): 0 = healthy plant without with curling yellow leaves; 1 = mild dwarfing (∼4/5 height of healthy plant) with curling yellow leaves (<20% on the top); 2 = dwarfing (∼2/3 height of healthy plant) with leaf curling and yellowing (<40% on the top); 3 = obvious dwarfing (∼1/2 to 2/3 height of healthy plant) with typical curling yellow leaves (60%); and 4 = severe dwarfing (<1/2 height of healthy plant) with typical curling yellow leaves (60% to all), or early withering. The disease index and control rate were calculated as follows.
Disease
index
(
%
)
=
∑
(
Disease
severity
×
corresponding
plant
number
)
4
×
Total
number
of
plants
×
100
Control
rate
(
%
)
=
Disease
index
of
control
group
-
disease
index
of
treatment
group
Disease
index
of
control
group
×
100
Plant height was measured 28 dpvi and leaf samples were taken for subsequent analyses. Three biological samples (repetitions) were collected from each group, and each biological repetition was a mixture of three leaf samples collected from the same part of three random plants in each block. The sample consisting of the first and second leaves from the apical growing point was used for virus quantification; the third leaf was used for whole-transcriptome RNA sequencing (RNA-seq), bisulfite sequencing, and small RNA-seq, as well as for quantitative validation of functional genes; and the fourth leaf was used for the analysis of plant hormone contents. The leaf samples destined for virus quantification and sequencing analyses were rinsed with distilled water and blotted dry, then immediately wrapped with tin foil and placed in liquid nitrogen. The samples for hormone analysis were weighed and then placed into liquid nitrogen. All samples were brought to the laboratory where they were kept at −80°C until used.
Virus genotyping and quantification
Total DNA was extracted from tomato plants using a modified hexadecyltrimethylammonium bromide (CTAB) method. PCR amplification was performed as described in the Supplementary Note 1, with previously reported primers for the full-length sequence of the complete TYLCV genome (Supplementary Table 1). The amplified products were sequenced to obtain the complete genome sequence. Then we determined the virus type by comparing the obtained sequence with available virus genotypes in the NCBI database. Gene sequence similarity analysis was performed using the MegAlign program in DNAstar (DNASTAR, Inc., Madison, WI, USA). The MEGA-X software (Kumar et al., 2018) and MUSCLE algorithm were used to align the nucleic acid sequences of multiple sequences, while the neighbor joining algorithm was applied to build phylogenetic trees for taxonomic analysis of the virus. There were 23 different TYLCV genotypes in the evolutionary tree (including the virus genotypes identified in this study). Their accession numbers in the NCBI database, origins, and the similarities to the genotype of the virus used in this study are listed in Supplementary Table 2. The search of virus open reading frame (ORF) was performed using the NCBI Open Reading Frame Finder.1
Quantification of TYLCV in tomato leaves was performed using the method of Sade et al. (2014). Total DNA was extracted from leaf samples using the M5 Plant Genomic DNA Kit (Mei5 Biotech Co., Ltd., Beijing, China). The concentration and purity of extracted DNA were measured using spectrophotometry with an ultra microplate (Epoch, BioTek, Winooski, VT, USA). The DNA integrity was detected using agarose gel electrophoresis. The DNA samples were subjected to real-time fluorescent quantitative PCR on an iQ5 PCR system (Bio-Rad, Hercules, CA, USA) as described in the Supplementary Note 1. The PCR amplification was performed using TYLCV-V1 as the primer (Supplementary Table 1) and the EF1 gene as the internal reference. There were three technical replicates for each sample. Data collection and analysis were performed using the Bio-Rad iQ5 software (v2.1.97). The relative expression level of the viral gene was calculated using the 2–ΔΔCT method (Schmittgen and Livak, 2008).
Whole-transcriptome RNA-seq and gene functional analysis
Frozen leaf samples were sent to BGI Biotech Co., Ltd. (Shenzhen, China) for transcriptome sequencing using the BGISEQ-500 platform (BGI Biotech). Total RNA was isolated from the leaf samples using Trizol reagent (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. Then mRNA libraries were constructed. Whole-transcriptome RNA-seq data processing was performed as described in Supplementary Note 2. The detection of significantly differentially expressed genes (DEGs) was performed following the method described by Wang et al. (2010). Genes with a | fold-change| > 2 and p-value ≤ 0.001 were defined as DEGs.
The DEGs were subjected to Gene Ontology (GO) and pathway enrichment analysis using the PANTHER classification system2 (Mi et al., 2005) and the Kyoto Encyclopedia of Genes and Genomes (KEGG),3 respectively (Kanehisa et al., 2008). The basic local alignment search tool (BLASTn) was adopted to identify homologous genes in the tomato genome.4 The raw sequences were deposited at the Sequence Read Archive in the National Center for Biotechnology Information (NCBI)5 under BioProject accession number: PRJNA553064.
Quantitative PCR validation
The validation of functional genes was performed by reverse transcription quantitative real-time PCR (RT-qPCR) using six genes involved in gene silencing, plant growth, and disease resistance. The PCR protocols are provided in Supplementary Note 3. There were three technical replicates for each sample, and the primer sequences used are provided in Supplementary Table 1. The relative expression levels of the selected genes and small RNAs were calculated as described for virus quantification.
ABA analysis
The content of ABA in tomato leaves was measured using a previously reported method (Pan et al., 2010) as described in Supplementary Data (Note S4). Briefly, the ABA was extracted from leaf samples with 80% methanol/water (v/v) followed by a mixture of 5% acetic acid/ethyl acetate and water (1:1, v/v). The levels of ABA in the extracts were measured using gas chromatography coupled to triple quadrupole mass spectrometry (Agilent, Santa Clara, CA, USA).
Viral genome methylation analysis by bisulfite sequencing
Gene-specific DNA methylation was assessed by a next generation sequencing-based BSP, according to previously published method (Gao et al., 2014). In brief, BSP primers were designed using the online MethPrimer software and listed in Supplementary Table 1. A total of 1 μg of genomic DNA was converted using the ZYMO EZ DNA Methylation-Gold Kit (Zymo Research, Irvine, CA, USA) and one twentieth of the elution products were used as templates for PCR amplification with 35 cycles using KAPA 2G Robust HotStart PCR Kit (Kapa Biosystems, Wilmington, MA, USA). For each sample, BSP products of multiple genes were pooled equally, 5′-phosphorylated, 3′-dA-tailed and ligated to barcoded adapter using T4 DNA ligase (NEB). Barcoded libraries from all samples were sequenced on Illumina platform. For the bisulfite sequencing reads of each sample, firstly, adapters and low-quality reads were removed using software Trimmomatic-0.36. After removing the adapter sequences and filtering out the low-quality reads, the clean sequencing reads were directly aligned to the target sequences using software Bsmap (v2.73) with the default parameters which combines genome hashing and bitwise masking to achieve fast and accurate bisulfite mapping. Methylation level of C base in TYLCV genome was calculated as follows.
M
e
t
h
y
l
a
t
i
o
n
r
a
t
e
o
f
C
s
i
t
e
(
%
)
=
methylated
sequences
methylated
sequences
+
demethylated
sequences
×
100
Small RNA-seq and identification of virus-derived siRNA
Small RNA-seq was used to detected the changes in the types and abundances of siRNA in tomato leaves, which have lengths mainly concentrated in the range of 19–24 nt (Seo et al., 2013). For the sample preparation and library construction, small RNA-seq was performed using the BGISEQ-500 platform (BGI Biotech). Sample selection and total RNA extraction followed the same methods that were used for the whole-transcriptome RNA-seq. Small RNA libraries were then constructed and small RNA-seq data were processed as described in Supplementary Note 5. The siRNAs with a p-value ≤ 0.05 were considered significantly differentially expressed. The raw sequences were deposited at the Sequence Read Archive in the NCBI (see text footnote 5) under BioProject accession number: PRJNA553309.
Given that siRNAs can bind to viral mRNA for transcriptional and post-transcriptional inhibition of viral replication, we aligned all the differentially expressed (p < 0.05) siRNAs against the genome of the TYLCV-SJ which identified in this study. The alignment of sequences was performed using the SOAP2.22 short-read alignment software package (Li et al., 2009), with no mismatches allowed. For siRNAs that completely matched to the viral genome sequence, their name, sequence, matched site, abundance in the two treatment groups (Ba. amyloliquefaciens Ba13 and control), and corresponding viral genotype were recorded.
Statistical analysis
Data were statistically analyzed using DPS 4.0.6 The comparison of group means was performed using the Student’s t-test, except that the disease index values obtained at two time points were compared using the least-significant difference (LSD) test. A p-value of less than 0.05 was considered significant. The relative expression levels (fold-changes) of genes and small RNAs were plotted using Excel 2010 (Microsoft Corp., Redmond, WA, USA), and the relative quantities (chromatographic peaks) of plant hormones were drawn using Origin 2022b (OriginLab Corp., Northampton, MA, USA). The enrichment of DEGs in the biotic stress pathway was analyzed using Mapman (Thimm et al., 2004).
Discussion
RNAi is an antiviral system which allow plants to directly act on viruses. The application of Ba. amyloliquefaciens Ba13 as a biocontrol strain against TYLCV necessitates a holistic understanding of how this strain modulates the RNAi-based antiviral system in its host plants. In the present study, we found that inoculation with Ba. amyloliquefaciens Ba13 increased the expression of RNAi-related genes (e.g., AGO3, AGO4, AGO5, AGO7) in tomato leaves, with improved levels of viral genome methylation and ABA content. These results indicate that the tomato plants might have initiated RNAi-related pathways after TYLCV invasion, whereas the application of Ba. amyloliquefaciens Ba13 could improve plant resistance by enhancing RNAi.
As an integral part of RNAi, DNA methylation-mediated TGS is considered to be an important pathway for plant defense against infection by geminiviruses (Mei et al., 2020). For DNA viruses, AGO4 is the main effector protein in the antiviral silencing machinery. Based on the whole-transcriptome RNA-seq data, we observed upregulation of AGO4 expression in tomato leaves upon Ba. amyloliquefaciens Ba13 inoculation. Taking into account the role of AGO4–siRNA complex in viral genome methylation, we expected that Ba. amyloliquefaciens Ba13 mediated increased methylation level of the TYLCV genome. This hypothesis was supported by the higher methylation levels at CpG sites in AV1 and AV2 following Ba. amyloliquefaciens Ba13 inoculation. Previous study has found that after infection with the tomato leaf curl New Delhi virus, the methylation level of the viral intergenic spacer in a resistant cultivar was substantially higher than in a sensitive cultivar after infection. Additionally, Brough et al. (1992) introduced the cytosine methylated DNA of tomato golden mosaic virus into tobacco protoplasts, leading to a reduction in viral replication to 1/20–1/5 of the original level. These observations indicate that Ba. amyloliquefaciens Ba13 inoculation increased the methylation level of viral DNA through upregulation of AGO4 expression, which in turn inhibited AV1 and AV2 expression and enhanced plant resistance to TYLCV. It has been shown that the methylation frequency of CpG sites in the plant genome markedly increases during virus invasion. CpG methylation frequency is a key factor in inhibiting geminivirus replication and reducing plant symptoms (Rodriguez-Negrete et al., 2013). Further study is needed to ascertain whether inoculation with Ba. amyloliquefaciens Ba13 can affect the frequency of plant genome methylation, thereby improving tomato resistance to TYLCV.
In addition to AGO4, other AGO family members (e.g., AGO3, AGO5, AGO7) play non-negligible roles in plant antiviral defense. AGO3 has been reported to take part in ABA-mediated antiviral defense processes (Alazem et al., 2014). Promoters of AGO genes contain several ABA responsive elements, and their expression levels are affected by exogenous ABA (Alazem et al., 2017). Therefore, ABA can enhance the expression of AGO3—which is necessary for plants to resist Bamboo mosaic virus (Alazem et al., 2017). In the case of TYLCV infection, Ba. amyloliquefaciens Ba13 inoculation increased AGO3 expression in tomato plants. To clarify whether the upregulation of AGO3 expression is related to ABA, we looked at ABA biosynthesis related genes and ABA content in tomato leaves. We observed a remarkable upregulation in the expression of ABF upon Ba. amyloliquefaciens Ba13 inoculation, along with an increase in the ABA content. These results provide strong evidence that the application of Ba. amyloliquefaciens Ba13 improved ABA biosynthesis and enhanced the expression of AGO3, thereby facilitating plant defense against TYLCV.
With regard to AGO5 and AGO7 expression, we observed a distinct upregulation in tomato leaves upon Ba. amyloliquefaciens Ba13 inoculation. AGO5 can bind to small RNAs derived from viruses or viroids. Potato virus X (PVX) infection induces AGO5 expression, which is essential for Arabidopsis to limit PVX infection (Brosseau and Moffett, 2015). AGO7 is also involved in PTGS and thus plays a role in plant defense against viruses (Garcia-Ruiz et al., 2015). For example, AGO7 has been shown to collaborate with AGO1 in clearing viral RNA with different levels of secondary structure (Qu et al., 2008). AGO7 additionally participates in plant defense against Turnip crinkle virus, and together with AGO3, it becomes the major contributor to virus clearance in the leaves (Zheng et al., 2019). The collective results allow us to conclude that the application of Ba. amyloliquefaciens Ba13 upregulated the expression of AGO5 and AGO7, which synergistically acted to improve the efficiency of clearing viral RNAs in tomato plants.
While the expression levels of multiple AGO genes were upregulated in tomato plants inoculated with Ba. amyloliquefaciens Ba13, plant resistance to TYLCV was enhanced. It has been reported that an increase in AGO gene expression enhances plant resistance to different viruses such as Bamboo mosaic virus, Potato virus X, Turnip crinkle virus (Qu et al., 2008; Brosseau and Moffett, 2015; Alazem et al., 2017). This leads us to posit that Ba. amyloliquefaciens Ba13 is likely to regulate plant defense against other viruses by upregulating AGO expression. We additionally found that Streptomyces pactum Act12 can increase tomato resistance to TYLCV by upregulating AGO gene expression (Li Y. L. et al., 2019). This suggests that other biocontrol bacteria may perform similar functions as observed in Ba. amyloliquefaciens Ba13. More experiments are required to verify whether the effects and mechanisms of Ba. amyloliquefaciens Ba13 are universal for other viruses and biocontrol strains.
RNAi plays a role in disease control via siRNAs, which target viral mRNAs with corresponding sequences at specific sites for degradation. We compared the obtained siRNA sequences against the TYLCV genotype (MN910280.1) and found 10 siRNAs completely matched to our obtained viral sequence. The siRNA matcher site in the ORF were corresponding to AV1, C1, and C2. AV1 and C1 are, respectively, involved in coding virus coat protein and replication-related proteins, which help TYLCV growth and spread. C2 can revert TGS and decrease DNA methylation in plants by interfering with the methyl cycle (Buchmann et al., 2009). The complete matching of these siRNAs to the functional gene regions of TYLCV indicates that the siRNAs may interfere with viral protein synthesis by complementary binding and degrading mRNA of gene transcription. This result illustrates that RNAi-related pathways were initiated in tomato plants in response to TYLCV infection.
We observed Ba. amyloliquefaciens Ba13-mediated downregulation of all 10 siRNAs matched to the TYLCV genotype in tomato leaves, consistent with the pattern observed in a previous study. Rodriguez-Negrete et al. (2009) investigated the recovery of plants (with less susceptible symptoms) after infection with pepper golden mosaic virus. They found a lower total abundance of virus-derived siRNAs but a higher DNA methylation level of the viral genome in the recovered tissue than in the diseased tissue of pepper plants. Additionally, we found that the disease symptoms of Ba. amyloliquefaciens Ba13-inoculated tomato plants were alleviated compared with those of non-inoculated controls. The decrease in siRNA abundances could be related to the decrease in virus quantity in tomato plants treated with Ba. amyloliquefaciens Ba13. Ba. amyloliquefaciens Ba13 inoculation enhanced plant defense against TYLCV, allowing faster recovery from viral infection.
Our results indicated that the expression levels of genes related to SA, ABA, and RNAi pathways were upregulated in tomato leaves after inoculation with Ba. amyloliquefaciens Ba13. In fact, the defense processes of plants against viruses involve the participation and collaboration of multiple mechanisms, such as RNAi, dominant resistance, and SAR (Nicaise, 2014). Thus, inoculation with Ba. amyloliquefaciens Ba13 is likely to induce multifaceted plant response through various defense systems. Although both SAR and RNAi play a role in plant defense against viruses, RNAi has higher specificity for viruses compared with the broad-spectrum resistance mechanism of SAR and can target degradation of viruses. While this study has focused on RNAi modulation by Ba. amyloliquefaciens Ba13, the relationship between RNAi and other defense mechanisms in inoculated plants should be clarified by further studies.
Furthermore, we found that the application of Ba. amyloliquefaciens Ba13 strongly affected the expression of genes related to photosystem and auxin response pathways. A total of 12 DEGs related to the biosynthesis of LHCII chlorophyll a/b binding protein 1 (LHCb1), LHCb3, and LHCb7 in the LHCII pathway (the most abundant light harvester in plants) were upregulated in tomato leaves upon Ba. amyloliquefaciens Ba13 inoculation. Additionally, 13 of the 14 DEGs related to auxin/indoleacetic acid and small auxin-up RNA, as well as two DEGs related to GH3, were upregulated in tomato leaves of inoculated plants. Montasser et al. (2012) found that the content of photosynthesis related chemical components (e.g., chlorophyll a, b) decreased in tomato leaves infected with TYLCV compared with normal leaves. Wu et al. (2012) showed that auxin/indoleacetic acid proteins were significantly expressed in tomato plants during root and stem development. Small auxin-up RNA plays a major role in modulating cell elongation and positively regulates cell expansion by modulating the transport of auxin (Chae et al., 2012). Furthermore, GH-like genes mediate the interaction between auxin and salicylic acid signaling pathways. Overexpression of GH3 can enhance the systemic disease resistance in plants (Zhang et al., 2007). Taken together, these results indicate that the presence of Ba. amyloliquefaciens Ba13 could attenuate the negative effects of viral infection on leaf photosynthetic efficiency, promote plant growth, and as such, enhance plant defense against TYLCV.