Gene network analysis of Arabidopsis thaliana flower development through dynamic gene perturbations

Floral Homeotic Factors: A Question of Specificity

MADS-domain transcription factors are involved in the control of a multitude of processes in eukaryotes, and in plants, they play particularly important roles during reproductive development. Among the members of this large family of regulatory proteins are the floral organ identity factors, which specify the identities of the different types of floral organs in a combinatorial manner. Much has been learned over the past three decades about the function of these master regulators. For example, it has been shown that they have similar DNA-binding activities and that their genome-wide binding patterns exhibit large overlaps. At the same time, it appears that only a minority of binding events lead to changes in gene expression and that the different floral organ identity factors have distinct sets of target genes. Thus, binding of these transcription factors to the promoters of target genes alone may not be sufficient for their regulation. How these master regulators achieve specificity in a developmental context is currently not well understood. Here, we review what is known about their activities and highlight open questions that need to be addressed to gain more detailed insights into the molecular mechanisms underlying their functions. We discuss evidence for the involvement of cofactors as well as the results from studies on transcription factors in animals that may be instructive for a better understanding of how the floral organ identity factors achieve regulatory specificity.

Floral Induction Systems for the Study of Arabidopsis Flower Development

Assessing the molecular changes that occur over the course of flower development is hampered by difficulties in isolating sufficient amounts of floral tissue at specific developmental stages. This is especially problematic when investigating molecular events at early stages of Arabidopsis flower development, as floral buds are minute and are initiated sequentially so that a single flower on an inflorescence is at a given developmental stage. Moreover, young floral buds are hidden by older flowers, which presents an additional challenge for dissection. To circumvent these issues, floral induction systems that allow the simultaneous induction of a large number of flowers on the inflorescence of a single plant were developed. To allow the plant community to avail of the full benefits of these systems, we address some common problems that can be encountered when growing these plants and collecting floral buds for analysis.

Multi-Omics Methods Applied to Flower Development

Developmental processes in multicellular organisms depend on the proficiency of cells to orchestrate different gene expression programs. Over the past years, several studies of reproductive organ development have considered genomic analyses of transcription factors and global gene expression changes, modeling complex gene regulatory networks. Nevertheless, the dynamic view of developmental processes requires, as well, the study of the proteome in its expression, complexity, and relationship with the transcriptome. In this chapter, we describe a dual extraction method-for protein and RNA-for the characterization of genome expression at proteome level and its correlation to transcript expression data. We also present a shotgun proteomic procedure (LC-MS/MS) followed by a pipeline for the imputation of missing values in mass spectrometry results.

Carbon-nitrogen trading in symbiotic nodules depends on magnesium import

Symbiotic nitrogen fixation provides large amounts of nitrogen for global agricultural systems with little environmental or economic costs. The basis of symbiosis is the nutrient exchange occurring between legumes and rhizobia, but key regulators controlling nutrient exchange are largely unknown. Here, we reveal that magnesium (Mg), an important nutrient factor that preferentially accumulates in inner cortical cells of soybean nodules, shows the most positive correlation with nodule carbon (C) import and nitrogen (N) export. We further identified a pair of Mg transporter genes, GmMGT4 and GmMGT5, that are specifically expressed in the nodule cortex, modulating both nodule Mg import and C-N transport processes. The GmMGT4&5-dependent Mg import activates the activity of a plasmodesmata-located β-1,3-glucanase GmBG2 and consequently keeps plasmodesmata permeable for C-N transport in nodule inner cortical cells. Our studies discovered an important regulating pathway for host plants fine-tuning nodule C-N trading to achieve optimal growth, which may be helpful for optimizing nutrient management for soybean production.

Spatiotemporal Gene Expression Profiling and Network Inference: A Roadmap for Analysis, Visualization, and Key Gene Identification

Gene expression data analysis and the prediction of causal relationships within gene regulatory networks (GRNs) have guided the identification of key regulatory factors and unraveled the dynamic properties of biological systems. However, drawing accurate and unbiased conclusions requires a comprehensive understanding of relevant tools, computational methods, and their workflows. The topics covered in this chapter encompass the entire workflow for GRN inference including: (1) experimental design; (2) RNA sequencing data processing; (3) differentially expressed gene (DEG) selection; (4) clustering prior to inference; (5) network inference techniques; and (6) network visualization and analysis. Moreover, this chapter aims to present a workflow feasible and accessible for plant biologists without a bioinformatics or computer science background. To address this need, TuxNet, a user-friendly graphical user interface that integrates RNA sequencing data analysis with GRN inference, is chosen for the purpose of providing a detailed tutorial.

The intervening domain is required for DNA-binding and functional identity of plant MADS transcription factors

The MADS transcription factors (TF) are an ancient eukaryotic protein family. In plants, the family is divided into two main lineages. Here, we demonstrate that DNA binding in both lineages absolutely requires a short amino acid sequence C-terminal to the MADS domain (M domain) called the Intervening domain (I domain) that was previously defined only in type II lineage MADS. Structural elucidation of the MI domains from the floral regulator, SEPALLATA3 (SEP3), shows a conserved fold with the I domain acting to stabilise the M domain. Using the floral organ identity MADS TFs, SEP3, APETALA1 (AP1) and AGAMOUS (AG), domain swapping demonstrate that the I domain alters genome-wide DNA-binding specificity and dimerisation specificity. Introducing AG carrying the I domain of AP1 in the Arabidopsis ap1 mutant resulted in strong complementation and restoration of first and second whorl organs. Taken together, these data demonstrate that the I domain acts as an integral part of the DNA-binding domain and significantly contributes to the functional identity of the MADS TF.

MADS transcription factors regulate multiple aspects of plant development. Here the authors show that the intervening I domain is conserved in both type I and type II plant MADS lineages and contributes to the functional identity of the protein by influencing both DNA binding activity and dimerisation specificity.

AINTEGUMENTA and AINTEGUMENTA-LIKE6 directly regulate floral homeotic, growth, and vascular development genes in young Arabidopsis flowers

Arabidopsis flower primordia give rise to organ primordia in stereotypical positions within four concentric whorls. Floral organ primordia in each whorl undergo distinct developmental programs to become one of four organ types (sepals, petals, stamens, and carpels). The Arabidopsis transcription factors AINTEGUMENTA (ANT) and AINTEGUMENTA-LIKE6 (AIL6) are required for correct positioning of floral organ initiation, contribute to the specification of floral organ identity, and regulate the growth and morphogenesis of developing floral organs. To gain insight into the molecular mechanisms by which ANT and AIL6 contribute to floral organogenesis, we identified the genome-wide binding sites of both ANT and AIL6 in stage 3 flower primordia, the developmental stage at which sepal primordia become visible and class B and C floral homeotic genes are first expressed. AIL6 binds to a subset of ANT sites, suggesting that AIL6 regulates some but not all of the same target genes as ANT. ANT- and AIL6-binding sites are associated with genes involved in many biological processes related to meristem and flower organ development. Comparison of genes associated with both ANT and AIL6 ChIP-Seq peaks and those differentially expressed after perturbation of ANT and/or AIL6 activity identified likely direct targets of ANT and AIL6 regulation. These include class B and C floral homeotic genes, growth regulatory genes, and genes involved in vascular development.

ANT and AIL6: masters of the master regulators during flower development

This article comments on: Krizek BA, Bantle AT, Heflin JM, Han H, Freese NH, Loraine AE. 2021. AINTEGUMENTA and AINTEGUMENTA-LIKE6 directly regulate floral homeotic, growth, and vascular development genes in young Arabidopsis flowers. Journal of Experimental Botany 72, 5478–5493.

Expression of KNUCKLES in the Stem Cell Domain Is Required for Its Function in the Control of Floral Meristem Activity in Arabidopsis

In the model plant Arabidopsis thaliana, the zinc-finger transcription factor KNUCKLES (KNU) plays an important role in the termination of floral meristem activity, a process that is crucial for preventing the overgrowth of flowers. The KNU gene is activated in floral meristems by the floral organ identity factor AGAMOUS (AG), and it has been shown that both AG and KNU act in floral meristem control by directly repressing the stem cell regulator WUSCHEL (WUS), which leads to a loss of stem cell activity. When we re-examined the expression pattern of KNU in floral meristems, we found that KNU is expressed throughout the center of floral meristems, which includes, but is considerably broader than the WUS expression domain. We therefore hypothesized that KNU may have additional functions in the control of floral meristem activity. To test this, we employed a gene perturbation approach and knocked down KNU activity at different times and in different domains of the floral meristem. In these experiments we found that early expression in the stem cell domain, which is characterized by the expression of the key meristem regulatory gene CLAVATA3 (CLV3), is crucial for the establishment of KNU expression. The results of additional genetic and molecular analyses suggest that KNU represses floral meristem activity to a large extent by acting on CLV3. Thus, KNU might need to suppress the expression of several meristem regulators to terminate floral meristem activity efficiently.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

Whole-transcriptome analysis of differentially expressed genes in the mutant and normal capitula of Chrysanthemum morifolium

BACKGROUND

Chrysanthemum morifolium is one of the most economically important and popular floricultural crops in the family Asteraceae. Chrysanthemum flowers vary considerably in terms of colors and shapes. However, the molecular mechanism controlling the development of chrysanthemum floral colors and shapes remains an enigma. We analyzed a cut-flower chrysanthemum variety that produces normal capitula composed of ray florets with normally developed pistils and purple corollas and mutant capitula comprising ray florets with green corollas and vegetative buds instead of pistils.

RESULTS

We conducted a whole-transcriptome analysis of the differentially expressed genes (DEGs) in the mutant and normal capitula using third-generation and second-generation sequencing techniques. We identified the DEGs between the mutant and normal capitula to reveal important regulators underlying the differential development. Many transcription factors and genes related to the photoperiod and GA pathways, floral organ identity, and the anthocyanin biosynthesis pathway were differentially expressed between the normal and mutant capitula. A qualitative analysis of the pigments in the florets of normal and mutant capitula indicated anthocyanins were synthesized and accumulated in the florets of normal capitula, but not in the florets of mutant capitula. These results provide clues regarding the molecular basis of the replacement of Chrysanthemum morifolium ray florets with normally developed pistils and purple corollas with mutant ray florets with green corollas and vegetative buds. Additionally, the study findings will help to elucidate the molecular mechanisms underlying floral organ development and contribute to the development of techniques for studying the regulation of flower shape and color, which may enhance chrysanthemum breeding.

CONCLUSIONS

The whole-transcriptome analysis of DEGs in mutant and normal C. morifolium capitula described herein indicates the anthocyanin deficiency of the mutant capitula may be related to the mutation that replaces ray floret pistils with vegetative buds. Moreover, pistils may be required for the anthocyanin biosynthesis in the corollas of chrysanthemum ray florets.

The Arabidopsis transcription factor AINTEGUMENTA orchestrates patterning genes and auxin signaling in the establishment of floral growth and form

Understanding how flowers form is an important problem in plant biology, as human food supply depends on flower and seed production. Flower development also provides an excellent model for understanding how cell division, expansion and differentiation are coordinated during organogenesis. In the model plant Arabidopsis thaliana, floral organogenesis requires AINTEGUMENTA (ANT) and AINTEGUMENTA-LIKE 6 (AIL6)/PLETHORA 3 (PLT3), two members of the Arabidopsis AINTEGUMENTA-LIKE/PLETHORA (AIL/PLT) transcription factor family. Together, ANT and AIL6/PLT3 regulate aspects of floral organogenesis, including floral organ initiation, growth, identity specification and patterning. Previously, we used RNA-Seq to identify thousands of genes with disrupted expression in ant ail6 mutant flowers, indicating that ANT and AIL6/PLT3 influence a vast transcriptional network. The immediate downstream targets of ANT and AIL6/PLT3 in flowers are unknown, however. To identify direct targets of ANT regulation, we performed an RNA-Seq time-course experiment in which we induced ANT activity in transgenic plants bearing an ANT-glucocorticoid receptor fusion construct. In addition, we performed a ChIP-Seq experiment that identified ANT binding sites in developing flowers. These experiments identified 200 potential ANT target genes based on their proximity to ANT binding sites and differential expression in response to ANT. These 200 candidate target genes were involved in functions such as polarity specification, floral organ development, meristem development and auxin signaling. In addition, we identified several genes associated with lateral organ growth that may mediate the role of ANT in organ size control. These results reveal new features of the ANT transcriptional network by linking ANT to previously unknown regulatory targets.

Gene Regulatory Network Inference: Connecting Plant Biology and Mathematical Modeling

Plant responses to environmental and intrinsic signals are tightly controlled by multiple transcription factors (TFs). These TFs and their regulatory connections form gene regulatory networks (GRNs), which provide a blueprint of the transcriptional regulations underlying plant development and environmental responses. This review provides examples of experimental methodologies commonly used to identify regulatory interactions and generate GRNs. Additionally, this review describes network inference techniques that leverage gene expression data to predict regulatory interactions. These computational and experimental methodologies yield complex networks that can identify new regulatory interactions, driving novel hypotheses. Biological properties that contribute to the complexity of GRNs are also described in this review. These include network topology, network size, transient binding of TFs to DNA, and competition between multiple upstream regulators. Finally, this review highlights the potential of machine learning approaches to leverage gene expression data to predict phenotypic outputs.

Photosynthetic activity of reproductive organs

During seed development, carbon is reallocated from maternal tissues to support germination and subsequent growth. As this pool of resources is depleted post-germination, the plant begins autotrophic growth through leaf photosynthesis. Photoassimilates derived from the leaf are used to sustain the plant and form new organs, including other vegetative leaves, stems, bracts, flowers, fruits, and seeds. In contrast to the view that reproductive tissues act only as resource sinks, many studies demonstrate that flowers, fruits, and seeds are photosynthetically active. The photosynthetic contribution to development is variable between these reproductive organs and between species. In addition, our understanding of the developmental control of photosynthetic activity in reproductive organs is vastly incomplete. A further complication is that reproductive organ photosynthesis (ROP) appears to be particularly important under suboptimal growth conditions. Therefore, the topic of ROP presents the community with a challenge to integrate the fields of photosynthesis, development, and stress responses. Here, we attempt to summarize our understanding of the contribution of ROP to development and the molecular mechanisms underlying its control.

Inducible Promoter Systems for Gene Perturbation Experiments in Arabidopsis

Assessing molecular changes that occur through altering a gene's activity is often hampered by difficulties that arise due to the typically static nature of the introduced perturbation. This is especially problematic when investigating molecular events at specific stages and/or in certain tissues or organs during Arabidopsis development. To circumvent these issues, we have employed chemically inducible artificial microRNAs (amiRNAs) for the specific knockdown of developmental regulators. For our own research, we have combined this gene perturbation approach with a floral induction system, which allows the simultaneous induction of a large number of flowers on the inflorescence of a single plant, and the ability to knock down a gene's activity at any given stage of development. To enable the plant community to avail of the full benefits of these systems, we describe, in this chapter, strategies for amiRNA-mediated gene perturbations and address some common problems that can be encountered when generating inducible amiRNA constructs, growing these plants, and collecting floral buds for analysis.

Floral homeotic proteins modulate the genetic program for leaf development to suppress trichome formation in flowers

As originally proposed by Goethe in 1790, floral organs are derived from leaf-like structures. The conversion of leaves into different types of floral organ is mediated by floral homeotic proteins, which, as described by the ABCE model of flower development, act in a combinatorial manner. However, how these transcription factors bring about this transformation process is not well understood. We have previously shown that floral homeotic proteins are involved in suppressing the formation of branched trichomes, a hallmark of leaf development, on reproductive floral organs of Arabidopsis Here, we present evidence that the activities of the C function gene AGAMOUS (AG) and the related SHATTERPROOF1/2 genes are superimposed onto the regulatory network that controls the distribution of trichome formation in an age-dependent manner. We show that AG regulates cytokinin responses and genetically interacts with the organ polarity gene KANADI1 to suppress trichome initiation on gynoecia. Thus, our results show that parts of the genetic program for leaf development remain active during flower formation but have been partially rewired through the activities of the floral homeotic proteins.

A Specific Knockdown of Transcription Factor Activities in Arabidopsis

Transcription factors are pivotal for the control of development and the response of organisms to changes in the environment. Therefore, a detailed understanding of their functions is of central importance for biology. Over the years, different experimental methods have been developed to study the activities of transcription factors in plants. These methods include perturbation assays, where the activity of a given transcription factor is disrupted and subsequently, the resulting effects are monitored using molecular, genomic, or physiological approaches. Perturbation assays can also be used to distinguish primary roles of transcription factors of interest from secondary effects. Thus, molecular genetic experiments after perturbation can be advantageous or even necessary for the precise understanding of transcription factor function at a certain stage of plant development or in a single tissue or organ type. In this chapter, we describe several commonly used techniques to knock down transcription factor activities and provide detailed information on how those techniques are employed in the model plant Arabidopsis thaliana.

Profiling Histone Modifications in Synchronized Floral Tissues for Quantitative Resolution of Chromatin and Transcriptome Dynamics

Covalent histone modifications and their effects on chromatin state and accessibility play a key role in the regulation of gene expression in eukaryotes. To gain insights into their functions during plant growth and development, the distribution of histone modifications can be analyzed at a genome-wide scale through chromatin immunoprecipitation assays followed by sequencing of the isolated genomic DNA. Here, we present a protocol for systematic analysis of the distribution and dynamic changes of selected histone modifications, during flower development in the model plant Arabidopsis thaliana. This protocol utilizes a previously established floral induction system to synchronize flower development, which allows the collection of sufficient plant material for analysis by genomic technologies. In this chapter, we describe how to use this system to study, from the same set of samples, chromatin and transcriptome dynamics during early stages of flower formation.

From Genes to Networks: Characterizing Gene-Regulatory Interactions in Plants

Plants, like other eukaryotes, have evolved complex mechanisms to coordinate gene expression during development, environmental response, and cellular homeostasis. Transcription factors (TFs), accompanied by basic cofactors and posttranscriptional regulators, are key players in gene-regulatory networks (GRNs). The coordinated control of gene activity is achieved by the interplay of these factors and by physical interactions between TFs and DNA. Here, we will briefly outline recent technological progress made to elucidate GRNs in plants. We will focus on techniques that allow us to characterize physical interactions in GRNs in plants and to analyze their regulatory consequences. Targeted manipulation allows us to test the relevance of specific gene-regulatory interactions. The combination of genome-wide experimental approaches with mathematical modeling allows us to get deeper insights into key-regulatory interactions and combinatorial control of important processes in plants.

Exploring Plant Co-Expression and Gene-Gene Interactions with CORNET 3.0

Selecting and filtering a reference expression and interaction dataset when studying specific pathways and regulatory interactions can be a very time-consuming and error-prone task. In order to reduce the duplicated efforts required to amass such datasets, we have created the CORNET (CORrelation NETworks) platform which allows for easy access to a wide variety of data types: coexpression data, protein-protein interactions, regulatory interactions, and functional annotations. The CORNET platform outputs its results in either text format or through the Cytoscape framework, which is automatically launched by the CORNET website.CORNET 3.0 is the third iteration of the web platform designed for the user exploration of the coexpression space of plant genomes, with a focus on the model species Arabidopsis thaliana. Here we describe the platform: the tools, data, and best practices when using the platform. We indicate how the platform can be used to infer networks from a set of input genes, such as upregulated genes from an expression experiment. By exploring the network, new target and regulator genes can be discovered, allowing for follow-up experiments and more in-depth study. We also indicate how to avoid common pitfalls when evaluating the networks and how to avoid over interpretation of the results.All CORNET versions are available at http://bioinformatics.psb.ugent.be/cornet/ .

Gene-regulatory networks controlling inflorescence and flower development in Arabidopsis thaliana

Reproductive development in plants is controlled by complex and intricate gene-regulatory networks of transcription factors. These networks integrate the information from endogenous, hormonal and environmental regulatory pathways. Many of the key players have been identified in Arabidopsis and other flowering plant species, and their interactions and molecular modes of action are being elucidated. An emerging theme is that there is extensive crosstalk between different pathways, which can be accomplished at the molecular level by modulation of transcription factor activity or of their downstream targets. In this review, we aim to summarize current knowledge on transcription factors and epigenetic regulators that control basic developmental programs during inflorescence and flower morphogenesis in the model plant Arabidopsis thaliana. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Patterns of gene expression during Arabidopsis flower development from the time of initiation to maturation

Background

The formation of flowers is one of the main model systems to elucidate the molecular mechanisms that control developmental processes in plants. Although several studies have explored gene expression during flower development in the model plant Arabidopsis thaliana on a genome-wide scale, a continuous series of expression data from the earliest floral stages until maturation has been lacking. Here, we used a floral induction system to close this information gap and to generate a reference dataset for stage-specific gene expression during flower formation.

Results

Using a floral induction system, we collected floral buds at 14 different stages from the time of initiation until maturation. Using whole-genome microarray analysis, we identified 7,405 genes that exhibit rapid expression changes during flower development. These genes comprise many known floral regulators and we found that the expression profiles for these regulators match their known expression patterns, thus validating the dataset. We analyzed groups of co-expressed genes for over-represented cellular and developmental functions through Gene Ontology analysis and found that they could be assigned specific patterns of activities, which are in agreement with the progression of flower development. Furthermore, by mapping binding sites of floral organ identity factors onto our dataset, we were able to identify gene groups that are likely predominantly under control of these transcriptional regulators. We further found that the distribution of paralogs among groups of co-expressed genes varies considerably, with genes expressed predominantly at early and intermediate stages of flower development showing the highest proportion of such genes.

Conclusions

Our results highlight and describe the dynamic expression changes undergone by a large number of genes during flower development. They further provide a comprehensive reference dataset for temporal gene expression during flower formation and we demonstrate that it can be used to integrate data from other genomics approaches such as genome-wide localization studies of transcription factor binding sites.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1699-6) contains supplementary material, which is available to authorized users.