Shadow enhancers as a source of evolutionary novelty

A complex regulatory landscape involved in the development of mammalian external genitals

Developmental genes are often controlled by large regulatory landscapes matching topologically associating domains (TADs). In various contexts, the associated chromatin backbone is modified by specific enhancer-enhancer and enhancer-promoter interactions. We used a TAD flanking the mouse HoxD cluster to study how these regulatory architectures are formed and deconstructed once their function achieved. We describe this TAD as a functional unit, with several regulatory sequences acting together to elicit a transcriptional response. With one exception, deletion of these sequences didn't modify the transcriptional outcome, a result at odds with a conventional view of enhancer function. The deletion and inversion of a CTCF site located near these regulatory sequences did not affect transcription of the target gene. Slight modifications were nevertheless observed, in agreement with the loop extrusion model. We discuss these unexpected results considering both conventional and alternative explanations relying on the accumulation of poorly specific factors within the TAD backbone.

Loss-of-function tolerance of enhancers in the human genome

Previous studies have surveyed the potential impact of loss-of-function (LoF) variants and identified LoF-tolerant protein-coding genes. However, the tolerance of human genomes to losing enhancers has not yet been evaluated. Here we present the catalog of LoF-tolerant enhancers using structural variants from whole-genome sequences. Using a conservative approach, we estimate that individual human genomes possess at least 28 LoF-tolerant enhancers on average. We assessed the properties of LoF-tolerant enhancers in a unified regulatory network constructed by integrating tissue-specific enhancers and gene-gene interactions. We find that LoF-tolerant enhancers tend to be more tissue-specific and regulate fewer and more dispensable genes relative to other enhancers. They are enriched in immune-related cells while enhancers with low LoF-tolerance are enriched in kidney and brain/neuronal stem cells. We developed a supervised learning approach to predict the LoF-tolerance of all enhancers, which achieved an area under the receiver operating characteristics curve (AUROC) of 98%. We predict 3,519 more enhancers would be likely tolerant to LoF and 129 enhancers that would have low LoF-tolerance. Our predictions are supported by a known set of disease enhancers and novel deletions from PacBio sequencing. The LoF-tolerance scores provided here will serve as an important reference for disease studies.

Enhancers are elements where transcription factors bind and regulate the expression of protein-coding genes. Although multiple previous studies have focused on which genes can tolerate loss-of-function (LoF), none has systematically evaluated the tolerance of all enhancers in the human genome to LoF. Individual studies have shown a broad range of phenotypic effects of enhancer LoF. The phenotypic effects of enhancer LoF likely fall into a spectrum where deletion of LoF-tolerant enhancers would not elicit substantial phenotypic impact, while some enhancers are likely to cause fitness defects when deleted. Here we report a systematic computational approach that uses machine learning and properties of enhancers in a unified human regulatory network with tissue-specific annotations to predict the LoF-tolerance of all enhancers identified in the human genome. The LoF-tolerance scores of enhancers provided in this study can significantly facilitate the interpretation and prioritization of non-coding sequence variants for disease and functional studies.

Three-dimensional genome: developmental technologies and applications in precision medicine

In the 20th century, our familiar structure of DNA was the double helix. Due to technical limitations, we do not have a good way to understand the finer structure of the genome, let alone its transcriptional regulation. Until the advent of 3C technologies, we were no longer blind to this one. Three-dimensional (3D) genomics is a new subject, which mainly studies the 3D structure and transcriptional regulation of eukaryotic genomes. Now, this field mainly has Hi-C series and CHIA-PET series technologies. Through 3D genomics, we can understand the basic structure of DNA, understand the growth and development of organisms and the occurrence of diseases, so as to promote human medical and health undertakings. The review introduces the main research techniques of 3D genomics and their characteristics, the latest development of 3D genome structure, the relationship between diseases and 3D genome structure, the applications of 3D genome in precision medicine, and the development of the 4D nucleome project.

Independent Transposon Exaptation Is a Widespread Mechanism of Redundant Enhancer Evolution in the Mammalian Genome

Many regulatory networks appear to involve partially redundant enhancers. Traditionally, such enhancers have been hypothesized to originate mainly by sequence duplication. An alternative model postulates that they arise independently, through convergent evolution. This mechanism appears to be counterintuitive to natural selection: Redundant sequences are expected to either diverge and acquire new functions or accumulate mutations and become nonfunctional. Nevertheless, we show that at least 31% of the redundant enhancer pairs in the human genome (and 17% in the mouse genome) indeed originated in this manner. Specifically, for virtually all transposon-derived redundant enhancer pairs, both enhancer partners have evolved independently, from the exaptation of two different transposons. In addition to conferring robustness to the system, redundant enhancers could provide an evolutionary advantage by fine-tuning gene expression. Consistent with this hypothesis, we observed that the target genes of redundant enhancers exhibit higher expression levels and tissue specificity as compared with other genes. Finally, we found that although enhancer redundancy appears to be an intrinsic property of certain mammalian regulatory networks, the corresponding enhancers are largely species-specific. In other words, the redundancy in these networks is most likely a result of convergent evolution.

The design and logic of terminal patterning in Drosophila

Terminal regions of the early Drosophila embryo are patterned by the highly conserved ERK cascade, giving rise to the nonsegmented terminal structures of the future larva. In less than an hour, this signaling event establishes several gene expression boundaries and sets in motion a sequence of elaborate morphogenetic events. Genetic studies of terminal patterning discovered signaling components and transcription factors that are involved in numerous developmental contexts and deregulated in human diseases. This review summarizes current understanding of signaling and morphogenesis during terminal patterning and discusses several open questions that can now be rigorously investigated using live imaging, omics, and optogenetic approaches. The anatomical simplicity of the terminal patterning system and its amenability to a broad range of increasingly sophisticated genetic perturbations will continue to make it a premier quantitative model for studying multiple aspects of tissue patterning by dynamically controlled cell signaling pathways.

How to study enhancers in non-traditional insect models

Transcriptional enhancers are central to the function and evolution of genes and gene regulation. At the organismal level, enhancers play a crucial role in coordinating tissue- and context-dependent gene expression. At the population level, changes in enhancers are thought to be a major driving force that facilitates evolution of diverse traits. An amazing array of diverse traits seen in insect morphology, physiology and behavior has been the subject of research for centuries. Although enhancer studies in insects outside of Drosophila have been limited, recent advances in functional genomic approaches have begun to make such studies possible in an increasing selection of insect species. Here, instead of comprehensively reviewing currently available technologies for enhancer studies in established model organisms such as Drosophila, we focus on a subset of computational and experimental approaches that are likely applicable to non-Drosophila insects, and discuss the pros and cons of each approach. We discuss the importance of validating enhancer function and evaluate several possible validation methods, such as reporter assays and genome editing. Key points and potential pitfalls when establishing a reporter assay system in non-traditional insect models are also discussed. We close with a discussion of how to advance enhancer studies in insects, both by improving computational approaches and by expanding the genetic toolbox in various insects. Through these discussions, this Review provides a conceptual framework for studying the function and evolution of enhancers in non-traditional insect models.

Enhancer Domains Predict Gene Pathogenicity and Inform Gene Discovery in Complex Disease

Non-coding transcriptional regulatory elements are critical for controlling the spatiotemporal expression of genes. Here, we demonstrate that the sizes and number of enhancers linked to a gene reflect its disease pathogenicity. Moreover, genes with redundant enhancer domains are depleted of cis-acting genetic variants that disrupt gene expression, and they are buffered against the effects of disruptive non-coding mutations. Our results demonstrate that dosage-sensitive genes have evolved a robustness to the disruptive effects of genetic variation by expanding their regulatory domains. This solves a puzzle about why genes associated with human disease are depleted of cis-eQTLs (cis-expression quantitative trait loci), suggesting that this relationship might complicate gene identification in causal genome-wide association studies (GWASs) using eQTL information, and establishes a framework for identifying non-coding regulatory variation with phenotypic consequences.

Enhancer jungles establish robust tissue-specific regulatory control in the human genome

An increasing number of studies suggest that functionally redundant enhancers safeguard development via buffering gene expression against environmental and genetic perturbations. Here, we identified over-represented clusters of enhancers (enhancer jungles or EJs) in human B lymphoblastoid cells. We found that EJs tend to associate with genes involved in the activation of the immune system response. Although spanning multiple genes, the enhancers within an EJ tend to collaborate with each other on regulating a single gene. The employment of homotypic transcription factor binding sites (TFBSs) in EJ enhancers and heterotypic TFBSs between constituent enhancers within an EJ may safeguard a robust transcriptional output of the target gene. EJ enhancers evolve under a weaker selective pressure compared to regular enhancers (REs), and approximately 35% of EJs do not have orthologues in the mouse genome. In GM12878, these human-specific EJs appear to regulate genes associated with the adaptive immune system response, while the conserved EJs are associated with innate immunity. Recently acquired human EJs are associated with the higher level of target gene expression compared with conserved EJs, thus facilitating the environmental adaptation of the organism during evolution. In short, the existence of EJs is a common regulatory architecture conferring a robust regulatory control for key lineage genes.

Two pair-rule responsive enhancers regulate wingless transcription in the Drosophila blastoderm embryo

BACKGROUND

While many developmentally relevant enhancers act in a modular fashion, there is growing evidence for nonadditive interactions between distinct cis-regulatory enhancers. We investigated if nonautonomous enhancer interactions underlie transcription regulation of the Drosophila segment polarity gene, wingless.

RESULTS

We identified two wg enhancers active at the blastoderm stage: wg 3613u, located from -3.6 to -1.3 kb upstream of the wg transcription start site (TSS) and 3046d, located in intron two of the wg gene, from 3.0 to 4.6 kb downstream of the TSS. Genetic experiments confirm that Even Skipped (Eve), Fushi-tarazu (Ftz), Runt, Odd-paired (Opa), Odd-skipped (Odd), and Paired (Prd) contribute to spatially regulated wg expression. Interestingly, there are enhancer specific differences in response to the gain or loss of function of pair-rule gene activity. Although each element recapitulates aspects of wg expression, a composite reporter containing both enhancers more faithfully recapitulates wg regulation than would be predicted from the sum of their individual responses.

CONCLUSION

These results suggest that the regulation of wg by pair-rule genes involves nonadditive interactions between distinct cis-regulatory enhancers.

Coacting enhancers can have complementary functions within gene regulatory networks and promote canalization

Developmental genes are often regulated by multiple enhancers exhibiting similar spatiotemporal outputs, which are generally considered redundantly acting though few have been studied functionally. Using CRISPR-Cas9, we created deletions of two enhancers, brk5' and brk3', that drive similar but not identical expression of the gene brinker (brk) in early Drosophila embryos. Utilizing both in situ hybridization and quantitative mRNA analysis, we investigated the changes in the gene network state caused by the removal of one or both of the early acting enhancers. brk5' deletion generally phenocopied the gene mutant, including expansion of the BMP ligand decapentaplegic (dpp) as well as inducing variability in amnioserosa tissue cell number suggesting a loss of canalization. In contrast, brk3' deletion presented unique phenotypes including dorsal expansion of several ventrally expressed genes and a decrease in amnioserosa cell number. Similarly, deletions were made for two enhancers associated with the gene short-gastrulation (sog), sog.int and sog.dist, demonstrating that they also exhibit distinct patterning phenotypes and affect canalization. In summary, this study shows that similar gene expression driven by coacting enhancers can support distinct, and sometimes complementary, functions within gene regulatory networks and, moreover, that phenotypes associated with individual enhancer deletion mutants can provide insight into new gene functions.

Combinatorial chromatin dynamics foster accurate cardiopharyngeal fate choices

During embryogenesis, chromatin accessibility profiles control lineage-specific gene expression by modulating transcription, thus impacting multipotent progenitor states and subsequent fate choices. Subsets of cardiac and pharyngeal/head muscles share a common origin in the cardiopharyngeal mesoderm, but the chromatin landscapes that govern multipotent progenitors competence and early fate choices remain largely elusive. Here, we leveraged the simplicity of the chordate model Ciona to profile chromatin accessibility through stereotyped transitions from naive Mesp+ mesoderm to distinct fate-restricted heart and pharyngeal muscle precursors. An FGF-Foxf pathway acts in multipotent progenitors to establish cardiopharyngeal-specific patterns of accessibility, which govern later heart vs. pharyngeal muscle-specific expression profiles, demonstrating extensive spatiotemporal decoupling between early cardiopharyngeal enhancer accessibility and late cell-type-specific activity. We found that multiple cis-regulatory elements, with distinct chromatin accessibility profiles and motif compositions, are required to activate Ebf and Tbx1/10, two key determinants of cardiopharyngeal fate choices. We propose that these 'combined enhancers' foster spatially and temporally accurate fate choices, by increasing the repertoire of regulatory inputs that control gene expression, through either accessibility and/or activity.

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.

Enhancer Dysfunction in 3D Genome and Disease

Spatiotemporal patterns of gene expression depend on enhancer elements and other factors during individual development and disease progression. The rapid progress of high-throughput techniques has led to well-defined enhancer chromatin properties. Various genome-wide methods have revealed a large number of enhancers and the discovery of three-dimensional (3D) genome architecture showing the distant interacting mechanisms of enhancers that loop to target gene promoters. Whole genome sequencing projects directed at cancer have led to the discovery of substantial enhancer dysfunction in misregulating gene expression and in tumor initiation and progression. Results from genome-wide association studies (GWAS) combined with functional genomics analyses have elucidated the functional impacts of many cancer risk-associated variants that are enriched within the enhancer regions of chromatin. Risk variants dysregulate the expression of enhancer variant-associated genes via 3D genomic interactions. Moreover, these enhancer variants often alter the chromatin binding affinity for cancer-relevant transcription factors, which in turn leads to aberrant expression of the genes associated with cancer susceptibility. In this review, we investigate the extent to which these genetic regulatory circuits affect cancer predisposition and how the recent development of genome-editing methods have enabled the determination of the impacts of genomic variation and alteration on cancer phenotype, which will eventually lead to better management plans and treatment responses to human cancer in the clinic.

Functional Consequences of Chromosomal Rearrangements on Gene Expression: Not So Deleterious After All?

Chromosomes are folded and organized into topologically associating domains (TADs) which provide a framework for the interaction of enhancers with the promoter of their target gene(s). Structural rearrangements observed during evolution or in disease contexts suggest that changes in genome organization strongly affect gene expression and can have drastic phenotypic effects. In this review, I will discuss how recent genomic engineering experiments reveal a more contrasted picture, suggesting that TADs are important but not always essential for gene expression regulation.

Sufficiency analysis of estrogen responsive enhancers using synthetic activators

Enhancers bound by synthetic activators can recreate a transcriptional response to estrogen, drive different levels of gene expression, and work independently to regulate transcription.

Multiple regulatory regions bound by the same transcription factor have been shown to simultaneously control a single gene’s expression. However, it remains unclear how these regulatory regions combine to regulate transcription. Here, we test the sufficiency of promoter-distal estrogen receptor α-binding sites (ERBSs) for activating gene expression by recruiting synthetic activators in the absence of estrogens. Targeting either dCas9-VP16(10x) or dCas9-p300(core) to ERBS induces H3K27ac and activates nearby expression in a manner similar to an estrogen induction, with dCas9-VP16(10x) acting as a stronger activator. The sufficiency of individual ERBSs is highly correlated with their necessity, indicating an inherent activation potential that is associated with the binding of RNA polymerase II and several transcription factors. By targeting ERBS combinations, we found that ERBSs work independently to control gene expression when bound by synthetic activators. The sufficiency results contrast necessity assays that show synergy between these ERBSs, suggesting that synergy occurs between ERBSs in terms of activator recruitment, whereas directly recruiting activators leads to independent effects on gene expression.

Transcriptional autoregulation of zebrafish tbx6 is required for somite segmentation

The presumptive somite boundary in the presomitic mesoderm (PSM) is defined by the anterior border of the expression domain of Tbx6 protein. During somite segmentation, the expression domain of Tbx6 is regressed by Ripply-meditated degradation of Tbx6 protein. Although the expression of zebrafish tbx6 remains restricted to the PSM, the transcriptional regulation of tbx6 remains poorly understood. Here, we show that the expression of zebrafish tbx6 is maintained by transcriptional autoregulation. We find that a proximal-located cis-regulatory module, TR1, which contains two putative T-box sites, is required for somite segmentation in the intermediate body and for proper expression of segmentation genes. Embryos with deletion of TR1 exhibit significant reduction of tbx6 expression at the 12-somite stage, although its expression is initially observed. Additionally, Tbx6 is associated with TR1 and activates its own expression in the anterior PSM. Furthermore, the anterior expansion of tbx6 expression in ripply gene mutants is suppressed in a TR1-dependent manner. The results suggest that the autoregulatory loop of zebrafish tbx6 facilitates immediate removal of Tbx6 protein through termination of its own transcription at the anterior PSM.

Chromatin-associated RNAs as facilitators of functional genomic interactions

Mammalian genomes are extensively transcribed, which produces a large number of both coding and non-coding transcripts. Various RNAs are physically associated with chromatin, through being either retained in cis at their site of transcription or recruited in trans to other genomic regions. Driven by recent technological innovations for detecting chromatin-associated RNAs, diverse roles are being revealed for these RNAs and associated RNA-binding proteins (RBPs) in gene regulation and genome function. Such functions include locus-specific roles in gene activation and silencing, as well as emerging roles in higher-order genome organization, such as involvement in long-range enhancer-promoter interactions, transcription hubs, heterochromatin, nuclear bodies and phase transitions.

Modular tissue-specific regulation of doublesex underpins sexually dimorphic development in Drosophila

The ability of a single genome to produce distinct and often dramatically different male and female forms is one of the wonders of animal development. In Drosophila melanogaster, most sexually dimorphic traits are controlled by sex-specific isoforms of the doublesex (dsx) transcription factor, and dsx expression is mostly limited to cells that give rise to sexually dimorphic traits. However, it is unknown how this mosaic of sexually dimorphic and monomorphic organs arises. Here, we characterize the cis-regulatory sequences that control dsx expression in the foreleg, which contains multiple types of sex-specific sensory organs. We find that separate modular enhancers are responsible for dsx expression in each sexually dimorphic organ. Expression of dsx in the sex comb is co-regulated by two enhancers with distinct spatial and temporal specificities that are separated by a genitalia-specific enhancer. The sex comb-specific enhancer from D. willistoni, a species that primitively lacks sex combs, is not active in the foreleg. Thus, the mosaic of sexually dimorphic and monomorphic organs depends on modular regulation of dsx transcription by dedicated cell type-specific enhancers.

Structural Redundancy in Supracellular Actomyosin Networks Enables Robust Tissue Folding

Tissue morphogenesis is strikingly robust. Yet, how tissues are sculpted under challenging conditions is unknown. Here, we combined network analysis, experimental perturbations, and computational modeling to determine how network connectivity between hundreds of contractile cells on the ventral side of the Drosophila embryo ensures robust tissue folding. We identified two network properties that mechanically promote robustness. First, redundant supracellular cytoskeletal network paths ensure global connectivity, even with network degradation. By forming many more connections than are required, morphogenesis is not disrupted by local network damage, analogous to the way redundancy guarantees the large-scale function of vasculature and transportation networks. Second, directional stiffening of edges oriented orthogonal to the folding axis promotes furrow formation at lower contractility levels. Structural redundancy and directional network stiffening ensure robust tissue folding with proper orientation.

Super-enhancers: critical roles and therapeutic targets in hematologic malignancies

Super-enhancers (SEs) in a broad range of human cell types are large clusters of enhancers with aberrant high levels of transcription factor binding, which are central to drive expression of genes in controlling cell identity and stimulating oncogenic transcription. Cancer cells acquire super-enhancers at oncogene and cancerous phenotype relies on these abnormal transcription propelled by SEs. Furthermore, specific inhibitors targeting SEs assembly and activation have offered potential targets for treating various tumors including hematological malignancies. Here, we first review the identification, functional significance of SEs. Next, we summarize recent findings of SEs and SE-driven gene regulation in normal hematopoiesis and hematologic malignancies. The importance and various modes of SE-mediated MYC oncogene amplification are illustrated. Finally, we highlight the progress of SEs as selective therapeutic targets in basic research and clinical trials. Some open questions regarding functional significance and future directions of targeting SEs in the clinic will be discussed too.

Multi-enhancer transcriptional hubs confer phenotypic robustness

We previously showed in Drosophila melanogaster embryos that low-affinity Ultrabithorax (Ubx)-responsive shavenbaby (svb) enhancers drive expression using localized transcriptional environments and that active svb enhancers on different chromosomes tended to colocalize (Tsai et al., 2017). Here, we test the hypothesis that these multi-enhancer 'hubs' improve phenotypic resilience to stress by buffering against decreases in transcription factor concentrations and transcriptional output. Deleting a redundant enhancer from the svb locus led to reduced trichome numbers in embryos raised at elevated temperatures. Using high-resolution fluorescence microscopy, we observed lower Ubx concentration and transcriptional output in this deletion allele. Transcription sites of the full svb cis-regulatory region inserted into a different chromosome colocalized with the svb locus, increasing Ubx concentration, the transcriptional output of svb, and partially rescuing the phenotype. Thus, multiple enhancers could reinforce a local transcriptional hub to buffer against environmental stresses and genetic perturbations, providing a mechanism for phenotypical robustness.

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

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

p63 cooperates with CTCF to modulate chromatin architecture in skin keratinocytes

The transcription factor p63 regulates epidermal genes and the enhancer landscape in skin keratinocytes. Its molecular function in controlling the chromatin structure is, however, not yet completely understood. Here, we integrated multi-omics profiles, including the transcriptome, transcription factor DNA-binding and chromatin accessibility, in skin keratinocytes isolated from EEC syndrome patients carrying p63 mutations, to examine the role of p63 in shaping the chromatin architecture. We found decreased chromatin accessibility in p63- and CTCF-bound open chromatin regions that potentially contributed to gene deregulation in mutant keratinocytes. Cooperation of p63 and CTCF seemed to assist chromatin interactions between p63-bound enhancers and gene promoters in skin keratinocytes. Our study suggests an intriguing model where cell type-specific transcription factors such as p63 cooperate with the genome organizer CTCF in the three-dimensional chromatin space to regulate the transcription program important for the proper cell identity.

Perturbing Enhancer Activity in Cancer Therapy

Tight regulation of gene transcription is essential for normal development, tissue homeostasis, and disease-free survival. Enhancers are distal regulatory elements in the genome that provide specificity to gene expression programs and are frequently misregulated in cancer. Recent studies examined various enhancer-driven malignant dependencies and identified different approaches to specifically target these programs. In this review, we describe numerous features that make enhancers good transcriptional targets in cancer therapy and discuss different approaches to overcome enhancer perturbation. Interestingly, a number of approved therapeutic agents, such as cyclosporine, steroid hormones, and thiazolidinediones, actually function by affecting enhancer landscapes by directly targeting very specific transcription factor programs. More recently, a broader approach to targeting deregulated enhancer programs has been achieved via Bromodomain and Extraterminal (BET) inhibition or perturbation of transcription-related cyclin-dependent kinases (CDK). One challenge to enhancer-targeted therapy is proper patient stratification. We suggest that monitoring of enhancer RNA (eRNA) expression may serve as a unique biomarker of enhancer activity that can help to predict and monitor responsiveness to enhancer-targeted therapies. A more thorough investigation of cancer-specific enhancers and the underlying mechanisms of deregulation will pave the road for an effective utilization of enhancer modulators in a precision oncology approach to cancer treatment.

The Drosophila Pioneer Factor Zelda Modulates the Nuclear Microenvironment of a Dorsal Target Enhancer to Potentiate Transcriptional Output

Connecting the developmental patterning of tissues to the mechanistic control of RNA polymerase II remains a long-term goal of developmental biology. Many key elements have been identified in the establishment of spatial-temporal control of transcription in the early Drosophila embryo, a model system for transcriptional regulation. The dorsal-ventral axis of the Drosophila embryo is determined by the graded distribution of Dorsal (Dl), a homolog of the nuclear factor κB (NF-κB) family of transcriptional activators found in humans [1, 2]. A second maternally deposited factor, Zelda (Zld), is uniformly distributed in the embryo and is thought to act as a pioneer factor, increasing enhancer accessibility for transcription factors, such as Dl [3-9]. Here, we utilized the MS2 live imaging system to evaluate the expression of the Dl target gene short gastrulation (sog) to better understand how a pioneer factor affects the kinetic parameters of transcription. Our experiments indicate that Zld modifies probability of activation, the timing of this activation, and the rate at which transcription occurs. Our results further show that this effective rate increase is due to an increased accumulation of Dl at the site of transcription, suggesting that transcription factor "hubs" induced by Zld [10] functionally regulate transcription.

Shaping the nebulous enhancer in the era of high-throughput assays and genome editing

Since the 1st discovery of transcriptional enhancers in 1981, their textbook definition has remained largely unchanged in the past 37 years. With the emergence of high-throughput assays and genome editing, which are switching the paradigm from bottom-up discovery and testing of individual enhancers to top-down profiling of enhancer activities genome-wide, it has become increasingly evidenced that this classical definition has left substantial gray areas in different aspects. Here we survey a representative set of recent research articles and report the definitions of enhancers they have adopted. The results reveal that a wide spectrum of definitions is used usually without the definition stated explicitly, which could lead to difficulties in data interpretation and downstream analyses. Based on these findings, we discuss the practical implications and suggestions for future studies.

Redundant and Cryptic Enhancer Activities of the Drosophila yellow Gene

Cis-regulatory sequences known as enhancers play a key role in regulating gene expression. Evolutionary changes in these DNA sequences contribute to phenotypic evolution. The Drosophila yellow gene, which is required for pigmentation, has emerged as a model system for understanding how cis-regulatory sequences evolve, providing some of the most detailed insights available into how activities of orthologous enhancers have diverged between species. Here, we examine the evolution of yellow cis-regulatory sequences on a broader scale, by comparing the distribution and function of yellow enhancer activities throughout the 5' intergenic and intronic sequences of Drosophila melanogaster, D. pseudoobscura, and D. willistoni We find that cis-regulatory sequences driving expression in a particular tissue are not as modular as previously described, but rather have many redundant and cryptic enhancer activities distributed throughout the regions surveyed. Interestingly, cryptic enhancer activities of sequences from one species often drove patterns of expression observed in other species, suggesting that the frequent evolutionary changes in yellow expression observed among Drosophila species may be facilitated by gaining and losing repression of preexisting cis-regulatory sequences.

Signal Integration by Shadow Enhancers and Enhancer Duplications Varies across the Drosophila Embryo

Transcription of developmental genes is controlled by multiple enhancers. Frequently, more than one enhancer can activate transcription from the same promoter in the same cells. How is regulatory information from multiple enhancers combined to determine the overall expression output? We measure nascent transcription driven by a pair of shadow enhancers, each enhancer of the pair separately, and each duplicated, using live imaging in Drosophila embryos. This set of constructs allows us to quantify the input-output function describing signal integration by two enhancers. We show that signal integration performed by these shadow enhancers and duplications varies across the expression pattern, implying that how their activities are combined depends on the transcriptional regulators bound to the enhancers in different parts of the embryo. Characterizing signal integration by multiple enhancers is a critical step in developing conceptual and computational models of gene expression at the locus level, where multiple enhancers control transcription together.

The cis-Regulatory Atlas of the Mouse Immune System

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

A Multiplexed Assay for Exon Recognition Reveals that an Unappreciated Fraction of Rare Genetic Variants Cause Large-Effect Splicing Disruptions

Mutations that lead to splicing defects can have severe consequences on gene function and cause disease. Here, we explore how human genetic variation affects exon recognition by developing a multiplexed functional assay of splicing using Sort-seq (MFASS). We assayed 27,733 variants in the Exome Aggregation Consortium (ExAC) within or adjacent to 2,198 human exons in the MFASS minigene reporter and found that 3.8% (1,050) of variants, most of which are extremely rare, led to large-effect splice-disrupting variants (SDVs). Importantly, we find that 83% of SDVs are located outside of canonical splice sites, are distributed evenly across distinct exonic and intronic regions, and are difficult to predict a priori. Our results indicate extant, rare genetic variants can have large functional effects on splicing at appreciable rates, even outside the context of disease, and MFASS enables their empirical assessment at scale.

Integrative view on how erythropoietin signaling controls transcription patterns in erythroid cells

PURPOSE OF REVIEW

Erythropoietin (EPO) is necessary and sufficient to trigger dynamic transcriptional patterns that drive the differentiation of erythroid precursor cells into mature, enucleated red cells. Because the molecular cloning and Food and Drug Administration approval for the therapeutic use of EPO over 30 years ago, a detailed understanding of how EPO works has advanced substantially. Yet, the precise epigenetic and transcriptional mechanisms by which EPO signaling controls erythroid expression patterns remains poorly understood. This review focuses on the current state of erythroid biology in regards to EPO signaling from human genetics and functional genomics perspectives.

RECENT FINDINGS

The goal of this review is to provide an integrative view of the gene regulatory underpinnings for erythroid expression patterns that are dynamically shaped during erythroid differentiation. Here, we highlight vignettes connecting recent insights into a genome-wide association study linking an EPO mutation to anemia, a study linking EPO-signaling to signal transducer and activator of transcription 5 (STAT5) chromatin occupancy and enhancers, and studies that examine the molecular mechanisms driving topological chromatin organization in erythroid cells.

SUMMARY

The genetic, epigenetic, and gene regulatory mechanisms underlying how hormone signal transduction influences erythroid gene expression remains only partly understood. A detailed understanding of these molecular pathways and how they intersect with one another will provide the basis for novel strategies to treat anemia and potentially other hematological diseases. As new regulators and signal transducers of EPO-signaling continue to emerge, new clinically relevant targets may be identified that improve the specificity and effectiveness of EPO therapy.

A Genome-wide Framework for Mapping Gene Regulation via Cellular Genetic Screens

Over one million candidate regulatory elements have been identified across the human genome, but nearly all are unvalidated and their target genes uncertain. Approaches based on human genetics are limited in scope to common variants and in resolution by linkage disequilibrium. We present a multiplex, expression quantitative trait locus (eQTL)-inspired framework for mapping enhancer-gene pairs by introducing random combinations of CRISPR/Cas9-mediated perturbations to each of many cells, followed by single-cell RNA sequencing (RNA-seq). Across two experiments, we used dCas9-KRAB to perturb 5,920 candidate enhancers with no strong a priori hypothesis as to their target gene(s), measuring effects by profiling 254,974 single-cell transcriptomes. We identified 664 (470 high-confidence) cis enhancer-gene pairs, which were enriched for specific transcription factors, non-housekeeping status, and genomic and 3D conformational proximity to their target genes. This framework will facilitate the large-scale mapping of enhancer-gene regulatory interactions, a critical yet largely uncharted component of the cis-regulatory landscape of the human genome.

Functional dissection of the ARGONAUTE7 promoter

ARGONAUTES are the central effector proteins of RNA silencing which bind target transcripts in a small RNA-guided manner. Arabidopsis thaliana has 10 ARGONAUTE (AGO) genes, with specialized roles in RNA-directed DNA methylation, post-transcriptional gene silencing, and antiviral defense. To better understand specialization among AGO genes at the level of transcriptional regulation we tested a library of 1497 transcription factors for binding to the promoters of AGO1,AGO10, and AGO7 using yeast 1-hybrid assays. A ranked list of candidate DNA-binding TFs revealed binding of the AGO7 promoter by a number of proteins in two families: the miR156-regulated SPL family and the miR319-regulated TCP family, both of which have roles in developmental timing and leaf morphology. Possible functions for SPL and TCP binding are unclear: we showed that these binding sites are not required for the polar expression pattern of AGO7, nor for the function of AGO7 in leaf shape. Normal AGO7 transcription levels and function appear to depend instead on an adjacent 124-bp region. Progress in understanding the structure of this promoter may aid efforts to understand how the conserved AGO7-triggered TAS3 pathway functions in timing and polarity.

Dichotomy in redundant enhancers points to presence of initiators of gene regulation

Background

The regulatory landscape of a gene locus often consists of several functionally redundant enhancers establishing phenotypic robustness and evolutionary stability of its regulatory program. However, it is unclear what mechanisms are employed by redundant enhancers to cooperatively orchestrate gene expression.

Results

By comparing redundant enhancers to single enhancers (enhancers present in a single copy in a gene locus), we observed that the DNA sequence encryption differs between these two classes of enhancers, suggesting a difference in their regulatory mechanisms. Initiator enhancers, which are a subset of redundant enhancers and show similar sequence encryption to single enhancers, differ from the rest of redundant enhancers in their sequence encryption, evolutionary conservation and proximity to target genes. Genes hosting initiator enhancers in their loci feature elevated levels of expression. Initiator enhancers show a high level of 3D chromatin contacts with both transcription start sites and regular enhancers, suggesting their roles as primary activators and intermediate catalysts of gene expression, through which the regulatory signals of redundant enhancers are propagated to the target genes. In addition, GWAS and eQTLs variants are significantly enriched in initiator enhancers compared to redundant enhancers, suggesting a key functional role these sequences play in gene regulation.

Conclusions

The specific characteristics and widespread abundance of initiator enhancers advocate for a possible universal hierarchical mechanism of tissue-specific gene regulation involving multiple redundant enhancers acting through initiator enhancers.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5335-0) contains supplementary material, which is available to authorized users.

Studying Transcriptional Enhancers: The Founder Fallacy, Validation Creep, and Other Biases

Transcriptional enhancers play a major role in regulating metazoan gene expression. Recent developments in genomics and next-generation sequencing have accelerated and revitalized the study of this important class of sequence elements. Increased interest and attention, however, has also led to troubling trends in the enhancer literature. In this Opinion, I describe some of these issues and show how they arise from shifting and nonuniform enhancer definitions, and genome-era biases. I discuss how they can lead to interpretative errors and an unduly narrow focus on certain aspects of enhancer biology to the potential exclusion of others.

Promoter-enhancer looping and shadow enhancers of the mouse αA-crystallin locus

Gene regulation by enhancers is important for precise temporal and spatial gene expression. Enhancers can drive gene expression regardless of their location, orientation or distance from the promoter. Changes in chromatin conformation and chromatin looping occur to bring the promoter and enhancers into close proximity. αA-crystallin ranks among one of the most abundantly expressed genes and proteins in the mammalian lens. The αA-crystallin locus is characterized by a 16 kb chromatin domain marked by two distal enhancers, 5′ DCR1 and 3′ DCR3. Here we used chromatin conformation capture (3C) analysis and transgenic approaches to analyze temporal control of the mouse αA-crystallin gene. We find that DCR1 is necessary, but not sufficient alone to drive expression at E10.5 in the mouse lens pit. Chromatin looping revealed interaction between the promoter and the region 3′ to DCR1, identifying a novel enhancer region in the αA-crystallin locus. We determined that this novel enhancer region, DCR1S, recapitulates the temporal control by DCR1. Acting as shadow enhancers, DCR1 and DCR1S are able to control expression in the lens vesicle at E11.5. It remains to be elucidated however, which region of the αA-crystallin locus is responsible for expression in the lens pit at E10.5.

Summary: The αA-crystallin ranks amongst the most highly expressed tissue-specific genes. It is an advantageous model system to probe both promoter-enhancer looping and to identify distal enhancers and their temporal/spatial activities.

Buffering and Amplifying Transcriptional Noise During Cell Fate Specification

The molecular processes that drive gene transcription are inherently noisy. This noise often manifests in the form of transcriptional bursts, producing fluctuations in gene activity over time. During cell fate specification, this noise is often buffered to ensure reproducible developmental outcomes. However, sometimes noise is utilized as a “bet-hedging” mechanism to diversify functional roles across a population of cells. Studies of bacteria, yeast, and cultured cells have provided insights into the nature and roles of noise in transcription, yet we are only beginning to understand the mechanisms by which noise influences the development of multicellular organisms. Here we discuss the sources of transcriptional noise and the mechanisms that either buffer noise to drive reproducible fate choices or amplify noise to randomly specify fates.

Heart enhancers with deeply conserved regulatory activity are established early in zebrafish development

During the phylotypic period, embryos from different genera show similar gene expression patterns, implying common regulatory mechanisms. Here we set out to identify enhancers involved in the initial events of cardiogenesis, which occurs during the phylotypic period. We isolate early cardiac progenitor cells from zebrafish embryos and characterize 3838 open chromatin regions specific to this cell population. Of these regions, 162 overlap with conserved non-coding elements (CNEs) that also map to open chromatin regions in human. Most of the zebrafish conserved open chromatin elements tested drive gene expression in the developing heart. Despite modest sequence identity, human orthologous open chromatin regions recapitulate the spatial temporal expression patterns of the zebrafish sequence, potentially providing a basis for phylotypic gene expression patterns. Genome-wide, we discover 5598 zebrafish-human conserved open chromatin regions, suggesting that a diverse repertoire of ancient enhancers is established prior to organogenesis and the phylotypic period.

During early embryogenesis, critical cardiac specification events occur. Here the authors isolate cardiac progenitor cells from early zebrafish embryos and characterize accessible chromatin regions specific to this cell population, finding that many of these regions overlap with conserved non-coding elements that are ortholgous to accessible chromatin regions in human.

Augmentation of a wound response element accompanies the origin of a Hox-regulated Drosophila abdominal pigmentation trait

A challenge for evolutionary research is to uncover how new morphological traits evolve the coordinated spatial and temporal expression patterns of genes that govern their formation during development. Detailed studies are often limited to characterizing how one or a few genes contributed to a trait's emergence, and thus our knowledge of how entire GRNs evolve their coordinated expression of each gene remains unresolved. The melanic color patterns decorating the male abdominal tergites of Drosophila (D.) melanogaster evolved in part by novel expression patterns for genes acting at the terminus of a pigment metabolic pathway, driven by cis-regulatory elements (CREs) with distinct mechanisms of Hox regulation. Here, we examined the expression and evolutionary histories of two important enzymes in this pathway, encoded by the pale and Ddc genes. We found that while both genes exhibit dynamic patterns of expression, a robust pattern of Ddc expression specifically evolved in the lineage of fruit flies with pronounced melanic abdomens. Derived Ddc expression requires the activity of a CRE previously shown to activate expression in response to epidermal wounding. We show that a binding site for the Grainy head transcription factor that promotes the ancestral wound healing function of this CRE is also required for abdominal activity. Together with previous findings in this system, our work shows how the GRN for a novel trait emerged by assembling unique yet similarly functioning CREs from heterogeneous starting points.

Diversification of transcription factor-DNA interactions and the evolution of gene regulatory networks

Sequence-specific transcription factors (TFs) bind short DNA sequences in the genome to regulate the expression of target genes. In the last decade, numerous technical advances have enabled the determination of the DNA-binding specificities of many of these factors. Large-scale screens of many TFs enabled the creation of databases of TF DNA-binding specificities, typically represented as position weight matrices (PWMs). Although great progress has been made in determining and predicting binding specificities systematically, there are still many surprises to be found when studying a particular TF's interactions with DNA in detail. Paralogous TFs' binding specificities can differ in subtle ways, in a manner that is not immediately apparent from looking at their PWMs. These differences affect gene regulatory outputs and enable TFs to rewire transcriptional networks over evolutionary time. This review discusses recent observations made in the study of TF-DNA interactions that highlight the importance of continued in-depth analysis of TF-DNA interactions and their inherent complexity. This article is categorized under: Biological Mechanisms > Regulatory Biology.

Enhancer and super-enhancer: Positive regulators in gene transcription

Enhancer is a positive regulator for spatiotemporal development in eukaryotes. As a cluster, super-enhancer is closely related to cell identity- and fate-determined processes. Both of them function tightly depending on their targeted transcription factors, cofactors, and genes through distal genomic interactions. They have been recognized as critical components and played positive roles in transcriptional regulatory network or factory. Recent advances of next-generation sequencing have dramatically expanded our ability and knowledge to interrogate the molecular mechanism of enhancer and super-enhancer for transcription. Here, we review the history, importance, advances and challenges on enhancer and super-enhancer field. This will benefit our understanding of their function mechanism for transcription underlying precise gene expression.

Diet-Induced Circadian Enhancer Remodeling Synchronizes Opposing Hepatic Lipid Metabolic Processes

Overnutrition disrupts circadian metabolic rhythms by mechanisms that are not well understood. Here, we show that diet-induced obesity (DIO) causes massive remodeling of circadian enhancer activity in mouse liver, triggering synchronous high-amplitude circadian rhythms of both fatty acid (FA) synthesis and oxidation. SREBP expression was rhythmically induced by DIO, leading to circadian FA synthesis and, surprisingly, FA oxidation (FAO). DIO similarly caused a high-amplitude circadian rhythm of PPARα, which was also required for FAO. Provision of a pharmacological activator of PPARα abrogated the requirement of SREBP for FAO (but not FA synthesis), suggesting that SREBP indirectly controls FAO via production of endogenous PPARα ligands. The high-amplitude rhythm of PPARα imparted time-of-day-dependent responsiveness to lipid-lowering drugs. Thus, acquisition of rhythmicity for non-core clock components PPARα and SREBP1 remodels metabolic gene transcription in response to overnutrition and enables a chronopharmacological approach to metabolic disorders.

Cooperation, cis-interactions, versatility and evolutionary plasticity of multiple cis-acting elements underlie krox20 hindbrain regulation

Cis-regulation plays an essential role in the control of gene expression, and is particularly complex and poorly understood for developmental genes, which are subject to multiple levels of modulation. In this study, we performed a global analysis of the cis-acting elements involved in the control of the zebrafish developmental gene krox20. krox20 encodes a transcription factor required for hindbrain segmentation and patterning, a morphogenetic process highly conserved during vertebrate evolution. Chromatin accessibility analysis reveals a cis-regulatory landscape that includes 6 elements participating in the control of initiation and autoregulatory aspects of krox20 hindbrain expression. Combining transgenic reporter analyses and CRISPR/Cas9-mediated mutagenesis, we assign precise functions to each of these 6 elements and provide a comprehensive view of krox20 cis-regulation. Three important features emerged. First, cooperation between multiple cis-elements plays a major role in the regulation. Cooperation can surprisingly combine synergy and redundancy, and is not restricted to transcriptional enhancer activity (for example, 4 distinct elements cooperate through different modes to maintain autoregulation). Second, several elements are unexpectedly versatile, which allows them to be involved in different aspects of control of gene expression. Third, comparative analysis of the elements and their activities in several vertebrate species reveals that this versatility is underlain by major plasticity across evolution, despite the high conservation of the gene expression pattern. These characteristics are likely to be of broad significance for developmental genes.

Enhancer Logic and Mechanics in Development and Disease

Enhancers are distally located genomic cis-regulatory elements that integrate spatiotemporal cues to coordinate gene expression in a tissue-specific manner during metazoan development. Enhancer function depends on a combination of bound transcription factors and cofactors that regulate local chromatin structure, as well as on the topological interactions that are necessary for their activity. Numerous genome-wide studies concur that the vast majority of disease-associated variations occur within non-coding genomic sequences, in other words the 'cis-regulome', and this underscores their relevance for human health. Advances in DNA sequencing and genome-editing technologies have dramatically expanded our ability to identify enhancers and investigate their properties in vivo, revealing an extraordinary level of interconnectivity underlying cis-regulatory networks. We discuss here these recently developed methodologies, as well as emerging trends and remaining questions in the field of enhancer biology, and how perturbation of enhancer activities/functions results in enhanceropathies.

A Distalless-responsive enhancer of the Hox gene Sex combs reduced is required for segment- and sex-specific sensory organ development in Drosophila

Hox genes are involved in the patterning of animal body parts at multiple levels of regulatory hierarchies. Early expression of Hox genes in different domains along the embryonic anterior-posterior (A/P) axis in insects, vertebrates, and other animals establishes segmental or regional identity. However, Hox gene function is also required later in development for the patterning and morphogenesis of limbs and other organs. In Drosophila, spatiotemporal modulation of Sex combs reduced (Scr) expression within the first thoracic (T1) leg underlies the generation of segment- and sex-specific sense organ patterns. High Scr expression in defined domains of the T1 leg is required for the development of T1-specific transverse bristle rows in both sexes and sex combs in males, implying that the patterning of segment-specific sense organs involves incorporation of Scr into the leg development and sex determination gene networks. We sought to gain insight into this process by identifying the cis-and trans-regulatory factors that direct Scr expression during leg development. We have identified two cis-regulatory elements that control spatially modulated Scr expression within T1 legs. One of these enhancers directs sexually dimorphic expression and is required for the formation of T1-specific bristle patterns. We show that the Distalless and Engrailed homeodomain transcription factors act through sequences in this enhancer to establish elevated Scr expression in spatially defined domains. This enhancer functions to integrate Scr into the intrasegmental gene regulatory network, such that Scr serves as a link between leg patterning, sex determination, and sensory organ development.

Disrupting the three-dimensional regulatory topology of the Pitx1 locus results in overtly normal development

Developmental gene expression patterns are orchestrated by thousands of distant-acting transcriptional enhancers. However, identifying enhancers essential for the expression of their target genes has proven challenging. Maps of long-range regulatory interactions may provide the means to identify enhancers crucial for developmental gene expression. To investigate this hypothesis, we used circular chromosome conformation capture coupled with interaction maps in the mouse limb to characterize the regulatory topology of Pitx1, which is essential for hindlimb development. We identified a robust hindlimb-specific interaction between Pitx1 and a putative hindlimb-specific enhancer. To interrogate the role of this interaction in Pitx1 regulation, we used genome editing to delete this enhancer in mouse. Although deletion of the enhancer completely disrupts the interaction, Pitx1 expression in the hindlimb is only mildly affected, without any detectable compensatory interactions between the Pitx1 promoter and potentially redundant enhancers. Pitx1 enhancer null mice did not exhibit any of the characteristic morphological defects of the Pitx1-/- mutant. Our results suggest that robust, tissue-specific physical interactions at essential developmental genes have limited predictive value for identifying enhancer mutations with strong loss-of-function phenotypes.

Enhancer identification and activity evaluation in the red flour beetle, Tribolium castaneum

Evolution of cis-regulatory elements (such as enhancers) plays an important role in the production of diverse morphology. However, a mechanistic understanding is often limited by the absence of methods for studying enhancers in species other than established model systems. Here, we sought to establish methods to identify and test enhancer activity in the red flour beetle, Tribolium castaneum To identify possible enhancer regions, we first obtained genome-wide chromatin profiles from various tissues and stages of Tribolium using FAIRE (formaldehyde-assisted isolation of regulatory elements)-sequencing. Comparison of these profiles revealed a distinct set of open chromatin regions in each tissue and at each stage. In addition, comparison of the FAIRE data with sets of computationally predicted (i.e. supervised cis-regulatory module-predicted) enhancers revealed a very high overlap between the two datasets. Second, using nubbin in the wing and hunchback in the embryo as case studies, we established the first universal reporter assay system that works in various contexts in Tribolium, and in a cross-species context. Together, these advances will facilitate investigation of cis-evolution and morphological diversity in Tribolium and other insects.

Evolutionary emergence of the rac3b/rfng/sgca regulatory cluster refined mechanisms for hindbrain boundaries formation

Developmental programs often rely on parallel morphogenetic mechanisms that guarantee precise tissue architecture. While redundancy constitutes an obvious selective advantage, little is known on how novel morphogenetic mechanisms emerge during evolution. In zebrafish, rhombomeric boundaries behave as an elastic barrier, preventing cell intermingling between adjacent compartments. Here, we identify the fundamental role of the small-GTPase Rac3b in actomyosin cable assembly at hindbrain boundaries. We show that the novel rac3b/rfng/sgca regulatory cluster, which is specifically expressed at the boundaries, emerged in the Ostariophysi superorder by chromosomal rearrangement that generated new cis-regulatory interactions. By combining 4C-seq, ATAC-seq, transgenesis, and CRISPR-induced deletions, we characterized this regulatory domain, identifying hindbrain boundary-specific cis-regulatory elements. Our results suggest that the capacity of boundaries to act as an elastic mesh for segregating rhombomeric cells evolved by cooption of critical genes to a novel regulatory block, refining the mechanisms for hindbrain segmentation.