Genome-Wide Identification and Validation of Reference Genes in Infected Tomato Leaves for Quantitative RT-PCR Analyses

The Gram-negative bacterium Xanthomonas campestris pv. vesicatoria (Xcv) causes bacterial spot disease of pepper and tomato by direct translocation of type III effector proteins into the plant cell cytosol. Once in the plant cell the effectors interfere with host cell processes and manipulate the plant transcriptome. Quantitative RT-PCR (qRT-PCR) is usually the method of choice to analyze transcriptional changes of selected plant genes. Reliable results depend, however, on measuring stably expressed reference genes that serve as internal normalization controls. We identified the most stably expressed tomato genes based on microarray analyses of Xcv-infected tomato leaves and evaluated the reliability of 11 genes for qRT-PCR studies in comparison to four traditionally employed reference genes. Three different statistical algorithms, geNorm, NormFinder and BestKeeper, concordantly determined the superiority of the newly identified reference genes. The most suitable reference genes encode proteins with homology to PHD finger family proteins and the U6 snRNA-associated protein LSm7. In addition, we identified pepper orthologs and validated several genes as reliable normalization controls for qRT-PCR analysis of Xcv-infected pepper plants. The newly identified reference genes will be beneficial for future qRT-PCR studies of the Xcv-tomato and Xcv-pepper pathosystems, as well as for the identification of suitable normalization controls for qRT-PCR studies of other plant-pathogen interactions, especially, if related plant species are used in combination with bacterial pathogens.


Introduction
The analysis of gene transcription profiles is a powerful tool to uncover the roles of specific genes in cellular processes and to place them into regulatory networks. Quantitative reverse transcription PCR (qRT-PCR), also termed real-time RT-PCR, is the method of choice to analyze changes in gene transcription because of its high sensitivity, large dynamic range and accuracy [1]. The reliability of results strongly depends on suitable reference genes for novel tomato reference genes. In addition, we identified the pepper orthologs of these genes and found several to be suitable normalization controls for qRT-PCR analyses in pepper during biotic stress.

Plant material and inoculations
Tomato (Solanum lycopersicum) plants of cultivar (cv.) MoneyMaker and pepper (Capsicum annuum) cv. ECW-30R plants were grown in the greenhouse under standard conditions (day and night temperatures of 23°C and 19°C, respectively, for tomato, and 25°C and 19°C for pepper, with 16 h light and 40 to 60% humidity). For qRT-PCR studies, tomato and pepper plants were transferred to a Percival growth chamber (Percival Scientific, Perry, USA) three days before inoculation. Mature leaves of seven-week-old tomato and pepper plants were inoculated with mock (10 mM MgCl 2 ) or Xcv (5×10 8 cfu/ml in 10 mM MgCl 2 ) using a needleless syringe.

Microarray analyses
For microarray studies, 12 tomato plants were inoculated per experiment. To minimize differences in gene expression due to leaf-to-leaf variability, Xcv strains and 10 mM MgCl 2 , respectively, were infiltrated into the same leaves. Four leaf discs (0.5 cm diameter) per inoculum and leaf were harvested, immediately frozen in liquid nitrogen and stored at -80°C. In the first study, Xcv 85-10 and 85-10ΔhrcN were inoculated; leaf material was harvested 45 min and 6, 10 and 24 hours post infiltration (hpi). Leaf material of four plants was pooled for each timepoint (16 leaf discs per sample, three technical replicates). In the second study, 85-10ΔhrcN and 10 mM MgCl 2 were infiltrated and leaf material was harvested at 0, 4, 8 and 16 hpi and pooled as above. In addition, four leaf discs per plant were harvested as control before treatment. This was performed three times independently with four plants each (biological replicates). The experimental setup is summarized in S1 Fig. Total RNA was extracted using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Hilden, Germany) and treated with DNase I (Roche, Mannheim, Germany) for 30 min. Approximately 1.5 μg total RNA was sent to Source BioScience (Berlin, Germany) for cDNA synthesis and microarray hybridizations.
For the tomato whole-genome chip (Source BioScience), oligonucleotides for 34,383 annotated tomato genes [according to the international tomato annotation group (ITAG, version 2. 3)] were spotted on Agilent custom arrays. Five 50-bp oligonucleotides per gene were tested on an Agilent custom array 4x180K, and a set of suitable oligos was chosen for the final chip. Due to space limitations (8x60K), 25,985 randomly chosen genes were represented twice with different oligonucleotides, whereas 8398 genes were represented by one oligonucleotide each. Finally, seven identical 8x60K chips were used for sample analysis. Different chips were hybridized with biological and technical replicates, respectively. cDNA synthesis, labelling, hybridization, washing, scanning and data collection was performed by Source BioScience according to Agilent standard protocols.

Data processing and statistical analyses
Microarray raw data (column "gProcessedSignal") were analyzed by the statistical software R [56]. All experiments of one study (treatments, time points and replicates) were normalized by quantile normalization on the probe level using the "preprocessCore" R package (version 1.26.1, http://www.bioconductor.org/packages/release/bioc/html/preprocessCore.html). For each gene, values for transcript accumulation were obtained as the arithmetic mean of the intensities of all probes representing the gene. The coefficient of variation (CV) was computed for each gene as the standard deviation of its transcript levels across all experiments divided by its mean transcript level. To evaluate the similarity of expression patterns in biological and technical replicates, normalized log-expression values of the individual experiments were clustered hierarchically using the R function hclust [56]. The distance between the expression vectors of experiments was determined as one minus the Pearson correlation of log-expression values using the R function cor.dist from the bioDist package of the Bioconductor suite [57]. Clustering was performed using complete linkage, which yields compact clusters with high intra-cluster correlations. Dendrograms were plotted using the specific plot function of the R class hclust [56].

Quantitative reverse transcription polymerase chain reaction (qRT-PCR)
Templates for qRT-PCR were produced as follows: three to four leaf discs (1.3 cm diameter) from different plants infiltrated with Xcv and MgCl 2 , respectively, were pooled for RNA isolation using the QIAGEN RNeasy Plant Mini Kit. Oligo-dT-and random hexamer-primed cDNA was synthesized with the Maxima First Strand cDNA Synthesis Kit (Thermo Scientific, Schwerte, Germany). qRT-PCR was performed on a CFX96 thermal cycler (Bio-Rad, Munich, Germany) using a SYBR Green-based PCR reaction mixture (Absolute Blue qPCR SYBR Green Fluorescein Mix; Thermo Scientific) and 8 ng template cDNA. Oligonucleotide sequences are listed in S1 Table. To compare Ct (cycle threshold) values measured on different plates using different reaction mixtures, automatically calculated thresholds of all plates were set manually to the highest threshold obtained. The efficiency of PCR reactions was determined for each primer pair using a dilution series of template plotted into a standard curve. To ensure amplification specificity, amplicons were subjected to melting curve analysis and analyzed on 1% agarose gels. Transcript levels were determined as technical duplicates of biological triplicates.

Selection of candidate reference genes for gene expression studies in tomato
To identify reference genes suitable for the analysis of Xcv-induced changes in the mRNA levels of tomato genes we evaluated the results of two whole-genome microarray screens. For the first screen, S. lycopersicum cv. MM plants were inoculated with the Xcv wild-type (WT) strain 85-10 and the T3S-deficient derivative 85-10ΔhrcN, respectively. Leaf material was harvested 45 min and 6, 10 and 24 hours post infiltration (hpi). In the second screen, S. lycopersicum cv. MM plants challenged by 85-10ΔhrcN inoculation were compared to mock-infiltrated tomato plants, and leaf material was harvested at 0, 4, 8 and 16 hpi. Transcriptional changes of 34,383 annotated tomato genes were analyzed using "Agilent custom arrays". Hierarchical cluster analysis illustrates similar expression patterns in biological and technical replicates confirming that the experimental treatments worked (S1 Fig). In the first screen, two samples ("85-10; 45 mpi; #2" and "85-10; 6 hpi; #2") showed aberrant gene expression patterns compared to the corresponding replicates resulting in separate clustering (S1A Fig). Both samples were excluded from further data evaluation.
The microarray analyses revealed a high variability in the expression patterns of housekeeping genes conventionally used as references in transcript studies [30,[61][62][63][64] (Fig 1). To identify the most stably transcribed genes, the coefficient of variation (CV) was determined for each gene, which is defined as the standard deviation of its expression levels across all experiments (treatments and time-points) divided by its mean expression level. Genes with a log 2 mean expression level below 7 or above 13 were excluded to account for the bigger influence of random noise on low expression values, and for saturation effects of microarrays at high mRNA levels, respectively. Genes with CV values 0.12 in both microarray studies were ranked by increasing CV in the second screen which delivers more reliable data compared to the first study (biological instead of technical replicates). The best 50 candidate reference genes are listed in Table 1. The tomato sequences were classified based on the functional categories of their A. thaliana orthologs which were identified by BLASTx [65] against "The Arabidopsis Information Resource" database (TAIR Blast 2.2.8; S2 Fig). Only predicted proteins that displayed minimum 40% amino-acid identity over at least 70% of the tomato sequence were taken into account. This allowed a functional classification of approximately three quarters of the sequences (74%), most of them possessing putative functions in protein expression (transcription and splicing) and turnover (ubiquitination/proteolysis; S2 Fig).
Evaluation of the expression stability of novel and traditional tomato reference genes qRT-PCR analyses of the 11 most stably expressed genes (Table 1) were performed to validate their expression stability in S. lycopersicum cv. MM infected with Xcv. The genes encode a TFIID subunit (TAF6), importin β (IMP-β), a PHD finger family protein (PHD), a cytochrome c oxidase subunit (COX), polyribonucleotide 5´-hydroxyl-kinase Clp1 (CLP1), a ubiquitin carboxyl-terminal hydrolase family protein (UCH), a polypyrimidine tract-binding protein-like protein (PTBL), U6 snRNA-associated Sm-like protein LSm7 (LSM7) and an acyl carrier protein (ACP), as well as two unknown proteins (UP1 and UP2; Table 1, S2 Table). For comparison, four housekeeping genes were analyzed that are widely used as references, namely actin (ACT), EF-1α, GAPDH and ubiquitin (UBI). First, suitability of oligonucleotides (S1 Table) and target sequences was confirmed. Melting curves and gel electrophoresis revealed unique amplicons for all oligonucleotide combinations used validating their specificity (S3 Fig). PCR efficiencies ranged between 80.48 and 99.71% (S1 Table). For expression analysis of the reference gene candidates, total RNA was isolated from tomato leaves 0, 6, 10 and 24 h after treatment with 10 mM MgCl 2 , Xcv 85-10, 85-10ΔhrcN and 85-10(pavrBs4), respectively. The latter strain induces the ETI, i.e., the HR in S. lycopersicum cv. MM due to the Bs4-dependent recognition of the avirulence protein AvrBs4, a member of the TAL effector family [66]. Technical duplicates of three biological replicates were subjected to qRT-PCR analysis. Average Ct (cycle threshold) values of the new reference gene candidates ranged from 27.1 (CLP1) to 31.1 (UP1; Fig 2). To select the optimal reference genes, we used three different algorithms to evaluate our qRT-PCR results: geNorm [58], NormFinder [59] and BestKeeper [60]. geNorm analysis. The geNorm software provides a ranking of the tested genes based on a stability value M which is calculated by average pairwise variation of each candidate gene combination [58]. The lower the M value, the higher the expression stability of the gene. Eventually, the algorithm selects an optimal pair of reference genes out of the candidate set analyzed. Considering a cutoff of M 0.5, the traditional references GAPDH, ACT and UBI proved unreliable for the normalization of qRT-PCR data under the experimental conditions chosen (Fig 3). By contrast, all newly identified tomato candidate genes and EF-1α represent suitable references, with IMP-β and PHD being optimal (Fig 3).
NormFinder analysis. Next, we analyzed the qRT-PCR data using NormFinder [59]. The stability value M calculated by this "model-based variance estimation approach" considers not only the "overall expression variation" measured in different samples, but additionally takes into account variations among and inside sample subgroups [59]. Thus, the algorithm avoids co-regulated reference genes which display systematic intergroup variation and would lead to erroneous conclusions. Since we are interested in changes of plant gene expression levels induced by different Xanthomonas strains but also in expression level changes over a certain time period, two separate NormFinder analyses were performed with sample subgroups defined based on treatment [MgCl 2 , Xcv 85-10, 85-10ΔhrcN and 85-10(pavrBs4)] and time- point of sampling (0, 6, 10 and 24 hpi), respectively. As shown in Fig 4, all tested genes fulfill the minimal requirement for suitable reference genes, i.e., possess an M value below 1.5. However, the traditionally employed reference genes ACT and GAPDH were considerably less stable than the other genes, whereas EF-1α and UBI seemed more suitable under the chosen experimental conditions. The top-ranked references, however, were among the newly identified candidate genes, namely COX > PHD > CLP1 > LSM7 with respect to the grouping by treatment (Fig 4A).
BestKeeper analysis. We compared the six most stable new reference genes according to NormFinder with the four classical reference genes using BestKeeper [60]. This tool evaluates the suitability of up to 10 reference genes based on the calculation of Ct value variations, performing pair-wise correlations of all candidate gene combinations. Extreme samples (x-fold over-/under-expression) are also considered. As shown in Table 2, expression of all genes except for ACT and GAPDH fluctuated in a range compatible with standard deviations (SD) [± Ct] < 1 and SD [± x-fold] < 2, which represents an acceptable overall variation [60]. Notably, BestKeeper evaluated all six new reference gene candidates as better suited than the four traditional housekeeping genes, with PHD > CLP1 > LSM7 > COX being the top four. Taken together, regardless of the ranking order, geNorm, NormFinder and BestKeeper evidenced the superior expression stability of the new tomato reference genes under the experimental conditions chosen.

Quantification of immunity marker genes in infected tomato leaves
We applied our findings to the analysis of two target genes previously reported to be induced during PTI and ETI, respectively, LRR22 [67] and an UDP-glucosyltransferase gene (UGT, Solyc09g092500 [68]). For this, total RNA was analyzed from tomato leaves six hours after treatment with 10 mM MgCl 2 , Xcv 85-10, 85-10ΔhrcN and 85-10(pavrBs4), respectively. To increase the accuracy of normalization we took into account two reference genes. We compared the two best reference genes identified by geNorm (IMP-β and PHD), NormFinder (PHD and COX) and BestKeeper (CLP1 and PHD) with the two least-stable genes, GAPDH and ACT, for their ability to provide reliable relative quantification of SlLRR22 and SlUGT by qRT-PCR. As shown in Fig 5, accumulation of SlLRR22 transcript was approximately two-fold higher in the leaves treated with 85-10ΔhrcN than in the mock control if compared to any of the new reference gene combinations. By contrast, comparison to the suboptimal references revealed an apparent five-fold induction of gene expression. In addition, referring to ACT and GAPDH suggested a more than two-fold upregulation of SlLRR22 by the Xcv WT strain 85-10 and by 85-10(pavrBs4), the latter induction being significant, which was not detectable with any of the superior reference genes. Notably, standard deviations between the different biological datasets were substantially lower if one of the new reference gene combinations was employed. The analysis of the ETI marker gene, SlUGT, did not show pronounced differences in the expression pattern depending on the reference genes chosen. In all cases, transcript abundance was significantly higher in the leaves treated with the avirulent strain 85-10 (pavrBs4) than in the mock-infiltrated leaves. However, a slight induction of SlUGT expression by both Xcv 85-10 and 85-10ΔhrcN was only detected when the traditional references were employed. A possible explanation for these results is downregulation of ACT and/or GAPDH by Xcv infection. To test this possibility, the expression of both genes was analyzed using the newly identified reference genes as normalization controls. As shown in S4 Fig, GAPDH transcript levels were indeed significantly lower in the leaf material inoculated with bacteria compared to the mock control, whereas ACT appeared not to be changed under these conditions.

Selection and validation of pepper reference genes based on tomato orthologs
Based on the tomato microarray data, pepper orthologs of the eleven most stably expressed genes (Table 1) were identified by BLASTx against the European Nucleotide Archive (http://www.ebi. ac.uk/ena). Oligonucleotides for qRT-PCR were derived (S1 Table), and melting curve analysis and gel electrophoresis confirmed specific products for nine candidate genes (S5 Fig). PCR efficiencies ranged between 72.09 and 99.32% (S1 Table). For expression analysis, pepper ECW-30R (Bs3) leaves were infiltrated with 10 mM MgCl 2 , Xcv 85-10, 85-10ΔhrcN and 85-10(pavrBs3), respectively, and leaf material was harvested at 0, 6, 10 and 24 hpi. 85-10(pavrBs3) translocates the effector AvrBs3 which induces the HR in Bs3 pepper plants. Technical duplicates of three biological replicates were subjected to qRT-PCR analysis. Average Ct values of the new reference gene candidates ranged from 27.4 (UCH) to 38.8 (TAF6; Fig 6). For comparison, the four classical reference genes EF-1α, GAPDH, ACT and β-tubulin (TUB) were also analyzed. The data were evaluated similarly to the analysis of the tomato reference genes described above. GeNorm analysis revealed that only three genes, UCH, LSM7 and PHD, match the cutoff-value for a reliable reference gene (M 0.5). In general, the pepper orthologs of the newly identified tomato reference genes were more stably expressed than the traditional pepper references (Fig 7).
Using NormFinder, the classical reference gene EF-1α matched the requirements of a suitable reference gene (M < 1.5) when the sample subgroups were defined by treatment (Fig 8A), but turned out to be completely unreliable when the classification was based on the time-point of sampling (Fig 8B). GAPDH and TUB matched the minimal requirements of a reliable reference gene but were considerably less stable than the other genes tested, while ACT appeared more suitable. Notably, all newly identified reference genes were evaluated as reliable normalization controls with UCH > PHD > UP2 > LSM7 as the top-four when the grouping was based on treatment ( Fig 8A). BestKeeper analysis of the six best new pepper reference genes according to NormFinder and the four classical references surprisingly revealed that only one gene, UCH, fulfilled both requirements for a suitable normalization control in qRT-PCR studies, i.e., SD [± Ct] < 1 and SD [± x-fold] < 2 ( Table 3). Most of the other genes matched at least the threshold for SD [± x-fold], whereas EF-1α appeared to be completely unreliable as reference gene (Table 3). SlUGT (b) in S. lycopersicum cv. MM leaves treated with 10 mM MgCl 2 (mock) or 5×10 8 cfu/ml of Xcv 85-10, 85-10ΔhrcN and 85-10(pavrBs4), respectively, 6 hpi. qRT-PCR data were normalized with different reference gene pairs. Values are mean-fold changes in mRNA levels in Xcv-infected relative to mockinoculated leaves for three biological replicates. Error bars indicate standard deviation (SD). Letters denote statistically significant differences (Student´s ttest, P < 0.05). doi:10.1371/journal.pone.0136499.g005

Quantification of PTI and ETI marker genes in infected pepper leaves
We compared the two best reference genes from pepper identified by geNorm, i.e., UCH and LSM7 with the best reference genes according to BestKeeper, UCH and ACT, and the traditional reference genes EF-1α and GAPDH for their ability to provide reliable relative quantification of the target genes LRR22 and TFT4, which are induced during PTI and ETI, respectively [67,69]. As shown in Fig 9, employment of different reference genes did not result in substantial differences in the expression patterns of CaLRR22 and CaTFT4. The T3S-deficient Xcv strain 85-10ΔhrcN led to significantly higher expression of CaLRR22 compared to the WT strain 85-10 and Xcv 85-10(pavrBs3). CaTFT4 was induced significantly during the incompatible interaction with Xcv 85-10(pavrBs3), similarly to the reported induction after recognition of the type III effector AvrBs2 [69]. However, the observed differences in target gene expression levels were only judged as significant when the newly identified reference genes were used, but not with the traditional combination EF-1α/GAPDH. Utilization of the newly identified normalization controls resulted in significantly lower standard deviations underlining the higher reproducibility of the results in different experiments (Fig 9).

Discussion
Correct normalization of gene transcripts depends on the choice of suitable reference genes. This is essential for reliable analyses of gene expression by qRT-PCR and has to be established for specific experimental conditions [4]. Based on microarray expression analyses of >34,000 genes, we identified and validated 11 novel tomato reference genes with superior expression stability under biotic stress conditions, i.e., challenge by the bacterial pathogen Xcv. Although Novel Tomato Reference Genes the new reference genes do not comprise "classical" housekeeping genes, homologies on the protein level indicate putative roles in basic cell functions, e.g., oxidation/reduction processes (COX), mRNA processing (LSM7, CLP1, PTBL), regulation of transcription/chromatin dynamics (PHD), nuclear import (IMP-β) and fatty acid biosynthesis (ACP; S2 Table). The three statistical programs we used for the evaluation of gene expression stability, geNorm, NormFinder and BestKeeper, slightly differed in the ranking of the reference gene candidates, which was also observed in previous studies and is probably due to different algorithms underlying the programs [23,30,31,70]. Importantly, the newly identified genes were usually evaluated as more stable than the traditional housekeeping genes we analyzed for comparison and, notably, always included the optimal normalization control identified by the respective program. Based on our results, we recommend the use of PHD and LSM7 as reference genes for normalization in future plant gene expression studies in the Xcv-tomato pathosystem.
To the best of our knowledge, previous studies of pepper and tomato comparing reference gene stabilities selected candidates solely based on homology. It was shown that different genes, often housekeeping genes, are preferable under different conditions [23,30,38,62,[71][72][73]. Notably, our microarray data revealed that the expression of classical tomato housekeeping genes varied considerably, confirmed by qRT-PCR studies of selected genes. In particular, GAPDH and ACT were attested a variability too high for a reliable reference gene by geNorm and BestKeeper, respectively. Therefore, we do not recommend the further employment of these genes as normalization controls in qRT-PCR analysis of tomato genes after pathogen infection, especially because we clearly showed an Xcv-dependent downregulation of GAPDH expression. Taken together, our results demonstrate the advantage of an unbiased, whole transcriptome-based approach to identify suitable reference genes. Concordantly, several whole-transcriptome analyses of different plant species and experimental setups identified other than traditional housekeeping genes as the most stably expressed genes [37][38][39][40][41][42][43].
It is, however, not feasible to perform microarray analyses for reference gene identification every time the experimental setup is changed. Therefore, one has to resort also to the homology-based selection of candidate genes. The identification of suitable candidates can be strongly improved by using orthologs of genes that were experimentally verified as appropriate Novel Tomato Reference Genes references in related organisms under similar experimental or developmental conditions [38,[74][75][76]. We used such an approach to identify the pepper orthologs of our new superior tomato reference genes and determined UCH and PHD as the most suitable references for normalization of plant gene expression in the Xcv-pepper pathosystem. Interestingly, one of the traditional reference genes, ACT, also turned out to be stably expressed in our experimental setup. This contradicts the results of Wan et al. who described ACT as relatively unstable under different abiotic stresses and hormonal treatments [72]. On the other hand, EF-1α turned out to be the most unstable pepper gene in our analyses although it was published as one of the least-variably expressed genes under abiotic stress conditions and hormone treatment [71]. This underpins the observation that a chosen gene can be stable under certain conditions but highly variable under others [3]. It should be noted that differences between the pepper lines used in the different studies might also play a role.
Although our selection of pepper orthologs of the new tomato reference genes surely represents an improvement compared to the selection of genes based on their known or suspected housekeeping roles, the ranking of our tomato reference genes and their pepper equivalents illustrates that the expression of gene orthologs can distinctly differ even between related plant species. In general, the M values calculated by NormFinder were lower for the tomato genes compared with their pepper orthologs. This difference appeared even more pronounced using geNorm which judged only three of the pepper genes tested as reliable reference genes. Similarly, using Bestkeeper, only one pepper gene, UCH, matched both requirements for a suitable reference gene. Therefore, we would like to emphasize that, even if our new pepper reference genes proved to be superior to most of the classical normalization controls we analyzed, a whole-transcriptome analysis of Xcv-challenged pepper plants might uncover even more suitable reference genes.
Taken together, the newly discovered tomato reference genes proved to be superior normalization controls for qRT-PCR studies of Xcv-infected tomato plants. In addition, they led to successful identification of the pepper orthologs as reliable reference genes in qRT-PCR analyses of the Xcv-pepper pathosystem. Similarly, these genes might be useful for the identification of suitable qRT-PCR normalization controls in other plant species for the analysis of plant gene expression during pathogen infection. Supporting Information Relative expression of PTI and ETI marker genes in Xcv-infected and mock-treated pepper leaves. Expression patterns of CaLRR22 (a) and CaTFT4 (b) in C. annuum ECW-30R leaves treated with 10 mM MgCl 2 (mock) or 5×10 8 cfu/ml of Xcv 85-10, 85-10ΔhrcN and 85-10(pavrBs3), respectively, six hpi. qRT-PCR data were normalized with different reference gene pairs. Values are mean fold changes in mRNA levels in Xcv-infected relative to mockinoculated leaves for three biological replicates. Error bars indicate SD. Letters denote statistically significant differences (Student´s t-test, P < 0.05).