The segmentation gene network in Drosophila embryo solves the fundamental problem of embryonic patterning: how to establish a periodic pattern of gene expression, which determines both the positions and the identities of body segments. The gap gene network constitutes the first zygotic regulatory tier in this process. Here we have applied the systems-level approach to investigate the regulatory effect of gap gene Kruppel (Kr) on segmentation gene expression. We acquired a large dataset on the expression of gap genes in Kr null mutants and demonstrated that the expression levels of these genes are significantly reduced in the second half of cycle 14A. To explain this novel biological result we applied the gene circuit method which extracts regulatory information from spatial gene expression data. Previous attempts to use this formalism to correctly and quantitatively reproduce gap gene expression in mutants for a trunk gap gene failed, therefore here we constructed a revised model and showed that it correctly reproduces the expression patterns of gap genes in Kr null mutants. We found that the remarkable alteration of gap gene expression patterns in Kr mutants can be explained by the dynamic decrease of activating effect of Cad on a target gene and exclusion of Kr gene from the complex network of gap gene interactions, that makes it possible for other interactions, in particular, between hb and gt, to come into effect. The successful modeling of the quantitative aspects of gap gene expression in mutant for the trunk gap gene Kr is a significant achievement of this work. This result also clearly indicates that the oversimplified representation of transcriptional regulation in the previous models is one of the reasons for unsuccessful attempts of mutant simulations.

Abstract Motivation: To construct an integrated map of Drosophila segmentation gene expression from partial data taken from individual embryos. Results: Spline and wavelet based registration techniques were developed to register Drosophila segmentation gene expression data. As ground control points for registration we used the locations of extrema on gene expression patterns, represented in 1D. The registration method was characterized by unprecedented high accuracy. A method for constructing the integrated pattern of gene expression at cellular resolution was designed. These patterns were constructed for 9 segmentation genes belonging to gap and pair-rule classes. Availability: Programs running on SPARC/SunOS 5.0 platform or Alpha/DU 4.0E systems may be obtained from the authors. Contact: [email protected] * To whom correspondence should be addressed.

A 1D wavelet-based registration technique was applied to Drosophila segmentation gene expression data. We apply the fast redundant dyadic wavelet transform with a continuous basis as a method for extracting the ground control points. This transform was chosen for its basic properties, such as noise reduction and good localization of characteristic features.