Soybean kinome: functional classification and gene expression patterns

Genome-wide characterization of the sunflower kinome: classification, evolutionary analysis and expression patterns under different stresses

Protein kinases play a significant role in plant responses to biotic and abiotic stresses, as well as in growth and development. While the kinome has been extensively investigated in crops such as Arabidopsis thaliana, soybean, common bean, and cotton, studies on protein kinases in sunflower remain limited. Our objective is to explore protein kinases in sunflower to bridge the research gap and enhance the understanding of their functions. We identified a total of 2,583 protein kinases from sunflower, which were classified into 22 families and 121 subfamilies. By comparing the subfamily members between sunflower and other species, we found that three subfamilies in sunflower-RLK-Pelle_CrRLK1L-1, RLK-Pelle_SD-2b, and RLK-Pelle_WAK-had undergone significant expansion. We then investigated the chromosomal distribution, molecular weight, isoelectric point, transmembrane domain, signal peptide, and structural and evolutionary diversity of the protein kinases. Through these studies, we have obtained a basic understanding of protein kinases in sunflower. To investigate the role of protein kinases in sunflower's response to biotic and abiotic stresses, we obtained 534 transcriptome datasets from various research groups, covering eight types of abiotic stress and two types of biotic stress. For the first time, we overcame the batch effects in the data and utilized a gene scoring system developed by our lab to perform a comprehensive analysis of multiple transcriptome datasets from different research groups. Ultimately, 73 key protein kinases were identified from numerous candidates, and functional annotation revealed that they are key members of signaling pathways such as ABA, MAPK, and SOS, actively participating in sunflower's response to biotic and abiotic stresses. In summary, through the exploration of protein kinases in sunflower, we have filled the gap in protein kinase research and provided a substantial amount of foundational data. By using the new scoring method to eliminate batch effects between transcriptome datasets, we achieved the first comprehensive analysis of large-scale transcriptome data. This method allows for a more thorough and detailed identification of key protein kinases that are widely regulated under various stress conditions, providing numerous candidate genes for sunflower stress resistance research.

The soybean immune receptor GmBIR1 regulates host transcriptome, spliceome, and immunity during cyst nematode infection

BAK1-INTERACTING RECEPTOR LIKE KINASE1 (BIR1) is a negative regulator of various aspects of disease resistance and immune responses. Here, we investigated the functional role of soybean (Glycine max) BIR1 (GmBIR1) during soybean interaction with soybean cyst nematode (SCN, Heterodera glycines) and the molecular mechanism through which GmBIR1 regulates plant immunity. Overexpression of wild-type variant of GmBIR1 (WT-GmBIR1) using transgenic soybean hairy roots significantly increased soybean susceptibility to SCN, whereas overexpression of kinase-dead variant (KD-GmBIR1) significantly increased plant resistance. Transcriptome analysis revealed that genes oppositely regulated in WT-GmBIR1 and KD-GmBIR1 upon SCN infection were enriched primarily in defense and immunity-related functions. Quantitative phosphoproteomic analysis identified 208 proteins as putative substrates of the GmBIR1 signaling pathway, 114 of which were differentially phosphorylated upon SCN infection. In addition, the phosphoproteomic data pointed to a role of the GmBIR1 signaling pathway in regulating alternative pre-mRNA splicing. Genome-wide analysis of splicing events provided compelling evidence supporting a role of the GmBIR1 signaling pathway in establishing alternative splicing during SCN infection. Our results provide novel mechanistic insights into the function of the GmBIR1 signaling pathway in regulating soybean transcriptome and spliceome via differential phosphorylation of splicing factors and regulation of splicing events of pre-mRNA decay- and spliceosome-related genes.

Functional annotation of proteins for signaling network inference in non-model species

Molecular biology aims to understand cellular responses and regulatory dynamics in complex biological systems. However, these studies remain challenging in non-model species due to poor functional annotation of regulatory proteins. To overcome this limitation, we develop a multi-layer neural network that determines protein functionality directly from the protein sequence. We annotate kinases and phosphatases in Glycine max. We use the functional annotations from our neural network, Bayesian inference principles, and high resolution phosphoproteomics to infer phosphorylation signaling cascades in soybean exposed to cold, and identify Glyma.10G173000 (TOI5) and Glyma.19G007300 (TOT3) as key temperature regulators. Importantly, the signaling cascade inference does not rely upon known kinase motifs or interaction data, enabling de novo identification of kinase-substrate interactions. Conclusively, our neural network shows generalization and scalability, as such we extend our predictions to Oryza sativa, Zea mays, Sorghum bicolor, and Triticum aestivum. Taken together, we develop a signaling inference approach for non-model species leveraging our predicted kinases and phosphatases.

Soybean gene co-expression network analysis identifies two co-regulated gene modules associated with nodule formation and development

Gene co-expression network analysis is an efficient systems biology approach for the discovery of novel gene functions and trait-associated gene modules. To identify clusters of functionally related genes involved in soybean nodule formation and development, we performed a weighted gene co-expression network analysis. Two nodule-specific modules (NSM-1 and NSM-2, containing 304 and 203 genes, respectively) were identified. The NSM-1 gene promoters were significantly enriched in cis-binding elements for ERF, MYB, and C2H2-type zinc transcription factors, whereas NSM-2 gene promoters were enriched in cis-binding elements for TCP, bZIP, and bHLH transcription factors, suggesting a role of these regulatory factors in the transcriptional activation of nodule co-expressed genes. The co-expressed gene modules included genes with potential novel roles in nodulation, including those involved in xylem development, transmembrane transport, the ethylene signalling pathway, cytoskeleton organization, cytokinesis and regulation of the cell cycle, regulation of meristem initiation and growth, transcriptional regulation, DNA methylation, and histone modifications. Functional analysis of two co-expressed genes using TILLING mutants provided novel insight into the involvement of unsaturated fatty acid biosynthesis and folate metabolism in nodule formation and development. The identified gene co-expression modules provide valuable resources for further functional genomics studies to dissect the genetic basis of nodule formation and development in soybean.

Genome-wide characterization of the common bean kinome: Catalog and insights into expression patterns and genetic organization

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and is the core regulator of cellular signaling. Even considering this substantial importance, the kinome of common bean (Phaseolus vulgaris) has not been profiled yet. Here, we identified and characterised the complete set of kinases of common bean, performing an in-depth investigation with phylogenetic analyses and measurements of gene distribution, structural organization, protein properties, and expression patterns over a large set of RNA-Sequencing data. Being composed of 1,203 PKs distributed across all P. vulgaris chromosomes, this set represents 3.25% of all predicted proteins for the species. These PKs could be classified into 20 groups and 119 subfamilies, with a more pronounced abundance of subfamilies belonging to the receptor-like kinase (RLK)-Pelle group. In addition to provide a vast and rich reservoir of data, our study supplied insights into the compositional similarities between PK subfamilies, their evolutionary divergences, highly variable functional profile, structural diversity, and expression patterns, modeled with coexpression networks for investigating putative interactions associated with stress response.

The rubber tree kinome: Genome-wide characterization and insights into coexpression patterns associated with abiotic stress responses

The protein kinase (PK) superfamily constitutes one of the largest and most conserved protein families in eukaryotic genomes, comprising core components of signaling pathways in cell regulation. Despite its remarkable relevance, only a few kinase families have been studied in Hevea brasiliensis. A comprehensive characterization and global expression analysis of the PK superfamily, however, is currently lacking. In this study, with the aim of providing novel inferences about the mechanisms associated with the stress response developed by PKs and retained throughout evolution, we identified and characterized the entire set of PKs, also known as the kinome, present in the Hevea genome. Different RNA-sequencing datasets were employed to identify tissue-specific expression patterns and potential correspondences between different rubber tree genotypes. In addition, coexpression networks under several abiotic stress conditions, such as cold, drought and latex overexploitation, were employed to elucidate associations between families and tissues/stresses. A total of 1,809 PK genes were identified using the current reference genome assembly at the scaffold level, and 1,379 PK genes were identified using the latest chromosome-level assembly and combined into a single set of 2,842 PKs. These proteins were further classified into 20 different groups and 122 families, exhibiting high compositional similarities among family members and with two phylogenetically close species Manihot esculenta and Ricinus communis. Through the joint investigation of tandemly duplicated kinases, transposable elements, gene expression patterns, and coexpression events, we provided insights into the understanding of the cell regulation mechanisms in response to several conditions, which can often lead to a significant reduction in rubber yield.

Genome-Wide Association Studies in Sunflower: Towards Sclerotinia sclerotiorum and Diaporthe/Phomopsis Resistance Breeding

Diseases caused by necrotrophic fungi, such as the cosmopolitan Sclerotinia sclerotiorum and the Diaporthe/Phomopsis complex, are among the most destructive diseases of sunflower worldwide. The lack of complete resistance combined with the inefficiency of chemical control makes assisted breeding the best strategy for disease control. In this work, we present an integrated genome-wide association (GWA) study investigating the response of a diverse panel of sunflower inbred lines to both pathogens. Phenotypic data for Sclerotinia head rot (SHR) consisted of five disease descriptors (disease incidence, DI; disease severity, DS; area under the disease progress curve for DI, AUDPCI, and DS, AUDPCS; and incubation period, IP). Two disease descriptors (DI and DS) were evaluated for two manifestations of Diaporthe/Phomopsis: Phomopsis stem canker (PSC) and Phomopsis head rot (PHR). In addition, a principal component (PC) analysis was used to derive transformed phenotypes as inputs to a univariate GWA (PC-GWA). Genotypic data comprised a panel of 4269 single nucleotide polymorphisms (SNP), generated via genotyping-by-sequencing. The GWA analysis revealed 24 unique marker-trait associations for SHR, 19 unique marker-trait associations for Diaporthe/Phomopsis diseases, and 7 markers associated with PC1 and PC2. No common markers were found for the response to the two pathogens. Nevertheless, epistatic interactions were identified between markers significantly associated with the response to S. sclerotiorum and Diaporthe/Phomopsis. This suggests that, while the main determinants of resistance may differ for the two pathogens, there could be an underlying common genetic basis. The exploration of regions physically close to the associated markers yielded 364 genes, of which 19 were predicted as putative disease resistance genes. This work presents the first simultaneous evaluation of two manifestations of Diaporthe/Phomopsis in sunflower, and undertakes a comprehensive GWA study by integrating PSC, PHR, and SHR data. The multiple regions identified, and their exploration to identify candidate genes, contribute not only to the understanding of the genetic basis of resistance, but also to the development of tools for assisted breeding.

Insights into the regulation of wild soybean tolerance to salt-alkaline stress

Soybean is an important grain and oil crop. In China, there is a great contradiction between soybean supply and demand. China has around 100 million ha of salt-alkaline soil, and at least 10 million could be potentially developed for cultivated land. Therefore, it is an effective way to improve soybean production by breeding salt-alkaline-tolerant soybean cultivars. Compared with wild soybean, cultivated soybean has lost a large number of important genes related to environmental adaptation during the long-term domestication and improvement process. Therefore, it is greatly important to identify the salt-alkaline tolerant genes in wild soybean, and investigate the molecular basis of wild soybean tolerance to salt-alkaline stress. In this review, we summarized the current research regarding the salt-alkaline stress response in wild soybean. The genes involved in the ion balance and ROS scavenging in wild soybean were summarized. Meanwhile, we also introduce key protein kinases and transcription factors that were reported to mediate the salt-alkaline stress response in wild soybean. The findings summarized here will facilitate the molecular breeding of salt-alkaline tolerant soybean cultivars.

Silveira S, Brommonschenkel SH, Fontes EPB. RLPredictiOme, a Machine Learning-Derived Method for High-Throughput Prediction of Plant Receptor-like Proteins, Reveals Novel Classes of Transmembrane Receptors

Cell surface receptors play essential roles in perceiving and processing external and internal signals at the cell surface of plants and animals. The receptor-like protein kinases (RLK) and receptor-like proteins (RLPs), two major classes of proteins with membrane receptor configuration, play a crucial role in plant development and disease defense. Although RLPs and RLKs share a similar single-pass transmembrane configuration, RLPs harbor short divergent C-terminal regions instead of the conserved kinase domain of RLKs. This RLP receptor structural design precludes sequence comparison algorithms from being used for high-throughput predictions of the RLP family in plant genomes, as has been extensively performed for RLK superfamily predictions. Here, we developed the RLPredictiOme, implemented with machine learning models in combination with Bayesian inference, capable of predicting RLP subfamilies in plant genomes. The ML models were simultaneously trained using six types of features, along with three stages to distinguish RLPs from non-RLPs (NRLPs), RLPs from RLKs, and classify new subfamilies of RLPs in plants. The ML models achieved high accuracy, precision, sensitivity, and specificity for predicting RLPs with relatively high probability ranging from 0.79 to 0.99. The prediction of the method was assessed with three datasets, two of which contained leucine-rich repeats (LRR)-RLPs from Arabidopsis and rice, and the last one consisted of the complete set of previously described Arabidopsis RLPs. In these validation tests, more than 90% of known RLPs were correctly predicted via RLPredictiOme. In addition to predicting previously characterized RLPs, RLPredictiOme uncovered new RLP subfamilies in the Arabidopsis genome. These include probable lipid transfer (PLT)-RLP, plastocyanin-like-RLP, ring finger-RLP, glycosyl-hydrolase-RLP, and glycerophosphoryldiester phosphodiesterase (GDPD, GDPDL)-RLP subfamilies, yet to be characterized. Compared to the only Arabidopsis GDPDL-RLK, molecular evolution studies confirmed that the ectodomain of GDPDL-RLPs might have undergone a purifying selection with a predominance of synonymous substitutions. Expression analyses revealed that predicted GDPGL-RLPs display a basal expression level and respond to developmental and biotic signals. The results of these biological assays indicate that these subfamily members have maintained functional domains during evolution and may play relevant roles in development and plant defense. Therefore, RLPredictiOme provides a framework for genome-wide surveys of the RLP superfamily as a foundation to rationalize functional studies of surface receptors and their relationships with different biological processes.

Genome wide study of cysteine rich receptor like proteins in Gossypium sp

Cysteine-rich receptor-like-kinases (CRKs), a transmembrane subfamily of receptor-like kinase, play crucial roles in plant adaptation. As such cotton is the major source of fiber for the textile industry, but environmental stresses are limiting its growth and production. Here, we have performed a deep computational analysis of CRKs in five Gossypium species, including G. arboreum (60 genes), G. raimondii (74 genes), G. herbaceum (65 genes), G. hirsutum (118 genes), and G. barbadense (120 genes). All identified CRKs were classified into 11 major classes and 43 subclasses with the finding of several novel CRK-associated domains including ALMT, FUSC_2, Cript, FYVE, and Pkinase. Of these, DUF26_DUF26_Pkinase_Tyr was common and had elevated expression under different biotic and abiotic stresses. Moreover, the 35 land plants comparison identified several new CRKs domain-architectures. Likewise, several SNPs and InDels were observed in CLCuD resistant G. hirsutum. The miRNA target side prediction and their expression profiling in different tissues predicted miR172 as a major CRK regulating miR. The expression profiling of CRKs identified multiple clusters with co-expression under certain stress conditions. The expression analysis under CLCuD highlighted the role of GhCRK057, GhCRK059, GhCRK058, and GhCRK081 in resistant accession. Overall, these results provided primary data for future potential functional analysis as well as a reference study for other agronomically important crops.

Kinase-dead mutation: A novel strategy for improving soybean resistance to soybean cyst nematode Heterodera glycines

Protein kinases phosphorylate proteins for functional changes and are involved in nearly all cellular processes, thereby regulating almost all aspects of plant growth and development, and responses to biotic and abiotic stresses. We generated two independent co-expression networks of soybean genes using control and stress response gene expression data and identified 392 differentially highly interconnected kinase hub genes among the two networks. Of these 392 kinases, 90 genes were identified as "syncytium highly connected hubs", potentially essential for activating kinase signalling pathways in the nematode feeding site. Overexpression of wild-type coding sequences of five syncytium highly connected kinase hub genes using transgenic soybean hairy roots enhanced plant susceptibility to soybean cyst nematode (SCN; Heterodera glycines) Hg Type 0 (race 3). In contrast, overexpression of kinase-dead variants of these five syncytium kinase hub genes significantly enhanced soybean resistance to SCN. Additionally, three of the five tested kinase hub genes enhanced soybean resistance to SCN Hg Type 1.2.5.7 (race 2), highlighting the potential of the kinase-dead approach to generate effective and durable resistance against a wide range of SCN Hg types. Subcellular localization analysis revealed that kinase-dead mutations do not alter protein cellular localization, confirming the structure-function of the kinase-inactive variants in producing loss-of-function phenotypes causing significant decrease in nematode susceptibility. Because many protein kinases are highly conserved and are involved in plant responses to various biotic and abiotic stresses, our approach of identifying kinase hub genes and their inactivation using kinase-dead mutation could be translated for biotic and abiotic stress tolerance.

Genome-wide analysis reveals widespread roles for RcREM genes in floral organ development in Rosa chinensis

Members of the REM (Reproductive Meristem) gene family are expressed primarily in reproductive meristems and floral organs. However, their evolution and their functional profiles in flower development remain poorly understood. Here, we performed genome-wide identification and evolutionary analysis of the REM gene family in Rosaceae. This family has been greatly expanded in rose (Rosa chinensis) compared to other species, primarily through tandem duplication. Expression analysis revealed that most RcREM genes are specifically expressed in reproductive organs and that their specific expression patterns are dramatically altered in rose plants with mutations affecting floral organs. Protein-protein interaction analysis indicated that RcREM14 interact with RcAP1 (one of the homology of A class genes in ABCDE model), highlighting the roles of RcREM genes in floral organ identity. Finally, co-expression network analysis indicated that RcREM genes are co-expressed with a high proportion of key genes that regulate flowering time, floral organ development, and cell proliferation and expansion in R. chinensis.

Pecan kinome: classification and expression analysis of all protein kinases in Carya illinoinensis

Protein kinases (PKs) are involved in plant growth and stress responses, and constitute one of the largest superfamilies due to numerous gene duplications. However, limited PKs have been functionally described in pecan, an economically important nut tree. Here, the comprehensive identification, annotation and classification of the entire pecan kinome are reported. A total of 967 PK genes were identified from the pecan genome, and further classified into 20 different groups and 121 subfamilies using the kinase domain sequences, which were verified by phylogenetic analysis. The receptor-like kinase (RLK) group contained 565 members, which constituted the largest group. Gene duplication contributed to the expansion of pecan kinome, 169 segmental duplication events including 285 PK genes were found, and the Ka/Ks ratio revealed they experienced strong negative selection. The RNA-Seq data of PK genes in pecan were further analyzed at the subfamily level, and different PK subfamilies performed various expression patterns across pecan embryo development or drought treatment, suggesting PK genes in pecan are involved in embryo development and drought stress response. Taken together, this study provides insight into the classification, expansion, evolution, and expression of pecan PKs. Our findings regarding expansion, expression and co-expression analyses lay a good foundation for future research to understand the roles of pecan PKs, and more efficiently determine the key candidate genes.

The Cowpea Kinome: Genomic and Transcriptomic Analysis Under Biotic and Abiotic Stresses

The present work represents a pioneering effort, being the first to analyze genomic and transcriptomic data from Vigna unguiculata (cowpea) kinases. We evaluated the cowpea kinome considering its genome-wide distribution and structural characteristics (at the gene and protein levels), sequence evolution, conservation among Viridiplantae species, and gene expression in three cowpea genotypes under different stress situations, including biotic (injury followed by virus inoculation-CABMV or CPSMV) and abiotic (root dehydration). The structural features of cowpea kinases (VuPKs) indicated that 1,293 bona fide VuPKs covered 20 groups and 118 different families. The RLK-Pelle was the largest group, with 908 members. Insights on the mechanisms of VuPK genomic expansion and conservation among Viridiplantae species indicated dispersed and tandem duplications as major forces for VuPKs' distribution pattern and high orthology indexes and synteny with other legume species, respectively. K _a /K _s ratios showed that almost all (91%) of the tandem duplication events were under purifying selection. Candidate cis-regulatory elements were associated with different transcription factors (TFs) in the promoter regions of the RLK-Pelle group. C2H2 TFs were closely associated with the promoter regions of almost all scrutinized families for the mentioned group. At the transcriptional level, it was suggested that VuPK up-regulation was stress, genotype, or tissue dependent (or a combination of them). The most prominent families in responding (up-regulation) to all the analyzed stresses were RLK-Pelle_DLSV and CAMK_CAMKL-CHK1. Concerning root dehydration, it was suggested that the up-regulated VuPKs are associated with ABA hormone signaling, auxin hormone transport, and potassium ion metabolism. Additionally, up-regulated VuPKs under root dehydration potentially assist in a critical physiological strategy of the studied cowpea genotype in this assay, with activation of defense mechanisms against biotic stress while responding to root dehydration. This study provides the foundation for further studies on the evolution and molecular function of VuPKs.

TILLING-by-Sequencing+ Reveals the Role of Novel Fatty Acid Desaturases (GmFAD2-2s) in Increasing Soybean Seed Oleic Acid Content

Soybean is the second largest source of oil worldwide. Developing soybean varieties with high levels of oleic acid is a primary goal of the soybean breeders and industry. Edible oils containing high level of oleic acid and low level of linoleic acid are considered with higher oxidative stability and can be used as a natural antioxidant in food stability. All developed high oleic acid soybeans carry two alleles; GmFAD2-1A and GmFAD2-1B. However, when planted in cold soil, a possible reduction in seed germination was reported when high seed oleic acid derived from GmFAD2-1 alleles were used. Besides the soybean fatty acid desaturase (GmFAD2-1) subfamily, the GmFAD2-2 subfamily is composed of five members, including GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E. Segmental duplication of GmFAD2-1A/GmFAD2-1B, GmFAD2-2A/GmFAD2-2C, GmFAD2-2A/GmFAD2-2D, and GmFAD2-2D/GmFAD2-2C have occurred about 10.65, 27.04, 100.81, and 106.55 Mya, respectively. Using TILLING-by-Sequencing+ technology, we successfully identified 12, 8, 10, 9, and 19 EMS mutants at the GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E genes, respectively. Functional analyses of newly identified mutants revealed unprecedented role of the five GmFAD2-2A, GmFAD2-2B, GmFAD2-2C, GmFAD2-2D, and GmFAD2-2E members in controlling the seed oleic acid content. Most importantly, unlike GmFAD2-1 members, subcellular localization revealed that members of the GmFAD2-2 subfamily showed a cytoplasmic localization, which may suggest the presence of an alternative fatty acid desaturase pathway in soybean for converting oleic acid content without substantially altering the traditional plastidial/ER fatty acid production.

Orthology and synteny analysis of receptor-like kinases "RLK" and receptor-like proteins "RLP" in legumes

BACKGROUND

Legume species are an important plant model because of their protein-rich physiology. The adaptability and productivity of legumes are limited by major biotic and abiotic stresses. Responses to these stresses directly involve plasma membrane receptor proteins known as receptor-like kinases and receptor-like proteins. Evaluating the homology relations among RLK and RLP for seven legume species, and exploring their presence among synteny blocks allow an increased understanding of evolutionary relations, physical position, and chromosomal distribution in related species and their shared roles in stress responses.

RESULTS

Typically, a high proportion of RLK and RLP legume proteins belong to orthologous clusters, which is confirmed in this study, where between 66 to 90% of the RLKs and RLPs per legume species were classified in orthologous clusters. One-third of the evaluated syntenic blocks had shared RLK/RLP genes among both legumes and non-legumes. Among the legumes, between 75 and 98% of the RLK/RLP were present in syntenic blocks. The distribution of chromosomal segments between Phaseolus vulgaris and Vigna unguiculata, two species that diverged ~ 8 mya, were highly similar. Among the RLK/RLP synteny clusters, seven experimentally validated resistance RLK/RLP genes were identified in syntenic blocks. The RLK resistant genes FLS2, BIR2, ERECTA, IOS1, and AtSERK1 from Arabidopsis and SLSERK1 from Solanum lycopersicum were present in different pairwise syntenic blocks among the legume species. Meanwhile, only the LYM1- RLP resistant gene from Arabidopsis shared a syntenic blocks with Glycine max.

CONCLUSIONS

The orthology analysis of the RLK and RLP suggests a dynamic evolution in the legume family, with between 66 to 85% of RLK and 83 to 88% of RLP belonging to orthologous clusters among the species evaluated. In fact, for the 10-species comparison, a lower number of singleton proteins were reported among RLP compared to RLK, suggesting that RLP positions are more physically conserved compared to RLK. The identification of RLK and RLP genes among the synteny blocks in legumes revealed multiple highly conserved syntenic blocks on multiple chromosomes. Additionally, the analysis suggests that P. vulgaris is an appropriate anchor species for comparative genomics among legumes.

Wild soybean SNARE proteins BET1s mediate the subcellular localization of the cytoplasmic receptor-like kinases CRCK1s to modulate salt stress responses

Plants have evolved numerous receptor-like kinases (RLKs) that modulate environmental stress responses. However, little is known regarding soybean (Glycine max) RLKs. We have previously identified that Glycine soja Ca2+ /CAM-binding RLK (GsCBRLK) is involved in salt tolerance. Here, we report that soluble NSF attachment protein receptor proteins BET1s mediate subcellular localization of calmodulin-binding receptor-like cytoplasmic kinases CRCK1s to modulate salt stress responses. Direct interaction between GsCBRLK and GsBET11a was initially identified via yeast two-hybrid and bimolecular fluorescence complementation assays. Further analysis demonstrated conserved interaction between BET1s and CRCK1s. GsCBRLK interacted with all BET1 proteins in wild soybean (Glycine soja) and Arabidopsis, and GsBET11a strongly associated with GsCRCK1a-1d, but slightly with AtCRCK1. In addition, GsBET11a interacted with GsCBRLK via its C-terminal transmembrane domain (TMD), where the entire TMD, not the sequence, was critical for the interaction. Moreover, the N-terminal variable domain (VD) of GsCBRLK was responsible for interacting with GsBET11a, and the intensity of interaction between GsCBRLK/AtCRCK1 and GsBET11a was dependent on VD. Furthermore, GsBET11a was able to mediate the GsCBRLK subcellular localization via direct interaction with VD. Additionally, knockout of AtBET11 or AtBET12 individually did not alter GsCBRLK localization, while GsBET11a expression caused partial internalization of GsCBRLK from the plasma membrane (PM). We further suggest the necessity of GsCBRLK VD for its PM localization via N-terminal truncation assays. Finally, GsBET11a was shown to confer enhanced salt stress tolerance when overexpressed in Arabidopsis and soybean. These results revealed the conserved and direct interaction between BET1s and CRCK1s, and suggested their involvement in salt stress responses.

The Wild Sugarcane and Sorghum Kinomes: Insights Into Expansion, Diversification, and Expression Patterns

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.

Soybean TILLING-by-Sequencing+ reveals the role of novel GmSACPD members in unsaturated fatty acid biosynthesis while maintaining healthy nodules

Developing soybean lines with high levels of stearic acid is a primary goal of the soybean industry. Most high-stearic-acid soybeans carry different GmSACPD-C mutated alleles. However, due to the dual role of GmSACPD-C in seeds and nodule development, all derived deleterious GmSACPD-C mutant alleles are of extremely poor agronomic value because of defective nodulation. The soybean stearoyl-acyl carrier protein desaturase (GmSACPD) gene family is composed of five members. Comparative genomics analysis indicated that SACPD genes were duplicated and derived from a common ancestor that is still present in chlorophytic algae. Synteny analysis showed the presence of segment duplications between GmSACPD-A/GmSACPD-B, and GmSACPD-C/GmSACPD-D. GmSACPD-E was not contained in any duplicated segment and may be the result of tandem duplication. We developed a TILLING by Target Capture Sequencing (Tilling-by-Sequencing+) technology, a versatile extension of the conventional TILLING by sequencing, and successfully identified 12, 14, and 18 ethyl methanesulfonate mutants at the GmSACPD-A, GmSACPD-B, and GmSACPD-D genes, respectively. Functional analysis of all identified mutants revealed an unprecedented role of GmSACPD-A, GmSACPD-B, and GmSACPD-D in unsaturated fatty acid biosynthesis without affecting nodule development and structure. This discovery will positively impact the development of high-stearic-acid lines to enhance soybean nutritional value without potential developmental tradeoffs.

The wild strawberry kinome: identification, classification and transcript profiling of protein kinases during development and in response to gray mold infection

Background

Protein kinases (PKs) play an important role in signaling cascades and are one of the largest and most conserved protein super families in plants. Despite their importance, the woodland strawberry (Fragaria vesca) kinome and expression patterns of PK genes remain to be characterized.

Results

Here, we report on the identification and classification of 954 Fragaria vesca PK genes, which were classified into nine groups and 124 gene families. These genes were distributed unevenly among the seven chromosomes, and the number of introns per gene varied from 0 to 47. Almost half of the putative PKs were predicted to localize to the nucleus and 24.6% were predicted to localize to the cell membrane. The expansion of the woodland strawberry PK gene family occurred via different duplication mechanisms and tandem duplicates occurred relatively late as compared to other duplication types. Moreover, we found that tandem and transposed duplicated PK gene pairs had undergone stronger diversifying selection and evolved relatively faster than WGD genes. The GO enrichment and transcriptome analysis implicates the involvement of strawberry PK genes in multiple biological processes and molecular functions in differential tissues, especially in pollens. Finally, 109 PKs, mostly the receptor-like kinases (RLKs), were found transcriptionally responsive to Botrytis cinerea infection.

Conclusions

The findings of this research expand the understanding of the evolutionary dynamics of PK genes in plant species and provide a potential link between cell signaling pathways and pathogen attack.

A pathogenesis-related protein GmPR08-Bet VI promotes a molecular interaction between the GmSHMT08 and GmSNAP18 in resistance to Heterodera glycines

Soybean cyst nematode (SCN, Heterodera glycines) is the most devastating pest affecting soybean production worldwide. SCN resistance requires both the GmSHMT08 and the GmSNAP18 in 'Peking'-type resistance. Here, we describe the molecular interaction between GmSHMT08 and GmSNAP18, which is potentiated by a pathogenesis-related protein GmPR08-Bet VI. Like GmSNAP18 and GmSHMT08, GmPR08-Bet VI expression was induced in response to SCN and its overexpression decreased SCN cysts by 65% in infected transgenic soybean roots. Overexpression of GmPR08-Bet VI did not have an effect on SCN resistance when the two cytokinin-binding sites in GmPR08-Bet VI were mutated, indicating a new role of GmPR08-Bet VI in SCN resistance. GmPR08-Bet VI was mapped to a QTL for resistance to SCN using different mapping populations. GmSHMT08, GmSNAP18 and GmPR08-Bet VI localize to the cytosol and plasma membrane. GmSNAP18 expression and localization hyper-accumulated at the plasma membrane and was specific to the root cells surrounding the nematode in SCN-resistant soybeans. Genes encoding key components of the salicylic acid signalling pathway were induced under SCN infection. GmSNAP18 and GmPR08-Bet VI were also induced under salicylic acid and cytokinin exogenous treatments, while GmSHMT08 was induced only when the resistant GmSNAP18 was present, pointing to the presence of a molecular crosstalk between SCN-resistant genes and defence genes. Expression analysis of GmSHMT08 and GmSNAP18 identified the need of a minimum expression requirement to trigger the SCN resistance reaction. These results provide insight into a new response mechanism towards plant nematode resistance involving haplotype compatibility, gene dosage and hormone signalling.

Genome reorganization of the GmSHMT gene family in soybean showed a lack of functional redundancy in resistance to soybean cyst nematode

In soybeans, eighteen members constitute the serine hydroxymethyltransferase (GmSHMT) gene family, of which the cytosolic-targeted GmSHMT08c member has been reported to mediate resistance to soybean cyst nematode (SCN). This work presents a comprehensive study of the SHMT gene family members, including synteny, phylogeny, subcellular localizations, haplotypes, protein homology modeling, mutational, and expression analyses. Phylogenetic analysis showed that SHMT genes are divided into four classes reflecting their subcellular distribution (cytosol, nucleus, mitochondrion, and chloroplast). Subcellular localization of selected GmSHMT members supports their in-silico predictions and phylogenetic distribution. Expression and functional analyses showed that GmSHMT genes display many overlapping, but some divergent responses during SCN infection. Furthermore, mutational analysis reveals that all isolated EMS mutants that lose their resistance to SCN carry missense and nonsense mutations at the GmSHMT08c, but none of the Gmshmt08c mutants carried mutations in the other GmSHMT genes. Haplotype clustering analysis using the whole genome resequencing data from a collection of 106 diverse soybean germplams (15X) was performed to identify allelic variants and haplotypes within the GmSHMT gene family. Interestingly, only the cytosolic-localized GmSHMT08c presented SNP clusters that were associated with SCN resistance, supporting our mutational analysis. Although eight GmSHMT members respond to the nematode infestation, functional and mutational analysis has shown the absence of functional redundancy in resistance to SCN. Structural analysis and protein homology modeling showed the presence of spontaneous mutations at important residues within the GmSHMT proteins, suggesting the presence of altered enzyme activities based on substrate affinities. Due to the accumulation of mutations during the evolution of the soybean genome, the other GmSHMT members have undergone neofunctionalization and subfunctionalization events.

Genome-wide analysis of the soybean CRK-family and transcriptional regulation by biotic stress signals triggering plant immunity

Cysteine-rich receptor-like kinases (CRKs) are transmembrane proteins that exhibit ectodomains containing the domain of unknown function 26 (DUF26). The CRKs form a large subfamily of receptor-like kinases in plants, and their possible functions remain to be elucidated. Several lines of evidence suggest that CRKs play important roles in plant defense responses to environmental stress, including plant immunity. We performed a genome-wide analysis of CRK encoding genes in soybean (Glycine max). We found 91 GmCRKs distributed in 16 chromosomes, and identified several tandem and segmental duplications, which influenced the expansion of this gene family. According to our phylogenetic analysis, GmCRKs are grouped in four clades. Furthermore, 12% of the members exhibited GmCRKs with a duplicated bi-modular organization of the ectodomains, containing four DUF26 domains. Expression analysis of GmCRKs was performed by exploring publicly available databases, and by RT-qPCR analysis of selected genes in soybean leaves responding to biotic stress signals. GmCRKs exhibited diverse expression patterns in leaves, stems, roots, and other tissues. Some of them were highly expressed in only one type of tissue, suggesting predominant roles in specific tissues. Furthermore, several GmCRKs were induced with PAMPs, DAMPs and the pathogens Phakopsora pachyrhizi and Phytophthora sojae. Expression profiles of several GmCRKs encoding highly similar proteins exhibited antagonist modes of regulation. The results suggest a fine-tuning control of GmCRKs transcriptional regulation in response to external stimuli, including PAMPs and DAMPs. This study offers a comprehensive view of the GmCRKs family in soybean, and provides a foundation for evolutionary and functional analysis of this family of plant proteins involved in the perception of pathogens and activation of plant immunity.

The kinome of pineapple: catalog and insights into functions in crassulacean acid metabolism plants

BACKGROUND

Crassulacean acid metabolism (CAM) plants use water 20-80% more efficiently by shifting stomata opening and primary CO₂ uptake and fixation to the nighttime. Protein kinases (PKs) play pivotal roles in this biological process. However, few PKs have been functionally analyzed precisely due to their abundance and potential functional redundancy (caused by numerous gene duplications).

RESULTS

In this study, we systematically identified a total of 758 predicted PK genes in the genome of a CAM plant, pineapple (Ananas comosus). The pineapple kinome was classified into 20 groups and 116 families based on the kinase domain sequences. The RLK was the largest group, containing 480 members, and over half of them were predicted to locate at the plasma membrane. Both segmental and tandem duplications make important contributions to the expansion of pineapple kinome based on the synteny analysis. Ka/Ks ratios showed all of the duplication events were under purifying selection. The global expression analysis revealed that pineapple PKs exhibit different tissue-specific and diurnal expression patterns. Forty PK genes in a cluster performed higher expression levels in green leaf tip than in white leaf base, and fourteen of them had strong differential expression patterns between the photosynthetic green leaf tip and the non-photosynthetic white leaf base tissues.

CONCLUSIONS

Our findings provide insights into the evolution and biological function of pineapple PKs and a foundation for further functional analysis of PKs in CAM plants. The gene duplication, expression, and coexpression analysis helped us to rapidly identify the key candidates in pineapple kinome, which may play roles in the carbon fixation process in pineapple and help engineering CAM pathway into C3 crops for improved drought tolerance.

Genome-wide classification, evolutionary analysis and gene expression patterns of the kinome in Gossypium

The protein kinase (PK, kinome) family is one of the largest families in plants and regulates almost all aspects of plant processes, including plant development and stress responses. Despite their important functions, comprehensive functional classification, evolutionary analysis and expression patterns of the cotton PK gene family has yet to be performed on PK genes. In this study, we identified the cotton kinomes in the Gossypium raimondii, Gossypium arboretum, Gossypium hirsutum and Gossypium barbadense genomes and classified them into 7 groups and 122-24 subfamilies using software HMMER v3.0 scanning and neighbor-joining (NJ) phylogenetic analysis. Some conserved exon-intron structures were identified not only in cotton species but also in primitive plants, ferns and moss, suggesting the significant function and ancient origination of these PK genes. Collinearity analysis revealed that 16.6 million years ago (Mya) cotton-specific whole genome duplication (WGD) events may have played a partial role in the expansion of the cotton kinomes, whereas tandem duplication (TD) events mainly contributed to the expansion of the cotton RLK group. Synteny analysis revealed that tetraploidization of G. hirsutum and G. barbadense contributed to the expansion of G. hirsutum and G. barbadense PKs. Global expression analysis of cotton PKs revealed stress-specific and fiber development-related expression patterns, suggesting that many cotton PKs might be involved in the regulation of the stress response and fiber development processes. This study provides foundational information for further studies on the evolution and molecular function of cotton PKs.

Ectopic Expression of GsSRK in Medicago sativa Reveals Its Involvement in Plant Architecture and Salt Stress Responses

Receptor-like kinases (RLK) play fundamental roles in plant growth and stress responses. Compared with other RLKs, little information is provided concerning the S-locus LecRLK subfamily, which is characterized by an extracellular G-type lectin domain and an S-locus-glycop domain. Until now, the function of the G-type lectin domain is still unknown. In a previous research, we identified a Glycine soja S-locus LecRLK gene GsSRK, which conferred increased salt stress tolerance in transgenic Arabidopsis. In this study, to investigate the role of the G-type lectin domain and to breed transgenic alfalfa with superior salt stress tolerance, we transformed the full-length GsSRK (GsSRK-f) and a truncated version of GsSRK (GsSRK-t) deleting the G-type lectin domain into alfalfa. Our results showed that overexpression of GsSRK-t, but not GsSRK-f, resulted in changes of plant architecture, as evidenced by more branches but shorter shoots of GsSRK-t transgenic alfalfa, indicating a potential role of the extracellular G-type lectin domain in regulating plant architecture. Furthermore, we also found that transgenic alfalfa overexpressing either GsSRK-f or GsSRK-t showed increased salt stress tolerance, and GsSRK-t transgenic alfalfa displayed better growth (more branches and higher fresh weight) than GsSRK-f lines under salt stress. In addition, our results suggested that both GsSRK-f and GsSRK-t were involved in ion homeostasis, ROS scavenging, and osmotic regulation. Under salt stress, the Na+ content in the transgenic lines was significantly lower, while the K+ content was slightly higher than that in WT. Moreover, the transgenic lines displayed reduced ion leakage and MDA content, but increased SOD activity and proline content than WT. Notably, no obvious difference in these physiological indices was observed between GsSRK-f and GsSRK-t transgenic lines, implying that deletion of the GsSRK G-type lectin domain does not affect its physiological function in salt stress responses. In conclusion, results in this research reveal the dual role of GsSRK in regulating both plant architecture and salt stress responses.

Phytopathogen-induced changes to plant methylomes

DNA methylation is a dynamic and reversible type of epigenetic mark that contributes to cellular physiology by affecting transcription activity, transposon mobility and genome stability. When plants are infected with pathogens, plant DNA methylation patterns can change, indicating an epigenetic interplay between plant host and pathogen. In most cases methylation can change susceptibility. While DNA hypomethylation appears to be a common phenomenon during the susceptible interaction, the levels and patterns of hypomethylation in transposable elements and genic regions may mediate distinct responses against various plant pathogens. The effect of DNA methylation on the plant immune response and other cellular activities and molecular functions is established by localized differential DNA methylation via cis-regulatory mechanisms as well as through trans-acting mechanisms. Understanding the epigenetic differences that control the phenotypic variations between susceptible and resistant interactions should facilitate the identification of new sources of resistance mediated by epigenetic mechanisms, which can be exploited to endow pathogen resistance to crops.

The grapevine kinome: annotation, classification and expression patterns in developmental processes and stress responses

Protein kinases (PKs) have evolved as the largest family of molecular switches that regulate protein activities associated with almost all essential cellular functions. Only a fraction of plant PKs, however, have been functionally characterized even in model plant species. In the present study, the entire grapevine kinome was identified and annotated using the most recent version of the grapevine genome. A total of 1168 PK-encoding genes were identified and classified into 20 groups and 121 families, with the RLK-Pelle group being the largest, with 872 members. The 1168 kinase genes were unevenly distributed over all 19 chromosomes, and both tandem and segmental duplications contributed to the expansion of the grapevine kinome, especially of the RLK-Pelle group. Ka/Ks values indicated that most of the tandem and segmental duplication events were under purifying selection. The grapevine kinome families exhibited different expression patterns during plant development and in response to various stress treatments, with many being coexpressed. The comprehensive annotation of grapevine kinase genes, their patterns of expression and coexpression, and the related information facilitate a more complete understanding of the roles of various grapevine kinases in growth and development, responses to abiotic stress, and evolutionary history.

Genome-wide characterization, evolution, and expression analysis of the leucine-rich repeat receptor-like protein kinase (LRR-RLK) gene family in Rosaceae genomes

BACKGROUND

Leucine-rich repeat receptor-like protein kinase (LRR-RLK) is the largest gene family of receptor-like protein kinases (RLKs) and actively participates in regulating the growth, development, signal transduction, immunity, and stress responses of plants. However, the patterns of LRR-RLK gene family evolution in the five main Rosaceae species for which genome sequences are available have not yet been reported. In this study, we performed a comprehensive analysis of LRR-RLK genes for five Rosaceae species: Fragaria vesca (strawberry), Malus domestica (apple), Pyrus bretschneideri (Chinese white pear), Prunus mume (mei), and Prunus persica (peach), which contained 201, 244, 427, 267, and 258 LRR-RLK genes, respectively.

RESULTS

All LRR-RLK genes were further grouped into 23 subfamilies based on the hidden Markov models approach. RLK-Pelle_LRR-XII-1, RLK-Pelle_LRR-XI-1, and RLK-Pelle_LRR-III were the three largest subfamilies. Synteny analysis indicated that there were 236 tandem duplicated genes in the five Rosaceae species, among which subfamilies XII-1 (82 genes) and XI-1 (80 genes) comprised 68.6%.

CONCLUSIONS

Our results indicate that tandem duplication made a large contribution to the expansion of the subfamilies. The gene expression, tissue-specific expression, and subcellular localization data revealed that LRR-RLK genes were differentially expressed in various organs and tissues, and the largest subfamily XI-1 was highly expressed in all five Rosaceae species, suggesting that LRR-RLKs play important roles in each stage of plant growth and development. Taken together, our results provide an overview of the LRR-RLK family in Rosaceae genomes and the basis for further functional studies.

Genome-wide identification, classification, evolutionary analysis and gene expression patterns of the protein kinase gene family in wheat and Aegilops tauschii

In this study we systematically identified and classified PKs in Triticum aestivum, Triticum urartu and Aegilops tauschii. Domain distribution and exon-intron structure analyses of PKs were performed, and we found conserved exon-intron structures within the exon phases in the kinase domain. Collinearity events were determined, and we identified various T. aestivum PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed that some PKs might participate in the signaling pathways of stress response and developmental processes. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of Fusarium graminearum to validate the prediction of microarray. The protein kinase (PK) gene superfamily is one of the largest families in plants and participates in various plant processes, including growth, development, and stress response. To better understand wheat PKs, we conducted genome-wide identification, classification, evolutionary analysis and expression profiles of wheat and Ae. tauschii PKs. We identified 3269, 1213 and 1448 typical PK genes in T. aestivum, T. urartu and Ae. tauschii, respectively, and classified them into major groups and subfamilies. Domain distributions and gene structures were analyzed and visualized. Some conserved intron-exon structures within the conserved kinase domain were found in T. aestivum, T. urartu and Ae. tauschii, as well as the primitive land plants Selaginella moellendorffii and Physcomitrella patens, revealing the important roles and conserved evolutionary history of these PKs. We analyzed the collinearity events of T. aestivum PKs and identified PKs from polyploidizations and tandem duplication events. Global expression pattern analysis of T. aestivum PKs revealed tissue-specific and stress-specific expression profiles, hinting that some wheat PKs may regulate abiotic and biotic stress response signaling pathways. QRT-PCR of 15 selected PKs were performed under drought treatment and with infection of F. graminearum to validate the prediction of microarray. Our results will provide the foundational information for further studies on the molecular functions of wheat PKs.

Characterization of Disease Resistance Loci in the USDA Soybean Germplasm Collection Using Genome-Wide Association Studies

Genetic resistance is a key strategy for disease management in soybean. Over the last 50 years, soybean germplasm has been phenotyped for resistance to many pathogens, resulting in the development of disease-resistant elite breeding lines and commercial cultivars. While biparental linkage mapping has been used to identify disease resistance loci, genome-wide association studies (GWAS) using high-density and high-quality markers such as single nucleotide polymorphisms (SNPs) has become a powerful tool to associate molecular markers and phenotypes. The objective of our study was to provide a comprehensive understanding of disease resistance in the United States Department of Agriculture Agricultural Research Service Soybean Germplasm Collection by using phenotypic data in the public Germplasm Resources Information Network and public SNP data (SoySNP50K). We identified SNPs significantly associated with disease ratings from one bacterial disease, five fungal diseases, two diseases caused by nematodes, and three viral diseases. We show that leucine-rich repeat (LRR) receptor-like kinases and nucleotide-binding site-LRR candidate resistance genes were enriched within the linkage disequilibrium regions of the significant SNPs. We review and present a global view of soybean resistance loci against multiple diseases and discuss the power and the challenges of using GWAS to discover disease resistance in soybean.

Systems medicine approaches for peptide array-based protein kinase profiling: progress and prospects

INTRODUCTION

Pharmacological manipulation of signalling pathways is becoming an increasingly important avenue for the rational clinical management of disease but is hampered by a lack of technologies that allow the generation of comprehensive descriptions of cellular signalling.

AREAS COVERED

Herein, the authors discuss the potential of peptide array-based kinome profiling for evaluating cellular signalling in the context of drug discovery. Expert commentary: Genomic and proteomic approaches have been of significant value to our elucidation of the molecular mechanisms that govern physiology. However, an equally, if not more important goal, is to define those proteins that participate in signalling pathways that ultimately control cell fate, especially kinases. Traditional genetic and biochemical approaches can certainly provide answers here, but for technical and practical reasons, are typically pursued one gene or pathway at a time. A more comprehensive approach is one in which peptide arrays of kinase-specific substrates are incubated with cell lysates and (33)P-ATP generating comprehensive descriptions, or where arrays are interrogated with phosphospecific antibodies. Both approaches allow analysis of cellular signalling without a priori assumptions to possibly influenced pathways.

Perception of pathogenic or beneficial bacteria and their evasion of host immunity: pattern recognition receptors in the frontline

Plants are continuously monitoring the presence of microorganisms to establish an adapted response. Plants commonly use pattern recognition receptors (PRRs) to perceive microbe- or pathogen-associated molecular patterns (MAMPs/PAMPs) which are microorganism molecular signatures. Located at the plant plasma membrane, the PRRs are generally receptor-like kinases (RLKs) or receptor-like proteins (RLPs). MAMP detection will lead to the establishment of a plant defense program called MAMP-triggered immunity (MTI). In this review, we overview the RLKs and RLPs that assure early recognition and control of pathogenic or beneficial bacteria. We also highlight the crucial function of PRRs during plant-microbe interactions, with a special emphasis on the receptors of the bacterial flagellin and peptidoglycan. In addition, we discuss the multiple strategies used by bacteria to evade PRR-mediated recognition.