Comparative transcriptome analysis of tomato (Solanum lycopersicum) in response to exogenous abscisic acid

Background Abscisic acid (ABA) can regulate the expressions of many stress-responsive genes in plants. However, in defense responses to pathogens, mounting evidence suggests that ABA plays variable roles. Little information exists about genome-wide gene expression in ABA responses in tomato (Solanum lycopersicum L.), a model fruit crop plant. Results Global transcriptome profiles of tomato leaf responses to exogenous ABA were generated using Illumina RNA-sequencing. More than 173 million base pair reads were mapped onto the tomato reference genome and the expression pattern differences between treated and control leaves were assessed. In total, 50,616 transcripts were generated. Among them, 42,583 were functionally annotated in the NCBI non-redundant database and 47,877 in the tomato genome reference. Additionally, 31,107 transcripts were categorized into 57 functional groups based on Gene Ontology terms, and 14,371 were assigned to 310 Kyoto Encyclopedia of Genes and Genomes pathways. In both the ABA treatment and control samples, 39,671 transcripts were available to analyze their expressions, of which 21,712 (54.73%) responded to exogenous ABA. Of these transcripts, 2,787 were significantly differently expressed genes (DEGs). Many known and novel ABA-induced and -repressed genes were found. Exogenous ABA can influence the ABA signaling pathway with PYR/PYL/RCARs-PP2Cs-SnRK2s as the center. Eighteen PYL genes were detected. A large number of genes related to various transcription factors, heat shock proteins, pathogen resistance, and the salicylic acid, jasmonic acid, and ethylene signaling pathways were up-regulated by exogenous ABA. Conclusions The results indicated that ABA has the potential to improve pathogen-resistance and abiotic stress tolerance in tomato. This study presents the global expression analysis of ABA-regulated transcripts in tomato and provides a robust database for investigating the functions of genes induced by ABA. Electronic supplementary material The online version of this article (doi:10.1186/1471-2164-14-841) contains supplementary material, which is available to authorized users.


Background
Tomato (Solanum lycopersicum L.) is an important fruit crop and a model system in plants. The basic chromosome number of tomato is 2n = 24, and wild forms range from diploids to hexaploids [1,2]. The genome of the tomato was recently published; it possesses about 35,000 genes, a rich resource to be studied by scientists [3].
The hormone abscisic acid (ABA) regulates numerous developmental and functional processes, including stomatal aperture and hydraulic conductivity, seed dormancy, a key phase transitions throughout the plant lifecycle [4][5][6]. An important role of ABA in a variety of plants is to increase tolerance to stresses such as drought, salinity, cold, and heat [7][8][9]. ABA was recognized as an important signaling molecule that can trigger signaling and regulating mechanisms to cope with adverse stresses. Under drought conditions, plants synthesize and redistribute ABA, leading to guard cell responses that close stomata, thereby reducing plant water loss [10]. Research in chickpea (Cicer arietinum) has shown that ABA can improve cold tolerance by, in part, improving the water status of leaves and anti-oxidative ability and can enhance heat tolerance via the accumulation of osmoprotectants [11,12]. Rice (Oryza sativa) seeds pretreated with ABA had improved salttolerance, with lower Na + and Cllevels and Na + /K + ratio as well as higher K + and Ca 2+ concentrations, proline accumulation, soluble sugar content, and grain yield [13].
Mounting evidence suggests that ABA plays contrasting roles during defense responses to pathogens [14,15]. It interfered at multiple levels with pathogen stress signaling and suppressed or promoted phenylalanine ammonia lyase (PAL) activity [16,17]. The different plant species assayed in these studies can account, in part, for the various observed effects of ABA on pathogen resistance. ABA is associated predominantly with pathogen susceptibility in tomato, but both negative and positive effects were reported in Arabidopsis [18]. Moreover, the synergistic or antagonistic roll of ABA in disease resistance is universally acknowledged to depend on the type of pathogen, its specific way of entering the host, the timing of the defense response, and the type of tissue affected, as well as the concentration of ABA [19]. ABA responses during biotic and abiotic stresses primarily affect levels of gene expression [7]. ABA-induced genes mainly encoded proteins associated with dehydrins and enzymes that detoxify reactive oxygen species, enzymes of compatible solute metabolism, a variety of transporters, regulatory proteins such as transcription factors (TFs), protein kinases and phosphatases, and enzymes involved in phospholipid signaling. ABA-repressed genes usually encode proteins that are associated with growth, including cell wall, ribosomal, plasma membrane, and chloroplast proteins [8].
To comprehensively understand the functions of ABA and related genes, high-throughput screening techniques have proven quite effective [20][21][22][23]. Matsui et al. (2008) completed an Arabidopsis transcriptome analysis after 2 and 10 h of drought, cold, high-salinity, and ABA-treatment conditions using a tiling array, and the results indicated that approximately 16% of the genes in the transcriptomic library were significantly regulated by ABA. In addition, there was greater crosstalk between ABA and drought and highsalinity stress signaling processes than between ABA-and cold-stress signaling. Similar results were discovered in expression profiles of rice genes under various stresses and ABA application using a cDNA microarray and RNA gelblot analyses [20]. The impact of ABA on gene expression far exceeds that of other plant hormones. About 3,000 genes detected by microarray were significantly regulated by ABA in Arabidopsis seedlings; this number was twice and 23 times as many genes as were influenced by methyl jasmonate and gibberellin, respectively [24].
With the rapid development of next-generation sequencing technology, RNA deep-sequencing (RNA-seq) becomes more efficient and less expensive [25]. The results of RNA-seq are highly reproducible, both technically and biologically [26,27]. In recent years, numerous transcriptome data sets have been produced, made publicly available, and reanalyzed by other researchers. However, so far, few RNA-seq analyses of the genome-wide responses of genes to ABA in plants have been reported.
We investigated ABA responses in tomato, and here report the results of a comparative transcriptome analysis of exogenous ABA-treated tomato leaves using RNA-seq technology. The goals were to (i) construct a tomato leaf transcriptome; (ii) compare and analyze the transcripts in control and ABA-treated plants; and (iii) gain insight into stress tolerance and pathogen-resistance induced by ABA in tomato. This study presents the transcriptome of tomato leaves responding to ABA and provides a genetic resource that can be used for crop improvement.

Global transcriptome analysis
Tomato seedling sprayed with 7.58 μmol L -1 ABA solution were chosen to study ABA responses. Twenty-four hours later, the young third leaves of randomly-selected plants from both the ABA-treated and control groups were collected, labeled as a1d (1 d after ABA-treatment) and c1d (1 d, control), respectively. The tomato ABA/control RNA samples were used for deep sequencing on an Illumina HiSeq 2000 platform. Each read in the Solexa paired-end (PE) sequencing was 101 bp in length. Sequencing generated 266.98 million reads, a total raw dataset of 26.96 Gb. After trimming, 212.78 million clean reads remained, corresponding to 20.95 Gb clean data. The dataset of each sample, including c1d and a1d, was represented by over 100 million reads, a read density sufficient for the subsequent quantitative analysis of genes. The sequence reads were aligned to the tomato reference genome using SOAPaligner/soap2 software (http://soap.genomics.org.cn), allowing two base mismatches [28]. Of the total reads, 81.97% matched either to a unique (36.53%) or to multiple (45.44%) genomic locations (Table 1).
After aligning and assembling, 27,597 genes and 35,051 transcripts were identified in the transcriptome. Of those spliced transcripts, 9,138 (26.07%) matched completely with the annotated tomato genome, 12,501 (35.67%) were potentially novel isoforms, 3,339 (9.53%) were unknown or intergenic transcripts, and 1,634 (4.66%) mapped to the complementary strand of an annotated gene. These results suggested that some transcripts were probably generated from alternative mRNA splicing or were new transcripts.
The genes and transcripts resulting from Illumina sequencing were merged with the annotated reference genome to generate 37,633 genes and 51,606 transcripts using Cuffmerge [29]. We removed transcripts with lengths less than 150 bp, and the remaining 50,770 transcripts (98.38% of all transcripts), corresponding to 37,093 genes, were used for functional annotation and expression profiling in each sample. Transcripts of length 150-200 bp accounted for 8.14% of the total, those 200-600 bp for 32.53%, those 600-1000 bp for 20.28%, those 1000-1600 bp for 20.72%, those 1600-2200 bp for 9.44%, those 2200-3000 bp for 5.41%, and those >3000 bp for 12.27%, as shown in Figure 1.

Annotation, functional classification, and KEGG analysis of all detected transcripts
Of the 50,770 transcripts, 45,704 (90.02%) possessed open reading frames (ORF), while ORFs in the other 5,066 (9.98%) were not predicted (Additional file 1: Table S1). Of all transcripts, 47,877 were described in the tomato genome, and 42,583 had homologs in the NCBI nonredundant (NR) protein database. Many of these genes were reported to respond to drought, cold, high salinity, or ABA, and included both regulatory and functional proteins. Transcription factors (1,799 transcripts) were the largest group of regulatory proteins (Additional file 1: Table S1). This result was consistent with the previous reports in Arabidopsis and rice using microarray and tiling arrays analysis [20,23].
For the global functional analysis, all identified tomato transcripts were classified into different functional categories using Blast2GO (version 2.3.5) [30]. A total of 31,107 transcripts could be annotated in Gene Ontology (GO) and were classified into 57 functional groups, including 23 groups in biological process, 19 in cellular component, and 15 in molecular function ( Figure 2). Within biological process, "metabolic process" (GO: 0008152) with 17,827 transcripts and "cellular process" (GO: 0009987) with 18,169 transcripts were predominant. In the category of cellular component, the three main groups were "cell" (GO: 0005623, 18,493 transcripts), "organelle" (GO: 0043226, 14,423 transcripts) and "cell part" (GO: 0044464, 18,493). The categories "binding" (GO: 0005488) and "catalytic activity" (GO: 0003824) were most common in molecular function, represented by 15,938 and 15,869 transcripts, respectively.

Analysis of expressed genes
In total, 39,671 (78.14% of 50,770) transcripts expressed in the tomato leaf transcriptome (Additional file 1: Table S3), including 38,626 in c1d and 37,989 in a1d, respectively. There were 1,682 transcripts expressed only in the ABAtreated group (a1d), 1,045 transcripts expressed only in the control group (c1d), and 36,944 transcripts expressed in both libraries, which indicated that ABA may activate or repress a fraction of unique transcripts.
Based on expression level, the transcripts were divided into five groups ( Table 2). The group with expression levels of 1-10 FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) was largest, representing 44.64% of all transcripts, while those with expression levels of <1 FPKM accounted for 28.87%, those with 10-100 FPKM 23.49%, those with 100-1000 FPKM 2.73%, and those with >1000 FPKM only 0.26%. Among all expressed transcripts, 17,959 showed no changes in expression levels (|log 2 fold-change (log 2 FC)| < 0.25). The expressions of 21,712 (54.73%) transcripts were altered by exogenous ABA (14,559 were up-regulated, 7,153 were down-regulated) (Additional file 1: Table S3). Of the transcripts with altered expression, 2,787 (12.84% of 21,712) were significantly changed (|log 2 FC| ≥ 1 and false discovery rate (FDR) < 0.05) and labeled differentially expressed genes (DEGs), including 1,952 that were up-regulated and 835 that were down-regulated ( Figure 4). Among the DEGs, there were more up-regulated transcripts than downregulated ones, indicating that many genes responded positively to ABA treatment. This result was similar to those reported in previous studies in Arabidopsis and rice [31][32][33].
Molecular functional classification based on GO was also performed on the 2,787 DEGs. The up-and downregulated DEGs annotated in GO were grouped into 54 and 38 groups on GO.level3, respectively (Additional file 1: Table S4) and into ten classes on GO.level2. The most common categories were "catalytic activity" (654 upregulated, 295 down-regulated) and "binding" (661 and 260, respectively), followed by "transmembrane transporter activity", "substrate-specific transporter activity", "transcription factor activity", "translation factor activity", and "signal transducer activity".
The discovery of PYR/PYL/RCAR receptors in these samples greatly advances the understanding of the ABA signaling pathway. Among 18 PYLs in the tomato transcriptome, the expression level of eight did not differ (|log 2 FC| < 0.25) in ABA treatment compared with control. The other 10 PYLs showed slight differences in their transcript abundance, from 0.55 to 1.43 fold, with five up-regulated and five down-regulated genes. Two up-regulated PYLs showed moderate transcript level differences (TCONS_00012484, change from 14.16 FPKM in c1d to 18.07 in a1d, TCONS_00036624, from 14.79 to 17.26 FPKM). The other eight differentially-expressed PYLs were found at low abundances (FPKM <10). The results indicated that exogenous ABA had a slight impact on PYL expression because the ABA-mediated signaling cascade is initiated by the perception of ABA receptors [35]. The slight change in upstream ABAreceptor gene expression could initiate responses in a series of downstream genes. RCARs/PYR1/PYLs were recently found to comprise 14 members in Arabidopsis thaliana [6,36]. All PYL family members are ABA-binding proteins and that regulate the target PP2Cs in a combinatorial manner [37]. However, we detected 18 PYL members in this tomato transcriptome, enriching the PYL receptor gene family. Among the 23 PP2Cs in this study, 13 transcripts changed their expression levels (|log 2 FC| ≥ 0.25), treated by ABA, including five that were up-regulated and eight that were down-regulated. Of them, only one transcript with low levels (TCONS_00014920) was significantly higher than the control by 2.05-fold (from 0.37 to 0.74 FPKM). Three PP2Cs showed significantly decreased expression in ABAtreatment (a1d) by 2.42-fold (TCONS_00014686, from Figure 4 Differential expression analysis of all transcripts in the control (c1d) and ABA-treatment (a1d) libraries. Transcripts that differed by less than 20% (|log 2 FC| < 0.25) were assumed to not change in expression level. Transcripts that satisfied the conditions of "FDR < 0.05" and "|log2 fold-change (log 2 FC)| ≥1" were considered significantly differentially-expressed genes (DEGs). Other transcripts are noted as "up-regulated slightly" or "down-regulated slightly".  [8]. In the absence of ABA, PYR/PYL/RCARs are not bound to PP2Cs, and PP2C activity is high, which prevents SnRK2 activation. In the presence of ABA, the combination of ABA and PYR/PYL/RCARs bind to and inhibit PP2Cs, which allows the accumulation of phosphorylated SnRK2s and subsequent phosphorylation of ABFs [8].
Among the 12 SnRK2 transcripts, five differed slightly in level compared with control. Of them, one showed increased expression (TCONS_00011840 from 4.41 to 5.34 FPKM) and four showed decreased expression (TCONS_00017720, TCONS_00017721, TCONS_00044604, TCONS_00046427). Although TCONS_00017721 had moderate transcript abundance (FPKM from 40.56 to 26.80 FPKM), the other three down-regulated SnRK2s were expressed at low levels (FPKM < 1.04). SnRK2s are categorized into subclasses I, II and III. Subclass I members are not activated by ABA, and subclass II and III members are weakly and strongly activated by ABA, respectively [38]. Among 10 SnRK2 members in Arabidopsis and rice, SAPK8, SAPK9, and SAPK10 were in subclass III, and SAPK1 to SAPK7 were in subclass I and II [38]. In this transcriptome, the only induced SnRK2 (TCONS_00011840) that was annotated as serine/threonine protein kinase SAPK8-like protein was of subclass III. The other four repressed SnRK2s were possibly SAPK1, SAPK3, or SAPK7 based on their NCBI blastx matches. Kobayashi et al. (2004) demonstrated that only three SnRK2 family members were activated by ABA but all were activated by hyperosmotic stress, and there were no members activated only by ABA. This study also suggested that a small percentage of SnRK2s participated in the signaling transduction under exogenous ABA stress.
In the 18 ABFs, seven were up-regulated and three were down-regulated. TCONS_00038922 (from 0.16 to 0.43 FPKM) and TCONS_00038921 (from 0.43 to 0.10 FPKM) were significantly altered. This result indicated that ABFs were mainly positive regulators of ABA response. The members of the ABF gene family belongs to a subfamily of bZIP TFs [6]. Most ABFs are known induced in vegetative tissue under drought stress and ABA treatment in Arabidopsis [9]. However, in the transcriptome, ABFs were expressed at low levels in both libraries, so transcription regulator factors may target a cascade action [39].
Overall, these results agreed with what is known of the ABA regulation pathway [6,8,40]. Thus exogenous ABA can activate the signaling pathway, and it affected PYLs, PP2Cs and SnRK2s, and ABFs. Moreover, members of these gene families had different responses to ABA, indicating that they probably played positive or negative regulatory roles in the ABA pathway. Furthermore, many genes in these families had no differences in expression abundance, implying that there exists functional redundancy in the ABA pathway or that these genes participate in an ABAindependent pathway [8].

Transcription factors induced by ABA
With respect to TFs, the results of this study supported those of previous research. The TF families of bZIP (111 transcripts), bHLH (146), MYB (212), AP2/ERF (84), NAC (63), and WRKY (81) were identified in this study. Moreover, 41 Heat shock factors (HSFs) and 61 MADS-boxes were also detected ( Table 3). In the bZIP family, ABA-responsive TF (four transcripts) and bZIP (107 transcripts) were two subgroups. Dehydration-responsive element-binding protein (DREB, 11) and bHLH (135) belonged to the bHLH family. The MYB family included MYB (193) and MYC (19). The AP2 family was divided into three subgroups, including AP2 (26 transcripts), AP2/ERF domain-containing TF (55), and RAV1 (3). Reportedly, the most common classes of regulatory sequences induced by ABA are the G-box ABREs recognized by members of the bZIP family, and many ABAregulated genes also contain binding sites for proteins of the MYB and MYC families. NAC, WRKY, bHLH, and Zn-finger classes were also found to participate in some ABA responses [8].
Although TFs were up-regulated after ABA treatment, they had low abundances in both libraries (FPKM < 10) ( Table 3, Additional file 1: Table S5). The TFs interacted with cis-elements in the promoter regions of several responsive genes and thus control the expressions of many downstream genes, triggering cascade reactions of many physiological processes and controlling biochemical reactions in plant cells [39]. Thus, a slight alteration in the transcript abundance of TFs can result in a substantial change in downstream gene expression and physiological responses. Therefore, although in this tomato transcriptome ABA affected TFs with low transcript abundances, the potential results cannot be neglected. Exceptionally, dehydration-responsive element binding proteins (DREBs) that bind to drought-response elements (DRE) cis-elements were substantially down-regulated (five down-regulated and two up-regulated). Reportedly, most DREBs are independent of ABA [41].
HSP90s and HSP70s are important members of the HSP protein family. The major role of HSP90 is to manage protein folding [43,44], and it also plays key roles in signal-transduction networks, cell-cycle control, protein degradation, and protein trafficking [45]. HSP70s function in preventing aggregation, assisting refolding, protein import and translocation, signal transduction, and transcriptional activation [42,46]. In this transcriptome, three up-regulated and one down-regulated HSP90s were weakly expressed. Two of the three up-regulated HSP90s had high transcript abundances (TCONS_00021180 from 95.04 to 110.50 FPKM, TCONS_00032910 from 868.34 to 966.19 FPKM), while the down-regulated HSP90 was expressed at low abundance. Among 36 HSP70s, 16 transcripts differed in expression level, with eight up-regulated and eight down-regulated, although the differences in expression were slight. The up-regulated TCONS_00037214 (from 989.52 to 1116.60 FPKM) and the down-regulated TCONS_00017664 (from 105.12 to 80.45 FPKM) showed high and moderate transcript abundances, respectively. The other changed HSP70s were expressed at low abundances.
Although many HSP70s and HSP90s were slightly upregulated, some showed high transcript abundances. For a highly-abundant transcript, a slightly changed ratio may substantially alter its abundance and induce physiological changes. This research indicated that HSP70 and HSP90 may play an important role in ABA responses.
Small HSPs are a large family present ubiquitously in all organisms. Plants have many more sHSPs than other eukaryotes [47]. In 24 sHSPs, 15 slightly changed their expression levels by ABA, including four up-regulated and 11 down-regulated transcripts. All up-regulated sHSPs had low transcript abundances in both the ABA treatment and control. Among the down-regulated sHSPs, there were five with moderate abundance (10 < FPKM < 100) in both libraries (TCONS_00036382, TCONS_00026922, TCONS_00028974,  TCONS_00035725, TCONS_00023776), and the other six had low abundances. sHSPs bind denatured proteins and prevent their irreversible aggregation, so they are referred to as 'paramedics of the cell' [48]. Therefore, they have a crucial and fundamental role in plant cell biology [47,48], and contribute to the molecular adaptation to a wide range of environmental stresses, including heat, cold, drought, salinity, and oxidative stress [42]. At present, studies of ABA responses to HSP are limited. HSP101 mRNAs in wheat leaves were reported to be induced by 2 h of dehydration and 50 μM ABA [49]. In an examination of nine rice HSP genes under abiotic stresses and ABA treatment, only OsHSP71.1 was induced by ABA while OsHSP24.1 was suppressed [50]. In this study, several HSP70s and HSP90s were up-regulated, while most sHSPs were down-regulated. These observations imply that the HSP genes may play different roles in plant development and abiotic stress responses. Furthermore, the relationship between ABA and HSPs will be interesting to investigate in the future.
In a stress environment, the first step in O 2 reduction produces diffusible hydroperoxyl (HO 2 − ) and superoxide (O 2 − ) radicals [52]. Meanwhile, SODs are induced to rapidly disproportionate O 2 − into oxygen and H 2 O 2 . However, of nine SODs, one with high abundance (TCONS_00025638, from 204.01 to 265.26 FPKM) was slightly down-regulated, and another (TCONS_00025639, from 0 to 0.012 FPKM) showed slightly increased expression level. The other six SODs showed no changes (|log 2 FC| < 0.25). It was also found that the expression levels of SODs showed slightly change after ABA treatment in Arabidopsis [23].
CAT reduces H 2 O 2 directly to water and oxygen. In this study, the most of the genes encoding CAT were up-regulated by ABA, suggesting the importance of CATs in scavenging ROS under ABA stress. Among the seven identified CAT genes (five up-regulated, two unchanged), the abundances of TCONS_00003149, TCONS_00003243, TCONS_00049411, TCONS_00051177, and TCONS_00003244 in ABA treatment were higher than in control by 3.33, 2.20, 1.52, 1.45, and 1.23 fold, respectively. In particular, TCONS_00049411 (from 115.88 to 176.01 FPKM) and TCONS_00051177 (from 91.12 to 131.86 FPKM) had high abundance, and TCONS_00003243 was expressed a moderate levels (from 9.43 to 20.70 FPKM).
There were 133 transcripts encoding Trxs (25 upregulated, 29 down-regulated) and 10 transcripts encoding PrxRs (one up-regulated, five down-regulated) in the PrxR/Trx pathway. All up-regulated Trxs were expressed at low abundance in both libraries (FPKM < 10), except TCONS_00023279 and TCONS_00024016. Among the down-regulated PrxR and Trx genes, more than half were moderately or highly abundant.
The results of this study showed that most transcripts of APX, MDAR, DHAR, GR, GPX, GLR, PrxR, and Trx were down-regulated or did not change in expression level. These enzymes reduce oxidative damage by H 2 O 2 in cells directly or indirectly, depending on the GSH redox state, which is catalyzed by GR [54][55][56]. Previous research demonstrated that H 2 O 2 and the GSH/GSSG couple may interact with other signaling pathways that regulate the expression of antioxidant genes during stress responses [57]. However, this study showed that exogenous ABA repressed the expression of the genes encoding enzymes related to the GSH/GSSG couple, suggesting that ABA signaling may be independent of GSH redox signaling.
Furthermore, more than half of the GSTs, which may catalyze the conjugation reaction of GSH and protect SH-containing enzymes against oxidation, were induced by ABA [58]. ABA probably also caused the downregulation of genes encoding enzymes that depend on GSH. POD catalyzes an early step in lignin synthesis in the phenylalanine metabolism pathway, which can increase the rigidity of cell walls to protect plants from pathogenic bacteria [59] and allowing them to resist wilting and mechanical weakness imparted by abiotic stresses [60]. This study found 33 ABA-induced POD genes that maybe participated in these processes.

Genes related to pathogen resistance
The phytohormone ABA plays a multifaceted and pivotal role in plant immunity. In the whole transcriptome, a large number of genes related to disease resistance were detected, including 16 PAL transcripts (six up-regulated, three down-regulated), five polyphenol oxidase (PPO) transcripts (four up-regulated, one unchanged), 62 beta-1, 3-glucanase (GLU) transcripts (14 up-regulated, 12 downregulated), and 36 chitinase transcripts (15 up-regulated and four down-regulated). There were more up-regulated isoforms than down-regulated ones, and the expression changes of up-regulated isoforms were greater than downregulated one. The results suggested that ABA has a greater positive than negative impact on tomato pathogen resistance.
In this transcriptome, many genes in the salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) signal transduction pathways were detected. These pathways play well-established roles in phytopathogen defenses [64]. TGACG/as-1 bindings (TGAs, 26 transcripts in this study) and pathogenesis-related proteins (PR1s, 45 transcripts) are critical component in the SA signaling pathway, while jasmonate resistant 1 (JAR1s, 15 transcripts) and jasmonate ZIM-domain proteins (JAZs, 18 transcripts) are crucial in the JA signaling pathway and ethylene receptors (ETRs, 14 transcripts) and ethylene response factor/ethylene-responsive TFs (ERF/EREBPs, 15 transcripts) are vital to the ET signaling pathway ( Table 4, Additional file 1: Table S5). Most of these genes were up-regulated by exogenous ABA and most were expressed at low abundance.
Salicylic acid, a key phytohormone in plant-pathogen interactions, is typically associated with systemic expression of pathogenesis-related protein genes against biotrophs and hemibiotrophs. This response often culminates in a localized cell death termed the hypersensitive response [65][66][67]. JA and ET predominantly protected Arabidopsis against biotrophic pathogens [68]. In tomato, the comparative transcriptome results indicated that ABA played a positive role in biotic stress resulting from the crosstalk between ABA and SA, JA, and ET.
In general, ABA is thought to play a negative role in plant resistance to both biotrophic and necrotrophic fungi and bacteria, and ABA interacts antagonistically with SA, JA, and ET in various ways [69]. For example, in tomato, the ABA-deficient mutant sitiens was much more resistant to Erwinia chrysanthemi, and exogenous ABA (100 μM) made both the wild type and sitiens more susceptible to E. chrysanthemi [18,70]. ABA (100 μM) decreased the disease resistance of rice to rice blast (Magnaporthe grisea) by suppressing some genes in the SA signaling pathway [71]. However, in this study, ABA largely up-regulated the genes related to biotic stresses, indicating that ABA has potential to promote pathogen resistance in tomato. Furthermore, the concentration of ABA applied here (7.58 μM) was much lower than that used in most previous reports (usually 100 μM). We hypothesize that the ABA concentration is very important in determining whether plant resistance to pathogens will be improved or suppressed.
ABA not only mediates tolerance to adverse environmental conditions, but also affects the resistance to biotic stresses. It plays critical roles in the interrelationship between biotic and abiotic stress responses. Further insights into the core network of ABA signaling in plants will be crucial for genetically engineering and breeding crop species with both improved abiotic stress tolerance and pathogen resistance [72]. ABA has also been recognized to crosstalk with the "developmental" hormones auxin, gibberellic acid, cytokinins, and brassinosteroids in plant immunity [14,15,69]. In this study, the expressions of many genes related to these hormones were also affected by ABA (Additional file 1: Table S1). These genes may directly or indirectly contribute to plant pathogen resistance.
Quantitative real-time-PCR validation of differentiallyexpressed transcripts from RNA-seq To confirm the accuracy and reproducibility of this Illumina RNA-seq result, 25 transcripts were chosen randomly for quantitative real-time (qRT) PCR. Those genes were involved in metabolism, information transfer, or were hypothetical proteins, and included up-regulated, downregulated, and unaffected transcripts. The primer sequences, gene functions, and FPKM and qRT-PCR values are listed in Additional file 1: Table S6. The qRT-PCR results generally agreed (84%) with the changes in transcript abundance determined by RNA-seq, suggesting the reliability of the RNA-seq data.

Conclusion
Genome-wide transcriptome analysis results indicated that exogenous ABA can influence the ABA signaling pathway with the core of PYR/PYL/RCARs-PP2Cs-SnRK2s. Exogenous ABA up-regulated many genes related to stress tolerance and pathogen resistance, including various TFs, HSP90s, HSP70s, CATs, GSTs, PODs, PALs, PPOs and chitinases, as well as the genes TGA, PR1, JAR1, JAZ, ERT and ERF/EREBP in the SA, JA, and ET signaling pathways. These results suggested that ABA has the potential to improve the abiotic and biotic tolerance of tomato. The study extends the knowledge of the complex molecular and cellular events during ABA signaling and of ABA-induced genes. Challenges and opportunities remain in exploring the complex interactions between ABA and defense in whole plants.

Plant material for RNA-Seq
Seeds of tomato cv. Hongtaiyang 903 were planted and grown in plastic pots filled with organic loam in April, 2012, and grown in a glasshouse in Chengdu, Sichuan Province, China. Sixty pots (four seedlings per pot) were used in this experiment. They were watered every other day. After 45 days when the plants had 5-7 leaves, the 60 pots were divided into two groups. The plants in 30 pots were sprayed with 400 mL of 7.58 μmol L -1 ABA solution, and the plants in the other 30 pots were sprayed with the same volume of purified water as a control. Twenty-four hours later, the young third leaves of ten randomly-selected plants from both the ABAtreated and control groups were collected. The leaves of each group were combined, immediately snap-frozen, and stored in liquid nitrogen. These samples were labeled a1d (1 d after ABA-treatment) and c1d (1 d, control).

Read trimming and optimization
For each set of sequencing reads, the sequencing adapters were trimmed using SeqPrep (https://github.com/jstjohn/ SeqPrep), and then low quality bases (Solexa/Illumina quality score < 25) of the 3′ ends were trimmed using inhouse Perl scripts. The remaining high-quality reads were submitted for mapping analysis against the reference genome sequences (ftp://ftp.solgenomics.net/tomato_genome/ annotation/ITAG2.3_release/) using Tophat [73]. The mapped reads were then assembled with Cufflink [29]. The assembled results and original genome annotations were merged and used for further annotation and differentialexpression analysis.
All of these RNA-Seq reads were deposited in Sequence Read Archive database (http://www.ncbi.nlm.nih.gov/Traces/ sra/) under the Accession number of SRR926182 and SRR926185.

Differential expression analysis
The Tophat (http://tophat.cbcb.umd.edu/) and Cufflinks (http://cufflinks.cbcb.umd.edu/) programs provide FPKM (Fragments Per Kilobase of exon model per Million mapped fragments) values within a 95% confidence interval. Differential expression was analyzed and calculated according to the count values of each transcript in the two libraries using edgeR (the Empirical analysis of Digital Gene Expression in R) software [77]. "FDR < 0.05" and "|log 2 fold-change (log 2 FC)| ≥1" were used as the thresholds for judging significant difference in transcript expression. Transcripts with |log 2 FC| < 0.25 were assumed have no change in expression levels.

qRT-PCR verification
qRT-PCR was performed to verify the expression patterns revealed by the RNA-seq study. The purified RNA samples were reverse-transcribed using the PrimeScript RT Reagent Kit with gDNA Eraser (Takara, Dalian, China) following the manufacturer's protocol. Twenty-five transcripts were selected randomly for the qRT-PCR assay. Gene specific qRT-PCR primers (18-20 bp) (Additional file 1: Table S6) were designed using Premier 5.0 software (Premier Biosoft International, Palo Alto, CA). qRT-PCR was performed using SybrGreen qRT-PCR Master Mix (Ruian Biotechnologies, Shanghai, China) in an ABI 7500 FAST Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). PCR conditions were 2 min at 95°C, followed by 40 cycles of heating at 95°C for 10 s and annealing at 60°C for 40 s. Three replicates were performed, and the amplicons were used for melting curve analysis to check the amplification specificity. The relative expression level of each gene was calculated as 2 -(ΔΔCt) and the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase gene (Accession No. U93208) from S. lycopersicum was used to normalize the amount of template cDNA added in each reaction.

Additional file
Additional file 1: Table S1. Annotation of transcripts according to the NCBI non-redundant, Gene Ontology, and STRING databases and the tomato genome. Table S2. KEGG pathway mapping. Table S3. Differently-expressed transcripts in the control (c1d) and ABA-treatment (a1d) libraries. Table S4. Molecular functional classification of the significantly differentially-expressed genes (DEGs) based on GO.level2 and GO.level3. Table S5. Genes related to ABA signaling transduction, transcription factors, heat shock proteins (HSPs), reactive oxygen scavenging enzymes, and pathogen resistance.