Replication dynamics identifies the folding principles of the inactive X chromosome.

Chromosome-wide late replication is an enigmatic hallmark of the inactive X chromosome (Xi). How it is established and what it represents remains obscure. By single-cell DNA replication sequencing, here we show that the entire Xi is reorganized to replicate rapidly and uniformly in late S-phase during X-chromosome inactivation (XCI), reflecting its relatively uniform structure revealed by 4C-seq. Despite this uniformity, only a subset of the Xi became earlier replicating in SmcHD1-mutant cells. In the mutant, these domains protruded out of the Xi core, contacted each other and became transcriptionally reactivated. 4C-seq suggested that they constituted the outermost layer of the Xi even before XCI and were rich in escape genes. We propose that this default positioning forms the basis for their inherent heterochromatin instability in cells lacking the Xi-binding protein SmcHD1 or exhibiting XCI escape. These observations underscore the importance of 3D genome organization for heterochromatin stability and gene regulation.

Epigenetic plasticity safeguards heterochromatin configuration in mammals.

Heterochromatin is a key architectural feature of eukaryotic chromosomes critical for cell type-specific gene expression and genome stability. In the mammalian nucleus, heterochromatin segregates from transcriptionally active genomic regions and exists in large, condensed, and inactive nuclear compartments. However, the mechanisms underlying the spatial organization of heterochromatin need to be better understood. Histone H3 lysine 9 trimethylation (H3K9me3) and lysine 27 trimethylation (H3K27me3) are two major epigenetic modifications that enrich constitutive and facultative heterochromatin, respectively. Mammals have at least five H3K9 methyltransferases (SUV39H1, SUV39H2, SETDB1, G9a and GLP) and two H3K27 methyltransferases (EZH1 and EZH2). In this study, we addressed the role of H3K9 and H3K27 methylation in heterochromatin organization using a combination of mutant cells for five H3K9 methyltransferases and an EZH1/2 dual inhibitor, DS3201. We showed that H3K27me3, which is normally segregated from H3K9me3, was redistributed to regions targeted by H3K9me3 after the loss of H3K9 methylation and that the loss of both H3K9 and H3K27 methylation resulted in impaired condensation and spatial organization of heterochromatin. Our data demonstrate that the H3K27me3 pathway safeguards heterochromatin organization after the loss of H3K9 methylation in mammalian cells.

Cell cycle dynamics and developmental dynamics of the 3D genome: toward linking the two timescales.

In the mammalian cell nucleus, chromosomes are folded differently in interphase and mitosis. Interphase chromosomes are relatively decondensed and display at least two unique layers of higher-order organization: topologically associating domains (TADs) and cell-type-specific A/B compartments, which correlate well with early/late DNA replication timing (RT). In mitosis, these structures rapidly disappear but are gradually reconstructed during G1 phase, coincident with the establishment of the RT program. However, these structures also change dynamically during cell differentiation and reprogramming, and yet we are surprisingly ignorant about the relationship between their cell cycle dynamics and developmental dynamics. In this review, we summarize the recent findings on this topic, discuss how these two processes might be coordinated with each other and its potential significance.

Highly rigid H3.1/H3.2-H3K9me3 domains set a barrier for cell fate reprogramming in trophoblast stem cells.

The placenta is a highly evolved, specialized organ in mammals. It differs from other organs in that it functions only for fetal maintenance during gestation. Therefore, there must be intrinsic mechanisms that guarantee its unique functions. To address this question, we comprehensively analyzed epigenomic features of mouse trophoblast stem cells (TSCs). Our genome-wide, high-throughput analyses revealed that the TSC genome contains large-scale (>1-Mb) rigid heterochromatin architectures with a high degree of histone H3.1/3.2-H3K9me3 accumulation, which we termed TSC-defined highly heterochromatinized domains (THDs). Importantly, depletion of THDs by knockdown of CAF1, an H3.1/3.2 chaperone, resulted in down-regulation of TSC markers, such as Cdx2 and Elf5, and up-regulation of the pluripotent marker Oct3/4, indicating that THDs maintain the trophoblastic nature of TSCs. Furthermore, our nuclear transfer technique revealed that THDs are highly resistant to genomic reprogramming. However, when H3K9me3 was removed, the TSC genome was fully reprogrammed, giving rise to the first TSC cloned offspring. Interestingly, THD-like domains are also present in mouse and human placental cells in vivo, but not in other cell types. Thus, THDs are genomic architectures uniquely developed in placental lineage cells, which serve to protect them from fate reprogramming to stably maintain placental function.

SAF-A promotes origin licensing and replication fork progression to ensure robust DNA replication.

The organisation of chromatin is closely intertwined with biological activities of chromosome domains, including transcription and DNA replication status. Scaffold-attachment factor A (SAF-A), also known as heterogeneous nuclear ribonucleoprotein U (HNRNPU), contributes to the formation of open chromatin structure. Here, we demonstrate that SAF-A promotes the normal progression of DNA replication and enables resumption of replication after inhibition. We report that cells depleted of SAF-A show reduced origin licensing in G1 phase and, consequently, reduced origin activation frequency in S phase. Replication forks also progress less consistently in cells depleted of SAF-A, contributing to reduced DNA synthesis rate. Single-cell replication timing analysis revealed two distinct effects of SAF-A depletion: first, the boundaries between early- and late-replicating domains become more blurred; and second, SAF-A depletion causes replication timing changes that tend to bring regions of discordant domain compartmentalisation and replication timing into concordance. Associated with these defects, SAF-A-depleted cells show elevated formation of phosphorylated histone H2AX (γ-H2AX) and tend to enter quiescence. Overall, we find that SAF-A protein promotes robust DNA replication to ensure continuing cell proliferation.

Regulation of mammalian 3D genome organization and histone H3K9 dimethylation by H3K9 methyltransferases.

Histone H3 lysine 9 dimethylation (H3K9me2) is a highly conserved silencing epigenetic mark. Chromatin marked with H3K9me2 forms large domains in mammalian cells and overlaps well with lamina-associated domains and the B compartment defined by Hi-C. However, the role of H3K9me2 in 3-dimensional (3D) genome organization remains unclear. Here, we investigated genome-wide H3K9me2 distribution, transcriptome, and 3D genome organization in mouse embryonic stem cells following the inhibition or depletion of H3K9 methyltransferases (MTases): G9a, GLP, SETDB1, SUV39H1, and SUV39H2. We show that H3K9me2 is regulated by all five MTases; however, H3K9me2 and transcription in the A and B compartments are regulated by different MTases. H3K9me2 in the A compartments is primarily regulated by G9a/GLP and SETDB1, while H3K9me2 in the B compartments is regulated by all five MTases. Furthermore, decreased H3K9me2 correlates with changes to more active compartmental state that accompanied transcriptional activation. Thus, H3K9me2 contributes to inactive compartment setting.

Mapping replication timing domains genome wide in single mammalian cells with single-cell DNA replication sequencing.

Replication timing (RT) domains are stable units of chromosome structure that are regulated in the context of development and disease. Conventional genome-wide RT mapping methods require many S-phase cells for either the effective enrichment of replicating DNA through bromodeoxyuridine (BrdU) immunoprecipitation or the determination of copy-number differences during S-phase, which precludes their application to non-abundant cell types and single cells. Here, we provide a simple, cost-effective, and robust protocol for single-cell DNA replication sequencing (scRepli-seq). The scRepli-seq methodology relies on whole-genome amplification (WGA) of genomic DNA (gDNA) from single S-phase cells and next-generation sequencing (NGS)-based determination of copy-number differences that arise between replicated and unreplicated DNA. Haplotype-resolved scRepli-seq, which distinguishes pairs of homologous chromosomes within a single cell, is feasible by using single-nucleotide polymorphism (SNP)/indel information. We also provide computational pipelines for quality control, normalization, and binarization of the scRepli-seq data. The experimental portion of this protocol (before sequencing) takes 3 d.

Microrheology for Hi-C Data Reveals the Spectrum of the Dynamic 3D Genome Organization.

The one-dimensional information of genomic DNA is hierarchically packed inside the eukaryotic cell nucleus and organized in a three-dimensional (3D) space. Genome-wide chromosome conformation capture (Hi-C) methods have uncovered the 3D genome organization and revealed multiscale chromatin domains of compartments and topologically associating domains (TADs). Moreover, single-nucleosome live-cell imaging experiments have revealed the dynamic organization of chromatin domains caused by stochastic thermal fluctuations. However, the mechanism underlying the dynamic regulation of such hierarchical and structural chromatin units within the microscale thermal medium remains unclear. Microrheology is a way to measure dynamic viscoelastic properties coupling between thermal microenvironment and mechanical response. Here, we propose a new, to our knowledge, microrheology for Hi-C data to analyze the dynamic compliance property as a measure of rigidness and flexibility of genomic regions along with the time evolution. Our method allows the conversion of an Hi-C matrix into the spectrum of the dynamic rheological property along the genomic coordinate of a single chromosome. To demonstrate the power of the technique, we analyzed Hi-C data during the neural differentiation of mouse embryonic stem cells. We found that TAD boundaries behave as more rigid nodes than the intra-TAD regions. The spectrum clearly shows the dynamic viscoelasticity of chromatin domain formation at different timescales. Furthermore, we characterized the appearance of synchronous and liquid-like intercompartment interactions in differentiated cells. Together, our microrheology data derived from Hi-C data provide physical insights into the dynamics of the 3D genome organization.

Multifaceted Hi-C benchmarking: what makes a difference in chromosome-scale genome scaffolding?

BACKGROUND: Hi-C is derived from chromosome conformation capture (3C) and targets chromatin contacts on a genomic scale. This method has also been used frequently in scaffolding nucleotide sequences obtained by de novo genome sequencing and assembly, in which the number of resultant sequences rarely converges to the chromosome number. Despite its prevalent use, the sample preparation methods for Hi-C have not been intensively discussed, especially from the standpoint of genome scaffolding. RESULTS: To gain insight into the best practice of Hi-C scaffolding, we performed a multifaceted methodological comparison using vertebrate samples and optimized various factors during sample preparation, sequencing, and computation. As a result, we identified several key factors that helped improve Hi-C scaffolding, including the choice and preparation of tissues, library preparation conditions, the choice of restriction enzyme(s), and the choice of scaffolding program and its usage. CONCLUSIONS: This study provides the first comparison of multiple sample preparation kits/protocols and computational programs for Hi-C scaffolding by an academic third party. We introduce a customized protocol designated "inexpensive and controllable Hi-C (iconHi-C) protocol," which incorporates the optimal conditions identified in this study, and demonstrate this technique on chromosome-scale genome sequences of the Chinese softshell turtle Pelodiscus sinensis.

The Eleanor ncRNAs activate the topological domain of the ESR1 locus to balance against apoptosis.

MCF7 cells acquire estrogen-independent proliferation after long-term estrogen deprivation (LTED), which recapitulates endocrine therapy resistance. LTED cells can become primed for apoptosis, but the underlying mechanism is largely unknown. We previously reported that Eleanor non-coding RNAs (ncRNAs) upregulate the ESR1 gene in LTED cells. Here, we show that Eleanors delineate the topologically associating domain (TAD) of the ESR1 locus in the active nuclear compartment of LTED cells. The TAD interacts with another transcriptionally active TAD, which is 42.9 Mb away from ESR1 and contains a gene encoding the apoptotic transcription factor FOXO3. Inhibition of a promoter-associated Eleanor suppresses all genes inside the Eleanor TAD and the long-range interaction between the two TADs, but keeps FOXO3 active to facilitate apoptosis in LTED cells. These data indicate a role of ncRNAs in chromatin domain regulation, which may underlie the apoptosis-prone nature of therapy-resistant breast cancer cells and could be good therapeutic targets.

Single-cell DNA replication profiling identifies spatiotemporal developmental dynamics of chromosome organization.

In mammalian cells, chromosomes are partitioned into megabase-sized topologically associating domains (TADs). TADs can be in either A (active) or B (inactive) subnuclear compartments, which exhibit early and late replication timing (RT), respectively. Here, we show that A/B compartments change coordinately with RT changes genome wide during mouse embryonic stem cell (mESC) differentiation. While A to B compartment changes and early to late RT changes were temporally inseparable, B to A changes clearly preceded late to early RT changes and transcriptional activation. Compartments changed primarily by boundary shifting, altering the compartmentalization of TADs facing the A/B compartment interface, which was conserved during reprogramming and confirmed in individual cells by single-cell Repli-seq. Differentiating mESCs altered single-cell Repli-seq profiles gradually but uniformly, transiently resembling RT profiles of epiblast-derived stem cells (EpiSCs), suggesting that A/B compartments might also change gradually but uniformly toward a primed pluripotent state. These results provide insights into how megabase-scale chromosome organization changes in individual cells during differentiation.

Genome-wide stability of the DNA replication program in single mammalian cells.

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a genome-wide methodology that measures copy number differences between replicated and unreplicated DNA. Using scRepli-seq, we demonstrate that replication-domain organization is conserved among individual mouse embryonic stem cells (mESCs). Differentiated mESCs exhibited distinct profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication profiles of homologous autosomes, while the inactive X chromosome was clearly replicated later than its active counterpart. However, a small degree of cell-to-cell replication-timing heterogeneity was present, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

Practical Analysis of Hi-C Data: Generating A/B Compartment Profiles.

Recent advances in next-generation sequencing (NGS) and chromosome conformation capture (3C) analysis have led to the development of Hi-C, a genome-wide version of the 3C method. Hi-C has identified new levels of chromosome organization such as A/B compartments, topologically associating domains (TADs) as well as large megadomains on the inactive X chromosome, while allowing the identification of chromatin loops at the genome scale. Despite its powerfulness, Hi-C data analysis is much more involved compared to conventional NGS applications such as RNA-seq or ChIP-seq and requires many more steps. This presents a significant hurdle for those who wish to implement Hi-C technology into their laboratory. On the other hand, genomics data repository sites sometimes contain processed Hi-C data sets, allowing researchers to perform further analysis without the need for high-spec workstations and servers. In this chapter, we provide a detailed description on how to calculate A/B compartment profiles from processed Hi-C data on the autosomes and the active/inactive X chromosomes.

Srf destabilizes cellular identity by suppressing cell-type-specific gene expression programs.

Multicellular organisms consist of multiple cell types. The identity of these cells is primarily maintained by cell-type-specific gene expression programs; however, mechanisms that suppress these programs are poorly defined. Here we show that serum response factor (Srf), a transcription factor that is activated by various extracellular stimuli, can repress cell-type-specific genes and promote cellular reprogramming to pluripotency. Manipulations that decrease ß-actin monomer quantity result in the nuclear accumulation of Mkl1 and the activation of Srf, which downregulate cell-type-specific genes and alter the epigenetics of regulatory regions and chromatin organization. Mice overexpressing Srf exhibit various pathologies including an ulcerative colitis-like symptom and a metaplasia-like phenotype in the pancreas. Our results demonstrate an unexpected function of Srf via a mechanism by which extracellular stimuli actively destabilize cell identity and suggest Srf involvement in a wide range of diseases.

Adapting CRISPR/Cas9 for functional genomics screens.

The use of CRISPR/Cas9 (clustered regularly interspaced short palindromic repeats/CRISPR-associated protein) for targeted genome editing has been widely adopted and is considered a "game changing" technology. The ease and rapidity by which this approach can be used to modify endogenous loci in a wide spectrum of cell types and organisms makes it a powerful tool for customizable genetic modifications as well as for large-scale functional genomics. The development of retrovirus-based expression platforms to simultaneously deliver the Cas9 nuclease and single guide (sg) RNAs provides unique opportunities by which to ensure stable and reproducible expression of the editing tools and a broad cell targeting spectrum, while remaining compatible with in vivo genetic screens. Here, we describe methods and highlight considerations for designing and generating sgRNA libraries in all-in-one retroviral vectors for such applications.

Protospacer adjacent motif (PAM)-distal sequences engage CRISPR Cas9 DNA target cleavage.

The clustered regularly interspaced short palindromic repeat (CRISPR)-associated enzyme Cas9 is an RNA-guided nuclease that has been widely adapted for genome editing in eukaryotic cells. However, the in vivo target specificity of Cas9 is poorly understood and most studies rely on in silico predictions to define the potential off-target editing spectrum. Using chromatin immunoprecipitation followed by sequencing (ChIP-seq), we delineate the genome-wide binding panorama of catalytically inactive Cas9 directed by two different single guide (sg) RNAs targeting the Trp53 locus. Cas9:sgRNA complexes are able to load onto multiple sites with short seed regions adjacent to (5')NGG(3') protospacer adjacent motifs (PAM). Yet among 43 ChIP-seq sites harboring seed regions analyzed for mutational status, we find editing only at the intended on-target locus and one off-target site. In vitro analysis of target site recognition revealed that interactions between the 5' end of the guide and PAM-distal target sequences are necessary to efficiently engage Cas9 nucleolytic activity, providing an explanation for why off-target editing is significantly lower than expected from ChIP-seq data.

Classifying leukemia types with chromatin conformation data.

BACKGROUND: Although genetic or epigenetic alterations have been shown to affect the three-dimensional organization of genomes, the utility of chromatin conformation in the classification of human disease has never been addressed. RESULTS: Here, we explore whether chromatin conformation can be used to classify human leukemia. We map the conformation of the HOXA gene cluster in a panel of cell lines with 5C chromosome conformation capture technology, and use the data to train and test a support vector machine classifier named 3D-SP. We show that 3D-SP is able to accurately distinguish leukemias expressing MLL-fusion proteins from those expressing only wild-type MLL, and that it can also classify leukemia subtypes according to MLL fusion partner, based solely on 5C data. CONCLUSIONS: Our study provides the first proof-of-principle demonstration that chromatin conformation contains the information value necessary for classification of leukemia subtypes.

Hox in motion: tracking HoxA cluster conformation during differentiation.

Three-dimensional genome organization is an important higher order transcription regulation mechanism that can be studied with the chromosome conformation capture techniques. Here, we combined chromatin organization analysis by chromosome conformation capture-carbon copy, computational modeling and epigenomics to achieve the first integrated view, through time, of a connection between chromatin state and its architecture. We used this approach to examine the chromatin dynamics of the HoxA cluster in a human myeloid leukemia cell line at various stages of differentiation. We found that cellular differentiation involves a transient activation of the 5'-end HoxA genes coinciding with a loss of contacts throughout the cluster, and by specific silencing at the 3'-end with H3K27 methylation. The 3D modeling of the data revealed an extensive reorganization of the cluster between the two previously reported topologically associated domains in differentiated cells. Our results support a model whereby silencing by polycomb group proteins and reconfiguration of CTCF interactions at a topologically associated domain boundary participate in changing the HoxA cluster topology, which compartmentalizes the genes following differentiation.

A Torrent of data: mapping chromatin organization using 5C and high-throughput sequencing.

The study of three-dimensional genome organization is an exciting research area, which has benefited from the rapid development of high-resolution molecular mapping techniques over the past decade. These methods are derived from the chromosome conformation capture (3C) technique and are each aimed at improving some aspect of 3C. All 3C technologies use formaldehyde fixation and proximity-based ligation to capture chromatin contacts in cell populations and consider in vivo spatial proximity more or less inversely proportional to the frequency of measured interactions. The 3C-carbon copy (5C) method is among the most quantitative of these approaches. 5C is extremely robust and can be used to study chromatin organization at various scales. Here, we present a modified 5C analysis protocol adapted for sequencing with an Ion Torrent Personal Genome Machine™ (PGM™). We explain how Torrent 5C libraries are produced and sequenced. We also describe the statistical and computational methods we developed to normalize and analyze raw Torrent 5C sequence data. The Torrent 5C protocol should facilitate the study of in vivo chromatin architecture at high resolution because it benefits from high accuracy, greater speed, low running costs, and the flexibility of in-house next-generation sequencing.

Discovering genome regulation with 3C and 3C-related technologies.

It has been known for some time that eukaryotic genomic DNA is packaged in the form of highly organized chromatin in vivo. This organization is important not only to reduce the length of chromosomes during interphase but also because it represents a type of higher-order genome regulation mechanism. Indeed, spatial chromatin architecture is known to be important for transcription, DNA replication and repair. Chromosome structure can be observed at different scales and studied with a variety of complementary techniques. For example, microscopy can provide single cell information while technologies such as the chromosome conformation capture (3C) method and its derivatives can yield higher-resolution data from cell populations. In this review, we report on the biological questions addressed with 3C and 3C-related techniques and what has been uncovered to date. We also explore what these methods may further reveal about the regulation of genomic DNA activities.