A Systems Approach Uncovers Restrictions for Signal Interactions Regulating Genome-wide Responses to Nutritional Cues in Arabidopsis

As sessile organisms, plants must cope with multiple and combined variations of signals in their environment. However, very few reports have studied the genome-wide effects of systematic signal combinations on gene expression. Here, we evaluate a high level of signal integration, by modeling genome-wide expression patterns under a factorial combination of carbon (C), light (L), and nitrogen (N) as binary factors in two organs (O), roots and leaves. Signal management is different between C, N, and L and in shoots and roots. For example, L is the major factor controlling gene expression in leaves. However, in roots there is no obvious prominent signal, and signal interaction is stronger. The major signal interaction events detected genome wide in Arabidopsis roots are deciphered and summarized in a comprehensive conceptual model. Surprisingly, global analysis of gene expression in response to C, N, L, and O revealed that the number of genes controlled by a signal is proportional to the magnitude of the gene expression changes elicited by the signal. These results uncovered a strong constraining structure in plant cell signaling pathways, which prompted us to propose the existence of a “code” of signal integration

Modeling the global effect of the basic-leucine zipper transcription factor 1 (bZIP1) on nitrogen and light regulation in Arabidopsis

Background: Nitrogen and light are two major regulators of plant metabolism and development. While genes involved in the control of each of these signals have begun to be identified, regulators that integrate gene responses to nitrogen and light signals have yet to be determined. Here, we evaluate the role of bZIP1, a transcription factor involved in light and nitrogen sensing, by exposing wild-type (WT) and bZIP1 T-DNA null mutant plants to a combinatorial space of nitrogen (N) and light (L) treatment conditions and performing transcriptome analysis. We use ANOVA analysis combined with clustering and Boolean modeling, to evaluate the role of bZIP1 in mediating L and N signaling genome-wide. Results: This transcriptome analysis demonstrates that a mutation in the bZIP1 gene can alter the L and/or N-regulation of several gene clusters. More surprisingly, the bZIP1 mutation can also trigger N and/or L regulation of genes that are not normally controlled by these signals in WT plants. This analysis also reveals that bZIP1 can, to a large extent, invert gene regulation (e. g., several genes induced by N in WT plants are repressed by N in the bZIP1 mutant). Conclusion: These findings demonstrate that the bZIP1 mutation triggers a genome-wide de-regulation in response to L and/or N signals that range from i) a reduction of the L signal effect, to ii) unlocking gene regulation in response to L and N combinations. This systems biology approach demonstrates that bZIP1 tunes L and N signaling relationships genome-wide, and can suppress regulatory mechanisms hypothesized to be needed at different developmental stages and/or environmental conditions

Cross-Species Network Analysis Uncovers Conserved Nitrogen-Regulated Network Modules in Rice

In this study, we used a cross-species network approach to uncover nitrogen-regulated network modules conserved across a model and a crop species. By translating gene “network knowledge” from the data-rich model Arabidopsis (Arabidopsis thaliana) to a crop (Oryza sativa), we identified evolutionarily conserved N-regulatory modules as targets for translational studies to improve N-use efficiency in transgenic plants. To uncover such conserved N-regulatory network modules, we first generated a N-regulatory network based solely on rice (O. sativa) transcriptome and gene interaction data. Next, we enhanced the “network knowledge” in the rice N-regulatory network using transcriptome and gene interaction data from Arabidopsis and new data from Arabidopsis and rice plants exposed to the same N-treatment conditions. This cross-species network analysis uncovered a set of N-regulated transcription factors (TFs) predicted to target the same genes and network modules in both species. Supernode analysis of the TFs and their targets in these conserved network modules uncovered genes directly related to nitrogen use (e.g. N-assimilation) and to other shared biological processes indirectly related to nitrogen. This cross-species network approach was validated with members of two TF families in the supernode network, bZIP-TGA and HRS1/HHO family, have recently been experimentally validated to mediate the N-response in Arabidopsis.Fil: Obertello, Mariana. University of New York; Estados Unidos. Consejo Nacional de Investigaciones Científicas y Técnicas. Instituto de Investigaciones en Ingeniería Genética y Biología Molecular ; ArgentinaFil: Shrivastava, Stuti. University of New York; Estados UnidosFil: Katari, Manpreet S.. University of New York; Estados UnidosFil: Coruzzi, Gloria M.. University of New York; Estados Unido

Hit-and-run transcriptional control by bZIP1 mediates rapid nutrient signaling in Arabidopsis

The dynamic nature of gene regulatory networks allows cells to rapidly respond to environmental change. However, the underlying temporal connections are missed, even in kinetic studies, as transcription factor (TF) binding within at least one time point is required to identify primary targets. The TF-regulated but unbound genes are dismissed as secondary targets. Instead, we report that these genes comprise transient TF-target interactions most relevant to rapid signal transduction. We temporally perturbed a master TF (Basic Leucine Zipper 1, bZIP1) and the nitrogen (N) signal it transduces and integrated TF regulation and binding data from the same cell samples. Our enabling approach could identify primary TF targets based solely on gene regulation, in the absence of TF binding. We uncovered three classes of primary TF targets: (i) poised (TF-bound but not TF-regulated), (ii) stable (TF-bound and TF-regulated), and (iii) transient (TF-regulated but not TF-bound), the largest class. Unexpectedly, the transient bZIP1 targets are uniquely relevant to rapid N signaling in planta, enriched in dynamic N-responsive genes, and regulated by TF and N signal interactions. These transient targets include early N responders nitrate transporter 2.1 and NIN-like protein 3, bound by bZIP1 at 1-5 min, but not at later time points following TF perturbation. Moreover, promoters of these transient targets are uniquely enriched with cis-regulatory motifs coinherited with bZIP1 binding sites, suggesting a recruitment role for bZIP1. This transient mode of TF action supports a classic, but forgotten, "hit-and-run" transcription model, which enables a "catalyst TF" to activate a large set of targets within minutes of signal perturbation

Correlation Network Analysis reveals a sequential reorganization of metabolic and transcriptional states during germination and gene-metabolite relationships in developing seedlings of Arabidopsis

<p>Abstract</p> <p>Background</p> <p>Holistic profiling and systems biology studies of nutrient availability are providing more and more insight into the mechanisms by which gene expression responds to diverse nutrients and metabolites. Less is known about the mechanisms by which gene expression is affected by endogenous metabolites, which can change dramatically during development. Multivariate statistics and correlation network analysis approaches were applied to non-targeted profiling data to investigate transcriptional and metabolic states and to identify metabolites potentially influencing gene expression during the heterotrophic to autotrophic transition of seedling establishment.</p> <p>Results</p> <p>Microarray-based transcript profiles were obtained from extracts of Arabidopsis seeds or seedlings harvested from imbibition to eight days-old. ¹H-NMR metabolite profiles were obtained for corresponding samples. Analysis of transcript data revealed high differential gene expression through seedling emergence followed by a period of less change. Differential gene expression increased gradually to day 8, and showed two days, 5 and 7, with a very high proportion of up-regulated genes, including transcription factor/signaling genes. Network cartography using spring embedding revealed two primary clusters of highly correlated metabolites, which appear to reflect temporally distinct metabolic states. Principle Component Analyses of both sets of profiling data produced a chronological spread of time points, which would be expected of a developmental series. The network cartography of the transcript data produced two distinct clusters comprising days 0 to 2 and days 3 to 8, whereas the corresponding analysis of metabolite data revealed a shift of day 2 into the day 3 to 8 group. A metabolite and transcript pair-wise correlation analysis encompassing all time points gave a set of 237 highly significant correlations. Of 129 genes correlated to sucrose, 44 of them were known to be sucrose responsive including a number of transcription factors.</p> <p>Conclusions</p> <p>Microarray analysis during germination and establishment revealed major transitions in transcriptional activity at time points potentially associated with developmental transitions. Network cartography using spring-embedding indicate that a shift in the state of nutritionally important metabolites precedes a major shift in the transcriptional state going from germination to seedling emergence. Pair-wise linear correlations of transcript and metabolite levels identified many genes known to be influenced by metabolites, and provided other targets to investigate metabolite regulation of gene expression during seedling establishment.</p

In Silico Evaluation of Predicted Regulatory Interactions in Arabidopsis thaliana

<p>Abstract</p> <p>Background</p> <p>Prediction of transcriptional regulatory mechanisms in <it>Arabidopsis </it>has become increasingly critical with the explosion of genomic data now available for both gene expression and gene sequence composition. We have shown in previous work <abbrgrp><abbr bid="B1">1</abbr></abbrgrp>, that a combination of correlation measurements and <it>cis</it>-regulatory element (CRE) detection methods are effective in predicting targets for candidate transcription factors for specific case studies which were validated. However, to date there has been no quantitative assessment as to which correlation measures or CRE detection methods used alone or in combination are most effective in predicting TF→target relationships on a genome-wide scale.</p> <p>Results</p> <p>We tested several widely used methods, based on correlation (Pearson and Spearman Rank correlation) and <it>cis-</it>regulatory element (CRE) detection (≥1 CRE or CRE over-representation), to determine which of these methods individually or in combination is the most effective by various measures for making regulatory predictions. To predict the regulatory targets of a transcription factor (TF) of interest, we applied these methods to microarray expression data for genes that were regulated over treatment and control conditions in wild type (WT) plants. Because the chosen data sets included identical experimental conditions used on TF over-expressor or T-DNA knockout plants, we were able to test the TF→target predictions made using microarray data from WT plants, with microarray data from mutant/transgenic plants. For each method, or combination of methods, we computed sensitivity, specificity, positive and negative predictive value and the F-measure of balance between sensitivity and positive predictive value (precision). This analysis revealed that the ≥1 CRE and Spearman correlation (used alone or in combination) were the most balanced CRE detection and correlation methods, respectively with regard to their power to accurately predict regulatory-target interactions.</p> <p>Conclusion</p> <p>These findings provide an approach and guidance for researchers interested in predicting transcriptional regulatory mechanisms using microarray data that they generate (or microarray data that is publically available) combined with CRE detection in promoter sequence data.</p

Predictive network modeling of the high-resolution dynamic plant transcriptome in response to nitrate

International audienceABSTRACT: BACKGROUND: Nitrate, acting as both a nitrogen source and a signaling molecule, controls many aspects of plant development. However, gene networks involved in plant adaptation to fluctuating nitrate environments have not yet been identified. RESULTS: Here we use time-series transcriptome data to decipher gene relationships and consequently to build core regulatory networks involved in Arabidopsis root adaptation to nitrate provision. The experimental approach has been to monitor genome-wide responses to nitrate at 3, 6, 9, 12, 15 and 20 minutes, using Affymetrix ATH1 gene chips. This high-resolution time course analysis demonstrated that the previously known primary nitrate response is actually preceded by a very fast gene expression modulation, involving genes and functions needed to prepare plants to use or reduce nitrate. A state-space model inferred from this microarray time-series data successfully predicts gene behavior in unlearnt conditions. CONCLUSIONS: The experiments and methods allow us to propose a temporal working model for nitrate-driven gene networks. This network model is tested both in silico and experimentally. For example, the over-expression of a predicted gene hub encoding a transcription factor induced early in the cascade indeed leads to the modification of the kinetic nitrate response of sentinel genes such as NIR, NIA2, and NRT1.1, and several other transcription factors. The potential nitrate /hormone connections implicated by this time-series data is also evaluated

UVR8 mediated spatial differences as a prerequisite for UV-B induced inflorescence phototropism

In Arabidopsis hypocotyls, phototropins are the dominant photoreceptors for the positive phototropism response towards unilateral ultraviolet-B (UV-B) radiation. We report a stark contrast of response mechanism with inflorescence stems with a central role for UV RESISTANCE LOCUS 8 (UVR8). The perception of UV-B occurs mainly in the epidermis and cortex with a lesser contribution of the endodermis. Unilateral UV-B exposure does not lead to a spatial difference in UVR8 protein levels but does cause differential UVR8 signal throughout the stem with at the irradiated side 1) increase of the transcription factor ELONGATED HYPOCOTYL 5 (HY5), 2) an associated strong activation of flavonoid biosynthesis genes and flavonoid accumulation, 3) increased GA2oxidase expression, diminished gibberellin1 levels and accumulation of DELLA protein REPRESSOR OF GA1 (RGA) and, 4) increased expression of the auxin transport regulator, PINOID, contributing to local diminished auxin signalling. Our molecular findings are in support of the Blaauw theory (1919), suggesting that differential growth occurs trough unilateral photomorphogenic growth inhibition. Together the data indicate phototropin independent inflorescence phototropism through multiple locally UVR8-regulated hormone pathways

Tree Peony Species Are a Novel Resource for Production of α-Linolenic Acid

Tree peony is known worldwide for its excellent ornamental and medical values, but recent reports that their seeds contain over 40% α-linolenic acid (ALA), an essential fatty acid for humans drew additional interest of biochemists. To understand the key factors that contribute to this rich accumulation of ALA, we carried out a comprehensive study of oil accumulation in developing seeds of nine wild tree peony species. The fatty acid content and composition was highly variable among the nine species; however, we selected a high- (P. rockii) and low-oil (P. lutea) accumulating species for a comparative transcriptome analysis. Similar to other oilseed transcriptomic studies, upregulation of select genes involved in plastidial fatty acid synthesis, and acyl editing, desaturation and triacylglycerol assembly in the endoplasmic reticulum was noted in seeds of P. rockii relative to P. lutea. Also, in association with the ALA content, transcript levels for fatty acid desaturases (SAD, FAD2 and FAD3), which encode for enzymes necessary for polyunsaturated fatty acid synthesis were higher in P. rockii compared to P. lutea. We further showed that the overexpression of PrFAD2 and PrFAD3 in Arabidopsis increased linoleic and α-linolenic acid content, respectively and modulated their final ratio in the seed oil. In conclusion, we identified the key steps that contribute to efficient ALA synthesis and validated the necessary desaturases in P. rockii that are responsible for not only increasing oil content but also modulating 18:2/18:3 ratio in seeds. Together, these results will aid to improve essential fatty acid content in seeds of tree peonies and other crops of agronomic interest