# New insights on leucine-rich repeats receptor-like kinase orthologous relationships in angiosperms

Jean-François Dufayard, Mathilde Bettembourg, Iris Fischer, Gaëtan Droc,

Emmanuel Guiderdoni, Christophe Perin, Nathalie Chantret, Anne Dievart

► To cite this version:

Jean-François Dufayard, Mathilde Bettembourg, Iris Fischer, Gaëtan Droc, Emmanuel Guiderdoni, et al. New insights on leucine-rich repeats receptor-like kinase orthologous relationships in angiosperms. Frontiers in Plant Science, 2017, 10.3389/fpls.2017.00381. hal-01605653

HAL Id: hal-01605653

https://hal.science/hal-01605653

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
# NEW INSIGHTS ON LEUCINE-RICH REPEATS RECEPTOR-LIKE KINASE ORTHOLOGOUS RELATIONSHIPS IN ANGIOSPERMS

Jean-François Dufayard¹, Mathilde Bettembourg¹, Iris Fischer², Gaetan Droc¹, Emmanuel Guiderdoni¹, Christophe Perin¹, Nathalie Chantret²*, Anne Dievart¹*

¹UMR AGAP, CIRAD, France, ²UMR AGAP, INRA, France

**Submitted to Journal:**
Frontiers in Plant Science

**Specialty Section:**
Plant Genetics and Genomics

**ISSN:**
1664-462X

**Article type:**
Original Research Article

**Received on:**
21 Nov 2016

**Accepted on:**
06 Mar 2017

**Provisional PDF published on:**
06 Mar 2017

**Frontiers website link:**
www.frontiersin.org

**Citation:**
Dufayard J, Bettembourg M, Fischer I, Droc G, Guiderdoni E, Perin C, Chantret N and Dievart A(2017) NEW INSIGHTS ON LEUCINE-RICH REPEATS RECEPTOR-LIKE KINASE ORTHOLOGOUS RELATIONSHIPS IN ANGIOSPERMS. *Front. Plant Sci.* 8:381. doi:10.3389/fpls.2017.00381

**Copyright statement:**
© 2017 Dufayard, Bettembourg, Fischer, Droc, Guiderdoni, Perin, Chantret and Dievart. This is an open-access article distributed under the terms of the [Creative Commons Attribution License (CC BY)](https://creativecommons.org/licenses/by/4.0/). The use, distribution and reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
1

2

3

4

5

6

7

8

9

10

11

12

13

14

15

16

NEW INSIGHTS ON LEUCINE-RICH REPEATS RECEPTOR-LIKE KINASE ORTHOLOGOUS RELATIONSHIPS IN ANGIOSPERMS

Jean-François Dufayard¹, Mathilde Bettembourg¹, Iris Fischer², Gaetan Droc¹, Emmanuel Guiderdoni¹, Christophe Périn¹, Nathalie Chantret²,* and Anne Diévart¹,*

¹CIRAD, UMR AGAP, F-34398 Montpellier, France

²INRA, UMR AGAP, F-34060 Montpellier, France

Correspondence:

Anne Diévart, CIRAD, UMR AGAP, Avenue Agropolis, TA A 108/03, 34 398 Montpellier Cedex 5, France, +33 4 67 61 75 70, anne.dievart@cirad.fr.

Nathalie Chantret, INRA, UMR AGAP, 2, Place P. Viala, 34060 Montpellier Cedex 1, France, +33 4 99 61 60 83, nathalie.chantret@supagro.inra.fr.
17 **Abbreviations**

18 RLK, Receptor-Like Kinase; LRR, Leucine-Rich Repeat; LRR-RLK, Leucine-Rich Repeats

19 Receptor-Like Kinase; TM, transmembrane domain; KD, kinase domain; SG, subgroup;

20 ECD, extracellular domain; ICD, intracellular domain; OG, orthologous group; Cys-pair,

21 cysteine-pair; MLD, malectin-like domain; SP, signal peptide; MD, monocots dicots.
**Abstract**

*Leucine-Rich Repeats Receptor-Like Kinase (LRR-RLK)* genes represent a large and complex gene family in plants, mainly involved in development and stress responses. These receptors are composed of an LRR-containing extracellular domain (ECD), a transmembrane domain (TM) and an intracellular kinase domain (KD). To provide new perspectives on functional analyses of these genes in model and non-model plant species, we performed a phylogenetic analysis on 8,360 LRR-RLK receptors in 31 angiosperm genomes (8 monocots and 23 dicots). We identified 101 orthologous groups (OGs) of genes being conserved among almost all monocot and dicot species analyzed. We observed that more than 10% of these OGs are absent in the Brassicaceae species studied. We show that the ECD structural features are not always conserved among orthologs, suggesting that functions may have diverged in some OG sets. Moreover, we looked at targets of positive selection footprints in 12 pairs of OGs and noticed that depending on the subgroups, positive selection occurred more frequently either in the ECDs or in the KDs. (167 words)

**Keywords**

LRR, receptor, kinase, angiosperms, phylogeny, orthologs.
# 1. Introduction

Receptor-like kinases (RLK) constitute one of the largest gene families in the plant kingdom. They are typically composed of an amino-terminal extracellular domain (ECD), a transmembrane domain (TM), and an intracellular domain containing the kinase domain (KD). Several phylogenetic studies of the RLK family were conducted, initially focusing on Arabidopsis and later including other plant species (Shiu and Bleecker, 2001b; 2003; Shiu et al., 2004; Lehti-Shiu et al., 2009; Liu et al., 2009; Sakamoto et al., 2012; Zan et al., 2013). Using a phylogeny inferred from their KD alignment the Arabidopsis RLK genes were classified into 44 subgroups (SGs) or subfamilies (Shiu and Bleecker, 2001a). Fifteen SGs have been described containing common motifs in their ECD. The ECD of the largest SG possesses Leucine-Rich Repeats (LRR) and this SG has therefore been named Leucine-Rich Repeats Receptor-Like Kinase (LRR-RLK) (Kobe and Deisenhofer, 1994; Kajava, 1998; Shiu and Bleecker, 2001a; Shiu et al., 2004; Lehti-Shiu et al., 2009). The first members of this large family were cloned in the 90's and their signaling pathways were extensively studied. Those members are ERECTA (ER), CLAVATA1 (CLV1), BRASSINOSTEROID INSENSITIVE 1 (BRI1), SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK), HAESA-RLK5 and Xa21 (Horn and Walker, 1994; Song et al., 1995; Torii et al., 1996; Clark et al., 1997; Li and Chory, 1997; Schmidt et al., 1997). To date, functions have been assigned to ~35% of the ~230 LRR-RLK members in *A. thaliana* and ~ to a lesser extent ~ other species (Wu et al., 2016). They are important mediators of cell-cell communication to relay developmental cues and environmental stimuli or to activate defense/resistance against pathogens (Mu et al., 1994; Muschietti et al., 1998; Antolin-Llovera et al., 2014a; Belkhadir et al., 2014; Jaouannet et al., 2014).

Functional analyses conducted on LRR-RLK genes over the last twenty years raveled the role of the domains located in the ECD of these receptors. The LRR domains are highly versatile in number allowing a whole range of protein-protein interactions. These include homo- or hetero-dimerization of receptors, in addition to ligand binding. Furthermore, some LRR-RLK receptors possess island domains - devoid of LRRs - located between LRR motifs (Li and Chory, 1997). They were identified in the BRI1 receptor as the binding site for the brassinosteroid hormone (Kinoshita et al., 2005; Hothorn et al., 2011; She et al., 2011). Few studies have also described the functions of other ECD domains. For example, two cysteine-pairs (Cys-pair) have been reported. The first one is located in the N-terminal part of the LRR-RLKs, approximately 60 AA from the start codon between the signal peptide (SP) and the first LRRs. The second one - if present - can be found between the last LRR and the TM domain (Dievart and Clark, 2003). Mutations in the Cys-pairs have been shown to affect the function of some LRR-RLKs, e.g. FLAGELLIN SENSING 2 (FLS2), a gene participating in the perception of the bacterial elicitor flagellin. However, there is also an example of a LRR receptor like protein (CLAVATA 2) for which mutations in Cys-pairs had no effect on the function of protein - at least in the meristem and roots (Noguchi et al., 1999; Song et al., 2010; Sun et al., 2012). In BRI1, a mutant harboring a mutation in Cys-pairs appears to be functional but seems to be retained in the endoplasmic reticulum and degraded. This suggests that this mutant protein does not pass the endoplasmic reticulum quality control (Hong et al., 2008). Although no general conclusions can be drawn so far on the importance of this motif, all the variations observed in Cys-pairs likely play a role in the folding, trafficking and/or the binding to other proteins. It was therefore suggested that this motif influences the signaling pathways activated downstream of the LRR-RLKs (Su et al., 2012). Another ECD, the malectin-like domain (MLD) lying in between the SP and the LRRs, has also been described in one LRR-RLK SG (Hok et al., 2011). In legumes and actinorhizal plants, the SYMBIOSIS RECEPTOR LIKE KINASE (SYMRK, also known as NORK or DMI2) receptor, involved in
phosphate-acquiring arbuscular mycorrhiza and in nitrogen-fixing root nodule symbiosis, possesses a malectin domain but the exact function of this receptor is still unclear (Antolin-Llovera et al., 2014a). It has been recently demonstrated that the SYMRK receptor is likely cleaved at the plasma membrane to release the N-glycosylated MLD (Antolin-Llovera et al., 2014b). Moreover, this cleavage would permit a physical interaction between the SYMRK and the LysM-type RLK NOD FACTOR RECEPTOR 5 and induces a rapid degradation of the SYMRK protein lacking its MLD. Thus, all the domains lying in the ECD with the LRRs play essential and complementary roles for specific LRR-RLK receptor functions.

Their central role in plant development and perception of environmental condition or stresses, their ubiquity in all angiosperms, and the complexity of their relationships make LRR-RLK genes an interesting candidate family to be studied in a phylogenetic framework (Shi et al., 2014). Such an analysis will be helpful to identify groups of orthologous genes and to compare functions between orthologs. However, inferring the phylogeny of such a large family raises several challenges. First, the vast number of sequences to be analyzed poses a problem of computational time and space. Second, the high rate of gene gains and losses during the evolution of the family, species-specific characteristics, and annotation errors result in complex orthologous relationships that are not always identified correctly by automatic gene annotation. For these reasons, large gene families - such as LRR-RLKs - are not well characterized on platforms like GreenPhylDB or Phytozome dedicated to automatic clustering (Conte et al., 2008; Rouard et al., 2011; Goodstein et al., 2012) and significant manual expertise is required to produce reliable results.

In the present article, we conducted a phylogenetic analysis of the LRR-RLK genes from 33 plant genomes with the objective to investigate the characteristics of genes belonging to the same orthologous groups, expected to be conserved among most monocot and dicot species analyzed. To do so, we first looked for and identified 101 orthologous groups (OGs) of genes present in most genomes analysed defined them as the LRR-RLK angiosperm "core" sets. We observed that ECD structural features were not always conserved in some OGs, suggesting that functions may have diverged among these orthologs. We also looked at selection footprints that led to the differentiations of pairs of OGs. This allowed us to investigate the putative role and function of uncharacterized genes in recently sequenced genomes from experimentally characterized LRR-RLK genes in model organisms.
## 2. Materials and methods

### 2.1. Plant genomes analysed

The 33 species analyzed represent a broad spectrum of land plants (embryophyta), with one bryophyta genome, *Physcomitrella patens* (PHYPA, moss) (Rensing et al., 2008), the *Selaginella moellendorffii* genome (SELML, spikemoss) (Banks et al., 2011), representative of the lycopodiopsida, and 31 species of magnoliophyta (angiosperms), divided into eight monocot species (*Phoenix dactylifera* (PHODC, date palm) (Al-Dous et al., 2011), *Musa acuminata* (MUSAC, banana) (D'Hont et al., 2012), two subspecies of *Oryza sativa* (rice), *Oryza sativa* ssp. *japonica* (ORYSJ) (Goff et al., 2002) and *Oryza sativa* ssp. *indica* (ORYSI) (Yu et al., 2002; Project, 2005), *Brachypodium distachyon* (BRADI, purple false brome) (Initiative, 2010), *Zea mays* (MAIZE, corn) (Schnable et al., 2009), *Sorghum bicolor* (SORBI) (Paterson et al., 2009) and *Setaria italica* (SETIT, Foxtail millet)) (Bennetzen et al., 2012; Zhang et al., 2012), and 23 dicot species (two *Solanum* species, *Solanum tuberosum* (SOLTU, potato) (Consortium, 2011) and *Solanum lycopersicum* (SOLLC, tomato) (Consortium, 2012), *Vitis vinifera* (VITVI, Grape Vine) (Characterization, 2007), *Lotus japonicus* (LOTJA) (Sato et al., 2008), *Medicago truncatula* (MEDTR, Barrel Medic) (Young et al., 2011), *Glycine max* (GLYMA, soybean) (Schmutz et al., 2010), *Cajanus cajan* (CAJCA, pigeon pea) (Varshney et al., 2012), *Prunus persica* (PRUPE, peach) (Ahmad et al., 2011), *Malus x domestica* (MALDO, apple) (Velasco et al., 2010; Jung et al., 2012), *Ricinus communis* (RICCO, castor oil plant) (Chan et al., 2010), *Jatropha curcas* (JATCU) (Sato et al., 2011), *Manihot esculenta* (MANES, Cassava) (Prochnik et al., 2012), *Populus trichocarpa* (POPTR, black cottonwood) (Tuskan et al., 2006), two *Cucumis* species, *Cucumis sativus* (CUCSA, cucumber) (Huang et al., 2009) and *Cucumis melo* (CUCME, Muskmelon) (Sebastian et al., 2010; Garcia-Mas et al., 2012), *Schrenkiella parvula* (SCHPA (formerly EUTPR)) (Dassanayake et al., 2011), *Eutrema salsugineum* (EUTSA (formerly THEHA)) (Oh et al., 2010), *Brassica rapa* (BRARA) (Wang et al., 2011), two *Arabidopsis* species, *Arabidopsis lyrata* (ARALY) and *Arabidopsis thaliana* (ARATH, thale cress) (Initiative, 2000; Hu et al., 2011), *Carica papaya* (CARPA, papaya) (Ming et al., 2008), *Gossypium raimondii* (GOSRA, cotton) (Wang et al., 2012) and *Theobroma cacao* (THECC, cacao tree) (Argout et al., 2011). Note that the genus *Thellungiella* is now known as *Eutrema*, so the species formerly known as *Thellungiella halophila* (THEHA) is now known as *Eutrema salsugineum* (EUTSA). In addition, the species of *Eutrema* sequenced at the JGI has been determined to be *salsugineum* (a close relative of *halophila*). Therefore, this genome is actually classified as *Eutrema salsugineum*. Since we started our analysis before this change, all sequences related to the *Eutrema salsugineum* genome are annotated "THEHA" or "EUTSA" in our paper. Also, in the course of our study, the name *Eutrema parvula* (EUTPR) has been changed for *Schrenkiella parvula* (SCHPA). So all *Schrenkiella parvula* (SCHPA) sequences are annotated "EUTPR" or "SCHPA". See **Suppl-Table 1** for bibliography and web links for download. Species tree representation (**Suppl Material online**) is based on several studies (Koch et al., 2000; Oh and Potter, 2005; Rensing et al., 2008; Forest and Chase, 2009; Magallcdn and Castillo, 2009; Wang et al., 2009; Oh et al., 2010; Sebastian et al., 2010; Arakaki et al., 2011; Reineke et al., 2011; Wang et al., 2011; Young et al., 2011; Consortium, 2012; Varshney et al., 2012; Wang et al., 2012; Zhang et al., 2012).

### 2.2. LRR-RLKs extraction, phylogeny and orthologous groups

On each of the 33 plant proteomes, the hmmsearch program was run to extract peptide sequences containing both LRR(s) and a KD (Eddy, 2009). Sequences containing both LRRs and KD were classified into SGs using a global phylogenetic analysis (Fischer et al., 2016). First, the KD of all these sequences was aligned using MAFFT with a progressive strategy
175 (Katoh et al., 2002). Then the alignment was cleaned with TrimAl configured to remove
176 every sites with more than 20% of gaps or with a similarity score lower than 0.001 (Capella-
177 Gutierrez et al., 2009). A similarity matrix was computed using ProtDist with a JTT model,
178 and then a global distance phylogeny was inferred using FastME configured with default
179 settings and SPR movements to optimize the tree topology (Felsenstein, 1989; Desper and
180 Gascuel, 2002). SGs were defined manually in the global phylogeny using the Arabidopsis
181 genes as reference (Shiu et al., 2004; Lehti-Shiu et al., 2009; Fischer et al., 2016). To extend
182 this dataset to receptor kinases devoid of LRRs in their ECD (sequences annotated
183 "No_LRR"), the BLASTP algorithm (default parameters) has been run SG per SG, using each
184 of the 7,767 KD sequences to search a database composed of the 33 proteomes as query
185 (Altschul et al., 1997). Blast outputs were parsed to keep only homolog sequences sharing
186 more than 90% identity with the query sequence. The new "No_LRR" sequences retrieved by
187 blast were assigned to the same SG as the query sequence. Then, phylogenies were inferred
188 for each of the 20 SGs. Each group of sequences was aligned using MAFFT with an iterative
189 strategy (maximum of 100 iterations) (Katoh et al., 2002). Alignments were cleaned using
190 TrimAl configured this time to remove sites with more than 80% of gaps (Capella-Gutierrez
191 et al., 2009). Then maximum likelihood phylogenies were inferred using PhyML 3.0,
192 configured with LG+gamma model, and the best of NNI and SPR topology optimization
193 (Guindon and Gascuel, 2003). Statistical branch supports were computed using the aLRT/SH-
194 like strategy (Guindon et al., 2010; Anisimova et al., 2011). Each of the 20 phylogenetic trees
195 has been reconciled with the species tree using RAP-Green (Dufayard et al., 2005)
196 (https://github.com/ SouthGreenPlatform/rap-green). By comparing the gene tree with the
197 species tree, this analysis allows us to root phylogenetic trees (Dufayard et al., 2005). We
198 tested this approach of rooting (by minimizing the number of inferred duplications and losses)
199 and compared it with rooting with outgroups (data not shown). The two methods provided
200 very close root locations that did not change the overall conclusions.
201

201

202 To define the MD OGs, monocots / dicots bifurcations (branch support threshold >
203 0.85) have been manually located in each of the 20 SG-specific trees. To be considered as
204 MD OGs, the minimum number of monocots and dicots species represented was 3 and 4
205 respectively, to avoid keeping groups of misannotated sequences as MD OGs. Thus, 101 MD
206 OGs have been defined. For each of them, the number of sequences in each of the 31 studied
207 angiosperm species was recorded. If no sequence was discovered in one species, it was
208 considered lost in this species.
209

209

210 **2.3. Structural features**

211 Number of motifs and positions (LRRs (PF00560.24) and KDs (PF00069.16)) are
212 outputs of the hmmsearch program (default parameters) (Eddy, 2009). Island domains were
213 determined based on predicted LRR positions in sequences. For SP and TM domains, the
214 TMHMM and TOPPRED softwares have been used (Claros and von Heijne, 1994; Krogh et
215 al., 2001). For the malectin domains in SG_I and SG_VIII-2, sequences of the domains were
216 extracted in SMART and aligned to build hmm motifs with the hmmbuild program (Eddy,
217 2009; Letunic et al., 2009). For Cys-pairs, hmm motifs were built based on subsets of
218 sequences known to possess these motifs (Eddy, 2009).
219

219

220 **2.4. Test of positive selection on ancestral branches**

221 Twelve sub-trees were considered: the OGs selected were those with a ‘simple’
222 organization, *i.e.* with a gene topology fitting approximately with the species tree. Three to
223 four sequences among those that are the most closely related to the OG were selected from the
224 whole SG tree as outgroup (**Suppl-Figure 4**). The sequences were re-aligned and the
225 alignment was cleaned as described previously (Fischer et al., 2014; Fischer et al., 2016). We
226 ran codeml branch/site models implemented in the PAML4 software (Yang, 2007). For each
227 OG, the following branch partition was defined: all branches but one were tagged as
228 ‘background’ branches and the branch between the duplication node and the node
229 corresponding to the split between monocot and dicot tagged as ‘foreground’ branch. Then
230 two models were compared: the null model (A₀), in which sites on the fore- and background
231 branches evolved under the same selective pressure (purifying or neutral), and a model
232 including positive selection (model A) in which some sites on the foreground branch evolved
233 under positive selection whereas sites on the background branches still evolved under
234 purifying selection or neutrality. The most likely model was inferred by a likelihood ratio test
235 (LRT). To take into account multiple testing, a Bonferroni correction was applied: the
236 significance threshold of 0.05 was divided by the number of tested branches (24). Sites
237 detected to be under positive selection at the codon level were manually curated for alignment
238 quality and reliability. In branches identified to have evolved under positive selection, Bayes
239 empirical Bayes was used to calculate the posterior probabilities at each codon and detect
240 those under positive selection (i.e. those with a posterior probability of ω>1 strictly above
241 95%).
242
# 3. Results and discussion

## 3.1. More than 200 LRR-RLK genes on average per angiosperm species

We conducted a phylogenetic analysis of the LRR-RLK gene family in 33 fully sequenced plant genomes to classify them into SGs and to highlight and describe general characteristics of these LRR-RLK gene sets (Fischer et al., 2016). Briefly, besides 31 angiosperms genomes - represented by eight monocots (including six poaceae) and 23 dicots - one bryophyte genome of *Physcomitrella patens* (PHYPA, moss) and one lycopodiopsida genome of *Selaginella moellendorffii* (SELML, spikemoss) were included (**Suppl Table 1**, see Materials and Methods for details and five-digit species code). As it has been done previously, we based our classification of LRR-RLK genes into SGs on the KD phylogeny (Shiu et al., 2004; Lehti-Shiu et al., 2009; Fischer et al., 2016). In our previous study (Fischer et al. 2016), the LRR-RLK dataset contained 7,767 sequences possessing at least one LRR in their ECD. Since we focus on structural features of ECDs and presence/absence of the genes in LRR-RLK OGs in this present study, we included LRR-RLK homologues for which LRRs were completely lost or were degenerated in this new dataset. Thus, 593 sequences (prefixed "No_LRR") were added to the original set of 7,767 sequences which lead to a total of 8,360 sequences (**Suppl Table 2**). Within each of the 20 SGs, KDs were aligned and SG-specific trees were obtained using a likelihood-based method (PHYML) (**Figure 1A** and **Suppl Material online**). Altogether, these LRR-RLK genes represent on average 0.71% and 0.66% of the monocot and dicot proteomes, respectively. Interestingly, in moss (PHYPA) and spikemoss (SELML), the proportions of LRR-RLKs per genome (0.41 % and 0.36%, respectively) are approximately half the ratio observed in angiosperms. Likewise, the average number of LRR-RLK genes in angiosperms is 263.6, with 260.7 LRR-RLK proteins (+/- 20.2 (SE)) in dicots and 268.8 (+/- 18.2 (SE)) in monocots. In PHYPA and SELML, 134 and 81 LRR-RLK genes have been retrieved, respectively. There is no significant difference in the average number of LRR-RLK genes between monocots and dicots but the number almost doubled in most angiosperms compared to PHYPA and SELML. It has to be noted that in some genomes (e.g. CARPA and LOTJA), the number of LRR-RLK genes is particularly low compared to other angiosperm genomes, suggesting that retention rates vary among genomes, and that many losses may have occured in some genomes (Fischer et al., 2016). Nevertheless, our results highlight the fact that, after the first wave of expansion in early land plants (Embryophyta), a second large amplification occured in angiosperm genomes which shaped the current LRR-RLK family size of more than 200 gene copies on average per genome.

Among the 20 SGs, SG_III, SG_XI, and SG_XIIa are the largest as they contain ~50% of the total number of LRR-RLK genes in the analyzed plant genomes (**Figure 1A**). The extensive expansions leading to their size do not follow the same amplification pattern (**Figure 1B**). This observation is highlighted in Figure 1B by the color code used for each species in the 20 SG-specific trees (with branches of monocots species in pink and red, and branches of dicots species in yellow, blue and green). Our results reveal that the high numbers of SG_XIIa genes is the consequence of many lineage-specific expansions (LSE) (See (Fischer et al., 2016) for details). These LSEs are relatively recent as they can be observed in phyla as well as species-specific lineages. On the contrary, in SG_III and SG_XI, expansions occured mainly before the early divergence of angiosperm lineages, even though LSEs can also be observed at different levels of resolution in the trees. Therefore, these numerous and diverse modes of expansions lead to complex paralogous and orthologous relationships.

## 3.2. 101 orthologous groups of monocot and dicot genes retained along angiosperm evolution
With the aim of transfering functional annotation from well studied genes from model
species to orthologous genes in other genomes, we first analysed in depth orthologous
relationships between monocots and dicots LRR-RLK genes in each SG. This analysis led us
to define what we named the "core set" of LRR-RLK genes in angiosperms: i.e. orthologous
genes which have not been completely lost in either monocots or dicots throughout the
angiosperms evolutionary history. To do so, the 20 SG-specific trees were scaned to locate
monocots/dicots bifurcations (Figure 2). Based on this analysis, 101 orthologous groups
(OGs) containing monocots and dicots sequences (named MD OGs) were characterized and
defined as the "core" set of LRR-RLK genes in angiosperms.

The SG analysis revealed that these 101 MD OGs are present in 19 of the 20 SGs,
with the majority of them in SG_III and SG_XI (Table 1). This highlights again the fact that
most of the LRR-RLK genes which underwent expansions before the monocots/dicots split
have been retained in these SGs. In order to go further in the description of orthologous
relationships, we qualified OGs as either “simple” or “complex”. In “simple” OGs, the
presence or absence of duplications within the monocot or dicot clades can clearly be inferred
from the phylogenetic tree. On the other hand, if several duplications occurred disorderly with
no obvious connection to the species tree we described these OGs as “complex” (OG_c).
Interestingly, these OG_c are over-represented in SG_IV, VIII-1, VIII-2, and XI. A total of
6739 genes are contained in the 101 OGs, representing 82.7% of the entire LRR-RLK gene
family. However, while 2956 genes are included in the 24 OG_c (average of 123.2 genes per
OG), 3783 genes belong to the 77 non complex OGs (average of 49.2 genes per OG),
highlighting differences in expansion/retention rates between these OGs. Moreover, looking
at the percentage of genes contained in OGs per SG reveals that, except for SG_VIII-2 and
XIIb, more than 70% of the LRR-RLK genes belong to OGs, with genes mainly in complex
OGs in SG_I, IV, VIII-1, VIII-2, XI and XIIa (Figure 3A). In some SGs, some OGs contain a
very large number of genes, such as 405 genes in one of SG_I OG (SG_I-3c), or 307 and 708
genes in two SG_XIIa OGs (Figure 3B and Suppl Table 2 for details). In these large OGs,
many species-specific duplications occured, and among the 10 OGs containing more than 100
genes, 8 are complex. In SG_I-3c for example, several genes have been studied in
Arabidopsis, such as IMPAIRED OOMYCETE SUSCEPTIBILITY 1 (IOS1), FLG22-
INDUCED RECEPTOR-LIKE KINASE 1 (FRK1), and light-repressible receptor protein
kinase (LRRPK) (Deeken and Kaldenhoff, 1997; Asai et al., 2002). All have been reported to
be involved in abiotic and biotic responses in dicots but no gene from the same OG have been
described so far in monocots. However, the fact that these genes are classified into the
complex mode of expansion suggests that in monocots too, these genes could be involved in
stress response.

**3.3. More than 10% of the LRR-RLK core sets are absent in Brassicaceae**

For each OG, we investigated whether some species were lacking members, focusing
particularly on the Brassicales, for which 6 species are included in our analysis (Figure 4).
Moreover, this clade contains the model plant Arabidopsis which is the reference for many
studies on LRR-RLK functions (Wu et al., 2016). Interestingly, among the 101 OGs, 14
(13.8%) have been completely lost in the Brassicaceae, and 3 of them are even absent in all
the Brassicales. This observation is an incentive to the extensive study of these receptors in
other plants than Arabidopsis, adding an argument to the fact that functions or interactions are
sometimes phylum-specific. For example, the SG_I-2 OG contains the SYMBIOSIS
RECEPTOR LIKE KINASE (SYMRK, also known as NORK or DMI2) receptor which is
involved in actinorhizal plants and legumes, respectively, in phosphate-acquiring arbuscular
mycorrhiza and nitrogen-fixing root nodule symbiosis (Antolin-Llovera et al., 2014a). For
this gene, we noticed that besides being absent in the Brassicaceae, which do not form mycorrhizal associations and root nodule symbiosis with rhizobia, some other characteristics of these receptors have been observed in monocots (see below).

## 3.4. Inference of functional information from experimentally characterized LRR-RLK genes to uncharacterized genes

The use of orthologous relationships to infer functional annotations relies on the fact that orthologs are expected to carry equivalent functions in different organisms. However, this can only be reliably inferred if, at least, structural characteristics and domain architecture are conserved between othologs. In our analysis, the phylogeny of the LRR-RLK proteins was computed on the well conserved KDs. However, LRR-RLK sequences are composed of several domains and, notably, of LRRs in their ECD. One could wonder whether the domains belonging to genes of the same OG are conserved. First, we took a detailed look at the predicted number of LRRs of all these receptors. The number of LRR motifs per protein is an important feature for homo- and hetero-complex formation between LRR-RLKs (Macho and Zipfel, 2014). Second, we investigated the presence of island domains in between LRRs. These domains have been described to be the binding site for the brassinosteroid hormone in some receptors (Kinoshita et al., 2005; Hothorn et al., 2011; She et al., 2011). Third, we analysed the presence of the MLD, a carbohydrate-binding domain, and the GDPC, a protein cleavage motif, which were shown to be located before the LRRs in some SG_I receptors. Fourth, we investigated the presence of Cys-pairs surrounding the LRRs in some SGs. The presence and organization of these domains is functionally important and has to be taken into account for transferring functional informations between orthologous genes. The description of structural features localized in the ECDs of these receptors allows a subclassification which, although reflecting the phylogeny of the KD of these receptors, also takes the structural differences of the ECD into account. Therefore, we subdivided the 20 SGs further according to these characteristics in their ECD (Figure 5 and Supplementary Table 3 for details).

### 3.4.1. Number of LRR motifs

First, for the moss and/or spikemoss genes which are orthologous to the angiosperm core sets of genes, we investigated if some sets of receptors varied in the number of LRR motifs. All those genes have a common ancestor, predating the divergence between moss and/or spikemoss and angiosperms and their KDs have all evolved in concert for ~450 MYA. Despite speciation events, the close phylogenetic relationship of all these LRR-RLK genes with moss and spikemoss orthologs suggests that signaling pathways downstream of these receptors could be conserved. In the OG containing the FLS2 receptor (SG_XIIa), we noticed that the number of LRRs in the PHYPA orthologs was lower than in angiosperms (Figure 6A). This peculiar differences noted in the PHYPA ECDs of the FLS2 orthologs could affect ligand binding or even suggest that ligands are not conserved. This would be in agreement with publications stating that the moss *Physcomitrella patens* does not carry an FLS2 ortholog and also shows no response to flg22 (Boller and Felix, 2009; Tanigaki et al., 2014). All other core gene sets for which *Physcomitrella/Selaginella* ECDs are conserved compared to angiosperms, provide interesting cases for which it would be worth to verify if the functions described for monocots and/or dicots members are entirely conserved in bryophytes and lycopsids.

Second, we focused on the number of predicted LRRs in the 7,767 LRR-containing sequences. Even if the number of LRRs per sequence is very variable, the distribution of the number of LRR per sequences shows three peaks at five, ten or 21 (Figure 6B). This
observation suggests that these numbers of LRRs per sequences may be optimal for the 3D conformation of these receptors and their interactions in homo- or heterocomplexes. To our knowledge, this observation has not been explicitly made in animal LRR-containing proteins but could also be true (Ng et al., 2011). In plants, one hetero-oligomeric protein complex has been described between the brassinosteroid (BR) receptor BRI1 (SG_Xb M4 in **Figure 5**) and BAK1/SERK3 or SERK1 receptors (SG_II B.3 in **Figure 5**) (Aker and de Vries, 2008;Chinchilla et al., 2009;Santiago et al., 2013;Sun et al., 2013). The complex crystal structure of SERK1 and BRI1 has revealed that the BRI1 C-terminal LRRs form a docking platform for the LRRs of the SERK1 co-receptor (Hothorn et al., 2011;Santiago et al., 2013). The SERK proteins have also been shown to serve various other BR-independent functions by forming heterocomplexes with SG_XIIa receptors (FLS2 and EFR, structure O in **Figure 5**) and the PEP1 RECEPTOR proteins (PEPR1, SG_XI N4 in **Figure 5**) (Chinchilla et al., 2007;Heese et al., 2007;Albrecht et al., 2008;Schulze et al., 2010;Roux et al., 2011;Koller and Bent, 2014). In rice, OsSERK2 (SG_II B.3 in **Figure 5**) forms a constitutive complex with the LRR-RLK Xa21 (SG_XIIa O in **Figure 5**) (Chen et al., 2014). Thus, these SERK co-receptors (4-5 LRRs) seem to play a central role in the regulation of multiple LRR-RLKs (>20 LRRs) by interacting directly with them (Aker and de Vries, 2008;Chinchilla et al., 2009;Kim et al., 2013;Santiago et al., 2013;Sun et al., 2013). Interestingly, SERK3/BAK1 has also been found in complex with the BAK1-INTERACTING RECEPTOR KINASE1 and 2 (BIR1 and BIR2) proteins, two receptors belonging to SG_Xa, another SG possessing five LRRs in its ECD (structure L in **Figure 5**) (Gao et al., 2009;Halter et al., 2014). The association between SERK3/BAK1 and BIR1/2 prevents the FLS2-BAK1 interaction before elicitation of the immune response. It is still unknown if the co-receptor status is restricted to the SERK genes subfamily of SG_II B.3, or if other LRR-RLK SGs could also be part of various signaling pathways in interaction with receptors possessing 20-25 LRRs. Indeed, other SGs such as SG_XI, SG_XIIIa, and SG_XIV corresponding to structures N2, Q and S1 respectively (**Figure 5**), also possess five LRRs in their ECD. In SG_XI-22 (structure N2 in **Figure 5**), the receptor SUPPRESSOR OF BIR1 1 (SOBIR1) also named EVERSHED (EVR) has been shown to be involved in floral organ shedding and in the regulation of several resistance signaling pathways with the BIR1, BAK1, and FLS2 receptors (Gao et al., 2009;Leslie et al., 2010). Moreover, the SOBIR1 receptor is also described as a coreceptor/adaptor forming complexes with many LRR-Receptor like proteins (LRR receptors devoid of a KD), suggesting that the SERK-type receptors (5-LRRs) could be considered as general adaptors important for functionality in complex with their receptor partners (Liebrand et al., 2013;Gust and Felix, 2014). In SG_XIIIa (structure Q), the FEIs receptors (FEI1 and FEI2, named after the Chinese word for fat), whose single mutants were indistinguishable from the wild type in development, work against a co-receptor function (Xu et al., 2008). However, for other SGs possessing five LRRs in their ECD, the question about the putative co-receptor function will remain unanswered until further molecular characterization is performed. The receptors belonging to SG_II B.3, contrary to SG_II B.1 and B.2, contain also a Pro-rich motif in their ECD. The question on the functionality of the Pro-rich domain in the ECD of the SERK proteins also remains to be answered. This domain of unknown function could provide a flexible hinge to the ECD (Schmidt et al., 1997;Kay et al., 2000;Baudino et al., 2001;Hecht et al., 2001;Chevalier et al., 2005). Interestingly, other SGs possess these kinds of motifs, e.g. SG_VI F3.2, for which no receptor has been studied yet.

439

### 3.4.2. Island domains

We positioned all the predicted LRRs on the 7,767 proteins and searched for islands between them (Suppl Table 4 LINK TO SUPPL TABLES). These islands are of particular

441

442
importance as they are the brassinolide hormone binding sites for the BRI1 and BRI1-like (BRL) receptors (structures M2-M4 in **Figure 5**) belonging to SG_Xb (Hothorn et al., 2011; She et al., 2011). We found that in SG_Xb, all sets of orthologs possess an island encompassing two (Structures M1 in **Figure 5**) or at least three (structures M2-M5) LRRs. One could ask if the island domain in genes of structure M5, for which no function has been described up to now and very similar to the M2-M4 structures, is also a binding site for the brassinosteroid hormone. The remaining OGs in SG_Xb (structure M1) contain the three Arabidopsis genes: PSKR1, PSKR2 and PLANT PEPTIDE CONTAINING SULFATED TYROSINE 1 RECEPTOR (PSY1R). These receptors have overlapping functions in promoting cellular proliferation, longevity and expansion (Matsubayashi et al., 2006; Amano et al., 2007; Hartmann et al., 2013). The PSKR subfamily is also required for PSK peptide signaling in sexual reproduction in plants (Stuhrwohldt et al., 2015). Moreover, these receptors play a role in modifying responses to biotic pathogens and wounding (Loivamaki et al., 2010; Mosher et al., 2013). In the three Arabidopsis PSKRs, one island of ~60 AA was detected in addition to other smaller ones specific to each receptor. These islands could be important for hormone binding like in BRI1 and BRL receptors (Hothorn et al., 2011; She et al., 2011). Indeed, it has been shown that the brassinosteroid hormone could play a role in the signaling pathways activated downstream of these receptors (Hartmann et al., 2013). In SG_IX (structure K) and SG_XV (structure T1), islands encompassing the size of at least two LRRs are also present. In SG_IX, the crystal structure of the Arabidopsis TRANSMEMBRANE KINASE 1 (TMK1, also known as BLK1 (BARK1-like Kinase 1)) suggests that the islands could be critical for structural integrity. In SG_XV, the crystal structure of RECEPTOR-LIKE PROTEIN KINASE 2 (RPK2, also known as TOAD2) suggests that the islands could be the site for ligand binding as in BRI1 (Liu et al., 2013; Song et al., 2014).

### 3.4.3. Additional domains

As mentioned previously, a MLD lying in between the SP and the LRRs has been described in some SG_I receptors (Figure 5) (Hok et al., 2011; Antolin-Llovera et al., 2014a). One of them is the SYMRK receptor (structure A1.1) involved in mycorrhizal associations and rhizobium-legumes symbiosis, but its exact function is still unclear (Antolin-Llovera et al., 2014a). Recently, it has been demonstrated that the SYMRK receptor is cleaved at a GDPC motif placed at the end of the MLD to release the N-glycosylated ectodomain in the absence of symbiotic stimulation (Antolin-Llovera et al., 2014b). Moreover, protein cleavage on this motif would permit a physical interaction with the LysM-type RLK NOD FACTOR RECEPTOR 5 (NFR5) and induces a rapid degradation of the SYMRK protein lacking its MLD. In this form, SYMRK could act as a co-receptor to initiate symbiotic signaling with NFR5 and would mirror the role played by the receptors of the SERK family. After MLD release, the structure of SYMRK could indeed resemble that of BAK1/SERK3. Our analysis of structural features of ECDs reveals that the monocot orthologs of the SYMRK receptor were much shorter than the dicot ones and that the monocot receptors are devoid of the MLD present in dicots (structure A1.2). Thus, the SYMRK activation mechanism will have to be further investigated in monocots to evaluate if it can fit into the dicots model, or if other receptors possessing a MLD (SG_I or others, see below) are involved in this process.

We therefore looked for malectin domains in all LRR-RLK sequences and found that SG_VIII-2 receptors also contain one. However, in SG_VIII-2, this domain is not located just after the SP but between the LRRs and the TM domain (Structure J in **Figure 5**). We also found the GDPC cleavage motif in most of the SG_IV, V and SG_VIII-1 receptors. However, contrary to SG_I, the GDPC site is located just before the first Cys-pairs in all other SGs. In
these SGs, the Cys of the GDPC motif is the first site of the Cys-pair. It is still unknown if
these receptors are also cleaved at this site and what the functional consequences would be. In
SG_VIII-2, which contains a malectin domain C-terminal of the LRRs, no GDPC sites are
present. This does not exclude the possibility that another cleavage site could be used to
truncate the protein. Thus, the function of the malectin domains in SG_VIII-2 will have to be
explored in the future to decipher their exact functional role and their potential involvement in
protein stabilization.

500

3.5. Positive selection in the divergence between ancestrally duplicated OGs

Twelve pairs of MD OGs present in almost all monocot and dicot species, and
harboring a gene topology fitting approximately the species tree, appear to be issued from
ancestral duplications predating the monocot/dicot divergence (Table 2). As these OG pairs
had a similar ECD structural organization and were kept in almost all species studied here, we
searched for potential positive selection footprints in the divergence leading to their
differentiation. Although these genes do not all have a known function, it is expected that they
underwent amino-acid changes leading to their sub- and/or neo-functionalization, and one can
wonder if and how positive selection could have driven these changes.

510

To answer this question, we tested whether some sites underwent positive selection on
the two branches starting from the ancestral duplication and ending at the monocot/dicot
divergence node of each OG. The detailed results of this analysis are presented in Table 2 and
Suppl-Table 5. Two pairs showed no signal on either of the two branches (SG_Xa-1/2,
SG_XIIa-2/3). For the pairs SG_II-3/4 and SG_XIIIa-1/2 a signal was detected for one branch
only but the signal on SG_XIIIa-1/2 may be a false positive or the sign of a lack of power,
since no sites appeared to be significant (see Materials and Methods for details). The eight
other pairs showed a signal of positive selection on each branch. Although the model
indicating positive selection performs significantly better than the null model, two pairs
(SG_III-3/4 and SG_XI-19/20) have no significant sites for one of the two tested branches.
This again indicates either a false positive or a lack of power. It is also possible that positive
selection acted on a large number of sites which results in none of them exceeding the
significance threshold. Finally, five pairs (SG_III-8/9, SG_XI-8/9, SG_XI-14/15, SG_XI-
17/18 and SG_XIIIb-1/2) have a strong signal with up to 26 sites validated after manual
curation. This result shows that in about half of the tested cases, several amino acid changes
fixed in the divergence between these genes are compatible with a signal of positive selection.

527

A total of 141 sites were manually validated as having experienced an episode of
positive selection during MD OG genes divergence. For the five pairs with a strong signal, the
repartition of these sites across the different domains of the LRR-RLK protein showed that
the LRRs and KDs are the most affected (Figure 7). More than half of the sites (78) fall in the
ECD, among which 68 are in the LRR domain; 51 sites fall in the intracellular domain (ICD),
among which 42 are in the KD (Table 2). Considering that these LRR and KD are the largest
domains, we normalized the number of positively selected sites by domain size. Kinase and
LRR appeared then to be affected equally (Chi-square test, p=0.26) by positive selection. This
result is very different from positive selection signatures observed in the recent paralogs, for
which LRR is the most strongly affected domain (Fischer et al., 2016). The number of sites
laying in the LRR domain allowed us to look for any specific distribution across the 24 amino
acids composing the motif. The repartition of the sites affected by positive selection within
the LRR is not homogeneous and the majority of them (67) fall in the 13 non canonical
positions (Suppl Table 5). However no notable pattern emerge from their distribution (data
not shown). Again, this result contrasts with what is observed in lineage-specific expanded
543
544
545

genes for which four positions are predominantly affected (Fischer et al., 2016). The
remaining domains are affected by a number of sites varying from 1 to 9.

546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565

This approach revealed a prevalence of sites targeted by positive selection in the ECD for the three couples of genes belonging to SG_XI as well as SG_XIIb. On the opposite, a tendency to target intracellular domains can be observed for SG_III. Two pairs of OGs for which positive selection footprints are detected in the ECD correspond to OGs whose duplication gave rise to the PXY and PXLs clades in SG_XI-17/18 (Fisher and Turner, 2007; Jung et al., 2015); and to the ERECTA (ER) and ERECTA-like (ERL) clades in SG_XIIb-1/2 (Torii et al., 1996; Sanchez-Rodriguez et al., 2009). Other pairs of OGs concern differentiation of the PEPR1 and 2 clade (SG_XI-9), of the STERILITY-REGULATING KINASE MEMBER 2 (SKM2) gene (SG_XI-14) or of the DspA/E-interacting protein of Malus x domestica Borkh 1 and 3 (DIPM1 and 3) clade (SG_III-9), from their respective sister clades, SG_XI-8, SG_XI-15 and SG_III-8 (Meng et al., 2006; Krol et al., 2010; Kang and Hardtke, 2016). In these clades, no gene has been described yet. The strong signal of positive selection detected for these five groups of genes indicates that the divergence between ancestral copies may have procured a selective advantageous: in the domain involved in ligands or partners binding for ECDs, or in the domain affecting downstream signaling pathways for ICDs. Indeed, during the early expansion of LRR-RLK that took place before angiosperm split, some duplicated LRR-RLK differentiated by fixation of a higher number of non-synonymous than synonymous mutations at some amino acid sites, indicating the emergence of probably new advantageous functions.

566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592

**3.6 Conclusions**

In this report, we provide a framework to aid in the classification and give new insights to new prospects for functional analysis of some plant LRR-RLKs. We have defined the "core set" of the large LRR-RLK gene family and classified these receptors based on their ECD features. These analyses reveal that even if the KDs of the LRR-RLKs are phylogenetically related, the ECDs may have been subjected to major (e.g. loss of LRRs revealed by the structural features characterization) or minor (e.g. point mutations revealed by the traces of positive selection analysis) modifications during the evolution of orthologs. These alterations could affect ligand recognition sites, dimerization with other receptors, and/or other processes involved in signal transduction. Indeed, the proper signal transduction via receptor kinases is not restricted to the binding of ligands to receptors located at the plasma membrane. Tightly regulated steps for proper folding of the proteins, trafficking from endomembranes to plasma membranes, and finally internalization and recycling of the receptors after ligand binding play essential roles in signal transduction (Shah et al., 2002; Robatzek et al., 2006; Salomon and Robatzek, 2006; Irani and Russinova, 2009; Beck et al., 2012; Di Rubbo et al., 2013; Offringa and Huang, 2013; Martins et al., 2015). Recently, an enthusiastic wave swept over the plant receptor kinases community concerning endoplasmic reticulum quality control since most of these steps take place in this cellular compartment (Saijo, 2010; Su et al., 2011; Huttner and Strasser, 2012; Tintor and Saijo, 2014). Newly synthetised membrane-resident proteins translocate first into the endoplasmic reticulum where they are subjected to folding and modifications like formation of disulfide bridges. It is also the place where nascent polypeptides are glycosylated - the most common post-traductional modification which is a crucial event during protein folding and quality control processes (Bieberich, 2014). The LRR-RLKs are part of the large family of plant proteins which are N-glycosylated and many N-glycosylation acceptor sequences are present in all ECDs. In some pattern recognition receptors and receptors involved in developmental processes, proteins with mutations at residues which will create misfolded proteins have been shown to be part of
593 endoplasmic reticulum protein complexes and directed to degradation (Hong et al., 2008;Hong et al., 2009;Li et al., 2009;Nekrasov et al., 2009;Lee et al., 2011;Su et al., 2011;Hong et al., 2012;Huttner and Strasser, 2012;Sun et al., 2012;Park et al., 2013;Robatzek and Wirthmueller, 2013;Huttner et al., 2014). The significance of all the structural feature 594 modifications which have been mentioned above are still mostly unknown but classic 595 biochemical and cell biological studies (e.g. domain swapping among orthologs and/or 596 targeted point mutations using CRISPR/Cas9) should help to explore their functions in details 597 and will provide many novel insights into the molecular characterization of LRR-RLKs. 598
600
601
602
603
604
605
606
607
608
609
610
611
612
613 **Author contributions:**

614 NC, CP, EG and AD designed the study; GD and AD performed the LRR-RLK extraction;

615 JFD and AD performed the phylogenetic clustering; NC and IF performed the selection

616 footprint analysis; IF, NC, JFD, MB and AD analyzed the data; JFD, AD, NC and IF wrote

617 the article.

618
619 **Funding:** This work was supported by the German Research Foundation (DFG) grant number

620 FI 1984/1-1 to IF; the Agropolis Resource Center for Crop Conservation, Adaptation and

621 Diversity (ARCAD); the Centre de coopération Internationale de Recherche en Agronomie

622 pour le Développement (CIRAD) PhD fellowship to MB; and the Agence Nationale de la

623 Recherche (ANR, France) ANR-08-GENM-021 to AD, CP and EG.

624
References

Ahmad, R., Parfitt, D.E., Fass, J., Ogundiwin, E., Dhingra, A., Gradziel, T.M., Lin, D., Joshi, N.A., Martinez-Garcia, P.J., and Crisosto, C.H. (2011). Whole genome sequencing of peach (Prunus persica L.) for SNP identification and selection. BMC Genomics 12, 569. doi: 10.1186/1471-2164-12-569.

Aker, J., and De Vries, S.C. (2008). Plasma membrane receptor complexes. Plant Physiol 147, 1560-1564. doi: 10.1104/pp.108.120501.

Al-Dous, E.K., George, B., Al-Mahmoud, M.E., Al-Jaber, M.Y., Wang, H., Salameh, Y.M., Al-Azwani, E.K., Chaluvadi, S., Pontaroli, A.C., Debarry, J., Arondel, V., Ohlrogge, J., Saie, I.J., Suliman-Elmeer, K.M., Bennetzen, J.L., Kruegger, R.R., and Malek, J.A. (2011). De novo genome sequencing and comparative genomics of date palm (Phoenix dactylifera). Nat Biotech 29, 521-527.

Albrecht, C., Russinova, E., Kemmerling, B., Kwaaitaal, M., and De Vries, S.C. (2008). Arabidopsis SOMATIC EMBRYOGENESIS RECEPTOR KINASE proteins serve brassinosteroid-dependent and -independent signaling pathways. Plant Physiol 148, 611-619. doi: 10.1104/pp.108.123216.

Altschul, S.F., Madden, T.L., Schäffer, A.A., Zhang, J., Zhang, Z., Miller, W., and Lipman, D. (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 25, 3389-3402.

Amano, Y., Tsubouchi, H., Shinohara, H., Ogawa, M., and Matsubayashi, Y. (2007). Tyrosine-sulfated glycopeptide involved in cellular proliferation and expansion in Arabidopsis. Proc Natl Acad Sci U S A 104, 18333-18338. doi: 10.1073/pnas.0706403104.

Anisimova, M., Gil, M., Dufayard, J.F., Dessimoz, C., and Gascuel, O. (2011). Survey of branch support methods demonstrates accuracy, power, and robustness of fast likelihood-based approximation schemes. Syst Biol 60, 685-699. doi: 10.1093/sysbio/syr041.

Antolin-Llovera, M., Petutsching, E.K., Ried, M.K., Lipka, V., Nurnberger, T., Robatzek, S., and Parniske, M. (2014a). Knowing your friends and foes--plant receptor-like kinases as initiators of symbiosis or defence. New Phytol 204, 791-802. doi: 10.1111/nph.13117.

Antolin-Llovera, M., Ried, M.K., and Parniske, M. (2014b). Cleavage of the SYMBIOSIS RECEPTOR-LIKE KINASE ectodomain promotes complex formation with Nod factor receptor 5. Curr Biol 24, 422-427. doi: 10.1016/j.cub.2013.12.053.

Arakaki, M., Christin, P.A., Nyffeler, R., Lendel, A., Eggli, U., Ogburn, R.M., Spriggs, E., Moore, M.J., and Edwards, E.J. (2011). Contemporaneous and recent radiations of the world's major succulent plant lineages. Proc Natl Acad Sci U S A 108, 8379-8384. doi: 10.1073/pnas.1100628108.

Argout, X., Salse, J., Aury, J.-M., Guiltinan, M.J., Droc, G., Gouzy, J., Allegre, M., Chaparro, C., Legavre, T., Maximova, S.N., Abrouk, M., Murat, F., Fouet, O., Poulain, J., Ruiz, M., Roguet, Y., Rodier-Goud, M.,
659 Barbosa-Neto, J.F., Sabot, F., Kudrna, D., Ammiraju, J.S.S., Schuster, S.C., Carlson, J.E., Sallet, E., Schiex, 660 T., Dievart, A., Kramer, M., Gelley, L., Shi, Z., Berard, A., Viot, C., Boccara, M., Risterucci, A.M., Guignon, 661 V., Sabau, X., Axtell, M.J., Ma, Z., Zhang, Y., Brown, S., Bourge, M., Golser, W., Song, X., Clement, D., 662 Rivallan, R., Tahi, M., Akaza, J.M., Pitollat, B., Gramacho, K., D'hont, A., Brunel, D., Infante, D., Kebe, I., 663 Costet, P., Wing, R., Mccombie, W.R., Guiderdoni, E., Quetier, F., Panaud, O., Wincker, P., Bocs, S., and 664 Lanaud, C. (2011). The genome of *Theobroma cacao*. *Nat Genet* 43, 101-108.

665 Asai, T., Tena, G., Plotnikova, J., Willmann, M.R., Chiu, W.L., Gomez-Gomez, L., Boller, T., Ausubel, F.M., 666 and Sheen, J. (2002). MAP kinase signalling cascade in Arabidopsis innate immunity. *Nature* 415, 977-983.

667 Banks, J.A., Nishiyama, T., Hasebe, M., Bowman, J.L., Gribskov, M., Depamphilis, C., Albert, V.A., Aono, N., 668 Aoyama, T., Ambrose, B.A., Ashton, N.W., Axtell, M.J., Barker, E., Barker, M.S., Bennetzen, J.L., Bonawitz, 669 N.D., Chapple, C., Cheng, C., Correa, L.G.G., Dacre, M., Debarry, J., Dreyer, I., Elias, M., Engstrom, E.M., 670 Estelle, M., Feng, L., Finet, C.D., Floyd, S.K., Frommer, W.B., Fujita, T., Gramzow, L., Gutensohn, M., 671 Harholt, J., Hattori, M., Heyl, A., Hirai, T., Hiwatashi, Y., Ishikawa, M., Iwata, M., Karol, K.G., Koehler, B., 672 Kolukisaoglu, U., Kubo, M., Kurata, T., Lalonde, S., Li, K., Li, Y., Litt, A., Lyons, E., Manning, G., 673 Maruyama, T., Michael, T.P., Mikami, K., Miyazaki, S., Morinaga, S.-I., Murata, T., Mueller- Roeber, B., 674 Nelson, D.R., Obara, M., Oguri, Y., Olmstead, R.G., Onodera, N., Petersen, B.L., Pils, B., Prigge, M., Rensing, 675 S.A., Riano-Pachon, D.M., Roberts, A.W., Sato, Y., Scheller, H.V., Schulz, B., Schulz, C., Shakirov, E.V., 676 Shibagaki, N., Shinohara, N., Shippen, D.E., Sorensen, I., Sotooka, R., Sugimoto, N., Sugita, M., Sumikawa, 677 N., Tanurdzic, M., Theissen, G., Ulvskov, P., Wakazuki, S., Weng, J.-K., Willats, W.W.G.T., Wipf, D., Wolf, 678 P.G., Yang, L., Zimmer, A.D., Zhu, Q., Mitros, T., Hellsten, U., Loqué, D., Otillar, R., Salamov, A., Schmutz, 679 J., Shapiro, H., Lindquist, E., et al. (2011). The Selaginella Genome Identifies Genetic Changes Associated 680 with the Evolution of Vascular Plants. *Science* 332, 960-963. doi: 10.1126/science.1203810.

681 Baudino, S., Hansen, S., Brettschneider, R., Hecht, V.F., Dresselhaus, T., Lorz, H., Dumas, C., and Rogowsky, 682 P.M. (2001). Molecular characterisation of two novel maize LRR receptor-like kinases, which belong to the 683 SERK gene family. *Planta* 213, 1-10.

684 Beck, M., Zhou, J., Faulkner, C., Maclean, D., and Robatzek, S. (2012). Spatio-temporal cellular dynamics of 685 the Arabidopsis flagellin receptor reveal activation status-dependent endosomal sorting. *Plant Cell* 24, 4205- 686 4219. doi: 10.1105/tpc.112.100263.

687 Belkhadir, Y., Yang, L., Hetzel, J., Dangl, J.L., and Chory, J. (2014). The growth-defense pivot: crisis 688 management in plants mediated by LRR-RK surface receptors. *Trends Biochem Sci* 39, 447-456. doi: 689 10.1016/j.tibs.2014.06.006.

690 Bennetzen, J.L., Schmutz, J., Wang, H., Percifield, R., Hawkins, J., Pontaroli, A.C., Estep, M., Feng, L., 691 Vaughn, J.N., Grimwood, J., Jenkins, J., Barry, K., Lindquist, E., Hellsten, U., Deshpande, S., Wang, X., Wu, 692 X., Mitros, T., Triplett, J., Yang, X., Ye, C.-Y., Mauro-Herrera, M., Wang, L., Li, P., Sharma, M., Sharma, R., 693 Ronald, P.C., Panaud, O., Kellogg, E.A., Brutnell, T.P., Doust, A.N., Tuskan, G.A., Rokhsar, D., and Devos, 694 K.M. (2012). Reference genome sequence of the model plant Setaria. *Nat Biotech* 30, 555-561.
695 Bieberich, E. (2014). Synthesis, Processing, and Function of N-glycans in N-glycoproteins. *Adv Neurobiol* 9, 47-
696 70. doi: 10.1007/978-1-4939-1154-7_3.

697 Boller, T., and Felix, G. (2009). A renaissance of elicitors: perception of microbe-associated molecular patterns
698 and danger signals by pattern-recognition receptors. *Annu Rev Plant Biol* 60, 379-406. doi:
699 10.1146/annurev.arplant.57.032905.105346.

700 Capella-Gutierrez, S., Silla-Martinez, J.M., and Gabaldon, T. (2009). trimAl: a tool for automated alignment
701 trimming in large-scale phylogenetic analyses. *Bioinformatics* 25, 1972-1973. doi:
702 10.1093/bioinformatics/btp348.

703 Chan, A.P., Crabtree, J., Zhao, Q., Lorenzi, H., Orvis, J., Puiu, D., Melake-Berhan, A., Jones, K.M., Redman, J.,
704 Chen, G., Cahoon, E.B., Gedil, M., Stanke, M., Haas, B.J., Wortman, J.R., Fraser-Liggett, C.M., Ravel, J., and
705 Rabinowicz, P.D. (2010). Draft genome sequence of the oilseed species *Ricinus communis*. *Nat Biotech* 28,
706 951-956.

707 Characterization, T.F.I.P.C.F.G.G. (2007). The grapevine genome sequence suggests ancestral hexaploidization
708 in major angiosperm phyla. *Nature* 449, 463-467.

709 Chen, X., Zuo, S., Schwessinger, B., Chern, M., Canlas, P.E., Ruan, D., Zhou, X., Wang, J., Daudi, A., Petzold,
710 C.J., Heazlewood, J.L., and Ronald, P.C. (2014). An XA21-associated kinase (OsSERK2) regulates immunity
711 mediated by the XA21 and XA3 immune receptors. *Mol Plant* 7, 874-892. doi: 10.1093/mp/ssu003.

712 Chevalier, D., Batoux, M., Fulton, L., Pfister, K., Yadav, R.K., Schellenberg, M., and Schneitz, K. (2005).
713 STRUBBELIG defines a receptor kinase-mediated signaling pathway regulating organ development in
714 Arabidopsis. *Proc Natl Acad Sci U S A* 102, 9074-9079.

715 Chinchilla, D., Shan, L., He, P., De Vries, S., and Kemmerling, B. (2009). One for all: the receptor-associated
716 kinase BAK1. *Trends Plant Sci* 14, 535-541. doi: 10.1016/j.tplants.2009.08.002.

717 Chinchilla, D., Zipfel, C., Robatzek, S., Kemmerling, B., Nurnberger, T., Jones, J.D., Felix, G., and Boller, T.
718 (2007). A flagellin-induced complex of the receptor FLS2 and BAK1 initiates plant defence. *Nature* 448, 497-
719 500. doi: 10.1038/nature05999.

720 Clark, S.E., Williams, R.W., and Meyerowitz, E.M. (1997). The CLAVATA1 gene encodes a putative receptor
721 kinase that controls shoot and floral meristem size in Arabidopsis. *Cell* 89, 575-585.

722 Claros, M.G., and Von Heijne, G. (1994). TopPred II: an improved software for membrane protein structure
723 predictions. *Comput Appl Biosci* 10, 685-686.

724 Consortium, T.P.G.S. (2011). Genome sequence and analysis of the tuber crop potato. *Nature* 475, 189-195.

725 Consortium, T.T.G. (2012). The tomato genome sequence provides insights into fleshy fruit evolution. *Nature*
726 485, 635-641.
727 Conte, M.G., Gaillard, S., Lanau, N., Rouard, M., and Perin, C. (2008). GreenPhylDB: a database for plant comparative genomics. *Nucleic Acids Res* 36, D991-998. doi: 10.1093/nar/gkm934.

729 D'hont, A., Denoeud, F., Aury, J.M., Baurens, F.C., Carreel, F., Garsmeur, O., Noel, B., Bocs, S., Droc, G., 730 Rouard, M., Da Silva, C., Jabbari, K., Cardi, C., Poulain, J., Souquet, M., Labadie, K., Jourda, C., Lengelle, J., 731 Rodier-Goud, M., Alberti, A., Bernard, M., Correa, M., Ayyampalayam, S., Mckain, M.R., Leebens-Mack, J., 732 Burgess, D., Freeling, M., Mbeguie, A.M.D., Chabannes, M., Wicker, T., Panaud, O., Barbosa, J., Hribova, E., 733 Heslop-Harrison, P., Habas, R., Rivallan, R., Francois, P., Poiron, C., Kilian, A., Burthia, D., Jenny, C., Bakry, 734 F., Brown, S., Guignon, V., Kema, G., Dita, M., Waalwijk, C., Joseph, S., Dievart, A., Jaillon, O., Leclercq, J., 735 Argout, X., Lyons, E., Almeida, A., Jeridi, M., Dolezel, J., Roux, N., Risterucci, A.M., Weissenbach, J., Ruiz, 736 M., Glaszmann, J.C., Quetier, F., Yahiaoui, N., and Wincker, P. (2012). The banana (*Musa acuminata*) 737 genome and the evolution of monocotyledonous plants. *Nature* 488, 213-217. doi: 10.1038/nature11241.

738 Dassanayake, M., Oh, D.-H., Haas, J.S., Hernandez, A., Hong, H., Ali, S., Yun, D.-J., Bressan, R.A., Zhu, J.-K., 739 Bohnert, H.J., and Cheeseman, J.M. (2011). The genome of the extremophile crucifer Thellungiella parvula. 740 *Nat Genet* 43, 913-918.

741 Deeken, R., and Kaldenhoff, R. (1997). Light-repressible receptor protein kinase: a novel photo-regulated gene 742 from *Arabidopsis thaliana*. *Planta* 202, 479-486.

743 Desper, R., and Gascuel, O. (2002). Fast and accurate phylogeny reconstruction algorithms based on the 744 minimum-evolution principle. *J Comput Biol* 9, 687-705. doi: 10.1089/106652702761034136.

745 Di Rubbo, S., Irani, N.G., Kim, S.Y., Xu, Z.Y., Gadeyne, A., Dejonghe, W., Vanhoutte, I., Persiau, G., 746 Eeckhout, D., Simon, S., Song, K., Kleine-Vehn, J., Friml, J., De Jaeger, G., Van Damme, D., Hwang, I., and 747 Russinova, E. (2013). The clathrin adaptor complex AP-2 mediates endocytosis of brassinosteroid insensitive1 748 in Arabidopsis. *Plant Cell* 25, 2986-2997. doi: 10.1105/tpc.113.114058.

749 Dievart, A., and Clark, S.E. (2003). Using mutant alleles to determine the structure and function of leucine-rich 750 repeat receptor-like kinases. *Curr Opin Plant Biol* 6, 507-516. doi: S136952660300089X [pii].

751 Dufayard, J.-F., Duret, L., Penel, S., Gouy, M., Rechenmann, F., and Perrière, G. (2005). Tree pattern matching 752 in phylogenetic trees: automatic search for orthologs or paralogs in homologous gene sequence databases. 753 *Bioinformatics* 21, 2596-2603. doi: 10.1093/bioinformatics/bti325.

754 Eddy, S.R. (2009). A new generation of homology search tools based on probabilistic inference. *Genome Inform* 755 23, 205-211. doi: 9781848165632_0019 [pii].

756 Felsenstein, J. (1989). PHYLIP - Phylogeny Inference Package. *Cladistics* 5, 164-165.

757 Fischer, I., Dainat, J., Ranwez, V., Glemin, S., Dufayard, J.F., and Chantret, N. (2014). Impact of recurrent gene 758 duplication on adaptation of plant genomes. *BMC Plant Biol* 14, 151. doi: 10.1186/1471-2229-14-151.
759 Fischer, I., Dievart, A., Droc, G., Dufayard, J.F., and Chantret, N. (2016). Evolutionary Dynamics of the Leucine-Rich Repeat Receptor-Like Kinase (LRR-RLK) Subfamily in Angiosperms. *Plant Physiol* 170, 1595-1610. doi: 10.1104/pp.15.01470.

762 Fisher, K., and Turner, S. (2007). PXY, a receptor-like kinase essential for maintaining polarity during plant vascular-tissue development. *Curr Biol* 17, 1061-1066. doi: 10.1016/j.cub.2007.05.049.

764 Forest, F., and Chase, M.W. (2009). "Eurosid I," in *The timetree of life*, ed. K.S. Hedges Sb. New York: Oxford University Press), 188-196.

766 Gao, M., Wang, X., Wang, D., Xu, F., Ding, X., Zhang, Z., Bi, D., Cheng, Y.T., Chen, S., Li, X., and Zhang, Y. (2009). Regulation of cell death and innate immunity by two receptor-like kinases in Arabidopsis. *Cell Host Microbe* 6, 34-44. doi: 10.1016/j.chom.2009.05.019.

769 Garcia-Mas, J., Benjak, A., Sanseverino, W., Bourgeois, M., Mir, G., Gonzã¡Lez, V.M., Hâ©Naff, E., Cã¢Mara, F., Cozzuto, L., Lowy, E., Alioto, T., Capella-Gutiã©Rrez, S., Blanca, J., Caã±Izares, J., Ziarsolo, P., Gonzalez-Ibeas, D., Rodrăguez-Moreno, L., Droege, M., Du, L., Alvarez-Tejado, M., Lorente-Galdos, B., Melâ©, M., Yang, L., Weng, Y., Navarro, A., Marques-Bonet, T., Aranda, M.A., Nuez, F., Picã³, B.N., Gabaldã³n, T., Roma, G., Guigã³, R., Casacuberta, J.M., Arãºs, P., and Puigdomã¨Nech, P. (2012). The genome of melon (Cucumis melo L.). *Proc of the Natl Acad Sci U S A* 109, 11872-11877. doi: 10.1073/pnas.1205415109.

776 Goff, S.A., Ricke, D., Lan, T.-H., Presting, G., Wang, R., Dunn, M., Glazebrook, J., Sessions, A., Oeller, P., Varma, H., Hadley, D., Hutchison, D., Martin, C., Katagiri, F., Lange, B.M., Moughamer, T., Xia, Y., Budworth, P., Zhong, J., Miguel, T., Paszkowski, U., Zhang, S., Colbert, M., Sun, W.-L., Chen, L., Cooper, B., Park, S., Wood, T.C., Mao, L., Quail, P., Wing, R., Dean, R., Yu, Y., Zharkikh, A., Shen, R., Sahasrabudhe, S., Thomas, A., Cannings, R., Gutin, A., Pruss, D., Reid, J., Tavtigian, S., Mitchell, J., Eldredge, G., Scholl, T., Miller, R.M., Bhatnagar, S., Adey, N., Rubano, T., Tusneem, N., Robinson, R., Feldhaus, J., Macalma, T., Oliphant, A., and Briggs, S. (2002). A Draft Sequence of the Rice Genome (*Oryza sativa* L. ssp. japonica). *Science* 296, 92-100. doi: 10.1126/science.1068275.

784 Goodstein, D., Shu, S., Howson, R., Neupane, R., Hayes, R., Fazo, J., Mitros, T., Dirks, W., Hellsten, U., and Putnam, N. (2012). Phytozome: a comparative platform for green plant genomics. *Nucleic Acids Res* 40, D1178 - D1186.

787 Guindon, S., Dufayard, J.F., Lefort, V., Anisimova, M., Hordijk, W., and Gascuel, O. (2010). New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. *Syst Biol* 59, 307-321. doi: 10.1093/sysbio/syq010.

790 Guindon, S., and Gascuel, O. (2003). A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. *Syst Biol* 52, 696-704. doi: 54QHX07WB5K5XCX4 [pii].
792
793

Gust, A.A., and Felix, G. (2014). Receptor like proteins associate with SOBIR1-type of adaptors to form bimolecular receptor kinases. *Curr Opin Plant Biol* 21, 104-111. doi: 10.1016/j.pbi.2014.07.007.

794
795
796
797

Halter, T., Imkampe, J., Mazzotta, S., Wierzba, M., Postel, S., Bucherl, C., Kiefer, C., Stahl, M., Chinchilla, D., Wang, X., Nurnberger, T., Zipfel, C., Clouse, S., Borst, J.W., Boeren, S., De Vries, S.C., Tax, F., and Kemmerling, B. (2014). The leucine-rich repeat receptor kinase BIR2 is a negative regulator of BAK1 in plant immunity. *Curr Biol* 24, 134-143. doi: 10.1016/j.cub.2013.11.047.

798
799
800

Hartmann, J., Stuhrwohldt, N., Dahlke, R.I., and Sauter, M. (2013). Phytosulfokine control of growth occurs in the epidermis, is likely to be non-cell autonomous and is dependent on brassinosteroids. *Plant J* 73, 579-590. doi: 10.1111/tpj.12056.

801
802
803

Hecht, V., Vielle-Calzada, J.P., Hartog, M.V., Schmidt, E.D., Boutilier, K., Grossniklaus, U., and De Vries, S.C. (2001). The Arabidopsis SOMATIC EMBRYOGENESIS RECEPTOR KINASE 1 gene is expressed in developing ovules and embryos and enhances embryogenic competence in culture. *Plant Physiol* 127, 803-816.

804
805
806

Heese, A., Hann, D.R., Gimenez-Ibanez, S., Jones, A.M.E., He, K., Li, J., Schroeder, J.I., Peck, S.C., and Rathjen, J.P. (2007). The receptor-like kinase SERK3/BAK1 is a central regulator of innate immunity in plants. *Proc Natl Acad Sci USA* 104, 12217-12222.

807
808
809

Hok, S., Danchin, E.G., Allasia, V., Panabieres, F., Attard, A., and Keller, H. (2011). An Arabidopsis (malectin-like) leucine-rich repeat receptor-like kinase contributes to downy mildew disease. *Plant Cell Environ* 34, 1944-1957. doi: 10.1111/j.1365-3040.2011.02390.x.

810
811
812

Hong, Z., Jin, H., Fitchette, A.C., Xia, Y., Monk, A.M., Faye, L., and Li, J. (2009). Mutations of an alpha1,6 mannosyltransferase inhibit endoplasmic reticulum-associated degradation of defective brassinosteroid receptors in Arabidopsis. *Plant Cell* 21, 3792-3802. doi: 10.1105/tpc.109.070284.

813
814
815

Hong, Z., Jin, H., Tzfira, T., and Li, J. (2008). Multiple mechanism-mediated retention of a defective brassinosteroid receptor in the endoplasmic reticulum of Arabidopsis. *Plant Cell* 20, 3418-3429. doi: 10.1105/tpc.108.061879.

816
817
818

Hong, Z., Kajiura, H., Su, W., Jin, H., Kimura, A., Fujiyama, K., and Li, J. (2012). Evolutionarily conserved glycan signal to degrade aberrant brassinosteroid receptors in Arabidopsis. *Proc Natl Acad Sci USA* 109, 11437-11442. doi: 10.1073/pnas.1119173109.

819
820

Horn, M.A., and Walker, J.C. (1994). Biochemical properties of the autophosphorylation of RLK5, a receptor-like protein kinase from *Arabidopsis thaliana*. *Biochim. Biophys. Acta* 1208, 65-74.

821
822

Hothorn, M., Belkhadir, Y., Dreux, M., Dabi, T., Noel, J.P., Wilson, I.A., and Chory, J. (2011). Structural basis of steroid hormone perception by the receptor kinase BRI1. *Nature* 474, 467-471. doi: 10.1038/nature10153.

823
824

Hu, T.T., Pattyn, P., Bakker, E.G., Cao, J., Cheng, J.-F., Clark, R.M., Fahlgren, N., Fawcett, J.A., Grimwood, J., Gundlach, H., Haberer, G., Hollister, J.D., Ossowski, S., Ottilar, R.P., Salamov, A.A., Schneeberger, K.,
825 Spannagl, M., Wang, X., Yang, L., Nasrallah, M.E., Bergelson, J., Carrington, J.C., Gaut, B.S., Schmutz, J.,
826 Mayer, K.F.X., Van De Peer, Y., Grigoriev, I.V., Nordborg, M., Weigel, D., and Guo, Y.-L. (2011). The
827 *Arabidopsis lyrata* genome sequence and the basis of rapid genome size change. *Nat Genet* 43, 476-481.

828 Huang, S., Li, R., Zhang, Z., Li, L., Gu, X., Fan, W., Lucas, W.J., Wang, X., Xie, B., Ni, P., Ren, Y., Zhu, H.,
829 Li, J., Lin, K., Jin, W., Fei, Z., Li, G., Staub, J., Kilian, A., Van Der Vossen, E.a.G., Wu, Y., Guo, J., He, J.,
830 Jia, Z., Ren, Y., Tian, G., Lu, Y., Ruan, J., Qian, W., Wang, M., Huang, Q., Li, B., Xuan, Z., Cao, J., Asan,
831 Wu, Z., Zhang, J., Cai, Q., Bai, Y., Zhao, B., Han, Y., Li, Y., Li, X., Wang, S., Shi, Q., Liu, S., Cho, W.K.,
832 Kim, J.-Y., Xu, Y., Heller-Uszynska, K., Miao, H., Cheng, Z., Zhang, S., Wu, J., Yang, Y., Kang, H., Li, M.,
833 Liang, H., Ren, X., Shi, Z., Wen, M., Jian, M., Yang, H., Zhang, G., Yang, Z., Chen, R., Liu, S., Li, J., Ma, L.,
834 Liu, H., Zhou, Y., Zhao, J., Fang, X., Li, G., Fang, L., Li, Y., Liu, D., Zheng, H., Zhang, Y., Qin, N., Li, Z.,
835 Yang, G., Yang, S., Bolund, L., Kristiansen, K., Zheng, H., Li, S., Zhang, X., Yang, H., Wang, J., Sun, R.,
836 Zhang, B., Jiang, S., Wang, J., Du, Y., and Li, S. (2009). The genome of the cucumber, *Cucumis sativus* L. *Nat
837 Genet* 41, 1275-1281.

838 Huttner, S., and Strasser, R. (2012). Endoplasmic reticulum-associated degradation of glycoproteins in plants.
839 *Front Plant Sci* 3, 67. doi: 10.3389/fpls.2012.00067.

840 Huttner, S., Veit, C., Vavra, U., Schoberer, J., Dicker, M., Maresch, D., Altmann, F., and Strasser, R. (2014). A
841 context-independent N-glycan signal targets the misfolded extracellular domain of Arabidopsis STRUBBELIG
842 to endoplasmic-reticulum-associated degradation. *Biochem J* 464, 401-411. doi: 10.1042/BJ20141057.

843 Initiative, T.a.G. (2000). Analysis of the genome sequence of the flowering plant *Arabidopsis thaliana*. *Nature*
844 408, 796-815.

845 Initiative, T.I.B. (2010). Genome sequencing and analysis of the model grass *Brachypodium distachyon*. *Nature*
846 463, 763-768.

847 Irani, N.G., and Russinova, E. (2009). Receptor endocytosis and signaling in plants. *Curr Opin Plant Biol* 12,
848 653-659. doi: 10.1016/j.pbi.2009.09.011.

849 Jaouannet, M., Rodriguez, P.A., Thorpe, P., Lenoir, C.J., Macleod, R., Escudero-Martinez, C., and Bos, J.I.
850 (2014). Plant immunity in plant-aphid interactions. *Front Plant Sci* 5, 663. doi: 10.3389/fpls.2014.00663.

851 Jung, C.G., Hwang, S.G., Park, Y.C., Park, H.M., Kim, D.S., Park, D.H., and Jang, C.S. (2015). Molecular
852 characterization of the cold- and heat-induced Arabidopsis PXL1 gene and its potential role in transduction
853 pathways under temperature fluctuations. *J Plant Physiol* 176, 138-146. doi: 10.1016/j.jplph.2015.01.001.

854 Jung, S., Cestaro, A., Troggio, M., Main, D., Zheng, P., Cho, I., Folta, K., Sosinski, B., Abbott, A., Celton, J.-
855 M., Arus, P., Shulaev, V., Verde, I., Morgante, M., Rokhsar, D., Velasco, R., and Sargent, D. (2012). Whole
856 genome comparisons of Fragaria, Prunus and Malus reveal different modes of evolution between Rosaceous
857 subfamilies. *BMC Genomics* 13, 129.
858

Kajava, A.V. (1998). Structural diversity of leucine-rich repeat proteins. *J Mol Biol* 277, 519-527.

859

860

Kang, Y.H., and Hardtke, C.S. (2016). Arabidopsis MAKR5 is a positive effector of BAM3-dependent CLE45 signaling. *EMBO Rep* 17, 1145-1154. doi: 10.15252/embr.201642450.

861

862

Katoh, K., Misawa, K., Kuma, K., and Miyata, T. (2002). MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. *Nucleic Acids Res* 30, 3059-3066.

863

864

Kay, B.K., Williamson, M.P., and Sudol, M. (2000). The importance of being proline: the interaction of proline-rich motifs in signaling proteins with their cognate domains. *FASEB J* 14, 231-241.

865

866

Kim, B.H., Kim, S.Y., and Nam, K.H. (2013). Assessing the diverse functions of BAK1 and its homologs in arabidopsis, beyond BR signaling and PTI responses. *Mol Cells* 35, 7-16. doi: 10.1007/s10059-013-2255-3.

867

868

Kinoshita, T., Cano-Delgado, A., Seto, H., Hiranuma, S., Fujioka, S., Yoshida, S., and Chory, J. (2005). Binding of brassinosteroids to the extracellular domain of plant receptor kinase BRI1. *Nature* 433, 167-171.

869

870

Kobe, B., and Deisenhofer, J. (1994). The leucine-rich repeat: a versatile binding motif. *Trends Biochem Sci* 19, 415-421.

871

872

873

Koch, M.A., Haubold, B., and Mitchell-Olds, T. (2000). Comparative evolutionary analysis of chalcone synthase and alcohol dehydrogenase loci in Arabidopsis, Arabis, and related genera (Brassicaceae). *Mol Biol Evol* 17, 1483-1498.

874

875

Koller, T., and Bent, A.F. (2014). FLS2-BAK1 extracellular domain interaction sites required for defense signaling activation. *PLoS One* 9, e111185. doi: 10.1371/journal.pone.0111185.

876

877

Krogh, A., Larsson, B., Von Heijne, G., and Sonnhammer, E. (2001). Predicting transmembrane protein topology with a hidden markov model: application to complete genomes. *J Mol Biol* 305, 567-580.

878

879

880

Krol, E., Mentzel, T., Chinchilla, D., Boller, T., Felix, G., Kemmerling, B., Postel, S., Arents, M., Jeworutzki, E., and Al-Rasheid, K. (2010). Perception of the arabidopsis danger signal peptide 1 involves the pattern recognition receptor AtPEPR1 and its close homologue AtPEPR2. *J Biol Chem* 285, 13471-13479.

881

882

883

Lee, H.Y., Bowen, C.H., Popescu, G.V., Kang, H.G., Kato, N., Ma, S., Dinesh-Kumar, S., Snyder, M., and Popescu, S.C. (2011). Arabidopsis RTNLB1 and RTNLB2 Reticulon-like proteins regulate intracellular trafficking and activity of the FLS2 immune receptor. *Plant Cell* 23, 3374-3391. doi: 10.1105/tpc.111.089656.

884

885

Lehti-Shiu, M.D., Zou, C., Hanada, K., and Shiu, S.H. (2009). Evolutionary history and stress regulation of plant receptor-like kinase/pelle genes. *Plant Physiol* 150, 12-26. doi: 10.1104/pp.108.134353.

886

887

Leslie, M.E., Lewis, M.W., Youn, J.Y., Daniels, M.J., and Liljegren, S.J. (2010). The EVERSHED receptor-like kinase modulates floral organ shedding in Arabidopsis. *Development* 137, 467-476. doi: 10.1242/dev.041335.
888 Letunic, I., Doerks, T., and Bork, P. (2009). SMART 6: recent updates and new developments. *Nucleic Acids Res* 37, D229-232. doi: 10.1093/nar/gkn808.

890 Li, J., and Chory, J. (1997). A putative leucine-rich repeat receptor kinase involved in brassinosteroid signal transduction. *Cell* 90, 929-938.

892 Li, J., Zhao-Hui, C., Batoux, M., Nekrasov, V., Roux, M., Chinchilla, D., Zipfel, C., and Jones, J.D. (2009). Specific ER quality control components required for biogenesis of the plant innate immune receptor EFR. *Proc Natl Acad Sci USA* 106, 15973-15978. doi: 10.1073/pnas.0905532106.

895 Liebrand, T.W., Van Den Berg, G.C., Zhang, Z., Smit, P., Cordewener, J.H., America, A.H., Sklenar, J., Jones, A.M., Tameling, W.I., Robatzek, S., Thomma, B.P., and Joosten, M.H. (2013). Receptor-like kinase SOBIR1/EVR interacts with receptor-like proteins in plant immunity against fungal infection. *Proc Natl Acad Sci USA* 110, 10010-10015. doi: 10.1073/pnas.1220015110.

899 Liu, P., Hu, Z., Zhou, B., Liu, S., and Chai, J. (2013). Crystal structure of an LRR protein with two solenoids. *Cell Res* 23, 303-305. doi: 10.1038/cr.2012.159.

901 Liu, P., Wei, W., Ouyang, S., Zhang, J.S., Chen, S.Y., and Zhang, W.K. (2009). Analysis of expressed receptor-like kinases (RLKs) in soybean. *Journal of genetics and genomics* 36, 611-619. doi: 10.1016/S1673-8527(08)60153-8.

904 Loivamaki, M., Stuhrwohldt, N., Deeken, R., Steffens, B., Roitsch, T., Hedrich, R., and Sauter, M. (2010). A role for PSK signaling in wounding and microbial interactions in Arabidopsis. *Physiol Plant* 139, 348-357. doi: 10.1111/j.1399-3054.2010.01371.x.

907 Macho, A.P., and Zipfel, C. (2014). Plant PRRs and the activation of innate immune signaling. *Mol Cell* 54, 263-272. doi: 10.1016/j.molcel.2014.03.028.

909 Magallcdn, S., and Castillo, A. (2009). Angiosperm diversification through time. *Am J Bot* 96, 349-365.

910 Martins, S., Dohmann, E.M., Cayrel, A., Johnson, A., Fischer, W., Pojer, F., Satiat-Jeunemaitre, B., Jaillais, Y., Chory, J., Geldner, N., and Vert, G. (2015). Internalization and vacuolar targeting of the brassinosteroid hormone receptor BRI1 are regulated by ubiquitination. *Nat Commun* 6, 6151. doi: 10.1038/ncomms7151.

913 Matsubayashi, Y., Ogawa, M., Kihara, H., Niwa, M., and Sakagami, Y. (2006). Disruption and overexpression of Arabidopsis phytosulfokine receptor gene affects cellular longevity and potential for growth. *Plant Physiol* 142, 45-53. doi: 10.1104/pp.106.081109.

916 Meng, X., Bonasera, J.M., Kim, J.F., Nissinen, R.M., and Beer, S.V. (2006). Apple Proteins that Interact with DspA/E, a Pathogenicity Effector of Erwinia amylovora, the Fire Blight Pathogen. *Molecular Plant-Microbe Interactions* 19, 53-61.
919 Ming, R., Hou, S., Feng, Y., Yu, Q., Dionne-Laporte, A., Saw, J.H., Senin, P., Wang, W., Ly, B.V., Lewis, 920 K.L.T., Salzberg, S.L., Feng, L., Jones, M.R., Skelton, R.L., Murray, J.E., Chen, C., Qian, W., Shen, J., Du, P., 921 Eustice, M., Tong, E., Tang, H., Lyons, E., Paull, R.E., Michael, T.P., Wall, K., Rice, D.W., Albert, H., Wang, 922 M.-L., Zhu, Y.J., Schatz, M., Nagarajan, N., Acob, R.A., Guan, P., Blas, A., Wai, C.M., Ackerman, C.M., Ren, 923 Y., Liu, C., Wang, J., Wang, J., Na, J.-K., Shakirov, E.V., Haas, B., Thimmapuram, J., Nelson, D., Wang, X., 924 Bowers, J.E., Gschwend, A.R., Delcher, A.L., Singh, R., Suzuki, J.Y., Tripathi, S., Neupane, K., Wei, H., 925 Irikura, B., Paidi, M., Jiang, N., Zhang, W., Presting, G., Windsor, A., Navajas-Perez, R., Torres, M.J., Feltus, 926 F.A., Porter, B., Li, Y., Burroughs, A.M., Luo, M.-C., Liu, L., Christopher, D.A., Mount, S.M., Moore, P.H., 927 Sugimura, T., Jiang, J., Schuler, M.A., Friedman, V., Mitchell-Olds, T., Shippen, D.E., Depamphilis, C.W., 928 Palmer, J.D., Freeling, M., Paterson, A.H., Gonsalves, D., Wang, L., and Alam, M. (2008). The draft genome 929 of the transgenic tropical fruit tree papaya (Carica papaya Linnaeus). *Nature* 452, 991-996.

930 Mosher, S., Seybold, H., Rodriguez, P., Stahl, M., Davies, K.A., Dayaratne, S., Morillo, S.A., Wierzba, M., 931 Favery, B., Keller, H., Tax, F.E., and Kemmerling, B. (2013). The tyrosine-sulfated peptide receptors PSKR1 932 and PSY1R modify the immunity of Arabidopsis to biotrophic and necrotrophic pathogens in an antagonistic 933 manner. *Plant J* 73, 469-482. doi: 10.1111/tpj.12050.

934 Mu, J.-H., Lee, H.-S., and Kao, T.-H. (1994). Characterization of a pollen-expressed receptor-like kinase gene of 935 *Petunia inflata* and the activity of its encoded kinase. *Plant Cell* 6, 709-721.

936 Muschietti, J., Eyal, Y., and Mccormick, S. (1998). Pollen tube localization implies a role in pollen-pistil 937 interactions for the tomato receptor-like protein kinases LePRK1 and LePRK2. *Plant Cell* 10, 319-330.

938 Nekrasov, V., Li, J., Batoux, M., Roux, M., Chu, Z.H., Lacombe, S., Rougon, A., Bittel, P., Kiss-Papp, M., 939 Chinchilla, D., Van Esse, H.P., Jorda, L., Schwessinger, B., Nicaise, V., Thomma, B.P., Molina, A., Jones, 940 J.D., and Zipfel, C. (2009). Control of the pattern-recognition receptor EFR by an ER protein complex in plant 941 immunity. *EMBO J* 28, 3428-3438. doi: 10.1038/emboj.2009.262.

942 Ng, A.C., Eisenberg, J.M., Heath, R.J., Huett, A., Robinson, C.M., Nau, G.J., and Xavier, R.J. (2011). Human 943 leucine-rich repeat proteins: a genome-wide bioinformatic categorization and functional analysis in innate 944 immunity. *Proc Natl Acad Sci U S A* 108 Suppl 1, 4631-4638. doi: 10.1073/pnas.1000093107.

945 Noguchi, T., Fujioka, S., Choe, S., Takatsuto, S., Yoshida, S., Yuan, H., Feldmann, K.A., and Tax, F.E. (1999). 946 Brassinosteroid-insensitive dwarf mutants of Arabidopsis accumulate brassinosteroids. *Plant Physiol* 121, 743- 947 752.

948 Offringa, R., and Huang, F. (2013). Phosphorylation-dependent trafficking of plasma membrane proteins in 949 animal and plant cells. *J Integr Plant Biol* 55, 789-808. doi: 10.1111/jipb.12096.

950 Oh, D.H., Dassanayake, M., Haas, J.S., Kropornika, A., Wright, C., D'urzo, M.P., Hong, H., Ali, S., Hernandez, 951 A., Lambert, G.M., Inan, G., Galbraith, D.W., Bressan, R.A., Yun, D.J., Zhu, J.K., Cheeseman, J.M., and 952 Bohnert, H.J. (2010). Genome structures and halophyte-specific gene expression of the extremophile
953

*Thellungiella parvula* in comparison with *Thellungiella salsuginea* (*Thellungiella halophila*) and *Arabidopsis*.

954

*Plant Physiol* 154, 1040-1052. doi: 10.1104/pp.110.163923.

955

Oh, S.H., and Potter, D. (2005). Molecular phylogenetic systematics and biogeography of tribe Neillieae

956

(Rosaceae) using DNA sequences of cpDNA, rDNA, and LEAFY. *Am J Bot* 92, 179-192.

957

Park, C.-J., Sharma, R., Lefebvre, B., Canlas, P.E., and Ronald, P.C. (2013). The endoplasmic reticulum-quality

958

control component SDF2 is essential for XA21-mediated immunity in rice. *Plant science* 210, 53-60.

959

Paterson, A.H., Bowers, J.E., Bruggmann, R., Dubchak, I., Grimwood, J., Gundlach, H., Haberer, G., Hellsten,

960

U., Mitros, T., Poliakov, A., Schmutz, J., Spannagl, M., Tang, H., Wang, X., Wicker, T., Bharti, A.K.,

961

Chapman, J., Feltus, F.A., Gowik, U., Grigoriev, I.V., Lyons, E., Maher, C.A., Martis, M., Narechania, A.,

962

Otillar, R.P., Penning, B.W., Salamov, A.A., Wang, Y., Zhang, L., Carpita, N.C., Freeling, M., Gingle, A.R.,

963

Hash, C.T., Keller, B., Klein, P., Kresovich, S., Mccann, M.C., Ming, R., Peterson, D.G., Mehboob Ur, R.,

964

Ware, D., Westhoff, P., Mayer, K.F.X., Messing, J., and Rokhsar, D.S. (2009). The *Sorghum bicolor* genome

965

and the diversification of grasses. *Nature* 457, 551-556.

966

Prochnik, S., Marri, P., Desany, B., Rabinowicz, P., Kodira, C., Mohiuddin, M., Rodriguez, F., Fauquet, C.,

967

Tohme, J., Harkins, T., Rokhsar, D., and Rounsley, S. (2012). The Cassava Genome: Current Progress, Future

968

Directions. *Tropical Plant Biology* 5, 88-94. doi: 10.1007/s12042-011-9088-z.

969

Project, I.R.G.S. (2005). The map-based sequence of the rice genome. *Nature* 436, 793-800. doi:

970

10.1038/nature03895.

971

Reineke, A.R., Bornberg-Bauer, E., and Gu, J. (2011). Evolutionary divergence and limits of conserved non-

972

coding sequence detection in plant genomes. *Nucleic Acids Res* 39, 6029-6043. doi: 10.1093/nar/gkr179.

973

Rensing, S.A., Lang, D., Zimmer, A.D., Terry, A., Salamov, A., Shapiro, H., Nishiyama, T., Perroud, P.-F.O.,

974

Lindquist, E.A., Kamisugi, Y., Tanahashi, T., Sakakibara, K., Fujita, T., Oishi, K., Shin-I, T., Kuroki, Y.,

975

Toyoda, A., Suzuki, Y., Hashimoto, S.-I., Yamaguchi, K., Sugano, S., Kohara, Y., Fujiyama, A., Anterola, A.,

976

Aoki, S., Ashton, N., Barbazuk, W.B., Barker, E., Bennetzen, J.L., Blankenship, R., Cho, S.H., Dutcher, S.K.,

977

Estelle, M., Fawcett, J.A., Gundlach, H., Hanada, K., Heyl, A., Hicks, K.A., Hughes, J., Lohr, M., Mayer, K.,

978

Melkozernov, A., Murata, T., Nelson, D.R., Pils, B., Prigge, M., Reiss, B., Renner, T., Rombauts, S., Rushton,

979

P.J., Sanderfoot, A., Schween, G., Shiu, S.-H., Stueber, K., Theodoulou, F.L., Tu, H., Van De Peer, Y.,

980

Verrier, P.J., Waters, E., Wood, A., Yang, L., Cove, D., Cuming, A.C., Hasebe, M., Lucas, S., Mishler, B.D.,

981

Reski, R., Grigoriev, I.V., Quatrano, R.S., and Boore, J.L. (2008). The Physcomitrella genome reveals

982

evolutionary insights into the conquest of land by plants. *Science* 319, 64-69. doi: 10.1126/science.1150646.

983

Robatzek, S., Chinchilla, D., and Boller, T. (2006). Ligand-induced endocytosis of the pattern recognition

984

receptor FLS2 in Arabidopsis. *Genes Dev* 20, 537-542. doi: 10.1101/gad.366506.

985

Robatzek, S., and Wirthmueller, L. (2013). Mapping FLS2 function to structure: LRRs, kinase and its working

986

bits. *Protoplasma* 250, 671-681. doi: 10.1007/s00709-012-0459-6.
987 Rouard, M., Guignon, V., Aluome, C., Laporte, M.A., Droc, G., Walde, C., Zmasek, C.M., Perin, C., and Conte, M.G. (2011). GreenPhylDB v2.0: comparative and functional genomics in plants. *Nucleic Acids Res* 39, D1095-1102. doi: 10.1093/nar/gkq811.

990 Roux, M., Schwessinger, B., Albrecht, C., Chinchilla, D., Jones, A., Holton, N., Malinovsky, F.G., Tor, M., De Vries, S., and Zipfel, C. (2011). The Arabidopsis leucine-rich repeat receptor-like kinases BAK1/SERK3 and BKK1/SERK4 are required for innate immunity to hemibiotrophic and biotrophic pathogens. *Plant Cell* 23, 2440-2455. doi: 10.1105/tpc.111.084301.

994 Saijo, Y. (2010). ER quality control of immune receptors and regulators in plants. *Cell Microbiol* 12, 716-724. doi: 10.1111/j.1462-5822.2010.01472.x.

996 Sakamoto, T., Deguchi, M., Brustolini, O., Santos, A., Silva, F., and Fontes, E.P. (2012). The tomato RLK superfamily: phylogeny and functional predictions about the role of the LRRII-RLK subfamily in antiviral defense. *BMC Plant Biology* 12, 229.

999 Salomon, S., and Robatzek, S. (2006). Induced Endocytosis of the Receptor Kinase FLS2. *Plant Signal Behav* 1, 293-295.

1001 Sanchez-Rodriguez, C., Estevez, J.M., Llorente, F., Hernandez-Blanco, C., Jorda, L., Pagan, I., Berrocal, M., Marco, Y., Somerville, S., and Molina, A. (2009). The ERECTA Receptor-Like Kinase Regulates Cell Wall-Mediated Resistance to Pathogens in Arabidopsis thaliana. *Mol Plant Microbe Interact* 22, 953-963. doi: 10.1094/MPMI-22-8-0953.

1005 Santiago, J., Henzler, C., and Hothorn, M. (2013). Molecular mechanism for plant steroid receptor activation by somatic embryogenesis co-receptor kinases. *Science* 341, 889-892. doi: 10.1126/science.1242468.

1007 Sato, S., Hirakawa, H., Isobe, S., Fukai, E., Watanabe, A., Kato, M., Kawashima, K., Minami, C., Muraki, A., Nakazaki, N., Takahashi, C., Nakayama, S., Kishida, Y., Kohara, M., Yamada, M., Tsuruoka, H., Sasamoto, S., Tabata, S., Aizu, T., Toyoda, A., Shin-I, T., Minakuchi, Y., Kohara, Y., Fujiyama, A., Tsuchimoto, S., Kajiyama, S.I., Makigano, E., Ohmido, N., Shibagaki, N., Cartagena, J.A., Wada, N., Kohinata, T., Atefeh, A., Yuasa, S., Matsunaga, S., and Fukui, K. (2011). Sequence Analysis of the Genome of an Oil-Bearing Tree, *Jatropha curcas* L. *DNA Research* 18, 65-76. doi: 10.1093/dnares/dsq030.

1013 Sato, S., Nakamura, Y., Kaneko, T., Asamizu, E., Kato, T., Nakao, M., Sasamoto, S., Watanabe, A., Ono, A., Kawashima, K., Fujishiro, T., Katoh, M., Kohara, M., Kishida, Y., Minami, C., Nakayama, S., Nakazaki, N., Shimizu, Y., Shinpo, S., Takahashi, C., Wada, T., Yamada, M., Ohmido, N., Hayashi, M., Fukui, K., Baba, T., Nakamichi, T., Mori, H., and Tabata, S. (2008). Genome Structure of the Legume, *Lotus japonicus*. *DNA Research* 15, 227-239. doi: 10.1093/dnares/dsn008.

1018 Schmidt, E.D., Guzzo, F., Toonen, M.A., and De Vries, S.C. (1997). A leucine-rich repeat containing receptor-like kinase marks somatic plant cells competent to form embryos. *Development* 124, 2049-2062.
1020 Schmutz, J., Cannon, S.B., Schlueter, J., Ma, J., Mitros, T., Nelson, W., Hyten, D.L., Song, Q., Thelen, J.J.,
1021 Cheng, J., Xu, D., Hellsten, U., May, G.D., Yu, Y., Sakurai, T., Umezawa, T., Bhattacharyya, M.K., Sandhu,
1022 D., Valliyodan, B., Lindquist, E., Peto, M., Grant, D., Shu, S., Goodstein, D., Barry, K., Futrell-Griggs, M.,
1023 Abernathy, B., Du, J., Tian, Z., Zhu, L., Gill, N., Joshi, T., Libault, M., Sethuraman, A., Zhang, X.-C.,
1024 Shinozaki, K., Nguyen, H.T., Wing, R.A., Cregan, P., Specht, J., Grimwood, J., Rokhsar, D., Stacey, G.,
1025 Shoemaker, R.C., and Jackson, S.A. (2010). Genome sequence of the palaeopolyploid soybean. *Nature* 463,
1026 178-183.

1027 Schnable, P.S., Ware, D., Fulton, R.S., Stein, J.C., Wei, F., Pasternak, S., Liang, C., Zhang, J., Fulton, L.,
1028 Graves, T.A., Minx, P., Reily, A.D., Courtney, L., Kruchowski, S.S., Tomlinson, C., Strong, C., Delehaunty,
1029 K., Fronick, C., Courtney, B., Rock, S.M., Belter, E., Du, F., Kim, K., Abbott, R.M., Cotton, M., Levy, A.,
1030 Marchetto, P., Ochoa, K., Jackson, S.M., Gillam, B., Chen, W., Yan, L., Higginbotham, J., Cardenas, M.,
1031 Waligorski, J., Applebaum, E., Phelps, L., Falcone, J., Kanchi, K., Thane, T., Scimone, A., Thane, N., Henke,
1032 J., Wang, T., Ruppert, J., Shah, N., Rotter, K., Hodges, J., Ingenthron, E., Cordes, M., Kohlberg, S., Sgro, J.,
1033 Delgado, B., Mead, K., Chinwalla, A., Leonard, S., Crouse, K., Collura, K., Kudrna, D., Currie, J., He, R.,
1034 Angelova, A., Rajasekar, S., Mueller, T., Lomeli, R., Scara, G., Ko, A., Delaney, K., Wissotski, M., Lopez, G.,
1035 Campos, D., Braidotti, M., Ashley, E., Golser, W., Kim, H., Lee, S., Lin, J., Dujmic, Z., Kim, W., Talag, J.,
1036 Zuccolo, A., Fan, C., Sebastian, A., Kramer, M., Spiegel, L., Nascimento, L., Zutavern, T., Miller, B.,
1037 Ambroise, C., Muller, S., Spooner, W., Narechania, A., Ren, L., Wei, S., Kumari, S., Faga, B., Levy, M.J.,
1038 Mcmahan, L., Van Buren, P., Vaughn, M.W., et al. (2009). The B73 Maize Genome: Complexity, Diversity,
1039 and Dynamics. *Science* 326, 1112-1115. doi: 10.1126/science.1178534.

1040 Schulze, B., Mentzel, T., Jehle, A.K., Mueller, K., Beeler, S., Boller, T., Felix, G., and Chinchilla, D. (2010).
1041 Rapid heteromerization and phosphorylation of ligand-activated plant transmembrane receptors and their
1042 associated kinase BAK1. *J Biol Chem* 285, 9444-9451.

1043 Sebastian, P., Schaefer, H., Telford, I.R., and Renner, S.S. (2010). Cucumber (*Cucumis sativus*) and melon (*C. melo*) have numerous wild relatives in Asia and Australia, and the sister species of melon is from Australia.
1044
1045 *Proc Natl Acad Sci U S A* 107, 14269-14273. doi: 10.1073/pnas.1005338107.

1046 Shah, K., Russinova, E., Gadella, T.W., Jr., Willemse, J., and De Vries, S.C. (2002). The Arabidopsis kinase-associated protein phosphatase controls internalization of the somatic embryogenesis receptor kinase 1. *Genes Dev* 16, 1707-1720. doi: 10.1101/gad.220402.

1049 She, J., Han, Z., Kim, T.W., Wang, J., Cheng, W., Chang, J., Shi, S., Wang, J., Yang, M., Wang, Z.Y., and Chai, J. (2011). Structural insight into brassinosteroid perception by BRI1. *Nature* 474, 472-476.

1051 Shi, T., Huang, H., Sanderson, M.J., and Tax, F.E. (2014). Evolutionary dynamics of leucine-rich repeat receptor-like kinases and related genes in plants: A phylogenomic approach. *J. Integr. Plant Biol.* 56, 1744-17909.

1054 Shiu, S.H., and Bleecker, A.B. (2001a). Plant receptor-like kinase gene family: diversity, function, and signaling.
1055 *Sci STKE* 2001, RE22.
1056 Shiu, S.H., and Bleecker, A.B. (2001b). Receptor-like kinases from Arabidopsis form a monophyletic gene family related to animal receptor kinases. *Proc Natl Acad Sci U S A* 98, 10763-10768. doi: 10.1073/pnas.181141598.

1059 Shiu, S.H., and Bleecker, A.B. (2003). Expansion of the receptor-like kinase/Pelle gene family and receptor-like proteins in Arabidopsis. *Plant Physiol* 132, 530-543.

1061 Shiu, S.H., Karlowski, W.M., Pan, R., Tzeng, Y.H., Mayer, K.F., and Li, W.H. (2004). Comparative analysis of the receptor-like kinase family in Arabidopsis and rice. *Plant Cell* 16, 1220-1234.

1063 Song, W., Han, Z., Sun, Y., and Chai, J. (2014). Crystal structure of a plant leucine rich repeat protein with two island domains. *Sci China Life Sci* 57, 137-144. doi: 10.1007/s11427-013-4586-x.

1065 Song, W.Y., Wang, G.L., Chen, L.L., Kim, H.S., Pi, L.Y., Holsten, T., Gardner, J., Wang, B., Zhai, W.X., Zhu, L.H., Fauquet, C., and Ronald, P. (1995). A receptor kinase-like protein encoded by the rice disease resistance gene, Xa21. *Science* 270, 1804-1806.

1068 Song, X., Guo, P., Li, C., and Liu, C.-M. (2010). The cysteine pairs in CLV2 are not necessary for sensing the CLV3 peptide in shoot and root meristems. *J Integr Plant Biol* 52, 774-781.

1070 Stuhrwohldt, N., Dahlke, R.I., Kutschmar, A., Peng, X., Sun, M.X., and Sauter, M. (2015). Phytosulfokine peptide signaling controls pollen tube growth and funicular pollen tube guidance in Arabidopsis thaliana. *Physiol Plant* 153, 643-653. doi: 10.1111/ppl.12270.

1073 Su, W., Liu, Y., Xia, Y., Hong, Z., and Li, J. (2011). Conserved endoplasmic reticulum-associated degradation system to eliminate mutated receptor-like kinases in Arabidopsis. *Proc Natl Acad Sci U S A* 108, 870-875. doi: 10.1073/pnas.1013251108.

1076 Su, W., Liu, Y., Xia, Y., Hong, Z., and Li, J. (2012). The Arabidopsis homolog of the mammalian OS-9 protein plays a key role in the endoplasmic reticulum-associated degradation of misfolded receptor-like kinases. *Mol Plant* 5, 929-940. doi: 10.1093/mp/sss042.

1079 Sun, W., Cao, Y., Jansen Labby, K., Bittel, P., Boller, T., and Bent, A.F. (2012). Probing the Arabidopsis flagellin receptor: FLS2-FLS2 association and the contributions of specific domains to signaling function. *Plant Cell* 24, 1096-1113. doi: 10.1105/tpc.112.095919.

1082 Sun, Y., Han, Z., Tang, J., Hu, Z., Chai, C., Zhou, B., and Chai, J. (2013). Structure reveals that BAK1 as a co-receptor recognizes the BRI1-bound brassinolide. *Cell Research* 23, 1326-1329. doi: doi:10.1038/cr.2013.131.

1084 Tanigaki, Y., Ito, K., Obuchi, Y., Kosaka, A., Yamato, K.T., Okanami, M., Lehtonen, M.T., Valkonen, J.P., and Akita, M. (2014). Physcomitrella patens has kinase-LRR R gene homologs and interacting proteins. *PLoS One* 9, e95118. doi: 10.1371/journal.pone.0095118.
1087 Tintor, N., and Saijo, Y. (2014). ER-mediated control for abundance, quality, and signaling of transmembrane immune receptors in plants. *Front Plant Sci* 5, 65. doi: 10.3389/fpls.2014.00065.

1089 Torii, K.U., Mitsukawa, N., Oosumi, T., Matsuura, Y., Yokoyama, R., Whittier, R.F., and Komeda, Y. (1996). The Arabidopsis ERECTA gene encodes a putative receptor protein kinase with extracellular leucine-rich repeats. *Plant Cell* 8, 735-746.

1092 Tuskan, G., Difazio, S., Jansson, S.,Bohlmann, J., Grigoriev, I., Hellsten, U., Putnam, N., Ralph, S., Rombauts, S., and Salamov, A. (2006). The genome of black cottonwood, *Populus trichocarpa* (torr. & Gray). *Science* 313, 1596-1604.

1095 Varshney, R.K., Chen, W., Li, Y., Bharti, A.K., Saxena, R.K., Schlueter, J.A., Donoghue, M.T.A., Azam, S., Fan, G., Whaley, A.M., Farmer, A.D., Sheridan, J., Iwata, A., Tuteja, R., Penmetsa, R.V., Wu, W., Upadhyaya, H.D., Yang, S.-P., Shah, T., Saxena, K.B., Michael, T., Mccombie, W.R., Yang, B., Zhang, G., Yang, H., Wang, J., Spillane, C., Cook, D.R., May, G.D., Xu, X., and Jackson, S.A. (2012). Draft genome sequence of pigeonpea (*Cajanus cajan*), an orphan legume crop of resource-poor farmers. *Nat Biotech* 30, 83-89.

1100 Velasco, R., Zharkikh, A., Affourtit, J., Dhingra, A., Cestaro, A., Kalyanaraman, A., Fontana, P., Bhatnagar, S.K., Troggio, M., Pruss, D., Salvi, S., Pindo, M., Baldi, P., Castelletti, S., Cavaiuolo, M., Coppola, G., Costa, F., Cova, V., Dal Ri, A., Goremykin, V., Komjanc, M., Longhi, S., Magnago, P., Malacarne, G., Malnoy, M., Micheletti, D., Moretto, M., Perazzolli, M., Si-Ammour, A., Vezzulli, S., Zini, E., Eldredge, G., Fitzgerald, L.M., Gutin, N., Lanchbury, J., Macalma, T., Mitchell, J.T., Reid, J., Wardell, B., Kodira, C., Chen, Z., Desany, B., Niazi, F., Palmer, M., Koepke, T., Jiwan, D., Schaeffer, S., Krishnan, V., Wu, C., Chu, V.T., King, S.T., Vick, J., Tao, Q., Mraz, A., Stormo, A., Stormo, K., Bogden, R., Ederle, D., Stella, A., Vecchietti, A., Kater, M.M., Masiero, S., Lasserre, P., Lespinasse, Y., Allan, A.C., Bus, V., Chagne, D., Crowhurst, R.N., Gleave, A.P., Lavezzo, E., Fawcett, J.A., Proost, S., Rouze, P., Sterck, L., Toppo, S., Lazzari, B., Hellens, R.P., Durel, C.-E., Gutin, A., Bumgarner, R.E., Gardiner, S.E., Skolnick, M., Egholm, M., Van De Peer, Y., Salamini, F., and Viola, R. (2010). The genome of the domesticated apple (*Malus x domestica* Borkh.). *Nat Genet* 42, 833-839.

1112 Wang, H., Moore, M.J., Soltis, P.S., Bell, C.D., Brockington, S.F., Alexandre, R., Davis, C.C., Latvis, M., Manchester, S.R., and Soltis, D.E. (2009). Rosid radiation and the rapid rise of angiosperm-dominated forests. *Proc Natl Acad Sci U S A* 106, 3853-3858. doi: 10.1073/pnas.0813376106.

1115 Wang, K., Wang, Z., Li, F., Ye, W., Wang, J., Song, G., Yue, Z., Cong, L., Shang, H., Zhu, S., Zou, C., Li, Q., Yuan, Y., Lu, C., Wei, H., Gou, C., Zheng, Z., Yin, Y., Zhang, X., Liu, K., Wang, B., Song, C., Shi, N., Kohel, R.J., Percy, R.G., Yu, J.Z., Zhu, Y.X., Wang, J., and Yu, S. (2012). The draft genome of a diploid cotton *Gossypium raimondii*. *Nat Genet* 44, 1098-1103. doi: 10.1038/ng.2371.

1119 Wang, X., Wang, H., Wang, J., Sun, R., Wu, J., Liu, S., Bai, Y., Mun, J.-H., Bancroft, I., Cheng, F., Huang, S., Li, X., Hua, W., Wang, J., Wang, X., Freeling, M., Pires, J.C., Paterson, A.H., Chalhoub, B., Wang, B., Hayward, A., Sharpe, A.G., Park, B.-S., Weisshaar, B., Liu, B., Li, B., Liu, B., Tong, C., Song, C., Duran, C., Peng, C., Geng, C., Koh, C., Lin, C., Edwards, D., Mu, D., Shen, D., Soumpourou, E., Li, F., Fraser, F.,
1123 Conant, G., Lassalle, G., King, G.J., Bonnema, G., Tang, H., Wang, H., Belcram, H., Zhou, H., Hirakawa, H., Abe, H., Guo, H., Wang, H., Jin, H., Parkin, I.a.P., Batley, J., Kim, J.-S., Just, J., Li, J., Xu, J., Deng, J., Kim, J.A., Li, J., Yu, J., Meng, J., Wang, J., Min, J., Poulain, J., Hatakeyama, K., Wu, K., Wang, L., Fang, L., Trick, M., Links, M.G., Zhao, M., Jin, M., Ramchiary, N., Drou, N., Berkman, P.J., Cai, Q., Huang, Q., Li, R., Tabata, S., Cheng, S., Zhang, S., Zhang, S., Huang, S., Sato, S., Sun, S., Kwon, S.-J., Choi, S.-R., Lee, T.-H., Fan, W., Zhao, X., Tan, X., Xu, X., Wang, Y., Qiu, Y., Yin, Y., Li, Y., Du, Y., et al. (2011). The genome of the mesopolyploid crop species *Brassica rapa*. *Nat Genet* 43, 1035-1039.

1130 Wu, Y., Xun, Q., Guo, Y., Zhang, J., Cheng, K., Shi, T., He, K., Hou, S., Gou, X., and Li, J. (2016). Genome-Wide Expression Pattern Analyses of the Arabidopsis Leucine-Rich Repeat Receptor-Like Kinases. *Mol Plant* 9, 289-300. doi: 10.1016/j.molp.2015.12.011.

1133 Xu, S.L., Rahman, A., Baskin, T.I., and Kieber, J.J. (2008). Two leucine-rich repeat receptor kinases mediate signaling, linking cell wall biosynthesis and ACC synthase in Arabidopsis. *Plant Cell* 20, 3065-3079. doi: 10.1105/tpc.108.063354.

1136 Yang, Z. (2007). PAML 4: phylogenetic analysis by maximum likelihood. *Mol Biol Evol* 24, 1586-1591. doi: 10.1093/molbev/msm088.

1138 Young, N.D., Debelle, F., Oldroyd, G.E.D., Geurts, R., Cannon, S.B., Udvardi, M.K., Benedito, V.A., Mayer, K.F.X., Gouzy, J., Schoof, H., Van De Peer, Y., Proost, S., Cook, D.R., Meyers, B.C., Spannagl, M., Cheung, F., De Mita, S., Krishnakumar, V., Gundlach, H., Zhou, S., Mudge, J., Bharti, A.K., Murray, J.D., Naoumkina, M.A., Rosen, B., Silverstein, K.a.T., Tang, H., Rombauts, S., Zhao, P.X., Zhou, P., Barbe, V., Bardou, P., Bechner, M., Bellec, A., Berger, A., Berges, H., Bidwell, S., Bisseling, T., Choisne, N., Couloux, A., Denny, R., Deshpande, S., Dai, X., Doyle, J.J., Dudez, A.-M., Farmer, A.D., Fouteau, S., Franken, C., Gibelin, C., Gish, J., Goldstein, S., Gonzalez, A.J., Green, P.J., Hallab, A., Hartog, M., Hua, A., Humphray, S.J., Jeong, D.-H., Jing, Y., Jocker, A., Kenton, S.M., Kim, D.-J., Klee, K., Lai, H., Lang, C., Lin, S., Macmil, S.L., Magdelenat, G., Matthews, L., Mccorrison, J., Monaghan, E.L., Mun, J.-H., Najar, F.Z., Nicholson, C., Noirot, C., O/Bleness, M., Paule, C.R., Poulain, J., Prion, F., Qin, B., Qu, C., Retzel, E.F., Riddle, C., Sallet, E., Samain, S., Samson, N., Sanders, I., Saurat, O., Scarpelli, C., Schiex, T., Segurens, B., Severin, A.J., Sherrier, D.J., Shi, R., Sims, S., Singer, S.R., Sinharoy, S., Sterck, L., Viollet, A., Wang, B.-B., et al. (2011). The Medicago genome provides insight into the evolution of rhizobial symbioses. *Nature* 480, 520-524.

1151 Yu, J., Hu, S., Wang, J., Wong, G.K., Li, S., Liu, B., Deng, Y., Dai, L., Zhou, Y., Zhang, X., Cao, M., Liu, J., Sun, J., Tang, J., Chen, Y., Huang, X., Lin, W., Ye, C., Tong, W., Cong, L., Geng, J., Han, Y., Li, L., Li, W., Hu, G., Li, J., Liu, Z., Qi, Q., Li, T., Wang, X., Lu, H., Wu, T., Zhu, M., Ni, P., Han, H., Dong, W., Ren, X., Feng, X., Cui, P., Li, X., Wang, H., Xu, X., Zhai, W., Xu, Z., Zhang, J., He, S., Xu, J., Zhang, K., Zheng, X., Dong, J., Zeng, W., Tao, L., Ye, J., Tan, J., Chen, X., He, J., Liu, D., Tian, W., Tian, C., Xia, H., Bao, Q., Li, G., Gao, H., Cao, T., Zhao, W., Li, P., Chen, W., Zhang, Y., Hu, J., Liu, S., Yang, J., Zhang, G., Xiong, Y., Li, Z., Mao, L., Zhou, C., Zhu, Z., Chen, R., Hao, B., Zheng, W., Chen, S., Guo, W., Tao, M., Zhu, L., Yuan, L., and Yang, H. (2002). A draft sequence of the rice genome (*Oryza sativa* L. ssp. indica). *Science* 296, 79-92.
1159 Zan, Y., Ji, Y., Zhang, Y., Yang, S., Song, Y., and Wang, J. (2013). Genome-wide identification, characterization and expression analysis of populus leucine-rich repeat receptor-like protein kinase genes. BMC Genomics 14, 318.

1162 Zhang, G., Liu, X., Quan, Z., Cheng, S., Xu, X., Pan, S., Xie, M., Zeng, P., Yue, Z., Wang, W., Tao, Y., Bian, C., Han, C., Xia, Q., Peng, X., Cao, R., Yang, X., Zhan, D., Hu, J., Zhang, Y., Li, H., Li, H., Li, N., Wang, J., Wang, C., Wang, R., Guo, T., Cai, Y., Liu, C., Xiang, H., Shi, Q., Huang, P., Chen, Q., Li, Y., Wang, J., Zhao, Z., and Wang, J. (2012). Genome sequence of foxtail millet (Setaria italica) provides insights into grass evolution and biofuel potential. Nat Biotech 30, 549-554.

1167

1168

1169
Tables

<table><caption>**Table 1.** Number of MD OGs per SG</caption><thead><tr><th>SG</th><th>Total number of MD OG</th><th>Number of OG_c</th></tr></thead><tbody><tr><td>SG_I</td><td>3</td><td>1</td></tr><tr><td>SG_II</td><td>5</td><td>1</td></tr><tr><td>SG_III</td><td>24</td><td>3</td></tr><tr><td>SG_IV</td><td>1</td><td>1</td></tr><tr><td>SG_V</td><td>4</td><td></td></tr><tr><td>SG_VI</td><td>5</td><td></td></tr><tr><td>SG_VIIa</td><td>3</td><td></td></tr><tr><td>SG_VIIb</td><td>2</td><td></td></tr><tr><td>SG_VIII-1</td><td>3</td><td>3</td></tr><tr><td>SG_VIII-2</td><td>3</td><td>3</td></tr><tr><td>SG_IX</td><td>3</td><td>1</td></tr><tr><td>SG_Xa</td><td>2</td><td></td></tr><tr><td>SG_Xb</td><td>7</td><td></td></tr><tr><td>SG_XI</td><td>22</td><td>9</td></tr><tr><td>SG_XIIa</td><td>6</td><td>2</td></tr><tr><td>SG_XIIb</td><td>1</td><td></td></tr><tr><td>SG_XIIla</td><td>2</td><td></td></tr><tr><td>SG_XIIlb</td><td>2</td><td></td></tr><tr><td>SG_XIV</td><td>0</td><td></td></tr><tr><td>SG_XV</td><td>3</td><td></td></tr><tr><td>Total</td><td>101</td><td>24</td></tr></tbody></table>
<table><caption>**Table 2.** Positive selection on branches of pairs of OGs</caption><thead><tr><th>SG</th><th>number of sequences</th><th>name of branch</th><th>structural SG</th><th>known genes</th><th>model</th><th>number of parameters</th><th>InL</th><th>P-value</th><th>number of sites</th><th>ECD</th><th>ICD</th></tr></thead><tbody><tr><td rowspan="4">SG_II</td><td rowspan="4">101</td><td>SG_II-3</td><td>B.1</td><td>DOCS1</td><td>A<sub>0</sub></td><td>203</td><td>-47756.49038</td><td>0.2960935</td><td>ns</td><td>-</td><td></td></tr><tr><td rowspan="3">SG_II-4</td><td rowspan="3">B.2</td><td rowspan="3">NIK3</td><td>A</td><td>204</td><td>-47755.94453</td><td rowspan="3">0.0012474</td><td rowspan="3">*</td><td rowspan="3">3</td><td rowspan="3">3</td></tr><tr><td>A<sub>0</sub></td><td>203</td><td>-47761.1425</td></tr><tr><td>A</td><td>204</td><td>-47755.93314</td></tr><tr><td rowspan="4">SG_III</td><td rowspan="4">92</td><td>SG_III-3</td><td>C2</td><td></td><td>A<sub>0</sub></td><td>185</td><td>-87752.00789</td><td>0.0013461</td><td>*</td><td>1</td><td>1</td></tr><tr><td rowspan="3">SG_III-4</td><td rowspan="3">C2</td><td></td><td>A</td><td>186</td><td>-87746.86878</td><td rowspan="3">0.0000735</td><td rowspan="3">**</td><td rowspan="3">0</td><td></td></tr><tr><td>A<sub>0</sub></td><td>185</td><td>-87749.51178</td></tr><tr><td>A</td><td>186</td><td>-87741.65272</td></tr><tr><td rowspan="4">SG_III</td><td rowspan="4">126</td><td>SG_III-9</td><td>C6.1</td><td>DIPM1, DIPM3</td><td>A<sub>0</sub></td><td>253</td><td>-58707.08864</td><td>0.000000</td><td>**</td><td>11</td><td>4</td><td>7</td></tr><tr><td rowspan="3">SG_III-8</td><td rowspan="3">C6.2</td><td></td><td>A</td><td>254</td><td>-58683.61908</td><td rowspan="3">0.000000</td><td rowspan="3">**</td><td rowspan="3">5</td><td rowspan="3">1</td><td rowspan="3">4</td></tr><tr><td>A<sub>0</sub></td><td>253</td><td>-58712.53117</td></tr><tr><td>A</td><td>254</td><td>-58695.21242</td></tr><tr><td rowspan="4">SG_Xa</td><td rowspan="4">145</td><td>SG_Xa-1</td><td>L</td><td>BIR1</td><td>A<sub>0</sub></td><td>291</td><td>-80168.41936</td><td>0.013983288</td><td>ns</td><td>-</td><td></td></tr><tr><td rowspan="3">SG_Xa-2</td><td rowspan="3">L</td><td rowspan="3">BIR2, BIR3, BIR4</td><td>A</td><td>292</td><td>-80165.39924</td><td rowspan="3">0.031737545</td><td rowspan="3">ns</td><td rowspan="3">-</td><td></td></tr><tr><td>A<sub>0</sub></td><td>291</td><td>-80166.7335</td></tr><tr><td>A</td><td>292</td><td>-80164.42719</td></tr><tr><td rowspan="4">SGXb</td><td rowspan="4">162</td><td>SG_Xb-5</td><td>M1.1</td><td>PSKR2</td><td>A<sub>0</sub></td><td>325</td><td>-148131.1793</td><td>0.000381394</td><td>**</td><td>1</td><td>1</td></tr><tr><td rowspan="3">SG_Xb-6</td><td rowspan="3">M1.1</td><td rowspan="3">PSKR1</td><td>A</td><td>326</td><td>-148124.8687</td><td rowspan="3">0.000004</td><td rowspan="3">**</td><td rowspan="3">1</td><td rowspan="3">1</td></tr><tr><td>A<sub>0</sub></td><td>325</td><td>-148131.1246</td></tr><tr><td>A</td><td>326</td><td>-148120.4449</td></tr><tr><td rowspan="3">SG_XI</td><td rowspan="3">99</td><td>SG_XI-20</td><td>N3.2</td><td>PXC3</td><td>A<sub>0</sub></td><td>199</td><td>-40840.3843</td><td>0.00030735</td><td>**</td><td>0</td><td></td></tr><tr><td rowspan="2">SG_XI-19</td><td rowspan="2">N3.2</td><td></td><td>A</td><td>200</td><td>-40833.87176</td><td rowspan="2">0.000838</td><td rowspan="2">*</td><td rowspan="2">1</td><td rowspan="2">1</td></tr><tr><td>A<sub>0</sub></td><td>199</td><td>-40838.95437</td></tr></tbody></table>
<table><thead><tr><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th><th></th></tr></thead><tbody><tr><td></td><td></td><td>SG_XI-8</td><td>N4</td><td>A<sub>0</sub></td><td>233</td><td>-80679.22227</td><td>0.000000</td><td>**</td><td>9</td><td>3</td><td>6</td></tr><tr><td></td><td></td><td></td><td></td><td>A</td><td>234</td><td>-80660.42399</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SG_XI</td><td rowspan="2">115</td><td>SG_XI-14</td><td>N6.2</td><td>SKM2</td><td>A<sub>0</sub></td><td>231</td><td>-126354.8399</td><td>3.48134E-19</td><td>**</td><td>20</td><td>12</td><td>6</td></tr><tr><td></td><td></td><td></td><td>A</td><td>232</td><td>-126314.768</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>SG_XI-15</td><td>N6.1</td><td></td><td>A<sub>0</sub></td><td>231</td><td>-126374.0677</td><td>0.000000</td><td>**</td><td>12</td><td>8</td><td>3</td></tr><tr><td></td><td></td><td></td><td></td><td>A</td><td>232</td><td>-126335.0153</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SG_XI</td><td rowspan="2">132</td><td>SG_XI-18</td><td>N6.1</td><td>PXY</td><td>A<sub>0</sub></td><td>265</td><td>-131270.9444</td><td>1.13374E-34</td><td>**</td><td>5</td><td>4</td><td>1</td></tr><tr><td></td><td></td><td></td><td>A</td><td>266</td><td>-131195.5224</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>SG_XI-17</td><td>N6.2</td><td>PXL1, PXL2</td><td>A<sub>0</sub></td><td>265</td><td>-131261.1148</td><td>2.96009E-40</td><td>**</td><td>23</td><td>14</td><td>9</td></tr><tr><td></td><td></td><td></td><td></td><td>A</td><td>266</td><td>-131172.9143</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SG_XIIa</td><td rowspan="2">101</td><td>SG_XIIa-2</td><td>O</td><td></td><td>A<sub>0</sub></td><td>203</td><td>-69490.32865</td><td>0.020062</td><td>ns</td><td>-</td><td></td><td></td></tr><tr><td></td><td></td><td></td><td>A</td><td>204</td><td>-69487.62542</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>SG_XIIa-3</td><td>O</td><td></td><td>A<sub>0</sub></td><td>203</td><td>-69490.53986</td><td>0.047369</td><td>ns</td><td>-</td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td>A</td><td>204</td><td>-69488.57373</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SG_XIIla</td><td rowspan="2">143</td><td>SG_XIIla-2</td><td>Q</td><td></td><td>A<sub>0</sub></td><td>287</td><td>-66690.11714</td><td>0.006790</td><td>ns</td><td>-</td><td></td><td></td></tr><tr><td></td><td></td><td></td><td>A</td><td>288</td><td>-66686.45331</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>SG_XIIla-1</td><td>Q</td><td>FEI1, FEI2</td><td>A<sub>0</sub></td><td>287</td><td>-66684.64365</td><td>0.000737</td><td>*</td><td>0</td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td>A</td><td>288</td><td>-66678.94678</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td rowspan="2">SG_XIIlb</td><td rowspan="2">143</td><td>SG_XIIlb-1</td><td>R.2</td><td>ERL1, ERL2</td><td>A<sub>0</sub></td><td>161</td><td>-79161.61197</td><td>3.21923E-09</td><td>**</td><td>26</td><td>14</td><td>8</td></tr><tr><td></td><td></td><td></td><td>A</td><td>162</td><td>-79144.08871</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td>SG_XIIlb-2</td><td>R.1</td><td>ER</td><td>A<sub>0</sub></td><td>161</td><td>-79156.50926</td><td>0.000012</td><td>**</td><td>11</td><td>8</td><td>3</td></tr><tr><td></td><td></td><td></td><td></td><td>A</td><td>162</td><td>-79146.91448</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

lnL: log-likelihood; P-values: ns, Likelihood ratio test not significant; *, significant at 5%-level; **, significant at 1%-level after the Bonferroni correction (see M&M).
## Figure legends

**Fig. 1. Number and expansion of LRR-RLK genes in 33 plant genomes**

The species tree is based on several studies (see Material and methods for details). Branches have been color coded to highlight monocots (red (Poaceae) and pink (others)) versus dicots (yellow (Solanaceae), green (Eurosids I) and blue (Eurosids II)) species. *Selaginella* and moss species are respectively in light and dark brown. **A.** Numbers of LRR-RLK genes found in each proteomes and each SG are indicated. The grey coloring, from light to dark, is proportional to the number of genes found. **B.** Unrooted phylogenetic trees of the three biggest SGs. Note the different expansion patterns in SG_III and SG_XI vs SG_XIIa.

**Fig. 2. Example of monocots/dicots OG characterization in SG_IX.** Branches are color coded according to Figure 1 with branches of monocots (M) species in pink and red, branches of dicots (D) species in yellow, blue and green, and branches of moss and spikemoss in light and dark brown respectively. OGs containing M and D genes are represented as MD OGs and numbered by SG. Within these OGs, the orthologous relationships can be “simple” or “complex” (OG_c). Note that the presence and number of paralogs after monocots and dicots divergence is not taken into account. Numbers at monocots/dicots bifurcations represent nodes with statistical branch supports aLRT/SH-like > 0.85.

**Fig. 3. Distribution of genes per OGs and per SGs.**

**A.** Percentage of genes contained in non complex (OG) and complex OGs (OG_c).

**B.** Number of genes per OGs. The width of boxes are proportional to the number of OGs. Note that SG_III and SG_XI are the largest.

**Fig. 4. Presence/absence of OGs in the Arabidopsis clade.** Presence/absence is represented by white/black boxes respectively.

**Fig. 5. LRR-RLK structural features.** Schematic representation of the subdivided SGs based on ECD characteristics with names of LRR-RLK genes discussed in the article.

**Fig. 6. LRR motifs.** **A.** Schematic representation of the FLS2 orthologs in moss (PHYPA), Selaginella (SELML) and monocots (M) and dicots (D). Note that the number of predicted LRR motifs in PHYPA is lower than in other phyla. **B.** Number of predicted LRRs per protein. The majority of the sequences contain around five, ten or 21 LRRs.

**Fig. 7. Number of sites under positive selection in SG_III, XI or XIIIb per domain.** Schematic representation of a LRR-RLK receptor containing around 20 LRRs.
**Supplementary data**

Supplementary tables, alignments and phylogenetic trees are available online:
http://phylogeny.southgreen.fr/kinase2/

**Suppl-Table 1.** List of species, abbreviation, link to the publication, name of the genome
 fasta file, location for download and version are listed

**Suppl-Table 2.** List of LRR-RLK genes

**Suppl-Table 3.** Detailed structural features of LRR-containing LRR-RLK genes

**Suppl-Table 4.** Positions of predicted LRRs in sequences and islands

**Suppl-Table 5.** Sites under positive selection

Provisional
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Number and expansion of LRR-RLK genes in 33 plant genomes</figcaption></figure>

The species tree is based on several studies (see Material and methods for details). Branches have been color coded to highlight monocots (red (Poaceae) and pink (others)) versus dicots (yellow (Solanaceae), green (Eurosids I) and blue (Eurosids II)) species. *Selaginella* and moss species are respectively in light and dark brown. **A.** Numbers of LRR-RLK genes found in each proteomes and each SG are indicated. The grey coloring, from light to dark, is proportional to the number of genes found. **B.** Unrooted phylogenetic trees of the three biggest SGs. Note the different expansion patterns in SG_III and SG_XI vs SG_XIIa.
<figure><img src="image_3.png" /></figure>

**Figure 2.** Example of monocots/dicots OG characterization in SG_IX. Branches are color coded according to Figure 1 with branches of monocots (M) species in pink and red, branches of dicots (D) species in yellow, blue and green, and branches of moss and spikemoss in light and dark brown respectively. OGs containing M and D genes are represented as MD OGs and numbered by SG. Within these OGs, the orthologous relationships can be “simple” or “complex” (OG_c). Note that the presence and number of paralogs after monocots and dicots divergence is not taken into account. Numbers at monocots/dicots bifurcations represent nodes with statistical branch supports aLRT/SH-like > 0.85.
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>Figure 3. Distribution of genes per OGs and per SGs.</figcaption></figure>

A. Percentage of genes contained in non complex (OG) and complex OGs (OG_c).

B. Number of genes per OGs. The width of boxes are proportional to the number of OGs. Note that SG_III and SG_XI are the largest.
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Presence/absence of OGs in the Arabidopsis clade. Presence/absence is represented by white/black boxes respectively.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** LRR-RLK structural features. Schematic representation of the subdivided SGs based on ECD characteristics with names of LRR-RLK genes discussed in the article.</figcaption></figure>
A

<figure><img src="image_8.png" /></figure>

B

<figure><img src="image_9.png" /></figure>

**Figure 6. LRR motifs.** A. Schematic representation of the FLS2 orthologs in moss (PHYPA), Selaginella (SELML) and monocots (M) and dicots (D). Note that the number of predicted LRR motifs in PHYPA is lower than in other phyla. B. Number of predicted LRRs per protein. The majority of the sequences contain around five, ten or 21 LRRs.
<figure><img src="image_10.png" /><figcaption>**Figure 7.** Number of sites under positive selection in SG_III, XI or XIIb per domain.
Schematic representation of a LRR-RLK receptor containing around 20 LRRs.</figcaption></figure>