Identification of Genes Differentially Expressed between Resistant and Susceptible Tomato Lines during Time-Course Interactions with Xanthomonas perforans Race T3

Bacterial spot caused by several Xanthomonas sp. is one of the most devastating diseases in tomato (Solanum lycopersicum L.). The genetics of hypersensitive resistance to X. perforans race T3 has been intensively investigated and regulatory genes during the infection of race T3 have been identified through transcriptional profiling. However, no work on isolating regulatory genes for field resistance has been reported. In this study, cDNA-amplified fragment length polymorphism technique was used to identify differentially expressed transcripts between resistant tomato accession PI 114490 and susceptible variety OH 88119 at 3, 4 and 5 days post-inoculation of the pathogen. Using 256 selective primer combinations, a total of 79 differentially expressed transcript-derived fragments (TDFs) representing 71 genes were obtained. Of which, 60 were up-regulated and 4 were down-regulated in both tomato lines, 4 were uniquely up-regulated and 2 were uniquely down-regulated in PI 114490, and 1 was specifically up-regulated in OH 88119. The expression patterns of 19 representative TDFs were further confirmed by semi-quantitative and/or quantitative real time RT-PCR. These results suggested that the two tomato lines activated partly similar defensive mechanism in response to race T3 infection. The data obtained here will provide some fundamental information for elucidating the molecular mechanism of response to race T3 infection in tomato plants with field resistance.


Introduction
Bacterial spot severely affects tomato (Solanum lycopersicum L.) fruit production and quality in both open field and protected area [1]. The causal agent is a complex of at least four species of Xanthomonas (X. euvescatoria, X. vescatoria, X. perforans, and X. gardneri) with five races designated T1 through T5 based on their virulence on tomato differential genotypes [2][3][4]. Due to the existence of multiple pathogen species and races, marginal efficacy of commonly applied chemicals, development of resistance to these chemicals in bacterial populations, and a lack of available disease resistance traits in commercial cultivars, control of the disease has not been effective once epidemics start [5,6]. Exploiting host resistance gene(s) combined with important defense response genes for developing cultivars with durable resistance is considered as an effective approach to manage the disease.
The resistance to race T3 of bacterial spot in tomato can be either qualitatively or quantitatively inherited. Two S. pimpinellifolium accessions PI 128216 and LA 1589 as well as one unimproved tomato breeding line Hawaii 7981 show both hypersensitive response (HR) and field resistance to race T3.
The HR is conditioned by single dominant genes of Xv3 in Hawaii 7981 [7], Rx4 in PI 128216 [8], and Rx LA1589 in LA 1589 [9]. These three genes have been mapped to the same region on chromosome 11 and might be the same gene or allelic genes [1]. Field resistance to race T3 in these lines shows partial or incomplete dominance, requiring interactions between Xv3 and some modifiers in Hawaii 7981 [10], or depending on gene dosage and genetic backgrounds in PI 128216 [11].
The S. lycopersicum var. cerasiforme accession PI 114490 has been considered as a durable source for resistance to bacterial spot due to its high level of field resistance to four races T1-T4 [5]. Classical genetic analysis of resistance to race T2 using F 2 and inbred backcross populations derived from PI 114490 suggests that genetic control is conferred by a minimum of two loci [12], and resistance to races T3 and T4 is conditioned by at least four quantitative trait loci (QTL) [13,14]. A couple of QTLs are common for resistance to races T2, T3, and T4. Although there is no or few lesions on PI 114490, the bacterial population of race T3 in its leaves is not significantly different from that in susceptible variety OH 88119 [15]. In addition, PI 114490 does not show HR to the pathogen of bacterial spot [1]. Therefore, the mechanism of resistance in PI 114490 seems complicated and different from that in those lines with HR to the pathogen.
Despite extensive efforts have been made on understanding the genetic basis of resistance, little is known about the defense regulation and mechanism underlying inducible response to race T3 of bacterial spot in tomato. To date, a total of 426 genes differentially expressed during the time-course of HR to race T3 in Hawaii 7981 [16] and 1345 genes triggered by recognition of the Xanthomonas type III effector AvrXv3 [17] have been identified. These genes have been predicted to participate in a complex molecular network of regulation including components of defense responses, stress transcriptional regulation factors, signal transduction components, and regulators of primary and secondary metabolisms. Furthermore, several genes involved in the disease resistance process have been isolated through this approach [18,19]. However, these studies are based on hypersensitive response to the pathogen of bacterial spot in tomato plants and may not cover all genes during the process of field resistance.
The complementary DNA-amplified fragment length polymorphism (cDNA-AFLP) is a powerful method to obtain a wide collection of differentially expressed transcript profiles even if rarely expressed during the process of response to abiotic or biotic stresses [20][21][22], enabling discovery of novel genes without any prior knowledge of gene sequences [23]. Moreover, compared with hybridization-based approaches including DNA chips and microarrays, cDNA-AFLP has a relatively low startup cost and enables distinguishing expression patterns of highly homologous gene family members [24]. The cDNA-AFLP technique has been successfully employed to identify differentially expressed genes in various plant-pathogen systems [22,23,[25][26][27]. In this study, we isolate 71 genes differentially expressed in resistant line PI 114490 and susceptible line OH 88119 during time-course of response to bacterial spot race T3 using the cDNA-AFLP technique. The expression patterns of 19 genes are validated using semiquantitative RT-PCR and/or quantitative real time RT-PCR (qRT-PCR). The data obtained here will provide some basic information for understanding the molecular mechanism of response to bacterial spot race T3 infection in tomato plants with field resistance.

Plant materials
Two tomato lines, the S. lycopersicum var. cerasiforme accession PI 114490 with field resistance and a variety OH 88119 susceptible to X. perforans race T3 [14], were used in this study. All seeds were surface sterilized with 4% NaClO for 5 min and germinated in 288 Square Plug Trays (Taizhou Longji Gardening Materials Co., Ltd, Zhejiang, China) in a sterilized mixture of peat and vermiculite (3:1). One-month old seedlings were transplanted into 10 cm (diameter) x 8 cm (height) pots filled with the same sterilized peat : vermiculite mixture and placed in a growth room with 28/25uC day/night cycle and 14 h photoperiod. Water and fertilizer were provided as needed.

Inoculum preparation and inoculation
X. perforans race T3 strain Xv829 was grown on yeast, dextrose, and calcium carbonate (YDC) agar medium [28] at 28uC for 48 to 72 h. Bacterial cells were collected and resuspended to A 600 = 0.2 (approximately 3610 8 colony forming units per ml) using sterile solution containing 10 mM MgSO 4 ?7H 2 O and 0.025% (v/v) Silwet L77. About five to six-week old seedlings were sprayinoculated with bacterial suspensions, whereas control plants (mock-treatment) were sprayed with the same sterile solution. Leaf samples were collected from the infected and mock-treated plants at 3, 4, and 5 days post-inoculation (DPI) for RNA isolation.

RNA isolation and cDNA synthesis
Total RNA was isolated from both infected and mock-treated leaves using Trizol reagent following the manufacturer's instructions (Invitrogen, CA, USA). The RNA integrity was monitored by 1.2% (w/v) agarose gel electrophoresis and the concentration was determined using a Nanodrop 2000 Spectrophotometer (Thermo Fisher Scientific, DE, USA). mRNA was isolated following the protocol of polyA Tract mRNA isolation system kit (Promega, WI, USA), and used as the template for double-stranded cDNA synthesis using M-MLV reverse transcriptase (TaKaRa, Dalian, Liaoning, China). The double-stranded cDNAs were purified using phenol/chloroform extraction and ethanol precipitation, and then stored at -20uC.

cDNA-AFLP analysis
The template for cDNA-AFLP was prepared using the method described in Bachem et al. [29]. Approximately 100 ng cDNA from each sample was first digested with restriction enzyme EcoR I at 37uC for 8 h, and then with Mse I at 65uC for 8 h. The doublestranded adapters were ligated to the digested products using T4 DNA ligase (Promega) at 16uC overnight. After ligation of the adapters, 1 ml diluted ligation mixture (1:10) was used as the template for pre-amplification in a 20 ml reaction mixture containing 10 ml Taq10 26Master Mix (Ausable Biotechnology, Beijing, China), 8 ml ddH 2 O, and 0.5 ml (10 mM) pre-amp primers. The PCR thermal cycling parameters were as follows: 94uC for 3 min, 26 cycles of 94uC for 30 s, 56uC for 1 min, and 72uC for 1min, and a final extension at 72uC for 5 min. The preamplified products were diluted (1:50) with ddH 2 O for selective amplification using 256 AFLP selective primer combinations under the following conditions: 94uC for 3 min, 12 cycles of 94uC for 30 s, 65uC for 30 s (decrease by 0.7 each cycle), and 72uC for 1 min, then by 23 cycles at 94uC for 30 s, 56uC for 30 s, and 72uC for 1min, and an extra extension at 72uC for 5 min for the last reaction. Following the selective amplification, 10 ml of each PCR product was mixed with 2 ml of formamide stop/loading buffer (95% formamide, 20 mM EDTA, pH 8.0, and bromophenol blue), denatured at 98uC for 10 min, and then chilled immediately on ice. The products were separated on 7% denaturing polyacrylamide gel at a constant power (95W) for 4.5 h using 20 cm645cm gel apparatus (Model: JY-CX1B, Beijing Junyi-Dongfang Electrophoresis Equipments Co., Ltd, Beijing, China ) and visualized using the silver-staining approach [30].
Isolation and sequence analysis of differentially expressed transcript-derived fragments (TDFs) AFLP fragments with persistent up-or down-regulated at 3, 4 and 5 DPI were excised from polyacrylamide gels with a clean razor blade, directly rinsed with 200 ml TE buffer (10 mM Tris and 1 mM EDTA, pH 8.0), and then incubated in 100 ml TE for 10 min using a boiling water bath. The fragments were reamplified using the same conditions of selective amplification protocol described above with the corresponding primers. The resulting PCR products were separated on 1.5% agarose gel and purified using Agarose Gel purification Kit (Biomed Co., Beijing, China).
The purified PCR products were ligated into the pMD19-T vector (Takara), transformed into E. coli DH5a competent cells, and plated on Luria-Bertani (LB) agar medium containing 50 mg/ml ampicillin. The recombinant clones were verified by PCR, and at least two positive clones were sequenced at Sunbiotech Company (Beijing, China). Each sequence was searched against the International Tomato Annotation Group (ITAG) Release 2.3 Predicted CDS (SL2.40) in the Sol Genomics Network (SGN) (http://solgenomics.net) database using BLASTX program to obtain the function annotation. Gene ontology analysis was performed using GoPipe [31] through BLASTP against Swiss-Prot and TrEMBL database. The WEGO software (http://wego. genomics.org.cn) was used for gene ontology (GO) classification and plotting [32].

RT-PCR and qRT-PCR analysis
To confirm the results of cDNA-AFLP analysis, 19 strongly differentially expressed TDFs were chosen for semi-quantitative RT-PCR and/or qRT-PCR analysis. Since the sequence lengths of TDFs were short, the sequences of corresponding CDS from SGN were used to design gene-specific primers. The specific primers for TDF792 were designed using the sequence of 39UTR region obtained using the SMARTer RACE cDNA Amplification Kit (Clontech Laboratories Inc., CA, USA). The specificity of each primer pair was confirmed by analyzing PCR products on agarose gel and melting curve during real-time PCR.
RT-PCR was performed using the method described in Pei et al. [8], and the products were separated on 1.2% agarose gel, stained with ethidium bromide, and photographed using a GIAS-4400 Gel Documentation System (Beony Science and Technology Co., Beijing, China). qRT-PCR was conducted in a 20 ml reaction volume consisting of 10 ml SYBR-Green supermix (TaKaRa), 0.2 ml ROXII, 0.2 ml each of forward and reverse primers (10 mM), and 9.4 ml cDNA templates diluted (1:40) from the first strand cDNA solution prepared for cDNA-AFLP. Reactions were heated at 95uC for 3 min followed by 40 cycles of 95uC for 15 s, 57uC for 30 s, and 72uC for 15 s using the ABI 7500 machine (Applied Biosystems, CA, USA). All samples were performed in triplicate with the control sample as the internal reference. The relative expression of the target genes was calculated by normalizing Ct values in each sample against the EF1-a using comparative C-T method [33].

Results
Responses to X. perforans race T3 in OH 88119 and PI 114490 The development of disease on plants was monitored each day after inoculation. No symptom was observed on both PI 114490 and OH 88119 plants prior to 4 DPI. Small lesions with yellow haloes developed on the leaves of OH 88119 on 5 DPI, and subsequently expanded to form large necrotic areas throughout the inoculated leaves in the late stage of the disease (Fig. 1A). In contrast, PI 114490 always showed few and small lesions restricted to portions of the leaves (Fig. 1B). This was consistent with our previous observations [15].
Based on the progress of disease development, the time-course transcriptional profiling of PI 114490 and OH 88119 were conducted to identify differentially expressed genes during the process of race T3 infection. Infected and mock-treated leaves were harvested at 3, 4 and 5 DPI for RNA isolation and cDNA-AFLP analysis.

Isolation and analysis of differentially expressed TDFs
A total of about 6400 TDFs with sizes of 100 bp to 500 bp were amplified using 256 AFLP selective primer combinations. Various types of differentially expressed TDFs between PI 114490 and OH 88119 at three time-points (3, 4, and 5 DPI) were observed, which mainly included commonly and persistently up-or down-regulated in both lines (Fig. 2a,b), specifically up-or down-regulated in either PI 114490 or OH 88119 (Fig. 2c,d), constitutively expressed only in PI 114490 or OH 88119 ( Fig. 2e-g). Considering the occurrence of fragments unique to PI 114490 or OH 88119 might be due to the mutation of restriction sites for EcoR I or Mse I, only TDFs commonly presented in both lines with persistently up-or downregulated at three time-points were excised from the polyacrylamide gel. Therefore, 122 TDFs were recovered from the gels, and 108 were successfully sequenced. Due to inconsistency between sequences derived from two clones, 29 TDFs were discarded and the remaining 79 TDFs were retained for further analysis (Table  1). Among the 79 TDFs, 68 (86.1%) showed up-regulated and four showed down-regulated expression in both PI 114490 and OH 88119, four were specifically up-regulated and two were specifically down-regulated in PI 114490, and one displayed specifically up-regulated in OH 88119 (Table 1).
Based on the ITAG annotation, the 79 TDFs represented 71 genes (Table 1). Four TDFs (TDF742, TDF747, TDF772, and TDF785) did not show sequence similarity to any gene in the

Gene ontology enrichment analysis of differentially expressed genes
Based on the functional annotation using GoPipe [31], 61 of the 71 differentially expression genes had GO terms. Gene ontology enrichment analyses assigned the 61 genes to 22 functional groups in three main categories including cellular component, molecular function, and biological process (Fig. 3). The dominant subcategories for each main category were cell, binding, and metabolic process, respectively. In molecular function category, genes involved in catalytic activity also well represented, while the cellular component category mainly consisted of cell part. It was noteworthy that 11 differentially expressed genes involved in ''response to stimulus'' category (GO:0050896) within biological process. With the exception of TDF352 that was specifically upregulated in OH 88119, the remaining 10 TDFs (30, 270, 702, 260, 318, 791, 352, 600, 87, and 255) were up-regulated in both PI 114490 and OH 88119. Functional annotation analysis revealed that some differentially expressed TDFs that were not classified into the 'response to stimulus' subcategory could also be induced during other biotic and abiotic stresses. The TDF781 without GO term showed approximately 91.5% of amino acid identity with PvPR4, one of the smallest wound-induced proteins that mRNA expression level increased 10-fold upon wounding over a period of 24 h in Phaseolus vulgaris [34]. Another TDF108 (So-lyc10g084320.1.1, Subtilisin-like protease) not assigned to the 'response to stimulus' subcategory was significantly up-regulated in tomato plants with response to Clavibacter michiganesis subsp. michiganesis at 4 and 8 DPI [35].

Validation of expression patterns for selected genes
To select a reliable reference gene for RT-PCR and qRT-PCR analysis, the expression of commonly used reference genes (EF1-a, 18S rRNA and Ubi) was analyzed with the same sample used for cDNA-AFLP analysis. When cDNA concentration of all the samples was adjusted to 0.4 mg/20 ml in qRT-PCR reaction system, the sample's Ct values of EF1-a and 18S rRNA were at    Combining the current results and internal reference genes used in literatures [36,37], the EF1-a was used to normalize the amount of cDNA of all samples. A total of 19 strongly modulated TDFs (Table 2) were used to validate the cDNA-AFLP results. Based on the functional annotation in literatures, majority of these TDFs were associated with massive defense response processes including jasmonic acid biosynthesis pathways (TDF711and TDF741), ethylene biosynthesis (TDF87 and TDF721), defense-related enzymes (TDF30,  TDF657, TDF673 TDF702, TDF108, TDF701, TDF712, and  TDF745), and cell death related (TDF792 and TDF781). The  remaining (TDF625, TDF626, TDF666, TDF667, and TDF794) did not match any known defense response genes.
The first strand cDNA of 3 or 5 DPI was used for the RT-PCR validation (Fig. 4). Except that TDF625 was down-regulated, other 18 TDFs showed up-regulation in both PI 114490 and OH 88119 upon race T3 infection, which were consistent with cDNA-AFLP data. qRT-PCR was used to confirm the fold change of expression levels for 10 TDFs at 3, 4, and 5 DPI. It was very clear that TDF625 was down-regulated and the other nine TDFs were upregulated at different time-points of infection in both PI 114490 and OH 88119, differing in their induction timing and the strength of expression patterns (Fig. 5). TDF30, DF702 and TDF721 displayed significantly high fold-induction in both tomato genotypes during three time-points of infection. At 3 DPI, the expressions of eight TDFs (30,108,702,721,711,657,673, and 792) with homologies to pathogenesis or defense-related genes were significantly induced to higher levels in resistant PI 114490 than susceptible OH88119. However, the up-regulated fold changes of TDF721, TDF711 and TDF657 were then reinforced significantly only in OH 88119 at 4 and 5 DPI. The expression level of TDF745 was significantly up-regulated in OH 88119 at 3 DPI, but fold change value was higher in PI 114490 at 4 and 5 DPI. These results were almost consistent with cDNA-AFLP analysis, indicating the reliability of cDNA-AFLP results.

Discussion
The progresses on understanding genetics of host resistance [1] and the emergency of the pathogen genome sequencing [38,39] have urged bacterial spot to be a suitable pathosystem for investigating the mechanism of tomato-Xanthomonas interactions. Although several loci conferring HR to races T1, T3, or T4 of bacterial spot have been mapped [1,5], and a candidate gene Rx4 conferring HR to race T3 in S. pimpinellifolium accession PI 128216 and genes responding to race T3 in tomato line Hawaii 7981 during the process of HR have been identified [8,16,17], none of these HR genes has been cloned. Therefore, the mechanism of resistance to bacterial spot in tomato remains unclear. Particularly, less work has been done on identification and characterization of genes during the disease process in tomato lines with only field resistance. Using the cDNA-AFLP approach, we here identified 79 TDFs differentially expressed during the disease process in resistant line PI 114490 and susceptible tomato line OH 88119. The expression patterns of some TDFs were further verified by RT-PCR and/or qRT-PCR, suggesting that these TDFs might contribute to host response to the pathogen of bacterial spot during the disease process in tomato. This will provide some fundamental information for investigating the mechanisms of resistance to bacterial spot in tomato.
The number (71) of differentially expressed genes identified in tomato lines PI 114490 and OH 88119 during the disease process was lower comparing to a total of 426 genes identified from tomato line Hawaii 7981 [16]. This could be due to several reasons. First, different techniques might contribute to the difference in numbers of genes identified in two studies. Suppression subtractive hybridization and microarray analysis were used in the previous study [16], while cDNA-AFLP technique was used in this study. Without parallel experiments using these techniques in one study, it was hard to tell the difference in terms of efficiency between two techniques. Second, different criteria might cause different numbers of genes identified in two studies. Gibly et al. [16] only used the resistant tomato line Hawaii 7981 to identify differentially expressed genes with RNAs isolated from pooled leaf tissues collected at various time-points after inoculation. Here we used two tomato lines including one resistance line PI 114490 and one susceptible line OH 88119, and only genes persistently up-or down-regulated at 3, 4, and 5 DPI in both lines were selected. This approach excluded a lot of genes with different expression patterns at three time-points in two tomato lines, which might result in loss of some genes involved in the disease process. Third, different responses to race T3 of tomato bacterial spot in two resistance lines and different inoculation method might also contribute to different numbers of differentially expressed genes identified in two studies. Gibly et al. [16] used infiltration inoculation to incite HR in Hawaii 7981 and we used spray inoculation to evaluate field resistance in PI 114490. Plants with HR should develop localized cell death at infection sites in a short time (,24 h) to restrict pathogen growth, while plants with field resistance might take up to several days to response. Fourth, the number (23.4) of TDFs per primer combination was lower comparing to a previous study of 48.8 TDFs per primer combination in the tomato-Cladosporium fulvum pathosystem [25], which could also contribute to the less genes identified in this study. However, the ratio of differentially expressed TDFs (1.2%) in this study was close to 1% in Gabriëls et al. [25]. Therefore, cDNA-AFLP technique could be used to identify genes differently expressed in tomato plants during the infection of bacterial spot pathogen.
The genetics of field resistance to bacterial spot race T3 in tomato line PI 114490 is complicated. At least five QTLs on chromosomes 1, 3, 8, 11, and 12 with small effect (6.5-11.7%) have been reported to date [14,15], suggesting that the high level of field resistance in PI 114490 requires interaction of some of these QTLs [14]. In the present study, 31 of 71 genes were on the five chromosomes having QTLs and 12 genes were at the QTL regions. Particularly, the TDF785 (Solyc11g051200.1.1, Cytochrome P450 like_TBP) locating on the QTL region of chromosome 11 was exclusively up-regulated in PI 114490, indicating that this TDF might be associated with the QTL. In addition, although most TDFs showed up-regulated in both PI 114490 and OH 88119, the expression patterns were different. For example, the up-regulation of TDF721 showed a reduction trend at 3, 4, and 5 DPI in PI 114490, but an increase trend in OH 88119 (Fig. 5). Therefore, the regulatory mechanism could also be complicated.
Three pathways, Microbe/Pathogen-Associated Molecular Patterns recognition, effector recognition, and phytohormone pathways, have been proposed for plant immune response. Each pathway involves in interactions among many genes. Phytohormones including salicylic acid (SA), gaseous ethylene (ET), and jasmonic acid (JA) can regulate local and systemic resistance to invasive pathogens [40]. The plant hormone ET is an important component of defense signaling [41,42]. Genes encoding AP2-like ethylene-responsive transcript factor (AP2/ERF) and 1-aminocyclopropane-1-carboxylate oxidase (ACC) play an important role in response to mechanical damage and symptom development by pathogen attack [35,43,44]. Tomato ethylene-responsive transcript factor (ERF) Pti4, Pti5 and Pti6 directly interact with Pto resistance gene [42] and play a role in activation of pathogenesisrelated genes, resulting in enhanced defense against certain bacterial and fungal pathogens [45][46][47]. Based on the annotation, the TDF87 identified in this study is a member of the ERF protein family. It showed different levels of up-regulation in both resistant tomato line PI 114490 and susceptible line OH 88119 during the X. perforans race T3 infection, suggesting that it might contribute to the defense against the pathogen. In addition to its contributing to resistance in some plant-pathogen interactions, ethylene is also involved in chlorosis symptom development. The ethyleneinsensitive Arabidopsis and soybean plants showed less chlorosis symptoms than wild-type plants [48,49]. The TDF721, a member of ACC family, was up-regulated during the pathogen infection in tomato lines PI 114490 and OH 88119. However, the resistant line PI 114490 showed less disease symptom on leaves than the susceptible line OH 88119, which might due to the differential upregulation of TDF721 in two tomato lines.
Many plant pathogens have evolved virulence strategies to regulate JA signaling pathway in their hosts to facilitate infection and production of disease symptoms [50,51]. We here identified four up-regulated genes (TDF281, TDF711, TDF741, and TDF744) that might contribute to the regulation of JA biosynthesis during tomato-X. perforans interactions. TDF281 and TDF711 were homologous to different Cytochrome P450 (Table 1). They were induced in both PI 114490 and OH 88119 but the expressions were higher in OH 88119 than in PI 114490 at 4 and 5 DPI (Fig. 5). Blast search revealed that both genes contain allene oxide synthase activity, which is the key enzyme in the JA biosynthesis pathway [52,53]. TDF741 and TDF744 were annotated as Jasmonate ZIM domain 2 (JAZ2) and Jasmonate ZIM-domain protein 3 (JAZ3), respectively ( Table 1). The differentially expressed JAZ genes were also observed in bean-Xanthomonas axonopodis pv. phaseoli interaction [54] and Arabidopsis -Pseudomonas syringae interactions [51]. Characterization of mutants encoding truncated JAZ proteins and RNAi line have demonstrated that some JAZ proteins are both negative regulators of JA signaling and disease symptom development [51,55,56]. Therefore, it is very likely that these two TDFs might play important roles in altering JA signaling to impact on symptom development in tomato-X. perforans race T3 interaction.
In summary, response to X. perforans in tomato plants involves a complicated network regulated by a large number of defense/ pathogenesis-related genes. Using the cDNA-AFLP technique, a total of 79 TDFs representing 71 genes with different expression patterns during the pathogen infection were identified in this study. Differential expression of several JA-and ET-related genes during the process of the pathogen infection indicates that coactivation of JA and ET signaling pathways may play key roles in response to X. perforans race T3 in tomato plants.