Architecture and dynamics of the abscisic acid gene regulatory network

The plant hormone abscisic acid (ABA) regulates essential processes in plant development and responsiveness to abiotic and biotic stresses. ABA perception triggers a post-translational signaling cascade that elicits the ABA gene regulatory network (GRN), encompassing hundreds of transcription factors (TFs) and thousands of transcribed genes. To further our knowledge of this GRN, we performed an RNA-seq time series experiment consisting of 14 time points in the 16 h following a one-time ABA treatment of 5-week-old Arabidopsis rosettes. During this time course, ABA rapidly changed transcription levels of 7151 genes, which were partitioned into 44 coexpressed modules that carry out diverse biological functions. We integrated our time-series data with publicly available TF-binding site data, motif data, and RNA-seq data of plants inhibited in translation, and predicted (i) which TFs regulate the different coexpression clusters, (ii) which TFs contribute the most to target gene amplitude, (iii) timing of engagement of different TFs in the ABA GRN, and (iv) hierarchical position of TFs and their targets in the multi-tiered ABA GRN. The ABA GRN was found to be highly interconnected and regulated at different amplitudes and timing by a wide variety of TFs, of which the bZIP family was most prominent, and upregulation of genes encompassed more TFs than downregulation. We validated our network models in silico with additional public TF-binding site data and transcription data of selected TF mutants. Finally, using a drought assay we found that the Trihelix TF GT3a is likely an ABA-induced positive regulator of drought tolerance.