BEST: a novel computational approach for comparing gene expression patterns from early stages of Drosophila melanogaster development

Embryonic gene expression patterns are an indispensable part of modern developmental biology. Currently, investigators must visually inspect numerous images containing embryonic expression patterns to identify spatially similar patterns for inferring potential genetic interactions. The lack of a computational approach to identify pattern similarities is an impediment to advancement in developmental biology research because of the rapidly increasing amount of available embryonic gene expression data. Therefore, we have developed computational approaches to automate the comparison of gene expression patterns contained in images of early stage Drosophila melanogaster embryos (prior to the beginning of germ-band elongation); similarities and differences in gene expression patterns in these early stages have extensive developmental effects. Here we describe a basic expression search tool (BEST) to retrieve best matching expression patterns for a given query expression pattern and a computational device for gene interaction inference using gene expression pattern images and information on the associated genotypes and probes. Analysis of a prototype collection of Drosophila gene expression pattern images is presented to demonstrate the utility of these methods in identifying biologically meaningful matches and inferring gene interactions by direct image content analysis. In particular, the use of BEST searches for gene expression patterns is akin to that of BLAST searches for finding similar sequences. These computational developmental biology methodologies are likely to make the great wealth of embryonic gene expression pattern data easily accessible and to accelerate the discovery of developmental networks.     

Analysis of the expression and function of Wnt-5a and Wnt-5b in developing and regenerating axolotl (Ambystoma mexicanum) limbs

Urodele amphibians are unique adult vertebrates because they are able to regenerate body parts after amputation. Studies of urodele limb regeneration, the key model system for vertebrate regeneration, have led to an understanding of the origin of blastema cells and the importance of positional interactions between blastema cells in the control of growth and pattern formation. Progress is now being made in the identification of the signaling pathways that regulate dedifferentiation, blastema morphogenesis, growth and pattern formation. Members of the Wnt family of secreted proteins are expressed in developing and regenerating limbs, and have the potential to control growth, pattern formation and differentiation. We have studied the expression of two non-canonical Wnt genes, Wnt-5a and Wnt-5b . We report that they are expressed in equivalent patterns during limb development and limb regeneration in the axolotl ( Ambystoma mexicanum ), and during limb development in other tetrapods, implying conservation of function. Our analysis of the effects of ectopic Wnt-5a expression is consistent with the hypothesis that canonical Wnt signaling functions during the early stages of regeneration to control the dedifferentiation of stump cells giving rise to the regeneration-competent cells of the blastema.     

The dominant hemimelia mutation uncouples epithelial-mesenchymal interactions and disrupts anterior mesenchyme formation in mouse hindlimbs.

Epithelial-mesenchymal interactions are essential for both limb outgrowth and pattern formation in the limb. Molecules capable of communication between these two tissues are known and include the signaling molecules SHH and FGF4, FGF8 and FGF10. Evidence suggests that the pattern and maintenance of expression of these genes are dependent on a number of factors including regulatory loops between genes expressed in the AER and those in the underlying mesenchyme. We show here that the mouse mutation dominant hemimelia (Dh) alters the pattern of gene expression in the AER such that Fgf4, which is normally expressed in a posterior domain, and Fgf8, which is expressed throughout are expressed in anterior patterns. We show that maintenance of Shh expression in the posterior mesenchyme is not dependent on either expression of Fgf4 or normal levels of Fgf8 in the overlying AER. Conversely, AER expression of Fgf4 is not directly dependent on Shh expression. Also the reciprocal regulatory loop proposed for Fgf8 in the AER and Fgf10 in the underlying mesenchyme is also uncoupled by this mutation. Early during the process of limb initiation, Dh is involved in regulating the width of the limb bud, the mutation resulting in selective loss of anterior mesenchyme. The Dh gene functions in the initial stages of limb development and we suggest that these initial roles are linked to mechanisms that pattern gene expression in the AER.     

The Dynamics of Turing Patterns for Morphogen-Regulated Growing Domains with Cellular Response Delays

Since its conception in 1952, the Turing paradigm for pattern formation has been the subject of numerous theoretical investigations. Experimentally, this mechanism was first demonstrated in chemical reactions over 20 years ago and, more recently, several examples of biological self-organisation have also been implicated as Turing systems. One criticism of the Turing model is its lack of robustness, not only with respect to fluctuations in the initial conditions, but also with respect to the inclusion of delays in critical feedback processes such as gene expression. In this work we investigate the possibilities for Turing patterns on growing domains where the morphogens additionally regulate domain growth, incorporating delays in the feedback between signalling and domain growth, as well as gene expression. We present results for the proto-typical Schnakenberg and Gierer–Meinhardt systems: exploring the dynamics of these systems suggests a reconsideration of the basic Turing mechanism for pattern formation on morphogen-regulated growing domains as well as highlighting when feedback delays on domain growth are important for pattern formation.     

Molecular anatomy of the developing limb in the coquí frog,Eleutherodactylus coqui

The vertebrate limb demonstrates remarkable similarity in basic organization across phylogenetically disparate groups. To gain further insight into how this morphological similarity is maintained in different developmental contexts, we explored the molecular anatomy of size‐reduced embryos of the Puerto Rican coquí frog, Eleutherodactylus coqui. This animal demonstrates direct development, a life‐history strategy marked by rapid progression from egg to adult and absence of a free‐living, aquatic larva. Nonetheless, coquí exhibits a basal anuran limb structure, with four toes on the forelimb and five toes on the hind limb. We investigated the extent to which coquí limb bud development conforms to the model of limb development derived from amniote studies. Toward this end, we characterized dynamic patterns of expression for 13 critical patterning genes across three principle stages of limb development. As expected, most genes demonstrate expression patterns that are essentially unchanged compared to amniote species. For example, we identified an EcFgf8‐expression domain within the apical ectodermal ridge (AER). This expression pattern defines a putatively functional AER signaling domain, despite the absence of a morphological ridge in coquí embryos. However, two genes, EcMeis2 and EcAlx4, demonstrate altered domains of expression, which imply a potential shift in gene function between coquí frogs and amniote model systems. Unexpectedly, several genes thought to be critical for limb patterning in other systems, including EcFgf4, EcWnt3a, EcWnt7a, and EcGremlin, demonstrated no evident expression pattern in the limb at the three stages we analyzed. The absence of EcFgf4 and EcWnt3a expression during limb patterning is perhaps not surprising, given that neither gene is critical for proper limb development in the mouse, based on knockout and expression analyses. In contrast, absence of EcWnt7a and EcGremlin is surprising, given that expression of these molecules appears to be absolutely essential in all other model systems so far examined. Although this analysis substantiates the existence of a core set of ancient limb‐patterning molecules, which likely mediate identical functions across highly diverse vertebrate forms, it also reveals remarkable evolutionary flexibility in the genetic control of a conserved morphological pattern across evolutionary time.     

Bioinformatics beyond sequence: mapping gene function in the embryo

The spatio-temporal expression pattern of a gene during development is a valuable piece of information. But there is no way to compare precisely the patterns of expression of different genes, or the way the patterns are changed in a mutant. One way to solve this problem is to construct digital reference images of development (a bioinformatics framework), to which expression patterns can be mapped and stored, then compared. Such frameworks are under active development in several model systems. They will form the basis of powerful and integrated gene expression databases, which facilitate comparisons between genes, tissues and species.     

Spatiotemporal mechanisms of morphogen gradient interpretation

Few mechanistic ideas from the pre-molecular era of biology have had as enduring an impact as the morphogen concept. In the classical view, cells in developing embryos obtain positional information by measuring morphogen concentrations and comparing them with fixed concentration thresholds; as a result, graded morphogen distributions map into discrete spatial arrangements of gene expression. Recent studies on Hedgehog and other morphogens suggest that establishing patterns of gene expression may be less a function of absolute morphogen concentrations, than of the dynamics of signal transduction, gene expression, and gradient formation. The data point away from any universal model of morphogen interpretation and suggest that organisms use multiple mechanisms for reading out developmental signals in order to accomplish specific patterning goals.     

A new computational growth model for sea urchin skeletons

A new computational model has been developed to simulate growth of regular sea urchin skeletons. The model incorporates the processes of plate addition and individual plate growth into a composite model of whole-body (somatic) growth. A simple developmental model based on hypothetical morphogens underlies the assumptions used to define the simulated growth processes. The data model is based on a Delaunay triangulation of plate growth center points, using the dual Voronoi polygons to define plate topologies. A spherical frame of reference is used for growth calculations, with affine deformation of the sphere (based on a Young-Laplace membrane model) to result in an urchin-like three-dimensional form. The model verifies that the patterns of coronal plates in general meet the criteria of Voronoi polygonalization, that a morphogen/threshold inhibition model for plate addition results in the alternating plate addition pattern characteristic of sea urchins, and that application of the Bertalanffy growth model to individual plates results in simulated somatic growth that approximates that seen in living urchins. The model suggests avenues of research that could explain some of the distinctions between modern sea urchins and the much more disparate groups of forms that characterized the Paleozoic Era.     

Initiation of dorso-ventral axis during chick limb development

We analysed spatio-temporal expression of dorso-ventral genes - Wnt-7a, En-1, Lmx-1 and Fgf-8 - during both normal and ectopic limb formation following fibroblast growth factor (FGF) application to the flank. We confirm that Wnt-7a is the first of these genes to be expressed in dorsal ectoderm in limb-forming regions. We also noticed patterns and kinetics of gene expression specific to chick that could account for differences observed in ridge formation between chick and mouse. We find that Wnt-7a expression, in dorsal ectoderm, is rapidly and locally induced by FGF application. In contrast, ectopic induction of Lmx-1 expression, in dorsal mesoderm, is much slower, occurs first at a distance from the FGF-2 bead and seems initially independent of direct Wnt-7a signalling during FGF-2 limb induction. Finally, we show that there is no contribution to extra-limb mesoderm from normal limb mesoderm and confirm that flank cells give rise to the extra limb. Furthermore, we suggest that an inhibitor present in the flank normally prevents Lmx-1 expression in this region and restricts its expression to limb-forming regions.     

A Mechanochemical Model for Embryonic Pattern Formation: Coupling Tissue Mechanics and Morphogen Expression

Motivated by recent experimental findings, we propose a novel mechanism of embryonic pattern formation based on coupling of tissue curvature with diffusive signaling by a chemical factor. We derive a new mathematical model using energy minimization approach and show that the model generates a variety of morphogen and curvature patterns agreeing with experimentally observed structures. The mechanism proposed transcends the classical Turing concept which requires interactions between two morphogens with a significantly different diffusivity. Our studies show how biomechanical forces may replace the elusive long-range inhibitor and lead to formation of stable spatially heterogeneous structures without existence of chemical prepatterns. We propose new experimental approaches to decisively test our central hypothesis that tissue curvature and morphogen expression are coupled in a positive feedback loop.     

Modes of morphogen cooperation for limb formation in vertebrates and insects.

Diffusing morphogens in cooperation can control gene expression in developing limbs. Additive cooperation corresponds to the Boolean operator OR and implies the equivalent action of the (suitably scaled) concentrations of two morphogens, either by their alternative binding to the same receptor or by another way of convergence of their effects during the signal transduction procedure. This cooperation can explain the spatial and temporal collinearities of the expression of hoxd genes in the vertebrate limb bud. A multiplicative cooperation of morphogens (corresponding to the Boolean operator AND), produced at the DPP and WG domains in the Drosophila leg imaginal disc, may account for the expression domains observed for Dll and dac. A molecular interpretation of the multiplicative morphogen cooperation is proposed. Some experiments are suggested for further testing of the model.     

The role and limits of a gradient based explanation of morphogenesis: a theoretical consideration

Development of an organism is a multi-dimensional process leading to the generation of complex species-specific structures. This specificity suggests machine-like organisation. The uneven distribution (gradient) of soluble substances (morphogens) and specific receptor-ligand interactions are known to cause differential gene expression. Therefore gradients of morphogens are used as a causal explanation of developmental processes. However each attempt to describe development causally should take into account both the local fine organisation and global robustness of morphogenesis. The classical view of the role of morphogens will be critically considered and possible alternative proposed. The core idea of my proposal is that the main function of morphogenetic substances could be a context dependent modification of cell behaviour. Both history and different features of morphogenetic fields create the framework for the activity of morphogenes.