Temporal control of late replication and coordination of origin firing by self-stabilizing Rif1-PP1 hubs in Drosophila.

In the metazoan S phase, coordinated firing of clusters of origins replicates different parts of the genome in a temporal program. Despite advances, neither the mechanism controlling timing nor that coordinating firing of multiple origins is fully understood. Rif1, an evolutionarily conserved inhibitor of DNA replication, recruits protein phosphatase 1 (PP1) and counteracts firing of origins by S-phase kinases. During the midblastula transition (MBT) in Drosophila embryos, Rif1 forms subnuclear hubs at each of the large blocks of satellite sequences and delays their replication. Each Rif1 hub disperses abruptly just prior to the replication of the associated satellite sequences. Here, we show that the level of activity of the S-phase kinase, DDK, accelerated this dispersal program, and that the level of Rif1-recruited PP1 retarded it. Further, Rif1-recruited PP1 supported chromatin association of nearby Rif1. This influence of nearby Rif1 can create a "community effect" counteracting kinase-induced dissociation such that an entire hub of Rif1 undergoes switch-like dispersal at characteristic times that shift in response to the balance of Rif1-PP1 and DDK activities. We propose a model in which the spatiotemporal program of late replication in the MBT embryo is controlled by self-stabilizing Rif1-PP1 hubs, whose abrupt dispersal synchronizes firing of associated late origins.

Dbf4 recruitment by forkhead transcription factors defines an upstream rate-limiting step in determining origin firing timing.

Initiation of eukaryotic chromosome replication follows a spatiotemporal program. The current model suggests that replication origins compete for a limited pool of initiation factors. However, it remains to be answered how these limiting factors are preferentially recruited to early origins. Here, we report that Dbf4 is enriched at early origins through its interaction with forkhead transcription factors Fkh1 and Fkh2. This interaction is mediated by the Dbf4 C terminus and was successfully reconstituted in vitro. An interaction-defective mutant, dbf4ΔC, phenocopies fkh alleles in terms of origin firing. Remarkably, genome-wide replication profiles reveal that the direct fusion of the DNA-binding domain (DBD) of Fkh1 to Dbf4 restores the Fkh-dependent origin firing but interferes specifically with the pericentromeric origin activation. Furthermore, Dbf4 interacts directly with Sld3 and promotes the recruitment of downstream limiting factors. These data suggest that Fkh1 targets Dbf4 to a subset of noncentromeric origins to promote early replication in a manner that is reminiscent of the recruitment of Dbf4 to pericentromeric origins by Ctf19.

DNA replication timing: Biochemical mechanisms and biological significance.

The regulation of DNA replication is a fascinating biological problem both from a mechanistic angle-How is replication timing regulated?-and from an evolutionary one-Why is replication timing regulated? Recent work has provided significant insight into the first question. Detailed biochemical understanding of the mechanism and regulation of replication initiation has made possible robust hypotheses for how replication timing is regulated. Moreover, technical progress, including high-throughput, single-molecule mapping of replication initiation and single-cell assays of replication timing, has allowed for direct testing of these hypotheses in mammalian cells. This work has consolidated the conclusion that differential replication timing is a consequence of the varying probability of replication origin initiation. The second question is more difficult to directly address experimentally. Nonetheless, plausible hypotheses can be made and one-that replication timing contributes to the regulation of chromatin structure-has received new experimental support.

Positive and Negative Regulation of DNA Replication Initiation.

Genomic DNA is replicated every cell cycle by the programmed activation of replication origins at specific times and chromosomal locations. The factors that define the locations of replication origins and their typical activation times in eukaryotic cells are poorly understood. Previous studies highlighted the role of activating factors and epigenetic modifications in regulating replication initiation. Here, we review the role that repressive pathways - and their alleviation - play in establishing the genomic landscape of replication initiation. Several factors mediate this repression, in particular, factors associated with inactive chromatin. Repression can support organized, yet stochastic, replication initiation, and its absence could explain instances of rapid and random replication or re-replication.

Mutation Rate Variability across Human Y-Chromosome Haplogroups.

A common assumption in dating patrilineal events using Y-chromosome sequencing data is that the Y-chromosome mutation rate is invariant across haplogroups. Previous studies revealed interhaplogroup heterogeneity in phylogenetic branch length. Whether this heterogeneity is caused by interhaplogroup mutation rate variation or nongenetic confounders remains unknown. Here, we analyzed whole-genome sequences from cultured cells derived from >1,700 males. We confirmed the presence of branch length heterogeneity. We demonstrate that sex-chromosome mutations that appear within cell lines, which likely occurred somatically or in vitro (and are thus not influenced by nongenetic confounders) are informative for germline mutational processes. Using within-cell-line mutations, we computed a relative Y-chromosome somatic mutation rate, and uncovered substantial variation (up to 83.3%) in this proxy for germline mutation rate among haplogroups. This rate positively correlates with phylogenetic branch length, indicating that interhaplogroup mutation rate variation is a likely cause of branch length heterogeneity.

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.

Conserved forkhead dimerization motif controls DNA replication timing and spatial organization of chromosomes in S. cerevisiae.

Forkhead Box (Fox) proteins share the Forkhead domain, a winged-helix DNA binding module, which is conserved among eukaryotes from yeast to humans. These sequence-specific DNA binding proteins have been primarily characterized as transcription factors regulating diverse cellular processes from cell cycle control to developmental fate, deregulation of which contributes to developmental defects, cancer, and aging. We recently identified Saccharomyces cerevisiae Forkhead 1 (Fkh1) and Forkhead 2 (Fkh2) as required for the clustering of a subset of replication origins in G1 phase and for the early initiation of these origins in the ensuing S phase, suggesting a mechanistic role linking the spatial organization of the origins and their activity. Here, we show that Fkh1 and Fkh2 share a unique structural feature of human FoxP proteins that enables FoxP2 and FoxP3 to form domain-swapped dimers capable of bridging two DNA molecules in vitro. Accordingly, Fkh1 self-associates in vitro and in vivo in a manner dependent on the conserved domain-swapping region, strongly suggestive of homodimer formation. Fkh1- and Fkh2-domain-swap-minus (dsm) mutations are functional as transcription factors yet are defective in replication origin timing control. Fkh1-dsm binds replication origins in vivo but fails to cluster them, supporting the conclusion that Fkh1 and Fkh2 dimers perform a structural role in the spatial organization of chromosomal elements with functional importance.

DNA replication timing alterations identify common markers between distinct progeroid diseases.

Progeroid syndromes are rare genetic disorders that phenotypically resemble natural aging. Different causal mutations have been identified, but no molecular alterations have been identified that are in common to these diseases. DNA replication timing (RT) is a robust cell type-specific epigenetic feature highly conserved in the same cell types from different individuals but altered in disease. Here, we characterized DNA RT program alterations in Hutchinson-Gilford progeria syndrome (HGPS) and Rothmund-Thomson syndrome (RTS) patients compared with natural aging and cellular senescence. Our results identified a progeroid-specific RT signature that is common to cells from three HGPS and three RTS patients and distinguishes them from healthy individuals across a wide range of ages. Among the RT abnormalities, we identified the tumor protein p63 gene (TP63) as a gene marker for progeroid syndromes. By using the redifferentiation of four patient-derived induced pluripotent stem cells as a model for the onset of progeroid syndromes, we tracked the progression of RT abnormalities during development, revealing altered RT of the TP63 gene as an early event in disease progression of both HGPS and RTS. Moreover, the RT abnormalities in progeroid patients were associated with altered isoform expression of TP63 Our findings demonstrate the value of RT studies to identify biomarkers not detected by other methods, reveal abnormal TP63 RT as an early event in progeroid disease progression, and suggest TP63 gene regulation as a potential therapeutic target.

ChECing out Rif1 action in freely cycling cells.

In buddying yeast, like all eukaryotes examined so far, DNA replication is under temporal control, such that some origins fire early and some late during S phase. This replication timing program is established in G1 phase, where chromatin states are thought to prevent binding of key-limiting initiation factors at late-firing origins. Although many factors are involved in replication initiation, a new player, Rif1, has recently entered the scene, with a spate of papers revealing a global role for the protein in the control of replication initiation timing from yeasts to humans. Since budding yeast Rif1 was known to bind only to telomeric and silent mating loci regions, it remained controversial whether Rif1 acts directly at replication origins or instead influences origin activity indirectly. In this perspective, we discuss our recent finding that Rif1 binds directly to the replication origins that it controls. In this study, we also found that Rif1's regulatory activity at origins is best revealed by an assay (sort-seq) that measures replication in unperturbed, freely cycling cultures, as opposed to commonly used protocols in which cells are first blocked in the G1 phase of the cell cycle by mating pheromone, then released into a synchronous S phase. Finally, we discuss how the sequestration of Rif1 at telomeres, through an interaction with the arrays of Rap1 molecules bound there, plays an important role in limiting Rif1's action primarily to telomere-proximal replication origins.

Rif1-Dependent Regulation of Genome Replication in Mammals.

Eukaryotic genomes are replicated starting from multiple origins of replication. Their usage is tightly regulated, and not all the potential origins are activated during a single cell cycle. In addition, the ones that are activated are activated in a sequential order. Why don't origins of replication normally all fire together? Is this important? And if so, why? Would any order of firing do, or does the specific sequence matter? How is this process regulated? These questions concern all eukaryotes but have proven extremely hard to address because replication timing is a process intricately connected with multiple aspects of nuclear function.

DNA replication timing: Coordinating genome stability with genome regulation on the X chromosome and beyond.

Recent studies based on next-generation DNA sequencing have revealed that the female inactive X chromosome is replicated in a rapid, unorganized manner, and undergoes increased rates of mutation. These observations link the organization of DNA replication timing to gene regulation on one hand, and to the generation of mutations on the other hand. More generally, the exceptional biology of the inactive X chromosome highlights general principles of genome replication. Cells may control replication timing by a combination of intrinsic replication origin properties, local chromatin states and global levels of replication factors, leading to a functional separation between the activity of genes and their mutation.

supraHex: an R/Bioconductor package for tabular omics data analysis using a supra-hexagonal map.

Biologists are increasingly confronted with the challenge of quickly understanding genome-wide biological data, which usually involve a large number of genomic coordinates (e.g. genes) but a much smaller number of samples. To meet the need for data of this shape, we present an open-source package called 'supraHex' for training, analysing and visualising omics data. This package devises a supra-hexagonal map to self-organise the input data, offers scalable functionalities for post-analysing the map, and more importantly, allows for overlaying additional data for multilayer omics data comparisons. Via applying to DNA replication timing data of mouse embryogenesis, we demonstrate that supraHex is capable of simultaneously carrying out gene clustering and sample correlation, providing intuitive visualisation at each step of the analysis. By overlaying CpG and expression data onto the trained replication-timing map, we also show that supraHex is able to intuitively capture an inherent relationship between late replication, low CpG density promoters and low expression levels. As part of the Bioconductor project, supraHex makes accessible to a wide community in a simple way, what would otherwise be a complex framework for the ultrafast understanding of any tabular omics data, both scientifically and artistically. This package can run on Windows, Mac and Linux, and is freely available together with many tutorials on featuring real examples at http://supfam.org/supraHex.

Incomplete reprogramming of DNA replication timing in induced pluripotent stem cells.

Induced pluripotent stem cells (iPSCs) are the foundation of cell therapy. Differences in gene expression, DNA methylation, and chromatin conformation, which could affect differentiation capacity, have been identified between iPSCs and embryonic stem cells (ESCs). Less is known about whether DNA replication timing, a process linked to both genome regulation and genome stability, is efficiently reprogrammed to the embryonic state. To answer this, we compare genome-wide replication timing between ESCs, iPSCs, and cells reprogrammed by somatic cell nuclear transfer (NT-ESCs). While NT-ESCs replicate their DNA in a manner indistinguishable from ESCs, a subset of iPSCs exhibits delayed replication at heterochromatic regions containing genes downregulated in iPSCs with incompletely reprogrammed DNA methylation. DNA replication delays are not the result of gene expression or DNA methylation aberrations and persist after cells differentiate to neuronal precursors. Thus, DNA replication timing can be resistant to reprogramming and influence the quality of iPSCs.

DNA replication in early mammalian embryos is patterned, predisposing lamina-associated regions to fragility.

DNA replication in differentiated cells follows a defined program, but when and how it is established during mammalian development is not known. Here we show using single-cell sequencing, that both bovine and mouse cleavage stage embryos progress through S-phase in a defined pattern. Late replicating regions are associated with the nuclear lamina from the first cell cycle after fertilization, and contain few active origins, and few but long genes. Chromosome breaks, which form spontaneously in bovine embryos at sites concordant with human embryos, preferentially locate to late replicating regions. In mice, late replicating regions show enhanced fragility due to a sparsity of dormant origins that can be activated under conditions of replication stress. This pattern predisposes regions with long neuronal genes to fragility and genetic change prior to segregation of soma and germ line. Our studies show that the formation of early and late replicating regions is among the first layers of epigenetic regulation established on the mammalian genome after fertilization.

Impaired Replication Timing Promotes Tissue-Specific Expression of Common Fragile Sites.

Common fragile sites (CFSs) are particularly vulnerable regions of the genome that become visible as breaks, gaps, or constrictions on metaphase chromosomes when cells are under replicative stress. Impairment in DNA replication, late replication timing, enrichment of A/T nucleotides that tend to form secondary structures, the paucity of active or inducible replication origins, the generation of R-loops, and the collision between replication and transcription machineries on particularly long genes are some of the reported characteristics of CFSs that may contribute to their tissue-specific fragility. Here, we validated the induction of two CFSs previously found in the human fetal lung fibroblast line, Medical Research Council cell strain 5 (MRC-5), in another cell line derived from the same fetal tissue, Institute for Medical Research-90 cells (IMR-90). After induction of CFSs through aphidicolin, we confirmed the expression of the CFS 1p31.1 on chromosome 1 and CFS 3q13.3 on chromosome 3 in both fetal lines. Interestingly, these sites were found to not be fragile in lymphocytes, suggesting a role for epigenetic or transcriptional programs for this tissue specificity. Both these sites contained late-replicating genes NEGR1 (neuronal growth regulator 1) at 1p31.1 and LSAMP (limbic system-associated membrane protein) at 3q13.3, which are much longer, 0.880 and 1.4 Mb, respectively, than the average gene length. Given the established connection between long genes and CFS, we compiled information from the literature on all previously identified CFSs expressed in fibroblasts and lymphocytes in response to aphidicolin, including the size of the genes contained in each fragile region. Our comprehensive analysis confirmed that the genes found within CFSs are longer than the average human gene; interestingly, the two longest genes in the human genome are found within CFSs: Contactin Associated Protein 2 gene (CNTNAP2) in a lymphocytes' CFS, and Duchenne muscular dystrophy gene (DMD) in a CFS expressed in both lymphocytes and fibroblasts. This indicates that the presence of very long genes is a unifying feature of all CFSs. We also obtained replication profiles of the 1p31.1 and 3q13.3 sites under both perturbed and unperturbed conditions using a combination of fluorescent in situ hybridization (FISH) and immunofluorescence against bromodeoxyuridine (BrdU) on interphase nuclei. Our analysis of the replication dynamics of these CFSs showed that, compared to lymphocytes where these regions are non-fragile, fibroblasts display incomplete replication of the fragile alleles, even in the absence of exogenous replication stress. Our data point to the existence of intrinsic features, in addition to the presence of long genes, which affect DNA replication of the CFSs in fibroblasts, thus promoting chromosomal instability in a tissue-specific manner.

Next-Generation Sequencing Enables Spatiotemporal Resolution of Human Centromere Replication Timing.

Centromeres serve a critical function in preserving genome integrity across sequential cell divisions, by mediating symmetric chromosome segregation. The repetitive, heterochromatic nature of centromeres is thought to be inhibitory to DNA replication, but has also led to their underrepresentation in human reference genome assemblies. Consequently, centromeres have been excluded from genomic replication timing analyses, leaving their time of replication unresolved. However, the most recent human reference genome, hg38, included models of centromere sequences. To establish the experimental requirements for achieving replication timing profiles for centromeres, we sequenced G1- and S-phase cells from five human cell lines, and aligned the sequence reads to hg38. We were able to infer DNA replication timing profiles for the centromeres in each of the five cell lines, which showed that centromere replication occurs in mid-to-late S phase. Furthermore, we found that replication timing was more variable between cell lines in the centromere regions than expected, given the distribution of variation in replication timing genome-wide. These results suggest the potential of these, and future, sequence models to enable high-resolution studies of replication in centromeres and other heterochromatic regions.

Germline Structural Variations Are Preferential Sites of DNA Replication Timing Plasticity during Development.

The DNA replication timing program is modulated throughout development and is also one of the main factors influencing the distribution of mutation rates across the genome. However, the relationship between the mutagenic influence of replication timing and its developmental plasticity remains unexplored. Here, we studied the distribution of copy number variations (CNVs) and single nucleotide polymorphisms across the zebrafish genome in relation to changes in DNA replication timing during embryonic development in this model vertebrate species. We show that CNV sites exhibit strong replication timing plasticity during development, replicating significantly early during early development but significantly late during more advanced developmental stages. Reciprocally, genomic regions that changed their replication timing during development contained a higher proportion of CNVs than developmentally constant regions. Developmentally plastic CNV sites, in particular those that become delayed in their replication timing, were enriched for the clustered protocadherins, a set of genes important for neuronal development that have undergone extensive genetic and epigenetic diversification during zebrafish evolution. In contrast, single nucleotide polymorphism sites replicated consistently early throughout embryonic development, highlighting a unique aspect of the zebrafish genome. Our results uncover a hitherto unrecognized interface between development and evolution.

DNA Replication Timing Enters the Single-Cell Era.

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.

Telomere-binding factors in the regulation of DNA replication.

Recent studies have indicated new roles for telomere-binding factors in the regulation of DNA replication, not only at the telomeres but also at the arm regions of the chromosome. Among these factors, Rif1, a conserved protein originally identified in yeasts as a telomere regulator, plays a major role in the spatiotemporal regulation of DNA replication during S phase. Its ability to interact with phosphatases and to create specific higher-order chromatin structures is central to the mechanism by which Rif1 exerts this function. In this review, we discuss recent progress in elucidating the roles of Rif1 and other telomere-binding factors in the regulation of chromosome events occurring at locations other than telomeres.

Shepherding DNA ends: Rif1 protects telomeres and chromosome breaks.

Cells have evolved conserved mechanisms to protect DNA ends, such as those at the termini of linear chromosomes, or those at DNA double-strand breaks (DSBs). In eukaryotes, DNA ends at chromosomal termini are packaged into proteinaceous structures called telomeres. Telomeres protect chromosome ends from erosion, inadvertent activation of the cellular DNA damage response (DDR), and telomere fusion. In contrast, cells must respond to damage-induced DNA ends at DSBs by harnessing the DDR to restore chromosome integrity, avoiding genome instability and disease. Intriguingly, Rif1 (Rap1-interacting factor 1) has been implicated in telomere homeostasis as well as DSB repair. The protein was first identified in Saccharomyces cerevisiae as being part of the proteinaceous telosome. In mammals, RIF1 is not associated with intact telomeres, but was found at chromosome breaks, where RIF1 has emerged as a key mediator of pathway choice between the two evolutionary conserved DSB repair pathways of non-homologous end-joining (NHEJ) and homologous recombination (HR). While this functional dichotomy has long been a puzzle, recent findings link yeast Rif1 not only to telomeres, but also to DSB repair, and mechanistic parallels likely exist. In this review, we will provide an overview of the actions of Rif1 at DNA ends and explore how exclusion of end-processing factors might be the underlying principle allowing Rif1 to fulfill diverse biological roles at telomeres and chromosome breaks.