Pathogenic variations in the sodium voltage-gated channel alpha subunit 1 (SCN1A) gene are responsible for multiple epilepsy phenotypes, including Dravet syndrome, febrile seizures (FS) and genetic epilepsy with FS plus. Phenotypic heterogeneity is a hallmark of SCN1A-related epilepsies, the causes of which are yet to be clarified. Genetic variation in the non-coding regulatory regions of SCN1A could be one potential causal factor. However, a comprehensive understanding of the SCN1A regulatory landscape is currently lacking. Here, we summarized the current state of knowledge of SCN1A regulation, providing details on its promoter and enhancer regions. We then integrated currently available data on SCN1A promoters by extracting information related to the SCN1A locus from genome-wide repositories and clearly defined the promoter and enhancer regions of SCN1A. Further, we explored the cellular specificity of differential SCN1A promoter usage. We also reviewed and integrated the available human brain-derived enhancer databases and mouse-derived data to provide a comprehensive computationally developed summary of SCN1A brain-active enhancers. By querying genome-wide data repositories, extracting SCN1A-specific data and integrating the different types of independent evidence, we created a comprehensive catalogue that better defines the regulatory landscape of SCN1A, which could be used to explore the role of SCN1A regulatory regions in disease.

Genetic changes that modify the function of transcriptional enhancers have been linked to the evolution of biological diversity across species. Multiple studies have focused on the role of nucleotide substitutions, transposition, and insertions and deletions in altering enhancer function. Here we show that turnover of CpG islands (CGIs), which contribute to enhancer activation, is broadly associated with changes in enhancer activity across mammals, including humans. We integrated maps of CGIs and enhancer activity-associated histone modifications obtained from multiple tissues in nine mammalian species and found that CGI content in enhancers was strongly associated with increased histone modification levels. CGIs showed widespread turnover across species and species-specific CGIs were strongly enriched for enhancers exhibiting species-specific activity across all tissues and species we examined. Genes associated with enhancers with species-specific CGIs showed concordant biases in their expression, supporting that CGI turnover contributes to gene regulatory innovation. Our results also implicate CGI turnover in the evolution of Human Gain Enhancers (HGEs), which show increased activity in human embryonic development and may have contributed to the evolution of uniquely human traits. Using a humanized mouse model, we show that a highly conserved HGE with a large CGI absent from the mouse ortholog shows increased activity at the human CGI in the humanized mouse diencephalon. Collectively, our results point to CGI turnover as a mechanism driving gene regulatory changes potentially underlying trait evolution in mammals.

Hox genes encode transcription factors whose roles in patterning animal body plans during embryonic development are well-documented. Multiple studies demonstrate that Hox genes continue to act in adult cells, in normal differentiation, in regenerative processes, and, with abnormal expression, in diverse types of cancers. However, surprisingly little is known about the regulatory mechanisms that govern Hox gene expression in specific cell types, as they differentiate during late embryonic development, and in the adult organism. The murine Hoxc8 gene determines the identity of multiple skeletal elements in the lower thoracic and lumbar region and continues to play a role in the proliferation and differentiation of cells in cartilage as the skeleton matures. This study was undertaken to identify regulatory elements in the Hoxc8 gene that control transcriptional activity, specifically in cartilage-producing chondrocytes. We report that an enhancer comprising two 416 and 224 bps long interacting DNA elements produces reporter gene activity when assayed on a heterologous transcriptional promoter in transgenic mice. This enhancer is distinct in spatial, temporal, and molecular regulation from previously identified regulatory sequences in the Hoxc8 gene that control its expression in early development. The identification of a tissue-specific Hox gene regulatory element now allows mechanistic investigations into Hox transcription factor expression and function in differentiating cell types and adult tissues and to specifically target these cells during repair processes and regeneration.

Eukaryotic gene expression regulation involves thousands of distal regulatory elements. Understanding the quantitative contribution of individual enhancers to gene expression is critical for assessing the role of disease-associated genetic risk variants. Yet, we lack the ability to accurately link genes with their distal regulatory elements. To address this, we used 3D enhancer-promoter (E-P) associations identified using split-pool recognition of interactions by tag extension (SPRITE) to build a predictive model of gene expression. Our model dramatically outperforms models using genomic proximity and can be used to determine the quantitative impact of enhancer loss on gene expression in different genetic backgrounds. We show that genes that form stable E-P hubs have less cell-to-cell variability in gene expression. Finally, we identified transcription factors that regulate stimulation-dependent E-P interactions. Together, our results provide a framework for understanding quantitative contributions of E-P interactions and associated genetic variants to gene expression.

Mitotically stable random monoallelic gene expression (RME) is documented for a small percentage of autosomal genes. We developed an in vivo genetic model to study the role of enhancers in RME using high-resolution single-cell analysis of natural killer (NK) cell receptor gene expression and enhancer deletions in the mouse germline. Enhancers of the RME NK receptor genes were accessible and enriched in H3K27ac on silent and active alleles alike in cells sorted according to allelic expression status, suggesting enhancer activation and gene expression status can be decoupled. In genes with multiple enhancers, enhancer deletion reduced gene expression frequency, in one instance converting the universally expressed gene encoding NKG2D into an RME gene, recapitulating all aspects of natural RME including mitotic stability of both the active and silent states. The results support the binary model of enhancer action, and suggest that RME is a consequence of general properties of gene regulation by enhancers rather than an RME-specific epigenetic program. Therefore, many and perhaps all genes may be subject to some degree of RME. Surprisingly, this was borne out by analysis of several genes that define different major hematopoietic lineages, that were previously thought to be universally expressed within those lineages: the genes encoding NKG2D, CD45, CD8α, and Thy-1. We propose that intrinsically probabilistic gene allele regulation is a general property of enhancer-controlled gene expression, with previously documented RME representing an extreme on a broad continuum.

Clinical studies in type 2 diabetes (T2D) primarily focused on the single nucleotide polymorphisms (SNPs) located in protein-coding regions. Recently, the SNPs located in noncoding regions have also been recognized to play an important role in disease susceptibility. The super enhancer is a cluster of transcriptional enhancers located in noncoding regions. It plays a critical role in cell-type specific gene expression. However, the exact mechanism of the super enhancer SNPs for T2D remains unclear. In this study, we integrated genome-wide association studies (GWASs) and T2D cell/tissue-specific histone modification ChIP-seq data to identify T2D-associated SNPs in super enhancer, followed by comprehensive bioinformatics analyses to further explore the functional importance of these SNPs. We identified several interesting T2D super enhancer SNPs. Interesting, most of them were clustered within the same or neighboring super enhancers. A number of SNPs are involved in chromatin interactive regulation and/or potentially influence the binding affinity of transcription factors. Gene Ontology (GO) analysis showed a significant enrichment in several well-known signaling pathways and regulatory process, e.g. WNT signaling pathway, which plays a key role in T2D metabolism. Our results highlighted the potential functional importance of T2D super enhancer SNPs, which may yield novel insights into the pathogenesis of T2D.

Hemangioblasts are known as the common precursors for primitive hematopoietic and endothelial lineages. Their existence has been supported mainly by the observation that both cell types develop in close proximity and by in vitro differentiation and genetic studies. However, more compelling evidence will arise from tracking their cell fates using a lineage-specific marker.

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the 4D arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find that alterations in chromatin mobility, not promoter-enhancer distance, is more informative about transcriptional status. Active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

Identifying transcriptional enhancers and their target genes is essential for understanding gene regulation and the impact of human genetic variation on disease1-6. Here we create and evaluate a resource of >13 million enhancer-gene regulatory interactions across 352 cell types and tissues, by integrating predictive models, measurements of chromatin state and 3D contacts, and largescale genetic perturbations generated by the ENCODE Consortium7. We first create a systematic benchmarking pipeline to compare predictive models, assembling a dataset of 10,411 elementgene pairs measured in CRISPR perturbation experiments, >30,000 fine-mapped eQTLs, and 569 fine-mapped GWAS variants linked to a likely causal gene. Using this framework, we develop a new predictive model, ENCODE-rE2G, that achieves state-of-the-art performance across multiple prediction tasks, demonstrating a strategy involving iterative perturbations and supervised machine learning to build increasingly accurate predictive models of enhancer regulation. Using the ENCODE-rE2G model, we build an encyclopedia of enhancer-gene regulatory interactions in the human genome, which reveals global properties of enhancer networks, identifies differences in the functions of genes that have more or less complex regulatory landscapes, and improves analyses to link noncoding variants to target genes and cell types for common, complex diseases. By interpreting the model, we find evidence that, beyond enhancer activity and 3D enhancer-promoter contacts, additional features guide enhancerpromoter communication including promoter class and enhancer-enhancer synergy. Altogether, these genome-wide maps of enhancer-gene regulatory interactions, benchmarking software, predictive models, and insights about enhancer function provide a valuable resource for future studies of gene regulation and human genetics.

The CCCTC-binding factor (CTCF) works together with the cohesin complex to drive the formation of chromatin loops and topologically associating domains, but its role in gene regulation has not been fully defined. Here, we investigated the effects of acute CTCF loss on chromatin architecture and transcriptional programs in mouse embryonic stem cells undergoing differentiation to neural precursor cells. We identified CTCF-dependent enhancer-promoter contacts genome-wide and found that they disproportionately affect genes that are bound by CTCF at the promoter and are dependent on long-distance enhancers. Disruption of promoter-proximal CTCF binding reduced both long-range enhancer-promoter contacts and transcription, which were restored by artificial tethering of CTCF to the promoter. Promoter-proximal CTCF binding is correlated with the transcription of over 2,000 genes across a diverse set of adult tissues. Taken together, the results of our study show that CTCF binding to promoters may promote long-distance enhancer-dependent transcription at specific genes in diverse cell types.

Enhancers are genetic elements that regulate spatiotemporal gene expression. Enhancer function requires transcription factor (TF) binding and correlates with histone modifications. However, the extent to which TF binding and histone modifications functionally define active enhancers remains unclear. Here, we combine chromatin immunoprecipitation with a massively parallel reporter assay (ChIP-STARR-seq) to identify functional enhancers in human embryonic stem cells (ESCs) genome-wide in a quantitative unbiased manner. Although active enhancers associate with TFs, only a minority of regions marked by NANOG, OCT4, H3K27ac, and H3K4me1 function as enhancers, with activity markedly changing under naive versus primed culture conditions. We identify an enhancer set associated with functions extending to non-ESC-specific processes. Moreover, although transposable elements associate with putative enhancers, only some exhibit activity. Similarly, within super-enhancers, large tracts are non-functional, with activity restricted to small sub-domains. This catalog of validated enhancers provides a valuable resource for further functional dissection of the regulatory genome.

Though numerous polymorphisms have been associated with risk of developing lymphoma, how these variants function to promote tumorigenesis is poorly understood. Here, we report that lymphoma risk SNPs, especially in the non-Hodgkin's lymphoma subtype chronic lymphocytic leukemia, are significantly enriched for co-localization with epigenetic marks of active gene regulation. These enrichments were seen in a lymphoid-specific manner for numerous ENCODE datasets, including DNase-hypersensitivity as well as multiple segmentation-defined enhancer regions. Furthermore, we identify putatively functional SNPs that are both in regulatory elements in lymphocytes and are associated with gene expression changes in blood. We developed an algorithm, UES, that uses a Monte Carlo simulation approach to calculate the enrichment of previously identified risk SNPs in various functional elements. This multiscale approach integrating multiple datasets helps disentangle the underlying biology of lymphoma, and more broadly, is generally applicable to GWAS results from other diseases as well.

Differences in cells' functions arise from differential activity of regulatory elements, including enhancers. Enhancers are cis-regulatory elements that cooperate with promoters through transcription factors to activate the expression of one or several genes by getting physically close to them in the 3D space of the nucleus. There is increasing evidence that genetic variants associated with common diseases are enriched in enhancers active in cell types relevant to these diseases. Identifying the enhancers associated with genes and conversely, the sets of genes activated by each enhancer (the so-called enhancer/gene or E/G relationships) across cell types, can help understanding the genetic mechanisms underlying human diseases. There are three broad approaches for the genome-wide identification of E/G relationships in a cell type: 1) genetic link methods or eQTL, 2) functional link methods based on 1D functional data such as open chromatin, histone mark or gene expression and 3) spatial link methods based on 3D data such as HiC. Since 1) and 3) are costly, the current strategy is to develop functional link methods and to use data from 1) and 3) as reference to evaluate them. However, there is still no consensus on the best functional link method to date, and method comparison remain seldom. Here, we compared the relative performances of three recent methods for the identification of enhancer-gene links, TargetFinder, Average-Rank, and the ABC model, using the three latest benchmarks from the field: a reference that combines 3D and eQTL data, called BENGI, and two genetic screening references, called CRiFF and CRiSPRi. Overall, none of the three methods performed best on the three references. CRiFF and CRISPRi reference sets are likely more reliable, but CRiFF is not genome-wide and CRiFF and CRISPRi are mostly available on the K562 cancer cell line. The BENGI reference set is genome-wide but likely contains many false positives. This study therefore calls for new reliable and genome-wide E/G reference data rather than new functional link E/G identification methods.

The spatiotemporal configuration of genes with distal regulatory elements, and the impact of chromatin mobility on transcription, remain unclear. Loop extrusion is an attractive model for bringing genetic elements together, but how this functionally interacts with transcription is also largely unknown. We combine live tracking of genomic loci and nascent transcripts with molecular dynamics simulations to assess the spatiotemporal arrangement of the Sox2 gene and its enhancer, in response to a battery of perturbations. We find a close link between chromatin mobility and transcriptional status: active elements display more constrained mobility, consistent with confinement within specialized nuclear sites, and alterations in enhancer mobility distinguish poised from transcribing alleles. Strikingly, we find that whereas loop extrusion and transcription factor-mediated clustering contribute to promoter-enhancer proximity, they have antagonistic effects on chromatin dynamics. This provides an experimental framework for the underappreciated role of chromatin dynamics in genome regulation.

Several heart malformations are associated with mutations in the regulatory regions of cardiac genes. Troponin I type 1b (tnni1b) is important for the formation of the atrioventricular canal in zebrafish hearts; however, the regulation of tnni1b is poorly understand. We aimed to identify a small but functional enhancer that is distal to tnni1b.

Genetic changes that altered the function of gene regulatory elements have been implicated in the evolution of human traits such as the expansion of the cerebral cortex. However, identifying the particular changes that modified regulatory activity during human evolution remain challenging. Here we used massively parallel enhancer assays in neural stem cells to quantify the functional impact of >32,000 human-specific substitutions in >4,300 human accelerated regions (HARs) and human gain enhancers (HGEs), which include enhancers with novel activities in humans. We found that >30% of active HARs and HGEs exhibited differential activity between human and chimpanzee. We isolated the effects of human-specific substitutions from background genetic variation to identify the effects of genetic changes most relevant to human evolution. We found that substitutions interacted in both additive and nonadditive ways to modify enhancer function. Substitutions within HARs, which are highly constrained compared to HGEs, showed smaller effects on enhancer activity, suggesting that the impact of human-specific substitutions is buffered in enhancers with constrained ancestral functions. Our findings yield insight into how human-specific genetic changes altered enhancer function and provide a rich set of candidates for studies of regulatory evolution in humans.

Enhancer trapping (ET) is a powerful approach to establish tissue- or cell-specific reporters and identify expression patterns of uncharacterized genes. Although a number of enhancer-trapping vectors have been developed and a large library of fish lines with distinct tissue- or cell-specific expression of reporter genes have been generated, the specificity and efficiency of trapping vectors need to be improved because of the bias interaction of minimal promoters with genomic enhancers. Accordingly, we generated an enhancer-trapping vector pTME that contains a minimal mouse metallothionein gene (mMTI) promoter upstream of EGFP reporter. In the first round of screening, twelve zebrafish lines that carry a single copy of ET cassettes were characterized to have tissue- or cell-specific EGFP expression. One of the highly conserved noncoding elements near an insertion site of trapping cassettes was characterized as an enhancer that can specifically regulate the expression of EGFP in cells of the central nervous system. In addition, the pTME vector contains a mutation-cassette that is able to effectively block the transcription of an endogenous gene in an ET line with ubiquitous EGFP expression. Thus, the pTME vector can be used as an alternative tool for both enhancer trapping and mutagenesis across a target genome.

Xenopus provides a simple and efficient model system to study nephrogenesis and explore the mechanisms causing renal developmental defects in human. Hnf1b (hepatocyte nuclear factor 1 homeobox b), a gene whose mutations are the most commonly identified genetic cause of developmental kidney disease, is required for the acquisition of a proximo-intermediate nephron segment in Xenopus as well as in mouse. Genetic networks involved in Hnf1b expression during kidney development remain poorly understood. We decided to explore the transcriptional regulation of Hnf1b in the developing Xenopus pronephros and mammalian renal cells. Using phylogenetic footprinting, we identified an evolutionary conserved sequence (CNS1) located several kilobases (kb) upstream the Hnf1b transcription start and harboring epigenomic marks characteristics of a distal enhancer in embryonic and adult renal cells in mammals. By means of functional expression assays in Xenopus and mammalian renal cell lines we showed that CNS1 displays enhancer activity in renal tissue. Using CRISPR/cas9 editing in Xenopus tropicalis, we demonstrated the in vivo functional relevance of CNS1 in driving hnf1b expression in the pronephros. We further showed the importance of Pax8-CNS1 interaction for CNS1 enhancer activity allowing us to conclude that Hnf1b is a direct target of Pax8. Our work identified for the first time a Hnf1b renal specific enhancer and may open important perspectives into the diagnosis for congenital kidney anomalies in human, as well as modeling HNF1B-related diseases.

Proinflammatory stimuli elicit rapid transcriptional responses via transduced signals to master regulatory transcription factors. To explore the role of chromatin-dependent signal transduction in the atherogenic inflammatory response, we characterized the dynamics, structure, and function of regulatory elements in the activated endothelial cell epigenome. Stimulation with tumor necrosis factor alpha prompted a dramatic and rapid global redistribution of chromatin activators to massive de novo clustered enhancer domains. Inflammatory super enhancers formed by nuclear factor-kappa B accumulate at the expense of immediately decommissioned, basal endothelial super enhancers, despite persistent histone hyperacetylation. Mass action of enhancer factor redistribution causes momentous swings in transcriptional initiation and elongation. A chemical genetic approach reveals a requirement for BET bromodomains in communicating enhancer remodeling to RNA Polymerase II and orchestrating the transition to the inflammatory cell state, demonstrated in activated endothelium and macrophages. BET bromodomain inhibition abrogates super enhancer-mediated inflammatory transcription, atherogenic endothelial responses, and atherosclerosis in vivo.

The bony skeleton is continuously renewed throughout adult life by the bone remodeling process, in which old or damaged bone is removed by osteoclasts via largely unknown mechanisms. Osteocytes regulate bone remodeling by producing the osteoclast differentiation factor RANKL (encoded by the TNFSF11 gene). However, the precise mechanisms underlying RANKL expression in osteocytes are still elusive. Here, we explored the epigenomic landscape of osteocytic cells and identified a hitherto-undescribed osteocytic cell-specific intronic enhancer in the TNFSF11 gene locus. Bioinformatics analyses showed that transcription factors involved in cell death and senescence act on this intronic enhancer region. Single-cell transcriptomic data analysis demonstrated that cell death signaling increased RANKL expression in osteocytic cells. Genetic deletion of the intronic enhancer led to a high-bone-mass phenotype with decreased levels of RANKL in osteocytic cells and osteoclastogenesis in the adult stage, while RANKL expression was not affected in osteoblasts or lymphocytes. These data suggest that osteocytes may utilize a specialized regulatory element to facilitate osteoclast formation at the bone surface to be resorbed by linking signals from cellular senescence/death and RANKL expression.