The genetic architecture of DNA replication timing in human pluripotent stem cells

DNA replication follows a strict spatiotemporal program that intersects with chromatin structure but has a poorly understood genetic basis. To systematically identify genetic regulators of replication timing, we exploited inter-individual variation in human pluripotent stem cells from 349 individuals. We show that the human genome’s replication program is broadly encoded in DNA and identify 1,617 cis-acting replication timing quantitative trait loci (rtQTLs) – sequence determinants of replication initiation. rtQTLs function individually, or in combinations of proximal and distal regulators, and are enriched at sites of histone H3 trimethylation of lysines 4, 9, and 36 together with histone hyperacetylation. H3 trimethylation marks are individually repressive yet synergistically associate with early replication. We identify pluripotency-related transcription factors and boundary elements as positive and negative regulators of replication timing, respectively. Taken together, human replication timing is controlled by a multi-layered mechanism with dozens of effectors working combinatorially and following principles analogous to transcription regulation.

Comprehensive analysis of DNA replication timing in genetic diseases and gene knockouts identifies MCM10 as a novel regulator of the replication program

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. MCM10 mutant cells demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in MCM10-mutant cells was predominantly comprised of replication initiation defects. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel modulator of DNA replication timing.

Asymmetrical deposition and modification of histone H3 variants are essential for zygote development

The pericentromeric heterochromatin of one-cell embryos forms a unique, ring-like structure around the nucleolar precursor body, which is absent in somatic cells. Here, we found that the histone H3 variants H3.1 and/or H3.2 (H3.1/H3.2) were localized asymmetrically between the male and female perinucleolar regions of the one-cell embryos; moreover, asymmetrical histone localization influenced DNA replication timing. The nuclear deposition of H3.1/3.2 in one-cell embryos was low relative to other preimplantation stages because of reduced H3.1/3.2 mRNA expression and incorporation efficiency. The forced incorporation of H3.1/3.2 into the pronuclei of one-cell embryos triggered a delay in DNA replication, leading to developmental failure. Methylation of lysine residue 27 (H3K27me3) of the deposited H3.1/3.2 in the paternal perinucleolar region caused this delay in DNA replication. These results suggest that reduced H3.1/3.2 in the paternal perinucleolar region is essential for controlled DNA replication and preimplantation development. The nuclear deposition of H3.1/3.2 is presumably maintained at a low level to avoid the detrimental effect of K27me3 methylation on DNA replication in the paternal perinucleolar region.

Low Replicative Stress Triggers Cell-Type Specific Inheritable Advanced Replication Timing

DNA replication timing (RT), reflecting the temporal order of origin activation, is known as a robust and conserved cell-type specific process. Upon low replication stress, the slowing of replication forks induces well-documented RT delays associated to genetic instability, but it can also generate RT advances that are still uncharacterized. In order to characterize these advanced initiation events, we monitored the whole genome RT from six independent human cell lines treated with low doses of aphidicolin. We report that RT advances are cell-type-specific and involve large heterochromatin domains. Importantly, we found that some major late to early RT advances can be inherited by the unstressed next-cellular generation, which is a unique process that correlates with enhanced chromatin accessibility, as well as modified replication origin landscape and gene expression in daughter cells. Collectively, this work highlights how low replication stress may impact cellular identity by RT advances events at a subset of chromosomal domains.

TIGER: inferring DNA replication timing from whole-genome sequence data

Motivation Genomic DNA replicates according to a reproducible spatiotemporal program, with some loci replicating early in S phase while others replicate late. Despite being a central cellular process, DNA replication timing studies have been limited in scale due to technical challenges. Results We present TIGER (Timing Inferred from Genome Replication), a computational approach for extracting DNA replication timing information from whole genome sequence data obtained from proliferating cell samples. The presence of replicating cells in a biological specimen leads to non-uniform representation of genomic DNA that depends on the timing of replication of different genomic loci. Replication dynamics can hence be observed in genome sequence data by analyzing DNA copy number along chromosomes while accounting for other sources of sequence coverage variation. TIGER is applicable to any species with a contiguous genome assembly and rivals the quality of experimental measurements of DNA replication timing. It provides a straightforward approach for measuring replication timing and can readily be applied at scale. Availability and Implementation TIGER is available at https://github.com/TheKorenLab/TIGER. Supplementary information Supplementary data are available at Bioinformatics online