Stepwise origin and evolution of a transcriptional activator and repressor system integrating nutrient signaling in plants

Plants possess a unique transcriptional regulatory system in which two related MYB-related transcription factors (TFs) coordinate gene expression according to phosphate (Pi) and nitrogen (N) availability. The Phosphorus Starvation Response (PSR) type TFs are transcriptional activators integrating the cellular Pi sensing machinery and gene regulation majorly under Pi starvation. The Hypersensitivity To Low Pi-Elicited Primary Root Shortening (HRS) type TFs are transcriptional repressors integrating the Pi and N availability signals through different feedback loops. They are highly connected through multiple signaling loops to finetune the transcriptional responses according to nutrient availability. Molecular functions of these TFs are fairly uncovered in model systems; however, how plants evolved this activator-repressor system is currently unknown. In this study, using sensitive evolutionary analysis, we identified a stepwise origin of the PSR-HRS regulatory system in plants. The PSR TFs were originated before the split of Prasinodermophyta and Chlorophyta. The HRS TFs were originated later in the Streptrophycean algae. We also identified the asymmetric expansion of this TF repertoire in land plants majorly shaped by genome duplication and triplication events. The phylogenetic reconstruction coupled with motif analysis revealed that the origin of the specific accessory motifs is a major contributing factor in the functional divergence which led to the evolution of different sub-families preceding the angiosperm radiation. The spatiotemporal gene expression analysis in different developmental stages and nutrient availability conditions in angiosperms identified a critical role of expression divergence in shaping the functions of these TF families which is essential for adaptive plasticity of plants.