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

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

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.

Geometric control of myosin II orientation during axis elongation

The actomyosin cytoskeleton is a crucial driver of morphogenesis. Yet how the behavior of large-scale cytoskeletal patterns in deforming tissues emerges from the interplay of geometry, genetics, and mechanics remains incompletely understood. Convergent extension in Drosophila melanogaster embryos provides the opportunity to establish a quantitative understanding of the dynamics of anisotropic non-muscle myosin II. Cell-scale analysis of protein localization in fixed embryos suggests that gene expression patterns govern myosin anisotropy via complex rules. However, technical limitations have impeded quantitative and dynamic studies of this process at the whole embryo level, leaving the role of geometry open. Here, we combine in toto live imaging with quantitative analysis of molecular dynamics to characterize the distribution of myosin anisotropy and the corresponding genetic patterning. We found pair rule gene expression continuously deformed, flowing with the tissue frame. In contrast, myosin anisotropy orientation remained approximately static and was only weakly deflected from the stationary dorsal-ventral axis of the embryo. We propose that myosin is recruited by a geometrically defined static source, potentially related to the embryo-scale epithelial tension, and account for transient deflections by cytoskeletal turnover and junction reorientation by flow. With only one parameter, this model quantitatively accounts for the time course of myosin anisotropy orientation in wild-type, twist, and even-skipped embryos, as well as embryos with perturbed egg geometry. Geometric patterning of the cytoskeleton suggests a simple physical strategy to ensure a robust flow and formation of shape.

The power of the (imperfect) palindrome: Sequence-specific roles of palindromic motifs in gene regulation

In human languages, a palindrome reads the same forward as backward (e.g., 'madam'). In regulatory DNA, a palindrome is an inverted sequence repeat that allows a transcription factor to bind as a homodimer or as a heterodimer with another type of transcription factor. Regulatory palindromes are typically imperfect, that is, the repeated sequences differ in at least one base pair, but the functional significance of this asymmetry remains poorly understood. Here, we review the use of imperfect palindromes in Drosophila photoreceptor differentiation and mammalian steroid receptor signaling. Moreover, we discuss mechanistic explanations for the predominance of imperfect palindromes over perfect palindromes in these two gene regulatory contexts. Lastly, we propose to elucidate whether specific imperfectly palindromic variants have specific regulatory functions in steroid receptor signaling and whether such variants can help predict transcriptional outcomes as well as the response of individual patients to drug treatments.

A combinatorial cis-regulatory logic restricts color-sensing Rhodopsins to specific photoreceptor subsets in Drosophila

Color vision in Drosophila melanogaster is based on the expression of five different color-sensing Rhodopsin proteins in distinct subtypes of photoreceptor neurons. Promoter regions of less than 300 base pairs are sufficient to reproduce the unique, photoreceptor subtype-specific rhodopsin expression patterns. The underlying cis-regulatory logic remains poorly understood, but it has been proposed that the rhodopsin promoters have a bipartite structure: the distal promoter region directs the highly restricted expression in a specific photoreceptor subtype, while the proximal core promoter region provides general activation in all photoreceptors. Here, we investigate whether the rhodopsin promoters exhibit a strict specialization of their distal (subtype specificity) and proximal (general activation) promoter regions, or if both promoter regions contribute to generating the photoreceptor subtype-specific expression pattern. To distinguish between these two models, we analyze the expression patterns of a set of hybrid promoters that combine the distal promoter region of one rhodopsin with the proximal core promoter region of another rhodopsin. We find that the function of the proximal core promoter regions extends beyond providing general activation: these regions play a previously underappreciated role in generating the non-overlapping expression patterns of the different rhodopsins. Therefore, cis-regulatory motifs in both the distal and the proximal core promoter regions recruit transcription factors that generate the unique rhodopsin patterns in a combinatorial manner. We compare this combinatorial regulatory logic to the regulatory logic of olfactory receptor genes and discuss potential implications for the evolution of rhodopsins.

Each type of sensory receptor neuron in our body expresses a specific sensory receptor protein, which allows us to detect and discriminate a variety of environmental stimuli. The regulatory logic that controls this spatially precise and highly restricted expression of sensory receptor proteins remains poorly understood. As a model system, we study the mechanisms that control the expression of different color-sensing Rhodopsin proteins in distinct subtypes of Drosophila photoreceptors, which is the basis for color vision. Compact promoter regions of less than 300 base pairs are sufficient to reproduce the non-overlapping rhodopsin patterns. However, the regulatory logic that underlies the combination (sometimes called ‘grammar’) of the cis-regulatory motifs (sometimes called ‘vocabulary’) within the rhodopsin promoters remains poorly understood. Here, we find that specific combinations of cis-regulatory motifs in the distal and the proximal core promoter regions of each rhodopsin direct its unique expression pattern.

Formation of polarized contractile interfaces by self-organized Toll-8/Cirl GPCR asymmetry

Interfaces between cells with distinct genetic identities elicit signals to organize local cell behaviors driving tissue morphogenesis. The Drosophila embryonic axis extension requires planar polarized enrichment of myosin-II powering oriented cell intercalations. Myosin-II levels are quantitatively controlled by GPCR signaling, whereas myosin-II polarity requires patterned expression of several Toll receptors. How Toll receptors polarize myosin-II and how this involves GPCRs remain unknown. Here, we report that differential expression of a single Toll receptor, Toll-8, polarizes myosin-II through binding to the adhesion GPCR Cirl/latrophilin. Asymmetric expression of Cirl is sufficient to enrich myosin-II, and Cirl localization is asymmetric at Toll-8 expression boundaries. Exploring the process dynamically, we reveal that Toll-8 and Cirl exhibit mutually dependent planar polarity in response to quantitative differences in Toll-8 expression between neighboring cells. Collectively, we propose that the cell surface protein complex Toll-8/Cirl self-organizes to generate local asymmetric interfaces essential for planar polarization of contractility.

•

Asymmetric expression of a single Toll receptor leads to Myo-II polarization

•

The adhesion GPCR Cirl binds to Toll-8 mediating Toll-8-induced Myo-II polarization

•

Toll-8 boundaries generate a Cirl interfacial asymmetry that can polarize Myo-II

•

Differences in Toll-8 levels lead to interdependent Toll-8 and Cirl planar polarity

Toll receptors remodel epithelia by directing planar-polarized Src and PI3K activity

Toll-like receptors are essential for animal development and survival, with conserved roles in innate immunity, tissue patterning, and cell behavior. The mechanisms by which Toll receptors signal to the nucleus are well characterized, but how Toll receptors generate rapid, localized signals at the cell membrane to produce acute changes in cell polarity and behavior is not known. We show that Drosophila Toll receptors direct epithelial convergent extension by inducing planar-polarized patterns of Src and PI3-kinase (PI3K) activity. Toll receptors target Src activity to specific sites at the membrane, and Src recruits PI3K to the Toll-2 complex through tyrosine phosphorylation of the Toll-2 cytoplasmic domain. Reducing Src or PI3K activity disrupts planar-polarized myosin assembly, cell intercalation, and convergent extension, whereas constitutive Src activity promotes ectopic PI3K and myosin cortical localization. These results demonstrate that Toll receptors direct cell polarity and behavior by locally mobilizing Src and PI3K activity.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Cooperative binding between distant transcription factors is a hallmark of active enhancers

Enhancers harbor binding motifs that recruit transcription factors (TFs) for gene activation. While cooperative binding of TFs at enhancers is known to be critical for transcriptional activation of a handful of developmental enhancers, the extent of TF cooperativity genome-wide is unknown. Here, we couple high-resolution nuclease footprinting with single-molecule methylation profiling to characterize TF cooperativity at active enhancers in the Drosophila genome. Enrichment of short micrococcal nuclease (MNase)-protected DNA segments indicates that the majority of enhancers harbor two or more TF-binding sites, and we uncover protected fragments that correspond to co-bound sites in thousands of enhancers. From the analysis of co-binding, we find that cooperativity dominates TF binding in vivo at the majority of active enhancers. Cooperativity is highest between sites spaced 50 bp apart, indicating that cooperativity occurs without apparent protein-protein interactions. Our findings suggest nucleosomes promoting cooperativity because co-binding may effectively clear nucleosomes and promote enhancer function.

Kinetic sculpting of the seven stripes of the Drosophila even-skipped gene

We used live imaging to visualize the transcriptional dynamics of the Drosophila melanogaster even-skipped gene at single-cell and high-temporal resolution as its seven stripe expression pattern forms, and developed tools to characterize and visualize how transcriptional bursting varies over time and space. We find that despite being created by the independent activity of five enhancers, even-skipped stripes are sculpted by the same kinetic phenomena: a coupled increase of burst frequency and amplitude. By tracking the position and activity of individual nuclei, we show that stripe movement is driven by the exchange of bursting nuclei from the posterior to anterior stripe flanks. Our work provides a conceptual, theoretical and computational framework for dissecting pattern formation in space and time, and reveals how the coordinated transcriptional activity of individual nuclei shapes complex developmental patterns.

A sensitive mNeonGreen reporter system to measure transcriptional dynamics in Drosophila development

The gene regulatory network governing anterior–posterior axis formation in Drosophila is a well-established paradigm to study transcription in developmental biology. The rapid temporal dynamics of gene expression during early stages of development, however, are difficult to track with standard techniques. We optimized the bright and fast-maturing fluorescent protein mNeonGreen as a real-time, quantitative reporter of enhancer expression. We derive enhancer activity from the reporter fluorescence dynamics with high spatial and temporal resolution, using a robust reconstruction algorithm. By comparing our results with data obtained with the established MS2-MCP system, we demonstrate the higher detection sensitivity of our reporter. We used the reporter to quantify the activity of variants of a simple synthetic enhancer, and observe increased activity upon reduction of enhancer–promoter distance or addition of binding sites for the pioneer transcription factor Zelda. Our reporter system constitutes a powerful tool to study spatio-temporal gene expression dynamics in live embryos.

Ceolin et al. optimise the fluorescent protein mNeonGreen to generate a reporter to measure enhancer expression in living Drosophila embryos. They also find that reduced enhancer–promoter distance or addition of binding sites for the pioneer transcription factor Zelda increases enhancer activity.

Advances in visualizing transcription factor - DNA interactions

At the heart of the transcription process is the specific interaction between transcription factors (TFs) and their target DNA sequences. Decades of molecular biology research have led to unprecedented insights into how TFs access the genome to regulate transcription. In the last 20 years, advances in microscopy have enabled scientists to add imaging as a powerful tool in probing two specific aspects of TF-DNA interactions: structure and dynamics. In this review, we examine how applications of diverse imaging technologies can provide structural and dynamic information that complements insights gained from molecular biology assays. As a case study, we discuss how applications of advanced imaging techniques have reshaped our understanding of TF behavior across the cell cycle, leading to a rethinking in the field of mitotic bookmarking.

A Mutation in the Drosophila melanogaster eve Stripe 2 Minimal Enhancer Is Buffered by Flanking Sequences

Enhancers are DNA sequences composed of transcription factor binding sites that drive complex patterns of gene expression in space and time. Until recently, studying enhancers in their genomic context was technically challenging. Therefore, minimal enhancers, the shortest pieces of DNA that can drive an expression pattern that resembles a gene's endogenous pattern, are often used to study features of enhancer function. However, evidence suggests that some enhancers require sequences outside the minimal enhancer to maintain function under environmental perturbations. We hypothesized that these additional sequences also prevent misexpression caused by a transcription factor binding site mutation within a minimal enhancer. Using the Drosophila melanogastereven-skipped stripe 2 enhancer as a case study, we tested the effect of a Giant binding site mutation (gt-2) on the expression patterns driven by minimal and extended enhancer reporter constructs. We found that, in contrast to the misexpression caused by the gt-2 binding site deletion in the minimal enhancer, the same gt-2 binding site deletion in the extended enhancer did not have an effect on expression. The buffering of expression levels, but not expression pattern, is partially explained by an additional Giant binding site outside the minimal enhancer. Deleting the gt-2 binding site in the endogenous locus had no significant effect on stripe 2 expression. Our results indicate that rules derived from mutating enhancer reporter constructs may not represent what occurs in the endogenous context.

Embryonic canalization and its limits-A view from temperature

Many animals are able to produce similar offspring over a range of environmental conditions. This property of the developmental process has been termed canalization-the channeling of developmental pathways to generate a stable outcome despite varying conditions. Temperature is one environmental parameter that has fundamental effects on cell physiology and biochemistry, yet developmental programs generally result in a stable phenotype under a range of temperatures. On the other hand, there are typically upper and lower temperature limits beyond which the developmental program is unable to produce normal offspring. This review summarizes data on how development is affected by temperature, particularly high temperature, in various animal species. It also brings together information on potential cell biological and developmental genetic factors that may be responsible for developmental stability in varying temperatures, and likely critical mechanisms that break down at high temperature. Also reviewed are possible means for studying temperature effects on embryogenesis and how to determine which factors are most critical at the high-temperature limits for normal development. Increased knowledge of these critical factors will point to the targets of selection under climate change, and more generally, how developmental robustness in varying environments is maintained.

Multimodal transcriptional control of pattern formation in embryonic development

Predicting how interactions between transcription factors and regulatory DNA sequence dictate rates of transcription and, ultimately, drive developmental outcomes remains an open challenge in physical biology. Using stripe 2 of the even-skipped gene in Drosophila embryos as a case study, we dissect the regulatory forces underpinning a key step along the developmental decision-making cascade: the generation of cytoplasmic mRNA patterns via the control of transcription in individual cells. Using live imaging and computational approaches, we found that the transcriptional burst frequency is modulated across the stripe to control the mRNA production rate. However, we discovered that bursting alone cannot quantitatively recapitulate the formation of the stripe and that control of the window of time over which each nucleus transcribes even-skipped plays a critical role in stripe formation. Theoretical modeling revealed that these regulatory strategies (bursting and the time window) respond in different ways to input transcription factor concentrations, suggesting that the stripe is shaped by the interplay of 2 distinct underlying molecular processes.

An in silico analysis of robust but fragile gene regulation links enhancer length to robustness

Organisms must ensure that expression of genes is directed to the appropriate tissues at the correct times, while simultaneously ensuring that these gene regulatory systems are robust to perturbation. This idea is captured by a mathematical concept called r-robustness, which says that a system is robust to a perturbation in up to r - 1 randomly chosen parameters. r-robustness implies that the biological system has a small number of sensitive parameters and that this number can be used as a robustness measure. In this work we use this idea to investigate the robustness of gene regulation using a sequence level model of the Drosophila melanogaster gene even-skipped. We consider robustness with respect to mutations of the enhancer sequence and with respect to changes of the transcription factor concentrations. We find that gene regulation is r-robust with respect to mutations in the enhancer sequence and identify a number of sensitive nucleotides. In both natural and in silico predicted enhancers, the number of nucleotides that are sensitive to mutation correlates negatively with the length of the sequence, meaning that longer sequences are more robust. The exact degree of robustness obtained is dependent not only on DNA sequence, but also on the local concentration of regulatory factors. We find that gene regulation can be remarkably sensitive to changes in transcription factor concentrations at the boundaries of expression features, while it is robust to perturbation elsewhere.

Evolution: How Many Phenotypes Do Regulatory Mutations Affect?

Mutations in gene regulatory regions are thought to play an important role in the evolution of morphological structures. This is largely due to their minimal pleiotropic effects, limiting their impact to one particular body part. A recent study finds that one such regulatory mutation may affect two particular morphological structures.

An LRR Receptor-Teneurin System Directs Planar Polarity at Compartment Boundaries

Epithelial cells dynamically self-organize in response to extracellular spatial cues relayed by cell-surface receptors. During convergent extension in Drosophila, Toll-related receptors direct planar polarized cell rearrangements that elongate the head-to-tail axis. However, many cells establish polarity in the absence of Toll receptor activity, indicating the presence of additional spatial cues. Here we demonstrate that the leucine-rich-repeat receptor Tartan and the teneurin Ten-m provide critical polarity signals at epithelial compartment boundaries. The Tartan and Ten-m extracellular domains interact in vitro, and Tartan promotes Ten-m localization to compartment boundaries in vivo. We show that Tartan and Ten-m are necessary for the planar polarity and organization of compartment boundary cells. Moreover, ectopic stripes of Tartan and Ten-m are sufficient to induce myosin accumulation at stripe boundaries. These results demonstrate that the Tartan/Ten-m and Toll receptor systems together create a high-resolution network of spatial cues that guides cell behavior during convergent extension.

Structure and cis-regulatory analysis of a Drosophila grainyhead neuroblast enhancer

Evolutionary analysis of cis-regulatory DNA reveals that enhancers consist of clusters of conserved sequence blocks (CSBs) that are made up of both unique and repeated sequence elements. This study seeks to address the basis for spatial and temporal regulation of neuroblas. A search for temporally restricted CNS NB enhancers identified one within the transcription factor grainyhead (grh) gene locus. The intronic enhancer, grh-15, contains two separable semi-autonomous activities, one that drives expression predominantly within the developing brain NBs and another in ventral cord NBs. To gain insight into the function of the CSBs constituting the brain-specific enhancer, we have systematically deleted each CSB and compared the activity of the altered enhancer to that of the full brain-specific enhancer. While our results indicate that information regulating enhancer activity is highly redundant, we have found that individual CSBs convey expression in subsets of larval lineages that are generated from either Type I or Type II NBs. These studies also highlight how evolutionary sequence conservation can be used as a guide the functional analysis of cis-regulatory DNA.

Expression of SULTR2;2, encoding a low-affinity sulphur transporter, in the Arabidopsis bundle sheath and vein cells is mediated by a positive regulator

The bundle sheath provides a conduit linking veins and mesophyll cells. In the C3 plant Arabidopsis thaliana, it also plays important roles in oxidative stress and sulphur metabolism. However, the mechanisms responsible for the patterns of gene expression that underpin these metabolic specializations are poorly understood. Here, we used the Arabidopsis SULTR2;2 gene as a model to better understand mechanisms that restrict expression to the bundle sheath. Deletion analysis indicated that the SULTR2;2 promoter contains a short region necessary for expression in the bundle sheath and veins. This sequence acts as a positive regulator and is tolerant to multiple consecutive deletions indicating considerable redundancy in the cis-elements involved. It is highly conserved in SULTR2;2 genes of the Brassicaceae and is functional in the distantly related C4 species Flaveria bidentis that belongs to the Asteraceae. We conclude that expression of SULTR2;2 in the bundle sheath and veins is underpinned by a highly redundant sequence that likely represents an ancient and conserved mechanism found in families as diverse as the Asteraceae and Brassicaceae.

Gene Regulatory Enhancers with Evolutionarily Conserved Activity Are More Pleiotropic than Those with Species-Specific Activity

Studies of regulatory activity and gene expression have revealed an intriguing dichotomy: There is substantial turnover in the regulatory activity of orthologous sequences between species; however, the expression level of orthologous genes is largely conserved. Understanding how distal regulatory elements, for example, enhancers, evolve and function is critical, as alterations in gene expression levels can drive the development of both complex disease and functional divergence between species. In this study, we investigated determinants of the conservation of regulatory enhancer activity for orthologous sequences across mammalian evolution. Using liver enhancers identified from genome-wide histone modification profiles in ten diverse mammalian species, we compared orthologous sequences that exhibited regulatory activity in all species (conserved-activity enhancers) to shared sequences active only in a single species (species-specific-activity enhancers). Conserved-activity enhancers have greater regulatory potential than species-specific-activity enhancers, as quantified by both the density and diversity of transcription factor binding motifs. Consistent with their greater regulatory potential, conserved-activity enhancers have greater regulatory activity in humans than species-specific-activity enhancers: They are active across more cellular contexts, and they regulate more genes than species-specific-activity enhancers. Furthermore, the genes regulated by conserved-activity enhancers are expressed in more tissues and are less tolerant of loss-of-function mutations than those targeted by species-specific-activity enhancers. These consistent results across various stages of gene regulation demonstrate that conserved-activity enhancers are more pleiotropic than their species-specific-activity counterparts. This suggests that pleiotropy is associated with the conservation of regulatory across mammalian evolution.

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

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

1st German Phage Symposium-Conference Report

In Germany, phage research and application can be traced back to the beginning of the 20th century. However, with the triumphal march of antibiotics around the world, the significance of bacteriophages faded in most countries, and respective research mainly focused on fundamental questions and niche applications. After a century, we pay tribute to the overuse of antibiotics that led to multidrug resistance and calls for new strategies to combat pathogenic microbes. Against this background, bacteriophages came into the spotlight of researchers and practitioners again resulting in a fast growing “phage community”. In October 2017, part of this community met at the 1st German Phage Symposium to share their knowledge and experiences. The participants discussed open questions and challenges related to phage therapy and the application of phages in general. This report summarizes the presentations given, highlights the main points of the round table discussion and concludes with an outlook for the different aspects of phage application.

Functional regulatory evolution outside of the minimal even-skipped stripe 2 enhancer

Transcriptional enhancers are regions of DNA that drive precise patterns of gene expression. Although many studies have elucidated how individual enhancers can evolve, most of this work has focused on what are called 'minimal' enhancers, the smallest DNA regions that drive expression that approximates an aspect of native gene expression. Here, we explore how the Drosophila erecta even-skipped (eve) locus has evolved by testing its activity in the divergent D. melanogaster genome. We found, as has been reported previously, that the D. erecta eve stripe 2 enhancer (eveS2) fails to drive appreciable expression in D. melanogaster However, we found that a large transgene carrying the entire D. erecta eve locus drives normal eve expression, including in stripe 2. We performed a functional dissection of the region upstream of the D. erecta eveS2 region and found multiple Zelda motifs that are required for normal expression. Our results illustrate how sequences outside of minimal enhancer regions can evolve functionally through mechanisms other than changes in transcription factor-binding sites that drive patterning.

Krox20 hindbrain regulation incorporates multiple modes of cooperation between cis-acting elements

Developmental genes can harbour multiple transcriptional enhancers that act simultaneously or in succession to achieve robust and precise spatiotemporal expression. However, the mechanisms underlying cooperation between cis-acting elements are poorly documented, notably in vertebrates. The mouse gene Krox20 encodes a transcription factor required for the specification of two segments (rhombomeres) of the developing hindbrain. In rhombomere 3, Krox20 is subject to direct positive feedback governed by an autoregulatory enhancer, element A. In contrast, a second enhancer, element C, distant by 70 kb, is active from the initiation of transcription independent of the presence of the KROX20 protein. Here, using both enhancer knock-outs and investigations of chromatin organisation, we show that element C possesses a dual activity: besides its classical enhancer function, it is also permanently required in cis to potentiate the autoregulatory activity of element A, by increasing its chromatin accessibility. This work uncovers a novel, asymmetrical, long-range mode of cooperation between cis-acting elements that might be essential to avoid promiscuous activation of positive autoregulatory elements.

The formation of multicellular organisms from the egg to the adult stage is largely under genetic control. The activation of specific genes is governed by regulatory DNA sequences present nearby on the chromosome. Most of these sequences promote activation and are called enhancers. In this paper, we study two enhancers governing the expression of a gene involved in the formation of the posterior brain in vertebrates. One of these enhancers is involved in a positive feedback loop: it is itself activated by the protein product of the gene that it regulates. The other enhancer was thought to be only involved in the initial accumulation of the protein, necessary for the subsequent activation of the feedback loop. Here we show that the second enhancer directly cooperates with the autoregulatory enhancer to increase its accessibility and its activity. Our work uncovers a novel, long-range mode of cooperation between enhancers that restricts the domain of action of autoregulatory enhancers within embryos and might be essential to avoid their inappropriate activation.

Kibra and Merlin Activate the Hippo Pathway Spatially Distinct from and Independent of Expanded

The Hippo pathway is emerging as a key evolutionarily conserved signaling mechanism that controls organ size. Three membrane-associated proteins, Kibra, Merlin, and Expanded, regulate pathway activity, but the precise molecular mechanism by which they function is still poorly understood. Here we provide evidence that Merlin and Kibra activate Hippo signaling in parallel to Expanded at a spatially distinct cellular domain, the medial apical cortex. Merlin and Kibra together recruit the adapter protein Salvador, which in turn recruits the core kinase Hippo. In addition, we show that Crumbs has a dual effect on Hippo signaling. Crumbs promotes the ability of Expanded to activate the pathway but also sequesters Kibra to downregulate Hippo signaling. Together, our findings elucidate the mechanism of Hippo pathway activation by Merlin and Kibra, identify a subcellular domain for Hippo pathway regulation, and demonstrate differential activity of upstream regulators in different subcellular domains.

SoxNeuro and Shavenbaby act cooperatively to shape denticles in the embryonic epidermis of Drosophila

During development, extracellular signals are integrated by cells to induce the transcriptional circuitry that controls morphogenesis. In the fly epidermis, Wingless (Wg)/Wnt signaling directs cells to produce either a distinctly shaped denticle or no denticle, resulting in a segmental pattern of denticle belts separated by smooth, or 'naked', cuticle. Naked cuticle results from Wg repression of shavenbaby (svb), which encodes a transcription factor required for denticle construction. We have discovered that although the svb promoter responds differentially to altered Wg levels, Svb alone cannot produce the morphological diversity of denticles found in wild-type belts. Instead, a second Wg-responsive transcription factor, SoxNeuro (SoxN), cooperates with Svb to shape the denticles. Co-expressing ectopic SoxN with svb rescued diverse denticle morphologies. Conversely, removing SoxN activity eliminated the residual denticles found in svb mutant embryos. Furthermore, several known Svb target genes are also activated by SoxN, and we have discovered two novel target genes of SoxN that are expressed in denticle-producing cells and that are regulated independently of Svb. We conclude that proper denticle morphogenesis requires transcriptional regulation by both SoxN and Svb.

Evolution of New cis-Regulatory Motifs Required for Cell-Specific Gene Expression in Caenorhabditis

Patterning of C. elegans vulval cell fates relies on inductive signaling. In this induction event, a single cell, the gonadal anchor cell, secretes LIN-3/EGF and induces three out of six competent precursor cells to acquire a vulval fate. We previously showed that this developmental system is robust to a four-fold variation in lin-3/EGF genetic dose. Here using single-molecule FISH, we find that the mean level of expression of lin-3 in the anchor cell is remarkably conserved. No change in lin-3 expression level could be detected among C. elegans wild isolates and only a low level of change—less than 30%—in the Caenorhabditis genus and in Oscheius tipulae. In C. elegans, lin-3 expression in the anchor cell is known to require three transcription factor binding sites, specifically two E-boxes and a nuclear-hormone-receptor (NHR) binding site. Mutation of any of these three elements in C. elegans results in a dramatic decrease in lin-3 expression. Yet only a single E-box is found in the Drosophilae supergroup of Caenorhabditis species, including C. angaria, while the NHR-binding site likely only evolved at the base of the Elegans group. We find that a transgene from C. angaria bearing a single E-box is sufficient for normal expression in C. elegans. Even a short 58 bp cis-regulatory fragment from C. angaria with this single E-box is able to replace the three transcription factor binding sites at the endogenous C. elegans lin-3 locus, resulting in the wild-type expression level. Thus, regulatory evolution occurring in cis within a 58 bp lin-3 fragment, results in a strict requirement for the NHR binding site and a second E-box in C. elegans. This single-cell, single-molecule, quantitative and functional evo-devo study demonstrates that conserved expression levels can hide extensive change in cis-regulatory site requirements and highlights the evolution of new cis-regulatory elements required for cell-specific gene expression.

Diversification of mechanisms regulating gene expression of key developmental factors is a major force in the evolution of development. However, in the past, comparisons of gene expression across different species have often been qualitative (i.e. ‘expression is on versus off’ in a certain cell) without precise quantification. New experimental methods now allow us to quantitatively compare the expression of gene homologs across species, with single cell resolution. Moreover, the development of genome editing tools enables the dissection of regulatory DNA sequences that drive gene expression. We use here a well-established “textbook” example of animal organogenesis in the microscopic nematode, Caenorhabditis elegans, focusing on the expression of lin-3, coding for the main inducer of the vulva, in a single cell called the anchor cell. We find that the lin-3 expression level is remarkably conserved, with 20–25 messenger RNAs per anchor cell, in species that are molecularly as distant as fish and mammals. This conservation occurs despite substantial changes and compensation in the regulatory elements required for cell-specific gene expression.

Early developmental gene enhancers affect subcortical volumes in the adult human brain

Genome-wide association screens aim to identify common genetic variants contributing to the phenotypic variability of complex traits, such as human height or brain morphology. The identified genetic variants are mostly within noncoding genomic regions and the biology of the genotype-phenotype association typically remains unclear. In this article, we propose a complementary targeted strategy to reveal the genetic underpinnings of variability in subcortical brain volumes, by specifically selecting genomic loci that are experimentally validated forebrain enhancers, active in early embryonic development. We hypothesized that genetic variation within these enhancers may affect the development and ultimately the structure of subcortical brain regions in adults. We tested whether variants in forebrain enhancer regions showed an overall enrichment of association with volumetric variation in subcortical structures of >13,000 healthy adults. We observed significant enrichment of genomic loci that affect the volume of the hippocampus within forebrain enhancers (empirical P = 0.0015), a finding which robustly passed the adjusted threshold for testing of multiple brain phenotypes (cutoff of P < 0.0083 at an alpha of 0.05). In analyses of individual single nucleotide polymorphisms (SNPs), we identified an association upstream of the ID2 gene with rs7588305 and variation in hippocampal volume. This SNP-based association survived multiple-testing correction for the number of SNPs analyzed but not for the number of subcortical structures. Targeting known regulatory regions offers a way to understand the underlying biology that connects genotypes to phenotypes, particularly in the context of neuroimaging genetics. This biology-driven approach generates testable hypotheses regarding the functional biology of identified associations. Hum Brain Mapp 37:1788-1800, 2016. © 2016 Wiley Periodicals, Inc.

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

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

Common binding by redundant group B Sox proteins is evolutionarily conserved in Drosophila

Background

Group B Sox proteins are a highly conserved group of transcription factors that act extensively to coordinate nervous system development in higher metazoans while showing both co-expression and functional redundancy across a broad group of taxa. In Drosophila melanogaster, the two group B Sox proteins Dichaete and SoxNeuro show widespread common binding across the genome. While some instances of functional compensation have been observed in Drosophila, the function of common binding and the extent of its evolutionary conservation is not known.

Results

We used DamID-seq to examine the genome-wide binding patterns of Dichaete and SoxNeuro in four species of Drosophila. Through a quantitative comparison of Dichaete binding, we evaluated the rate of binding site turnover across the genome as well as at specific functional sites. We also examined the presence of Sox motifs within binding intervals and the correlation between sequence conservation and binding conservation. To determine whether common binding between Dichaete and SoxNeuro is conserved, we performed a detailed analysis of the binding patterns of both factors in two species.

Conclusion

We find that, while the regulatory networks driven by Dichaete and SoxNeuro are largely conserved across the drosophilids studied, binding site turnover is widespread and correlated with phylogenetic distance. Nonetheless, binding is preferentially conserved at known cis-regulatory modules and core, independently verified binding sites. We observed the strongest binding conservation at sites that are commonly bound by Dichaete and SoxNeuro, suggesting that these sites are functionally important. Our analysis provides insights into the evolution of group B Sox function, highlighting the specific conservation of shared binding sites and suggesting alternative sources of neofunctionalisation between paralogous family members.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1495-3) contains supplementary material, which is available to authorized users.

Rethinking gene regulatory networks in light of alternative splicing, intrinsically disordered protein domains, and post-translational modifications

Models for genetic regulation and cell fate specification characteristically assume that gene regulatory networks (GRNs) are essentially deterministic and exhibit multiple stable states specifying alternative, but pre-figured cell fates. Mounting evidence shows, however, that most eukaryotic precursor RNAs undergo alternative splicing (AS) and that the majority of transcription factors contain intrinsically disordered protein (IDP) domains whose functionalities are context dependent as well as subject to post-translational modification (PTM). Consequently, many transcription factors do not have fixed cis-acting regulatory targets, and developmental determination by GRNs alone is untenable. Modeling these phenomena requires a multi-scale approach to explain how GRNs operationally interact with the intra- and intercellular environments. Evidence shows that AS, IDP, and PTM complicate gene expression and act synergistically to facilitate and promote time- and cell-specific protein modifications involved in cell signaling and cell fate specification and thereby disrupt a strict deterministic GRN-phenotype mapping. The combined effects of AS, IDP, and PTM give proteomes physiological plasticity, adaptive responsiveness, and developmental versatility without inefficiently expanding genome size. They also help us understand how protein functionalities can undergo major evolutionary changes by buffering mutational consequences.

Evolutionary meandering of intermolecular interactions along the drift barrier

Many cellular functions depend on highly specific intermolecular interactions, for example transcription factors and their DNA binding sites, microRNAs and their RNA binding sites, the interfaces between heterodimeric protein molecules, the stems in RNA molecules, and kinases and their response regulators in signal-transduction systems. Despite the need for complementarity between interacting partners, such pairwise systems seem to be capable of high levels of evolutionary divergence, even when subject to strong selection. Such behavior is a consequence of the diminishing advantages of increasing binding affinity between partners, the multiplicity of evolutionary pathways between selectively equivalent alternatives, and the stochastic nature of evolutionary processes. Because mutation pressure toward reduced affinity conflicts with selective pressure for greater interaction, situations can arise in which the expected distribution of the degree of matching between interacting partners is bimodal, even in the face of constant selection. Although biomolecules with larger numbers of interacting partners are subject to increased levels of evolutionary conservation, their more numerous partners need not converge on a single sequence motif or be increasingly constrained in more complex systems. These results suggest that most phylogenetic differences in the sequences of binding interfaces are not the result of adaptive fine tuning but a simple consequence of random genetic drift.

A positional Toll receptor code directs convergent extension in Drosophila

Elongation of the head-to-tail body axis by convergent extension is a conserved developmental process throughout metazoans. In Drosophila, patterns of transcription factor expression provide spatial cues that induce systematically oriented cell movements and promote tissue elongation. However, the mechanisms by which patterned transcriptional inputs control cell polarity and behaviour have long been elusive. We demonstrate that three Toll family receptors, Toll-2, Toll-6 and Toll-8, are expressed in overlapping transverse stripes along the anterior-posterior axis and act in combination to direct planar polarity and polarized cell rearrangements during convergent extension. Simultaneous disruption of all three receptors strongly reduces actomyosin-driven junctional remodelling and axis elongation, and an ectopic stripe of Toll receptor expression is sufficient to induce planar polarized actomyosin contractility. These results demonstrate that tissue-level patterns of Toll receptor expression provide spatial signals that link positional information from the anterior-posterior patterning system to the essential cell behaviours that drive convergent extension.

The embryo as a laboratory: quantifying transcription in Drosophila

Transcriptional regulation of gene expression is fundamental to most cellular processes, including determination of cellular fates. Quantitative studies of transcription in cultured cells have led to significant advances in identifying mechanisms underlying transcriptional control. Recent progress allowed implementation of these same quantitative methods in multicellular organisms to ask how transcriptional regulation unfolds both in vivo and at the single molecule level in the context of embryonic development. Here we review some of these advances in early Drosophila development, which bring the embryo on par with its single celled counterparts. In particular, we discuss progress in methods to measure mRNA and protein distributions in fixed and living embryos, and we highlight some initial applications that lead to fundamental new insights about molecular transcription processes. We end with an outlook on how to further exploit the unique advantages that come with investigating transcriptional control in the multicellular context of development.

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.

Natural Populations of Drosophila melanogaster Reveal Features of an Uncharacterized Circadian Property: The Lower Temperature Limit of Rhythmicity

Most cyclic biological processes are under control of a circadian molecular timing system that synchronizes these phenomena to the 24-h day. One generic property of circadian-controlled processes is that they operate within a specific temperature range, below which the manifestation of rhythm ceases. Little is known about the evolutionary relevance of the lower temperature limit of rhythmicity or about the mechanism underlying the loss of overt circadian behavior below this lower limit, especially in one model organism of chronobiology, Drosophila melanogaster. Natural populations of Drosophila are evolving under divergent selection pressures and so provide a source of diversity necessary to address these issues. Using lines derived from African populations, we find that there is natural variation in the expression of rhythmic behavior under low-temperature conditions. We found evidence that this variability is evolutionarily relevant at extremely low temperature (12 °C) because high-altitude populations exhibit selection for locally adapted genomes that contribute to rhythmic behavior. Lines resistant to 15 °C show an additional layer of diversity in their response to temperature extremes because some lines are resistant to low temperature (15 °C) only, whereas others are cross-resistant to high and low temperature (15 °C and 30 °C). Genetic analysis of one cold-resistant circadian line at 15 °C reveals that the phenotype maps to the X-chromosome but not to the core clock genes, per and sgg. Analysis of the central clock cells of this line reveals that maintenance of rhythm is associated with robust clock function, which is compromised in a standard laboratory strain. These data indicate that the cold-resistant circadian phenotype is clock based. This study highlights the importance of using natural populations to inform us of the basic features of circadian traits, especially those that might be under temperature-based selection.

Development of five digits is controlled by a bipartite long-range cis-regulator

Conservation within intergenic DNA often highlights regulatory elements that control gene expression from a long range. How conservation within a single element relates to regulatory information and how internal composition relates to function is unknown. Here, we examine the structural features of the highly conserved ZRS (also called MFCS1) cis-regulator responsible for the spatiotemporal control of Shh in the limb bud. By systematically dissecting the ZRS, both in transgenic assays and within in the endogenous locus, we show that the ZRS is, in effect, composed of two distinct domains of activity: one domain directs spatiotemporal activity but functions predominantly from a short range, whereas a second domain is required to promote long-range activity. We show further that these two domains encode activities that are highly integrated and that the second domain is crucial in promoting the chromosomal conformational changes correlated with gene activity. During limb bud development, these activities encoded by the ZRS are interpreted differently by the fore limbs and the hind limbs; in the absence of the second domain there is no Shh activity in the fore limb, and in the hind limb low levels of Shh lead to a variant digit pattern ranging from two to four digits. Hence, in the embryo, the second domain stabilises the developmental programme providing a buffer for SHH morphogen activity and this ensures that five digits form in both sets of limbs.

Subtle changes in motif positioning cause tissue-specific effects on robustness of an enhancer's activity

Deciphering the specific contribution of individual motifs within cis-regulatory modules (CRMs) is crucial to understanding how gene expression is regulated and how this process is affected by sequence variation. But despite vast improvements in the ability to identify where transcription factors (TFs) bind throughout the genome, we are limited in our ability to relate information on motif occupancy to function from sequence alone. Here, we engineered 63 synthetic CRMs to systematically assess the relationship between variation in the content and spacing of motifs within CRMs to CRM activity during development using Drosophila transgenic embryos. In over half the cases, very simple elements containing only one or two types of TF binding motifs were capable of driving specific spatio-temporal patterns during development. Different motif organizations provide different degrees of robustness to enhancer activity, ranging from binary on-off responses to more subtle effects including embryo-to-embryo and within-embryo variation. By quantifying the effects of subtle changes in motif organization, we were able to model biophysical rules that explain CRM behavior and may contribute to the spatial positioning of CRM activity in vivo. For the same enhancer, the effects of small differences in motif positions varied in developmentally related tissues, suggesting that gene expression may be more susceptible to sequence variation in one tissue compared to another. This result has important implications for human eQTL studies in which many associated mutations are found in cis-regulatory regions, though the mechanism for how they affect tissue-specific gene expression is often not understood.

3D chromatin interactions organize Yan chromatin occupancy and repression at the even-skipped locus

Long-range integration of transcriptional inputs is critical for gene expression, yet the mechanisms are poorly understood. Here, Webber et al. investigate how cis-regulatory elements at the even-skipped (eve) locus coordinate heart-specific gene expression. Targeted deletions of regions bound by the repressor Yan define two novel enhancers that contribute repressive inputs to stabilize tissue-specific output from a third enhancer. This study provides a new paradigm for chromatin-level integration of repressive inputs with specific patterning information to achieve robust gene expression.

Long-range integration of transcriptional inputs is critical for gene expression, yet the mechanisms remain poorly understood. We investigated the molecular determinants that confer fidelity to expression of the heart identity gene even-skipped (eve). Targeted deletion of regions bound by the repressor Yan defined two novel enhancers that contribute repressive inputs to stabilize tissue-specific output from a third enhancer. Deletion of any individual enhancer reduced Yan occupancy at the other elements, impacting eve expression, cell fate specification, and cardiac function. These long-range interactions may be stabilized by three-dimensional chromatin contacts that we detected between the elements. Our work provides a new paradigm for chromatin-level integration of general repressive inputs with specific patterning information to achieve robust gene expression.

Assessing constraints on the path of regulatory sequence evolution

Structural and functional constraints are known to play a major role in restricting the path of evolution of protein activities. However, constraints acting on evolving transcriptional regulatory sequences, e.g. enhancers, are largely unknown. Recently, we elucidated how a novel expression pattern of the Neprilysin-1 (Nep1) gene in the optic lobe of Drosophila santomea evolved via co-option of existing enhancer activities. Drosophila santomea, which has diverged from Drosophila yakuba by approximately 400 000 years has accumulated four fixed mutations that each contribute to the full activity of this enhancer. Recreating and testing the optic lobe enhancer of the ancestor of D. santomea and D. yakuba revealed that the strong D. santomea enhancer activity evolved from a weak ancestral activity. Because each mutation on the path from the D. yakuba/santomea ancestor to modern-day D. santomea contributes to the newly derived optic lobe enhancer activity, we sought here to use this system to study the path of evolution of enhancer sequences. We inferred likely paths of evolution of this enhancer by observing the transcriptional output of all possible intermediate steps between the ancestral D. yakuba/santomea enhancer and the modern D. santomea enhancer. Many possible paths had epistatic and cooperative effects. Furthermore, we found that several paths significantly increased ectopic transcriptional activity or affected existing enhancer activities from which the novel activity was co-opted. We suggest that these attributes highlight constraints that guide the path of evolution of enhancers.

TALE-mediated modulation of transcriptional enhancers in vivo

We tested whether Transcription Activator-Like Effectors (TALEs) can mediate repression and activation of endogenous enhancers in the Drosophila genome. TALE-repressors (TALERs) targeting each of the five even-skipped (eve) “stripe” enhancers generated repression specifically of the focal stripes. TALE-activators (TALEAs) targeting the eve promoter or eve enhancers caused increased expression primarily in cells normally activated by the promoter or targeted enhancer, respectfully. The phenotypic effects of TALER and TALEA expression in larvae and adults are consistent with the observed modulations of eve expression. In these assays, the Hairy repression domain did not exhibit previously described long-range transcriptional repression activity. The precise effects of the TALEAs support the view that repression acts in a dominant fashion on transcriptional activators and that the activity state of an enhancer influences TALE binding or the ability of VP16 to enhance transcription. TALEs thus provide a novel tool for detection and functional modulation of transcriptional enhancers in their native genomic context.

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.

Accurate measurements of dynamics and reproducibility in small genetic networks

Quantification of gene expression has become a central tool for understanding genetic networks. In many systems, the only viable way to measure protein levels is by immunofluorescence, which is notorious for its limited accuracy. Using the early Drosophila embryo as an example, we show that careful identification and control of experimental error allows for highly accurate gene expression measurements. We generated antibodies in different host species, allowing for simultaneous staining of four Drosophila gap genes in individual embryos. Careful error analysis of hundreds of expression profiles reveals that less than ∼20% of the observed embryo-to-embryo fluctuations stem from experimental error. These measurements make it possible to extract not only very accurate mean gene expression profiles but also their naturally occurring fluctuations of biological origin and corresponding cross-correlations. We use this analysis to extract gap gene profile dynamics with ∼1 min accuracy. The combination of these new measurements and analysis techniques reveals a twofold increase in profile reproducibility owing to a collective network dynamics that relays positional accuracy from the maternal gradients to the pair-rule genes.

The relationship between long-range chromatin occupancy and polymerization of the Drosophila ETS family transcriptional repressor Yan

ETS family transcription factors are evolutionarily conserved downstream effectors of Ras/MAPK signaling with critical roles in development and cancer. In Drosophila, the ETS repressor Yan regulates cell proliferation and differentiation in a variety of tissues; however, the mechanisms of Yan-mediated repression are not well understood and only a few direct target genes have been identified. Yan, like its human ortholog TEL1, self-associates through an N-terminal sterile α-motif (SAM), leading to speculation that Yan/TEL1 polymers may spread along chromatin to form large repressive domains. To test this hypothesis, we created a monomeric form of Yan by recombineering a point mutation that blocks SAM-mediated self-association into the yan genomic locus and compared its genome-wide chromatin occupancy profile to that of endogenous wild-type Yan. Consistent with the spreading model predictions, wild-type Yan-bound regions span multiple kilobases. Extended occupancy patterns appear most prominent at genes encoding crucial developmental regulators and signaling molecules and are highly conserved between Drosophila melanogaster and D. virilis, suggesting functional relevance. Surprisingly, although occupancy is reduced, the Yan monomer still makes extensive multikilobase contacts with chromatin, with an overall pattern similar to that of wild-type Yan. Despite its near-normal chromatin recruitment, the repressive function of the Yan monomer is significantly impaired, as evidenced by elevated target gene expression and failure to rescue a yan null mutation. Together our data argue that SAM-mediated polymerization contributes to the functional output of the active Yan repressive complexes that assemble across extended stretches of chromatin, but does not directly mediate recruitment to DNA or chromatin spreading.

Deciphering the transcriptional cis-regulatory code

Information about developmental gene expression resides in defined regulatory elements, called enhancers, in the non-coding part of the genome. Although cells reliably utilize enhancers to orchestrate gene expression, a cis-regulatory code that would allow their interpretation has remained one of the greatest challenges of modern biology. In this review, we summarize studies from the past three decades that describe progress towards revealing the properties of enhancers and discuss how recent approaches are providing unprecedented insights into regulatory elements in animal genomes. Over the next years, we believe that the functional characterization of regulatory sequences in entire genomes, combined with recent computational methods, will provide a comprehensive view of genomic regulatory elements and their building blocks and will enable researchers to begin to understand the sequence basis of the cis-regulatory code.

Disentangling the many layers of eukaryotic transcriptional regulation

Regulation of gene expression in eukaryotes is an extremely complex process. In this review, we break down several critical steps, emphasizing new data and techniques that have expanded current gene regulatory models. We begin at the level of DNA sequence where cis-regulatory modules (CRMs) provide important regulatory information in the form of transcription factor (TF) binding sites. In this respect, CRMs function as instructional platforms for the assembly of gene regulatory complexes. We discuss multiple mechanisms controlling complex assembly, including cooperative DNA binding, combinatorial codes, and CRM architecture. The second section of this review places CRM assembly in the context of nucleosomes and condensed chromatin. We discuss how DNA accessibility and histone modifications contribute to TF function. Lastly, new advances in chromosomal mapping techniques have provided increased understanding of intra- and interchromosomal interactions. We discuss how these topological maps influence gene regulatory models.