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.

Here we characterize the expression of the full system of genes which control the segmentation morphogenetic field of Drosophila at the protein level in one dimension. The data used for this characterization are quantitative with cellular resolution in space and about 6 min in time. We present the full quantitative profiles of all 14 segmentation genes which act before the onset of gastrulation. The expression patterns of these genes are first characterized in terms of their average or typical behavior. At this level, the expression of all of the genes has been integrated into a single atlas of gene expression in which the expression levels of all genes in each cell are specified. We show that expression domains do not arise synchronously, but rather each domain has its own specific dynamics of formation. Moreover, we show that the expression domains shift position in the direction of the cephalic furrow, such that domains in the anlage of the segmented germ band shift anteriorly while those in the presumptive head shift posteriorly. The expression atlas of integrated data is very close to the expression profiles of individual embryos during the latter part of the blastoderm stage. At earlier times gap gene domains show considerable variation in amplitude, and significant positional variability. Nevertheless, an average early gap domain is close to that of a median individual. In contrast, we show that there is a diversity of developmental trajectories among pair-rule genes at a variety of levels, including the order of domain formation and positional accuracy. We further show that this variation is dynamically reduced, or canalized, over time. As the first quantitatively characterized morphogenetic field, this system and its behavior constitute an extraordinarily rich set of materials for the study of canalization and embryonic regulation at the molecular level.

Genetic studies have revealed that segment determination in Drosophila melanogaster is based on hierarchical regulatory interactions among maternal coordinate and zygotic segmentation genes. The gap gene system constitutes the most upstream zygotic layer of this regulatory hierarchy, responsible for the initial interpretation of positional information encoded by maternal gradients. We present a detailed analysis of regulatory interactions involved in gap gene regulation based on gap gene circuits, which are mathematical gene network models used to infer regulatory interactions from quantitative gene expression data. Our models reproduce gap gene expression at high accuracy and temporal resolution. Regulatory interactions found in gap gene circuits provide consistent and sufficient mechanisms for gap gene expression, which largely agree with mechanisms previously inferred from qualitative studies of mutant gene expression patterns. Our models predict activation of Kr by Cad and clarify several other regulatory interactions. Our analysis suggests a central role for repressive feedback loops between complementary gap genes. We observe that repressive interactions among overlapping gap genes show anteroposterior asymmetry with posterior dominance. Finally, our models suggest a correlation between timing of gap domain boundary formation and regulatory contributions from the terminal maternal system.

rhomboid (rho) encodes a putative transmembrane receptor that is required for the differentiation of the ventral epidermis. It is initially expressed before the completion of cellularization in lateral stripes within the presumptive neuroectoderm. Here, we present evidence that the maternal morphogen dorsal (dl) acts in concert with basic helix-loop-helix (b-HLH) proteins, possibly including twist (twi), to activate rho in both lateral and ventral regions. Expression is blocked in ventral regions (the presumptive mesoderm) by sna/l (sna), which is also a direct target of the d/morphogen. A 300-bp region of the rho promoter (the NEE), which is sufficient for neuroectoderm expression, contains a cluster of d/and b-HLH activator sites that are closely linked to sna repressor sites. Mutations in these binding sites cause genetically predicted changes in the levels and limits of rho expression. In particular, the disruption of sna-binding sites causes a derepression of the pattern throughout ventral regions, providing evidence that sna is directly responsible for establishing the mesoderm/neuroectoderm boundary before gastrulation. The tight linkage of activator and repressor sites in the rho NEE is similar to the arrangement of binding sites observed in the even-skipped stripe 2 element, which is regulated by bicoid (bcal. This suggests that the d/and bcd morphogens use a similar mechanism to make stripes in the Drosophila embryo.

We present a fluorescence-based, multiplex in situ hybridization method that permits the simultaneous detection of five differently labeled antisense RNA probes and up to seven differ-ent transcripts in a single Drosophila embryo. We also show that it should be possible to increase the number of detected transcripts substantially with nascent transcript multiplex fluorescent in situ hybridization. These multiplex methods fill a current technological gap between high-resolution in situ hybridization with one or two fluorescently labeled probes and low-resolution but genome-wide microarray RNA profiling and should be of great utility in establishing gene networks.