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

A dual role for DNA binding by Runt in activation and repression of sloppy paired transcription

This work investigates the role of DNA binding by Runt in regulating the sloppy paired 1 (slp1) gene and in particular two distinct cis-regulatory elements that mediate regulation by Runt and other pair-rule transcription factors during Drosophila segmentation. We find that a DNA-binding-defective form of Runt is ineffective at repressing both the distal (DESE) and proximal (PESE) early stripe elements of slp1 and is also compromised for DESE-dependent activation. The function of Runt-binding sites in DESE is further investigated using site-specific transgenesis and quantitative imaging techniques. When DESE is tested as an autonomous enhancer, mutagenesis of the Runt sites results in a clear loss of Runt-dependent repression but has little to no effect on Runt-dependent activation. Notably, mutagenesis of these same sites in the context of a reporter gene construct that also contains the PESE enhancer results in a significant reduction of DESE-dependent activation as well as the loss of repression observed for the autonomous mutant DESE enhancer. These results provide strong evidence that DNA binding by Runt directly contributes to the regulatory interplay of interactions between these two enhancers in the early embryo.

Fully interpretable deep learning model of transcriptional control

Motivation

The universal expressibility assumption of Deep Neural Networks (DNNs) is the key motivation behind recent worksin the systems biology community to employDNNs to solve important problems in functional genomics and moleculargenetics. Typically, such investigations have taken a ‘black box’ approach in which the internal structure of themodel used is set purely by machine learning considerations with little consideration of representing the internalstructure of the biological system by the mathematical structure of the DNN. DNNs have not yet been applied to thedetailed modeling of transcriptional control in which mRNA production is controlled by the binding of specific transcriptionfactors to DNA, in part because such models are in part formulated in terms of specific chemical equationsthat appear different in form from those used in neural networks.

Results

In this paper, we give an example of a DNN whichcan model the detailed control of transcription in a precise and predictive manner. Its internal structure is fully interpretableand is faithful to underlying chemistry of transcription factor binding to DNA. We derive our DNN from asystems biology model that was not previously recognized as having a DNN structure. Although we apply our DNNto data from the early embryo of the fruit fly Drosophila, this system serves as a test bed for analysis of much larger datasets obtained by systems biology studies on a genomic scale. .

Availability and implementation

The implementation and data for the models used in this paper are in a zip file in the supplementary material.

Supplementary information

Supplementary data are available at Bioinformatics online.

Classification-Based Inference of Dynamical Models of Gene Regulatory Networks

Cell-fate decisions during development are controlled by densely interconnected gene regulatory networks (GRNs) consisting of many genes. Inferring and predictively modeling these GRNs is crucial for understanding development and other physiological processes. Gene circuits, coupled differential equations that represent gene product synthesis with a switch-like function, provide a biologically realistic framework for modeling the time evolution of gene expression. However, their use has been limited to smaller networks due to the computational expense of inferring model parameters from gene expression data using global non-linear optimization. Here we show that the switch-like nature of gene regulation can be exploited to break the gene circuit inference problem into two simpler optimization problems that are amenable to computationally efficient supervised learning techniques. We present FIGR (Fast Inference of Gene Regulation), a novel classification-based inference approach to determining gene circuit parameters. We demonstrate FIGR’s effectiveness on synthetic data generated from random gene circuits of up to 50 genes as well as experimental data from the gap gene system of Drosophila melanogaster, a benchmark for inferring dynamical GRN models. FIGR is faster than global non-linear optimization by a factor of 600 and its computational complexity scales much better with GRN size. On a practical level, FIGR can accurately infer the biologically realistic gap gene network in under a minute on desktop-class hardware instead of requiring hours of parallel computing. We anticipate that FIGR would enable the inference of much larger biologically realistic GRNs than was possible before.

Synthetic enhancer design by in silico compensatory evolution reveals flexibility and constraint in cis-regulation

BACKGROUND

Models that incorporate specific chemical mechanisms have been successful in describing the activity of Drosophila developmental enhancers as a function of underlying transcription factor binding motifs. Despite this, the minimum set of mechanisms required to reconstruct an enhancer from its constituent parts is not known. Synthetic biology offers the potential to test the sufficiency of known mechanisms to describe the activity of enhancers, as well as to uncover constraints on the number, order, and spacing of motifs.

RESULTS

Using a functional model and in silico compensatory evolution, we generated putative synthetic even-skipped stripe 2 enhancers with varying degrees of similarity to the natural enhancer. These elements represent the evolutionary trajectories of the natural stripe 2 enhancer towards two synthetic enhancers designed ab initio. In the first trajectory, spatially regulated expression was maintained, even after more than a third of binding sites were lost. In the second, sequences with high similarity to the natural element did not drive expression, but a highly diverged sequence about half the length of the minimal stripe 2 enhancer drove ten times greater expression. Additionally, homotypic clusters of Zelda or Stat92E motifs, but not Bicoid, drove expression in developing embryos.

CONCLUSIONS

Here, we present a functional model of gene regulation to test the degree to which the known transcription factors and their interactions explain the activity of the Drosophila even-skipped stripe 2 enhancer. Initial success in the first trajectory showed that the gene regulation model explains much of the function of the stripe 2 enhancer. Cases where expression deviated from prediction indicates that undescribed factors likely act to modulate expression. We also showed that activation driven Bicoid and Hunchback is highly sensitive to spatial organization of binding motifs. In contrast, Zelda and Stat92E drive expression from simple homotypic clusters, suggesting that activation driven by these factors is less constrained. Collectively, the 40 sequences generated in this work provides a powerful training set for building future models of gene regulation.

A sequence level model of an intact locus predicts the location and function of nonadditive enhancers

Metazoan gene expression is controlled through the action of long stretches of noncoding DNA that contain enhancers-shorter sequences responsible for controlling a single aspect of a gene's expression pattern. Models built on thermodynamics have shown how enhancers interpret protein concentration in order to determine specific levels of gene expression, but the emergent regulatory logic of a complete regulatory locus shows qualitative and quantitative differences from isolated enhancers. Such differences may arise from steric competition limiting the quantity of DNA that can simultaneously influence the transcription machinery. We incorporated this competition into a mechanistic model of gene regulation, generated efficient algorithms for this computation, and applied it to the regulation of Drosophila even-skipped (eve). This model finds the location of enhancers and identifies which factors control the boundaries of eve expression. This model predicts a new enhancer that, when assayed in vivo, drives expression in a non-eve pattern. Incorporation of chromatin accessibility eliminates this inconsistency.

Quantitative perturbation-based analysis of gene expression predicts enhancer activity in early Drosophila embryo

Enhancers constitute one of the major components of regulatory machinery of metazoans. Although several genome-wide studies have focused on finding and locating enhancers in the genomes, the fundamental principles governing their internal architecture and cis-regulatory grammar remain elusive. Here, we describe an extensive, quantitative perturbation analysis targeting the dorsal-ventral patterning gene regulatory network (GRN) controlled by Drosophila NF-κB homolog Dorsal. To understand transcription factor interactions on enhancers, we employed an ensemble of mathematical models, testing effects of cooperativity, repression, and factor potency. Models trained on the dataset correctly predict activity of evolutionarily divergent regulatory regions, providing insights into spatial relationships between repressor and activator binding sites. Importantly, the collective predictions of sets of models were effective at novel enhancer identification and characterization. Our study demonstrates how experimental dataset and modeling can be effectively combined to provide quantitative insights into cis-regulatory information on a genome-wide scale.

DOI:http://dx.doi.org/10.7554/eLife.08445.001

DNA contains regions known as genes, which may be “transcribed” to produce the RNA molecules that act as templates for building proteins and regulate cell activity. Proteins called transcription factors can bind to specific sequences of DNA to influence whether nearby genes are transcribed. For example, so-called enhancer regions of DNA contain several binding sites for transcription factors, and this binding activates gene transcription. Little is known about how the transcription factor binding sites are organized in enhancer regions, which makes it difficult to use DNA sequence information alone to predict the regulation of genes.

A transcription factor called Dorsal controls the activity of a network of genes that plays a crucial role in the development of fruit fly embryos. Dorsal binds to the enhancer region of a gene called rhomboid, which has been well studied and is known to be a fairly typical example of an enhancer region.

To understand the regulatory information encoded in the DNA sequences of enhancers, Sayal, Dresch et al. have now used a technique called perturbation analysis to investigate the interactions that are likely to occur between Dorsal and other transcription factors as they bind to the rhomboid enhancer. This technique involves systematically mutating the enhancer to remove different combinations of transcription factor binding sites and quantitatively investigating the effect this has on gene activity. A large set of mathematical models were then trained using this data and shown to correctly predict the activity of a range of other gene regulatory regions. The collective predictions of the models identified new enhancer regions and revealed details about how different types of transcription factor binding sites are arranged within enhancers.

As we enter an era where the DNA sequences of entire human populations are increasingly accessible, we would like to know the functional significance of changes in gene regulatory regions. Sayal, Dresch et al. show that the regulatory properties of specific control proteins are accessible by employing quantitative experiments and mathematical models. Similar studies will be required to learn how mutations found across the genome may alter gene expression, leading to better diagnosis and treatment of disease.

DOI:http://dx.doi.org/10.7554/eLife.08445.002

Natural variation of the expression pattern of the segmentation gene even-skipped in melanogaster

The evolution of canalized traits is a central question in evolutionary biology. Natural variation in highly conserved traits can provide clues about their evolutionary potential. Here we investigate natural variation in a conserved trait-even-skipped (eve) expression at the cellular blastoderm stage of embryonic development in Drosophila melanogaster. Expression of the pair-rule gene eve was quantitatively measured in three inbred lines derived from a natural population of D. melanogaster. One line showed marked differences in the spacing, amplitude and timing of formation of the characteristic seven-striped pattern over a 50-min period prior to the onset of gastrulation. Stripe 5 amplitude and the width of the interstripe between stripes 4 and 5 were both reduced in this line, while the interstripe distance between stripes 3 and 4 was increased. Engrailed expression in stage 10 embryos revealed a statistically significant increase in the length of parasegment 6 and a decrease in the length of parasegments 8 and 9. These changes are larger than those previously reported between D. melanogaster and D. pseudoobscura, two species that are thought to have diverged from a common ancestor over 25 million years ago. This line harbors a rare 448 bp deletion in the first intron of knirps (kni). This finding suggested that reduced Kni levels caused the deviant eve expression, and indeed we observed lower levels of Kni protein at early cycle 14A in L2 compared to the other two lines. A second of the three lines displayed an approximately 20% greater level of expression for all seven eve stripes. The three lines are each viable and fertile, and none display a segmentation defect as adults, suggesting that early-acting variation in eve expression is ameliorated by developmental buffering mechanisms acting later in development. Canalization of the segmentation pathway may reduce the fitness consequences of genetic variation, thus allowing the persistence of mutations with unexpectedly strong gene expression phenotypes.

A quantitative atlas of Even-skipped and Hunchback expression in Clogmia albipunctata (Diptera: Psychodidae) blastoderm embryos

Background

Comparative studies of developmental processes are one of the main approaches to evolutionary developmental biology (evo-devo). Over recent years, there has been a shift of focus from the comparative study of particular regulatory genes to the level of whole gene networks. Reverse-engineering methods can be used to computationally reconstitute and analyze the function and dynamics of such networks. These methods require quantitative spatio-temporal expression data for model fitting. Obtaining such data in non-model organisms remains a major technical challenge, impeding the wider application of data-driven mathematical modeling to evo-devo.

Results

We have raised antibodies against four segmentation gene products in the moth midge Clogmia albipunctata, a non-drosophilid dipteran species. We have used these antibodies to create a quantitative atlas of protein expression patterns for the gap gene hunchback (hb), and the pair-rule gene even-skipped (eve). Our data reveal differences in the dynamics of Hb boundary positioning and Eve stripe formation between C. albipunctata and Drosophila melanogaster. Despite these differences, the overall relative spatial arrangement of Hb and Eve domains is remarkably conserved between these two distantly related dipteran species.

Conclusions

We provide a proof of principle that it is possible to acquire quantitative gene expression data at high accuracy and spatio-temporal resolution in non-model organisms. Our quantitative data extend earlier qualitative studies of segmentation gene expression in C. albipunctata, and provide a starting point for comparative reverse-engineering studies of the evolutionary and developmental dynamics of the segmentation gene system.

Reverse-engineering post-transcriptional regulation of gap genes in Drosophila melanogaster

Systems biology proceeds through repeated cycles of experiment and modeling. One way to implement this is reverse engineering, where models are fit to data to infer and analyse regulatory mechanisms. This requires rigorous methods to determine whether model parameters can be properly identified. Applying such methods in a complex biological context remains challenging. We use reverse engineering to study post-transcriptional regulation in pattern formation. As a case study, we analyse expression of the gap genes Krüppel, knirps, and giant in Drosophila melanogaster. We use detailed, quantitative datasets of gap gene mRNA and protein expression to solve and fit a model of post-transcriptional regulation, and establish its structural and practical identifiability. Our results demonstrate that post-transcriptional regulation is not required for patterning in this system, but is necessary for proper control of protein levels. Our work demonstrates that the uniqueness and specificity of a fitted model can be rigorously determined in the context of spatio-temporal pattern formation. This greatly increases the potential of reverse engineering for the study of development and other, similarly complex, biological processes.

The analysis of pattern-forming gene networks is largely focussed on transcriptional regulation. However, post-transcriptional events, such as translation and regulation of protein stability also play important roles in the establishment of protein expression patterns and levels. In this study, we use a reverse-engineering approach—fitting mathematical models to quantitative expression data—to analyse post-transcriptional regulation of the Drosophila gap genes Krüppel, knirps and giant, involved in segment determination during early embryogenesis. Rigorous fitting requires us to establish whether our models provide a robust and unique solution. We demonstrate, for the first time, that this can be done in the context of a complex spatio-temporal regulatory system. This is an important methodological advance for reverse-engineering developmental processes. Our results indicate that post-transcriptional regulation is not required for pattern formation, but is necessary for proper regulation of gap protein levels. Specifically, we predict that translation rates must be tuned for rapid early accumulation, and protein stability must be increased for persistence of high protein levels at late stages of gap gene expression.

Quantitative imaging of gene expression in Drosophila embryos

Quantitative measurements derived using sophisticated microscopy techniques are essential for understanding the basic principles that control the behavior of biological systems. Here we describe a data pipeline developed to extract quantitative data on segmentation gene expression from confocal images of gene expression patterns in Drosophila. The pipeline consists of image segmentation, background removal, temporal characterization of an embryo, data registration, and data averaging. This pipeline has been successfully applied to obtain quantitative gene expression data at cellular resolution in space and at 6.5-min resolution in time. It has also enabled the construction of a spatiotemporal atlas of segmentation gene expression. We describe the software used to construct a workflow for extracting quantitative data on segmentation gene expression and the BREReA package, which implements the methods for background removal and registration of segmentation gene expression patterns.

A primer on thermodynamic-based models for deciphering transcriptional regulatory logic

A rigorous analysis of transcriptional regulation at the DNA level is crucial to the understanding of many biological systems. Mathematical modeling has offered researchers a new approach to understanding this central process. In particular, thermodynamic-based modeling represents the most biophysically informed approach aimed at connecting DNA level regulatory sequences to the expression of specific genes. The goal of this review is to give biologists a thorough description of the steps involved in building, analyzing, and implementing a thermodynamic-based model of transcriptional regulation. The data requirements for this modeling approach are described, the derivation for a specific regulatory region is shown, and the challenges and future directions for the quantitative modeling of gene regulation are discussed.

Rearrangements of 2.5 kilobases of noncoding DNA from the Drosophila even-skipped locus define predictive rules of genomic cis-regulatory logic

Rearrangements of about 2.5 kilobases of regulatory DNA located 5' of the transcription start site of the Drosophila even-skipped locus generate large-scale changes in the expression of even-skipped stripes 2, 3, and 7. The most radical effects are generated by juxtaposing the minimal stripe enhancers MSE2 and MSE3 for stripes 2 and 3 with and without small "spacer" segments less than 360 bp in length. We placed these fusion constructs in a targeted transformation site and obtained quantitative expression data for these transformants together with their controlling transcription factors at cellular resolution. These data demonstrated that the rearrangements can alter expression levels in stripe 2 and the 2-3 interstripe by a factor of more than 10. We reasoned that this behavior would place tight constraints on possible rules of genomic cis-regulatory logic. To find these constraints, we confronted our new expression data together with previously obtained data on other constructs with a computational model. The model contained representations of thermodynamic protein-DNA interactions including steric interference and cooperative binding, short-range repression, direct repression, activation, and coactivation. The model was highly constrained by the training data, which it described within the limits of experimental error. The model, so constrained, was able to correctly predict expression patterns driven by enhancers for other Drosophila genes; even-skipped enhancers not included in the training set; stripe 2, 3, and 7 enhancers from various Drosophilid and Sepsid species; and long segments of even-skipped regulatory DNA that contain multiple enhancers. The model further demonstrated that elevated expression driven by a fusion of MSE2 and MSE3 was a consequence of the recruitment of a portion of MSE3 to become a functional component of MSE2, demonstrating that cis-regulatory "elements" are not elementary objects.

Sex-specific pattern formation during early Drosophila development

The deleterious effects of different X-chromosome dosage in males and females are buffered by a process called dosage compensation, which in Drosophila is achieved through a doubling of X-linked transcription in males. The male-specific lethal complex mediates this process, but is known to act only after gastrulation. Recent work has shown that the transcription of X-linked genes is also upregulated in males prior to gastrulation; whether it results in functional dosage compensation is not known. Absent or partial early dosage compensation raises the possibility of sex-biased expression of key developmental genes, such as the segmentation genes controlling anteroposterior patterning. We assess the functional output of early dosage compensation by measuring the expression of even-skipped (eve) with high spatiotemporal resolution in male and female embryos. We show that eve has a sexually dimorphic pattern, suggesting an interaction with either X-chromosome dose or the sex determination system. By manipulating the gene copy number of an X-linked transcription factor, giant (gt), we traced sex-biased eve patterning to gt dose, indicating that early dosage compensation is functionally incomplete. Despite sex-biased eve expression, the gene networks downstream of eve are able to produce sex-independent segmentation, a point that we establish by measuring the proportions of segments in elongated germ-band embryos. Finally, we use a whole-locus eve transgene with modified cis regulation to demonstrate that segment proportions have a sex-dependent sensitivity to subtle changes in Eve expression. The sex independence of downstream segmentation despite this sensitivity to Eve expression implies that additional autosomal gene- or pathway-specific mechanisms are required to ameliorate the effects of partial early dosage compensation.

Lack of tailless leads to an increase in expression variability in Drosophila embryos

Developmental processes are robust, or canalised: dynamic patterns of gene expression across space and time are regulated reliably and precisely in the presence of genetic and environmental perturbations. It remains unclear whether canalisation relies on specific regulatory factors (such as heat-shock proteins), or whether it is based on more general redundancy and distributed robustness at the network level. The latter explanation implies that mutations in many regulatory factors should exhibit loss of canalisation. Here, we present a quantitative characterisation of segmentation gene expression patterns in mutants of the terminal gap gene tailless (tll) in Drosophila melanogaster. Our analysis provides new insights into the dynamic mechanisms underlying gap gene regulation, and reveals significantly increased variability of gene expression in the mutant compared to the wild-type background. We show that both position and timing of posterior segmentation gene expression domains vary strongly from embryo-to-embryo in tll mutants. This variability must be caused by a vulnerability in the regulatory system which is hidden or buffered in the wild-type, but becomes uncovered by the deletion of tll. Our analysis provides evidence that loss of canalisation in mutants could be more widespread than previously thought.

Quantitative dynamics and increased variability of segmentation gene expression in the Drosophila Krüppel and knirps mutants

Here we characterize the response of the Drosophila segmentation system to mutations in two gap genes, Kr and kni, in the form of single or double homozygotes and single heterozygotes. Segmentation gene expression in these genotypes was quantitatively monitored with cellular resolution in space and 6.5 to 13min resolution in time. As is the case with wild type, we found that gene expression domains in the posterior portion of the embryo shift to the anterior over time. In certain cases, such as the gt posterior domain in Kr mutants, the shifts are significantly larger than is seen in wild type embryos. We also investigated the effects of Kr and kni on the variability of gene expression. Mutations often produce variable phenotypes, and it is well known that the cuticular phenotype of Kr mutants is variable. We sought to understand the molecular basis of this effect. We find that throughout cycle 14A the relative levels of eve and ftz expression in stripes 2 and 3 are variable among individual embryos. Moreover, in Kr and kni mutants, unlike wild type, the variability in positioning of the posterior Hb domain and eve stripe 7 is not decreased or filtered with time. The posterior Gt domain in Kr mutants is highly variable at early times, but this variability decreases when this domain shifts in the anterior direction to the position of the neighboring Kni domain. In contrast to these findings, positional variability throughout the embryo does not decrease over time in double Kr;kni mutants. In heterozygotes the early expression patterns of segmentation genes resemble patterns seen in homozygous mutants but by the onset of gastrulation they become similar to the wild type patterns. Finally, we note that gene expression levels are reduced in Kr and kni mutant embryos and have a tendency to decrease over time. This is a surprising result in view of the role that mutual repression is thought to play in the gap gene system.

Spatial gene expression quantification: a tool for analysis of in situ hybridizations in sea anemone Nematostella vectensis

BACKGROUND

Spatial gene expression quantification is required for modeling gene regulation in developing organisms. The fruit fly Drosophila melanogaster is the model system most widely applied for spatial gene expression analysis due to its unique embryonic properties: the shape does not change significantly during its early cleavage cycles and most genes are differentially expressed along a straight axis. This system of development is quite exceptional in the animal kingdom.In the sea anemone Nematostella vectensis the embryo changes its shape during early development; there are cell divisions and cell movement, like in most other metazoans. Nematostella is an attractive case study for spatial gene expression since its transparent body wall makes it accessible to various imaging techniques.

FINDINGS

Our new quantification method produces standardized gene expression profiles from raw or annotated Nematostella in situ hybridizations by measuring the expression intensity along its cell layer. The procedure is based on digital morphologies derived from high-resolution fluorescence pictures. Additionally, complete descriptions of nonsymmetric expression patterns have been constructed by transforming the gene expression images into a three-dimensional representation.

CONCLUSIONS

We created a standard format for gene expression data, which enables quantitative analysis of in situ hybridizations from embryos with various shapes in different developmental stages. The obtained expression profiles are suitable as input for optimization of gene regulatory network models, and for correlation analysis of genes from dissimilar Nematostella morphologies. This approach is potentially applicable to many other metazoan model organisms and may also be suitable for processing data from three-dimensional imaging techniques.

Medium-throughput processing of whole mount in situ hybridisation experiments into gene expression domains

Understanding the function and evolution of developmental regulatory networks requires the characterisation and quantification of spatio-temporal gene expression patterns across a range of systems and species. However, most high-throughput methods to measure the dynamics of gene expression do not preserve the detailed spatial information needed in this context. For this reason, quantification methods based on image bioinformatics have become increasingly important over the past few years. Most available approaches in this field either focus on the detailed and accurate quantification of a small set of gene expression patterns, or attempt high-throughput analysis of spatial expression through binary pattern extraction and large-scale analysis of the resulting datasets. Here we present a robust, “medium-throughput” pipeline to process in situ hybridisation patterns from embryos of different species of flies. It bridges the gap between high-resolution, and high-throughput image processing methods, enabling us to quantify graded expression patterns along the antero-posterior axis of the embryo in an efficient and straightforward manner. Our method is based on a robust enzymatic (colorimetric) in situ hybridisation protocol and rapid data acquisition through wide-field microscopy. Data processing consists of image segmentation, profile extraction, and determination of expression domain boundary positions using a spline approximation. It results in sets of measured boundaries sorted by gene and developmental time point, which are analysed in terms of expression variability or spatio-temporal dynamics. Our method yields integrated time series of spatial gene expression, which can be used to reverse-engineer developmental gene regulatory networks across species. It is easily adaptable to other processes and species, enabling the in silico reconstitution of gene regulatory networks in a wide range of developmental contexts.

Modeling of gap gene expression in Drosophila Kruppel mutants

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.

Systems biology is aimed to develop an understanding of biological function or process as a system of interacting components. Here we apply the systems-level approach to understand how the blueprints for segments in the fruit fly Drosophila embryo arise. We obtain gene expression data and use the gene circuits method which allow us to reconstruct the segment determination process in the computer. To understand the system we need not only to describe it in detail, but also to comprehend what happens when certain stimuli or disruptions occur. Previous attempts to model segmentation gene expression patterns in a mutant for a trunk gap gene were unsuccessful. Here we describe the extension of the model that allows us to solve this problem in the context of Kruppel (Kr) gene. We show that remarkable alteration of gap gene expression patterns in Kr mutants can be explained by dynamic decrease of the activating effect of Cad on a target gene and exclusion of Kr 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 role of Bicoid cooperative binding in the patterning of sharp borders in Drosophila melanogaster

In Drosophila embryonic development, the Bicoid (Bcd) protein establishes positional information of downstream developmental genes like hunchback (hb), which has a strong anterior expression and a sharp on-off boundary in the mid-embryo. The role of Bcd cooperative binding in the positioning of the Hb pattern has been previously demonstrated. However, there are discrepancies in the reported results about the role of this mechanism in the sharp Hb border. Here, we determined the Hill coefficient (nH) required for Bcd to generate the sharp border of Hb in wild-type (WT) embryos. We found that an n(H) of approximately 6.3 (s.d. 1.4) and 10.8 (s.d. 4.0) is required to account for Hb sharpness at early and late cycle 14A, respectively. Additional mechanisms are possibly required because the high nH is likely unachievable for Bcd binding to the hb promoter. To test this idea, we determined the nH required to pattern the Hb profile of 15 embryos expressing an hb14F allele that is defective in self-activation and found nH to be 3.0 (s.d. 1.0). This result indicates that in WT embryos, the hb self-activation is important for Hb sharpness. Corroborating our results, we also found a progressive increase in the required value of n(H) spanning from 4.0 to 9.2 by determining this coefficient from averaged profiles of eight temporal classes at cycle 14A (T1 to T8). Our results indicate that there is a transition in the mechanisms responsible for the sharp Hb border during cycle 14A: in early stages of this cycle, Bcd cooperative binding is primarily responsible for Hb sharpness; in late cycle 14A, hb self-activation becomes the dominant mechanism.

Efficient reverse-engineering of a developmental gene regulatory network

Understanding the complex regulatory networks underlying development and evolution of multi-cellular organisms is a major problem in biology. Computational models can be used as tools to extract the regulatory structure and dynamics of such networks from gene expression data. This approach is called reverse engineering. It has been successfully applied to many gene networks in various biological systems. However, to reconstitute the structure and non-linear dynamics of a developmental gene network in its spatial context remains a considerable challenge. Here, we address this challenge using a case study: the gap gene network involved in segment determination during early development of Drosophila melanogaster. A major problem for reverse-engineering pattern-forming networks is the significant amount of time and effort required to acquire and quantify spatial gene expression data. We have developed a simplified data processing pipeline that considerably increases the throughput of the method, but results in data of reduced accuracy compared to those previously used for gap gene network inference. We demonstrate that we can infer the correct network structure using our reduced data set, and investigate minimal data requirements for successful reverse engineering. Our results show that timing and position of expression domain boundaries are the crucial features for determining regulatory network structure from data, while it is less important to precisely measure expression levels. Based on this, we define minimal data requirements for gap gene network inference. Our results demonstrate the feasibility of reverse-engineering with much reduced experimental effort. This enables more widespread use of the method in different developmental contexts and organisms. Such systematic application of data-driven models to real-world networks has enormous potential. Only the quantitative investigation of a large number of developmental gene regulatory networks will allow us to discover whether there are rules or regularities governing development and evolution of complex multi-cellular organisms.

Measuring gene expression noise in early Drosophila embryos: nucleus-to-nucleus variability

In recent years the analysis of noise in gene expression has widely attracted the attention of experimentalists and theoreticians. Experimentally, the approaches based on in vivo fluorescent reporters in single cells appear to be straightforward and effective tools for bacteria and yeast. However, transferring these approaches to multicellular organisms presents many methodological problems. Here we describe our approach to measure between-nucleus variability (noise) in the primary morphogenetic gradient of Bicoid (Bcd) in the precellular blastoderm stage of fruit fly (Drosophila) embryos. The approach is based on the comparison of results for fixed immunostained embryos with observations of live embryos carrying fluorescent Bcd (Bcd-GFP). We measure the noise using two-dimensional Singular Spectrum Analysis (2D SSA). We have found that the nucleus-to-nucleus noise in Bcd intensity, both for live (Bcd-GFP) and for fixed immunstained embryos, tends to be signal-independent. In addition, the character of the noise is sensitive to the nuclear masking technique used to extract quantitative intensities. Further, the method of decomposing the raw quantitative expression data into a signal (expression surface) and residual noise affects the character of the residual noise. We find that careful masking of confocal images and use of appropriate computational tools to decompose raw expression data into trend and noise makes it possible to extract and study the biological noise of gene expression.

Mathematics and biology: a Kantian view on the history of pattern formation theory

Driesch's statement, made around 1900, that the physics and chemistry of his day were unable to explain self-regulation during embryogenesis was correct and could be extended until the year 1972. The emergence of theories of self-organisation required progress in several areas including chemistry, physics, computing and cybernetics. Two parallel lines of development can be distinguished which both culminated in the early 1970s. Firstly, physicochemical theories of self-organisation arose from theoretical (Lotka 1910-1920) and experimental work (Bray 1920; Belousov 1951) on chemical oscillations. However, this research area gained broader acceptance only after thermodynamics was extended to systems far from equilibrium (1922-1967) and the mechanism of the prime example for a chemical oscillator, the Belousov-Zhabotinski reaction, was deciphered in the early 1970s. Secondly, biological theories of self-organisation were rooted in the intellectual environment of artificial intelligence and cybernetics. Turing wrote his The chemical basis of morphogenesis (1952) after working on the construction of one of the first electronic computers. Likewise, Gierer and Meinhardt's theory of local activation and lateral inhibition (1972) was influenced by ideas from cybernetics. The Gierer-Meinhardt theory provided an explanation for the first time of both spontaneous formation of spatial order and of self-regulation that proved to be extremely successful in elucidating a wide range of patterning processes. With the advent of developmental genetics in the 1980s, detailed molecular and functional data became available for complex developmental processes, allowing a new generation of data-driven theoretical approaches. Three examples of such approaches will be discussed. The successes and limitations of mathematical pattern formation theory throughout its history suggest a picture of the organism, which has structural similarity to views of the organic world held by the philosopher Immanuel Kant at the end of the eighteenth century.

Consequences of eukaryotic enhancer architecture for gene expression dynamics, development, and fitness

The regulatory logic of time- and tissue-specific gene expression has mostly been dissected in the context of the smallest DNA fragments that, when isolated, recapitulate native expression in reporter assays. It is not known if the genomic sequences surrounding such fragments, often evolutionarily conserved, have any biological function or not. Using an enhancer of the even-skipped gene of Drosophila as a model, we investigate the functional significance of the genomic sequences surrounding empirically identified enhancers. A 480 bp long "minimal stripe element" is able to drive even-skipped expression in the second of seven stripes but is embedded in a larger region of 800 bp containing evolutionarily conserved binding sites for required transcription factors. To assess the overall fitness contribution made by these binding sites in the native genomic context, we employed a gene-replacement strategy in which whole-locus transgenes, capable of rescuing even-skipped(-) lethality to adulthood, were substituted for the native gene. The molecular phenotypes were characterized by tagging Even-skipped with a fluorescent protein and monitoring gene expression dynamics in living embryos. We used recombineering to excise the sequences surrounding the minimal enhancer and site-specific transgenesis to create co-isogenic strains differing only in their stripe 2 sequences. Remarkably, the flanking sequences were dispensable for viability, proving the sufficiency of the minimal element for biological function under normal conditions. These sequences are required for robustness to genetic and environmental perturbation instead. The mutant enhancers had measurable sex- and dose-dependent effects on viability. At the molecular level, the mutants showed a destabilization of stripe placement and improper activation of downstream genes. Finally, we demonstrate through live measurements that the peripheral sequences are required for temperature compensation. These results imply that seemingly redundant regulatory sequences beyond the minimal enhancer are necessary for robust gene expression and that "robustness" itself must be an evolved characteristic of the wild-type enhancer.

A computational statistics approach for estimating the spatial range of morphogen gradients

A crucial issue in studies of morphogen gradients relates to their range: the distance over which they can act as direct regulators of cell signaling, gene expression and cell differentiation. To address this, we present a straightforward statistical framework that can be used in multiple developmental systems. We illustrate the developed approach by providing a point estimate and confidence interval for the spatial range of the graded distribution of nuclear Dorsal, a transcription factor that controls the dorsoventral pattern of the Drosophila embryo.

A regression system for estimation of errors introduced by confocal imaging into gene expression data in situ

Background

Accuracy of the data extracted from two-dimensional confocal images is limited due to experimental errors that arise in course of confocal scanning. The common way to reduce the noise in images is sequential scanning of the same specimen several times with the subsequent averaging of multiple frames. Attempts to increase the dynamical range of an image by setting too high values of microscope PMT parameters may cause clipping of single frames and introduce errors into the data extracted from the averaged images. For the estimation and correction of this kind of errors a method based on censoring technique (Myasnikova et al., 2009) is used. However, the method requires the availability of all the confocal scans along with the averaged image, which is normally not provided by the standard scanning procedure.

Results

To predict error size in the data extracted from the averaged image we developed a regression system. The system is trained on the learning sample composed of images obtained from three different microscopes at different combinations of PMT parameters, and for each image all the scans are saved. The system demonstrates high prediction accuracy and was applied for correction of errors in the data on segmentation gene expression in Drosophila blastoderm stored in the FlyEx database (http://urchin.spbcas.ru/flyex/, http://flyex.uchicago.edu/flyex/). The prediction method is realized as a software tool CorrectPattern freely available at http://urchin.spbcas.ru/asp/2011/emm/.

Conclusions

We created a regression system and software to predict the magnitude of errors in the data obtained from a confocal image based on information about microscope parameters used for the image acquisition. An important advantage of the developed prediction system is the possibility to accurately correct the errors in data obtained from strongly clipped images, thereby allowing to obtain images of the higher dynamical range and thus to extract more detailed quantitative information from them.

Gene expression noise in spatial patterning: hunchback promoter structure affects noise amplitude and distribution in Drosophila segmentation

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.

Non-additive interactions involving two distinct elements mediate sloppy-paired regulation by pair-rule transcription factors

The relatively simple combinatorial rules responsible for establishing the initial metameric expression of sloppy-paired-1 (slp1) in the Drosophila blastoderm embryo make this system an attractive model for investigating the mechanism of regulation by pair-rule transcription factors. This investigation of slp1 cis-regulatory architecture identifies two distinct elements, a proximal early stripe element (PESE) and a distal early stripe element (DESE) located from -3.1kb to -2.5kb and from -8.1kb to -7.1kb upstream of the slp1 promoter, respectively, that mediate this early regulation. The proximal element expresses only even-numbered stripes and mediates repression by Even-skipped (Eve) as well as by the combination of Runt and Fushi-tarazu (Ftz). A 272 basepair sub-element of PESE retains an Eve-dependent repression, but is expressed throughout the even-numbered parasegments due to the loss of repression by Runt and Ftz. In contrast, the distal element expresses both odd and even-numbered stripes and also drives inappropriate expression in the anterior half of the odd-numbered parasegments due to an inability to respond to repression by Eve. Importantly, a composite reporter gene containing both early stripe elements recapitulates pair-rule gene-dependent regulation in a manner beyond what is expected from combining their individual patterns. These results indicate that interactions involving distinct cis-elements contribute to the proper integration of pair-rule regulatory information. A model fully accounting for these results proposes that metameric slp1 expression is achieved through the Runt-dependent regulation of interactions between these two pair-rule response elements and the slp1 promoter.

Extraction and comparison of gene expression patterns from 2D RNA in situ hybridization images

MOTIVATION

Recent advancements in high-throughput imaging have created new large datasets with tens of thousands of gene expression images. Methods for capturing these spatial and/or temporal expression patterns include in situ hybridization or fluorescent reporter constructs or tags, and results are still frequently assessed by subjective qualitative comparisons. In order to deal with available large datasets, fully automated analysis methods must be developed to properly normalize and model spatial expression patterns.

RESULTS

We have developed image segmentation and registration methods to identify and extract spatial gene expression patterns from RNA in situ hybridization experiments of Drosophila embryos. These methods allow us to normalize and extract expression information for 78,621 images from 3724 genes across six time stages. The similarity between gene expression patterns is computed using four scoring metrics: mean squared error, Haar wavelet distance, mutual information and spatial mutual information (SMI). We additionally propose a strategy to calculate the significance of the similarity between two expression images, by generating surrogate datasets with similar spatial expression patterns using a Monte Carlo swap sampler. On data from an early development time stage, we show that SMI provides the most biologically relevant metric of comparison, and that our significance testing generalizes metrics to achieve similar performance. We exemplify the application of spatial metrics on the well-known Drosophila segmentation network.

AVAILABILITY

A Java webstart application to register and compare patterns, as well as all source code, are available from: http://tools.genome.duke.edu/generegulation/image_analysis/insitu

CONTACT

uwe.ohler@duke.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Statistical lower bounds on protein copy number from fluorescence expression images

MOTIVATION

Fluorescence imaging has become a commonplace for quantitatively measuring mRNA or protein expression in cells and tissues. However, such expression data are usually relative-absolute concentrations or molecular copy numbers are typically not known. While this is satisfactory for many applications, for certain kinds of quantitative network modeling and analysis of expression noise, absolute measures of expression are necessary.

RESULTS

We propose two methods for estimating molecular copy numbers from single uncalibrated expression images of tissues. These methods rely on expression variability between cells, due either to steady-state fluctuations or unequal distribution of molecules during cell division, to make their estimates. We apply these methods to 152 protein fluorescence expression images of Drosophila melanogaster embryos during early development, generating copy number estimates for 14 genes in the segmentation network. We also analyze the effects of noise on our estimators and compare with empirical findings. Finally, we confirm an observation of Bar-Even et al., made in the much different setting of Saccharomyces cerevisiae, that steady-state expression variance tends to scale with mean expression.

AVAILABILITY

The data are all drawn from FlyEx (explained within), and is available at http://flyex.ams.sunysb.edu/FlyEx/.

Gene circuit analysis of the terminal gap gene huckebein

The early embryo of Drosophila melanogaster provides a powerful model system to study the role of genes in pattern formation. The gap gene network constitutes the first zygotic regulatory tier in the hierarchy of the segmentation genes involved in specifying the position of body segments. Here, we use an integrative, systems-level approach to investigate the regulatory effect of the terminal gap gene huckebein (hkb) on gap gene expression. We present quantitative expression data for the Hkb protein, which enable us to include hkb in gap gene circuit models. Gap gene circuits are mathematical models of gene networks used as computational tools to extract regulatory information from spatial expression data. This is achieved by fitting the model to gap gene expression patterns, in order to obtain estimates for regulatory parameters which predict a specific network topology. We show how considering variability in the data combined with analysis of parameter determinability significantly improves the biological relevance and consistency of the approach. Our models are in agreement with earlier results, which they extend in two important respects: First, we show that Hkb is involved in the regulation of the posterior hunchback (hb) domain, but does not have any other essential function. Specifically, Hkb is required for the anterior shift in the posterior border of this domain, which is now reproduced correctly in our models. Second, gap gene circuits presented here are able to reproduce mutants of terminal gap genes, while previously published models were unable to reproduce any null mutants correctly. As a consequence, our models now capture the expression dynamics of all posterior gap genes and some variational properties of the system correctly. This is an important step towards a better, quantitative understanding of the developmental and evolutionary dynamics of the gap gene network.

GCPReg package for registration of the segmentation gene expression data in Drosophila

In modern functional genomics registration techniques areused to construct reference gene expression patterns and createa spatiotemporal atlas of the expression of all the genes in anetwork. In this paper we present a software package calledGCPReg, which can be used to register the expression patterns ofsegmentation genes in the early Drosophila embryo. The key task,which this package performs, is the extraction of spatially localizedcharacteristic features of expression patterns. To facilitatethis task, we have developed an easy-to-use interactive graphicalinterface. We describe GCPReg usage and demonstrate how thispackage can be applied to register gene expression patterns inwild type and mutants. GCPReg has been designed to operate ona UNIX platform and is freely available via the Internet at http://urchin.spbcas.ru/downloads/GCPReg/GCPReg.htm.

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.

Animals have an astonishing ability to develop reliably in spite of variable conditions during embryogenesis. More than 60 years ago, it was proposed that this property of development, called canalization, results from genetic interactions that adjust biochemical reactions so as to bring about reliable outcomes. Since then, a great deal of progress has been made in understanding the buffering of genotypic and environmental variation, and individual mutations that reveal variation have been identified. However, the mechanisms by which genetic interactions produce canalization are not yet well understood, because this requires molecular data on multiple developmental determinants and models that correctly predict complex interactions. We make use of gene expression data at both high spatial and temporal resolution for the gap genes involved in the segmentation of Drosophila. We also apply a mathematical model to show that cross regulation among the gap genes is responsible for canalization in this system. Furthermore, the model predicted specific interactions that cause canalization, and the prediction was validated experimentally. Our results show that groups of genes can act on one another to reduce variation and highlights the importance of genetic networks in generating robust development.

DuringDrosophila development, the expression patterns of gap genes are much less variable than the Bicoid morphogen gradient. Modeling and experiments show that this specific instance of canalization or developmental robustness occurs by gap gene cross regulation.

Estimation of errors introduced by confocal imaging into the data on segmentation gene expression in Drosophila

MOTIVATION

Currently the confocal scanning microscopy of fluorescently tagged molecules is extensively employed to acquire quantitative data on gene expression at cellular resolution. Following this approach, we generated a large dataset on the expression of segmentation genes in the Drosophila blastoderm, that is widely used in systems biology studies. As data accuracy is of critical importance for the success of studies in this field, we took a shot to evaluate possible errors introduced in the data by acquisition and processing methods. This article deals with errors introduced by confocal microscope.

RESULTS

In confocal imaging, the inevitable photon noise is commonly reduced by the averaging of multiple frames. The averaging may introduce errors into the data, if single frames are clipped by microscope hardware. A method based on censoring technique is used to estimate and correct this type of errors. Additional source of errors is the quantification of blurred images. To estimate and correct these errors, the Richardson-Lucy deconvolution method was modified to provide the higher accuracy of data read off from blurred images of the Drosophila blastoderm. We have found that the sizes of errors introduced by confocal imaging make up approximately 5-7% of the mean intensity values and do not disguise the dynamic behavior and characteristic features of gene expression patterns. We also defined a range of microscope parameters for the acquisition of sufficiently accurate data.

AVAILABILITY

http://urchin.spbcas.ru/downloads/step/step.htm

Image processing and analysis for quantifying gene expression from early Drosophila embryos

Correlation of quantities of transcriptional activators and repressors with the mRNA output of target genes is a central issue for modeling gene regulation. In multicellular organisms, both spatial and temporal differences in gene expression must be taken into account; this can be achieved by use of in situ hybridization followed by confocal laser scanning microscopy (CLSM). Here we present a method to correlate the protein levels of the short-range repressor Giant with lacZ mRNA produced by reporter genes using images of Drosophila blastoderm embryos taken by CLSM. The image stacks from CLSM are processed using a semiautomatic algorithm to produce correlations between the repressor levels and lacZ mRNA reporter genes. We show that signals derived from CLSM are proportional to actual mRNA levels. Our analysis reveals that a suggested parabolic form of the background fluorescence in confocal images of early Drosophila embryos is evident most prominently in flattened specimens, with intact embryos exhibiting a more linear background. The data extraction described in this paper is primarily conceived for analysis of synthetic reporter genes that are designed to decipher cis-regulatory grammar, but the techniques are generalizable for quantitative analysis of other engineered or endogenous genes in embryos.

FlyEx, the quantitative atlas on segmentation gene expression at cellular resolution

The datasets on gene expression are the valuable source of information about the functional state of an organism. Recently, we have acquired the large dataset on expression of segmentation genes in the Drosophila blastoderm. To provide efficient access to the data, we have developed the FlyEx database (http://urchin.spbcas.ru/flyex). FlyEx contains 4716 images of 14 segmentation gene expression patterns obtained from 1579 embryos and 9,500,000 quantitative data records. Reference data are available for all segmentation genes in cycles 11-13 and all temporal classes of cycle 14A. FlyEx supports operations on images of gene expression patterns. The database can be used to examine the quality of data, analyze the dynamics of formation of segmentation gene expression domains, as well as to estimate the variability of gene expression patterns. Currently, a user is able to monitor and analyze the dynamics of formation of segmentation gene expression domains over the whole period of segment determination, that amounts to 1.5 h of development. FlyEx supports the data downloads and construction of personal reference datasets, that makes it possible to more effectively use and analyze data.

Spatial bistability generates hunchback expression sharpness in the Drosophila embryo

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.

Methods for Acquisition of Quantitative Data from Confocal Images of Gene Expression in situ

In this review, we summarize original methods for the extraction of quantitative information from confocal images of gene-expression patterns. These methods include image segmentation, the extraction of quantitative numerical data on gene expression, and the removal of background signal and spatial registration. Finally, it is possible to construct a spatiotemporal atlas of gene expression from individual images recorded at each developmental stage. Initially all methods were developed to extract quantitative numerical information from confocal images of segmentation gene expression in Drosophila melanogaster. The application of these methods to Drosophila images makes it possible to reveal new mechanisms in the formation of segmentation gene expression domains, as well as to construct a quantitative atlas of segmentation gene expression. Most image processing procedures can be easily adapted to process a wide range of biological images.

Pipeline for acquisition of quantitative data on segmentation gene expression from confocal images

We describe a data pipeline developed to extract the quantitative data on segmentation gene expression from confocal images of gene expression patterns in Drosophila. The pipeline consists of five steps: image segmentation, background removal, temporal characterization of an embryo, data registration and data averaging. This pipeline was successfully applied to obtain quantitative gene expression data at cellular resolution in space and at the 6.5-minute resolution in time, as well as to construct a spatiotemporal atlas of segmentation gene expression. Each data pipeline step can be easily adapted to process a wide range of images of gene expression patterns.

Characterization of the Drosophila segment determination morphome

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.

Master equation simulation analysis of immunostained Bicoid morphogen gradient

Background

The concentration gradient of Bicoid protein which determines the developmental pathways in early Drosophila embryo is the best characterized morphogen gradient at the molecular level. Because different developmental fates can be elicited by different concentrations of Bicoid, it is important to probe the limits of this specification by analyzing intrinsic fluctuations of the Bicoid gradient arising from small molecular number. Stochastic simulations can be applied to further the understanding of the dynamics of Bicoid morphogen gradient formation at the molecular number level, and determine the source of the nucleus-to-nucleus expression variation (noise) observed in the Bicoid gradient.

Results

We compared quantitative observations of Bicoid levels in immunostained Drosophila embryos with a spatially extended Master Equation model which represents diffusion, decay, and anterior synthesis. We show that the intrinsic noise of an autonomous reaction-diffusion gradient is Poisson distributed. We demonstrate how experimental noise can be identified in the logarithm domain from single embryo analysis, and then separated from intrinsic noise in the normalized variance domain of an ensemble statistical analysis. We show how measurement sensitivity affects our observations, and how small amounts of rescaling noise can perturb the noise strength (Fano factor) observed. We demonstrate that the biological noise level in data can serve as a physical constraint for restricting the model's parameter space, and for predicting the Bicoid molecular number and variation range. An estimate based on a low variance ensemble of embryos suggests that the steady-state Bicoid molecular number in a nucleus should be larger than 300 in the middle of the embryo, and hence the gradient should extend to the posterior end of the embryo, beyond the previously assumed background limit. We exhibit the predicted molecular number gradient together with measurement effects, and make a comparison between conditions of higher and lower variance respectively.

Conclusion

Quantitative comparison of Master Equation simulations with immunostained data enabled us to determine narrow ranges for key biophysical parameters, which for this system can be independently validated. Intrinsic noise is clearly detectable as well, although the staining process introduces certain limits in resolution.

Three-dimensional morphology and gene expression in the Drosophila blastoderm at cellular resolution I: data acquisition pipeline

A suite of methods that provide the first quantitative three-dimensional description of gene expression and morphology with cellular resolution in whole Drosophila embryos is described.

Background

To model and thoroughly understand animal transcription networks, it is essential to derive accurate spatial and temporal descriptions of developing gene expression patterns with cellular resolution.

Results

Here we describe a suite of methods that provide the first quantitative three-dimensional description of gene expression and morphology at cellular resolution in whole embryos. A database containing information derived from 1,282 embryos is released that describes the mRNA expression of 22 genes at multiple time points in the Drosophila blastoderm. We demonstrate that our methods are sufficiently accurate to detect previously undescribed features of morphology and gene expression. The cellular blastoderm is shown to have an intricate morphology of nuclear density patterns and apical/basal displacements that correlate with later well-known morphological features. Pair rule gene expression stripes, generally considered to specify patterning only along the anterior/posterior body axis, are shown to have complex changes in stripe location, stripe curvature, and expression level along the dorsal/ventral axis. Pair rule genes are also found to not always maintain the same register to each other.

Conclusion

The application of these quantitative methods to other developmental systems will likely reveal many other previously unknown features and provide a more rigorous understanding of developmental regulatory networks.

On the dynamic nature of positional information

Morphogenetic fields are among the most fundamental concepts of embryology. However, they are also among the most ill-defined, since they consist of dynamic regulatory processes whose exact nature remains elusive. In order to achieve a more rigorous definition of a developmental field, Lewis Wolpert introduced the concept of positional information illustrated by his French Flag model. Here we argue that Wolpert's positional information - a static coordinate system defining a field - lacks essential properties of the original field concept. We show how data-driven mathematical modeling approaches now enable us to study regulatory processes in a way that is qualitatively different from our previous level of understanding. As an example, we review our recent analysis of segmentation gene expression in the blastoderm embryo of the fruit fly Drosophila melanogaster. Based on this analysis, we propose a revised French Flag, which incorporates the dynamic, feedback-driven nature of pattern formation in the Drosophila blastoderm.

Known maternal gradients are not sufficient for the establishment of gap domains in Drosophila melanogaster

Gap genes are among the first transcriptional targets of maternal morphogen gradients in the early Drosophila embryo. However, it remains unclear whether these gradients are indeed sufficient to establish the boundaries of localized gap gene expression patterns. In this study, we address this question using gap gene circuits, which are data-driven mathematical tools for extracting regulatory information from quantitative wild-type gene expression data. We present new, quantitative data on the mRNA expression patterns for the gap genes Krüppel (Kr), knirps (kni) and giant (gt) during the early blastoderm stage of Drosophila development. This data set shows significant differences in timing of gap gene expression compared to previous observations, and reveals that early gap gene expression is highly variable in both space and time. Gene circuit models fit to this data set were used for a detailed regulatory analysis of early gap gene expression. Our analysis shows that the proper balance of maternal repression and activation is essential for the correct positioning of gap domains, and that such balance can only be achieved for early expression of kni. In contrast, our results suggest that early expression of gt requires local neutralization of repressive input in the anterior region of the embryo, and that known maternal gradients are completely insufficient to position the boundaries of the early central Kr domain, or in fact any Kr-like domain in the central region of the blastoderm embryo. Based on this, we propose that unknown additional regulators must be involved in early gap gene regulation.

Analysis of pattern precision shows that Drosophila segmentation develops substantial independence from gradients of maternal gene products

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.

Quantitative and predictive model of transcriptional control of the Drosophila melanogaster even skipped gene

Here we present a quantitative and predictive model of the transcriptional readout of the proximal 1.7 kb of the control region of the Drosophila melanogaster gene even skipped (eve). The model is based on the positions and sequence of individual binding sites on the DNA and quantitative, time-resolved expression data at cellular resolution. These data demonstrated new expression features, first reported here. The model correctly predicts the expression patterns of mutations in trans, as well as point mutations, insertions and deletions in cis. It also shows that the nonclassical expression of stripe 7 driven by this fragment is activated by the protein Caudal (Cad), and repressed by the proteins Tailless (Tll) and Giant (Gt).