DNA replication timing

Superresolution imaging reveals spatiotemporal propagation of human replication foci mediated by CTCF-organized chromatin structures

Mammalian DNA replication is initiated at numerous replication origins, which are clustered into thousands of replication domains (RDs) across the genome. However, it remains unclear whether the replication origins within each RD are activated stochastically or preferentially near certain chromatin features. To understand how DNA replication in single human cells is regulated at the sub-RD level, we directly visualized and quantitatively characterized the spatiotemporal organization, morphology, and in situ epigenetic signatures of individual replication foci (RFi) across S-phase at superresolution using stochastic optical reconstruction microscopy. Importantly, we revealed a hierarchical radial pattern of RFi propagation dynamics that reverses directionality from early to late S-phase and is diminished upon caffeine treatment or CTCF knockdown. Together with simulation and bioinformatic analyses, our findings point to a "CTCF-organized REplication Propagation" (CoREP) model, which suggests a nonrandom selection mechanism for replication activation at the sub-RD level during early S-phase, mediated by CTCF-organized chromatin structures. Collectively, these findings offer critical insights into the key involvement of local epigenetic environment in coordinating DNA replication across the genome and have broad implications for our conceptualization of the role of multiscale chromatin architecture in regulating diverse cell nuclear dynamics in space and time.

Active and Repressed Chromatin Domains Exhibit Distinct Nucleosome Segregation during DNA Replication

Chromatin domains and their associated structures must be faithfully inherited through cellular division to maintain cellular identity. However, accessing the localized strategies preserving chromatin domain inheritance, specifically the transfer of parental, pre-existing nucleosomes with their associated post-translational modifications (PTMs) during DNA replication, is challenging in living cells. We devised an inducible, proximity-dependent labeling system to irreversibly mark replication-dependent H3.1 and H3.2 histone-containing nucleosomes at desired loci in mouse embryonic stem cells so that their fate after DNA replication could be followed. Strikingly, repressed chromatin domains are preserved through local re-deposition of parental nucleosomes. In contrast, nucleosomes decorating active chromatin domains do not exhibit such preservation. Notably, altering cell fate leads to an adjustment of the positional inheritance of parental nucleosomes that reflects the corresponding changes in chromatin structure. These findings point to important mechanisms that contribute to parental nucleosome segregation to preserve cellular identity.

Antagonistic control of DDK binding to licensed replication origins by Mcm2 and Rad53.. [DOI: 10.1101/2020.05.04.077628] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Eukaryotic replication origins are licensed by the loading of the replicative DNA helicase, Mcm2-7, in inactive double hexameric form around DNA. Subsequent origin activation is under control of multiple protein kinases that either promote or inhibit origin activation, which is important for genome maintenance. Using the reconstituted budding yeast DNA replication system, we find that the flexible N-terminal tail of Mcm2 promotes the stable recruitment of Dbf4-dependent kinase (DDK) to Mcm2-7 double hexamers, which in turn promotes DDK phosphorylation of Mcm4 and -6 and subsequent origin activation. Conversely, we demonstrate that the checkpoint kinase, Rad53, inhibits DDK binding to Mcm2-7 double hexamers. Unexpectedly, this function is not dependent on Rad53 kinase activity, but requires Rad53 activation by trans-autophosphorylation, suggesting steric inhibition of DDK by activated Rad53. These findings identify critical determinants of the origin activation reaction and uncover a novel mechanism for checkpoint-dependent origin inhibition.

Adefovir dipivoxil induces DNA replication stress and augments ATR inhibitor-related cytotoxicity

Replication stress is a common feature of cancer cells. Ataxia telangiectasia-mutated (ATM) and Rad3-related (ATR) signalling, a DNA damage repair (DDR) pathway, is activated by regions of single-stranded DNA (ssDNA) that can arise during replication stress. ATR delays cell cycle progression and prevents DNA replication fork collapse, which prohibits cell death and promotes proliferation. Several ATR inhibitors have been developed in order to restrain this protective mechanism in tumours. It is known, however, that despite other effective anticancer chemotherapy treatments targeting DDR pathways, resistance occurs. This begets the need to identify combination treatments to overcome resistance and prevent tumour cell growth. We conducted a drug screen to identify potential synergistic combination treatments by screening an ATR inhibitor (VE822) together with compounds from a bioactive small molecule library. The screen identified adefovir dipivoxil, a reverse transcriptase inhibitor and nucleoside analogue, as a compound that has increased cytotoxicity in the presence of ATR, but not ATM or DNA-dependant protein kinase (DNA-PK) inhibition. Here we demonstrate that adefovir dipivoxil induces DNA replication stress, activates ATR signalling and stalls cells in S phase. This simultaneous induction of replication stress and inhibition of ATR signalling lead to a marked increase in pan-nuclear γH2AX-positive cells, ssDNA accumulation and cell death, indicative of replication catastrophe.

Chromatin domains in space and their functional implications

Genome organization displays functional compartmentalization. Many factors, including epigenetic modifications, transcription factors, chromatin remodelers, and RNAs, shape chromatin domains and the three-dimensional genome organization. Various types of chromatin domains with distinct epigenetic and spatial features exhibit different transcriptional activities. As part of the efforts to better understand plant functional genomics, over the past a few years, spatial distribution patterns of plant chromatin domains have been brought to light. In this review, we discuss chromatin domains associated with the nuclear periphery and the nucleolus, as well as chromatin domains staying in proximity and showing physical interactions. The functional implication of these domains is discussed, with a particular focus on the transcriptional regulation and replication timing. Finally, from a biophysical point of view, we discuss potential roles of liquid-liquid phase separation in plant nuclei in the genesis and maintenance of spatial chromatin domains.

High-resolution Repli-Seq defines the temporal choreography of initiation, elongation and termination of replication in mammalian cells

BACKGROUND

DNA replication in mammalian cells occurs in a defined temporal order during S phase, known as the replication timing (RT) programme. Replication timing is developmentally regulated and correlated with chromatin conformation and local transcriptional potential. Here, we present RT profiles of unprecedented temporal resolution in two human embryonic stem cell lines, human colon carcinoma line HCT116, and mouse embryonic stem cells and their neural progenitor derivatives.

RESULTS

Fine temporal windows revealed a remarkable degree of cell-to-cell conservation in RT, particularly at the very beginning and ends of S phase, and identified 5 temporal patterns of replication in all cell types, consistent with varying degrees of initiation efficiency. Zones of replication initiation (IZs) were detected throughout S phase and interacted in 3D space preferentially with other IZs of similar firing time. Temporal transition regions were resolved into segments of uni-directional replication punctuated at specific sites by small, inefficient IZs. Sites of convergent replication were divided into sites of termination or large constant timing regions consisting of many synchronous IZs in tandem. Developmental transitions in RT occured mainly by activating or inactivating individual IZs or occasionally by altering IZ firing time, demonstrating that IZs, rather than individual origins, are the units of developmental regulation. Finally, haplotype phasing revealed numerous regions of allele-specific and allele-independent asynchronous replication. Allele-independent asynchronous replication was correlated with the presence of previously mapped common fragile sites.

CONCLUSIONS

Altogether, these data provide a detailed temporal choreography of DNA replication in mammalian cells.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

DNA copy-number measurement of genome replication dynamics by high-throughput sequencing: the sort-seq, sync-seq and MFA-seq family

Genome replication follows a defined temporal programme that can change during cellular differentiation and disease onset. DNA replication results in an increase in DNA copy number that can be measured by high-throughput sequencing. Here we present a protocol to determine genome replication dynamics using DNA copy-number measurements. Cell populations can be obtained in three variants of the method. First, sort-seq reveals the average replication dynamics across S phase in an unperturbed cell population; FACS is used to isolate replicating and non-replicating subpopulations from asynchronous cells. Second, sync-seq measures absolute replication time at specific points during S phase using a synchronized cell population. Third, marker frequency analysis can be used to reveal the average replication dynamics using copy-number analysis in any proliferating asynchronous cell culture. These approaches have been used to reveal genome replication dynamics in prokaryotes, archaea and a wide range of eukaryotes, including yeasts and mammalian cells. We have found this approach straightforward to apply to other organisms and highlight example studies from across the three domains of life. Here we present a Saccharomyces cerevisiae version of the protocol that can be performed in 7-10 d. It requires basic molecular and cellular biology skills, as well as a basic understanding of Unix and R.

The Emerging Roles of Fox Family Transcription Factors in Chromosome Replication, Organization, and Genome Stability

The forkhead box (Fox) transcription factors (TFs) are widespread from yeast to humans. Their mutations and dysregulation have been linked to a broad spectrum of malignant neoplasias. They are known as critical players in DNA repair, metabolism, cell cycle control, differentiation, and aging. Recent studies, especially those from the simple model eukaryotes, revealed unexpected contributions of Fox TFs in chromosome replication and organization. More importantly, besides functioning as a canonical TF in cell signaling cascades and gene expression, Fox TFs can directly participate in DNA replication and determine the global replication timing program in a transcription-independent mechanism. Yeast Fox TFs preferentially recruit the limiting replication factors to a subset of early origins on chromosome arms. Attributed to their dimerization capability and distinct DNA binding modes, Fkh1 and Fkh2 also promote the origin clustering and assemblage of replication elements (replication factories). They can mediate long-range intrachromosomal and interchromosomal interactions and thus regulate the four-dimensional chromosome organization. The novel aspects of Fox TFs reviewed here expand their roles in maintaining genome integrity and coordinating the multiple essential chromosome events. These will inevitably be translated to our knowledge and new treatment strategies of Fox TF-associated human diseases including cancer.

Resolution of Complex Issues in Genome Regulation and Cancer Requires Non-Linear and Network-Based Thermodynamics

The apparent lack of success in curing cancer that was evidenced in the last four decades of molecular medicine indicates the need for a global re-thinking both its nature and the biological approaches that we are taking in its solution. The reductionist, one gene/one protein method that has served us well until now, and that still dominates in biomedicine, requires complementation with a more systemic/holistic approach, to address the huge problem of cross-talk between more than 20,000 protein-coding genes, about 100,000 protein types, and the multiple layers of biological organization. In this perspective, the relationship between the chromatin network organization and gene expression regulation plays a fundamental role. The elucidation of such a relationship requires a non-linear thermodynamics approach to these biological systems. This change of perspective is a necessary step for developing successful 'tumour-reversion' therapeutic strategies.

Checkpoint inhibition of origin firing prevents DNA topological stress

A universal feature of DNA damage and replication stress in eukaryotes is the activation of a checkpoint-kinase response. In S-phase, the checkpoint inhibits replication initiation, yet the function of this global block to origin firing remains unknown. To establish the physiological roles of this arm of the checkpoint, we analyzed separation of function mutants in the budding yeast Saccharomyces cerevisiae that allow global origin firing upon replication stress, despite an otherwise normal checkpoint response. Using genetic screens, we show that lack of the checkpoint-block to origin firing results in a dependence on pathways required for the resolution of topological problems. Failure to inhibit replication initiation indeed causes increased DNA catenation, resulting in DNA damage and chromosome loss. We further show that such topological stress is not only a consequence of a failed checkpoint response but also occurs in an unperturbed S-phase when too many origins fire simultaneously. Together we reveal that the role of limiting the number of replication initiation events is to prevent DNA topological problems, which may be relevant for the treatment of cancer with both topoisomerase and checkpoint inhibitors.

The effects of manipulating levels of replication initiation factors on origin firing efficiency in yeast

Chromosome replication in Saccharomyces cerevisiae is initiated from ~300 origins that are regulated by DNA sequence and by the limited abundance of six trans-acting initiation proteins (Sld2, Sld3, Dpb11, Dbf4, Sld7 and Cdc45). We set out to determine how the levels of individual factors contribute to time of origin activation and/or origin efficiency using induced depletion of single factors and overexpression of sets of multiple factors. Depletion of Sld2 or Sld3 slows growth and S phase progression, decreases origin efficiency across the genome and impairs viability as a result of incomplete replication of the rDNA. We find that the most efficient early origins are relatively unaffected by depletion of either Sld2 or Sld3. However, Sld3 levels, and to a lesser extent Sld2 levels, are critical for firing of the less efficient early origins. Overexpression of Sld3 simultaneously with Sld2, Dpb11 and Dbf4 preserves the relative efficiency of origins. Only when Cdc45 and Sld7 are also overexpressed is origin efficiency equalized between early- and late-firing origins. Our data support a model in which Sld3 together with Cdc45 (and/or Sld7) is responsible for the differential efficiencies of origins across the yeast genome.

Eukaryotic chromosome duplication begins at sites called origins of replication along the chromosomal DNA. A conserved property of eukaryotic origins is that they vary in efficiency—the proportion of cells in a population in which they “fire”—and in the average time of activation within S phase, but the molecular details underlying this variation are not well understood. Previous work has shown that limiting concentrations of a set of conserved replication initiation proteins referred to as “SSDDCS” (Sld2, Sld3, Dbf4, Dpb11, Cdc45, and Sld7) are rate limiting for origin activation in budding yeast and other eukaryotes; combined overexpression of these proteins increases and/or advances origin firing. However, it remained possible that different factors affect different aspects of origin activation (e.g., timing vs. efficiency). Here, by depleting individual factors or by overexpressing sets of factors in budding yeast, we demonstrate that it is levels of Sld3, Cdc45 and/or Sld7 levels are primarily responsible for modulating the differences in relative origin efficiency and timing. This work gives further insights into what shapes the landscape of genome duplication.

Peripheral re-localization of constitutive heterochromatin advances its replication timing and impairs maintenance of silencing marks

The replication of the genome is a highly organized process, both spatially and temporally. Although a lot is known on the composition of the basic replication machinery, how its activity is regulated is mostly unknown. Several chromatin properties have been proposed as regulators, but a potential role of the nuclear DNA position remains unclear. We made use of the prominent structure and well-defined heterochromatic landscape of mouse pericentric chromosome domains as a well-studied example of late replicating constitutive heterochromatin. We established a method to manipulate its nuclear position and evaluated the effect on replication timing, DNA compaction and epigenetic composition. Using time-lapse microscopy, we observed that constitutive heterochromatin, known to replicate during late S-phase, was replicated in mid S-phase when repositioned to the nuclear periphery. Out-of-schedule replication resulted in deficient post-replicative maintenance of chromatin modifications, namely silencing marks. We propose that repositioned constitutive heterochromatin was activated in trans according to the domino model of origin firing by nearby (mid S) firing origins. In summary, our data provide, on the one hand, a novel approach to manipulate nuclear DNA position and, on the other hand, establish nuclear DNA position as a novel mechanism regulating DNA replication timing and epigenetic maintenance.

Comparing 3D Genome Organization in Multiple Species Using Phylo-HMRF

Recent whole-genome mapping approaches for the chromatin interactome have offered new insights into 3D genome organization. However, our knowledge of the evolutionary patterns of 3D genome in mammals remains limited. In particular, there are no existing phylogenetic-model-based methods to analyze chromatin interactions as continuous features. Here, we develop phylogenetic hidden Markov random field (Phylo-HMRF) to identify evolutionary patterns of 3D genome based on multi-species Hi-C data by jointly utilizing spatial constraints among genomic loci and continuous-trait evolutionary models. We used Phylo-HMRF to uncover cross-species 3D genome patterns based on Hi-C data from the same cell type in four primate species (human, chimpanzee, bonobo, and gorilla). The identified evolutionary patterns of 3D genome correlate with features of genome structure and function. This work provides a new framework to analyze multi-species continuous genomic features with spatial constraints and has the potential to help reveal the evolutionary principles of 3D genome organization.

Nascent chromatin occupancy profiling reveals locus- and factor-specific chromatin maturation dynamics behind the DNA replication fork

Proper regulation and maintenance of the epigenome is necessary to preserve genome function. However, in every cell division, the epigenetic state is disassembled and then reassembled in the wake of the DNA replication fork. Chromatin restoration on nascent DNA is a complex and regulated process that includes nucleosome assembly and remodeling, deposition of histone variants, and the re-establishment of transcription factor binding. To study the genome-wide dynamics of chromatin restoration behind the DNA replication fork, we developed nascent chromatin occupancy profiles (NCOPs) to comprehensively profile nascent and mature chromatin at nucleotide resolution. Although nascent chromatin is inherently less organized than mature chromatin, we identified locus-specific differences in the kinetics of chromatin maturation that were predicted by the epigenetic landscape, including the histone variant H2AZ, which marked loci with rapid maturation kinetics. The chromatin maturation at origins of DNA replication was dependent on whether the origin underwent initiation or was passively replicated from distal-originating replication forks, suggesting distinct chromatin assembly mechanisms surrounding activated and disassembled prereplicative complexes. Finally, we identified sites that were only occupied transiently by DNA-binding factors following passage of the replication fork, which may provide a mechanism for perturbations of the DNA replication program to shape the regulatory landscape of the genome.

Metazoan DNA replication origins

DNA replication starts with the opening of DNA at sites called DNA replication origins. From the single sequence-specific DNA replication origin of the small Escherichia coli genome, up to thousands of origins that are necessary to replicate the large human genome, strict sequence specificity has been lost. Nevertheless, genome-wide analyses performed in the recent years, using different mapping methods, demonstrated that there are precise locations along the metazoan genome from which replication initiates. These sites contain relaxed sequence consensus and epigenetic features. There is flexibility in the choice of origins to be used during a given cell cycle, probably imposed by evolution and developmental constraints. Here, we will briefly describe their main features.

Monitoring early S-phase origin firing and replication fork movement by sequencing nascent DNA from synchronized cells

A better understanding of DNA replication initiation in human cells and how this process is altered upon DNA replication stress requires the ability to study origin firing genome wide. Previously described methods of mapping DNA replication origins in higher eukaryotes rely principally on fractionation of DNA fragments based on their size and, optionally, on the presence of ribonucleotides at their 5' end. Here, we describe a protocol for EdUseq-HU, a method for mapping early S-phase replication origins. Cells, synchronized by mitotic shake-off, are released in medium containing 5-ethynyl-2'-deoxyuridine (EdU; to label nascent DNA) and hydroxyurea (HU; to limit fork progression after origin firing). After using click chemistry to tag the EdU label with a biotin conjugate that is cleavable under mild conditions, the nascent DNA is captured on streptavidin beads. One variant of EdUseq-HU allows mapping of DNA replication origins on the genome at a resolution of 10 kb, and a second variant monitors progression of replication forks. Using EdUseq-HU, the spatiotemporal program of DNA replication in human cell lines can be interrogated in <2 weeks. The protocol requires basic cell culture and molecular biology skills, as well as familiarity with the Perl programming language and the Linux operating system.

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

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

Outer kinetochore protein Dam1 promotes centromere clustering in parallel with Slk19 in budding yeast

A higher order organization of the centromeres in the form of clustering of these DNA loci has been observed in many organisms. While centromere clustering is biologically significant to achieve faithful chromosome segregation, the underlying molecular mechanism is yet to be fully understood. In budding yeast, a kinetochore-associated protein Slk19 is shown to have a role in clustering in association with the microtubules whereas removal of either Slk19 or microtubules alone does not have any effect on the centromere clustering. Furthermore, Slk19 is non-essential for growth and becomes cleaved during anaphase whereas clustering being an essential event occurs throughout the cell cycle. Hence, we searched for an additional factor involved in the clustering and since the integrity of the kinetochore complex is shown to be crucial for centromere clustering, we restricted our search within the complex. We observed that the outermost kinetochore protein Dam1 promotes centromere clustering through stabilization of the kinetochore integrity. While in the absence of Dam1 we failed to detect Slk19 at the centromere, on the other hand, we found almost no Dam1 at the centromere in the absence of Slk19 and microtubules suggesting interdependency between these two pathways. Strikingly, we observed that overexpression of Dam1 or Slk19 could restore the centromere clustering largely in the cells devoid of Slk19 and microtubules or Dam1, respectively. Thus, we propose that in budding yeast, centromere clustering is achieved at least by two parallel pathways, through Dam1 and another via Slk19, in concert with the microtubules suggesting that having a dual mechanism may be crucial for ensuring microtubule capture by the point centromeres where each attaches to only one microtubule.

Systematic Evaluation of Statistical Methods for Identifying Looping Interactions in 5C Data

Chromosome-Conformation-Capture-Carbon-Copy (5C) is a molecular technology based on proximity ligation that enables high-resolution and high-coverage inquiry of long-range looping interactions. Computational pipelines for analyzing 5C data involve a series of interdependent normalization procedures and statistical methods that markedly influence downstream biological results. A detailed analysis of the trade-offs inherent to all stages of 5C data analysis has not been reported. Here, we provide a comparative assessment of method performance at each step in the 5C analysis pipeline, including sequencing depth and library complexity correction, bias mitigation, spatial noise reduction, distance-dependent expected and variance estimation, statistical modeling, and loop detection. We discuss methodological advantages and disadvantages at each step and provide a full suite of algorithms, lib5C, to allow investigators to test the range of approaches on their own 5C data. Principles learned from our comparative analyses can be applied to protein-independent proximity ligation-based data, including Hi-C, 4C, and Capture-C.

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.

Higher order genomic organization and epigenetic control maintain cellular identity and prevent breast cancer

Cells establish and sustain structural and functional integrity of the genome to support cellular identity and prevent malignant transformation. In this review, we present a strategic overview of epigenetic regulatory mechanisms including histone modifications and higher order chromatin organization (HCO) that are perturbed in breast cancer onset and progression. Implications for dysfunctions that occur in hormone regulation, cell cycle control, and mitotic bookmarking in breast cancer are considered, with an emphasis on epithelial-to-mesenchymal transition and cancer stem cell activities. The architectural organization of regulatory machinery is addressed within the contexts of translating cancer-compromised genomic organization to advances in breast cancer risk assessment, diagnosis, prognosis, and identification of novel therapeutic targets with high specificity and minimal off target effects.

DNA replication dynamics of vole genome and its epigenetic regulation

Background

The genome of some vole rodents exhibit large blocks of heterochromatin coupled to their sex chromosomes. The DNA composition and transcriptional activity of these heterochromatin blocks have been studied, but little is known about their DNA replication dynamics and epigenetic composition.

Results

Here, we show prominent epigenetic marks of the heterochromatic blocks in the giant sex chromosomes of female Microtus cabrerae cells. While the X chromosomes are hypoacetylated and cytosine hypomethylated, they are either enriched for macroH2A and H3K27me3 typical for facultative heterochromatin or for H3K9me3 and HP1 beta typical for constitutive heterochromatin. Using pulse-chase replication labeling and time-lapse microscopy, we found that the heterochromatic block enriched for macroH2A/H3K27me3 of the X chromosome is replicated during mid-S-phase, prior to the heterochromatic block enriched for H3K9me3/HP1 beta, which is replicated during late S-phase. To test whether histone acetylation level regulates its replication dynamics, we induced either global hyperacetylation by pharmacological inhibition or by targeting a histone acetyltransferase to the heterochromatic region of the X chromosomes. Our data reveal that histone acetylation level affects DNA replication dynamics of the sex chromosomes’ heterochromatin and leads to a global reduction in replication fork rate genome wide.

Conclusions

In conclusion, we mapped major epigenetic modifications controlling the structure of the sex chromosome-associated heterochromatin and demonstrated the occurrence of differences in the molecular mechanisms controlling the replication timing of the heterochromatic blocks at the sex chromosomes in female Microtus cabrerae cells. Furthermore, we highlighted a conserved role of histone acetylation level on replication dynamics across mammalian species.

Electronic supplementary material

The online version of this article (10.1186/s13072-019-0262-0) contains supplementary material, which is available to authorized users.

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

Here, we report a single-cell DNA replication sequencing method, scRepli-seq, a 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 profiles, which were also conserved among cells. Haplotype-resolved scRepli-seq revealed similar replication 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, which was smallest at the beginning and the end of S phase. In addition, developmentally regulated domains were found to deviate from others and showed a higher degree of heterogeneity, thus suggesting a link to developmental plasticity. Moreover, allelic expression imbalance was found to strongly associate with replication-timing asynchrony. Our results form a foundation for single-cell-level understanding of DNA replication regulation and provide insights into three-dimensional genome organization.

Dynamics of DNA replication in a eukaryotic cell

Each genomic locus in a eukaryotic cell has a distinct average time of replication during S phase that depends on the spatial and temporal pattern of replication initiation events. Replication timing can affect genomic integrity because late replication is associated with an increased mutation rate. For most eukaryotes, the features of the genome that specify the location and timing of initiation events are unknown. To investigate these features for the fission yeast, Schizosaccharomyces pombe, we developed an integrative model to analyze large single-molecule and global genomic datasets. The model provides an accurate description of the complex dynamics of S. pombe DNA replication at high resolution. We present evidence that there are many more potential initiation sites in the S. pombe genome than previously identified and that the distribution of these sites is primarily determined by two factors: the sequence preferences of the origin recognition complex (ORC), and the interference of transcription with the assembly or stability of prereplication complexes (pre-RCs). We suggest that in addition to directly interfering with initiation, transcription has driven the evolution of the binding properties of ORC in S. pombe and other eukaryotic species to target pre-RC assembly to regions of the genome that are less likely to be transcribed.

On the Interplay of the DNA Replication Program and the Intra-S Phase Checkpoint Pathway

DNA replication in eukaryotes is achieved by the activation of multiple replication origins which needs to be precisely coordinated in space and time. This spatio-temporal replication program is regulated by many factors to maintain genome stability, which is frequently threatened through stresses of exogenous or endogenous origin. Intra-S phase checkpoints monitor the integrity of DNA synthesis and are activated when replication forks are stalled. Their activation leads to the stabilization of forks, to the delay of the replication program by the inhibition of late firing origins, and the delay of G2/M phase entry. In some cell cycles during early development these mechanisms are less efficient in order to allow rapid cell divisions. In this article, we will review our current knowledge of how the intra-S phase checkpoint regulates the replication program in budding yeast and metazoan models, including early embryos with rapid S phases. We sum up current models on how the checkpoint can inhibit origin firing in some genomic regions, but allow dormant origin activation in other regions. Finally, we discuss how numerical and theoretical models can be used to connect the multiple different actors into a global process and to extract general rules.

DNA fragility in the parallel evolution of pelvic reduction in stickleback fish

Evolution generates a remarkable breadth of living forms, but many traits evolve repeatedly, by mechanisms that are still poorly understood. A classic example of repeated evolution is the loss of pelvic hindfins in stickleback fish (Gasterosteus aculeatus). Repeated pelvic loss maps to recurrent deletions of a pelvic enhancer of the Pitx1 gene. Here, we identify molecular features contributing to these recurrent deletions. Pitx1 enhancer sequences form alternative DNA structures in vitro and increase double-strand breaks and deletions in vivo. Enhancer mutability depends on DNA replication direction and is caused by TG-dinucleotide repeats. Modeling shows that elevated mutation rates can influence evolution under demographic conditions relevant for sticklebacks and humans. DNA fragility may thus help explain why the same loci are often used repeatedly during parallel adaptive evolution.

Resolving Roadblocks to Telomere Replication

The maintenance of genome stability in eukaryotic cells relies on accurate and efficient replication along each chromosome following every cell division. The terminal position, repetitive sequence, and structural complexities of the telomeric DNA make the telomere an inherently difficult region to replicate within the genome. Thus, despite functioning to protect genome stability mammalian telomeres are also a source of replication stress and have been recognized as common fragile sites within the genome. Telomere fragility is exacerbated at telomeres that rely on the Alternative Lengthening of Telomeres (ALT) pathway. Like common fragile sites, ALT telomeres are prone to chromosome breaks and are frequent sites of recombination suggesting that ALT telomeres are subjected to chronic replication stress. Here, we will review the features of telomeric DNA that challenge the replication machinery and also how the cell overcomes these challenges to maintain telomere stability and ensure the faithful duplication of the human genome.

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.

Evaluation of mechanisms that may generate DNA lesions triggering antigenic variation in African trypanosomes

Antigenic variation by variant surface glycoprotein (VSG) coat switching in African trypanosomes is one of the most elaborate immune evasion strategies found among pathogens. Changes in the identity of the transcribed VSG gene, which is always flanked by 70-bp and telomeric repeats, can be achieved either by transcriptional or DNA recombination mechanisms. The major route of VSG switching is DNA recombination, which occurs in the bloodstream VSG expression site (ES), a multigenic site transcribed by RNA polymerase I. Recombinogenic VSG switching is frequently catalyzed by homologous recombination (HR), a reaction normally triggered by DNA breaks. However, a clear understanding of how such breaks arise-including whether there is a dedicated and ES-focused mechanism-is lacking. Here, we synthesize data emerging from recent studies that have proposed a range of mechanisms that could generate these breaks: action of a nuclease or nucleases; repetitive DNA, most notably the 70-bp repeats, providing an intra-ES source of instability; DNA breaks derived from the VSG-adjacent telomere; DNA breaks arising from high transcription levels at the active ES; and DNA lesions arising from replication-transcription conflicts in the ES. We discuss the evidence that underpins these switch-initiation models and consider what features and mechanisms might be shared or might allow the models to be tested further. Evaluation of all these models highlights that we still have much to learn about the earliest acting step in VSG switching, which may have the greatest potential for therapeutic intervention in order to undermine the key reaction used by trypanosomes for their survival and propagation in the mammalian host.

Rif1 inhibits replication fork progression and controls DNA copy number in Drosophila

Control of DNA copy number is essential to maintain genome stability and ensure proper cell and tissue function. In Drosophila polyploid cells, the SNF2-domain-containing SUUR protein inhibits replication fork progression within specific regions of the genome to promote DNA underreplication. While dissecting the function of SUUR's SNF2 domain, we identified an interaction between SUUR and Rif1. Rif1 has many roles in DNA metabolism and regulates the replication timing program. We demonstrate that repression of DNA replication is dependent on Rif1. Rif1 localizes to active replication forks in a partially SUUR-dependent manner and directly regulates replication fork progression. Importantly, SUUR associates with replication forks in the absence of Rif1, indicating that Rif1 acts downstream of SUUR to inhibit fork progression. Our findings uncover an unrecognized function of the Rif1 protein as a regulator of replication fork progression.

Chromatin conformation and transcriptional activity are permissive regulators of DNA replication initiation in Drosophila

Chromatin structure has emerged as a key contributor to spatial and temporal control over the initiation of DNA replication. However, despite genome-wide correlations between early replication of gene-rich, accessible euchromatin and late replication of gene-poor, inaccessible heterochromatin, a causal relationship between chromatin structure and replication initiation remains elusive. Here, we combined histone gene engineering and whole-genome sequencing in Drosophila to determine how perturbing chromatin structure affects replication initiation. We found that most pericentric heterochromatin remains late replicating in H3K9R mutants, even though H3K9R pericentric heterochromatin is depleted of HP1a, more accessible, and transcriptionally active. These data indicate that HP1a loss, increased chromatin accessibility, and elevated transcription do not result in early replication of heterochromatin. Nevertheless, a small amount of pericentric heterochromatin with increased accessibility replicates earlier in H3K9R mutants. Transcription is de-repressed in these regions of advanced replication but not in those regions of the H3K9R mutant genome that replicate later, suggesting that transcriptional repression may contribute to late replication. We also explored relationships among chromatin, transcription, and replication in euchromatin by analyzing H4K16R mutants. In Drosophila, the X Chromosome gene expression is up-regulated twofold and replicates earlier in XY males than it does in XX females. We found that H4K16R mutation prevents normal male development and abrogates hyperexpression and earlier replication of the male X, consistent with previously established genome-wide correlations between transcription and early replication. In contrast, H4K16R females are viable and fertile, indicating that H4K16 modification is dispensable for genome replication and gene expression.

Integrative detection and analysis of structural variation in cancer genomes

Structural variants (SVs) can contribute to oncogenesis through a variety of mechanisms. Despite their importance, the identification of SVs in cancer genomes remains challenging. Here, we present a framework that integrates optical mapping, high-throughput chromosome conformation capture (Hi-C), and whole-genome sequencing to systematically detect SVs in a variety of normal or cancer samples and cell lines. We identify the unique strengths of each method and demonstrate that only integrative approaches can comprehensively identify SVs in the genome. By combining Hi-C and optical mapping, we resolve complex SVs and phase multiple SV events to a single haplotype. Furthermore, we observe widespread structural variation events affecting the functions of noncoding sequences, including the deletion of distal regulatory sequences, alteration of DNA replication timing, and the creation of novel three-dimensional chromatin structural domains. Our results indicate that noncoding SVs may be underappreciated mutational drivers in cancer genomes.

Recent advances in understanding DNA replication: cell type-specific adaptation of the DNA replication program

DNA replication is an essential process occurring prior to cell division. Cell division coupled to proliferation ensures the growth and renewal of a large variety of specialized cell types generated during embryonic development. Changes in the DNA replication program occur during development. Embryonic undifferentiated cells show a high replication rate and fast proliferation, whereas more differentiated cells are characterized by reduced DNA synthesis and a low proliferation rate. Hence, the DNA replication program must adapt to the specific features of cells committed to different fates. Recent findings on DNA synthesis regulation in different cell types open new perspectives for developing efficient and more adapted therapies to treat various diseases such as genetic diseases and cancer. This review will put the emphasis on recent progress made in this field.

Dynamic changes in ORC localization and replication fork progression during tissue differentiation

Background

Genomic regions repressed for DNA replication, resulting in either delayed replication in S phase or underreplication in polyploid cells, are thought to be controlled by inhibition of replication origin activation. Studies in Drosophila polytene cells, however, raised the possibility that impeding replication fork progression also plays a major role.

Results

We exploited genomic regions underreplicated (URs) with tissue specificity in Drosophila polytene cells to analyze mechanisms of replication repression. By localizing the Origin Recognition Complex (ORC) in the genome of the larval fat body and comparing this to ORC binding in the salivary gland, we found that sites of ORC binding show extensive tissue specificity. In contrast, there are common domains nearly devoid of ORC in the salivary gland and fat body that also have reduced density of ORC binding sites in diploid cells. Strikingly, domains lacking ORC can still be replicated in some polytene tissues, showing absence of ORC and origins is insufficient to repress replication. Analysis of the width and location of the URs with respect to ORC position indicates that whether or not a genomic region lacking ORC is replicated is controlled by whether replication forks formed outside the region are inhibited.

Conclusions

These studies demonstrate that inhibition of replication fork progression can block replication across genomic regions that constitutively lack ORC. Replication fork progression can be inhibited in both tissue-specific and genome region-specific ways. Consequently, when evaluating sources of genome instability it is important to consider altered control of replication forks in response to differentiation.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4992-3) contains supplementary material, which is available to authorized users.

Continuous-Trait Probabilistic Model for Comparing Multi-species Functional Genomic Data

A large amount of multi-species functional genomic data from high-throughput assays are becoming available to help understand the molecular mechanisms for phenotypic diversity across species. However, continuous-trait probabilistic models, which are key to such comparative analysis, remain under-explored. Here we develop a new model, called phylogenetic hidden Markov Gaussian processes (Phylo-HMGP), to simultaneously infer heterogeneous evolutionary states of functional genomic features in a genome-wide manner. Both simulation studies and real data application demonstrate the effectiveness of Phylo-HMGP. Importantly, we applied Phylo-HMGP to analyze a new cross-species DNA replication timing (RT) dataset from the same cell type in five primate species (human, chimpanzee, orangutan, gibbon, and green monkey). We demonstrate that our Phylo-HMGP model enables discovery of genomic regions with distinct evolutionary patterns of RT. Our method provides a generic framework for comparative analysis of multi-species continuous functional genomic signals to help reveal regions with conserved or lineage-specific regulatory roles.

The organization of genome duplication is a critical determinant of the landscape of genome maintenance

Genome duplication is essential for cell proliferation, and the mechanisms regulating its execution are highly conserved. These processes give rise to a spatiotemporal organization of replication initiation across the genome, referred to as the replication program. Despite the identification of such programs in diverse eukaryotic organisms, their biological importance for cellular physiology remains largely unexplored. We address this fundamental question in the context of genome maintenance, taking advantage of the inappropriate origin firing that occurs when fission yeast cells lacking the Rad3/ATR checkpoint kinase are subjected to replication stress. Using this model, we demonstrate that the replication program quantitatively dictates the extent of origin de-regulation and the clustered localization of these events. Furthermore, our results uncover an accumulation of abnormal levels of single-stranded DNA (ssDNA) and the Rad52 repair protein at de-regulated origins. We show that these loci constitute a defining source of the overall ssDNA and Rad52 hotspots in the genome, generating a signature pattern of instability along the chromosomes. We then induce a genome-wide reprogramming of origin usage and evaluate its consequences in our experimental system. This leads to a complete redistribution of the sites of both inappropriate initiation and associated Rad52 recruitment. We therefore conclude that the organization of genome duplication governs the checkpoint control of origin-associated hotspots of instability and plays an integral role in shaping the landscape of genome maintenance.

Rif1 prolongs the embryonic S phase at the Drosophila mid-blastula transition

In preparation for dramatic morphogenetic events of gastrulation, rapid embryonic cell cycles slow at the mid-blastula transition (MBT). In Drosophila melanogaster embryos, down-regulation of cyclin-dependent kinase 1 (Cdk1) activity initiates this slowing by delaying replication of heterochromatic satellite sequences and extending S phase. We found that Cdk1 activity inhibited the chromatin association of Rap1 interacting factor 1 (Rif1), a candidate repressor of replication. Furthermore, Rif1 bound selectively to satellite sequences following Cdk1 down-regulation at the MBT. In the next S phase, Rif1 dissociated from different satellites in an orderly schedule that anticipated their replication. Rif1 lacking potential phosphorylation sites failed to dissociate and dominantly prevented completion of replication. Loss of Rif1 in mutant embryos shortened the post-MBT S phase and rescued embryonic cell cycles disrupted by depletion of the S phase–promoting kinase, cell division cycle 7 (Cdc7). Our work shows that Rif1 and S phase kinases compose a replication timer controlling first the developmental onset of late replication and then the precise schedule of replication within S phase. In addition, we describe how onset of late replication fits into the progressive maturation of heterochromatin during development.

Cells divide rapidly in the early embryos of most animals. However, during a conserved period of development known as the mid-blastula transition (MBT), the cell cycle slows down dramatically. In Drosophila embryos, genome duplication abruptly slows to initiate this cell cycle prolongation. This is achieved through the onset of late replication, a well-recognized phenomenon in which specific sequences of the genome await replication until long after other sequences have finished. Even though this temporal program of replication is a major determinant of the duration of S phase, the factors involved in this process remain unknown. Here, we use genetics and real-time microscopy to visualize replication in developing fly embryos and show that the protein Rap1 interacting factor 1 (Rif1) mediates the introduction of late replication at the MBT. We find that at this stage, Rif1 binds to and selectively delays the replication of large blocks of repetitive DNA known as satellite sequences. During the rapid cell cycles before the MBT, we show that the cyclin-dependent kinase 1 (Cdk1) prevents Rif1 from slowing down DNA replication by driving its removal from the chromatin. The developmental down-regulation of Cdk1 at the MBT allows Rif1 to associate with the satellite sequences and initiate cell cycle slowing. Our work provides new insights into the temporal programming of S phase and into the embryonic origin of late replication.

Lgr5+ intestinal stem cells reside in an unlicensed G1 phase

During late mitosis and the early G1 phase, the origins of replication are licensed by binding to double hexamers of MCM2-7. In this study, we investigated how licensing and proliferative commitment are coupled in the epithelium of the small intestine. We developed a method for identifying cells in intact tissue containing DNA-bound MCM2-7. Interphase cells above the transit-amplifying compartment had no DNA-bound MCM2-7, but still expressed the MCM2-7 protein, suggesting that licensing is inhibited immediately upon differentiation. Strikingly, we found most proliferative Lgr5+ stem cells are in an unlicensed state. This suggests that the elongated cell-cycle of intestinal stem cells is caused by an increased G1 length, characterized by dormant periods with unlicensed origins. Significantly, the unlicensed state is lost in Apc-mutant epithelium, which lacks a functional restriction point, causing licensing immediately upon G1 entry. We propose that the unlicensed G1 phase of intestinal stem cells creates a temporal window when proliferative fate decisions can be made.

Eukaryotic DNA replication: Orchestrated action of multi-subunit protein complexes

Genome duplication is an essential process to preserve genetic information between generations. The eukaryotic cell cycle is composed of functionally distinct phases: G1, S, G2, and M. One of the key replicative proteins that participate at every stage of DNA replication is the Mcm2-7 complex, a replicative helicase. In the G1 phase, inactive Mcm2-7 complexes are loaded on the replication origins by replication-initiator proteins, ORC and Cdc6. Two kinases, S-CDK and DDK, convert the inactive origin-loaded Mcm2-7 complex to an active helicase, the CMG complex in the S phase. The activated CMG complex begins DNA unwinding and recruits enzymes essential for DNA synthesis to assemble a replisome at the replication fork. After completion of DNA synthesis, the inactive CMG complex on the replicated DNA is removed from chromatin to terminate DNA replication. In this review, we will discuss the structure, function, and regulation of the molecular machines involved in each step of DNA replication.

Profiling DNA Replication Timing Using Zebrafish as an In Vivo Model System

DNA replication timing is an important cellular characteristic, exhibiting significant relationships with chromatin structure, transcription, and DNA mutation rates. Changes in replication timing occur during development and in cancer, but the role replication timing plays in development and disease is not known. Zebrafish were recently established as an in vivo model system to study replication timing. Here is detailed the protocols for using the zebrafish to determine DNA replication timing. After sorting cells from embryos and adult zebrafish, high-resolution genome-wide DNA replication timing patterns can be constructed by determining changes in DNA copy number through analysis of next generation sequencing data. The zebrafish model system allows for evaluation of the replication timing changes that occur in vivo throughout development, and can also be used to assess changes in individual cell types, disease models, or mutant lines. These methods will enable studies investigating the mechanisms and determinants of replication timing establishment and maintenance during development, the role replication timing plays in mutations and tumorigenesis, and the effects of perturbing replication timing on development and disease.

A site specific model and analysis of the neutral somatic mutation rate in whole-genome cancer data

Background

Detailed modelling of the neutral mutational process in cancer cells is crucial for identifying driver mutations and understanding the mutational mechanisms that act during cancer development. The neutral mutational process is very complex: whole-genome analyses have revealed that the mutation rate differs between cancer types, between patients and along the genome depending on the genetic and epigenetic context. Therefore, methods that predict the number of different types of mutations in regions or specific genomic elements must consider local genomic explanatory variables. A major drawback of most methods is the need to average the explanatory variables across the entire region or genomic element. This procedure is particularly problematic if the explanatory variable varies dramatically in the element under consideration.

Results

To take into account the fine scale of the explanatory variables, we model the probabilities of different types of mutations for each position in the genome by multinomial logistic regression. We analyse 505 cancer genomes from 14 different cancer types and compare the performance in predicting mutation rate for both regional based models and site-specific models. We show that for 1000 randomly selected genomic positions, the site-specific model predicts the mutation rate much better than regional based models.

We use a forward selection procedure to identify the most important explanatory variables. The procedure identifies site-specific conservation (phyloP), replication timing, and expression level as the best predictors for the mutation rate. Finally, our model confirms and quantifies certain well-known mutational signatures.

Conclusion

We find that our site-specific multinomial regression model outperforms the regional based models. The possibility of including genomic variables on different scales and patient specific variables makes it a versatile framework for studying different mutational mechanisms. Our model can serve as the neutral null model for the mutational process; regions that deviate from the null model are candidates for elements that drive cancer development.

Electronic supplementary material

The online version of this article (10.1186/s12859-018-2141-2) contains supplementary material, which is available to authorized users.

Similarity in replication timing between polytene and diploid cells is associated with the organization of the Drosophila genome

Morphologically, polytene chromosomes of Drosophila melanogaster consist of compact “black” bands alternating with less compact “grey” bands and interbands. We developed a comprehensive approach that combines cytological mapping data of FlyBase-annotated genes and novel tools for predicting cytogenetic features of chromosomes on the basis of their protein composition and determined the genomic coordinates for all black bands of polytene chromosome 2R. By a PCNA immunostaining assay, we obtained the replication timetable for all the bands mapped. The results allowed us to compare replication timing between polytene chromosomes in salivary glands and chromosomes from cultured diploid cell lines and to observe a substantial similarity in the global replication patterns at the band resolution level. In both kinds of chromosomes, the intervals between black bands correspond to early replication initiation zones. Black bands are depleted of replication initiation events and are characterized by a gradient of replication timing; therefore, the time of replication completion correlates with the band length. The bands are characterized by low gene density, contain predominantly tissue-specific genes, and are represented by silent chromatin types in various tissues. The borders of black bands correspond well to the borders of topological domains as well as to the borders of the zones showing H3K27me3, SUUR, and LAMIN enrichment. In conclusion, the characteristic pattern of polytene chromosomes reflects partitioning of the Drosophila genome into two global types of domains with contrasting properties. This partitioning is conserved in different tissues and determines replication timing in Drosophila.

DNA replication through a chromatin environment

Compaction of the genome into the nuclear space is achieved by wrapping DNA around octameric assemblies of histone proteins to form nucleosomes, the fundamental repeating unit of chromatin. Aside from providing a means by which to fit larger genomes into the cell, chromatinization of DNA is a crucial means by which the cell regulates access to the genome. While the complex role that chromatin plays in gene transcription has been appreciated for a long time, it is now also apparent that crucial aspects of DNA replication are linked to the biology of chromatin. This review will focus on recent advances in our understanding of how the chromatin environment influences key aspects of DNA replication.This article is part of the themed issue 'Chromatin modifiers and remodellers in DNA repair and signalling'.

A mechanism for epigenetic control of DNA replication

DNA replication origins in hyperacetylated euchromatin fire preferentially during early S phase. However, how acetylation controls DNA replication timing is unknown. TICRR/TRESLIN is an essential protein required for the initiation of DNA replication. Here, we report that TICRR physically interacts with the acetyl-histone binding bromodomain (BRD) and extraterminal (BET) proteins BRD2 and BRD4. Abrogation of this interaction impairs TICRR binding to acetylated chromatin and disrupts normal S-phase progression. Our data reveal a novel function for BET proteins and establish the TICRR-BET interaction as a potential mechanism for epigenetic control of DNA replication.

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

Fang et al. show that Dbf4 is enriched at early origins through its interaction with forkhead transcription factors Fkh1 and Fkh2. Dbf4 interacts directly with Sld3 and promotes the recruitment of downstream limiting factors.

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.