Regulation of DNA replication within the immunoglobulin heavy-chain locus during B cell commitment

Visualizing DNA replication by single-molecule analysis of replicated DNA

Single-molecule analysis of replicated DNA (SMARD) is a unique technique that enables visualization of DNA replication at specific genomic regions at single-molecule resolution. Here, we present a protocol for visualizing DNA replication by SMARD. We describe steps for pulse labeling DNA, followed by isolating and stretching of genomic DNA. We then detail the detection of the replication at chromosomal regions through immunostaining and fluorescence in situ hybridization. Using SMARD, we can visualize replication initiation, progression, termination, and fork stalling. For complete details on the use and execution of this protocol, please refer to Norio et al. (2001) and Gerhardt et al. (2014).1,2.

Where and when to start: Regulating DNA replication origin activity in eukaryotic genomes

In eukaryotic genomes, hundreds to thousands of potential start sites of DNA replication named origins are dispersed across each of the linear chromosomes. During S-phase, only a subset of origins is selected in a stochastic manner to assemble bidirectional replication forks and initiate DNA synthesis. Despite substantial progress in our understanding of this complex process, a comprehensive 'identity code' that defines origins based on specific nucleotide sequences, DNA structural features, the local chromatin environment, or 3D genome architecture is still missing. In this article, we review the genetic and epigenetic features of replication origins in yeast and metazoan chromosomes and highlight recent insights into how this flexibility in origin usage contributes to nuclear organization, cell growth, differentiation, and genome stability.

Genome-wide measurement of DNA replication fork directionality and quantification of DNA replication initiation and termination with Okazaki fragment sequencing

Studying the dynamics of genome replication in mammalian cells has been historically challenging. To reveal the location of replication initiation and termination in the human genome, we developed Okazaki fragment sequencing (OK-seq), a quantitative approach based on the isolation and strand-specific sequencing of Okazaki fragments, the lagging strand replication intermediates. OK-seq quantitates the proportion of leftward- and rightward-oriented forks at every genomic locus and reveals the location and efficiency of replication initiation and termination events. Here we provide the detailed experimental procedures for performing OK-seq in unperturbed cultured human cells and budding yeast and the bioinformatics pipelines for data processing and computation of replication fork directionality. Furthermore, we present the analytical approach based on a hidden Markov model, which allows automated detection of ascending, descending and flat replication fork directionality segments revealing the zones of replication initiation, termination and unidirectional fork movement across the entire genome. These tools are essential for the accurate interpretation of human and yeast replication programs. The experiments and the data processing can be accomplished within six days. Besides revealing the genome replication program in fine detail, OK-seq has been instrumental in numerous studies unravelling mechanisms of genome stability, epigenome maintenance and genome evolution.

Rif1-Dependent Control of Replication Timing

Successful duplication of the genome requires the accurate replication of billions of base pairs of DNA within a relatively short time frame. Failure to accurately replicate the genome results in genomic instability and a host of diseases. To faithfully and rapidly replicate the genome, DNA replication must be tightly regulated and coordinated with many other nuclear processes. These regulations, however, must also be flexible as replication kinetics can change through development and differentiation. Exactly how DNA replication is regulated and how this regulation changes through development is an active field of research. One aspect of genome duplication where much remains to be discovered is replication timing (RT), which dictates when each segment of the genome is replicated during S phase. All organisms display some level of RT, yet the precise mechanisms that govern RT remain are not fully understood. The study of Rif1, a protein that actively regulates RT from yeast to humans, provides a key to unlock the underlying molecular mechanisms controlling RT. The paradigm for Rif1 function is to delay helicase activation within certain regions of the genome, causing these regions to replicate late in S phase. Many questions, however, remain about the intricacies of Rif1 function. Here, we review the current models for the activity of Rif1 with the goal of trying to understand how Rif1 functions to establish the RT program.

Super-resolution microscopy reveals stochastic initiation of replication in Drosophila polytene chromosomes

Studying the probability distribution of replication initiation along a chromosome is a huge challenge. Drosophila polytene chromosomes in combination with super-resolution microscopy provide a unique opportunity for analyzing the probabilistic nature of replication initiation at the ultrastructural level. Here, we developed a method for synchronizing S-phase induction among salivary gland cells. An analysis of the replication label distribution in the first minutes of S phase and in the following hours after the induction revealed the dynamics of replication initiation. Spatial super-resolution structured illumination microscopy allowed identifying multiple discrete replication signals and to investigate the behavior of replication signals in the first minutes of the S phase at the ultrastructural level. We identified replication initiation zones where initiation occurs stochastically. These zones differ significantly in the probability of replication initiation per time unit. There are zones in which initiation occurs on most strands of the polytene chromosome in a few minutes. In other zones, the initiation on all strands takes several hours. Compact bands are free of replication initiation events, and the replication runs from outer edges to the middle, where band shapes may alter.

Stochastic initiation of DNA replication across the human genome

Wang et al. (2021) comprehensively map DNA replication initiation events across the human genome using single-molecule optical resolution mapping and find that initiation events are randomly distributed across broad initiation zones that are only utilized in a stochastic fashion across a population of cells.

DNA molecular combing-based replication fork directionality profiling

The replication strategy of metazoan genomes is still unclear, mainly because definitive maps of replication origins are missing. High-throughput methods are based on population average and thus may exclusively identify efficient initiation sites, whereas inefficient origins go undetected. Single-molecule analyses of specific loci can detect both common and rare initiation events along the targeted regions. However, these usually concentrate on positioning individual events, which only gives an overview of the replication dynamics. Here, we computed the replication fork directionality (RFD) profiles of two large genes in different transcriptional states in chicken DT40 cells, namely untranscribed and transcribed DMD and CCSER1 expressed at WT levels or overexpressed, by aggregating hundreds of oriented replication tracks detected on individual DNA fibres stretched by molecular combing. These profiles reconstituted RFD domains composed of zones of initiation flanking a zone of termination originally observed in mammalian genomes and were highly consistent with independent population-averaging profiles generated by Okazaki fragment sequencing. Importantly, we demonstrate that inefficient origins do not appear as detectable RFD shifts, explaining why dispersed initiation has remained invisible to population-based assays. Our method can both generate quantitative profiles and identify discrete events, thereby constituting a comprehensive approach to study metazoan genome replication.

Mammalian DNA Replication Timing

Immediately following the discovery of the structure of DNA and the semi-conservative replication of the parental DNA sequence into two new DNA strands, it became apparent that DNA replication is organized in a temporal and spatial fashion during the S phase of the cell cycle, correlated with the large-scale organization of chromatin in the nucleus. After many decades of limited progress, technological advances in genomics, genome engineering, and imaging have finally positioned the field to tackle mechanisms underpinning the temporal and spatial regulation of DNA replication and the causal relationships between DNA replication and other features of large-scale chromosome structure and function. In this review, we discuss these major recent discoveries as well as expectations for the coming decade.

Genome-wide mapping of human DNA replication by optical replication mapping supports a stochastic model of eukaryotic replication

The heterogeneous nature of eukaryotic replication kinetics and the low efficiency of individual initiation sites make mapping the location and timing of replication initiation in human cells difficult. To address this challenge, we have developed optical replication mapping (ORM), a high-throughput single-molecule approach, and used it to map early-initiation events in human cells. The single-molecule nature of our data and a total of >2,500-fold coverage of the human genome on 27 million fibers averaging ∼300 kb in length allow us to identify initiation sites and their firing probability with high confidence. We find that the distribution of human replication initiation is consistent with inefficient, stochastic activation of heterogeneously distributed potential initiation complexes enriched in accessible chromatin. These observations are consistent with stochastic models of initiation-timing regulation and suggest that stochastic regulation of replication kinetics is a fundamental feature of eukaryotic replication, conserved from yeast to humans.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

In silico analysis of DNA re-replication across a complete genome reveals cell-to-cell heterogeneity and genome plasticity

DNA replication is a complex and remarkably robust process: despite its inherent uncertainty, manifested through stochastic replication timing at a single-cell level, multiple control mechanisms ensure its accurate and timely completion across a population. Disruptions in these mechanisms lead to DNA re-replication, closely connected to genomic instability and oncogenesis. Here, we present a stochastic hybrid model of DNA re-replication that accurately portrays the interplay between discrete dynamics, continuous dynamics and uncertainty. Using experimental data on the fission yeast genome, model simulations show how different regions respond to re-replication and permit insight into the key mechanisms affecting re-replication dynamics. Simulated and experimental population-level profiles exhibit a good correlation along the genome, robust to model parameters, validating our approach. At a single-cell level, copy numbers of individual loci are affected by intrinsic properties of each locus, in cis effects from adjoining loci and in trans effects from distant loci. In silico analysis and single-cell imaging reveal that cell-to-cell heterogeneity is inherent in re-replication and can lead to genome plasticity and a plethora of genotypic variations.

Chromatin and Nuclear Architecture: Shaping DNA Replication in 3D

In eukaryotes, DNA replication progresses through a finely orchestrated temporal and spatial program. The 3D genome structure and nuclear architecture have recently emerged as fundamental determinants of the replication program. Factors with established roles in replication have been recognized as genome organization regulators. Exploiting paradigms from yeasts and mammals, we discuss how DNA replication is regulated in time and space through DNA-associated trans-acting factors, diffusible limiting replication initiation factors, higher-order chromatin folding, dynamic origin localization, and specific nuclear microenvironments. We present an integrated model for the regulation of DNA replication in 3D and highlight the importance of accurate spatio-temporal regulation of DNA replication in physiology and disease.

Preimplantation Genetic Testing for Chromosomal Abnormalities: Aneuploidy, Mosaicism, and Structural Rearrangements

There is a high incidence of chromosomal abnormalities in early human embryos, whether they are generated by natural conception or by assisted reproductive technologies (ART). Cells with chromosomal copy number deviations or chromosome structural rearrangements can compromise the viability of embryos; much of the naturally low human fecundity as well as low success rates of ART can be ascribed to these cytogenetic defects. Chromosomal anomalies are also responsible for a large proportion of miscarriages and congenital disorders. There is therefore tremendous value in methods that identify embryos containing chromosomal abnormalities before intrauterine transfer to a patient being treated for infertility—the goal being the exclusion of affected embryos in order to improve clinical outcomes. This is the rationale behind preimplantation genetic testing for aneuploidy (PGT-A) and structural rearrangements (-SR). Contemporary methods are capable of much more than detecting whole chromosome abnormalities (e.g., monosomy/trisomy). Technical enhancements and increased resolution and sensitivity permit the identification of chromosomal mosaicism (embryos containing a mix of normal and abnormal cells), as well as the detection of sub-chromosomal abnormalities such as segmental deletions and duplications. Earlier approaches to screening for chromosomal abnormalities yielded a binary result of normal versus abnormal, but the new refinements in the system call for new categories, each with specific clinical outcomes and nuances for clinical management. This review intends to give an overview of PGT-A and -SR, emphasizing recent advances and areas of active development.

FORK-seq: replication landscape of the Saccharomyces cerevisiae genome by nanopore sequencing

Genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200-nucleotide resolution. By quantifying BrdU incorporation along pulse-chased replication intermediates from Saccharomyces cerevisiae, we orient 58,651 replication tracks reproducing population-based replication directionality profiles and map 4964 and 4485 individual initiation and termination events, respectively. Although most events cluster at known origins and fork merging zones, 9% and 18% of initiation and termination events, respectively, occur at many locations previously missed. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

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.

Capturing the dynamics of genome replication on individual ultra-long nanopore sequence reads

Replication of eukaryotic genomes is highly stochastic, making it difficult to determine the replication dynamics of individual molecules with existing methods. We report a sequencing method for the measurement of replication fork movement on single molecules by detecting nucleotide analog signal currents on extremely long nanopore traces (D-NAscent). Using this method, we detect 5-bromodeoxyuridine (BrdU) incorporated by Saccharomyces cerevisiae to reveal, at a genomic scale and on single molecules, the DNA sequences replicated during a pulse-labeling period. Under conditions of limiting BrdU concentration, D-NAscent detects the differences in BrdU incorporation frequency across individual molecules to reveal the location of active replication origins, fork direction, termination sites, and fork pausing/stalling events. We used sequencing reads of 20-160 kilobases to generate a whole-genome single-molecule map of DNA replication dynamics and discover a class of low-frequency stochastic origins in budding yeast. The D-NAscent software is available at https://github.com/MBoemo/DNAscent.git .

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.

Mutating for Good: DNA Damage Responses During Somatic Hypermutation

Somatic hypermutation (SHM) of immunoglobulin (Ig) genes plays a key role in antibody mediated immunity. SHM in B cells provides the molecular basis for affinity maturation of antibodies. In this way SHM is key in optimizing antibody dependent immune responses. SHM is initiated by targeting the Activation-Induced Cytidine Deaminase (AID) to rearranged V(D)J and switch regions of Ig genes. The mutation rate of this programmed mutagenesis is ~10-3 base pairs per generation, a million-fold higher than the non-AID targeted genome of B cells. AID is a processive enzyme that binds single-stranded DNA and deaminates cytosines in DNA. Cytosine deamination generates highly mutagenic deoxy-uracil (U) in the DNA of both strands of the Ig loci. Mutagenic processing of the U by the DNA damage response generates the entire spectrum of base substitutions characterizing SHM at and around the initial U lesion. Starting from the U as a primary lesion, currently five mutagenic DNA damage response pathways have been identified in generating a well-defined SHM spectrum of C/G transitions, C/G transversions, and A/T mutations around this initial lesion. These pathways include (1) replication opposite template U generates transitions at C/G, (2) UNG2-dependent translesion synthesis (TLS) generates transversions at C/G, (3) a hybrid pathway comprising non-canonical mismatch repair (ncMMR) and UNG2-dependent TLS generates transversions at C/G, (4) ncMMR generates mutations at A/T, and (5) UNG2- and PCNA Ubiquitination (PCNA-Ub)-dependent mutations at A/T. Furthermore, specific strand-biases of SHM spectra arise as a consequence of a biased AID targeting, ncMMR, and anti-mutagenic repriming. Here, we review mammalian SHM with special focus on the mutagenic DNA damage response pathways involved in processing AID induced Us, the origin of characteristic strand biases, and relevance of the cell cycle.

hDNA2 nuclease/helicase promotes centromeric DNA replication and genome stability

DNA2 is a nuclease/helicase that is involved in Okazaki fragment maturation, replication fork processing, and end resection of DNA double‐strand breaks. Similar such helicase activity for resolving secondary structures and structure‐specific nuclease activity are needed during DNA replication to process the chromosome‐specific higher order repeat units present in the centromeres of human chromosomes. Here, we show that DNA2 binds preferentially to centromeric DNA. The nuclease and helicase activities of DNA2 are both essential for resolution of DNA structural obstacles to facilitate DNA replication fork movement. Loss of DNA2‐mediated clean‐up mechanisms impairs centromeric DNA replication and CENP‐A deposition, leading to activation of the ATR DNA damage checkpoints at centromeric DNA regions and late‐S/G2 cell cycle arrest. Cells that escape arrest show impaired metaphase plate formation and abnormal chromosomal segregation. Furthermore, the DNA2 inhibitor C5 mimics DNA2 knockout and synergistically kills cancer cells when combined with an ATR inhibitor. These findings provide mechanistic insights into how DNA2 supports replication of centromeric DNA and give further insights into new therapeutic strategies.

DNA replication origins-where do we begin?

In this review, Prioleau and MacAlpine summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

For more than three decades, investigators have sought to identify the precise locations where DNA replication initiates in mammalian genomes. The development of molecular and biochemical approaches to identify start sites of DNA replication (origins) based on the presence of defining and characteristic replication intermediates at specific loci led to the identification of only a handful of mammalian replication origins. The limited number of identified origins prevented a comprehensive and exhaustive search for conserved genomic features that were capable of specifying origins of DNA replication. More recently, the adaptation of origin-mapping assays to genome-wide approaches has led to the identification of tens of thousands of replication origins throughout mammalian genomes, providing an unprecedented opportunity to identify both genetic and epigenetic features that define and regulate their distribution and utilization. Here we summarize recent advances in our understanding of how primary sequence, chromatin environment, and nuclear architecture contribute to the dynamic selection and activation of replication origins across diverse cell types and developmental stages.

Order from clutter: selective interactions at mammalian replication origins

Mammalian chromosome duplication progresses in a precise order and is subject to constraints that are often relaxed in developmental disorders and malignancies. Molecular information about the regulation of DNA replication at the chromatin level is lacking because protein complexes that initiate replication seem to bind chromatin indiscriminately. High-throughput sequencing and mathematical modelling have yielded detailed genome-wide replication initiation maps. Combining these maps and models with functional genetic analyses suggests that distinct DNA-protein interactions at subgroups of replication initiation sites (replication origins) modulate the ubiquitous replication machinery and supports an emerging model that delineates how indiscriminate DNA-binding patterns translate into a consistent, organized replication programme.

Insights into specificity, redundancy and new cellular functions of C/EBPa and C/EBPb transcription factors through interactome network analysis

BACKGROUND

C/EBPa and C/EBPb are transcription factors with tissue specific expression regulating several important cellular processes. They work by recruiting protein complexes to a common DNA recognition motif and both are able to compensate each other's absence in many cell types, thus showing functional redundancy. They also play distinct roles in specific cellular pathways and their abnormal functioning gives raise to different human pathologies.

METHODS

To investigate the molecular basis of C/EBPa and C/EBPb specificity and redundancy we characterized their in vivo protein-protein interaction networks by Tandem Affinity Purification (TAP) and Mass Spectrometry (MS). To unravel the functional features of C/EBPa and C/EBPb proteomes we studied the statistical enrichment of binding partners related to Gene Ontology (GO) terms and KEGG pathways.

RESULTS

Our data confirmed that the C/EBPa and C/EBPb regulate biological processes like cell proliferation, apoptosis and transformation. We found that both C/EBPa and C/EBPb are involved in other cellular pathways such as RNA maturation, RNA splicing and DNA repair. Specific interactions of C/EBPa with MRE11, RUVBL1 and RUVBL2 components of DNA repair system were confirmed by co-immunoprecipitation assays.

CONCLUSIONS

Our comparative analysis of the C/EBPa and C/EBPb proteomes provides an insight for understanding both their redundant and specific roles in cells indicating their involvement in new pathways. Such novel predicted functions are relevant to normal cellular processes and disease phenotypes controlled by these transcription factors.

GENERAL SIGNIFICANCE

Functional characterization of C/EBPa and C/EBPb proteomes suggests they can regulate novel pathways and indicate potential molecular targets for therapeutic intervention.

PrimPol prevents APOBEC/AID family mediated DNA mutagenesis

PrimPol is a DNA damage tolerant polymerase displaying both translesion synthesis (TLS) and (re)-priming properties. This led us to study the consequences of a PrimPol deficiency in tolerating mutagenic lesions induced by members of the APOBEC/AID family of cytosine deaminases. Interestingly, during somatic hypermutation, PrimPol counteracts the generation of C>G transversions on the leading strand. Independently, mutation analyses in human invasive breast cancer confirmed a pro-mutagenic activity of APOBEC3B and revealed a genome-wide anti-mutagenic activity of PRIMPOL as well as most Y-family TLS polymerases. PRIMPOL especially prevents APOBEC3B targeted cytosine mutations within TpC dinucleotides. As C transversions induced by APOBEC/AID family members depend on the formation of AP-sites, we propose that PrimPol reprimes preferentially downstream of AP-sites on the leading strand, to prohibit error-prone TLS and simultaneously stimulate error-free homology directed repair. These in vivo studies are the first demonstrating a critical anti-mutagenic activity of PrimPol in genome maintenance.

Replication landscape of the human genome

Despite intense investigation, human replication origins and termini remain elusive. Existing data have shown strong discrepancies. Here we sequenced highly purified Okazaki fragments from two cell types and, for the first time, quantitated replication fork directionality and delineated initiation and termination zones genome-wide. Replication initiates stochastically, primarily within non-transcribed, broad (up to 150 kb) zones that often abut transcribed genes, and terminates dispersively between them. Replication fork progression is significantly co-oriented with the transcription. Initiation and termination zones are frequently contiguous, sometimes separated by regions of unidirectional replication. Initiation zones are enriched in open chromatin and enhancer marks, even when not flanked by genes, and often border ‘topologically associating domains' (TADs). Initiation zones are enriched in origin recognition complex (ORC)-binding sites and better align to origins previously mapped using bubble-trap than λ-exonuclease. This novel panorama of replication reveals how chromatin and transcription modulate the initiation process to create cell-type-specific replication programs.

The physical origin and termination sites of DNA replication in human cells have remained elusive. Here the authors use Okazaki fragment sequencing to reveal global replication patterns and show how chromatin and transcription modulate the process.

The DNA repair endonuclease Mus81 facilitates fast DNA replication in the absence of exogenous damage

The Mus81 endonuclease resolves recombination intermediates and mediates cellular responses to exogenous replicative stress. Here, we show that Mus81 also regulates the rate of DNA replication during normal growth by promoting replication fork progression while reducing the frequency of replication initiation events. In the absence of Mus81 endonuclease activity, DNA synthesis is slowed and replication initiation events are more frequent. In addition, Mus81-deficient cells fail to recover from exposure to low doses of replication inhibitors and cell viability is dependent on the XPF endonuclease. Despite an increase in replication initiation frequency, cells lacking Mus81 use the same pool of replication origins as Mus81-expressing cells. Therefore, decelerated DNA replication in Mus81-deficient cells does not initiate from cryptic or latent origins not used during normal growth. These results indicate that Mus81 plays a key role in determining the rate of DNA replication without activating a novel group of replication origins.

The spatial and temporal organization of origin firing during the S-phase of fission yeast

Eukaryotes duplicate their genomes using multiple replication origins, but the organization of origin firing along chromosomes and during S-phase is not well understood. Using fission yeast, we report the first genome-wide analysis of the spatial and temporal organization of replication origin firing, analyzed using single DNA molecules that can approach the full length of chromosomes. At S-phase onset, origins fire randomly and sparsely throughout the chromosomes. Later in S-phase, clusters of fired origins appear embedded in the sparser regions, which form the basis of nuclear replication foci. The formation of clusters requires proper histone methylation and acetylation, and their locations are not inherited between cell cycles. The rate of origin firing increases gradually, peaking just before mid S-phase. Toward the end of S-phase, nearly all the available origins within the unreplicated regions are fired, contributing to the timely completion of genome replication. We propose that the majority of origins do not fire as a part of a deterministic program. Instead, origin firing, both individually and as clusters, should be viewed as being mostly stochastic.

Peaks cloaked in the mist: the landscape of mammalian replication origins

Replication of mammalian genomes starts at sites termed replication origins, which historically have been difficult to locate as a result of large genome sizes, limited power of genetic identification schemes, and rareness and fragility of initiation intermediates. However, origins are now mapped by the thousands using microarrays and sequencing techniques. Independent studies show modest concordance, suggesting that mammalian origins can form at any DNA sequence but are suppressed by read-through transcription or that they can overlap the 5' end or even the entire gene. These results require a critical reevaluation of whether origins form at specific DNA elements and/or epigenetic signals or require no such determinants.

Best practices for mapping replication origins in eukaryotic chromosomes

Understanding the regulatory principles ensuring complete DNA replication in each cell division is critical for deciphering the mechanisms that maintain genomic stability. Recent advances in genome sequencing technology facilitated complete mapping of DNA replication sites and helped move the field from observing replication patterns at a handful of single loci to analyzing replication patterns genome-wide. These advances address issues, such as the relationship between replication initiation events, transcription, and chromatin modifications, and identify potential replication origin consensus sequences. This unit summarizes the technological and fundamental aspects of replication profiling and briefly discusses novel insights emerging from mining large datasets, published in the last 3 years, and also describes DNA replication dynamics on a whole-genome scale.

Chromatin structure and replication origins: determinants of chromosome replication and nuclear organization

The DNA replication program is, in part, determined by the epigenetic landscape that governs local chromosome architecture and directs chromosome duplication. Replication must coordinate with other biochemical processes occurring concomitantly on chromatin, such as transcription and remodeling, to insure accurate duplication of both genetic and epigenetic features and to preserve genomic stability. The importance of genome architecture and chromatin looping in coordinating cellular processes on chromatin is illustrated by two recent sets of discoveries. First, chromatin-associated proteins that are not part of the core replication machinery were shown to affect the timing of DNA replication. These chromatin-associated proteins could be working in concert, or perhaps in competition, with the transcriptional machinery and with chromatin modifiers to determine the spatial and temporal organization of replication initiation events. Second, epigenetic interactions are mediated by DNA sequences that determine chromosomal replication. In this review, we summarize recent findings and current models linking spatial and temporal regulation of the replication program with epigenetic signaling. We discuss these issues in the context of the genome's three-dimensional structure with an emphasis on events occurring during the initiation of DNA replication.

Allele-specific genome-wide profiling in human primary erythroblasts reveal replication program organization

We have developed a new approach to characterize allele-specific timing of DNA replication genome-wide in human primary basophilic erythroblasts. We show that the two chromosome homologs replicate at the same time in about 88% of the genome and that large structural variants are preferentially associated with asynchronous replication. We identified about 600 megabase-sized asynchronously replicated domains in two tested individuals. The longest asynchronously replicated domains are enriched in imprinted genes suggesting that structural variants and parental imprinting are two causes of replication asynchrony in the human genome. Biased chromosome X inactivation in one of the two individuals tested was another source of detectable replication asynchrony. Analysis of high-resolution TimEX profiles revealed small variations termed timing ripples, which were undetected in previous, lower resolution analyses. Timing ripples reflect highly reproducible, variations of the timing of replication in the 100 kb-range that exist within the well-characterized megabase-sized replication timing domains. These ripples correspond to clusters of origins of replication that we detected using novel nascent strands DNA profiling methods. Analysis of the distribution of replication origins revealed dramatic differences in initiation of replication frequencies during S phase and a strong association, in both synchronous and asynchronous regions, between origins of replication and three genomic features: G-quadruplexes, CpG Islands and transcription start sites. The frequency of initiation in asynchronous regions was similar in the two homologs. Asynchronous regions were richer in origins of replication than synchronous regions.

DNA replication in mammalian cells proceeds according to a distinct order. Genes that are expressed tend to replicate before genes that are not expressed. We report here that we have developed a method to measure the timing of replication of the maternal and paternal chromosomes separately. We found that the paternal and maternal chromosomes replicate at exactly the same time in the large majority of the genome and that the 12% of the genome that replicated asynchronously was enriched in imprinted genes and in structural variants. Previous experiments have shown that chromosomes could be divided into replication timing domains that are a few hundred thousand to a few megabases in size. We show here that these domains can be divided into sub-domains defined by ripples in the timing profile. These ripples corresponded to clusters of origins of replication. Finally, we show that the frequency of initiation in asynchronous regions was similar in the two homologs.

Epigenetic Regulation of Individual Modules of the immunoglobulin heavy chain locus 3' Regulatory Region

The Igh locus undergoes an amazing array of DNA rearrangements and modifications during B cell development. During early stages, the variable region gene is constructed from constituent variable (V), diversity (D), and joining (J) segments (VDJ joining). B cells that successfully express an antibody can be activated, leading to somatic hypermutation (SHM) focused on the variable region, and class switch recombination (CSR), which substitutes downstream constant region genes for the originally used Cμ constant region gene. Many investigators, ourselves included, have sought to understand how these processes specifically target the Igh locus and avoid other loci and potential deleterious consequences of malignant transformation. Our laboratory has concentrated on a complex regulatory region (RR) that is located downstream of Cα, the most 3′ of the Igh constant region genes. The ~40 kb 3′ RR, which is predicted to serve as a downstream major regulator of the Igh locus, contains two distinct segments: an ~28 kb region comprising four enhancers, and an adjacent ~12 kb region containing multiple CTCF and Pax5 binding sites. Analysis of targeted mutations in mice by a number of investigators has concluded that the entire 3′ RR enhancer region is essential for SHM and CSR (but not for VDJ joining) and for high levels of expression of multiple isotypes. The CTCF/Pax5 binding region is a candidate for influencing VDJ joining early in B cell development and serving as a potential insulator of the Igh locus. Components of the 3′ RR are subject to a variety of epigenetic changes during B cell development, i.e., DNAse I hypersensitivity, histone modifications, and DNA methylation, in association with transcription factor binding. I propose that these changes provide a foundation by which regulatory elements in modules of the 3′ RR function by interacting with each other and with target sequences of the Igh locus.

A new light on DNA replication from the inactive X chromosome

While large portions of the mammalian genome are known to replicate sequentially in a distinct, tissue-specific order, recent studies suggest that the inactive X chromosome is duplicated rapidly via random, synchronous DNA synthesis at numerous adjacent regions. The rapid duplication of the inactive X chromosome was observed in high-resolution studies visualizing DNA replication patterns in the nucleus, and by allele-specific DNA sequencing studies measuring the extent of DNA synthesis. These studies conclude that inactive X chromosomes complete replication earlier than previously thought and suggest that the strict order of DNA replication detected in the majority of genomic regions is not preserved in non-transcribed, "silent" chromatin. These observations alter current concepts about the regulation of DNA replication in non-transcribed portions of the genome in general and in the inactive X-chromosome in particular.

Inferring the spatiotemporal DNA replication program from noisy data

We generalize a stochastic model of DNA replication to the case where replication-origin-initiation rates vary locally along the genome and with time. Using this generalized model, we address the inverse problem of inferring initiation rates from experimental data concerning replication in cell populations. Previous work based on curve fitting depended on arbitrarily chosen functional forms for the initiation rate, with free parameters that were constrained by the data. We introduce a nonparametric method of inference that is based on Gaussian process regression. The method replaces specific assumptions about the functional form of the initiation rate with more general prior expectations about the smoothness of variation of this rate, along the genome and in time. Using this inference method, we recover, with high precision, simulated replication schemes from noisy data that are typical of current experiments.

Dissection of the beta-globin replication-initiation region reveals specific requirements for replicator elements during gene amplification

Gene amplification plays a pivotal role in malignant transformation of human cells. A plasmid with both a mammalian replication-initiation region (IR)/origin/replicator and a nuclear matrix-attachment region (MAR) is spontaneously amplified in transfected cells by a mechanism that involves amplification at the extrachromosomal site, followed by amplification at the chromosomal arm, ultimately generating a long homogeneously staining region (HSR). Several observations suggest that replication initiation from IR sequences might mediate amplification. To test this idea, we previously dissected c-myc and DHFR IRs to identify the minimum sequence required to support amplification. In this study, we applied an improved analysis that discriminates between two amplification steps to the ß-globin RepP IR, which contains separate elements already known to be essential for initiation on the chromosome arm. The IR sequence was required at least for the extrachromosomal amplification step. In addition to the vector-encoded MAR, amplification also required an AT-rich region and a MAR-like element, consistent with the results regarding replicator activity on the chromosome. However, amplification did not require the AG-rich tract necessary for replicator activity, but instead required a novel sequence containing another AG-rich tract. The differential sequence requirement might be a consequence of extrachromosomal replication.

DNA replication timing

Patterns of replication within eukaryotic genomes correlate with gene expression, chromatin structure, and genome evolution. Recent advances in genome-scale mapping of replication kinetics have allowed these correlations to be explored in many species, cell types, and growth conditions, and these large data sets have allowed quantitative and computational analyses. One striking new correlation to emerge from these analyses is between replication timing and the three-dimensional structure of chromosomes. This correlation, which is significantly stronger than with any single histone modification or chromosome-binding protein, suggests that replication timing is controlled at the level of chromosomal domains. This conclusion dovetails with parallel work on the heterogeneity of origin firing and the competition between origins for limiting activators to suggest a model in which the stochastic probability of individual origin firing is modulated by chromosomal domain structure to produce patterns of replication. Whether these patterns have inherent biological functions or simply reflect higher-order genome structure is an open question.

The replication domain model: regulating replicon firing in the context of large-scale chromosome architecture

The "Replicon Theory" of Jacob, Brenner, and Cuzin has reliably served as the paradigm for regulating the sites where individual replicons initiate replication. Concurrent with the replicon model was Taylor's demonstration that plant and animal chromosomes replicate segmentally in a defined temporal sequence, via cytologically defined units too large to be accounted for by a single replicon. Instead, there seemed to be a program to choreograph when chromosome units replicate during S phase, executed by initiation at clusters of individual replicons within each segment. Here, we summarize recent molecular evidence for the existence of such units, now known as "replication domains", and discuss how the organization of large chromosomes into structural units has added additional layers of regulation to the original replicon model.

Genetic and epigenetic determinants of DNA replication origins, position and activation

In the genome of eukaryotic cells, DNA synthesis is initiated at multiple sites called origins of DNA replication. Origins must fire only once per cell cycle and how this is achieved is now well understood. However, little is known about the mechanisms that determine when and where replication initiates in a given cell. A large body of evidence indicates that origins are not equal in terms of efficiency and timing of activation. Origin usage also changes concomitantly with the different cell differentiation programs. As DNA replication occurs in the context of chromatin, initiation could be influenced by multiple parameters, such as nucleosome positioning, histone modifications, and three-dimensional (3D) organization of the nucleus. This view is supported by recent genome-wide studies showing that DNA replication profiles are shaped by genetic and epigenetic processes that act both at the local and global levels to regulate origin function in eukaryotic cells.