Tomato Genome-Wide Transcriptional Responses to Fusarium Wilt and Tomato Mosaic Virus

Since gene expression approaches constitute a starting point for investigating plant–pathogen systems, we performed a transcriptional analysis to identify a set of genes of interest in tomato plants infected with F. oxysporum f. sp. lycopersici (Fol) and Tomato Mosaic Virus (ToMV). Differentially expressed tomato genes upon inoculation with Fol and ToMV were identified at two days post-inoculation. A large overlap was found in differentially expressed genes throughout the two incompatible interactions. However, Gene Ontology enrichment analysis evidenced specific categories in both interactions. Response to ToMV seems more multifaceted, since more than 70 specific categories were enriched versus the 30 detected in Fol interaction. In particular, the virus stimulated the production of an invertase enzyme that is able to redirect the flux of carbohydrates, whereas Fol induced a homeostatic response to prevent the fungus from killing cells. Genomic mapping of transcripts suggested that specific genomic regions are involved in resistance response to pathogen. Coordinated machinery could play an important role in prompting the response, since 60% of pathogen receptor genes (NB-ARC-LRR, RLP, RLK) were differentially regulated during both interactions. Assessment of genomic gene expression patterns could help in building up models of mediated resistance responses.


Introduction
Disease resistance in plants depends on the ability of the host to recognize pathogens and initiate defense mechanisms that limit infection. A basal type of immunity in plants is conferred by the recognition of conserved microorganism-associated molecular patterns (MAMPs) by specific pattern-recognition receptors (PRRs) that protect hosts against non-specialized pathogens. Plants are also capable of establishing immune responses by using pathogen receptors (mainly Nucleotide binding-ARC/leucine-rich repeat proteins NB-ARC-LRR but also Receptor Like Proteins RLP; Receptor-Like Kinase RLK) known as resistance (R) proteins able to recognize the presence of pathogen effector molecules and to activate effector-triggered immunity (ETI). This response is typically associated with programmed cell death of the infected cells and the production of antimicrobial molecules [1]. Modulation of hormone pathways is required to restrict pathogen invasion, re-allocate resources, control cell death and modify plant architecture [2]. In addition, systemic-acquired resistance, which immunizes against subsequent infections, could also occur [1].
Tomato (Solanum lycopersicum) has served as an important model system for studying the genetics and molecular basis of resistance mechanisms in plants. The breadth of pathogen classes affecting tomato underscores the importance of tomato pathosystems as amenable models for studying the plant immune system [3]. Plants build up appropriate defense responses without draining energy resources to unsustainable levels through the cross-talk and finetuning of different defense pathways [4]. Identification of host genes involved in resistance responses could be helpful both for understanding plant defence mechanisms against pathogens and building up a biological model of plant-pathogen interaction.
Fusarium oxysporum f. sp. lycopersici (Fol) is one of the main diseases affecting tomato, causing vascular wilt disease by colonizing the xylem vessels of roots and stems. Three physiological races, termed 1, 2 and 3, have been identified. The FOL AVR genes (AVR1, AVR2 and AVR3) are carried in different combinations in different FOL races. AVR1 ( = SIX4) is unique to race 1 whereas AVR2 ( = SIX3) is found in races 1 and 2. AVR3 ( = SIX1), which exists in all races I, I2 and I3 are known confer resistance to FOL in tomato. Race 1resistant cultivars (I i2 i3), races 1 and 2-resistant cultivars (I I2 i3), and races 1, 2 and 3-resistant cultivars (I I2 I3) have been developed [5]. To date only I-2 gene, a coiled-coil/nucleotide binding-ARC/leucine-rich repeat (CC-NB-ARC-LRR) protein, that confers resistance to race 2 of Fusarium oxysporum f. sp. lycopersici has been cloned [6;7]. Responses mainly involves the callose deposition, the accumulation of phenolics and the formation of tyloses (outgrowths of xylem contact cells) and gels in the infected vessels [8]. Interaction between tomato and Fusarium oxysporum f. sp.lycopersici has been investigated in depth, becoming a model system for disease resistance response [9].
Tomato mosaic virus (ToMV) is a positive-sense ssRNA virus belonging to the Tobamovirus genus. ToMV infects tomato plants systemically, causing mosaic symptoms, which are characterized by intermingled light and dark green regions [10]. The Tm-2 gene of tomato and its allelic form, Tm-2 (2), confer resistance to Tomato mosaic virus (ToMV) and encode members of the CC-NB-ARC-LRR protein class [11].
These two pathogens have different infection strategies: ToMV is a biotrophic leaf pathogen that causes tomato leaf mold, whereas Fol is a saprotrophic soilborne vascular pathogen that challenges several Solanaceae. The aim of this study was to compare tomato global transcriptional profiles in response to host attack by ToMV and Fol in order to identify genomic differences and similarities in incompatible interactions between a foliar and a vascular pathogen, making use of the recently sequenced tomato genome. First we examined global tomato transcriptional profiling during both incompatible interactions to compare transcriptional changes occurring in tomato plants. Then we performed in-depth gene annotation, highlighting enriched GO categories and any metabolic perturbations arising. Finally, we explored the genome arrangement of expressed genes along the chromosome in order to connect genomic and transcriptional events.

Plant Material and Inoculation Protocol
For our experiments we used the Rosso Delta tomato (S.lycopersicum) variety, resistant to Fol races 1 and 2 and ToMV race 0. Tomato plants were inoculated with Fol strain ATCC 16605 containing the Avr1 and Avr2, (Plant Research International, NL). Plantlets were infected at the stage of expanded cotyledons: roots were cut and dipped in a suspension at a concentration of 1610 6 conidia/ml, according to CPVO technical protocol TP/044/3 [12] (Experiment I). Inoculation with Tomato Mosaic Virus was carried out with the strain GM6s of ToMV (Plant Research International, NL), according to the method reported in the CPVO technical protocol TP/044/3 [12] (Experiment II). Plants with expanded cotyledons were inoculated with the sap obtained from infected desiccated tomato leaves. Viral transmission was ensured by mechanical tissue abrasion caused by diatomaceous powder.

Sample Collection
Two days after treatment, infected and non-infected tomato leaf samples of two independent replicas of Experiment I and II were collected. Roots and leaves were removed from the plants, weighed and immediately frozen in liquid nitrogen and stored at 280uC. RNA was isolated from whole plants using the Rneasy Plant Kit according to the manual instructions (Quiagen Valencia, USA). RNA sample concentration was determined using a Nanodrop photometer. An Agilent 2100 Bioanalyzer was used to check the quality of RNA.

Chip Design and Microarray Hybridization
Transcriptomic analysis was performed using a 90 K Toma-tArray 2.0 microarray synthesized using the CombiMatrix platform at the Plant Functional Genomics Center of the University of Verona. The chip (TomatoArray2.0) carries 25,789 non redundant probes (23,282 unique probes and 2,507 probes with more than one target) randomly distributed in triplicate across the array, each comprising a 35-40-mer oligonucleotide. The source of sequence information included tentative consensus sequences (TCs) derived from the DFCI Tomato Gene Index Release 12.0 and expressed sequence tags. Eight bacterial oligonucleotide sequences provided by CombiMatrix, eight probes designed on eight Ambion spikes and 40 probes based on Bacillus anthracis, Haemophilus ducreyi and Alteromonas phage sequences were used as negative controls. Complete description of the chip is available at the Gene Expression Omnibus under the platform accession GPL13934. Microarray analysis was used to investigate tomato gene expression profiles two days after the infection with Fol and ToMV, comparing with the uninfected control profile. After checking the quantity and quality by spectrophotometry using NanoDrop 1000 (Thermo Scientific), total RNA (2 mg) was amplified and labeled using the RNA ampULSe kit (Kreatech). For array hybridization, according to the manufacturer's recommendations, 4 mg of labeled aRNA was employed. Cy5 labeled aRNA was hybridized with the microarray at 45uC for 16 h; then arrays were washed with hybridization Solution at 45uC for 5 minutes and 36SS, 0.56SSPET, PBST wash and PBS wash at room temperature. Pre-hybridization, hybridization, washing and imaging steps were carried out according to the manufacturer's instructions. Microarrays were stripped for reuse (up to three reuses per chip) with the CombiMatrix CustomArrayTM Stripping Kit according to the manufacturer's instructions (product number 610049). After hybridization and washing, the microarray was scanned using a Perkin Elmer Scan Array 4000XL (software ScanArray Express Microarray Analysis System v4.0).

Data Analysis
Gene expression levels corresponding to eight microarrays were processed using R software (R Core Team 2013 [13], URL http://www.R-project.org/.) and the limma package [14]. Raw data were investigated for quality assessment, preprocessed and normalized using a suitable subset of control probes available in the CombiMatrix array and a quantile normalization technique. Three technical replicates within each array and two biological replicates were employed to assess differential expression for each experiment (Experiments 1 and 2) to compare the different experimental conditions (non-inoculated vs inoculated) using a linear model for microarray [15]. The significance of the differential expression was assessed taking into account the multiple testing setting and controlling the False Discovery Rate (FDR) at FDR = 0.01 and the empirical Bayes moderated tstatistics available in the limma package. Linear model fitting was used to assess differential expression of the two pairs of contrasted experimental conditions (non-inoculated vs inoculated), taking into account the presence of technical replicates as well as biological replicates. Differential expression data are deposited in Gene Expression Omnibus (GEO) database repository with accession number (GSE52336).

Annotation of the Gene Chip Probes
An in-house pipeline was developed to annotate tomato tentative consensus sequences (TCs) used to develop CombiMatrix CustomArrayTM probes. The queried tomato genes were identified by mapping TC sequences to the tomato CDS sequence using BLASTn (E-value 1e-3). The latest version of the tomato gff3 annotation files was parsed to extract the cds sequences of genes probed. Blast2GO http://blast2go.bioinfo.cipf.es/was used to provide automatic high-throughput annotation, gene ontology mapping and categorization of tomato proteins identified. An expectation value threshold of 1e-6 in BLASTp analysis was performed. Blast2GO was used for the statistical analysis of GOterm frequency differences. The enrichment analysis of the GO terms was based on Fisher's exact test and corrects for multiple testing. A cut-off FDR value of 0.05 was used. An interactive graph of the GO category enrichment REVIGO htt://revigo.org/ and a graphic representation Cytoscape htt://cytoscape.org// were also performed.

Identification of Differentially Expressed Genes Induced by F. oxysporum f. sp. lycopersici (Fol) and ToMV Inoculation
The transcriptional responses of resistant tomato seedlings, inoculated with Fol and ToMV, were evaluated by querying 15,734 tomato genes. This analysis allowed us to compare changes occurring during both incompatible interactions with a vascular fungus and with a virus. Differentially expressed tomato genes (FDR-adjusted p-value ,0.01) were identified at 2 DPI, comparing inoculated and non-inoculated plants. The time point was selected in order to identify genes involved in the initial stages of the defense against the pathogens.
In Fol incompatible interaction, transcriptome variation resulted in 3,753 differentially induced genes. In particular, 2,392 genes (about 64%) were up-regulated (Fig. 1A), indicating considerable gene activation during infection. As for the vascular pathogens, 3,501 transcriptional changes were monitored in tomato upon inoculation with ToMV, of which about 2000 (52%) were overexpressed. When the total number of differentially regulated genes between the two incompatible interactions were compared, roughly half (2,205 genes) overlapped, of which a small subset of 131 genes had the opposite expression direction during the two interactions. As shown in Fig. 1B, a marked host gene induction was evident, since more than half of the differentially expressed genes appeared to be induced in both interactions. Transcriptional changes during the two tomato-pathogen interactions were further investigated in order to find the best strategy for building a biological model for a plant-pathogen system.

Gene Set Enrichment Analysis
We performed a GO-term annotation analysis of all transcripts identified by gene chip probe matching. Through this analysis, we were able to assign functional annotations to 15,378 chip transcripts. In order to facilitate the formulation of biological hypotheses, a categorization of gene expression data was performed. Indeed, the Gene Ontology research effort was extremely useful for structuring data description.
In the tomato-ToMV interaction 226 enriched GO terms were detected, 70% of which belonged to a biological process (P) category, while in the fungus interaction fewer enriched categories (185) were identified, 78% of which belonged to a biological process (P) category ( Fig. 2A). Comparing the total number of GO categories between the two incompatible interactions, roughly 60% overlapped.
A heatmap obtained through hierarchical analysis evidenced four specific GO term clusters. Clusters 1 and 4 included GO ontology terms that did not show significant differences in two experiments. Cluster 2 contained about 90% of the enriched GO terms. Subcluster 2b consisted of specific GO terms of tomato-Fol interaction, whilst subclusters 2a and 2d as well as cluster 3 included tomato-ToMV interaction specific ontology terms (Fig. 2B).
To reduce redundancy of functional categories between the enriched GO terms, we used the semantic similarity approach. The gene ontology interactive graph-based network produced ( Fig. 3) encapsulated functional homology between genes of the two tomato-pathogen interactions. In the tomato-Fol interaction, the network comprised 53 nodes and 310 edges, while for the tomato-ToMV interaction the network was more complex (58 nodes and 383 edges).
The network of tomato-Fol interaction showed that the biotic stimulus node (GO:0009607) correlated with the five internal nodes with gene ontology terms related to external stimuli response (GO:0009637, GO:0009746, GO:0009628, GO:0042742, GO:0042744). The GO:0009607 category includes a large number of genes that play a role in signal transduction and regulation of gene expression in the defense response, representing about 10% of differentially expressed genes in tomato-Fol interaction. The external stimuli sub-network was also connected with hydrogen peroxide catabolism node (GO:0042744). It contained 28 genes, 70% of which were overexpressed. In particular, four haem peroxidases (Solyc01g006290.2.1; So-lyc01g006300.2.1; Solyc06g050440.2.1; Solyc11g018800.1.1), which play a role in host defense by inhibiting the hyphal extension of invading pathogens, were included in this node. Fig. 3B shows the network obtained with significant GO-term changes in tomato-ToMV interaction. The biotic stimuli GO:0009607 node correlated with four GO-term nodes (GO:0009637, GO:0006950, GO:0010038, GO:0009628), associated to 264 genes. Interestingly, no external link was evidenced for this interaction. These data indicate that incompatible interactions induced both well documented stress-responsive genes as well as unknown genes that might play a role in multi-stress responses. Investigation of single GO categories could evidence GO terms related to each assessed specific response.

Investigation of Perturbed Biological Processes
To obtain an overview of processes involved in fungus and virus incompatible reactions in-depth GO (gene ontology) term analysis was performed. Significant differences were observed in specific regulated biological processes between the two incompatible interactions. Notably, 30 specific enriched GO-term categories for tomato-Fol interaction, and 71 for tomato-ToMV interaction were identified (Table S1and S2 in File S1).
Tomato-Fol interaction revealed changes in cell structure, light perception, iron transport, aromatic amino acid synthesis and ROS response. In particular, we found six specific GO terms associated to homeostatic process (GO:0006873 cellular ion homeostasis; GO:0055080 cation homeostasis; GO:0050801 ion homeostasis; GO:0019725 cellular homeostasis; GO:0048878 chemical homeostasis and GO:0042592 homeostatic process) that included 104 genes. This strongly suggests that homeostasis plays an important role in tomato defense to Fol. Interestingly, the master gene of inflammation, NF-kB (Solyc02g094530.1.1) was also up-regulated in tomato-Fol interaction. This gene is a key player in anti-apoptotic signaling and is able to prevent apoptotic signaling by inhibiting map-kinases. Indeed, tomato map-kinases (Solyc12g019460.1.1; Solyc11g072630.1.1) were found to be heavily down-regulated during this interaction. Table 1 shows 14 auxin-binding genes and two abscisic acid receptors upregulated in tomato-Fol interaction. Photosynthetic and carbohydrate derivative processes, nitrogen biosynthesis, therpene and carotenoid metabolism, and cadmium/ copper function were challenged in tomato-ToMV interaction. Three specific GO terms associated to photosynthetic processes (GO:0009853 photorespiration; GO:0019685 photosynthesis dark reaction and GO:0009773 photosynthetic electron transport in photosystem I) were identified and were investigated further. Table 2 shows the enriched categories associated to photosynthetic processes containing 191, genes over 86% of which were downregulated. Interestingly, photosynthetic gene transcription repres-sion was inversely correlated with the pathogenesis-related gene induction (Table 3). Moreover, a close linkage clearly emerged between the plant carbohydrate status and the resulting plant pathogen interaction. We found 17 carbon fixation enzymes differentially expressed during tomato-ToMV interaction. Furthermore, several genes involved directly in the Calvin cycle were down-expressed (Table 4), and b-fructofuranosidase (So-lyc10g085640.1.1), a key invertase involved in sink source portioning, was overexpressed. We also found up-regulated six gibberellin-modulated genes and two genes involved in SA

Relationship between Genomics and Transcriptional Changes
In order to investigate the chromosomal arrangement of differential expressed sequences identified by tomato chip, a genomic expression map was constructed (Fig. 4). Our results indicated that differentially expressed genes (about 33% of all genes analyzed) during the two interactions were arranged along the chromosomes with different levels of expression. The fungus and virus distribution showed slight differences. Indeed, the number of expressed genes identified for single chromosomes was always greater in tomato-ToMV interaction than in tomato-Fol interaction, with the exception of chromosome 8. Chromosomal areas (on chromosomes 1, 2, 8, 11 and 12) showed a higher density of up-regulated genes in the fungus-interaction. By contrast, on chromosome 8 we found a region with a higher density of downregulated genes in the virus-interaction. In both interactions, the genomic regions rich in differentially expressed genes were located in the vicinity of telomeres, inter alia the richest regions of genes were the long arm of chromosomes 10 and 12.
The genome-wide distribution of differentially expressed genes, based both on chromosome size and gene density was not uniformly distributed along the genome (chromosome size: fungusinteraction x2 = 107; virus-interaction x2 = 117. Gene density: x2 value included between 27 and 170). Fig. 5 and 6 showed the genomic distribution of under-expressed and over-expressed genes during the two interactions. If we use the genome location of the coding genes spotted on the microarray as a reference we could highlight a significant difference in the distribution of overexpressed genes as well as under-expressed ones (Kolmogorov-Smirnov and Cramer von Mises tests in Table S3 and S4 in File S1). Indeed, differentially expressed genes (p-value ,0.10) were evidenced on all chromosomes, except for chromosomes 2, 4 and 7. Stronger evidence of peculiar chromosomal displacement of over-expressed genes may be found on chromosomes 3, 5, 6 and 1 (p-value ,0.05). Significantly differential displacement of underexpressed genes (p-value ,0.05) was detected mainly in chromosomes 5 and 9. This suggested the presence of a dynamic genomic landscapes involved in pathogen response.
For each interaction studied, distribution of expressed specific genes is shown in Fig. 7. Most were located in the same regions, but clear differences could be highlighted. The 151 common differentially expressed genes with opposite direction were mainly arranged on chromosomes 7, 9 11 and 12. On chromosomes 9 and 11, in the vicinity of genes Tm2(2) conferring resistance to Tomato mosaic virus and I and I2 conferring resistance to Fusarium oxysporum, we highlighted differences in gene expression between the two experiments. Eleven genes were specifically expressed in ToMV interaction in this area of chromosome 9 (including: pentatricopeptide repeat, ATP dependent RNA helicase, cytochrome P450 and pectinesterase) and 19 of the I2 region of chromosome 11 in Fol interaction (including: BEL1-like homeodomain protein 8, Calmodulin-like protein and MYB transcription factor and two proteins of unknown function).whereas around of I locus in the 1 Mb region delimitated from TG523 and CT168 [16], 8 specific genes were identified (including a Pumilio-like IPR011989 protein and Serine-threonine protein kinase). Fig. 8 shows differential expressed pathogen recognition genes (NB-ARC-LRR, RLP, RLK) both in tomato-Fol and tomato-ToMV interactions. Pronounced activation of all pathogen recognition gene classes assessed was evidenced in tomato-Fol interaction since 70% of all the differentially expressed genes were up-regulated. Instead, in the tomato-ToMV interaction the number of upregulated and down-regulated pathogen recognition genes was the same: 45% of all the differentially regulated pathogen recognition genes overlapped. Of these genes, 10 had the opposite expression in the two experiments (Table S5 and S6 in File S2).

Discussion
Plant-expressed molecules during pathogen challenging can give insights into the underlying defense mechanisms. Plants have evolved a defense system against microbial pathogens that involves the regulation of gene expression, cascade signaling activation, hormone balancing and synthesis of defensive metabolites [17]. Several microarray experiments have been successfully carried out in the sphere of molecular plant-microbe interaction, shedding light on the mechanisms controlling plant disease resistance and the crosstalk among the signaling pathways involved [4]. More than 3500 genes, out of 15,734 assessed, were activated in tomato Fol and ToMV interactions, more than 60% of which overlapped. Pathogens from different kingdoms deploy independently evolved virulence proteins that interact with a limited set of highly connected plant cellular hubs [18]. Network analysis suggests that plant response to Fol and ToMV infection involved activation or repression of genes implicated in pathways shared by most responses to environmental stimuli. However, GO enrichment analysis evidenced specific enriched categories in both interactions. Response to ToMV seemed more multifaceted, since more than 70 specific categories were enriched versus the 30 detected in Fol interaction. Biotrophic plant pathogens are generally accepted to have a more intricate biological interaction with their host plant than saprophytic plant pathogenic fungi [19]. Once an attack is perceived, plant metabolism must balance potentially competing demands for resources to support defense versus requirements for cellular maintenance, growth and reproduction [20;21;22;23]. In tomato-Fol interaction our investigation evidenced a number of overexpressed genes associated to maintenance of cellular structures and cellular homeostasis. These are very important metabolic activities required by plants to survive fungus-inflicted stresses. Fusarium oxysporum is a necrotrophic fungus [24] that kills host cells prior to the infection through the predicted deployment of toxins and enzymes, which induces cell death. Hence the expression of anti-apoptosis genes could confer resistance to Fol [25]. The NF-kB signaling enhancement found in this work can inhibit apoptosis and combat the effects of oxidative stress by switching down the anti-apoptotic signals mediated by the mapkinase cascade. These events lead us to hypothesize that the resistant plant challenged by Fol induces a homeostatic response to prevent the fungus from killing the cells to obtain its nourishment. Fourteen auxin-binding genes were also activated in this interaction. Auxin signaling is an important phytoregulator for resistance to saprophytic/necrotrophic fungi [26]. The degradation of AUX/IAA proteins allows activation of auxin response factors (ARFs) and the expression of auxin-responsive genes [27]. These protein families can have an enormous number of interactions, capable of fine-tuning specific responses within the auxin signaling pathway [28]. Moreover, the buildup of IAA near the sites of pathogen ingress constitutes one of the main host factors that determine plant resistance to Fusarium wilt [29]. Essential regulators of Jasmonate-ZIM domain proteins (JAZs) in tomato-Fol interaction were activated, suggesting attenuated JA  signaling to reduce the senescence process. Indeed, also a pair of abscisic acid receptors were up-regulated in Fol interaction.
Interactions between multiple components of ABA and the JAethylene signaling pathways modulate defense and stress-responsive gene expression in response to Fol, confirming that these hormones mediate resistance in necrothophic interaction [30;31]. The SA pathway proved to be switched down. Indeed, silencing a tomato gene encoding SA methyltransferase was found to enhance resistance to F. oxysporum f. sp. lycopersici biotic and abiotic stresses [32].    In tomato-ToMV interaction, we detected a reduced photosynthetic activity and elevated carbohydrate catabolism. Pathogen attack may have caused metabolic alterations that induced a hormone response for reprogramming the cellular metabolism. The photosynthesis and assimilatory metabolism were switched off to initiate respiration and other processes required for defense. This occurrence reflected the reallocation of plant metabolites from normal growth processes to defensive functions after the elicitation of induced plant responses by virus infection [33]. A close linkage between the plant carbohydrate-status with the resulting plant pathogen interaction was also evident. During interaction with the virus, the host produced salicylic acid and increased the production of an invertase enzyme that is able to redirect the flux of carbohydrate acquisition. The infection of ToMV could have led to the alteration of the concentration gradient of sucrose in phloem. The plant overexpressed bfructofuranosidase (Solyc10g085640.1.1) to ensure the sucrose transport of source cells to the sink cells. Overexpression of invertase in tomato-virus interaction could be elicited by the expression of pathogenesis-related proteins and by salicylic acid in order to increase resistance to virus infection [34]. In order to ensure the success of defense, this mechanism appears to be crucial. The phenomenon of ''high sugar resistance'' was described long ago [35] and the finding that various pathogenesis-related genes were sugar-inducible [36] supports this hypothesis. We also found six gibberellin, one SAM dependent carboxyl methyltransferase (Solyc01g014320.2.1) and one jasmonate-zim domain (Solyc07g042170.2.1) genes up-regulated, suggesting that GAs activated tomato immune responses to ToMV by modulating the levels of salicylic acid and/or jasmonic acid [37].
Plants rely on the innate immunity and on systemic signals emanating from metabolic alterations [38]. The ETI system seems to prompt response in the right direction thanks to metabolic clues and hormone signaling. Salicylic acid (SA) primarily triggers resistance against biotrophic pathogens, whereas a combination of jasmonic acid (JA) and ethylene (ET) signaling activates resistance to necrotrophic pathogens [39;40]. The metabolic changes associated with defense response in two incompatible tomatopathogen interactions suggested that the response to the specific metabolic alteration (photosynthesis/carbohydrate metabolism and homeostasis activity) in tomato was pathogen-specific and contributed substantially to monogenetic gene resistance. Our results confirm that resistance to pathogens depends on a sophisticated interplay among different biological pathways and that hormonal directionality is critical to the outcome of a response (carbohydrate biosynthesis and homeostatic activity). Several hormones involved in pathogen perception, activation of defense products restricting pathogen invasion, have been identified [2]. How they trigger ETI defense is currently unclear. In the absence of a pathogen, NB-ARC genes are in an autoinhibition state, which is relieved upon pathogen perception [41]. Although different forms of plant immunity share the same signaling mechanisms, they use the same mechanisms in very different ways [42]. The mechanisms that lead to rapid metabolism switch and connection among all the defense pathways is still not clear. Signaling able to fine-tune the defense mechanism could be activated by genes in proximity of the main pathogen receptors.
Genome organization of functional gene networks to tolerate alterations can result in plasticity. In two specific interactions the chromosome distribution of expressed genes showed wide overlapping regions, except for the region holding genes I2 and Tm2 (2). Indeed, genome regions can be enriched in genes with the specific function for fine-tuning gene expression in a compensatory way. In Fol interaction on chromosome 11 region harboring I2 gene a calmodulin gene, a myb factor and a BEL-l like homeodomain protein 8 were specifically expressed. Calmodulin (CaM) plays an important role in sensing and transducing changes in cellular Ca2 + concentration in response to several biotic and abiotic stresses, In ToMV interaction Pentatricopeptide-repeat, ATP-dependent and RNA helicase could be involved in the process of silencing and virus replication.
A high level of NB-ARC genes activation and of other gene classes potentially involved in pathogen recognition was also found, suggesting a host-coordinated reaction of defense machinery to monitor integrity of cellular proteins. The NB-ARC system minimizes the cost of defense for the plant, as multiple NB-ARC-LRR proteins can be maintained at a low level in the absence of a pathogen and rapidly induced under pathogen attack through miRNA regulation [43]. Pathogen-encoded suppressors of RNAsilencing mechanisms might result in the induction of multiple NB-ARC-LRR defense proteins [44]. Investigation of differentially regulated pathogen recognition genes could lead to the identification of specific modulation patterns.
Albeit still fragmented, our depiction provides a global view of the I-I2 and Tm2(2) mediated resistance process, in which the gene expression network may be considered the starting point to construct a genomic model for R-mediated response in the two pathosystems investigated.

Supporting Information
File S1 Contains the following files: Table S1. Specific enriched GO-term categories in tomato-Fol interaction; Table S2. Specific enriched GO-term categories in tomato-ToMV interaction; Table S3. Results of Kolmogorov-Smirnov and Cramer von Mises analysis for testing the differences between the position distribution of over-expressed genes and coding genes spotted on the microarray chip; Table S4. Results of Kolmogorov-Smirnov and Cramer von Mises analysis for testing the differences between the position distribution of under-expressed genes and coding genes spotted on the microarray chip. (DOCX) File S2 Contains the following files: Table S5. Differential expressed pathogen recognition gene (NB-ARC-LRR, RLP, RLK) in tomato-Fol interaction; Table S6. Differential expressed pathogen recognition gene (NB-ARC-LRR, RLP, RLK) in tomato-ToMV interaction. (XLSX)

Author Contributions
Conceived and designed the experiments: MRE. Performed the experiments: FF AF LS. Analyzed the data: GA LT. Contributed reagents/ materials/analysis tools: GA AF. Wrote the paper: GA MRE LF.