C-TALE, a new cost-effective method for targeted enrichment of Hi-C/3C-seq libraries

Exchange of subtelomeric regions between chromosomes 4q and 10q reverts the FSHD genotype and phenotype

The most common form of facioscapulohumeral dystrophy (FSHD1) is caused by a partial loss of the D4Z4 macrosatellite repeat array in the subtelomeric region of chromosome 4. Patients with FSHD1 typically carry 1 to 10 D4Z4 repeats, whereas nonaffected individuals have 11 to 150 repeats. The ~150-kilobyte subtelomeric region of the chromosome 10q exhibits a ~99% sequence identity to the 4q, including the D4Z4 array. Nevertheless, contractions of the chr10 array do not cause FSHD or any known disease, as in most people D4Z4 array on chr10 is flanked by the nonfunctional polyadenylation signal, not permitting the DUX4 expression. Here, we attempted to correct the FSHD genotype by a CRISPR-Cas9-induced exchange of the chr4 and chr10 subtelomeric regions. We demonstrated that the induced t(4;10) translocation can generate recombinant genotypes translated into improved FSHD phenotype. FSHD myoblasts with the t(4;10) exhibited reduced expression of the DUX4 targets, restored PAX7 target expression, reduced sensitivity to oxidative stress, and improved differentiation capacity.

Transcription decouples estrogen-dependent changes in enhancer-promoter contact frequencies and spatial proximity

How enhancers regulate their target genes in the context of 3D chromatin organization is extensively studied and models which do not require direct enhancer-promoter contact have recently emerged. Here, we use the activation of estrogen receptor-dependent enhancers in a breast cancer cell line to study enhancer-promoter communication at two loci. This allows high temporal resolution tracking of molecular events from hormone stimulation to efficient gene activation. We examine how both enhancer-promoter spatial proximity assayed by DNA fluorescence in situ hybridization, and contact frequencies resulting from chromatin in situ fragmentation and proximity ligation, change dynamically during enhancer-driven gene activation. These orthogonal methods produce seemingly paradoxical results: upon enhancer activation enhancer-promoter contact frequencies increase while spatial proximity decreases. We explore this apparent discrepancy using different estrogen receptor ligands and transcription inhibitors. Our data demonstrate that enhancer-promoter contact frequencies are transcription independent whereas altered enhancer-promoter proximity depends on transcription. Our results emphasize that the relationship between contact frequencies and physical distance in the nucleus, especially over short genomic distances, is not always a simple one.

The long non-coding RNA Meg3 mediates imprinted gene expression during stem cell differentiation

The imprinted Dlk1-Dio3 domain comprises the developmental genes Dlk1 and Rtl1, which are silenced on the maternal chromosome in different cell types. On this parental chromosome, the domain's imprinting control region activates a polycistron that produces the lncRNA Meg3 and many miRNAs (Mirg) and C/D-box snoRNAs (Rian). Although Meg3 lncRNA is nuclear and associates with the maternal chromosome, it is unknown whether it controls gene repression in cis. We created mouse embryonic stem cells (mESCs) that carry an ectopic poly(A) signal, reducing RNA levels along the polycistron, and generated Rian-/- mESCs as well. Upon ESC differentiation, we found that Meg3 lncRNA (but not Rian) is required for Dlk1 repression on the maternal chromosome. Biallelic Meg3 expression acquired through CRISPR-mediated demethylation of the paternal Meg3 promoter led to biallelic Dlk1 repression, and to loss of Rtl1 expression. lncRNA expression also correlated with DNA hypomethylation and CTCF binding at the 5'-side of Meg3. Using Capture Hi-C, we found that this creates a Topologically Associating Domain (TAD) organization that brings Meg3 close to Dlk1 on the maternal chromosome. The requirement of Meg3 for gene repression and TAD structure may explain how aberrant MEG3 expression at the human DLK1-DIO3 locus associates with imprinting disorders.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

3DGenBench: a web-server to benchmark computational models for 3D Genomics

Modeling 3D genome organisation has been booming in the last years thanks to the availability of experimental datasets of genomic contacts. However, the field is currently missing the standardisation of methods and metrics to compare predictions and experiments. We present 3DGenBench, a web server available at https://inc-cost.eu/benchmarking/, that allows benchmarking computational models of 3D Genomics. The benchmark is performed using a manually curated dataset of 39 capture Hi-C profiles in wild type and genome-edited mouse cells, and five genome-wide Hi-C profiles in human, mouse, and Drosophila cells. 3DGenBench performs two kinds of analysis, each supplied with a specific scoring module that compares predictions of a computational method to experimental data using several metrics. With 3DGenBench, the user obtains model performance scores, allowing an unbiased comparison with other models. 3DGenBench aims to become a reference web server to test new 3D genomics models and is conceived as an evolving platform where new types of analysis will be implemented in the future.

Towards the Genetic Architecture of Complex Gene Expression Traits: Challenges and Prospects for eQTL Mapping in Humans

The discovery of expression quantitative trait loci (eQTLs) and their target genes (eGenes) has not only compensated for the limitations of genome-wide association studies for complex phenotypes but has also provided a basis for predicting gene expression. Efforts have been made to develop analytical methods in statistical genetics, a key discipline in eQTL analysis. In particular, mixed model– and deep learning–based analytical methods have been extremely beneficial in mapping eQTLs and predicting gene expression. Nevertheless, we still face many challenges associated with eQTL discovery. Here, we discuss two key aspects of these challenges: 1, the complexity of eTraits with various factors such as polygenicity and epistasis and 2, the voluminous work required for various types of eQTL profiles. The properties and prospects of statistical methods, including the mixed model method, Bayesian inference, the deep learning method, and the integration method, are presented as future directions for eQTL discovery. This review will help expedite the design and use of efficient methods for eQTL discovery and eTrait prediction.

Assessment of Multiway Interactions with Tri-C

Tri-C is a chromosome conformation capture (3C) approach that can efficiently identify multiway chromatin interactions with viewpoints of interest. As opposed to pair-wise interactions identified in methods such as Hi-C, 4C, and Capture-C, the detection of multiway interactions allows researchers to investigate how multiple cis-regulatory elements interact together in higher-order structures in single nuclei and address questions regarding structural cooperation between these elements. Here, we describe the procedure for designing and performing a Tri-C experiment.

Single-cell Hi-C data analysis: safety in numbers

Over the past decade, genome-wide assays for chromatin interactions in single cells have enabled the study of individual nuclei at unprecedented resolution and throughput. Current chromosome conformation capture techniques survey contacts for up to tens of thousands of individual cells, improving our understanding of genome function in 3D. However, these methods recover a small fraction of all contacts in single cells, requiring specialised processing of sparse interactome data. In this review, we highlight recent advances in methods for the interpretation of single-cell genomic contacts. After discussing the strengths and limitations of these methods, we outline frontiers for future development in this rapidly moving field.

Investigation of the Basic Steps in the Chromosome Conformation Capture Procedure

A constellation of chromosome conformation capture methods (С-methods) are an important tool for biochemical analysis of the spatial interactions between DNA regions that are separated in the primary sequence. All these methods are based on the long sequence of basic steps of treating cells, nuclei, chromatin, and finally DNA, thus representing a significant technical challenge. Here, we present an in-depth study of the basic steps in the chromatin conformation capture procedure (3С), which was performed using Drosophila Schneider 2 cells as a model. We investigated the steps of cell lysis, nuclei washing, nucleoplasm extraction, chromatin treatment with SDS/Triton X-100, restriction enzyme digestion, chromatin ligation, reversion of cross-links, DNA extraction, treatment of a 3C library with RNases, and purification of the 3C library. Several options were studied, and optimal conditions were found. Our work contributes to the understanding of the 3C basic steps and provides a useful guide to the 3C procedure.

Sensitivity of cohesin-chromatin association to high-salt treatment corroborates non-topological mode of loop extrusion

Cohesin is a key organizer of chromatin folding in eukaryotic cells. The two main activities of this ring-shaped protein complex are the maintenance of sister chromatid cohesion and the establishment of long-range DNA-DNA interactions through the process of loop extrusion. Although the basic principles of both cohesion and loop extrusion have been described, we still do not understand several crucial mechanistic details. One of such unresolved issues is the question of whether a cohesin ring topologically embraces DNA string(s) during loop extrusion. Here, we show that cohesin complexes residing on CTCF-occupied genomic sites in mammalian cells do not interact with DNA topologically. We assessed the stability of cohesin-dependent loops and cohesin association with chromatin in high-ionic-strength conditions in G1-synchronized HeLa cells. We found that increased salt concentration completely displaces cohesin from those genomic regions that correspond to CTCF-defined loop anchors. Unsurprisingly, CTCF-anchored cohesin loops also dissipate in these conditions. Because topologically engaged cohesin is considered to be salt resistant, our data corroborate a non-topological model of loop extrusion. We also propose a model of cohesin activity throughout the interphase, which essentially equates the termination of non-topological loop extrusion with topological loading of cohesin. This theoretical framework enables a parsimonious explanation of various seemingly contradictory experimental findings.

Three-dimensional genome organization in the central nervous system, implications for neuropsychological disorders

Chromosomes in eukaryotic cell nuclei are highly compacted and finely organized into hierarchical three-dimensional (3D) configuration. In recent years, scientists have gained deeper understandings of 3D genome structures and revealed novel evidence linking 3D genome organization to various important cell events on the molecular level. Most importantly, alteration of 3D genome architecture has emerged as an intriguing higher order mechanism that connects disease-related genetic variants in multiple heterogenous and polygenic neuropsychological disorders, delivering novel insights into the etiology. In this review, we provide a brief overview of the hierarchical structures of 3D genome and two proposed regulatory models, loop extrusion and phase separation. We then focus on recent Hi-C data in the central nervous system and discuss 3D genome alterations during normal brain development and in mature neurons. Most importantly, we make a comprehensive review on current knowledge and discuss the role of 3D genome in multiple neuropsychological disorders, including schizophrenia, repeat expansion disorders, 22q11 deletion syndrome, and others.

Subunits of the PBAP Chromatin Remodeler Are Capable of Mediating Enhancer-Driven Transcription in Drosophila

The chromatin remodeler SWI/SNF is an important participant in gene activation, functioning predominantly by opening the chromatin structure on promoters and enhancers. Here, we describe its novel mode of action in which SWI/SNF factors mediate the targeted action of an enhancer. We studied the functions of two signature subunits of PBAP subfamily, BAP170 and SAYP, in Drosophila. These subunits were stably tethered to a transgene reporter carrying the hsp70 core promoter. The tethered subunits mediate transcription of the reporter in a pattern that is generated by enhancers close to the insertion site in multiple loci throughout the genome. Both tethered SAYP and BAP170 recruit the whole PBAP complex to the reporter promoter. However, we found that BAP170-dependent transcription is more resistant to the depletion of other PBAP subunits, suggesting that BAP170 may play a more critical role in establishing enhancer-dependent transcription.

A modified protocol of Capture-C allows affordable and flexible high-resolution promoter interactome analysis

Large-scale epigenomic projects have mapped hundreds of thousands of potential regulatory sites in the human genome, but only a small proportion of these elements are proximal to transcription start sites. It is believed that the majority of these sequences are remote promoter-activating genomic sites scattered within several hundreds of kilobases from their cognate promoters and referred to as enhancers. It is still unclear what principles, aside from relative closeness in the linear genome, determine which promoter(s) is controlled by a given enhancer; however, this understanding is of great fundamental and clinical relevance. In recent years, C-methods (chromosome conformation capture-based methods) have become a powerful tool for the identification of enhancer-promoter spatial contacts that, in most cases, reflect their functional link. Here, we describe a new hybridisation-based promoter Capture-C protocol that makes use of biotinylated dsDNA probes generated by PCR from a custom pool of long oligonucleotides. The described protocol allows high-resolution promoter interactome description, providing a flexible and cost-effective alternative to the existing promoter Capture-C modifications. Based on the obtained data, we propose several tips on probe design that could potentially improve the results of future experiments.

Studying RNA-DNA interactome by Red-C identifies noncoding RNAs associated with various chromatin types and reveals transcription dynamics

Non-coding RNAs (ncRNAs) participate in various biological processes, including regulating transcription and sustaining genome 3D organization. Here, we present a method termed Red-C that exploits proximity ligation to identify contacts with the genome for all RNA molecules present in the nucleus. Using Red-C, we uncovered the RNA-DNA interactome of human K562 cells and identified hundreds of ncRNAs enriched in active or repressed chromatin, including previously undescribed RNAs. Analysis of the RNA-DNA interactome also allowed us to trace the kinetics of messenger RNA production. Our data support the model of co-transcriptional intron splicing, but not the hypothesis of the circularization of actively transcribed genes.