Replication initiation patterns in the beta-globin loci of totipotent and differentiated murine cells: evidence for multiple initiation regions

Large-Scale Chromatin Structure-Function Relationships during the Cell Cycle and Development: Insights from Replication Timing

Chromosome architecture has received a lot of attention since the recent development of genome-scale methods to measure chromatin interactions (Hi-C), enabling the first sequence-based models of chromosome tertiary structure. A view has emerged of chromosomes as a string of structural units (topologically associating domains; TADs) whose boundaries persist through the cell cycle and development. TADs with similar chromatin states tend to aggregate, forming spatially segregated chromatin compartments. However, high-resolution Hi-C has revealed substructure within TADs (subTADs) that poses a challenge for models that attribute significance to structural units at any given scale. More than 20 years ago, the DNA replication field independently identified stable structural (and functional) units of chromosomes (replication foci) as well as spatially segregated chromatin compartments (early and late foci), but lacked the means to link these units to genomic map units. Genome-wide studies of replication timing (RT) have now merged these two disciplines by identifying individual units of replication regulation (replication domains; RDs) that correspond to TADs and are arranged in 3D to form spatiotemporally segregated subnuclear compartments. Furthermore, classifying RDs/TADs by their constitutive versus developmentally regulated RT has revealed distinct classes of chromatin organization, providing unexpected insight into the relationship between large-scale chromosome structure and function.

Mapping of Replication Origins in the X Inactivation Center of Vole Microtus levis Reveals Extended Replication Initiation Zone

DNA replication initiates at specific positions termed replication origins. Genome-wide studies of human replication origins have shown that origins are organized into replication initiation zones. However, only few replication initiation zones have been described so far. Moreover, few origins were mapped in other mammalian species besides human and mouse. Here we analyzed pattern of short nascent strands in the X inactivation center (XIC) of vole Microtus levis in fibroblasts, trophoblast stem cells, and extraembryonic endoderm stem cells and confirmed origins locations by ChIP approach. We found that replication could be initiated in a significant part of XIC. We also analyzed state of XIC chromatin in these cell types. We compared origin localization in the mouse and vole XIC. Interestingly, origins associated with gene promoters are conserved in these species. The data obtained allow us to suggest that the X inactivation center of M. levis is one extended replication initiation zone.

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.

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.

Specification of DNA replication origins and genomic base composition in fission yeasts

In the "Replicon Theory", Jacob, Brenner and Cuzin proposed the existence of replicators and initiators as the two major actors in DNA replication. Over the years, many protein components of initiators have been shown to be conserved in different organisms during evolution. By contrast, replicator DNA sequences (often referred to as replication origins) have diverged beyond possible comparison between eukaryotic genomes. Replication origins in the fission yeast Schizosaccharomyces pombe are made up of A+T-rich sequences that do not share any consensus elements. The information encoded in these replicators is interpreted by the Orc4 subunit of the ORC (origin recognition complex), which is unique among eukaryotes in that it contains a large domain harboring nine AT-hook subdomains that target ORC to a great variety of A+T-rich sequences along the chromosomes. Recently, the genomes of other Schizosaccharomyces species have been sequenced and the regions encompassing their replication origins have been identified. DNA sequence analysis and comparison of the organization of their Orc4 proteins have revealed species-specific differences that contribute to our understanding of how the specification of replication origins has evolved during the phylogenetic divergence of fission yeasts.

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.

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

The temporal order of replication of mammalian chromosomes appears to be linked to their functional organization, but the process that establishes and modifies this order during cell differentiation remains largely unknown. Here, we studied how the replication of the Igh locus initiates, progresses, and terminates in bone marrow pro-B cells undergoing B cell commitment. We show that many aspects of DNA replication can be quantitatively explained by a mechanism involving the stochastic firing of origins (across the S phase and the Igh locus) and extensive variations in their firing rate (along the locus). The firing rate of origins shows a high degree of coordination across Igh domains that span tens to hundreds of kilobases, a phenomenon not observed in simple eukaryotes. Differences in domain sizes and firing rates determine the temporal order of replication. During B cell commitment, the expression of the B-cell-specific factor Pax5 sharply alters the temporal order of replication by modifying the rate of origin firing within various Igh domains (particularly those containing Pax5 binding sites). We propose that, within the Igh C_H-3′RR domain, Pax5 is responsible for both establishing and maintaining high rates of origin firing, mostly by controlling events downstream of the assembly of pre-replication complexes.

Each time a mammalian cell duplicates its genome in preparation for cell division it activates thousands of so called “DNA origins of replication.” The timely and complete duplication of the genome depends on careful orchestration of origin activation, which is modified when cells differentiate to perform a specific function. We currently lack a universally accepted model of origin regulation that can explain the replication dynamics in complex eukaryotes. Here, we studied the mouse immunoglobulin heavy-chain locus, one of the antibody-encoding portions of the genome, where origins change activity when antibody-producing B cells differentiate in the bone marrow. We show that multiple aspects of DNA replication initiation, progression, and termination can be explained mathematically by the interplay between randomly firing origins and two independent variables: the speed of progression of replication forks and the firing rate of origins along the locus. The rate of origin firing varies extensively along the locus during B cell differentiation and, thus, is a dominant factor in establishing the temporal order of replication. A differentiation factor called Pax5 can alter the temporal order of replication by modifying the rate of origin firing across various parts of the locus.

Unraveling cell type-specific and reprogrammable human replication origin signatures associated with G-quadruplex consensus motifs

DNA replication is highly regulated, ensuring faithful inheritance of genetic information through each cell cycle. In metazoans, this process is initiated at many thousands of DNA replication origins whose cell type-specific distribution and usage are poorly understood. We exhaustively mapped the genome-wide location of replication origins in human cells using deep sequencing of short nascent strands and identified ten times more origin positions than we expected; most of these positions were conserved in four different human cell lines. Furthermore, we identified a consensus G-quadruplex-forming DNA motif that can predict the position of DNA replication origins in human cells, accounting for their distribution, usage efficiency and timing. Finally, we discovered a cell type-specific reprogrammable signature of cell identity that was revealed by specific efficiencies of conserved origin positions and not by the selection of cell type-specific subsets of origins.

Prevention of transcriptional silencing by a replicator-binding complex consisting of SWI/SNF, MeCP1, and hnRNP C1/C2

Transcriptional silencing selectively impedes gene expression. Silencing is often accompanied by replication delay and can be prevented by replicator sequences. Here we report a replicator-binding protein complex involved in the prevention of transcriptional silencing. The protein complex interacts with an essential asymmetric region within the human β-globin Rep-P replicator and includes hnRNP C1/C2, SWI/SNF complex, and MeCP1, which are members of the locus control region (LCR)-associated remodeling complex (LARC). Interaction between LARC and Rep-P prevented transcriptional silencing and replication delay. Transgenes that did not contain the asymmetric LARC-binding region of Rep-P replicated late and exhibited stable silencing that could not be affected by a DNA methylation inhibitor. In contrast, transgenes that contain a mutation of the asymmetric region of Rep-P that could not bind LARC exhibited a silent state that could transiently be reactivated by DNA demethylation. The effect of DNA demethylation was transient, and prolonged exposure to a methylation inhibitor induced distinct, stable, methylation-independent silencing. These observations suggest that the interaction of LARC complex with replicators plays a role in preventing gene silencing and provides support for a novel, epigenetic mechanism of resistance to methylation inhibitors.

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.

A key commitment step in erythropoiesis is synchronized with the cell cycle clock through mutual inhibition between PU.1 and S-phase progression

Hematopoietic progenitors undergo differentiation while navigating several cell division cycles, but it is unknown whether these two processes are coupled. We addressed this question by studying erythropoiesis in mouse fetal liver in vivo. We found that the initial upregulation of cell surface CD71 identifies developmentally matched erythroblasts that are tightly synchronized in S-phase. We show that DNA replication within this but not subsequent cycles is required for a differentiation switch comprising rapid and simultaneous committal transitions whose precise timing was previously unknown. These include the onset of erythropoietin dependence, activation of the erythroid master transcriptional regulator GATA-1, and a switch to an active chromatin conformation at the β-globin locus. Specifically, S-phase progression is required for the formation of DNase I hypersensitive sites and for DNA demethylation at this locus. Mechanistically, we show that S-phase progression during this key committal step is dependent on downregulation of the cyclin-dependent kinase p57(KIP2) and in turn causes the downregulation of PU.1, an antagonist of GATA-1 function. These findings therefore highlight a novel role for a cyclin-dependent kinase inhibitor in differentiation, distinct to their known function in cell cycle exit. Furthermore, we show that a novel, mutual inhibition between PU.1 expression and S-phase progression provides a "synchromesh" mechanism that "locks" the erythroid differentiation program to the cell cycle clock, ensuring precise coordination of critical differentiation events.

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 multiprotein complex necessary for both transcription and DNA replication at the β-globin locus

DNA replication, repair, transcription and chromatin structure are intricately associated nuclear processes, but the molecular links between these events are often obscure. In this study, we have surveyed the protein complexes that bind at β-globin locus control region, and purified and characterized the function of one such multiprotein complex from human erythroleukemic K562 cells. We further validated the existence of this complex in human CD34+ cell-derived normal erythroid cells. This complex contains ILF2/ILF3 transcription factors, p300 acetyltransferase and proteins associated with DNA replication, transcription and repair. RNAi knockdown of ILF2, a DNA-binding component of this complex, abrogates the recruitment of the complex to its cognate DNA sequence and inhibits transcription, histone acetylation and usage of the origin of DNA replication at the β-globin locus. These results imply a direct link between mammalian DNA replication, transcription and histone acetylation mediated by a single multiprotein complex.

Single-molecule analysis reveals changes in the DNA replication program for the POU5F1 locus upon human embryonic stem cell differentiation

Human embryonic stem cells (hESCs), due to their pluripotent nature, represent a particularly relevant model system to study the relationship between the replication program and differentiation state. Here, we define the basic properties of the replication program in hESCs and compare them to the programs of hESC-derived multipotent cells (neural rosette cells) and primary differentiated cells (microvascular endothelial cells [MECs]). We characterized three genomic loci: two pluripotency regulatory genes, POU5F1 (OCT4) and NANOG, and the IGH locus, a locus that is transcriptionally active specifically in B-lineage cells. We applied a high-resolution approach to capture images of individual replicated DNA molecules. We demonstrate that for the loci studied, several basic properties of replication, including the average speed of replication forks and the average density of initiation sites, were conserved among the cells analyzed. We also demonstrate, for the first time, the presence of initiation zones in hESCs. However, significant differences were evident in other aspects of replication for the DNA segment containing the POU5F1 gene. Specifically, the locations of centers of initiation zones and the direction of replication fork progression through the POU5F1 gene were conserved in two independent hESC lines but were different in hESC-derived multipotent cells and MECs. Thus, our data identify features of the replication program characteristic of hESCs and define specific changes in replication during hESC differentiation.

Cell cycle-dependent accumulation of histone H3.3 and euchromatic histone modifications in pericentromeric heterochromatin in response to a decrease in DNA methylation levels

In mammals, DNA methylation is an important epigenetic mark that is associated with gene silencing, particularly in constitutive heterochromatin. However, the effect of DNA methylation on other epigenetic properties of chromatin is controversial. In this study, we show that inhibition of DNA methylation in mouse fibroblast cells affects histone modification and the subnuclear localization of histone H3.3 in a cell cycle-dependent manner. Using a DNA methyltransferase (Dnmt) inhibitor 5-aza-2'-deoxycytidine (5-aza-dC), we found that reduced levels of DNA methylation were associated with the activation of transcription from centromeric and pericentromeric satellite repeats. The de-repressed pericentromeric chromatin was enriched in euchromatic histone modifications such as acetylation of histone H4, and di- and tri-methylation of lysine 4 on histone H3. Spatio-temporal analysis showed that the accumulation of these euchromatic histone modifications occurred during the second S phase following 5-aza-dC treatment, corresponding precisely with a shift in replication timing of the pericentromeric satellite repeats from middle/late S phase to early S phase. Moreover, we found that histone H3.3 was deposited on the pericentromeric heterochromatin prior to the accumulation of the euchromatic histone modifications. These results suggest that DNA CpG methylation is essential for the proper organization of pericentromeric heterochromatin in differentiated mouse cells.

Domain-wide regulation of DNA replication timing during mammalian development

Studies of replication timing provide a handle into previously impenetrable higher-order levels of chromosome organization and their plasticity during development. Although mechanisms regulating replication timing are not clear, novel genome-wide studies provide a thorough survey of the extent to which replication timing is regulated during most of the early cell fate transitions in mammals, revealing coordinated changes of a defined set of 400-800 kb chromosomal segments that involve at least half the genome. Furthermore, changes in replication time are linked to changes in sub-nuclear organization and domain-wide transcriptional potential, and tissue-specific replication timing profiles are conserved from mouse to human, suggesting that the program has developmental significance. Hence, these studies have provided a solid foundation for linking megabase level chromosome structure to function, and suggest a central role for replication in domain-level genome organization.

Structural diversity and dynamics of genomic replication origins in Schizosaccharomyces pombe

DNA replication origins (ORI) in Schizosaccharomyces pombe colocalize with adenine and thymine (A+T)-rich regions, and earlier analyses have established a size from 0.5 to over 3 kb for a DNA fragment to drive replication in plasmid assays. We have asked what are the requirements for ORI function in the chromosomal context. By designing artificial ORIs, we have found that A+T-rich fragments as short as 100 bp without homology to S. pombe DNA are able to initiate replication in the genome. On the other hand, functional dissection of endogenous ORIs has revealed that some of them span a few kilobases and include several modules that may be as short as 25-30 contiguous A+Ts capable of initiating replication from ectopic chromosome positions. The search for elements with these characteristics across the genome has uncovered an earlier unnoticed class of low-efficiency ORIs that fire late during S phase. These results indicate that ORI specification and dynamics varies widely in S. pombe, ranging from very short elements to large regions reminiscent of replication initiation zones in mammals.

Cyclin A-Cdk1 regulates the origin firing program in mammalian cells

Somatic mammalian cells possess well-established S-phase programs with specific regions of the genome replicated at precise times. The ATR-Chk1 pathway plays a central role in these programs, but the mechanism for how Chk1 regulates origin firing remains unknown. We demonstrate here the essential role of cyclin A2-Cdk1 in the regulation of late origin firing. Activity of cyclin A2-Cdk1 was hardly detected at the onset of S phase, but it was obvious at middle to late S phase under unperturbed condition. Chk1 depletion resulted in increased expression of Cdc25A, subsequent hyperactivation of cyclin A2-Cdk1, and abnormal replication at early S phase. Hence, the ectopic expression of cyclin A2-Cdk1AF (constitutively active mutant) fusion constructs resulted in abnormal origin firing, causing the premature appearance of DNA replication at late origins at early S phase. Intriguingly, inactivation of Cdk1 in temperature-sensitive Cdk1 mutant cell lines (FT210) resulted in a prolonged S phase and inefficient activation of late origin firing even at late S phase. Our results thus suggest that cyclin A2-Cdk1 is a key regulator of S-phase programs.

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.

Cooperative sequence modules determine replication initiation sites at the human beta-globin locus

The human beta globin locus contains two adjacent replicators, each capable of initiating DNA replication when transferred from its native locus to ectopic sites. Here, we report a detailed analysis of the sequence requirements for replication initiation from these replicators. In both replicators, initiation required a combination of an asymmetric purine:pyrimidine sequence and several AT-rich stretches. Modules from the two replicators could combine to initiate replication. AT-rich sequences were essential for replicator activity: a low frequency of initiation was observed in DNA fragments that included a short stretch of AT-rich sequences, whereas inclusion of additional AT-rich stretches increased initiation efficiency. By contrast, replication initiated at a low level without the asymmetric purine:pyrimidine modules but they were required in synergy to achieve efficient initiation. These data support a combinatorial model for replicator activity and suggest that the initiation of DNA replication requires interaction between at least two distinct sequence modules.

HoxB domain induction silences DNA replication origins in the locus and specifies a single origin at its boundary

In multicellular organisms, changes in the DNA replication programme could act to integrate differentiation with cell division in various developmental and transcriptional contexts. Here, we have addressed the use of DNA replication origins during differentiation in the HoxB domain-a cluster of nine genes developmentally regulated in a collinear manner. In undifferentiated mouse P19 cells, we detected several DNA replication origins in the 100 kb HoxB locus, indicating a relaxed origin use when the locus is transcriptionally silent. By contrast, in retinoic-acid-induced differentiated cells, when HoxB transcription is activated, a general silencing of DNA replication origins occurs in the locus except one located downstream of Hoxb1, at the 3' boundary of the HoxB domain. Silencing of the replication origins is associated with histone hyperacetylation, whereas the active Hoxb1 origin persists as a hypoacetylated island. These findings provide direct evidence for the differentiated use of origins in HoxB genes, and we suggest that this regulation might contribute to the regulated expression of HoxB genes during development.

Broadening of DNA replication origin usage during metazoan cell differentiation

We have examined whether replication of the chicken beta-globin locus changes during differentiation of primary erythroid progenitors into erythrocytes. In undifferentiated progenitors, four principal initiation sites and a replication fork pausing region (RFP) were observed. Forty-eight hours after induction of differentiation, the principal sites were maintained, even in the activated beta(A)-globin gene, some minor sites were enhanced, three new sites appeared and the RFP disappeared. One of the activated origins showed increased histone H3 K9K14 diacetylation, but the others did not. These results demonstrate a broadening of DNA replication origin usage during differentiation of untransformed metazoan cells and indicate that histone H3 diacetylation, other histone modifications so far reported and transcription are not crucial determinants of origin selection in this system.

Mapping of DNA replication origins to noncoding genes of the X-inactivation center

In mammals, few DNA replication origins have been identified. Although there appears to be an association between origins and epigenetic regulation, their underlying link to monoallelic gene expression remains unclear. Here, we identify novel origins of DNA replication (ORIs) within the X-inactivation center (Xic). We analyze 86 kb of the Xic using an unbiased approach and find an unexpectedly large number of functional ORIs. Although there has been a tight correlation between ORIs and CpG islands, we find that ORIs are not restricted to CpG islands and there is no dependence on transcriptional activity. Interestingly, these ORIs colocalize to important genetic elements or genes involved in X-chromosome inactivation. One prominent ORI maps to the imprinting center and to a domain within Tsix known to be required for X-chromosome counting and choice. Location and/or activity of ORIs appear to be modulated by removal of specific Xic elements. These data provide a foundation for testing potential relationships between DNA replication and epigenetic regulation in future studies.

Replication initiation from a novel origin identified in the Th2 cytokine cluster locus requires a distant conserved noncoding sequence

Lineage commitment of Th cells is associated with the establishment of specific transcriptional programs of cytokines. However, how Th cell differentiation affects the program of DNA replication has not been addressed. To gain insight into interplays between differentiation-induced transcription regulation and initiation of DNA replication, we took advantage of an in vitro differentiation system of naive T cells, in which one can manipulate their differentiation into Th1 or Th2 cells. We searched for replication origins in the murine IL-4/IL-13 locus and compared their profiles in the two Th cell lineages which were derived in vitro from the same precursor T cells. We identified a replication origin (ori(IL-13)) downstream from exon 4 of IL-13 and showed that this origin functions in both Th2 and Th1 cells. A distant regulatory element called CNS-1 (conserved noncoding sequence 1) in the IL-4/IL-13 intergenic region coincides with a Th2-specific DNase I-hypersensitive site and is required for efficient, coordinated expression of Th2 cytokines. Replication initiation from ori(IL-13) is significantly reduced in Th1 and Th2 cells derived from CNS-1-deficient mice. However, the replication timing of this locus is consistently early during S phase in both Th1 and Th2 cells under either the wild-type or CNS-1 deletion background. Thus, the conserved noncoding element in the intergenic region regulates replication initiation from a distant replication origin in a manner independent from its effect on lineage-specific transcription but not the replication timing of the segment surrounding this origin.

Preventing gene silencing with human replicators

Transcriptional silencing, one of the major impediments to gene therapy in humans, is often accompanied by replication during late S-phase. We report that transcriptional silencing and late replication were prevented by DNA sequences that can initiate DNA replication (replicators). When replicators were included in silencing-prone transgenes, they did not undergo transcriptional silencing, replicated early and maintained histone acetylation patterns characteristic of euchromatin. A mutant replicator, which could not initiate replication, could not prevent gene silencing and replicated late when included in identical transgenes and inserted at identical locations. These observations suggest that replicators introduce epigenetic chromatin changes that facilitate initiation of DNA replication and affect gene silencing. Inclusion of functional replicators in gene therapy vectors may provide a tool for stabilizing gene expression patterns.

Recruitment of ORC or CDC6 to DNA is sufficient to create an artificial origin of replication in mammalian cells

Origins of replication are expected to recruit initiation proteins like origin recognition complex (ORC) and Cdc6 in eukaryotes and provide a platform for unwinding DNA. Here we test whether localization of initiation proteins onto DNA is sufficient for origin function. Different components of the ORC complex and Cdc6 stimulated prereplicative complex (pre-RC) formation and replication initiation when fused to the GAL4 DNA-binding domain and recruited to plasmid DNA containing a tandem array of GAL4-binding sites. Replication occurred once per cell cycle and was inhibited by Geminin, indicating that the plasmid was properly licensed during the cell cycle. The GAL4 fusion protein recruits other polypeptides of the ORC-Cdc6 complex, and nascent strand abundance was highest near the GAL4-binding sites. Therefore, the artificial origin recapitulates many of the regulatory features of physiological origins and is valuable for studies on replication initiation in mammalian cells. We demonstrated the utility of this system by showing the functional importance of the ATPase domains of human Cdc6 and Orc1 and the dispensability of the N-terminal segments of Orc1 and Orc2 in this assay. Artificial recruitment of a eukaryotic cellular replication initiation factor to a DNA sequence can create a functional origin of replication, providing a robust genetic assay for these factors and a novel approach to generating episomal vectors for gene therapy.

Initiation of DNA replication at the human beta-globin 3' enhancer

The origin of DNA replication in the human β-globin gene contains an initiation region (IR) and two flanking auxiliary elements. Two replicator modules are located within the upstream auxiliary sequence and the IR core, but the functional sequences in the downstream auxiliary element are unknown. Here, we use a combination of benzoylated-naphthoylated DEAE (BND) cellulose purification and nascent strand abundance assays to show that replication initiation occurs at the β-globin 3′ enhancer on human chromosome 11 in the Hu11 hybrid murine erythroleukemia (MEL) cell line. To examine replicator function, 3′ enhancer fragments were inserted into an ectopic site in MEL cells via an optimized FRT/EGFP-FLP integration system. These experiments demonstrate that the 1.6 kb downstream auxiliary element is a third replicator module called bGRep-E in erythroid cells. The minimal 260 bp 3′ enhancer is required but not sufficient to initiate efficient replication, suggesting cooperation with adjacent sequences. The minimal 3′ enhancer also cooperates with elements in an expressing HS3β/γ-globin construct to initiate replication. These data indicate that the β-globin replicator has multiple initiation sites in three closely spaced replicator modules. We conclude that a mammalian enhancer can cooperate with adjacent sequences to create an efficient replicator module.

The replicon revisited: an old model learns new tricks in metazoan chromosomes

The origins of DNA replication were proposed in the replicon model to be specified genetically by replicator elements that coordinate the initiation of DNA synthesis with gene expression and cell growth. Recent studies have identified DNA sequences in mammalian cells that fulfil the genetic criteria for replicators and are beginning to uncover the sequence requirements for the initiation of DNA replication. Mammalian replicators are com- posed of non-redundant modules that cooperate to direct initiation to specific chromosomal sites. Conversely, replicators do not show strong sequence similarity, and their ability to initiate replication depends on the chromosomal context and epigenetic factors, as well as their primary sequence. Here, we review the properties of metazoan replicators, and discuss the genetic and epigenetic factors that determine where and when DNA replication is initiated.

Complementation of replication origin function in mouse embryonic stem cells by human DNA sequences

A functional origin of replication was mapped to the transcriptional promoter and exon 1 of the hypoxanthine-guanine phosphoribosyltransferase (HPRT) gene in the mouse and human genomes. This origin was lost in mouse embryonic stem (ES) cells with a spontaneous deletion (approximately 36 kb) at the 5' end of the HPRT locus. Restoration of HPRT activity by homologous recombination with human/mouse chimeric sequences reconstituted replication origin activity in two independent ES cell lines. Quantitative PCR analyses of abundance of genetic markers in size-fractionated nascent DNA indicated that initiation of DNA replication coincided with the site of insertion in the mouse genome of the 335 bp of human DNA containing the HPRT exon 1 and a truncated promoter. The genetic information contained in the human sequence and surrounding mouse DNA was analyzed for cis-acting elements that might contribute to selection and functional activation of a mammalian origin of DNA replication.

The histone deacetylase inhibitor trichostatin A alters the pattern of DNA replication origin activity in human cells

Eukaryotic chromatin structure limits the initiation of DNA replication spatially to chromosomal origin zones and temporally to the ordered firing of origins during S phase. Here, we show that the level of histone H4 acetylation correlates with the frequency of replication initiation as measured by the abundance of short nascent DNA strands within the human c-myc and lamin B2 origins, but less well with the frequency of initiation across the β-globin locus. Treatment of HeLa cells with trichostatin A (TSA) reversibly increased the acetylation level of histone H4 globally and at these initiation sites. At all three origins, TSA treatment transiently promoted a more dispersive pattern of initiations, decreasing the abundance of nascent DNA at previously preferred initiation sites while increasing the nascent strand abundance at lower frequency genomic initiation sites. When cells arrested in late G₁ were released into TSA, they completed S phase more rapidly than untreated cells, possibly due to the earlier initiation from late-firing origins, as exemplified by the β-globin origin. Thus, TSA may modulate replication origin activity through its effects on chromatin structure, by changing the selection of initiation sites, and by advancing the time at which DNA synthesis can begin at some initiation sites.

Bidirectional replication initiates at sites throughout the mitochondrial genome of birds

Analysis of mitochondrial replication intermediates of Gallus gallus on fork-direction gels indicates that replication occurs in both directions around circular mitochondrial DNA. This finding was corroborated by a study of chick mitochondrial DNA on standard neutral two-dimensional agarose gels, which yielded archetypal initiation arcs in fragments covering the entire genome. There was, however, considerable variation in initiation arc intensity. The majority of initiation events map to regions flanking the major non-coding region, in particular the NADH dehydrogenase subunit 6 (ND6) gene. Initiation point mapping of the ND6 gene identified prominent free 5' ends of DNA, which are candidate start sites for DNA synthesis. Therefore we propose that the initiation zone of G. gallus mitochondrial DNA encompasses most, if not all, of the genome, with preferred initiation sites in regions flanking the major non-coding region. Comparison with mammals suggests a common mechanism of initiation of mitochondrial DNA replication in higher vertebrates.

Differentiation-induced replication-timing changes are restricted to AT-rich/long interspersed nuclear element (LINE)-rich isochores

The replication timing of some genes is developmentally regulated, but the significance of replication timing to cellular differentiation has been difficult to substantiate. Studies have largely been restricted to the comparison of a few genes in established cell lines derived from different tissues, and most of these genes do not change replication timing. Hence, it has not been possible to predict how many or what types of genes might be subject to such control. Here, we have evaluated the replication timing of 54 tissue-specific genes in mouse embryonic stem cells before and after differentiation to neural precursors. Strikingly, genes residing within isochores rich in GC and poor in long interspersed nuclear elements (LINEs) did not change their replication timing, whereas half of genes within isochores rich in AT and long interspersed nuclear elements displayed programmed changes in replication timing that accompanied changes in gene expression. Our results provide direct evidence that differentiation-induced autosomal replication-timing changes are a significant part of mammalian development, provide a means to predict genes subject to such regulation, and suggest that replication timing may be more related to the evolution of metazoan genomes than to gene function or expression pattern.

Initiation sites for human DNA replication at a putative ribulose-5-phosphate 3-epimerase gene

Replication of the human genome requires the activation of thousands of replicons distributed along each one of the chromosomes. Each replicon contains an initiation, or origin, site, at which DNA synthesis begins. However, very little information is known about the nature and positioning of these initiation sites along human chromosomes. We have recently focused our attention to a 1.1 kb region of human chromosome 2 which functioned as an episomal origin in the yeast Saccharomyces cerevisiae. This region corresponded to the largest exon of a putative ribulose-5-phosphate-3-epimerase gene (RPE). In the present study we have used a real-time PCR-based nascent strand DNA abundance assay to map initiation sites for DNA replication in in vivo human chromosomes around a 13.4 kb region encompassing the putative RPE gene. By applying this analysis to a 1-1.4 kb nascent strand DNA fraction isolated from both normal skin fibroblasts, and the breast cell line MCF10; we have identified five initiation sites within the 13.4 kb region of chromosome 2. The initiation sites appear to map to similar positions in both cell lines and occur outside the coding regions of the putative RPE gene.

DNA replication initiates at domains overlapping with nuclear matrix attachment regions in the xenopus and mouse c-myc promoter

Only a very few origins have been mapped in different multicellular organisms, and they do not share detectable consensus sequence elements. Moreover, it is not clear if origins are localized at similar positions in the corresponding locus in genomes of different organisms. Here, we have mapped DNA replication origins in the c-myc locus both in Xenopus and mouse, allowing a comparison of the corresponding sites in three different animal species (Xenopus, mouse, human). An origin of DNA replication is present in the three homologous c-myc loci. In Xenopus, a main DNA replication origin was located 3 kilobases (kb) upstream of the active c-myc promoter, whereas, in mouse, we detected an origin 1 kb upstream of the promoter, as previously mapped in human c-myc. We also identified a nuclear matrix attachment region in both Xenopus and mouse, which is localized to two different regions of the c-myc promoter region. However, in both cases, the nuclear matrix attachment sites are close to the DNA replication origin mapped in the locus. These data suggest that global features of chromatin organization in different organisms may contribute to DNA replication origin localization.

The human beta-globin replication initiation region consists of two modular independent replicators

Previous studies have shown that mammalian cells contain replicator sequences, which can determine where DNA replication initiates. However, the specific sequences that confer replicator activity were not identified. Here we report a detailed analysis of replicator sequences that dictate initiation of DNA replication from the human beta-globin locus. This analysis suggests that the beta-globin replication initiation region contains two adjacent, redundant replicators. Each replicator was capable of initiating DNA replication independently at ectopic sites. Within each of these two replicators, we identified short, discrete, nonredundant sequences, which cooperatively determine replicator activity. Experiments with somatic cell hybrids further demonstrated that the requirements for initiation at ectopic sites were similar to the requirements for initiation within native human chromosomes. The replicator clustering and redundancy exemplified in the human beta-globin locus may account for the extreme difficulty in identifying replicator sequences in mammalian cells and suggest that mammalian replication initiation sites may be determined by cooperative sequence modules.

Identification of novel initiation sites for human DNA replication around ARSH1, a previously characterized yeast replicator

Replication of mammalian chromosomes depends on the activation of a large number of origins of DNA replication distributed along the chromosomes. We have focused our attention on a human DNA region, named ARSH1, localized to chromosome 2, that had been previously shown to act as an episomal origin in the yeast Saccharomyces cerevisiae. In the present study we have used a nascent strand DNA abundance assay to map initiation sites for DNA replication in in vivo human chromosomes around a 5 kb region encompassing ARSH1. This analysis applied to a 1-1.4 kb nascent DNA strand fraction isolated from normal skin fibroblasts revealed the presence of two major initiations sites surrounding the ARSH1 region. With an equivalent DNA fraction obtained from HeLa cells, in addition to these sites, a broad initiation profile was observed which included the ARSH1 region. This DNA region however was not sufficient to support episomal replication of an ARSH1-containing plasmid transfected into HeLa cells.

Mammalian mitochondrial DNA replicates bidirectionally from an initiation zone

Previous data from our laboratory suggested that replication of mammalian mitochondrial DNA initiates exclusively at or near to the formerly designated origin of heavy strand replication, OH, and proceeds unidirectionally from that locus. New results obtained using two-dimensional agarose gel electrophoresis of replication intermediates demonstrate that replication of mitochondrial DNA initiates from multiple origins across a broad zone. After fork arrest near OH, replication is restricted to one direction only. The initiation zone of bidirectional replication includes the genes for cytochrome b and NADH dehydrogenase subunits 5 and 6.

Dynamics of DNA replication in mammalian somatic cells: nucleotide pool modulates origin choice and interorigin spacing

Selection of active origins and regulation of interorigin spacing are poorly understood in mammalian cells. Using tricolor analysis of combed DNA molecules, we studied an amplified locus containing the known origin, oriGNAI3. We visualized replication firing events at this and other discrete regions and established a strict correlation between AT richness and initiation sites. We found that oriGNAI3 is the prominent origin of the domain, the firing of which correlates with silencing of neighboring sites and establishes large interorigin distances. We demonstrate that cells reversibly respond to a reduction in nucleotide availability by slowing the rate of replication fork progression; in addition, the efficiency of initiation at oriGNAI3 is lowered while other normally dormant origins in the region are activated, which results in an overall increase in the density of initiation events. Thus, nucleotide pools are involved in the specification of active origins, which in turn defines their density along chromosomes.

Dynamic alterations of replication timing in mammalian cells

BACKGROUND

The eukaryotic genome is divided into distinct replication timing domains, which are activated during S phase in a strictly conserved order. Cellular differentiation can alter replication timing in some loci, but recent experiments yielded conflicting data regarding the relationship between gene expression and replication timing. The genetic and epigenetic determinants of replication timing in mammalian cells have yet to be elucidated.

RESULTS

We developed a mammalian experimental system in which the timing of DNA replication can be altered in a controlled manner. This system utilizes sequences from the human beta-globin locus that exhibit orientation-dependent transcriptional silencing when inserted into the murine genome. We found that before insertion, the murine target site replicated late during S phase. After insertion, replication timing depended on the orientation of the transgene. In a transcription-permissive orientation, the transgene and flanking sequences replicated early. In the reverse (silencing-prone) orientation, these sequences replicated late. Early replication correlated with histone modifications of the transgene chromatin but could be observed in the absence of the beta-globin promoter. Importantly, the replication timing switch did not require a replication origin within the transgene.

CONCLUSIONS

Transgene insertions into mammalian heterochromatin can alter the timing of DNA replication at the insertion site. This differentiation-independent replication timing switch did not necessitate insertion of an active promoter or a replication origin. These observations suggest that the timing of DNA replication can be manipulated by changes in DNA sequence, but that the determinants of replication timing are distinct from the sequences that specify replication initiation sites.

Replication of the chicken beta-globin locus: early-firing origins at the 5' HS4 insulator and the rho- and betaA-globin genes show opposite epigenetic modifications

Chromatin structure is believed to exert a strong effect on replication origin function. We have studied the replication of the chicken beta-globin locus, whose chromatin structure has been extensively characterized. This locus is delimited by hypersensitive sites (HSs) that mark the position of insulator elements. A stretch of condensed chromatin and another HS separate the beta-globin domain from an adjacent folate receptor (FR) gene. We demonstrate here that in erythroid cells that express the FR but not the globin genes, replication initiates at four sites within the beta-globin domain, one at the 5' HS4 insulator and the other three near the rho- and beta(A)-globin genes. Three origins consist of G+C-rich sequences enriched in CpG dinucleotides. The fourth origin is A+T rich. Together with previous work, these data reveal that the insulator origin has unmethylated CpGs, hyperacetylated histones H3 and H4, and lysine 4-methylated histone H3. In contrast, opposite modifications are observed at the other G+C-rich origins. We also show that the whole region, including the stretch of condensed chromatin, replicates early in S phase in these cells. Therefore, different early-firing origins within the same locus may have opposite patterns of epigenetic modifications. The role of insulator elements in DNA replication is discussed.

Multiple functional elements comprise a Mammalian chromosomal replicator

The structure of replication origins in metazoans is only nominally similar to that in model organisms, such as Saccharomyces cerevisiae. By contrast to the compact origins of budding yeast, in metazoans multiple elements act as replication start sites or control replication efficiency. We first reported that replication forks diverge from an origin 5' to the human c-myc gene and that a 2.4-kb core fragment of the origin displays autonomous replicating sequence activity in plasmids and replicator activity at an ectopic chromosomal site. Here we have used clonal HeLa cell lines containing mutated c-myc origin constructs integrated at the same chromosomal location to identify elements important for DNA replication. Replication activity was measured before or after integration of the wild-type or mutated origins using PCR-based nascent DNA abundance assays. We find that deletions of several segments of the c-myc origin, including the DNA unwinding element and transcription factor binding sites, substantially reduced replicator activity, whereas deletion of the c-myc promoter P1 had only a modest effect. Substitution mutagenesis indicated that the sequence of the DNA unwinding element, rather than the spacing of flanking sequences, is critical. These results identify multiple functional elements essential for c-myc replicator activity.

Replicating by the clock

The eukaryotic genome is divided into well-defined DNA regions that are programmed to replicate at different times during S phase. Active genes are generally associated with early replication, whereas inactive genes replicate late. This expression pattern might be facilitated by the differential restructuring of chromatin at the time of replication in early or late S phase.

The origin of a developmentally regulated Igh replicon is located near the border of regulatory domains for Igh replication and expression

The 3' Ig heavy chain locus (Igh) regulatory region is the most downstream known element of the murine Igh gene cluster. We report here that the nearest non-Igh genes-Crip, Crp2, and Mta1-are located approximately 70 kb further downstream and are beyond the end of the domain of Igh transcriptional regulation. We have localized an origin of replication in MEL cells to a 3-kb segment located between the 3' Igh regulatory region and Crip. Sequences downstream of this origin are replicated by forks that move in both directions. Sequences upstream of this origin (Igh-C, -D, and -J) are replicated in a single direction through a 500-kb segment in which no active bidirectional origins can be detected. We propose that this origin may lie at or near the end of the Igh regulation domain.

Site-specific ORC binding, pre-replication complex assembly and DNA synthesis at Schizosaccharomyces pombe replication origins

Previous studies have shown that the Schizo saccharomyces pombe Orc4 subunit is solely responsible for in vitro binding of origin recognition complex (ORC) to specific AT-rich sites within S.pombe replication origins. Using ARS3001, a S.pombe replication origin consisting of four genetically required sites, we show that, in situ as well as in vitro, Orc4 binds strongly to the Delta3 site, weakly to the Delta6 site and not at all to the remaining sequences. In situ, the footprint over Delta3 is extended during G(1) phase, but only when Cdc18 is present and Mcm proteins are bound to chromatin. Moreover, this footprint extends into the adjacent Delta2 site, where leading strand DNA synthesis begins. Therefore, we conclude that ARS3001 consists of a single primary ORC binding site that assembles a pre-replication complex and initiates DNA synthesis, plus an additional novel origin element (Delta9) that neither binds ORC nor functions as a centromere, but does bind an as yet unidentified protein throughout the cell cycle. Schizosaccharomyces pombe may be an appropriate paradigm for the complex origins found in the metazoa.