Brassinosteroids (BRs) are plant steroid hormones that regulate cell division and stress response. Here we use a systems biology approach to integrate multi-omic datasets and unravel the molecular signaling events of BR response in Arabidopsis. We profile the levels of 26,669 transcripts, 9,533 protein groups, and 26,617 phosphorylation sites from Arabidopsis seedlings treated with brassinolide (BL) for six different lengths of time. We then construct a network inference pipeline called Spatiotemporal Clustering and Inference of Omics Networks (SC-ION) to integrate these data. We use our network predictions to identify putative phosphorylation sites on BES1 and experimentally validate their importance. Additionally, we identify BRONTOSAURUS (BRON) as a transcription factor that regulates cell division, and we show that BRON expression is modulated by BR-responsive kinases and transcription factors. This work demonstrates the power of integrative network analysis applied to multi-omic data and provides fundamental insights into the molecular signaling events occurring during BR response.

Brassinosteroids are plant steroid hormones that regulate diverse processes, such as cell division and cell elongation, through gene regulatory networks that vary in space and time. By using time series single-cell RNA sequencing to profile brassinosteroid-responsive gene expression specific to different cell types and developmental stages of the Arabidopsis root, we identified the elongating cortex as a site where brassinosteroids trigger a shift from proliferation to elongation associated with increased expression of cell wall–related genes. Our analysis revealed HOMEOBOX FROM ARABIDOPSIS THALIANA 7 ( HAT7 ) and GT-2-LIKE 1 ( GTL1 ) as brassinosteroid-responsive transcription factors that regulate cortex cell elongation. These results establish the cortex as a site of brassinosteroid-mediated growth and unveil a brassinosteroid signaling network regulating the transition from proliferation to elongation, which illuminates aspects of spatiotemporal hormone responses.

The receptor kinase FERONIA (FER) is a versatile regulator of plant growth and development, biotic and abiotic stress responses, and reproduction. To gain new insights into the molecular interplay of these processes and to identify new FER functions, we carried out quantitative transcriptome, proteome, and phosphoproteome profiling of Arabidopsis (Arabidopsis thaliana) wild-type and fer-4 loss-of-function mutant plants. Gene ontology terms for phytohormone signaling, abiotic stress, and biotic stress were significantly enriched among differentially expressed transcripts, differentially abundant proteins, and/or misphosphorylated proteins, in agreement with the known roles for FER in these processes. Analysis of multiomics data and subsequent experimental evidence revealed previously unknown functions for FER in endoplasmic reticulum (ER) body formation and glucosinolate biosynthesis. FER functions through the transcription factor NAI1 to mediate ER body formation. FER also negatively regulates indole glucosinolate biosynthesis, partially through NAI1. Furthermore, we found that a group of abscisic acid (ABA)-induced transcription factors is hypophosphorylated in the fer-4 mutant and demonstrated that FER acts through the transcription factor ABA INSENSITIVE5 (ABI5) to negatively regulate the ABA response during cotyledon greening. Our integrated omics study therefore reveals novel functions for FER and provides new insights into the underlying mechanisms of FER function.

Brassinosteroids (BRs) regulate plant growth, development, and stress responses by activating the core transcription factor BRI1-EMS-SUPPRESSOR1 (BES1), whose degradation occurs through the proteasome and autophagy pathways. The E3 ubiquitin ligase(s) that modify BES1 for autophagy-mediated degradation remain to be fully defined. Here, we identified an F-box family E3 ubiquitin ligase named BES1-ASSOCIATED F-BOX1 (BAF1) in Arabidopsis thaliana. BAF1 interacts with BES1 and mediates its ubiquitination and degradation. Our genetic data demonstrated that BAF1 inhibits BR signaling in a BES1-dependent manner. Moreover, BAF1 targets BES1 for autophagic degradation in a selective manner. BAF1-triggered selective autophagy of BES1 depends on the ubiquitin binding receptor DOMINANT SUPPRESSOR OF KAR2 (DSK2). Sucrose starvation-induced selective autophagy of BES1, but not bulk autophagy, was significantly compromised in baf1 mutant and BAF1-ΔF (BAF1 F-box decoy) overexpression plants, but clearly increased by BAF1 overexpression. The baf1 and BAF1-ΔF overexpression plants had increased BR-regulated growth but were sensitive to long-term sucrose starvation, while BAF1 overexpression plants had decreased BR-regulated growth but were highly tolerant of sucrose starvation. Our results not only established BAF1 as an E3 ubiquitin ligase that targets BES1 for degradation through selective autophagy pathway, but also revealed a mechanism for plants to reduce growth during sucrose starvation.

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 BRASSI-NOSTEROID 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 geneproducts 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.

Autophagy is a conserved recycling process with important functions in plant growth, development, and stress responses. Phytohormones also play key roles in the regulation of some of the same processes. Increasing evidence indicates that a close relationship exists between autophagy and phytohormone signaling pathways, and the mechanisms of interaction between these pathways have begun to be revealed. Here, we review recent advances in our understanding of how autophagy regulates hormone signaling and, conversely, how hormones regulate the activity of autophagy, both in plant growth and development and in environmental stress responses. We highlight in particular recent mechanistic insights into the coordination between autophagy and signaling events controlled by the stress hormone abscisic acid and by the growth hormones brassinosteroid and cytokinin and briefly discuss potential connections between autophagy and other phytohormones.

Plant architecture is one of the most important factors that determines crop yield potential and productivity. In apple (Malus domestica), genetic improvement of tree architecture has been challenging due to a long juvenile phase and growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is one of the four paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis (Arabidopsis thaliana). The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T > C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that co-localizes with Region III, one of the five conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Over-expressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response, but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.

Plant architecture is one of the most important factors that determines crop yield potential and productivity. In apple (Malus), genetic improvement of tree architecture has been challenging due to a long juvenile phase and their growth as complex trees composed of a distinct scion and a rootstock. To better understand the genetic control of apple tree architecture, the dominant weeping growth phenotype was investigated. We report the identification of MdLAZY1A (MD13G1122400) as the genetic determinant underpinning the Weeping (W) locus that largely controls weeping growth in Malus. MdLAZY1A is one of the four paralogs in apple that are most closely related to AtLAZY1 involved in gravitropism in Arabidopsis. The weeping allele (MdLAZY1A-W) contains a single nucleotide mutation c.584T>C that leads to a leucine to proline (L195P) substitution within a predicted transmembrane domain that co-localizes with Region III, one of the five conserved regions in LAZY1-like proteins. Subcellular localization revealed that MdLAZY1A localizes to the plasma membrane and nucleus in plant cells. Over-expressing the weeping allele in apple cultivar Royal Gala (RG) with standard growth habit impaired its gravitropic response and altered the growth to weeping-like. Suppressing the standard allele (MdLAZY1A-S) by RNA interference (RNAi) in RG similarly changed the branch growth direction to downward. Overall, the L195P mutation in MdLAZY1A is genetically causal for weeping growth, underscoring not only the crucial roles of residue L195 and Region III in MdLAZY1A-mediated gravitropic response, but also a potential DNA base editing target for tree architecture improvement in Malus and other crops.