Search

This protocol is an extension to: Nat. Protoc. 6, 870-895 (2014); doi:10.1038/nprot.2011.328; published online 02 June 2011Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early- and late-replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and subnuclear position. Moreover, RT is regulated during development and is altered in diseases. Here, we describe E/L Repli-seq, an extension of our Repli-chip protocol. E/L Repli-seq is a rapid, robust and relatively inexpensive protocol for analyzing RT by next-generation sequencing (NGS), allowing genome-wide assessment of how cellular processes are linked to RT. Briefly, cells are pulse-labeled with BrdU, and early and late S-phase fractions are sorted by flow cytometry. Labeled nascent DNA is immunoprecipitated from both fractions and sequenced. Data processing leads to a single bedGraph file containing the ratio of nascent DNA from early versus late S-phase fractions. The results are comparable to those of Repli-chip, with the additional benefits of genome-wide sequence information and an increased dynamic range. We also provide computational pipelines for downstream analyses, for parsing phased genomes using single-nucleotide polymorphisms (SNPs) to analyze RT allelic asynchrony, and for direct comparison to Repli-chip data. This protocol can be performed in up to 3 d before sequencing, and requires basic cellular and molecular biology skills, as well as a basic understanding of Unix and R.

The role of DNA sequence in determining replication timing (RT) and chromatin higher order organization remains elusive. To address this question, we have developed an extra-chromosomal replication system (E-BACs) consisting of ∼200 kb human bacterial artificial chromosomes (BACs) modified with Epstein-Barr virus (EBV) stable segregation elements. E-BACs were stably maintained as autonomous mini-chromosomes in EBNA1-expressing HeLa or human induced pluripotent stem cells (hiPSCs) and established distinct RT patterns. An E-BAC harboring an early replicating chromosomal region replicated early during S phase, while E-BACs derived from RT transition regions (TTRs) and late replicating regions replicated in mid to late S phase. Analysis of E-BAC interactions with cellular chromatin (4C-seq) revealed that the early replicating E-BAC interacted broadly throughout the genome and preferentially with the early replicating compartment of the nucleus. In contrast, mid- to late-replicating E-BACs interacted with more specific late replicating chromosomal segments, some of which were shared between different E-BACs. Together, we describe a versatile system in which to study the structure and function of chromosomal segments that are stably maintained separately from the influence of cellular chromosome context.

Transcriptionally inactive genes are often positioned at the nuclear lamina (NL), as part of large lamina‐associated domains (LADs). Activation of such genes is often accompanied by repositioning toward the nuclear interior. How this process works and how it impacts flanking chromosomal regions are poorly understood. We addressed these questions by systematic activation or inactivation of individual genes, followed by detailed genome‐wide analysis of NL interactions, replication timing, and transcription patterns. Gene activation inside LADs typically causes NL detachment of the entire transcription unit, but rarely more than 50–100 kb of flanking DNA, even when multiple neighboring genes are activated. The degree of detachment depends on the expression level and the length of the activated gene. Loss of NL interactions coincides with a switch from late to early replication timing, but the latter can involve longer stretches of DNA. Inactivation of active genes can lead to increased NL contacts. These extensive datasets are a resource for the analysis of LAD rewiring by transcription and reveal a remarkable flexibility of interphase chromosomes., Genomic maps acquired after transcriptional activation or inactivation of genes reveals how transcription affects their interaction with the nuclear lamina.

The 3D organization of mammalian chromatin was described more than 30 years ago by visualizing sites of DNA synthesis at different times during the S phase of the cell cycle. These early cytogenetic studies revealed structurally stable chromosome domains organized into subnuclear compartments. Active-gene-rich domains in the nuclear interior replicate early, whereas more condensed chromatin domains that are largely at the nuclear and nucleolar periphery replicate later. During the past decade, this spatiotemporal DNA replication programme has been mapped along the genome and found to correlate with epigenetic marks, transcriptional activity and features of 3D genome architecture such as chromosome compartments and topologically associated domains. But the causal relationship between these features and DNA replication timing and the regulatory mechanisms involved have remained an enigma. The recent identification of cis-acting elements regulating the replication time and 3D architecture of individual replication domains and of long non-coding RNAs that coordinate whole chromosome replication provide insights into such mechanisms. Different genomic regions are replicated at different times during the S phase of the cell cycle, forming early- and late-replicating domains that occupy different locations in the nucleus. The recent identification of specific DNA sequences and long non-coding RNAs that regulate DNA replication timing is providing key insights into the roles of replication timing and into timing and 3D organization.

Transcriptionally inactive genes are often positioned at the nuclear lamina (NL), as part of large lamina-associated domains (LADs). Activation of such genes is often accompanied by repositioning towards the nuclear interior. How this process works and how it impacts flanking chromosomal regions is poorly understood. We addressed these questions by systematic manipulation of gene activity and detailed analysis of NL interactions. Activation of genes inside LADs typically causes detachment of the entire transcription unit but rarely more than 50-100 kb of flanking DNA, even when multiple neighboring genes are activated. The degree of detachment depends on the expression level and the length of the activated gene. Loss of NL interactions coincides with a switch from late to early replication timing, but the latter can involve longer stretches of DNA. These findings show how NL interactions can be shaped locally by transcription and point to a remarkable flexibility of interphase chromosomes.

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.

SUMMARYThe temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, butcis-elements regulating spatio-temporal control of replication have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discretecis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.Highlightscis-elements (ERCEs) regulate large scale chromosome structure and functionMultiple ERCEs cooperatively control domain-wide replicationERCEs harbor prominent active chromatin features and form CTCF-independent loopsERCEs enable genetic dissection of large-scale chromosome structure-function.

The temporal order of DNA replication (replication timing [RT]) is highly coupled with genome architecture, but cis-elements regulating either remain elusive. We created a series of CRISPR-mediated deletions and inversions of a pluripotency-associated topologically associating domain (TAD) in mouse ESCs. CTCF-associated domain boundaries were dispensable for RT. CTCF protein depletion weakened most TAD boundaries but had no effect on RT or A/B compartmentalization genome-wide. By contrast, deletion of three intra-TAD CTCF-independent 3D contact sites caused a domain-wide early-to-late RT shift, an A-to-B compartment switch, weakening of TAD architecture, and loss of transcription. The dispensability of TAD boundaries and thenecessity of these "early replication control elements" (ERCEs) was validated by deletions and inversions at additional domains. Our results demonstrate that discrete cis-regulatory elements orchestrate domain-wide RT, A/B compartmentalization, TAD architecture, and transcription, revealing fundamental principles linking genome structure and function.

It is now well accepted that defined architectural compartments within the cell nucleus can regulate the transcriptional activity of chromosomal domains within their vicinity. However, it is generally unclear how these compartments are formed. The nuclear periphery has received a great deal of attention as a repressive compartment that is implicated in many cellular functions during development and disease. The inner nuclear membrane, the nuclear lamina, and associated proteins compose the nuclear periphery and together they interact with proximal chromatin creating a repressive environment. A new study by Harr et al. identifies specific protein:DNA interactions and epigenetic states necessary to re-position chromatin to the nuclear periphery in a cell-type specific manner. Here, we review concepts in gene positioning within the nucleus and current accepted models of dynamic gene repositioning within the nucleus during differentiation. This study highlights that myriad pathways lead to nuclear organization.

Mammalian genomes are partitioned into domains that replicate in a defined temporal order. These domains can replicate at similar times in all cell types (constitutive) or at cell type-specific times (developmental). Genome-wide chromatin conformation capture (Hi-C) has revealed sub-megabase topologically associating domains (TADs), which are the structural counterparts of replication domains. Hi-C also segregates inter-TAD contacts into defined 3D spatial compartments that align precisely to genome-wide replication timing profiles. Determinants of the replication-timing program are re-established during early G1 phase of each cell cycle and lost in G2 phase, but it is not known when TAD structure and inter-TAD contacts are re-established after their elimination during mitosis. Here, we use multiplexed 4C-seq to study dynamic changes in chromatin organization during early G1. We find that both establishment of TADs and their compartmentalization occur during early G1, within the same time frame as establishment of the replication-timing program. Once established, this 3D organization is preserved either after withdrawal into quiescence or for the remainder of interphase including G2 phase, implying 3D structure is not sufficient to maintain replication timing. Finally, we find that developmental domains are less well compartmentalized than constitutive domains and display chromatin properties that distinguish them from early and late constitutive domains. Overall, this study uncovers a strong connection between chromatin re-organization during G1, establishment of replication timing, and its developmental control.

The temporal order of DNA replication (replication timing, RT) is highly coupled with genome architecture, but cis-elements regulating spatio-temporal control of replicatio have remained elusive. We performed an extensive series of CRISPR mediated deletions and inversions and high-resolution capture Hi-C of a pluripotency associated domain (DppA2/4) in mouse embryonic stem cells. Whereas CTCF mediated loops and chromatin domain boundaries were dispensable, deletion of three intra-domain prominent CTCF-independent 3D contact sites caused a domain-wide delay in RT, shift in sub-nuclear chromatin compartment and loss of transcriptional activity, These “early replication control elements” (ERCEs) display prominent chromatin features resembling enhancers/promoters and individual and pair-wise deletions of the ERCEs confirmed their partial redundancy and interdependency in controlling domain-wide RT and transcription. Our results demonstrate that discrete cis-regulatory elements mediate domain-wide RT, chromatin compartmentalization, and transcription, representing a major advance in dissecting the relationship between genome structure and function.

International audience; Progeroid syndromes are rare genetic disorders that phenotypically resemble natural aging. Different causal mutations have been identified, but no molecular alterations have been identified that are in common to these diseases. DNA replication timing (RT) is a robust cell type-specific epigenetic feature highly conserved in the same cell types from different individuals but altered in disease. Here, we characterized DNA RT program alterations in Hutchinson–Gilford progeria syndrome (HGPS) and Rothmund–Thomson syndrome (RTS) patients compared with natural aging and cellular senescence. Our results identified a progeroid-specific RT signature that is common to cells from three HGPS and three RTS patients and distinguishes them from healthy individuals across a wide range of ages. Among the RT abnormalities, we identified the tumor protein p63 gene (TP63) as a gene marker for progeroid syndromes. By using the redifferentiation of four patient-derived induced pluripotent stem cells as a model for the onset of progeroid syndromes, we tracked the progression of RT abnormalities during development, revealing altered RT of the TP63 gene as an early event in disease progression of both HGPS and RTS. Moreover, the RT abnormalities in progeroid patients were associated with altered isoform expression of TP63. Our findings demonstrate the value of RT studies to identify biomarkers not detected by other methods, reveal abnormal TP63 RT as an early event in progeroid disease progression, and suggest TP63 gene regulation as a potential therapeutic target.

Cycling cells duplicate their DNA content during S phase, following a defined program called replication timing (RT). Early and late replicating regions differ in terms of mutation rates, transcriptional activity, chromatin marks and sub-nuclear position. Moreover, RT is regulated during development and is altered in disease. Exploring mechanisms linking RT to other cellular processes in normal and diseased cells will be facilitated by rapid and robust methods with which to measure RT genome wide. Here, we describe a rapid, robust and relatively inexpensive protocol to analyze genome-wide RT by next-generation sequencing (NGS). This protocol yields highly reproducible results across laboratories and platforms. We also provide computational pipelines for analysis, parsing phased genomes using single nucleotide polymorphisms (SNP) for analyzing RT allelic asynchrony, and for direct comparison to Repli-chip data obtained by analyzing nascent DNA by microarrays.

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.

A recent flurry of reports correlates replication timing (RT) with mutation rates during both evolution and cancer. Specifically, point mutations and copy number losses correlate with late replication, while copy number gains and other rearrangements correlate with early replication. In some cases, plausible mechanisms have been proposed. Point mutation rates may reflect temporal variation in repair mechanisms. Transcription-induced double-strand breaks are expected to occur in transcriptionally active early replicating chromatin. Fusion partners are generally in close proximity, and chromatin in close proximity replicates at similar times. However, temporal enrichment of copy number gains and losses remains an enigma. Moreover, many conclusions are compromised by a lack of matched RT and sequence datasets, the filtering out of developmental variation in RT, and the use of somatic cell lines to make inferences about germline evolution.