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

Visualizing DNA replication by single-molecule analysis of replicated DNA

Single-molecule analysis of replicated DNA (SMARD) is a unique technique that enables visualization of DNA replication at specific genomic regions at single-molecule resolution. Here, we present a protocol for visualizing DNA replication by SMARD. We describe steps for pulse labeling DNA, followed by isolating and stretching of genomic DNA. We then detail the detection of the replication at chromosomal regions through immunostaining and fluorescence in situ hybridization. Using SMARD, we can visualize replication initiation, progression, termination, and fork stalling. For complete details on the use and execution of this protocol, please refer to Norio et al. (2001) and Gerhardt et al. (2014).1,2.

DNA polymerase zeta contributes to heterochromatin replication to prevent genome instability

The DNA polymerase zeta (Polζ) plays a critical role in bypassing DNA damage. REV3L, the catalytic subunit of Polζ, is also essential in mouse embryonic development and cell proliferation for reasons that remain incompletely understood. In this study, we reveal that REV3L protein interacts with heterochromatin components including repressive histone marks and localizes in pericentromeric regions through direct interaction with HP1 dimer. We demonstrate that Polζ/REV3L ensures progression of replication forks through difficult‐to‐replicate pericentromeric heterochromatin, thereby preventing spontaneous chromosome break formation. We also find that Rev3l‐deficient cells are compromised in the repair of heterochromatin‐associated double‐stranded breaks, eliciting deletions in late‐replicating regions. Lack of REV3L leads to further consequences that may be ascribed to heterochromatin replication and repair‐associated functions of Polζ, with a disruption of the temporal replication program at specific loci. This is correlated with changes in epigenetic landscape and transcriptional control of developmentally regulated genes. These results reveal a new function of Polζ in preventing chromosome instability during replication of heterochromatic regions.

In Vivo Analysis of mtDNA Replication at the Single Molecule Level and with High Resolution

Single molecule analysis of replicating DNA (SMARD) is a powerful methodology that allows in vivo analysis of replicating DNA; identification of origins of replication, assessment of fork directionality, and measurement of replication fork speed. SMARD, which has been extensively used to study replication of nuclear DNA, involves incorporation of thymidine analogs to nascent DNA chains and their subsequent visualization through immune detection. Here, we adapt and fine-tune the SMARD technique to the specifics of human and mouse mitochondrial DNA. The mito-SMARD protocol allows researchers to gain in vivo insight into mitochondrial DNA (mtDNA) replication at the single molecule level and with high resolution.

Pluripotent stem cells with low differentiation potential contain incompletely reprogrammed DNA replication

Reprogrammed pluripotent stem cells (PSCs) are valuable for research and potentially for cell replacement therapy. However, only a fraction of reprogrammed PSCs are developmentally competent. Genomic stability and accurate DNA synthesis are fundamental for cell development and critical for safety. We analyzed whether defects in DNA replication contribute to genomic instability and the diverse differentiation potentials of reprogrammed PSCs. Using a unique single-molecule approach, we visualized DNA replication in isogenic PSCs generated by different reprogramming approaches, either somatic cell nuclear transfer (NT-hESCs) or with defined factors (iPSCs). In PSCs with lower differentiation potential, DNA replication was incompletely reprogrammed, and genomic instability increased during replicative stress. Reprogramming of DNA replication did not correlate with DNA methylation. Instead, fewer replication origins and a higher frequency of DNA breaks in PSCs with incompletely reprogrammed DNA replication were found. Given the impact of error-free DNA synthesis on the genomic integrity and differentiation proficiency of PSCs, analyzing DNA replication may be a useful quality control tool.

The Replication of Frataxin Gene Is Assured by Activation of Dormant Origins in the Presence of a GAA-Repeat Expansion

It is well known that DNA replication affects the stability of several trinucleotide repeats, but whether replication profiles of human loci carrying an expanded repeat differ from those of normal alleles is poorly understood in the endogenous context. We investigated this issue using cell lines from Friedreich's ataxia patients, homozygous for a GAA-repeat expansion in intron 1 of the Frataxin gene. By interphase, FISH we found that in comparison to the normal Frataxin sequence the replication of expanded alleles is slowed or delayed. According to molecular combing, origins never fired within the normal Frataxin allele. In contrast, in mutant alleles dormant origins are recruited within the gene, causing a switch of the prevalent fork direction through the expanded repeat. Furthermore, a global modification of the replication profile, involving origin choice and a differential distribution of unidirectional forks, was observed in the surrounding 850 kb region. These data provide a wide-view of the interplay of events occurring during replication of genes carrying an expanded repeat.

Replication timing and transcriptional control: beyond cause and effect-part III

DNA replication is essential for faithful transmission of genetic information and is intimately tied to chromosome structure and function. Genome duplication occurs in a defined temporal order known as the replication-timing (RT) program, which is regulated during the cell cycle and development in discrete units referred to as replication domains (RDs). RDs correspond to topologically-associating domains (TADs) and are spatio-temporally compartmentalized in the nucleus. While improvements in experimental tools have begun to reveal glimpses of causality, they have also unveiled complex context-dependent relationships that challenge long recognized correlations of RT to chromatin organization and gene regulation. In particular, RDs/TADs that switch RT during development march to the beat of a different drummer.

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.

Epigenetic modifications in human fragile X pluripotent stem cells; Implications in fragile X syndrome modeling

Patients with fragile X syndrome (FXS) exhibit moderate to severe intellectual disabilities. In addition, one-third of FXS patients show characteristics of autism spectrum disorder. FXS is caused by a trinucleotide repeat expansion, which leads to silencing of the fragile X mental retardation (FMR1) gene. The absence of the FMR1 gene product, FMRP, is the reason for the disease symptoms. It has been suggested that repeat instability and transcription of the FMR1 gene occur during early embryonic development, while after cell differentiation repeats become stable and the FMR1 gene is silent. Epigenetic marks, such as DNA methylation, are associated with gene silencing and repeat stability at the FMR1 locus. However, the mechanisms leading to gene silencing and repeat expansion are still ambiguous, because studies at the human genomic locus were limited until now. The FXS pluripotent stem cells, recently derived from FXS adult cells and FXS blastocysts, are new useful tools to examine these mechanisms at the human endogenous FMR1 locus. This review summarizes the epigenetic features and experimental studies of FXS human embryonic and FXS induced pluripotent stem cells, generated so far. This article is part of a Special Issue entitled SI: Exploiting human neurons.

Dynamic changes in replication timing and gene expression during lineage specification of human pluripotent stem cells

Duplication of the genome in mammalian cells occurs in a defined temporal order referred to as its replication-timing (RT) program. RT changes dynamically during development, regulated in units of 400-800 kb referred to as replication domains (RDs). Changes in RT are generally coordinated with transcriptional competence and changes in subnuclear position. We generated genome-wide RT profiles for 26 distinct human cell types, including embryonic stem cell (hESC)-derived, primary cells and established cell lines representing intermediate stages of endoderm, mesoderm, ectoderm, and neural crest (NC) development. We identified clusters of RDs that replicate at unique times in each stage (RT signatures) and confirmed global consolidation of the genome into larger synchronously replicating segments during differentiation. Surprisingly, transcriptome data revealed that the well-accepted correlation between early replication and transcriptional activity was restricted to RT-constitutive genes, whereas two-thirds of the genes that switched RT during differentiation were strongly expressed when late replicating in one or more cell types. Closer inspection revealed that transcription of this class of genes was frequently restricted to the lineage in which the RT switch occurred, but was induced prior to a late-to-early RT switch and/or down-regulated after an early-to-late RT switch. Analysis of transcriptional regulatory networks showed that this class of genes contains strong regulators of genes that were only expressed when early replicating. These results provide intriguing new insight into the complex relationship between transcription and RT regulation during human development.

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.

Epigenetic landscape for initiation of DNA replication

The key genetic process of DNA replication is initiated at specific sites referred to as replication origins. In eukaryotes, origins of DNA replication are not specified by a defined nucleotide sequence. Recent studies have shown that the structural context and topology of DNA sequence, chromatin features, and its transcriptional activity play an important role in origin choice. During differentiation and development, significant changes in chromatin organization and transcription occur, influencing origin activity and choice. In the last few years, a number of different genome-wide studies have broadened the understanding of replication origin regulation. In this review, we discuss the epigenetic factors and mechanisms that modulate origin choice and firing.

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.

Genome-wide copy number profiling of single cells in S-phase reveals DNA-replication domains

Single-cell genomics is revolutionizing basic genome research and clinical genetic diagnosis. However, none of the current research or clinical methods for single-cell analysis distinguishes between the analysis of a cell in G1-, S- or G2/M-phase of the cell cycle. Here, we demonstrate by means of array comparative genomic hybridization that charting the DNA copy number landscape of a cell in S-phase requires conceptually different approaches to that of a cell in G1- or G2/M-phase. Remarkably, despite single-cell whole-genome amplification artifacts, the log2 intensity ratios of single S-phase cells oscillate according to early and late replication domains, which in turn leads to the detection of significantly more DNA imbalances when compared with a cell in G1- or G2/M-phase. Although these DNA imbalances may, on the one hand, be falsely interpreted as genuine structural aberrations in the S-phase cell’s copy number profile and hence lead to misdiagnosis, on the other hand, the ability to detect replication domains genome wide in one cell has important applications in DNA-replication research. Genome-wide cell-type-specific early and late replicating domains have been identified by analyses of DNA from populations of cells, but cell-to-cell differences in DNA replication may be important in genome stability, disease aetiology and various other cellular processes.

Regulation of DNA replication timing on human chromosome by a cell-type specific DNA binding protein SATB1

BACKGROUND

Replication timing of metazoan DNA during S-phase may be determined by many factors including chromosome structures, nuclear positioning, patterns of histone modifications, and transcriptional activity. It may be determined by Mb-domain structures, termed as "replication domains", and recent findings indicate that replication timing is under developmental and cell type-specific regulation.

METHODOLOGY/PRINCIPAL FINDINGS

We examined replication timing on the human 5q23/31 3.5-Mb segment in T cells and non-T cells. We used two independent methods to determine replication timing. One is quantification of nascent replicating DNA in cell cycle-fractionated stage-specific S phase populations. The other is FISH analyses of replication foci. Although the locations of early- and late-replicating domains were common between the two cell lines, the timing transition region (TTR) between early and late domains were offset by 200-kb. We show that Special AT-rich sequence Binding protein 1 (SATB1), specifically expressed in T-cells, binds to the early domain immediately adjacent to TTR and delays the replication timing of the TTR. Measurement of the chromosome copy number along the TTR during synchronized S phase suggests that the fork movement may be slowed down by SATB1.

CONCLUSIONS

Our results reveal a novel role of SATB1 in cell type-specific regulation of replication timing along the chromosome.

Human telomeres replicate using chromosome-specific, rather than universal, replication programs

Human telomere replication initiates either from within telomere repeats or from within the subtelomere using a chromosome-specific replication program that appears conserved between different cell types.

Telomeric and adjacent subtelomeric heterochromatin pose significant challenges to the DNA replication machinery. Little is known about how replication progresses through these regions in human cells. Using single molecule analysis of replicated DNA (SMARD), we delineate the replication programs—i.e., origin distribution, termination site location, and fork rate and direction—of specific telomeres/subtelomeres of individual human chromosomes in two embryonic stem (ES) cell lines and two primary somatic cell types. We observe that replication can initiate within human telomere repeats but was most frequently accomplished by replisomes originating in the subtelomere. No major delay or pausing in fork progression was detected that might lead to telomere/subtelomere fragility. In addition, telomeres from different chromosomes from the same cell type displayed chromosome-specific replication programs rather than a universal program. Importantly, although there was some variation in the replication program of the same telomere in different cell types, the basic features of the program of a specific chromosome end appear to be conserved.

Global analysis reveals multiple pathways for unique regulation of mRNA decay in induced pluripotent stem cells

Pluripotency is a unique state in which cells can self-renew indefinitely but also retain the ability to differentiate into other cell types upon receipt of extracellular cues. Although it is clear that stem cells have a distinct transcriptional program, little is known about how alterations in post-transcriptional mechanisms, such as mRNA turnover, contribute to the achievement and maintenance of pluripotency. Here we have assessed the rates of decay for the majority of mRNAs expressed in induced pluripotent stem (iPS) cells and the fully differentiated human foreskin fibroblasts (HFFs) they were derived from. Comparison of decay rates in the two cell types led to the discovery of three independent regulatory mechanisms that allow coordinated turnover of specific groups of mRNAs. One mechanism results in increased stability of many histone mRNAs in iPS cells. A second pathway stabilizes a large set of zinc finger protein mRNAs, potentially through reduced levels of miRNAs that target them. Finally, a group of transcripts bearing 3' UTR C-rich sequence elements, many of which encode transcription factors, are significantly less stable in iPS cells. Intriguingly, two poly(C)-binding proteins that recognize this type of element are reciprocally expressed in iPS and HFF cells. Overall, our results highlight the importance of post-transcriptional control in pluripotent cells and identify miRNAs and RNA-binding proteins whose activity may coordinately control expression of a wide range of genes in iPS cells.

Regulation of DNA replication during development

As development unfolds, DNA replication is not only coordinated with cell proliferation, but is regulated uniquely in specific cell types and organs. This differential regulation of DNA synthesis requires crosstalk between DNA replication and differentiation. This dynamic aspect of DNA replication is highlighted by the finding that the distribution of replication origins varies between differentiated cell types and changes with differentiation. Moreover, differential DNA replication in some cell types can lead to increases or decreases in gene copy number along chromosomes. This review highlights the recent advances and technologies that have provided us with new insights into the developmental regulation of DNA replication.

DNA replication fading as proliferating cells advance in their commitment to terminal differentiation

Terminal differentiation is the process by which cycling cells stop proliferating to start new specific functions. It involves dramatic changes in chromatin organization as well as gene expression. In the present report we used cell flow cytometry and genome wide DNA combing to investigate DNA replication during murine erythroleukemia-induced terminal cell differentiation. The results obtained indicated that the rate of replication fork movement slows down and the inter-origin distance becomes shorter during the precommitment and commitment periods before cells stop proliferating and accumulate in G1. We propose this is a general feature caused by the progressive heterochromatinization that characterizes terminal cell differentiation.

Evidence for sequential and increasing activation of replication origins along replication timing gradients in the human genome

Genome-wide replication timing studies have suggested that mammalian chromosomes consist of megabase-scale domains of coordinated origin firing separated by large originless transition regions. Here, we report a quantitative genome-wide analysis of DNA replication kinetics in several human cell types that contradicts this view. DNA combing in HeLa cells sorted into four temporal compartments of S phase shows that replication origins are spaced at 40 kb intervals and fire as small clusters whose synchrony increases during S phase and that replication fork velocity (mean 0.7 kb/min, maximum 2.0 kb/min) remains constant and narrowly distributed through S phase. However, multi-scale analysis of a genome-wide replication timing profile shows a broad distribution of replication timing gradients with practically no regions larger than 100 kb replicating at less than 2 kb/min. Therefore, HeLa cells lack large regions of unidirectional fork progression. Temporal transition regions are replicated by sequential activation of origins at a rate that increases during S phase and replication timing gradients are set by the delay and the spacing between successive origin firings rather than by the velocity of single forks. Activation of internal origins in a specific temporal transition region is directly demonstrated by DNA combing of the IGH locus in HeLa cells. Analysis of published origin maps in HeLa cells and published replication timing and DNA combing data in several other cell types corroborate these findings, with the interesting exception of embryonic stem cells where regions of unidirectional fork progression seem more abundant. These results can be explained if origins fire independently of each other but under the control of long-range chromatin structure, or if replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA. These findings shed a new light on the replication timing program of mammalian genomes and provide a general model for their replication kinetics.

Eukaryotic chromosomes replicate from multiple replication origins that fire at different times in S phase. The mechanisms that specify origin position and firing time and coordinate origins to ensure complete genome duplication are unclear. Previous studies proposed either that origins are arranged in temporally coordinated groups or fire independently of each other in a stochastic manner. Here, we have performed a quantitative analysis of human genome replication kinetics using a combination of DNA combing, which reveals local patterns of origin firing and replication fork progression on single DNA molecules, and massive sequencing of newly replicated DNA, which reveals the population-averaged replication timing profile of the entire genome. We show that origins are activated synchronously in large regions of uniform replication timing but more gradually in temporal transition regions and that the rate of origin firing increases as replication progresses. Large regions of unidirectional fork progression are abundant in embryonic stem cells but rare in differentiated cells. We propose a model in which replication forks progressing from early origins stimulate initiation in nearby unreplicated DNA in a manner that explains the shape of the replication timing profile. These results provide a fundamental insight into the temporal regulation of mammalian genome replication.

Single molecule analysis of replicated DNA reveals the usage of multiple KSHV genome regions for latent replication

Kaposi's sarcoma associated herpesvirus (KSHV), an etiologic agent of Kaposi's sarcoma, Body Cavity Based Lymphoma and Multicentric Castleman's Disease, establishes lifelong latency in infected cells. The KSHV genome tethers to the host chromosome with the help of a latency associated nuclear antigen (LANA). Additionally, LANA supports replication of the latent origins within the terminal repeats by recruiting cellular factors. Our previous studies identified and characterized another latent origin, which supported the replication of plasmids ex-vivo without LANA expression in trans. Therefore identification of an additional origin site prompted us to analyze the entire KSHV genome for replication initiation sites using single molecule analysis of replicated DNA (SMARD). Our results showed that replication of DNA can initiate throughout the KSHV genome and the usage of these regions is not conserved in two different KSHV strains investigated. SMARD also showed that the utilization of multiple replication initiation sites occurs across large regions of the genome rather than a specified sequence. The replication origin of the terminal repeats showed only a slight preference for their usage indicating that LANA dependent origin at the terminal repeats (TR) plays only a limited role in genome duplication. Furthermore, we performed chromatin immunoprecipitation for ORC2 and MCM3, which are part of the pre-replication initiation complex to determine the genomic sites where these proteins accumulate, to provide further characterization of potential replication initiation sites on the KSHV genome. The ChIP data confirmed accumulation of these pre-RC proteins at multiple genomic sites in a cell cycle dependent manner. Our data also show that both the frequency and the sites of replication initiation vary within the two KSHV genomes studied here, suggesting that initiation of replication is likely to be affected by the genomic context rather than the DNA sequences.

Genomic approaches to the initiation of DNA replication and chromatin structure reveal a complex relationship

The mechanisms regulating the coordinate activation of tens of thousands of replication origins in multicellular organisms remain poorly explored. Recent advances in genomics have provided valuable information about the sites at which DNA replication is initiated and the selection mechanisms of specific sites in both yeast and vertebrates. Studies in yeast have advanced to the point that it is now possible to develop convincing models for origin selection. A general model has emerged, but yeast data have also revealed an unsuspected diversity of strategies for origin positioning. We focus here on the ways in which chromatin structure may affect the formation of pre-replication complexes, a prerequisite for origin activation. We also discuss the need to exercise caution when trying to extrapolate yeast models directly to more complex vertebrate genomes.