KRSA: An R package and R Shiny web application for an end-to-end upstream kinase analysis of kinome array data.

Phosphorylation by serine-threonine and tyrosine kinases is critical for determining protein function. Array-based platforms for measuring reporter peptide signal levels allow for differential phosphorylation analysis between conditions for distinct active kinases. Peptide array technologies like the PamStation12 from PamGene allow for generating high-throughput, multi-dimensional, and complex functional proteomics data. As the adoption rate of such technologies increases, there is an imperative need for software tools that streamline the process of analyzing such data. We present Kinome Random Sampling Analyzer (KRSA), an R package and R Shiny web-application for analyzing kinome array data to help users better understand the patterns of functional proteomics in complex biological systems. KRSA is an All-In-One tool that reads, formats, fits models, analyzes, and visualizes PamStation12 kinome data. While the underlying algorithm has been experimentally validated in previous publications, we demonstrate KRSA workflow on dorsolateral prefrontal cortex (DLPFC) in male (n = 3) and female (n = 3) subjects to identify differential phosphorylation signatures and upstream kinase activity. Kinase activity differences between males and females were compared to a previously published kinome dataset (11 female and 7 male subjects) which showed similar global phosphorylation signals patterns.

Leveraging breeding programs and genomic data in Norway spruce (Picea abies L. Karst) for GWAS analysis.

BACKGROUND: Genome-wide association studies (GWAS) identify loci underlying the variation of complex traits. One of the main limitations of GWAS is the availability of reliable phenotypic data, particularly for long-lived tree species. Although an extensive amount of phenotypic data already exists in breeding programs, accounting for its high heterogeneity is a great challenge. We combine spatial and factor-analytics analyses to standardize the heterogeneous data from 120 field experiments of 483,424 progenies of Norway spruce to implement the largest reported GWAS for trees using 134 605 SNPs from exome sequencing of 5056 parental trees. RESULTS: We identify 55 novel quantitative trait loci (QTLs) that are associated with phenotypic variation. The largest number of QTLs is associated with the budburst stage, followed by diameter at breast height, wood quality, and frost damage. Two QTLs with the largest effect have a pleiotropic effect for budburst stage, frost damage, and diameter and are associated with MAP3K genes. Genotype data called from exome capture, recently developed SNP array and gene expression data indirectly support this discovery. CONCLUSION: Several important QTLs associated with growth and frost damage have been verified in several southern and northern progeny plantations, indicating that these loci can be used in QTL-assisted genomic selection. Our study also demonstrates that existing heterogeneous phenotypic data from breeding programs, collected over several decades, is an important source for GWAS and that such integration into GWAS should be a major area of inquiry in the future.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6.

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.

Characterization of testis-specific serine/threonine kinase 1-like (TSSK1-like) gene and expression patterns in diploid and triploid Pacific abalone (Haliotis discus hannai; Gastropoda; Mollusca) males.

Testis-specific serine/threonine kinase 1-like (TSSK1-like), which plays important roles in late-phase spermatogenesis and male fertility, was characterized in Pacific abalone Haliotis discus hannai, an important commercial marine gastropod. Further, its expression patterns were assessed in diploid and induced triploid males showing differential degrees of testis maturation. Abalone TSSK1-like shared conserved structural features with mammalian TSSK1s and other potential metazoan orthologs, especially regarding the catalytic STKc domain. Phylogenetically, abalone TSSK1-like displayed a genetic affiliation with its molluscan TSSK1-like orthologs and human TSSK1. Additionally, abalone TSSK1-like gene showed a tetrapartite exon-intron organization, unlike the intronless structure of most amniotic tetrapodian TSSK1s. Molecular phylogenetic analysis in the metazoan lineage suggested a possible revision in the origin of the earliest ancestral TSSK1. Further, abalone TSSK1-like showed testis-predominant expression, which was significantly influenced by both age and seasonal reproductive cycles. Comparative expression analyses between diploid and triploid abalone males suggested that robust TSSK1-like expression occurred primarily at the post-meiotic stage. Additionally, RT-PCR assay indicates that mature abalone sperms retain TSSK1-like transcripts after release. Taken together, this study provides useful insights for further studies to assess male reproduction and sterility and/or partial fertility of induced male triploidy in abalone species.

CaM Kinase: Still Inspiring at 40.

The Ca2+/calmodulin (CaM)-dependent protein kinase II (CaMKII) was touted as a memory molecule, even before its involvement in long-term potentiation (LTP) was shown. The enzyme has not disappointed, with subsequent demonstrations of remarkable structural and regulatory properties. Its neuronal functions now extend to long-term depression (LTD), and last year saw the first direct evidence for memory storage by CaMKII. Although CaMKII may have taken the spotlight, it is a member of a large family of diverse and interesting CaM kinases. Our aim is to place CaMKII in context of the other CaM kinases and then review certain aspects of this kinase that are of current interest.

The crystal structure of the protein kinase HIPK2 reveals a unique architecture of its CMGC-insert region.

The homeodomain-interacting protein kinase (HIPK) family is comprised of four nuclear protein kinases, HIPK1-4. HIPK proteins phosphorylate a diverse range of transcription factors involved in cell proliferation, differentiation, and apoptosis. HIPK2, thus far the best-characterized member of this largely understudied family of protein kinases, plays a role in the activation of p53 in response to DNA damage. Despite this tumor-suppressor function, HIPK2 is also found overexpressed in several cancers, and its hyperactivation causes chronic fibrosis. There are currently no structures of HIPK2 or of any other HIPK kinase. Here, we report the crystal structure of HIPK2's kinase domain bound to CX-4945, a casein kinase 2α (CK2α) inhibitor currently in clinical trials against several cancers. The structure, determined at 2.2 Å resolution, revealed that CX-4945 engages the HIPK2 active site in a hybrid binding mode between that seen in structures of CK2α and Pim1 kinases. The HIPK2 kinase domain crystallized in the active conformation, which was stabilized by phosphorylation of the activation loop. We noted that the overall kinase domain fold of HIPK2 closely resembles that of evolutionarily related dual-specificity tyrosine-regulated kinases (DYRKs). Most significant structural differences between HIPK2 and DYRKs included an absence of the regulatory N-terminal domain and a unique conformation of the CMGC-insert region and of a newly defined insert segment in the αC-ß4 loop. This first crystal structure of HIPK2 paves the way for characterizing the understudied members of the HIPK family and for developing HIPK2-directed therapies for managing cancer and fibrosis.

MOB (Mps one Binder) Proteins in the Hippo Pathway and Cancer.

The family of MOBs (monopolar spindle-one-binder proteins) is highly conserved in the eukaryotic kingdom. MOBs represent globular scaffold proteins without any known enzymatic activities. They can act as signal transducers in essential intracellular pathways. MOBs have diverse cancer-associated cellular functions through regulatory interactions with members of the NDR/LATS kinase family. By forming additional complexes with serine/threonine protein kinases of the germinal centre kinase families, other enzymes and scaffolding factors, MOBs appear to be linked to an even broader disease spectrum. Here, we review our current understanding of this emerging protein family, with emphases on post-translational modifications, protein-protein interactions, and cellular processes that are possibly linked to cancer and other diseases. In particular, we summarise the roles of MOBs as core components of the Hippo tissue growth and regeneration pathway.

Lotus SHAGGY-like kinase 1 is required to suppress nodulation in Lotus japonicus.

Glycogen synthase kinase/SHAGGY-like kinases (SKs) are a highly conserved family of signaling proteins that participate in many developmental, cell-differentiation, and metabolic signaling pathways in plants and animals. Here, we investigate the involvement of SKs in legume nodulation, a process requiring the integration of multiple signaling pathways. We describe a group of SKs in the model legume Lotus japonicus (LSKs), two of which respond to inoculation with the symbiotic nitrogen-fixing bacterium Mesorhizobium loti. RNAi knock-down plants and an insertion mutant for one of these genes, LSK1, display increased nodulation. Ηairy-root lines overexpressing LSK1 form only marginally fewer mature nodules compared with controls. The expression levels of genes involved in the autoregulation of nodulation (AON) mechanism are affected in LSK1 knock-down plants at low nitrate levels, both at early and late stages of nodulation. At higher levels of nitrate, these same plants show the opposite expression pattern of AON-related genes and lose the hypernodulation phenotype. Our findings reveal an additional role for the versatile SK gene family in integrating the signaling pathways governing legume nodulation, and pave the way for further study of their functions in legumes.

Leucine-rich repeat receptor-like gene screen reveals that Nicotiana RXEG1 regulates glycoside hydrolase 12 MAMP detection.

Activation of innate immunity by membrane-localized receptors is conserved across eukaryotes. Plant genomes contain hundreds of such receptor-like genes and those encoding proteins with an extracellular leucine-rich repeat (LRR) domain represent the largest family. Here, we develop a high-throughput approach to study LRR receptor-like genes on a genome-wide scale. In total, 257 tobacco rattle virus-based constructs are generated to target 386 of the 403 identified LRR receptor-like genes in Nicotiana benthamiana for silencing. Using this toolkit, we identify the LRR receptor-like protein Response to XEG1 (RXEG1) that specifically recognizes the glycoside hydrolase 12 protein XEG1. RXEG1 associates with XEG1 via the LRR domain in the apoplast and forms a complex with the LRR receptor-like kinases BAK1 and SOBIR1 to transduce the XEG1-induced defense signal. Thus, this genome-wide silencing assay is demonstrated to be an efficient toolkit to pinpoint new immune receptors, which will contribute to developing durable disease resistance.

Identification and sequence analysis of putative Sulla species nod factor receptor.

In legumes, LysM domains of receptors-like kinases (RLKs) mediate rhizobial NFs perception; which are required for infection and establishment of symbiosis without triggering the host immune response. In this study, we identify the LysM extracellular domain sequences of putative Sulla species Nod factor receptors (S. pallida, S. capitata and S. coronaria). The Blast search displayed high identity percentages with genes encoding LjNFR5-like of several legumes. Phylogenetic trees were built using the partial nod factor receptor and predicted amino acid sequences, which grouped Sulla in a separate clade. The multiple alignments of the LysM2 domains revealed that amino acids found to be important in other legume species are not conserved in Sulla species. Further examination of the predicted proteins sequences (LysM2 domain) showed that the three species were different in the two crucial sites for Nod factor perception.

Receptor-Like Kinase LYK9 in Pisum sativum L. Is the CERK1-Like Receptor that Controls Both Plant Immunity and AM Symbiosis Development.

Plants are able to discriminate and respond to structurally related chitooligosaccharide (CO) signals from pathogenic and symbiotic fungi. In model plants Arabidopsis thaliana and Oryza sativa LysM-receptor like kinases (LysM-RLK) AtCERK1 and OsCERK1 (chitin elicitor receptor kinase 1) were shown to be involved in response to CO signals. Based on phylogenetic analysis, the pea Pisum sativum L. LysM-RLK PsLYK9 was chosen as a possible candidate given its role on the CERK1-like receptor. The knockdown regulation of the PsLyk9 gene by RNA interference led to increased susceptibility to fungal pathogen Fusarium culmorum. Transcript levels of PsPAL2, PsPR10 defense-response genes were significantly reduced in PsLyk9 RNAi roots. PsLYK9's involvement in recognizing short-chain COs as most numerous signals of arbuscular mycorrhizal (AM) fungi, was also evaluated. In transgenic roots with PsLyk9 knockdown treated with short-chain CO5, downregulation of AM symbiosis marker genes (PsDELLA3, PsNSP2, PsDWARF27) was observed. These results clearly indicate that PsLYK9 appears to be involved in the perception of COs and subsequent signal transduction in pea roots. It allows us to conclude that PsLYK9 is the most likely CERK1-like receptor in pea to be involved in the control of plant immunity and AM symbiosis formation.

Evolution of protein phosphorylation across 18 fungal species.

Living organisms have evolved protein phosphorylation, a rapid and versatile mechanism that drives signaling and regulates protein function. We report the phosphoproteomes of 18 fungal species and a phylogenetic-based approach to study phosphosite evolution. We observe rapid divergence, with only a small fraction of phosphosites conserved over hundreds of millions of years. Relative to recently acquired phosphosites, ancient sites are enriched at protein interfaces and are more likely to be functionally important, as we show for sites on H2A1 and eIF4E. We also observe a change in phosphorylation motif frequencies and kinase activities that coincides with the whole-genome duplication event. Our results provide an evolutionary history for phosphosites and suggest that rapid evolution of phosphorylation can contribute strongly to phenotypic diversity.

Atg1-independent induction of autophagy by the Drosophila Ulk3 homolog, ADUK.

Although canonical autophagy regulation requires a multi-protein complex centered on the Ser/Thr-kinase Atg1 (mammalian Ulk1/2), alternative signals can induce autophagy independent of Atg1 through unknown mechanisms. Here we identify the Drosophila Ulk3 ortholog, another Drosophila Unc-51-like kinase (ADUK), as an Atg1-independent autophagy inducer. ADUK interacts with Atg1 complex members Atg13 and 200 kDa FAK family kinase-interacting protein, and requires Atg13 but not Atg1 for autophagy induction. Loss of ADUK shortens adult lifespan and reduces the autophagic response to a chemical stressor, dimethyl sulfoxide. However, ADUK is not required for autophagy induction by Atg1-dependent nutrient or developmental cues. Atg1 and ADUK/Ulk3 thus represent alternative catalytic components of a shared autophagy kinase complex.

LRR-RLK family from two Citrus species: genome-wide identification and evolutionary aspects.

BACKGROUND: Leucine-rich repeat receptor-like kinases (LRR-RLKs) represent the largest subfamily of plant RLKs. The functions of most LRR-RLKs have remained undiscovered, and a few that have been experimentally characterized have been shown to have important roles in growth and development as well as in defense responses. Although RLK subfamilies have been previously studied in many plants, no comprehensive study has been performed on this gene family in Citrus species, which have high economic importance and are frequent targets for emerging pathogens. In this study, we performed in silico analysis to identify and classify LRR-RLK homologues in the predicted proteomes of Citrus clementina (clementine) and Citrus sinensis (sweet orange). In addition, we used large-scale phylogenetic approaches to elucidate the evolutionary relationships of the LRR-RLKs and further narrowed the analysis to the LRR-XII group, which contains several previously described cell surface immune receptors. RESULTS: We built integrative protein signature databases for Citrus clementina and Citrus sinensis using all predicted protein sequences obtained from whole genomes. A total of 300 and 297 proteins were identified as LRR-RLKs in C. clementina and C. sinensis, respectively. Maximum-likelihood phylogenetic trees were estimated using Arabidopsis LRR-RLK as a template and they allowed us to classify Citrus LRR-RLKs into 16 groups. The LRR-XII group showed a remarkable expansion, containing approximately 150 paralogs encoded in each Citrus genome. Phylogenetic analysis also demonstrated the existence of two distinct LRR-XII clades, each one constituted mainly by RD and non-RD kinases. We identified 68 orthologous pairs from the C. clementina and C. sinensis LRR-XII genes. In addition, among the paralogs, we identified a subset of 78 and 62 clustered genes probably derived from tandem duplication events in the genomes of C. clementina and C. sinensis, respectively. CONCLUSIONS: This work provided the first comprehensive evolutionary analysis of the LRR-RLKs in Citrus. A large expansion of LRR-XII in Citrus genomes suggests that it might play a key role in adaptive responses in host-pathogen co-evolution, related to the perennial life cycle and domestication of the citrus crop species.

Comprehensive phylogeny, biogeography and new classification of the diverse bee tribe Megachilini: Can we use DNA barcodes in phylogenies of large genera?

Classification and evolutionary studies of particularly speciose clades pose important challenges, as phylogenetic analyses typically sample a small proportion of the existing diversity. We examine here one of the largest bee genera, the genus Megachile - the dauber and leafcutting bees. Besides presenting a phylogeny based on five nuclear genes (5480 aligned nucleotide positions), we attempt to use the phylogenetic signal of mitochondrial DNA barcodes, which are rapidly accumulating and already include a substantial proportion of the known species diversity in the genus. We used barcodes in two ways: first, to identify particularly divergent lineages and thus to guide taxon sampling in our nuclear phylogeny; second, to augment taxon sampling by combining nuclear markers (as backbone for ancient divergences) with DNA barcodes. Our results indicate that DNA barcodes bear phylogenetic signal limited to very recent divergences (3-4 my before present). Sampling within clades of very closely related species may be augmented using this technique, but our results also suggest statistically supported, but incongruent placements of some taxa. However, the addition of one single nuclear gene (LW-rhodopsin) to the DNA barcode data was enough to recover meaningful placement with high clade support values for nodes up to 15 million years old. We discuss different proposals for the generic classification of the tribe Megachilini. Finding a classification that is both in agreement with our phylogenetic hypotheses and practical in terms of diagnosability is particularly challenging as our analyses recover several well-supported clades that include morphologically heterogeneous lineages. We favour a classification that recognizes seven morphologically well-delimited genera in Megachilini: Coelioxys, Gronoceras, Heriadopsis, Matangapis, Megachile, Noteriades and Radoszkowskiana. Our results also lead to the following classification changes: the groups known as Dinavis, Neglectella, Eurymella and Phaenosarus are reestablished as valid subgenera of the genus Megachile, while the subgenus Alocanthedon is placed in synonymy with M. (Callomegachile), the subgenera Parachalicodoma and Largella with M. (Pseudomegachile), Anodonteutricharaea with M. (Paracella), Platysta with M. (Eurymella), and Grosapis and Eumegachile with M. (Megachile) (new synonymies). In addition, we use maximum likelihood reconstructions of ancestral geographic ranges to infer the origin of the tribe and reconstruct the main dispersal routes explaining the current, cosmopolitan distribution of this genus.

MAGED1 is a novel regulator of a select subset of bHLH PAS transcription factors.

Transcription factors of the basic helix-loop-helix (bHLH) PER-ARNT-SIM (PAS) family generally have critical and nonredundant biological roles, but some bHLH PAS proteins compete for common cofactors or recognise similar DNA elements. Identifying factors that regulate function of bHLH PAS proteins, particularly in cells where multiple family members are coexpressed, is important for understanding bHLH PAS factor biology. This study identifies and characterises a novel interaction between melanoma-associated antigen D1 (MAGED1) and select members of the bHLH PAS transcription factor family. MAGED1 binds and positively regulates the transcriptional activity of family members SIM1, SIM2, NPAS4 and ARNT2, but does not interact with AhR, HIF1α and ARNT. This interaction is mediated by PAS repeat regions which also form the interface for bHLH PAS dimerisation, and accordingly MAGED1 is not found in complex with bHLH PAS dimers. We show that MAGED1 does not affect bHLH PAS protein levels and cannot be acting as a coactivator of transcriptionally active heterodimers, but rather appears to interact with nascent bHLH PAS proteins in the cytoplasm to enhance their function prior to nuclear import. As a selective regulator, MAGED1 may play an important role in the biology of these specific factors and in general bHLH PAS protein dynamics.

The evolutionary origin of CIPK16: A gene involved in enhanced salt tolerance.

Calcineurin B-like protein interacting protein kinases (CIPKs) are key regulators of pre-transcriptional and post-translational responses to abiotic stress. Arabidopsis thaliana CIPK16 (AtCIPK16) was identified from a forward genetic screen as a gene that mediates lower shoot salt accumulation and improved salinity tolerance in Arabidopsis and transgenic barley. Here, we aimed to gain an understanding of the evolution of AtCIPK16, and orthologues of CIPK16 in other plant species including barley, by conducting a phylogenetic analysis of terrestrial plant species. The resulting protein sequence based phylogenetic trees revealed a single clade that included AtCIPK16 along with two segmentally duplicated CIPKs, AtCIPK5 and AtCIPK25. No monocots had proteins that fell into this clade; instead the most closely related monocot proteins formed a group basal to the entire CIPK16, 5 and 25 clade. We also found that AtCIPK16 contains a core Brassicales specific indel and a putative nuclear localisation signal, which are synapomorphic characters of CIPK16 genes. In addition, we present a model that proposes the evolution of CIPK16, 5 and 25 clade.

MAP4K Family Kinases in Immunity and Inflammation.

MAP kinase kinase kinase kinases (MAP4Ks) belong to the mammalian Ste20-like family of serine/threonine kinases. MAP4Ks including MAP4K1/HPK1, MAP4K2/GCK, MAP4K3/GLK, MAP4K4/HGK, MAP4K5/KHS, and MAP4K6/MINK have been reported to induce JNK activation through activating the MAP3K-MAP2K cascade. The physiological roles of MAP4Ks in immunity and inflammation are largely unknown until recent studies using biochemical approaches and knockout mice. Surprisingly, JNK is not the major target of MAP4Ks in immune cells; MAP4Ks regulate immune responses through novel targets. HPK1 inhibits T-cell receptor (TCR) signaling and B-cell receptor signaling via inducing phosphorylation/ubiquitination of SLP-76 and BLNK, respectively. GLK activates TCR signaling through phosphorylating/activating PKCθ. T-cell-mediated immune responses and Th17-mediated experimental autoimmune diseases are enhanced in HPK1 knockout mice but ameliorated in GLK knockout mice. Consistently, HPK1 levels are decreased in peripheral blood mononuclear cells and T cells from patients with psoriatic arthritis and systemic lupus erythematosus (SLE), respectively. Moreover, GLK levels are increased in T cells from patients with SLE, rheumatoid arthritis, or adult-onset Still's disease; the percentages of GLK-overexpression T cells are correlated with the disease activity. In addition, HGK phosphorylates and induces TRAF2 protein degradation, leading to negative regulation of IL-6 production in resting T cells. Loss of HGK in T cells results in spontaneous systemic inflammation and type 2 diabetes in mice. HGK is also involved in cancer cell migration. To date, the phenotypes of knockout mice for GCK, KHS, and MINK have not been reported; the roles of these three MAP4Ks in immune cell signaling are discussed in this review. Taken together, MAP4K family kinases play diverse roles in immune cell signaling, immune responses, and inflammation.

MdSOS2L1 phosphorylates MdVHA-B1 to modulate malate accumulation in response to salinity in apple.

KEY MESSAGE: Salt-induced phosphorylation of MdVHA-B1 protein was mediated by MdSOS2L1 protein kinase, and thereby increasing malate content in apple. Salinity is an important environmental factor that influences malate accumulation in apple. However, the molecular mechanism by which salinity regulates this process is poorly understood. In this work, we found that MdSOS2L1, a novel AtSOS2-LIKE protein kinase, interacts with V-ATPase subunit MdVHA-B1. Furthermore, MdSOS2L1 directly phosphorylates MdVHA-B1 at Ser(396) site to modulate malate accumulation in response to salt stress. Meanwhile, a series of transgenic analyses in apple calli showed that the MdSOS2L1-MdVHAB1 pathway was involved in the regulation of malate accumulation. Finally, a viral vector-based transformation approach demonstrated that the MdSOS2L1-MdVHAB1 pathway also modulated malate accumulation in apple fruits with or without salt stress. Collectively, our findings provide a new insight into the mechanism by which MdSOS2L1 phosphorylates MdVHA-B1 to modulate malate accumulation in response to salinity in apple.

The phenome analysis of mutant alleles in Leucine-Rich Repeat Receptor-Like Kinase genes in rice reveals new potential targets for stress tolerant cereals.

Plants are constantly exposed to a variety of biotic and abiotic stresses that reduce their fitness and performance. At the molecular level, the perception of extracellular stimuli and the subsequent activation of defense responses require a complex interplay of signaling cascades, in which protein phosphorylation plays a central role. Several studies have shown that some members of the Leucine-Rich Repeat Receptor-Like Kinase (LRR-RLK) family are involved in stress and developmental pathways. We report here a systematic analysis of the role of the members of this gene family by mutant phenotyping in the monocotyledon model plant rice, Oryza sativa. We have then targeted 176 of the ∼320 LRR-RLK genes (55.7%) and genotyped 288 mutant lines. Position of the insertion was confirmed in 128 lines corresponding to 100 LRR-RLK genes (31.6% of the entire family). All mutant lines harboring homozygous insertions have been screened for phenotypes under normal conditions and under various abiotic stresses. Mutant plants have been observed at several stages of growth, from seedlings in Petri dishes to flowering and grain filling under greenhouse conditions. Our results show that 37 of the LRR-RLK rice genes are potential targets for improvement especially in the generation of abiotic stress tolerant cereals.