Global organization of replication time zones of the mouse genome

Organization of DNA replication

The discovery of the DNA double helix structure half a century ago immediately suggested a mechanism for its duplication by semi-conservative copying of the nucleotide sequence into two DNA daughter strands. Shortly after, a second fundamental step toward the elucidation of the mechanism of DNA replication was taken with the isolation of the first enzyme able to polymerize DNA from a template. In the subsequent years, the basic mechanism of DNA replication and its enzymatic machinery components were elucidated, mostly through genetic approaches and in vitro biochemistry. Most recently, the spatial and temporal organization of the DNA replication process in vivo within the context of chromatin and inside the intact cell are finally beginning to be elucidated. On the one hand, recent advances in genome-wide high throughput techniques are providing a new wave of information on the progression of genome replication at high spatial resolution. On the other hand, novel super-resolution microscopy techniques are just starting to give us the first glimpses of how DNA replication is organized within the context of single intact cells with high spatial resolution. The integration of these data with time lapse microscopy analysis will give us the ability to film and dissect the replication of the genome in situ and in real time.

Impact of replication timing on non-CpG and CpG substitution rates in mammalian genomes

Neutral nucleotide substitutions occur at varying rates along genomes, and it remains a major issue to unravel the mechanisms that cause these variations and to analyze their evolutionary consequences. Here, we study the role of replication in the neutral substitution pattern. We obtained a high-resolution replication timing profile of the whole human genome by massively parallel sequencing of nascent BrdU-labeled replicating DNA. These data were compared to the neutral substitution rates along the human genome, obtained by aligning human and chimpanzee genomes using macaque and orangutan as outgroups. All substitution rates increase monotonously with replication timing even after controlling for local or regional nucleotide composition, crossover rate, distance to telomeres, and chromatin compaction. The increase in non-CpG substitution rates might result from several mechanisms including the increase in mutation-prone activities or the decrease in efficiency of DNA repair during the S phase. In contrast, the rate of C --> T transitions in CpG dinucleotides increases in later-replicating regions due to increasing DNA methylation level that reflects a negative correlation between timing and gene expression. Similar results are observed in the mouse, which indicates that replication timing is a main factor affecting nucleotide substitution dynamics at non-CpG sites and constitutes a major neutral process driving mammalian genome evolution.

Replication factory activation can be decoupled from the replication timing program by modulating Cdk levels

Cdk activity can differentially regulate the number of active replication factories, replication rates, and the rate of progression through the timing program during S phase.

In the metazoan replication timing program, clusters of replication origins located in different subchromosomal domains fire at different times during S phase. We have used Xenopus laevis egg extracts to drive an accelerated replication timing program in mammalian nuclei. Although replicative stress caused checkpoint-induced slowing of the timing program, inhibition of checkpoint kinases in an unperturbed S phase did not accelerate it. Lowering cyclin-dependent kinase (Cdk) activity slowed both replication rate and progression through the timing program, whereas raising Cdk activity increased them. Surprisingly, modest alteration of Cdk activity changed the amount of DNA synthesized during different stages of the timing program. This was associated with a change in the number of active replication factories, whereas the distribution of origins within active factories remained relatively normal. The ability of Cdks to differentially effect replication initiation, factory activation, and progression through the timing program provides new insights into the way that chromosomal DNA replication is organized during S phase.

Decreased replication origin activity in temporal transition regions

Experimental attempts to activate replication origins within the temporal transition region in the IgH locus in mouse embryonic stem cells were not successful, and thus, why and how they become activated in B cells remains unclear.

In the mammalian genome, early- and late-replicating domains are often separated by temporal transition regions (TTRs) with novel properties and unknown functions. We identified a TTR in the mouse immunoglobulin heavy chain (Igh) locus, which contains replication origins that are silent in embryonic stem cells but activated during B cell development. To investigate which factors contribute to origin activation during B cell development, we systematically modified the genetic and epigenetic status of the endogenous Igh TTR and used a single-molecule approach to analyze DNA replication. Introduction of a transcription unit into the Igh TTR, activation of gene transcription, and enhancement of local histone modifications characteristic of active chromatin did not lead to origin activation. Moreover, very few replication initiation events were observed when two ectopic replication origin sequences were inserted into the TTR. These findings indicate that the Igh TTR represents a repressive compartment that inhibits replication initiation, thus maintaining the boundaries between early and late replication domains.

Timing of replication is a determinant of neutral substitution rates but does not explain slow Y chromosome evolution in rodents

Mutation rates, assayed as substitution rates of putatively neutral sites, are highly variable around mammalian genomes: There is heterogeneity between genes, between autosomes, and between X, Y, and autosomes. The differences between X, Y, and autosomes are typically assumed to reflect the greater number of cell divisions in the male germ-line. Such an effect can neither account for within-autosome differences nor does it predict the differences between X, Y, and autosome observed in rodents. It has recently been proposed that in primates, the time during S-phase when a gene is replicated is an important determinant of neutral rates of evolution. Here we ask 1) whether we can replicate this result in rodents, 2) whether different autosomes replicate on average at different times, and 3) whether this might explain differences in their substitution rates. Finally we ask 4) whether X, Y, and autosome replicate at different times and 5) whether any difference might explain why the number of replication events alone cannot explain their substitution rates. We find that, as in primates, autosomal intronic rates of evolution increase significantly during S-phase. Different autosomes do have different average replication times, and together with rearrangement, this is a significant predictor of between-autosome differences in substitution rate. Although we find that autosomal, X-, and Y-linked genes replicate at different times, it is paradoxical that the Y-linked genes replicate latest, and replicate more often, but are not especially fast evolving. These results support the hypothesis that replication timing is an important source of substitution rate heterogeneity.

Sequencing newly replicated DNA reveals widespread plasticity in human replication timing

Faithful transmission of genetic material to daughter cells involves a characteristic temporal order of DNA replication, which may play a significant role in the inheritance of epigenetic states. We developed a genome-scale approach--Repli Seq--to map temporally ordered replicating DNA using massively parallel sequencing and applied it to study regional variation in human DNA replication time across multiple human cell types. The method requires as few as 8,000 cytometry-fractionated cells for a single analysis, and provides high-resolution DNA replication patterns with respect to both cell-cycle time and genomic position. We find that different cell types exhibit characteristic replication signatures that reveal striking plasticity in regional replication time patterns covering at least 50% of the human genome. We also identified autosomal regions with marked biphasic replication timing that include known regions of monoallelic expression as well as many previously uncharacterized domains. Comparison with high-resolution genome-wide profiles of DNaseI sensitivity revealed that DNA replication typically initiates within foci of accessible chromatin comprising clustered DNaseI hypersensitive sites, and that replication time is better correlated with chromatin accessibility than with gene expression. The data collectively provide a unique, genome-wide picture of the epigenetic compartmentalization of the human genome and suggest that cell-lineage specification involves extensive reprogramming of replication timing patterns.

Genome-wide dynamics of replication timing revealed by in vitro models of mouse embryogenesis

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.

Transcriptional co-activator LEDGF interacts with Cdc7-activator of S-phase kinase (ASK) and stimulates its enzymatic activity

Lens epithelium-derived growth factor (LEDGF) is an important co-factor of human immunodeficiency virus DNA integration; however, its cellular functions are poorly characterized. We now report identification of the Cdc7-activator of S-phase kinase (ASK) heterodimer as a novel interactor of LEDGF. Both kinase subunits co-immunoprecipitated with endogenous LEDGF from human cell extracts. Truncation analyses identified the integrase-binding domain of LEDGF as essential and minimally sufficient for the interaction with Cdc7-ASK. Reciprocally, the interaction required autophosphorylation of the kinase and the presence of 50 C-terminal residues of ASK. The kinase phosphorylated LEDGF in vitro, with Ser-206 being the major target, and LEDGF phosphorylated at this residue could be detected during S phase of the cell cycle. LEDGF potently stimulated the enzymatic activity of Cdc7-ASK, increasing phosphorylation of MCM2 in vitro by more than 10-fold. This enzymatic stimulation as well as phosphorylation of LEDGF depended on the protein-protein interaction. Intriguingly, removing the C-terminal region of ASK, involved in the interaction with LEDGF, resulted in a hyperactive kinase. Our results indicate that the interaction with LEDGF relieves autoinhibition of Cdc7-ASK kinase, imposed by the C terminus of ASK.

Increased rate of human mutations where DNA and RNA polymerases collide

Gene density and orientation of genes in eukaryotes seem to be correlated with the replication origin and the mutation rate is greater in late replicating regions; however, the reason for these patterns is unknown. Here, we investigate predicted replication origins in the human genome and find that levels of polymorphism as well as divergence from the chimpanzee genome are greater in genes transcribed on the lagging strand than those on the leading strand. This might be caused by interference between RNA and DNA polymerases, and avoidance of collisions between these enzymes might be an evolutionary force shaping gene orientation and density surrounding replication start sites. Physical constraints might have a larger influence on genome evolution than previously thought.

Temporal regulation of DNA replication in mammalian cells

Eukaryotic cells follow a temporal program to duplicate their genomes. Chromosomes are divided into domains with a specific DNA replication timing (RT), not dictated by DNA sequence alone, which is conserved from one cell cycle to the next. Timing of replication correlates with gene density, transcriptional activity, chromatin structure and nuclear position, making it an intriguing epigenetic mark. The differentiation from embryonic stem cells to specialized cell types is accompanied by global changes in the RT program. This review covers our current understanding of the mechanisms that determine RT in mammalian cells, its possible biological significance and how unscheduled alterations of the RT program may predispose to human disease.

Predictable dynamic program of timing of DNA replication in human cells

The organization of mammalian DNA replication is poorly understood. We have produced high-resolution dynamic maps of the timing of replication in human erythroid, mesenchymal, and embryonic stem (ES) cells using TimEX, a method that relies on gaussian convolution of massive, highly redundant determinations of DNA copy-number variations during S phase to produce replication timing profiles. We first obtained timing maps of 3% of the genome using high-density oligonucleotide tiling arrays and then extended the TimEX method genome-wide using massively parallel sequencing. We show that in untransformed human cells, timing of replication is highly regulated and highly synchronous, and that many genomic segments are replicated in temporal transition regions devoid of initiation, where replication forks progress unidirectionally from origins that can be hundreds of kilobases away. Absence of initiation in one transition region is shown at the molecular level by single molecule analysis of replicated DNA (SMARD). Comparison of ES and erythroid cells replication patterns revealed that these cells replicate about 20% of their genome in different quarters of S phase. Importantly, we detected a strong inverse relationship between timing of replication and distance to the closest expressed gene. This relationship can be used to predict tissue-specific timing of replication profiles from expression data and genomic annotations. We also provide evidence that early origins of replication are preferentially located near highly expressed genes, that mid-firing origins are located near moderately expressed genes, and that late-firing origins are located far from genes.

Open chromatin encoded in DNA sequence is the signature of 'master' replication origins in human cells

For years, progress in elucidating the mechanisms underlying replication initiation and its coupling to transcriptional activities and to local chromatin structure has been hampered by the small number (approximately 30) of well-established origins in the human genome and more generally in mammalian genomes. Recent in silico studies of compositional strand asymmetries revealed a high level of organization of human genes around 1000 putative replication origins. Here, by comparing with recently experimentally identified replication origins, we provide further support that these putative origins are active in vivo. We show that regions approximately 300-kb wide surrounding most of these putative replication origins that replicate early in the S phase are hypersensitive to DNase I cleavage, hypomethylated and present a significant enrichment in genomic energy barriers that impair nucleosome formation (nucleosome-free regions). This suggests that these putative replication origins are specified by an open chromatin structure favored by the DNA sequence. We discuss how this distinctive attribute makes these origins, further qualified as 'master' replication origins, priviledged loci for future research to decipher the human spatio-temporal replication program. Finally, we argue that these 'master' origins are likely to play a key role in genome dynamics during evolution and in pathological situations.

Shifts in replication timing actively affect histone acetylation during nucleosome reassembly

The entire genome is replicated in a programmed manner, with specific regions undergoing DNA synthesis at different times in S phase. Active genes generally replicate in early S phase, while repressed genes replicate late, and for some loci this process is developmentally regulated. Using a nuclear microinjection system, we demonstrate that DNA sequences originally packaged into nucleosomes containing deacetylated histones during late S become reassembled with acetylated histones after undergoing replication in early S. Conversely, a change from early to late replication timing is accompanied by repackaging into nucleosomes containing deacetylated histones. This is carried out by differential cell-cycle-controlled acetylation and deacetylation of histones H3 and H4. These studies provide strong evidence that switches in replication timing may play a role in the regulation of nucleosome structure during development.

Transcription initiation activity sets replication origin efficiency in mammalian cells

Genomic mapping of DNA replication origins (ORIs) in mammals provides a powerful means for understanding the regulatory complexity of our genome. Here we combine a genome-wide approach to identify preferential sites of DNA replication initiation at 0.4% of the mouse genome with detailed molecular analysis at distinct classes of ORIs according to their location relative to the genes. Our study reveals that 85% of the replication initiation sites in mouse embryonic stem (ES) cells are associated with transcriptional units. Nearly half of the identified ORIs map at promoter regions and, interestingly, ORI density strongly correlates with promoter density, reflecting the coordinated organisation of replication and transcription in the mouse genome. Detailed analysis of ORI activity showed that CpG island promoter-ORIs are the most efficient ORIs in ES cells and both ORI specification and firing efficiency are maintained across cell types. Remarkably, the distribution of replication initiation sites at promoter-ORIs exactly parallels that of transcription start sites (TSS), suggesting a co-evolution of the regulatory regions driving replication and transcription. Moreover, we found that promoter-ORIs are significantly enriched in CAGE tags derived from early embryos relative to all promoters. This association implies that transcription initiation early in development sets the probability of ORI activation, unveiling a new hallmark in ORI efficiency regulation in mammalian cells.

The duplication of the genetic information of a cell starts from specific sites on the chromosomes called DNA replication origins. Their number varies from a few hundred in yeast cells to several thousands in human cells, distributed along the genome at comparable distances in both systems. An important question in the field is to understand how origins of replication are specified and regulated in the mammalian genome, as neither their location nor their activity can be directly inferred from the DNA sequence. Previous studies at individual origins and, more recently, at large scale across 1% of the human genome, have revealed that most origins overlap with transcriptional regulatory elements, and specifically with gene promoters. To gain insight into the nature of the relationship between active transcription and origin specification we have combined a genomic mapping of origins at 0.4% of the mouse genome with detailed studies of activation efficiency. The data identify two types of origins with distinct regulatory properties: highly efficient origins map at CpG island-promoters and low efficient origins locate elsewhere in association with transcriptional units. We also find a remarkable parallel organisation of the replication initiation sites and transcription start sites at efficient promoter-origins that suggests a prominent role of transcription initiation in setting the efficiency of replication origin activation.

Replication timing and epigenetic reprogramming of gene expression: a two-way relationship?

An overall link between the potential for gene transcription and the timing of replication in S phase is now well established in metazoans. Here we discuss emerging evidence that highlights the possibility that replication timing is causally linked with epigenetic reprogramming. In particular, we bring together conclusions from a range of studies to propose a model in which reprogramming factors determine the timing of replication and the implementation of reprogramming events requires passage through S phase. These considerations have implications for our understanding of development, evolution and diseases such as cancer.

NA-Seq: a discovery tool for the analysis of chromatin structure and dynamics during differentiation

It is well established that epigenetic modulation of genome accessibility in chromatin occurs during biological processes. Here we describe a method based on restriction enzymes and next-generation sequencing for identifying accessible DNA elements using a small amount of starting material, and use it to examine myeloid differentiation of primary human CD34+ cells. The accessibility of several classes of cis-regulatory elements was a predictive marker of in vivo DNA binding by transcription factors, and was associated with distinct patterns of histone posttranslational modifications. We also mapped large chromosomal domains with differential accessibility in progenitors and maturing cells. Accessibility became restricted during differentiation, correlating with a decreased number of expressed genes and loss of regulatory potential. Our data suggest that a permissive chromatin structure in multipotent cells is progressively and selectively closed during differentiation, and illustrate the use of our method for the identification of functional cis-regulatory elements.

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.

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

The relationship between DNA replication and human genome organization

Assessment of the impact of DNA replication on genome architecture in Eukaryotes has long been hampered by the scarcity of experimental data. Recent work, relying on computational predictions of origins of replication, suggested that replication might be a major determinant of gene organization in human (Huvet et al. 2007. Human gene organization driven by the coordination of replication and transcription. Genome Res. 17:1278-1285). Here, we address this question by analyzing the first large-scale data set of experimentally determined origins of replication in human: 283 origins identified in HeLa cells, in 1% of the genome covered by ENCODE regions (Cadoret et al. 2008. Genome-wide studies highlight indirect links between human replication origins and gene regulation. Proc Natl Acad Sci USA. 105:15837-15842). We show that origins of replication are not randomly distributed as they display significant overlap with promoter regions and CpG islands. The hypothesis of a selective pressure to avoid frontal collisions between replication and transcription polymerases is not supported by experimental data as we find no evidence for gene orientation bias in the proximity of origins of replication. The lack of a significant orientation bias remains manifest even when considering only genes expressed at a high rate, or in a wide number of tissues, and is not affected by the regional replication timing. Gene expression breadth does not appear to be correlated with the distance from the origins of replication. We conclude that the impact of DNA replication on human genome organization is considerably weaker than previously proposed.

ReplicationDomain: a visualization tool and comparative database for genome-wide replication timing data

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.