Protein sequences were gathered from the Joint Genome Institute (JGI) Phytozome 12 database (https://phytozome-next.jgi.doe.gov) unless specified otherwise (Goodstein et al., 2012). Sequences from Antheros agrestis, Antheros angustus, and Antheros punctatus were gathered from the hornwort sequence database (https://www.hornworts.uzh.ch/en/download.html) of the Szovenyi group. Sequences from Salvinia cucullate were gathered from the Fernbase (https://www.fernbase.org/). Sequences from Picea abies were gathered from the gymnosperm PLAZA plant comparative genomics database (https://bioinformatics.psb.ugent.be/plaza) (Proost et al., 2009).

RLCK XI candidates were identified using the hmmsearch program from the HMMER3 package (Eddy, 2009). The hmmsearch was conducted on 60 proteomes ( Supplementary Table 1 ) across the plant lineage against an RLCK XI profile Hidden Markov Model (HMM). The RLCK XI profile HMM comprised the four known members of RLCK XI from A. thaliana. The resulting 152 sequences that passed an E-value threshold of 0.01 were classified as RLCK XI candidates ( Supplementary Table 2 ).

RLCK XI candidates were annotated for domains using the Conserved Domain Search Tool (https://www.ncbi.nlm.nih.gov/Structure/bwrpsb/bwrpsb.cgi) (Lu et al., 2020). Nuclear localization signals within RLCK XI sequences were identified using the Eukaryotic Linear Motif (ELM) resource (Kumar et al., 2022). Conserved motifs were mapped using the Multiple Em for Motif Elicitation (MEME) online toolbox (http://meme-suite.org/) (Bailey et al., 2015). Transmembrane regions were annotated using Deep Transmembrane Hidden Markov Model (DeepTMHMM) (https://dtu.biolib.com/app/DeepTMHMM/) (Hallgren et al., 2022). Intrinsically disordered regions (IDRs) were analyzed using the IUPred2A program configured to the default (long) setting (Erdos and Dosztanyi, 2020). Sequence logos were built using WebLogo (https://weblogo.berkeley.edu/logo.cgi) (Crooks et al., 2004).

For protein alignment of AthRLCK XI and representative RLKs in Figure 1 , we aligned members of RLCK XI from A. thaliana together with well-studied RLK AthCLV1 from LRR-RLK XI-1 subfamily, AthBIK1 from RLCK VIIa-2 subfamily, and the distant RLK homolog HsaIRAK1 from humans. For Supplementary Figure 1 , we used the RLK dataset from Dievart et al. (2020), which included grouped subfamilies of RLKs from the iTAK database of protein kinases (http://itak.feilab.net/cgi-bin/itak/index.cgi). To identify suitable RLK representatives, each RLK subfamily was aligned with AthRLCK XI-2. The sequence with the highest pairwise comparison score with AthRLCK XI-2 was chosen to represent their subfamily. AthRLCK XI-2 was selected because it has the highest average pairwise similarity score among all AthRLCK XI members. Each RLK representative selected for alignment contained all eleven conserved protein kinase subdomains. RLK subfamilies whose members lacked one or more conserved protein kinase subdomains were excluded from the alignment. To determine the evolutionary conservation of the RLCK XI subfamily in plants, full RLCK XI protein sequences were aligned using the L-INS-i strategy from MAFFT version 7 (Katoh and Standley, 2013). The resulting alignment was used to construct a phylogenetic tree using the maximum likelihood method from IQ-TREE (Nguyen et al., 2015). The Phylogenetic tree was illustrated using iTOL (https://itol.embl.de) (Letunic and Bork, 2021).

The KID, M2, M3, and NLS sequences found in AthRLCK XI-2 were used as queries for BLAST search against the NCBI non-redundant (nr) protein database (Altschul et al., 2005). The BLAST search was performed using the BLASTP algorithm. The E-value threshold was set to 1 to broaden the coverage of the database search due to the lack of prior knowledge involving sequences containing the RLCK XI-KID.

Full-length protein-coding sequences of AthRLCK XI members and different individual domains were PCR amplified, subcloned into pCR8/GW/TOPO vector, and subsequently recombined into the binary vector by LR reaction following the manufacturer’s instructions (constructs are listed in Supplementary Table 3 ). Next, the constructs were transfected into Arabidopsis protoplasts by polyethylene glycol for transient expression (Yoo et al., 2007). Sixteen hours later, transfected cells were imaged using a Zeiss LSM 780 confocal microscope (Germany). Image analysis was performed with LSM software ZEN.

The full-length protein-coding sequences of AthRLCK XI-2 and AthRLCK XI-3 with the GFP tag at the C-terminal end were cloned into pMDC7 vector (Curtis and Grossniklaus, 2003) and transformed into A. thaliana Columbia (Col-0) accession for inducible expression by estradiol. Transgenic plants were generated through an Agrobacterium tumefaciens (GV3101) dipping procedure (Clough and Bent, 1998). The successful transgenic T2 lines were grown on agar plates with half-strength modified Hoagland nutrient solution (Chiou et al., 2006). Ten-day-old seedlings were transferred to the same medium containing 10 μM β-estradiol. The expression of GFP-tagged AthRLCK XI-2 and AthRLCK XI-3 in the root was observed two days after β-estradiol treatment. The confocal microscopy images of roots were taken using an Andor Dragonfly spinning disk confocal microscope with an objective lens PL APO 40X/1.10 W CORR. Excitation/emission wavelengths were 488 nm/500–550 nm for GFP.

Publicly available RNA-seq data were obtained from the Arabidopsis RNA-seq database (ARS) (Zhang et al., 2020) (http://ipf.sustech.edu.cn/pub/athrna/). Data from Arabidopsis tissues and biotic stress response were collected for RLCK XI members from A. thaliana.

Receptor-like kinases possess a well-defined protein kinase domain that has a length of 250-300 amino acids and includes eleven conserved subdomains (I-XI) (Hanks and Hunter, 1995; Shiu and Bleecker, 2001; Lehti-Shiu and Shiu, 2012). There are over 600 RLKs found in the model organism Arabidopsis thaliana that can be classified into more than 50 subfamilies (Shiu and Bleecker, 2001; Shiu and Bleecker, 2003; Dievart et al., 2020). However, only a handful of these RLKs have been studied (Fan et al., 2018; Dievart et al., 2020). Among the RLKs found in A. thaliana, we noticed that the four members of RLCK XI, namely, AthRLCK XI-1 (AT1G80870), AthRLCK XI-2 (AT2G45590), AthRLCK XI-3 (AT4G25390), and AthRLCK XI-4 (AT5G51770), had a longer protein sequence compared to other RLCKs by an average of 280 amino acids. We then investigated the origin of the difference in sequence length by aligning the protein kinase domain of AthRLCK XI members with well-studied RLKs AthCLV1, AthBIK1, and HsaIRAK1 to represent the different structural classifications of RLKs (Shiu and Bleecker, 2001; Dievart et al., 2020). The alignment revealed an unannotated region of low similarity that only aligned within members of RLCK XI ( Figure 1 ). Closer inspection showed that this region resided in the activation segment between the conserved subdomain VII and VIII of the protein kinase domain. The unannotated region within the protein kinase domain of RLCK XI-1 to 4 had an average length of 276 amino acids, compared to the 20-35 amino acids in the activation segment of AthCLV1, BIK1, and HsaIRAK1 in the same region. We refer to this region as the RLCK XI kinase insert domain (RLCK XI-KID). To reinforce the uniqueness of the RLCK XI-KID among RLKs, we performed additional multiple sequence alignment with representatives from several RLK and RLCK subfamilies from A. thaliana. The resulting alignment revealed that none of these other RLCK or RLK subfamilies in A. thaliana had sequences similar to the KID ( Supplementary Figure 1 ). These results suggest that members of AthRLCK XI possess a unique kinase insert domain located between conserved subdomain VII and VIII within their protein kinase domain.

To understand whether the members of AthRLCK XI are conserved, we built a profile Hidden Markov Model (HMM) using AthRLCK XI members to search the Phytozome 12 database for potential homologs in the plant lineage ( Supplementary Table 1 ). A total of 152 RLCK XI candidate sequences were retrieved using the hmmsearch program of the HMMER3 package ( Supplementary Table 2 ). Similar to the naming convention of AthRLCKs, these RLCK XI candidates were assigned names based on their physical locations on the genome of their respective plant species. Annotation of individual protein sequences revealed that all RLCK XI candidates possessed the same domain architecture as those seen in AthRLCK XI members. These results show that RLCK XI is evolutionarily conserved in land plants ( Figure 2 ).

The tree was divided into two major clades. Clade 1 consisted of members from ferns, gymnosperms, basal angiosperms, monocots, and dicots. On the other hand, clade 2 consisted of bryophytes, lycophytes, a basal angiosperm, monocots, and dicots. We identified eight RLCK XI members from the bryophytes Anthoceros agrestis, Anthoceros angustus, Anthoceros punctatus, Marchantia polymorpha, Physcomitrella patens, and Sphagnum fallax. The number of RLCK XI members increased from a single member in bryophytes to two members in basal angiosperms. In turn, the RLCK XI subfamily increased in number in monocots and dicots, with 2-3 members in monocots, and up to 8 members in the genome of dicot Kalanchoe laxiflora ( Figure 2 ). Interestingly, AthRLCK XI-1 was grouped separately from the other AthRLCK XI members in clade 2. RLCK XI members from the same monocot and dicot species could be found in both clades, but RLCK XI members from gymnosperms, ferns, lycophytes, and bryophytes from the same species were exclusively found in one clade but not both. We did not detect RLCK XI members in the chlorophytes, such as Chlamydomonas reinhardtii, Dunaliella salina, Volvox carteri, Coccomyxa subellipsoidea, Micromonas pusilla, and Ostreococcus lucimarinus. Our results show that the RLCK XI subfamily is conserved in land plants. The presence of RLCK XI among land plants suggest that they may play a crucial role in the evolution of the land plant lineage.

To test whether the RLCK XI-KID influenced phylogenetic relationships, we performed the same procedure using only their KID sequences ( Supplementary Figure 2 ). The constructed RLCK XI-KID tree demonstrated similarities to the full-length tree. The RLCK-KID tree remained split into two major clades where lineages and plant species were grouped similarly to the full-length tree. RLCK XI members from monocots and dicots clustered similarly in the RLCK XI-KID tree. However, RLCK XI from M. polymorpha, P. patens, and S. fallax became a sister of ferns, gymnosperms, and angiosperms in clade 1 ( Supplementary Figure 2 ). These results suggest that the RLCK XI-KID is a conserved region among members of RLCK XI.

Additionally, we performed a BLAST search using the KID of AthRLCK XI-2 to examine whether the RLCK XI-KID is present in sequences other than RLCK XI. A total of 1346 hits that included both hypothetical and predicted proteins were identified and all of them could be traced as members of the RLCK XI subfamily based on their sequence length and coverage of the AthRLCK XI-2 KID ( Supplementary Table 4 ). These results show that the RLCK XI-KID is exclusively found in members of RLCK XI within the plant lineage.

We observed that the length of the RLCK XI-KID has been dynamic in the course of land plant evolution, ranging from 199-527 amino acids, with 280 amino acids on average. The longest KID belonged to MpoRLCK XI-1, which is the only RLCK XI homolog found in M. polymorpha ( Supplementary Table 2 ).

Previous studies have shown that the activation segments of protein kinases function as flexible contact sites with high intrinsic disorder relative to the surrounding protein kinase core (Gogl et al., 2019). We thus examined whether the RLCK XI-KID also exhibits this property. The full-length sequences of RLCK XI members were used to predict the distribution of IDRs through IUPred2A (Dosztanyi et al., 2005). Based on the propensity for disordered regions, the kinase domain-containing conserved subdomains I-VII and VII-XI have a low propensity for being disordered. In contrast, the KIDs clearly show a high tendency for disorder among most RLCK XI members ( Figure 3 ). This result reveals that, like the activation segment of protein kinases, RLCK XI-KIDs also presented qualities found in IDRs.

Next, we wanted to know whether conserved motifs existed within the RLCK XI-KID. We observed that the 152 RLCK XI-KID sequences we collected had an average sequence identity of 40%. We examined these sequences for ungapped sequence motifs using the Multiple Em for Motif Elicitation (MEME) suite (http://meme-suite.org/) (Bailey et al., 2015) as a means to determine the origin of the observed sequence similarity. We found twelve conserved motifs with E-values of less than 0.01 across all RLCK XI-KID sequences ( Figure 4A ). These motifs were numbered from lowest E-value (M1) to highest (M12). Interestingly, sub-motif D/EW and D/EWW variations were observed in motifs M2, M3, and M6, while motifs with polylysine residues were detected in M4 and M11. We used the D/EWW-containing M2 and M3 sequences of AthRLCK XI-2 to perform a BLAST search against the NCBI non-redundant protein database to determine whether these motifs could be found outside of RLCK XI. We only found M2 or M3-containing sequences belonging to members of RLCK XI as determined by sequence homology search ( Supplementary Tables 5 , 6 ). These results suggest that the conserved motifs within the KID are unique to members of RLCK XI.

In our motif search, we observed that each of the 12 motifs identified only occurred once in all 152 RLCK XI-KID sequences. We wanted to investigate whether these motifs were positioned at random or arranged in a discernable pattern along the RLCK XI-KID. We divided the RLCK XI-KID into twelve positions (P1-P12) based on the single occurrence of each of the twelve motifs in all RLCK XI-KID sequences. If a discernable pattern exists in the KID, then each position would reflect a high frequency of specific motifs in that position. When we explored the position or arrangement of these motifs ( Figure 4B ), we found that M8 and M9 were positioned most frequently at the N- and C-terminal ends of the RLCK XI-KID. Interestingly, M3 and M4 had the highest frequency in multiple positions. M3 had the highest frequency in P3, P4, and P5. Similarly, M4 had the highest frequency in P9, P10, and P11. We noticed a prevalence of charged residues in M1, M2, M3, M4, and M11. On the other hand, M10 was composed of mostly serine residues. M10, M11, and M12 were found across multiple positions on the KID. We hypothesize that the prevalence of M10, M11, and M12 was in part due to the short sequence that made up their respective motifs. These results show that conserved motifs within the RLCK XI-KID have distinct sequence and spatial properties. The conservation of these motifs implies the involvement of potential functions.

Disordered regions have been shown to contain a high frequency of interfaces for protein-protein interaction. Some of these sites exhibit both regulatory and signaling functions and are referred to as eukaryotic linear motifs (ELMs) (Kumar et al., 2022). ELMs are evolutionarily malleable sequences and interact with relatively low affinity due to the limited number of residues that make direct contact with the binding partner (Kumar et al., 2022). These features confer ELMs with the ability to mediate transient interactions to maintain robust cell signaling (Davey et al., 2012; Kumar et al., 2022).

To explore the functional relevance of RLCK XI-KID, we searched the ELM resources (Kumar et al., 2022) to examine whether the KIDs of four Arabidopsis RLCK XI members contain any ELMs. We identified several ELMs, annotated as “TRG_NLS_Bipartite_1”, “TRG_NLS_MonoCore_2”, “TRG_NLS_MonoExtC_3” or “TRG_NLS_MonoExtN_4” ( Supplementary Table 7 ), which belong to bipartite variants of the classical basic and charged nuclear localization signal (NLS). Interestingly, these NLSs repeatedly occur in the KID ( Figure 5A ) and overlap with the conserved motifs M2 and M3, suggesting that the conserved motifs found in RLCK XI members could be functionally relevant. Similarly, the BLAST search for sequences containing the NLS that overlapped M2 and M3 were members of RLCK XI in plants ( Supplementary Table 8 ). In addition to the KID, NLS could be identified in the N and/or C termini of several AthRLCK XI members except for AthRLCK XI-1 ( Figure 5A ). Similar to Arabidopsis, two RLCK XI members from rice also contain NLSs in their KIDs.

To examine whether the observed NLSs in the KID of AthRLCK XI members are functional, we tagged the protein with a green fluorescent protein (GFP) at the C-terminus. We performed a transient expression assay in Arabidopsis protoplasts to observe their subcellular localization using confocal microscopy. While the free GFP was distributed in the cytosol and nucleus, we found the GFP-tagged KIDs of all four AthRLCK XI members and one OsaRLCK XI member were strictly confined in the nucleus ( Figure 5 ). These results indicate that the predicted NLSs in the KID of RLCK XI are functional in mediating nuclear import.

The nuclear localization of the KID in AthRLCK XI members was unexpected due to the presence of a single predicted transmembrane domain in all four members. To validate whether RLCK XI members are membrane-bound RLKs, we transiently expressed the full-length AthRLCK XI tagged with GFP at the C-terminus and observed their subcellular localization in Arabidopsis protoplasts. In the following experiments, we only focused on AthRLCK XI-2 and -3 because we were unable to clone the full-length transcripts of AthRLCK XI-1 and -4 due to unknown reasons. Transient expression of full-length AthRLCK XI-2-GFP and -3-GFP revealed their localization on the plasma membrane and in the nucleus with a lesser signal in the cytosol ( Figures 6B, D ). On the other hand, protoplasts expressing AthRLCK XI-2 and -3 without a transmembrane region (ΔTM-GFP) showed localization signals exclusively in the nucleus, similar to the GFP-tagged KID sequence. However, expression of AthRLCK XI-2 and -3 without the KID (ΔKID-GFP) resulted in signals on the plasma membrane and cytosol, where the signals often formed aggregates ( Figures 6B, D ). Intriguingly, expression of AthRLCK XI-2 and -3 that lacked a transmembrane region and KID (ΔTMΔKID-GFP) resulted in punctate signals in the nucleus (AthRLCK XI-2) or the cytosol with less uniform nuclear signals (AthRLCK XI-3) ( Figures 6B, D ). We presume this nuclear localization may result from the predicted NLS sequences in the C-terminal region of AthRLCK XI-2 and -3. Indeed, further removal of the C-terminus (ΔTMΔKIDΔC-GFP) resulted in punctate signals in the cytosol and none in the nucleus ( Figure 6A ). These results suggest KID is essential for the proper nuclear localization of AthRLCK XI-2 and -3. Consistent with the results observed in the protoplasts, full-length AthRLCK XI-2 and -3 mainly localized on the plasma membrane and in the nucleus of transgenic Arabidopsis root tips when induced by β-estradiol ( Figures 6C, E ).

The observation of dual localization of full-length AthRLCK XI-2 and -3 on the plasma membrane and nucleus raises a possibility that the plasma membrane-anchored N-terminal portions of AthRLCK XI-2 and -3 may be proteolytically cleaved, triggered by an unknown signal. As a result, the transmembrane domain-truncated C-terminal portion of AthRLCK XI is redirected to the nucleus via the NLS signals in the KID to transmit and relay the signal. It is worth noting that although NLS sequences could be found in the KID and the C-terminal end of AthRLCK XI-2 and -3, the C-terminal NLS sequences alone were insufficient to mediate the nuclear localization observed in full-length AthRLCK XI-2 and -3 ( Figure 6A ).

Members of the RLCK XI subfamily from A. thaliana have been classified in previous literature as part of large-scale phylogenetic studies aiming to place RLKs as a whole in the protein kinase superfamily (Shiu and Bleecker, 2001; Shiu and Bleecker, 2003; Lehti-Shiu et al., 2009). However, to our knowledge, there has been no specific report of a KID in any of them. Full-length protein alignment of representative protein kinases alongside members of RLCK XI clearly shows the presence of the KID with an average length of 280 amino acids that splits the conserved subdomains I-VII from subdomains VIII-XI ( Figure 1 and Supplementary Figure 1 ). In protein kinases, this region corresponds to the activation segment, which is 20-35 amino acids long and is less conserved within the kinase domain (Nolen et al., 2004). The activation segment, usually exhibiting qualities similar to IDRs, functions as a site for domain interaction and for regulating kinase activity (Taylor and Kornev, 2011; Gogl et al., 2019). KIDs were previously described as anomalous sequences exhibiting qualities present in IDRs with little sequence similarity between protein kinase domains of RTKs (Locascio and Donoghue, 2013; Koike et al., 2020). Although RLCK XI KIDs exhibited properties similar to IDRs ( Figure 3 ), the presence of several conserved motifs that overlapped with predicted NLS within this region suggests that the RLCK XI KID has a role in RLCK XI function ( Figure 4 and Supplementary Table S8 ). Indeed, classical NLS have often been found in IDRs of cargo proteins (Wubben et al., 2020). Furthermore, KIDs found in RTKs contained phosphorylation sites that regulated protein function (Locascio and Donoghue, 2013).

It would be worth exploring whether the KID affects RLCK XI protein kinase activity. To conduct the in vitro kinase activity assay, we attempted to produce recombinant proteins of AthRLCK XI-2 and AthRLCK XI-3 without the N-terminal transmembrane domain (ΔTM) by expressing them in E. coli. Unfortunately, the recombinant proteins were not expressed for unknown reasons. We suspect the intrinsically disordered property of the KID may hinder their expression in E. coli. Nevertheless, we could glean some insight from previous studies. First, we’ve observed that members of AthRLCK XI have all the subdomains (from subdomains I to XI) found in functional protein kinases. In contrast, pseudokinases have been shown to lack one or more of these subdomains in their protein sequence (Kannan and Taylor, 2008; McClendon et al., 2014; Kwon et al., 2019). Second, the location of the KID corresponds to the activation segment, which contains regulatory regions and phosphorylation sites that regulate protein kinase activity (Nolen et al., 2004; Taylor et al., 2012). We hypothesize that the RLCK XI-KID has the potential to influence protein kinase activity due to being in the same region as the activation segment in other protein kinases. Interestingly, members of RLCK XI are also classified as non-RD protein kinases. This property refers to protein kinase subdomain VI, which follows the HxD motif wherein x could be R or another amino acid (Kornev et al., 2006). RD protein kinases have been shown to require autophosphorylation on the activation loop for regulation (Adams, 2003; Steichen et al., 2010). In contrast, non-RD protein kinases can be regulated through mechanisms other than autophosphorylation (Dardick et al., 2012; Bender et al., 2021). The presence of the KID and the non-RD subdomain VI in RLCK XI may indicate that these protein kinases are regulated through an alternate mechanism from phosphorylation within the activation segment.

The RLK family of protein kinases has expanded in land plants, which resulted in the various subfamilies observed in present-day plant species (Shiu and Bleecker, 2003; Lehti-Shiu and Shiu, 2012; Dievart et al., 2020). Recent studies have identified sequences belonging to the RLK family in glaucophytes (Gong and Han, 2021). However, we have only traced members of RLCK XI up until bryophytes. The conservation and small number of the RLCK XI subfamily is similar to previously described RLK subfamilies found in plants (Lehti-Shiu and Shiu, 2012). However, none of these other RLK subfamilies have the KID ( Supplementary Figure 1 ). Having identified the unique domain architecture in members of RLCK XI, we speculate that during an event in the course of land plant evolution, likely aquatic-to-land plant transition, this RLK acquired the KID we observe today. How the ancestral RLK acquired the KID requires further study. Nevertheless, their conservation among land plants suggests they may have played a role in land plant evolution.

A single-pass transmembrane region in the AthRLCK XI sequences suggested their localization on the plasma membrane ( Figure 5A ). However, detecting multiple NLSs within the RLCK XI-KID, the juxtamembrane region, and their C-terminal sequence also suggested their nuclear localization ( Figure 5A ). Different truncated variants of RLCK XI proteins ( Figure 6A ) proved that the NLS were functional in directing their nuclear localization. The dual localization of full-length AthRLCK XI-2 and -3 on the plasma membrane and in the nucleus further verified the functionality of both the transmembrane region and NLS ( Figures 6B–E ). We attempted to tag either the N-terminus or C-terminus of full-length AthRLCK XI-2 and AthRLCK XI-3. However, the N-terminal GFP-tagged protein resulted in irregular aggregates in the cytosol of protoplasts. We speculate that N-terminal GFP tagging of AthRLCK XI disrupt their targeting on the plasma membrane. We, therefore, only showed the result of the C-terminal GFP-tagged proteins. In the future, adding a different fluorescent protein tag downstream of the transmembrane domain but upstream of the kinase domain of the C-terminal GFP-tagged AthRLCK XI may allow us to track changes in their subcellular localization.

The dual localization of AthRLCK XI-2 and -3 on the plasma membrane and in the nucleus suggests that these proteins may undergo proteolytic cleavage before translocating to the nucleus. Indeed, proteolytic cleavage of plant RLKs to release the intracellular fragment has been reported. In rice, the LRR-RLK Xa21 confers broad-spectrum resistance to Xanthomonas oryzae pv. Oryzae (Song et al., 1995). Xa21 undergoes proteolytic cleavage, leading to the translocation of its intracellular kinase domain to the nucleus. Xa21 has been shown to bind WRKY62 in the nucleus of rice protoplasts (Peng et al., 2008; Park and Ronald, 2012). Several RLCK members were also shown to undergo translocation from plasma membranes to the nucleus as a part of their signaling cascade. BIK1 is translocated to the nucleus and interacts with WRKY transcription factors that regulate Jasmonic acid and salicylic acid levels (Lal et al., 2018). Similarly, chitin perception-induced PBS1-LIKE 19 (PBL19) can be translocated to the nucleus, where it interacts with WRKY8 to regulate its own transcription (Li et al., 2022). In another case, the RLCK Puccinia striiformis-Induced Protein Kinase 1 (PsIPK1) found in wheat interacts with the fungal effector PsSpg1 from Puccinia striiformis to promote its nuclear localization to phosphorylate CCAAT-binding transcription factor 1d (TaCBF1d) and enhance susceptibility to the pathogen (Wang et al., 2022). Since RLKs are widely recognized as receptors and mediators for plant development and stress response, we speculate that upon being triggered by a specific signal, RLCK XI may be cleaved by a protease and undergo translocation from the plasma membrane to the nucleus. RLCK XI may function to transmit the extracellular signal to the nucleus for a signaling cascade. The signal and protease responsible for inducing this translocation mechanism have yet to be determined.

It would be fascinating to know the physiological role of RLCK XI members. Unfortunately, to the best of our knowledge, there have been no prior studies on them. According to the publicly available expression profiles (Arabidopsis RNA-seq database, http://ipf.sustech.edu.cn) (Zhang et al., 2020), all the members of AthRLCK XI were expressed ubiquitously except in the pollen ( Supplementary Figures 3A–D ). Members of AthRLCK XI were also expressed under different biotic stresses with some variations ( Supplementary Figures 3E–H ). We attempted to characterize the T-DNA insertion Arabidopsis mutants of all four AthRLCK XI members individually, but none of the single mutants showed developmental defects. Higher-order mutants or specific stress conditions may be employed to examine their function in the future. On the other hand, the transgenic Arabidopsis lines with inducible expression of AthRLCK XI ( Figures 6C, E ) will be excellent materials for functional characterization in the future.

In summary, although the function of the RLCK XI members remains elusive, our characterization of RLCK XI in this study has expanded on the structural diversity known for the protein kinase domain and the assortment of evolutionary changes undergone by receptor-like kinases in plants.