The even-skipped locus is contained in a 16-kb chromatin domain

The Density of Regulatory Information Is a Major Determinant of Evolutionary Constraint on Noncoding DNA in Drosophila

Evolutionary analyses have estimated that ∼60% of nucleotides in intergenic regions of the Drosophila melanogaster genome are functionally relevant, suggesting that regulatory information may be encoded more densely in intergenic regions than has been revealed by most functional dissections of regulatory DNA. Here, we approached this issue through a functional dissection of the regulatory region of the gene shavenbaby (svb). Most of the ∼90 kb of this large regulatory region is highly conserved in the genus Drosophila, though characterized enhancers occupy a small fraction of this region. By analyzing the regulation of svb in different contexts of Drosophila development, we found that the regulatory information that drives svb expression in the abdominal pupal epidermis is organized in a different way than the elements that drive svb expression in the embryonic epidermis. While in the embryonic epidermis svb is activated by compact enhancers separated by large inactive DNA regions, svb expression in the pupal epidermis is driven by regulatory information distributed over broader regions of svb cis-regulatory DNA. In the same vein, we observed that other developmental genes also display a dense distribution of putative regulatory elements in their regulatory regions. Furthermore, we found that a large percentage of conserved noncoding DNA of the Drosophila genome is contained within regions of open chromatin. These results suggest that part of the evolutionary constraint on noncoding DNA of Drosophila is explained by the density of regulatory information, which may be greater than previously appreciated.

Mechanisms of Interaction between Enhancers and Promoters in Three Drosophila Model Systems

In higher eukaryotes, the regulation of developmental gene expression is determined by enhancers, which are often located at a large distance from the promoters they regulate. Therefore, the architecture of chromosomes and the mechanisms that determine the functional interaction between enhancers and promoters are of decisive importance in the development of organisms. Mammals and the model animal Drosophila have homologous key architectural proteins and similar mechanisms in the organization of chromosome architecture. This review describes the current progress in understanding the mechanisms of the formation and regulation of long-range interactions between enhancers and promoters at three well-studied key regulatory loci in Drosophila.

A timer gene network is spatially regulated by the terminal system in the Drosophila embryo

In insect embryos, anteroposterior patterning is coordinated by the sequential expression of the 'timer' genes caudal, Dichaete, and odd-paired, whose expression dynamics correlate with the mode of segmentation. In Drosophila, the timer genes are expressed broadly across much of the blastoderm, which segments simultaneously, but their expression is delayed in a small 'tail' region, just anterior to the hindgut, which segments during germband extension. Specification of the tail and the hindgut depends on the terminal gap gene tailless, but beyond this the regulation of the timer genes is poorly understood. We used a combination of multiplexed imaging, mutant analysis, and gene network modelling to resolve the regulation of the timer genes, identifying 11 new regulatory interactions and clarifying the mechanism of posterior terminal patterning. We propose that a dynamic Tailless expression gradient modulates the intrinsic dynamics of a timer gene cross-regulatory module, delineating the tail region and delaying its developmental maturation.

An insulator blocks access to enhancers by an illegitimate promoter, preventing repression by transcriptional interference

Several distinct activities and functions have been described for chromatin insulators, which separate genes along chromosomes into functional units. Here, we describe a novel mechanism of functional separation whereby an insulator prevents gene repression. When the homie insulator is deleted from the end of a Drosophila even skipped (eve) locus, a flanking P-element promoter is activated in a partial eve pattern, causing expression driven by enhancers in the 3’ region to be repressed. The mechanism involves transcriptional read-through from the flanking promoter. This conclusion is based on the following. Read-through driven by a heterologous enhancer is sufficient to repress, even when homie is in place. Furthermore, when the flanking promoter is turned around, repression is minimal. Transcriptional read-through that does not produce anti-sense RNA can still repress expression, ruling out RNAi as the mechanism in this case. Thus, transcriptional interference, caused by enhancer capture and read-through when the insulator is removed, represses eve promoter-driven expression. We also show that enhancer-promoter specificity and processivity of transcription can have decisive effects on the consequences of insulator removal. First, a core heat shock 70 promoter that is not activated well by eve enhancers did not cause read-through sufficient to repress the eve promoter. Second, these transcripts are less processive than those initiated at the P-promoter, measured by how far they extend through the eve locus, and so are less disruptive. These results highlight the importance of considering transcriptional read-through when assessing the effects of insulators on gene expression.

Several distinct activities and functions have been described for chromatin insulators, which are regulatory DNA elements that separate genes along chromosomes into functional units. Here, we describe how insulators can prevent repression of one gene by preventing inappropriate transcription of another gene, without blocking read-through of transcription per se. When the insulator homie is deleted from the end of a transgenic eve locus, a flanking transposable element promoter is activated by eve enhancers, causing repression of the eve promoter. The mechanism involves transcriptional read-through from the flanking promoter, which disrupts normal eve enhancer-promoter activities. When the flanking promoter is turned around, repression of eve is minimal. Thus, transcriptional interference, caused by enhancer capture and read-through when the insulator is removed, represses the eve promoter. These results show a novel role for transcriptional read-through in the effects of insulators on gene expression.

Rationally-engineered reproductive barriers using CRISPR & CRISPRa: an evaluation of the synthetic species concept in Drosophila melanogaster

The ability to erect rationally-engineered reproductive barriers in animal or plant species promises to enable a number of biotechnological applications such as the creation of genetic firewalls, the containment of gene drives or novel population replacement and suppression strategies for genetic control. However, to date no experimental data exist that explores this concept in a multicellular organism. Here we examine the requirements for building artificial reproductive barriers in the metazoan model Drosophila melanogaster by combining CRISPR-based genome editing and transcriptional transactivation (CRISPRa) of the same loci. We directed 13 single guide RNAs (sgRNAs) to the promoters of 7 evolutionary conserved genes and used 11 drivers to conduct a misactivation screen. We identify dominant-lethal activators of the eve locus and find that they disrupt development by strongly activating eve outside its native spatio-temporal context. We employ the same set of sgRNAs to isolate, by genome editing, protective INDELs that render these loci resistant to transactivation without interfering with target gene function. When these sets of genetic components are combined we find that complete synthetic lethality, a prerequisite for most applications, is achievable using this approach. However, our results suggest a steep trade-off between the level and scope of dCas9 expression, the degree of genetic isolation achievable and the resulting impact on fly fitness. The genetic engineering strategy we present here allows the creation of single or multiple reproductive barriers and could be applied to other multicellular organisms such as disease vectors or transgenic organisms of economic importance.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

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

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

Odd-paired controls frequency doubling in Drosophila segmentation by altering the pair-rule gene regulatory network

The Drosophila embryo transiently exhibits a double-segment periodicity, defined by the expression of seven 'pair-rule' genes, each in a pattern of seven stripes. At gastrulation, interactions between the pair-rule genes lead to frequency doubling and the patterning of 14 parasegment boundaries. In contrast to earlier stages of Drosophila anteroposterior patterning, this transition is not well understood. By carefully analysing the spatiotemporal dynamics of pair-rule gene expression, we demonstrate that frequency-doubling is precipitated by multiple coordinated changes to the network of regulatory interactions between the pair-rule genes. We identify the broadly expressed but temporally patterned transcription factor, Odd-paired (Opa/Zic), as the cause of these changes, and show that the patterning of the even-numbered parasegment boundaries relies on Opa-dependent regulatory interactions. Our findings indicate that the pair-rule gene regulatory network has a temporally modulated topology, permitting the pair-rule genes to play stage-specific patterning roles.

Determinants of Chromosome Architecture: Insulator Pairing in cis and in trans

The chromosomes of multicellular animals are organized into a series of topologically independent looped domains. This domain organization is critical for the proper utilization and propagation of the genetic information encoded by the chromosome. A special set of architectural elements, called boundaries or insulators, are responsible both for subdividing the chromatin into discrete domains and for determining the topological organization of these domains. Central to the architectural functions of insulators are homologous and heterologous insulator:insulator pairing interactions. The former (pairing between copies of the same insulator) dictates the process of homolog alignment and pairing in trans, while the latter (pairing between different insulators) defines the topology of looped domains in cis. To elucidate the principles governing these architectural functions, we use two insulators, Homie and Nhomie, that flank the Drosophila even skipped locus. We show that homologous insulator interactions in trans, between Homie on one homolog and Homie on the other, or between Nhomie on one homolog and Nhomie on the other, mediate transvection. Critically, these homologous insulator:insulator interactions are orientation-dependent. Consistent with a role in the alignment and pairing of homologs, self-pairing in trans is head-to-head. Head-to-head self-interactions in cis have been reported for other fly insulators, suggesting that this is a general principle of self-pairing. Homie and Nhomie not only pair with themselves, but with each other. Heterologous Homie-Nhomie interactions occur in cis, and we show that they serve to delimit a looped chromosomal domain that contains the even skipped transcription unit and its associated enhancers. The topology of this loop is defined by the heterologous pairing properties of Homie and Nhomie. Instead of being head-to-head, which would generate a circular loop, Homie-Nhomie pairing is head-to-tail. Head-to-tail pairing in cis generates a stem-loop, a configuration much like that observed in classical lampbrush chromosomes. These pairing principles provide a mechanistic underpinning for the observed topologies within and between chromosomes.

The chromosomes of multicellular animals are organized into a series of topologically independent looped domains. This domain organization is critical for the proper utilization and propagation of the genetic information encoded by the chromosome. A special set of architectural elements, called boundaries or insulators, are responsible for both subdividing the chromatin fiber into discrete domains, and determining the topological organization of these domains. Central to the architectural functions of insulators are heterologous and homologous insulator:insulator pairing interactions. In Drosophila, the former defines the topology of individual looped domains in cis, while the latter dictates the process of homolog alignment and pairing in trans. Here we use two insulators from the even skipped locus to elucidate the principles governing these two architectural functions. These principles align with several longstanding observations, and resolve a number of conundrums regarding chromosome topology and function.

The Drosophila eve insulator Homie promotes eve expression and protects the adjacent gene from repression by polycomb spreading

Insulators can block the action of enhancers on promoters and the spreading of repressive chromatin, as well as facilitating specific enhancer-promoter interactions. However, recent studies have called into question whether the activities ascribed to insulators in model transgene assays actually reflect their functions in the genome. The Drosophila even skipped (eve) gene is a Polycomb (Pc) domain with a Pc-group response element (PRE) at one end, flanked by an insulator, an arrangement also seen in other genes. Here, we show that this insulator has three major functions. It blocks the spreading of the eve Pc domain, preventing repression of the adjacent gene, TER94. It prevents activation of TER94 by eve regulatory DNA. It also facilitates normal eve expression. When Homie is deleted in the context of a large transgene that mimics both eve and TER94 regulation, TER94 is repressed. This repression depends on the eve PRE. Ubiquitous TER94 expression is “replaced” by expression in an eve pattern when Homie is deleted, and this effect is reversed when the PRE is also removed. Repression of TER94 is attributable to spreading of the eve Pc domain into the TER94 locus, accompanied by an increase in histone H3 trimethylation at lysine 27. Other PREs can functionally replace the eve PRE, and other insulators can block PRE-dependent repression in this context. The full activity of the eve promoter is also dependent on Homie, and other insulators can promote normal eve enhancer-promoter communication. Our data suggest that this is not due to preventing promoter competition, but is likely the result of the insulator organizing a chromosomal conformation favorable to normal enhancer-promoter interactions. Thus, insulator activities in a native context include enhancer blocking and enhancer-promoter facilitation, as well as preventing the spread of repressive chromatin.

Insulators are specialized DNA elements that can separate the genome into functional units. Most of the current thinking about these elements comes from studies done with model transgenes. Studies of insulators within the specialized Hox gene complexes have suggested that model transgenes can reflect the normal functions of these elements in their native context. However, recent genome-wide studies have called this into question. This work analyzes the native function of an insulator that resides between the Drosophila genes eve and TER94, which are expressed in very different patterns. Also, the eve gene is a Polycomb (Pc) domain, a specialized type of chromatin that is found in many places throughout the genome. We show that this insulator has three major functions. It blocks the spreading of the eve Pc domain, preventing repression of TER94. It prevents activation of TER94 by eve regulatory DNA. It also facilitates normal eve expression. Each of these activities are consistent with those seen with model transgenes, and other known insulators can provide these functions in this context. This work provides a novel and convincing example of the normal role of insulators in regulating the eukaryotic genome, as well as providing insights into their mechanisms of action.

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.

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

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

The regulatory region controlling the nitrate-responsive expression of a nitrate reductase gene, NIA1, in Arabidopsis

Nitrate reductase (NR) is the enzyme that catalyzes the first step of nitrate assimilation. It is well known that the expression of NR genes is rapidly induced in various plants by nitrate. Previously, the activity of a tobacco NR gene promoter was reported to be high in tobacco plants grown on medium containing ammonium as the sole nitrogen source, but low in tobacco plants grown on nitrate-containing medium. This cast some doubt on the role of the NR gene promoter in the nitrate-inducible expression of this gene. Furthermore, in previous studies, transformation with genomic fragments containing NR loci restored the reduced NR activity in NR mutants to a limited extent, suggesting a complex regulation of NR gene expression. Here, we show that although the 1.9 kb promoter of an NR gene in Arabidopsis, NIA1, is not activated by nitrate, the expression of a GUS (β-glucuronidase) reporter gene inserted between the 5'- and 3'-flanking sequences of the NIA1 coding region is strongly induced by nitrate. When the 3'-flanking sequence was fused downstream of the GUS gene under the control of the 35S minimal promoter, its expression was also strongly induced by nitrate. Furthermore, dissection analysis of the 3'-flanking region revealed that the sequence downstream of the transcriptional terminator rather than the 3'-untranslated region plays a role in nitrate-inducible expression, indicating a requirement for the 3'-flanking sequence for the nitrate-inducible transcription of NIA1. We also show that the 2.7 kb promoter sequence of NIA2, another NR gene of Arabidopsis, cannot direct nitrate-inducible expression.

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

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

A chromatin insulator mediates transgene homing and very long-range enhancer-promoter communication

Insulator sequences help to organize the genome into discrete functional regions by preventing inappropriate cross-regulation. This is thought to be mediated in part through associations with other insulators located elsewhere in the genome. Enhancers that normally drive Drosophila even skipped (eve) expression are located closer to the TER94 transcription start site than to that of eve. We discovered that the region between these genes has enhancer-blocking activity, and that this insulator region also mediates homing of P-element transgenes to the eve-TER94 genomic neighborhood. Localization of these activities to within 0.6 kb failed to separate them. Importantly, homed transgenic promoters respond to endogenous eve enhancers from great distances, and this long-range communication depends on the homing/insulator region, which we call Homie. We also find that the eve promoter contributes to long-distance communication. However, even the basal hsp70 promoter can communicate with eve enhancers across distances of several megabases, when the communication is mediated by Homie. These studies show that, while Homie blocks enhancer-promoter communication at short range, it facilitates long-range communication between distant genomic regions, possibly by organizing a large chromosomal loop between endogenous and transgenic Homies.

Effects of random mutations in the human immunodeficiency virus type 1 transcriptional promoter on viral fitness in different host cell environments

A mutation's effect on fitness or phenotype may in part depend on the interaction of the mutation with the environment. The resulting phenotype or fitness is important, since it determines the adaptive potential of a species. To date, most studies have focused on alterations to protein-coding regions of the genome and their consequential fitness effects. Non-protein-coding regulatory regions have been largely neglected, although they make up a large and important part of an organism's genome. Here, we use human immunodeficiency virus type 1 as a model system to investigate fitness effects of random mutations in noncoding DNA-binding sites of the transcriptional promoter. We determined 242 fitness values for 35 viral promoter mutants with one, two, or three mutations across seven distinct cellular environments and identified that (i) all mutants have an effect in at least one cellular environment; (ii) fitness effects are highly dependent on the cellular environment; (iii) disadvantageous and advantageous mutations occur at high and similar frequencies; and (iv) epistatic effects of multiple mutations are rare. Our results underline the evolutionary potential of regulatory regions and indicate that DNA-binding sites evolve under strong selection, while at the same time, they are very plastic to environmental change.

Embryonic even skipped-dependent muscle and heart cell fates are required for normal adult activity, heart function, and lifespan

The Drosophila pair-rule gene even skipped (eve) is required for embryonic segmentation and later in specific cell lineages in both the nervous system and the mesoderm. We previously generated eve mesoderm-specific mutants by combining an eve null mutant with a rescuing transgene that includes the entire locus, but with the mesodermal enhancer removed. This allowed us to analyze in detail the defects that result from a precisely targeted elimination of mesodermal eve expression in the context of an otherwise normal embryo. Absence of mesodermal eve causes a highly selective loss of the entire eve-expressing lineage in this germ layer, including those progeny that do not continue to express eve, suggesting that mesodermal eve precursor specification is not implemented. Despite the resulting absence of a subset of muscles and pericardial cells, mesoderm-specific eve mutants survive to fertile adulthood, providing an opportunity to examine the effects of these developmental abnormalities on adult fitness and heart function. We find that in these mutants, flying ability, myocardial performance under normal and stressed conditions, and lifespan are severely reduced. These data imply a nonautonomous role of the affected pericardial cells and body wall muscles in developing and/or maintaining cardiac performance and possibly other functions contributing to normal lifespan. Given the similarities of molecular-genetic control between Drosophila and vertebrates, these findings suggest that peri/epicardial influences may well be important for proper myocardial function.

cis-regulatory logic of short-range transcriptional repression in Drosophila melanogaster

Bioinformatics analysis of transcriptional control is guided by knowledge of the characteristics of cis-regulatory regions or enhancers. Features such as clustering of binding sites and co-occurrence of binding sites have aided enhancer identification, but quantitative predictions of enhancer function are not yet generally feasible. To facilitate the analysis of regulatory sequences in Drosophila melanogaster, we identified quantitative parameters that affect the activity of short-range transcriptional repressors, proteins that play key roles in development. In addition to the previously noted distance dependence, repression is strongly influenced by the stoichiometry, affinity, spacing, and arrangement of activator binding sites. Repression is insensitive to the type of activation domain, suggesting that short-range repression may primarily affect activators at the level of DNA binding. The activity of several short-range, but not long-range, repressors is circumscribed by the same quantitative parameters. This cis-regulatory "grammar" may aid the identification of enhancers regulated by short-range repressors and facilitate bioinformatic prediction of the functional output of transcriptional regulatory sequences.

Functional evolution of a cis-regulatory module

Lack of knowledge about how regulatory regions evolve in relation to their structure-function may limit the utility of comparative sequence analysis in deciphering cis-regulatory sequences. To address this we applied reverse genetics to carry out a functional genetic complementation analysis of a eukaryotic cis-regulatory module-the even-skipped stripe 2 enhancer-from four Drosophila species. The evolution of this enhancer is non-clock-like, with important functional differences between closely related species and functional convergence between distantly related species. Functional divergence is attributable to differences in activation levels rather than spatiotemporal control of gene expression. Our findings have implications for understanding enhancer structure-function, mechanisms of speciation and computational identification of regulatory modules.

Identifying spatially similar gene expression patterns in early stage fruit fly embryo images: binary feature versus invariant moment digital representations

Background

Modern developmental biology relies heavily on the analysis of embryonic gene expression patterns. Investigators manually inspect hundreds or thousands of expression patterns to identify those that are spatially similar and to ultimately infer potential gene interactions. However, the rapid accumulation of gene expression pattern data over the last two decades, facilitated by high-throughput techniques, has produced a need for the development of efficient approaches for direct comparison of images, rather than their textual descriptions, to identify spatially similar expression patterns.

Results

The effectiveness of the Binary Feature Vector (BFV) and Invariant Moment Vector (IMV) based digital representations of the gene expression patterns in finding biologically meaningful patterns was compared for a small (226 images) and a large (1819 images) dataset. For each dataset, an ordered list of images, with respect to a query image, was generated to identify overlapping and similar gene expression patterns, in a manner comparable to what a developmental biologist might do. The results showed that the BFV representation consistently outperforms the IMV representation in finding biologically meaningful matches when spatial overlap of the gene expression pattern and the genes involved are considered. Furthermore, we explored the value of conducting image-content based searches in a dataset where individual expression components (or domains) of multi-domain expression patterns were also included separately. We found that this technique improves performance of both IMV and BFV based searches.

Conclusions

We conclude that the BFV representation consistently produces a more extensive and better list of biologically useful patterns than the IMV representation. The high quality of results obtained scales well as the search database becomes larger, which encourages efforts to build automated image query and retrieval systems for spatial gene expression patterns.

Computational identification of developmental enhancers: conservation and function of transcription factor binding-site clusters in Drosophila melanogaster and Drosophila pseudoobscura

27 predicted gene-regulatory regions in the Drosophila melanogaster genome were analyzed in vivo, confirming 15 active enhancer regions. A comparison with Drosophila pseudoobscura sequences revealed that conservation of binding-site clusters accurately discriminates functional regions from non-functional ones.

Background

The identification of sequences that control transcription in metazoans is a major goal of genome analysis. In a previous study, we demonstrated that searching for clusters of predicted transcription factor binding sites could discover active regulatory sequences, and identified 37 regions of the Drosophila melanogaster genome with high densities of predicted binding sites for five transcription factors involved in anterior-posterior embryonic patterning. Nine of these clusters overlapped known enhancers. Here, we report the results of in vivo functional analysis of 27 remaining clusters.

Results

We generated transgenic flies carrying each cluster attached to a basal promoter and reporter gene, and assayed embryos for reporter gene expression. Six clusters are enhancers of adjacent genes: giant, fushi tarazu, odd-skipped, nubbin, squeeze and pdm2; three drive expression in patterns unrelated to those of neighboring genes; the remaining 18 do not appear to have enhancer activity. We used the Drosophila pseudoobscura genome to compare patterns of evolution in and around the 15 positive and 18 false-positive predictions. Although conservation of primary sequence cannot distinguish true from false positives, conservation of binding-site clustering accurately discriminates functional binding-site clusters from those with no function. We incorporated conservation of binding-site clustering into a new genome-wide enhancer screen, and predict several hundred new regulatory sequences, including 85 adjacent to genes with embryonic patterns.

Conclusions

Measuring conservation of sequence features closely linked to function - such as binding-site clustering - makes better use of comparative sequence data than commonly used methods that examine only sequence identity.

The regulatory content of intergenic DNA shapes genome architecture

The relationship between regulatory complexity and gene spacing was examined in Caenorhabditis elegans and Drosophila melanogaster. Intergenic distance, and hence genome architecture, is shaped by regulatory information contained in noncoding DNA.

Background

Factors affecting the organization and spacing of functionally unrelated genes in metazoan genomes are not well understood. Because of the vast size of a typical metazoan genome compared to known regulatory and protein-coding regions, functional DNA is generally considered to have a negligible impact on gene spacing and genome organization. In particular, it has been impossible to estimate the global impact, if any, of regulatory elements on genome architecture.

Results

To investigate this, we examined the relationship between regulatory complexity and gene spacing in Caenorhabditis elegans and Drosophila melanogaster. We found that gene density directly reflects local regulatory complexity, such that the amount of noncoding DNA between a gene and its nearest neighbors correlates positively with that gene's regulatory complexity. Genes with complex functions are flanked by significantly more noncoding DNA than genes with simple or housekeeping functions. Genes of low regulatory complexity are associated with approximately the same amount of noncoding DNA in D. melanogaster and C. elegans, while loci of high regulatory complexity are significantly larger in the more complex animal. Complex genes in C. elegans have larger 5' than 3' noncoding intervals, whereas those in D. melanogaster have roughly equivalent 5' and 3' noncoding intervals.

Conclusions

Intergenic distance, and hence genome architecture, is highly nonrandom. Rather, it is shaped by regulatory information contained in noncoding DNA. Our findings suggest that in compact genomes, the species-specific loss of nonfunctional DNA reveals a landscape of regulatory information by leaving a profile of functional DNA in its wake.

Even-skipped, acting as a repressor, regulates axonal projections in Drosophila

Nervous system-specific eve mutants were created by removing regulatory elements from a 16 kb transgene capable of complete rescue of normal eve function. When transgenes lacking the regulatory element for either RP2+a/pCC, EL or U/CQ neurons were placed in an eve-null background, eve expression was completely eliminated in the corresponding neurons, without affecting other aspects of eve expression. Many of these transgenic flies were able to survive to fertile adulthood. In the RP2+a/pCC mutant flies: (1) both RP2 and aCC showed abnormal axonal projection patterns, failing to innervate their normal target muscles; (2) the cell bodies of these neurons were positioned abnormally; and (3) in contrast to the wild type, pCC axons often crossed the midline. The Eve HD alone was able to provide a weak, partial rescue of the mutant phenotype, while both the Groucho-dependent and -independent repressor domains contributed equally to full rescue of each aspect of the mutant phenotype. Complete rescue was also obtained with a chimeric protein containing the Eve HD and the Engrailed repressor domain. Consistent with the apparent sufficiency of repressor function, a fusion protein between the Gal4 DNA-binding domain and Eve repressor domains was capable of actively repressing UAS target genes in these neurons. A key target of the repressor function of Eve was Drosophila Hb9, the derepression of which correlated with the mutant phenotype in individual eve-mutant neurons. Finally, homologues of Eve from diverse species were able to rescue the eve mutant phenotype, indicating conservation of both targeting and repression functions in the nervous system.

Homotypic regulatory clusters in Drosophila

Cis-regulatory modules (CRMs) are transcription regulatory DNA segments (approximately 1 Kb range) that control the expression of developmental genes in higher eukaryotes. We analyzed clustering of known binding motifs for transcription factors (TFs) in over 60 known CRMs from 20 Drosophila developmental genes, and we present evidence that each type of recognition motif forms significant clusters within the regulatory regions regulated by the corresponding TF. We demonstrate how a search with a single binding motif can be applied to explore gene regulatory networks and to discover coregulated genes in the genome. We also discuss the potential of the clustering method in interpreting the differential response of genes to various levels of transcriptional regulators.

Mutants highlight the modular control of butterfly eyespot patterns

The eyespots on butterfly wings are thought to be serially homologous pattern elements. Yet eyespots differ greatly in number, shape, color, and size, within and among species. To what extent do these serially homologues have separate developmental identities, upon which selection acts to create diversity? We examined x-ray-induced mutations for the eyespots of the nymphalid butterfly Bicyclus anynana that highlight the modular control of these serially homologous wing pattern elements. These mutations reduce or eliminate individual eyespots, or groups of eyespots, with no further effect on the wing color pattern. The collection of mutants highlights a greater potential developmental repertoire than that observed across the genus Bicyclus. We studied in detail one such mutation, of codominant effect, that causes the elimination of two adjacent eyespots on the ventral hindwing. By analyzing the expression of genes known to be involved in eyespot formation, we found an alteration in the differentiation of the "organizing" cells at the eyespot's center. No such cells differentiate in the wing subdivisions lacking the two eyespots in the mutants. We propose several developmental models, based on wing compartmentalization in Drosophila, that provide the first framework for thinking about the molecular evolution of butterfly wing pattern modularity.

The repressor activity of Even-skipped is highly conserved, and is sufficient to activate engrailed and to regulate both the spacing and stability of parasegment boundaries

During segmentation of the Drosophila embryo, even skipped is required to activate engrailed stripes and to organize odd-numbered parasegments. A 16 kb transgene containing the even skipped coding region can rescue normal engrailed expression, as well as all other aspects of segmentation, in even skipped null mutants. To better understand its mechanism of action, we functionally dissected the Even-skipped protein in the context of this transgene. We found that Even-skipped utilizes two repressor domains to carry out its function. Each of these domains can function autonomously in embryos when fused with the Gal4 DNA-binding domain. A chimeric protein consisting only of the Engrailed repressor domain and the Even-skipped homeodomain, but not the homeodomain alone, was able to restore function, indicating that the repression of target genes is sufficient for even skipped function at the blastoderm stage, while the homeodomain is sufficient to recognize those target genes. When Drosophila Even skipped was replaced by its homologs from other species, including a mouse homolog, they could provide substantial function, indicating that these proteins can recognize similar target sites and also provide repressor activity. Using this rescue system, we show that broad, early even skipped stripes are sufficient for activation of both odd- and even-numbered engrailed stripes. Furthermore, these ‘unrefined’ stripes organize odd-numbered parasegments in a dose-dependent manner, while the refined, late stripes, which coincide cell-for-cell with parasegment boundaries, are required to ensure the stability of the boundaries.

Cis elements of the villin gene control expression in restricted domains of the vertical (crypt) and horizontal (duodenum, cecum) axes of the intestine

Villin, an actin bundling protein found in the apical brush border of absorptive tissues, is one of the first structural genes to be transcriptionally activated in the embryonic intestinal endoderm. In the adult, villin is broadly expressed in every cell of the intestinal epithelium on both the vertical axis (crypt to villus tip) and the horizontal axis (duodenum through colon) of the intestine. Here, we document that a 12.4-kilobase region of the mouse villin gene drives high level expression of two different reporter genes (LacZ and Cre recombinase) within the entire intestinal epithelium of transgenic mice. Deletion of a portion of this transgene results in reduction of beta-galactosidase activity in restricted domains of the small intestine (duodenum) and large intestine (cecum). In addition, expression is reduced in the crypt compartment throughout the intestine. Thus, the global expression pattern of villin in the intestine is apparently the consequence of an amalgam of distinct and individual domain-specific control processes. That is, expression of villin in the duodenum and cecum requires different regulatory sequences than the rest of the intestine, and the expression of villin in crypts is regulated by different circuitry than expression of villin on villus tips.

Exploiting transcription factor binding site clustering to identify cis-regulatory modules involved in pattern formation in the Drosophila genome

A major challenge in interpreting genome sequences is understanding how the genome encodes the information that specifies when and where a gene will be expressed. The first step in this process is the identification of regions of the genome that contain regulatory information. In higher eukaryotes, this cis-regulatory information is organized into modular units [cis-regulatory modules (CRMs)] of a few hundred base pairs. A common feature of these cis-regulatory modules is the presence of multiple binding sites for multiple transcription factors. Here, we evaluate the extent to which the tendency for transcription factor binding sites to be clustered can be used as the basis for the computational identification of cis-regulatory modules. By using published DNA binding specificity data for five transcription factors active in the early Drosophila embryo, we identified genomic regions containing unusually high concentrations of predicted binding sites for these factors. A significant fraction of these binding site clusters overlap known CRMs that are regulated by these factors. In addition, many of the remaining clusters are adjacent to genes expressed in a pattern characteristic of genes regulated by these factors. We tested one of the newly identified clusters, mapping upstream of the gap gene giant (gt), and show that it acts as an enhancer that recapitulates the posterior expression pattern of gt.

Assessing the impact of comparative genomic sequence data on the functional annotation of the Drosophila genome

BACKGROUND

It is widely accepted that comparative sequence data can aid the functional annotation of genome sequences; however, the most informative species and features of genome evolution for comparison remain to be determined.

RESULTS

We analyzed conservation in eight genomic regions (apterous, even-skipped, fushi tarazu, twist, and Rhodopsins 1, 2, 3 and 4) from four Drosophila species (D. erecta, D. pseudoobscura, D. willistoni, and D. littoralis) covering more than 500 kb of the D. melanogaster genome. All D. melanogaster genes (and 78-82% of coding exons) identified in divergent species such as D. pseudoobscura show evidence of functional constraint. Addition of a third species can reveal functional constraint in otherwise non-significant pairwise exon comparisons. Microsynteny is largely conserved, with rearrangement breakpoints, novel transposable element insertions, and gene transpositions occurring in similar numbers. Rates of amino-acid substitution are higher in uncharacterized genes relative to genes that have previously been studied. Conserved non-coding sequences (CNCSs) tend to be spatially clustered with conserved spacing between CNCSs, and clusters of CNCSs can be used to predict enhancer sequences.

CONCLUSIONS

Our results provide the basis for choosing species whose genome sequences would be most useful in aiding the functional annotation of coding and cis-regulatory sequences in Drosophila. Furthermore, this work shows how decoding the spatial organization of conserved sequences, such as the clustering of CNCSs, can complement efforts to annotate eukaryotic genomes on the basis of sequence conservation alone.

Groucho augments the repression of multiple Even skipped target genes in establishing parasegment boundaries

Groucho acts as a co-repressor for several Drosophila DNA binding transcriptional repressors. Several of these proteins have been found to contain both Groucho-dependent and -independent repression domains, but the extent to which this distinction has functional consequences for the regulation of different target genes is not known. The product of the pair-rule gene even skipped has previously been shown to contain a Groucho-independent repression activity. In the Even skipped protein, outside the Groucho-independent repression domain, we have identified a conserved C-terminal motif (LFKPY), similar to motifs that mediate Groucho interaction in Hairy, Runt and Huckebein. Even skipped interacts with Groucho in yeast and in vitro, and groucho and even skipped genetically interact in vivo. Even skipped with a mutated Groucho interaction motif, which abolished binding to Groucho, showed a significantly reduced ability to rescue the even skipped null phenotype when driven by the complete even skipped regulatory region. Replacing this motif with a heterologous Groucho interaction motif restored the rescuing function of Even skipped in segmentation. Further functional assays demonstrated that the Even skipped C terminus acts as a Groucho-dependent repression domain in early Drosophila embryos. This novel repression domain was active on two target genes that are normally repressed by Even skipped at different concentrations, paired and sloppy paired. When the Groucho interaction motif is mutated, repression of each target gene is reduced to a similar extent, with some activity remaining. Thus, the ability of Even skipped to repress different target genes at different concentrations does not appear to involve differential recruitment or function of Groucho. The accumulation of multiple domains of similar function within a single protein may be a common evolutionary mechanism that fine-tunes the level of activity for different regulatory functions.

Ras pathway specificity is determined by the integration of multiple signal-activated and tissue-restricted transcription factors

Ras signaling elicits diverse outputs, yet how Ras specificity is generated remains incompletely understood. We demonstrate that Wingless (Wg) and Decapentaplegic (Dpp) confer competence for receptor tyrosine kinase-mediated induction of a subset of Drosophila muscle and cardiac progenitors by acting both upstream of and in parallel to Ras. In addition to regulating the expression of proximal Ras pathway components, Wg and Dpp coordinate the direct effects of three signal-activated (dTCF, Mad, and Pointed-functioning in the Wg, Dpp, and Ras/MAPK pathways, respectively) and two tissue-restricted (Twist and Tinman) transcription factors on a progenitor identity gene enhancer. The integration of Pointed with the combinatorial effects of dTCF, Mad, Twist, and Tinman determines inductive Ras signaling specificity in muscle and heart development.

Neural cell fate in rca1 and cycA mutants: the roles of intrinsic and extrinsic factors in asymmetric division in the Drosophila central nervous system

In the central nervous system (CNS) of Drosophila embryos lacking regulator of cyclin A (rca1) or cyclin A, we observe that several ganglion mother cells (GMCs) fail to divide. Whereas GMCs normally produce two sibling neurons that acquire different fates ('A/B'), non-dividing GMCs differentiate exclusively in the manner of one of their progeny ('B'). In zygotic numb mutants, sibling neuron fate alterations ('A/B' to 'A/A') occur infrequently or do not occur in some sibling pairs; we have determined that depletion of both maternal and zygotic numb causes sibling neurons to acquire equalized fates ('A/A') with near-complete expressivity. In rca1, numb mutant embryos, we observe binary cell fate changes ('B' to 'A') in several GMCs as well. Finally, we have demonstrated that expression of Delta in the mesoderm is sufficient to attain both sibling fates. Our results indicate that the intrinsic determinant Numb is absolutely required to attain differential sibling neuron fates. While the extrinsic factors Notch and Delta are also required to attain both fates, our results indicate that Delta signal can be received from outside the sibling pair.