High-resolution temporal transcript profiling during Arabidopsis thaliana gynoecium morphogenesis uncovers the chronology of gene regulatory network activity and reveals novel developmental regulators

The gynoecium is the most complex organ formed by the flowering plants. It encloses the ovules, provides a surface for pollen contact and self-incompatibility reactions, allows pollen tube growth and, post fertilization, and develops into the fruit. Consequently, the regulation of gynoecium morphogenesis is complex and appropriate timing of this process in part determines reproductive success. However, little is known about the global control of gynoecium development, even though many regulatory genes have been characterized. Here, we characterized dynamic gene expression changes using laser-microdissected gynoecium tissue from four developmental stages in Arabidopsis. We provide a high-resolution map of global expression dynamics during gynoecium morphogenesis and link these to the gynoecium interactome. We reveal groups of genes acting together early and others acting late in morphogenesis. Clustering of co-expressed genes enables comparisons between the leaf, shoot apex, and gynoecium transcriptomes allowing the dissection of common and distinct regulators. Furthermore, our results lead to the discovery of the LESSER FERTILITY1-4 (LEF1-4) genes, which, when mutated, lead to impaired gynoecium expansion, illustrating that global transcriptome analyses reveal yet unknown developmental regulators. Our data show that highly interacting proteins, such as SEPALLATA3, AGAMOUS, and TOPLESS are expressed more evenly during development, but switch interactors in time, whereas stage-specific proteins have only few interactors. Our analysis connects specific transcriptional regulator activities, protein interactions, and underlying metabolic processes towards the development of a dynamic network model for gynoecium development.