Flowers are icons in developmental studies of complex structures. The vast majority of 250,000 angiosperm plant species have flowers with a conserved organ plan bearing sepals, petals, stamens, and carpels in the center. The combinatorial model for the activity of the so-called ABC homeotic floral genes has guided extensive experimental studies in Arabidopsis thaliana and many other plant species. However, a mechanistic and dynamical explanation for the ABC model and prevalence among flowering plants is lacking. Here, we put forward a simple discrete model that postulates logical rules that formally summarize published ABC and non-ABC gene interaction data for Arabidopsis floral organ cell fate determination and integrates this data into a dynamic network model. This model shows that all possible initial conditions converge to few steady gene activity states that match gene expression profiles observed experimentally in primordial floral organ cells of wild-type and mutant plants. Therefore, the network proposed here provides a dynamical explanation for the ABC model and shows that precise signaling pathways are not required to restrain cell types to those found in Arabidopsis, but these are rather determined by the overall gene network dynamics. Furthermore, we performed robustness analyses that clearly show that the cell types recovered depend on the network architecture rather than on specific values of the model's gene interaction parameters. These results support the hypothesis that such a network constitutes a developmental module, and hence provide a possible explanation for the overall conservation of the ABC model and overall floral plan among angiosperms. In addition, we have been able to predict the effects of differences in network architecture between Arabidopsis and Petunia hybrida.

Organismal development and many cell biological processes are organized in a modular fashion, where regulatory molecules form groups with many interactions within a group and few interactions between groups. Thus, the activity of elements within a module depends little on elements outside of it. Modularity facilitates the production of heritable variation and of evolutionary innovations. There is no consensus on how modularity might evolve, especially for modules in development. We show that modularity can increase in gene regulatory networks as a byproduct of specialization in gene activity. Such specialization occurs after gene regulatory networks are selected to produce new gene activity patterns that appear in a specific body structure or under a specific environmental condition. Modules that arise after specialization in gene activity comprise genes that show concerted changes in gene activities. This and other observations suggest that modularity evolves because it decreases interference between different groups of genes. Our work can explain the appearance and maintenance of modularity through a mechanism that is not contingent on environmental change. We also show how modularity can facilitate co-option, the utilization of existing gene activity to build new gene activity patterns, a frequent feature of evolutionary innovations.

BackgroundMany important evolutionary adaptations originate in the modification of gene regulatory circuits to produce new gene activity phenotypes. How do evolving populations sift through an astronomical number of circuits to find circuits with new adaptive phenotypes? The answer may often involve phenotypic plasticity. Phenotypic plasticity allows a genotype to produce different - alternative - phenotypes after non-genetic perturbations that include gene expression noise, environmental change, or epigenetic modification.ResultsWe here analyze a well-studied model of gene regulatory circuits. A circuit's genotype encodes the regulatory interactions among circuit genes, and its phenotype corresponds to a stable gene activity pattern the circuit forms. For this model, we study how genotypes are arranged in genotype space, where the distance between two genotypes reflects the number of regulatory mutations that set those genotypes apart. Specifically, we address whether this arrangement favors adaptive evolution mediated by plasticity. We find that plasticity facilitates the origin of genotypes that produce a new phenotype in response to non-genetic perturbations. We also find that selection can then stabilize the new phenotype genetically, allowing it to become a circuit's dominant gene expression phenotype. These are generic properties of the circuits we study here.ConclusionsTaken together, our observations suggest that phenotypic plasticity frequently facilitates the evolution of novel beneficial gene activity patterns in gene regulatory circuits.

In contrast to the classical view of development as a preprogrammed and deterministic process, recent studies have demonstrated that stochastic perturbations of highly non-linear systems may underlie the emergence and stability of biological patterns. Herein, we address the question of whether noise contributes to the generation of the stereotypical temporal pattern in gene expression during flower development. We modeled the regulatory network of organ identity genes in the Arabidopsis thaliana flower as a stochastic system. This network has previously been shown to converge to ten fixed-point attractors, each with gene expression arrays that characterize inflorescence cells and primordial cells of sepals, petals, stamens, and carpels. The network used is binary, and the logical rules that govern its dynamics are grounded in experimental evidence. We introduced different levels of uncertainty in the updating rules of the network. Interestingly, for a level of noise of around 0.5–10%, the system exhibited a sequence of transitions among attractors that mimics the sequence of gene activation configurations observed in real flowers. We also implemented the gene regulatory network as a continuous system using the Glass model of differential equations, that can be considered as a first approximation of kinetic-reaction equations, but which are not necessarily equivalent to the Boolean model. Interestingly, the Glass dynamics recover a temporal sequence of attractors, that is qualitatively similar, although not identical, to that obtained using the Boolean model. Thus, time ordering in the emergence of cell-fate patterns is not an artifact of synchronous updating in the Boolean model. Therefore, our model provides a novel explanation for the emergence and robustness of the ubiquitous temporal pattern of floral organ specification. It also constitutes a new approach to understanding morphogenesis, providing predictions on the population dynamics of cells with different genetic configurations during development.

Background: Dynamical models are instrumental for exploring the way information required to generate robust developmental patterns arises from complex interactions among genetic and non-genetic factors. We address this fundamental issue of developmental biology studying the leaf and root epidermis of Arabidopsis. We propose an experimentally-grounded model of gene regulatory networks (GRNs) that are coupled by protein diffusion and comprise a meta-GRN implemented on cellularised domains.

Genes and proteins form complex dynamical systems or gene regulatory networks (GRN) that can reach several steady states (attractors). These may be associated with distinct cell types. In plants, the ABC combinatorial model establishes the necessary gene combinations for floral organ cell specification. We have developed dynamic gene regulatory network (GRN) models to understand how the combinatorial selection of gene activity is established during floral organ primordia specification as a result of the concerted action of ABC and non-ABC genes. Our analyses have shown that the floral organ specification GRN reaches six attractors with gene configurations observed in primordial cell types during early stages of flower development and four that correspond to regions of the inflorescence meristem. This suggests that it is the overall GRN dynamics rather than precise signals that underlie the ABC model. Furthermore, our analyses suggest that the steady states of the GRN are robust to random alterations of the logical functions that define the gene interactions. Here we have updated the GRN model and have systematically altered the outputs of all the logical functions and addressed in which cases the original attractors are recovered. We then reduced the original threestate GRN to a two-state (Boolean) GRN and performed the same systematic perturbation analysis. Interestingly, the Boolean GRN reaches the same number and type of attractors as reached by the three-state GRN, and it responds to perturbations in a qualitatively identical manner as the original GRN. These results suggest that a Boolean model is sufficient to capture the dynamical features of the floral network and provide additional support for the robustness of the floral GRN. These findings further support that the GRN model provides a dynamical explanation for the ABC model and that the floral GRN robustness could be behind the widespread conservation of the floral plan among eudicotyledoneous plants. Other aspects of evolution of flower organ arrangement and ABC gene expression patterns are discussed in the context of the approach proposed here.Á lvaro Chaos, Max Aldana and Elena Alvarez-Buylla contributed equally to this work.

In Arabidopsis thaliana, leaf and root epidermis hairs exhibit contrasting spatial arrangements even though the genetic networks regulating their respective cell-fate determination have very similar structures and components. We integrated available experimental data for leaf and root hair patterning in dynamic network models which may be reduced to activator-inhibitor models. This integration yielded expected results for these kinds of dynamic models, including striped and dotted cell patterns which are characteristic of root and leaf epidermis, respectively. However, these formal tools have led us to novel insights on current data and to put forward precise hypotheses which can be addressed experimentally. In particular, despite subtle differences in the root and leaf networks, these have equivalent dynamical behaviors. Our simulations also suggest that only when a biasing signal positively affects an activator in the network, the system recovers striped cellular patterns similar to those of root epidermis. We also postulate that cell shape may affect pattern stability in the root. Our results thus support the idea that in this and other cases, contrasting spatial cell patterns and other evolutionary morphogenetic novelties originate from conserved genetic network modules subject to divergent contextual traits.

Nongenetic perturbations, such as environmental change or developmental noise, can induce novel phenotypes. If an induced phenotype appears recurrently and confers a fitness advantage, selection may promote its genetic stabilization. Nongenetic perturbations can thus initiate evolutionary innovation. Genetic variation that is not usually phenotypically visible may play an important role in this process. Populations under stabilizing selection on a phenotype that is robust to mutations can accumulate such variation. After nongenetic perturbations, this variation can produce new phenotypes. We here study the relationship between a phenotype's mutational robustness and a population's potential to generate novel phenotypic variation. To this end, we use a well‐studied model of transcriptional regulation circuits that are important in many evolutionary innovations. We find that phenotypic robustness promotes phenotypic variability in response to nongenetic perturbations, but not in response to mutation. Our work suggests that nongenetic perturbations may initiate innovation more frequently in mutationally robust gene expression traits.