Search

Inferring phenotypic outcomes from genomic features is both a promise and challenge for systems biology. Using gene expression data to predict phenotypic outcomes, and functionally validating the genes with predictive powers are two challenges we address in this study. We applied an evolutionarily informed machine learning approach to predict phenotypes based on transcriptome responses shared both within and across species. Specifically, we exploited the phenotypic diversity in nitrogen use efficiency and evolutionarily conserved transcriptome responses to nitrogen treatments across Arabidopsis accessions and maize varieties. We demonstrate that using evolutionarily conserved nitrogen responsive genes is a biologically principled approach to reduce the feature dimensionality in machine learning that ultimately improved the predictive power of our gene-to-trait models. Further, we functionally validated seven candidate transcription factors with predictive power for NUE outcomes in Arabidopsis and one in maize. Moreover, application of our evolutionarily informed pipeline to other species including rice and mice models underscores its potential to uncover genes affecting any physiological or clinical traits of interest across biology, agriculture, or medicine., Predicting complex phenotypes from genomic information is still a challenge. Here, the authors use an evolutionarily informed machine learning approach within and across species to predict genes affecting nitrogen utilization in crops, and show their approach is also useful in mammalian systems.

The TARGET system allows for the rapid identification of direct regulated gene targets of transcription factors (TFs). It employs the transient transformation of plant protoplasts with inducible nuclear entry of the TF and subsequent transcriptomic and/or ChIP-seq analysis. The ability to separate direct TF-target gene regulatory interactions from indirect downstream responses and the significantly shorter amount of time required to perform the assay, compared to the generation of transgenics, make this plant cell-based approach a valuable tool for a higher throughput approach to identify the genome-wide targets of multiple TFs, to build validated transcriptional networks in plants. Here, we describe the use of the TARGET system in Arabidopsis seedling root protoplasts to map the gene regulatory network downstream of transcription factors-of-interest.

The TARGET system allows for the rapid identification of direct regulated gene targets of transcription factors (TFs). It employs the transient transformation of plant protoplasts with inducible nuclear entry of the TF and subsequent transcriptomic and/or ChIP-seq analysis. The ability to separate direct TF-target gene regulatory interactions from indirect downstream responses and the significantly shorter amount of time required to perform the assay, compared to the generation of transgenics, makes this plant cell-based approach a valuable tool for a higher through-put approach to identify the genome-wide targets of multiple TFs, to build validated transcriptional networks in plants. Here, we describe the use of the TARGET system in Arabidopsis seedling root protoplasts to map the gene regulatory network downstream of transcription factors-of-interest. NB. The original uploaded pdf contained a typo; PEG solution contains 0.4 M mannitol (not 0.4 mM).

The ability to accurately predict the causal relationships from transcription factors to genes would greatly enhance our understanding of transcriptional dynamics. This could lead to applications in which one or more transcription factors could be manipulated to effect a change in genes leading to the enhancement of some desired trait. Here we present a method called OutPredict that constructs a model for each gene based on time series (and other) data and that predicts gene's expression in a previously unseen subsequent time point. The model also infers causal relationships based on the most important transcription factors for each gene model, some of which have been validated from previous physical experiments. The method benefits from known network edges and steady-state data to enhance predictive accuracy. Our results across B. subtilis, Arabidopsis, E.coli, Drosophila and the DREAM4 simulated in silico dataset show improved predictive accuracy ranging from 40% to 60% over other state-of-the-art methods. We find that gene expression models can benefit from the addition of steady-state data to predict expression values of time series. Finally, we validate, based on limited available data, that the influential edges we infer correspond to known relationships significantly more than expected by chance or by state-of-the-art methods.

We discuss how both nitrogen and water availability combine to impact plant biological responses at the molecular, transcriptomic, and physiological level., Nitrogen (N) and water (W) are crucial inputs for plant survival as well as costly resources for agriculture. Given their importance, the molecular mechanisms that plants rely on to signal changes in either N or W status have been under intense scrutiny. However, how plants sense and respond to the combination of N and W signals at the molecular level has received scant attention. The purpose of this review is to shed light on what is currently known about how plant responses to N are impacted by W status. We review classic studies which detail how N and W combinations have both synergistic and antagonistic effects on key plant traits, such as root architecture and stomatal aperture. Recent molecular studies of N and W interactions show that mutations in genes involved in N metabolism affect drought responses, and vice versa. Specifically, perturbing key N signaling genes may lead to changes in drought-responsive gene expression programs, which is supported by a meta-analysis we conduct on available transcriptomic data. Additionally, we cite studies that show how combinatorial transcriptional responses to N and W status might drive crop phenotypes. Through these insights, we suggest research strategies that could help to develop crops adapted to marginal soils depleted in both N and W, an important task in the face of climate change.

As sessile organisms, plants are finely tuned to respond dynamically to developmental, circadian, and environmental cues. Genome-wide studies investigating these types of cues have uncovered the intrinsically different ways they can impact gene expression over time. Recent advances in single cell sequencing and time-based bioinformatic algorithms are now beginning to reveal the dynamics of these time-based responses within individual cells and plant tissues. Here, we review what these techniques have revealed about the spatiotemporal nature of gene regulation, paying particular attention to the three distinct ways in which plant tissues are time sensitive. (i) First, we discuss how studying plant cell identity can reveal developmental trajectories hidden in pseudotime. (ii) Next, we present evidence which indicates that plant cell types keep their own local time through tissue specific regulation of the circadian clock. (iii) Finally, we review what determines the speed of environmental signaling responses, and how they can be contingent on developmental and circadian time. By these means, this review sheds light on how these different scales of time-based responses can act with tissue and cell-type specificity to elicit changes in whole plant systems.

The Atacama Desert in Chile—hyperarid and with high–ultraviolet irradiance levels—is one of the harshest environments on Earth. Yet, dozens of species grow there, including Atacama-endemic plants. Herein, we establish the Talabre–Lejia transect (TLT) in the Atacama as an unparalleled natural laboratory to study plant adaptation to extreme environmental conditions. We characterized climate, soil, plant, and soil–microbe diversity at 22 sites (every 100 m of altitude) along the TLT over a 10-y period. We quantified drought, nutrient deficiencies, large diurnal temperature oscillations, and pH gradients that define three distinct vegetational belts along the altitudinal cline. We deep-sequenced transcriptomes of 32 dominant plant species spanning the major plant clades, and assessed soil microbes by metabarcoding sequencing. The top-expressed genes in the 32 Atacama species are enriched in stress responses, metabolism, and energy production. Moreover, their root-associated soils are enriched in growth-promoting bacteria, including nitrogen fixers. To identify genes associated with plant adaptation to harsh environments, we compared 32 Atacama species with the 32 closest sequenced species, comprising 70 taxa and 1,686,950 proteins. To perform phylogenomic reconstruction, we concatenated 15,972 ortholog groups into a supermatrix of 8,599,764 amino acids. Using two codon-based methods, we identified 265 candidate positively selected genes (PSGs) in the Atacama plants, 64% of which are located in Pfam domains, supporting their functional relevance. For 59/184 PSGs with an Arabidopsis ortholog, we uncovered functional evidence linking them to plant resilience. As some Atacama plants are closely related to staple crops, these candidate PSGs are a “genetic goldmine” to engineer crop resilience to face climate change.

Plant responses to multiple environmental stimuli must be integrated to enable them to adapt their metabolism and development. Light and nitrogen (N) are two such stimuli whose downstream signaling pathways must be intimately connected to each other to control plant energy status. Here, we describe the functional role of the WRKY1 transcription factor in controlling genome-wide transcriptional reprogramming of Arabidopsis (Arabidopsis thaliana) leaves in response to individual and combined light and N signals. This includes a cross-regulatory network consisting of 724 genes regulated by WRKY1 and involved in both N and light signaling pathways. The loss of WRKY1 gene function has marked effects on the light and N response of genes involved in N uptake and assimilation (primary metabolism) as well as stress response pathways (secondary metabolism). Our results at the transcriptome and at the metabolite analysis level support a model in which WRKY1 enables plants to activate genes involved in the recycling of cellular carbon resources when light is limiting but N is abundant and upregulate amino acid metabolism when both light and N are limiting. In this potential energy conservation mechanism, WRKY1 integrates information about cellular N and light energy resources to trigger changes in plant metabolism.

International audience

Charting a temporal path in gene networks requires linking early transcription factor (TF)-triggered events to downstream effects. We scale-up a cell-based TF-perturbation assay to identify direct regulated targets of 33 nitrogen (N)-early response TFs encompassing 88% of N-responsive Arabidopsis genes. We uncover a duality where each TF is an inducer and repressor, and in vitro cis-motifs are typically specific to regulation directionality. Validated TF-targets (71,836) are used to refine precision of a time-inferred root network, connecting 145 N-responsive TFs and 311 targets. These data are used to chart network paths from direct TF1-regulated targets identified in cells to indirect targets responding only in planta via Network Walking. We uncover network paths from TGA1 and CRF4 to direct TF2 targets, which in turn regulate 76% and 87% of TF1 indirect targets in planta, respectively. These results have implications for N-use and the approach can reveal temporal networks for any biological system., Temporal control of transcriptional networks enables organisms to adapt to changing environment. Here, the authors use a scaled-up cell-based assay to identify direct targets of nitrogen-early responsive transcription factors and validate a network path mediating dynamic nitrogen signaling in Arabidopsis.

Plants need to cope with strong variations of nitrogen availability in the soil. Although many molecular players are being discovered concerning how plants perceive NO3− provision, it is less clear how plants recognize a lack of nitrogen. Following nitrogen removal, plants activate their nitrogen starvation response (NSR), which is characterized by the activation of very high-affinity nitrate transport systems (NRT2.4 and NRT2.5) and other sentinel genes involved in N remobilization such as GDH3. Using a combination of functional genomics via transcription factor perturbation and molecular physiology studies, we show that the transcription factors belonging to the HHO subfamily are important regulators of NSR through two potential mechanisms. First, HHOs directly repress the high-affinity nitrate transporters, NRT2.4 and NRT2.5. hho mutants display increased high-affinity nitrate transport activity, opening up promising perspectives for biotechnological applications. Second, we show that reactive oxygen species (ROS) are important to control NSR in wild-type plants and that HRS1 and HHO1 overexpressors and mutants are affected in their ROS content, defining a potential feed-forward branch of the signaling pathway. Taken together, our results define the relationships of two types of molecular players controlling the NSR, namely ROS and the HHO transcription factors. This work (i) up opens perspectives on a poorly understood nutrient-related signaling pathway and (ii) defines targets for molecular breeding of plants with enhanced NO3− uptake.

The TARGET system allows for the rapid identification of direct regulated gene targets of transcription factors (TFs). It employs the transient transformation of plant protoplasts with inducible nuclear entry of the TF and subsequent transcriptomic and/or ChIP-seq analysis. The ability to separate direct TF-target gene regulatory interactions from indirect downstream responses and the significantly shorter amount of time required to perform the assay, compared to the generation of transgenics, makes this plant cell-based approach a valuable tool for a higher through-put approach to identify the genome-wide targets of multiple TFs, to build validated transcriptional networks in plants. Here, we describe the use of the TARGET system in Arabidopsis seedling root protoplasts to map the gene regulatory network downstream of transcription factors-of-interest.

All aspects of transcription and its regulation involve dynamic events. However, capturing these dynamic events in gene regulatory networks (GRNs) offers both a promise and a challenge. The promise is that capturing and modeling the dynamic changes in GRNs will allow us to understand how organisms adapt to a changing environment. The ability to mount a rapid transcriptional response to environmental changes is especially important in nonmotile organisms such as plants. The challenge is to capture these dynamic, genome-wide events and model them in GRNs. In this review, we cover recent progress in capturing dynamic interactions of transcription factors with their targets—at both the local and genome-wide levels—and how they are used to learn how GRNs operate as a function of time. We also discuss recent advances that employ time-based machine learning approaches to forecast gene expression at future time points, a key goal of systems biology.

Deciphering gene regulatory networks (GRNs) is both a promise and challenge of systems biology. The promise lies in identifying key transcription factors (TFs) that enable an organism to react to changes in its environment. The challenge lies in validating GRNs that involve hundreds of TFs with hundreds of thousands of interactions with their genome-wide targets experimentally determined by high-throughput sequencing. To address this challenge, we developed ConnecTF, a species-independent, web-based platform that integrates genome-wide studies of TF–target binding, TF–target regulation, and other TF-centric omic datasets and uses these to build and refine validated or inferred GRNs. We demonstrate the functionality of ConnecTF by showing how integration within and across TF–target datasets uncovers biological insights. Case study 1 uses integration of TF–target gene regulation and binding datasets to uncover TF mode-of-action and identify potential TF partners for 14 TFs in abscisic acid signaling. Case study 2 demonstrates how genome-wide TF–target data and automated functions in ConnecTF are used in precision/recall analysis and pruning of an inferred GRN for nitrogen signaling. Case study 3 uses ConnecTF to chart a network path from NLP7, a master TF in nitrogen signaling, to direct secondary TF2s and to its indirect targets in a Network Walking approach. The public version of ConnecTF (https://ConnecTF.org) contains 3,738,278 TF–target interactions for 423 TFs in Arabidopsis, 839,210 TF–target interactions for 139 TFs in maize (Zea mays), and 293,094 TF–target interactions for 26 TFs in rice (Oryza sativa). The database and tools in ConnecTF will advance the exploration of GRNs in plant systems biology applications for model and crop species.

Deciphering gene regulatory networks (GRNs) is both a promise and challenge of systems biology. The promise is identifying key transcription factors (TFs) that enable an organism to react to changes in its environment. The challenge is constructing GRNs that involve hundreds of TFs and hundreds of thousands of interactions with their genome-wide target genes validated by high-throughput sequencing. To address this challenge, we developed ConnecTF, a species-independent web-based platform for constructing validated GRNs and to refine inferred GRNs via combined analysis of genome-wide studies of TF-target gene binding, TF-target regulation and other TF-centric omic data. We demonstrate the functionality of ConnecTF in three case studies, showing how integration within and across TF-target datasets uncovers biological insights. Case study 1 uses integration of TF-target gene regulation and binding datasets to uncover mode-of-action and identify potential TF partners for 14 TFs in abscisic acid signaling. Case study 2 demonstrates how genome-wide TF-target data and automated functions in ConnecTF are used to conduct precision/recall analysis and pruning of an inferred GRN for nitrogen signaling. In case study 3, we use ConnecTF to chart a network path from NLP7, a master TF in nitrogen signaling, to direct secondary TF2s, to its indirect targets, in an approach called Network Walking. The public version of ConnecTF (https://ConnecTF.org) contains 3,738,278 TF-target interactions for 423 TFs in Arabidopsis, and 839,210 TF-target interactions for 139 TFs in maize. The database and tools in ConnecTF should advance the exploration of GRNs in plant systems biology applications for models and crops.

International audience; Nitrogen (N) is an essential macronutrient for plants, and a major limiting factor for plant growth and crop production. Nitrate is the main source of N available for plants in agricultural soils and in many natural environments. Sustaining agricultural productivity is of paramount importance in the current scenario of increasing world population, diversification of crop uses, and climate change. Plant productivity for major crops around the world is still supported by excess application of N-based fertilizers with detrimental economic and environmental impacts. Thus, understanding how plants regulate nitrate uptake and metabolism is key for developing new crops with enhanced N use efficiency (NUE) and to cope with future world food demands. The study of plant responses to nitrate has gained considerable interest over the last thirty years. This review provides an overview of key findings in nitrate research, spanning biochemistry, molecular genetics, genomics and systems biology. We discuss how we reached our current understanding of nitrate transport, local and systemic nitrate sensing/signaling, and the regulatory networks underlying nitrate-controlled outputs in plants. We hope this summary serves not only as a time-line and information repository, but also as a base to outline important open questions for future research.

An increase in nutrient dose leads to proportional increases in crop biomass and agricultural yield. However, the molecular underpinnings of this nutrient dose–response are largely unknown. To investigate, we assayed changes in the Arabidopsis root transcriptome to different doses of nitrogen (N)—a key plant nutrient—as a function of time. By these means, we found that rate changes of genome-wide transcript levels in response to N-dose could be explained by a simple kinetic principle: the Michaelis–Menten (MM) model. Fitting the MM model allowed us to estimate the maximum rate of transcript change ( V max ), as well as the N-dose at which one-half of V max was achieved ( K m ) for 1,153 N-dose–responsive genes. Since transcription factors (TFs) can act in part as the catalytic agents that determine the rates of transcript change, we investigated their role in regulating N-dose–responsive MM-modeled genes. We found that altering the abundance of TGA1, an early N-responsive TF, perturbed the maximum rates of N-dose transcriptomic responses ( V max ), K m , as well as the rate of N-dose–responsive plant growth. We experimentally validated that MM-modeled N-dose–responsive genes included both direct and indirect TGA1 targets, using a root cell TF assay to detect TF binding and/or TF regulation genome-wide. Taken together, our results support a molecular mechanism of transcriptional control that allows an increase in N-dose to lead to a proportional change in the rate of genome-wide expression and plant growth.

Chromatin modification has gained increased attention for its role in the regulation of plant responses to environmental changes, but the specific mechanisms and molecular players remain elusive. Here, we show that the Arabidopsis (Arabidopsis thaliana) histone methyltransferase SET DOMAIN GROUP8 (SDG8) mediates genome-wide changes in H3K36 methylation at specific genomic loci functionally relevant to nitrate treatments. Moreover, we show that the specific H3K36 methyltransferase encoded by SDG8 is required for canonical RNA processing, and that RNA isoform switching is more prominent in the sdg8-5 deletion mutant than in the wild type. To demonstrate that SDG8-mediated regulation of RNA isoform expression is functionally relevant, we examined a putative regulatory gene, CONSTANS, CO-like, and TOC1 101 (CCT101), whose nitrogen-responsive isoform-specific RNA expression is mediated by SDG8. We show by functional expression in shoot cells that the different RNA isoforms of CCT101 encode distinct regulatory proteins with different effects on genome-wide transcription. We conclude that SDG8 is involved in plant responses to environmental nitrogen supply, affecting multiple gene regulatory processes including genome-wide histone modification, transcriptional regulation, and RNA processing, and thereby mediating developmental and metabolic processes related to nitrogen use.

Dynamic reprogramming of gene regulatory networks (GRNs) enables organisms to rapidly respond to environmental perturbation. However, the underlying transient interactions between transcription factors (TFs) and genome-wide targets typically elude biochemical detection. Here, we capture both stable and transient TF-target interactions genome-wide within minutes after controlled TF nuclear import using time-series chromatin immunoprecipitation (ChIP-seq) and/or DNA adenine methyltransferase identification (DamID-seq). The transient TF-target interactions captured uncover the early mode-of-action of NIN-LIKE PROTEIN 7 (NLP7), a master regulator of the nitrogen signaling pathway in plants. These transient NLP7 targets captured in root cells using temporal TF perturbation account for 50% of NLP7-regulated genes not detectably bound by NLP7 in planta. Rapid and transient NLP7 binding activates early nitrogen response TFs, which we validate to amplify the NLP7-initiated transcriptional cascade. Our approaches to capture transient TF-target interactions genome-wide can be applied to validate dynamic GRN models for any pathway or organism of interest., Conventional methods cannot reveal transient transcription factors (TFs) and targets interactions. Here, Alvarez et al. capture both stable and transient TF-target interactions by time-series ChIP-seq and/or DamID-seq in a cell-based TF perturbation system and show NLP7 as a master TF to initiate a rapid nitrogen-response cascade.

Significance Our study exploits time—the relatively unexplored fourth dimension of gene regulatory networks (GRNs)—to learn the temporal transcriptional logic underlying dynamic nitrogen (N) signaling in plants. We introduce several conceptual innovations to the analysis of time-series data in the area of predictive GRNs. Our resulting network now provides the “transcriptional logic” for transcription factor perturbations aimed at improving N-use efficiency, an important issue for global food production in marginal soils and for sustainable agriculture. More broadly, the combination of the time-based approaches we develop and deploy can be applied to uncover the temporal “transcriptional logic” for any response system in biology, agriculture, or medicine., This study exploits time, the relatively unexplored fourth dimension of gene regulatory networks (GRNs), to learn the temporal transcriptional logic underlying dynamic nitrogen (N) signaling in plants. Our “just-in-time” analysis of time-series transcriptome data uncovered a temporal cascade of cis elements underlying dynamic N signaling. To infer transcription factor (TF)-target edges in a GRN, we applied a time-based machine learning method to 2,174 dynamic N-responsive genes. We experimentally determined a network precision cutoff, using TF-regulated genome-wide targets of three TF hubs (CRF4, SNZ, and CDF1), used to “prune” the network to 155 TFs and 608 targets. This network precision was reconfirmed using genome-wide TF-target regulation data for four additional TFs (TGA1, HHO5/6, and PHL1) not used in network pruning. These higher-confidence edges in the GRN were further filtered by independent TF-target binding data, used to calculate a TF “N-specificity” index. This refined GRN identifies the temporal relationship of known/validated regulators of N signaling (NLP7/8, TGA1/4, NAC4, HRS1, and LBD37/38/39) and 146 additional regulators. Six TFs—CRF4, SNZ, CDF1, HHO5/6, and PHL1—validated herein regulate a significant number of genes in the dynamic N response, targeting 54% of N-uptake/assimilation pathway genes. Phenotypically, inducible overexpression of CRF4 in planta regulates genes resulting in altered biomass, root development, and 15NO3− uptake, specifically under low-N conditions. This dynamic N-signaling GRN now provides the temporal “transcriptional logic” for 155 candidate TFs to improve nitrogen use efficiency with potential agricultural applications. Broadly, these time-based approaches can uncover the temporal transcriptional logic for any biological response system in biology, agriculture, or medicine.

Plant responses to multiple stimuli must be integrated to trigger transcriptional cascades that lead to changes in plant metabolism and development. Light (L) and nitrogen (N) are two signaling pathways that are intimately connected to each other and to plant energy status. Here, we describe the functional role of the WRKY1 transcription factor in mediating the regulation between L and N signaling pathways in Arabidopsis thaliana. WRKY1 participates in genome-wide transcriptional reprogramming in leaves in response to individual and combined L and N signals. A regulatory network was identified, consisting of 724 genes regulated by WRKY1 and involved in both N and L signaling pathways. The loss of WRKY1 gene function has marked effects on the L and N response of genes involved in N uptake and assimilation (primary metabolism) as well as stress response pathways (secondary metabolism). Our results support a model in which WRKY1 enables plants to activate genes involved in the recycling of cellular carbon resources when L is limiting but N is abundant, and up-regulate amino acid metabolism genes when both L and N are limiting. In this potential energy conservation mechanism, WRKY1 integrates responses to N and light-energy status to trigger changes in plant metabolism.SummaryBased on transcriptome analysis, the WRKY1 transcription factor mediates regulation of nitrogen and light signaling pathways in a potential energy conservation mechanism.

Changes in nutrient dose have dramatic effects on gene expression and development. One outstanding question is whether organisms respond to changes in absolute nutrient amount (moles) vs. its concentration in water (molarity). This question is particularly relevant to plants, as soil drying can alter nutrient concentration, without changing its absolute amount. To compare the effects of amount vs. concentration, we expose rice to a factorial matrix varying the dose of nitrogen (N) and water (W) over a range of combinations, and quantify transcriptome and phenotype responses. Using linear models, we identify distinct dose responses to either N-moles, W-volume, N-molarity (N/W), or their synergistic interaction (N×W). Importantly, genes whose expression patterns are best explained by N-dose and W interactions (N/W or N×W) in seedlings are associated with crop outcomes in replicated field trials. Such N-by-W responsive genes may assist future efforts to develop crops resilient to increasingly arid, low nutrient soils., Scarcity of water and nutrients limit crop yields. Here Swift et al. show that gene expression in rice responds differently to changes in the absolute amount of nitrogen available compared to nitrogen concentration and identify expression profiles associated with crop performance in arid, low-nutrient soils.

The long-distance signaling network allowing a plant to properly develop its root system is crucial to optimize root foraging in areas where nutrients are available. Cytokinin is an essential element of the systemic signaling network leading to the enhancement of lateral root proliferation in areas where nitrate is available. Here, we explore more precisely: (i) which particular traits of lateral root growth (density and length of emerged lateral roots) are the targets of systemic signaling in a context of heterogeneous nitrate supply; and (ii) if the systemic signaling depends only on cytokinin or on a combination of several signalings.

The ability to accurately predict the causal relationships from transcription factors to genes would greatly enhance our understanding of transcriptional dynamics. This could lead to applications in which one or more transcription factors could be manipulated to effect a change in genes leading to the enhancement of some desired trait. Here we present a method called OutPredict that constructs a model for each gene based on time series (and other) data and that predicts gene's expression in a previously unseen subsequent time point. The model also infers causal relationships based on the most important transcription factors for each gene model, some of which have been validated from previous physical experiments. The method benefits from known network edges and steady-state data to enhance predictive accuracy. Our results across B. subtilis, Arabidopsis, E.coli, Drosophila and the DREAM4 simulated in silico dataset show improved predictive accuracy ranging from 40% to 60% over other state-of-the-art methods. We find that gene expression models can benefit from the addition of steady-state data to predict expression values of time series. Finally, we validate, based on limited available data, that the influential edges we infer correspond to known relationships significantly more than expected by chance or by state-of-the-art methods.

Plants need to cope with strong variations in the nitrogen content of the soil solution. Although many molecular actors are being discovered concerning how plants perceive NO3- provision, it is less clear how plants recognize a lack of Nitrogen. Indeed, following N removal plants activate their Nitrogen Starvation Response (NSR) being characterized in particular by the activation of very high affinity nitrate transport systems (NRT2.4, NRT2.5) and other sentinel genes such as GDH3. Here we show using a combination of functional genomics (via TF perturbation) and molecular physiology studies, that the GARP Transcription Factors (TFs) belonging the HHO sub-family are important regulators of the NSR through two potential mechanisms. First, HHOs directly repress NRT2.4 and NRT2.5 high-affinity nitrate transporters. Genotypes affected in HHO genes (mutants and overexpressors) display modified high-affinity nitrate transport activities opening interesting perspectives in biotechnology applications. Second, we show that Reactive Oxygen Species (ROS) are important to control NSR in wild type plants and that HRS1 and HHO1 overexpressors are affected in their ROS content, defining a potential feedforward branch of the signaling pathway. Taken together our results define two new classes of molecular actors in the control of NSR including ROS and the first transcription factors to date. This work (i) opens perspectives on a poorly understood nutrient related signaling pathway, and (ii) defines targets for molecular breeding of plants with enhanced NO3- uptake.

Chromatin immunoprecipitation (ChIP) is a widely used method to map the position of DNA-binding proteins such as histones and transcription factors (TFs) upon their interaction with particular regions of the genome. To examine the genomic distribution of a TF in specific cell types in response to a change in nitrogen concentration, we developed a micro-ChIP (μChIP) protocol that requires only ~5000 Arabidopsis cells transiently expressing the Arabidopsis TF Basic Leucine Zipper 1 (bZIP1) fused to the glucocorticoid receptor (GR) domain that mediates nuclear import in the presence of dexamethasone. The DNA fragments obtained from the immunoprecipitation of bZIP1-DNA complexes were analyzed by next-generation sequencing (ChIP-seq), which helped uncover genome-wide associations between a bZIP1 and its targets in plant cells upon fluctuations in nitrogen availability.

Chromatin immunoprecipitation (ChIP) is a widely used method to map the position of DNA-binding proteins such as histones and transcription factors (TFs) upon their interaction with particular regions of the genome. To examine the genomic distribution of a TF in specific cell types in response to a change in nitrogen concentration, we developed a micro-ChIP (μChIP) protocol that requires only ~5000 Arabidopsis cells transiently expressing the Arabidopsis TF Basic Leucine Zipper 1 (bZIP1) fused to the glucocorticoid receptor (GR) domain that mediates nuclear import in the presence of dexamethasone. The DNA fragments obtained from the immunoprecipitation of bZIP1-DNA complexes were analyzed by next-generation sequencing (ChIP-seq), which helped uncover genome-wide associations between a bZIP1 and its targets in plant cells upon fluctuations in nitrogen availability.

Modeling dynamic gene regulatory networks (GRNs) is a new frontier in systems biology. It has special implications for plants, whose survival requires rapid deployment of GRNs in response to environmental changes. However, capturing and dissecting transient interactions of transcription factors (TFs) and their targets in GRNs remains a considerable experimental challenge. Here we review recent progress in understanding GRNs as a function of time and discuss the relevance of these findings in plants to studies in other eukaryotes. We cover progress in profiling and modeling time-course transcriptome changes across plant species and the insights they have provided into the regulatory mechanisms underlying these temporal transcriptional responses, with a focus on the dynamic behavior of TFs. Lastly, we review state-of-the-art techniques to monitor the single-molecule dynamics of TFs in vivo. Together, these advances have helped develop new models for dynamic transcriptional control with relevance across eukaryotes.

Understanding how transcription factor (TF) binding is related to gene regulation is a moving target. We recently uncovered genome-wide evidence for a "Hit-and-Run" model of transcription. In this model, a master TF "hits" a target promoter to initiate a rapid response to a signal. As the "hit" is transient, the model invokes recruitment of partner TFs to sustain transcription over time. Following the "run", the master TF "hits" other targets to propagate the response genome-wide. As such, a TF may act as a "catalyst" to mount a broad and acute response in cells that first sense the signal, while the recruited TF partners promote long-term adaptive behavior in the whole organism. This "Hit-and-Run" model likely has broad relevance, as TF perturbation studies across eukaryotes show small overlaps between TF-regulated and TF-bound genes, implicating transient TF-target binding. Here, we explore this "Hit-and-Run" model to suggest molecular mechanisms and its biological relevance.

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 ; Argentina Fil: Shrivastava, Stuti. University of New York; Estados Unidos Fil: Katari, Manpreet S.. University of New York; Estados Unidos Fil: Coruzzi, Gloria M.. University of New York; Estados Unidos

Shaping of root architecture is a quintessential developmental response that involves the concerted action of many different cell types, is highly dynamic, and underpins root plasticity. To determine to what extent the environmental regulation of lateral root development is a product of cell-type preferential activities, we tracked transcriptomic responses to two different treatments that both change root development in Arabidopsis thaliana at an unprecedented level of temporal detail. We found that individual transcripts are expressed with a very high degree of temporal and spatial specificity, yet biological processes are commonly regulated, in a mechanism we term response nonredundancy. Using causative gene network inference to compare the genes regulated in different cell types and during responses to nitrogen and a biotic interaction, we found that common transcriptional modules often regulate the same gene families but control different individual members of these families, specific to response and cell type. This reinforces that the activity of a gene cannot be defined simply as molecular function; rather, it is a consequence of spatial location, expression timing, and environmental responsiveness.

The architecture of plant roots affects essential functions including nutrient and water uptake, soil anchorage, and symbiotic interactions. Root architecture comprises many features that arise from the growth of the primary and lateral roots. These root features are dictated by the genetic background but are also highly responsive to the environment. Thus, root system architecture (RSA) represents an important and complex trait that is highly variable, affected by genotype × environment interactions, and relevant to survival/performance. Quantification of RSA in Arabidopsis (Arabidopsis thaliana) using plate-based tissue culture is a very common and relatively rapid assay, but quantifying RSA represents an experimental bottleneck when it comes to medium- or high-throughput approaches used in mutant or genotype screens. Here, we present RootScape, a landmark-based allometric method for rapid phenotyping of RSA using Arabidopsis as a case study. Using the software AAMToolbox, we created a 20-point landmark model that captures RSA as one integrated trait and used this model to quantify changes in the RSA of Arabidopsis (Columbia) wild-type plants grown under different hormone treatments. Principal component analysis was used to compare RootScape with conventional methods designed to measure root architecture. This analysis showed that RootScape efficiently captured nearly all the variation in root architecture detected by measuring individual root traits and is 5 to 10 times faster than conventional scoring. We validated RootScape by quantifying the plasticity of RSA in several mutant lines affected in hormone signaling. The RootScape analysis recapitulated previous results that described complex phenotypes in the mutants and identified novel gene × environment interactions.

Plants form the basis of the food webs that sustain animal life. Exogenous factors, such as nutrients and sunlight, and endogenous factors, such as hormones, cooperate to control both the growth and the development of plants. We assessed how Arabidopsis thaliana integrated nutrient and hormone signaling pathways to control root growth and development by investigating the effects of combinatorial treatment with the nutrients nitrate and ammonium; the hormones auxin, cytokinin, and abscisic acid; and all binary combinations of these factors. We monitored and integrated short-term genome-wide changes in gene expression over hours and longterm effects on root development and architecture over several days. Our analysis revealed trends in nutrient and hormonal signal cross-talk and feedback, including responses that exhibited logic gate behavior, which means that they were triggered only when specific combinations of signals were present. From the data, we developed a multivariate network model comprising the signaling molecules, the early gene expression modulation, and the subsequent changes in root phenotypes. This multivariate network model pinpoints several genes that play key roles in the control of root development and may help understand how eukaryotes manage multifactorial signaling inputs.

Nitrate (NO(3)(-)) is a key nutrient as well as a signaling molecule that impacts both metabolism and development of plants. Understanding the complexity of the regulatory networks that control nitrate uptake, metabolism, and associated responses has the potential to provide solutions that address the major issues of nitrate pollution and toxicity that threaten agricultural and ecological sustainability and human health. Recently, major advances have been made in cataloguing the nitrate transcriptome and in identifying key components that mediate nitrate signaling. In this perspective, we describe the genes involved in nitrate regulation and how they influence nitrate transport and assimilation, and we discuss the role of systems biology approaches in elucidating the gene networks involved in NO(3)(-) signaling adaptation to fluctuating environments.

One of the most striking examples of plant developmental plasticity to changing environmental conditions is the modulation of root system architecture (RSA) in response to nitrate supply. Despite the fundamental and applied significance of understanding this process, the molecular mechanisms behind nitrate-regulated changes in developmental programs are still largely unknown. Small RNAs (sRNAs) have emerged as master regulators of gene expression in plants and other organisms. To evaluate the role of sRNAs in the nitrate response, we sequenced sRNAs from control and nitrate-treated Arabidopsis seedlings using the 454 sequencing technology. miR393 was induced by nitrate in these experiments. miR393 targets transcripts that code for a basic helix-loop-helix (bHLH) transcription factor and for the auxin receptors TIR1, AFB1, AFB2, and AFB3. However, only AFB3 was regulated by nitrate in roots under our experimental conditions. Analysis of the expression of this miR393/ AFB3 module, revealed an incoherent feed-forward mechanism that is induced by nitrate and repressed by N metabolites generated by nitrate reduction and assimilation. To understand the functional role of this N-regulatory module for plant development, we analyzed the RSA response to nitrate in AFB3 insertional mutant plants and in miR393 overexpressors. RSA analysis in these plants revealed that both primary and lateral root growth responses to nitrate were altered. Interestingly, regulation of RSA by nitrate was specifically mediated by AFB3 , indicating that miR393/ AFB3 is a unique N-responsive module that controls root system architecture in response to external and internal N availability in Arabidopsis .

We use measures of congruence on a combined expressed sequenced tag genome phylogeny to identify proteins that have potential significance in the evolution of seed plants. Relevant proteins are identified based on the direction of partitioned branch and hidden support on the hypothesis obtained on a 16-species tree, constructed from 2,557 concatenated orthologous genes. We provide a general method for detecting genes or groups of genes that may be under selection in directions that are in agreement with the phylogenetic pattern. Gene partitioning methods and estimates of the degree and direction of support of individual gene partitions to the overall data set are used. Using this approach, we correlate positive branch support of specific genes for key branches in the seed plant phylogeny. In addition to basic metabolic functions, such as photosynthesis or hormones, genes involved in posttranscriptional regulation by small RNAs were significantly overrepresented in key nodes of the phylogeny of seed plants. Two genes in our matrix are of critical importance as they are involved in RNA-dependent regulation, essential during embryo and leaf development. These are Argonaute and the RNA-dependent RNA polymerase 6 found to be overrepresented in the angiosperm clade. We use these genes as examples of our phylogenomics approach and show that identifying partitions or genes in this way provides a platform to explain some of the more interesting organismal differences among species, and in particular, in the evolution of plants.

BMAA is a cycad-derived glutamate receptor agonist that causes a two- to three-fold increase in hypocotyl elongation on Arabidopsis seedlings grown in the light. To probe the role of plant glutamate receptors and their downstream mediators, we utilized a previously described genetic screen to identify a novel, BMAA insensitive morphology (bim) mutant, bim409. The normal BMAA-induced hypocotyl elongation response observed on wild-type seedlings grown in the light is impaired in the bim409 mutant. This BMAA-induced phenotype is light-specific, as the bim409 mutant exhibits normal hypocotyl elongation in etiolated (dark grown) plants (+ or − BMAA). The mutation in bim409 was identified to be in a gene encoding the Proteosomal Regulatory Particle AAA-ATPase-3 (RPT3). Possible roles of the proteosome in Glu-mediated signaling in plants is discussed.