Transcriptomic analysis reveals tomato genes whose expression is induced specifically during effector-triggered immunity and identifies the Epk1 protein kinase which is required for the host response to three bacterial effector proteins
© Pombo et al.; licensee BioMed Central Ltd. 2014
Received: 2 July 2014
Accepted: 9 October 2014
Published: 17 October 2014
Plants have two related immune systems to defend themselves against pathogen attack. Initially, pattern-triggered immunity is activated upon recognition of microbe-associated molecular patterns by pattern recognition receptors. Pathogenic bacteria deliver effector proteins into the plant cell that interfere with this immune response and promote disease. However, some plants express resistance proteins that detect the presence of specific effectors leading to a robust defense response referred to as effector-triggered immunity. The interaction of tomato with Pseudomonas syringae pv. tomato is an established model system for understanding the molecular basis of these plant immune responses.
We apply high-throughput RNA sequencing to this pathosystem to identify genes whose expression changes specifically during pattern-triggered or effector-triggered immunity. We then develop reporter genes for each of these responses that will enable characterization of the host response to the large collection of P. s. pv. tomato strains that express different combinations of effectors. Virus-induced gene silencing of 30 of the effector-triggered immunity-specific genes identifies Epk1 which encodes a predicted protein kinase from a family previously unknown to be involved in immunity. Knocked-down expression of Epk1 compromises effector-triggered immunity triggered by three bacterial effectors but not by effectors from non-bacterial pathogens. Epistasis experiments indicate that Epk1 acts upstream of effector-triggered immunity-associated MAP kinase signaling.
Using RNA-seq technology we identify genes involved in specific immune responses. A functional genomics screen led to the discovery of Epk1, a novel predicted protein kinase required for plant defense activation upon recognition of three different bacterial effectors.
Plants are vulnerable to attack by many pathogenic microorganisms. To respond to these assaults, plants have evolved two interlinked layers of immunity. Plants initially use pattern recognition receptors to recognize microorganisms by detecting certain conserved features referred to as microbe- or pathogen-associated molecular patterns (MAMPs or PAMPs) ,. Such pattern-triggered immunity (PTI) leads to production of reactive oxygen species, activation of mitogen-activated protein kinase (MAPK) cascades, changes in the intracellular calcium concentration and transcriptional reprogramming -. However, pathogens such as Pseudomonas syringae undermine PTI by delivering virulence proteins (effectors) into the plant cell using a type III secretion system . In a further evolutionary step, some plants acquired intracellular proteins that detect, either directly or indirectly, the presence of specific effectors. This layer of defense, termed effector-triggered immunity (ETI), is often associated with localized programmed cell death (PCD) called the hypersensitive response that may limit pathogen spread ,,.
The interaction between tomato (Solanum lycopersicum) and Pseudomonas syringae pv. tomato (Pst), the causative agent of bacterial speck disease, has been extensively used to study the molecular basis of host responses to bacterial infection ,. Among the bacterial MAMPs perceived by tomato, the best characterized derive from the flagellin protein encoded by the fliC gene. This protein, which forms the flagellum and therefore plays a key role in motility, possesses two MAMPs that can be detected by tomato: flg22 perceived by the FLS2 receptor ,; and flgII-28, which is recognized by an unidentified receptor referred to as FLS3 . We have recently reported that flagellin-derived MAMPs in Pst are the primary elicitors of PTI in tomato, resulting in extensive transcriptional changes . Pst strain DC3000 translocates approximately 30 effectors into plant cells and two of these, AvrPto and AvrPtoB, act early in the tomato-Pst interaction by interfering with pattern recognition receptor functions and thereby suppressing PTI and promoting bacterial virulence -.
Certain wild relatives of tomato have evolved a specific ETI mechanism to recognize and respond to the presence of AvrPto or AvrPtoB in the plant cell. This mechanism involves members of the Pto kinase family which physically interact with these effectors and act with the nucleotide binding-leucine-rich repeat (NB-LRR) protein Prf to activate ETI ,,. Prf is embedded within the Pto family gene cluster on chromosome 5 and this region has been introgressed into many tomato cultivars to confer resistance to bacterial speck disease ,. Changes in gene expression that occur during Pto/Prf-mediated ETI in response to AvrPto were previously analyzed using GeneCalling, an mRNA profiling technology . This study was limited by the lack of a tomato genome sequence and gene annotation, but nevertheless identified 432 ETI-induced genes, including members of 11 transcription factor gene families. The experimental design used in this study did not allow the determination of whether any of these genes were also induced during PTI.
Several recent reports have used microarrays and a series of Arabidopsis mutants to analyze changes that occur during ETI and PTI ,. Although PTI involves recognition of MAMPs and ETI recognition of effectors, these studies found that a majority of genes whose expression changes during the defense response are affected by both ETI and PTI -. Analysis of transcriptome changes also indicated that the PTI response was transient and more vulnerable to being undermined by the pathogen, whereas the ETI response was more prolonged and robust . These differences appeared to be due to how ETI and PTI use the same signaling networks rather than to their use of distinct signaling networks . In the case of PTI, gene expression changes were characterized as being synergistic and for this reason more vulnerable to pathogen suppression. ETI, however, uses the same signaling network in a compensatory way, making it difficult for the pathogen to interfere with this response . Overall, these studies have provided important insights into the dynamics and temporal aspects of plant immunity.
RNA-Seq is a powerful high-throughput technology that is being broadly used for transcriptome studies in different cells and treatments -. The RNA-Seq approach has recently been applied to the analysis of PTI in the tomato-Pst system, where it allowed the identification of a subset of genes whose expression is induced by flagellin-derived MAMPs but reduced by the activity of AvrPto and AvrPtoB effectors, referred to as FIRE (flagellin-induced, repressed by effectors) genes . These FIRE genes were screened by a virus-induced gene silencing (VIGS) approach which identified a cell wall-associated kinase, SlWAK1, that is required for an effective PTI response .
Here we describe the use of the tomato-Pst system in an experimental design that allowed the identification of host genes whose expression changes specifically in ETI or PTI. We developed reporter genes for each of these responses which will be useful for analyzing the response of tomato to Pst mutants that lack combinations of effectors, other virulence factors, or MAMPs. We then focused on a subset of 30 genes whose expression was induced specifically during ETI and used these in a VIGS screen to determine if they play a demonstrable role in ETI against Pst. This approach identified a predicted protein kinase, Epk1, that has not previously been implicated in the plant immune response and which may act in a pathway unique to the immune response triggered by bacterial effectors.
Analysis of transcriptome modifications during Pto/Prf-mediated effector-triggered immunity in tomato
In order to study the transcriptome changes in tomato during Pto/Prf-mediated ETI, we infiltrated tomato Rio Grande (RG)-PtoR resistant plants (plants that have a functional Pto/Prf signaling pathway, Pto/Pto, Prf/Prf) and two different susceptible plants: RG-prf3 and RG-prf19 (Pto/Pto, prf/prf), with Pseudomonas syringae pv. tomato DC3000 (DC3000) (Figure S1A in Additional file 1). The susceptible lines have a non-functional Prf gene due to a 1.1 kb deletion or a G-insertion at position 2,584 (which causes a frameshift), respectively . We collected leaf tissue at 4 and 6 h after inoculation (hai) to assess early changes in host gene expression after translocation of DC3000 effectors AvrPto and AvrPtoB, which occur at about 3 hai. The plants were then maintained in the same conditions to observe signs of disease. As expected, RG-PtoR plants did not develop speck disease whereas the RG-prf plants did (Figure S1B in Additional file 1).
For the data analysis, we considered a gene to be `expressed’ if it had three or more RPKMs (reads per kilobase per million of mapped reads) in at least one of the treatments analyzed. The cutoff used for the comparisons was P <0.05 and ≥2-fold expression change . We compared the transcriptome changes observed in resistant (RG-PtoR) and susceptible (RG-prf3 or RG-prf19) plants and found that the number of differentially expressed genes increased from 4 to 6 h (Figure S1C in Additional file 1). We considered a gene to be induced during ETI when its expression was higher in RG-PtoR than in RG-prf3 or RG-prf19 plants. On the other hand, a gene was considered to have reduced expression during ETI when it was expressed higher in RG-prf3 or RG-prf19 than in RG-PtoR plants after DC3000 inoculation. The number of genes with reduced expression was higher than those with induced expression in both PtoR/prf3 and PtoR/prf19 combinations. As expected, the overall gene expression differences of the two susceptible lines were very small, with only 27 differentially expressed genes between RG-prf3 and RG-prf19.
Taking just the ETI-induced genes from this experiment at 6 h, we compared them with PTI-induced genes reported previously . For this purpose we considered as PTI-induced genes those increased by flgII-28, DC3000 ΔhrcQ-U, Pseudomonas fluorescens or the Pseudomonas putida treatments . This comparison revealed that essentially the same number of genes was induced only in ETI, only in PTI, or in both responses (Figure S1D in Additional file 1). Thus, although there was overlap in the ETI- and PTI-associated transcriptomes, there also appeared to be unique gene expression changes associated with each of these immune responses.
Transcriptome changes associated with Pto/Prf-mediated ETI and flagellin-activated PTI revealed a set of genes specifically induced in each immune response
Summary of the experiments performed for the RNA-Seq analysis
5 × 106Δcfu/ml
DC3000 ΔavrPto ΔavrPtoBc
DC3000 ΔavrPto ΔavrPtoB ΔfliCd
We performed a bibliographic search to identify genes in the three categories that have been described previously as being implicated in plant-pathogen interactions (Additional file 3). Approximately 25% of the genes had at least one publication associating them with plant immunity (Additional file 3). In the ETI-specific category, we found 10 genes that had been previously confirmed to be involved in ETI mediated by the Pto kinase using a loss of function approach .
Gene Ontology term analysis supports the central importance of protein kinase activity in plant immunity
We performed a Gene Ontology (GO) term analysis using the tomato genome sequence as a reference (Figure 2A; Additional file 5). The `Process’ GO terms of `defense’ or `response to stimulus’ or `stress’ were the most common for the PTI-specific and the shared PTI/ETI genes, whereas `transport’ and `cellular localization’ had the highest enrichment for the ETI-specific genes. Additional Process GO terms associated with ETI were `response to stress’, `response to stimulus’, `signal transduction’ and `regulation of cell death’ (Additional file 5). In the `Function’ category of GO terms all three groups of genes showed similar results with `protein kinase activity’ and `lipid-binding’, `carbohydrate-binding’, or `protein-binding` being most prominent. Among the ETI-specific genes, 176 of them are kinases representing 27% of the total expressed kinases in our analysis (657 expressed kinases out of 1,150 total predicted in tomato ), supporting the relevance of protein kinases in plant immunity.
We categorized the protein kinase-encoding genes whose expression was induced by our treatments (Figure 2B). Of the 64 protein kinase families present in tomato , 46 (72%) had a member that was induced in at least one of the categories (ETI-specific, both, or PTI-specific). PTI-specific kinases occurred in just five families and all of them also contained ETI-specific kinases. A total of 20 kinase families have genes induced during both PTI and ETI. The remaining 26 are families that possess only ETI-specific induced kinases and many of these have low numbers of induced kinase genes (from 1 to 5).
We also categorized induced transcription factor genes and found that a total of 249 occurred in all three groups, with 190 being induced specifically during ETI (Figure S4A in Additional file 6). Analysis of the transcription factor families indicated that those with the largest number of induced genes (AP2-EREBP, C2H2, MYB, NAC and WRKY) contain genes from all the three categories (ETI-specific, common genes and PTI-specific). However, most of the families have only ETI-specific induced transcription factors (Figure S4B in Additional file 6).
Development of Pto/Prf-mediated ETI and flagellin-activated PTI reporter genes
RG-PtoR plants were again infiltrated with different DC3000 strains to assess the response of the reporter genes to ETI and PTI, but in this case we took leaf samples at 3, 6, 9 and 12 hai to determine the temporal pattern of their induction. By using qRT-PCR we confirmed each gene to be ETI- or PTI-specific at 6 hai as found from the RNA-Seq data (Figure 3B,C). In the case of the ETI-induced genes, the two UDP-glucosyltransferase genes were shown to be good reporters at 6, 9, and 12 hai. The laccase gene was ETI-specific only at 6 h after treatment (Figure 3B). All three PTI-specific genes were induced at each time point only during this immune response starting at 6 hai (Figure 3C). None of the six genes proved to be effective at distinguishing the two immune responses at 3 hai, which is likely before translocation of the effectors into the plant cell.
Identification of a novel kinase involved in plant immunity using a virus-induced gene silencing screen
To determine if the ETI-specific genes play a demonstrable role in this immune response we selected a subset of 30 genes that encode predicted protein kinases and transcription factors and performed a loss-of-function screen using VIGS in Nicotiana benthamiana (Additional file 3). Leaves of silenced plants were infiltrated with a mix of Agrobacterium tumefaciens strains carrying 35S:Pto or 35S:avrPto. Three genes with known roles in ETI were included as positive controls (SlMAPKKKα, SlMEK2 and SlPrf; Additional file 7) ,,. Negative controls were a fragment of the green fluorescent protein (GFP) gene and an Escherichia coli-derived DNA fragment (Ec1) (Additional file 7).
Silencing effectiveness of the SlEpkvirus-induced gene silencing construct in N. benthamiana
A VIGS construct using a 260 bp fragment of the tomato Epk1 gene was used for silencing N. benthamiana plants (Additional file 7). To analyze the specificity of the silencing, we first identified genes in N. benthamiana that have 100% nucleotide identity to the Epk1 fragment over a stretch of ≥17 nucleotides. Twelve such genes were identified, although six were not considered further due to low or no expression (based on N. benthamiana RNA-Seq data) or because they appeared to derive from a sequence encoding a non-functional protein (Additional file 10). We performed qRT-PCR using Ec1- and SlEpk1-silenced N. benthamiana plants on the six remaining genes (Additional file 10). Of these six genes, transcript abundance for four was reduced in the SlEpk1-silenced plants, with the degree ranging from 20 to 80% (Figure S5C in Additional file 8). It was notable that the more similar the N. benthamiana gene sequence is to SlEpk1, the better it was silenced (Figure S5C in Additional file 8; Additional file 11). All four of the silenced genes are in the same clade with SlEpk1 (Additional file 11, clade A); the two unaffected genes belong to another clade (Additional file 11, clade B). These results indicate that one or more genes in N. benthamiana are Epk1-related genes and one or more of them could contribute to ETI.
Silencing of SlEpkcompromises resistance to Pseudomonas syringae pv. tabaci (AvrPto)
SlEpk1 also plays a role in effector-triggered immunity activated by the bacterial effector HopQ1-1
SlEpkdoes not contribute to ETI activated by several non-bacterial effectors
SlEpkappears to act upstream of effector-triggered immunity-associated MAPK signaling
To investigate the position at which SlEpk1 might act in the signaling pathway activated upon AvrPto recognition by Pto, we performed an epistasis experiment that involved silencing and transiently expressing genes known to be involved in ETI . N. benthamiana plants were silenced for SlEpk1, MAPKKKα, MEK2 or SlPrf, with Ec1 as a negative control. MAPKKKα and MEK2DD, both of which cause ETI-associated PCD, were then transiently expressed under an estradiol inducible promoter in the silenced leaves. If a silenced gene acts downstream or directly with MAPKKKα or MEK2, we expect to see compromised PCD in these experiments.
We discovered in two independent experiments that there are plant genes whose expression is induced specifically during ETI or PTI. There are previous reports describing an overlap in transcriptomic changes during ETI/PTI , and we also observed that a majority of the genes whose expression is altered at 6 h after activation of PTI are also altered at the same time point after activation of ETI. However, we identified a high percentage of ETI-induced genes that were not affected during PTI. The significant differences we observed between ETI and PTI could be due to more robust signaling in the case of ETI and because the delivery of effectors into the plant cell suppresses flagellin-activated PTI . Previous studies have not focused on immune response-specific gene expression changes as we have done here, but have instead looked for shared signaling networks between ETI and PTI -. The discovery of ETI- and PTI-specific genes raises the possibility that these distinct classes are enriched for genes that play a role in these specific immune responses and our bibliographic search supported this hypothesis. In agreement with this, our screening of just 30 ETI-specific genes successfully identified a predicted protein kinase from a family previously unknown to have a role in the plant immune response.
Our GO term analysis revealed characteristics about the ETI- and PTI-specific gene sets that distinguish them from each other. For example, the PTI-specific group is enriched in secondary metabolite biosynthetic processes mainly related to the phenylpropanoid pathway known to be important in plant defense. Although it is not exclusive to PTI, membrane receptor kinase activity seems to have a primary role in this form of plant immunity. As expected, only in the ETI-specific gene group is there an enrichment in processes related to the hypersensitive response ,, such as programmed cell death, ion homeostasis, regulation of apoptotic process, response to hydrogen peroxide and other ROS, membrane organization and vesicle-mediated transport. A number of genes described previously as involved in incompatible plant-pathogen interactions were found to be specifically induced during ETI, supporting the idea that gene expression differences are a useful way to distinguish the two immune responses .
Several previous studies have used reporter genes as a readout for activation of PTI and its subsequent suppression by pathogen effectors -. Reporter genes for PTI have also been developed for tomato and N. benthamiana,,. However, these reporters have typically not been shown to be specific for PTI and may also be induced during ETI. In fact, using our RNA-Seq data we found that three tomato PTI reporter genes developed previously  are also induced during ETI. In addition, the tomato genes most closely related to FRK1, a commonly used PTI reporter gene in Arabidopsis, are not expressed or not induced during ETI or PTI, highlighting the importance of plant genera- or species-specific reporter genes. Plant-pathogen interactions can result in host cell death due to ETI or disease that is associated with activation and subsequent suppression of PTI. Reporter genes that distinguish between these two immune responses will be a valuable tool to characterize the plant response to Pst strains that have mutations affecting specific effectors, other virulence factors or MAMPs. Our RNA-Seq data and subsequent qRT-PCR analyses identified six such reporter genes for tomato that can be used as specific plant defense markers in a biologically relevant time frame after bacterial inoculation (from 6 to 12 hai). In the future, development of specific reporter genes for N. benthamiana will be useful for understanding the contribution of SlEpk1 to ETI.
Our loss-of-function approach to screening a set of 30 ETI-specific genes identified a novel predicted protein kinase, SlEpk1. This protein is most closely related to the GmPK6/AtMRK1-like protein kinase family. In Arabidopsis, this family of kinases is in the MAPKKK group of proteins . However, in tomato it is considered as an independent family with around 51 members . In both Arabidopsis and tomato, the genes in this family are known only from genomic analysis without any reports about their role in a specific biological process.
There are many examples of proteins that play a role in ETI against diverse pathogens. For example, RAR1 has been shown to act in an ETI pathway triggered by many coiled-coil NB-LRR proteins and EDS1 acts in pathways triggered by many different TIR-NB-LRR R proteins ,. Similarly, silencing of the genes encoding MAPKKKα, MEK2 or TFT7, a 14-3-3 protein, compromises ETI elicited by R proteins directed at diverse pathogens ,,,. In contrast, our experiments with five different R protein/effector pairs indicate that SlEpk1 may act in a pathway that responds specifically to Pst effectors. If future experiments confirm this specificity, it suggests that Epk1 may act early in these ETI pathways. Consistent with the apparent specificity of SlEpk1 for bacterial ETI, we discovered that this kinase likely acts upstream of MAPKKKα. It is possible that Epk1 participates in a complex directly with Pto/Prf and in a protein complex involved in HopQ1-1 recognition. Alternatively, SlEpk1 might act downstream of these recognition complexes at a juncture shared by bacterial resistance pathways, possibly even interacting with MAPKKKα.
Further experiments are needed to understand the function of SlEpk1 in plant immunity. For example, it will be important to determine if SlEpk1 is an active kinase and if its kinase activity is required for its function in plant defense triggered by bacterial effectors. If SlEpk1 is an active kinase, it will be important to determine its possible substrates. SlEpk1 and its closest tobacco and Arabidopsis genes are described as serine/threonine tyrosine protein kinases ,. In the case of tobacco, DSK1 was confirmed to be phosphorylated on all three of these amino acids . Lately, tyrosine phosphorylation of EF-TU RECEPTOR (EFR) upon MAMP recognition was shown to play an important role in innate immunity . Consequently, if Epk1 is an active kinase, it will be interesting to determine which amino acids in its kinase domain are phosphorylated. An Epk1-like protein has been reported to localize to the chloroplast in tobacco and to the plasma membrane in Arabidopsis,. Determining the subcellular localization of the tomato SlEpk1 protein may also provide insights into the role of SlEpk1 in plant immunity.
Using RNA-Seq technology, we identified genes whose transcript abundance is increased specifically during PTI or ETI. Highlighting the relevance of these sets of induced genes, 25% of them were implicated previously as contributing to plant-pathogen interactions. We used these data to develop specific reporter genes that will be useful for future studies of plant responses to different bacterial pathogens. A functional genomics screen identified a predicted protein kinase, Epk1, as playing a role in ETI response to three different bacterial effectors. The ETI- and PTI-specific gene sets provide a unique resource for further dissection of the plant immune responses.
Materials and methods
Bacterial strains and growth conditions
Pseudomonas strains were grown on King’s B medium at 30°C. A. tumefaciens and E. coli were grown on Luria-Bertani (LB) medium at 30°C and 37°C, respectively. Antibiotics used were: ampicillin (100 μg/ml), kanamycin (50 μg/ml), rifampicin (10 μg/ml), spectinomycin (50 μg/ml) and gentamycin (10 μg/ml). Bacterial strains are listed in Additional file 12.
Plant material and bacterial infiltrations
To analyze changes in transcript abundance during ETI, 4-week old resistant Rio Grande-PtoR (RG-PtoR), susceptible Rio Grande-prf3 (RG-prf3) or Rio Grande-prf19 (RG-prf19) plants were vacuum-infiltrated with a suspension of 2 × 107cfu/ml of Pseudomonas syringae pv. tomato (Pst) DC3000. The experiment was repeated in three successive weeks (three biological replicates) and leaf samples were collected in each experiment at 4 and 6 hai, frozen in liquid N2 and stored at -80°C until processed.
For comparisons between ETI and PTI, RG-PtoR plants were vacuum-infiltrated with 5 × 106cfu/ml of different DC3000 mutants (DC3000 ΔhopQ1-1 ΔfliC, DC3000 ΔhopQ1-1 ΔavrPto ΔavrPtoB and DC3000 ΔhopQ1-1 ΔavrPto ΔavrPtoB ΔfliC) (Table 1). Leaf samples from infiltrations performed in three successive weeks (three biological replicates) were taken at 6 hai and manipulated as described above. For all experiments, after the samples were collected plants were maintained in the same conditions to observe disease symptoms (3 to 4 days after infiltration).
RNA-Seq library preparation and analysis
Total RNA was isolated with the Plant RNA isolation reagent (Life Technologies, Grand Island, NY, USA) according to the manufacturer’s instructions. Libraries for sequencing were constructed as described in  except that barcode sequences were in the reverse PCR primer. Barcoded libraries were multiplexed by 12 or 15 in each lane and sequenced on an Illumina HiSeq 2000 system using the single-end mode. The length of the reads was around 51 bp. Detailed information about the quality of reads in each replicate is provided in Additional file 13. Analysis of the RNA-Seq data was performed as described in . Genes were considered to be induced when the expression was higher in RG-PtoR than in RG-prf3 or RG-prf19 plants infiltrated with DC3000; or when the expression was higher in RG-PtoR plants inoculated with DC3000 ΔfliC or DC3000 ΔavrPto ΔavrPtoB compared with DC3000 ΔavrPto ΔavrPtoB ΔfliC. In the opposite situation the genes were considered to be suppressed. Sequence reads have been deposited in the NCBI Sequence Read Archive (SRA) under accession number SRP043126 and SRP043127. Processed data are available from the Tomato Functional Genomics Database  (accessions D010 and D011).
Development of effector-triggered and pattern-triggered immunity reporter genes
RG-PtoR plants were vacuum-infiltrated with the DC3000 strains as shown in Table 1. Leaf tissue was collected from four biological replicates at 3, 6, 9 and 12 hai. Total RNA was isolated using the Plant RNA isolation reagent (Life Technologies) following the manufacturer’s instructions. Total RNA (8 μg) was processed with TURBO DNA-free kit (Life Technologies) for 60 minutes at 37°C. After DNase treatment, 4 μg RNA was used to prepare cDNA using SuperScript III First-Strand Synthesis System (Life Technologies) with oligo(dT)20. qRT-PCR was performed as described previously . The sequences of the primers used were: Solyc09g092500 F: 5′-TTGGACAGATCAAGGGACTAATG -3′, R: 5′-CACTCTCAACCACACCATCTT-3′; Solyc10g085880 F: 5′-CCTGGATTGTTCGACAAGAT-3′, R: 5′-CTCCTCCGCTTTCTTCATTT-3′; Solyc04g072280 F: 5′-AACGTCCCGATCGTAGAA-3′, R: 5′-GGATGATCAACTCCACCTAATAA-3′; Solyc02g069960 F: 5′-AGCCAACAAAGCTCAGGAA-3′, R: 5′-CATCCCAGTTGCCATGTTCTA-3′; Solyc11g044390 F: 5′- TCCTAATGACTTGTCCGGATTT-3′, R: 5′-AGTATCACTAGGGCAAGCAAATA-3′; Solyc04g077180 F: 5′-CAGCATTCTGTGGGCTATAC-3′, R: 5′- CCGAAGAAGAAGAGGTTTCC-3′. Data were normalized using: SlATPase (Solyc04g081090) F: 5′-TTGCTGAAGCCTTGGCTCTTTACG-3′, R: 5′-ACCAGCGCGAGAAGAAAGGATGAT-3′; SlEF1α (Solyc06g005060) F: 5′-TCCAAAGATGGTCAGACCCGTGAA-3′, R: 5′-ATACCTAGCCTTGGAGTACTTGGG-3′ and SlCBL1 (Solyc12g015870) F: 5′-CCATCCAAATGCTCCGATCGATGA-3′, R: 5′-TGCCTCTCAATGAAGCCTTGTTGC-3′. Cycling conditions during qRT-PCR were 50°C for 2 minutes, 95°C for 10 minutes, and 40 cycles of 95°C for 30 s, 55°C for 30 s and 72°C for 30 s. A Tukey’s HSD test (P <0.01) was used for statistical analysis of the results.
Virus-induced gene silencing
Tomato genes for which a clear ortholog could be identified in N. benthamiana were selected for designing VIGS constructs. Fragments of 280 to 350 bp were chosen using a VIGS tool . Primers were designed inside this region using Primer3 . PCR amplification was performed using cDNA obtained from RG-PtoR leaves tissue infiltrated with 2 × 107cfu/ml DC3000 or N. benthamiana 35S:Pto plants infiltrated with 2 × 107cfu/ml DC3000 ΔhopQ1-1. PCR products were cloned into pCR8/GW/TOPO vector (Life Technologies) and recombined into a Gateway compatible TRV2 vector . After sequence confirmation, plasmids were transformed into A. tumefaciens GV2260 and gene silencing was performed as described previously .
Screening of effector-triggered immunity-specific genes using programmed cell death assays
A. tumefaciens GV2260 strains carrying pBTEX:Pto or pBTEX:avrPto were grown overnight on solid LB medium (rifampicin and kanamycin). The following day the strains were incubated for 5 h in induction medium (0.05 M MES, 0.5% D-glucose, 0.0272% NaH2PO4, pH 5.6) supplemented with 20× AB salts (2% NH4Cl, 0.6% MgSO4.7H2O, 0.3% KCl, 0.005% FeSO4.7H2O, 0.2% CaCl2.2H2O) and 200 μM acetosyringone. Induced cultures were washed with 10 mM MES, 10 mM MgCl2 and combined in a 1:1 ratio with a final OD600 = 0.15 and 200 μM acetosyringone. The resulting suspensions were syringe-infiltrated into leaves of the silenced plants. Scoring was performed visually during several days as follows: PCD = more than 20% of the infiltrated area (within a marked circle) showed cell death; no PCD = less than 20% of the infiltrated area showed cell death. The number of infiltrated areas in each category was then used to calculate the PCD percentage. Plants silenced for GFP and an E. coli gene-based fragment (Ec1, which contains a 56% GC and not a single ≥17 bp-long 100% identical stretch in N. benthamiana) were used as negative controls. MAPKKKα, MEK2 and Prf were used as positive controls. The gene fragments used for silencing are provided in Additional file 7.
Measurement of reactive oxygen species production
Discs from young leaves of Ec1-, SlEpk1- and NbFls2-silenced plants were excised with a 4-mm-diameter cork borer. Leaf disks were floated adaxial side up in a 96-well black plate (Greiner Bio-one, Kremsmuenster, Austria) containing 180 μl of water per well, and left at room temperature overnight. The next day, the water was removed, and 100 μl of a solution containing the following was added: 500 nM flg22 (GenScript, Piscataway, NJ USA), luminol at 34 μg/ml, and horseradish peroxidase at 20 μg/ml (Sigma, St Louis, MO, USA) in water. Luminescence was measured using the Synergy HT plate reader (Biotek, Winooski, VT, USA). Four leaf disks per plant were taken, and six plants silenced for each gene were considered in each experiment.
Disease assays using Pseudomonas syringae pv. tabaci
N. benthamiana (Nb) or Nb-35S:Pto VIGS-silenced plants were syringe- or vacuum-infiltrated with 5 × 104cfu/ml P. s. tabaci expressing HopQ1-1, AvrPto or carrying an empty vector. The plants were scored as described in the PCD assay above for the presence of disease-associated cell-death and used for bacterial growth assays. Fisher’s exact test (P <0.05) was used to determine significant differences.
Bacterial population assays
Seven-week-old VIGS-silenced plants were vacuum-infiltrated with a suspension of 5 × 104cfu/ml P. s. tabaci (HopQ1-1) or 6 × 104cfu/ml P. s. tabaci (AvrPto) in 10 mM MgCl2 and 0.002% Silwet L-77. To measure bacterial populations, three 0.43 cm2 disks were taken from the oldest expanding leaves and processed twice in a Tissue Lyser (Qiagen, Germantown, MD, USA) for 30 s at 25/s frequency with 0.25 ml of 10 mM MgCl2. The volume was adjusted to 1 ml and serial dilutions were plated on solid LB medium with antibiotics. In each experiment, six biological replicates per construct were used. A Student’s t-test (P <0.05) was used to determine significant differences.
Agrobacterium-mediated transient assays
ETI elicitors were transiently expressed in VIGS-silenced N. benthamiana plants to induce PCD. Pto, AvrPto1-387, Rx2/CP, RPP13/ATR13′41 and Gpa2/RBP-1 were expressed under the 35S promoter (pBTEX), SlMAPKKKα full-length, SlMAPKKKα KD-, MEK2DD and MEK2 WT were expressed using an estradiol-inducible system (pER8), and Bax was expressed using a dexamethasone-inducible system (pTA7002). Expression was induced 48 hai with 2 μM estradiol or 10 μM dexamethasone mixed with 0.02% Tween 20.
Total RNA isolation, DNAse treatment, cDNA synthesis and qRT-PCR were performed as described above, using N. benthamiana leaf tissue from Ec1- and Epk1-silenced plants. The genes chosen for the analysis are shown in Additional file 10. Primers used for each gene analyzed were: Nb00051202g0009.1 F: 5′-TGTTGGGTCAAATGATTCTCACA-3′, R: 5′- GCTCAACCCATTAGAAACTCTGA-3′; Nb00020954g0005.1 F: 5′-TGTCGGGTCAAATGATTCTC-3′, R: 5′-CAACCCGTTAGAAACTCTCC-3′; Nb00014536g0001.1 F: 5′-AGGCTACCTAATGATGATGAAA-3′, R: 5′-GGTGGAGTTGAGACAATAGAG-3′; Nb00003176g0019.1 F: 5′-ACCATGAACATTATGGACTGT-3′, R: 5′-CAACAGAACCTCCACCATT-3′; Nb00029791g0013.1 F: 5′-AGAGGCTCCAAAGTTCGCA-3′, R: 5′-AACCGAACCTCCACCAAT-3′; Nb00042373g0002.1 F: 5′-ATGATAAGACTACCTGATGATGAC-3′, R: 5′-GGTGGAGTTGAGATAAGAGAATAA-3′. Data were normalized using NbPP2a and NbEF1α . Cycling conditions during qRT-PCR were 50°C for 2 minutes, 95°C for 10 minutes, and 40 cycles of 95°C for 30 s, 58°C for 30 s and 72°C for 30 s. A pairwise Student’s t-test (P <0.01) was used to determine significance differences.
Gene Ontology term analysis
Genes with expression levels ≥3 RPKM in at least one treatment, >2-fold induction and P <0.05 were used for the GO term analysis using the GO::Term Finder module .
The SeaView program  was used to perform the phylogeny analysis with the default parameters for the GTR model (nucleotides) and the JTT model (proteins). PhyML with one hundred bootstraps was used for each analysis. The tree figure was created using FigTree .
Bibliographic search associated with effector-triggered and pattern-triggered immunity-induced genes
SwissProt and TrEMBL protein databases , and EMBL-EBI plant EST and STD nucleotides databases  were downloaded and filtered using Ruby custom scripts to obtain plant sequences with PUBMED publications . These sequences were compared with the tomato gene models (ITAG2.3) using BLAST  to find putative orthologs to the tomato genes. Custom Ruby and BioRuby  scripts, available upon request, were used to filter the bibliography to obtain publications related to plant immunity associated with these genes.
The RNA-Seq data are available from the NCBI (accessions SRP043126 and SRP043127) and the Tomato Functional Genomics Database  (accessions D010 and D011).
MAP and GBM designed the research, MAP and DMD performed the research, MAP, YZ, NFP, ZF and GBM analyzed the data, MAP and GBM wrote the paper. All authors read and approved the manuscript.
green fluorescent protein
hours after inoculation
microbe-associated molecular pattern
mitogen-activated protein kinase
nucleotide binding-leucine-rich repeat
programmed cell death
- Pst Pseudomonas syringae :
quantitative reverse transcription-polymerase chain reaction
reactive oxygen species
reads per kilobase of exon model per million mapped reads
virus-induced gene silencing
We thank Paige Reeves for her assistance in the greenhouse, Hailei Wei and Alan Collmer for providing advice and bacterial strains, and Hernan Rosli and Simon Schwizer for helpful comments on the manuscript. This research was supported by National Science Foundation Plant Genome Research Program grant IOS-1025642 (GBM, Alan Collmer, PI) and by Israel-USA Binational Agriculture Research and Development Fund grant IS-4510-12C (GBM).
- Schwessinger B, Ronald PC: Plant innate immunity: perception of conserved microbial signatures. Annu Rev Plant Biol. 2012, 63: 451-482. 10.1146/annurev-arplant-042811-105518.PubMedView ArticleGoogle Scholar
- Boller T, He SY: Innate immunity in plants: an arms race between pattern recognition receptors in plants and effectors in microbial pathogens. Science. 2009, 324: 742-744. 10.1126/science.1171647.PubMedPubMed CentralView ArticleGoogle Scholar
- Dodds PN, Rathjen JP: Plant immunity: towards an integrated view of plant-pathogen interactions. Nat Rev Genet. 2010, 11: 539-548. 10.1038/nrg2812.PubMedView ArticleGoogle Scholar
- Segonzac C, Zipfel C: Activation of plant pattern-recognition receptors by bacteria. Curr Opin Microbiol. 2011, 14: 54-61. 10.1016/j.mib.2010.12.005.PubMedView ArticleGoogle Scholar
- Nguyen HP, Chakravarthy S, Velásquez AC, McLane HS, Zeng L, Park D-W, Collmer A, Martin GB: Methods to study PAMP-triggered immunity using tomato and Nicotiana benthamiana. Mol Plant Microbe Interact. 2010, 23: 991-999. 10.1094/MPMI-23-8-0991.PubMedView ArticleGoogle Scholar
- Maekawa T, Kufer TA, Schulze-Lefert P: NLR functions in plant and animal immune systems: so far and yet so close. Nat Immunol. 2011, 12: 817-826. 10.1038/ni.2083.PubMedView ArticleGoogle Scholar
- Jones JD, Dangl JL: The plant immune system. Nature. 2006, 444: 323-329. 10.1038/nature05286.PubMedView ArticleGoogle Scholar
- Pedley KF, Martin GB: Molecular basis of Pto-mediated resistance to bacterial speck disease in tomato. Annu Rev Phytopathol. 2003, 41: 215-243. 10.1146/annurev.phyto.41.121602.143032.PubMedView ArticleGoogle Scholar
- Oh CS, Martin GB: Effector-triggered immunity mediated by the Pto kinase. Trends Plant Sci. 2011, 16: 132-140. 10.1016/j.tplants.2010.11.001.PubMedView ArticleGoogle Scholar
- Gomez-Gomez L, Boller T: FLS2: an LRR receptor-like kinase involved in the perception of the bacterial elicitor flagellin in Arabidopsis. Mol Cell. 2000, 5: 1003-1011. 10.1016/S1097-2765(00)80265-8.PubMedView ArticleGoogle Scholar
- Zeng W, He SY: A prominent role of the flagellin receptor FLAGELLIN-SENSING2 in mediating stomatal response to Pseudomonas syringae pv tomato DC3000 in Arabidopsis. Plant Physiol. 2010, 153: 1188-1198. 10.1104/pp.110.157016.PubMedPubMed CentralView ArticleGoogle Scholar
- Clarke CR, Chinchilla D, Hind SR, Taguchi F, Miki R, Ichinose Y, Martin GB, Leman S, Felix G, Vinatzer BA: Allelic variation in two distinct Pseudomonas syringae flagellin epitopes modulates the strength of plant immune responses but not bacterial motility. New Phytol. 2013, 200: 847-860. 10.1111/nph.12408.PubMedPubMed CentralView ArticleGoogle Scholar
- Rosli H, Zheng Y, Pombo M, Zhong S, Bombarely A, Fei Z, Collmer A, Martin G: Transcriptomics-based screen for genes induced by flagellin and repressed by pathogen effectors identifies a cell wall-associated kinase involved in plant immunity. Genome Biol. 2013, 14: R139-10.1186/gb-2013-14-12-r139.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin GB: Suppression and activation of the plant immune system by Pseudomonas syringae effectors AvrPto and AvrPtoB. Effectors in Plant-Microbe Interactions. Edited by: Martin F, Kamoun S. 2012, Wiley-Blackwell, John Wiley and Sons, London, 123-154.Google Scholar
- Cunnac S, Chakravarthy S, Kvitko BH, Russell AB, Martin GB, Collmer A: Genetic disassembly and combinatorial reassembly identify a minimal functional repertoire of type III effectors in Pseudomonas syringae. Proc Natl Acad Sci U S A. 2011, 108: 2975-2980. 10.1073/pnas.1013031108.PubMedPubMed CentralView ArticleGoogle Scholar
- Martin GB, Brommonschenkel SH, Chunwongse J, Frary A, Ganal MW, Spivey R, Wu T, Earle ED, Tanksley SD: Map-based cloning of a protein kinase gene conferring disease resistance in tomato. Science. 1993, 262: 1432-1436. 10.1126/science.7902614.PubMedView ArticleGoogle Scholar
- Tang X, Frederick RD, Zhou J, Halterman DA, Jia Y, Martin GB: Initiation of plant disease resistance by physical interaction of AvrPto and Pto kinase. Science. 1996, 274: 2060-2063. 10.1126/science.274.5295.2060.PubMedView ArticleGoogle Scholar
- Salmeron JM, Oldroyd GED, Rommens CMT, Scofield SR, Kim HS, Lavelle DT, Dahlbeck D, Staskawicz BJ: Tomato Prf is a member of the leucine-rich repeat class of plant disease resistance genes and lies embedded within the Pto kinase gene cluster. Cell. 1996, 86: 123-133. 10.1016/S0092-8674(00)80083-5.PubMedView ArticleGoogle Scholar
- Mysore KS, Crasta OR, Tuori RP, Folkerts O, Swirsky PB, Martin GB: Comprehensive transcript profiling of Pto- and Prf-mediated host defense responses to infection by Pseudomonas syringae pv. tomato. Plant J. 2002, 32: 299-315. 10.1046/j.1365-313X.2002.01424.x.PubMedView ArticleGoogle Scholar
- Tsuda K, Katagiri F: Comparing signaling mechanisms engaged in pattern-triggered and effector-triggered immunity. Curr Opin Plant Biol. 2010, 13: 459-465. 10.1016/j.pbi.2010.04.006.PubMedView ArticleGoogle Scholar
- Katagiri F, Tsuda K: Understanding the plant immune system. Mol Plant Microbe Interact. 2010, 23: 1531-1536. 10.1094/MPMI-04-10-0099.PubMedView ArticleGoogle Scholar
- Tao Y, Xie Z, Chen W, Glazebrook J, Chang HS, Han B, Zhu T, Zou G, Katagiri F: Quantitative nature of Arabidopsis responses during compatible and incompatible interactions with the bacterial pathogen Pseudomonas syringae. Plant Cell. 2003, 15: 317-330. 10.1105/tpc.007591.PubMedPubMed CentralView ArticleGoogle Scholar
- Navarro L, Zipfel C, Rowland O, Keller I, Robatzek S, Boller T, Jones JD: The transcriptional innate immune response to flg22. Interplay and overlap with Avr gene-dependent defense responses and bacterial pathogenesis. Plant Physiol. 2004, 135: 1113-1128. 10.1104/pp.103.036749.PubMedPubMed CentralView ArticleGoogle Scholar
- Tsuda K, Sato M, Glazebrook J, Cohen JD, Katagiri F: Interplay between MAMP-triggered and SA-mediated defense responses. Plant J. 2008, 53: 763-775. 10.1111/j.1365-313X.2007.03369.x.PubMedView ArticleGoogle Scholar
- Tsuda K, Sato M, Stoddard T, Glazebrook J, Katagiri F: Network properties of robust immunity in plants. PLoS Genet. 2009, 5: e1000772-10.1371/journal.pgen.1000772.PubMedPubMed CentralView ArticleGoogle Scholar
- Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev. 2009, 10: 57-63. 10.1038/nrg2484.View ArticleGoogle Scholar
- Ansorge WJ: Next generation DNA sequencing techniques. Nat Biotechnol. 2009, 25: 195-203.Google Scholar
- Haas BJ, Zody MC: Advancing RNA-Seq analysis. Nat Biotechnol. 2010, 28: 421-423. 10.1038/nbt0510-421.PubMedView ArticleGoogle Scholar
- Wang L, Li P, Brutnell TP: Exploring plant transcriptome using ultra high-throughput sequencing. Brief Funct Genomics. 2010, 9: 118-128. 10.1093/bfgp/elp057.PubMedView ArticleGoogle Scholar
- Salmeron JM, Barker SJ, Carland FM, Mehta AY, Staskawicz BJ: Tomato mutants altered in bacterial disease resistance provide evidence for a new controlling pathogen recognition. Plant Cell. 1994, 6: 511-520. 10.1105/tpc.6.4.511.PubMedPubMed CentralView ArticleGoogle Scholar
- iTAK - Plant Transcription factor and Protein Kinase Identifier and Classifier , [http://bioinfo.bti.cornell.edu/cgi-bin/itak/index.cgi]
- del Pozo O, Pedley KF, Martin GB: MAPKKKa is a positive regulator of cell death associated with both plant immunity and disease. EMBO J. 2004, 23: 3072-3082. 10.1038/sj.emboj.7600283.PubMedView ArticleGoogle Scholar
- Ekengren SK, Liu Y, Schiff M, Dinesh-Kumar SP, Martin GB: Two MAPK cascades, NPR1, and TGA transcription factors play a role in Pto-mediated disease resistance in tomato. Plant J. 2003, 36: 905-917. 10.1046/j.1365-313X.2003.01944.x.PubMedView ArticleGoogle Scholar
- Wulf J, Pascuzzi PE, Fahmy A, Martin GB, Nicholson LK: The solution structure of type III effector protein AvrPto reveals conformational and dynamic features important for plant pathogenesis. Structure. 2004, 12: 1257-1268. 10.1016/j.str.2004.04.017.PubMedView ArticleGoogle Scholar
- Wei CF, Kvitko BH, Shimizu R, Crabill E, Alfano JR, Lin NC, Martin GB, Huang HC, Collmer A: A Pseudomonas syringae pv. tomato DC3000 mutant lacking the type III effector HopQ1-1 is able to cause disease in the model plant Nicotiana benthamiana. Plant J. 2007, 51: 32-46. 10.1111/j.1365-313X.2007.03126.x.PubMedView ArticleGoogle Scholar
- Rosebrock TR, Zeng L, Brady JJ, Abramovitch RB, Xiao F, Martin GB: A bacterial E3 ubiquitin ligase targets a host protein kinase to disrupt plant immunity. Nature. 2007, 448: 370-374. 10.1038/nature05966.PubMedPubMed CentralView ArticleGoogle Scholar
- Bendahmane A, Querci M, Kanyuka K, Baulcombe DC: Agrobacterium transient expression system as a tool for the isolation of disease resistance genes: application to the Rx2 locus in potato. Plant J. 2000, 21: 73-81. 10.1046/j.1365-313x.2000.00654.x.PubMedView ArticleGoogle Scholar
- Rentel MC, Leonelli L, Dahlbeck D, Zhao B, Staskawicz BJ: Recognition of the Hyaloperonospora parasitica effector ATR13 triggers resistance against oomycete, bacterial, and viral pathogens. Proc Natl Acad Sci U S A. 2008, 105: 1091-1096. 10.1073/pnas.0711215105.PubMedPubMed CentralView ArticleGoogle Scholar
- Sacco MA, Koropacka K, Grenier E, Jaubert MJ, Blanchard A, Goverse A, Smant G, Moffett P: The cyst nematode SPRYSEC protein RBP-1 elicits Gpa2- and RanGAP2-dependent plant cell death. PLoS Pathog. 2009, 5: e1000564-10.1371/journal.ppat.1000564.PubMedPubMed CentralView ArticleGoogle Scholar
- Lacomme C, Santa Cruz S: Bax-induced cell death in tobacco is similar to the hypersensitive response. Proc Natl Acad Sci U S A. 1999, 96: 7956-7961. 10.1073/pnas.96.14.7956.PubMedPubMed CentralView ArticleGoogle Scholar
- Xin X-F, He SY: Pseudomonas syringae pv. tomato DC3000: a model pathogen for probing disease susceptibility and hormone signaling in plants. Annu Rev Phytopathol. 2013, 51: 473-498. 10.1146/annurev-phyto-082712-102321.PubMedView ArticleGoogle Scholar
- Mur LAJ, Kenton P, Lloyd AJ, Ougham H, Prats E: The hypersensitive response; the centenary is upon us but how much do we know?. J Exp Bot. 2008, 59: 501-520. 10.1093/jxb/erm239.PubMedView ArticleGoogle Scholar
- Coll NS, Epple P, Dangl JL: Programmed cell death in the plant immune system. Cell Death Differ. 2011, 18: 1247-1256. 10.1038/cdd.2011.37.PubMedPubMed CentralView ArticleGoogle Scholar
- Kang L, Li J, Zhao T, Xiao F, Tang X, Thilmony R, He S, Zhou J-M: Interplay of the Arabidopsis nonhost resistance gene NHO1 with bacterial virulence. Proc Natl Acad Sci U S A. 2003, 100: 3519-3524. 10.1073/pnas.0637377100.PubMedPubMed CentralView ArticleGoogle Scholar
- Li X, Lin H, Zhang W, Zou Y, Zhang J, Tang X, Zhou J-M: Flagellin induces innate immunity in nonhost interactions that is suppressed by Pseudomonas syringae pv. tomato effectors. Proc Natl Acad Sci U S A. 2005, 102: 12990-12995. 10.1073/pnas.0502425102.PubMedPubMed CentralView ArticleGoogle Scholar
- He P, Shan L, Lin NC, Martin GB, Kemmerling B, Nurnberger T, Sheen J: Specific bacterial suppressors of MAMP signaling upstream of MAPKKK in Arabidopsis innate immunity. Cell. 2006, 125: 563-575. 10.1016/j.cell.2006.02.047.PubMedView ArticleGoogle Scholar
- Schulze S, Kay S, Buttner D, Egler M, Eschen-Lippold L, Hause G, Kruger A, Lee J, Muller O, Scheel D, Szczesny R, Thieme F, Bonas U: Analysis of new type III effectors from Xanthomonas uncovers XopB and XopS as suppressors of plant immunity. New Phytol. 2012, 195: 894-911. 10.1111/j.1469-8137.2012.04210.x.PubMedView ArticleGoogle Scholar
- Kim JG, Li X, Roden JA, Taylor KW, Aakre CD, Su B, Lalonde S, Kirik A, Chen Y, Baranage G, McLane H, Martin GB, Mudgett MB: Xanthomonas T3S effector XopN suppresses PAMP-triggered immunity and interacts with a tomato atypical receptor-like kinase and TFT1. Plant Cell. 2009, 21: 1305-1323. 10.1105/tpc.108.063123.PubMedPubMed CentralView ArticleGoogle Scholar
- Zheng X, McLellan H, Fraiture M, Liu X, Boevink PC, Gilroy EM, Chen Y, Kandel K, Sessa G, Birch PR, Brunner F: Functionally redundant RXLR effectors from Phytophthora infestans act at different steps to suppress early flg22-triggered immunity. PLoS Pathog. 2014, 10: e1004057-10.1371/journal.ppat.1004057.PubMedPubMed CentralView ArticleGoogle Scholar
- Ichimura K, Shinozaki K, Tena G, Sheen J, Henry Y, Champion A, Kreis M, Zhang SQ, Hirt H, Wilson C, Heberle-Bors E, Ellis BE, Morris PC, Innes RW, Ecker JR, Scheel D, Klessig DF, Machida Y, Mundy J, Ohashi Y, Walker JC, Grp M: Mitogen-activated protein kinase cascades in plants: a new nomenclature. Trends Plant Sci. 2002, 7: 301-308. 10.1016/S1360-1385(02)02302-6.View ArticleGoogle Scholar
- Martin GB, Bogdanove AJ, Sessa G: Understanding the functions of plant disease resistance proteins. Annu Rev Plant Biol. 2003, 54: 23-61. 10.1146/annurev.arplant.54.031902.135035.PubMedView ArticleGoogle Scholar
- Bernoux M, Ellis JG, Dodds PN: New insights in plant immunity signaling activation. Curr Opin Plant Biol. 2011, 14: 512-518. 10.1016/j.pbi.2011.05.005.PubMedPubMed CentralView ArticleGoogle Scholar
- Jin H, Liu Y, Yang KY, Kim CY, Baker B, Zhang S: Function of a mitogen-activated protein kinase pathway in N gene-mediated resistance in tobacco. Plant J. 2003, 33: 719-731. 10.1046/j.1365-313X.2003.01664.x.PubMedView ArticleGoogle Scholar
- Oh CS, Pedley KF, Martin GB: Tomato 14-3-3 protein 7 positively regulates immunity-associated programmed cell death by enhancing protein abundance and signaling ability of MAPKKK alpha. Plant Cell. 2010, 22: 260-272. 10.1105/tpc.109.070664.PubMedPubMed CentralView ArticleGoogle Scholar
- Cho HS, Yoon GM, Lee SS, Kim YA, Hwang I, Choi D, Pai H-S: A novel dual-specificity protein kinase targeted to the chloroplast in tobacco. FEBS Lett. 2001, 497: 124-130. 10.1016/S0014-5793(01)02439-5.PubMedView ArticleGoogle Scholar
- Rudrabhatla P, Reddy MM, Rajasekharan R: Genome-wide analysis and experimentation of plant serine/ threonine/tyrosine-specific protein kinases. Plant Mol Biol. 2006, 60: 293-319. 10.1007/s11103-005-4109-7.PubMedView ArticleGoogle Scholar
- Macho AP, Schwessinger B, Ntoukakis V, Brutus A, Segonzac C, Roy S, Kadota Y, Oh MH, Sklenar J, Derbyshire P, Lozano-Duran R, Malinovsky FG, Monaghan J, Menke FL, Huber SC, He SY, Zipfel C: A bacterial tyrosine phosphatase inhibits plant pattern recognition receptor activation. Science. 2014, 343: 1509-1512. 10.1126/science.1248849.PubMedView ArticleGoogle Scholar
- Bayer RG, Stael S, Rocha AG, Mair A, Vothknecht UC, Teige M: Chloroplast-localized protein kinases: a step forward towards a complete inventory. J Exp Bot. 2012, 63: 1713-1723. 10.1093/jxb/err377.PubMedPubMed CentralView ArticleGoogle Scholar
- Zhong S, Joung JG, Zheng Y, Chen YR, Liu B, Shao Y, Xiang JZ, Fei Z, Giovannoni JJ: High-throughput illumina strand-specific RNA sequencing library preparation. Cold Spring Harb Protoc. 2011, 2011: 940-949. 10.1101/pdb.prot5652.PubMedView ArticleGoogle Scholar
- Tomato Functional Genomics Database , [http://ted.bti.cornell.edu/]
- VIGS Tool , [http://solgenomics.net/tools/vigs]
- Untergasser A, Cutcutache I, Koressaar T, Ye J, Faircloth BC, Remm M, Rozen SG: Primer3–new capabilities and interfaces. Nucleic Acids Res. 2012, 40: e115-10.1093/nar/gks596.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu Y, Schiff M, Dinesh-Kumar SP: Virus-induced gene silencing in tomato. Plant J. 2002, 31: 777-786. 10.1046/j.1365-313X.2002.01394.x.PubMedView ArticleGoogle Scholar
- Velasquez AC, Chakravarthy S, Martin GB: Virus-induced gene silencing (VIGS) in Nicotiana benthamiana and tomato. J Vis Exp. 2009, 28: 1292-1299.PubMedGoogle Scholar
- Liu D, Shi L, Han C, Yu J, Li D, Zhang Y: Validation of reference genes for gene expression studies in virus-infected Nicotiana benthamiana using quantitative real-time PCR. PLoS One. 2012, 7: e46451-10.1371/journal.pone.0046451.PubMedPubMed CentralView ArticleGoogle Scholar
- Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, Sherlock G: GO::TermFinder–open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004, 20: 3710-3715. 10.1093/bioinformatics/bth456.PubMedPubMed CentralView ArticleGoogle Scholar
- Gouy M, Guindon S, Gascuel O: SeaView version 4: a multiplatform graphical user interface for sequence alignment and phylogenetic tree building. Mol Biol Evol. 2010, 27: 221-224. 10.1093/molbev/msp259.PubMedView ArticleGoogle Scholar
- FigTree , [http://tree.bio.ed.ac.uk/software/figtree]
- Consortium TU: Activities at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2014, 42: D191-D198. 10.1093/nar/gkt1140.View ArticleGoogle Scholar
- Stoesser G, Baker W, van den Broek A, Camon E, Garcia-Pastor M, Kanz C, Kulikova T, Lombard V, Lopez R, Parkinson H, Redaschi N, Sterk P, Stoehr P, Tuli MA: The EMBL nucleotide sequence database. Nucleic Acids Res. 2001, 29: 17-21. 10.1093/nar/29.1.17.PubMedPubMed CentralView ArticleGoogle Scholar
- PUBMED , [http://www.ncbi.nlm.nih.gov/pubmed/]
- Camacho C, Coulouris G, Avagyan V, Ma N, Papadopoulos J, Bealer K, Madden T: BLAST+: architecture and applications. BMC Bioinformatics. 2009, 10: 421-10.1186/1471-2105-10-421.PubMedPubMed CentralView ArticleGoogle Scholar
- Goto N, Prins P, Nakao M, Bonnal R, Aerts J, Katayama T: BioRuby: bioinformatics software for the Ruby programming language. Bioinformatics. 2010, 26: 2617-2619. 10.1093/bioinformatics/btq475.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.