Search

Abstract Background Eukaryotic DNA replication is regulated at the level of large chromosomal domains (0.5–5 megabases in mammals) within which replicons are activated relatively synchronously. These domains replicate in a specific temporal order during S-phase and our genome-wide analyses of replication timing have demonstrated that this temporal order of domain replication is a stable property of specific cell types. Results We have developed ReplicationDomain http://www.replicationdomain.org as a web-based database for analysis of genome-wide replication timing maps (replication profiles) from various cell lines and species. This database also provides comparative information of transcriptional expression and is configured to display any genome-wide property (for instance, ChIP-Chip or ChIP-Seq data) via an interactive web interface. Our published microarray data sets are publicly available. Users may graphically display these data sets for a selected genomic region and download the data displayed as text files, or alternatively, download complete genome-wide data sets. Furthermore, we have implemented a user registration system that allows registered users to upload their own data sets. Upon uploading, registered users may choose to: (1) view their data sets privately without sharing; (2) share with other registered users; or (3) make their published or "in press" data sets publicly available, which can fulfill journal and funding agencies' requirements for data sharing. Conclusion ReplicationDomain is a novel and powerful tool to facilitate the comparative visualization of replication timing in various cell types as well as other genome-wide chromatin features and is considerably faster and more convenient than existing browsers when viewing multi-megabase segments of chromosomes. Furthermore, the data upload function with the option of private viewing or sharing of data sets between registered users should be a valuable resource for the scientific community.

Mammalian chromosomes form a hierarchical structure within the cell nucleus, from chromatin loops, megabase (Mb)-sized topologically associating domains (TADs) to larger-scale A/B compartments. The molecular basis of the structures of loops and TADs has been actively studied. However, the A and B compartments, which correspond to early-replicating euchromatin and late-replicating heterochromatin, respectively, are still relatively unexplored. In this review, we focus on the A/B compartments, discuss their close relationship to DNA replication timing (RT), and introduce recent findings on the features of subcompartments revealed by detailed classification of the A/B compartments. In doing so, we speculate on the structure, potential function, and developmental dynamics of A/B compartments and subcompartments in mammalian cells., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Author(s). Published by Elsevier Ltd.. All rights reserved.)

Faithful DNA replication is essential for genome integrity 1-4 . Under-replicated DNA leads to defects in chromosome segregation, which are common during embryogenesis 5-8 . However, the regulation of DNA replication remains poorly understood in early mammalian embryos. Here we constructed a single-cell genome-wide DNA replication atlas of pre-implantation mouse embryos and identified an abrupt replication program switch accompanied by a transient period of genomic instability. In 1- and 2-cell embryos, we observed the complete absence of a replication timing program, and the entire genome replicated gradually and uniformly using extremely slow-moving replication forks. In 4-cell embryos, a somatic-cell-like replication timing program commenced abruptly. However, the fork speed was still slow, S phase was extended, and markers of replication stress, DNA damage and repair increased. This was followed by an increase in break-type chromosome segregation errors specifically during the 4-to-8-cell division with breakpoints enriched in late-replicating regions. These errors were rescued by nucleoside supplementation, which accelerated fork speed and reduced the replication stress. By the 8-cell stage, forks gained speed, S phase was no longer extended and chromosome aberrations decreased. Thus, a transient period of genomic instability exists during normal mouse development, preceded by an S phase lacking coordination between replisome-level regulation and megabase-scale replication timing regulation, implicating a link between their coordination and genome stability., (© 2024. The Author(s).)

Over a decade has passed since the development of the Hi-C method for genome-wide analysis of 3D genome organization. Hi-C utilizes next-generation sequencing (NGS) technology to generate large-scale chromatin interaction data, which has accumulated across a diverse range of species and cell types, particularly in eukaryotes. There is thus a growing need to streamline the process of Hi-C data analysis to utilize these data sets effectively. Hi-C generates data that are much larger compared to other NGS techniques such as chromatin immunoprecipitation sequencing (ChIP-seq) or RNA-seq, making the data reanalysis process computationally expensive. In an effort to bridge this resource gap, the 4D Nucleome (4DN) Data Portal has reanalyzed approximately 600 Hi-C data sets, allowing users to access and utilize the analyzed data. In this chapter, we provide detailed instructions for the implementation of the common workflow language (CWL)-based Hi-C analysis pipeline adopted by the 4DN Data Portal ecosystem. This reproducible and portable pipeline generates standard Hi-C contact matrices in formats such as .hic or .mcool from FASTQ files. It enables users to output their own Hi-C data in the same format as those registered in the 4DN Data portal, facilitating comparative analysis using data registered in the portal. Our custom-made scripts are available on GitHub at https://github.com/kuzobuta/4dn_cwl_pipeline ., (© 2025. The Author(s), under exclusive license to Springer Science+Business Media, LLC, part of Springer Nature.)

Cell death and proliferation are at a glance dichotomic events, but occasionally coupled. Caspases, traditionally known to execute apoptosis, play non-apoptotic roles, but their exact mechanism remains elusive. Here, using Drosophila intestinal stem cells (ISCs), we discovered that activation of caspases induces massive cell proliferation rather than cell death. We elucidate that a positive feedback circuit exists between caspases and JNK, which can simultaneously drive cell proliferation and cell death. In ISCs, signalling from JNK to caspases is defective, which skews the balance towards proliferation. Mechanistically, two-tiered regulation of the DIAP1 inhibitor rpr , through its transcription and its protein localization, exists. This work provides a conceptual framework that explains how caspases perform apoptotic and non-apoptotic functions in vivo and how ISCs accomplish their resistance to cell death., (© 2024 Sulekh et al.)

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., (© 2023. The Author(s).)

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., (© The Author(s) 2023. Published by Oxford University Press on behalf of Nucleic Acids Research.)

Epigenetic abnormalities are extremely widespread in cancer. Some of them are mere consequences of transformation, but some actively contribute to cancer initiation and progression; they provide powerful new biological markers, as well as new targets for therapies. In this review, we examine the recent literature and focus on one particular aspect of epigenome deregulation: large-scale chromatin changes, causing global changes of DNA methylation or histone modifications. After a brief overview of the one-dimension (1D) and three-dimension (3D) epigenome in healthy cells and of its homeostasis mechanisms, we use selected examples to describe how many different events (mutations, changes in metabolism, and infections) can cause profound changes to the epigenome and fuel cancer. We then present the consequences for therapies and briefly discuss the role of single-cell approaches for the future progress of the field.

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., (Copyright © 2022 The Authors. Published by Elsevier Ltd.. All rights reserved.)

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., Competing Interests: Competing interests The authors declare no competing or financial interests., (© 2022. Published by The Company of Biologists Ltd.)

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., (© 2022 Hada et al.; Published by Cold Spring Harbor Laboratory Press.)

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.

The fine-scale dynamics from euchromatin (EC) to facultative heterochromatin (fHC) has remained largely unclear. Here, we focus on Xist and its silencing initiator Tsix as a paradigm of transcription-mediated conversion from EC to fHC. In mouse epiblast stem cells, induction of Tsix recapitulates the conversion at the Xist promoter. Investigating the dynamics reveals that the conversion proceeds in a stepwise manner. Initially, a transient opened chromatin structure is observed. In the second step, gene silencing is initiated and dependent on Tsix, which is reversible and accompanied by simultaneous changes in multiple histone modifications. At the last step, maintenance of silencing becomes independent of Tsix and irreversible, which correlates with occupation of the -1 position of the transcription start site by a nucleosome and initiation of DNA methylation introduction. This study highlights the hierarchy of multiple chromatin events upon stepwise gene silencing establishment., Competing Interests: Declaration of interests The authors declare no competing interests., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Multiple epigenetic pathways underlie the temporal order of DNA replication (replication timing) in the contexts of development and disease. DNA methylation by DNA methyltransferases (Dnmts) and downstream chromatin reorganization and transcriptional changes are thought to impact DNA replication, yet this remains to be comprehensively tested. Using cell-based and genome-wide approaches to measure replication timing, we identified a number of genomic regions undergoing subtle but reproducible replication timing changes in various Dnmt -mutant mouse embryonic stem (ES) cell lines that included a cell line with a drug-inducible Dnmt3a2 expression system. Replication timing within pericentromeric heterochromatin (PH) was shown to be correlated with redistribution of H3K27me3 induced by DNA hypomethylation: Later replicating PH coincided with H3K27me3-enriched regions. In contrast, this relationship with H3K27me3 was not evident within chromosomal arm regions undergoing either early-to-late (EtoL) or late-to-early (LtoE) switching of replication timing upon loss of the Dnmts. Interestingly, Dnmt-sensitive transcriptional up- and downregulation frequently coincided with earlier and later shifts in replication timing of the chromosomal arm regions, respectively. Our study revealed the previously unrecognized complex and diverse effects of the Dnmts loss on the mammalian DNA replication landscape.

During embryogenesis in mammals, the 3-dimensional (3D) genome organization changes globally in parallel with transcription changes in a cell-type specific manner. This involves the progressive formation of heterochromatin, the best example of which is the inactive X chromosome (Xi) in females, originally discovered as a compact 3D structure at the nuclear periphery known as the Barr body. The heterochromatin formation on the autosomes and the Xi is tightly associated with the differentiation state and the developmental potential of cells, making it an ideal readout of the cellular epigenetic state. At a glance, the heterochromatin appears to be uniform. However, recent studies are beginning to reveal a more complex picture, with multiple hierarchical levels co-existing within the heterochromatin compartment. Such hierarchical levels appear to exist in the heterochromatin compartment on autosomes as well as on the Xi. Here, we review recent progress in our understanding of the 3D genome organization changes during the period of differentiation surrounding pluripotency in vivo and in vitro, with a focus on the heterochromatin compartment. We first look at the whole genome, then focus on the Xi, and discuss their differences and similarities. Finally, we present a unified view of how the heterochromatin compartment is formed and regulated during early development. In particular, we emphasize that there are multiple layers within the heterochromatic compartment on both the autosomes and the Xi, with regulatory mechanisms common and specific to each layer., (© 2021 Japanese Society of Developmental Biologists.)

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.

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., (Copyright © 2020 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

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 Author(s) 2020. Published by Oxford University Press.)

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.

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.

In mammalian cells, DNA replication timing is controlled at the level of megabase (Mb)-sized chromosomal domains and correlates well with transcription, chromatin structure, and three-dimensional (3D) genome organization. Because of these properties, DNA replication timing is an excellent entry point to explore genome regulation at various levels and a variety of studies have been carried out over the years. However, DNA replication timing studies traditionally required at least tens of thousands of cells, and it was unclear whether the replication domains detected by cell population analyses were preserved at the single-cell level. Recently, single-cell DNA replication profiling methods became available, which revealed that the Mb-sized replication domains detected by cell population analyses were actually well preserved in individual cells. In this article, we provide a brief overview of our current knowledge on DNA replication timing regulation in mammals based on cell population studies, outline the findings from single-cell DNA replication profiling, and discuss future directions and challenges.

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.

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.

It has been 8 years since the concept of naïve and primed pluripotent stem cell states was first proposed. Both are states of pluripotency, but exhibit slightly different properties. The naïve state represents the cellular state of the preimplantation mouse blastocyst inner cell mass, while the primed state is representative of the post-implantation epiblast cells. These two cell types exhibit clearly distinct developmental potential, as evidenced by the fact that naïve cells are able to contribute to blastocyst chimeras, while primed cells cannot. However, the epigenetic differences that underlie the distinct developmental potential of these cell types remain unclear, which is rather surprising given the large amount of active investigation over the years. Elucidating such epigenetic differences should lead to a better understanding of the fundamental properties of these states of pluripotency and the means by which the naïve-to-primed transition occurs, which may provide insights into the essence of stem cell commitment.

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.

Chromatin DNA must be read out for various cellular functions, and copied for the next cell division. These processes are targets of many anticancer agents. Platinum-based drugs, such as cisplatin, have been used extensively in cancer chemotherapy. The drug-DNA interaction causes DNA crosslinks and subsequent cytotoxicity. Recently, it was reported that an azolato-bridged dinuclear platinum(II) complex, 5-H-Y, exhibits a different anticancer spectrum from cisplatin. Here, using an interdisciplinary approach, we reveal that the cytotoxic mechanism of 5-H-Y is distinct from that of cisplatin. 5-H-Y inhibits DNA replication and also RNA transcription, arresting cells in the S/G2 phase, and are effective against cisplatin-resistant cancer cells. Moreover, it causes much less DNA crosslinking than cisplatin, and induces chromatin folding. 5-H-Y will expand the clinical applications for the treatment of chemotherapy-insensitive cancers.

Many types of epigenetic profiling have been used to classify stem cells, stages of cellular differentiation, and cancer subtypes. Existing methods focus on local chromatin features such as DNA methylation and histone modifications that require extensive analysis for genome-wide coverage. Replication timing has emerged as a highly stable cell type-specific epigenetic feature that is regulated at the megabase-level and is easily and comprehensively analyzed genome-wide. Here, we describe a cell classification method using 67 individual replication profiles from 34 mouse and human cell lines and stem cell-derived tissues, including new data for mesendoderm, definitive endoderm, mesoderm and smooth muscle. Using a Monte-Carlo approach for selecting features of replication profiles conserved in each cell type, we identify "replication timing fingerprints" unique to each cell type and apply a k nearest neighbor approach to predict known and unknown cell types. Our method correctly classifies 67/67 independent replication-timing profiles, including those derived from closely related intermediate stages. We also apply this method to derive fingerprints for pluripotency in human and mouse cells. Interestingly, the mouse pluripotency fingerprint overlaps almost completely with previously identified genomic segments that switch from early to late replication as pluripotency is lost. Thereafter, replication timing and transcription within these regions become difficult to reprogram back to pluripotency, suggesting these regions highlight an epigenetic barrier to reprogramming. In addition, the major histone cluster Hist1 consistently becomes later replicating in committed cell types, and several histone H1 genes in this cluster are downregulated during differentiation, suggesting a possible instrument for the chromatin compaction observed during differentiation. Finally, we demonstrate that unknown samples can be classified independently using site-specific PCR against fingerprint regions. In sum, replication fingerprints provide a comprehensive means for cell characterization and are a promising tool for identifying regions with cell type-specific organization.

Replication timing profiles are cell type-specific and reflect genome organization changes during differentiation. In this protocol, we describe how to analyze genome-wide replication timing (RT) in mammalian cells. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU) and sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated, amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray, which is currently more cost effective than high-throughput sequencing and equally capable of resolving features at the biologically relevant level of tens to hundreds of kilobases. We also present a guide to analyzing the resulting data sets based on methods we use routinely. Subjects include normalization, scaling and data quality measures, LOESS (local polynomial) smoothing of RT values, segmentation of data into domains and assignment of timing values to gene promoters. Finally, we cover clustering methods and means to relate changes in the replication program to gene expression and other genetic and epigenetic data sets. Some experience with R or similar programming languages is assumed. All together, the protocol takes ∼3 weeks per batch of samples.

Facioscapulohumeral muscular dystrophy (FSHD) is linked to contraction of an array of tandem 3.3-kb repeats (D4Z4) at 4q35.2 from 11-100 copies to 1-10 copies. The extent to which D4Z4 contraction at 4q35.2 affects overall 4q35.2 chromatin organization remains unclear. Because DNA replication timing is highly predictive of long-range chromatin interactions, we generated genome-wide replication-timing profiles for FSHD and control myogenic precursor cells. We compared non-immortalized myoblasts from four FSHD patients and three control individuals to each other and to a variety of other human cell types. This study also represents the first genome-wide comparison of replication timing profiles in non-immortalized human cell cultures. Myoblasts from both control and FSHD individuals all shared a myoblast-specific replication profile. In contrast, male and female individuals were readily distinguished by monoallelic differences in replication timing at DXZ4 and other regions across the X chromosome affected by X inactivation. We conclude that replication timing is a robust cell-type specific feature that is unaffected by FSHD-related D4Z4 contraction.

To identify evolutionarily conserved features of replication timing and their relationship to epigenetic properties, we profiled replication timing genome-wide in four human embryonic stem cell (hESC) lines, hESC-derived neural precursor cells (NPCs), lymphoblastoid cells, and two human induced pluripotent stem cell lines (hiPSCs), and compared them with related mouse cell types. Results confirm the conservation of coordinately replicated megabase-sized "replication domains" punctuated by origin-suppressed regions. Differentiation-induced replication timing changes in both species occur in 400- to 800-kb units and are similarly coordinated with transcription changes. A surprising degree of cell-type-specific conservation in replication timing was observed across regions of conserved synteny, despite considerable species variation in the alignment of replication timing to isochore GC/LINE-1 content. Notably, hESC replication timing profiles were significantly more aligned to mouse epiblast-derived stem cells (mEpiSCs) than to mouse ESCs. Comparison with epigenetic marks revealed a signature of chromatin modifications at the boundaries of early replicating domains and a remarkably strong link between replication timing and spatial proximity of chromatin as measured by Hi-C analysis. Thus, early and late initiation of replication occurs in spatially separate nuclear compartments, but rarely within the intervening chromatin. Moreover, cell-type-specific conservation of the replication program implies conserved developmental changes in spatial organization of chromatin. Together, our results reveal evolutionarily conserved aspects of developmentally regulated replication programs in mammals, demonstrate the power of replication profiling to distinguish closely related cell types, and strongly support the hypothesis that replication timing domains are spatially compartmentalized structural and functional units of three-dimensional chromosomal architecture.

Differentiation of mouse embryonic stem cells (mESCs) is accompanied by changes in replication timing. To explore the relationship between replication timing and cell fate transitions, we constructed genome-wide replication-timing profiles of 22 independent mouse cell lines representing 10 stages of early mouse development, and transcription profiles for seven of these stages. Replication profiles were cell-type specific, with 45% of the genome exhibiting significant changes at some point during development that were generally coordinated with changes in transcription. Comparison of early and late epiblast cell culture models revealed a set of early-to-late replication switches completed at a stage equivalent to the post-implantation epiblast, prior to germ layer specification and down-regulation of key pluripotency transcription factors [POU5F1 (also known as OCT4)/NANOG/SOX2] and coinciding with the emergence of compact chromatin near the nuclear periphery. These changes were maintained in all subsequent lineages (lineage-independent) and involved a group of irreversibly down-regulated genes, at least some of which were repositioned closer to the nuclear periphery. Importantly, many genomic regions of partially reprogrammed induced pluripotent stem cells (piPSCs) failed to re-establish ESC-specific replication-timing and transcription programs. These regions were enriched for lineage-independent early-to-late changes, which in female cells included the inactive X chromosome. Together, these results constitute a comprehensive "fate map" of replication-timing changes during early mouse development. Moreover, they support a model in which a distinct set of replication domains undergoes a form of "autosomal Lyonization" in the epiblast that is difficult to reprogram and coincides with an epigenetic commitment to differentiation prior to germ layer specification.

Studies of replication timing provide a handle into previously impenetrable higher-order levels of chromosome organization and their plasticity during development. Although mechanisms regulating replication timing are not clear, novel genome-wide studies provide a thorough survey of the extent to which replication timing is regulated during most of the early cell fate transitions in mammals, revealing coordinated changes of a defined set of 400-800 kb chromosomal segments that involve at least half the genome. Furthermore, changes in replication time are linked to changes in sub-nuclear organization and domain-wide transcriptional potential, and tissue-specific replication timing profiles are conserved from mouse to human, suggesting that the program has developmental significance. Hence, these studies have provided a solid foundation for linking megabase level chromosome structure to function, and suggest a central role for replication in domain-level genome organization.

It has been exactly 50 years since it was discovered that duplication of the eukaryotic genome follows a defined temporal order as cells progress through S-phase. While the mechanism of this replication-timing program still remains a mystery, various correlations of this program with both static and dynamic properties of chromatin render it an attractive forum to explore previously impenetrable higher-order organization of chromosomes. Indeed, studies of DNA replication have provided a simple and straightforward approach to address physical organization of the genome, both along the length of the chromosome as well as in the context of the 3-dimensional space in the cell nucleus. In this chapter, we summarize the 50-years history of the pursuit for understanding the replication-timing program and its developmental regulation, primarily in mammalian cells. We begin with the discovery of the replication-timing program, discuss developmental regulation of this program during X-inactivation in females as well as on autosomes and then describe the recent findings from genome-wide dissection of this program, with special reference to what takes place during mouse embryonic stem cell differentiation. We make an attempt to interpret what these findings might represent and discuss their potential relevance to embryonic development. In doing so, we revive an old concept of "autosomal Lyonization" to describe "facultative heterochromatinization" and irreversible silencing of individual replication domains on autosomes reminiscent of the stable silencing of the inactive X chromosome, which takes place at a stage equivalent to the postimplantation epiblast in mice.

We have investigated the role of the histone methyltransferase G9a in the establishment of silent nuclear compartments. Following conditional knockout of the G9a methyltransferase in mouse ESCs, 167 genes were significantly up-regulated, and no genes were strongly down-regulated. A partially overlapping set of 119 genes were up-regulated after differentiation of G9a-depleted cells to neural precursors. Promoters of these G9a-repressed genes were AT rich and H3K9me2 enriched but H3K4me3 depleted and were not highly DNA methylated. Representative genes were found to be close to the nuclear periphery, which was significantly enriched for G9a-dependent H3K9me2. Strikingly, although 73% of total genes were early replicating, more than 71% of G9a-repressed genes were late replicating, and a strong correlation was found between H3K9me2 and late replication. However, G9a loss did not significantly affect subnuclear position or replication timing of any non-pericentric regions of the genome, nor did it affect programmed changes in replication timing that accompany differentiation. We conclude that G9a is a gatekeeper for a specific set of genes localized within the late replicating nuclear periphery.

Replication timing is frequently discussed superficially in terms of its relationship to transcriptional activity via chromatin structure. However, so little is known about what regulates where and when replication initiates that it has been impossible to identify mechanistic and causal relationships. Moreover, much of our knowledge base has been anecdotal, derived from analyses of a few genes in unrelated cell lines. Recent studies have revisited long-standing hypotheses using genome-wide approaches. In particular, the foundation of this field was recently shored up with incontrovertible evidence that cellular differentiation is accompanied by coordinated changes in replication timing and transcription. These changes accompany subnuclear repositioning, and take place at the level of megabase-sized domains that transcend localized changes in chromatin structure or transcription. Inferring from these results, we propose that there exists a key transition during the middle of S-phase and that changes in replication timing traversing this period are associated with subnuclear repositioning and changes in the activity of certain classes of promoters.

Although early replication has long been associated with accessible chromatin, replication timing is not included in most discussions of epigenetic marks. This is partly due to a lack of understanding of the mechanisms behind this association but the issue has also been confounded by studies concluding that there are very few changes in replication timing during development. Recently, the first genome-wide study of replication timing during the course of differentiation revealed extensive changes that were strongly associated with changes in transcriptional activity and subnuclear organization. Domains of temporally coordinate replication delineate discrete units of chromosome structure and function that are characteristic of particular differentiation states. Hence, although we are still a long way from understanding the functional significance of replication timing, it is clear that replication timing is a distinct epigenetic signature of cell differentiation state.

Background: Eukaryotic DNA replication is regulated at the level of large chromosomal domains (0.5-5 megabases in mammals) within which replicons are activated relatively synchronously. These domains replicate in a specific temporal order during S-phase and our genome-wide analyses of replication timing have demonstrated that this temporal order of domain replication is a stable property of specific cell types., Results: We have developed ReplicationDomain http://www.replicationdomain.org as a web-based database for analysis of genome-wide replication timing maps (replication profiles) from various cell lines and species. This database also provides comparative information of transcriptional expression and is configured to display any genome-wide property (for instance, ChIP-Chip or ChIP-Seq data) via an interactive web interface. Our published microarray data sets are publicly available. Users may graphically display these data sets for a selected genomic region and download the data displayed as text files, or alternatively, download complete genome-wide data sets. Furthermore, we have implemented a user registration system that allows registered users to upload their own data sets. Upon uploading, registered users may choose to: (1) view their data sets privately without sharing; (2) share with other registered users; or (3) make their published or "in press" data sets publicly available, which can fulfill journal and funding agencies' requirements for data sharing., Conclusion: ReplicationDomain is a novel and powerful tool to facilitate the comparative visualization of replication timing in various cell types as well as other genome-wide chromatin features and is considerably faster and more convenient than existing browsers when viewing multi-megabase segments of chromosomes. Furthermore, the data upload function with the option of private viewing or sharing of data sets between registered users should be a valuable resource for the scientific community.

DNA replication in mammals is regulated via the coordinate firing of clusters of replicons that duplicate megabase-sized chromosome segments at specific times during S-phase. Cytogenetic studies show that these "replicon clusters" coalesce as subchromosomal units that persist through multiple cell generations, but the molecular boundaries of such units have remained elusive. Moreover, the extent to which changes in replication timing occur during differentiation and their relationship to transcription changes has not been rigorously investigated. We have constructed high-resolution replication-timing profiles in mouse embryonic stem cells (mESCs) before and after differentiation to neural precursor cells. We demonstrate that chromosomes can be segmented into multimegabase domains of coordinate replication, which we call "replication domains," separated by transition regions whose replication kinetics are consistent with large originless segments. The molecular boundaries of replication domains are remarkably well conserved between distantly related ESC lines and induced pluripotent stem cells. Unexpectedly, ESC differentiation was accompanied by the consolidation of smaller differentially replicating domains into larger coordinately replicated units whose replication time was more aligned to isochore GC content and the density of LINE-1 transposable elements, but not gene density. Replication-timing changes were coordinated with transcription changes for weak promoters more than strong promoters, and were accompanied by rearrangements in subnuclear position. We conclude that replication profiles are cell-type specific, and changes in these profiles reveal chromosome segments that undergo large changes in organization during differentiation. Moreover, smaller replication domains and a higher density of timing transition regions that interrupt isochore replication timing define a novel characteristic of the pluripotent state.

The replication timing of some genes is developmentally regulated, but the significance of replication timing to cellular differentiation has been difficult to substantiate. Studies have largely been restricted to the comparison of a few genes in established cell lines derived from different tissues, and most of these genes do not change replication timing. Hence, it has not been possible to predict how many or what types of genes might be subject to such control. Here, we have evaluated the replication timing of 54 tissue-specific genes in mouse embryonic stem cells before and after differentiation to neural precursors. Strikingly, genes residing within isochores rich in GC and poor in long interspersed nuclear elements (LINEs) did not change their replication timing, whereas half of genes within isochores rich in AT and long interspersed nuclear elements displayed programmed changes in replication timing that accompanied changes in gene expression. Our results provide direct evidence that differentiation-induced autosomal replication-timing changes are a significant part of mammalian development, provide a means to predict genes subject to such regulation, and suggest that replication timing may be more related to the evolution of metazoan genomes than to gene function or expression pattern.

The Xenopus LIM homeodomain (LIM-HD) protein, Xlim-1, is expressed in the Spemann organizer and cooperates with its positive regulator, Ldb1, to activate organizer gene expression. While this activation is presumably mediated through Xlim-1/Ldb1 tetramer formation, the mechanisms regulating proper Xlim-1/Ldb1 stoichiometry remains largely unknown. We isolated the Xenopus ortholog (XRnf12) of the RING finger protein Rnf12/RLIM and explored its functional interactions with Xlim-1 and Ldb1. Although XRnf12 functions as a E3 ubiquitin ligase for Ldb1 and causes proteasome-dependent degradation of Ldb1, we found that co-expression of a high level of Xlim-1 suppresses Ldb1 degradation by XRnf12. This suppression requires both the LIM domains of Xlim-1 and the LIM interaction domain of Ldb1, suggesting that Ldb1, when bound to Xlim-1, escapes degradation by XRnf12. We further show that a high level of Ldb1 suppresses the organizer activity of Xlim-1/Ldb1, suggesting that excess Ldb1 molecules disturb Xlim-1/Ldb1 stoichiometry. Consistent with this, Ldb1 overexpression in the dorsal marginal zone suppresses expression of several organizer genes including postulated Xlim-1 targets, and importantly, this suppression is rescued by co-expression of XRnf12. These data suggest that XRnf12 confers proper Ldb1 protein levels and Xlim-1/Ldb1 stoichiometry for their functions in the organizer. Together with the similarity in the expression pattern of Ldb1 and XRnf12 throughout early embryogenesis, we propose Rnf12/RLIM as a specific regulator of Ldb1 to ensure its proper interactions with LIM-HD proteins and possibly other Ldb1-interacting proteins in the organizer as well as in other tissues.

The Spemann organizer plays a central role in neural induction, patterning of the neuroectoderm and mesoderm, and morphogenetic movements during early embryogenesis. By seeking genes whose expression is activated by the organizer-specific LIM homeobox gene Xlim-1 in Xenopus animal caps, we isolated the receptor tyrosine kinase Xror2. Xror2 is expressed initially in the dorsal marginal zone, then in the notochord and the neuroectoderm posterior to the midbrain-hindbrain boundary. mRNA injection experiments revealed that overexpression of Xror2 inhibits convergent extension of the dorsal mesoderm and neuroectoderm in whole embryos, as well as the elongation of animal caps treated with activin, whereas it does not appear to affect cell differentiation of neural tissue and notochord. Interestingly, mutant constructs in which the kinase domain was point-mutated or deleted (named Xror2-TM) also inhibited convergent extension, and did not counteract the wild-type, suggesting that the ectodomain of Xror2 per se has activities that may be modulated by the intracellular domain. In relation to Wnt signaling for planar cell polarity, we observed: (1) the Frizzled-like domain in the ectodomain is required for the activity of wild-type Xror2 and Xror2-TM; (2) co-expression of Xror2 with Xwnt11, Xfz7, or both, synergistically inhibits convergent extension in embryos; (3) inhibition of elongation by Xror2 in activin-treated animal caps is reversed by co-expression of a dominant negative form of Cdc42 that has been suggested to mediate the planar cell polarity pathway of Wnt; and (4) the ectodomain of Xror2 interacts with Xwnts in co-immunoprecipitation experiments. These results suggest that Xror2 cooperates with Wnts to regulate convergent extension of the axial mesoderm and neuroectoderm by modulating the planar cell polarity pathway of Wnt.

The Xenopus LIM homeodomain protein Xlim-1 is specifically expressed in the Spemann organizer region and assumed to play a role in the establishment of the body axis as a transcriptional activator. To further elucidate the mechanism underlying the regulation of its transcriptional activity, we focused on the region C-terminal to the homeodomain of Xlim-1 (CT239-403) and divided it into five regions, CCR1-5 (C-terminal conserved regions), based on similarity between Xlim-1 and its paralog, Xlim-5. The role of Xlim-1 CT239-403 in the Spemann organizer was analyzed by assaying the axis-forming ability of a series of CCR-mutated constructs in Xenopus embryos. We show that high doses of Xlim-1 constructs deleted of CCR1 or CCR2 initiate secondary axis formation in the absence of its coactivator Ldb1 (LIM-domain-binding protein 1), suggesting that CCR1 and CCR2 are involved in negative regulation of Xlim-1. In contrast, while Xlim-1 is capable of initiating secondary axis formation at low doses in the presence of Ldb1, deletion of CCR2 (aa 275-295) or substitution of five conserved tyrosines in CCR2 with alanines (CCR2-5YA) abolished the activity. In addition, UAS-GAL4 one-hybrid reporter assays in Xenopus showed that CCR2, but not CCR2-5YA, with its flanking regions (aa 261-315) functions as a transactivation domain when fused to the GAL4 DNA-binding domain. Finally, we show that none of the known transcriptional coactivators tested (CBP, SRC-1, and TIF2) interacts with the Xlim-1 transactivation domain (aa 261-315). Thus, Xlim-1 not only contains a unique tyrosine-rich activation domain but also contains a negative regulatory domain in CT239-403, suggesting a complex regulatory mechanism underlying the transcriptional activity of Xlim-1 in the organizer.