fourSig: a method for determining chromosomal interactions in 4C-Seq data

Gene-to-gene coordinated regulation of transcription and alternative splicing by 3D chromatin remodeling upon NF-κB activation

The NF-κB protein p65/RelA plays a pivotal role in coordinating gene expression in response to diverse stimuli, including viral infections. At the chromatin level, p65/RelA regulates gene transcription and alternative splicing through promoter enrichment and genomic exon occupancy, respectively. The intricate ways in which p65/RelA simultaneously governs these functions across various genes remain to be fully elucidated. In this study, we employed the HTLV-1 Tax oncoprotein, a potent activator of NF-κB, to investigate its influence on the three-dimensional organization of the genome, a key factor in gene regulation. We discovered that Tax restructures the 3D genomic landscape, bringing together genes based on their regulation and splicing patterns. Notably, we found that the Tax-induced gene-gene contact between the two master genes NFKBIA and RELA is associated with their respective changes in gene expression and alternative splicing. Through dCas9-mediated approaches, we demonstrated that NFKBIA-RELA interaction is required for alternative splicing regulation and is caused by an intragenic enrichment of p65/RelA on RELA. Our findings shed light on new regulatory mechanisms upon HTLV-1 Tax and underscore the integral role of p65/RelA in coordinated regulation of NF-κB-responsive genes at both transcriptional and splicing levels in the context of the 3D genome.

Tissue-specific Grb10/Ddc insulator drives allelic architecture for cardiac development

Allele-specific expression of imprinted gene clusters is governed by gametic DNA methylation at master regulators called imprinting control regions (ICRs). Non-gametic or secondary differentially methylated regions (DMRs) at promoters and exonic regions reinforce monoallelic expression but do not control an entire cluster. Here, we unveil an unconventional secondary DMR that is indispensable for tissue-specific imprinting of two previously unlinked genes, Grb10 and Ddc. Using polymorphic mice, we mapped an intronic secondary DMR at Grb10 with paternal-specific CTCF binding (CBR2.3) that forms contacts with Ddc. Deletion of paternal CBR2.3 removed a critical insulator, resulting in substantial shifting of chromatin looping and ectopic enhancer-promoter contacts. Destabilized gene architecture precipitated abnormal Grb10-Ddc expression with developmental consequences in the heart and muscle. Thus, we redefine the Grb10-Ddc imprinting domain by uncovering an unconventional intronic secondary DMR that functions as an insulator to instruct the tissue-specific, monoallelic expression of multiple genes-a feature previously ICR exclusive.

seRNA PAM controls skeletal muscle satellite cell proliferation and aging through trans regulation of Timp2 expression synergistically with Ddx5

Muscle satellite cells (SCs) are responsible for muscle homeostasis and regeneration and lncRNAs play important roles in regulating SC activities. Here, in this study, we identify PAM (Pax7 Associated Muscle lncRNA) that is induced in activated/proliferating SCs upon injury to promote SC proliferation as myoblast cells. PAM is generated from a myoblast-specific super-enhancer (SE); as a seRNA it binds with a number of target genomic loci predominantly in trans. Further studies demonstrate that it interacts with Ddx5 to tether PAM SE to its inter-chromosomal targets Timp2 and Vim to activate the gene expression. Lastly, we show that PAM expression is increased in aging SCs, which leads to enhanced inter-chromosomal interaction and target genes upregulation. Altogether, our findings identify PAM as a previously unknown lncRNA that regulates both SC proliferation and aging through its trans gene regulatory activity.

Allele-specific aberration of imprinted domain chromosome architecture associates with large offspring syndrome

Large offspring syndrome (LOS) and Beckwith-Wiedemann syndrome are similar epigenetic congenital overgrowth conditions in ruminants and humans, respectively. We have reported global loss-of-imprinting, methylome epimutations, and gene misregulation in LOS. However, less than 4% of gene misregulation can be explained with short range (<20kb) alterations in DNA methylation. Therefore, we hypothesized that methylome epimutations in LOS affect chromosome architecture which results in misregulation of genes located at distances >20kb in cis and in trans (other chromosomes). Our analyses focused on two imprinted domains that frequently reveal misregulation in these syndromes, namely KvDMR1 and IGF2R. Using bovine fetal fibroblasts, we identified CTCF binding at IGF2R imprinting control region but not KvDMR1, and allele-specific chromosome architecture of these domains in controls. In LOS, analyses identified erroneous long-range contacts and clustering tendency in the direction of expression of misregulated genes. In conclusion, altered chromosome architecture is associated with LOS.

•

IGF2R imprinted domain has allele-specific chromosome architecture in bovines

•

In bovines, CTCF binds at IGF2R imprinting control region but not at KvDMR1

•

Bovine large offspring syndrome (LOS) shows altered chromosome architecture at IGF2R

•

Misregulated genes in LOS exhibit genomic location-based clustering tendency

Detection of Allele-Specific 3D Chromatin Interactions Using High-Resolution In-Nucleus 4C-seq

Chromosome conformation capture techniques are a set of methods used to determine 3D genome organization through the capture and identification of physical contacts between pairs of genomic loci. Among them, 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing) allows for the identification and quantification of the sequences interacting with a preselected locus of interest. 4C-seq has been widely used in the literature, mainly to study chromatin loops between enhancers and promoters or between CTCF binding sites and to identify chromatin domain boundaries. As 3D-contacts may be established in an allele-specific manner, we describe an up-to-date allele-specific 4C-seq protocol, starting from the selection of allele-specific viewpoints to Illumina sequencing. This protocol has mainly been optimized for cultured mammalian cells, but can be adapted for other cell types with relatively minor changes in initial steps.

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.

The role of CTCF in the organization of the centromeric 11p15 imprinted domain interactome

DNA methylation, chromatin-binding proteins, and DNA looping are common components regulating genomic imprinting which leads to parent-specific monoallelic gene expression. Loss of methylation (LOM) at the human imprinting center 2 (IC2) on chromosome 11p15 is the most common cause of the imprinting overgrowth disorder Beckwith-Wiedemann Syndrome (BWS). Here, we report a familial transmission of a 7.6 kB deletion that ablates the core promoter of KCNQ1. This structural alteration leads to IC2 LOM and causes recurrent BWS. We find that occupancy of the chromatin organizer CTCF is disrupted proximal to the deletion, which causes chromatin architecture changes both in cis and in trans. We also profile the chromatin architecture of IC2 in patients with sporadic BWS caused by isolated LOM to identify conserved features of IC2 regulatory disruption. A strong interaction between CTCF sites around KCNQ1 and CDKN1C likely drive their expression on the maternal allele, while a weaker interaction involving the imprinting control region element may impede this connection and mediate gene silencing on the paternal allele. We present an imprinting model in which KCNQ1 transcription is necessary for appropriate CTCF binding and a novel chromatin conformation to drive allele-specific gene expression.

A non-coding genetic variant associated with abdominal aortic aneurysm alters ERG gene regulation

Abdominal aortic aneurysm (AAA) is a major cause of sudden death in the elderly. While AAA has some overlapping genetic and environmental risk factors with atherosclerosis, there are substantial differences, and AAA-specific medication is lacking. A recent meta-analysis of genome-wide association studies has identified four novel single-nucleotide polymorphisms (SNPs) specifically associated with AAA. Here, we investigated the gene regulatory function for one of four non-coding SNPs associated with AAA, rs2836411, which is located in an intron of the ERG gene. Rs2836411 resides within a >70 kb super-enhancer that has high levels of H3K27ac and H3K4me1 in vascular endothelial and haematopoietic cell types. Enhancer luciferase assays in cell lines showed that the risk allele significantly alters enhancer activity. The risk allele also correlates with reduced ERG expression in aortic and other vascular tissues. To identify whether rs2836411 directly contacts the promoters of ERG and/or of genes further away, we performed allele-specific circular chromosome conformation capture sequencing. In vascular endothelial cells, which express ERG, the SNP region interacts highly within the super-enhancer, while in vascular smooth muscle cells, which do not express ERG, the interactions are distributed across a wider region that includes neighbouring genes. Furthermore, the risk allele has fewer interactions within the super-enhancer compared to the protective allele. In conclusion, our results indicate that rs2836411 likely affects ERG expression by altering enhancer activity and changing local chromatin interactions. ERG is involved in vascular development, angiogenesis, and inflammation in atherosclerosis; therefore mechanistically, rs2836411 could contribute to AAA by modulating ERG levels.

Molecular and computational approaches to map regulatory elements in 3D chromatin structure

Epigenetic marks do not change the sequence of DNA but affect gene expression in a cell-type specific manner by altering the activities of regulatory elements. Development of new molecular biology assays, sequencing technologies, and computational approaches enables us to profile the human epigenome in three-dimensional structure genome-wide. Here we describe various molecular biology techniques and bioinformatic tools that have been developed to measure the activities of regulatory elements and their chromatin interactions. Moreover, we list currently available three-dimensional epigenomic data sets that are generated in various human cell types and tissues to assist in the design and analysis of research projects.

Benchmarking of 4C-seq pipelines based on real and simulated data

Motivation

With its capacity for high-resolution data output in one region of interest, chromosome conformation capture combined with high-throughput sequencing (4C-seq) is a state-of-the-art next-generation sequencing technique that provides epigenetic insights, and regularly advances current medical research. However, 4C-seq data are complex and prone to biases, and while specialized programs exist, an unbiased, extensive benchmarking is still lacking. Furthermore, neither substantial datasets with fully characterized ground truth, nor simulation programs for realistic 4C-seq data have been published.

Results

We conducted a benchmarking study on 66 4C-seq samples from 20 datasets, and developed a novel 4C-seq simulation software, Basic4CSim, to allow for detailed comparisons of 4C-seq algorithms on 50 simulated datasets with 10–120 samples each. Simulations and benchmarking were adapted to address different characteristics of 4C-seq data. Simulated data were compared with published samples to validate simulation settings. We identified differences between 4C-seq algorithms in terms of precision, recall, interaction structure, and run time, and observed general trends. Novel differential pipeline versions of single-sample based 4C-seq algorithms were included in the benchmarking. While no single tool was optimally suited for both near-cis and far-cis, and both single-sample and differential analyses, choosing a high-performing algorithm variant did improve results considerably. For near-cis scenarios, r3Cseq, peakC and FourCSeq offered high precision, while fourSig demonstrated high overall F₁ scores in far-cis analyses. Finally, 4C-seq simulations may aid in the development of improved analysis algorithms.

Availability and implementation

Basic4CSim is available at https://github.com/walter–ca/Basic4CSim.

Supplementary information

Supplementary data are available at Bioinformatics online.

Application of Chromosome Conformation Capture Method for Detection MYC/TRD Chromosomal Translocation in Leukemia Cell Line

Background: Chromosomal breakpoints are the most common cause of hereditary diseases and cancers. Today, many standard clinical methods such as cytogenetic and PCR based techniques are used which have limitation regarding detection resolution. Chromosome conformation capture is a method for detecting gene proximity and chromosomal rearrangements. Materials and Methods: In this study, SKW3 cell line was used for detecting t(8;14)(q24;q11) using a 3C-based technique. SKW3 cell line was used for 3C library preparation. For Inverse PCR, two regions were selected in upstream and downstream of the viewpoint locus on chromosome 8-MYC gene based on EcoRI restriction sites. The captured sequence with intra-chromosomal interaction between chr8-c-MYC and chr14-TRD was selected for the translocation PCR primer design. Results: The DNA fragment captured in 3C PCR showed a specific TRD sequence translocated downstream of the MYC gene. Translocation PCR demonstrated the existence of (8; 14) (q24; q11) MYC /TRD in both library and genomic DNA. Conclusion: This result demonstrated 3C- based method could be used as a useful low-cost easy operating technique in chromosomal rearrangements detection. In this study, the integration of whole genome library monitoring and PCR method was used as a high- through put method in chromosomal breakpoints detection.

Analysis of 4C-seq data: A comparison of methods

The circular chromosome conformation capture technique followed by sequencing (4C-seq) has been used in a number of studies to investigate chromosomal interactions between DNA fragments. Computational pipelines have been developed and published that offer various possibilities of 4C-seq data processing and statistical analysis. Here, we present an overview of four of such pipelines (fourSig, FourCSeq, 4C-ker and w4Cseq) taking into account the most important stages of computations. We provide comparisons of the methods and discuss their advantages and possible weaknesses. We illustrate the results with the use of data obtained for two different species, in a study devoted to vernalization control in Arabidopsis thaliana by the FLOWERING LOCUS C (FLC) gene and to long-range chromatin interactions in mouse embryonic stem cells.

4See: A Flexible Browser to Explore 4C Data

It is established that transcription of many metazoan genes is regulated by distal regulatory sequences beyond the promoter. Enhancers have been identified at up to megabase distances from their regulated genes, and/or proximal to or within the introns of unregulated genes. The unambiguous identification of the target genes of newly identified regulatory elements can thus be challenging. Well-studied enhancers have been found to come into direct physical proximity with regulated genes, presumably by the formation of chromatin loops. Chromosome conformation capture (3C) derivatives that assess the frequency of proximity between different genetic elements is thus a popular method for exploring gene regulation by distal regulatory elements. For studies of chromatin loops and promoter-enhancer communication, 4C (circular chromosome conformation capture) is one of the methods of choice, optimizing cost (required sequencing depth), throughput, and resolution. For ease of visual inspection of 4C data we present 4See, a versatile and user-friendly browser. 4See allows 4C profiles from the same bait to be flexibly plotted together, allowing biological replicates to either be compared, or pooled for comparisons between different cell types or experimental conditions. 4C profiles can be integrated with gene tracks, linear epigenomic profiles, and annotated regions of interest, such as called significant interactions, allowing rapid data exploration with limited computational resources or bioinformatics expertise.

Distinct IL-1α-responsive enhancers promote acute and coordinated changes in chromatin topology in a hierarchical manner

How cytokine-driven changes in chromatin topology are converted into gene regulatory circuits during inflammation still remains unclear. Here, we show that interleukin (IL)-1α induces acute and widespread changes in chromatin accessibility via the TAK1 kinase and NF-κB at regions that are highly enriched for inflammatory disease-relevant SNPs. Two enhancers in the extended chemokine locus on human chromosome 4 regulate the IL-1α-inducible IL8 and CXCL1-3 genes. Both enhancers engage in dynamic spatial interactions with gene promoters in an IL-1α/TAK1-inducible manner. Microdeletions of p65-binding sites in either of the two enhancers impair NF-κB recruitment, suppress activation and biallelic transcription of the IL8/CXCL2 genes, and reshuffle higher-order chromatin interactions as judged by i4C interactome profiles. Notably, these findings support a dominant role of the IL8 "master" enhancer in the regulation of sustained IL-1α signaling, as well as for IL-8 and IL-6 secretion. CRISPR-guided transactivation of the IL8 locus or cross-TAD regulation by TNFα-responsive enhancers in a different model locus supports the existence of complex enhancer hierarchies in response to cytokine stimulation that prime and orchestrate proinflammatory chromatin responses downstream of NF-κB.

Enhancer variants associated with Alzheimer's disease affect gene expression via chromatin looping

BACKGROUND

Genome-wide association studies (GWASs) have identified single-nucleotide polymorphisms (SNPs) that may be genetic factors underlying Alzheimer's disease (AD). However, how these AD-associated SNPs (AD SNPs) contribute to the pathogenesis of this disease is poorly understood because most of them are located in non-coding regions, such as introns and intergenic regions. Previous studies reported that some disease-associated SNPs affect regulatory elements including enhancers. We hypothesized that non-coding AD SNPs are located in enhancers and affect gene expression levels via chromatin loops.

METHODS

To characterize AD SNPs within non-coding regions, we extracted 406 AD SNPs with GWAS p-values of less than 1.00 × 10- 6 from the GWAS catalog database. Of these, we selected 392 SNPs within non-coding regions. Next, we checked whether those non-coding AD SNPs were located in enhancers that typically regulate gene expression levels using publicly available data for enhancers that were predicted in 127 human tissues or cell types. We sought expression quantitative trait locus (eQTL) genes affected by non-coding AD SNPs within enhancers because enhancers are regulatory elements that influence the gene expression levels. To elucidate how the non-coding AD SNPs within enhancers affect the gene expression levels, we identified chromatin-chromatin interactions by Hi-C experiments.

RESULTS

We report the following findings: (1) nearly 30% of non-coding AD SNPs are located in enhancers; (2) eQTL genes affected by non-coding AD SNPs within enhancers are associated with amyloid beta clearance, synaptic transmission, and immune responses; (3) 95% of the AD SNPs located in enhancers co-localize with their eQTL genes in topologically associating domains suggesting that regulation may occur through chromatin higher-order structures; (4) rs1476679 spatially contacts the promoters of eQTL genes via CTCF-CTCF interactions; (5) the effect of other AD SNPs such as rs7364180 is likely to be, at least in part, indirect through regulation of transcription factors that in turn regulate AD associated genes.

CONCLUSION

Our results suggest that non-coding AD SNPs may affect the function of enhancers thereby influencing the expression levels of surrounding or distant genes via chromatin loops. This result may explain how some non-coding AD SNPs contribute to AD pathogenesis.

Male-Specific Long Noncoding RNA TTTY15 Inhibits Non-Small Cell Lung Cancer Proliferation and Metastasis via TBX4

Gender affects cancer susceptibility. Currently, there are only a few studies on Y chromosome-linked long noncoding RNAs (lncRNAs), and the potential association between lncRNAs and cancers in males has not been fully elucidated. Here, we examined the expression of testis-specific transcript Y-linked 15 (TTTY15) in 37 males with non-small cell lung cancer (NSCLC), and performed circular chromosome conformation capture with next-generation sequencing to determine the genomic interaction regions of the TTTY15 gene. Our results showed that the expression levels of TTTY15 were lower in NSCLC tissues. Lower TTTY15 expression levels were associated with Tumor-Node-Metastasis (TNM) stage. A TTTY15 knockdown promoted malignant transformation of NSCLC cells. Based on the bioinformatics analysis of circular chromosome conformation capture data, we found that T-box transcription factor 4 (TBX4) may be a potential target gene of TTTY15. The RNA immunoprecipitation and chromatin immunoprecipitation results showed that TTTY15 may interact with DNA (cytosine-5)-methyltransferase 3A (DNMT3A), and the TTTY15 knockdown increased the binding of DNMT3A to the TBX4 promoter. We concluded that low TTTY15 expression correlates with worse prognosis among patients with NSCLC. TTTY15 promotes TBX4 expression via DNMT3A-mediated regulation. The identification of lncRNAs encoded by male-specific genes may help to identify potential targets for NSCLC therapy.

GemC1 governs multiciliogenesis through direct interaction with and transcriptional regulation of p73

A distinct combination of transcription factors elicits the acquisition of a specific fate and the initiation of a differentiation program. Multiciliated cells (MCCs) are a specialized type of epithelial cells that possess dozens of motile cilia on their apical surface. Defects in cilia function have been associated with ciliopathies that affect many organs, including brain and airway epithelium. Here we show that the geminin coiled-coil domain-containing protein 1 GemC1 (also known as Lynkeas) regulates the transcriptional activation of p73, a transcription factor central to multiciliogenesis. Moreover, we show that GemC1 acts in a trimeric complex with transcription factor E2F5 and tumor protein p73 (officially known as TP73), and that this complex is important for the activation of the p73 promoter. We also provide in vivo evidence that GemC1 is necessary for p73 expression in different multiciliated epithelia. We further show that GemC1 regulates multiciliogenesis through the control of chromatin organization, and the epigenetic marks/tags of p73 and Foxj 1. Our results highlight novel signaling cues involved in the commitment program of MCCs across species and tissues.This article has an associated First Person interview with the first author of the paper.

High-throughput Identification of Gene Regulatory Sequences Using Next-generation Sequencing of Circular Chromosome Conformation Capture (4C-seq)

The identification of regulatory elements for a given target gene poses a significant technical challenge owing to the variability in the positioning and effect sizes of regulatory elements to a target gene. Some progress has been made with the bioinformatic prediction of the existence and function of proximal epigenetic modifications associated with activated gene expression using conserved transcription factor binding sites. Chromatin conformation capture studies have revolutionized our ability to discover physical chromatin contacts between sequences and even within an entire genome. Circular chromatin conformation capture coupled with next-generation sequencing (4C-seq), in particular, is designed to discover all possible physical chromatin interactions for a given sequence of interest (viewpoint), such as a target gene or a regulatory enhancer. Current 4C-seq strategies directly sequence from within the viewpoint but require numerous and diverse viewpoints to be simultaneously sequenced to avoid the technical challenges of uniform base calling (imaging) with next generation sequencing platforms. This volume of experiments may not be practical for many laboratories. Here, we report a modified approach to the 4C-seq protocol that incorporates both an additional restriction enzyme digest and qPCR-based amplification steps that are designed to facilitate a greater capture of diverse sequence reads and mitigate the potential for PCR bias, respectively. Our modified 4C method is amenable to the standard molecular biology lab for assessing chromatin architecture.

4C-seq characterization of Drosophila BEAF binding regions provides evidence for highly variable long-distance interactions between active chromatin

Chromatin organization is crucial for nuclear functions such as gene regulation, DNA replication and DNA repair. Insulator binding proteins, such as the Drosophila Boundary Element-Associated Factor (BEAF), are involved in chromatin organization. To further understand the role of BEAF, we detected cis- and trans-interaction partners of four BEAF binding regions (viewpoints) using 4C (circular chromosome conformation capture) and analyzed their association with different genomic features. Previous genome-wide mapping found that BEAF usually binds near transcription start sites, often of housekeeping genes, so our viewpoints were selected to reflect this. Our 4C data show the interaction partners of our viewpoints are highly variable and generally enriched for active chromatin marks. The most consistent association was with housekeeping genes, a feature in common with our viewpoints. Fluorescence in situ hybridization indicated that the long-distance interactions occur even in the absence of BEAF. These data are most consistent with a model in which BEAF is redundant with other factors found at active promoters. Our results point to principles of long-distance interactions made by active chromatin, supporting a previously proposed model in which condensed chromatin is sticky and associates into topologically associating domains (TADs) separated by active chromatin. We propose that the highly variable long-distance interactions we detect are driven by redundant factors that open chromatin to promote transcription, combined with active chromatin filling spaces between TADs while packing of TADs relative to each other varies from cell to cell.

Genome-wide analyses reveal a role of Polycomb in promoting hypomethylation of DNA methylation valleys

BACKGROUND

Previous studies showed that the majority of developmental genes are devoid of DNA methylation at promoters even when they are repressed. Such hypomethylated regions at developmental genes are unusually large and extend well beyond proximal promoters, forming DNA methylation valleys (DMVs) or DNA methylation canyons. However, it remains elusive how most developmental genes can evade DNA methylation regardless of their transcriptional states.

RESULTS

We show that DMVs are hypomethylated in development and are highly conserved across vertebrates. Importantly, DMVs are hotspots of regulatory regions for key developmental genes and show low levels of deamination mutation rates. By analyzing a panel of DNA methylomes from mouse tissues, we identify a subset of DMVs that are dynamically methylated. These DMVs are strongly enriched for Polycomb-deposited H3K27me3 when the associated genes are silenced, and surprisingly show elevated DNA methylation upon gene activation. 4C-seq analyses indicates that Polycomb-bound DMVs form insulated and self-interacting chromatin domains. Further investigations show that DNA hypomethylation is better correlated with the binding of Polycomb than with H3K27me3. In support of a role of Polycomb in DMV hypomethylation, we observe aberrant methylation in DMVs in mouse embryonic stem cells deficient in the EED protein. Finally, we show that Polycomb regulates hypomethylation of DMVs likely through ten-eleven translocation (TET) proteins.

CONCLUSIONS

We show that Polycomb promotes the hypomethylation of DMVs near key developmental genes. These data reveal a delicate interplay between histone modifiers and DNA methylation, which contributes to their division at distinct gene targets, allowing lineage-specifying genes to largely maintain DNA methylation-free at regulatory elements.

A Dementia-Associated Risk Variant near TMEM106B Alters Chromatin Architecture and Gene Expression

Neurodegenerative diseases pose an extraordinary threat to the world's aging population, yet no disease-modifying therapies are available. Although genome-wide association studies (GWASs) have identified hundreds of risk loci for neurodegeneration, the mechanisms by which these loci influence disease risk are largely unknown. Here, we investigated the association between common genetic variants at the 7p21 locus and risk of the neurodegenerative disease frontotemporal lobar degeneration. We showed that variants associated with disease risk correlate with increased expression of the 7p21 gene TMEM106B and no other genes; co-localization analyses implicated a common causal variant underlying both association with disease and association with TMEM106B expression in lymphoblastoid cell lines and human brain. Furthermore, increases in the amount of TMEM106B resulted in increases in abnormal lysosomal phenotypes and cell toxicity in both immortalized cell lines and neurons. We then combined fine-mapping, bioinformatics, and bench-based approaches to functionally characterize all candidate causal variants at this locus. This approach identified a noncoding variant, rs1990620, that differentially recruits CTCF in lymphoblastoid cell lines and human brain to influence CTCF-mediated long-range chromatin-looping interactions between multiple cis-regulatory elements, including the TMEM106B promoter. Our findings thus provide an in-depth analysis of the 7p21 locus linked by GWASs to frontotemporal lobar degeneration, nominating a causal variant and causal mechanism for allele-specific expression and disease association at this locus. Finally, we show that genetic variants associated with risk of neurodegenerative diseases beyond frontotemporal lobar degeneration are enriched in CTCF-binding sites found in brain-relevant tissues, implicating CTCF-mediated gene regulation in risk of neurodegeneration more generally.

A DNA Contact Map for the Mouse Runx1 Gene Identifies Novel Haematopoietic Enhancers

The transcription factor Runx1 is essential for definitive haematopoiesis, and the RUNX1 gene is frequently translocated or mutated in leukaemia. Runx1 is transcribed from two promoters, P1 and P2, to give rise to different protein isoforms. Although the expression of Runx1 must be tightly regulated for normal blood development, the mechanisms that regulate Runx1 isoform expression during haematopoiesis remain poorly understood. Gene regulatory elements located in non-coding DNA are likely to be important for Runx1 transcription. Here we use circular chromosome conformation capture sequencing to identify DNA interactions with the P1 and P2 promoters of Runx1, and the previously identified +24 enhancer, in the mouse multipotent haematopoietic progenitor cell line HPC-7. The active promoter, P1, interacts with nine non-coding regions that are occupied by transcription factors within a 1 Mb topologically associated domain. Eight of nine regions function as blood-specific enhancers in zebrafish, of which two were previously shown to harbour blood-specific enhancer activity in mice. Interestingly, the +24 enhancer interacted with multiple distant regions on chromosome 16, suggesting it may regulate the expression of additional genes. The Runx1 DNA contact map identifies connections with multiple novel and known haematopoietic enhancers that are likely to be involved in regulating Runx1 expression in haematopoietic progenitor cells.

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.

Exploiting native forces to capture chromosome conformation in mammalian cell nuclei

Mammalian interphase chromosomes fold into a multitude of loops to fit the confines of cell nuclei, and looping is tightly linked to regulated function. Chromosome conformation capture (3C) technology has significantly advanced our understanding of this structure‐to‐function relationship. However, all 3C‐based methods rely on chemical cross‐linking to stabilize spatial interactions. This step remains a “black box” as regards the biases it may introduce, and some discrepancies between microscopy and 3C studies have now been reported. To address these concerns, we developed “i3C”, a novel approach for capturing spatial interactions without a need for cross‐linking. We apply i3C to intact nuclei of living cells and exploit native forces that stabilize chromatin folding. Using different cell types and loci, computational modeling, and a methylation‐based orthogonal validation method, “TALE‐iD”, we show that native interactions resemble cross‐linked ones, but display improved signal‐to‐noise ratios and are more focal on regulatory elements and CTCF sites, while strictly abiding to topologically associating domain restrictions.

Long noncoding RNA TUG1 is downregulated in non-small cell lung cancer and can regulate CELF1 on binding to PRC2

Background

Long noncoding RNAs (lncRNAs) play crucial roles in tumorigenesis, and lncRNA taurine-upregulated gene 1 (TUG1) has been proven to be associated with several human cancers. However, the mechanisms of TUG1-involved regulation remain largely unknown.

Methods

We examined the expressions of TUG1 in a cohort of 89 patients with non-small cell lung cancer (NSCLC) to determine the association between TUG1 expression and clinical parameters. We used circular chromosome conformation capture (4C) coupled with next-generation sequencing to explore the genome regions that interact with TUG1 and the TUG1-mediated regulation.

Results

TUG1 was significantly downregulated, and the TUG1 downregulation correlated with sex (p = 0.006), smoking status (p = 0.016), and tumor differentiation grade (p = 0.001). Knockdown of TUG1 significantly promoted the proliferation of NSCLC cells. According to the bioinformatic analysis result of TUG1 4C sequencing data, 83 candidate genes and their interaction regions were identified. Among these candidate genes, CUGBP and Elav-like family member 1 (CELF1) are potential targets of TUG1 in-trans regulation. To confirm the interaction between TUG1 and CELF1, relative expressions of CELF1 were examined in TUG1 knockdown H520 cells; results showed that CELF1 was significantly upregulated in TUG1 knockdown H520 cells. RNA immunoprecipitation was then performed to examine whether TUG1 RNA was bound to PRC2, a TUG1-involved regulation mechanism reported in previous studies. The results demonstrated that TUG1 RNA was bound to enhancer of zeste protein 2/embryonic ectoderm development (EZH2/EED), which is essential for PRC2. Finally, our designed ChIP assay revealed that the EZH2/EED was bound to the promotor region of CELF1 within 992 bp upstream of the transcript start site.

Conclusion

TUG1 is downregulated in NSCLC. Using TUG1 4C sequencing and bioinformatic analysis, we found CELF1 to be a potential target of TUG1 RNA in in-trans regulation. Moreover, subsequent experiments showed that TUG1 RNA could bind to PRC2 in the promotor region of CELF1 and negatively regulate CELF1 expressions in H520 cells. Our results may facilitate developing new treatment modalities targeting TUG1/PRC2/CELF1 interactions in patients with NSCLC.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2569-6) contains supplementary material, which is available to authorized users.

Long-Range Chromosome Interactions Mediated by Cohesin Shape Circadian Gene Expression

Mammalian circadian rhythm is established by the negative feedback loops consisting of a set of clock genes, which lead to the circadian expression of thousands of downstream genes in vivo. As genome-wide transcription is organized under the high-order chromosome structure, it is largely uncharted how circadian gene expression is influenced by chromosome architecture. We focus on the function of chromatin structure proteins cohesin as well as CTCF (CCCTC-binding factor) in circadian rhythm. Using circular chromosome conformation capture sequencing, we systematically examined the interacting loci of a Bmal1-bound super-enhancer upstream of a clock gene Nr1d1 in mouse liver. These interactions are largely stable in the circadian cycle and cohesin binding sites are enriched in the interactome. Global analysis showed that cohesin-CTCF co-binding sites tend to insulate the phases of circadian oscillating genes while cohesin-non-CTCF sites are associated with high circadian rhythmicity of transcription. A model integrating the effects of cohesin and CTCF markedly improved the mechanistic understanding of circadian gene expression. Further experiments in cohesin knockout cells demonstrated that cohesin is required at least in part for driving the circadian gene expression by facilitating the enhancer-promoter looping. This study provided a novel insight into the relationship between circadian transcriptome and the high-order chromosome structure.

Chromatin interactions and candidate genes at ten prostate cancer risk loci

Genome-wide association studies have identified more than 100 common single nucleotide polymorphisms (SNPs) that are associated with prostate cancer risk. However, the vast majority of these SNPs lie in noncoding regions of the genome. To test whether these risk SNPs regulate their target genes through long-range chromatin interactions, we applied capture-based 3C sequencing technology to investigate possible cis-interactions at ten prostate cancer risk loci in six cell lines. We identified significant physical interactions between risk regions and their potential target genes including CAPG at 2p11.2, C2orf43 at 2p24.1, RFX6 at 6q22.1, NFASC at 1q32.1, MYC at 8q24.1 and AGAP7P at 10q11.23. Most of the interaction peaks were co-localized to regions of active histone modification and transcription factor binding sites. Expression quantitative trait locus (eQTL) analysis showed suggestive eQTL signals at rs1446669, rs699664 and rs1078004 for CAPG (p < 0.004), rs13394027 for C2orf43 (p = 2.25E-27), rs10993994 and rs4631830 for AGAP7P (p < 8.02E-5). Further analysis revealed an enhancer activity at genomic region surrounding rs4631830 which was expected to disrupt HOXB-like DNA binding affinity. This study identifies a set of candidate genes and their potential regulatory variants, and provides additional evidence showing the role of long-range chromatin interactions in prostate cancer etiology.

4C-ker: A Method to Reproducibly Identify Genome-Wide Interactions Captured by 4C-Seq Experiments

4C-Seq has proven to be a powerful technique to identify genome-wide interactions with a single locus of interest (or "bait") that can be important for gene regulation. However, analysis of 4C-Seq data is complicated by the many biases inherent to the technique. An important consideration when dealing with 4C-Seq data is the differences in resolution of signal across the genome that result from differences in 3D distance separation from the bait. This leads to the highest signal in the region immediately surrounding the bait and increasingly lower signals in far-cis and trans. Another important aspect of 4C-Seq experiments is the resolution, which is greatly influenced by the choice of restriction enzyme and the frequency at which it can cut the genome. Thus, it is important that a 4C-Seq analysis method is flexible enough to analyze data generated using different enzymes and to identify interactions across the entire genome. Current methods for 4C-Seq analysis only identify interactions in regions near the bait or in regions located in far-cis and trans, but no method comprehensively analyzes 4C signals of different length scales. In addition, some methods also fail in experiments where chromatin fragments are generated using frequent cutter restriction enzymes. Here, we describe 4C-ker, a Hidden-Markov Model based pipeline that identifies regions throughout the genome that interact with the 4C bait locus. In addition, we incorporate methods for the identification of differential interactions in multiple 4C-seq datasets collected from different genotypes or experimental conditions. Adaptive window sizes are used to correct for differences in signal coverage in near-bait regions, far-cis and trans chromosomes. Using several datasets, we demonstrate that 4C-ker outperforms all existing 4C-Seq pipelines in its ability to reproducibly identify interaction domains at all genomic ranges with different resolution enzymes.

Chromatin topology is coupled to Polycomb group protein subnuclear organization

The genomes of metazoa are organized at multiple scales. Many proteins that regulate genome architecture, including Polycomb group (PcG) proteins, form subnuclear structures. Deciphering mechanistic links between protein organization and chromatin architecture requires precise description and mechanistic perturbations of both. Using super-resolution microscopy, here we show that PcG proteins are organized into hundreds of nanoscale protein clusters. We manipulated PcG clusters by disrupting the polymerization activity of the sterile alpha motif (SAM) of the PcG protein Polyhomeotic (Ph) or by increasing Ph levels. Ph with mutant SAM disrupts clustering of endogenous PcG complexes and chromatin interactions while elevating Ph level increases cluster number and chromatin interactions. These effects can be captured by molecular simulations based on a previously described chromatin polymer model. Both perturbations also alter gene expression. Organization of PcG proteins into small, abundant clusters on chromatin through Ph SAM polymerization activity may shape genome architecture through chromatin interactions.

In the loop: promoter-enhancer interactions and bioinformatics

Enhancer-promoter regulation is a fundamental mechanism underlying differential transcriptional regulation. Spatial chromatin organization brings remote enhancers in contact with target promoters in cis to regulate gene expression. There is considerable evidence for promoter-enhancer interactions (PEIs). In the recent years, genome-wide analyses have identified signatures and mapped novel enhancers; however, being able to precisely identify their target gene(s) requires massive biological and bioinformatics efforts. In this review, we give a short overview of the chromatin landscape and transcriptional regulation. We discuss some key concepts and problems related to chromatin interaction detection technologies, and emerging knowledge from genome-wide chromatin interaction data sets. Then, we critically review different types of bioinformatics analysis methods and tools related to representation and visualization of PEI data, raw data processing and PEI prediction. Lastly, we provide specific examples of how PEIs have been used to elucidate a functional role of non-coding single-nucleotide polymorphisms. The topic is at the forefront of epigenetic research, and by highlighting some future bioinformatics challenges in the field, this review provides a comprehensive background for future PEI studies.

Basic4Cseq: an R/Bioconductor package for analyzing 4C-seq data

SUMMARY

Basic4Cseq is an R/Bioconductor package for basic filtering, analysis and subsequent near-cis visualization of 4C-seq data. The package processes aligned 4C-seq raw data stored in binary alignment/map (BAM) format and maps the short reads to a corresponding virtual fragment library. Functions are included to create virtual fragment libraries providing chromosome position and further information on 4C-seq fragments (length and uniqueness of the fragment ends, and blindness of a fragment) for any BSGenome package. An optional filter is included for BAM files to remove invalid 4C-seq reads, and further filter functions are offered for 4C-seq fragments. Additionally, basic quality controls based on the read distribution are included. Fragment data in the vicinity of the experiment's viewpoint are visualized as coverage plot based on a running median approach and a multi-scale contact profile. Wig files or csv files of the fragment data can be exported for further analyses and visualizations of interactions with other programs.

AVAILABILITY AND IMPLEMENTATION

Basic4Cseq is implemented in R and available at http://www.bioconductor.org/. A vignette with detailed descriptions of the functions is included in the package.

CONTACT

Carolin.Walter@uni-muenster.de

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Evidence for local regulatory control of escape from imprinted X chromosome inactivation

X chromosome inactivation (XCI) is an epigenetic process that almost completely inactivates one of two X chromosomes in somatic cells of mammalian females. A few genes are known to escape XCI and the mechanism for this escape remains unclear. Here, using mouse trophoblast stem (TS) cells, we address whether particular chromosomal interactions facilitate escape from imprinted XCI. We demonstrate that promoters of genes escaping XCI do not congregate to any particular region of the genome in TS cells. Further, the escape status of a gene was uncorrelated with the types of genomic features and gene activity located in contacted regions. Our results suggest that genes escaping imprinted XCI do so by using the same regulatory sequences as their expressed alleles on the active X chromosome. We suggest a model where regulatory control of escape from imprinted XCI is mediated by genomic elements located in close linear proximity to escaping genes.