Genome-wide analysis of replication timing in mammalian cells: Troubleshooting problems encountered when comparing different cell types

DNA is replicated in a defined temporal order that is developmentally regulated and constitutes a unique and stable fingerprint of a given cell type. Recently, we developed a robust assay to profile replication timing genome wide that can be applied to essentially any proliferating cell population. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU). The cells are sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated (BrdU IP), amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray (or sequenced). Since the basic steps of this protocol have been detailed elsewhere, here we focus on problems encountered when adapting this protocol to different cell types or tissue sources and modifications that have been successfully applied to troubleshoot these problems. There is an increasing demand for such studies to address how replication is regulated during development, its relationship to chromatin architecture and other chromosome functions, and the relevance of cell culture models to regulation in the native organismal niche.     

Cellular senescence induces replication stress with almost no affect on DNA replication timing

Organismal aging entails a gradual decline of normal physiological functions and a major contributor to this decline is withdrawal of the cell cycle, known as senescence. Senescence can result from telomere diminution leading to a finite number of population doublings, known as replicative senescence (RS), or from oncogene overexpression, as a protective mechanism against cancer. Senescence is associated with large-scale chromatin re-organization and changes in gene expression. Replication stress is a complex phenomenon, defined as the slowing or stalling of replication fork progression and/or DNA synthesis, which has serious implications for genome stability, and consequently in human diseases. Aberrant replication fork structures activate the replication stress response leading to the activation of dormant origins, which is thought to be a safeguard mechanism to complete DNA replication on time. However, the relationship between replicative stress and the changes in the spatiotemporal program of DNA replication in senescence progression remains unclear.

Here, we studied the DNA replication program during senescence progression in proliferative and pre-senescent cells from donors of various ages by single DNA fiber combing of replicated DNA, origin mapping by sequencing short nascent strands and genome-wide profiling of replication timing (TRT).

We demonstrate that, progression into RS leads to reduced replication fork rates and activation of dormant origins, which are the hallmarks of replication stress. However, with the exception of a delay in RT of the CREB5 gene in all pre-senescent cells, RT was globally unaffected by replication stress during entry into either oncogene-induced or RS. Consequently, we conclude that RT alterations associated with physiological and accelerated aging, do not result from senescence progression. Our results clarify the interplay between senescence, aging and replication programs and demonstrate that RT is largely resistant to replication stress.     

Replication timing analysis in polyploid cells reveals Rif1 uses multiple mechanisms to promote underreplication in Drosophila

Regulation of DNA replication and copy number is necessary to promote genome stability and maintain cell and tissue function. DNA replication is regulated temporally in a process known as replication timing (RT). Rap1-interacting factor 1 (Rif1) is a key regulator of RT and has a critical function in copy number control in polyploid cells. Previously, we demonstrated that Rif1 functions with SUUR to inhibit replication fork progression and promote underreplication (UR) of specific genomic regions. How Rif1-dependent control of RT factors into its ability to promote UR is unknown. By applying a computational approach to measure RT in Drosophila polyploid cells, we show that SUUR and Rif1 have differential roles in controlling UR and RT. Our findings reveal that Rif1 acts to promote late replication, which is necessary for SUUR-dependent underreplication. Our work provides new insight into the process of UR and its links to RT.     

Spatio-temporal organization of DNA replication in murine embryonic stem, primary, and immortalized cells

The extent to which chromosomal domains are reorganized within the nucleus during differentiation is central to our understanding of how cells become committed to specific developmental lineages. Spatio-temporal patterns of DNA replication are a reflection of this organization. Here, we demonstrate that the temporal order and relative duration of these replication patterns during S-phase are similar in murine pluripotent embryonic stem (ES) cells, primary adult myoblasts, and an immortalized fibroblast line. The observed patterns were independent of fixation and denaturation techniques. Importantly, the same patterns were detected when fluorescent nucleotides were introduced into living cells, demonstrating their physiological relevance. These data suggest that heritable gene silencing during commitment to specific cell lineages is not mediated by global changes in the sub-nuclear organization and replication timing of chromosome domains.     

Determination of DNA replication kinetics in synchronized human cells using a PCR-based assay.

Studies on the temporal order of DNA replication are difficult due to the lack of sensitivity of methods available for replication kinetic analysis. To overcome problems associated with the current techniques, we propose a PCR-based assay to determine the replication time of any single-copy DNA sequence in complex genomes. Human cells labeled with 5-bromodeoxyuridine (BrdU) were flow sorted, according to their DNA content, at different times after synchronous release from the G1/S phase boundary. The selective removal of newly-replicated BrdU-substituted DNA was achieved by UV light irradiation followed by S1 nuclease treatment. The timing of replication of selected DNA sequences (housekeeping, tissue-specific, and non-coding loci) was determined by polymerase chain reaction (PCR) amplification using appropriate primers. DNA sequences localized in inactive replication units allowed amplification whereas those that have replicated will not be amplified by PCR. Using this sensitive and quantitative assay the replication kinetic analysis of a number of different DNA sequences can be performed from a single sorting experiment.     

Regulation of the replication of the murine immunoglobulin heavy chain gene locus: evaluation of the role of the 3' regulatory region.

DNA replication in mammalian cells is a precisely controlled physical and temporal process, likely involving cis-acting elements that control the region(s) from which replication initiates. In B cells, previous studies showed replication timing to be early throughout the immunoglobulin heavy chain (Igh) locus. The implication from replication timing studies in the B-cell line MPC11 was that early replication of the Igh locus was regulated by sequences downstream of the C alpha gene. A potential candidate for these replication control sequences was the 3' regulatory region of the Igh locus. Our results demonstrate, however, that the Igh locus maintains early replication in a B-cell line in which the 3' regulatory region has been deleted from one allele, thus indicating that replication timing of the locus is independent of this region. In non-B cells (murine erythroleukemia cells [MEL]), previous studies of segments within the mouse Igh locus demonstrated that DNA replication likely initiated downstream of the Igh gene cluster. Here we use recently cloned DNA to demonstrate that segments located sequentially downstream of the Igh 3' regulatory region continue to replicate progressively earlier in S phase in MEL. Furthermore, analysis by two-dimensional gel electrophoresis indicates that replication forks proceed exclusively in the 3'-to-5' direction through the region 3' of the Igh locus. Extrapolation from these data predicts that initiation of DNA replication occurs in MEL at one or more sites within a 90-kb interval located between 40 and 130 kb downstream of the 3' regulatory region.     

Reversal of DNA methylation with 5-azacytidine alters chromosome replication patterns in human lymphocyte and fibroblast cultures

Prior studies demonstrated that developmental or induced methylation of DNA can inactivate associated gene loci. Such DNA methylation can be reversed and specific genes reactivated by treatment with 5-azacytidine (5- azaC ). The present cytogenetic studies using replication banding methods show that 5- azaC treatment also results in an increase or decrease in replication staining at one or more band locations in human lymphocyte and fibroblast chromosomes. New replication band locations are not formed. These changes in replication staining, which reflect changes in timing of replication, are different between these two tissues. However, in both tissues, the delayed onset of replication in the heterocyclic, inactive X is shortened by 5- azaC . A correlation is thus suggested between the induced temporal change to earlier DNA replication, and induced hypomethylation and gene activation. The temporal effect on chromosome replication in 5- azaC -treated cells depends on the portion of the S-period studied. Toward the beginning of S, early-replication patterns are increased in both lymphocytes and fibroblasts. Toward the end of S, late-replication patterns are increased only in lymphocytes, suggesting a differential effect of 5- azaC in: (1) early-vs. late-S, and (2) lymphocytes vs. fibroblasts. Generally, 5- azaC has its greatest effect on the inactive chromosome regions that are typically late-replicating prior to 5- azaC treatment. These observed changes in replication band staining suggest that DNA methylation may modify regional groups of genes in concert.     

Visualization of DNA replication in the vertebrate model system DT40 using the DNA fiber technique

Maintenance of replication fork stability is of utmost importance for dividing cells to preserve viability and prevent disease. The processes involved not only ensure faithful genome duplication in the face of endogenous and exogenous DNA damage but also prevent genomic instability, a recognized causative factor in tumor development. Here, we describe a simple and cost-effective fluorescence microscopy-based method to visualize DNA replication in the avian B-cell line DT40. This cell line provides a powerful tool to investigate protein function in vivo by reverse genetics in vertebrate cells(1). DNA fiber fluorography in DT40 cells lacking a specific gene allows one to elucidate the function of this gene product in DNA replication and genome stability. Traditional methods to analyze replication fork dynamics in vertebrate cells rely on measuring the overall rate of DNA synthesis in a population of pulse-labeled cells. This is a quantitative approach and does not allow for qualitative analysis of parameters that influence DNA synthesis. In contrast, the rate of movement of active forks can be followed directly when using the DNA fiber technique(2-4). In this approach, nascent DNA is labeled in vivo by incorporation of halogenated nucleotides (Fig 1A). Subsequently, individual fibers are stretched onto a microscope slide, and the labeled DNA replication tracts are stained with specific antibodies and visualized by fluorescence microscopy (Fig 1B). Initiation of replication as well as fork directionality is determined by the consecutive use of two differently modified analogues. Furthermore, the dual-labeling approach allows for quantitative analysis of parameters that influence DNA synthesis during the S-phase, i.e. replication structures such as ongoing and stalled forks, replication origin density as well as fork terminations. Finally, the experimental procedure can be accomplished within a day, and requires only general laboratory equipment and a fluorescence microscope.     

Chromatin assembly at kinetochores is uncoupled from DNA replication

The specification of metazoan centromeres does not depend strictly on centromeric DNA sequences, but also requires epigenetic factors. The mechanistic basis for establishing a centromeric "state" on the DNA remains unclear. In this work, we have directly examined replication timing of the prekinetochore domain of human chromosomes. Kinetochores were labeled by expression of epitope-tagged CENP-A, which stably marks prekinetochore domains in human cells. By immunoprecipitating CENP-A mononucleosomes from synchronized cells pulsed with [(3)H]thymidine we demonstrate that CENP-A-associated DNA is replicated in mid-to-late S phase. Cytological analysis of DNA replication further demonstrated that centromeres replicate asynchronously in parallel with numerous other genomic regions. In contrast, quantitative Western blot analysis demonstrates that CENP-A protein synthesis occurs later, in G2. Quantitative fluorescence microscopy and transient transfection in the presence of aphidicolin, an inhibitor of DNA replication, show that CENP-A can assemble into centromeres in the absence of DNA replication. Thus, unlike most genomic chromatin, histone synthesis and assembly are uncoupled from DNA replication at the kinetochore. Uncoupling DNA replication from CENP-A synthesis suggests that regulated chromatin assembly or remodeling could play a role in epigenetic centromere propagation.     

The replication timing of the amplified dihydrofolate reductase genes in the Chinese hamster ovary cell line CHOC 400

We have examined the timing of replication of the amplified dihydrofolate reductase genes in the methotrexate-resistant Chinese hamster ovary cell line CHOC 400 using two synchronization procedures. DNA replicated in the presence of 5-bromodeoxyuridine was collected from cells of various times during the DNA synthesis phase and the extent of replication for defined sequences was determined by Southern blotting analysis of CsCl density gradient fractions. We report that under these conditions the DHFR gene replicates throughout the course of S phase in a mode similar to the bulk of the replicated genomic DNA. This contrasts with previous data that shows the non-amplified DHFR gene replicates during the first quarter of S phase. Therefore, we conclude that gene amplification alters the replication timing of the DHFR gene in CHOC 400 cells.     

Rif1 regulates the replication timing domains on the human genome

DNA replication is spatially and temporally regulated during S-phase. DNA replication timing is established in early-G1-phase at a point referred to as timing decision point. However, how the genome-wide replication timing domains are established is unknown. Here, we show that Rif1 (Rap1-interacting-factor-1), originally identified as a telomere-binding factor in yeast, is a critical determinant of the replication timing programme in human cells. Depletion of Rif1 results in specific loss of mid-S replication foci profiles, stimulation of initiation events in early-S-phase and changes in long-range replication timing domain structures. Analyses of replication timing show replication of sequences normally replicating early is delayed, whereas that normally replicating late is advanced, suggesting that replication timing regulation is abrogated in the absence of Rif1. Rif1 tightly binds to nuclear-insoluble structures at late-M-to-early-G1 and regulates chromatin-loop sizes. Furthermore, Rif1 colocalizes specifically with the mid-S replication foci. Thus, Rif1 establishes the mid-S replication domains that are restrained from being activated at early-S-phase. Our results indicate that Rif1 plays crucial roles in determining the replication timing domain structures in human cells through regulating higher-order chromatin architecture.