ORIGINAL RESEARCH article

Front. Plant Sci., 13 August 2021

Sec. Plant Systematics and Evolution

Volume 12 - 2021 | https://doi.org/10.3389/fpls.2021.650171

A Comprehensive Phylogenetic Analysis of the MAP4K Family in the Green Lineage

  • 1. Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium

  • 2. VIB Center for Plant Systems Biology, Ghent, Belgium

Article metrics

View details

9

Citations

3,9k

Views

1,4k

Downloads

Abstract

The kinase-mediated phosphorylation impacts every basic cellular process. While mitogen-activated protein kinase technology kinase kinases (MAP4Ks) are evolutionarily conserved, there is no comprehensive overview of the MAP4K family in the green lineage (Viridiplantae). In this study, we identified putative MAP4K members from representative species of the two core groups in the green lineage: Chlorophyta, which is a diverse group of green algae, and Streptophyta, which is mostly freshwater green algae and land plants. From that, we inferred the evolutionary relationships of MAP4K proteins through a phylogenetic reconstruction. Furthermore, we provided a classification of the MAP4Ks in the green lineage into three distinct.

Introduction

Post-translational modifications (PTMs) are vital for plants to sense and respond to environmental changes (Hashiguchi and Komatsu, 2016; Vu et al., 2018). Among these, reversible protein phosphorylation is one of the most widespread and pivotal PTMs that affect every basic cellular process (Ardito et al., 2017; Vu et al., 2018). Phosphorylation is catalyzed by protein kinases, which predominantly phosphorylate substrate proteins on serine, threonine, or tyrosine residues, and protein phosphatases, which mediate the reverse reaction. One of the well-explored, highly conserved kinase subfamilies are components of the mitogen-activated protein kinase (MAPK) cascades, which act through linear sequential serine/threonine and/or tyrosine phosphorylation (Xu and Zhang, 2015). An MAPK signaling module is composed of an MAPK kinase kinase (MAP3K), a MAPK kinase (MAP2K), and a MAPK (Jagodzik et al., 2018). Typically, a MAPK cascade regulates signal transduction through a MAP3K that phosphorylates a MAP2K, which then activates a MAPK by phosphorylation (Xu and Zhang, 2015; Krysan and Colcombet, 2018). Plant MAPK cascades act downstream of receptor-like protein kinases or G-proteins and play important roles in immunity, abiotic stress, and plant growth and development (Cho et al., 2008; Cheng et al., 2015; Xu and Zhang, 2015; Liu and Zhou, 2018; Yan et al., 2018; Zhu et al., 2019).

However, the yeast Ste20 acts as a MAPK kinase kinase kinase (MAP4K) that activates Ste11, a MAP3K, which, in turn, activates downstream components in the mating signaling pathway (Leberer et al., 1992). Interestingly, MAP4Ks are conserved not only in yeast but also in mammals and plants (Dan et al., 2001; Champion et al., 2004; Chuang et al., 2016). The MAP4Ks in plants were first identified in Brassica (Leprince et al., 1999), but, so far, plant MAP4Ks have not been extensively characterized (Pan and De Smet, 2020). In Arabidopsis thaliana, 10 MAP4Ks were described (Jonak et al., 2002), but not all MAP4Ks have been functionally characterized: BLUE LIGHT SIGNALING 1 (BLUS1)/MAP4K10 is essential for stomatal opening in response to blue light (Takemiya et al., 2013; Schnabel et al., 2018), SALT INDUCIBLE KINASE 1 (SIK1)/MAP4K3 is critical for cell proliferation and expansion during organ growth and development (Xiong et al., 2016), and both SIK1 and MAP4K4 regulate the flg22-triggered immunity response (Jiang et al., 2019). In addition, a TARGET OF TEMPERATURE 3 (TOT3)/MAP4K4 interacts with TOT3-INTERACTING PROTEIN 4 (TOI4)/MAP4K6 and TOI5/MAP5K5 and controls thermomorphogenesis in Arabidopsis and wheat (Vu et al., 2021).

Although that MAP4Ks are so conserved throughout multiple clades, there is no comprehensive overview of the MAP4K family in the green lineage (Viridiplantae). Viridiplantae is a monophyletic clade of photosynthetic eukaryotic organisms that play important roles in both terrestrial and aquatic ecologies (Leliaert et al., 2012). The group separated around 1 billion years ago into Chlorophyta and Streptophyta, the latter including Embryophyta (the clade of the land plants) (Morris et al., 2018; Fürst-Jansen et al., 2020). One major evolutionary milestone was the conquest of land by the Embryophyta, which is from a freshwater streptophyte algal common ancestor, followed by an expansion in biodiversity (Becker and Marin, 2009). The habitat shift was accompanied by the modification of a series of traits that granted success in land colonization (e.g., the development of specialized organs and tissues such as stomata, roots, leaves, seeds, and flowers), but also a set of molecular innovations (Pires and Dolan, 2012; Ishizaki, 2017).

In this study, we identify putative MAP4K members from representative species of Chlorophyta and Streptophyta, the two clades in the Viridiplantae, and infer the phylogenetic relationships in the green lineage underpinning MAP4K evolution and diversification.

Materials and Methods

Genome-Wide Identification of (Putative) MAP4Ks

To identify and reconstruct the evolution of putative MAP4Ks throughout the green lineage (Viridiplantae), 6 Chlorophytes, and 29 Streptophytes (including 24 Embryophytes) were chosen depending on the genome availability and annotation status (Supplementary Table 1). The MAP4K sequences were identified using a hidden Markov model (HMM) profile generated and calibrated using the HMMER software v3.3 (Eddy, 1998) that was trained based on an alignment of the 10 MAP4Ks previously identified in A. thaliana (Supplementary Figure 1; Jonak et al., 2002). Local searches were performed on primary protein libraries (Leebens-Mack et al., 2019). To maximize the balance between sensitivity and selectivity, only hits with both per-sequence and per-domain E-values lower than 2e-100 and primary protein isoforms were considered for downstream analysis. The sequences retrieved were then further confirmed for the presence of the protein kinases domain (PF00069) in the Pfam database (El-Gebali et al., 2019) and InterPro Scan v5 (Jones et al., 2014) and inspected considering the annotation status for some of the libraries; the 5% longest and shortest were reviewed using a homology-based method (Korf et al., 2001). For this process, the canonical genomic structure [(e.g., evolutionarily conserved exons, introns, and untranslated regions (UTRs)] of the best reciprocal blast hits among A. thaliana (Streptophyta-dicotyledon), Oryza sativa (Streptophyta-monocotyledon), and Chlamydomonas reinhardtii (Chlorophyta) were compared for validation (Altschul et al., 1990, 1997; Moreno-Hagelsieb and Latimer, 2008).

Alignment and Phylogenetic Analyses

To build a phylogenetic reconstruction for the MAP4K group, a MAPK family member (MPK3, AT3G45640.1) was first added to the list of putative MAP4Ks as a possible outgroup. Following this, full-length sequences were aligned using the GISMO software, which applies a top-down sequence alignment and is optimized for protein sampling (Neuwald and Altschul, 2016) (Supplementary Data Sheet 1). Subsequently, we inferred the maximum likelihood phylogenetic tree using the concatenated alignment under the best fitting model JTT + I + G4 in IQ-TREE v2.0.6 (Neuwald and Altschul, 2016), executing 1,000 ultrafast replicates. The resulting Newick consensus tree was visualized using the Interactive Tree Of Live version 4 (iTOL) (Letunic and Bork, 2019). The taxonomy tree was generated with the phyloT online tool (phylot.biobyte.de). In addition, a Bayesian Markov chain Monte Carlo (MCMC) sampling method for phylogenetic reconstruction was used. The alignment was submitted to the best fitting model JTT + G4 implemented in PhyloBayes with 5,000 generations and, after the convergence of the two chains, a burn-in of 25% was applied. The resulting consensus tree (data not shown) was inspected and compared with the previous tree, revealing a high conservation of the three clades (100%, 100%, and 99.41% for clade I, II, and II, respectively) and high posterior probability support (1,0.83, and 1 for Clade I, II, and III, respectively), but several unresolved polytomies.

Feature Search and Characterization

To gain insight into the conserved motifs among plant MAP4K sequences, the MEME suite 5.3.31 (Bailey and Gribskov, 1998) was used. The number of discoverable motifs was set to 15 and the width of motifs was set to an interval of a minimum of 3 and maximum of 50 amino acids. To create sequence logos, WEBLOGO2 was used with the default settings (Crooks et al., 2004). The analysis of disorder was performed by the web-based IUPred2A tool3 (Mészáros et al., 2018) utilizing the A. thaliana MAP4K protein sequences. The resulting disorder probability was used to categorize each residue as either ordered (<0.4), intermediate (0.4–0.6), or disordered (>0.6). Analyses of the biotic and abiotic influences on A. thaliana MAP4K expression was done through analyses of available transcriptome data on the Genevestigator platform4 (Zimmermann et al., 2004, 2005; Hruz et al., 2008) on August 25, 2020. The eFP browser5 (Hayes et al., 2008) with standard settings was used for the cell, tissue, and organ MAP4K expression analyses.

Results and Discussion

The MAP4K Family Members in A. thaliana

In A. thaliana, the MAP4K family consists of 10 members, namely, MAP4K1 (At1G53165), MAP4K2 (At3G15220), MAP4K3/SIK1 (At1G69220), MAP4K4/TOT3 (At5G14720), MAP4K5/TOI5 (At4G24100), MAP4K6/TOI4 (At4G10730), MAP4K7 (At1G70430), MAP4K8 (At1G79640), MAP4K9 (At1G23700), and MAP4K10/BLUS1 (At4G14480) (Pan and De Smet, 2020; Vu et al., 2021). Most A. thaliana MAP4K family members share common features, such as a MAP4K signature motif and several conserved residues in the kinase domain, and display a C-terminal half from the end of the kinase domain to the stop codon that is largely intrinsically disordered (Supplementary Figure 1). However, some differences can be observed. For example, SIK1 contains a long intrinsically disordered N-terminal part before the kinase domain, which is important for interaction with MOB1 to control cell expansion (Xiong et al., 2016; Supplementary Figure 1). Furthermore, MAP4K9 lacks the HRD and DFG motifs in the kinase domain that are vital for catalytic activity, does not contain an intrinsically disordered C-terminal part, and does not display the conserved MAP4K signature and GXGXXG/A motifs (Supplementary Figure 1). In Homo sapiens, MAP4Ks contain a characteristic, highly conserved C-terminal citron-homology domain (CNH) in all seven known members of the group (Chuang et al., 2016; Seo et al., 2020), which is also present in MAP4K orthologs in Drosophila melanogaster and Caenorhabditis elegans (Fiedler et al., 2014), which is important for protein–protein interactions (Poinat et al., 2002; Machida et al., 2004; Taira et al., 2004). Strikingly, a fully conserved CNH domain is absent in A. thaliana MAP4Ks.

To gain insight into the expression patterns of the A. thaliana MAP4K family during development, the MAP4K expression data were compiled and visualized from online repositories, namely, eFP browser (Winter et al., 2007) and Genevestigator v3 (Hruz et al., 2008). These in silico expression patterns suggested that most MAP4Ks are broadly expressed throughout the plant, with some MAP4Ks expressed stronger or specifically in one organ, such as in pollen (Supplementary Figure 2). This indicates distinct roles for some of these kinases. Furthermore, the MAP4K expression was hardly affected by hormones or nutrients and was mainly regulated by abiotic and biotic stresses (Supplementary Table 2).

Identification of Putative MAP4K Family Members in the Green Lineage

While there is some information on the MAP4K family members in a few plant species, such as A. thaliana, Brassica napus, Zea mays, Solanum chacoense, and Vitis vinifera (Leprince et al., 1999; Llompart et al., 2003; Major et al., 2009; Cakir and Kılıçkaya, 2015; Pan and De Smet, 2020), a comprehensive overview is lacking. Therefore, to reconstruct the evolutionary relationships of the MAP4K family in the green lineage, we analyzed 34 representative species with available genome or protein information (Tuskan et al., 2006; Jaillon et al., 2007; Merchant et al., 2007; Ouyang et al., 2007; Palenik et al., 2007; Worden et al., 2009; Blanc et al., 2010; Prochnik et al., 2010; Schmutz et al., 2010; Banks et al., 2011; Hu et al., 2011; Young et al., 2011; Lamesch et al., 2012; Sato et al., 2012; Albert et al., 2013; Motamayor et al., 2013; Nystedt et al., 2013; Sharma et al., 2013; Blanc-Mathieu et al., 2014; Hori et al., 2014; Liu et al., 2014; Hirsch et al., 2016; Beier et al., 2017; Bowman et al., 2017; Appels et al., 2018; De Clerck et al., 2018; Lang et al., 2018; Li et al., 2018; McCormick et al., 2018; Nishiyama et al., 2018; Wan et al., 2018; Cheng et al., 2019; Wang et al., 2020; Zhang et al., 2020; Supplementary Table 1). In this analysis, we included representative green algae from Chlorophyta [Chlamydomonas reinhardtii and Volvox carteri (Chlorophyceae), Micromonas pusilla (Mamiellophyceae), Ostreococcus tauri and Chlorella variabilis (Trebouxiophyceae), and Ulva mutabilis (Ulvophyceae)] and Streptophyta [Klebsormidium nitens (Klebsormidiophyceae), Chlorokybus atmophyticus (Chlorokybophyceae), Mesostigma viride (Mesostigmatophyceae), Mesotaenium endlicherianum (Zygnematophyceae), and Chara braunii (Charophyceae)]. In addition, several land plants from the Streptophyta were included to trace the expansion of MAP4Ks in land plants: the bryophyte Physcomitrium patens, the liverwort Marchantia polymorpha, the hornwort Anthoceros agrestis, the lycophyte Selaginella moellendorffii, the ferns Azolla filiculoides and Salvinia cucullata, the gymnosperms Picea abies and Gnetum montanum, and 16 angiosperms (including Amborella trichopoda, 10 eudicots, and 5 monocots). Through in-depth analyses of protein-containing databases (see section “Materials and Methods”), we generated a set of 249 putative MAP4K family members with a minimum and maximum amino acid sequence length of 225 and 2,251, respectively (Supplementary Table 3 and Supplementary Data Sheet 2). These results showed that putative members of the MAP4K family are present in both Chlorophyta and Streptophyta and all representative species used in this analysis (Figure 1 and Supplementary Tables 1, 3). The resulting phylogenetic consensus tree revealed three major clades (Figure 1). With respect to clade III, four subclades could be distinguished in angiosperms: one subclade with AtMAP4K8 (Subclade IIIA), a second subclade with AtMAP4K4, AtMAP4K7, and AtMAP4K9 (Subclade IIIB), the third subclade with AtMAP4K5 and AtMAP4K6 (Subclade IIIC), and the fourth clade with AtMAP4K10 (Subclade IIID) (Figure 1 and Supplementary Table 3).

FIGURE 1

FIGURE 1

Combined phylogenetic tree and motif locations of the green lineage MAP4Ks. (A) Inferred phylogenetic tree of putative green lineage mitogen-activated protein kinase MAP4Ks among representative species. The maximum likelihood tree was inferred from the multiple alignments of 239 putative MAP4K sequences and 10 previously described Arabidopsis thaliana MAP4Ks sequences across 34 Viridiplantae primary proteomes using IQ-Tree. Three main clades are indicated: I, II, and III. Clade III contains additional subclades. IIIA, IIIB, IIIC, and IIID. Bootstrap values > 24% are indicated in the branches. Red dots indicate Arabidopsis MAP4Ks. The Arabidopsis MAP4Ks included in the (sub)clades are MAP4K1 (I), MAP4K2 (I), MAP4K3/SIK1 (II), MAP4K4 (IIIB), MAP4K5 (IIIC), MAP4K6 (IIIC), MAP4K7 (IIIB), MAP4K8 (IIIA), MAP4K9 (IIIB), and MAP4K10/BLUS1 (IIID). (B) Motif distribution along the putative MAP4K sequences identified by MEME tool. The height of motif “block” is proportional to -log (p-value), the taller the block, and the lower the probability of a wrong match. The motif sequence is indicated in the figure.

Conserved Features of the MAP4Ks in Plants

As proteins are generally composed of one or more functional regions or domains that can provide insight into their function and evolutionary relationships (Lee et al., 2007), we explored to what level the key domains and motifs are conserved in the putative MAP4K family members (Figures 1, 2A, Supplementary Figure 3, and Supplementary Table 3).

FIGURE 2

FIGURE 2

MAP4K features. (A) Sequence logos of main conserved motifs. (B) Length of C-terminal half for all MAP4Ks and indicated clades.

First, we assessed to what extent the identified sequences have a likely functional kinase domain based on some key features, such as a glycine-rich loop motif (GxGxxG) that is essential for nucleotide binding and is an integral part of the ATP-binding site. In general, in about 20% of kinases, the third G is substituted by A or S, which disrupts kinase activity (Hanks and Quinn, 1991; Hemmer et al., 1997; Torkamani et al., 2008; Chmielecki et al., 2010; Steinberg, 2018; Zhang et al., 2018). We also observed this in several MAP4Ks (Supplementary Table 3 and Supplementary Figure 3). In addition, we analyzed HRD and DFG motifs and an invariant lysine (K), all of which are pivotal for catalytic activity (Hanks et al., 1988; Carrera et al., 1993; Hanks and Hunter, 1995; Xu et al., 2000; Roskoski, 2004; Kannan and Neuwald, 2005; Zhang et al., 2015; Modi and Dunbrack, 2019; Pan and De Smet, 2020; Supplementary Table 3 and Supplementary Figure 3). Due to changes in these amino acid triads, the putative MAP4Ks that lack all or some of the above-mentioned kinase features are possibly kinase-inactive MAP4Ks (Supplementary Table 3; Paul and Srinivasan, 2020).

Second, we evaluated the presence of the previously proposed signature motif for (mammalian) MAP4Ks (GTPyWMAPEv, with Y and V being less conserved) located in the kinase subdomain VIII (Sells and Chernoff, 1997; Dan et al., 2001). The majority (94%) of the putative MAP4Ks in the green lineage share this well-conserved signature motif (Supplementary Table 3 and Supplementary Figure 3). Interestingly, two A. thaliana and A. lyrata MAP4K9s from subclade IIIB do not contain this conserved signature motif and show severe alterations in the HRD, DFG, and glycine-rich loop motifs (GxGxxG) (Supplementary Table 3), in addition, which raises the question if these should still be considered as true (functional) MAP4Ks.

Third, we investigated the length of the C-terminal half of the MAP4Ks. Although the disordered C-terminal half is essential for the interaction with substrates, it is largely not conserved across sequences (Lee et al., 2017). For almost 70% of the putative MAP4Ks, the length of the C-terminal half was between 300 and 600 amino acids (Figure 2B, Supplementary Figure 4, and Supplementary Table 3). Putative MAP4Ks that (partially) lack such a C-terminal half, such as SMO169G0133.1 in Selaginella, TraesCS6A02G353400.1 in wheat, and Bo9g089620.1 in Brassica, will likely be affected in the interaction with substrates (Zhang et al., 2018; Jiang et al., 2019).

Strikingly, none out of the 249 sequences contained a fully conserved CNH domain, suggesting that this characteristic element for several non-plant MAP4Ks might not have been present in the ancestor of the Viridiplantae. However, we identified conserved but clade-specific motifs at a similar position where the CNH domain is found in other eukaryotes (Figure 1, Supplementary Figures 15, and Supplementary Table 3). Furthermore, an experimental analysis will be required to reveal whether these motifs are involved in similar functions as the CNH domain, such as protein–protein interaction.

Finally, the phototropin-mediated phosphorylation of BLUS1 Ser-348 alleviates the autoinhibitory activity of the C-terminal part on its kinase activity (Hosotani et al., 2021). This serine residue belongs to a conserved RRISGWNF consensus motif (Takemiya et al., 2013). Therefore, we checked for the presence of such a motif in all the putative MAP4Ks. This revealed that an S(G/A)WNF motif is absent in Clade I and Clade II of MAP4Ks, but present in 64% of Clade III MAP4Ks (Supplementary Table 3 and Supplementary Figure 3).

Evolutionary Insight in the MAP4K Family in the Green Lineage

Based on the tree topology of the inferred phylogenetic reconstruction of the MAP4K family, we divided the group into three major clades (I, II, and III). In addition, based on the protein motif structure and expansion in flowering plants, we defined four subclades in clade III (IIIA, IIIB, IIIC, and IIID) (Figure 1). While MAP4Ks are present in all the plant species we investigated, the absolute number of group members increased in land plants, especially in Clade III, and some species seemed to lack members for a particular clade (Figure 3). For example, O. tauri (Chlorophyta), Mesostigma viride (Streptophyta), Picea abies, and Solanum tuberosum lacked a clade II member, and clade I was absent in C. atmophyticus (Streptophyta), Salvinia cucullate (Fern), P. abies, S. tuberosum, and O. sativa (Figures 1, 3). Furthermore, Subclade IIIA (which contained angiosperm-exclusive MAP4Ks) appeared as a sister group of Subclade IIIB (which was exclusive to seed plant MAP4Ks) and orthologs for both clades were identified in ferns. Likely, the increase in the number of MAP4Ks in vascular plants (ferns, gymnosperms, and angiosperms) is associated with the increased rate of whole-genome duplications or polyploidization largely observed in land plants; thus, it is often invoked as one of the main causal agents for diversification and land colonization (Bowers et al., 2003; Qiao et al., 2019; Zhang et al., 2019). In addition, Subclade IIID (mainly defined by the presence of BLUS1) contained only members of gymnosperms and angiosperms (Figure 1), extending this beyond the angiosperm lineage (Harris et al., 2020). However, the single Gnetum montanum (TnS000821935t04) MAP4K sequence identified as belonging to this clade has retained only some of the described angiosperm motifs and lacks the LQQRMISLTEELQKEKLKNVQ motif in the disordered domain, which is likely essential for protein–protein interactions (Figure 1). Therefore, whether the gymnosperm MAP4K in subclade IIID retains a similar function as the other members of the clade remains unknown.

FIGURE 3

FIGURE 3

Summary of MAP4K family occurrence in the green lineage. Table based on the results of the inferred MAP4K phylogenetic tree within selected species. The presence (√) (and the number of MAP4Ks) or absence (x) of a likely functional ortholog are indicated. It should be noted that for those with missing clades, this can also be due to an incomplete genome.

Conclusion

In this study, we identified putative MAP4K members from representative plant species and inferred the phylogenetic relationships in the green lineage underpinning MAP4K evolution and diversification. As a whole, the presence and absence of typical sequence motifs in (putative) MAP4K sequences likely resulted in a functional diversification within the MAPK4 family. The next important and necessary step is thus to evaluate the functionality and activity of these MAP4Ks, not only in A. thaliana, but also in other members of the green lineage. Importantly, this also relates to some of the motifs that were identified and the large intrinsically disordered C-terminal half. Some of these analyses will be facilitated by the fewer MAP4K family members identified in this study in, for example, Marchantia, which is a valuable model for genetic studies (Ishizaki et al., 2016; Poveda, 2020).

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Statements

Data availability statement

The original contributions presented in the study are included in the article and Supplementary Material. Further inquiries can be directed to the corresponding author.

Author contributions

LP, CF, and LDV performed analyses. LP, CF, LDV, and ID interpreted the results and wrote the manuscript. ID coordinated the study. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the Research Foundation, Flanders (FWO.OPR.2019.0009.01). LP is a recipient of a Ph.D. grant from the Chinese Scholarship Council. LDV was supported by the post-doctoral fellowship of the Special Research Fund (Bijzonder Onderzoeksfonds, Ghent University).

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2021.650171/full#supplementary-material

References

  • 1

    AlbertV. A.BarbazukW. B.DepamphilisC. W.DerJ. P.Leebens-MackJ.MaH.et al (2013). The Amborella genome and the evolution of flowering plants.Science342:1241089. 10.1126/science.1241089

  • 2

    AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool.J. Mol. Biol.215403410. 10.1016/S0022-2836(05)80360-2

  • 3

    AltschulS. F.MaddenT. L.SchäfferA. A.ZhangJ.ZhangZ.MillerW.et al (1997). Gapped BLAST and PSI-BLAST: a new generation of protein database search programs.Nucleic Acids Res.2533893402. 10.1093/nar/25.17.3389

  • 4

    AppelsR.EversoleK.SteinN.FeuilletC.KellerB.RogersJ.et al (2018). Shifting the limits in wheat research and breeding using a fully annotated reference genome.Science361:eaar7191. 10.1126/science.aar7191

  • 5

    ArditoF.GiulianiM.PerroneD.TroianoG.Lo MuzioL. (2017). The crucial role of protein phosphorylation in cell signaling and its use as targeted therapy.Int. J. Mol. Med.40271280. 10.3892/ijmm.2017.3036

  • 6

    BaileyT. L.GribskovM. (1998). Combining evidence using p-values: application to sequence homology searches.Bioinformatics (Oxford, England)144854. 10.1093/bioinformatics/14.1.48

  • 7

    BanksJ. A.NishiyamaT.HasebeM.BowmanJ. L.GribskovM.DePamphilisC.et al (2011). The Selaginella genome identifies genetic changes associated with the evolution of vascular plants.Science332960963. 10.1126/science.1203810

  • 8

    BeckerB.MarinB. (2009). Streptophyte algae and the origin of embryophytes.Ann. Bot.1039991004. 10.1093/aob/mcp044

  • 9

    BeierS.HimmelbachA.ColmseeC.ZhangX.-Q.BarreroR. A.ZhangQ.et al (2017). Construction of a map-based reference genome sequence for barley, Hordeum vulgare L.Sci. Data4124. 10.1038/sdata.2017.44

  • 10

    BlancG.DuncanG.AgarkovaI.BorodovskyM.GurnonJ.KuoA.et al (2010). The Chlorella variabilis NC64A genome reveals adaptation to photosymbiosis, coevolution with viruses, and cryptic sex.Plant Cell.2229432955. 10.1105/tpc.110.076406

  • 11

    Blanc-MathieuR.VerhelstB.DerelleE.RombautsS.BougetF.-Y.CarréI.et al (2014). An improved genome of the model marine alga Ostreococcus tauri unfolds by assessing Illumina de novo assemblies.BMC Genomics15:1103. 10.1186/1471-2164-15-1103

  • 12

    BowersJ. E.ChapmanB. A.RongJ.PatersonA. H. (2003). Unravelling angiosperm genome evolution by phylogenetic analysis of chromosomal duplication events.Nature422433438. 10.1038/nature01521

  • 13

    BowmanJ. L.KohchiT.YamatoK. T.JenkinsJ.ShuS.IshizakiK.et al (2017). Insights into land plant evolution parnered from the Marchantia polymorpha genome.Cell171287.e5304.e5. 10.1016/j.cell.2017.09.030

  • 14

    CakirB.KılıçkayaO. (2015). Mitogen-activated protein kinase cascades in Vitis vinifera.Front. Plant Sci.6:556. 10.3389/fpls.2015.00556

  • 15

    CarreraA. C.AlexandrovK.RobertsT. M. (1993). The conserved lysine of the catalytic domain of protein kinases is actively involved in the phosphotransfer reaction and not required for anchoring ATP.Proc. Natl. Acad. Sci. U.S.A.90442446. 10.1073/pnas.90.2.442

  • 16

    ChampionA.PicaudA.HenryY. (2004). Reassessing the MAP3K and MAP4K relationships.Trends Plant Sci.9123129. 10.1016/j.tplants.2004.01.005

  • 17

    ChengS.XianW.FuY.MarinB.KellerJ.WuT.et al (2019). Genomes of subaerial Zygnematophyceae provide insights into land plant evolution.Cell1791057.e41067.e4. 10.1016/j.cell.2019.10.019

  • 18

    ChengZ.LiJ.-F.NiuY.ZhangX.-C.WoodyO. Z.XiongY.et al (2015). Pathogen-secreted proteases activate a novel plant immune pathway.Nature521213216. 10.1038/nature14243

  • 19

    ChmieleckiJ.PeiferM.JiaP.SocciN. D.HutchinsonK.VialeA.et al (2010). Targeted next-generation sequencing of DNA regions proximal to a conserved GXGXXG signaling motif enables systematic discovery of tyrosine kinase fusions in cancer.Nucleic Acids Res.3869856996. 10.1093/nar/gkq579

  • 20

    ChoS. K.LarueC. T.ChevalierD.WangH.JinnT.-L.ZhangS.et al (2008). Regulation of floral organ abscission in Arabidopsis thaliana.Proc. Natl. Acad. Sci. U.S.A.1051562915634. 10.1073/pnas.0805539105

  • 21

    ChuangH. C.WangX.TanT. H. (2016). MAP4K family kinases in immunity and inflammation.Adv. Immunol.129277314. 10.1016/bs.ai.2015.09.006

  • 22

    CrooksG. E.HonG.ChandoniaJ.-M.BrennerS. E. (2004). WebLogo: a sequence logo generator.Genome Res.1411881190. 10.1101/gr.849004

  • 23

    DanI.WatanabeN. M.KusumiA. (2001). The Ste20 group kinases as regulators of MAP kinase cascades.Trends Cell Biol.11220230. 10.1016/s0962-8924(01)01980-8

  • 24

    De ClerckO.KaoS. M.BogaertK. A.BlommeJ.FoflonkerF.KwantesM.et al (2018). Insights into the evolution of multicellularity from the Sea Lettuce genome.Curr. Biol. : CB282921.e52933.e5. 10.1016/j.cub.2018.08.015

  • 25

    EddyS. R. (1998). Profile hidden Markov models.Bioinformatics (Oxford, England)14755763. 10.1093/bioinformatics/14.9.755

  • 26

    El-GebaliS.MistryJ.BatemanA.EddyS. R.LucianiA.PotterS. C.et al (2019). The Pfam protein families database in 2019.Nucleic Acids Res.47D427D432. 10.1093/nar/gky995

  • 27

    FiedlerL. R.MaifoshieE.SchneiderM. D. (2014). Mouse models of heart failure: cell signaling and cell survival.Curr. Topics Dev. Biol.109171247. 10.1016/b978-0-12-397920-9.00002-0

  • 28

    Fürst-JansenJ. M.de VriesS.de VriesJ. (2020). Evo-physio: on stress responses and the earliest land plants.J. Exp. Bot.7132543269. 10.1093/jxb/eraa007

  • 29

    HanksS. K.HunterT. (1995). The eukaryotic protein kinase superfamily: kinase (catalytic) domain structure and classification 1.FASEB J.9576596. 10.1096/fasebj.9.8.7768349

  • 30

    HanksS. K.QuinnA. M. (1991). Protein kinase catalytic domain sequence database: identification of conserved features of primary structure and classification of family members.Methods Enzymol.2003862. 10.1016/0076-6879(91)00126-h

  • 31

    HanksS. K.QuinnA. M.HunterT. (1988). The protein kinase family: conserved features and deduced phylogeny of the catalytic domains.Science2414252. 10.1126/science.3291115

  • 32

    HarrisB. J.HarrisonC. J.HetheringtonA. M.WilliamsT. A. (2020). Phylogenomic evidence for the monophyly of Bryophytes and the reductive evolution of stomata.Curr. Biol. : CB302001.e22012.e2. 10.1016/j.cub.2020.03.048

  • 33

    HashiguchiA.KomatsuS. (2016). Impact of post-translational modifications of crop proteins under abiotic stress.Proteomes.4:42. 10.3390/proteomes4040042

  • 34

    HayesJ.Jr.UnderwoodR. M.PettisJ. (2008). A survey of honey bee colony losses in the US, fall 2007 to spring 2008.PLoS One3:e4071. 10.1371/journal.pone.0004071

  • 35

    HemmerW.McGloneM.TsigelnyI.TaylorS. S. (1997). Role of the glycine triad in the ATP-binding site of cAMP-dependent protein kinase.J. Biol. Chem.2721694616954. 10.1074/jbc.272.27.16946

  • 36

    HirschC. N.HirschC. D.BrohammerA. B.BowmanM. J.SoiferI.BaradO.et al (2016). Draft assembly of elite inbred line PH207 provides insights into genomic and transcriptome diversity in maize.Plant Cell2827002714. 10.1105/tpc.16.00353

  • 37

    HoriK.MaruyamaF.FujisawaT.TogashiT.YamamotoN.SeoM.et al (2014). Klebsormidium flaccidum genome reveals primary factors for plant terrestrial adaptation.Nat. Commun.5:3978. 10.1038/ncomms4978

  • 38

    HosotaniS.YamauchiS.KobayashiH.FujiS.KoyaS.ShimazakiK.-I.et al (2021). A BLUS1 kinase signal and a decrease in intercellular CO2 concentration are necessary for stomatal opening in response to blue light.Plant Cell3318131827. 10.1093/plcell/koab067

  • 39

    HruzT.LauleO.SzaboG.WessendorpF.BleulerS.OertleL.et al (2008). Genevestigator V3: a reference expression database for the meta-analysis of transcriptomes.Adv. Bioinform.2008:420747. 10.1155/2008/420747

  • 40

    HuT. T.PattynP.BakkerE. G.CaoJ.ChengJ.-F.ClarkR. M.et al (2011). The Arabidopsis lyrata genome sequence and the basis of rapid genome size change.Nat. Genet.43476481. 10.1038/ng.807

  • 41

    IshizakiK. (2017). Evolution of land plants: insights from molecular studies on basal lineages.Biosci. Biotechnol. Biochem.817380. 10.1080/09168451.2016.1224641

  • 42

    IshizakiK.NishihamaR.YamatoK. T.KohchiT. (2016). Molecular genetic tools and techniques for Marchantia polymorpha Research.Plant Cell Physiol.57262270. 10.1093/pcp/pcv097

  • 43

    JagodzikP.Tajdel-ZielinskaM.CieslaA.MarczakM.LudwikowA. (2018). Mitogen-activated protein kinase cascades in plant hormone signaling.Front. Plant Sci.9:1387. 10.3389/fpls.2018.01387

  • 44

    JaillonO.AuryJ.-M.NoelB.PolicritiA.ClepetC.CasagrandeA.et al (2007). The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla.Nature449463467. 10.1038/nature06148

  • 45

    JiangY.HanB.ZhangH.MariappanK. G.BigeardJ.ColcombetJ.et al (2019). MAP4K4 associates with BIK1 to regulate plant innate immunity.EMBO Rep.20:e47965. 10.15252/embr.201947965

  • 46

    JonakC.ÖkrészL.BögreL.HirtH. (2002). Complexity, cross talk and integration of plant MAP kinase signalling.Curr. Opin. Plant Biol.5415424. 10.1016/s1369-5266(02)00285-6

  • 47

    JonesP.BinnsD.ChangH.-Y.FraserM.LiW.McAnullaC.et al (2014). InterProScan 5: genome-scale protein function classification.Bioinformatics3012361240. 10.1093/bioinformatics/btu031

  • 48

    KannanN.NeuwaldA. F. (2005). Did protein kinase regulatory mechanisms evolve through elaboration of a simple structural component?J. Mol. Biol.351956972. 10.1016/j.jmb.2005.06.057

  • 49

    KorfI.FlicekP.DuanD.BrentM. R. (2001). Integrating genomic homology into gene structure prediction.Bioinformatics17S140S148. 10.1093/bioinformatics/17.suppl_1.S140

  • 50

    KrysanP. J.ColcombetJ. (2018). Cellular complexity in MAPK signaling in plants: questions and emerging tools to answer them.Front. Plant Sci.9:1674. 10.3389/fpls.2018.01674

  • 51

    LameschP.BerardiniT.LiD.SwarbreckD.WilksC.SasidharanR.et al (2012). A gene expression map of Arabidopsis thaliana development.Nucleic Acids Res.40D1202D1210. 10.1038/ng1543

  • 52

    LangD.UllrichK. K.MuratF.FuchsJ.JenkinsJ.HaasF. B.et al (2018). The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.Plant J.93515533. 10.1111/tpj.13801

  • 53

    LebererE.DignardD.HarcusD.ThomasD. Y.WhitewayM. (1992). The protein kinase homologue Ste20p is required to link the yeast pheromone response G-protein beta gamma subunits to downstream signalling components.EMBO J.1148154824. 10.1002/j.1460-2075.1992.tb05587.x

  • 54

    LeeD.RedfernO.OrengoC. (2007). Predicting protein function from sequence and structure.Nat. Rev. Mol. Cell Biol.89951005. 10.1038/nrm2281

  • 55

    LeeS.DevamaniT.SongH. D.SandhuM.LarsenA.SommeseR.et al (2017). Distinct structural mechanisms determine substrate affinity and kinase activity of protein kinase Cα.J. Biol. Chem.2921630016309. 10.1074/jbc.M117.804781

  • 56

    Leebens-MackJ. H.BarkerM. S.CarpenterE. J.DeyholosM. K.GitzendannerM. A.GrahamS. W.et al (2019). One thousand plant transcriptomes and the phylogenomics of green plants.Nature574679685. 10.1038/s41586-019-1693-2

  • 57

    LeliaertF.SmithD. R.MoreauH.HerronM. D.VerbruggenH.DelwicheC. F.et al (2012). Phylogeny and molecular evolution of the green algae.Crit. Rev. Plant Sci.31146. 10.1080/07352689.2011.615705

  • 58

    LeprinceA.-S.JouannicS.HamalA.KreisM.HenryY. (1999). Molecular characterisation of plant cDNAs BnMAP4Kα1 and BnMAP4Kα2 belonging to the GCK/SPS1 subfamily of MAP kinase kinase kinase kinase.Biochim. Biophys. Acta (BBA)-Gene Struct. Exp.1444113. 10.1016/S0167-4781(98)00246-2

  • 59

    LetunicI.BorkP. (2019). Interactive Tree Of Life (iTOL) v4: recent updates and new developments.Nucleic Acids Res.47W256W259. 10.1093/nar/gkz239

  • 60

    LiF.-W.BrouwerP.Carretero-PauletL.ChengS.De VriesJ.DelauxP.-M.et al (2018). Fern genomes elucidate land plant evolution and cyanobacterial symbioses.Nat. Plants4460472. 10.1038/s41477-018-0188-8

  • 61

    LiuS.LiuY.YangX.TongC.EdwardsD.ParkinI. A.et al (2014). The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes.Nat. Commun.5111. 10.1038/ncomms4930

  • 62

    LiuY.ZhouJ. (2018). MAPping kinase regulation of ICE1 in freezing tolerance.Trends Plant Sci.239193. 10.1016/j.tplants.2017.12.002

  • 63

    LlompartB.CastellsE.RíoA.RocaR.FerrandoA.StiefelV.et al (2003). The direct activation of MIK, a germinal center kinase (GCK)-like kinase, by MARK, a maize atypical receptor kinase, suggests a new mechanism for signaling through kinase-dead receptors.J. Biol. Chem.2784810548111. 10.1074/jbc.M307482200

  • 64

    MachidaN.UmikawaM.TakeiK.SakimaN.MyagmarB.-E.TairaK.et al (2004). Mitogen-activated protein kinase kinase kinase kinase 4 as a putative effector of Rap2 to activate the c-Jun N-terminal kinase.J. Biol. Chem.2791571115714. 10.1074/jbc.C300542200

  • 65

    MajorG.DaigleC.Stafford-RichardT.TebbjiF.LafleurE.CaronS.et al (2009). Characterization of ScMAP4K1, a MAP kinase kinase kinase kinase involved in ovule, seed and fruit development in Solanum chacoense Bitt.Plant Biol.102746.

  • 66

    McCormickR. F.TruongS. K.SreedasyamA.JenkinsJ.ShuS.SimsD.et al (2018). The Sorghum bicolor reference genome: improved assembly, gene annotations, a transcriptome atlas, and signatures of genome organization.Plant J.93338354. 10.1111/tpj.13781

  • 67

    MerchantS. S.ProchnikS. E.VallonO.HarrisE. H.KarpowiczS. J.WitmanG. B.et al (2007). The Chlamydomonas genome reveals the evolution of key animal and plant functions.Science (New York, NY)318245250. 10.1126/science.1143609

  • 68

    MészárosB.ErdõsG.DosztányiZ. (2018). IUPred2A: context-dependent prediction of protein disorder as a function of redox state and protein binding.Nucleic Acids Res.46W329W337. 10.1093/nar/gky384

  • 69

    ModiV.DunbrackR. L. (2019). Defining a new nomenclature for the structures of active and inactive kinases.Proc. Natl. Acad. Sci. U.S.A.11668186827. 10.1073/pnas.1814279116

  • 70

    Moreno-HagelsiebG.LatimerK. (2008). Choosing BLAST options for better detection of orthologs as reciprocal best hits.Bioinformatics24319324. 10.1093/bioinformatics/btm585

  • 71

    MorrisJ. L.PuttickM. N.ClarkJ. W.EdwardsD.KenrickP.PresselS.et al (2018). The timescale of early land plant evolution.Proc. Natl. Acad. Sci. U.S.A.115E2274E2283. 10.1073/pnas.1719588115

  • 72

    MotamayorJ. C.MockaitisK.SchmutzJ.HaiminenN.LivingstoneD.CornejoO.et al (2013). The genome sequence of the most widely cultivated cacao type and its use to identify candidate genes regulating pod color.Genome Biol.14125. 10.1186/gb-2013-14-6-r53

  • 73

    NeuwaldA. F.AltschulS. F. (2016). Bayesian top-down protein sequence alignment with inferred position-specific gap penalties.PLoS Comput. Biol.12:e1004936. 10.1371/journal.pcbi.1004936

  • 74

    NishiyamaT.SakayamaH.de VriesJ.BuschmannH.Saint-MarcouxD.UllrichK. K.et al (2018). The Chara genome: secondary complexity and implications for plant terrestrialization.Cell174448.e4464.e4. 10.1016/j.cell.2018.06.033

  • 75

    NystedtB.StreetN. R.WetterbomA.ZuccoloA.LinY.-C.ScofieldD. G.et al (2013). The Norway spruce genome sequence and conifer genome evolution.Nature497579584. 10.1038/nature12211

  • 76

    OuyangS.ZhuW.HamiltonJ.LinH.CampbellM.ChildsK.et al (2007). The TIGR rice genome annotation resource: improvements and new features.Nucleic Acids Res.35D883D887. 10.1093/nar/gkl976

  • 77

    PalenikB.GrimwoodJ.AertsA.RouzéP.SalamovA.PutnamN.et al (2007). The tiny eukaryote Ostreococcus provides genomic insights into the paradox of plankton speciation.Proc. Natl. Acad. Sci. U.S.A.104:7705. 10.1073/pnas.0611046104

  • 78

    PanL.De SmetI. (2020). Expanding the mitogen-activated protein kinase (MAPK) universe: an update on MAP4Ks.Front. Plant Sci.11:1220. 10.3389/fpls.2020.01220

  • 79

    PaulA.SrinivasanN. (2020). Genome-wide and structural analyses of pseudokinases encoded in the genome of Arabidopsis thaliana provide functional insights.Proteins: Struct. Funct. Bioinform.8816201638. 10.1002/prot.25981

  • 80

    PiresN. D.DolanL. (2012). Morphological evolution in land plants: new designs with old genes.Philosop. Trans. R. Soc. B: Biol.Sci.367508518. 10.1098/rstb.2011.0252

  • 81

    PoinatP.De ArcangelisA.SookhareeaS.ZhuX.HedgecockE. M.LabouesseM.et al (2002). A conserved interaction between β1 integrin/PAT-3 and Nck-interacting kinase/MIG-15 that mediates commissural axon navigation in C. elegans.Curr. Biol.12622631. 10.1016/s0960-9822(02)00764-9

  • 82

    PovedaJ. (2020). Marchantia polymorpha as a model plant in the evolutionary study of plant-microorganism interactions.Curr. Plant Biol.23:100152. 10.1016/j.cpb.2020.100152

  • 83

    ProchnikS. E.UmenJ.NedelcuA. M.HallmannA.MillerS. M.NishiiI.et al (2010). Genomic analysis of organismal complexity in the multicellular green alga Volvox carteri.Science329:223. 10.1126/science.1188800

  • 84

    QiaoX.LiQ.YinH.QiK.LiL.WangR.et al (2019). Gene duplication and evolution in recurring polyploidization–diploidization cycles in plants.Genome Biol.20:38. 10.1186/s13059-019-1650-2

  • 85

    RoskoskiR.Jr. (2004). Src protein–tyrosine kinase structure and regulation.Biochem. Biophys. Res. Commun.32411551164. 10.1016/j.bbrc.2004.09.171

  • 86

    SatoS.TabataS.HirakawaH.AsamizuE.ShirasawaK.IsobeS.et al (2012). The tomato genome sequence provides insights into fleshy fruit evolution.Nature485635641. 10.1038/nature11119

  • 87

    SchmutzJ.CannonS. B.SchlueterJ.MaJ.MitrosT.NelsonW.et al (2010). Genome sequence of the palaeopolyploid soybean.Nature463178183. 10.1038/nature08670

  • 88

    SchnabelJ.HombachP.WaksmanT.GiurianiG.PetersenJ.ChristieJ. M. (2018). A chemical genetic approach to engineer phototropin kinases for substrate labeling.J. Biol. Chem.29356135623. 10.1074/jbc.RA118.001834

  • 89

    SellsM. A.ChernoffJ. (1997). Emerging from the Pak: the p21-activated protein kinase family.Trends Cell Biol.7162167. 10.1016/S0962-8924(97)01003-9

  • 90

    SeoG.HanH.VargasR. E.YangB.LiX.WangW. (2020). MAP4K interactome reveals STRN4 as a key STRIPAK complex component in Hippo pathway regulation.Cell Rep.32:107860. 10.1016/j.celrep.2020.107860

  • 91

    SharmaS. K.BolserD.de BoerJ.SønderkærM.AmorosW.CarboniM. F.et al (2013). Construction of reference chromosome-scale pseudomolecules for potato: integrating the potato genome with genetic and physical maps.G3 Genes Genomes Genet. 320312047. 10.1534/g3.113.007153

  • 92

    SteinbergS. F. (2018). Post-translational modifications at the ATP-positioning G-loop that regulate protein kinase activity.Pharmacol. Res.135181187. 10.1016/j.phrs.2018.07.009

  • 93

    TairaK.UmikawaM.TakeiK.MyagmarB.-E.ShinzatoM.MachidaN.et al (2004). The Traf2-and Nck-interacting kinase as a putative effector of Rap2 to regulate actin cytoskeleton.J. Biol. Chem.2794948849496. 10.1074/jbc.M406370200

  • 94

    TakemiyaA.SugiyamaN.FujimotoH.TsutsumiT.YamauchiS.HiyamaA.et al (2013). Phosphorylation of BLUS1 kinase by phototropins is a primary step in stomatal opening.Nat. Commun.418. 10.1038/ncomms3094

  • 95

    TorkamaniA.KannanN.TaylorS. S.SchorkN. J. (2008). Congenital disease SNPs target lineage specific structural elements in protein kinases.Proc. Natl. Acad. Sci. U.S.A.10590119016. 10.1073/pnas.0802403105

  • 96

    TuskanG. A.DiFazioS.JanssonS.BohlmannJ.GrigorievI.HellstenU.et al (2006). The genome of black cottonwood, Populus Trichocarpa (Torr. &Gray).Science313:1596. 10.1126/science.1128691

  • 97

    VuL. D.GevaertK.De SmetI. (2018). Protein language: post-translational modifications talking to each other.Trends Plant Sci.2310681080. 10.1016/j.tplants.2018.09.004

  • 98

    VuL. D.XuX.ZhuT.PanL.Van ZantenM.De JongD.et al (2021). The membrane-localized protein kinase MAP4K4/TOT3 regulates thermomorphogenesis.Nat. Commun.12114. 10.1038/s41467-021-23112-0

  • 99

    WanT.LiuZ.-M.LiL.-F.LeitchA. R.LeitchI. J.LohausR.et al (2018). A genome for gnetophytes and early evolution of seed plants.Nat. Plants48289. 10.1038/s41477-017-0097-2

  • 100

    WangS.LiL.LiH.SahuS. K.WangH.XuY.et al (2020). Genomes of early-diverging streptophyte algae shed light on plant terrestrialization.Nat. Plants695106. 10.1038/s41477-019-0560-3

  • 101

    WinterD.VinegarB.NahalH.AmmarR.WilsonG. V.ProvartN. J. (2007). An “Electronic Fluorescent Pictograph” browser for exploring and analyzing large-scale biological data sets.PLoS One2:e718. 10.1371/journal.pone.0000718

  • 102

    WordenA. Z.LeeJ.-H.MockT.RouzéP.SimmonsM. P.AertsA. L.et al (2009). Green evolution and dynamic adaptations revealed by genomes of the marine picoeukaryotes Micromonas.Science324:268. 10.1126/science.1167222

  • 103

    XiongJ.CuiX.YuanX.YuX.SunJ.GongQ. (2016). The Hippo/STE20 homolog SIK1 interacts with MOB1 to regulate cell proliferation and cell expansion in Arabidopsis.J. Exp. Bot.6714611475. 10.1093/jxb/erv538

  • 104

    XuB.-E.EnglishJ. M.WilsbacherJ. L.StippecS.GoldsmithE. J.CobbM. H. (2000). WNK1, a novel mammalian serine/threonine protein kinase lacking the catalytic lysine in subdomain II.J. Biol. Chem.2751679516801. 10.1074/jbc.275.22.16795

  • 105

    XuJ.ZhangS. (2015). Mitogen-activated protein kinase cascades in signaling plant growth and development.Trends Plant Sci.205664. 10.1016/j.tplants.2014.10.001

  • 106

    YanH.ZhaoY.ShiH.LiJ.WangY.TangD. (2018). BRASSINOSTEROID-SIGNALING KINASE1 phosphorylates MAPKKK5 to regulate immunity in Arabidopsis.Plant Physiol.17629913002. 10.1104/pp.17.01757

  • 107

    YoungN. D.DebelléF.OldroydG. E. D.GeurtsR.CannonS. B.UdvardiM. K.et al (2011). The Medicago genome provides insight into the evolution of rhizobial symbioses.Nature480520524. 10.1038/nature10625

  • 108

    ZhangJ.FuX.-X.LiR.-Q.ZhaoX.LiuY.LiM.-H.et al (2020). The hornwort genome and early land plant evolution.Nat. Plants6107118. 10.1038/s41477-019-0588-4

  • 109

    ZhangK.WangX.ChengF. (2019). Plant polyploidy: origin, evolution, and Its influence on crop domestication.Horticult. Plant J.5231239. 10.1016/j.hpj.2019.11.003

  • 110

    ZhangL.WangJ.-C.HouL.CaoP.-R.WuL.ZhangQ.-S.et al (2015). Functional role of histidine in the conserved His-x-Asp motif in the catalytic core of protein kinases.Sci. Rep.5:10115. 10.1016/j.hpj.2019.11.003

  • 111

    ZhangM.ChiangY.-H.ToruñoT. Y.LeeD.MaM.LiangX.et al (2018). The MAP4 kinase SIK1 ensures robust extracellular ROS burst and antibacterial immunity in plants.Cell Host Microbe.24379.e5391.e5. 10.1016/j.chom.2018.08.007

  • 112

    ZhuQ.ShaoY.GeS.ZhangM.ZhangT.HuX.et al (2019). A MAPK cascade downstream of IDA–HAE/HSL2 ligand–receptor pair in lateral root emergence.Nat. Plants5414423. 10.1038/s41477-019-0396-x

  • 113

    ZimmermannP.HennigL.GruissemW. (2005). Gene-expression analysis and network discovery using Genevestigator.Trends Plant Sci.10407409. 10.1016/j.tplants.2005.07.003

  • 114

    ZimmermannP.Hirsch-HoffmannM.HennigL.GruissemW. (2004). GENEVESTIGATOR. Arabidopsis microarray database and analysis toolbox.Plant Physiol.13626212632. 10.1104/pp.104.046367

Summary

Keywords

signaling, MAP4K, phylogenetic analysis, evolution, motifs

Citation

Pan L, Fonseca De Lima CF, Vu LD and De Smet I (2021) A Comprehensive Phylogenetic Analysis of the MAP4K Family in the Green Lineage. Front. Plant Sci. 12:650171. doi: 10.3389/fpls.2021.650171

Received

27 January 2021

Accepted

22 July 2021

Published

13 August 2021

Volume

12 - 2021

Edited by

Angelica Cibrian-Jaramillo, Centro de Investigaciones y Estudios Avanzados, Instituto Politécnico Nacional de México (CINVESTAV), Mexico

Reviewed by

Luis Delaye, Centro de Investigaciones y Estudios Avanzados, Instituto Politécnico Nacional de México (CINVESTAV), Mexico; Jan de Vries, University of Göttingen, Germany; Adrian Reyes-Prieto, University of New Brunswick Fredericton, Canada

Updates

Copyright

*Correspondence: Ive De Smet,

These authors have contributed equally to this work

This article was submitted to Plant Systematics and Evolution, a section of the journal Frontiers in Plant Science

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics