Information-incorporated gene network construction with FDR control.

MOTIVATION: Large-scale gene expression studies allow gene network construction to uncover associations among genes. To study direct associations among genes, partial correlation-based networks are preferred over marginal correlations. However, FDR control for partial correlation-based network construction is not well-studied. In addition, currently available partial correlation-based methods cannot take existing biological knowledge to help network construction while controlling FDR. RESULTS: In this paper, we propose a method called Partial Correlation Graph with Information Incorporation (PCGII). PCGII estimates partial correlations between each pair of genes by regularized node-wise regression that can incorporate prior knowledge while controlling the effects of all other genes. It handles high-dimensional data where the number of genes can be much larger than the sample size and controls FDR at the same time. We compare PCGII with several existing approaches through extensive simulation studies and demonstrate that PCGII has better FDR control and higher power. We apply PCGII to a plant gene expression dataset where it recovers confirmed regulatory relationships and a hub node, as well as several direct associations that shed light on potential functional relationships in the system. We also introduce a method to supplement observed data with a pseudogene to apply PCGII when no prior information is available, which also allows checking FDR control and power for real data analysis. AVAILABILITY AND IMPLEMENTATION: R package is freely available for download at https://cran.r-project.org/package=PCGII.

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.

FERONIA functions through Target of Rapamycin (TOR) to negatively regulate autophagy.

FERONIA (FER) receptor kinase plays versatile roles in plant growth and development, biotic and abiotic stress responses, and reproduction. Autophagy is a conserved cellular recycling process that is critical for balancing plant growth and stress responses. Target of Rapamycin (TOR) has been shown to be a master regulator of autophagy. Our previous multi-omics analysis with loss-of-function fer-4 mutant implicated that FER functions in the autophagy pathway. We further demonstrated here that the fer-4 mutant displayed constitutive autophagy, and FER is required for TOR kinase activity measured by S6K1 phosphorylation and by root growth inhibition assay to TOR kinase inhibitor AZD8055. Taken together, our study provides a previously unknown mechanism by which FER functions through TOR to negatively regulate autophagy.

The deubiquitinating enzymes UBP12 and UBP13 positively regulate recovery after carbon starvation by modulating BES1 stability in Arabidopsis thaliana.

BRI1-EMS-SUPPRESSOR1 (BES1), a core transcription factor in the brassinosteroid (BR) signaling pathway, primarily regulates plant growth and development by influencing BR-regulated gene expression. Several E3 ubiquitin (Ub) ligases regulate BES1 stability, but little is known about BES1 deubiquitination, which antagonizes E3 ligase-mediated ubiquitination to maintain BES1 homeostasis. Here, we report that two Arabidopsis thaliana deubiquitinating enzymes, Ub-SPECIFIC PROTEASE (UBP) 12 and UBP13, interact with BES1. UBP12 and UBP13 removed Ub from polyubiquitinated BES1 to stabilize both phosphorylated and dephosphorylated forms of BES1. A double mutant, ubp12-2w ubp13-3, lacking UBP12 and UBP13 function showed both BR-deficient and BR-insensitive phenotypes, whereas transgenic plants overexpressing UBP12 or UBP13 exhibited an increased BR response. Expression of UBP12 and UPB13 was induced during recovery after carbon starvation, which led to BES1 accumulation and quick recovery of stressed plants. Our work thus establishes a mechanism by which UBP12 and UBP13 regulate BES1 protein abundance to enhance BR-regulated growth during recovery after carbon starvation.

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.

Integrated omics reveal novel functions and underlying mechanisms of the receptor kinase FERONIA in Arabidopsis thaliana.

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.

A Heterogeneous RISC-V Processor for Efficient DNN Application in Smart Sensing System.

Extracting features from sensing data on edge devices is a challenging application for which deep neural networks (DNN) have shown promising results. Unfortunately, the general micro-controller-class processors which are widely used in sensing system fail to achieve real-time inference. Accelerating the compute-intensive DNN inference is, therefore, of utmost importance. As the physical limitation of sensing devices, the design of processor needs to meet the balanced performance metrics, including low power consumption, low latency, and flexible configuration. In this paper, we proposed a lightweight pipeline integrated deep learning architecture, which is compatible with open-source RISC-V instructions. The dataflow of DNN is organized by the very long instruction word (VLIW) pipeline. It combines with the proposed special intelligent enhanced instructions and the single instruction multiple data (SIMD) parallel processing unit. Experimental results show that total power consumption is about 411 mw and the power efficiency is about 320.7 GOPS/W.

Integrated omics networks reveal the temporal signaling events of brassinosteroid response in Arabidopsis.

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.

The F-box E3 ubiquitin ligase BAF1 mediates the degradation of the brassinosteroid-activated transcription factor BES1 through selective autophagy in Arabidopsis.

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.

A BIN2-GLK1 Signaling Module Integrates Brassinosteroid and Light Signaling to Repress Chloroplast Development in the Dark.

Arabidopsis GLYCOGEN SYNTHASE KINASE 3 (GSK3)-like kinases play various roles in plant development, including chloroplast development, but the underlying molecular mechanism is not well defined. Here, we demonstrate that transcription factors GLK1 and GLK2 interact with and are phosphorylated by the BRASSINOSTEROID insensitive2 (BIN2). The loss-of-function mutant of BIN2 and its homologs, bin2-3 bil1 bil2, displays abnormal chloroplast development, whereas the gain-of-function mutant, bin2-1, exhibits insensitivity to BR-induced de-greening and reduced numbers of thylakoids per granum, suggesting that BIN2 positively regulates chloroplast development. Furthermore, BIN2 phosphorylates GLK1 at T175, T238, T248, and T256, and mutations of these phosphorylation sites alter GLK1 protein stability and DNA binding and impair plant responses to BRs/darkness. On the other hand, BRs and darkness repress the BIN2-GLK module to enhance BR/dark-mediated de-greening and impair the formation of the photosynthetic apparatus. Our results thus provide a mechanism by which BRs modulate photomorphogenesis and chloroplast development.

Measuring Protein Half-life in Arabidopsis thaliana.

Post-translational modifications play important roles in controlling protein function and can lead to altered protein stability. Protein stability can be determined after treatment with the protein synthesis inhibitor Cycloheximide. Cycloheximide is a translational inhibitor that inhibits protein synthesis via cytoplasmic ribosomes. Here we describe how to measure the stability of MYC2 in the context of regulation by FERONIA receptor kinase. First, we describe how to measure MYC2 stability in wild-type and feronia mutant; then we describe similar assays in transgenic plants expressing MYC2-FLAG and MYC2A12-FLAG (12 FERONIA phosphorylation sites are mutated to Alanine and the mutant protein is stabilized). MYC2 can be induced by mechanical touch, which can be a confounding factor in protein level measurement. In this protocol, we take that into consideration and try to achieve more accurate measurement.

FERONIA mutation induces high levels of chloroplast-localized Arabidopsides which are involved in root growth.

The FERONIA (FER) signaling pathway is known to have diverse roles in Arabidopsis thaliana, such as growth, reproduction, and defense, but how this receptor kinase is involved in various biological processes is not well established. In this work, we applied multiple mass spectrometry techniques to identify metabolites involved in the FER signaling pathway and to understand their biological roles. A direct infusion Fourier transform ion cyclotron resonance (FT-ICR)-MS approach was used for initial screening of wild-type and feronia (fer) mutant plant extracts, and Arabidopsides were found to be significantly enriched in the mutant. As Arabidopsides are known to be induced by wounding, further experiments on wounded and non-wounded leaf samples were carried out to investigate these oxylipins as well as related phytohormones using a quadrupole-time-of-flight (Q-TOF) MS by direct injection and LC-MS/MS. In a root growth bioassay with Arabidopside A isolated from fer mutants, the wild-type showed significant root growth inhibition compared with the fer mutant. Our results therefore implicated Arabidopsides, and Arabidopside A specifically, in FER functions and/or signaling. Finally, matrix-assisted laser desorption/ionization MS imaging (MALDI-MSI) was used to visualize the localization of Arabidopsides, and we confirmed that Arabidopsides are highly abundant at wounding sites in both wild-type and fer mutant leaves. More significantly, five micron high-spatial resolution MALDI-MSI revealed that Arabidopsides are localized to the chloroplasts where many stress signaling molecules are made.

FERONIA Receptor Kinase Contributes to Plant Immunity by Suppressing Jasmonic Acid Signaling in Arabidopsis thaliana.

Bacterial pathogens use effectors and phytotoxins to facilitate infection of host plants. Coronatine (COR) is one of the phytotoxins produced in bacterial pathogens, such as Pseudomonas syringae pv. tomato DC3000 (pst DC3000). COR structurally and functionally mimics the active form of the plant hormone jasmonic acid (JA), JA-isoleucine (JA-Ile), and can hijack the host JA-signaling pathway to achieve host disease susceptibility [1]. COR utilizes the transcription factor MYC2, a master regulator of JA signaling, to activate NAC transcription factors, which functions to inhibit accumulation of salicylic acid (SA) and thus compromise host immunity [2]. It has been demonstrated that SA can antagonize JA signaling through NONEXPRESSOR of PATHOGENESIS-RELATED GENE1 (NPR1) [3] and downstream transcription factors TGAs [4] and WRKYs [5, 6]. However, the detailed mechanism by which host plants counteract COR-mediated susceptibility is largely unknown. Here, we show that the receptor kinase FERONIA (FER) functions to inhibit JA and COR signaling by phosphorylating and destabilizing MYC2, thereby positively regulating immunity. Conversely, the peptide ligand RALF23 acts through FER to stabilize MYC2 and elevate JA signaling, negatively contributing to plant immunity. Our results establish the RALF23-FER-MYC2 signaling module and provide a previously unknown mechanism by which host plants utilize FER signaling to counteract COR-mediated host disease susceptibility.

RD26 mediates crosstalk between drought and brassinosteroid signalling pathways.

Brassinosteroids (BRs) regulate plant growth and stress responses via the BES1/BZR1 family of transcription factors, which regulate the expression of thousands of downstream genes. BRs are involved in the response to drought, however the mechanistic understanding of interactions between BR signalling and drought response remains to be established. Here we show that transcription factor RD26 mediates crosstalk between drought and BR signalling. When overexpressed, BES1 target gene RD26 can inhibit BR-regulated growth. Global gene expression studies suggest that RD26 can act antagonistically to BR to regulate the expression of a subset of BES1-regulated genes, thereby inhibiting BR function. We show that RD26 can interact with BES1 protein and antagonize BES1 transcriptional activity on BR-regulated genes and that BR signalling can also repress expression of RD26 and its homologues and inhibit drought responses. Our results thus reveal a mechanism coordinating plant growth and drought tolerance.

Identification of Brassinosteroid Target Genes by Chromatin Immunoprecipitation Followed by High-Throughput Sequencing (ChIP-seq) and RNA-Sequencing.

Brassinosteroids (BRs) play important roles in many growth and developmental processes. BRs signal to regulate BR-INSENSITIVE1-ETHYL METHANESULFONATE-SUPPRESSOR1 (BES1) and BRASSINAZOLE-RESISTANT1 (BZR1) transcription factors (TFs), which, in turn, regulate several hundreds of transcription factors (termed BES1/BZR1-targeted TFs or BTFs) and thousands of genes to mediate various BR responses. Chromatin Immunoprecipitation followed by high-throughput sequencing (ChIP-seq) with BES1/BZR1 and BTFs is an important approach to identify BR target genes. In combination with RNA-sequencing experiments, these genomic methods have become powerful tools to detect BR target genes and reveal transcriptional networks underlying BR-regulated processes.

Role of brassinosteroid signaling in modulating Tobacco mosaic virus resistance in Nicotiana benthamiana.

Plant steroid hormones, brassinosteroids (BRs), play essential roles in plant growth, development and stress responses. However, mechanisms by which BRs interfere with plant resistance to virus remain largely unclear. In this study, we used pharmacological and genetic approaches in combination with infection experiments to investigate the role of BRs in plant defense against Tobacco Mosaic Virus (TMV) in Nicotiana benthamiana. Exogenous applied BRs enhanced plant resistance to virus infection, while application of Bikinin (inhibitor of glycogen synthase kinase-3), which activated BR signaling, increased virus susceptibility. Silencing of NbBRI1 and NbBSK1 blocked BR-induced TMV resistance, and silencing of NbBES1/BZR1 blocked Bikinin-reduced TMV resistance. Silencing of NbMEK2, NbSIPK and NbRBOHB all compromised BR-induced virus resistance and defense-associated genes expression. Furthermore, we found MEK2-SIPK cascade activated while BES1/BZR1 inhibited RBOHB-dependent ROS production, defense gene expression and virus resistance induced by BRs. Thus, our results revealed BR signaling had two opposite effects on viral defense response. On the one hand, BRs enhanced virus resistance through MEK2-SIPK cascade and RBOHB-dependent ROS burst. On the other hand, BES1/BZR1 inhibited RBOHB-dependent ROS production and acted as an important mediator of the trade-off between growth and immunity in BR signaling.

Light intensity affects chlorophyll synthesis during greening process by metabolite signal from mitochondrial alternative oxidase in Arabidopsis.

Although mitochondrial alternative oxidase (AOX) has been proposed to play essential roles in high light stress tolerance, the effects of AOX on chlorophyll synthesis are unclear. Previous studies indicated that during greening, chlorophyll accumulation was largely delayed in plants whose mitochondrial cyanide-resistant respiration was inhibited by knocking out nuclear encoded AOX gene. Here, we showed that this delay of chlorophyll accumulation was more significant under high light condition. Inhibition of cyanide-resistant respiration was also accompanied by the increase of plastid NADPH/NADP(+) ratio, especially under high light treatment which subsequently blocked the import of multiple plastidial proteins, such as some components of the photosynthetic electron transport chain, the Calvin-Benson cycle enzymes and malate/oxaloacetate shuttle components. Overexpression of AOX1a rescued the aox1a mutant phenotype, including the chlorophyll accumulation during greening and plastidial protein import. It thus suggests that light intensity affects chlorophyll synthesis during greening process by a metabolic signal, the AOX-derived plastidial NADPH/NADP(+) ratio change. Further, our results thus revealed a molecular mechanism of chloroplast-mitochondria interactions.

Histone lysine methyltransferase SDG8 is involved in brassinosteroid-regulated gene expression in Arabidopsis thaliana.

The plant steroid hormones, brassinosteroids (BRs), play important roles in plant growth, development, and responses to environmental stresses. BRs signal through receptors localized to the plasma membrane and other signaling components to regulate the BES1/BZR1 family of transcription factors, which modulates the expression of thousands of genes. How BES1/BZR1 and their interacting proteins function to regulate the large number of genes are not completely understood. Here we report that histone lysine methyltransferase SDG8, implicated in histone 3 lysine 36 di- and trimethylation (H3K36me2 and me3), is involved in BR-regulated gene expression. BES1 interacts with SDG8, directly or indirectly through IWS1, a transcription elongation factor involved in BR-regulated gene expression. The knockout mutant sdg8 displays a reduced growth phenotype with compromised BR responses. Global gene expression studies demonstrated that, while BR regulates about 5000 genes in wild-type plants, the hormone regulates fewer than 700 genes in sdg8 mutant. In addition, more than half of BR-regulated genes are differentially affected in sdg8 mutant. A Chromatin Immunoprecipitation (ChIP) experiment showed that H3K36me3 is reduced in BR-regulated genes in the sdg8 mutant. Based on these results, we propose that SDG8 plays an essential role in mediating BR-regulated gene expression. Our results thus reveal a major mechanism by which histone modifications dictate hormonal regulation of gene expression.

Transcription factors involved in brassinosteroid repressed gene expression and their regulation by BIN2 kinase.

Brassinosteroids (BRs) play important roles in plant growth, development and responses to environmental cues. BRs signal through plasma membrane-localized receptors to control BES1/BZR1 family transcription factors, which mediate the expression of thousands of genes. BRs activate and repress approximately equal numbers of genes. BES1/BZR1 interact with other transcription factors, histone-modifying enzymes, and transcription elongation factors to activate BR-induced genes. However, the mechanisms by which BES1/BZR1 mediate the BR-repressed gene expression are not well understood. Recent studies revealed that 2 BR-repressed transcription factors, MYBL2 and HAT1, cooperate with BES1 to downregulate BR-repressed genes expression. Moreover, BIN2 kinase, a well-established negative regulator in the BR signaling pathway, phosphorylates MYBL2 and HAT1. While BIN2 phosphorylates and destabilizes BES1/BZR1, BIN2 phosphorylated MYBL2 and HAT1 appear to be stabilized. These results not only extended our understanding of BR-repressed gene expression, but also revealed multiple inputs of BR signaling into BR transcriptional networks.