Ichiro Hiratani

Replication timing and transcriptional control: beyond cause and effect―part II

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.

G9a selectively represses a class of late-replicating genes at the nuclear periphery

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 as an epigenetic mark

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 sub-nuclear 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.

Space and time in the nucleus: developmental control of replication timing and chromosome architecture.

All eukaryotic cells replicate segments of their genomes in a defined temporal sequence. In multicellular organisms, at least half of the genome is subject to changes in this temporal sequence during development. We now know that this temporal sequence and its developmentally regulated changes are conserved across distantly related species, suggesting that it either represents or reflects something biologically important. However, both the mechanism and the significance of this program remain unknown. We recently demonstrated a remarkably strong genome-wide correlation between replication timing and chromatin interaction maps, stronger than any other chromosomal property analyzed to date, indicating that sequences localized close to one another replicate at similar times. This provides molecular confirmation of long-standing cytogenetic evidence for spatial compartmentalization of early- and late-replicating DNA and supports our earlier model that replication timing is reestablished in each G(1) phase, coincident with the anchorage of chromosomal segments at specific locations within the nucleus (timing decision point [TDP]). Here, we review the evidence linking the replication program to the three-dimensional architecture of chromatin in the nucleus and discuss what such a link might mean for the mechanism and significance of a developmentally regulated replication program.

Genome-scale analysis of replication timing: from bench to bioinformatics

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.

Replication Timing: A Fingerprint for Cell Identity and Pluripotency

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.

DNA Replication Timing Is Maintained Genome-Wide in Primary Human Myoblasts Independent of D4Z4 Contraction in FSH Muscular Dystrophy

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.

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

Here, we report the establishment of a single-cell DNA replication sequencing method, scRepli-seq, which is a simple 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 replication profiles, which were conserved from cell to cell. Haplotype-resolved scRepli-seq revealed similar replication timing 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, and we discovered that developmentally regulated domains are a source of such variability, suggesting a link between cell-to-cell heterogeneity and developmental plasticity. Together, our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into 3D genome organization.