Spatiotemporal analysis identifies ABF2 and ABF3 as key hubs of endodermal response to nitrate.

Nitrate is a nutrient and a potent signal that impacts global gene expression in plants. However, the regulatory factors controlling temporal and cell type-specific nitrate responses remain largely unknown. We assayed nitrate-responsive transcriptome changes in five major root cell types of the Arabidopsis thaliana root as a function of time. We found that gene-expression response to nitrate is dynamic and highly localized and predicted cell type-specific transcription factor (TF)-target interactions. Among cell types, the endodermis stands out as having the largest and most connected nitrate-regulatory gene network. ABF2 and ABF3 are major hubs for transcriptional responses in the endodermis cell layer. We experimentally validated TF-target interactions for ABF2 and ABF3 by chromatin immunoprecipitation followed by sequencing and a cell-based system to detect TF regulation genome-wide. Validated targets of ABF2 and ABF3 account for more than 50% of the nitrate-responsive transcriptome in the endodermis. Moreover, ABF2 and ABF3 are involved in nitrate-induced lateral root growth. Our approach offers an unprecedented spatiotemporal resolution of the root response to nitrate and identifies important components of cell-specific gene regulatory networks.

ConnecTF: A platform to integrate transcription factor-gene interactions and validate regulatory networks.

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.

Transient genome-wide interactions of the master transcription factor NLP7 initiate a rapid nitrogen-response cascade.

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.

SDG8-Mediated Histone Methylation and RNA Processing Function in the Response to Nitrate Signaling.

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.

Network Walking charts transcriptional dynamics of nitrogen signaling by integrating validated and predicted genome-wide interactions.

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 transcriptional logic of dynamic regulatory networks underlying nitrogen signaling and use in plants.

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.