Conservation of epigenetic regulation, ORC binding and developmental timing of DNA replication origins in the genus Drosophila

Finishing the egg

Gamete development is a fundamental process that is highly conserved from early eukaryotes to mammals. As germ cells develop, they must coordinate a dynamic series of cellular processes that support growth, cell specification, patterning, the loading of maternal factors (RNAs, proteins, and nutrients), differentiation of structures to enable fertilization and ensure embryonic survival, and other processes that make a functional oocyte. To achieve these goals, germ cells integrate a complex milieu of environmental and developmental signals to produce fertilizable eggs. Over the past 50 years, Drosophila oogenesis has risen to the forefront as a system to interrogate the sophisticated mechanisms that drive oocyte development. Studies in Drosophila have defined mechanisms in germ cells that control meiosis, protect genome integrity, facilitate mRNA trafficking, and support the maternal loading of nutrients. Work in this system has provided key insights into the mechanisms that establish egg chamber polarity and patterning as well as the mechanisms that drive ovulation and egg activation. Using the power of Drosophila genetics, the field has begun to define the molecular mechanisms that coordinate environmental stresses and nutrient availability with oocyte development. Importantly, the majority of these reproductive mechanisms are highly conserved throughout evolution, and many play critical roles in the development of somatic tissues as well. In this chapter, we summarize the recent progress in several key areas that impact egg chamber development and ovulation. First, we discuss the mechanisms that drive nutrient storage and trafficking during oocyte maturation and vitellogenesis. Second, we examine the processes that regulate follicle cell patterning and how that patterning impacts the construction of the egg shell and the establishment of embryonic polarity. Finally, we examine regulatory factors that control ovulation, egg activation, and successful fertilization.

Anatomy and evolution of a DNA replication origin

DNA amplification occurs at the DNA puff II/9A locus in the fungus fly Sciara coprophila. As a foundation to study the molecular mechanism for the initiating events of II/9A DNA re-replication, we have sequenced 14 kb spanning a DNase hypersensitive site (DHS) upstream of the 1 kb amplification origin and through transcription units II/9-1 and II/9-2 downstream of the origin. These elements are annotated as well as the ORC binding site at the origin and the transition point (TP) between continuous and discontinuous DNA syntheses that marks the origin of bidirectional replication at the nucleotide level. A 9 bp motif found at the TP is repeated near the other end of the 1 kb ORI and may identify a putative second TP. The steroid hormone ecdysone induces DNA amplification as well as transcription and puffing at locus II/9A. Within the 14 kb, several matches to the ecdysone response element (EcRE) consensus sequence were identified, including some in the amplification origin region. EcRE O-P is at a central axis of a remarkable symmetry, equidistant to the TPs that are themselves equidistant to EcRE O-1 and EcRE O-2. DNA sequence alterations have occurred throughout the II/9A region in a newly discovered polymorphism (#2). Polymorphism #2 is not specific to developmental stage, sex, or tissue, and it does not impair DNA amplification. The DHS, both 9 bp TP sequences, and EcREs O-1, O-P, and O-2 are conserved between the polymorphism #1 and #2 sequences, suggesting their functional importance and retention during evolutionary selection. Moreover, a 72 bp sequence in the Sciara DHS at DNA puff II/9A is conserved in DNA puff C-3 of Rhynchosciara americana. Comparisons are discussed between the Sciara II/9A amplicon and the chorion locus amplicon on the third chromosome of Drosophila.

Prospect of reprogramming replication licensing for cancer drug development

Eukaryotic chromosomal DNA replication is preceded by replication licensing which involves the identification of the origin of replication by origin recognition complex (ORC). The ORC loads pre-replication complexes (pre-RCs) through a series of tightly regulated mechanisms where the ORC interacts with Cdc6 to recruit cdt1-MCM2-7 complexes to the origin of replication. In eukaryotes, adherence to regulatory mechanisms of the replication program is required to ensure that all daughter cells carry the exact copy of genetic material as the parent cell. Failure of which leads to the development of genome instability syndromes like cancer, diabetes, etc. In an event of such occurrence, preventing cells from carrying the defaulted genetic material and passing it to other cells hinges on the regulation of chromosomal DNA replication. Thus, understanding the mechanisms underpinning chromosomal DNA replication and particularly replication licensing can expose druggable enzymes, effector molecules, and secondary messengers that can be targeted for diagnosis and therapeutic purposes. Effectively drugging these molecular markers to reprogram pre-replication events can be used to control the fate of chromosomal DNA replication for the treatment of genome instability disorders and in this case, cancer. This review discusses available knowledge of replication licensing in the contest of molecular drug discovery for the treatment of cancer.

Cryo-EM of dynamic protein complexes in eukaryotic DNA replication

DNA replication in Eukaryotes is a highly dynamic process that involves several dozens of proteins. Some of these proteins form stable complexes that are amenable to high-resolution structure determination by cryo-EM, thanks to the recent advent of the direct electron detector and powerful image analysis algorithm. But many of these proteins associate only transiently and flexibly, precluding traditional biochemical purification. We found that direct mixing of the component proteins followed by 2D and 3D image sorting can capture some very weakly interacting complexes. Even at 2D average level and at low resolution, EM images of these flexible complexes can provide important biological insights. It is often necessary to positively identify the feature-of-interest in a low resolution EM structure. We found that systematically fusing or inserting maltose binding protein (MBP) to selected proteins is highly effective in these situations. In this chapter, we describe the EM studies of several protein complexes involved in the eukaryotic DNA replication over the past decade or so. We suggest that some of the approaches used in these studies may be applicable to structural analysis of other biological systems.

DNA sequence templates adjacent nucleosome and ORC sites at gene amplification origins in Drosophila

Eukaryotic origins of DNA replication are bound by the origin recognition complex (ORC), which scaffolds assembly of a pre-replicative complex (pre-RC) that is then activated to initiate replication. Both pre-RC assembly and activation are strongly influenced by developmental changes to the epigenome, but molecular mechanisms remain incompletely defined. We have been examining the activation of origins responsible for developmental gene amplification in Drosophila. At a specific time in oogenesis, somatic follicle cells transition from genomic replication to a locus-specific replication from six amplicon origins. Previous evidence indicated that these amplicon origins are activated by nucleosome acetylation, but how this affects origin chromatin is unknown. Here, we examine nucleosome position in follicle cells using micrococcal nuclease digestion with Ilumina sequencing. The results indicate that ORC binding sites and other essential origin sequences are nucleosome-depleted regions (NDRs). Nucleosome position at the amplicons was highly similar among developmental stages during which ORC is or is not bound, indicating that being an NDR is not sufficient to specify ORC binding. Importantly, the data suggest that nucleosomes and ORC have opposite preferences for DNA sequence and structure. We propose that nucleosome hyperacetylation promotes pre-RC assembly onto adjacent DNA sequences that are disfavored by nucleosomes but favored by ORC.

Induction of endocycles represses apoptosis independently of differentiation and predisposes cells to genome instability

The endocycle is a common developmental cell cycle variation wherein cells become polyploid through repeated genome duplication without mitosis. We previously showed that Drosophila endocycling cells repress the apoptotic cell death response to genotoxic stress. Here, we investigate whether it is differentiation or endocycle remodeling that promotes apoptotic repression. We find that when nurse and follicle cells switch into endocycles during oogenesis they repress the apoptotic response to DNA damage caused by ionizing radiation, and that this repression has been conserved in the genus Drosophila over 40 million years of evolution. Follicle cells defective for Notch signaling failed to switch into endocycles or differentiate and remained apoptotic competent. However, genetic ablation of mitosis by knockdown of Cyclin A or overexpression of fzr/Cdh1 induced follicle cell endocycles and repressed apoptosis independently of Notch signaling and differentiation. Cells recovering from these induced endocycles regained apoptotic competence, showing that repression is reversible. Recovery from fzr/Cdh1 overexpression also resulted in an error-prone mitosis with amplified centrosomes and high levels of chromosome loss and fragmentation. Our results reveal an unanticipated link between endocycles and the repression of apoptosis, with broader implications for how endocycles may contribute to genome instability and oncogenesis.

RNAP-II molecules participate in the anchoring of the ORC to rDNA replication origins

The replication of genomic DNA is limited to a single round per cell cycle. The first component, which recognises and remains bound to origins from recognition until activation and replication elongation, is the origin recognition complex. How origin recognition complex (ORC) proteins remain associated with chromatin throughout the cell cycle is not yet completely understood. Several genome-wide studies have undoubtedly demonstrated that RNA polymerase II (RNAP-II) binding sites overlap with replication origins and with the binding sites of the replication components. RNAP-II is no longer merely associated with transcription elongation. Several reports have demonstrated that RNAP-II molecules affect chromatin structure, transcription, mRNA processing, recombination and DNA repair, among others. Most of these activities have been reported to directly depend on the interaction of proteins with the C-terminal domain (CTD) of RNAP-II. Two-dimensional gels results and ChIP analysis presented herein suggest that stalled RNAP-II molecules bound to the rDNA chromatin participate in the anchoring of ORC proteins to origins during the G1 and S-phases. The results show that in the absence of RNAP-II, Orc1p, Orc2p and Cdc6p do not bind to origins. Moreover, co-immunoprecipitation experiments suggest that Ser2P-CTD and hypophosphorylated RNAP-II interact with Orc1p. In the context of rDNA, cryptic transcription by RNAP-II did not negatively interfere with DNA replication. However, the results indicate that RNAP-II is not necessary to maintain the binding of ORCs to the origins during metaphase. These findings highlight for the first time the potential importance of stalled RNAP-II in the regulation of DNA replication.

Emerging players in the initiation of eukaryotic DNA replication

Faithful duplication of the genome in eukaryotes requires ordered assembly of a multi-protein complex called the pre-replicative complex (pre-RC) prior to S phase; transition to the pre-initiation complex (pre-IC) at the beginning of DNA replication; coordinated progression of the replisome during S phase; and well-controlled regulation of replication licensing to prevent re-replication. These events are achieved by the formation of distinct protein complexes that form in a cell cycle-dependent manner. Several components of the pre-RC and pre-IC are highly conserved across all examined eukaryotic species. Many of these proteins, in addition to their bona fide roles in DNA replication are also required for other cell cycle events including heterochromatin organization, chromosome segregation and centrosome biology. As the complexity of the genome increases dramatically from yeast to human, additional proteins have been identified in higher eukaryotes that dictate replication initiation, progression and licensing. In this review, we discuss the newly discovered components and their roles in cell cycle progression.

Late replication domains in polytene and non-polytene cells of Drosophila melanogaster

In D. melanogaster polytene chromosomes, intercalary heterochromatin (IH) appears as large dense bands scattered in euchromatin and comprises clusters of repressed genes. IH displays distinctly low gene density, indicative of their particular regulation. Genes embedded in IH replicate late in the S phase and become underreplicated. We asked whether localization and organization of these late-replicating domains is conserved in a distinct cell type. Using published comprehensive genome-wide chromatin annotation datasets (modENCODE and others), we compared IH organization in salivary gland cells and in a Kc cell line. We first established the borders of 60 IH regions on a molecular map, these regions containing underreplicated material and encompassing ∼12% of Drosophila genome. We showed that in Kc cells repressed chromatin constituted 97% of the sequences that corresponded to IH bands. This chromatin is depleted for ORC-2 binding and largely replicates late. Differences in replication timing between the cell types analyzed are local and affect only sub-regions but never whole IH bands. As a rule such differentially replicating sub-regions display open chromatin organization, which apparently results from cell-type specific gene expression of underlying genes. We conclude that repressed chromatin organization of IH is generally conserved in polytene and non-polytene cells. Yet, IH domains do not function as transcription- and replication-regulatory units, because differences in transcription and replication between cell types are not domain-wide, rather they are restricted to small “islands” embedded in these domains. IH regions can thus be defined as a special class of domains with low gene density, which have narrow temporal expression patterns, and so displaying relatively conserved organization.

The origin recognition complex: a biochemical and structural view

The origin recognition complex (ORC) was first discovered in the baker's yeast in 1992. Identification of ORC opened up a path for subsequent molecular level investigations on how eukaryotic cells initiate and control genome duplication each cell cycle. Twenty years after the first biochemical isolation, ORC is now taking on a three-dimensional shape, although a very blurry shape at the moment, thanks to the recent electron microscopy and image reconstruction efforts. In this chapter, we outline the current biochemical knowledge about ORC from several eukaryotic systems, with emphasis on the most recent structural and biochemical studies. Despite many species-specific properties, an emerging consensus is that ORC is an ATP-dependent machine that recruits other key proteins to form pre-replicative complexes (pre-RCs) at many origins of DNA replication, enabling the subsequent initiation of DNA replication in S phase.

Analysis of model replication origins in Drosophila reveals new aspects of the chromatin landscape and its relationship to origin activity and the prereplicative complex

A study of model DNA replication origins in Drosophila reveals a codependence between histone acetylation and pre-RC assembly and leads to a chromatin switch model for the coordination of origin and promoter activity during development.

Epigenetic regulation exerts a major influence on origins of DNA replication during development. The mechanisms for this regulation, however, are poorly defined. We showed previously that acetylation of nucleosomes regulates the origins that mediate developmental gene amplification during Drosophila oogenesis. Here we show that developmental activation of these origins is associated with acetylation of multiple histone lysines. Although these modifications are not unique to origin loci, we find that the level of acetylation is higher at the active origins and quantitatively correlated with the number of times these origins initiate replication. All of these acetylation marks were developmentally dynamic, rapidly increasing with origin activation and rapidly declining when the origins shut off and neighboring promoters turn on. Fine-scale analysis of the origins revealed that both hyperacetylation of nucleosomes and binding of the origin recognition complex (ORC) occur in a broad domain and that acetylation is highest on nucleosomes adjacent to one side of the major site of replication initiation. It was surprising to find that acetylation of some lysines depends on binding of ORC to the origin, suggesting that multiple histone acetyltransferases may be recruited during origin licensing. Our results reveal new insights into the origin epigenetic landscape and lead us to propose a chromatin switch model to explain the coordination of origin and promoter activity during development.

Preferential re-replication of Drosophila heterochromatin in the absence of geminin

To ensure genomic integrity, the genome must be duplicated exactly once per cell cycle. Disruption of replication licensing mechanisms may lead to re-replication and genomic instability. Cdt1, also known as Double-parked (Dup) in Drosophila, is a key regulator of the assembly of the pre-replicative complex (pre-RC) and its activity is strictly limited to G1 by multiple mechanisms including Cul4-Ddb1 mediated proteolysis and inhibition by geminin. We assayed the genomic consequences of disregulating the replication licensing mechanisms by RNAi depletion of geminin. We found that not all origins of replication were sensitive to geminin depletion and that heterochromatic sequences were preferentially re-replicated in the absence of licensing mechanisms. The preferential re-activation of heterochromatic origins of replication was unexpected because these are typically the last sequences to be duplicated in a normal cell cycle. We found that the re-replication of heterochromatin was regulated not at the level of pre-RC activation, but rather by the formation of the pre-RC. Unlike the global assembly of the pre-RC that occurs throughout the genome in G1, in the absence of geminin, limited pre-RC assembly was restricted to the heterochromatin by elevated cyclin A-CDK activity. These results suggest that there are chromatin and cell cycle specific controls that regulate the re-assembly of the pre-RC outside of G1.

Catastrophic consequences may occur if the cell fails to either completely copy the genome or if it duplicates some regions of the genome more than once in a cell cycle. The cell must coordinate thousands of DNA replication start sites (origins) to ensure that the entire genome is copied and that no replication origin is activated more than once in a cell cycle. The cell accomplishes this coordination by confining the selection and activation of replication origins to discrete phases of the cell cycle. Start sites can only be selected or ‘licensed’ for DNA replication in G1 and similarly, they can only be activated for the initiation of DNA replication in S phase. Disruption of the mechanisms that regulate this ‘licensing’ process have been shown to result in extensive re-replication, genomic instability and tumorigenesis in a variety of eukaryotic systems. Here we use genomic approaches in Drosophila to identify which origins of replication are susceptible to re-initiation of DNA replication in the absence of replication licensing controls. Unexpectedly, we find that sequences in the heterochromatin, which were thought to contain only inefficient origins of replication, are preferentially re-replicated. These results provide insights into how origins of replication are selected and regulated in distinct chromatin environments to maintain genomic stability.

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.

Common fragile sites are characterized by histone hypoacetylation

Common fragile sites (CFSs) represent large, highly unstable regions of the human genome. CFS sequences are sensitive to perturbation of replication; however, the molecular basis for the instability at CFSs is poorly understood. We hypothesized that a unique epigenetic pattern may underlie the unusual sensitivity of CFSs to replication interference. To examine this hypothesis, we analyzed chromatin modification patterns within the six human CFSs with the highest levels of breakage, and their surrounding non-fragile regions (NCFSs). Chromatin at most of the CFSs analyzed has significantly less histone acetylation than that of their surrounding NCFSs. Trichostatin A and/or 5-azadeoxycytidine treatment reduced chromosome breakage at CFSs. Furthermore, chromatin at the most commonly expressed CFS, the FRA3B, is more resistant to micrococcal nuclease than that of the flanking non-fragile sequences. These results demonstrate that histone hypoacetylation is a characteristic epigenetic pattern of CFSs, and chromatin within CFSs might be relatively more compact than that of the NCFSs, indicating a role for chromatin conformation in genomic instability at CFSs. Moreover, lack of histone acetylation at CFSs may contribute to the defective response to replication stress characteristic of CFSs, leading to the genetic instability characteristic of this regions.

The origin recognition complex protein family

The proteins of the origin recognition complex are found throughout all eukaryotes and have roles beyond that of DNA replication.

Origin recognition complex (ORC) proteins were first discovered as a six-subunit assemblage in budding yeast that promotes the initiation of DNA replication. Orc1-5 appear to be present in all eukaryotes, and include both AAA+ and winged-helix motifs. A sixth protein, Orc6, shows no structural similarity to the other ORC proteins, and is poorly conserved between budding yeast and most other eukaryotic species. The replication factor Cdc6 has extensive sequence similarity with Orc1 and phylogenetic analysis suggests the genes that encode them may be paralogs. ORC proteins have also been found in the archaea, and the bacterial DnaA replication protein has ORC-like functional domains. In budding yeast, Orc1-6 are bound to origins of DNA replication throughout the cell cycle. Following association with Cdc6 in G1 phase, the sequential hydrolysis of Cdc6 - then ORC-bound ATP loads the Mcm2-7 helicase complex onto DNA. Localization of ORC subunits to the kinetochore and centrosome during mitosis and to the cleavage furrow during cytokinesis has been observed in metazoan cells and, along with phenotypes observed following knockdown with short interfering RNAs, point to additional roles at these cell-cycle stages. In addition, ORC proteins function in epigenetic gene silencing through interactions with heterochromatin factors such as Sir1 in budding yeast and HP1 in higher eukaryotes. Current avenues of research have identified roles for ORC proteins in the development of neuronal and muscle tissue, and are probing their relationship to genome integrity.

Chromatin state marks cell-type- and gender-specific replication of the Drosophila genome

Duplication of eukaryotic genomes during S phase is coordinated in space and time. In order to identify zones of initiation and cell-type- as well as gender-specific plasticity of DNA replication, we profiled replication timing, histone acetylation, and transcription throughout the Drosophila genome. We observed two waves of replication initiation with many distinct zones firing in early-S phase and multiple, less defined peaks at the end of S phase, suggesting that initiation becomes more promiscuous in late-S phase. A comparison of different cell types revealed widespread plasticity of replication timing on autosomes. Most occur in large regions, but only half coincide with local differences in transcription. In contrast to confined autosomal differences, a global shift in replication timing occurs throughout the single male X chromosome. Unlike in females, the dosage-compensated X chromosome replicates almost exclusively early. This difference occurs at sites that are not transcriptionally hyperactivated, but show increased acetylation of Lys 16 of histone H4 (H4K16ac). This suggests a transcription-independent, yet chromosome-wide process related to chromatin. Importantly, H4K16ac is also enriched at initiation zones as well as early replicating regions on autosomes during S phase. Together, our study reveals novel organizational principles of DNA replication of the Drosophila genome and suggests that H4K16ac is more closely correlated with replication timing than is transcription.

Intrinsically bent DNA sites in the Drosophila melanogaster third chromosome amplified domain

Bent DNA sites promote the curvature of DNA in both eukaryotic and prokaryotic chromosomes. Here, we investigate the localization and structure of intrinsically bent DNA sites in the extensively characterized Drosophila melanogaster third chromosome DAFC-66D segment (Drosophila amplicon in the follicle cells). This region contains the amplification control element ACE3, which is a replication enhancer that acts in cis to activate the major replication origin ori-beta. Through both electrophoretic and in silico analysis, we have identified three major bent DNA sites in DAFC-66D. The bent DNA site (b1) is localized in the ACE3 element, whereas the other two bent DNA sites (b2 and b3) are localized in the ori-beta region. Four additional bent DNA sites were identified in the intron of the S18 gene and near the TATA box of the S15, S19, and S16 genes. The identification of DNA bent sites in genomic regions previously characterized as functionally relevant for DNA amplification further supports a function for DNA bent sites in DNA replication in eukaryotes.

Detection of replication origins using comparative genomics and recombinational ARS assay

Effective experimental techniques are available to identify replication origin regions in eukaryotic cells. Genome-wide identification of the precise sequence elements that direct origin activity is however still not straightforward, even in the yeast Saccharomyces cerevisiae which has the best characterised eukaryotic replication origins. The availability of genome sequences for a series of closely related (sensu stricto) budding yeasts has allowed us to take a 'comparative genomics' approach to this problem. Since they represent functional protein-binding sites, origin sequences are conserved better than the surrounding intergenic sequence within the genomes of closely related yeasts. We describe here how phylogenetic comparison data can be used to identify candidate replication origin sequences in the S. cerevisiae genome, and how large numbers of such candidate sites can simultaneously be assayed for ability to initiate replication. Similar approaches could potentially be used to identify protein-binding sequence elements having other functions, as well as replication origin sites in other species.

Beyond heterochromatin: SIR2 inhibits the initiation of DNA replication

Over the last decade, data have accumulated that support a role for chromatin structure in regulating the initiation of DNA replication and its timing during S-phase. However, the mechanisms underlying how chromatin structure influences replication initiation are not always understood. For example, in Drosophila histone acetylation at the ACE3 and Ori-beta sequences near one of the amplified chorion loci is correlated with ORC (origin recognition complex) binding and re-replication of this locus. Whether histone acetylation promotes ORC binding or some later step in replication is not known. In yeast, hypo-acetylated heterochromatin and telomeric regions replicate late in S-phase but the mechanisms that restrict the initiation of replication at these loci are not fully understood. Nonetheless, it seems likely that histone acetylation and other types of histone modification will significantly impact DNA replication. A recent study published in Molecular Cell reveals a role for the conserved NAD(+)-dependent histone deacetylase, Sir2, in inhibiting the assembly of the multiprotein complex necessary for the selection and activation of yeast replication origins. Here, we highlight key conclusions from this study, place them in perspective with earlier work, and outline important future questions.

Genome-wide studies highlight indirect links between human replication origins and gene regulation

To get insights into the regulation of replication initiation, we systematically mapped replication origins along 1% of the human genome in HeLa cells. We identified 283 origins, 10 times more than previously known. Origin density is strongly correlated with genomic landscapes, with clusters of closely spaced origins in GC-rich regions and no origins in large GC-poor regions. Origin sequences are evolutionarily conserved, and half of them map within or near CpG islands. Most of the origins overlap transcriptional regulatory elements, providing further evidence of a connection with gene regulation. Moreover, we identify c-JUN and c-FOS as important regulators of origin selection. Half of the identified replication initiation sites do not have an open chromatin configuration, showing the absence of a direct link with gene regulation. Replication timing analyses coupled with our origin mapping suggest that a relatively strict origin-timing program regulates the replication of the human genome.