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.

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. Microscopy studies have suggested that chromosomal DNA is composed of multiple, megabase-sized segments, each replicated at different times during S-phase of the cell cycle. However, a molecular definition of these coordinately replicated sequences and the stability of the boundaries between them has not been established. We constructed genome-wide replication-timing maps in mouse embryonic stem cells, identifying multimegabase coordinately replicated chromosome segments—“replication domains”—separated by remarkably distinct temporal boundaries. These domain boundaries were shared between several unrelated embryonic stem cell lines, including somatic cells reprogrammed to pluripotency (so-called induced pluripotent stem cells). However, upon differentiation to neural precursor cells, domains encompassing approximately 20% of the genome changed their replication timing, temporally consolidating into fewer, larger replication domains that were conserved between different neural precursor cell lines. Domains that changed replication timing showed a unique sequence composition, a strongly biased directionality for changes in resident gene expression, and altered radial positioning within the three-dimensional space in the cell nucleus, suggesting that changes in replication timing are related to the reorganization of higher-order chromosome structure and function during differentiation. Moreover, the property of smaller discordantly replicating domains may define a novel characteristic of pluripotency.

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.

Centromeres are specified by sequence-independent epigenetic mechanisms in most organisms. Rarely, centromere repositioning results in neocentromere formation at ectopic sites. However, the mechanisms governing how and where neocentromeres form are unknown. Here, we established a chromosome-engineering system in chicken DT40 cells that allowed us to efficiently isolate neocentromere-containing chromosomes. Neocentromeres appear to be structurally and functionally equivalent to native centromeres. Chromatin immunoprecipitation sequencing (ChIP-seq) analysis with 18 neocentromeres revealed that the centromere-specific histone H3 variant CENP-A occupies an ∼40 kb region at each neocentromere, which has no preference for specific DNA sequence motifs. Furthermore, we found that neocentromeres were not associated with histone modifications H3K9me3, H3K4me2, and H3K36me3 or with early replication timing. Importantly, low but significant levels of CENP-A are detected around endogenous centromeres, which are capable of seeding neocentromere assembly if the centromere core is removed. In summary, our experimental system provides valuable insights for understanding how neocentromeres form.

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.