Phenotypic robustness conferred by apparently redundant transcriptional enhancers

Regulatory genome annotation of 33 insect species

Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules-e.g., enhancers and silencers-that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

Evolutionary Innovations in Conserved Regulatory Elements Associate With Developmental Genes in Mammals

Transcriptional enhancers orchestrate cell type- and time point-specific gene expression programs. Genetic variation within enhancer sequences is an important contributor to phenotypic variation including evolutionary adaptations and human disease. Certain genes and pathways may be more prone to regulatory evolution than others, with different patterns across diverse organisms, but whether such patterns exist has not been investigated at a sufficient scale. To address this question, we identified signatures of accelerated sequence evolution in conserved enhancer elements throughout the mammalian phylogeny at an unprecedented scale. While different genes and pathways were enriched for regulatory evolution in different parts of the tree, we found a striking overall pattern of pleiotropic genes involved in gene regulatory and developmental processes being enriched for accelerated enhancer evolution. These genes were connected to more enhancers than other genes, which was the basis for having an increased amount of sequence acceleration over all their enhancers combined. We provide evidence that sequence acceleration is associated with turnover of regulatory function. Detailed study of one acceleration event in an enhancer of HES1 revealed that sequence evolution led to a new activity domain in the developing limb that emerged concurrently with the evolution of digit reduction in hoofed mammals. Our results provide evidence that enhancer evolution has been a frequent contributor to regulatory innovation at conserved developmental signaling genes in mammals.

Neuronal enhancers fine-tune adaptive circuit plasticity

Neuronal activity-regulated gene expression plays a crucial role in sculpting neural circuits that underpin adaptive brain function. Transcriptional enhancers are now recognized as key components of gene regulation that orchestrate spatiotemporally precise patterns of gene transcription. We propose that the dynamics of enhancer activation uniquely position these genomic elements to finely tune activity-dependent cellular plasticity. Enhancer specificity and modularity can be exploited to gain selective genetic access to specific cell states, and the precise modulation of target gene expression within restricted cellular contexts enabled by targeted enhancer manipulation allows for fine-grained evaluation of gene function. Mounting evidence also suggests that enduring stimulus-induced changes in enhancer states can modify target gene activation upon restimulation, thereby contributing to a form of cell-wide metaplasticity. We advocate for focused exploration of activity-dependent enhancer function to gain new insight into the mechanisms underlying brain plasticity and cognitive dysfunction.

Super-silencer perturbation by EZH2 and REST inhibition leads to large loss of chromatin interactions and reduction in cancer growth

Human silencers have been shown to regulate developmental gene expression. However, the functional importance of human silencers needs to be elucidated, such as whether they can form 'super-silencers' and whether they are linked to cancer progression. Here, we show two silencer components of the FGF18 gene can cooperate through compensatory chromatin interactions to form a super-silencer. Double knockout of two silencers exhibited synergistic upregulation of FGF18 expression and changes in cell identity. To perturb the super-silencers, we applied combinational treatment of an enhancer of zeste homolog 2 inhibitor GSK343, and a repressor element 1-silencing transcription factor inhibitor, X5050 ('GR'). Interestingly, GR led to severe loss of topologically associated domains and loops, which were associated with reduced CTCF and TOP2A mRNA levels. Moreover, GR synergistically upregulated super-silencer-controlled genes related to cell cycle, apoptosis and DNA damage, leading to anticancer effects in vivo. Overall, our data demonstrated a super-silencer example and showed that GR can disrupt super-silencers, potentially leading to cancer ablation.

Interpreting cis-regulatory interactions from large-scale deep neural networks

The rise of large-scale, sequence-based deep neural networks (DNNs) for predicting gene expression has introduced challenges in their evaluation and interpretation. Current evaluations align DNN predictions with orthogonal experimental data, providing insights into generalization but offering limited insights into their decision-making process. Existing model explainability tools focus mainly on motif analysis, which becomes complex when interpreting longer sequences. Here we present cis-regulatory element model explanations (CREME), an in silico perturbation toolkit that interprets the rules of gene regulation learned by a genomic DNN. Applying CREME to Enformer, a state-of-the-art DNN, we identify cis-regulatory elements that enhance or silence gene expression and characterize their complex interactions. CREME can provide interpretations across multiple scales of genomic organization, from cis-regulatory elements to fine-mapped functional sequence elements within them, offering high-resolution insights into the regulatory architecture of the genome. CREME provides a powerful toolkit for translating the predictions of genomic DNNs into mechanistic insights of gene regulation.

Morphogens in the evolution of size, shape and patterning

Much of the striking diversity of life on Earth has arisen from variations in the way that the same molecules and networks operate during development to shape and pattern tissues and organs into different morphologies. However, we still understand very little about the potential for diversification exhibited by different, highly conserved mechanisms during evolution, or, conversely, the constraints that they place on evolution. With the aim of steering the field in new directions, we focus on morphogen-mediated patterning and growth as a case study to demonstrate how conserved developmental mechanisms can adapt during evolution to drive morphological diversification and optimise functionality, and to illustrate how evolution algorithms and computational tools can be used alongside experiments to provide insights into how these conserved mechanisms can evolve. We first introduce key conserved properties of morphogen-driven patterning mechanisms, before summarising comparative studies that exemplify how changes in the spatiotemporal expression and signalling levels of morphogens impact the diversification of organ size, shape and patterning in nature. Finally, we detail how theoretical frameworks can be used in conjunction with experiments to probe the role of morphogen-driven patterning mechanisms in evolution. We conclude that morphogen-mediated patterning is an excellent model system and offers a generally applicable framework to investigate the evolution of developmental mechanisms.

Enhancer-promoter hubs organize transcriptional networks promoting oncogenesis and drug resistance

Recent advances in high-resolution mapping of spatial interactions among regulatory elements support the existence of complex topological assemblies of enhancers and promoters known as enhancer-promoter hubs or cliques. Yet, organization principles of these multi-interacting enhancer-promoter hubs and their potential role in regulating gene expression in cancer remain unclear. Here, we systematically identify enhancer-promoter hubs in breast cancer, lymphoma, and leukemia. We find that highly interacting enhancer-promoter hubs form at key oncogenes and lineage-associated transcription factors potentially promoting oncogenesis of these diverse cancer types. Genomic and optical mapping of interactions among enhancer and promoter elements further show that topological alterations in hubs coincide with transcriptional changes underlying acquired resistance to targeted therapy in T cell leukemia and B cell lymphoma. Together, our findings suggest that enhancer-promoter hubs are dynamic and heterogeneous topological assemblies with the potential to control gene expression circuits promoting oncogenesis and drug resistance.

Beyond modular enhancers: new questions in cis-regulatory evolution

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

The evolution of developmental biology through conceptual and technological revolutions

Developmental biology-the study of the processes by which cells, tissues, and organisms develop and change over time-has entered a new golden age. After the molecular genetics revolution in the 80s and 90s and the diversification of the field in the early 21st century, we have entered a phase when powerful technologies provide new approaches and open unexplored avenues. Progress in the field has been accelerated by advances in genomics, imaging, engineering, and computational biology and by emerging model systems ranging from tardigrades to organoids. We summarize how revolutionary technologies have led to remarkable progress in understanding animal development. We describe how classic questions in gene regulation, pattern formation, morphogenesis, organogenesis, and stem cell biology are being revisited. We discuss the connections of development with evolution, self-organization, metabolism, time, and ecology. We speculate how developmental biology might evolve in an era of synthetic biology, artificial intelligence, and human engineering.

Enhancer-promoter hubs organize transcriptional networks promoting oncogenesis and drug resistance

Recent advances in high-resolution mapping of spatial interactions among regulatory elements support the existence of complex topological assemblies of enhancers and promoters known as enhancer-promoter hubs or cliques. Yet, organization principles of these multi-interacting enhancer-promoter hubs and their potential role in regulating gene expression in cancer remains unclear. Here, we systematically identified enhancer-promoter hubs in breast cancer, lymphoma, and leukemia. We found that highly interacting enhancer-promoter hubs form at key oncogenes and lineage-associated transcription factors potentially promoting oncogenesis of these diverse cancer types. Genomic and optical mapping of interactions among enhancer and promoter elements further showed that topological alterations in hubs coincide with transcriptional changes underlying acquired resistance to targeted therapy in T cell leukemia and B cell lymphoma. Together, our findings suggest that enhancer-promoter hubs are dynamic and heterogeneous topological assemblies with the potential to control gene expression circuits promoting oncogenesis and drug resistance.

Accelerated evolution in the human lineage led to gain and loss of transcriptional enhancers in the RBFOX1 locus

A long-standing goal of evolutionary biology is to decode how changes in gene regulatory networks contribute to human-specific traits. Human accelerated regions (HARs) are prime candidates for driving gene regulatory modifications in human development. The RBFOX1 locus is densely populated with HARs, providing a set of potential regulatory elements that could have changed its expression in the human lineage. Here, we examined the role of RBFOX1-HARs using transgenic zebrafish reporter assays and identified 15 transcriptional enhancers that are active in the developing nervous system, 9 of which displayed differential activity between the human and chimpanzee sequences. The engineered loss of two selected RBFOX1-HARs in knockout mouse models modified Rbfox1 expression at specific developmental stages and tissues in the brain, influencing the expression and splicing of a high number of Rbfox1 target genes. Our results provided insight into the spatial and temporal changes in gene expression driven by RBFOX1-HARs.

Genomic context sensitizes regulatory elements to genetic disruption

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

The single-cell chromatin landscape in gonadal cell lineage specification

Gonad development includes sex determination and divergent maturation of the testes and ovaries. Recent advances in measuring gene expression in single cells are providing new insights into this complex process. However, the underlying epigenetic regulatory mechanisms remain unclear. Here, we profiled chromatin accessibility in mouse gonadal cells of both sexes from embryonic day 11.5 to 14.5 using single-cell assay for transposase accessible chromatin by sequencing (scATAC-seq). Our results showed that individual cell types can be inferred by the chromatin landscape, and that cells can be temporally ordered along developmental trajectories. Integrative analysis of transcriptomic and chromatin-accessibility maps identified multiple putative regulatory elements proximal to key gonadal genes Nr5a1, Sox9 and Wt1. We also uncover cell type-specific regulatory factors underlying cell type specification. Overall, our results provide a better understanding of the epigenetic landscape associated with the progressive restriction of cell fates in the gonad.

Gene regulatory network co-option is sufficient to induce a morphological novelty in Drosophila

Identifying the molecular origins by which new morphological structures evolve is one of the long standing problems in evolutionary biology. To date, vanishingly few examples provide a compelling account of how new morphologies were initially formed, thereby limiting our understanding of how diverse forms of life derived their complex features. Here, we provide evidence that the large projections on the Drosophila eugracilis phallus that are implicated in sexual conflict have evolved through co-option of the trichome genetic network. These unicellular apical projections on the phallus postgonal sheath are reminiscent of trichomes that cover the Drosophila body but are up to 20-fold larger in size. During their development, they express the transcription factor Shavenbaby, the master regulator of the trichome network. Consistent with the co-option of the Shavenbaby network during the evolution of the D. eugracilis projections, somatic mosaic CRISPR/Cas9 mutagenesis shows that shavenbaby is necessary for their proper length. Moreover, mis-expression of Shavenbaby in the sheath of D. melanogaster , a naïve species that lacks these extensions, is sufficient to induce small trichomes. These induced extensions rely on a genetic network that is shared to a large extent with the D. eugracilis projections, indicating its co-option but also some genetic rewiring. Thus, by leveraging a genetically tractable evolutionarily novelty, our work shows that the trichome-forming network is flexible enough that it can be co-opted in a new context, and subsequently refined to produce unique apical projections that are barely recognizable compared to their simpler ancestral beginnings.

Interpreting Cis-Regulatory Interactions from Large-Scale Deep Neural Networks for Genomics

The rise of large-scale, sequence-based deep neural networks (DNNs) for predicting gene expression has introduced challenges in their evaluation and interpretation. Current evaluations align DNN predictions with experimental perturbation assays, which provides insights into the generalization capabilities within the studied loci but offers a limited perspective of what drives their predictions. Moreover, existing model explainability tools focus mainly on motif analysis, which becomes complex when interpreting longer sequences. Here we introduce CREME, an in silico perturbation toolkit that interrogates large-scale DNNs to uncover rules of gene regulation that it learns. Using CREME, we investigate Enformer, a prominent DNN in gene expression prediction, revealing cis-regulatory elements (CREs) that directly enhance or silence target genes. We explore the intricate complexity of higher-order CRE interactions, the relationship between CRE distance from transcription start sites on gene expression, as well as the biochemical features of enhancers and silencers learned by Enformer. Moreover, we demonstrate the flexibility of CREME to efficiently uncover a higher-resolution view of functional sequence elements within CREs. This work demonstrates how CREME can be employed to translate the powerful predictions of large-scale DNNs to study open questions in gene regulation.

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

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

Deciphering the origin of developmental stability: The role of intracellular expression variability in evolutionary conservation

Progress in evolutionary developmental biology (evo-devo) has deepened our understanding of how intrinsic properties of embryogenesis, along with natural selection and population genetics, shape phenotypic diversity. A focal point of recent empirical and theoretical research is the idea that highly developmentally stable phenotypes are more conserved in evolution. Previously, we demonstrated that in Japanese medaka (Oryzias latipes), embryonic stages and genes with high stability, estimated through whole-embryo RNA-seq, are highly conserved in subsequent generations. However, the precise origin of the stability of gene expression levels evaluated at the whole-embryo level remained unclear. Such stability could be attributed to two distinct sources: stable intracellular expression levels or spatially stable expression patterns. Here we demonstrate that stability observed in whole-embryo RNA-seq can be attributed to stability at the cellular level (low variability in gene expression at the cellular levels). We quantified the intercellular variations in expression levels and spatial gene expression patterns for seven key genes involved in patterning dorsoventral and rostrocaudal regions during early development in medaka. We evaluated intracellular variability by counting transcripts and found its significant correlation with variation observed in whole-embryo RNA-seq data. Conversely, variation in spatial gene expression patterns, assessed through intraindividual left-right asymmetry, showed no correlation. Given the previously reported correlation between stability and conservation of expression levels throughout embryogenesis, our findings suggest a potential general trend: the stability or instability of developmental systems-and the consequent evolutionary diversity-may be primarily anchored in intrinsic fundamental elements such as the variability of intracellular states.

Extreme restructuring of cis-regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis-regulatory sequences. It remains unclear how such drastic cis-regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis-regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ~125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 (CLV3) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis-regulatory sequences. Our results demonstrate remarkable malleability in cis-regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis-regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis-regulation to effectively engineer trait variation from conserved productivity genes in crops.

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.

Deep molecular learning of transcriptional control of a synthetic CRE enhancer and its variants

Massively parallel reporter assay measures transcriptional activities of various cis-regulatory modules (CRMs) in a single experiment. We developed a thermodynamic computational model framework that calculates quantitative levels of gene expression directly from regulatory DNA sequences. Using the framework, we investigated the molecular mechanisms of cis-regulatory mutations of a synthetic enhancer that cause abnormal gene expression. We found that, in a human cell line, competitive binding between family transcription factors (TFs) with slightly different binding preferences significantly increases the accuracy of recapitulating the transcriptional effects of thousands of single- or multi-mutations. We also discovered that even if various harmful mutations occurred in an activator binding site, CRM could stably maintain or even increase gene expression through a certain form of competitive binding between family TFs. These findings enhance understanding the effect of SNPs and indels on CRMs and would help building robust custom-designed CRMs for biologics production and gene therapy.

Extreme restructuring of cis -regulatory regions controlling a deeply conserved plant stem cell regulator

A striking paradox is that genes with conserved protein sequence, function and expression pattern over deep time often exhibit extremely divergent cis -regulatory sequences. It remains unclear how such drastic cis -regulatory evolution across species allows preservation of gene function, and to what extent these differences influence how cis- regulatory variation arising within species impacts phenotypic change. Here, we investigated these questions using a plant stem cell regulator conserved in expression pattern and function over ∼125 million years. Using in-vivo genome editing in two distantly related models, Arabidopsis thaliana (Arabidopsis) and Solanum lycopersicum (tomato), we generated over 70 deletion alleles in the upstream and downstream regions of the stem cell repressor gene CLAVATA3 ( CLV3 ) and compared their individual and combined effects on a shared phenotype, the number of carpels that make fruits. We found that sequences upstream of tomato CLV3 are highly sensitive to even small perturbations compared to its downstream region. In contrast, Arabidopsis CLV3 function is tolerant to severe disruptions both upstream and downstream of the coding sequence. Combining upstream and downstream deletions also revealed a different regulatory outcome. Whereas phenotypic enhancement from adding downstream mutations was predominantly weak and additive in tomato, mutating both regions of Arabidopsis CLV3 caused substantial and synergistic effects, demonstrating distinct distribution and redundancy of functional cis -regulatory sequences. Our results demonstrate remarkable malleability in cis -regulatory structural organization of a deeply conserved plant stem cell regulator and suggest that major reconfiguration of cis -regulatory sequence space is a common yet cryptic evolutionary force altering genotype-to-phenotype relationships from regulatory variation in conserved genes. Finally, our findings underscore the need for lineage-specific dissection of the spatial architecture of cis -regulation to effectively engineer trait variation from conserved productivity genes in crops.

Author summary

We investigated the evolution of cis -regulatory elements (CREs) and their interactions in the regulation of a plant stem cell regulator gene, CLAVATA3 (CLV3) , in Arabidopsis and tomato. Despite diverging ∼125 million years ago, the function and expression of CLV3 is conserved in these species; however, cis -regulatory sequences upstream and downstream have drastically diverged, preventing identification of conserved non-coding sequences between them. We used CRISPR-Cas9 to engineer dozens of mutations within the cis -regulatory regions of Arabidopsis and tomato CLV3. In tomato, our results show that tomato CLV3 function primarily relies on interactions among CREs in the 5' non-coding region, unlike Arabidopsis CLV3 , which depends on a more balanced distribution of functional CREs between the 5' and 3' regions. Therefore, despite a high degree of functional conservation, our study demonstrates divergent regulatory strategies between two distantly related CLV3 orthologs, with substantial alterations in regulatory sequences, their spatial arrangement, and their relative effects on CLV3 regulation. These results suggest that regulatory regions are not only extremely robust to mutagenesis, but also that the sequences underlying this robustness can be lineage-specific for conserved genes, due to the complex and often redundant interactions among CREs that ensure proper gene function amidst large-scale sequence turnover.

Functional analysis of the Drosophila eve locus in response to non-canonical combinations of gap gene expression levels

Transcription factor combinations play a key role in shaping cellular identity. However, the precise relationship between specific combinations and downstream effects remains elusive. Here, we investigate this relationship within the context of the Drosophila eve locus, which is controlled by gap genes. We measure spatiotemporal levels of four gap genes in heterozygous and homozygous gap mutant embryos and correlate them with the striped eve activity pattern. Although changes in gap gene expression extend beyond the manipulated gene, the spatial patterns of Eve expression closely mirror canonical activation levels in wild type. Interestingly, some combinations deviate from the wild-type repertoire but still drive eve activation. Although in homozygous mutants some Eve stripes exhibit partial penetrance, stripes consistently emerge at reproducible positions, even with varying gap gene levels. Our findings suggest a robust molecular canalization of cell fates in gap mutants and provide insights into the regulatory constraints governing multi-enhancer gene loci.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

Enduring questions in regenerative biology and the search for answers

The potential for basic research to uncover the inner workings of regenerative processes and produce meaningful medical therapies has inspired scientists, clinicians, and patients for hundreds of years. Decades of studies using a handful of highly regenerative model organisms have significantly advanced our knowledge of key cell types and molecular pathways involved in regeneration. However, many questions remain about how regenerative processes unfold in regeneration-competent species, how they are curtailed in non-regenerative organisms, and how they might be induced (or restored) in humans. Recent technological advances in genomics, molecular biology, computer science, bioengineering, and stem cell research hold promise to collectively provide new experimental evidence for how different organisms accomplish the process of regeneration. In theory, this new evidence should inform the design of new clinical approaches for regenerative medicine. A deeper understanding of how tissues and organs regenerate will also undoubtedly impact many adjacent scientific fields. To best apply and adapt these new technologies in ways that break long-standing barriers and answer critical questions about regeneration, we must combine the deep knowledge of developmental and evolutionary biologists with the hard-earned expertise of scientists in mechanistic and technical fields. To this end, this perspective is based on conversations from a workshop we organized at the Banbury Center, during which a diverse cross-section of the regeneration research community and experts in various technologies discussed enduring questions in regenerative biology. Here, we share the questions this group identified as significant and unanswered, i.e., known unknowns. We also describe the obstacles limiting our progress in answering these questions and how expanding the number and diversity of organisms used in regeneration research is essential for deepening our understanding of regenerative capacity. Finally, we propose that investigating these problems collaboratively across a diverse network of researchers has the potential to advance our field and produce unexpected insights into important questions in related areas of biology and medicine.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal discs; in contrast, small IDE-reporter transgenes are expressed mainly in the anterior compartment. We found that En binds to the IDEs and suggest that it may directly repress IDE function and modulate En expression levels. We identified two en IDEs, O and S. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from an endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data suggest that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic memory set by the embryonic enhancers. This epigenetic regulation is similar to that of the Ultrabithorax IDEs and we suggest that the activity of late-acting enhancers in other genes may be similarly regulated.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

Deep multiple-instance learning accurately predicts gene haploinsufficiency and deletion pathogenicity

Copy number losses (deletions) are a major contributor to the etiology of severe genetic disorders. Although haploinsufficient genes play a critical role in deletion pathogenicity, current methods for deletion pathogenicity prediction fail to integrate multiple lines of evidence for haploinsufficiency at the gene level, limiting their power to pinpoint deleterious deletions associated with genetic disorders. Here we introduce DosaCNV, a deep multiple-instance learning framework that, for the first time, models deletion pathogenicity jointly with gene haploinsufficiency. By integrating over 30 gene-level features potentially predictive of haploinsufficiency, DosaCNV shows unmatched performance in prioritizing pathogenic deletions associated with a broad spectrum of genetic disorders. Furthermore, DosaCNV outperforms existing methods in predicting gene haploinsufficiency even though it is not trained on known haploinsufficient genes. Finally, DosaCNV leverages a state-of-the-art technique to quantify the contributions of individual gene-level features to haploinsufficiency, allowing for human-understandable explanations of model predictions. Altogether, DosaCNV is a powerful computational tool for both fundamental and translational research.

The interaction between enhancer variants and environmental factors as an overlooked aetiological paradigm in human complex disease

The interactions between genetic and environmental risk factors contribute to the aetiology of complex human diseases. Genome-wide association studies (GWAS) have revealed that most of the genetic variants associated with complex diseases are located in the non-coding part of the genome, preferentially within enhancers. Enhancers are distal cis-regulatory elements composed of clusters of transcription factors binding sites that positively regulate the expression of their target genes. The generation of genome-wide maps for histone marks (e.g., H3K27ac), chromatin accessibility and transcription factor and coactivator (e.g., p300) binding profiles have enabled the identification of enhancers across many human cell types and tissues. Nonetheless, the functional and pathological consequences of the majority of disease-associated genetic variants located within enhancers seem to be rather minor under normal conditions, thus questioning their medical relevance. Here we propose that, due to the prevalence of enhancer redundancy, the pathological effects of many disease-associated non-coding genetic variants might be preferentially (or even only) manifested under environmental stress.

Regulatory changes associated with the head to trunk developmental transition

BACKGROUND

Development of vertebrate embryos is characterized by early formation of the anterior tissues followed by the sequential extension of the axis at their posterior end to build the trunk and tail structures, first by the activity of the primitive streak and then of the tail bud. Embryological, molecular and genetic data indicate that head and trunk development are significantly different, suggesting that the transition into the trunk formation stage involves major changes in regulatory gene networks.

RESULTS

We explored those regulatory changes by generating differential interaction networks and chromatin accessibility profiles from the posterior epiblast region of mouse embryos at embryonic day (E)7.5 and E8.5. We observed changes in various cell processes, including several signaling pathways, ubiquitination machinery, ion dynamics and metabolic processes involving lipids that could contribute to the functional switch in the progenitor region of the embryo. We further explored the functional impact of changes observed in Wnt signaling associated processes, revealing a switch in the functional relevance of Wnt molecule palmitoleoylation, essential during gastrulation but becoming differentially required for the control of axial extension and progenitor differentiation processes during trunk formation. We also found substantial changes in chromatin accessibility at the two developmental stages, mostly mapping to intergenic regions and presenting differential footprinting profiles to several key transcription factors, indicating a significant switch in the regulatory elements controlling head or trunk development. Those chromatin changes are largely independent of retinoic acid, despite the key role of this factor in the transition to trunk development. We also tested the functional relevance of potential enhancers identified in the accessibility assays that reproduced the expression profiles of genes involved in the transition. Deletion of these regions by genome editing had limited effect on the expression of those genes, suggesting the existence of redundant enhancers that guarantee robust expression patterns.

CONCLUSIONS

This work provides a global view of the regulatory changes controlling the switch into the axial extension phase of vertebrate embryonic development. It also revealed mechanisms by which the cellular context influences the activity of regulatory factors, channeling them to implement one of several possible biological outputs.

Dynamic network-guided CRISPRi screen identifies CTCF-loop-constrained nonlinear enhancer gene regulatory activity during cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those contributing to complex diseases, have weak effects on gene expression. Our gene regulatory network modeling identified that nonlinear enhancer gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Using human embryonic stem cell definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. We discovered a comprehensive set of enhancers for each of the core endoderm-specifying transcription factors. Many enhancers had strong effects mid-transition but weak effects post-transition, consistent with the nonlinear temporal responses to enhancer perturbation predicted by the modeling. Integrating three-dimensional genomic information, we were able to develop a CTCF-loop-constrained Interaction Activity model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and systematic enhancer discovery in both normal and pathological cell state transitions.

Precisely timed regulation of enhancer activity defines the binary expression pattern of Fushi tarazu in the Drosophila embryo

The genes that drive development each typically have many different enhancers. Properly coordinating the action of these different enhancers is crucial to correctly specifying cell-fate decisions, yet it remains poorly understood how their activity is choregraphed in time. To shed light on this question, we used recently developed single-cell live imaging tools to dissect the regulation of Fushi tarazu (Ftz) in Drosophila melanogaster embryos. Ftz is a transcription factor that is expressed in asymmetric stripes by two distinct enhancers: autoregulatory and zebra. The anterior edge of each stripe needs to be sharply defined to specify essential lineage boundaries. Here, we tracked how boundary cells commit to either a high-Ftz or low-Ftz fate by measuring Ftz protein traces in real time and simultaneously quantifying transcription from the endogenous locus and individual enhancers. This revealed that the autoregulatory enhancer does not establish this fate choice. Instead, it perpetuates the decision defined by zebra. This is contrary to the prevailing view that autoregulation drives the fate decision by causing bi-stable Ftz expression. Furthermore, we showed that the autoregulatory enhancer is not activated based on a Ftz-concentration threshold but through a timing-based mechanism. We hypothesize that this is regulated by several ubiquitously expressed pioneer-like transcription factors, which have recently been shown to act as timers in the embryo. Our work provides new insight into how precisely timed enhancer activity can directly regulate the dynamics of gene regulatory networks, which may be a general mechanism for ensuring that embryogenesis runs like clockwork.

Speciation and development

Understanding general principles about the origin of species remains one of the foundational challenges in evolutionary biology. The genomic divergence between groups of individuals can spawn hybrid inviability and hybrid sterility, which presents a tantalizing developmental problem. Divergent developmental programs may yield either conserved or divergent phenotypes relative to ancestral traits, both of which can be responsible for reproductive isolation during the speciation process. The genetic mechanisms of developmental evolution involve cis- and trans-acting gene regulatory change, protein-protein interactions, genetic network structures, dosage, and epigenetic regulation, all of which also have roots in population genetic and molecular evolutionary processes. Toward the goal of demystifying Darwin's "mystery of mysteries," this review integrates microevolutionary concepts of genetic change with principles of organismal development, establishing explicit links between population genetic process and developmental mechanisms in the production of macroevolutionary pattern. This integration aims to establish a more unified view of speciation that binds process and mechanism.

The activity of engrailed imaginal disc enhancers is modulated epigenetically by chromatin and autoregulation

engrailed (en) encodes a homeodomain transcription factor crucial for the proper development of Drosophila embryos and adults. Like many developmental transcription factors, en expression is regulated by many enhancers, some of overlapping function, that drive expression in spatially and temporally restricted patterns. The en embryonic enhancers are located in discrete DNA fragments that can function correctly in small reporter transgenes. In contrast, the en imaginal disc enhancers (IDEs) do not function correctly in small reporter transgenes. En is expressed in the posterior compartment of wing imaginal disks; small IDE-reporter transgenes are expressed in the anterior compartment, the opposite of what is expected. Our data show that the En protein binds to en IDEs, and we suggest that En directly represses IDE function. We identified two en IDEs, 'O' and 'S'. Deletion of either of these IDEs from a 79kb HA-en rescue transgene (HAen79) caused a loss-of-function en phenotype when the HAen79 transgene was the sole source of En. In contrast, flies with a deletion of the same IDEs from the endogenous en gene had no phenotype, suggesting a resiliency not seen in the HAen79 rescue transgene. Inserting a gypsy insulator in HAen79 between en regulatory DNA and flanking sequences strengthened the activity of HAen79, giving better function in both the ON and OFF transcriptional states. Altogether our data show that the en IDEs stimulate expression in the entire imaginal disc, and that the ON/OFF state is set by epigenetic regulators. Further, the endogenous locus imparts a stability to en function not seen even in a large transgene, reflecting the importance of both positive and negative epigenetic influences that act over relatively large distances in chromatin.

Esearch3D: propagating gene expression in chromatin networks to illuminate active enhancers

Most cell type-specific genes are regulated by the interaction of enhancers with their promoters. The identification of enhancers is not trivial as enhancers are diverse in their characteristics and dynamic in their interaction partners. We present Esearch3D, a new method that exploits network theory approaches to identify active enhancers. Our work is based on the fact that enhancers act as a source of regulatory information to increase the rate of transcription of their target genes and that the flow of this information is mediated by the folding of chromatin in the three-dimensional (3D) nuclear space between the enhancer and the target gene promoter. Esearch3D reverse engineers this flow of information to calculate the likelihood of enhancer activity in intergenic regions by propagating the transcription levels of genes across 3D genome networks. Regions predicted to have high enhancer activity are shown to be enriched in annotations indicative of enhancer activity. These include: enhancer-associated histone marks, bidirectional CAGE-seq, STARR-seq, P300, RNA polymerase II and expression quantitative trait loci (eQTLs). Esearch3D leverages the relationship between chromatin architecture and transcription, allowing the prediction of active enhancers and an understanding of the complex underpinnings of regulatory networks. The method is available at: https://github.com/InfOmics/Esearch3D and https://doi.org/10.5281/zenodo.7737123.

Comparative three-dimensional genome architectures of adipose tissues provide insight into human-specific regulation of metabolic homeostasis

Elucidating the regulatory mechanisms of human adipose tissues (ATs) evolution is essential for understanding human-specific metabolic regulation, but the functional importance and evolutionary dynamics of three-dimensional (3D) genome organizations of ATs are not well defined. Here, we compared the 3D genome architectures of anatomically distinct ATs from humans and six representative mammalian models. We recognized evolutionarily conserved and human-specific chromatin conformation in ATs at multiple scales, including compartmentalization, topologically associating domain (TAD), and promoter-enhancer interactions (PEI), which have not been described previously. We found PEI are much more evolutionarily dynamic with respect to compartmentalization and topologically associating domain. Compared to conserved PEIs, human-specific PEIs are enriched for human-specific sequence, and the binding motifs of their potential mediators (transcription factors) are less conserved. Our data also demonstrated that genes involved in the evolutionary dynamics of chromatin organization have weaker transcriptional conservation than those associated with conserved chromatin organization. Furthermore, the genes involved in energy metabolism and the maintenance of metabolic homeostasis are enriched in human-specific chromatin organization, while housekeeping genes, health-related genes, and genetic variations are enriched in evolutionarily conserved compared to human-specific chromatin organization. Finally, we showed extensively divergent human-specific 3D genome organizations among one subcutaneous and three visceral ATs. Together, these findings provide a global overview of 3D genome architecture dynamics between ATs from human and mammalian models and new insights into understanding the regulatory evolution of human ATs.

Mutational scans reveal differential evolvability of Drosophila promoters and enhancers

Rapid enhancer and slow promoter evolution have been demonstrated through comparative genomics. However, it is not clear how this information is encoded genetically and if this can be used to place evolution in a predictive context. Part of the challenge is that our understanding of the potential for regulatory evolution is biased primarily toward natural variation or limited experimental perturbations. Here, to explore the evolutionary capacity of promoter variation, we surveyed an unbiased mutation library for three promoters in Drosophila melanogaster. We found that mutations in promoters had limited to no effect on spatial patterns of gene expression. Compared to developmental enhancers, promoters are more robust to mutations and have more access to mutations that can increase gene expression, suggesting that their low activity might be a result of selection. Consistent with these observations, increasing the promoter activity at the endogenous locus of shavenbaby led to increased transcription yet limited phenotypic changes. Taken together, developmental promoters may encode robust transcriptional outputs allowing evolvability through the integration of diverse developmental enhancers.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Shadow enhancers mediate trade-offs between transcriptional noise and fidelity

Enhancers are stretches of regulatory DNA that bind transcription factors (TFs) and regulate the expression of a target gene. Shadow enhancers are two or more enhancers that regulate the same target gene in space and time and are associated with most animal developmental genes. These multi-enhancer systems can drive more consistent transcription than single enhancer systems. Nevertheless, it remains unclear why shadow enhancer TF binding sites are distributed across multiple enhancers rather than within a single large enhancer. Here, we use a computational approach to study systems with varying numbers of TF binding sites and enhancers. We employ chemical reaction networks with stochastic dynamics to determine the trends in transcriptional noise and fidelity, two key performance objectives of enhancers. This reveals that while additive shadow enhancers do not differ in noise and fidelity from their single enhancer counterparts, sub- and superadditive shadow enhancers have noise and fidelity trade-offs not available to single enhancers. We also use our computational approach to compare the duplication and splitting of a single enhancer as mechanisms for the generation of shadow enhancers and find that the duplication of enhancers can decrease noise and increase fidelity, although at the metabolic cost of increased RNA production. A saturation mechanism for enhancer interactions similarly improves on both of these metrics. Taken together, this work highlights that shadow enhancer systems may exist for several reasons: genetic drift or the tuning of key functions of enhancers, including transcription fidelity, noise and output.

Decoding enhancer complexity with machine learning and high-throughput discovery

Enhancers are genomic DNA elements controlling spatiotemporal gene expression. Their flexible organization and functional redundancies make deciphering their sequence-function relationships challenging. This article provides an overview of the current understanding of enhancer organization and evolution, with an emphasis on factors that influence these relationships. Technological advancements, particularly in machine learning and synthetic biology, are discussed in light of how they provide new ways to understand this complexity. Exciting opportunities lie ahead as we continue to unravel the intricacies of enhancer function.

Selective Enhancer Dependencies in MYC -Intact and MYC -Rearranged Germinal Center B-cell Diffuse Large B-cell Lymphoma

High expression of MYC and its target genes define a subset of germinal center B-cell diffuse large B-cell lymphoma (GCB-DLBCL) associated with poor outcomes. Half of these high-grade cases show chromosomal rearrangements between the MYC locus and heterologous enhancer-bearing loci, while focal deletions of the adjacent non-coding gene PVT1 are enriched in MYC -intact cases. To identify genomic drivers of MYC activation, we used high-throughput CRISPR-interference (CRISPRi) profiling of candidate enhancers in the MYC locus and rearrangement partner loci in GCB-DLBCL cell lines and mantle cell lymphoma (MCL) comparators that lacked common rearrangements between MYC and immunoglobulin (Ig) loci. Rearrangements between MYC and non-Ig loci were associated with unique dependencies on specific enhancer subunits within those partner loci. Notably, fitness dependency on enhancer modules within the BCL6 super-enhancer ( BCL6 -SE) cluster regulated by a transcription factor complex of MEF2B, POU2F2, and POU2AF1 was higher in cell lines bearing a recurrent MYC::BCL6 -SE rearrangement. In contrast, GCB-DLBCL cell lines without MYC rearrangement were highly dependent on a previously uncharacterized 3' enhancer within the MYC locus itself (GCBME-1), that is regulated in part by the same triad of factors. GCBME-1 is evolutionarily conserved and active in normal germinal center B cells in humans and mice, suggesting a key role in normal germinal center B cell biology. Finally, we show that the PVT1 promoter limits MYC activation by either native or heterologous enhancers and demonstrate that this limitation is bypassed by 3' rearrangements that remove PVT1 from its position in cis with the rearranged MYC gene.

Key points

CRISPR-interference screens identify a conserved germinal center B cell MYC enhancer that is essential for GCB-DLBCL lacking MYC rearrangements. Functional profiling of MYC partner loci reveals principles of MYC enhancer-hijacking activation by non-immunoglobulin rearrangements.

Precise modulation of transcription factor levels identifies features underlying dosage sensitivity

Transcriptional regulation exhibits extensive robustness, but human genetics indicates sensitivity to transcription factor (TF) dosage. Reconciling such observations requires quantitative studies of TF dosage effects at trait-relevant ranges, largely lacking so far. TFs play central roles in both normal-range and disease-associated variation in craniofacial morphology; we therefore developed an approach to precisely modulate TF levels in human facial progenitor cells and applied it to SOX9, a TF associated with craniofacial variation and disease (Pierre Robin sequence (PRS)). Most SOX9-dependent regulatory elements (REs) are buffered against small decreases in SOX9 dosage, but REs directly and primarily regulated by SOX9 show heightened sensitivity to SOX9 dosage; these RE responses partially predict gene expression responses. Sensitive REs and genes preferentially affect functional chondrogenesis and PRS-like craniofacial shape variation. We propose that such REs and genes underlie the sensitivity of specific phenotypes to TF dosage, while buffering of other genes leads to robust, nonlinear dosage-to-phenotype relationships.

3D enhancer-promoter interactions and multi-connected hubs: Organizational principles and functional roles

The spatiotemporal control of gene expression is dependent on the activity of cis-acting regulatory sequences, called enhancers, which regulate target genes over variable genomic distances and, often, by skipping intermediate promoters, suggesting mechanisms that control enhancer-promoter communication. Recent genomics and imaging technologies have revealed highly complex enhancer-promoter interaction networks, whereas advanced functional studies have started interrogating the forces behind the physical and functional communication among multiple enhancers and promoters. In this review, we first summarize our current understanding of the factors involved in enhancer-promoter communication, with a particular focus on recent papers that have revealed new layers of complexities to old questions. In the second part of the review, we focus on a subset of highly connected enhancer-promoter "hubs" and discuss their potential functions in signal integration and gene regulation, as well as the putative factors that might determine their dynamics and assembly.

Massively Parallel Reporter Assays for High-Throughput In Vivo Analysis of Cis-Regulatory Elements

The rapid improvement of descriptive genomic technologies has fueled a dramatic increase in hypothesized connections between cardiovascular gene expression and phenotypes. However, in vivo testing of these hypotheses has predominantly been relegated to slow, expensive, and linear generation of genetically modified mice. In the study of genomic cis-regulatory elements, generation of mice featuring transgenic reporters or cis-regulatory element knockout remains the standard approach. While the data obtained is of high quality, the approach is insufficient to keep pace with candidate identification and therefore results in biases introduced during the selection of candidates for validation. However, recent advances across a range of disciplines are converging to enable functional genomic assays that can be conducted in a high-throughput manner. Here, we review one such method, massively parallel reporter assays (MPRAs), in which the activities of thousands of candidate genomic regulatory elements are simultaneously assessed via the next-generation sequencing of a barcoded reporter transcript. We discuss best practices for MPRA design and use, with a focus on practical considerations, and review how this emerging technology has been successfully deployed in vivo. Finally, we discuss how MPRAs are likely to evolve and be used in future cardiovascular research.

Dynamic network-guided CRISPRi screen reveals CTCF loop-constrained nonlinear enhancer-gene regulatory activity in cell state transitions

Comprehensive enhancer discovery is challenging because most enhancers, especially those affected in complex diseases, have weak effects on gene expression. Our network modeling revealed that nonlinear enhancer-gene regulation during cell state transitions can be leveraged to improve the sensitivity of enhancer discovery. Utilizing hESC definitive endoderm differentiation as a dynamic transition system, we conducted a mid-transition CRISPRi-based enhancer screen. The screen discovered a comprehensive set of enhancers (4 to 9 per locus) for each of the core endoderm lineage-specifying transcription factors, and many enhancers had strong effects mid-transition but weak effects post-transition. Through integrating enhancer activity measurements and three-dimensional enhancer-promoter interaction information, we were able to develop a CTCF loop-constrained Interaction Activity (CIA) model that can better predict functional enhancers compared to models that rely on Hi-C-based enhancer-promoter contact frequency. Our study provides generalizable strategies for sensitive and more comprehensive enhancer discovery in both normal and pathological cell state transitions.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Evolution of modular and pleiotropic enhancers

Cis-regulatory elements (CREs), or enhancers, are segments of noncoding DNA that regulate the spatial and temporal expression of nearby genes. Sometimes, genes are expressed in more than one tissue, and this can be driven by two main types of CREs: tissue-specific "modular" CREs, where different CREs drive expression of the gene in the different tissues, or by "pleiotropic" CREs, where the same CRE drives expression in the different tissues. In this perspective, we will discuss some of the ways (i) modular and pleiotropic CREs might originate; (ii) propose that modular CREs might derive from pleiotropic CREs via a process of duplication, degeneration, and complementation (the CRE-DDC model); and (iii) propose that hotspot loci of evolution are associated with the origin of modular CREs belonging to any gene in a regulatory network.

Widespread cis- and trans-regulatory evolution underlies the origin, diversification, and loss of a sexually dimorphic fruit fly pigmentation trait

Changes in gene expression are a prominent feature of morphological evolution. These changes occur to hierarchical gene regulatory networks (GRNs) of transcription factor genes that regulate the expression of trait-building differentiation genes. While changes in the expression of differentiation genes are essential to phenotypic evolution, they can be caused by mutations within cis-regulatory elements (CREs) that drive their expression (cis-evolution) or within genes for CRE-interacting transcription factors (trans-evolution). Locating these mutations remains a challenge, especially when experiments are limited to one species that possesses the ancestral or derived phenotype. We investigated CREs that control the expression of the differentiation genes tan and yellow, the expression of which evolved during the gain, modification, and loss of dimorphic pigmentation among Sophophora fruit flies. We show these CREs to be necessary components of a pigmentation GRN, as deletion from Drosophila melanogaster (derived dimorphic phenotype) resulted in lost expression and lost male-specific pigmentation. We evaluated the ability of orthologous CRE sequences to drive reporter gene expression in species with modified (Drosophila auraria), secondarily lost (Drosophila ananassae), and ancestrally absent (Drosophila willistoni) pigmentation. We show that the transgene host frequently determines CRE activity, implicating trans-evolution as a significant factor for this trait's diversity. We validated the gain of dimorphic Bab transcription factor expression as a trans-change contributing to the dimorphic trait. Our findings suggest an amenability to change for the landscape of trans-regulators and begs for an explanation as to why this is so common compared to the evolution of differentiation gene CREs.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

A novel role for trithorax in the gene regulatory network for a rapidly evolving fruit fly pigmentation trait

Animal traits develop through the expression and action of numerous regulatory and realizator genes that comprise a gene regulatory network (GRN). For each GRN, its underlying patterns of gene expression are controlled by cis-regulatory elements (CREs) that bind activating and repressing transcription factors. These interactions drive cell-type and developmental stage-specific transcriptional activation or repression. Most GRNs remain incompletely mapped, and a major barrier to this daunting task is CRE identification. Here, we used an in silico method to identify predicted CREs (pCREs) that comprise the GRN which governs sex-specific pigmentation of Drosophila melanogaster. Through in vivo assays, we demonstrate that many pCREs activate expression in the correct cell-type and developmental stage. We employed genome editing to demonstrate that two CREs control the pupal abdomen expression of trithorax, whose function is required for the dimorphic phenotype. Surprisingly, trithorax had no detectable effect on this GRN's key trans-regulators, but shapes the sex-specific expression of two realizator genes. Comparison of sequences orthologous to these CREs supports an evolutionary scenario where these trithorax CREs predated the origin of the dimorphic trait. Collectively, this study demonstrates how in silico approaches can shed novel insights on the GRN basis for a trait's development and evolution.