Identification of the phosphorylation targets of symbiotic receptor-like kinases using a high-throughput multiplexed assay for kinase specificity

A rare non-canonical splice site in Trema orientalis SYMRK does not affect its dual symbiotic functioning in endomycorrhiza and rhizobium nodulation

BACKGROUND

Nitrogen-fixing nodules occur in ten related taxonomic lineages interspersed with lineages of non-nodulating plant species. Nodules result from an endosymbiosis between plants and diazotrophic bacteria; rhizobia in the case of legumes and Parasponia and Frankia in the case of actinorhizal species. Nodulating plants share a conserved set of symbiosis genes, whereas related non-nodulating sister species show pseudogenization of several key nodulation-specific genes. Signalling and cellular mechanisms critical for nodulation have been co-opted from the more ancient plant-fungal arbuscular endomycorrhizal symbiosis. Studies in legumes and actinorhizal plants uncovered a key component in symbiotic signalling, the LRR-type SYMBIOSIS RECEPTOR KINASE (SYMRK). SYMRK is essential for nodulation and arbuscular endomycorrhizal symbiosis. To our surprise, however, despite its arbuscular endomycorrhizal symbiosis capacities, we observed a seemingly critical mutation in a donor splice site in the SYMRK gene of Trema orientalis, the non-nodulating sister species of Parasponia. This led us to investigate the symbiotic functioning of SYMRK in the Trema-Parasponia lineage and to address the question of to what extent a single nucleotide polymorphism in a donor splice site affects the symbiotic functioning of SYMRK.

RESULTS

We show that SYMRK is essential for nodulation and endomycorrhization in Parasponia andersonii. Subsequently, it is revealed that the 5'-intron donor splice site of SYMRK intron 12 is variable and, in most dicotyledon species, doesn't contain the canonical dinucleotide 'GT' signature but the much less common motif 'GC'. Strikingly, in T. orientalis, this motif is converted into a rare non-canonical 5'-intron donor splice site 'GA'. This SYMRK allele, however, is fully functional and spreads in the T. orientalis population of Malaysian Borneo. A further investigation into the occurrence of the non-canonical GA-AG splice sites confirmed that these are extremely rare.

CONCLUSION

SYMRK functioning is highly conserved in legumes, actinorhizal plants, and Parasponia. The gene possesses a non-common 5'-intron GC donor splice site in intron 12, which is converted into a GA in T. orientalis accessions of Malaysian Borneo. The discovery of this functional GA-AG splice site in SYMRK highlights a gap in our understanding of splice donor sites.

Brassinosteroids modulate autophagy through phosphorylation of RAPTOR1B by the GSK3-like kinase BIN2 in Arabidopsis

Macroautophagy/autophagy is a conserved recycling process that maintains cellular homeostasis during environmental stress. Autophagy is negatively regulated by TOR (target of rapamycin), a nutrient-regulated protein kinase that in plants is activated by several phytohormones, leading to increased growth. However, the detailed molecular mechanisms by which TOR integrates autophagy and hormone signaling are poorly understood. Here, we show that TOR modulates brassinosteroid (BR)-regulated plant growth and stress-response pathways. Active TOR was required for full BR-mediated growth in Arabidopsis thaliana. Autophagy was constitutively up-regulated upon blocking BR biosynthesis or signaling, and down-regulated by increasing the activity of the BR pathway. BIN2 (brassinosteroid-insensitive 2) kinase, a GSK3-like kinase functioning as a negative regulator in BR signaling, directly phosphorylated RAPTOR1B (regulatory-associated protein of TOR 1B), a substrate-recruiting subunit in the TOR complex, at a conserved serine residue within a typical BIN2 phosphorylation motif. Mutation of RAPTOR1B serine 916 to alanine, to block phosphorylation by BIN2, repressed autophagy and increased phosphorylation of the TOR substrate ATG13a (autophagy-related protein 13a). By contrast, this mutation had only a limited effect on growth. We present a model in which RAPTOR1B is phosphorylated and inhibited by BIN2 when BRs are absent, activating the autophagy pathway. When BRs signal and inhibit BIN2, RAPTOR1B is thus less inhibited by BIN2 phosphorylation. This leads to increased TOR activity and ATG13a phosphorylation, and decreased autophagy activity. Our studies define a new mechanism by which coordination between BR and TOR signaling pathways helps to maintain the balance between plant growth and stress responses.

Integration of multi-omics data reveals interplay between brassinosteroid and Target of Rapamycin Complex signaling in Arabidopsis

Brassinosteroids (BRs) and Target of Rapamycin Complex (TORC) are two major actors coordinating plant growth and stress responses. Brassinosteroids function through a signaling pathway to extensively regulate gene expression and TORC is known to regulate translation and autophagy. Recent studies have revealed connections between these two pathways, but a system-wide view of their interplay is still missing. We quantified the level of 23 975 transcripts, 11 183 proteins, and 27 887 phosphorylation sites in wild-type Arabidopsis thaliana and in mutants with altered levels of either BRASSINOSTEROID INSENSITIVE 2 (BIN2) or REGULATORY ASSOCIATED PROTEIN OF TOR 1B (RAPTOR1B), two key players in BR and TORC signaling, respectively. We found that perturbation of BIN2 or RAPTOR1B levels affects a common set of gene-products involved in growth and stress responses. Furthermore, we used the multi-omic data to reconstruct an integrated signaling network. We screened 41 candidate genes identified from the reconstructed network and found that loss of function mutants of many of these proteins led to an altered BR response and/or modulated autophagy activity. Altogether, these results establish a predictive network that defines different layers of molecular interactions between BR- or TORC-regulated growth and autophagy.

NIPK, a protein pseudokinase that interacts with the C subunit of the transcription factor NF-Y, is involved in rhizobial infection and nodule organogenesis

Heterotrimeric Nuclear Factor Y (NF-Y) transcription factors are key regulators of the symbiotic program that controls rhizobial infection and nodule organogenesis. Using a yeast two-hybrid screening, we identified a putative protein kinase of Phaseolus vulgaris that interacts with the C subunit of the NF-Y complex. Physical interaction between NF-YC1 Interacting Protein Kinase (NIPK) and NF-YC1 occurs in the cytoplasm and the plasma membrane. Only one of the three canonical amino acids predicted to be required for catalytic activity is conserved in NIPK and its putative homologs from lycophytes to angiosperms, indicating that NIPK is an evolutionary conserved pseudokinase. Post-transcriptional silencing on NIPK affected infection and nodule organogenesis, suggesting NIPK is a positive regulator of the NF-Y transcriptional complex. In addition, NIPK is required for activation of cell cycle genes and early symbiotic genes in response to rhizobia, including NF-YA1 and NF-YC1. However, strain preference in co-inoculation experiments was not affected by NIPK silencing, suggesting that some functions of the NF-Y complex are independent of NIPK. Our work adds a new component associated with the NF-Y transcriptional regulators in the context of nitrogen-fixing symbiosis.

Advances in functional proteomics to study plant-pathogen interactions

Pathogen infection triggers complex signaling networks in plant cells that ultimately result in either susceptibility or resistance. We have made substantial progress in dissecting many of these signaling events, and it is becoming clear that changes in proteome composition and protein activity are major drivers of plant-microbe interactions. Here, we highlight different approaches to analyze the functional proteomes of hosts and pathogens and discuss how they have been used to further our understanding of plant disease. Global proteome profiling can quantify the dynamics of proteins, posttranslational modifications, and biological pathways that contribute to immune-related outcomes. In addition, emerging techniques such as enzyme activity-based profiling, proximity labeling, and kinase-substrate profiling are being used to dissect biochemical events that operate during infection. Finally, we discuss how these functional approaches can be integrated with other profiling data to gain a mechanistic, systems-level view of plant and pathogen signaling.

Tomato SD1, encoding a kinase-interacting protein, is a major locus controlling stem development

Stems serve as key determinants of plant development by connecting and supporting parts of the plant body, transporting nutrients important for long-distance communication that affect crop yield, and producing new organs. Nonetheless, studies on the regulation of stem development in crops are rather limited. Here, we found a significant correlation (P<0.001) between stem diameter (SD) and fruit size in tomato (Solanum lycopersicum). We performed a genome-wide association study and identified a novel quantitative trait locus (QTL), SDR9 (stem diameter regulator on CHROMOSOME 9), that co-localized with a gene encoding a kinase-interacting family protein (KIP), which is the most likely candidate gene related to SD (hereafter referred to as SD1). Overexpression of SD1 in thin-stem accessions resulted in increased SD. In contrast, suppressed expression of SD1 in thick-stem accessions using RNA interference exhibited the opposite effect. Further microscopic analyses showed that SD1 affected the stem diameter by controlling the size and number of secondary phloem cells. An 11-bp indel in the promoter region of SD1 that disrupts a gibberellin-responsive cis-element was linked to SD. Expression analysis revealed that SD1 was mainly expressed at the cambium of the stem and positively regulates stem development. Evolutionary analysis revealed that the thick-stem allele of SD1 was selected during the recent process of tomato improvement. Our results provide novel genetic and molecular insight into natural variation of SD in tomato and may accelerate the breeding of high yield tomato.

Identification of nodulation-related genes in Medicago truncatula using genome-wide association studies and co-expression networks

Genome-wide association studies (GWAS) have proven to be a valuable approach for identifying genetic intervals associated with phenotypic variation in Medicago truncatula. These intervals can vary in size, depending on the historical local recombination. Typically, significant intervals span numerous gene models, limiting the ability to resolve high-confidence candidate genes underlying the trait of interest. Additional genomic data, including gene co-expression networks, can be combined with the genetic mapping information to successfully identify candidate genes. Co-expression network analysis provides information about the functional relationships of each gene through its similarity of expression patterns to other well-defined clusters of genes. In this study, we integrated data from GWAS and co-expression networks to pinpoint candidate genes that may be associated with nodule-related phenotypes in M. truncatula. We further investigated a subset of these genes and confirmed that several had existing evidence linking them nodulation, including MEDTR2G101090 (PEN3-like), a previously validated gene associated with nodule number.