Flower Development as an Interplay between Dynamical Physical Fields and Genetic Networks

In this paper we propose a model to describe the mechanisms by which undifferentiated cells attain gene configurations underlying cell fate determination during morphogenesis. Despite the complicated mechanisms that surely intervene in this process, it is clear that the fundamental fact is that cells obtain spatial and temporal information that bias their destiny. Our main hypothesis assumes that there is at least one macroscopic field that breaks the symmetry of space at a given time. This field provides the information required for the process of cell differentiation to occur by being dynamically coupled to a signal transduction mechanism that, in turn, acts directly upon the gene regulatory network (GRN) underlying cell-fate decisions within cells. We illustrate and test our proposal with a GRN model grounded on experimental data for cell fate specification during organ formation in early Arabidopsis thaliana flower development. We show that our model is able to recover the multigene configurations characteristic of sepal, petal, stamen and carpel primordial cells arranged in concentric rings, in a similar pattern to that observed during actual floral organ determination. Such pattern is robust to alterations of the model parameters and simulated failures predict altered spatio-temporal patterns that mimic those described for several mutants. Furthermore, simulated alterations in the physical fields predict a pattern equivalent to that found in Lacandonia schismatica, the only flowering species with central stamens surrounded by carpels

Modeling Stochasticity and Variability in Gene Regulatory Networks

Modeling stochasticity in gene regulatory networks is an important and complex problem in molecular systems biology. To elucidate intrinsic noise, several modeling strategies such as the Gillespie algorithm have been used successfully. This paper contributes an approach as an alternative to these classical settings. Within the discrete paradigm, where genes, proteins, and other molecular components of gene regulatory networks are modeled as discrete variables and are assigned as logical rules describing their regulation through interactions with other components. Stochasticity is modeled at the biological function level under the assumption that even if the expression levels of the input nodes of an update rule guarantee activation or degradation there is a probability that the process will not occur due to stochastic effects. This approach allows a finer analysis of discrete models and provides a natural setup for cell population simulations to study cell-to-cell variability. We applied our methods to two of the most studied regulatory networks, the outcome of lambda phage infection of bacteria and the p53-mdm2 complex.Comment: 23 pages, 8 figure

Positive selection on the K domain of the AGAMOUS protein in the Zingiberales suggests a mechanism for the evolution of androecial morphology

BACKGROUND: The ABC model of flower development describes the molecular basis for specification of floral organ identity in model eudicots such as Arabidopsis and Antirrhinum. According to this model, expression of C-class genes is linked to stamen and gynoecium organ identity. The Zingiberales is an order of tropical monocots in which the evolution of floral morphology is characterized by a marked increase in petaloidy in the androecium. Petaloidy is a derived characteristic of the ginger families and seems to have arisen in the common ancestor of the ginger clade. We hypothesize that duplication of the C-class AGAMOUS (AG) gene followed by divergence of the duplicated AG copies during the diversification of the ginger clade lineages explains the evolution of petaloidy in the androecium. In order to address this hypothesis, we carried out phylogenetic analyses of the AG gene family across the Zingiberales and investigated patterns of gene expression within the androecium. RESULTS: Phylogenetic analysis supports a scenario in which Zingiberales-specific AG genes have undergone at least one round of duplication. Gene duplication was immediately followed by divergence of the retained copies. In particular, we detect positive selection in the third alpha-helix of the K domain of Zingiberales AGAMOUS copy 1 (ZinAG-1). A single fixed amino acid change is observed in ZinAG-1 within the ginger clade when compared to the banana grade. Expression analyses of AG and APETALA1/FRUITFULL (AP1/FUL) in Musa basjoo is similar to A- and C-class gene expressions in the Arabidopsis thaliana model, while Costus spicatus exhibits simultaneous expression of AG and AP1/FUL in most floral organs. We propose that this novel expression pattern could be correlated with the evolution of androecial petaloidy within the Zingiberales. CONCLUSIONS: Our results present an intricate story in which duplication of the AG lineage has lead to the retention of at least two diverged Zingiberales-specific copies, ZinAG-1 and Zingiberales AGAMOUS copy 2 (ZinAG-2). Positive selection on ZinAG-1 residues suggests a mechanism by which AG gene divergence may explain observed morphological changes in Zingiberales flowers. Expression data provides preliminary support for the proposed mechanism, although further studies are required to fully test this hypothesis