The Chinese hamster dihydrofolate reductase replication origin decision point follows activation of transcription and suppresses initiation of replication within transcription units

The double life of mammalian DNA replication origins

Mammalian DNA replication origins have been historically difficult to identify and their determinants are still unresolved. Here, we first review methods developed over the last decades to map replication initiation sites either directly via initiation intermediates or indirectly via determining replication fork directionality profiles. We also discuss the factors that may specify these sites as replication initiation sites. Second, we address the controversy that has emerged from these results over whether origins are narrowly defined and localized to specific sites or are more dispersed and organized into broad zones. Ample evidence in favor of both scenarios currently creates an impression of unresolved confusion in the field. We attempt to formulate a synthesis of both models and to reconcile discrepant findings. It is evident that not only one approach is sufficient in isolation but that the combination of several is instrumental toward understanding initiation sites in mammalian genomes. We argue that an aggregation of several individual and often inefficient initiation sites into larger initiation zones and the existence of efficient unidirectional initiation sites and fork stalling at the borders of initiation zones can reconcile the different observations.

Master transcription factor binding sites constitute the core of early replication control elements

Eukaryotic genomes replicate in a defined temporal order called the replication timing (RT) program. RT is developmentally regulated with potential to drive cell fate transitions, but mechanisms controlling RT remain elusive. We previously identified "Early Replication Control Elements" (ERCEs) necessary for early RT, domain-wide transcription, 3D chromatin architecture and compartmentalization in mouse embryonic stem cells (mESCs) but, deletions identifying ERCEs were large and encompassed many putative regulatory elements. Here, we show that ERCEs are compound elements whose RT activity can largely be accounted for by multiple sites of diverse master transcription factor binding (subERCEs), distinguished from other such sites by their long-range interactions. While deletion of subERCEs had large effects on both transcription and RT, deleting transcription start sites eliminated nearly all transcription with moderate effects on RT. Our results suggest a model in which subERCEs respond to diverse master transcription factors by functioning both as transcription enhancers and as elements that organize chromatin domains structurally and support early RT, potentially providing a feed-forward loop to drive robust epigenomic change during cell fate transitions.

Transcription can be sufficient, but is not necessary, to advance replication timing

DNA replication timing (RT) is correlated with transcription during cell fate changes but there are many exceptions and our understanding of this relationship suffers from a paucity of reductionist approaches. Here, we manipulated length and strength of transcription in hybrid-genome mouse embryonic stem cells (mESCs) at a single locus upstream of the silent, late replicating, Pleiotrophin (Ptn) gene, directly comparing RT to nascent transcription rates at engineered vs. wild-type alleles. First, we inserted four reporter genes that differ only in their promoter. Two promoters transcribed the reporter gene at high rates and advanced RT. The other two transcribed at lower rates and did not advance RT. Since these promoters may prove useful in applications where effects on RT are undesirable, we confirmed the inability of one of them to advance RT at numerous ectopic sites. We next juxtaposed these same four promoters upstream of the Ptn transcription start site where they all transcribed the 96kb Ptn gene and advanced RT to different extents correlated with transcription rates. Indeed, a doxycycline-responsive promoter, which could not advance RT when induced as a small reporter gene, elicited a rapid and reversible RT advance proportional to the rate of transcription, providing direct evidence that transcription itself can advance RT. However, deletion of the Ptn promoter and enhancer, followed by directed differentiation to neural precursors, eliminated induction of transcription throughout the entire Ptn replication domain, without preventing the switch to early replication. Our results provide a solid empirical base with which to re-evaluate many decades of literature, demonstrating that length and strength of transcription is sufficient but not necessary to advance RT. Our results also provide a robust system in which to rapidly effect an RT change, permitting mechanistic studies of the role of transcription in RT and the consequences of RT changes to epigenomic remodeling.

Human DNA topoisomerase I poisoning causes R loop-mediated genome instability attenuated by transcription factor IIS

DNA topoisomerase I can contribute to cancer genome instability. During catalytic activity, topoisomerase I forms a transient intermediate, topoisomerase I-DNA cleavage complex (Top1cc) to allow strand rotation and duplex relaxation, which can lead to elevated levels of DNA-RNA hybrids and micronuclei. To comprehend the underlying mechanisms, we have integrated genomic data of Top1cc-triggered hybrids and DNA double-strand breaks (DSBs) shortly after Top1cc induction, revealing that Top1ccs increase hybrid levels with different mechanisms. DSBs are at highly transcribed genes in early replicating initiation zones and overlap with hybrids downstream of accumulated RNA polymerase II (RNAPII) at gene 5'-ends. A transcription factor IIS mutant impairing transcription elongation further increased RNAPII accumulation likely due to backtracking. Moreover, Top1ccs can trigger micronuclei when occurring during late G1 or early/mid S, but not during late S. As micronuclei and transcription-replication conflicts are attenuated by transcription factor IIS, our results support a role of RNAPII arrest in Top1cc-induced transcription-replication conflicts leading to DSBs and micronuclei.

Origins of DNA replication in eukaryotes

Errors occurring during DNA replication can result in inaccurate replication, incomplete replication, or re-replication, resulting in genome instability that can lead to diseases such as cancer or disorders such as autism. A great deal of progress has been made toward understanding the entire process of DNA replication in eukaryotes, including the mechanism of initiation and its control. This review focuses on the current understanding of how the origin recognition complex (ORC) contributes to determining the location of replication initiation in the multiple chromosomes within eukaryotic cells, as well as methods for mapping the location and temporal patterning of DNA replication. Origin specification and configuration vary substantially between eukaryotic species and in some cases co-evolved with gene-silencing mechanisms. We discuss the possibility that centromeres and origins of DNA replication were originally derived from a common element and later separated during evolution.

Cohesin-mediated loop anchors confine the locations of human replication origins

DNA replication occurs through an intricately regulated series of molecular events and is fundamental for genome stability1,2. At present, it is unknown how the locations of replication origins are determined in the human genome. Here we dissect the role of topologically associating domains (TADs)3-6, subTADs7 and loops8 in the positioning of replication initiation zones (IZs). We stratify TADs and subTADs by the presence of corner-dots indicative of loops and the orientation of CTCF motifs. We find that high-efficiency, early replicating IZs localize to boundaries between adjacent corner-dot TADs anchored by high-density arrays of divergently and convergently oriented CTCF motifs. By contrast, low-efficiency IZs localize to weaker dotless boundaries. Following ablation of cohesin-mediated loop extrusion during G1, high-efficiency IZs become diffuse and delocalized at boundaries with complex CTCF motif orientations. Moreover, G1 knockdown of the cohesin unloading factor WAPL results in gained long-range loops and narrowed localization of IZs at the same boundaries. Finally, targeted deletion or insertion of specific boundaries causes local replication timing shifts consistent with IZ loss or gain, respectively. Our data support a model in which cohesin-mediated loop extrusion and stalling at a subset of genetically encoded TAD and subTAD boundaries is an essential determinant of the locations of replication origins in human S phase.

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.

Consequences and Resolution of Transcription-Replication Conflicts

Transcription–replication conflicts occur when the two critical cellular machineries responsible for gene expression and genome duplication collide with each other on the same genomic location. Although both prokaryotic and eukaryotic cells have evolved multiple mechanisms to coordinate these processes on individual chromosomes, it is now clear that conflicts can arise due to aberrant transcription regulation and premature proliferation, leading to DNA replication stress and genomic instability. As both are considered hallmarks of aging and human diseases such as cancer, understanding the cellular consequences of conflicts is of paramount importance. In this article, we summarize our current knowledge on where and when collisions occur and how these encounters affect the genome and chromatin landscape of cells. Finally, we conclude with the different cellular pathways and multiple mechanisms that cells have put in place at conflict sites to ensure the resolution of conflicts and accurate genome duplication.

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.

Control of DNA replication timing in the 3D genome

The 3D organization of mammalian chromatin was described more than 30 years ago by visualizing sites of DNA synthesis at different times during the S phase of the cell cycle. These early cytogenetic studies revealed structurally stable chromosome domains organized into subnuclear compartments. Active-gene-rich domains in the nuclear interior replicate early, whereas more condensed chromatin domains that are largely at the nuclear and nucleolar periphery replicate later. During the past decade, this spatiotemporal DNA replication programme has been mapped along the genome and found to correlate with epigenetic marks, transcriptional activity and features of 3D genome architecture such as chromosome compartments and topologically associated domains. But the causal relationship between these features and DNA replication timing and the regulatory mechanisms involved have remained an enigma. The recent identification of cis-acting elements regulating the replication time and 3D architecture of individual replication domains and of long non-coding RNAs that coordinate whole chromosome replication provide insights into such mechanisms.

Increased growth rate and productivity following stable depletion of miR-7 in a mAb producing CHO cell line causes an increase in proteins associated with the Akt pathway and ribosome biogenesis

Cell line engineering using microRNAs represents a desirable route for improving the efficiency of recombinant protein production by CHO cells. In this study we generated stable CHO DP12 cells expressing a miR-7 sponge transcript which sequesters miR-7 from its endogenous targets. Depletion of miR-7 results in a 65% increase in cell growth and >3-fold increase in yield of secreted IgG protein. Quantitative labelfree LC-MS/MS proteomic profiling was carried out to identify the targets of miR-7 and understand the functional drivers of the improved CHO cell phenotypes. Subcellular enrichment and total proteome analysis identified more than 3000 proteins per fraction resulting in over 5000 unique proteins identified per timepoint analysed. Early stage culture analysis identified 117 proteins overexpressed in miR-7 depleted cells. A subset of these proteins are involved in the Akt pathway which could be the underlying route for cell density improvement and may be exploited more specifically in the future. Late stage culture identified 160 proteins overexpressed in miR-7 depleted cells with some of these involved in ribosome biogenesis which may be causing the increased productivity through improved translational efficiency. This is the first in-depth proteomic profiling of the IgG producing CHO DP12 cell line stably depleted of miR-7. SIGNIFICANCE: Chinese hamster ovary (CHO) cells are the mammalian cell expression system of choice for production of recombinant therapeutic proteins. There is much research ongoing to characterise CHO cell factories through the application of systems biology approaches that will enable a fundamental understanding of CHO cell physiology, and as a result, a better knowledge and understanding of recombinant protein production. This study profiles the proteomic effects of microRNA-7 depletion on the IgG producing CHO DP12 cell line. This is one of the very few studies that attempts to identify the functioning proteins driving improved CHO cell phenotypes resulting from microRNA manipulation. Using subcellular enrichment and total proteome analysis we identified over 5000 unique proteins in miR-7 depleted CHO cells. This work has identified a cohort of proteins involved in the Akt pathway and ribosome biogenesis. These proteins may drive improved CHO cell phenotypes and are of great interest for future work.

Intragenic origins due to short G1 phases underlie oncogene-induced DNA replication stress

Oncogene-induced DNA replication stress contributes critically to the genomic instability that is present in cancer. However, elucidating how oncogenes deregulate DNA replication has been impeded by difficulty in mapping replication initiation sites on the human genome. Here, using a sensitive assay to monitor nascent DNA synthesis in early S phase, we identified thousands of replication initiation sites in cells before and after induction of the oncogenes CCNE1 and MYC. Remarkably, both oncogenes induced firing of a novel set of DNA replication origins that mapped within highly transcribed genes. These ectopic origins were normally suppressed by transcription during G1, but precocious entry into S phase, before all genic regions had been transcribed, allowed firing of origins within genes in cells with activated oncogenes. Forks from oncogene-induced origins were prone to collapse, as a result of conflicts between replication and transcription, and were associated with DNA double-stranded break formation and chromosomal rearrangement breakpoints both in our experimental system and in a large cohort of human cancers. Thus, firing of intragenic origins caused by premature S phase entry represents a mechanism of oncogene-induced DNA replication stress that is relevant for genomic instability in human cancer.

Replicating Large Genomes: Divide and Conquer

Complete duplication of large metazoan chromosomes requires thousands of potential initiation sites, only a small fraction of which are selected in each cell cycle. Assembly of the replication machinery is highly conserved and tightly regulated during the cell cycle, but the sites of initiation are highly flexible, and their temporal order of firing is regulated at the level of large-scale multi-replicon domains. Importantly, the number of replication forks must be quickly adjusted in response to replication stress to prevent genome instability. Here we argue that large genomes are divided into domains for exactly this reason. Once established, domain structure abrogates the need for precise initiation sites and creates a scaffold for the evolution of other chromosome functions.

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.

DNA replication stress as a hallmark of cancer

Human cancers share properties referred to as hallmarks, among which sustained proliferation, escape from apoptosis, and genomic instability are the most pervasive. The sustained proliferation hallmark can be explained by mutations in oncogenes and tumor suppressors that regulate cell growth, whereas the escape from apoptosis hallmark can be explained by mutations in the TP53, ATM, or MDM2 genes. A model to explain the presence of the three hallmarks listed above, as well as the patterns of genomic instability observed in human cancers, proposes that the genes driving cell proliferation induce DNA replication stress, which, in turn, generates genomic instability and selects for escape from apoptosis. Here, we review the data that support this model, as well as the mechanisms by which oncogenes induce replication stress. Further, we argue that DNA replication stress should be considered as a hallmark of cancer because it likely drives cancer development and is very prevalent.

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.

MYC and the control of DNA replication

The MYC oncogene is a multifunctional protein that is aberrantly expressed in a significant fraction of tumors from diverse tissue origins. Because of its multifunctional nature, it has been difficult to delineate the exact contributions of MYC's diverse roles to tumorigenesis. Here, we review the normal role of MYC in regulating DNA replication as well as its ability to generate DNA replication stress when overexpressed. Finally, we discuss the possible mechanisms by which replication stress induced by aberrant MYC expression could contribute to genomic instability and cancer.

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.

Control over DNA replication in time and space

DNA replication is precisely regulated in time and space, thereby safeguarding genomic integrity. In eukaryotes, replication initiates from multiple sites along the genome, termed origins of replication, and propagates bidirectionally. Dynamic origin bound complexes dictate where and when replication should initiate. During late mitosis and G1 phase, putative origins are recognized and become "licensed" through the assembly of pre-replicative complexes (pre-RCs) that include the MCM2-7 helicases. Subsequently, at the G1/S phase transition, a fraction of pre-RCs are activated giving rise to the establishment of replication forks. Origin location is influenced by chromatin and nuclear organization and origin selection exhibits stochastic features. The regulatory mechanisms that govern these cell cycle events rely on the periodic fluctuation of cyclin dependent kinase (CDK) activity through the cell cycle.

Evidence for a mammalian late-G1 phase inhibitor of replication licensing distinct from geminin or Cdk activity

Pre-replication complexes (pre-RCs) are assembled onto DNA during late mitosis and G1 to license replication origins for use in S phase. In order to prevent re-replication of DNA, licensing must be completely shutdown prior to entry into S phase. While mechanisms preventing re-replication during S phase and mitosis have been elucidated, the means by which cells first prevent licensing during late G1 phase are poorly understood. We have employed a hybrid mammalian / Xenopus egg extract replication system to dissect activities that inhibit replication licensing at different stages of the cell cycle in Chinese Hamster Ovary (CHO) cells. We find that soluble extracts from mitotic cells inhibit licensing through a combination of geminin and Cdk activities, while extracts from S-phase cells inhibit licensing predominantly through geminin alone. Surprisingly however, geminin did not accumulate until after cells enter S phase. Unlike extracts from cells in early G1 phase, extracts from late G1 phase and early S phase cells contained an inhibitor of licensing that could not be accounted for by either geminin or Cdk. Moreover, inhibiting cyclin and geminin protein synthesis or inhibiting Cdk activity early in G1 phase did not prevent the appearance of inhibitory activity. These results suggest that a soluble inhibitor of replication licensing appears prior to entry into S phase that is distinct from either geminin or Cdk activity. Our hybrid system should permit the identification of this and other novel cell cycle regulatory activities.

Genome-scale analysis of replication timing: from bench to bioinformatics

Replication timing profiles are cell type-specific and reflect genome organization changes during differentiation. In this protocol, we describe how to analyze genome-wide replication timing (RT) in mammalian cells. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU) and sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated, amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray, which is currently more cost effective than high-throughput sequencing and equally capable of resolving features at the biologically relevant level of tens to hundreds of kilobases. We also present a guide to analyzing the resulting data sets based on methods we use routinely. Subjects include normalization, scaling and data quality measures, LOESS (local polynomial) smoothing of RT values, segmentation of data into domains and assignment of timing values to gene promoters. Finally, we cover clustering methods and means to relate changes in the replication program to gene expression and other genetic and epigenetic data sets. Some experience with R or similar programming languages is assumed. All together, the protocol takes ∼3 weeks per batch of samples.

Highly stable loading of Mcm proteins onto chromatin in living cells requires replication to unload

Components of the minichromosome maintenance complex (Mcm2-7) remain indefinitely bound to chromatin during G1 phase and replication arrest.

The heterohexameric minichromosome maintenance protein complex (Mcm2-7) functions as the eukaryotic helicase during DNA replication. Mcm2-7 loads onto chromatin during early G1 phase but is not converted into an active helicase until much later during S phase. Hence, inactive Mcm complexes are presumed to remain stably bound from early G1 through the completion of S phase. Here, we investigated Mcm protein dynamics in live mammalian cells. We demonstrate that Mcm proteins are irreversibly loaded onto chromatin cumulatively throughout G1 phase, showing no detectable exchange with a gradually diminishing soluble pool. Eviction of Mcm requires replication; during replication arrest, Mcm proteins remained bound indefinitely. Moreover, the density of immobile Mcms is reduced together with chromatin decondensation within sites of active replication, which provides an explanation for the lack of colocalization of Mcm with replication fork proteins. These results provide in vivo evidence for an exceptionally stable lockdown mechanism to retain all loaded Mcm proteins on chromatin throughout prolonged cell cycles.

Pre-replication complex proteins assemble at regions of low nucleosome occupancy within the Chinese hamster dihydrofolate reductase initiation zone

Genome-scale mapping of pre-replication complex proteins has not been reported in mammalian cells. Poor enrichment of these proteins at specific sites may be due to dispersed binding, poor epitope availability or cell cycle stage-specific binding. Here, we have mapped sites of biotin-tagged ORC and MCM protein binding in G1-synchronized populations of Chinese hamster cells harboring amplified copies of the dihydrofolate reductase (DHFR) locus, using avidin-affinity purification of biotinylated chromatin followed by high-density microarray analysis across the DHFR locus. We have identified several sites of significant enrichment for both complexes distributed throughout the previously identified initiation zone. Analysis of the frequency of initiations across stretched DNA fibers from the DHFR locus confirmed a broad zone of de-localized initiation activity surrounding the sites of ORC and MCM enrichment. Mapping positions of mononucleosomal DNA empirically and computing nucleosome-positioning information in silico revealed that ORC and MCM map to regions of low measured and predicted nucleosome occupancy. Our results demonstrate that specific sites of ORC and MCM enrichment can be detected within a mammalian intitiation zone, and suggest that initiation zones may be regions of generally low nucleosome occupancy where flexible nucleosome positioning permits flexible pre-RC assembly sites.

Evaluating genome-scale approaches to eukaryotic DNA replication

Mechanisms regulating where and when eukaryotic DNA replication initiates remain a mystery. Recently, genome-scale methods have been brought to bear on this problem. The identification of replication origins and their associated proteins in yeasts is a well-integrated investigative tool, but corresponding data sets from multicellular organisms are scarce. By contrast, standardized protocols for evaluating replication timing have generated informative data sets for most eukaryotic systems. Here, I summarize the genome-scale methods that are most frequently used to analyse replication in eukaryotes, the kinds of questions each method can address and the technical hurdles that must be overcome to gain a complete understanding of the nature of eukaryotic replication origins.

A winding road to origin discovery

Studies in our laboratory over the last three decades have shown that the Chinese hamster dihydrofolate reductase (DHFR) origin of replication corresponds to a broad zone of inefficient initiation sites distributed throughout the spacer between the convergently transcribed DHFR and 2BE2121 genes. It is clear from mutational analysis that none of these sites is genetically required for controlling origin activity. However, the integrity of the promoter of the DHFR gene is needed to activate the downstream origin, while the 3' processing signals prevent invasion and inactivation of the downstream origin by transcription forks. Several other origins in metazoans have been shown to correspond to zones of inefficient sites, while a different subset appears to be similar to the fixed replicators that characterize origins in S. cerevisiae and lower organisms. These observations have led us to suggest a model in which the mammalian genome is dotted with a hierarchy of degenerate, redundant, and inefficient replicators at intervals of a kilobase or less, some of which may have evolved to be highly circumscribed and efficient. The activities of initiation sites are proposed to be largely regulated by local transcription and chromatin architecture. Recently, we and others have devised strategies for identifying active origins on a genome-wide scale in order to define their distributions between fixed and dispersive origin types and to detect relationships among origins, genes, and epigenetic markers. The global pictures emerging are suggestive but far from complete and appear to be plagued by some of the same uncertainties that have led to conflicting views of individual origins in the past (particularly DHFR). In this paper, we will trace the history of origin discovery in mammalian genomes, primarily using the well-studied DHFR origin as a model, because it has been analyzed by nearly every available origin mapping technique in several different laboratories, while many origins have been identified by only one. We will address the strengths and shortcomings of the various methods utilized to identify and characterize origins in complex genomes and will point out how we and others were sometimes led astray by false assumptions and biases, as well as insufficient information. The goal is to help guide future experiments that will provide a truly comprehensive and accurate portrait of origins and their regulation. After all, in the words of George Santayana, "Those who do not learn from history are doomed to repeat it."

Heterochromatin protein 1 (HP1) modulates replication timing of the Drosophila genome

The replication of a chromosomal region during S phase can be highly dynamic between cell types that differ in transcriptome and epigenome. Early replication timing has been positively correlated with several histone modifications that occur at active genes, while repressive histone modifications mark late replicating regions. This raises the question if chromatin modulates the initiating events of replication. To gain insights into this question, we have studied the function of heterochromatin protein 1 (HP1), which is a reader of repressive methylation at histone H3 lysine 9, in genome-wide organization of replication. Cells with reduced levels of HP1 show an advanced replication timing of centromeric repeats in agreement with the model that repressive chromatin mediates the very late replication of large clusters of constitutive heterochromatin. Surprisingly, however, regions with high levels of interspersed repeats on the chromosomal arms, in particular on chromosome 4 and in pericentromeric regions of chromosome 2, behave differently. Here, loss of HP1 results in delayed replication. The fact that these regions are bound by HP1 suggests a direct effect. Thus while HP1 mediates very late replication of centromeric DNA, it is also required for early replication of euchromatic regions with high levels of repeats. This observation of opposing functions of HP1 suggests a model where HP1-mediated repeat inactivation or replication complex loading on the chromosome arms is required for proper activation of origins of replication that fire early. At the same time, HP1-mediated repression at constitutive heterochromatin is required to ensure replication of centromeric repeats at the end of S phase.

Cytometry of chromatin bound Mcm6 and PCNA identifies two states in G1 that are separated functionally by the G1 restriction point

BACKGROUND

Cytometric measurements of DNA content and chromatin-bound Mcm2 have demonstrated bimodal patterns of expression in G1. These patterns, the replication licensing function of Mcm proteins, and a correlation between Mcm loading and cell cycle commitment for cells re-entering the cell cycle, led us to test the idea that cells expressing a defined high level of chromatin-bound Mcm6 in G1 are committed--i.e., past the G1 restriction point. We developed a cell-based assay for tightly-bound PCNA (PCNA*) and Mcm6 (Mcm6*), DNA content, and a mitotic marker to clearly define G1, S, G2, and M phases of the cell cycle. hTERT-BJ1, hTERT-RPE-1, and Molt4 cells were extracted with Triton X-100 followed by methanol fixation, stained with antibodies and DAPI, then measured by cytometry.

RESULTS

Bivariate analysis of cytometric data demonstrated complex patterns with distinct clustering for all combinations of the 4 variables. In G1, cells clustered in two groups characterized by low and high Mcm6* expression. Serum starvation and release experiments showed that residence in the high group was in late G1, just prior to S phase. Kinetic experiments, employing serum withdrawal, and stathmokinetic analysis with aphidicolin, mimosine or nocodazole demonstrated that cells with high levels of Mcm6* cycled with the committed phases of the cell cycle (S, G2, and M).

CONCLUSIONS

A multivariate assay for Mcm6*, PCNA*, DNA content, and a mitotic marker provides analysis capable of estimating the fraction of pre and post-restriction point G1 cells and supports the idea that there are at least two states in G1 defined by levels of chromatin bound Mcm proteins.

Identification and analysis of specific chromosomal region adjacent to exogenous Dhfr-amplified region in Chinese hamster ovary cell genome

Chinese hamster ovary (CHO) cells are widely used for the stable production of recombinant proteins. Gene amplification techniques are frequently used to improve of protein production, and the dihydrofolate reductase (DHFR) gene amplification system is most widely used in the CHO cell line. We previously constructed a CHO genomic bacterial artificial chromosome (BAC) library from a mouse Dhfr-amplified CHO DR1000L-4N cell line and one BAC clone (Cg0031N14) containing the CHO genomic DNA sequence adjacent to Dhfr was selected. To identify the specific chromosomal region adjacent to the exogenous Dhfr-amplified region in the CHO cell genome, we performed further screening of BAC clones to obtain other Dhfr-amplified regions in the CHO genome. From the screening by high-density replica filter hybridization using a digoxigenin-labeled pSV2-dhfr/hGM-CSF probe, we obtained 8 new BAC clones containing a Dhfr-amplified region. To define the structures of the 8 BAC clones, Southern blot analysis, BAC end sequencing and fluorescence in situ hybridization (FISH) were performed. These results revealed that all the selected BAC clones contained a large palindrome structure with a small inverted repeat in the junction region. This suggests that the obtained amplicon structure in the Dhfr-amplified region in the CHO genome plays an important role in exogenous gene amplification.

[Establishment of spatial and temporal program for mammalian chromosome replication]

It has been 55 years since the elucidation of the structure of DNA, suggesting an elegantly simple means for its self-replication. Who would have dreamed in 1953 that it would take longer for us to understand DNA replication than it would for us to uncover the basic rules of animal development? Without question, the mechanisms regulating where and when DNA replication initiates in the cells of our own body is the greatest remaining fundamental mystery in molecular biology. Cis-acting sequences that function as replication origins in mammalian cells have not been identified and the mechanisms that regulate where and when origins will fire during S-phase remain elusive. Indeed, the problem has been so difficult that most researchers move on to more lucrative fields. In this essay, I will summarize my laboratory's humble attempts to make some progress in this area. In doing so, I hope that I can inspire a few young scientists to breath fresh energy into this challenging field.

Cell cycle regulation of DNA replication

Eukaryotic DNA replication is regulated to ensure all chromosomes replicate once and only once per cell cycle. Replication begins at many origins scattered along each chromosome. Except for budding yeast, origins are not defined DNA sequences and probably are inherited by epigenetic mechanisms. Initiation at origins occurs throughout the S phase according to a temporal program that is important in regulating gene expression during development. Most replication proteins are conserved in evolution in eukaryotes and archaea, but not in bacteria. However, the mechanism of initiation is conserved and consists of origin recognition, assembly of prereplication (pre-RC) initiative complexes, helicase activation, and replisome loading. Cell cycle regulation by protein phosphorylation ensures that pre-RC assembly can only occur in G1 phase, whereas helicase activation and loading can only occur in S phase. Checkpoint regulation maintains high fidelity by stabilizing replication forks and preventing cell cycle progression during replication stress or damage.

Proliferation-dependent and cell cycle regulated transcription of mouse pericentric heterochromatin

Pericentric heterochromatin transcription has been implicated in Schizosaccharomyces pombe heterochromatin assembly and maintenance. However, in mammalian systems, evidence for such transcription is inconsistent. We identify two populations of RNA polymerase II–dependent mouse γ satellite repeat sequence–derived transcripts from pericentric heterochromatin that accumulate at different times during the cell cycle. A small RNA species was synthesized exclusively during mitosis and rapidly eliminated during mitotic exit. A more abundant population of large, heterogeneous transcripts was induced late in G1 phase and their synthesis decreased during mid S phase, which is coincident with pericentric heterochromatin replication. In cells that lack the Suv39h1,2 methyltransferases responsible for H3K9 trimethylation, transcription occurs from more sites but is still cell cycle regulated. Transcription is not detected in quiescent cells and induction during G1 phase is sensitive to serum deprivation or the cyclin-dependent kinase inhibitor roscovatine. We demonstrate that mammalian pericentric heterochromatin transcription is linked to cellular proliferation. Our data also provide an explanation for inconsistencies in the detection of such transcripts in different systems.

Replication in context: dynamic regulation of DNA replication patterns in metazoans

Replication in eukaryotes initiates from discrete genomic regions according to a strict, often tissue-specific temporal programme. However, the locations of initiation events within initiation regions vary, show sequence disparity and are affected by interactions with distal elements. Increasing evidence suggests that specification of replication sites and the timing of replication are dynamic processes that are regulated by tissue-specific and developmental cues, and are responsive to epigenetic modifications. Dynamic specification of replication patterns might serve to prevent or resolve possible spatial and/or temporal conflicts between replication, transcription and chromatin assembly, and facilitate subtle or extensive changes of gene expression during differentiation and development.

Replication origin plasticity, Taylor-made: inhibition vs recruitment of origins under conditions of replication stress

Among his many contributions to the field of chromosome structure and dynamics, J. Herbert Taylor showed that eukaryotic cells have many more potential replication origins than they use, which they can recruit when replication forks are slowed to complete S-phase in a timely fashion. Thirty years later, his findings raise an important but largely overlooked paradox. Although new data have confirmed his results, a larger body of data has revealed that slowing replication forks activates an S-phase checkpoint cascade that inhibits initiation from unfired origins until the stress is relieved. In this paper, in celebration of Taylor's work published in Chromosoma 30 years ago, I draw attention to this paradox and offer some plausible models to explain how replication stress can both inhibit and recruit new origins. I hope that this essay will stimulate further experimentation into the basis of Taylor's original findings.

Genome-wide mapping of ORC and Mcm2p binding sites on tiling arrays and identification of essential ARS consensus sequences in S. cerevisiae

BACKGROUND

Eukaryotic replication origins exhibit different initiation efficiencies and activation times within S-phase. Although local chromatin structure and function influences origin activity, the exact mechanisms remain poorly understood. A key to understanding the exact features of chromatin that impinge on replication origin function is to define the precise locations of the DNA sequences that control origin function. In S. cerevisiae, Autonomously Replicating Sequences (ARSs) contain a consensus sequence (ACS) that binds the Origin Recognition Complex (ORC) and is essential for origin function. However, an ACS is not sufficient for origin function and the majority of ACS matches do not function as ORC binding sites, complicating the specific identification of these sites.

RESULTS

To identify essential origin sequences genome-wide, we utilized a tiled oligonucleotide array (NimbleGen) to map the ORC and Mcm2p binding sites at high resolution. These binding sites define a set of potential Autonomously Replicating Sequences (ARSs), which we term nimARSs. The nimARS set comprises 529 ORC and/or Mcm2p binding sites, which includes 95% of known ARSs, and experimental verification demonstrates that 94% are functional. The resolution of the analysis facilitated identification of potential ACSs (nimACSs) within 370 nimARSs. Cross-validation shows that the nimACS predictions include 58% of known ACSs, and experimental verification indicates that 82% are essential for ARS activity.

CONCLUSION

These findings provide the most comprehensive, accurate, and detailed mapping of ORC binding sites to date, adding to the emerging picture of the chromatin organization of the budding yeast genome.

Discrete functional elements required for initiation activity of the Chinese hamster dihydrofolate reductase origin beta at ectopic chromosomal sites

The Chinese hamster dihydrofolate reductase (DHFR) DNA replication initiation region, the 5.8 kb ori-beta, can function as a DNA replicator at random ectopic chromosomal sites in hamster cells. We report a detailed genetic analysis of the DiNucleotide Repeat (DNR) element, one of several sequence elements necessary for ectopic ori-beta activity. Deletions within ori-beta identified a 132 bp core region within the DNR element, consisting mainly of dinucleotide repeats, and a downstream region that are required for ori-beta initiation activity at non-specific ectopic sites in hamster cells. Replacement of the DNR element with Xenopus or mouse transcriptional elements from rDNA genes restored full levels of initiation activity, but replacement with a nucleosome positioning element or a viral intron sequence did not. The requirement for the DNR element and three other ori-beta sequence elements was conserved when ori-beta activity was tested at either random sites or at a single specific ectopic chromosomal site in human cells. These results confirm the importance of specific cis-acting elements in directing the initiation of DNA replication in mammalian cells, and provide new evidence that transcriptional elements can functionally substitute for one of these elements in ori-beta.

Binding of AlF-C, an Orc1-binding transcriptional regulator, enhances replicator activity of the rat aldolase B origin

A region encompassing the rat aldolase B gene (aldB) promoter acts as a chromosomal origin of DNA replication (origin) in rat aldolase B-nonexpressing hepatoma cells. To examine replicator function of the aldB origin, we constructed recombinant mouse cell lines in which the rat aldB origin and the mutant derivatives were inserted into the same position at the mouse chromosome 8 by cre-mediated recombination. Nascent strand abundance assays revealed that the rat origin acts as a replicator at the ectopic mouse locus. Mutation of site C in the rat origin, which binds an Orc1-binding protein AlF-C in vitro, resulted in a significant reduction of the replicator activity in the mouse cells. Chromatin immunoprecipitation (ChIP) assays indicated that the reduction of replicator activity was paralleled with the reduced binding of AlF-C and Orc1, suggesting that sequence-specific binding of AlF-C to the ectopic rat origin leads to enhanced replicator activity in cooperation with Orc1. Involvement of AlF-C in replication in vivo was further examined for the aldB origin at its original rat locus and for a different rat origin identified in the present study, which contained an AlF-C-binding site. ChIP assays revealed that both replication origins bind AlF-C and Orc1. We think that the results presented here may represent one mode of origin recognition in mammalian cells.