Alteration of genome folding via contact domain boundary insertion

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Transcription shapes 3D chromatin organization by interacting with loop extrusion

Cohesin folds mammalian interphase chromosomes by extruding the chromatin fiber into numerous loops. "Loop extrusion" can be impeded by chromatin-bound factors, such as CTCF, which generates characteristic and functional chromatin organization patterns. It has been proposed that transcription relocalizes or interferes with cohesin and that active promoters are cohesin loading sites. However, the effects of transcription on cohesin have not been reconciled with observations of active extrusion by cohesin. To determine how transcription modulates extrusion, we studied mouse cells in which we could alter cohesin abundance, dynamics, and localization by genetic "knockouts" of the cohesin regulators CTCF and Wapl. Through Hi-C experiments, we discovered intricate, cohesin-dependent contact patterns near active genes. Chromatin organization around active genes exhibited hallmarks of interactions between transcribing RNA polymerases (RNAPs) and extruding cohesins. These observations could be reproduced by polymer simulations in which RNAPs were moving barriers to extrusion that obstructed, slowed, and pushed cohesins. The simulations predicted that preferential loading of cohesin at promoters is inconsistent with our experimental data. Additional ChIP-seq experiments showed that the putative cohesin loader Nipbl is not predominantly enriched at promoters. Therefore, we propose that cohesin is not preferentially loaded at promoters and that the barrier function of RNAP accounts for cohesin accumulation at active promoters. Altogether, we find that RNAP is an extrusion barrier that is not stationary, but rather, translocates and relocalizes cohesin. Loop extrusion and transcription might interact to dynamically generate and maintain gene interactions with regulatory elements and shape functional genomic organization.

Active enhancers strengthen insulation by RNA-mediated CTCF binding at chromatin domain boundaries

Vertebrate genomes are partitioned into chromatin domains or topologically associating domains (TADs), which are typically bound by head-to-head pairs of CTCF binding sites. Transcription at domain boundaries correlates with better insulation; however, it is not known whether the boundary transcripts themselves contribute to boundary function. Here we characterize boundary-associated RNAs genome-wide, focusing on the disease-relevant INK4a/ARF and MYC TAD. Using CTCF site deletions and boundary-associated RNA knockdowns, we observe that boundary-associated RNAs facilitate recruitment and clustering of CTCF at TAD borders. The resulting CTCF enrichment enhances TAD insulation, enhancer-promoter interactions, and TAD gene expression. Importantly, knockdown of boundary-associated RNAs results in loss of boundary insulation function. Using enhancer deletions and CRISPRi of promoters, we show that active TAD enhancers, but not promoters, induce boundary-associated RNA transcription, thus defining a novel class of regulatory enhancer RNAs.

PGG.SV: a whole-genome-sequencing-based structural variant resource and data analysis platform

Structural variations (SVs) play important roles in human evolution and diseases, but there is a lack of data resources concerning representative samples, especially for East Asians. Taking advantage of both next-generation sequencing and third-generation sequencing data at the whole-genome level, we developed the database PGG.SV to provide a practical platform for both regionally and globally representative structural variants. In its current version, PGG.SV archives 584 277 SVs obtained from whole-genome sequencing data of 6048 samples, including 1030 long-read sequencing genomes representing 177 global populations. PGG.SV provides (i) high-quality SVs with fine-scale and precise genomic locations in both GRCh37 and GRCh38, covering underrepresented SVs in existing sequencing and microarray data; (ii) hierarchical estimation of SV prevalence in geographical populations; (iii) informative annotations of SV-related genes, potential functions and clinical effects; (iv) an analysis platform to facilitate SV-based case-control association studies and (v) various visualization tools for understanding the SV structures in the human genome. Taken together, PGG.SV provides a user-friendly online interface, easy-to-use analysis tools and a detailed presentation of results. PGG.SV is freely accessible via https://www.biosino.org/pggsv.

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.

Constricted migration is associated with stable 3D genome structure differences in cancer cells

To spread from a localized tumor, metastatic cancer cells must squeeze through constrictions that cause major nuclear deformations. Since chromosome structure affects nucleus stiffness, gene regulation, and DNA repair, here, we investigate the relationship between 3D genome structure and constricted migration in cancer cells. Using melanoma (A375) cells, we identify phenotypic differences in cells that have undergone multiple rounds of constricted migration. These cells display a stably higher migration efficiency, elongated morphology, and differences in the distribution of Lamin A/C and heterochromatin. Hi-C experiments reveal differences in chromosome spatial compartmentalization specific to cells that have passed through constrictions and related alterations in expression of genes associated with migration and metastasis. Certain features of the 3D genome structure changes, such as a loss of B compartment interaction strength, are consistently observed after constricted migration in clonal populations of A375 cells and in MDA-MB-231 breast cancer cells. Our observations suggest that consistent types of chromosome structure changes are induced or selected by passage through constrictions and that these may epigenetically encode stable differences in gene expression and cellular migration phenotype.

Cohesin-mediated loop anchors confine the locations of human replication origins

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

3D or Not 3D: Shaping the Genome during Development

One of the most fundamental questions in developmental biology is how one fertilized cell can give rise to a fully mature organism and how gene regulation governs this process. Precise spatiotemporal gene expression is required for development and is believed to be achieved through a complex interplay of sequence-specific information, epigenetic modifications, trans-acting factors, and chromatin folding. Here we review the role of chromatin folding during development, the mechanisms governing 3D genome organization, and how it is established in the embryo. Furthermore, we discuss recent advances and debated questions regarding the contribution of the 3D genome to gene regulation during organogenesis. Finally, we describe the mechanisms that can reshape the 3D genome, including disease-causing structural variations and the emerging view that transposable elements contribute to chromatin organization.

Sequence-based modeling of three-dimensional genome architecture from kilobase to chromosome scale

To learn how genomic sequence influences multiscale three-dimensional (3D) genome architecture, this manuscript presents a sequence-based deep-learning approach, Orca, that predicts directly from sequence the 3D genome architecture from kilobase to whole-chromosome scale. Orca captures the sequence dependencies of structures including chromatin compartments and topologically associating domains, as well as diverse types of interactions from CTCF-mediated to enhancer-promoter interactions and Polycomb-mediated interactions with cell-type specificity. Orca enables various applications including predicting structural variant effects on multiscale genome organization and it recapitulated effects of experimentally studied variants at varying sizes (300 bp to 90 Mb). Moreover, Orca enables in silico virtual screens to probe the sequence basis of 3D genome organization at different scales. At the submegabase scale, it predicted specific transcription factor motifs underlying cell-type-specific genome interactions. At the compartment scale, virtual screens of sequence activities suggest a model for the sequence basis of chromatin compartments with a prominent role of transcription start sites.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

Genomic context sensitivity of insulator function

The specificity of interactions between genomic regulatory elements and potential target genes is influenced by the binding of insulator proteins such as CTCF, which can act as potent enhancer blockers when interposed between an enhancer and a promoter in a reporter assay. But not all CTCF sites genome-wide function as insulator elements, depending on cellular and genomic context. To dissect the influence of genomic context on enhancer blocker activity, we integrated reporter constructs with promoter-only, promoter and enhancer, and enhancer blocker configurations at hundreds of thousands of genomic sites using the Sleeping Beauty transposase. Deconvolution of reporter activity by genomic position reveals distinct expression patterns subject to genomic context, including a compartment of enhancer blocker reporter integrations with robust expression. The high density of integration sites permits quantitative delineation of characteristic genomic context sensitivity profiles and their decomposition into sensitivity to both local and distant DNase I hypersensitive sites. Furthermore, using a single-cell expression approach to test the effect of integrated reporters for differential expression of nearby endogenous genes reveals that CTCF insulator elements do not completely abrogate reporter effects on endogenous gene expression. Collectively, our results lend new insight into genomic regulatory compartmentalization and its influence on the determinants of promoter–enhancer specificity.

A chromosomal loop anchor mediates bacterial genome organization

Nucleoprotein complexes play an integral role in genome organization of both eukaryotes and prokaryotes. Apart from their role in locally structuring and compacting DNA, several complexes are known to influence global organization by mediating long-range anchored chromosomal loop formation leading to spatial segregation of large sections of DNA. Such megabase-range interactions are ubiquitous in eukaryotes, but have not been demonstrated in prokaryotes. Here, using a genome-wide sedimentation-based approach, we found that a transcription factor, Rok, forms large nucleoprotein complexes in the bacterium Bacillus subtilis. Using chromosome conformation capture and live-imaging of DNA loci, we show that these complexes robustly interact with each other over large distances. Importantly, these Rok-dependent long-range interactions lead to anchored chromosomal loop formation, thereby spatially isolating large sections of DNA, as previously observed for insulator proteins in eukaryotes.

HPV16-LINC00393 Integration Alters Local 3D Genome Architecture in Cervical Cancer Cells

High-risk human papillomavirus (hrHPV) infection and integration were considered as essential onset factors for the development of cervical cancer. However, the mechanism on how hrHPV integration influences the host genome structure remains not fully understood. In this study, we performed in situ high-throughput chromosome conformation capture (Hi-C) sequencing, chromatin immunoprecipitation and sequencing (ChIP-seq), and RNA-sequencing (RNA-seq) in two cervical cells, 1) NHEK normal human epidermal keratinocyte; and 2) HPV16-integrated SiHa tumorigenic cervical cancer cells. Our results reveal that the HPV-LINC00393 integrated chromosome 13 exhibited significant genomic variation and differential gene expression, which was verified by calibrated CTCF and H3K27ac ChIP-Seq chromatin restructuring. Importantly, HPV16 integration led to differential responses in topologically associated domain (TAD) boundaries, with a decrease in the tumor suppressor KLF12 expression downstream of LINC00393. Overall, this study provides significant insight into the understanding of HPV16 integration induced 3D structural changes and their contributions on tumorigenesis, which supplements the theory basis for the cervical carcinogenic mechanism of HPV16 integration.

Targeted Gene Delivery: Where to Land

Genome-editing technologies have the potential to correct most genetic defects involved in blood disorders. In contrast to mutation-specific editing, targeted gene insertion can correct most of the mutations affecting the same gene with a single therapeutic strategy (gene replacement) or provide novel functions to edited cells (gene addition). Targeting a selected genomic harbor can reduce insertional mutagenesis risk, while enabling the exploitation of endogenous promoters, or selected chromatin contexts, to achieve specific transgene expression levels/patterns and the modulation of disease-modifier genes. In this review, we will discuss targeted gene insertion and the advantages and limitations of different genomic harbors currently under investigation for various gene therapy applications.

CTCF and transcription influence chromatin structure re-configuration after mitosis

During mitosis, transcription is globally attenuated and chromatin architecture is dramatically reconfigured. We exploited the M- to G1-phase progression to interrogate the contributions of the architectural factor CTCF and the process of transcription to genome re-sculpting in newborn nuclei. Depletion of CTCF during the M- to G1-phase transition alters short-range compartmentalization after mitosis. Chromatin domain boundary re-formation is impaired upon CTCF loss, but a subset of boundaries, characterized by transitions in chromatin states, is established normally. Without CTCF, structural loops fail to form, leading to illegitimate contacts between cis-regulatory elements (CREs). Transient CRE contacts that are normally resolved after telophase persist deeply into G1-phase in CTCF-depleted cells. CTCF loss-associated gains in transcription are often linked to increased, normally illegitimate enhancer-promoter contacts. In contrast, at genes whose expression declines upon CTCF loss, CTCF seems to function as a conventional transcription activator, independent of its architectural role. CTCF-anchored structural loops facilitate formation of CRE loops nested within them, especially those involving weak CREs. Transcription inhibition does not significantly affect global architecture or transcription start site-associated boundaries. However, ongoing transcription contributes considerably to the formation of gene domains, regions of enriched contacts along gene bodies. Notably, gene domains emerge in ana/telophase prior to completion of the first round of transcription, suggesting that epigenetic features in gene bodies contribute to genome reconfiguration prior to transcription. The focus on the de novo formation of nuclear architecture during G1 entry yields insights into the contributions of CTCF and transcription to chromatin architecture dynamics during the mitosis to G1-phase progression.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

Engineering three-dimensional genome folding

Animal genomes are partitioned and folded at various scales that contribute distinctly to nuclear processes. While structural features have been disrupted either globally or at select loci in loss-of-function studies, gain-of-function studies that probe the role of genome architecture have lagged behind. Here we examine recent advances in experimentally creating chromatin loops, contact domains, boundaries and compartments. Furthermore, we explore parallels between this emerging theme and natural evolution of mammalian genomes with increasing architectural complexity. Finally, we provide a perspective on how insights arising from recent gain-of-function studies may inform future endeavors toward engineering the three-dimensional genome.

H3K27ac bookmarking promotes rapid post-mitotic activation of the pluripotent stem cell program without impacting 3D chromatin reorganization

During self-renewal, cell-type-defining features are drastically perturbed in mitosis and must be faithfully reestablished upon G1 entry, a process that remains largely elusive. Here, we characterized at a genome-wide scale the dynamic transcriptional and architectural resetting of mouse pluripotent stem cells (PSCs) upon mitotic exit. We captured distinct waves of transcriptional reactivation with rapid induction of stem cell genes and transient activation of lineage-specific genes. Topological reorganization at different hierarchical levels also occurred in an asynchronous manner and showed partial coordination with transcriptional resetting. Globally, rapid transcriptional and architectural resetting associated with mitotic retention of H3K27 acetylation, supporting a bookmarking function. Indeed, mitotic depletion of H3K27ac impaired the early reactivation of bookmarked, stem-cell-associated genes. However, 3D chromatin reorganization remained largely unaffected, suggesting that these processes are driven by distinct forces upon mitotic exit. This study uncovers principles and mediators of PSC molecular resetting during self-renewal.

Regulation of 3D chromatin organization by CTCF

Studies of nuclear architecture using chromosome conformation capture methods have provided a detailed view of how chromatin folds in the 3D nuclear space. New variants of this technology now afford unprecedented resolution and allow the identification of ever smaller folding domains that offer new insights into the mechanisms by which this organization is established and maintained. Here we review recent results in this rapidly evolving field with an emphasis on CTCF function, with the goal of gaining a mechanistic understanding of the principles by which chromatin is folded in the eukaryotic nucleus.

Engineering 3D genome organization

Cancers and developmental disorders are associated with alterations in the 3D genome architecture in space and time (the fourth dimension). Mammalian 3D genome organization is complex and dynamic and plays an essential role in regulating gene expression and cellular function. To study the causal relationship between genome function and its spatio-temporal organization in the nucleus, new technologies for engineering and manipulating the 3D organization of the genome have been developed. In particular, CRISPR-Cas technologies allow programmable manipulation at specific genomic loci, enabling unparalleled opportunities in this emerging field of 3D genome engineering. We review advances in mammalian 3D genome engineering with a focus on recent manipulative technologies using CRISPR-Cas and related technologies.

High-resolution targeted 3C interrogation of cis-regulatory element organization at genome-wide scale

Chromosome conformation capture (3C) provides an adaptable tool for studying diverse biological questions. Current 3C methods generally provide either low-resolution interaction profiles across the entire genome, or high-resolution interaction profiles at limited numbers of loci. Due to technical limitations, generation of reproducible high-resolution interaction profiles has not been achieved at genome-wide scale. Here, to overcome this barrier, we systematically test each step of 3C and report two improvements over current methods. We show that up to 30% of reporter events generated using the popular in situ 3C method arise from ligations between two individual nuclei, but this noise can be almost entirely eliminated by isolating intact nuclei after ligation. Using Nuclear-Titrated Capture-C, we generate reproducible high-resolution genome-wide 3C interaction profiles by targeting 8055 gene promoters in erythroid cells. By pairing high-resolution 3C interaction calls with nascent gene expression we interrogate the role of promoter hubs and super-enhancers in gene regulation.

Three-dimensional folding dynamics of the Xenopus tropicalis genome

Animal interphase chromosomes are organized into topologically associating domains (TADs). How TADs are formed is not fully understood. Here, we combined high-throughput chromosome conformation capture and gene silencing to obtain insights into TAD dynamics in Xenopus tropicalis embryos. First, TAD establishment in X. tropicalis is similar to that in mice and flies and does not depend on zygotic genome transcriptional activation. This process is followed by further refinements in active and repressive chromatin compartments and the appearance of loops and stripes. Second, within TADs, higher self-interaction frequencies at one end of the boundary are associated with higher DNA occupancy of the architectural proteins CTCF and Rad21. Third, the chromatin remodeling factor ISWI is required for de novo TAD formation. Finally, TAD structures are variable in different tissues. Our work shows that X. tropicalis is a powerful model for chromosome architecture analysis and suggests that chromatin remodeling plays an essential role in de novo TAD establishment.

Dynamic 3D Chromatin Reorganization during Establishment and Maintenance of Pluripotency

Higher-order chromatin structure is tightly linked to gene expression and therefore cell identity. In recent years, the chromatin landscape of pluripotent stem cells has become better characterized, and unique features at various architectural levels have been revealed. However, the mechanisms that govern establishment and maintenance of these topological characteristics and the temporal and functional relationships with transcriptional or epigenetic features are still areas of intense study. Here, we will discuss progress and limitations of our current understanding regarding how the 3D chromatin topology of pluripotent stem cells is established during somatic cell reprogramming and maintained during cell division. We will also discuss evidence and theories about the driving forces of topological reorganization and the functional links with key features and properties of pluripotent stem cell identity.