Genomic screens identify a new phytobacterial microbe-associated molecular pattern and the cognate Arabidopsis receptor-like kinase that mediates its immune elicitation
- G. Adam Mott1,
- Shalabh Thakur1,
- Elwira Smakowska2,
- Pauline W. Wang3,
- Youssef Belkhadir2,
- Darrell Desveaux†1, 3Email author and
- David S. Guttman†1, 3Email author
© Mott et al. 2016
Received: 29 June 2015
Accepted: 14 April 2016
Published: 9 May 2016
The recognition of microbe-associated molecular patterns during infection is central to the mounting of an effective immune response. In spite of their importance, it remains difficult to identify these molecules and the host receptors required for their perception, ultimately limiting our understanding of the role of these molecules in the evolution of host-pathogen relationships.
We employ a comparative genomics screen to identify six new immune eliciting peptides from the phytopathogenic bacterium Pseudomonas syringae. We then perform a reverse genetic screen to identify Arabidopsis thaliana leucine-rich repeat receptor-like kinases required for the recognition of these elicitors. We test the six elicitors on 187 receptor-like kinase knock-down insertion lines using a high-throughput peroxidase-based immune assay and identify multiple lines that show decreased immune responses to specific peptides. From this primary screen data, we focused on the interaction between the xup25 peptide from a bacterial xanthine/uracil permease and the Arabidopsis receptor-like kinase xanthine/uracil permease sensing 1; a family XII protein closely related to two well-characterized receptor-like kinases. We show that xup25 treatment increases pathogenesis-related gene induction, callose deposition, seedling growth inhibition, and resistance to virulent bacteria, all in a xanthine/uracil permease sensing 1-dependent manner. Finally, we show that this kinase-like receptor can bind the xup25 peptide directly. These results identify xup25 as a P. syringae microbe-associated molecular pattern and xanthine/uracil permease sensing 1 as a receptor-like kinase that detects the xup25 epitope to activate immune responses.
The present study demonstrates an efficient method to identify immune elicitors and the plant receptors responsible for their perception. Further exploration of these molecules will increase our understanding of plant-pathogen interactions and the basis for host specificity.
Effective immunity in plants relies upon a multi-tiered innate immune recognition system to successfully identify and appropriately respond to microbial invaders [1, 2]. This response requires the ability to quickly detect the presence of potential pathogens, effective mechanisms to disseminate that information through the organism, and finally appropriate physiological responses capable of controlling and clearing infection. The initial threat detection is largely accomplished through the recognition of microbe-associated molecular patterns (MAMPs), which are highly conserved immune elicitors derived from invading microbes. In the case of proteinaceous MAMPs, the genes encoding these microbial signatures are likely required for survival and, therefore, under strong negative selection as a whole . If this were not the case, they would be modified or eliminated through natural selection in order to subvert host recognition. Nevertheless, while the genes encoding MAMPs are under strong negative selection overall, individual residues can show signals of positive selection for diversity . This variation in what are otherwise conserved proteins may help the microbe avoid or dampen host recognition. In fact, MAMP sequence diversity has been shown to be associated with variation in the intensity of the immune response elicited by MAMP peptides [4–6].
The recognition of MAMPs is mediated through direct binding to a pattern-recognition receptor (PRR) on the plant cell surface. Characterized plant PRRs of peptide MAMPs contain an extracellular domain with a number of leucine-rich repeat (LRR) domains responsible for MAMP binding and can be divided based on their intracellular moieties into receptor-like kinases (RLK) and the closely related receptor-like proteins (RLP), which lack an intracellular kinase domain [7, 8]. These two gene families in Arabidopsis thaliana have 223 and 57 LRR-containing members, respectively [9–11]; a level of genetic diversification suggestive of the key role these proteins play in plant evolutionary success. Immune signaling requires PRRs to work in complex with regulatory co-receptors, such as the Brassinosteroid Insensitive1-associated receptor kinase1 (BAK1) that regulates LRR-containing PRR activation [12, 13]. This protein is an LRR-RLK that interacts with the PRR and appears to help shape the receptor pocket that binds directly to the MAMP [12, 13]. While this co-receptor is not required for PRR-MAMP binding, it does participate in MAMP binding and both of these LRR-RLKs are required for full activation of immune-associated receptor-like cytoplasmic kinases that transduce the immune activation signal.
The importance of MAMPs to plant immunity has led to considerable interest in their identification, yet surprisingly few MAMPs have been characterized. In fact, the majority of work on proteinaceous MAMPs to date has focused on flagellin and elongation factor Tu (EF-Tu) and their respective eliciting peptides flg22 and elf18 [5, 6]. This is due in no small part to the technical challenges involved in the identification of both MAMPs and the plant components required for their perception using traditional means. More recently, with the proliferation of high-throughput sequencing and advances in bioinformatic techniques, predictive methods for MAMP identification have arisen. We have previously described such an approach based on the unique evolutionary signature of MAMP peptides that arises due to strong negative pressure to maintain protein function coupled to positive pressure to mutate sites in order to avoid plant perception . While this early demonstration successfully identified peptides that elicited an immune response, it was not aimed at identifying the cognate host receptors required for recognition of these elicitors and consequently only addressed half of the issue. Another predictive approach has also recently been used to identify novel damage-associated molecular patterns (DAMPs) using sequence homology to known DAMPs . Such approaches will become even more powerful as more genomes are sequenced and more immune elicitors are identified. Nevertheless, the host signaling components required for perception of these novel elicitors have yet to be identified.
The present study couples our computational pipeline for elicitor identification to a reverse genetic screen for cognate receptors using a high-throughput peroxidase assay. We show here both the identification of a novel Pseudomonas syringae MAMP that stimulates defence gene expression, induces callose deposition in planta, causes seedling growth inhibition, and functionally protects the plant from subsequent pathogen challenge, as well as the identification of an Arabidopsis thaliana LRR-RLK that binds the peptide and is required for the perception of this novel MAMP.
Computational identification of immune elicitors
The screen for novel peptide elicitors was performed by examining P. syringae genomes for genes that show an overall pattern of strong negative selection for maintenance of required protein function, coupled to localized strong positive selection for avoidance of plant perception. We have used this approach previously to identify immune elicitors in P. syringae , but elected to generate a new set of candidates for this study due to the availability of vastly more genome sequences, which provided much more predictive power (see “Discussion” for more details). We examined 54 P. syringae isolates (Additional file 1: Table S1) to identify 3157 ortholog core gene families, as defined by being present in at least 90 % of the isolates. A PAML selection analysis of this core gene set identified 172 genes that had at least one positively selected residue with ω ratio significantly greater than 1.0, and which rejected the neutral and negative selection null model (M7) in favor of the positive selection model (M8). Our previous study found that many putative immune elicitors had clusters of positively selected residues , therefore we used the Bayes Empirical Bayes estimates from the M8 model to narrow the list to those with clusters of at least three positively selected sites in a window of 25 residues, producing a candidate pool of 61 genes (Additional file 1: Table S2). This candidate pool included the well-established MAMPs EF-Tu and flagellin, ranking fifth and ninth, respectively, based on the number of predicted positively selected sites.
Predicted peptide MAMPs
NCBI Gene ID
Xanthine/uracil permease family protein
Conserved hypothetical protein
ABC transporter, permease protein
Outer membrane protein
Membrane protein, putative
High-throughput peroxidase immune elicitation assay
We developed a sensitive, high-throughput, microtiter-plate based screening assay based on innate immunity-induced peroxidase (POX) activity as a quantitative measure of immune function in order to efficiently screen large numbers of candidate elicitor peptides and LRR-RLK T-DNA insertion lines. The apoplastic POX enzymes play a role in the plant response to both wounding and pathogen infection [15–17]. These two responses, however, are temporally distinct, a fact that we used to avoid confounding our assay with wound-induced POX activity caused during tissue harvest. While POX activity has been used previously to study the effects of MAMP treatment , this is its first adaptation to a high-throughput screen using intact plant tissue. Importantly, the assay also showed excellent responses at a variety of doses of flg22, with sensitivity similar to the traditional ROS assay (Additional file 2: Figure S1).
Reverse genetic screen to identify MAMP receptors/co-receptors
We performed a primary screen on the 187 T-DNA LRR-RLK insertion lines for their ability to respond to our six candidate elicitors using the high-throughput POX assay. Each assay involved testing six replicate plants per genotype, with ten leaf cores per plant. All controls and peptides were assayed in parallel. Overall, the 187 T-DNA lines were assayed 3513 times (excluding negative controls), with a median of 496 assays per elicitor, 16.5 assays per T-DNA line, and 2.7 assays per elicitor per line.
The XPS1 LRR-RLK is required for xup25-induced immunity
XPS1 specifically binds xup25
In spite of the important role that MAMPs play in plant immunity, only a small number have been identified, and an even smaller number of proteinaceous MAMPs have been linked to specific LRR-containing receptors. The traditional approaches for MAMP identification rely upon refining complex mixtures of molecules from pathogens, which becomes increasingly difficult when working with MAMPs that have weak biological effects under specific test conditions. In addition, our experimental models can only crudely mimic the complex milieu of microbial commensals, mutualists, and parasites that constantly interact with plants in the real environment. Interactions that we perceive as weak under rarefied laboratory conditions may in fact be extremely important under specific ecological conditions. While these ecological and evolutionary processes normally limit progress, our computational approach to MAMP identification instead takes advantage of these events by mining for the genetic footprints left by these interactions.
This study used comparative genomics analyses to identify peptides carrying genetic signatures consistent with expectations for MAMPs. Six of these immune-eliciting peptides were used in a reverse genetic screen against 187 A. thaliana T-DNA insertion lines of putative PRRs. While most of the knockout lines retained their ability to mount an immune response when challenged with the peptides, one line, xps1, was unable to mount any measurable immune response to the xup25 peptide. Importantly, the xps1 line retained its ability to respond to the other peptides, including flg22. Finally, the XPS1 ectodomain is shown to be involved in the direct binding of xup25.
These results indicate that XPS1 can bind to xup25 and is specifically required for perception and, therefore, xup25 and XPS1 can be considered a MAMP and PRR pair in the strictest sense.
The xup25 peptide is derived from a xanthine/uracil permease family protein (conservation is shown in Additional file 9: Figure S8). There is no clear reason why this protein should be a source of a peptide MAMP beyond the fact that it is highly conserved. In fact, all six candidate elicitors identified in this study are broadly conserved among the bacteria, while some are even found in the archaea (data not shown). Four of the candidates, including the protein encoding xup25, are membrane proteins, which while not a requirement for a MAMP (q.v. EF-Tu), may increase the likelihood that the peptides come in contact with host PRRs.
The amplitude of the immune response we observe with the newly identified peptide elicitors is lower than that observed for flg22. Recent work has shown that natural variation within a MAMP sequence can have large effects on the resulting immune response [4, 29], raising the question of whether we would have observed stronger immune responses if we had designed our peptides based on protein sequences from more divergent bacteria with non-host interactions, rather than P. syringae. As there has not yet been a broad survey of the relationship between MAMP origin and the strength of immune elicitation, we can only speculate on the matter based on the limited data available. In the best studied case, the flg22 peptide from Pseudomonas aeruginosa elicits a stronger immune response in A. thaliana than the corresponding peptide from P. syringae . This finding may seem logical, as P. aeruginosa is a non-host species rarely found on foliar surfaces, suggesting flg22 may be at least partially responsible for this niche restriction. Upon closer examination, this theory is contradicted by the fact that the flg22 peptide does not elicit an efficient immune response in all non-host plant species studied (e.g. celery ), nor do all natural accessions of A. thaliana carry the FLS2 PRR that recognizes flg22 . In fact, the opposite argument could also be made, namely that MAMPs from pathogens of a specific plant host should elicit a strong immune response, since these are ecologically relevant interactions. This argument is supported by the observation that FLS2 restricts the growth of even highly virulent pathogens such as Pto DC3000 . While type III secreted effectors also strongly influence disease and immunity, these proteins are highly variable even among closely related strains [32–34] and, consequently, likely mediate a relatively small number of all plant-microbe interactions. We elected to test candidate elicitor peptides derived from the host strain Pto DC3000, as well as a number of homologous peptides derived from the non-host strain P. syringae pv. phaseolicola 1448a, but found no significant differential response (data not shown). Future work examining homologous peptides from a wider range of species will help inform which of these two positions is more generally correct.
Previous work in our group has used a similar bioinformatics approach to identifying immune eliciting peptides , but the current report represents a significant improvement in several respects. Most notably, the previous work was performed when only three P. syringae isolates had been sequenced, necessitating the inclusion of three strains of Xanthomonas campestris to achieve sufficient sample size for the selection analyses. The inclusion of these divergent genomes strongly influenced the prediction of the set of core genes. Additionally, the use of two quite divergent species groups negatively influenced the quality of the alignments and ultimately the predictions of positive and negative selection as gap positions are ignored in PAML selection analyses. In the earlier report, we characterized the immune eliciting properties of eight bioinformatically predicted peptides. The proteins from which these peptides were derived were also analyzed as part of this report, however the use of expanded genome data from 54 P. syringae strains resulted in our decision to select new candidates for a variety of reasons. Three of the previously identified proteins (PSPTO2757, PSPTO3620, and PSPTO4436) were not present in a sufficient number of strains to be considered part of the core genome in this report. This observation may be true or simply an artefact due to the extensive use of incomplete draft genomes in the current study. Of the remaining five candidates from the previous report, four had positively selected sites (PSPTO3468, PSPTO0175, PSPTO5086, and PSPTO0624), while the last did not show a significant signal of positive selection in the new analysis (PSPTO1253). The inclusion of PSPTO1253 was unique in the previous report as it was based on the presence of a single PSS, making it the least robust prediction of the set. We elected not to use the four previously identified peptides that were also present in this new study since the expanded dataset and subsequent refined analysis identified other candidates with stronger signals based on the number and clustering of positively selected sites. The current analysis of 54 P. syringae genomes provides better sequence alignments, a much more reliable prediction of the species-specific core genome, greatly enhanced statistical power to detect signatures of positive selection, and ultimately a more robust prediction of MAMP candidates.
XPS1 is notable as it forms a part of the LRR-RLK family XII along with FLS2 and EFR, the receptors that recognize the peptide MAMPs flg22 and elf18, respectively [19, 20]. This clade also contains the A. thaliana ortholog of the rice LRR-RLK XA21 (based on reciprocal BLAST analysis), which recognizes a sulphated peptide from the RaxX protein to confer resistance to Xanthomonas oryzae pv. oryzae [21, 22, 24]. In addition to the phylogenetic proximity, the presence of a predicted 23 LRR domains (data not shown) in the XPS1 extracellular domain suggests the protein acts as a receptor. It has previously been suggested  that the number of LRR domains present in a given LRR-RLK can be used to predict whether that protein functions as a co-receptor, such as BAK1 with only five LRR domains [36, 37], or receptor, such as FLS2 or EFR which have 28 and 21 LRR domains, respectively [19, 20]. These observations suggest that the clade of LRR-RLKs containing FLS2, EFR, and XPS1 may be involved generally in MAMP perception and provides several excellent candidates for further inquiry.
There are few previous examples of studies that quantitate the binding affinity between a PRR and its peptide MAMP. While the techniques used in those reports differ from those used here, they generally report IC50 values in the low nM range [18, 20, 38]. While not directly comparable, the apparent EC50 of XPS1 for xup25 is in the high nM range. Our attempts to fit our data using a dissociation constant (KD) model showed highly significant variation between the two assays we performed (using independent protein preparations) even at different MST powers. Therefore, we were not able to calculate an accurate overall KD for the XPS1-xup25 binding. We attribute this variation in KD to non-optimized binding events driven by the heterogeneous N-glycosylation of the XPS1 ectodomain in our expression system (as observed in coomassie stained samples of the purified protein (Fig. 8a)). Alternatively, the variations in binding constants might reflect the use of a non-optimized xup25 peptide sequence in our binding assays. In this respect, we anticipate that interaction assays involving shorter or longer variants of xup25 will be very informative. Additionally, while we definitively show that the XPS1 binds to xup25 and is required for its biological activity, this does not eliminate the possibility that the peptide binds to other cell surface molecules including co-receptors. A missing required co-receptor in these in vitro binding assays would also explain the observed variability and current efforts in the lab are focused on discovery of other members of the binding complex.
To date, the identification of PRRs has relied primarily on forward genetic screens, with limited examples of successful reverse genetics approaches [14, 20, 39]. While this forward genetics approach has been very successful, it is limited when applied to MAMPs with weak biological effects or less robust phenotypes. We instead employed a reverse genetics approach to screen the predicted elicitors against a collection of predicted A. thaliana immune receptors. This method required a collection of LRR-RLK knockout lines and the use of a high-throughput assay to measure immune elicitation (the POX assay), but enabled us to screen a wide diversity of candidate immune receptors. While most LRR-RLK T-DNA lines showed a wild-type response to all tested peptides, some displayed general hyper- or hypo-sensitivity to the suite of elicitors, while a smaller group (including xps1) showed specific loss of sensitivity to a single peptide; thereby providing a wealth of leads for the dissection of plant immune signaling pathways.
As genomic and transcriptomic data become more available, these bioinformatics approaches to MAMP and PRR identification will become even more powerful. The continued elucidation of the true range of MAMP epitopes and their cognate receptors will not only benefit plant health, but also increase our understanding of the evolution of the plant-pathogen interaction.
Here we expand on our previous work, utilizing a predictive bioinformatics approach to identify potential MAMPs based on their unique evolutionary signatures, and then use these immune elicitors in a reverse genetic screen to identify their putative receptors or co-receptors. We show that the xup25 peptide derived from a xanthine/uracil permease protein is a novel MAMP with the associated biological effects, including plant growth inhibition, induction of defence gene expression and callose deposition, and the stimulation of a functional immune response sufficient to suppress the growth of a pathogen. Most significantly, xup25 specifically requires the function of the LRR-RLK XPS1 to elicit all of these responses and directly binds to this plant receptor. While we have focused on one particular interaction to validate the methodology, both the computational MAMP screen and the subsequent reverse genetic receptor screen have provided numerous candidates for future study. Of particular interest are those LRR-RLKs that show consistent patterns of hypo- and hyper-sensitivity to multiple elicitors. These proteins may play key roles in regulating the A. thaliana basal immune response, and therefore, are potentially important targets for the development of resistant crops.
Plant materials and growth conditions
All plants were grown with a 12-h photoperiod (100–150 μmol/m2s) at 21 °C, followed by a 12-h dark period at 20 °C in a Conviron growth room (Conviron). All T-DNA lines were selected using the T-DNA Express database of the Salk Institute Genomic Analysis Laboratory (SIGnAL, http://signal.salk.edu/cgi-bin/tdnaexpress). Preference was given to lines with insertions in coding sequences. The genotypes of lines of interest were confirmed by PCR using primers designed using the T-DNA primer design tool (http://signal.salk.edu/tdnaprimers.2.html).
The analysis was performed using protein sequences from the genomes of 54 P. syringae strains. The dataset included genomes from 28 strains available in public databases and 26 strains sequenced by the University of Toronto Centre for Analysis of Genome Evolution and Function (Additional file 1: Table S1). The ortholog prediction was performed using the OrthoMCL program with protein blast e-value 1e-10 . The protein sequences were predicted as ortholog pairs if they shared at least 70 % identity over 70 % sequence length. OrthoMCL was used to cluster ortholog pairs into the ortholog protein families using a markov chain clustering algorithm. The ortholog families that were present in at least 90 % of the strains were defined as core gene families.
The codon-based multiple alignments were constructed with translated protein sequences using translatorX and the MUSCLE multiple alignment program [41, 42]. The protein-based phylogenetic tree was built for each core gene family using the FastTree program . The patterns of natural selection were detected in the predicted core gene families by using the codeml program in the PAML application suite . The prediction of positive selection in the core genes was performed by applying codeml random-site models M7 (neutral selection) and M8 (positive selection) that allow the ω ratio dN/dS to vary among sites. A likelihood ratio test was used to reject the neutral and negative selection null model (M7) in favor of positive selection model (M8) with the chi-square P value <0.01 for all candidate core genes. Empirical Bayes estimates from the model M8 giving a posterior probability >0.5 were used to identify positively selected codons.
The predicted peptide MAMPs were designed from the P. syringae Pto DC3000 genome sequence as 25 amino acid sequences. The flg22 sequence is taken from the P. aeruginosa genome sequence as a 22 amino acid sequence. All peptides were purchased from the Sheldon Biotechnology Centre at McGill University and GenScript.
Single leaves were excised from six plants per genotype of mature A. thaliana plants. From each leaf ten size one leaf cores were taken and washed for 1 h in 1 mL of 1X MS solution with agitation. After washing, leaves were transferred to individual wells of a clear 96-well assay plate avoiding the use of the edge wells to minimize evaporation effects. Each well received 50 μL of 1X MS buffer alone or supplemented with 1 μM of each peptide. Thus each leaf was tested with each treatment, allowing for paired statistical testing. Plates were sealed with parafilm and incubated for 20 h with agitation. The leaf disks were removed and each well received 50 μL of a 1 mg/mL solution of 5-aminosalicylic acid (A79809, Sigma-Aldrich) pH 6.0 with 0.01 % hydrogen peroxide. The reaction was allowed to proceed for 1–3 min and stopped by the addition of 20 μM 2 N NaOH prior to reading the OD600 on a POLARstar OPTIMA microplate reader (BMG Labtech).
LRR-RLK prediction and alignment
RLKs were identified by identifying all proteins with both an LRR domain similar to the LRR domain of either FLS2 or EFR and a kinase domain similar to the kinase domain of FLS2. Similarity searches were performed using both Delta BLAST  and PSI-BLAST  with the respective conserved domain. The A. thaliana RefSeq database was queried using a minimum E-value threshold of 1e-08, and the searches were iterated until no new sequences were added. RLK candidate proteins were selected for further analysis only if they had a match to both an LRR domain and a kinase domain, resulting in 227 candidates. Global multiple sequence alignment was performed with the MAFFT E-INS-i algorithm  and a maximum likelihood phylogeny was constructed in MEGA v5  using partial deletion and 25 % site coverage cutoff, the WAG with Freqs. (+F) substitution mode, gamma distribution with five rate categories, and 100 bootstraps.
Confirmation of T-DNA insertion lines
The SALK lines SALK_101647 (xps1-1) and SALK_101668 (xps1-2) were obtained from the Arabidopsis Biological Resource Center (ABRC) . RT-PCR analyses were performed to verify alterations in transcript levels. Plant RNA was extracted from leaves, which were excised, pooled, and flash frozen in liquid nitrogen prior to homogenization by mortar and pestle under liquid nitrogen. Total RNA was extracted using TRI Reagent as per manufacturer’s instructions (Sigma-Aldrich), followed by DNaseI treatment (Thermo Scientific). A total of 2 μg of RNA was used for cDNA synthesis according to the manufacturer’s instructions with the SuperScript II reverse transcriptase (Life Technologies). Amplification of UBQ10 served as an internal control in RT-PCR assays. The PCR primers used were: XPS1 LRR domain 5’- TCTTCACTAATATTCCTG-3’ and 5’- TGATAGTTTATTGTATGA-3’, UBQ10 (At4g05320) 5’-CACACTCCACTTGGTCTTGCGT-3’ and 5’-TGGTCTTTCCGGTGAGAGTCTTCA-3’. The PCR included 35 cycles for XPS1 LRR domain, and 32 cycles for UBQ10.
Real-time PCR assay
A single leaf on two individual plants was pressure infiltrated with water alone or supplemented with 10 μM peptide. After 24 h of treatment the leaves were excised, pooled, and flash frozen in liquid nitrogen prior to homogenization by mortar and pestle under liquid nitrogen. Total RNA was extracted using TRI Reagent as per manufacturer’s instructions (Sigma-Aldrich), followed by DNaseI treatment (Thermo Scientific). A total of 2 μg of RNA was used for cDNA synthesis according to the manufacturer’s instructions with the SuperScript II reverse transcriptase (Life Technologies). The PCR was carried out in a final reaction volume of 20 μL using Maxima SYBR Green/ROX qPCR Master Mix (Thermo Scientific) using the recommended conditions on a BioRad C1000 thermal cycler (BioRad). The PCR primers used were: PR1 (At2g14610) 5’-CGGAGCTACGCAGAACAACT-3’ and 5’-CTCGCTAACCCACATGTTCA-3’, UBQ10 (At4g05320) 5’-CACACTCCACTTGGTCTTGCGT-3’ and 5’-TGGTCTTTCCGGTGAGAGTCTTCA-3’. The PCR program used for all reactions was: 95 °C for 10 min, followed by 40 cycles consisting of 95 °C for 15 s, 55 °C for 30 s, and 72 °C for 30 s. Each sample was assayed in triplicate and data were evaluated using BioRad CFX manager software, version 1.6 (BioRad). The relative PR1 expression was normalized using the UBQ10 control and reported as fold increase in expression relative to the water treated control using the ΔΔCt method .
A. thaliana plants were pressure infiltrated with water alone, or water containing 10 μM or 100 nM of peptide. The leaves were harvested after 24 h of treatment and cleared with lactophenol. The callose deposits were stained with aniline blue and imaged by epifluorescent microscopy. The images were analyzed with ImageJ  software and the percentage of callose coverage in each image was calculated excluding staining of trichomes and the vasculature. A minimum of two images were taken from each of two leaves from two plants per treatment. The experiment was replicated four times with similar results.
Seedling growth inhibition assay
Seedlings of A. thaliana ecotype Col-0, xps1-1, or xps1-2 were grown for 5 days on MS-Agar plates prior to transfer of two seedlings to each well of a 24-well plate containing 400 μL of 0.5× MS medium with 1 % sucrose. The seedlings were treated with water or the noted dose of peptide and grown for a further 10–14 days. The seedlings were removed, briefly dried, and weighed (fresh weight). A minimum of four plants per genotype per treatment were included, and the experiment was replicated three times with similar results; all data are presented. For each replicate the mean fresh weight of the water control seedlings was calculated and used to normalize the weights of all seedlings of that genotype as a percentage of the control value.
Virulence suppression assay
A. thaliana plants were pressure infiltrated with water alone or water containing 10 μM of peptide 24 h prior to infection with Pto DC3000 (0.0005 OD600). Bacterial growth was measured after 24 h. Briefly, surface sterilized leaf disks were homogenized and the resulting culture was serially diluted and plated on solid selective media. The resulting colony counts were used to determine the number of colony forming units per surface area of the leaf. Each treatment was conducted on four leaves per plant of six individual plants. The experiment was repeated three times with similar results.
XPS1 extracellular domain expression and purification
For XPS1 expression in Hi5 cells, the ectodomain of XPS1 (36-650) was inserted into the baculovirus transfer vector pMelBac B (Invitrogen). The extracellular domain of XPS1 (amino acids 36-650) was cloned by ligation independent cloning between the existing Honey bee melettin signal sequence and the C-terminal Strep II-9x His tag. All clones were verified by Sanger Sequencing. Primers (Forward: GGT CGT ATA CAT TTC TTA CAT CTA TGC GAC GGA GGA GAC TGA TAA ACA AGC ATT GC; Reverse: GCA CCC TGG AAG TAC AGG TTC TCT GAT GAA TGC CTT CTT GGT AAT TCA ACA GAG C) were designed to have an amplification part homologous to the desired boundaries of XPS1 ectodomain and extensions for RecA-mediated SLIC strategy. The StrepII-9x His fused XPS1 ectodomain was produced by secreted expression in baculovirus-infected insect cells, harvested 72 h post infection and purified by Ni-NTA affinity chromatography (Qiagen). The samples were processed twice with a Superdex 200 16/60 column (GE Healthcare) pre-equilibrated with 10 mM Bis-Tris pH 6.0, 150 mM NaCl. Protein purity was checked by SDS-PAGE. The identity of the XPS1 ectodomain was further confirmed by anti-His immunoblots and mass spectrometry.
MST binding assays
The XPS1 ectodomain was kept at a constant concentration (0.5 μM) in a buffer containing 10 mM Bis-Tris pH6, 150 mM NaCl, 5 % glycerol and 0.01 % Tween, whereas varying peptide concentrations were added. Approximately 4–6 μL of each sample was loaded in a fused silica capillary (NanoTemper Technologies). Measurements were performed at room temperature in a Monolith NT.label free instrument at a constant LED power of 20 and MST power of 60 %. Measurements were performed repeatedly on two independent protein preparations to ensure reproducibility. The data were analyzed by plotting xup25 concentrations against percent changes of normalized fluorescence (ΔFnorm [%] y-axis). Curve fitting was performed by using the NT affinity analysis software from Nanotemper.
Luminol-based ROS generation assay
Measurement of ROS production was performed using a modified version of the luminol-based assay as described . Ten leaf discs (diameter 4 mm) were taken from the leaves of 4-week-old Col-0 A. thaliana. Leaf disks were washed with 200 μL of sterile water for 20 h. The water was removed and replaced with 100 mM Tris-HCl (pH 8.0) containing 20 μg/mL horseradish peroxidase and 34 μg/mL luminol. Wells were treated with water control or the noted dose of peptide. The luminescence was summed over a 2-s interval every 2 min for a period of 60 min on an Infinite M1000Pro microplate reader (TECAN Group Ltd.). Each treatment was performed on six wells.
Sequence alignment and sequence logo
The protein sequence of the xup25 containing xanthine uracil permease gene from each species used was extracted and aligned using MUSCLE . The alignment was uploaded to the Gene Slider module of the Bio-Analytical Resource for Plant Biology (bar.utoronto.ca) to produce the sequence logo.
No ethics approval was necessary for this study.
All of the genomic data produced for this study have been submitted to NCBI; the BioProject Accession numbers for these genomes along with all the publicly available genomes used can be found in Additional file 1: Table S1. The RefSeq ID numbers for all LRR-RLK sequences used can be found in Additional file 1: Table S5.
We would like to thank the members of the DSG and DD labs for helpful discussion, and in particular Sarah Cao and Alina Guna for assistance in plant maintenance. We thank the ABRC for T-DNA insertion lines used in this study. We thank the VBCF Protein Technologies Facility for excellent assistance.
This work was supported by Natural Sciences and Engineering Research Council of Canada Discovery Awards to DSG and DD; a Canada Research Chair in Plant-Microbe Systems Biology (DD) or Comparative Genomics (DSG); the Centre for the Analysis of Genome Evolution and Function (DD and DSG). This work was also supported by funds from the Austrian Academy of Science through the Gregor Mendel Institute (YB).
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Dodds PN, Rathjen JP. Plant immunity: towards an integrated view of plant-pathogen interactions. Nat Rev Genet. 2010;11:539–48.View ArticlePubMedGoogle Scholar
- Mott GA, Middleton MA, Desveaux D, Guttman DS. Peptides and small molecules of the plant-pathogen apoplastic arena. Front Plant Sci. 2014;5:677.View ArticlePubMedPubMed CentralGoogle Scholar
- McCann HC, Nahal H, Thakur S, Guttman DS. Identification of innate immunity elicitors using molecular signatures of natural selection. Proc Natl Acad Sci U S A. 2012;109:4215–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun W, Dunning FM, Pfund C, Weingarten R, Bent AF. Within-species flagellin polymorphism in Xanthomonas campestris pv campestris and its impact on elicitation of Arabidopsis FLAGELLIN SENSING2-dependent defenses. Plant Cell. 2006;18:764–79.View ArticlePubMedPubMed CentralGoogle Scholar
- Felix G, Duran JD, Volko S, Boller T. Plants have a sensitive perception system for the most conserved domain of bacterial flagellin. Plant J. 1999;18:265–76.View ArticlePubMedGoogle Scholar
- Kunze G, Zipfel C, Robatzek S, Niehaus K, Boller T, Felix G. The N terminus of bacterial elongation factor Tu elicits innate immunity in Arabidopsis plants. Plant Cell. 2004;16:3496–507.View ArticlePubMedPubMed CentralGoogle Scholar
- Monaghan J, Zipfel C. Plant pattern recognition receptor complexes at the plasma membrane. Curr Opin Plant Biol. 2012;15:349–57.View ArticlePubMedGoogle Scholar
- Greeff C, Roux M, Mundy J, Petersen M. Receptor-like kinase complexes in plant innate immunity. Front Plant Sci. 2012;3:209.PubMedPubMed CentralGoogle Scholar
- Shiu S-H, Bleecker AB. Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. Proc Natl Acad Sci U S A. 2001;98:10763–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Wang G, Ellendorff U, Kemp B, Mansfield JW, Forsyth A, Mitchell K, et al. A genome-wide functional investigation into the roles of receptor-like proteins in Arabidopsis. Plant Physiol. 2008;147:503–17.View ArticlePubMedPubMed CentralGoogle Scholar
- Gou X, He K, Yang H, Yuan T, Lin H, Clouse SD, et al. Genome-wide cloning and sequence analysis of leucine-rich repeat receptor-like protein kinase genes in Arabidopsis thaliana. BMC Genomics. 2010;11:19.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun Y, Li L, Macho AP, Han Z, Hu Z, Zipfel C, et al. Structural basis for flg22-induced activation of the Arabidopsis FLS2-BAK1 immune complex. Science. 2013;342:624–8.View ArticlePubMedGoogle Scholar
- Schulze B, Mentzel T, Jehle AK, Mueller K, Beeler S, Boller T, et al. Rapid heteromerization and phosphorylation of ligand-activated plant transmembrane receptors and their associated kinase BAK1. J Biol Chem. 2010;285:9444–51.View ArticlePubMedPubMed CentralGoogle Scholar
- Hou S, Wang X, Chen D, Yang X, Wang M, Turrà D, et al. The secreted peptide PIP1 amplifies immunity through receptor-like kinase 7. PLoS Pathog. 2014;10:e1004331.View ArticlePubMedPubMed CentralGoogle Scholar
- Minibayeva F, Beckett RP, Kranner I. Roles of apoplastic peroxidases in plant response to wounding. Phytochemistry. 2015;112:122–9.View ArticlePubMedGoogle Scholar
- Daudi A, Cheng Z, O’Brien JA, Mammarella N, Khan S, Ausubel FM, et al. The apoplastic oxidative burst peroxidase in Arabidopsis is a major component of pattern-triggered immunity. Plant Cell. 2012;24:275–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Mammarella ND, Cheng Z, Fu ZQ, Daudi A, Bolwell GP, Dong X, et al. Apoplastic peroxidases are required for salicylic acid-mediated defense against Pseudomonas syringae. Phytochemistry. 2015;112:110–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Chinchilla D, Bauer Z, Regenass M, Boller T, Felix G. The Arabidopsis receptor kinase FLS2 binds flg22 and determines the specificity of flagellin perception. Plant Cell. 2006;18:465–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Gómez-Gómez 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–11.View ArticlePubMedGoogle Scholar
- Zipfel C, Kunze G, Chinchilla D, Caniard A, Jones JDG, Boller T, et al. Perception of the bacterial PAMP EF-Tu by the receptor EFR restricts Agrobacterium-mediated transformation. Cell. 2006;125:749–60.View ArticlePubMedGoogle Scholar
- Song WY, Wang GL, Chen LL, Kim HS, Pi LY, Holsten T, et al. A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. Science. 1995;270:1804–6.View ArticlePubMedGoogle Scholar
- Holton N, Nekrasov V, Ronald PC, Zipfel C. The phylogenetically-related pattern recognition receptors EFR and XA21 recruit similar immune signaling components in monocots and dicots. PLoS Pathog. 2015;11:e1004602.View ArticlePubMedPubMed CentralGoogle Scholar
- Bahar O, Pruitt R, Luu DD, Schwessinger B, Daudi A, Liu F, et al. The Xanthomonas Ax21 protein is processed by the general secretory system and is secreted in association with outer membrane vesicles. PeerJ. 2014;2:e242.View ArticlePubMedPubMed CentralGoogle Scholar
- Pruitt RN, Schwessinger B, Joe A, Thomas N, Liu F, Albert M, et al. The rice immune receptor XA21 recognizes a tyrosine-sulfated protein from a Gram-negative bacterium. Sci Adv. 2015;1:e1500245–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Gómez-Gómez L, Felix G, Boller T. A single locus determines sensitivity to bacterial flagellin in Arabidopsis thaliana. Plant J. 1999;18:277–84.View ArticlePubMedGoogle Scholar
- Asai T, Tena G, Plotnikova J, Willmann MR, Chiu W-L, Gomez-Gomez L, et al. MAP kinase signalling cascade in Arabidopsis innate immunity. Nature. 2002;415:977–83.View ArticlePubMedGoogle Scholar
- Zipfel C, Robatzek S, Navarro L, Oakeley EJ, Jones JDG, Felix G, et al. Bacterial disease resistance in Arabidopsis through flagellin perception. Nature. 2004;428:764–7.View ArticlePubMedGoogle Scholar
- Seidel SAI, Dijkman PM, Lea WA, van den Bogaart G, Jerabek-Willemsen M, Lazic A, et al. Microscale thermophoresis quantifies biomolecular interactions under previously challenging conditions. Methods. 2013;59:301–15.View ArticlePubMedPubMed CentralGoogle Scholar
- Clarke CR, Chinchilla D, Hind SR, Taguchi F, Miki R, Ichinose Y, et al. 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–60.View ArticlePubMedPubMed CentralGoogle Scholar
- Naito K, Taguchi F, Suzuki T, Inagaki Y, Toyoda K, Shiraishi T, et al. Amino acid sequence of bacterial microbe-associated molecular pattern flg22 is required for virulence. Mol Plant Microbe Interact. 2008;21:1165–74.View ArticlePubMedGoogle Scholar
- Nekrasov V, Li J, Batoux M, Roux M, Chu Z-H, Lacombe S, et al. Control of the pattern-recognition receptor EFR by an ER protein complex in plant immunity. EMBO J. 2009;28:3428–38.View ArticlePubMedPubMed CentralGoogle Scholar
- Cai R, Lewis J, Yan S, Liu H, Clarke CR, Campanile F, et al. The plant pathogen Pseudomonas syringae pv. tomato is genetically monomorphic and under strong selection to evade tomato immunity. PLoS Pathog. 2011;7:e1002130.View ArticlePubMedPubMed CentralGoogle Scholar
- O’Brien HE, Thakur S, Guttman DS. Evolution of plant pathogenesis in pseudomonas syringae: a genomics perspective. Annu Rev Phytopathol. 2011;49:269–89.View ArticlePubMedGoogle Scholar
- O’Brien HE, Thakur S, Gong Y, Fung P, Zhang J, Yuan L, et al. Extensive remodeling of the Pseudomonas syringae pv. avellanae type III secretome associated with two independent host shifts onto hazelnut. BMC Microbiol. 2012;12:141.View ArticlePubMedPubMed CentralGoogle Scholar
- Belkhadir Y, Yang L, Hetzel J, Dangl JL, Chory J. The growth-defense pivot: crisis management in plants mediated by LRR-RK surface receptors. Trends Biochem Sci. 2014;39:447–56.View ArticlePubMedPubMed CentralGoogle Scholar
- Li J, Wen J, Lease KA, Doke JT, Tax FE, Walker JC. BAK1, an Arabidopsis LRR receptor-like protein kinase, interacts with BRI1 and modulates Brassinosteroid signaling. Cell. 2002;110:213–22.View ArticlePubMedGoogle Scholar
- Nam KH, Li J. BRI1/BAK1, a receptor kinase pair mediating Brassinosteroid signaling. Cell. 2002;110:203–12.View ArticlePubMedGoogle Scholar
- Bauer Z, Gómez-Gómez L, Boller T, Felix G. Sensitivity of different ecotypes and mutants of Arabidopsis thaliana toward the bacterial elicitor flagellin correlates with the presence of receptor-binding sites. J Biol Chem. 2001;276:45669–76.View ArticlePubMedGoogle Scholar
- Danna CH, Millet YA, Koller T, Han S, Bent AF, Ronald PC, et al. The Arabidopsis flagellin receptor FLS2 mediates the perception of Xanthomonas Ax21 secreted peptides. Proc Natl Acad Sci U S A. 2011;108:9286–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Li L, Stoeckert CJ, Roos DS. OrthoMCL: identification of ortholog groups for eukaryotic genomes. Genome Res. 2003;13:2178–89.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. MUSCLE: a multiple sequence alignment method with reduced time and space complexity. BMC Bioinformatics. 2004;5:113.View ArticlePubMedPubMed CentralGoogle Scholar
- Abascal F, Zardoya R, Telford MJ. TranslatorX: multiple alignment of nucleotide sequences guided by amino acid translations. Nucleic Acids Res. 2010;38:W7–13.View ArticlePubMedPubMed CentralGoogle Scholar
- Price MN, Dehal PS, Arkin AP. FastTree: computing large minimum evolution trees with profiles instead of a distance matrix. Mol Biol Evol. 2009;26:1641–50.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007;24:1586–91.View ArticlePubMedGoogle Scholar
- Boratyn GM, Schäffer AA, Agarwala R, Altschul SF, Lipman DJ, Madden TL. Domain enhanced lookup time accelerated BLAST. Biol Direct. 2012;7:12.View ArticlePubMedPubMed CentralGoogle Scholar
- Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.View ArticlePubMedPubMed CentralGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33:511–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Alonso JM, Stepanova AN, Leisse TJ, Kim CJ, Chen H, Shinn P, et al. Genome-wide insertional mutagenesis of Arabidopsis thaliana. Science. 2003;301:653–7.View ArticlePubMedGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.View ArticlePubMedGoogle Scholar