Transcriptional regulation of a pair-rule stripe in Drosophila

Syntax compensates for poor binding sites to encode tissue specificity of developmental enhancers

Transcriptional enhancers are short segments of DNA that switch genes on and off in response to a variety of intrinsic and extrinsic signals. Despite the discovery of the first enhancer more than 30 y ago, the relationship between primary DNA sequence and enhancer activity remains obscure. In particular, the importance of "syntax" (the order, orientation, and spacing of binding sites) is unclear. A high-throughput screen identified synthetic notochord enhancers that are activated by the combination of ZicL and ETS transcription factors in Ciona embryos. Manipulation of these enhancers elucidated a "regulatory code" of sequence and syntax features for notochord-specific expression. This code enabled in silico discovery of bona fide notochord enhancers, including those containing low-affinity binding sites that would be excluded by standard motif identification methods. One of the newly identified enhancers maps upstream of the known enhancer that regulates Brachyury (Ci-Bra), a key determinant of notochord specification. This newly identified Ci-Bra shadow enhancer contains binding sites with very low affinity, but optimal syntax, and therefore mediates surprisingly strong expression in the notochord. Weak binding sites are compensated by optimal syntax, whereas enhancers containing high-affinity binding affinities possess suboptimal syntax. We suggest this balance has obscured the importance of regulatory syntax, as noncanonical binding motifs are typically disregarded by enhancer detection methods. As a result, enhancers with low binding affinities but optimal syntax may be a vastly underappreciated feature of the regulatory genome.

The Wnt and Delta-Notch signalling pathways interact to direct pair-rule gene expression via caudal during segment addition in the spider Parasteatoda tepidariorum

In short germ arthropods, posterior segments are added sequentially from a growth zone or segment addition zone (SAZ) during embryogenesis. Studies in spiders such as the common house spider, Parasteatoda tepidariorum, have provided insights into the gene regulatory network (GRN) that underlies the development of the SAZ, and revealed the involvement of two important signalling pathways. It was shown that Wnt8 maintains a pool of undifferentiated cells in the SAZ, but this ligand is also required for dynamic Delta (Dl) expression associated with the formation of new segments. However, it remains unclear how these pathways interact during SAZ formation and subsequently regulate segment addition. Here we show that Delta-Notch signalling is required for Wnt8 expression in posterior SAZ cells, but represses the expression of this Wnt gene in anterior SAZ cells. We also found that these two signalling pathways are required for the expression of the spider orthologues of the segmentation genes even-skipped (eve) and runt-1 (run-1), at least in part via the transcription factor encoded by caudal (cad). Moreover, it appears that dynamic expression of eve in this spider does not require a feedback loop with run-1, as is found in the pair-rule circuit of the beetle Tribolium. Taken together, our results suggest that the development of posterior segments in Parasteatoda is directed by dynamic interactions between Wnt8 and Delta-Notch signalling that are read out by cad, which is necessary but not sufficient to regulate the expression of the pair-rule genes eve and run-1. Our study therefore provides new insights towards better understanding the evolution and developmental regulation of segmentation in other arthropods including insects.

Tissue patterning and cellular mechanics

In development, cells organize into biological tissues through cell growth, migration, and differentiation. Globally, this process is dictated by a genetically encoded program in which secreted morphogens and cell-cell interactions prompt the adoption of unique cell fates. Yet, at its lowest level, development is achieved through the modification of cell-cell adhesion and actomyosin-based contractility, which set the level of tension within cells and dictate how they pack together into tissues. The regulation of tension within individual cells and across large groups of cells is a major driving force of tissue organization and the basis of all cell shape change and cell movement in development.

Suboptimization of developmental enhancers

Transcriptional enhancers direct precise on-off patterns of gene expression during development. To explore the basis for this precision, we conducted a high-throughput analysis of the Otx-a enhancer, which mediates expression in the neural plate of Ciona embryos in response to fibroblast growth factor (FGF) signaling and a localized GATA determinant. We provide evidence that enhancer specificity depends on submaximal recognition motifs having reduced binding affinities ("suboptimization"). Native GATA and ETS (FGF) binding sites contain imperfect matches to consensus motifs. Perfect matches mediate robust but ectopic patterns of gene expression. The native sites are not arranged at optimal intervals, and subtle changes in their spacing alter enhancer activity. Multiple tiers of enhancer suboptimization produce specific, but weak, patterns of expression, and we suggest that clusters of weak enhancers, including certain "superenhancers," circumvent this trade-off in specificity and activity.

A gene expression atlas of a bicoid-depleted Drosophila embryo reveals early canalization of cell fate

In developing embryos, gene regulatory networks drive cells towards discrete terminal fates, a process called canalization. We studied the behavior of the anterior-posterior segmentation network in Drosophila melanogaster embryos by depleting a key maternal input, bicoid (bcd), and measuring gene expression patterns of the network at cellular resolution. This method results in a gene expression atlas containing the levels of mRNA or protein expression of 13 core patterning genes over six time points for every cell of the blastoderm embryo. This is the first cellular resolution dataset of a genetically perturbed Drosophila embryo that captures all cells in 3D. We describe the technical developments required to build this atlas and how the method can be employed and extended by others. We also analyze this novel dataset to characterize the degree and timing of cell fate canalization in the segmentation network. We find that in two layers of this gene regulatory network, following depletion of bcd, individual cells rapidly canalize towards normal cell fates. This result supports the hypothesis that the segmentation network directly canalizes cell fate, rather than an alternative hypothesis whereby cells are initially mis-specified and later eliminated by apoptosis. Our gene expression atlas provides a high resolution picture of a classic perturbation and will enable further computational modeling of canalization and gene regulation in this transcriptional network.

Shadow enhancers enable Hunchback bifunctionality in the Drosophila embryo

Hunchback (Hb) is a bifunctional transcription factor that activates and represses distinct enhancers. Here, we investigate the hypothesis that Hb can activate and repress the same enhancer. Computational models predicted that Hb bifunctionally regulates the even-skipped (eve) stripe 3+7 enhancer (eve3+7) in Drosophila blastoderm embryos. We measured and modeled eve expression at cellular resolution under multiple genetic perturbations and found that the eve3+7 enhancer could not explain endogenous eve stripe 7 behavior. Instead, we found that eve stripe 7 is controlled by two enhancers: the canonical eve3+7 and a sequence encompassing the minimal eve stripe 2 enhancer (eve2+7). Hb bifunctionally regulates eve stripe 7, but it executes these two activities on different pieces of regulatory DNA--it activates the eve2+7 enhancer and represses the eve3+7 enhancer. These two "shadow enhancers" use different regulatory logic to create the same pattern.

Estimating binding properties of transcription factors from genome-wide binding profiles

The binding of transcription factors (TFs) is essential for gene expression. One important characteristic is the actual occupancy of a putative binding site in the genome. In this study, we propose an analytical model to predict genomic occupancy that incorporates the preferred target sequence of a TF in the form of a position weight matrix (PWM), DNA accessibility data (in the case of eukaryotes), the number of TF molecules expected to be bound specifically to the DNA and a parameter that modulates the specificity of the TF. Given actual occupancy data in the form of ChIP-seq profiles, we backwards inferred copy number and specificity for five Drosophila TFs during early embryonic development: Bicoid, Caudal, Giant, Hunchback and Kruppel. Our results suggest that these TFs display thousands of molecules that are specifically bound to the DNA and that whilst Bicoid and Caudal display a higher specificity, the other three TFs (Giant, Hunchback and Kruppel) display lower specificity in their binding (despite having PWMs with higher information content). This study gives further weight to earlier investigations into TF copy numbers that suggest a significant proportion of molecules are not bound specifically to the DNA.

Deciphering the combinatorial architecture of a Drosophila homeotic gene enhancer

In Drosophila, the 330 kb bithorax complex regulates cellular differentiation along the anterior–posterior axis during development in the thorax and abdomen and is comprised of three homeotic genes: Ultrabithorax, abdominal-A, and Abdominal-B. The expression of each of these genes is in turn controlled through interactions between transcription factors and a number of cis-regulatory modules in the neighboring intergenic regions. In this study, we examine how the sequence architecture of transcription factor binding sites mediates the functional activity of one of these cis-regulatory modules. Using computational, mathematical modeling and experimental molecular genetic approaches we investigate the IAB7b enhancer, which regulates Abdominal-B expression specifically in the presumptive seventh and ninth abdominal segments of the early embryo. A cross-species comparison of the IAB7b enhancer reveals an evolutionarily conserved signature motif containing two FUSHI-TARAZU activator transcription factor binding sites. We find that the transcriptional repressors KNIRPS, KRUPPEL and GIANT are able to restrict reporter gene expression to the posterior abdominal segments, using different molecular mechanisms including short-range repression and competitive binding. Additionally, we show the functional importance of the spacing between the two FUSHI-TARAZU binding sites and discuss the potential importance of cooperativity for transcriptional activation. Our results demonstrate that the transcriptional output of the IAB7b cis-regulatory module relies on a complex set of combinatorial inputs mediated by specific transcription factor binding and that the sequence architecture at this enhancer is critical to maintain robust regulatory function.

Dynamic regulation of eve stripe 2 expression reveals transcriptional bursts in living Drosophila embryos

We present the use of recently developed live imaging methods to examine the dynamic regulation of even-skipped (eve) stripe 2 expression in the precellular Drosophila embryo. Nascent transcripts were visualized via MS2 RNA stem loops. The eve stripe 2 transgene exhibits a highly dynamic pattern of de novo transcription, beginning with a broad domain of expression during nuclear cycle 12 (nc12), and progressive refinement during nc13 and nc14. The mature stripe 2 pattern is surprisingly transient, constituting just ∼15 min of the ∼90-min period of expression. Nonetheless, this dynamic transcription profile faithfully predicts the limits of the mature stripe visualized by conventional in situ detection methods. Analysis of individual transcription foci reveals intermittent bursts of de novo transcription, with duration cycles of 4-10 min. We discuss a multistate model of transcription regulation and speculate on its role in the dynamic repression of the eve stripe 2 expression pattern during development.

Naturally occurring deletions of hunchback binding sites in the even-skipped stripe 3+7 enhancer

Changes in regulatory DNA contribute to phenotypic differences within and between taxa. Comparative studies show that many transcription factor binding sites (TFBS) are conserved between species whereas functional studies reveal that some mutations segregating within species alter TFBS function. Consistently, in this analysis of 13 regulatory elements in Drosophila melanogaster populations, single base and insertion/deletion polymorphism are rare in characterized regulatory elements. Experimentally defined TFBS are nearly devoid of segregating mutations and, as has been shown before, are quite conserved. For instance 8 of 11 Hunchback binding sites in the stripe 3+7 enhancer of even-skipped are conserved between D. melanogaster and Drosophila virilis. Oddly, we found a 72 bp deletion that removes one of these binding sites (Hb8), segregating within D. melanogaster. Furthermore, a 45 bp deletion polymorphism in the spacer between the stripe 3+7 and stripe 2 enhancers, removes another predicted Hunchback site. These two deletions are separated by ∼250 bp, sit on distinct haplotypes, and segregate at appreciable frequency. The Hb8Δ is at 5 to 35% frequency in the new world, but also shows cosmopolitan distribution. There is depletion of sequence variation on the Hb8Δ-carrying haplotype. Quantitative genetic tests indicate that Hb8Δ affects developmental time, but not viability of offspring. The Eve expression pattern differs between inbred lines, but the stripe 3 and 7 boundaries seem unaffected by Hb8Δ. The data reveal segregating variation in regulatory elements, which may reflect evolutionary turnover of characterized TFBS due to drift or co-evolution.

Ancestral resurrection of the Drosophila S2E enhancer reveals accessible evolutionary paths through compensatory change

Upstream regulatory sequences that control gene expression evolve rapidly, yet the expression patterns and functions of most genes are typically conserved. To address this paradox, we have reconstructed computationally and resurrected in vivo the cis-regulatory regions of the ancestral Drosophila eve stripe 2 element and evaluated its evolution using a mathematical model of promoter function. Our feed-forward transcriptional model predicts gene expression patterns directly from enhancer sequence. We used this functional model along with phylogenetics to generate a set of possible ancestral eve stripe 2 sequences for the common ancestors of 1) D. simulans and D. sechellia; 2) D. melanogaster, D. simulans, and D. sechellia; and 3) D. erecta and D. yakuba. These ancestral sequences were synthesized and resurrected in vivo. Using a combination of quantitative and computational analysis, we find clear support for functional compensation between the binding sites for Bicoid, Giant, and Krüppel over the course of 40-60 My of Drosophila evolution. We show that this compensation is driven by a coupling interaction between Bicoid activation and repression at the anterior and posterior border necessary for proper placement of the anterior stripe 2 border. A multiplicity of mechanisms for binding site turnover exemplified by Bicoid, Giant, and Krüppel sites, explains how rapid sequence change may occur while maintaining the function of the cis-regulatory element.

Cellular resolution models for even skipped regulation in the entire Drosophila embryo

Transcriptional control ensures genes are expressed in the right amounts at the correct times and locations. Understanding quantitatively how regulatory systems convert input signals to appropriate outputs remains a challenge. For the first time, we successfully model even skipped (eve) stripes 2 and 3+7 across the entire fly embryo at cellular resolution. A straightforward statistical relationship explains how transcription factor (TF) concentrations define eve's complex spatial expression, without the need for pairwise interactions or cross-regulatory dynamics. Simulating thousands of TF combinations, we recover known regulators and suggest new candidates. Finally, we accurately predict the intricate effects of perturbations including TF mutations and misexpression. Our approach imposes minimal assumptions about regulatory function; instead we infer underlying mechanisms from models that best fit the data, like the lack of TF-specific thresholds and the positional value of homotypic interactions. Our study provides a general and quantitative method for elucidating the regulation of diverse biological systems. DOI:http://dx.doi.org/10.7554/eLife.00522.001.

Homeogenetic inductive mechanism of segmentation in polychaete tail regeneration

Segmentation is a body-patterning strategy in which new segments are specified from a segment-addition zone containing uncommitted cells. However, the cell-recruitment process is poorly understood. Here we investigated in detail the segmentation in a polychaete annelid, Perinereis nuntia (Lophotrochozoa), in which new segments emerge at the boundary between the posterior end of the segmented region and the terminal pygidium. Cells at this border synchronously remodel their chromatin, enter the cell cycle, and undergo oriented cell division, before being added to new segments. wingless is expressed at the posterior edge of the pre-existing segment, abutted by hedgehog in the first row of the new segment. Overstimulation of Wingless signaling caused excess cells to enter the cell cycle, prolonging segmentation and widening the new segment. Thus, segment addition may occur by a homeogenetic mechanism, in which Wingless expressed in the differentiated segment coordinates the stepwise recruitment of undifferentiated cells from the segment/pygidium boundary.

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.

Metazoan genes, including those of humans, contain large noncoding regions that are required for viability. Sequence variations in these regions are statistically associated with human disease, but the mechanisms underlying these associations are not well understood. These regions regulate transcription and are frequently larger than the gene's transcript by an order of magnitude. In this paper we attempt to elucidate the regulatory code of these noncoding segments of DNA by means of quantitative spatially resolved gene expression data and a computational model. The expression data comes from the early embryo of the fruit fly Drosophila melanogaster. We chose a family of DNA constructs to analyze that drive very different patterns of expression when very small changes in DNA sequence are made, reasoning that this sensitivity would reveal important properties of the regulatory code. The model reproduced the training data with precision greater than the expected accuracy of the training data itself. It was able to correctly predict from DNA sequence the expression of 44 segments of DNA from many genes and species.

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.

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.

HOT DNAs: a novel class of developmental enhancers

Enhancers mediate localized patterns of gene expression during development. A common feature of "traditional" enhancers is the presence of clustered binding motifs for sequence-specific transcription factors (TFs). In this issue of Genes & Development, Kvon and colleagues (pp. 908-913) present new evidence that HOT (highly occupied transcription) DNAs direct specific patterns of gene expression, despite being depleted for TF-binding motifs.

A system of repressor gradients spatially organizes the boundaries of Bicoid-dependent target genes

The homeodomain (HD) protein Bicoid (Bcd) is thought to function as a gradient morphogen that positions boundaries of target genes via threshold-dependent activation mechanisms. Here, we analyze 66 Bcd-dependent regulatory elements and show that their boundaries are positioned primarily by repressive gradients that antagonize Bcd-mediated activation. A major repressor is the pair-rule protein Runt (Run), which is expressed in an opposing gradient and is necessary and sufficient for limiting Bcd-dependent activation. Evidence is presented that Run functions with the maternal repressor Capicua and the gap protein Kruppel as the principal components of a repression system that correctly orders boundaries throughout the anterior half of the embryo. These results put conceptual limits on the Bcd morphogen hypothesis and demonstrate how the Bcd gradient functions within the gene network that patterns the embryo.

Mechanical force alters morphogenetic movements and segmental gene expression patterns during Drosophila embryogenesis

The development of an organism is accompanied by various cellular morphogenetic movements, changes in cellular as well as nuclear morphology and transcription programs. Recent evidence suggests that intra and inter-cellular connections mediated by various adhesion proteins contribute to defining nuclear morphology. In addition, three dimensional organization of the cell nucleus regulate the transcription programs. However the link between cellular morphogenetic movements and its coupling to nuclear function in a developmental context is poorly understood. In this paper we use a point perturbation by tissue level laser ablation and sheet perturbation by application of force using magnetic tweezers to alter cellular morphogenetic movements and probe its impact on nuclear morphology and segmental gene expression patterns. Mechanical perturbations during blastoderm stage in a developing Drosophila embryo resulted in localized alterations in nuclear morphology and cellular movement. In addition, global defects in germ-band (GB) extension and retraction are observed when external force is applied during morphogenetic movements, suggesting a long-range physical coupling within the GB layer of cells. Further local application of force resulted in redistribution of non muscle myosin-II in the GB layer. Finally these perturbations lead to altered segmental gene (engrailed) expression patterns later during the development. Our observations suggest that there exists a tight regulation between nuclear morphology and cellular adhesive connections during morphogenetic movement of cells in the embryo. The observed spatial changes in patterning genes, with perturbation, highlight the importance of nuclear integrity to cellular movement in establishing gene expression program in a developmental system.

Estimation of pairwise sequence similarity of mammalian enhancers with word neighbourhood counts

MOTIVATION

The identity of cells and tissues is to a large degree governed by transcriptional regulation. A major part is accomplished by the combinatorial binding of transcription factors at regulatory sequences, such as enhancers. Even though binding of transcription factors is sequence-specific, estimating the sequence similarity of two functionally similar enhancers is very difficult. However, a similarity measure for regulatory sequences is crucial to detect and understand functional similarities between two enhancers and will facilitate large-scale analyses like clustering, prediction and classification of genome-wide datasets.

RESULTS

We present the standardized alignment-free sequence similarity measure N2, a flexible framework that is defined for word neighbourhoods. We explore the usefulness of adding reverse complement words as well as words including mismatches into the neighbourhood. On simulated enhancer sequences as well as functional enhancers in mouse development, N2 is shown to outperform previous alignment-free measures. N2 is flexible, faster than competing methods and less susceptible to single sequence noise and the occurrence of repetitive sequences. Experiments on the mouse enhancers reveal that enhancers active in different tissues can be separated by pairwise comparison using N2.

CONCLUSION

N2 represents an improvement over previous alignment-free similarity measures without compromising speed, which makes it a good candidate for large-scale sequence comparison of regulatory sequences.

AVAILABILITY

The software is part of the open-source C++ library SeqAn (www.seqan.de) and a compiled version can be downloaded at http://www.seqan.de/projects/alf.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Comparisons of the embryonic development of Drosophila, Nasonia, and Tribolium

Studying the embryogenesis of diverse insect species is crucial to understanding insect evolution. Here, we review current advances in understanding the development of two emerging model organisms: the wasp Nasonia vitripennis and the beetle Tribolium castaneum in comparison with the well-studied fruit fly Drosophila melanogaster. Although Nasonia represents the most basally branching order of holometabolous insects, it employs a derived long germband mode of embryogenesis, more like that of Drosophila, whereas Tribolium undergoes an intermediate germband mode of embryogenesis, which is more similar to the ancestral mechanism. Comparing the embryonic development and genetic regulation of early patterning events in these three insects has given invaluable insights into insect evolution. The similar mode of embryogenesis of Drosophila and Nasonia is reflected in their reliance on maternal morphogenetic gradients. However, they employ different genes as maternal factors, reflecting the evolutionary distance separating them. Tribolium, on the other hand, relies heavily on self-regulatory mechanisms other than maternal cues, reflecting its sequential nature of segmentation and the need for reiterated patterning.

Transcription factor binding site redundancy in embryonic enhancers of the Drosophila bithorax complex

The molecular control of gene expression in development is mediated through the activity of embryonic enhancer cis-regulatory modules. This activity is determined by the combination of repressor and activator transcription factors that bind at specific DNA sequences in the enhancer. A proposed mechanism to ensure a high fidelity of transcriptional output is functional redundancy between closely spaced binding sites within an enhancer. Here I show that at the bithorax complex in Drosophila there is selective redundancy for both repressor and activator factor binding sites in vivo. The absence of compensatory binding sites is responsible for two rare gain-of-function mutations in the complex.

Analysis of gap gene regulation in a 3D organism-scale model of the Drosophila melanogaster embryo

The axial bodyplan of Drosophila melanogaster is determined during a process called morphogenesis. Shortly after fertilization, maternal bicoid mRNA is translated into Bicoid (Bcd). This protein establishes a spatially graded morphogen distribution along the anterior-posterior (AP) axis of the embryo. Bcd initiates AP axis determination by triggering expression of gap genes that subsequently regulate each other's expression to form a precisely controlled spatial distribution of gene products. Reaction-diffusion models of gap gene expression on a 1D domain have previously been used to infer complex genetic regulatory network (GRN) interactions by optimizing model parameters with respect to 1D gap gene expression data. Here we construct a finite element reaction-diffusion model with a realistic 3D geometry fit to full 3D gap gene expression data. Though gap gene products exhibit dorsal-ventral asymmetries, we discover that previously inferred gap GRNs yield qualitatively correct AP distributions on the 3D domain only when DV-symmetric initial conditions are employed. Model patterning loses qualitative agreement with experimental data when we incorporate a realistic DV-asymmetric distribution of Bcd. Further, we find that geometry alone is insufficient to account for DV-asymmetries in the final gap gene distribution. Additional GRN optimization confirms that the 3D model remains sensitive to GRN parameter perturbations. Finally, we find that incorporation of 3D data in simulation and optimization does not constrain the search space or improve optimization results.

Molecular dissection of cis-regulatory modules at the Drosophila bithorax complex reveals critical transcription factor signature motifs

At the Drosophila melanogaster bithorax complex (BX-C) over 330kb of intergenic DNA is responsible for directing the transcription of just three homeotic (Hox) genes during embryonic development. A number of distinct enhancer cis-regulatory modules (CRMs) are responsible for controlling the specific expression patterns of the Hox genes in the BX-C. While it has proven possible to identify orthologs of known BX-C CRMs in different Drosophila species using overall sequence conservation, this approach has not proven sufficiently effective for identifying novel CRMs or defining the key functional sequences within enhancer CRMs. Here we demonstrate that the specific spatial clustering of transcription factor (TF) binding sites is important for BX-C enhancer activity. A bioinformatic search for combinations of putative TF binding sites in the BX-C suggests that simple clustering of binding sites is frequently not indicative of enhancer activity. However, through molecular dissection and evolutionary comparison across the Drosophila genus we discovered that specific TF binding site clustering patterns are an important feature of three known BX-C enhancers. Sub-regions of the defined IAB5 and IAB7b enhancers were both found to contain an evolutionarily conserved signature motif of clustered TF binding sites which is critical for the functional activity of the enhancers. Together, these results indicate that the spatial organization of specific activator and repressor binding sites within BX-C enhancers is of greater importance than overall sequence conservation and is indicative of enhancer functional activity.

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.

Huckebein is part of a combinatorial repression code in the anterior blastoderm

The hierarchy of the segmentation cascade responsible for establishing the Drosophila body plan is composed by gap, pair-rule and segment polarity genes. However, no pair-rule stripes are formed in the anterior regions of the embryo. This lack of stripe formation, as well as other evidence from the literature that is further investigated here, led us to the hypothesis that anterior gap genes might be involved in a combinatorial mechanism responsible for repressing the cis-regulatory modules (CRMs) of hairy (h), even-skipped (eve), runt (run), and fushi-tarazu (ftz) anterior-most stripes. In this study, we investigated huckebein (hkb), which has a gap expression domain at the anterior tip of the embryo. Using genetic methods we were able to detect deviations from the wild-type patterns of the anterior-most pair-rule stripes in different genetic backgrounds, which were consistent with Hkb-mediated repression. Moreover, we developed an image processing tool that, for the most part, confirmed our assumptions. Using an hkb misexpression system, we further detected specific repression on anterior stripes. Furthermore, bioinformatics analysis predicted an increased significance of binding site clusters in the CRMs of h 1, eve 1, run 1 and ftz 1when Hkb was incorporated in the analysis, indicating that Hkb plays a direct role in these CRMs. We further discuss that Hkb and Slp1, which is the other previously identified common repressor of anterior stripes, might participate in a combinatorial repression mechanism controlling stripe CRMs in the anterior parts of the embryo and define the borders of these anterior stripes.

Gene expression in time and space: additive vs hierarchical organization of cis-regulatory regions

In higher eukaryotes, individual genes are often intermingled with other genes and spread out across tens to hundreds of kilobases, even though only small portions of their sequence are devoted to protein coding. Yet, in this seemingly extended and tangled mess, the cell is able to precisely regulate gene expression in both time and space. Over the past few decades, numerous elements, like enhancers, silencers and insulators have been found that shed some light on how the precise control of gene expression is achieved. Through these discoveries, an additive model of gene expression was envisioned, where the addition of the patterning details imparted by regulatory elements would create the final pattern of gene expression. Although many genes can be described using this model, recent work in the Drosophila bithorax complex suggests that this model may be somewhat simplistic and, in fact, regulatory elements sometimes seem to communicate with each other to form a functional hierarchy that is far from additive.

Initiator elements function to determine the activity state of BX-C enhancers

A >300 kb cis-regulatory region is required for the proper expression of the three bithorax complex (BX-C) homeotic genes. Based on genetic and transgenic analysis, a model has been proposed in which the numerous BX-C cis-regulatory elements are spatially restricted through the activation or repression of parasegment-specific chromatin domains. Particular early embryonic enhancers, called initiators, have been proposed to control this complex process. Here, in order to better understand the process of domain activation, we have undertaken a systematic in situ dissection of the iab-6 cis-regulatory domain using a new method, called InSIRT. Using this method, we create and genetically characterize mutations affecting iab-6 function, including mutations specifically modifying the iab-6 initiator. Through our mutagenesis of the iab-6 initiator, we provide strong evidence that initiators function not to directly control homeotic gene expression but rather as domain control centers to determine the activity state of the enhancers and silencers within a cis-regulatory domain.

Understanding how genes become activated is one of the primary areas of research in modern biology. In order to decipher the DNA components required for this process, scientists have traditionally turned to transgenic reporter assays, where DNA elements are removed from their native environment and placed next to a simplified reporter gene to monitor transcriptional activation. Although this approach is powerful, it can result in artifacts stemming from the channelization of regulatory element activities into predetermined classes. In this manuscript, we investigate the biological role of elements from the Drosophila bithorax complex, called initiators. In transgenic assays, these elements have been categorized as enhancers. However, genetic analysis suggests that, in situ, these elements perform a far more complex function. Here, using a new method to repeatedly target a genetic locus for mutagenesis, we show that initiators function as control elements that coordinate the activity of nearby enhancers and silencers. Overall, our study highlights how gene expression can be controlled through a hierarchical arrangement of cis-regulatory elements.

Fates-shifted is an F-box protein that targets Bicoid for degradation and regulates developmental fate determination in Drosophila embryos

Bicoid (Bcd) is a morphogenetic protein that instructs patterning along the anterior-posterior (A-P) axis in Drosophila embryos. Despite extensive studies, what controls the formation of a normal concentration gradient of Bcd remains an unresolved and controversial question. In this report we show that Bcd protein degradation is mediated by the ubiquitin-proteasome pathway. We identify a novel F-box protein, encoded by fates-shifted (fsd), that plays an important role in Bcd protein degradation by targeting it for ubiquitination. Embryos from females lacking fsd have an altered Bcd gradient profile, resulting in a shift of the fatemap along the A-P axis. Our study represents a first experimental demonstration that, contrary to an alternative hypothesis, Bcd protein degradation is required for normal gradient formation and developmental fate determination.

Dynamic evolution of precise regulatory encodings creates the clustered site signature of enhancers

Concentration gradients of morphogenic proteins pattern the embryonic axes of Drosophila by activating different genes at different concentrations. The neurogenic ectoderm enhancers (NEEs) activate different genes at different threshold levels of the Dorsal (Dl) morphogen, which patterns the dorsal/ventral axis. NEEs share a unique arrangement of highly constrained DNA-binding sites for Dl, Twist (Twi), Snail (Sna) and Suppressor of Hairless (Su(H)), and encode the threshold variable in the precise length of DNA that separates one well-defined Dl element from a Twi element. However, NEEs also possess dense clusters of variant Dl sites. Here, we show that these increasingly variant sites are eclipsed relic elements, which were superseded by more recently evolved threshold encodings. Given the divergence in egg size during Drosophila lineage evolution, the observed characteristic clusters of divergent sites indicate a history of frequent selection for changes in threshold responses to the Dl morphogen gradient and confirm the NEE structure/function model.

giant is a bona fide gap gene in the intermediate germband insect, Oncopeltus fasciatus

Drosophila undergoes a form of development termed long germ segmentation, where all segments are specified nearly simultaneously so that by the blastoderm stage, the entire body plan has been determined. This mode of segmentation is evolutionarily derived. Most insects undergo short or intermediate germ segmentation, where only anterior segments are specified early, and posterior segments are sequentially specified during germband elongation. These embryological differences imply that anterior and posterior segments might rely upon different molecular mechanisms. In Drosophila, embryos mutant for giant show a gap in the anterior as well fusions of several abdominal segments. In Tribolium, a short germ beetle, giant is required for segmental identity, but not formation, in gnathal segments and also for segmentation of the entire abdomen. This raises the possibility that giant might not act as a gap gene in short and intermediate germ insects. Oncopeltus fasciatus is an intermediate germ insect that is an outgroup to the clade containing Drosophila and Tribolium. We cloned the Oncopeltus homolog of giant and determined its expression and function during segmentation. We find that Oncopeltus giant is a canonical gap gene in the maxillary and labial segments and also plays a gap-like role in the first four abdominal segments. Our results suggest that giant was a bona fide gap gene in the ancestor of these insects with this role being lost in the lineage leading towards Tribolium. This highlights the conservation of anterior patterning and evolutionary plasticity of the genetic regulation controlling posterior segmentation, even in short and intermediate germ insects.

Evolutionary mirages: selection on binding site composition creates the illusion of conserved grammars in Drosophila enhancers

The clustering of transcription factor binding sites in developmental enhancers and the apparent preferential conservation of clustered sites have been widely interpreted as proof that spatially constrained physical interactions between transcription factors are required for regulatory function. However, we show here that selection on the composition of enhancers alone, and not their internal structure, leads to the accumulation of clustered sites with evolutionary dynamics that suggest they are preferentially conserved. We simulated the evolution of idealized enhancers from Drosophila melanogaster constrained to contain only a minimum number of binding sites for one or more factors. Under this constraint, mutations that destroy an existing binding site are tolerated only if a compensating site has emerged elsewhere in the enhancer. Overlapping sites, such as those frequently observed for the activator Bicoid and repressor Krüppel, had significantly longer evolutionary half-lives than isolated sites for the same factors. This leads to a substantially higher density of overlapping sites than expected by chance and the appearance that such sites are preferentially conserved. Because D. melanogaster (like many other species) has a bias for deletions over insertions, sites tended to become closer together over time, leading to an overall clustering of sites in the absence of any selection for clustered sites. Since this effect is strongest for the oldest sites, clustered sites also incorrectly appear to be preferentially conserved. Following speciation, sites tend to be closer together in all descendent species than in their common ancestors, violating the common assumption that shared features of species' genomes reflect their ancestral state. Finally, we show that selection on binding site composition alone recapitulates the observed number of overlapping and closely neighboring sites in real D. melanogaster enhancers. Thus, this study calls into question the common practice of inferring "cis-regulatory grammars" from the organization and evolutionary dynamics of developmental enhancers.

Dissecting the regulatory switches of development: lessons from enhancer evolution in Drosophila

Cis-regulatory modules are non-protein-coding regions of DNA essential for the control of gene expression. One class of regulatory modules is embryonic enhancers, which drive gene expression during development as a result of transcription factor protein binding at the enhancer sequences. Recent comparative studies have begun to investigate the evolution of the sequence architecture within enhancers. These analyses are illuminating the way that developmental biologists think about enhancers by revealing their molecular mechanism of function.

Statistical tests for natural selection on regulatory regions based on the strength of transcription factor binding sites

Background

Although cis-regulatory changes play an important role in evolution, it remains difficult to establish the contribution of natural selection to regulatory differences between species. For protein coding regions, powerful tests of natural selection have been developed based on comparisons of synonymous and non-synonymous substitutions, and analogous tests for regulatory regions would be of great utility.

Results

Here, tests for natural selection on regulatory regions are proposed based on nucleotide substitutions that occur in characterized transcription factor binding sites (an important type functional element within regulatory regions). In the absence of selection, these substitutions will tend to reduce the strength of existing binding sites. On the other hand, purifying selection will act to preserve the binding sites in regulatory regions, while positive selection can act to create or destroy binding sites, as well as change their strength. Using standard models of binding site strength and molecular evolution in the absence of selection, this intuition can be used to develop statistical tests for natural selection. Application of these tests to two well-characterized regulatory regions in Drosophila provides evidence for purifying selection.

Conclusion

This demonstrates that it is possible to develop tests for selection on regulatory regions based on the specific functional constrains on these sequences.

Organization of developmental enhancers in the Drosophila embryo

Most cell-specific enhancers are thought to lack an inherent organization, with critical binding sites distributed in a more or less random fashion. However, there are examples of fixed arrangements of binding sites, such as helical phasing, that promote the formation of higher-order protein complexes on the enhancer DNA template. Here, we investigate the regulatory ‘grammar’ of nearly 100 characterized enhancers for developmental control genes active in the early Drosophila embryo. The conservation of grammar is examined in seven divergent Drosophila genomes. Linked binding sites are observed for particular combinations of binding motifs, including Bicoid–Bicoid, Hunchback–Hunchback, Bicoid–Dorsal, Bicoid–Caudal and Dorsal–Twist. Direct evidence is presented for the importance of Bicoid–Dorsal linkage in the integration of the anterior–posterior and dorsal–ventral patterning systems. Hunchback–Hunchback interactions help explain unresolved aspects of segmentation, including the differential regulation of the eve stripe 3 + 7 and stripe 4 + 6 enhancers. We also present evidence that there is an under-representation of nucleosome positioning sequences in many enhancers, raising the possibility for a subtle higher-order structure extending across certain enhancers. We conclude that grammar of gene control regions is pervasively used in the patterning of the Drosophila embryo.

Transcription of the C. elegans let-7 microRNA is temporally regulated by one of its targets, hbl-1

The let-7 family of microRNAs (miRNAs) are important regulators of developmental timing and cell differentiation and are often misexpressed in human cancer. In C. elegans, let-7 controls cell fate transitions from larval stage 4 (L4) to adulthood by post-transcriptionally down-regulating lineage-abnormal 41 (lin-41) and hunchback-like 1 (hbl-1). Primary let-7 (pri-let-7) transcripts are up-regulated in the L3, yet little is known about what controls this transcriptional up-regulation. We sought factors that either turn on let-7 transcription or keep it repressed until the correct time. Here we report that one of let-7's targets, the transcription factor Hunchback-like 1 (HBL-1), is responsible for inhibiting the transcription of let-7 in specific tissues until the L3. hbl-1 is a known developmental timing regulator and inhibits adult development in larval stages. Therefore, one important function of HBL-1 in maintaining larval stage fates is inhibition of let-7. Indeed, our results reveal let-7 as the first known target of the HBL-1 transcription factor in C. elegans and suggest a negative feedback loop mechanism for let-7 and HBL-1 regulation.

Cytoplasmic Prep1 interacts with 4EHP inhibiting Hoxb4 translation

Background

Homeobox genes are essential for embryonic patterning and cell fate determination. They are regulated mostly at the transcriptional level. In particular, Prep1 regulates Hox transcription in association with Pbx proteins. Despite its nuclear role as a transcription factor, Prep1 is located in the cytosol of mouse oocytes from primary to antral follicles. The homeodomain factor Bicoid (Bcd) has been shown to interact with 4EHP (eukaryotic translation initiation factor 4E homolog protein) to repress translation of Caudal mRNA and to drive Drosophila embryo development. Interestingly, Prep1 contains a putative binding motif for 4EHP, which may reflect a novel unknown function.

Methodology/Principal Findings

In this paper we show by confocal microscopy and deconvolution analysis that Prep1 and 4EHP co-localize in the cytosol of growing mouse oocytes, demonstrating their interaction by co-immunoprecipitation and pull-down experiments. A functional 4EHP-binding motif present in Prep1 has been also identified by mutagenesis analysis. Moreover, Prep1 inhibits (>95%) the in vitro translation of a luciferase reporter mRNA fused to the Hoxb4 3′UTR, in the presence of 4EHP. RNA electrophoretic mobility shift assay was used to demonstrate that Prep1 binds the Hoxb4 3′UTR. Furthermore, conventional histology and immunohistochemistry has shown a dramatic oocyte growth failure in hypomorphic mouse Prep1^i/i females, accompanied by an increased production of Hoxb4. Finally, Hoxb4 overexpression in mouse zygotes showed a slow in vitro development effect.

Conclusions

Prep1 has a novel cytoplasmic, 4EHP-dependent, function in the regulation of translation. Mechanistically, the Prep1-4EHP interaction might bridge the 3′UTR of Hoxb4 mRNA to the 5′ cap structure. This is the first demonstration that a mammalian homeodomain transcription factor regulates translation, and that this function can be possibly essential for the development of female germ cells and involved in mammalian zygote development.

Evolution acts on enhancer organization to fine-tune gradient threshold readouts

The elucidation of principles governing evolution of gene regulatory sequence is critical to the study of metazoan diversification. We are therefore exploring the structure and organizational constraints of regulatory sequences by studying functionally equivalent cis-regulatory modules (CRMs) that have been evolving in parallel across several loci. Such an independent dataset allows a multi-locus study that is not hampered by nonfunctional or constrained homology. The neurogenic ectoderm enhancers (NEEs) of Drosophila melanogaster are one such class of coordinately regulated CRMs. The NEEs share a common organization of binding sites and as a set would be useful to study the relationship between CRM organization and CRM activity across evolving lineages. We used the D. melanogaster transgenic system to screen for functional adaptations in the NEEs from divergent drosophilid species. We show that the individual NEE modules across a genome in any one lineage have independently evolved adaptations to compensate for lineage-specific developmental and/or genomic changes. Specifically, we show that both the site composition and the site organization of NEEs have been finely tuned by distinct, lineage-specific selection pressures in each of the three divergent species that we have examined: D. melanogaster, D. pseudoobscura, and D. virilis. Furthermore, by precisely altering the organization of NEEs with different morphogen gradient threshold readouts, we show that CRM organizational evolution is sufficient for explaining changes in enhancer activity. Thus, evolution can act on CRM organization to fine-tune morphogen gradient threshold readouts over a wide dynamic range. Our study demonstrates that equivalence classes of CRMs are powerful tools for detecting lineage-specific adaptations by gene regulatory sequences.

The regulatory control of genes allows an organism to generate a diversity of cell types throughout its body. Gene regulation involves specialized DNA sequences called transcriptional enhancers that increase the expression of genes in specific places and times. Enhancers contain clusters of specific DNA sequences that are uniquely recognized by DNA binding proteins, whose activities are also regulated in space and time. The critical role that DNA enhancers play in generating the diversity of cell types within a single organism suggests that changes in these DNA sequences may also underlie the diversity of organismal forms produced by evolution. However, few examples linking specific changes in enhancer sequences to functional adaptations have been documented. We studied a group of neuro-embryonic enhancers that turn on a certain group of genes in different fruit fly species that have been diverging from each other for ∼50 million years. Each species has experienced unique changes in its protein-coding sequences, gene regulatory sequences, egg morphology, and developmental timing. We found that the organizational spacing between the protein binding sites in these enhancers has evolved in a manner that is consistent with functional adaptations compensating for the dynamic and idiosyncratic evolutionary history of each species.

Characterizing neuroectodermal enhancers located throughout the genomes of threeDrosophila species reveals that each enhancer evolved functional adaptations that compensate for the evolutionary changes specific to each lineage.

One thousand and one ways of making functionally similar transcriptional enhancers

Expression of most genes is regulated by the interaction of multiple transcription factors with cis-regulatory sequences. Many studies have focused on how changes in promoters and enhancers alter gene expression and phenotype. Recently, Hare et al., using elegant wet and computational approaches uncovered a series of enhancers driving the expression of the even-skipped gene in scavenger flies (Sepsidae). Despite the strong sequence divergence between the enhancers in sepsids and drosophilids, they lead to remarkably similar patterns of gene expression in transgenic Drosophila embryos. This can be explained by the existence of intra-enhancer compensatory mutations and the presence of overlapping/near binding sites for activators and repressors.