Positional information in developing embryos is specified by spatial gradients of transcriptional regulators. One of the classic systems for studying this is the activation of the hunchback (hb) gene in early fruit fly (Drosophila) segmentation by the maternally-derived gradient of the Bicoid (Bcd) protein. Gene regulation is subject to intrinsic noise which can produce variable expression. This variability must be constrained in the highly reproducible and coordinated events of development. We identify means by which noise is controlled during gene expression by characterizing the dependence of hb mRNA and protein output noise on hb promoter structure and transcriptional dynamics. We use a stochastic model of the hb promoter in which the number and strength of Bcd and Hb (self-regulatory) binding sites can be varied. Model parameters are fit to data from WT embryos, the self-regulation mutant hb^14F, and lacZ reporter constructs using different portions of the hb promoter. We have corroborated model noise predictions experimentally. The results indicate that WT (self-regulatory) Hb output noise is predominantly dependent on the transcription and translation dynamics of its own expression, rather than on Bcd fluctuations. The constructs and mutant, which lack self-regulation, indicate that the multiple Bcd binding sites in the hb promoter (and their strengths) also play a role in buffering noise. The model is robust to the variation in Bcd binding site number across a number of fly species. This study identifies particular ways in which promoter structure and regulatory dynamics reduce hb output noise. Insofar as many of these are common features of genes (e.g. multiple regulatory sites, cooperativity, self-feedback), the current results contribute to the general understanding of the reproducibility and determinacy of spatial patterning in early development. Noise is an intrinsic part of biochemical systems such as gene regulation networks. Noisy gene expression has been well documented in populations of single cells, and is likely a key mechanism in evolutionary change. But in developing embryos, cells within a tissue must overcome such variability in order to provide the uniformity required to coordinate multiple events. Reproducibility and determinacy of the spatial protein patterns preceding tissue differentiation is a critical aspect of development. In this study, we use anterior-posterior (AP) segmentation in the fruit fly (Drosophila) to understand how gene regulation dynamics control noise. One of the earliest AP patterning events is the anterior activation of the hunchback (hb) gene by the maternally-derived Bicoid (Bcd) protein gradient. This interaction has been very well characterized, providing the tools for us to develop a stochastic model of hb gene regulation to make predictions about expression noise, and to corroborate these experimentally. For hb, we show that self-regulation is a critical part of controlling noise, and the multiple Bcd binding sites in the hb promoter also enhance pattern reproducibility. To the degree that such features are shared by other genes, these noise-reducing mechanisms may be common to many pattern forming events.

During embryonic development, the positional information provided by concentration gradients of maternal factors directs pattern formation by providing spatially dependent cues for gene expression. In the fruit fly, Drosophila melanogaster, a classic example of this is the sharp on–off activation of the hunchback (hb) gene at midembryo, in response to local concentrations of the smooth anterior–posterior Bicoid (Bcd) gradient. The regulatory region for hb contains multiple binding sites for the Bcd protein as well as multiple binding sites for the Hb protein. Some previous studies have suggested that Bcd is sufficient for properly sharpened Hb expression, yet other evidence suggests a need for additional regulation. We experimentally quantified the dynamics of hb gene expression in flies that were wild-type, were mutant for hb self-regulation or Bcd binding, or contained an artificial promoter construct consisting of six Bcd and two Hb sites. In addition to these experiments, we developed a reaction–diffusion model of hb transcription, with Bcd cooperative binding and hb self-regulation, and used Zero Eigenvalue Analysis to look for multiple stationary states in the reaction network. Our model reproduces the hb developmental dynamics and correctly predicts the mutant patterns. Analysis of our model indicates that the Hb sharpness can be produced by spatial bistability, in which hb self-regulation produces two stable levels of expression. In the absence of self-regulation, the bistable behavior vanishes and Hb sharpness is disrupted. Bcd cooperative binding affects the position where bistability occurs but is not itself sufficient for a sharp Hb pattern. Our results show that the control of Hb sharpness and positioning, by hb self-regulation and Bcd cooperativity, respectively, are separate processes that can be altered independently. Our model, which matches the changes in Hb position and sharpness observed in different experiments, provides a theoretical framework for understanding the data and in particular indicates that spatial bistability can play a central role in threshold-dependent reading mechanisms of positional information. Pattern formation during embryonic development, or morphogenesis, is one of the most intriguing problems in biology, entailing the sequence of processes by which a relatively simple system, the fertilized egg, becomes a mature organism. In these processes, the genetic information, stored at the molecular scale in the DNA, is translated into the macroscopic spatial expression patterns that precede the tissue–organ scale of body organization. It can also be understood as a flux of information from the genetic to the organ–system level. In the fruit fly Drosophila melanogaster, one of the early processes during its embryonic development is the formation of the sharp Hunchback protein pattern. To generate this pattern, the hunchback gene interprets the position-dependent information in the shallow maternal Bicoid gradient and converts it into the sharp Hunchback protein pattern. We propose that bistability in the dynamics of hunchback gene regulation can account for this information reading process, and we show that this bistable mechanism can be produced by the ability of this gene to regulate its own expression. The solution of this problem offers new approaches to understand the phenomenon of morphogenesis.

We analyze the relation between maternal gradients and segmentation in Drosophila, by quantifying spatial precision in protein patterns. Segmentation is first seen in the striped expression patterns of the pair‐rule genes, such as even‐skipped (eve). We compare positional precision between Eve and the maternal gradients of Bicoid (Bcd) and Caudal (Cad) proteins, showing that Eve position could be initially specified by the maternal protein concentrations but that these do not have the precision to specify the mature striped pattern of Eve. By using spatial trends, we avoid possible complications in measuring single boundary precision (e.g., gap gene patterns) and can follow how precision changes in time. During nuclear cleavage cycles 13 and 14, we find that Eve becomes increasingly correlated with egg length, whereas Bcd does not. This finding suggests that the change in precision is part of a separation of segmentation from an absolute spatial measure, established by the maternal gradients, to one precise in relative (percent egg length) units. Developmental Dynamics 235:2949–2960, 2006.

A study is made by computation of the interplay between the pattern formation of growth catalysts on a plant surface and the expansion of the surface to generate organismal shape. Consideration is made of the localization of morphogenetically active regions, and the occurrence within them of symmetry-breaking processes such as branching from an initially dome-shaped tip or meristem. Representation of a changing and growing three-dimensional shape is necessary, as two-dimensional work cannot distinguish, for example, formation of an annulus from dichotomous branching. For the formation of patterns of chemical concentrations, the Brusselator reaction-diffusion model is used, applied on a hemispherical shell and generating patterns that initiate as surface spherical harmonics. The initial shape is hemispherical, represented as a mesh of triangles. These are combined into finite elements, each made up of all the triangles surrounding each node. Chemical pattern is converted into shape change by moving nodes outwards according to the concentration of growth catalyst at each, to relieve misfits caused by area increase of the finite element. New triangles are added to restore the refinement of the mesh in rapidly growing regions. The postulated mechanism successfully generates: tip growth (or stalk extension by an apical meristem) to ten times original hemisphere height; tip flattening and resumption of apical advance; and dichotomous branching and higher-order branching to make whorled structures. Control of the branching plane in successive dichotomous branchings is tackled with partial success and clarification of the issues. The representation of a growing plant surface in computations by an expanding mesh that has no artefacts constraining changes of shape and symmetry has been achieved. It is shown that one type of pattern-forming mechanism, Turing-type reaction-diffusion, acting within a surface to pattern a growth catalyst, can generate some of the most important types of morphogenesis in plant development.

Reaction–diffusion theory for pattern formation is considered in relation to processes of biological development in which there is continuous growth and shape change as each new pattern forms. This is particularly common in the plant kingdom, for both unicellular and multicellular organisms. In addition to the feedbacks in the chemical dynamics, there is then another loop linking size and shape changes with the reaction–diffusion patterning of growth controllers in the growing region. In studies by computation, the codes must incorporate, alongside the usual solvers of the partial differential dynamic equations, a versatile growth code, to express any kind of shape change. We have found that regulation of shape change in particular ways (e.g. to make narrow-angle branchings) demands new features in our chemical mechanisms. Our growth algorithm is for a surface growing tangentially, but moving outward and changing shape to accommodate the extra area. This is potentially applicable both to the tunica layer of multicellular plant meristems and to the growing tip of the cell surface, e.g. in the morphogenesis of single-celled chlorophyte algae which display branching processes: whorl formation in Acetabularia (Dasycladales) and repeated dichotomous branching in Micrasterias (Desmidiaceae). For computational studies, a hemispherical shell is a reasonable idealization of the initial shape. We describe results of two types of study: (1) Pattern formation by three reaction–diffusion models, with contrasted nonlinearities, on the hemispherical shell, particularly to find conditions for robust formation of annular pattern or pattern for dichotomous branching, both of which are common in plants. (2) Sequential dichotomous branchings in a system growing and changing in shape from the hemispherical start.