A cis-regulatory element regulates ERAP2 expression through autoimmune disease risk SNPs.

Single-nucleotide polymorphisms (SNPs) near the ERAP2 gene are associated with various autoimmune conditions, as well as protection against lethal infections. Due to high linkage disequilibrium, numerous trait-associated SNPs are correlated with ERAP2 expression; however, their functional mechanisms remain unidentified. We show by reciprocal allelic replacement that ERAP2 expression is directly controlled by the splice region variant rs2248374. However, disease-associated variants in the downstream LNPEP gene promoter are independently associated with ERAP2 expression. Allele-specific conformation capture assays revealed long-range chromatin contacts between the gene promoters of LNPEP and ERAP2 and showed that interactions were stronger in patients carrying the alleles that increase susceptibility to autoimmune diseases. Replacing the SNPs in the LNPEP promoter by reference sequences lowered ERAP2 expression. These findings show that multiple SNPs act in concert to regulate ERAP2 expression and that disease-associated variants can convert a gene promoter region into a potent enhancer of a distal gene.

Formalin-Fixed, Paraffin-Embedded-Targeted Locus Capture: A Next-Generation Sequencing Technology for Accurate DNA-Based Gene Fusion Detection in Bone and Soft Tissue Tumors.

Chromosomal rearrangements are important drivers in cancer, and their robust detection is essential for diagnosis, prognosis, and treatment selection, particularly for bone and soft tissue tumors. Current diagnostic methods are hindered by limitations, including difficulties with multiplexing targets and poor quality of RNA. A novel targeted DNA-based next-generation sequencing method, formalin-fixed, paraffin-embedded-targeted locus capture (FFPE-TLC), has shown advantages over current diagnostic methods when applied on FFPE lymphomas, including the ability to detect novel rearrangements. We evaluated the utility of FFPE-TLC in bone and soft tissue tumor diagnostics. FFPE-TLC sequencing was successfully applied on noncalcified and decalcified FFPE samples (n = 44) and control samples (n = 19). In total, 58 rearrangements were identified in 40 FFPE tumor samples, including three previously negative samples, and none was identified in the FFPE control samples. In all five discordant cases, FFPE-TLC could identify gene fusions where other methods had failed due to either detection limits or poor sample quality. FFPE-TLC achieved a high specificity and sensitivity (no false positives and negatives). These results indicate that FFPE-TLC is applicable in cancer diagnostics to simultaneously analyze many genes for their involvement in gene fusions. Similar to the observation in lymphomas, FFPE-TLC is a good DNA-based alternative to the conventional methods for detection of rearrangements in bone and soft tissue tumors.

Elevated enhancer-oncogene contacts and higher oncogene expression levels by recurrent CTCF inactivating mutations in acute T cell leukemia.

Monoallelic inactivation of CCCTC-binding factor (CTCF) in human cancer drives altered methylated genomic states, altered CTCF occupancy at promoter and enhancer regions, and deregulated global gene expression. In patients with T cell acute lymphoblastic leukemia (T-ALL), we find that acquired monoallelic CTCF-inactivating events drive subtle and local genomic effects in nearly half of t(5; 14) (q35; q32.2) rearranged patients, especially when CTCF-binding sites are preserved in between the BCL11B enhancer and the TLX3 oncogene. These solitary intervening sites insulate TLX3 from the enhancer by inducing competitive looping to multiple binding sites near the TLX3 promoter. Reduced CTCF levels or deletion of the intervening CTCF site abrogates enhancer insulation by weakening competitive looping while favoring TLX3 promoter to BCL11B enhancer looping, which elevates oncogene expression levels and leukemia burden.

CAF-1 deposits newly synthesized histones during DNA replication using distinct mechanisms on the leading and lagging strands.

During every cell cycle, both the genome and the associated chromatin must be accurately replicated. Chromatin Assembly Factor-1 (CAF-1) is a key regulator of chromatin replication, but how CAF-1 functions in relation to the DNA replication machinery is unknown. Here, we reveal that this crosstalk differs between the leading and lagging strand at replication forks. Using biochemical reconstitutions, we show that DNA and histones promote CAF-1 recruitment to its binding partner PCNA and reveal that two CAF-1 complexes are required for efficient nucleosome assembly under these conditions. Remarkably, in the context of the replisome, CAF-1 competes with the leading strand DNA polymerase epsilon (Polϵ) for PCNA binding. However, CAF-1 does not affect the activity of the lagging strand DNA polymerase Delta (Polδ). Yet, in cells, CAF-1 deposits newly synthesized histones equally on both daughter strands. Thus, on the leading strand, chromatin assembly by CAF-1 cannot occur simultaneously to DNA synthesis, while on the lagging strand these processes may be coupled. We propose that these differences may facilitate distinct parental histone recycling mechanisms and accommodate the inherent asymmetry of DNA replication.

Current challenges in understanding the role of enhancers in disease.

Enhancers play a central role in the spatiotemporal control of gene expression and tend to work in a cell-type-specific manner. In addition, they are suggested to be major contributors to phenotypic variation, evolution and disease. There is growing evidence that enhancer dysfunction due to genetic, structural or epigenetic mechanisms contributes to a broad range of human diseases referred to as enhanceropathies. Such mechanisms often underlie the susceptibility to common diseases, but can also play a direct causal role in cancer or Mendelian diseases. Despite the recent gain of insights into enhancer biology and function, we still have a limited ability to predict how enhancer dysfunction impacts gene expression. Here we discuss the major challenges that need to be overcome when studying the role of enhancers in disease etiology and highlight opportunities and directions for future studies, aiming to disentangle the molecular basis of enhanceropathies.

Building regulatory landscapes reveals that an enhancer can recruit cohesin to create contact domains, engage CTCF sites and activate distant genes.

Developmental gene expression is often controlled by distal regulatory DNA elements called enhancers. Distant enhancer action is restricted to structural chromosomal domains that are flanked by CTCF-associated boundaries and formed through cohesin chromatin loop extrusion. To better understand how enhancers, genes and CTCF boundaries together form structural domains and control expression, we used a bottom-up approach, building series of active regulatory landscapes in inactive chromatin. We demonstrate here that gene transcription levels and activity over time reduce with increased enhancer distance. The enhancer recruits cohesin to stimulate domain formation and engage flanking CTCF sites in loop formation. It requires cohesin exclusively for the activation of distant genes, not of proximal genes, with nearby CTCF boundaries supporting efficient long-range enhancer action. Our work supports a dual activity model for enhancers: its classic role of stimulating transcription initiation and elongation from target gene promoters and a role of recruiting cohesin for the creation of chromosomal domains, the engagement of CTCF sites in chromatin looping and the activation of distal target genes.

Role of the cellular factor CTCF in the regulation of bovine leukemia virus latency and three-dimensional chromatin organization.

Bovine leukemia virus (BLV)-induced tumoral development is a multifactorial phenomenon that remains incompletely understood. Here, we highlight the critical role of the cellular CCCTC-binding factor (CTCF) both in the regulation of BLV transcriptional activities and in the deregulation of the three-dimensional (3D) chromatin architecture surrounding the BLV integration site. We demonstrated the in vivo recruitment of CTCF to three conserved CTCF binding motifs along the provirus. Next, we showed that CTCF localized to regions of transitions in the histone modifications profile along the BLV genome and that it is implicated in the repression of the 5'Long Terminal Repeat (LTR) promoter activity, thereby contributing to viral latency, while favoring the 3'LTR promoter activity. Finally, we demonstrated that BLV integration deregulated the host cellular 3D chromatin organization through the formation of viral/host chromatin loops. Altogether, our results highlight CTCF as a new critical effector of BLV transcriptional regulation and BLV-induced physiopathology.

Interplay between CTCF boundaries and a super enhancer controls cohesin extrusion trajectories and gene expression.

To understand how chromatin domains coordinate gene expression, we dissected select genetic elements organizing topology and transcription around the Prdm14 super enhancer in mouse embryonic stem cells. Taking advantage of allelic polymorphisms, we developed methods to sensitively analyze changes in chromatin topology, gene expression, and protein recruitment. We show that enhancer insulation does not rely strictly on loop formation between its flanking boundaries, that the enhancer activates the Slco5a1 gene beyond its prominent domain boundary, and that it recruits cohesin for loop extrusion. Upon boundary inversion, we find that oppositely oriented CTCF terminates extrusion trajectories but does not stall cohesin, while deleted or mutated CTCF sites allow cohesin to extend its trajectory. Enhancer-mediated gene activation occurs independent of paused loop extrusion near the gene promoter. We expand upon the loop extrusion model to propose that cohesin loading and extrusion trajectories originating at an enhancer contribute to gene activation.

Robust detection of translocations in lymphoma FFPE samples using targeted locus capture-based sequencing.

In routine diagnostic pathology, cancer biopsies are preserved by formalin-fixed, paraffin-embedding (FFPE) procedures for examination of (intra-) cellular morphology. Such procedures inadvertently induce DNA fragmentation, which compromises sequencing-based analyses of chromosomal rearrangements. Yet, rearrangements drive many types of hematolymphoid malignancies and solid tumors, and their manifestation is instructive for diagnosis, prognosis, and treatment. Here, we present FFPE-targeted locus capture (FFPE-TLC) for targeted sequencing of proximity-ligation products formed in FFPE tissue blocks, and PLIER, a computational framework that allows automated identification and characterization of rearrangements involving selected, clinically relevant, loci. FFPE-TLC, blindly applied to 149 lymphoma and control FFPE samples, identifies the known and previously uncharacterized rearrangement partners. It outperforms fluorescence in situ hybridization (FISH) in sensitivity and specificity, and shows clear advantages over standard capture-NGS methods, finding rearrangements involving repetitive sequences which they typically miss. FFPE-TLC is therefore a powerful clinical diagnostics tool for accurate targeted rearrangement detection in FFPE specimens.

Novel orthogonal methods to uncover the complexity and diversity of nuclear architecture.

Recent years have seen a vast expansion of knowledge on three-dimensional (3D) genome organization. The majority of studies on chromosome topology consists of pairwise interaction data of bulk populations of cells and therefore conceals heterogenic and more complex folding patterns. Here, we discuss novel methodologies to study the variation in genome topologies between different cells and techniques that allow analysis of complex, multi-way interactions. These technologies will aid the interpretation of genome-wide chromosome conformation data and provide strategies to further dissect the interplay between genome architecture and transcription regulation.

Genetic Dissection of a Super Enhancer Controlling the Nppa-Nppb Cluster in the Heart.

RATIONALE: ANP (atrial natriuretic peptide) and BNP (B-type natriuretic peptide), encoded by the clustered genes Nppa and Nppb, are important prognostic, diagnostic, and therapeutic proteins in cardiac disease. The spatiotemporal expression pattern and stress-induction of the Nppa and Nppb are tightly regulated, possibly involving their coregulation by an evolutionary conserved enhancer cluster. OBJECTIVE: To explore the physiological functions of the enhancer cluster and elucidate the genomic mechanism underlying Nppa-Nppb coregulation in vivo. METHODS AND RESULTS: By analyzing epigenetic data we uncovered an enhancer cluster with super enhancer characteristics upstream of Nppb. Using CRISPR/Cas9 genome editing, the enhancer cluster or parts thereof, Nppb and flanking regions or the entire genomic block spanning Nppa-Nppb, respectively, were deleted from the mouse genome. The impact on gene regulation and phenotype of the respective mouse lines was investigated by transcriptomic, epigenomic, and phenotypic analyses. The enhancer cluster was essential for prenatal and postnatal ventricular expression of Nppa and Nppb but not of any other gene. Enhancer cluster-deficient mice showed enlarged hearts before and after birth, similar to Nppa-Nppb compound knockout mice we generated. Analysis of the other deletion alleles indicated the enhancer cluster engages the promoters of Nppa and Nppb in a competitive rather than a cooperative mode, resulting in increased Nppa expression when Nppb and flanking sequences were deleted. The enhancer cluster maintained its active epigenetic state and selectivity when its target genes are absent. In enhancer cluster-deficient animals, Nppa was induced but remained low in the postmyocardial infarction border zone and in the hypertrophic ventricle, involving regulatory sequences proximal to Nppa. CONCLUSIONS: Coordinated ventricular expression of Nppa and Nppb is controlled in a competitive manner by a shared super enhancer, which is also required to augment stress-induced expression and to prevent premature hypertrophy.

Epigenetic Analyses of Human Left Atrial Tissue Identifies Gene Networks Underlying Atrial Fibrillation.

BACKGROUND: Atrial fibrillation (AF) often arises from structural abnormalities in the left atria (LA). Annotation of the noncoding genome in human LA is limited, as are effects on gene expression and chromatin architecture. Many AF-associated genetic variants reside in noncoding regions; this knowledge gap impairs efforts to understand the molecular mechanisms of AF and cardiac conduction phenotypes. METHODS: We generated a model of the LA noncoding genome by profiling 7 histone post-translational modifications (active: H3K4me3, H3K4me2, H3K4me1, H3K27ac, H3K36me3; repressive: H3K27me3, H3K9me3), CTCF binding, and gene expression in samples from 5 individuals without structural heart disease or AF. We used MACS2 to identify peak regions (P<0.01), applied a Markov model to classify regulatory elements, and annotated this model with matched gene expression data. We intersected chromatin states with expression quantitative trait locus, DNA methylation, and HiC chromatin interaction data from LA and left ventricle. Finally, we integrated genome-wide association data for AF and electrocardiographic traits to link disease-related variants to genes. RESULTS: Our model identified 21 epigenetic states, encompassing regulatory motifs, such as promoters, enhancers, and repressed regions. Genes were regulated by proximal chromatin states; repressive states were associated with a significant reduction in gene expression (P<2×10-16). Chromatin states were differentially methylated, promoters were less methylated than repressed regions (P<2×10-16). We identified over 15 000 LA-specific enhancers, defined by homeobox family motifs, and annotated several cardiovascular disease susceptibility loci. Intersecting AF and PR genome-wide association studies loci with long-range chromatin conformation data identified a gene interaction network dominated by NKX2-5, TBX3, ZFHX3, and SYNPO2L. CONCLUSIONS: Profiling the noncoding genome provides new insights into the gene expression and chromatin regulation in human LA tissue. These findings enabled identification of a gene network underlying AF; our experimental and analytic approach can be extended to identify molecular mechanisms for other cardiac diseases and traits.

How chromosome topologies get their shape: views from proximity ligation and microscopy methods.

The 3D organization of our genome is an important determinant for the transcriptional output of a gene in (patho)physiological contexts. The spatial organization of linear chromosomes within nucleus is dominantly inferred using two distinct approaches, chromosome conformation capture (3C) and DNA fluorescent in situ hybridization (DNA-FISH). While 3C and its derivatives score genomic interaction frequencies based on proximity ligation events, DNA-FISH methods measure physical distances between genomic loci. Despite these approaches probe different characteristics of chromosomal topologies, they provide a coherent picture of how chromosomes are organized in higher-order structures encompassing chromosome territories, compartments, and topologically associating domains. Yet, at the finer topological level of promoter-enhancer communication, the imaging-centered and the 3C methods give more divergent and sometimes seemingly paradoxical results. Here, we compare and contrast observations made applying visual DNA-FISH and molecular 3C approaches. We emphasize that the 3C approach, due to its inherently competitive ligation step, measures only 'relative' proximities. A 3C interaction enriched between loci, therefore does not necessarily translates into a decrease in absolute spatial distance. Hence, we advocate caution when modeling chromosome conformations.

Epigenetic and Transcriptional Networks Underlying Atrial Fibrillation.

Genome-wide association studies have uncovered over a 100 genetic loci associated with atrial fibrillation (AF), the most common arrhythmia. Many of the top AF-associated loci harbor key cardiac transcription factors, including PITX2, TBX5, PRRX1, and ZFHX3. Moreover, the vast majority of the AF-associated variants lie within noncoding regions of the genome where causal variants affect gene expression by altering the activity of transcription factors and the epigenetic state of chromatin. In this review, we discuss a transcriptional regulatory network model for AF defined by effector genes in Genome-wide association studies loci. We describe the current state of the field regarding the identification and function of AF-relevant gene regulatory networks, including variant regulatory elements, dose-sensitive transcription factor functionality, target genes, and epigenetic states. We illustrate how altered transcriptional networks may impact cardiomyocyte function and ionic currents that impact AF risk. Last, we identify the need for improved tools to identify and functionally test transcriptional components to define the links between genetic variation, epigenetic gene regulation, and atrial function.

Recently Evolved Enhancers Emerge with High Interindividual Variability and Less Frequently Associate with Disease.

Mutations in non-coding regulatory DNA such as enhancers underlie a wide variety of diseases including developmental disorders and cancer. As enhancers rapidly evolve, understanding their function and configuration in non-human disease models can have important clinical applications. Here, we analyze enhancer configurations in tissues isolated from the common marmoset, a widely used primate model for human disease. Integrating these data with human and mouse data, we find that enhancers containing trait-associated variants are preferentially conserved. In contrast, most human-specific enhancers are highly variable between individuals, with a subset failing to contact promoters. These are located further away from genes and more often reside in inactive B-compartments. Our data show that enhancers typically emerge as instable elements with minimal biological impact prior to their integration in a transcriptional program. Furthermore, our data provide insight into which trait variations in enhancers can be faithfully modeled using the common marmoset.

Hominin-specific regulatory elements selectively emerged in oligodendrocytes and are disrupted in autism patients.

Speciation is associated with substantial rewiring of the regulatory circuitry underlying the expression of genes. Determining which changes are relevant and underlie the emergence of the human brain or its unique susceptibility to neural disease has been challenging. Here we annotate changes to gene regulatory elements (GREs) at cell type resolution in the brains of multiple primate species spanning most of primate evolution. We identify a unique set of regulatory elements that emerged in hominins prior to the separation of humans and chimpanzees. We demonstrate that these hominin gains perferentially affect oligodendrocyte function postnatally and are preferentially affected in the brains of autism patients. This preference is also observed for human-specific GREs suggesting this system is under continued selective pressure. Our data provide a roadmap of regulatory rewiring across primate evolution providing insight into the genomic changes that underlie the emergence of the brain and its susceptibility to neural disease.

Multi-contact 4C: long-molecule sequencing of complex proximity ligation products to uncover local cooperative and competitive chromatin topologies.

We present the experimental protocol and data analysis toolbox for multi-contact 4C (MC-4C), a new proximity ligation method tailored to study the higher-order chromatin contact patterns of selected genomic sites. Conventional chromatin conformation capture (3C) methods fragment proximity ligation products for efficient analysis of pairwise DNA contacts. By contrast, MC-4C is designed to preserve and collect large concatemers of proximity ligated fragments for long-molecule sequencing on an Oxford Nanopore or Pacific Biosciences platform. Each concatemer of proximity ligation products represents a snapshot topology of a different individual allele, revealing its multi-way chromatin interactions. By inverse PCR with primers specific for a fragment of interest (the viewpoint) and DNA size selection, sequencing is selectively targeted to thousands of different complex interactions containing this viewpoint. A tailored statistical analysis toolbox is able to generate background models and three-way interaction profiles from the same dataset. These profiles can be used to distinguish whether contacts between more than two regulatory sequences are mutually exclusive or, conversely, simultaneously occurring at chromatin hubs. The entire procedure can be completed in 2 w, and requires standard molecular biology and data analysis skills and equipment, plus access to a third-generation sequencing platform.

4C-seq from beginning to end: A detailed protocol for sample preparation and data analysis.

Chromosome conformation capture (3C) methods measure DNA contact frequencies based on nuclear proximity ligation, to uncover in vivo genomic folding patterns. 4C-seq is a derivative 3C method, designed to search the genome for sequences contacting a selected genomic site of interest. 4C-seq employs inverse PCR and next generation sequencing to amplify, identify and quantify its proximity ligated DNA fragments. It generates high-resolution contact profiles for selected genomic sites based on limited amounts of sequencing reads. 4C-seq can be used to study multiple aspects of genome organization. It primarily serves to identify specific long-range DNA contacts between individual regulatory DNA modules, forming for example regulatory chromatin loops between enhancers and promoters, or architectural chromatin loops between cohesin- and CTCF- associated domain boundaries. Additionally, 4C-seq contact profiles can reveal the contours of contact domains and can identify the structural domains that co-occupy the same nuclear compartment. Here, we present an improved step-by-step protocol for sample preparation and the generation of 4C-seq sequencing libraries, including an optimized PCR and 4C template purification strategy. In addition, a data processing pipeline is provided which processes multiplexed 4C-seq reads directly from FASTQ files and generates files compatible with standard genome browsers for visualization and further statistical analysis of the data such as peak calling using peakC. The protocols and the pipeline presented should readily allow anyone to generate, visualize and interpret their own high resolution 4C contact datasets.

Identification and Characterization of a Transcribed Distal Enhancer Involved in Cardiac Kcnh2 Regulation.

The human ether-a-go-go-related gene KCNH2 encodes the voltage-gated potassium channel underlying IKr, a current critical for the repolarization phase of the cardiac action potential. Mutations in KCNH2 that cause a reduction of the repolarizing current can result in cardiac arrhythmias associated with long-QT syndrome. Here, we investigate the regulation of KCNH2 and identify multiple active enhancers. A transcribed enhancer ∼85 kbp downstream of Kcnh2 physically contacts the promoters of two Kcnh2 isoforms in a cardiac-specific manner in vivo. Knockdown of its ncRNA transcript results in reduced expression of Kcnh2b and two neighboring mRNAs, Nos3 and Abcb8, in vitro. Genomic deletion of the enhancer, including the ncRNA transcription start site, from the mouse genome causes a modest downregulation of both Kcnh2a and Kcnh2b in the ventricles. These findings establish that the regulation of Kcnh2a and Kcnh2b is governed by a complex regulatory landscape that involves multiple partially redundantly acting enhancers.