Precision of hunchback expression in the Drosophila embryo

Gene activity fully predicts transcriptional bursting dynamics

Transcription commonly occurs in bursts, with alternating productive (ON) and quiescent (OFF) periods, governing mRNA production rates. Yet, how transcription is regulated through bursting dynamics remains unresolved. Here, we conduct real-time measurements of endogenous transcriptional bursting with single-mRNA sensitivity. Leveraging the diverse transcriptional activities in early fly embryos, we uncover stringent relationships between bursting parameters. Specifically, we find that the durations of ON and OFF periods are linked. Regardless of the developmental stage or body-axis position, gene activity levels predict individual alleles' average ON and OFF periods. Lowly transcribing alleles predominantly modulate OFF periods (burst frequency), while highly transcribing alleles primarily tune ON periods (burst size). These relationships persist even under perturbations of cis-regulatory elements or trans-factors and account for bursting dynamics measured in other species. Our results suggest a novel mechanistic constraint governing bursting dynamics rather than a modular control of distinct parameters by distinct regulatory processes.

Multifaceted effects on even-skipped transcriptional dynamics upon Krüppel dosage changes

Although fluctuations in transcription factor (TF) dosage are often well tolerated, TF dosage modulation can change the target gene expression dynamics and result in significant non-lethal developmental phenotypes. Using MS2/MCP-mediated quantitative live imaging in early Drosophila embryos, we analyzed how changing levels of the gap gene Krüppel (Kr) affects transcriptional dynamics of the pair-rule gene even-skipped (eve). Halving the Kr dosage leads to a transient posterior expansion of the eve stripe 2 and an anterior shift of stripe 5. Surprisingly, the most significant changes are observed in eve stripes 3 and 4, the enhancers of which do not contain Kr-binding sites. In Kr heterozygous embryos, both stripes 3 and 4 display narrower widths, anteriorly shifted boundaries and reduced mRNA production levels. We show that Kr dosage indirectly affects stripe 3 and 4 dynamics by modulating other gap gene dynamics. We quantitatively correlate moderate body segment phenotypes of Kr heterozygotes with spatiotemporal changes in eve expression. Our results indicate that nonlinear relationships between TF dosage and phenotypes underlie direct TF-DNA and indirect TF-TF interactions.

LucFlow: A method to measure Luciferase reporter expression in single cells

Reporter assays, in which the expression of an inert protein is driven by gene regulatory elements such as promoters and enhancers, are a workhorse for investigating gene regulation. Techniques for measuring reporter gene expression vary from single-cell or single-molecule approaches having low throughput to bulk Luciferase assays that have high throughput. We developed a Luciferase Reporter Assay using Flow-Cytometry (LucFlow), which measures reporter expression in single cells immunostained for Luciferase. We optimized and tested LucFlow with a murine cell line that can be differentiated into neutrophils, into which promoter-reporter and enhancer-promoter-reporter constructs have been integrated in a site-specific manner. The single-cell measurements are comparable to bulk ones but we found that dead cells have no detectable Luciferase protein, so that bulk assays underestimate reporter expression. LucFlow is able to achieve a higher accuracy than bulk methods by excluding dead cells during flow cytometry. Prior to fixation and staining, the samples are spiked with stained cells that can be discriminated during flow cytometry and control for tube-to-tube variation in experimental conditions. Computing fold change relative to control cells allows LucFlow to achieve a high level of precision. LucFlow, therefore, enables the accurate and precise measurement of reporter expression in a high throughput manner.

Setting the stage for development: the maternal-to-zygotic transition in Drosophila

The zygote has a daunting task ahead of itself; it must develop from a single cell (fertilized egg) into a fully functioning adult with a multitude of different cell types. In the beginning, the zygote has help from its mother, in the form of gene products deposited into the egg, but eventually, it must rely on its own resources to proceed through development. The transfer of developmental control from the mother to the embryo is called the maternal-to-zygotic transition (MZT). All animals undergo this transition, which is defined by two main processes-the degradation of maternal RNAs and the synthesis of new RNAs from the zygote's own genome. Here, we review the regulation of the MZT in Drosophila, but given the broad conservation of this essential process, much of the regulation is shared among metazoans.

Shaping the scaling characteristics of gap gene expression patterns in Drosophila

How patterns are formed to scale with tissue size remains an unresolved problem. Here we investigate embryonic patterns of gap gene expression along the anterior-posterior (AP) axis in Drosophila. We use embryos that greatly differ in length and, importantly, possess distinct length-scaling characteristics of the Bicoid (Bcd) gradient. We systematically analyze the dynamic movements of gap gene expression boundaries in relation to both embryo length and Bcd input as a function of time. We document the process through which such dynamic movements drive both an emergence of a global scaling landscape and evolution of boundary-specific scaling characteristics. We show that, despite initial differences in pattern scaling characteristics that mimic those of Bcd in the anterior, such characteristics of final patterns converge. Our study thus partitions the contributions of Bcd input and regulatory dynamics inherent to the AP patterning network in shaping embryonic pattern's scaling characteristics.

Differential regulation of alternative promoters emerges from unified kinetics of enhancer-promoter interaction

Many eukaryotic genes contain alternative promoters with distinct expression patterns. How these promoters are differentially regulated remains elusive. Here, we apply single-molecule imaging to quantify the transcriptional regulation of two alternative promoters (P1 and P2) of the Bicoid (Bcd) target gene hunchback in syncytial blastoderm Drosophila embryos. Contrary to the previous notion that Bcd only activates P2, we find that Bcd activates both promoters via the same two enhancers. P1 activation is less frequent and requires binding of more Bcd molecules than P2 activation. Using a theoretical model to relate promoter activity to enhancer states, we show that the two promoters follow common transcription kinetics driven by sequential Bcd binding at the two enhancers. Bcd binding at either enhancer primarily activates P2, while P1 activation relies more on Bcd binding at both enhancers. These results provide a quantitative framework for understanding the kinetic mechanisms of complex eukaryotic gene regulation.

Alternative promoters differ in their expression patterns, whose mechanisms are not well understood. Here the authors show that alternative promoters of a Drosophila embryonic gene hunchback are regulated by different action modes of two enhancers.

Application of Domain- and Genotype-Specific Models to Infer Post-Transcriptional Regulation of Segmentation Gene Expression in Drosophila

Unlike transcriptional regulation, the post-transcriptional mechanisms underlying zygotic segmentation gene expression in early Drosophila embryo have been insufficiently investigated. Condition-specific post-transcriptional regulation plays an important role in the development of many organisms. Our recent study revealed the domain- and genotype-specific differences between mRNA and the protein expression of Drosophila hb, gt, and eve genes in cleavage cycle 14A. Here, we use this dataset and the dynamic mathematical model to recapitulate protein expression from the corresponding mRNA patterns. The condition-specific nonuniformity in parameter values is further interpreted in terms of possible post-transcriptional modifications. For hb expression in wild-type embryos, our results predict the position-specific differences in protein production. The protein synthesis rate parameter is significantly higher in hb anterior domain compared to the posterior domain. The parameter sets describing Gt protein dynamics in wild-type embryos and Kr mutants are genotype-specific. The spatial discrepancy between gt mRNA and protein posterior expression in Kr mutants is well reproduced by the whole axis model, thus rejecting the involvement of post-transcriptional mechanisms. Our models fail to describe the full dynamics of eve expression, presumably due to its complex shape and the variable time delays between mRNA and protein patterns, which likely require a more complex model. Overall, our modeling approach enables the prediction of regulatory scenarios underlying the condition-specific differences between mRNA and protein expression in early embryo.

Dynamics of hunchback translation in real-time and at single-mRNA resolution in the Drosophila embryo

The Hunchback (Hb) transcription factor is crucial for anterior-posterior patterning of the Drosophila embryo. The maternal hb mRNA acts as a paradigm for translational regulation due to its repression in the posterior of the embryo. However, little is known about the translatability of zygotically transcribed hb mRNAs. Here, we adapt the SunTag system, developed for imaging translation at single-mRNA resolution in tissue culture cells, to the Drosophila embryo to study the translation dynamics of zygotic hb mRNAs. Using single-molecule imaging in fixed and live embryos, we provide evidence for translational repression of zygotic SunTag-hb mRNAs. Whereas the proportion of SunTag-hb mRNAs translated is initially uniform, translation declines from the anterior over time until it becomes restricted to a posterior band in the expression domain. We discuss how regulated hb mRNA translation may help establish the sharp Hb expression boundary, which is a model for precision and noise during developmental patterning. Overall, our data show how use of the SunTag method on fixed and live embryos is a powerful combination for elucidating spatiotemporal regulation of mRNA translation in Drosophila.

Cost-precision trade-off relation determines the optimal morphogen gradient for accurate biological pattern formation

Spatial boundaries formed during animal development originate from the pre-patterning of tissues by signaling molecules, called morphogens. The accuracy of boundary location is limited by the fluctuations of morphogen concentration that thresholds the expression level of target gene. Producing more morphogen molecules, which gives rise to smaller relative fluctuations, would better serve to shape more precise target boundaries; however, it incurs more thermodynamic cost. In the classical diffusion-depletion model of morphogen profile formation, the morphogen molecules synthesized from a local source display an exponentially decaying concentration profile with a characteristic length λ. Our theory suggests that in order to attain a precise profile with the minimal cost, λ should be roughly half the distance to the target boundary position from the source. Remarkably, we find that the profiles of morphogens that pattern the Drosophila embryo and wing imaginal disk are formed with nearly optimal λ. Our finding underscores the cost-effectiveness of precise morphogen profile formation in Drosophila development.

Shikonin Inhibits Tumor Growth of ESCC by suppressing PKM2 mediated Aerobic Glycolysis and STAT3 Phosphorylation

Background: Shikonin, a small molecule inhibitor of pyruvate kinase 2 (PKM2), has been demonstrated to play the antitumor effect in various cancers. However, the specific effects and related regulatory mechanism of Shikonin in esophageal squamous cell carcinoma (ESCC) have not been clearly declared.

Materials and methods: Cell viability was valued through 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay. Glucose consumption, lactate production, glycolytic intermediates and pyruvate kinase enzymatic activity were measured using corresponding assay kits. Patient-derived xenograft (PDX) models were constructed to observe the anti-ESCC effect of Shikonin in vivo. PKM2, p-PKM2, signal transducer and activator of transcription 3 (STAT3), p-STAT3, glucose transporter 1 (GLUT1) and hexokinase 2 (HK2) in ESCC tissues were assessed by western blot. The expression of PKM2, p-PKM2, p-STAT3, GLUT1 and HK2 was assessed by immunohistochemistry (IHC) in ESCC tissue based on PDXs.

Results: Shikonin effectively inhibited cell proliferation in dose-dependent and time-dependent manner compared with the control group. The detection of glycolysis showed that Shikonin suppressed the glucose consumption, lactate production, glycolytic intermediates and pyruvate kinase enzymatic activity. Furthermore, Shikonin not only inhibited the growth of ESCC, but also decreased the expression of p-PKM2 and p-STAT3 in vivo. Finally, Shikonin suppressed the expression of GLUT1 and HK2 proteins which are related to glycolysis.

Conclusion: Shikonin has a significant antitumor effect in the ESCC by suppressing PKM2 mediated aerobic glycolysis and regulating PKM2/STAT3 signal pathway.

Dynamic regulation of anterior-posterior patterning genes in living Drosophila embryos

Expression of the gap and pair-rule genes plays an essential role in body segmentation during Drosophila embryogenesis.^1-5 However, it remains unclear how precise expression patterns of these key developmental genes arise from stochastic transcriptional activation at the single-cell level. Here, I employed genome-editing and live-imaging approaches to comprehensively visualize regulation of the gap and pair-rule genes at the endogenous loci. Quantitative image analysis revealed that the total duration of active transcription (transcription period) is a major determinant of spatial patterning of gene expression in early embryos. The length of the transcription period is determined by the continuity of bursting activities in individual nuclei, with the core expression domain producing more bursts than boundary regions. Each gene exhibits a distinct rate of nascent RNA production during transcriptional bursting, which contributes to gene-to-gene variability in the total output. I also provide evidence for "enhancer interference," wherein a distal weak enhancer interferes with transcriptional activation by a strong proximal enhancer to downregulate the length of the transcription period without changing the transcription rate. Analysis of the endogenous hunchback (hb) locus revealed that the removal of the distal shadow enhancer induces strong ectopic transcriptional activation, which suppresses refinement of the initial broad expression domain into narrower stripe patterns at the anterior part of embryos. This study provides key insights into the link between transcriptional bursting, enhancer-promoter interaction, and spatiotemporal patterning of gene expression during animal development.

Real-time single-cell characterization of the eukaryotic transcription cycle reveals correlations between RNA initiation, elongation, and cleavage

The eukaryotic transcription cycle consists of three main steps: initiation, elongation, and cleavage of the nascent RNA transcript. Although each of these steps can be regulated as well as coupled with each other, their in vivo dissection has remained challenging because available experimental readouts lack sufficient spatiotemporal resolution to separate the contributions from each of these steps. Here, we describe a novel application of Bayesian inference techniques to simultaneously infer the effective parameters of the transcription cycle in real time and at the single-cell level using a two-color MS2/PP7 reporter gene and the developing fruit fly embryo as a case study. Our method enables detailed investigations into cell-to-cell variability in transcription-cycle parameters as well as single-cell correlations between these parameters. These measurements, combined with theoretical modeling, suggest a substantial variability in the elongation rate of individual RNA polymerase molecules. We further illustrate the power of this technique by uncovering a novel mechanistic connection between RNA polymerase density and nascent RNA cleavage efficiency. Thus, our approach makes it possible to shed light on the regulatory mechanisms in play during each step of the transcription cycle in individual, living cells at high spatiotemporal resolution.

Dynamic patterning by morphogens illuminated by cis-regulatory studies

Morphogen concentration changes in space as well as over time during development. However, how these dynamics are interpreted by cells to specify fate is not well understood. Here, we focus on two morphogens: the maternal transcription factors Bicoid and Dorsal, which directly regulate target genes to pattern Drosophila embryos. The actions of these factors at enhancers has been thoroughly dissected and provides a rich platform for understanding direct input by morphogens and their changing roles over time. Importantly, Bicoid and Dorsal do not work alone; we also discuss additional inputs that work with morphogens to control spatiotemporal gene expression in embryos.

Quantitative dissection of transcription in development yields evidence for transcription-factor-driven chromatin accessibility

Thermodynamic models of gene regulation can predict transcriptional regulation in bacteria, but in eukaryotes, chromatin accessibility and energy expenditure may call for a different framework. Here, we systematically tested the predictive power of models of DNA accessibility based on the Monod-Wyman-Changeux (MWC) model of allostery, which posits that chromatin fluctuates between accessible and inaccessible states. We dissected the regulatory dynamics of hunchback by the activator Bicoid and the pioneer-like transcription factor Zelda in living Drosophila embryos and showed that no thermodynamic or non-equilibrium MWC model can recapitulate hunchback transcription. Therefore, we explored a model where DNA accessibility is not the result of thermal fluctuations but is catalyzed by Bicoid and Zelda, possibly through histone acetylation, and found that this model can predict hunchback dynamics. Thus, our theory-experiment dialogue uncovered potential molecular mechanisms of transcriptional regulatory dynamics, a key step toward reaching a predictive understanding of developmental decision-making.

Transcriptional Enhancers in Drosophila

Key discoveries in Drosophila have shaped our understanding of cellular "enhancers." With a special focus on the fly, this chapter surveys properties of these adaptable cis-regulatory elements, whose actions are critical for the complex spatial/temporal transcriptional regulation of gene expression in metazoa. The powerful combination of genetics, molecular biology, and genomics available in Drosophila has provided an arena in which the developmental role of enhancers can be explored. Enhancers are characterized by diverse low- or high-throughput assays, which are challenging to interpret, as not all of these methods of identifying enhancers produce concordant results. As a model metazoan, the fly offers important advantages to comprehensive analysis of the central functions that enhancers play in gene expression, and their critical role in mediating the production of phenotypes from genotype and environmental inputs. A major challenge moving forward will be obtaining a quantitative understanding of how these cis-regulatory elements operate in development and disease.

Distinct Roles of Broadly Expressed Repressors Support Dynamic Enhancer Action and Change in Time

How broadly expressed repressors regulate gene expression is incompletely understood. To gain insight, we investigated how Suppressor of Hairless-Su(H)-and Runt regulate expression of bone morphogenetic protein (BMP) antagonist short-gastrulation via the sog_Distal enhancer. A live imaging protocol was optimized to capture this enhancer's spatiotemporal output throughout the early Drosophila embryo, finding in this context that Runt regulates transcription initiation, Su(H) regulates transcription rate, and both factors control spatial expression. Furthermore, whereas Su(H) functions as a dedicated repressor, Runt temporally switches from repressor to activator. Our results demonstrate that broad repressors play temporally distinct roles and contribute to dynamic gene expression. Both Run and Su(H)'s ability to influence the spatiotemporal domains of gene expression may serve to counterbalance activators and function in this manner as important regulators of the maternal-to-zygotic transition in early embryos.

Shadow enhancers can suppress input transcription factor noise through distinct regulatory logic

Shadow enhancers, groups of seemingly redundant enhancers, are found in a wide range of organisms and are critical for robust developmental patterning. However, their mechanism of action is unknown. We hypothesized that shadow enhancers drive consistent expression levels by buffering upstream noise through a separation of transcription factor (TF) inputs at the individual enhancers. By measuring the transcriptional dynamics of several Kruppel shadow enhancer configurations in live Drosophila embryos, we showed that individual member enhancers act largely independently. We found that TF fluctuations are an appreciable source of noise that the shadow enhancer pair can better buffer than duplicated enhancers. The shadow enhancer pair is also uniquely able to maintain low levels of expression noise across a wide range of temperatures. A stochastic model demonstrated the separation of TF inputs is sufficient to explain these findings. Our results suggest the widespread use of shadow enhancers is partially due to their noise suppressing ability.

In all higher organisms, life begins with a single cell. During the early stages of development, this single cell grows and divides multiple times to develop into the many different kinds of cells that make up an organism. This is a highly regulated process during which cells receive instructions telling them what kind of cell to become.

These instructions are relayed via genes, and a particular combination of activated genes determines the cell’s fate. Specific pieces of DNA, known as enhancers, act as switches that control when and where genes are active, while so-called shadow enhancers are found in groups and work together to turn on the same gene in a similar way.

Shadow enhancers are often active during the early stages of life to direct the formation of specialized cells in different parts of the body. But so far, it has been unclear why it is beneficial to the divide the role of activating genes across several shadow enhancers rather than a single one.

Here, Waymack et al. examined shadow enhancers around a gene called Kruppel in embryos of the fruit fly Drosophila melanogaster. Manipulating the shadow enhancers showed that they help to make gene activity more resistant to changes. Factors such as fluctuations in temperature have different effects on each shadow enhancer. Having several shadow enhancers working together ensures that, whatever happens, the right genes still get activated. For genes like Kruppel, which are key for healthy development, the ability to withstand unexpected changes is a valuable evolutionary benefit.

The study of Waymack et al. reveals why shadow enhancers are involved in the regulation of many genes, which may help to better understand developmental defects. Many conditions caused by such defects are influenced by both genetics and the environment. Genetic illnesses can vary in severity, which may be related to the roles of shadow enhancers. As such, studying shadow enhancers could lead to new approaches for treating genetic diseases.

A mechanism for hunchback promoters to readout morphogenetic positional information in less than a minute

Cell fate decisions in the fly embryo are rapid: hunchback genes decide in minutes whether nuclei follow the anterior/posterior developmental blueprint by reading out positional information in the Bicoid morphogen. This developmental system is a prototype of regulatory decision processes that combine speed and accuracy. Traditional arguments based on fixed-time sampling of Bicoid concentration indicate that an accurate readout is impossible within the experimental times. This raises the general issue of how speed-accuracy tradeoffs are achieved. Here, we compare fixed-time to on-the-fly decisions, based on comparing the likelihoods of anterior/posterior locations. We found that these more efficient schemes complete reliable cell fate decisions within the short embryological timescales. We discuss the influence of promoter architectures on decision times and error rates, present concrete examples that rapidly readout the morphogen, and predictions for new experiments. Lastly, we suggest a simple mechanism for RNA production and degradation that approximates the log-likelihood function.

Odd-paired is a pioneer-like factor that coordinates with Zelda to control gene expression in embryos

Pioneer factors such as Zelda (Zld) help initiate zygotic transcription in Drosophila early embryos, but whether other factors support this dynamic process is unclear. Odd-paired (Opa), a zinc-finger transcription factor expressed at cellularization, controls the transition of genes from pair-rule to segmental patterns along the anterior-posterior axis. Finding that Opa also regulates expression through enhancer sog_Distal along the dorso-ventral axis, we hypothesized Opa's role is more general. Chromatin-immunoprecipitation (ChIP-seq) confirmed its in vivo binding to sog_Distal but also identified widespread binding throughout the genome, comparable to Zld. Furthermore, chromatin assays (ATAC-seq) demonstrate that Opa, like Zld, influences chromatin accessibility genome-wide at cellularization, suggesting both are pioneer factors with common as well as distinct targets. Lastly, embryos lacking opa exhibit widespread, late patterning defects spanning both axes. Collectively, these data suggest Opa is a general timing factor and likely late-acting pioneer factor that drives a secondary wave of zygotic gene expression.

Tunable self-cleaving ribozymes for modulating gene expression in eukaryotic systems

Advancements in the field of synthetic biology have been possible due to the development of genetic tools that are able to regulate gene expression. However, the current toolbox of gene regulatory tools for eukaryotic systems have been outpaced by those developed for simple, single-celled systems. Here, we engineered a set of gene regulatory tools by combining self-cleaving ribozymes with various upstream competing sequences that were designed to disrupt ribozyme self-cleavage. As a proof-of-concept, we were able to modulate GFP expression in mammalian cells, and then showed the feasibility of these tools in Drosophila embryos. For each system, the fold-reduction of gene expression was influenced by the location of the self-cleaving ribozyme/upstream competing sequence (i.e. 5' vs. 3' untranslated region) and the competing sequence used. Together, this work provides a set of genetic tools that can be used to tune gene expression across various eukaryotic systems.

Constraints and limitations on the transcriptional response downstream of the Bicoid morphogen gradient

The regulation of the hunchback promoter expression by the maternal Bicoid gradient has been studied as a model system in development for many years. Yet, at the level of quantitative agreement between data and theoretical models, even the first step of this regulation, transcription, continues to be challenging. This situation is slowly progressing, thanks to quantitative live-imaging techniques coupled to advanced statistical data analysis and modeling. Here, we outline the current state of our knowledge of this apparently "simple" step, highlighting the newly appreciated role of bursty transcription dynamics and its regulation.

Lighting up the central dogma for predictive developmental biology

Although the last 30years have witnessed the mapping of the wiring diagrams of the gene regulatory networks that dictate cell fate and animal body plans, specific understanding building on such network diagrams that shows how DNA regulatory regions control gene expression lags far behind. These networks have yet to yield the predictive power necessary to, for example, calculate how the concentration dynamics of input transcription factors and DNA regulatory sequence prescribes output patterns of gene expression that, in turn, determine body plans themselves. Here, we argue that reaching a predictive understanding of developmental decision-making calls for an interplay between theory and experiment aimed at revealing how the regulation of the processes of the central dogma dictate network connections and how network topology guides cells toward their ultimate developmental fate. To make this possible, it is crucial to break free from the snapshot-based understanding of embryonic development facilitated by fixed-tissue approaches and embrace new technologies that capture the dynamics of developmental decision-making at the single cell level, in living embryos.

Dissecting the sharp response of a canonical developmental enhancer reveals multiple sources of cooperativity

Developmental enhancers integrate graded concentrations of transcription factors (TFs) to create sharp gene expression boundaries. Here we examine the hunchback P2 (HbP2) enhancer which drives a sharp expression pattern in the Drosophila blastoderm embryo in response to the transcriptional activator Bicoid (Bcd). We systematically interrogate cis and trans factors that influence the shape and position of expression driven by HbP2, and find that the prevailing model, based on pairwise cooperative binding of Bcd to HbP2 is not adequate. We demonstrate that other proteins, such as pioneer factors, Mediator and histone modifiers influence the shape and position of the HbP2 expression pattern. Comparing our results to theory reveals how higher-order cooperativity and energy expenditure impact boundary location and sharpness. Our results emphasize that the bacterial view of transcription regulation, where pairwise interactions between regulatory proteins dominate, must be reexamined in animals, where multiple molecular mechanisms collaborate to shape the gene regulatory function.

The regulatory control of Cebpa enhancers and silencers in the myeloid and red-blood cell lineages

Cebpa encodes a transcription factor (TF) that plays an instructive role in the development of multiple myeloid lineages. The expression of Cebpa itself is finely modulated, as Cebpa is expressed at high and intermediate levels in neutrophils and macrophages respectively and downregulated in non-myeloid lineages. The cis-regulatory logic underlying the lineage-specific modulation of Cebpa's expression level is yet to be fully characterized. Previously, we had identified 6 new cis-regulatory modules (CRMs) in a 78kb region surrounding Cebpa. We had also inferred the TFs that regulate each CRM by fitting a sequence-based thermodynamic model to a comprehensive reporter activity dataset. Here, we report the cis-regulatory logic of Cebpa CRMs at the resolution of individual binding sites. We tested the binding sites and functional roles of inferred TFs by designing and constructing mutated CRMs and comparing theoretical predictions of their activity against empirical measurements in a myeloid cell line. The enhancers were confirmed to be activated by combinations of PU.1, C/EBP family TFs, Egr1, and Gfi1 as predicted by the model. We show that silencers repress the activity of the proximal promoter in a dominant manner in G1ME cells, which are derived from the red-blood cell lineage. Dominant repression in G1ME cells can be traced to binding sites for GATA and Myb, a motif shared by all of the silencers. Finally, we demonstrate that GATA and Myb act redundantly to silence the proximal promoter. These results indicate that dominant repression is a novel mechanism for resolving hematopoietic lineages. Furthermore, Cebpa has a fail-safe cis-regulatory architecture, featuring several functionally similar CRMs, each of which contains redundant binding sites for multiple TFs. Lastly, by experimentally demonstrating the predictive ability of our sequence-based thermodynamic model, this work highlights the utility of this computational approach for understanding mammalian gene regulation.

Precision in a rush: Trade-offs between reproducibility and steepness of the hunchback expression pattern

Fly development amazes us by the precision and reproducibility of gene expression, especially since the initial expression patterns are established during very short nuclear cycles. Recent live imaging of hunchback promoter dynamics shows a stable steep binary expression pattern established within the three minute interphase of nuclear cycle 11. Considering expression models of different complexity, we explore the trade-off between the ability of a regulatory system to produce a steep boundary and minimize expression variability between different nuclei. We show how a limited readout time imposed by short developmental cycles affects the gene’s ability to read positional information along the embryo’s anterior posterior axis and express reliably. Comparing our theoretical results to real-time monitoring of the hunchback transcription dynamics in live flies, we discuss possible regulatory strategies, suggesting an important role for additional binding sites, gradients or non-equilibrium binding and modified transcription factor search strategies.

Despite very limited time, organisms develop in reproducible ways. In the early stages of fly development the information about maternal signals is read out in a few minutes to produce steep and precise gene expression patterns. Motivated by recent live imaging experiments in fly embryos, we explore the consequences of the trade-off between a rushed but reproducible readout and a steep expression pattern on the regulatory modules of gene expression. We show that the current view of one anterior gradient morphogen binding to six binding sites is quantitatively inconsistent with the experimental data given the short readout time, suggesting other regulatory features.

3 minutes to precisely measure morphogen concentration

Morphogen gradients provide concentration-dependent positional information along polarity axes. Although the dynamics of the establishment of these gradients is well described, precision and noise in the downstream activation processes remain elusive. A simple paradigm to address these questions is the Bicoid morphogen gradient that elicits a rapid step-like transcriptional response in young fruit fly embryos. Focusing on the expression of the major Bicoid target, hunchback (hb), at the onset of zygotic transcription, we used the MS2-MCP approach which combines fluorescent labeling of nascent mRNA with live imaging at high spatial and temporal resolution. Removing 36 putative Zelda binding sites unexpectedly present in the original MS2 reporter, we show that the 750 bp of the hb promoter are sufficient to recapitulate endogenous expression at the onset of zygotic transcription. After each mitosis, in the anterior, expression is turned on to rapidly reach a plateau with all nuclei expressing the reporter. Consistent with a Bicoid dose-dependent activation process, the time period required to reach the plateau increases with the distance to the anterior pole. Despite the challenge imposed by frequent mitoses and high nuclei-to-nuclei variability in transcription kinetics, it only takes 3 minutes at each interphase for the MS2 reporter loci to distinguish subtle differences in Bicoid concentration and establish a steadily positioned and steep (Hill coefficient ~ 7) expression boundary. Modeling based on the cooperativity between the 6 known Bicoid binding sites in the hb promoter region, assuming rate limiting concentrations of the Bicoid transcription factor at the boundary, is able to capture the observed dynamics of pattern establishment but not the steepness of the boundary. This suggests that a simple model based only on the cooperative binding of Bicoid is not sufficient to describe the spatiotemporal dynamics of early hb expression.

Heterosis Is a Systemic Property Emerging From Non-linear Genotype-Phenotype Relationships: Evidence From in Vitro Genetics and Computer Simulations

Heterosis, the superiority of hybrids over their parents for quantitative traits, represents a crucial issue in plant and animal breeding as well as evolutionary biology. Heterosis has given rise to countless genetic, genomic and molecular studies, but has rarely been investigated from the point of view of systems biology. We hypothesized that heterosis is an emergent property of living systems resulting from frequent concave relationships between genotypic variables and phenotypes, or between different phenotypic levels. We chose the enzyme-flux relationship as a model of the concave genotype-phenotype (GP) relationship, and showed that heterosis can be easily created in the laboratory. First, we reconstituted in vitro the upper part of glycolysis. We simulated genetic variability of enzyme activity by varying enzyme concentrations in test tubes. Mixing the content of "parental" tubes resulted in "hybrids," whose fluxes were compared to the parental fluxes. Frequent heterotic fluxes were observed, under conditions that were determined analytically and confirmed by computer simulation. Second, to test this model in a more realistic situation, we modeled the glycolysis/fermentation network in yeast by considering one input flux, glucose, and two output fluxes, glycerol and acetaldehyde. We simulated genetic variability by randomly drawing parental enzyme concentrations under various conditions, and computed the parental and hybrid fluxes using a system of differential equations. Again we found that a majority of hybrids exhibited positive heterosis for metabolic fluxes. Cases of negative heterosis were due to local convexity between certain enzyme concentrations and fluxes. In both approaches, heterosis was maximized when the parents were phenotypically close and when the distributions of parental enzyme concentrations were contrasted and constrained. These conclusions are not restricted to metabolic systems: they only depend on the concavity of the GP relationship, which is commonly observed at various levels of the phenotypic hierarchy, and could account for the pervasiveness of heterosis.

LlamaTags: A Versatile Tool to Image Transcription Factor Dynamics in Live Embryos

Embryonic cell fates are defined by transcription factors that are rapidly deployed, yet attempts to visualize these factors in vivo often fail because of slow fluorescent protein maturation. Here, we pioneer a protein tag, LlamaTag, which circumvents this maturation limit by binding mature fluorescent proteins, making it possible to visualize transcription factor concentration dynamics in live embryos. Implementing this approach in the fruit fly Drosophila melanogaster, we discovered stochastic bursts in the concentration of transcription factors that are correlated with bursts in transcription. We further used LlamaTags to show that the concentration of protein in a given nucleus heavily depends on transcription of that gene in neighboring nuclei; we speculate that this inter-nuclear signaling is an important mechanism for coordinating gene expression to delineate straight and sharp boundaries of gene expression. Thus, LlamaTags now make it possible to visualize the flow of information along the central dogma in live embryos.

Transcriptional precision and accuracy in development: from measurements to models and mechanisms

During development, genes are transcribed at specific times, locations and levels. In recent years, the emergence of quantitative tools has significantly advanced our ability to measure transcription with high spatiotemporal resolution in vivo. Here, we highlight recent studies that have used these tools to characterize transcription during development, and discuss the mechanisms that contribute to the precision and accuracy of the timing, location and level of transcription. We attempt to disentangle the discrepancies in how physicists and biologists use the term ‘precision' to facilitate interactions using a common language. We also highlight selected examples in which the coupling of mathematical modeling with experimental approaches has provided important mechanistic insights, and call for a more expansive use of mathematical modeling to exploit the wealth of quantitative data and advance our understanding of animal transcription.

Summary: This Review highlights how high-resolution quantitative tools and theoretical models have formed our current view of the mechanisms determining precision and accuracy in the timing, location and level of transcription in the Drosophila embryo.

On the importance of protein diffusion in biological systems: The example of the Bicoid morphogen gradient

Morphogens are proteins that form concentration gradients in embryos and developing tissues, where they act as postal codes, providing cells with positional information and allowing them to behave accordingly. Bicoid was the first discovered morphogen, and remains one of the most studied. It regulates segmentation in flies, forming a striking exponential gradient along the anterior-posterior axis of early Drosophila embryos, and activating the transcription of multiple target genes in a concentration-dependent manner. In this review, the work done by us and by others to characterize the mobility of Bicoid in D. melanogaster embryos is presented. The central role played by the diffusion of Bicoid in both the establishment of the gradient and the activation of target genes is discussed, and placed in the context of the need for these processes to be all at once rapid, precise and robust. The Bicoid system, and morphogen gradients in general, remain amongst the most amazing examples of the coexistence, often observed in living systems, of small-scale disorder and large-scale spatial order. This article is part of a Special Issue entitled: Biophysics in Canada, edited by Lewis Kay, John Baenziger, Albert Berghuis and Peter Tieleman.

Dense Bicoid hubs accentuate binding along the morphogen gradient

Morphogen gradients direct the spatial patterning of developing embryos; however, the mechanisms by which these gradients are interpreted remain elusive. Here we used lattice light-sheet microscopy to perform in vivo single-molecule imaging in early Drosophila melanogaster embryos of the transcription factor Bicoid that forms a gradient and initiates patterning along the anteroposterior axis. In contrast to canonical models, we observed that Bicoid binds to DNA with a rapid off rate throughout the embryo such that its average occupancy at target loci is on-rate-dependent. We further observed Bicoid forming transient "hubs" of locally high density that facilitate binding as factor levels drop, including in the posterior, where we observed Bicoid binding despite vanishingly low protein levels. We propose that localized modulation of transcription factor on rates via clustering provides a general mechanism to facilitate binding to low-affinity targets and that this may be a prevalent feature of other developmental transcription factors.

Composition and dosage of a multipartite enhancer cluster control developmental expression of Ihh (Indian hedgehog)

Copy number variations (CNVs) often include non-coding sequence and putative enhancers but how these rearrangements induce disease is poorly understood. Here we investigate CNVs involving the regulatory landscape of Indian hedgehog (IHH), causing multiple, highly localised phenotypes including craniosynostosis and synpolydactyly^1,2. We show through transgenic reporter and genome editing studies in mice that Ihh is regulated by a constellation of at least 9 enhancers with individual tissue specificities in the digit anlagen, growth plates, skull sutures and fingertips. Consecutive deletions show that they function in an additive manner resulting in growth defects of the skull and long bones. Duplications, in contrast, cause not only dose-dependent upregulation but also misexpression of Ihh, leading to abnormal phalanges, fusion of sutures and syndactyly. Thus, precise spatio-temporal control of developmental gene expression is achieved by complex multipartite enhancer ensembles. Alterations in the composition of such clusters can result in gene misexpression and disease.

Broadly expressed repressors integrate patterning across orthogonal axes in embryos

The role of spatially localized repressors in supporting embryonic patterning is well appreciated, but, alternatively, the role ubiquitously expressed repressors play in this process is not well understood. We investigated the function of two broadly expressed repressors, Runt (Run) and Suppressor of Hairless [Su(H)], in patterning the Drosophila embryo. Previous studies have shown that Run and Su(H) regulate gene expression along anterior-posterior (AP) or dorsal-ventral (DV) axes, respectively, by spatially limiting activator action, but here we characterize a different role. Our data show that broadly expressed repressors silence particular enhancers within cis-regulatory systems, blocking their expression throughout the embryo fully but transiently, and, in this manner, regulate spatiotemporal outputs along both axes. Our results suggest that Run and Su(H) regulate the temporal action of enhancers and are not dedicated regulators of one axis but, instead, act coordinately to pattern both axes, AP and DV.

Decoding temporal interpretation of the morphogen Bicoid in the early Drosophila embryo

Morphogen gradients provide essential spatial information during development. Not only the local concentration but also duration of morphogen exposure is critical for correct cell fate decisions. Yet, how and when cells temporally integrate signals from a morphogen remains unclear. Here, we use optogenetic manipulation to switch off Bicoid-dependent transcription in the early Drosophila embryo with high temporal resolution, allowing time-specific and reversible manipulation of morphogen signalling. We find that Bicoid transcriptional activity is dispensable for embryonic viability in the first hour after fertilization, but persistently required throughout the rest of the blastoderm stage. Short interruptions of Bicoid activity alter the most anterior cell fate decisions, while prolonged inactivation expands patterning defects from anterior to posterior. Such anterior susceptibility correlates with high reliance of anterior gap gene expression on Bicoid. Therefore, cell fates exposed to higher Bicoid concentration require input for longer duration, demonstrating a previously unknown aspect of Bicoid decoding.

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

Regulatory logic driving stable levels of defective proventriculus expression during terminal photoreceptor specification in flies

How differential levels of gene expression are controlled in post-mitotic neurons is poorly understood. In the Drosophila retina, expression of the transcription factor Defective Proventriculus (Dve) at distinct cell type-specific levels is required for terminal differentiation of color- and motion-detecting photoreceptors. Here, we find that the activities of two cis-regulatory enhancers are coordinated to drive dve expression in the fly eye. Three transcription factors act on these enhancers to determine cell-type specificity. Negative autoregulation by Dve maintains expression from each enhancer at distinct homeostatic levels. One enhancer acts as an inducible backup ('dark' shadow enhancer) that is normally repressed but becomes active in the absence of the other enhancer. Thus, two enhancers integrate combinatorial transcription factor input, feedback and redundancy to generate cell type-specific levels of dve expression and stable photoreceptor fate. This regulatory logic may represent a general paradigm for how precise levels of gene expression are established and maintained in post-mitotic neurons.

Gene expression boundary scaling and organ size regulation in the Drosophila embryo

How the shape and size of tissues and organs is regulated during development is a major question in developmental biology. Such regulation relies upon both intrinsic cues (such as signaling networks) and extrinsic inputs (such as from neighboring tissues). Here, we focus on pattern formation and organ development during Drosophila embryogenesis. In particular, we outline the importance of both biochemical and mechanical tissue-tissue interactions in size regulation. We describe how the Drosophila embryo can potentially provide novel insights into how shape and size are regulated during development. We focus on gene expression boundary scaling in the early embryo and how size is regulated in three organs (hindgut, trachea, and ventral nerve cord) later in development, with particular focus on the role of tissue-tissue interactions. Overall, we demonstrate that Drosophila embryogenesis provides a suitable model system for studying spatial and temporal scaling and size control in vivo.

Information Integration and Energy Expenditure in Gene Regulation

The quantitative concepts used to reason about gene regulation largely derive from bacterial studies. We show that this bacterial paradigm cannot explain the sharp expression of a canonical developmental gene in response to a regulating transcription factor (TF). In the absence of energy expenditure, with regulatory DNA at thermodynamic equilibrium, information integration across multiple TF binding sites can generate the required sharpness, but with strong constraints on the resultant "higher-order cooperativities." Even with such integration, there is a "Hopfield barrier" to sharpness; for n TF binding sites, this barrier is represented by the Hill function with the Hill coefficient n. If, however, energy is expended to maintain regulatory DNA away from thermodynamic equilibrium, as in kinetic proofreading, this barrier can be breached and greater sharpness achieved. Our approach is grounded in fundamental physics, leads to testable experimental predictions, and suggests how a quantitative paradigm for eukaryotic gene regulation can be formulated.

Forces directing germ-band extension in Drosophila embryos

Body axis elongation by convergent extension is a conserved developmental process found in all metazoans. Drosophila embryonic germ-band extension is an important morphogenetic process during embryogenesis, by which the length of the germ-band is more than doubled along the anterior-posterior axis. This lengthening is achieved by typical convergent extension, i.e. narrowing the lateral epidermis along the dorsal-ventral axis and simultaneous extension along the anterior-posterior axis. Germ-band extension is largely driven by cell intercalation, whose directionality is determined by the planar polarity of the tissue and ultimately by the anterior-posterior patterning system. In addition, extrinsic tensile forces originating from the invaginating endoderm induce cell shape changes, which transiently contribute to germ-band extension. Here, we review recent progress in understanding of the role of mechanical forces in germ-band extension.

Precision of Readout at the hunchback Gene: Analyzing Short Transcription Time Traces in Living Fly Embryos

The simultaneous expression of the hunchback gene in the numerous nuclei of the developing fly embryo gives us a unique opportunity to study how transcription is regulated in living organisms. A recently developed MS2-MCP technique for imaging nascent messenger RNA in living Drosophila embryos allows us to quantify the dynamics of the developmental transcription process. The initial measurement of the morphogens by the hunchback promoter takes place during very short cell cycles, not only giving each nucleus little time for a precise readout, but also resulting in short time traces of transcription. Additionally, the relationship between the measured signal and the promoter state depends on the molecular design of the reporting probe. We develop an analysis approach based on tailor made autocorrelation functions that overcomes the short trace problems and quantifies the dynamics of transcription initiation. Based on live imaging data, we identify signatures of bursty transcription initiation from the hunchback promoter. We show that the precision of the expression of the hunchback gene to measure its position along the anterior-posterior axis is low both at the boundary and in the anterior even at cycle 13, suggesting additional post-transcriptional averaging mechanisms to provide the precision observed in fixed embryos.

The fly embryo provides a natural laboratory to study the dynamics of transcription and its implications for the developing organism. Using live imaging experiments we investigate the nature of transcription regulation of the hunchback gene—the first to read out the maternal Bicoid gradient. While traditional time trace analysis methods based on OFF time distributions or autocorrelation functions fail for short signals, our tailored autocorrelation function overcomes these limitations revealing bursty dynamics that is reproducible between cell cycles and embryos. The inferred rates result in a lot of variability in the readout of nuclei sensing similar Bicoid concentrations, suggesting additional readout mechanisms than a one-to-one mapping of the input onto the output.

The Gene Regulatory Network of Lens Induction Is Wired through Meis-Dependent Shadow Enhancers of Pax6

Lens induction is a classical developmental model allowing investigation of cell specification, spatiotemporal control of gene expression, as well as how transcription factors are integrated into highly complex gene regulatory networks (GRNs). Pax6 represents a key node in the gene regulatory network governing mammalian lens induction. Meis1 and Meis2 homeoproteins are considered as essential upstream regulators of Pax6 during lens morphogenesis based on their interaction with the ectoderm enhancer (EE) located upstream of Pax6 transcription start site. Despite this generally accepted regulatory pathway, Meis1-, Meis2- and EE-deficient mice have surprisingly mild eye phenotypes at placodal stage of lens development. Here, we show that simultaneous deletion of Meis1 and Meis2 in presumptive lens ectoderm results in arrested lens development in the pre-placodal stage, and neither lens placode nor lens is formed. We found that in the presumptive lens ectoderm of Meis1/Meis2 deficient embryos Pax6 expression is absent. We demonstrate using chromatin immunoprecipitation (ChIP) that in addition to EE, Meis homeoproteins bind to a remote, ultraconserved SIMO enhancer of Pax6. We further show, using in vivo gene reporter analyses, that the lens-specific activity of SIMO enhancer is dependent on the presence of three Meis binding sites, phylogenetically conserved from man to zebrafish. Genetic ablation of EE and SIMO enhancers demostrates their requirement for lens induction and uncovers an apparent redundancy at early stages of lens development. These findings identify a genetic requirement for Meis1 and Meis2 during the early steps of mammalian eye development. Moreover, they reveal an apparent robustness in the gene regulatory mechanism whereby two independent "shadow enhancers" maintain critical levels of a dosage-sensitive gene, Pax6, during lens induction.

Maternal AP determinants in the Drosophila oocyte and embryo

An animal embryo cannot initiate its journey of forming a new life on its own. It must rely on maternally provided resources and inputs to kick-start its developmental process. In Drosophila, the initial polarities of the embryo along both the anterior-posterior (AP) and dorsal-ventral (DV) axes are also specified by maternal determinants. Over the past several decades, genetic and molecular studies have identified and characterized such determinants, as well as the zygotic genetic regulatory networks that control patterning in the early embryo. Extensive studies of oogenesis have also led to a detailed knowledge of the cellular and molecular interactions that control the formation of a mature egg. Despite these efforts, oogenesis and embryogenesis have been studied largely as separate problems, except for qualitative aspects with regard to maternal regulation of the asymmetric localization of maternal determinants. Can oogenesis and embryogenesis be viewed from a unified perspective at a quantitative level, and can that improve our understanding of how robust embryonic patterning is achieved? Here, we discuss the basic knowledge of the regulatory mechanisms controlling oogenesis and embryonic patterning along the AP axis. We explore properties of the maternal Bicoid gradient in relation to embryo size in search for a unified framework for robust AP patterning. WIREs Dev Biol 2016, 5:562-581. doi: 10.1002/wdev.235 For further resources related to this article, please visit the WIREs website.

Time to move on: Modeling transcription dynamics during an embryonic transition away from maternal control

In a recent study, we investigated the regulation of hunchback (hb) transcription dynamics in Drosophila embryos. Our results suggest that shutdown of hb transcription at early nuclear cycle (nc) 14 is an event associated with the global changes taking place during the mid-blastula transition (MBT). Here we have developed a simple model of hb transcription dynamics during this transition time. With kinetic parameters estimated from our published experimental data, the model describes the dynamical processes of hb gene transcription and hb mRNA accumulation. With two steps, transcription onset upon exiting the previous mitosis followed by a sudden impact that blocks gene activation, the model recapitulates the observed dynamics of hb transcription during the nc14 interphase. The timing of gene inactivation is essential, as its alterations lead to changes in both hb transcription dynamics and hb mRNA levels. Our model provides a clear dynamical picture of hb transcription regulation as one of the many, actively regulated events concurrently taking place during the MBT.

Sequential construction of a model for modular gene expression control, applied to spatial patterning of theDrosophilagenehunchback

Gene network simulations are increasingly used to quantify mutual gene regulation in biological tissues. These are generally based on linear interactions between single-entity regulatory and target genes. Biological genes, by contrast, commonly have multiple, partially independent, cis-regulatory modules (CRMs) for regulator binding, and can produce variant transcription and translation products. We present a modeling framework to address some of the gene regulatory dynamics implied by this biological complexity. Spatial patterning of the hunchback (hb) gene in Drosophila development involves control by three CRMs producing two distinct mRNA transcripts. We use this example to develop a differential equations model for transcription which takes into account the cis-regulatory architecture of the gene. Potential regulatory interactions are screened by a genetic algorithms (GAs) approach and compared to biological expression data.

Hill function-based models of transcriptional switches: impact of specific, nonspecific, functional and nonfunctional binding

We explore minimalist models of transcription in which we take into account that a cis-regulatory sequence is embedded in, and interacts with, a complex genome. The classical Hill equation is the simplest way to represent a transcriptional response. However, it may overlook the fact that a transcription factor (TF) establishes specific and nonspecific nonfunctional interactions with chromatin. Classical papers have shown that nonfunctional binding (not leading to transcription) may influence gene expression. We examine how the presence of additional binding sites for a TF, besides those on the gene(s) of interest, affect the shape and parameters of the transcriptional response. We consider two conditions: at equilibrium and at steady-state. In many cases the TF level is determined by the position of the cell within a spatial or temporal gradient. We show that such gradients can be adjusted by evolutionary selection to compensate for the alteration of the gene transcription response by the presence of nonfunctional binding sites. Finally, we analyse how the transcriptional response is affected by a decrease in TF concentration, as in cases of haploinsufficiency. We show that the nonlinearity of the transcriptional response as a function of [TF] exacerbates the effect of a decrease in the latter, at least for weakly expressed TFs. Although decades of work on TFs have led to the impression that almost everything is known about the control of gene expression, we show that even the simplest models of transcription control have not delivered all their secrets yet.

New methods to image transcription in living fly embryos: the insights so far, and the prospects

The regulation of transcription is a fundamental process underlying the determination of cell identity and its maintenance during development. In the last decades, most of the transcription factors, which have to be expressed at the right place and at the right time for the proper development of the fly embryo, have been identified. However, mostly because of the lack of methods to visualize transcription as the embryo develops, their coordinated spatiotemporal dynamics remains largely unexplored. Efforts have been made to decipher the transcription process with single molecule resolution at the single cell level. Recently, the fluorescent labeling of nascent RNA in developing fly embryos allowed the direct visualization of ongoing transcription at single loci within each nucleus. Together with powerful imaging and quantitative data analysis, these new methods provide unprecedented insights into the temporal dynamics of the transcription process and its intrinsic noise. Focusing on the Drosophila embryo, we discuss how the detection of single RNA molecules enhanced our comprehension of the transcription process and we outline the potential next steps made possible by these new imaging tools. In combination with genetics and theoretical analysis, these new imaging methods will aid the search for the mechanisms responsible for the robustness of development. For further resources related to this article, please visit the WIREs website.

Gap Gene Regulatory Dynamics Evolve along a Genotype Network

Developmental gene networks implement the dynamic regulatory mechanisms that pattern and shape the organism. Over evolutionary time, the wiring of these networks changes, yet the patterning outcome is often preserved, a phenomenon known as “system drift.” System drift is illustrated by the gap gene network—involved in segmental patterning—in dipteran insects. In the classic model organism Drosophila melanogaster and the nonmodel scuttle fly Megaselia abdita, early activation and placement of gap gene expression domains show significant quantitative differences, yet the final patterning output of the system is essentially identical in both species. In this detailed modeling analysis of system drift, we use gene circuits which are fit to quantitative gap gene expression data in M. abdita and compare them with an equivalent set of models from D. melanogaster. The results of this comparative analysis show precisely how compensatory regulatory mechanisms achieve equivalent final patterns in both species. We discuss the larger implications of the work in terms of “genotype networks” and the ways in which the structure of regulatory networks can influence patterns of evolutionary change (evolvability).