High-throughput mapping of origins of replication in human cells

The Protective Role of Dormant Origins in Response to Replicative Stress

Genome stability requires tight regulation of DNA replication to ensure that the entire genome of the cell is duplicated once and only once per cell cycle. In mammalian cells, origin activation is controlled in space and time by a cell-specific and robust program called replication timing. About 100,000 potential replication origins form on the chromatin in the gap 1 (G1) phase but only 20⁻30% of them are active during the DNA replication of a given cell in the synthesis (S) phase. When the progress of replication forks is slowed by exogenous or endogenous impediments, the cell must activate some of the inactive or "dormant" origins to complete replication on time. Thus, the many origins that may be activated are probably key to protect the genome against replication stress. This review aims to discuss the role of these dormant origins as safeguards of the human genome during replicative stress.

The replicative helicase MCM recruits cohesin acetyltransferase ESCO2 to mediate centromeric sister chromatid cohesion

Chromosome segregation depends on sister chromatid cohesion which is established by cohesin during DNA replication. Cohesive cohesin complexes become acetylated to prevent their precocious release by WAPL before cells have reached mitosis. To obtain insight into how DNA replication, cohesion establishment and cohesin acetylation are coordinated, we analysed the interaction partners of 55 human proteins implicated in these processes by mass spectrometry. This proteomic screen revealed that on chromatin the cohesin acetyltransferase ESCO2 associates with the MCM2-7 subcomplex of the replicative Cdc45-MCM-GINS helicase. The analysis of ESCO2 mutants defective in MCM binding indicates that these interactions are required for proper recruitment of ESCO2 to chromatin, cohesin acetylation during DNA replication, and centromeric cohesion. We propose that MCM binding enables ESCO2 to travel with replisomes to acetylate cohesive cohesin complexes in the vicinity of replication forks so that these complexes can be protected from precocious release by WAPL Our results also indicate that ESCO1 and ESCO2 have distinct functions in maintaining cohesion between chromosome arms and centromeres, respectively.

NASP antagonize chromatin accessibility through maintaining histone H3K9me1 in hepatocellular carcinoma

The regulation of histone deposits mediated by multi-chaperone complexes under physiological conditions remains to be further investigated. Here, we studied the function of nuclear autoantigenic sperm protein (NASP) in the regulation of liver cancer. We found that NASP levels in liver tumors were generally higher than in normal liver tissues and NASP down-regulation inhibited liver cancer cells from forming tumors. We further analyzed cellular responses and epigenetic mechanisms of the histone H3-H4 shortage induced by NASP knockdown in liver cancer cells. The results showed that the major effects of NASP knockdown were globally enhanced chromatin accessibility, which facilitates transcription release, and failure of replication initiation. Furthermore, we demonstrated that NASP depletion led to a global decrease of histone H3K9me1 modification associated with newly H3 processing, which occurred directly at the promoters of up-regulated anti-tumor genes BACH2 and RunX1T1. This also resulted in a synergistic effect on enhanced apoptosis with Myc and p53 decreases. Overall, our work provides new insights into the roles of NASP in tumorigenesis and cancer prevention.

FANCD2 Facilitates Replication through Common Fragile Sites

Common fragile sites (CFSs) are genomic regions that are unstable under conditions of replicative stress. Although the characteristics of CFSs that render them vulnerable to stress are associated mainly with replication, the cellular pathways that protect CFSs during replication remain unclear. Here, we identify and describe a role for FANCD2 as a trans-acting facilitator of CFS replication, in the absence of exogenous replicative stress. In the absence of FANCD2, replication forks stall within the AT-rich fragility core of CFS, leading to dormant origin activation. Furthermore, FANCD2 deficiency is associated with DNA:RNA hybrid formation at CFS-FRA16D, and inhibition of DNA:RNA hybrid formation suppresses replication perturbation. In addition, we also found that FANCD2 reduces the number of potential sites of replication initiation. Our data demonstrate that FANCD2 protein is required to ensure efficient CFS replication and provide mechanistic insight into how FANCD2 regulates CFS stability.

A journey through the microscopic ages of DNA replication

Scientific discoveries and technological advancements are inseparable but not always take place in a coherent chronological manner. In the next, we will provide a seemingly unconnected and serendipitous series of scientific facts that, in the whole, converged to unveil DNA and its duplication. We will not cover here the many and fundamental contributions from microbial genetics and in vitro biochemistry. Rather, in this journey, we will emphasize the interplay between microscopy development culminating on super resolution fluorescence microscopy (i.e., nanoscopy) and digital image analysis and its impact on our understanding of DNA duplication. We will interlace the journey with landmark concepts and experiments that have brought the cellular DNA replication field to its present state.

De novo identification of replication-timing domains in the human genome by deep learning

Motivation: The de novo identification of the initiation and termination zones—regions that replicate earlier or later than their upstream and downstream neighbours, respectively—remains a key challenge in DNA replication.

Results: Building on advances in deep learning, we developed a novel hybrid architecture combining a pre-trained, deep neural network and a hidden Markov model (DNN-HMM) for the de novo identification of replication domains using replication timing profiles. Our results demonstrate that DNN-HMM can significantly outperform strong, discriminatively trained Gaussian mixture model–HMM (GMM-HMM) systems and other six reported methods that can be applied to this challenge. We applied our trained DNN-HMM to identify distinct replication domain types, namely the early replication domain (ERD), the down transition zone (DTZ), the late replication domain (LRD) and the up transition zone (UTZ), using newly replicated DNA sequencing (Repli-Seq) data across 15 human cells. A subsequent integrative analysis revealed that these replication domains harbour unique genomic and epigenetic patterns, transcriptional activity and higher-order chromosomal structure. Our findings support the ‘replication-domain’ model, which states (1) that ERDs and LRDs, connected by UTZs and DTZs, are spatially compartmentalized structural and functional units of higher-order chromosomal structure, (2) that the adjacent DTZ-UTZ pairs form chromatin loops and (3) that intra-interactions within ERDs and LRDs tend to be short-range and long-range, respectively. Our model reveals an important chromatin organizational principle of the human genome and represents a critical step towards understanding the mechanisms regulating replication timing.

Availability and implementation: Our DNN-HMM method and three additional algorithms can be freely accessed at https://github.com/wenjiegroup/DNN-HMM. The replication domain regions identified in this study are available in GEO under the accession ID GSE53984.

Contact:shuwj@bmi.ac.cn or boxc@bmi.ac.cn

Supplementary information:Supplementary data are available at Bioinformatics online.

The hunt for origins of DNA replication in multicellular eukaryotes

Origins of DNA replication (ORIs) occur at defined regions in the genome. Although DNA sequence defines the position of ORIs in budding yeast, the factors for ORI specification remain elusive in metazoa. Several methods have been used recently to map ORIs in metazoan genomes with the hope that features for ORI specification might emerge. These methods are reviewed here with analysis of their advantages and shortcomings. The various factors that may influence ORI selection for initiation of DNA replication are discussed.

Ubiquitous human 'master' origins of replication are encoded in the DNA sequence via a local enrichment in nucleosome excluding energy barriers

As the elementary building block of eukaryotic chromatin, the nucleosome is at the heart of the compromise between the necessity of compacting DNA in the cell nucleus and the required accessibility to regulatory proteins. The recent availability of genome-wide experimental maps of nucleosome positions for many different organisms and cell types has provided an unprecedented opportunity to elucidate to what extent the DNA sequence conditions the primary structure of chromatin and in turn participates in the chromatin-mediated regulation of nuclear functions, such as gene expression and DNA replication. In this study, we use in vivo and in vitro genome-wide nucleosome occupancy data together with the set of nucleosome-free regions (NFRs) predicted by a physical model of nucleosome formation based on sequence-dependent bending properties of the DNA double-helix, to investigate the role of intrinsic nucleosome occupancy in the regulation of the replication spatio-temporal programme in human. We focus our analysis on the so-called replication U/N-domains that were shown to cover about half of the human genome in the germline (skew-N domains) as well as in embryonic stem cells, somatic and HeLa cells (mean replication timing U-domains). The 'master' origins of replication (MaOris) that border these megabase-sized U/N-domains were found to be specified by a few hundred kb wide regions that are hyper-sensitive to DNase I cleavage, hypomethylated, and enriched in epigenetic marks involved in transcription regulation, the hallmarks of localized open chromatin structures. Here we show that replication U/N-domain borders that are conserved in all considered cell lines have an environment highly enriched in nucleosome-excluding-energy barriers, suggesting that these ubiquitous MaOris have been selected during evolution. In contrast, MaOris that are cell-type-specific are mainly regulated epigenetically and are no longer favoured by a local abundance of intrinsic NFRs encoded in the DNA sequence. At the smaller few hundred bp scale of gene promoters, CpG-rich promoters of housekeeping genes found nearby ubiquitous MaOris as well as CpG-poor promoters of tissue-specific genes found nearby cell-type-specific MaOris, both correspond to in vivo NFRs that are not coded as nucleosome-excluding-energy barriers. Whereas the former promoters are likely to correspond to high occupancy transcription factor binding regions, the latter are an illustration that gene regulation in human is typically cell-type-specific.

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.

The role of transcription in the activation of a Drosophila amplification origin

The mechanisms that underlie metazoan DNA replication initiation, especially the connection between transcription and replication origin activation, are not well understood. To probe the role of transcription in origin activation, we exploited a specific replication origin in Drosophila melanogaster follicle cells, ori62, which coincides with the yellow-g2 transcription unit and exhibits transcription-dependent origin firing. Within a 10-kb genomic fragment that contains ori62 and is sufficient for amplification, RNA-sequencing analysis revealed that all detected RNAs mapped solely to the yellow-g2 gene. To determine whether transcription is required in cis for ori62 firing, we generated a set of tagged yellow-g2 transgenes in which we could prevent local transcription across ori62 by deletions in the yellow-g2 promoter. Surprisingly, inhibition of yellow-g2 transcription by promoter deletions did not affect ori62 firing. Our results reveal that transcription in cis is not required for ori62 firing, raising the possibility that a trans-acting factor is required specifically for the activation of ori62. This finding illustrates that a diversity of mechanisms can be used in the regulation of metazoan DNA replication initiation.

Targeted tandem duplication of a large chromosomal segment in Aspergillus oryzae

We describe here the first successful construction of a targeted tandem duplication of a large chromosomal segment in Aspergillus oryzae. The targeted tandem chromosomal duplication was achieved by using strains that had a 5'-deleted pyrG upstream of the region targeted for tandem chromosomal duplication and a 3'-deleted pyrG downstream of the target region. Consequently,strains bearing a 210-kb targeted tandem chromosomal duplication near the centromeric region of chromosome 8 and strains bearing a targeted tandem chromosomal duplication of a 700-kb region of chromosome 2 were successfully constructed. The strains bearing the tandem chromosomal duplication were efficiently obtained from the regenerated protoplast of the parental strains. However, the generation of the chromosomal duplication did not depend on the introduction of double-stranded breaks(DSBs) by I-SceI. The chromosomal duplications of these strains were stably maintained after five generations of culture under nonselective conditions. The strains bearing the tandem chromosomal duplication in the 700-kb region of chromosome 2 showed highly increased protease activity in solid-state culture, indicating that the duplication of large chromosomal segments could be a useful new breeding technology and gene analysis method.

Stable DNA methylation boundaries and expanded trinucleotide repeats: role of DNA insertions

The human genome segment upstream of the FMR1 (fragile X mental retardation 1) gene (Xq27.3) contains several genetic signals, among them is a DNA methylation boundary that is located 65-70 CpGs upstream of the CGG repeat. In fragile X syndrome (FXS), the boundary is lost, and the promoter is inactivated by methylation spreading. Here we document boundary stability in spite of critical expansions of the CGG trinucleotide repeat in male or female premutation carriers and in high functioning males (HFMs). HFMs carry a full CGG repeat expansion but exhibit an unmethylated promoter and lack the FXS phenotype. The boundary is also stable in Turner (45, X) females. A CTCF-binding site is located slightly upstream of the methylation boundary and carries a unique G-to-A polymorphism (single nucleotide polymorphism), which occurs 3.6 times more frequently in genomes with CGG expansions. The increased frequency of this single nucleotide polymorphism might have functional significance. In CGG expansions, the CTCF region does not harbor additional mutations. In FXS individuals and often in cells transgenomic for EBV (Epstein Barr Virus) DNA or for the telomerase gene, the large number of normally methylated CpGs in the far-upstream region of the boundary is decreased about 4-fold. A methylation boundary is also present in the human genome segment upstream of the HTT (huntingtin) promoter (4p16.3) and is stable both in normal and Huntington disease chromosomes. Hence, the vicinity of an expanded repeat does not per se compromise methylation boundaries. Methylation boundaries exert an important function as promoter safeguards.

The spatiotemporal program of DNA replication is associated with specific combinations of chromatin marks in human cells

The duplication of mammalian genomes is under the control of a spatiotemporal program that orchestrates the positioning and the timing of firing of replication origins. The molecular mechanisms coordinating the activation of about predicted origins remain poorly understood, partly due to the intrinsic rarity of replication bubbles, making it difficult to purify short nascent strands (SNS). The precise identification of origins based on the high-throughput sequencing of SNS constitutes a new methodological challenge. We propose a new statistical method with a controlled resolution, adapted to the detection of replication origins from SNS data. We detected an average of 80,000 replication origins in different cell lines. To evaluate the consistency between different protocols, we compared SNS detections with bubble trapping detections. This comparison demonstrated a good agreement between genome-wide methods, with 65% of SNS-detected origins validated by bubble trapping, and 44% of bubble trapping origins validated by SNS origins, when compared at the same resolution. We investigated the interplay between the spatial and the temporal programs of replication at fine scales. We show that most of the origins detected in regions replicated in early S phase are shared by all the cell lines investigated whereas cell-type-specific origins tend to be replicated in late S phase. We shed a new light on the key role of CpG islands, by showing that 80% of the origins associated with CGIs are constitutive. Our results further show that at least 76% of CGIs are origins of replication. The analysis of associations with chromatin marks at different timing of cell division revealed new potential epigenetic regulators driving the spatiotemporal activity of replication origins. We highlight the potential role of H4K20me1 and H3K27me3, the coupling of which is correlated with increased efficiency of replication origins, clearly identifying those marks as potential key regulators of replication origins.

Replication is the mechanism by which genomes are duplicated into two exact copies. Genomic stability is under the control of a spatiotemporal program that orchestrates both the positioning and the timing of firing of about 50,000 replication starting points, also called replication origins. Replication bubbles found at origins have been very difficult to map due to their short lifespan. Moreover, with the flood of data characterizing new sequencing technologies, the precise statistical analysis of replication data has become an additional challenge. We propose a new method to map replication origins on the human genome, and we assess the reliability of our finding using experimental validation and comparison with origins maps obtained by bubble trapping. This fine mapping then allowed us to identify potential regulators of the replication dynamics. Our study highlights the key role of CpG Islands and identifies new potential epigenetic regulators (methylation of lysine 4 on histone H4, and tri-methylation of lysine 27 on histone H3) whose coupling is correlated with an increase in the efficiency of replication origins, suggesting those marks as potential key regulators of replication. Overall, our study defines new potentially important pathways that might regulate the sequential firing of origins during genome duplication.

Human genome replication proceeds through four chromatin states

Advances in genomic studies have led to significant progress in understanding the epigenetically controlled interplay between chromatin structure and nuclear functions. Epigenetic modifications were shown to play a key role in transcription regulation and genome activity during development and differentiation or in response to the environment. Paradoxically, the molecular mechanisms that regulate the initiation and the maintenance of the spatio-temporal replication program in higher eukaryotes, and in particular their links to epigenetic modifications, still remain elusive. By integrative analysis of the genome-wide distributions of thirteen epigenetic marks in the human cell line K562, at the 100 kb resolution of corresponding mean replication timing (MRT) data, we identify four major groups of chromatin marks with shared features. These states have different MRT, namely from early to late replicating, replication proceeds though a transcriptionally active euchromatin state (C1), a repressive type of chromatin (C2) associated with polycomb complexes, a silent state (C3) not enriched in any available marks, and a gene poor HP1-associated heterochromatin state (C4). When mapping these chromatin states inside the megabase-sized U-domains (U-shaped MRT profile) covering about 50% of the human genome, we reveal that the associated replication fork polarity gradient corresponds to a directional path across the four chromatin states, from C1 at U-domains borders followed by C2, C3 and C4 at centers. Analysis of the other genome half is consistent with early and late replication loci occurring in separate compartments, the former correspond to gene-rich, high-GC domains of intermingled chromatin states C1 and C2, whereas the latter correspond to gene-poor, low-GC domains of alternating chromatin states C3 and C4 or long C4 domains. This new segmentation sheds a new light on the epigenetic regulation of the spatio-temporal replication program in human and provides a framework for further studies in different cell types, in both health and disease.

Previous studies revealed spatially coherent and biological-meaningful chromatin mark combinations in human cells. Here, we analyze thirteen epigenetic mark maps in the human cell line K562 at 100 kb resolution of MRT data. The complexity of epigenetic data is reduced to four chromatin states that display remarkable similarities with those reported in fly, worm and plants. These states have different MRT: (C1) is transcriptionally active, early replicating, enriched in CTCF; (C2) is Polycomb repressed, mid-S replicating; (C3) lacks of marks and replicates late and (C4) is a late-replicating gene-poor HP1 repressed heterochromatin state. When mapping these states inside the 876 replication U-domains of K562, the replication fork polarity gradient observed in these U-domains comes along with a remarkable epigenetic organization from C1 at U-domain borders to C2, C3 and ultimately C4 at centers. The remaining genome half displays early replicating, gene rich and high GC domains of intermingled C1 and C2 states segregating from late replicating, gene poor and low GC domains of concatenated C3 and/or C4 states. This constitutes the first evidence of epigenetic compartmentalization of the human genome into replication domains likely corresponding to autonomous units in the 3D chromatin architecture.

Multiscale analysis of genome-wide replication timing profiles using a wavelet-based signal-processing algorithm

In this protocol, we describe the use of the LastWave open-source signal-processing command language (http://perso.ens-lyon.fr/benjamin.audit/LastWave/) for analyzing cellular DNA replication timing profiles. LastWave makes use of a multiscale, wavelet-based signal-processing algorithm that is based on a rigorous theoretical analysis linking timing profiles to fundamental features of the cell's DNA replication program, such as the average replication fork polarity and the difference between replication origin density and termination site density. We describe the flow of signal-processing operations to obtain interactive visual analyses of DNA replication timing profiles. We focus on procedures for exploring the space-scale map of apparent replication speeds to detect peaks in the replication timing profiles that represent preferential replication initiation zones, and for delimiting U-shaped domains in the replication timing profile. In comparison with the generally adopted approach that involves genome segmentation into regions of constant timing separated by timing transition regions, the present protocol enables the recognition of more complex patterns of the spatio-temporal replication program and has a broader range of applications. Completing the full procedure should not take more than 1 h, although learning the basics of the program can take a few hours and achieving full proficiency in the use of the software may take days.

Measures of compositional strand bias related to replication machinery and its applications

The compositional asymmetry of complementary bases in nucleotide sequences implies the existence of a mutational or selectional bias in the two strands of the DNA duplex, which is commonly shaped by strand-specific mechanisms in transcription or replication. Such strand bias in genomes, frequently visualized by GC skew graphs, is used for the computational prediction of transcription start sites and replication origins, as well as for comparative evolutionary genomics studies. The use of measures of compositional strand bias in order to quantify the degree of strand asymmetry is crucial, as it is the basis for determining the applicability of compositional analysis and comparing the strength of the mutational bias in different biological machineries in various species. Here, we review the measures of strand bias that have been proposed to date, including the ∆GC skew, the B₁ index, the predictability score of linear discriminant analysis for gene orientation, the signal-to-noise ratio of the oligonucleotide bias, and the GC skew index. These measures have been predominantly designed for and applied to the analysis of replication-related mutational processes in prokaryotes, but we also give research examples in eukaryotes.

Kinetics of DNA replication and telomerase reaction at a single-seeded telomere in human cells

In most cancer cells, telomerase is activated to elongate telomere DNA, thereby ensuring numerous rounds of cell divisions. It is thus important to understand how telomerase and the replication fork react with telomeres in human cells. However, the highly polymorphic and repetitive nature of the nucleotide sequences in human subtelomeric regions hampers the precise analysis of sequential events taking place at telomeres in S phase. Here, we have established HeLa cells harboring a single-seeded telomere abutted by a unique subtelomere DNA sequence, which has enabled us to specifically focus on the seeded telomere. We have also developed a modified chromatin immunoprecipitation (ChIP) method that uses restriction digestion instead of sonication to fragment chromatin DNA (RES-ChIP), and a method for immunoprecipitating 5-bromo-2'-deoxyuridine (BrdU)-labeled single-stranded DNA by incubating DNA with anti-BrdU antibody in the nondenaturing condition. We have shown that DNA replication of the seeded telomere takes place during a relatively narrow time window in S phase, and telomerase synthesizes telomere DNA after the replication. Moreover, we have demonstrated that the telomerase catalytic subunit TERT associates with telomeres before telomere DNA replication. These results provide a temporal and spatial framework for understanding DNA replication and telomerase reaction at human telomeres.

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.

Genome-wide localization of replication factors

Chromatin Immunoprecipitation (ChIP) is a powerful tool for the identification and characterization of protein-DNA interactions in vivo. ChIP has been utilized to study diverse nuclear processes such as transcription regulation, chromatin modification, DNA recombination and DNA replication at specific loci and, more recently, across the entire genome. Advances in genomic approaches, and whole genome sequencing in particular, have made it possible and affordable to comprehensively identify specific protein binding sites throughout the genomes of higher eukaryotes. The dynamic nature of the DNA replication program and the transient occupancy of many replication factors throughout the cell cycle present additional challenges that may not pertain to the mapping of site specific transcription factors. Here we discuss the specific considerations that need to be addressed in the application of ChIP to the genome-wide location analysis of protein factors involved in DNA replication.

Initiation of DNA Replication in the Human Genome

Replication of the human genome relies on the presence of thousands of origins distributed along each of the chromosomes. The activation of these origins occurs in a highly regulated manner to ensure that chromosomes are faithfully duplicated only once during each cell cycle. Failure in this regulation can lead to abnormal cell proliferation, or/and genomic instability, the hallmarks of cancer cells. The mechanisms determining how, when, and where origins are activated remains still a mystery. However recent technological advances have facilitated the study of DNA replication in a genome-wide scale, and have provided a wealth of information on several features of this process. Here we present an overview of the current progress on our understanding of the initiation step of DNA replication in human cells, and its relationship to abnormal cell proliferation.

Developmental control of gene copy number by repression of replication initiation and fork progression

Precise DNA replication is crucial for genome maintenance, yet this process has been inherently difficult to study on a genome-wide level in untransformed differentiated metazoan cells. To determine how metazoan DNA replication can be repressed, we examined regions selectively under-replicated in Drosophila polytene salivary glands, and found they are transcriptionally silent and enriched for the repressive H3K27me3 mark. In the first genome-wide analysis of binding of the origin recognition complex (ORC) in a differentiated metazoan tissue, we find that ORC binding is dramatically reduced within these large domains, suggesting reduced initiation as one mechanism leading to under-replication. Inhibition of replication fork progression by the chromatin protein SUUR is an additional repression mechanism to reduce copy number. Although repressive histone marks are removed when SUUR is mutated and copy number restored, neither transcription nor ORC binding is reinstated. Tethering of the SUUR protein to a specific site is insufficient to block replication, however. These results establish that developmental control of DNA replication, at both the initiation and elongation stages, is a mechanism to change gene copy number during differentiation.

Common fragile sites: mechanisms of instability revisited

Common fragile sites (CFSs) are large chromosomal regions prone to breakage upon replication stress that are considered a driving force of oncogenesis. CFSs were long believed to contain sequences blocking fork progression, thus impeding replication completion and leading to DNA breaks upon chromosome condensation. However, recent studies show that delayed completion of DNA replication instead depends on a regional paucity in initiation events. Because the distribution and the timing of these events are cell type dependent, different chromosomal regions can be committed to fragility in different cell types. These new data reveal the epigenetic nature of CFSs and open the way to a reevaluation of the role played by these sites in the formation of chromosome rearrangements found in tumors from different tissues.

Interaction of HTLV-1 Tax with minichromosome maintenance proteins accelerates the replication timing program

The Tax oncoprotein encoded by the human T-cell leukemia virus type 1 plays a pivotal role in viral persistence and pathogenesis. Human T-cell leukemia virus type 1-infected cells proliferate faster than normal lymphocytes, expand through mitotic division, and accumulate genomic lesions. Here, we show that Tax associates with the minichromosome maintenance MCM2-7 helicase complex and localizes to origins of replication. Tax modulates the spatiotemporal program of origin activation and fires supplementary origins at the onset of S phase. Thereby, Tax increases the DNA replication rate, accelerates S phase progression, but also generates a replicative stress characterized by the presence of genomic lesions. Mechanistically, Tax favors p300 recruitment and histone hyperacetylation at late replication domains, advancing their replication timing in early S phase.

Chromatin remodeler sucrose nonfermenting 2 homolog (SNF2H) is recruited onto DNA replication origins through interaction with Cdc10 protein-dependent transcript 1 (Cdt1) and promotes pre-replication complex formation

From late mitosis to the G(1) phase of the cell cycle, ORC, CDC6, and Cdt1 form the machinery necessary to load MCM2-7 complexes onto DNA. Here, we show that SNF2H, a member of the ATP-dependent chromatin-remodeling complex, is recruited onto DNA replication origins in human cells in a Cdt1-dependent manner and positively regulates MCM loading. SNF2H physically interacted with Cdt1. ChIP assays indicated that SNF2H associates with replication origins specifically during the G(1) phase. Binding of SNF2H at origins was decreased by Cdt1 silencing and, conversely, enhanced by Cdt1 overexpression. Furthermore, SNF2H silencing prevented MCM loading at origins and moderately inhibited S phase progression. Although neither SNF2H overexpression nor SNF2H silencing appeared to impact rereplication induced by Cdt1 overexpression, Cdt1-induced checkpoint activation was inhibited by SNF2H silencing. Collectively, these data suggest that SNF2H may promote MCM loading at DNA replication origins via interaction with Cdt1 in human cells. Because efficient loading of excess MCM complexes is thought to be required for cells to tolerate replication stress, Cdt1- and SNF2H-mediated promotion of MCM loading may be biologically relevant for the regulation of DNA replication.

Analysis of a Drosophila amplicon in follicle cells highlights the diversity of metazoan replication origins

To investigate the properties of metazoan replication origins, recent studies in cell culture have adopted the strategy of identifying origins using genome-wide approaches and assessing correlations with such features as transcription and histone modifications. Drosophila amplicon in follicle cells (DAFCs), genomic regions that undergo repeated rounds of DNA replication to increase DNA copy number, serve as powerful in vivo model replicons. Because there are six DAFCs, compared with thousands of origins activated in the typical S phase, close molecular characterization of all DAFCs is possible. To determine the extent to which the six DAFCs are different or similar, we investigated the developmental and replication properties of the newly identified DAFC-34B. DAFC-34B contains two genes expressed in follicle cells, although the timing and spatial patterns of expression suggest that amplification is not a strategy to promote high expression at this locus. Like the previously characterized DAFC-62D, DAFC-34B displays origin activation at two separate stages of development. However, unlike DAFC-62D, amplification at the later stage is not transcription-dependent. We mapped the DAFC-34B amplification origin to 1 kb by nascent strand analysis and delineated cis requirements for origin activation, finding that a 6-kb region, but not the 1-kb origin alone, is sufficient for amplification. We analyzed the developmental localization of the origin recognition complex (ORC) and the minichromosome maintenance (MCM)2-7 complex, the replicative helicase. Intriguingly, the final round of origin activation at DAFC-34B occurs in the absence of detectable ORC, although MCMs are present, suggesting a new amplification initiation mechanism.

And-1 is required for the stability of histone acetyltransferase Gcn5

Histone acetyltransferases (HATs) play a central role in the modification of chromatin as well as in pathogenesis of a broad set of diseases including cancers. Gcn5 is the first identified transcription-related histone acetyltransferase (HAT) that has been implicated in the regulation of diverse cellular functions. However, how Gcn5 proteins are regulated remains largely unknown. Here we show that And-1 (a HMG domain-containing protein) has remarkable capability to regulate the stability of Gcn5 proteins and thereby histone H3 acetylation. We find that And-1 forms a complex with both histone H3 and Gcn5. Downregulation of And-1 results in Gcn5 degradation, leading to the reduction of H3K9 and H3K56 acetylation. And-1 overexpression stabilizes Gcn5 through protein-protein interactions in vivo. Furthermore, And-1 expression is increased in cancer cells in a manner correlating with increased Gcn5 and H3K9Ac and H3K56Ac. Thus, our data reveal not only a functional link between Gcn5 and And-1 that is essential to regulate Gcn5 protein stability and histone H3 acetylation, but also a potential role of And-1 in cancer.

Maps of open chromatin guide the functional follow-up of genome-wide association signals: application to hematological traits

Turning genetic discoveries identified in genome-wide association (GWA) studies into biological mechanisms is an important challenge in human genetics. Many GWA signals map outside exons, suggesting that the associated variants may lie within regulatory regions. We applied the formaldehyde-assisted isolation of regulatory elements (FAIRE) method in a megakaryocytic and an erythroblastoid cell line to map active regulatory elements at known loci associated with hematological quantitative traits, coronary artery disease, and myocardial infarction. We showed that the two cell types exhibit distinct patterns of open chromatin and that cell-specific open chromatin can guide the finding of functional variants. We identified an open chromatin region at chromosome 7q22.3 in megakaryocytes but not erythroblasts, which harbors the common non-coding sequence variant rs342293 known to be associated with platelet volume and function. Resequencing of this open chromatin region in 643 individuals provided strong evidence that rs342293 is the only putative causative variant in this region. We demonstrated that the C- and G-alleles differentially bind the transcription factor EVI1 affecting PIK3CG gene expression in platelets and macrophages. A protein–protein interaction network including up- and down-regulated genes in Pik3cg knockout mice indicated that PIK3CG is associated with gene pathways with an established role in platelet membrane biogenesis and thrombus formation. Thus, rs342293 is the functional common variant at this locus; to the best of our knowledge this is the first such variant to be elucidated among the known platelet quantitative trait loci (QTLs). Our data suggested a molecular mechanism by which a non-coding GWA index SNP modulates platelet phenotype.

Genome-wide scans have revealed multiple genetic regions underlying complex traits. However, the transition from an initial association signal to identifying the functional DNA change(s) has proved challenging. Many of the DNA changes discovered are located outside protein-coding regions and may exert their effects through gene regulation. We screened genetic regions associated with hematological traits in erythroblasts (red blood cells) and megakaryocytes (platelet-producing cells) and mapped sites of open chromatin, which harbor active gene regulatory elements. We investigated a DNA sequence change located within a site of open chromatin at chromosome 7 in megakaryocytes, but not erythroblasts, known to be associated with platelet volume. We showed that this DNA change is functional due to alteration of the binding site of a transcription factor, which regulates the expression of a gene that affects platelet characteristics. Mice lacking this gene revealed significant differences in expression of several important platelet genes compared to wild-type mice. The approach described here can be applied in different cell types to functionally follow-up association signals with many other biological traits by identification of the causative base change and how it affects gene function, thus paving the way to clinical benefit.

Preferential localization of human origins of DNA replication at the 5'-ends of expressed genes and at evolutionarily conserved DNA sequences

BACKGROUND

Replication of mammalian genomes requires the activation of thousands of origins which are both spatially and temporally regulated by as yet unknown mechanisms. At the most fundamental level, our knowledge about the distribution pattern of origins in each of the chromosomes, among different cell types, and whether the physiological state of the cells alters this distribution is at present very limited.

METHODOLOGY/PRINCIPAL FINDINGS

We have used standard λ-exonuclease resistant nascent DNA preparations in the size range of 0.7-1.5 kb obtained from the breast cancer cell line MCF-7 hybridized to a custom tiling array containing 50-60 nt probes evenly distributed among genic and non-genic regions covering about 1% of the human genome. A similar DNA preparation was used for high-throughput DNA sequencing. Array experiments were also performed with DNA obtained from BT-474 and H520 cell lines. By determining the sites showing nascent DNA enrichment, we have localized several thousand origins of DNA replication. Our major findings are: (a) both array and DNA sequencing assay methods produced essentially the same origin distribution profile; (b) origin distribution is largely conserved (>70%) in all cell lines tested; (c) origins are enriched at the 5'ends of expressed genes and at evolutionarily conserved intergenic sequences; and (d) ChIP on chip experiments in MCF-7 showed an enrichment of H3K4Me3 and RNA Polymerase II chromatin binding sites at origins of DNA replication.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that the program for origin activation is largely conserved among different cell types. Also, our work supports recent studies connecting transcription initiation with replication, and in addition suggests that evolutionarily conserved intergenic sequences have the potential to participate in origin selection. Overall, our observations suggest that replication origin selection is a stochastic process significantly dependent upon local accessibility to replication factors.

Tiling array data analysis: a multiscale approach using wavelets

BACKGROUND

Tiling array data is hard to interpret due to noise. The wavelet transformation is a widely used technique in signal processing for elucidating the true signal from noisy data. Consequently, we attempted to denoise representative tiling array datasets for ChIP-chip experiments using wavelets. In doing this, we used specific wavelet basis functions, Coiflets, since their triangular shape closely resembles the expected profiles of true ChIP-chip peaks.

RESULTS

In our wavelet-transformed data, we observed that noise tends to be confined to small scales while the useful signal-of-interest spans multiple large scales. We were also able to show that wavelet coefficients due to non-specific cross-hybridization follow a log-normal distribution, and we used this fact in developing a thresholding procedure. In particular, wavelets allow one to set an unambiguous, absolute threshold, which has been hard to define in ChIP-chip experiments. One can set this threshold by requiring a similar confidence level at different length-scales of the transformed signal. We applied our algorithm to a number of representative ChIP-chip data sets, including those of Pol II and histone modifications, which have a diverse distribution of length-scales of biochemical activity, including some broad peaks.

CONCLUSIONS

Finally, we benchmarked our method in comparison to other approaches for scoring ChIP-chip data using spike-ins on the ENCODE Nimblegen tiling array. This comparison demonstrated excellent performance, with wavelets getting the best overall score.

Genome-wide mapping of Arabidopsis thaliana origins of DNA replication and their associated epigenetic marks

Genome integrity requires faithful chromosome duplication. Origins of replication, the genomic sites at which DNA replication initiates, are scattered throughout the genome. Their mapping at a genomic scale in multicellular organisms has been challenging. In this study we profiled origins in Arabidopsis thaliana by high-throughput sequencing of newly synthesized DNA and identified ~1,500 putative origins genome-wide. This was supported by chromatin immunoprecipitation and microarray (ChIP-chip) experiments to identify ORC1- and CDC6-binding sites. We validated origin activity independently by measuring the abundance of nascent DNA strands. The midpoints of most A. thaliana origin regions are preferentially located within the 5' half of genes, enriched in G+C, histone H2A.Z, H3K4me2, H3K4me3 and H4K5ac, and depleted in H3K4me1 and H3K9me2. Our data help clarify the epigenetic specification of DNA replication origins in A. thaliana and have implications for other eukaryotes.

Large duplications at reciprocal translocation breakpoints that might be the counterpart of large deletions and could arise from stalled replication bubbles

Reciprocal chromosome translocations are often not exactly reciprocal. Most familiar are deletions at the breakpoints, up to megabases in extent. We describe here the opposite phenomenon-duplication of tens or hundreds of kilobases at the breakpoint junction, so that the same sequence is present on both products of a translocation. When the products of the translocation are mapped on the genome, they overlap. We report several of these "overlapping-breakpoint" duplications in breast cancer cell lines HCC1187, HCC1806, and DU4475. These lines also had deletions and essentially balanced translocations. In HCC1187 and HCC1806, we identified five cases of duplication ranging between 46 kb and 200 kb, with the partner chromosome showing deletions between 29 bp and 31 Mb. DU4475 had a duplication of at least 200 kb. Breakpoints were mapped using array painting, i.e., hybridization of chromosomes isolated by flow cytometry to custom oligonucleotide microarrays. Duplications were verified by fluorescent in situ hybridization (FISH), PCR on isolated chromosomes, and cloning of breakpoints. We propose that these duplications are the counterpart of deletions and that they are produced at a replication bubble, comprising two replication forks with the duplicated sequence in between. Both copies of the duplicated sequence would go to one daughter cell, on different products of the translocation, while the other daughter cell would show deletion. These duplications may have been overlooked because they may be missed by FISH and array-CGH and may be interpreted as insertions by paired-end sequencing. Such duplications may therefore be quite frequent.

Bubble-chip analysis of human origin distributions demonstrates on a genomic scale significant clustering into zones and significant association with transcription

We have used a novel bubble-trapping procedure to construct nearly pure and comprehensive human origin libraries from early S- and log-phase HeLa cells, and from log-phase GM06990, a karyotypically normal lymphoblastoid cell line. When hybridized to ENCODE tiling arrays, these libraries illuminated 15.3%, 16.4%, and 21.8% of the genome in the ENCODE regions, respectively. Approximately half of the origin fragments cluster into zones, and their signals are generally higher than those of isolated fragments. Interestingly, initiation events are distributed about equally between genic and intergenic template sequences. While only 13.2% and 14.0% of genes within the ENCODE regions are actually transcribed in HeLa and GM06990 cells, 54.5% and 25.6% of zonal origin fragments overlap transcribed genes, most with activating chromatin marks in their promoters. Our data suggest that cell synchronization activates a significant number of inchoate origins. In addition, HeLa and GM06990 cells activate remarkably different origin populations. Finally, there is only moderate concordance between the log-phase HeLa bubble map and published maps of small nascent strands for this cell line.

H2A.Z demarcates intergenic regions of the plasmodium falciparum epigenome that are dynamically marked by H3K9ac and H3K4me3

Epigenetic regulatory mechanisms and their enzymes are promising targets for malaria therapeutic intervention; however, the epigenetic component of gene expression in P. falciparum is poorly understood. Dynamic or stable association of epigenetic marks with genomic features provides important clues about their function and helps to understand how histone variants/modifications are used for indexing the Plasmodium epigenome. We describe a novel, linear amplification method for next-generation sequencing (NGS) that allows unbiased analysis of the extremely AT-rich Plasmodium genome. We used this method for high resolution, genome-wide analysis of a histone H2A variant, H2A.Z and two histone H3 marks throughout parasite intraerythrocytic development. Unlike in other organisms, H2A.Z is a constant, ubiquitous feature of euchromatic intergenic regions throughout the intraerythrocytic cycle. The almost perfect colocalisation of H2A.Z with H3K9ac and H3K4me3 suggests that these marks are preferentially deposited on H2A.Z-containing nucleosomes. By performing RNA-seq on 8 time-points, we show that acetylation of H3K9 at promoter regions correlates very well with the transcriptional status whereas H3K4me3 appears to have stage-specific regulation, being low at early stages, peaking at trophozoite stage, but does not closely follow changes in gene expression. Our improved NGS library preparation procedure provides a foundation to exploit the malaria epigenome in detail. Furthermore, our findings place H2A.Z at the cradle of P. falciparum epigenetic regulation by stably defining intergenic regions and providing a platform for dynamic assembly of epigenetic and other transcription related complexes.

Plasmodium falciparum is a unicellular pathogen that is responsible for the most severe form of malaria. Similar to other eukaryotic organisms, its genome is organized into chromosomes by proteins called histones. Modification or replacement of these histones has marked effects on the packaging grade of DNA and instructs the recruitment of protein complexes, thereby regulating essential cellular processes such as gene expression and replication. Here we unveil the genome-wide localization of two histone H3 modifications (K9ac/K4me3) and a histone variant, H2A.Z, during development of the parasite in the human red blood cells. We find that all three epigenetic features are predominantly present in intergenic regions of the P. falciparum genome, suggesting an interconnecting role in regulation of gene expression. H2A.Z levels appear to be largely invariable throughout intraerythrocytic development while placement/removal of the histone marks is dynamic with H3K9ac and H3K4me3 being transcription-coupled and stage-specific, respectively. These observations support a model in which H2A.Z-containing nucleosomes serve to demarcate regulatory regions in the parasite's genome and promote transcription initiation by guiding chromatin modifying and transcription initiating complexes. The findings and methodological developments presented in this paper provide a cornerstone for future epigenome research in eukaryotic pathogens and vital information to understand and to interfere with parasite development and survival.

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

Genome-wide mapping of nuclear mitochondrial DNA sequences links DNA replication origins to chromosomal double-strand break formation in Schizosaccharomyces pombe

Chromosomal double-strand breaks (DSBs) threaten genome integrity and repair of these lesions is often mutagenic. How and where DSBs are formed is a major question conveniently addressed in simple model organisms like yeast. NUMTs, nuclear DNA sequences of mitochondrial origin, are present in most eukaryotic genomes and probably result from the capture of mitochondrial DNA (mtDNA) fragments into chromosomal breaks. NUMT formation is ongoing and was reported to cause de novo human genetic diseases. Study of NUMTs is likely to contribute to the understanding of naturally occurring chromosomal breaks. We show that Schizosaccharomyces pombe NUMTs are exclusively located in noncoding regions with no preference for gene promoters and, when located into promoters, do not affect gene transcription level. Strikingly, most noncoding regions comprising NUMTs are also associated with a DNA replication origin (ORI). Chromatin immunoprecipitation experiments revealed that chromosomal NUMTs are probably not acting as ORI on their own but that mtDNA insertions occurred directly next to ORIs, suggesting that these loci may be prone to DSB formation. Accordingly, induction of excessive DNA replication origin firing, a phenomenon often associated with human tumor formation, resulted in frequent nucleotide deletion events within ORI3001 subtelomeric chromosomal locus, illustrating a novel aspect of DNA replication-driven genomic instability. How mtDNA is fragmented is another important issue that we addressed by sequencing experimentally induced NUMTs. This highlighted regions of S. pombe mtDNA prone to breaking. Together with an analysis of human NUMTs, we propose that these fragile sites in mtDNA may correspond to replication pause sites.

Molecular analysis of the replication program in unicellular model organisms

Eukaryotes have long been reported to show temporal programs of replication, different portions of the genome being replicated at different times in S phase, with the added possibility of developmentally regulated changes in this pattern depending on species and cell type. Unicellular model organisms, primarily the budding yeast Saccharomyces cerevisiae, have been central to our current understanding of the mechanisms underlying the regulation of replication origins and the temporal program of replication in particular. But what exactly is a temporal program of replication, and how might it arise? In this article, we explore this question, drawing again on the wealth of experimental information in unicellular model organisms.

Genome-wide analysis of the replication program in mammals

Microarray technology has facilitated the research of eukaryotic DNA replication on a genome-wide scale. Recent studies have shed light on the association between time of replication and chromosome structure, on the organization principles of the replication program, and on the correlation between replication timing and transcription. In this review, we summarize various genomic measurement approaches and the biological insights achieved through applying them in the study of the mammalian replication program.

Why are we where we are? Understanding replication origins and initiation sites in eukaryotes using ChIP-approaches

DNA replication initiates from origins of replication following a strict sequential activation programme and a conserved temporal order of activation. The number of replication initiation sites varies between species, according to the complexity of the genomes, with an average spacing of 100,000 bp. In contrast to yeast genomes, the location and definition of origins in mammalian genomes has been elusive. Historically, mammalian replication initiation sites have been mapped in situ by systematically searching specific genomic loci for sites that preferentially initiated DNA replication, potential origins by start-site mapping and autonomously replicating sequence experiments, and potential ORC and pre-replicative complex (pre-RC) sites by chromatin immunoprecipitation (ChIP) using antibodies for pre-RC proteins. In the past decade, ChIP has become an important method for analyzing protein/DNA interactions. Classically, ChIP is combined with Southern blotting or PCR. Recently, whole genome-ChIP methods have been very successful in unicellular eukaryotes to understand molecular mechanisms coordinating replication initiation and its flexibility in response to environmental changes. However, in mammalian systems, ChIP with pre-RC antibodies has often been challenging and genome-wide studies are scarce. In this review, we will appraise the progress that has been made in understanding replication origin organization using immunoprecipitation of the ORC and Mcm2-7 complexes. A special focus will be on the advantages and disadvantages of genome-wide ChIP-technologies and their potential impact on understanding metazoan replicators.

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.

Organization of DNA replication

The discovery of the DNA double helix structure half a century ago immediately suggested a mechanism for its duplication by semi-conservative copying of the nucleotide sequence into two DNA daughter strands. Shortly after, a second fundamental step toward the elucidation of the mechanism of DNA replication was taken with the isolation of the first enzyme able to polymerize DNA from a template. In the subsequent years, the basic mechanism of DNA replication and its enzymatic machinery components were elucidated, mostly through genetic approaches and in vitro biochemistry. Most recently, the spatial and temporal organization of the DNA replication process in vivo within the context of chromatin and inside the intact cell are finally beginning to be elucidated. On the one hand, recent advances in genome-wide high throughput techniques are providing a new wave of information on the progression of genome replication at high spatial resolution. On the other hand, novel super-resolution microscopy techniques are just starting to give us the first glimpses of how DNA replication is organized within the context of single intact cells with high spatial resolution. The integration of these data with time lapse microscopy analysis will give us the ability to film and dissect the replication of the genome in situ and in real time.

Impaired replication dynamics at the FRA3B common fragile site

Chromosomal common fragile sites (CFSs) are genetically unstable regions of the genome that are induced by conditions that impair DNA replication. In this report, we show that treatment with the DNA polymerase inhibitor, aphidicolin (APH), slows the replication rate throughout S phase. To investigate the unusual sensitivity of CFSs to APH-induced replication stress, we examined replication dynamics within a 50 kb region of the most frequently expressed CFS, FRA3B. We mapped four origins of replication, ori 1-4, using two independent methods. In untreated cells, we detected significantly less newly replicated DNA at FRA3B ori 1-3, as compared with three control origins located within non-fragile regions (NCFSs). In APH-treated cells, all FRA3B and control origins tested were active; however, there was a significant increase of nascent strand DNA at the control origins and, to a lesser extent, at the FRA3B ori 1-3. On the basis of these observations and the theoretical modeling of the nascent strand abundance assay developed in this study, we hypothesize that CFS origins may be less efficient, and that APH treatment slows replication fork movement near these origins to a greater extent, resulting in impaired DNA replication and, ultimately, leading to the genetic instability characteristic of CFSs.

HBO1 histone acetylase activity is essential for DNA replication licensing and inhibited by Geminin

HBO1, an H4-specific histone acetylase, is a coactivator of the DNA replication licensing factor Cdt1. HBO1 acetylase activity is required for licensing, because a histone acetylase (HAT)-defective mutant of HBO1 bound at origins is unable to load the MCM complex. H4 acetylation at origins is cell-cycle regulated, with maximal activity at the G1/S transition, and coexpression of HBO1 and Jade-1 increases histone acetylation and MCM complex loading. Overexpression of the Set8 histone H4 tail-binding domain specifically inhibits MCM loading, suggesting that histones are a physiologically relevant target for licensing. Lastly, Geminin inhibits HBO1 acetylase activity in the context of a Cdt1-HBO1 complex, and it associates with origins and inhibits H4 acetylation and licensing in vivo. Thus, H4 acetylation at origins by HBO1 is critical for replication licensing by Cdt1, and negative regulation of licensing by Geminin is likely to involve inhibition of HBO1 histone acetylase activity.

Drosophila ORC localizes to open chromatin and marks sites of cohesin complex loading

The origin recognition complex (ORC) is an essential DNA replication initiation factor conserved in all eukaryotes. In Saccharomyces cerevisiae, ORC binds to specific DNA elements; however, in higher eukaryotes, ORC exhibits little sequence specificity in vitro or in vivo. We investigated the genome-wide distribution of ORC in Drosophila and found that ORC localizes to specific chromosomal locations in the absence of any discernible simple motif. Although no clear sequence motif emerged, we were able to use machine learning approaches to accurately discriminate between ORC-associated sequences and ORC-free sequences based solely on primary sequence. The complex sequence features that define ORC binding sites are highly correlated with nucleosome positioning signals and likely represent a preferred nucleosomal landscape for ORC association. Open chromatin appears to be the underlying feature that is deterministic for ORC binding. ORC-associated sequences are enriched for the histone variant, H3.3, often at transcription start sites, and depleted for bulk nucleosomes. The density of ORC binding along the chromosome is reflected in the time at which a sequence replicates, with early replicating sequences having a high density of ORC binding. Finally, we found a high concordance between sites of ORC binding and cohesin loading, suggesting that, in addition to DNA replication, ORC may be required for the loading of cohesin on DNA in Drosophila.

Genomic study of replication initiation in human chromosomes reveals the influence of transcription regulation and chromatin structure on origin selection

DNA replication in metazoans initiates from multiple chromosomal loci called origins. This study identifies 150 new origins of replication that were confirmed by two methods of nascent strand purification. We discern the role of transcription initiation and regulation, as well as chromatin signatures in determining origin selection in human genome.

DNA replication in metazoans initiates from multiple chromosomal loci called origins. Currently, there are two methods to purify origin-centered nascent strands: lambda exonuclease digestion and anti-bromodeoxyuridine immunoprecipitation. Because both methods have unique strengths and limitations, we purified nascent strands by both methods, hybridized them independently to tiling arrays (1% genome) and compared the data to have an accurate view of genome-wide origin distribution. By this criterion, we identified 150 new origins that were reproducible across the methods. Examination of a subset of these origins by chromatin immunoprecipitation against origin recognition complex (ORC) subunits 2 and 3 showed 93% of initiation peaks to localize at/within 1 kb of ORC binding sites. Correlation of origins with functional elements of the genome revealed origin activity to be significantly enriched around transcription start sites (TSSs). Consistent with proximity to TSSs, we found a third of initiation events to occur at or near the RNA polymerase II binding sites. Interestingly, ∼50% of the early origin activity was localized within 5 kb of transcription regulatory factor binding region clusters. The chromatin signatures around the origins were enriched in H3K4-(di- and tri)-methylation and H3 acetylation modifications on histones. Affinity of origins for open chromatin was also reiterated by their proximity to DNAse I-hypersensitive sites. Replication initiation peaks were AT rich, and >50% of the origins mapped to evolutionarily conserved regions of the genome. In summary, these findings indicate that replication initiation is influenced by transcription initiation and regulation as well as chromatin structure.

A distinct DNA-methylation boundary in the 5'- upstream sequence of the FMR1 promoter binds nuclear proteins and is lost in fragile X syndrome

We have discovered a distinct DNA-methylation boundary at a site between 650 and 800 nucleotides upstream of the CGG repeat in the first exon of the human FMR1 gene. This boundary, identified by bisulfite sequencing, is present in all human cell lines and cell types, irrespective of age, gender, and developmental stage. The same boundary is found also in different mouse tissues, although sequence homology between human and mouse in this region is only 46.7%. This boundary sequence, in both the unmethylated and the CpG-methylated modes, binds specifically to nuclear proteins from human cells. We interpret this boundary as carrying a specific chromatin structure that delineates a hypermethylated area in the genome from the unmethylated FMR1 promoter and protecting it from the spreading of DNA methylation. In individuals with the fragile X syndrome (FRAXA), the methylation boundary is lost; methylation has penetrated into the FMR1 promoter and inactivated the FMR1 gene. In one FRAXA genome, the upstream terminus of the methylation boundary region exhibits decreased methylation as compared to that of healthy individuals. This finding suggests changes in nucleotide sequence and chromatin structure in the boundary region of this FRAXA individual. In the completely de novo methylated FMR1 promoter, there are isolated unmethylated CpG dinucleotides that are, however, not found when the FMR1 promoter and upstream sequences are methylated in vitro with the bacterial M-SssI DNA methyltransferase. They may arise during de novo methylation only in DNA that is organized in chromatin and be due to the binding of specific proteins.

Y RNA functions at the initiation step of mammalian chromosomal DNA replication

Non-coding Y RNAs have recently been identified as essential novel factors for chromosomal DNA replication in mammalian cell nuclei, but mechanistic details of their function have not been defined. Here, we identify the execution point for Y RNA function during chromosomal DNA replication in a mammalian cell-free system. We determined the effect of degradation of Y3 RNA on replication origin activation and on fork progression rates at single-molecule resolution by DNA combing and nascent-strand analysis. Degradation of Y3 RNA inhibits the establishment of new DNA replication forks at the G1- to S-phase transition and during S phase. This inhibition is negated by addition of exogenous Y1 RNA. By contrast, progression rates of DNA replication forks are not affected by degradation of Y3 RNA or supplementation with exogenous Y1 RNA. These data indicate that Y RNAs are required for the establishment, but not for the elongation, of chromosomal DNA replication forks in mammalian cell nuclei. We conclude that the execution point for non-coding Y RNA function is the activation of chromosomal DNA replication origins.