The impact of genomic distance on enhancer-promoter interactions at the CFTR locus.

We identified and characterized multiple cell-type selective enhancers of the CFTR gene promoter in previous work and demonstrated active looping of these elements to the promoter. Here we address the impact of genomic spacing on these enhancer:promoter interactions and on CFTR gene expression. Using CRISPR/Cas9, we generated clonal cell lines with deletions between the -35 kb airway enhancer and the CFTR promoter in the 16HBE14o- airway cell line, or between the intron 1 (185 + 10 kb) intestinal enhancer and the promoter in the Caco2 intestinal cell line. The effect of these deletions on CFTR transcript abundance, as well as the 3D looping structure of the locus was investigated in triplicate clones of each modification. Our results indicate that both small and larger deletions upstream of the promoter can perturb CFTR expression and -35 kb enhancer:promoter interactions in the airway cells, though the larger deletions are more impactful. In contrast, the small intronic deletions have no effect on CFTR expression and intron 1 enhancer:promoter interactions in the intestinal cells, whereas larger deletions do. Clonal variation following a specific CFTR modification is a confounding factor particularly in 16HBE14o- cells.

4C Analysis of EBV-Host DNA Interactome.

EBV persist as multicopy episomes in latently infected cells and alter transcriptional program of host systems. Knowledge of EBV tethering site helps us understand how EBV attaches to and regulates the host chromosome. Here, we introduce a step-by-step protocol for 4C-seq analysis, including cell fixation, 4C-DNA construction, and sequencing library preparation performed with EBV-positive Burkitt's lymphoma cells. The method can be applied in a variety of studies and cell-types to identify target loci associated with bait positions, such as viral episomes.

The pZRS non-coding regulatory mutation resulting in triphalangeal thumb-polysyndactyly syndrome changes the pattern of local interactions.

Herein, we report on a large Polish family presenting with a classical triphalangeal thumb-polysyndactyly syndrome (TPT-PS). This rare congenital limb anomaly is generally caused by microduplications encompassing the Sonic Hedgehog (SHH) limb enhancer, termed the zone of polarizing activity (ZPA) regulatory sequence (ZRS). Recently, a pathogenic variant in the pre-ZRS (pZRS), a conserved sequence located near the ZRS, has been described in a TPT-PS Dutch family. We performed targeted ZRS sequencing, array comparative genomic hybridization, and whole-exome sequencing. Next, we sequenced the recently described pZRS region. Finally, we performed a circular chromatin conformation capture-sequencing (4C-seq) assay on skin fibroblasts of one affected family member and control samples to examine potential alterations in the SHH regulatory domain and functionally characterize the identified variant. We found that all affected individuals shared a recently identified pathogenic point mutation in the pZRS region: NC_000007.14:g.156792782C>G (GRCh38/hg38), which is the same as in the Dutch family. The results of 4C-seq experiments revealed increased interactions within the whole SHH regulatory domain (SHH-LMBR1 TAD) in the patient compared to controls. Our study expands the number of TPT-PS families carrying a pathogenic alteration of the pZRS and underlines the importance of routine pZRS sequencing in the genetic diagnostics of patients with TPT-PS or similar phenotypes. The pathogenic mutation causative for TPT-PS in our patient gave rise to increased interactions within the SHH regulatory domain in yet unknown mechanism.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle.

BACKGROUND: Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined. RESULTS: In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers. CONCLUSIONS: This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Exploring high-resolution chromatin interaction changes and functional enhancers of myogenic marker genes during myogenic differentiation.

Skeletal muscle differentiation (myogenesis) is a complex and highly coordinated biological process regulated by a series of myogenic marker genes. Chromatin interactions between gene's promoters and their enhancers have an important role in transcriptional control. However, the high-resolution chromatin interactions of myogenic genes and their functional enhancers during myogenesis remain largely unclear. Here, we used circularized chromosome conformation capture coupled with next generation sequencing (4C-seq) to investigate eight myogenic marker genes in C2C12 myoblasts (C2C12-MBs) and C2C12 myotubes (C2C12-MTs). We revealed dynamic chromatin interactions of these marker genes during differentiation and identified 163 and 314 significant interaction sites (SISs) in C2C12-MBs and C2C12-MTs, respectively. The interacting genes of SISs in C2C12-MTs were mainly involved in muscle development, and histone modifications of the SISs changed during differentiation. Through functional genomic screening, we also identified 25 and 41 putative active enhancers in C2C12-MBs and C2C12-MTs, respectively. Using luciferase reporter assays for putative enhancers of Myog and Myh3, we identified eight activating enhancers. Furthermore, dCas9-KRAB epigenome editing and RNA-Seq revealed a role for Myog enhancers in the regulation of Myog expression and myogenic differentiation in the native genomic context. Taken together, this study lays the groundwork for understanding 3D chromatin interaction changes of myogenic genes during myogenesis and provides insights that contribute to our understanding of the role of enhancers in regulating myogenesis.

Circular Chromosome Conformation Capture Sequencing (4C-Seq ) in Primary Adherent Cells.

The three-dimensional structure of the genome is highly organized and is an important aspect of gene regulation. Chromatin interactions can be identified using chromosome conformation capture-based techniques, which rely on proximity ligation. Of these techniques, circular chromosome conformation capture sequencing (4C-seq) is used to identify all chromatin interactions occurring with a single chromosomal location (one versus all). Here we describe a 4C-seq protocol that has been optimized for primary adherent cells, for which the first digestion step is inefficient using standard 4C-seq protocols. It can, however, also be applied to other cell or tissue types. This protocol utilizes a standard DNA library preparation method using a commercial kit, and includes a description of the data processing steps.

Multimodal regulatory elements within a hormone-specific super enhancer control a heterogeneous transcriptional response.

The hormone-stimulated glucocorticoid receptor (GR) modulates transcription by interacting with thousands of enhancers and GR binding sites (GBSs) throughout the genome. Here, we examined the effects of GR binding on enhancer dynamics and investigated the contributions of individual GBSs to the hormone response. Hormone treatment resulted in genome-wide reorganization of the enhancer landscape in breast cancer cells. Upstream of the DDIT4 oncogene, GR bound to four sites constituting a hormone-dependent super enhancer. Three GBSs were required as hormone-dependent enhancers that differentially promoted histone acetylation, transcription frequency, and burst size. Conversely, the fourth site suppressed transcription and hormone treatment alleviated this suppression. GR binding within the super enhancer promoted a loop-switching mechanism that allowed interaction of the DDIT4 TSS with the active GBSs. The unique functions of each GR binding site contribute to hormone-induced transcriptional heterogeneity and demonstrate the potential for targeted modulation of oncogene expression.

Genetics and Genomics of SOST: Functional Analysis of Variants and Genomic Regulation in Osteoblasts.

SOST encodes the sclerostin protein, which acts as a key extracellular inhibitor of the canonical Wnt pathway in bone, playing a crucial role in skeletal development and bone homeostasis. The objective of this work was to assess the functionality of two variants previously identified (the rare variant rs570754792 and the missense variant p.Val10Ile) and to investigate the physical interactors of the SOST proximal promoter region in bone cells. Through a promoter luciferase reporter assay we show that the minor allele of rs570754792, a variant located in the extended TATA box motif, displays a significant decrease in promoter activity. Likewise, through western blot studies of extracellular and intracellular sclerostin, we observe a reduced expression of the p.Val10Ile mutant protein. Finally, using a circular chromosome conformation capture assay (4C-seq) in 3 bone cell types (MSC, hFOB, Saos-2), we have detected physical interactions between the SOST proximal promoter and the ECR5 enhancer, several additional enhancers located between EVT4 and MEOX1 and a distant region containing exon 18 of DHX8. In conclusion, SOST presents functional regulatory and missense variants that affect its expression and displays physical contacts with far reaching genomic sequences, which may play a role in its regulation within bone cells.

Epstein-Barr Virus Episome Physically Interacts with Active Regions of the Host Genome in Lymphoblastoid Cells.

The Epstein-Barr virus (EBV) episome is known to interact with the three-dimensional structure of the human genome in infected cells. However, the exact locations of these interactions and their potential functional consequences remain unclear. Recently, high-resolution chromatin conformation capture (Hi-C) assays in lymphoblastoid cells have become available, enabling us to precisely map the contacts between the EBV episome(s) and the human host genome. Using available Hi-C data at a 10-kb resolution, we have identified 15,000 reproducible contacts between EBV episome(s) and the human genome. These contacts are highly enriched in chromatin regions denoted by typical or super enhancers and active markers, including histone H3K27ac and H3K4me1. Additionally, these contacts are highly enriched at loci bound by host transcription factors that regulate B cell growth (e.g., IKZF1 and RUNX3), factors that enhance cell proliferation (e.g., HDGF), or factors that promote viral replication (e.g., NBS1 and NFIC). EBV contacts show nearly 2-fold enrichment in host regions bound by EBV nuclear antigen 2 (EBNA2) and EBNA3 transcription factors. Circular chromosome conformation capture followed by sequencing (4C-seq) using the EBV origin of plasmid replication (oriP) as a "bait" in lymphoblastoid cells further confirmed contacts with active chromatin regions. Collectively, our analysis supports interactions between EBV episome(s) and active regions of the human genome in lymphoblastoid cells.IMPORTANCE EBV is associated with ∼200,000 cancers each year. In vitro, EBV can transform primary human B lymphocytes into immortalized cell lines. EBV-encoded proteins, along with noncoding RNAs and microRNAs, hijack cellular proteins and pathways to control cell growth. EBV nuclear proteins usurp normal transcriptional programs to activate the expression of key oncogenes, including MYC, to provide a proliferation signal. EBV nuclear antigens also repress CDKN2A to suppress senescence. EBV membrane protein activates NF-κB to provide survival signals. EBV genomes are maintained by EBNA1, which tethers EBV episomes to the host chromosomes during mitosis. However, little is known about where EBV episomes are located in interphase cells. In interphase cells, EBV promoters drive the expression of latency genes, while oriP functions as an enhancer for these promoters. In this study, integrative analyses of published lymphoblastoid cell line (LCL) Hi-C data and our 4C-seq experiments position EBV episomes to host genomes with active epigenetic marks. These contact points were significantly enriched for super enhancers. The close proximity of EBV episomes and the super enhancers that are enriched for transcription cofactors or mediators in lymphoblasts may benefit EBV gene expression, suggesting a novel mechanism of transcriptional activation.

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.

Insulin promoter in human pancreatic ß cells contacts diabetes susceptibility loci and regulates genes affecting insulin metabolism.

Both type 1 and type 2 diabetes involve a complex interplay between genetic, epigenetic, and environmental factors. Our laboratory has been interested in the physical interactions, in nuclei of human pancreatic ß cells, between the insulin (INS) gene and other genes that are involved in insulin metabolism. We have identified, using Circularized Chromosome Conformation Capture (4C), many physical contacts in a human pancreatic ß cell line between the INS promoter on chromosome 11 and sites on most other chromosomes. Many of these contacts are associated with type 1 or type 2 diabetes susceptibility loci. To determine whether physical contact is correlated with an ability of the INS locus to affect expression of these genes, we knock down INS expression by targeting the promoter; 259 genes are either up or down-regulated. Of these, 46 make physical contact with INS We analyze a subset of the contacted genes and show that all are associated with acetylation of histone H3 lysine 27, a marker of actively expressed genes. To demonstrate the usefulness of this approach in revealing regulatory pathways, we identify from among the contacted sites the previously uncharacterized gene SSTR5-AS1 and show that it plays an important role in controlling the effect of somatostatin-28 on insulin secretion. These results are consistent with models in which clustering of genes supports transcriptional activity. This may be a particularly important mechanism in pancreatic ß cells and in other cells where a small subset of genes is expressed at high levels.

Regulatory landscape fusion in rhabdomyosarcoma through interactions between the PAX3 promoter and FOXO1 regulatory elements.

BACKGROUND: The organisation of vertebrate genomes into topologically associating domains (TADs) is believed to facilitate the regulation of the genes located within them. A remaining question is whether TAD organisation is achieved through the interactions of the regulatory elements within them or if these interactions are favoured by the pre-existence of TADs. If the latter is true, the fusion of two independent TADs should result in the rewiring of the transcriptional landscape and the generation of ectopic contacts. RESULTS: We show that interactions within the PAX3 and FOXO1 domains are restricted to their respective TADs in normal conditions, while in a patient-derived alveolar rhabdomyosarcoma cell line, harbouring the diagnostic t(2;13)(q35;q14) translocation that brings together the PAX3 and FOXO1 genes, the PAX3 promoter interacts ectopically with FOXO1 sequences. Using a combination of 4C-seq datasets, we have modelled the three-dimensional organisation of the fused landscape in alveolar rhabdomyosarcoma. CONCLUSIONS: The chromosomal translocation that leads to alveolar rhabdomyosarcoma development generates a novel TAD that is likely to favour ectopic PAX3:FOXO1 oncogene activation in non-PAX3 territories. Rhabdomyosarcomas may therefore arise from cells which do not normally express PAX3. The borders of this novel TAD correspond to the original 5'- and 3'- borders of the PAX3 and FOXO1 TADs, respectively, suggesting that TAD organisation precedes the formation of regulatory long-range interactions. Our results demonstrate that, upon translocation, novel regulatory landscapes are formed allowing new intra-TAD interactions between the original loci involved.

Unbiased Interrogation of 3D Genome Topology Using Chromosome Conformation Capture Coupled to High-Throughput Sequencing (4C-Seq).

The development and widespread implementation of chromosome conformation capture (3C) technology has allowed unprecedented new insight into how chromosomes are folded in three-dimensional (3D) space. 3C and its derivatives have contributed tremendously to the now widely accepted view that genome topology plays an important role in many major cellular processes, at a chromosome-wide scale, but certainly also at the level of individual genetic loci. A particularly popular application of 3C technology is to study transcriptional regulation, allowing researchers to draw maps of gene regulatory connections beyond the linear genome through addition of the third dimension. In this chapter, we provide a highly detailed protocol describing 3C coupled to high-throughput sequencing (referred to as 3C-Seq or more commonly 4C-Seq), allowing the unbiased interrogation of genome-wide chromatin interactions with specific genomic regions of interest. Interactions between spatially clustered DNA fragments are revealed by crosslinking the cells with formaldehyde, digesting the genome with a restriction endonuclease and performing a proximity ligation step to link interacting genomic fragments. Next, interactions with a selected DNA fragment are extracted from the 3C library through a second round of digestion and ligation followed by an inverse PCR. The generated products are immediately compatible with high-throughput sequencing, and amplicons from different PCR reactions can easily be multiplexed to dramatically increase throughput. Finally, we provide suggestions for data analysis and visualization.

Determination of High-Resolution 3D Chromatin Organization Using Circular Chromosome Conformation Capture (4C-seq).

3D chromatin organization is essential for many aspects of transcriptional regulation. Circular Chromosome Conformation Capture followed by Illumina sequencing (4C-seq) is among the most powerful techniques to determine 3D chromatin organization. 4C-seq, like other modifications of the original 3C technique, uses the principle of "proximity ligation" to identify and quantify ten thousands of genomic interactions at a kilobase scale in a single experiment for predefined loci in the genome.In this chapter we focus on the experimental steps in the 4C-seq protocol, providing detailed descriptions on the preparation of cells, the construction of the circularized 3C library and the generation of the Illumina high throughput sequencing library. This protocol is particularly suited for the use of mammalian tissue samples, but can be used with minimal changes on circulating cells and cell lines from other sources as well. In the final section of this chapter, we provide a brief overview of data analysis approaches, accompanied by links to publicly available analysis tools.

Assay for transposase-accessible chromatin and circularized chromosome conformation capture, two methods to explore the regulatory landscapes of genes in zebrafish.

Accurate transcriptional control of genes is fundamental for the correct functioning of organs and developmental processes. This control depends on the interplay between the promoter of genes and other noncoding sequences, whose interaction is mediated by 3D chromatin arrangements. Thus, the detailed description of transcriptional regulatory landscapes is essential to understand the mechanisms of transcriptional regulation. However, to achieve that, two important challenges have to be faced: (1) the identification of the noncoding sequences that contribute to gene transcription and (2) the association of these sequences to the respective genes they control. In this chapter, we describe two protocols that allow overcoming these important challenges: the assay for transposase-accessible chromatin using sequencing (ATAC-seq) and circularized chromosome conformation capture (4C-seq). ATAC-seq is a very efficient technique that, using a very low number of cells as starting material, allows the identification of active chromatin regions genome wide, whereas 4C-seq detects the subset of sequences that interact specifically with the promoter of a given gene. When combined, both techniques provide a comprehensive snapshot of the regulatory landscapes of developmental genes. The protocols we present here have been optimized for teleost fish samples, zebrafish and medaka, allowing the in-depth study of transcriptional regulation in these two emerging animal models. Given the amenability and easy genetic manipulation of these two experimental systems, we anticipate that they will be important in revealing general principles of the vertebrate regulatory genome.