Carlos E. Vanario-Alonso

Dynamic control of positional information in the early Drosophila embryo

Morphogen gradients contribute to pattern formation by determining positional information in morphogenetic fields. Interpretation of positional information is thought to rely on direct, concentration-threshold-dependent mechanisms for establishing multiple differential domains of target gene expression. In Drosophila, maternal gradients establish the initial position of boundaries for zygotic gap gene expression, which in turn convey positional information to pair-rule and segment-polarity genes, the latter forming a segmental pre-pattern by the onset of gastrulation. Here we report, on the basis of quantitative gene expression data, substantial anterior shifts in the position of gap domains after their initial establishment. Using a data-driven mathematical modelling approach, we show that these shifts are based on a regulatory mechanism that relies on asymmetric gap–gap cross-repression and does not require the diffusion of gap proteins. Our analysis implies that the threshold-dependent interpretation of maternal morphogen concentration is not sufficient to determine shifting gap domain boundary positions, and suggests that establishing and interpreting positional information are not independent processes in the Drosophila blastoderm.

Canalization of gene expression in the Drosophila blastoderm by gap gene cross regulation.

Developing embryos exhibit a robust capability to reduce phenotypic variations that occur naturally or as a result of experimental manipulation. This reduction in variation occurs by an epigenetic mechanism called canalization, a phenomenon which has resisted understanding because of a lack of necessary molecular data and of appropriate gene regulation models. In recent years, quantitative gene expression data have become available for the segment determination process in the Drosophila blastoderm, revealing a specific instance of canalization. These data show that the variation of the zygotic segmentation gene expression patterns is markedly reduced compared to earlier levels by the time gastrulation begins, and this variation is significantly lower than the variation of the maternal protein gradient Bicoid. We used a predictive dynamical model of gene regulation to study the effect of Bicoid variation on the downstream gap genes. The model correctly predicts the reduced variation of the gap gene expression patterns and allows the characterization of the canalizing mechanism. We show that the canalization is the result of specific regulatory interactions among the zygotic gap genes. We demonstrate the validity of this explanation by showing that variation is increased in embryos mutant for two gap genes, Krüppel and knirps, disproving competing proposals that canalization is due to an undiscovered morphogen, or that it does not take place at all. In an accompanying article in PLoS Computational Biology (doi:10.1371/journal.pcbi.1000303), we show that cross regulation between the gap genes causes their expression to approach dynamical attractors, reducing initial variation and providing a robust output. These results demonstrate that the Bicoid gradient is not sufficient to produce gap gene borders having the low variance observed, and instead this low variance is generated by gap gene cross regulation. More generally, we show that the complex multigenic phenomenon of canalization can be understood at a quantitative and predictive level by the application of a precise dynamical model.

Canalization of gene expression and domain shifts in the Drosophila blastoderm by dynamical attractors.

The variation in the expression patterns of the gap genes in the blastoderm of the fruit fly Drosophila melanogaster reduces over time as a result of cross regulation between these genes, a fact that we have demonstrated in an accompanying article in PLoS Biology (see Manu et al., doi:10.1371/journal.pbio.1000049). This biologically essential process is an example of the phenomenon known as canalization. It has been suggested that the developmental trajectory of a wild-type organism is inherently stable, and that canalization is a manifestation of this property. Although the role of gap genes in the canalization process was established by correctly predicting the response of the system to particular perturbations, the stability of the developmental trajectory remains to be investigated. For many years, it has been speculated that stability against perturbations during development can be described by dynamical systems having attracting sets that drive reductions of volume in phase space. In this paper, we show that both the reduction in variability of gap gene expression as well as shifts in the position of posterior gap gene domains are the result of the actions of attractors in the gap gene dynamical system. Two biologically distinct dynamical regions exist in the early embryo, separated by a bifurcation at 53% egg length. In the anterior region, reduction in variation occurs because of stability induced by point attractors, while in the posterior, the stability of the developmental trajectory arises from a one-dimensional attracting manifold. This manifold also controls a previously characterized anterior shift of posterior region gap domains. Our analysis shows that the complex phenomena of canalization and pattern formation in the Drosophila blastoderm can be understood in terms of the qualitative features of the dynamical system. The result confirms the idea that attractors are important for developmental stability and shows a richer variety of dynamical attractors in developmental systems than has been previously recognized.

Stripe forming architecture of the gap gene system.

In this report, we show that gap genes encode exactly one set of pair-rule stripes, which occur in the native even-skipped position. The core of this work is a detailed analysis that shows how this conclusion follows from the arrangement of gap domains in the embryo. This analysis shows that: (1) pattern forming information is transmitted from gap to pair-rule genes by means of a nonredundant set of morphogenetic gradients, and (2) the stripe forming capability of the gap genes is constrained by the arrangement of these gradients and by the fact that each gap domain consists of a pair of correlated gradients. We also show that in the blastoderm, the regulatory sign of a transcriptional regulator is unlikely to change in a concentration dependent manner. The principal analytic tool used to establish these results is the gene circuit method. Here, this method is applied to examine hybrid data sets consisting of real gene expression data for four gap genes and hypothetical pair-rule expression data generated by translating native even-skipped data along the anterior-posterior axis. In this way, we are able to investigate the stripe forming capabilities of the gap gene system in the complete absence of pair-rule cross regulation. We close with an inference about evolutionary development. It is argued that the constraints on gap gene architecture identified here are a consequence of selective pressures that minimize the number of gap genes required to determine segments in long-germ band insects.

A high-throughput method for quantifying gene expression data from early Drosophila embryos

We describe an automated high-throughput method to measure protein levels in single nuclei in blastoderm embryos of Drosophila melanogaster by means of immunofluorescence. The method consists of a chain of specific algorithms assembled into an image processing pipeline. This pipeline transforms a confocal scan of an embryo stained with fluorescently tagged antibodies into a text file. This text file contains a numerical identifier for each nucleus, the coordinates of its centroid, and the average concentrations of three proteins in that nucleus. The central algorithmic component of the method is the automatic identification of nuclei by edge detection with the use of watersheds as an error-correction step. This method provides high-throughput quantification at cellular resolution.

Sharp borders from fuzzy gradients

Critical boundaries in the early Drosophila embryo are set by morphogenetic gradients. A new quantitative study shows that the placement of one such boundary is more accurate than the gradient thought to set it. Genetic analysis of the accuracy of the process implicates a gene not previously thought to be involved.

Specialized features of Rhynchosciara americana embryogenesis

Insect embryo development is a complex process which requires nuclear and cellular division, cell shape alteration, and cell movement. This process needs to be orchestrated in a specific spatial and temporal fashion. Different insect species, despite similarities, present distinct morphogenetic pathways. We used the dipteran R. americana as a comparative model for embryo morphogenesis studies, following embryo development with different histochemical and immunohistochemical procedures. Despite the phylogenetic proximity with D. melanogaster, R. americana presents a peculiar morphogenesis. We show that at the initial phases of development, from egg fertilization to blastoderm formation, R. americana is similar to Drosophila. The first cleavages are nuclear and cellularization only begins after nuclei spread throughout the egg’s cortex. However after this stage a series of cell movements establishes a short compact germ band anlage, which gastrulates in a pattern quite different from Drosophila. After gastrulation the germ band elongates anterior–posteriorly and segmentation occurs simultaneously along the embryo. Embryo development from egg fertilization to larva hatching takes about 12 days. Our results show that R. americana presents a different morphogenetic pathway which does not fit in the current short, intermediate or long germ band classification.

Expression of dorsal-ventral genes during early development of Rhynchosciara americana embryos

The establishment of dorsal-ventral polarity in Drosophila is a complex process which involves the action of maternal and zygotically expressed genes. Interspecific differences in the expression pattern of some of these genes have been described in other species. Here we present the expression of dorsal-ventral genes during early embryogenesis in the lower dipteran Rhynchosciara americana. The expression of four genes, the ventralizing genes snail (sna) and twist (twi) and the dorsalizing genes decapentaplegic (dpp) and zerknullt (zen), was investigated by whole-mount in situ hybridization. Sense and antisense mRNA were transcribed in vitro using UTP-digoxigenin and hybridized at 55 degrees C with dechorionated fixed embryos. Staining was obtained with anti-digoxigenin alkaline phosphatase-conjugated antibody revealed with NBT-BCIP solution. The results showed that, in general, the spatial-temporal expression of R. americana dorsal-ventral genes is similar to that observed in Drosophila, where twi and sna are restricted to the ventral region, while dpp and zen are expressed in the dorsal side. The differences encountered were subtle and probably represent a particular aspect of dorsal-ventral axis determination in R. americana. In this lower dipteran sna is expressed slightly later than twi and dpp expression is expanded over the lateral ectoderm during cellular blastoderm stage. These data suggest that the establishment of dorsal-ventral polarity in R. americana embryos follows a program similar to that observed in Drosophila melanogaster.

Embryonic expression of the engrailed homologue of Rhynchosciara americana

The segment polarity gene engrailed is involved in the determination of segment posterior identity in Drosophila. engrailed has been largely used for comparative developmental studies due to its evolutionary conservation from nematodes to humans. By in situ hybridization of an engrailed cDNA probe from Drosophila to polytene chromosomes of fourth instar larvae of Rhynchosciara americana we have shown that engrailed-like sequences must be localized in band 6 of chromosome A in this species. The pattern of engrailed protein expression during R. americana embryo development is diffuse at first evolving into a nuclear striped pattern after quite a length of time. In addition, our results suggest a possible developmentally regulated molecular modification of engrailed protein in R. americana embryos.

Experimental determination of intrinsic drosophila embryo coordinates by evolutionary computation

Early fruit fly embryo development begins with the formation of a chemical blueprint that guides cellular movements and the development of organs and tissues. This blueprint sets the intrinsic spatial coordinates of the embryo. The coordinates are curvilinear from the start, becoming more curvilinear as cells start coherent movements several hours into development. This dynamic aspect of the curvature is an important characteristic of early embryogenesis: characterizing it is crucial for quantitative analysis and dynamic modeling of development. This presents a number of methodological problems for the elastic deformation of 3D and 4D data from confocal microscopy, to standardize images and follow temporal changes. The parameter searches for these deformations present hard optimization problems. Here we describe our evolutionary computation approaches to these problems. We outline some of the immediate applications of these techniques to crucial problems in Drosophila developmental biology.