Evolutionarily conserved replication timing profiles predict long-range chromatin interactions and distinguish closely related cell types

DNA replication timing and long-range DNA interactions predict mutational landscapes of cancer genomes

Somatic copy-number alterations (SCNA) are a hallmark of many cancer types, but the mechanistic basis underlying their genome-wide patterns remains incompletely understood. Here we integrate data on DNA replication timing, long-range interactions between genomic material, and 331,724 SCNAs from 2,792 cancer samples classified into 26 cancer types. We report that genomic regions of similar replication timing are clustered spatially in the nucleus, that the two boundaries of SCNAs tend to be found in such regions, and that regions replicated early and late display distinct patterns of frequencies of SCNA boundaries, SCNA size and a preference for deletions over insertions. We show that long-range interaction and replication timing data alone can identify a significant proportion of SCNAs in an independent test data set. We propose a model for the generation of SCNAs in cancer, suggesting that data on spatial proximity of regions replicating at the same time can be used to predict the mutational landscapes of cancer genomes.

Unravelling global genome organization by 3C-seq

Eukaryotic genomes exist in the cell nucleus as an elaborate three-dimensional structure which reflects various nuclear processes such as transcription, DNA replication and repair. Next-generation sequencing (NGS) combined with chromosome conformation capture (3C), referred to as 3C-seq in this article, has recently been applied to the yeast and human genomes, revealing genome-wide views of functional associations among genes and their regulatory elements. Here, we compare the latest genomic approaches such as 3C-seq and ChIA-PET, and provide a condensed overview of how eukaryotic genomes are functionally organized in the nucleus.

On the opportunistic nature of transcription and replication initiation in the metazoan genome

Cellular identity and its response to external or internal signalling variations are encoded in a cell's genome as regulatory information. The genomic regions that specify this type of information are highly variable and degenerated in their sequence determinants, as it is becoming increasingly evident through the application of genome-scale methods to study gene expression. Here, we speculate that the same scenario applies to the regulatory regions controlling where DNA replication starts in the metazoan genome. We propose that replication origins cannot be defined as unique genomic features, but rather that DNA synthesis initiates opportunistically from accessible DNA sites, making cells highly robust and adaptable to environmental or developmental changes.

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.

DNA replication timing is maintained genome-wide in primary human myoblasts independent of D4Z4 contraction in FSH muscular dystrophy

Facioscapulohumeral muscular dystrophy (FSHD) is linked to contraction of an array of tandem 3.3-kb repeats (D4Z4) at 4q35.2 from 11-100 copies to 1-10 copies. The extent to which D4Z4 contraction at 4q35.2 affects overall 4q35.2 chromatin organization remains unclear. Because DNA replication timing is highly predictive of long-range chromatin interactions, we generated genome-wide replication-timing profiles for FSHD and control myogenic precursor cells. We compared non-immortalized myoblasts from four FSHD patients and three control individuals to each other and to a variety of other human cell types. This study also represents the first genome-wide comparison of replication timing profiles in non-immortalized human cell cultures. Myoblasts from both control and FSHD individuals all shared a myoblast-specific replication profile. In contrast, male and female individuals were readily distinguished by monoallelic differences in replication timing at DXZ4 and other regions across the X chromosome affected by X inactivation. We conclude that replication timing is a robust cell-type specific feature that is unaffected by FSHD-related D4Z4 contraction.

Replication timing: a fingerprint for cell identity and pluripotency

Many types of epigenetic profiling have been used to classify stem cells, stages of cellular differentiation, and cancer subtypes. Existing methods focus on local chromatin features such as DNA methylation and histone modifications that require extensive analysis for genome-wide coverage. Replication timing has emerged as a highly stable cell type-specific epigenetic feature that is regulated at the megabase-level and is easily and comprehensively analyzed genome-wide. Here, we describe a cell classification method using 67 individual replication profiles from 34 mouse and human cell lines and stem cell-derived tissues, including new data for mesendoderm, definitive endoderm, mesoderm and smooth muscle. Using a Monte-Carlo approach for selecting features of replication profiles conserved in each cell type, we identify “replication timing fingerprints” unique to each cell type and apply a k nearest neighbor approach to predict known and unknown cell types. Our method correctly classifies 67/67 independent replication-timing profiles, including those derived from closely related intermediate stages. We also apply this method to derive fingerprints for pluripotency in human and mouse cells. Interestingly, the mouse pluripotency fingerprint overlaps almost completely with previously identified genomic segments that switch from early to late replication as pluripotency is lost. Thereafter, replication timing and transcription within these regions become difficult to reprogram back to pluripotency, suggesting these regions highlight an epigenetic barrier to reprogramming. In addition, the major histone cluster Hist1 consistently becomes later replicating in committed cell types, and several histone H1 genes in this cluster are downregulated during differentiation, suggesting a possible instrument for the chromatin compaction observed during differentiation. Finally, we demonstrate that unknown samples can be classified independently using site-specific PCR against fingerprint regions. In sum, replication fingerprints provide a comprehensive means for cell characterization and are a promising tool for identifying regions with cell type-specific organization.

While continued advances in stem cell and cancer biology have uncovered a growing list of clinical applications for stem cell technology, errors in indentifying cell lines have undermined a number of recent studies, highlighting a growing need for improvements in cell typing methods for both basic biological and clinical applications of stem cells. Induced pluripotent stem cells (iPSCs)—adult cells reprogrammed to a pluripotent state—show great promise for patient-specific stem cell treatments, but more efficient derivation of iPSCs depends on a more comprehensive understanding of pluripotency. Here, we describe a method to identify sets of regions that replicate at unique times in any given cell type (replication timing fingerprints) using pluripotent stem cells as an example, and show that genes in the pluripotency fingerprint belong to a class previously shown to be resistant to reprogramming in iPSCs, identifying potential new target genes for more efficient iPSC production. We propose that the order in which DNA is replicated (replication timing) provides a novel means for classifying cell types, and can reveal cell type specific features of genome organization.

Gene expression: sizing it all up

Genomic architecture appears to be a largely unexplored component of gene expression. That architecture can be related to chromatin domains, transposable element neighborhoods, epigenetic modifications of the genome, and more. Although surely not the end of the story, we are learning that when it comes to gene expression, size is also important. We have been surprised to find that certain patterns of expression, tissue specific versus constitutive, or high expression versus low expression, are often associated with physical attributes of the gene and genome. Multiple studies have shown an inverse relationship between gene expression patterns and various physical parameters of the genome such as intron size, exon size, intron number, and size of intergenic regions. An increase in expression level and breadth often correlates with a decrease in the size of physical attributes of the gene. Three models have been proposed to explain these relationships. Contradictory results were found in several organisms when expression level and expression breadth were analyzed independently. However, when both factors were combined in a single study a novel relationship was revealed. At low levels of expression, an increase in expression breadth correlated with an increase in genic, intergenic, and intragenic sizes. Contrastingly, at high levels of expression, an increase in expression breadth inversely correlated with the size of the gene. In this article we explore the several hypotheses regarding genome physical parameters and gene expression.

Genetic variation and DNA replication timing, or why is there late replicating DNA?

Mutation rates vary significantly within the genome and across species. Recent studies revealed a long suspected replication-timing effect on mutation rate, but the mechanisms that regulate the increase in mutation rate as the genome is replicated remain unclear. Evidence is emerging, however, that DNA repair systems, in general, are less efficient in late replicating heterochromatic regions compared to early replicating euchromatic regions of the genome. At the same time, mutation rates in both vertebrates and invertebrates have been shown to vary with generation time (GT). GT is correlated with genome size, which suggests a possible nucleotypic effect on species-specific mutation rates. These and other observations all converge on a role for DNA replication checkpoints in modulating generation times and mutation rates during the DNA synthetic phase (S phase) of the cell cycle. The following will examine the potential role of the intra-S checkpoint in regulating cell cycle times (GT) and mutation rates in eukaryotes. This article was published online on August 5, 2011. An error was subsequently identified. This notice is included in the online and print versions to indicate that both have been corrected October 4, 2011.

Genome-wide depletion of replication initiation events in highly transcribed regions

This report investigates the mechanisms by which mammalian cells coordinate DNA replication with transcription and chromatin assembly. In yeast, DNA replication initiates within nucleosome-free regions, but studies in mammalian cells have not revealed a similar relationship. Here, we have used genome-wide massively parallel sequencing to map replication initiation events, thereby creating a database of all replication initiation sites within nonrepetitive DNA in two human cell lines. Mining this database revealed that genomic regions transcribed at moderate levels were generally associated with high replication initiation frequency. In genomic regions with high rates of transcription, very few replication initiation events were detected. High-resolution mapping of replication initiation sites showed that replication initiation events were absent from transcription start sites but were highly enriched in adjacent, downstream sequences. Methylation of CpG sequences strongly affected the location of replication initiation events, whereas histone modifications had minimal effects. These observations suggest that high levels of transcription interfere with formation of pre-replication protein complexes. Data presented here identify replication initiation sites throughout the genome, providing a foundation for further analyses of DNA-replication dynamics and cell-cycle progression.

Defining the replication program through the chromatin landscape

DNA replication is an essential cell cycle event required for the accurate and timely duplication of the chromosomes. It is essential that the genome is replicated accurately and completely within the confines of S-phase. Failure to completely copy the genome has the potential to result in catastrophic genomic instability. Replication initiates in a coordinated manner from multiple locations, termed origins of replication, distributed across each of the chromosomes. The selection of these origins of replication is a dynamic process responding to both developmental and tissue-specific signals. In this review, we explore the role of the local chromatin environment in regulating the DNA replication program at the level of origin selection and activation. Finally, there is increasing molecular evidence that the DNA replication program itself affects the chromatin landscape, suggesting that DNA replication is critical for both genetic and epigenetic inheritance.

Nuclear visions enhanced: chromatin structure, organization and dynamics

The EMBO Workshop on 'Chromatin Structure, Organization and Dynamics' took place in April 2011 in Prague, Czech Republic. Scientists from the life sciences, chemistry and biophysics presented their latest data on the generation of three-dimensional and, eventually, four-dimensional models of the genome, working to correlate changes in the organization of chromatin with the functional state of the genome.

On emerging nuclear order

Although the nonrandom nature of interphase chromosome arrangement is widely accepted, how nuclear organization relates to genomic function remains unclear. Nuclear subcompartments may play a role by offering rich microenvironments that regulate chromatin state and ensure optimal transcriptional efficiency. Technological advances now provide genome-wide and four-dimensional analyses, permitting global characterizations of nuclear order. These approaches will help uncover how seemingly separate nuclear processes may be coupled and aid in the effort to understand the role of nuclear organization in development and disease.

Genome-scale analysis of replication timing: from bench to bioinformatics

Replication timing profiles are cell type-specific and reflect genome organization changes during differentiation. In this protocol, we describe how to analyze genome-wide replication timing (RT) in mammalian cells. Asynchronously cycling cells are pulse labeled with the nucleotide analog 5-bromo-2-deoxyuridine (BrdU) and sorted into S-phase fractions on the basis of DNA content using flow cytometry. BrdU-labeled DNA from each fraction is immunoprecipitated, amplified, differentially labeled and co-hybridized to a whole-genome comparative genomic hybridization microarray, which is currently more cost effective than high-throughput sequencing and equally capable of resolving features at the biologically relevant level of tens to hundreds of kilobases. We also present a guide to analyzing the resulting data sets based on methods we use routinely. Subjects include normalization, scaling and data quality measures, LOESS (local polynomial) smoothing of RT values, segmentation of data into domains and assignment of timing values to gene promoters. Finally, we cover clustering methods and means to relate changes in the replication program to gene expression and other genetic and epigenetic data sets. Some experience with R or similar programming languages is assumed. All together, the protocol takes ∼3 weeks per batch of samples.

Emerging patterns of epigenomic variation

Fuelled by new sequencing technologies, epigenome mapping projects are revealing epigenomic variation at all levels of biological complexity, from species to cells. Comparisons of methylation profiles among species reveal evolutionary conservation of gene body methylation patterns, pointing to the fundamental role of epigenomes in gene regulation. At the human population level, epigenomic changes provide footprints of the effects of genomic variants within the vast nonprotein-coding fraction of the genome, and comparisons of the epigenomes of parents and their offspring point to quantitative epigenomic parent-of-origin effects confounding classical Mendelian genetics. At the organismal level, comparisons of epigenomes from diverse cell types provide insights into cellular differentiation. Finally, comparisons of epigenomes from monozygotic twins help dissect genetic and environmental influences on human phenotypes and longitudinal comparisons reveal aging-associated epigenomic drift. The development of new bioinformatic frameworks for comparative epigenome analysis is putting epigenome maps within the reach of researchers across a wide spectrum of biological disciplines.

Chromatin folding – from biology to polymer models and back

There is rapidly growing evidence that folding of the chromatin fibre inside the interphase nucleus has an important role in the regulation of gene expression. In particular, the formation of loops mediated by the interaction between specific regulatory elements, for instance enhancers and promoters, is crucial in gene control. Biochemical studies that were based on the chromosome conformation capture (3C) technology have confirmed that eukaryotic genomes are highly looped. Insight into the underlying principles comes from polymer models that explore the properties of the chromatin fibre inside the nucleus. Recent models indicate that chromatin looping can explain various properties of interphase chromatin, including chromatin compaction and compartmentalisation of chromosomes. Entropic effects have a key role in these models. In this Commentary, we give an overview of the recent conjunction of ideas regarding chromatin looping in the fields of biology and polymer physics. Starting from simple linear polymer models, we explain how specific folding properties emerge upon introducing loops and how this explains a variety of experimental observations. We also discuss different polymer models that describe chromatin folding and compare them to experimental data. Experimentally testing the predictions of such polymer models and their subsequent improvement on the basis of measurements provides a solid framework to begin to understand how our genome is folded and how folding relates to function.

Substitution patterns are under different influences in primates and rodents

There are large-scale variations of the GC-content along mammalian chromosomes that have been called isochore structures. Primates and rodents have different isochore structures, which suggests that these lineages exhibit different modes of GC-content evolution. It has been shown that, in the human lineage, GC-biased gene conversion (gBGC), a neutral process associated with meiotic recombination, acts on GC-content evolution by influencing A or T to G or C substitution rates. We computed genome-wide substitution patterns in the mouse lineage from multiple alignments and compared them with substitution patterns in the human lineage. We found that in the mouse lineage, gBGC is active but weaker than in the human lineage and that male-specific recombination better predicts GC-content evolution than female-specific recombination. Furthermore, we were able to show that G or C to A or T substitution rates are predicted by a combination of different factors in both lineages. A or T to G or C substitution rates are most strongly predicted by meiotic recombination in the human lineage but by CpG odds ratio (the observed CpG frequency normalized by the expected CpG frequency) in the mouse lineage, suggesting that substitution patterns are under different influences in primates and rodents.

Reprogramming to pluripotency: stepwise resetting of the epigenetic landscape

In 2006, the "wall came down" that limited the experimental conversion of differentiated cells into the pluripotent state. In a landmark report, Shinya Yamanaka's group described that a handful of transcription factors (Oct4, Sox2, Klf4 and c-Myc) can convert a differentiated cell back to pluripotency over the course of a few weeks, thus reprograming them into induced pluripotent stem (iPS) cells. The birth of iPS cells started off a rush among researchers to increase the efficiency of the reprogramming process, to reveal the underlying mechanistic events, and allowed the generation of patient- and disease-specific human iPS cells, which have the potential to be converted into relevant specialized cell types for replacement therapies and disease modeling. This review addresses the steps involved in resetting the epigenetic landscape during reprogramming. Apparently, defined events occur during the course of the reprogramming process. Immediately, upon expression of the reprogramming factors, some cells start to divide faster and quickly begin to lose their differentiated cell characteristics with robust downregulation of somatic genes. Only a subset of cells continue to upregulate the embryonic expression program, and finally, pluripotency genes are upregulated establishing an embryonic stem cell-like transcriptome and epigenome with pluripotent capabilities. Understanding reprogramming to pluripotency will inform mechanistic studies of lineage switching, in which differentiated cells from one lineage can be directly reprogrammed into another without going through a pluripotent intermediate.

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.

Charting histone modifications and the functional organization of mammalian genomes

A succession of technological advances over the past decade have enabled researchers to chart maps of histone modifications and related chromatin structures with increasing accuracy, comprehensiveness and throughput. The resulting data sets highlight the interplay between chromatin and genome function, dynamic variations in chromatin structure across cellular conditions, and emerging roles for large-scale domains and higher-ordered chromatin organization. Here we review a selection of recent studies that have probed histone modifications and successive layers of chromatin structure in mammalian genomes, the patterns that have been identified and future directions for research.

Cell fate transitions and the replication timing decision point

Recent findings suggest that large-scale remodeling of three dimensional (3D) chromatin architecture occurs during a brief period in early G1 phase termed the replication timing decision point (TDP). In this speculative article, I suggest that the TDP may represent an as yet unappreciated window of opportunity for extracellular cues to influence 3D architecture during stem cell fate decisions. I also describe several testable predictions of this hypothesis.

Genomics tools for unraveling chromosome architecture

The spatial organization of chromosomes inside the cell nucleus is still poorly understood. This organization is guided by intra- and interchromosomal contacts and by interactions of specific chromosomal loci with relatively fixed nuclear 'landmarks' such as the nuclear envelope and the nucleolus. Researchers have begun to use new molecular genome-wide mapping techniques to uncover both types of molecular interactions, providing insights into the fundamental principles of interphase chromosome folding.

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.

Clusters, factories and domains: The complex structure of S-phase comes into focus

During S-phase of the cell cycle, chromosomal DNA is replicated according to a complex replication timing program, with megabase-sized domains replicating at different times. DNA fibre analysis reveals that clusters of adjacent replication origins fire near-synchronously. Analysis of replicating cells by light microscopy shows that DNA synthesis occurs in discrete foci or factories. The relationship between timing domains, origin clusters and replication foci is currently unclear. Recent work, using a hybrid Xenopus/hamster replication system, has shown that when CDK levels are manipulated during S-phase the activation of replication factories can be uncoupled from progression through the replication timing program. Here, we use data from this hybrid system to investigate potential relationships between timing domains, origin clusters and replication foci. We suggest that each timing domain typically comprises several replicon clusters, which are usually processed sequentially by replication factories. We discuss how replication might be regulated at different levels to create this complex organisation and the potential involvement of CDKs in this process.

Comparative analysis of DNA replication timing reveals conserved large-scale chromosomal architecture

Recent evidence suggests that the timing of DNA replication is coordinated across megabase-scale domains in metazoan genomes, yet the importance of this aspect of genome organization is unclear. Here we show that replication timing is remarkably conserved between human and mouse, uncovering large regions that may have been governed by similar replication dynamics since these species have diverged. This conservation is both tissue-specific and independent of the genomic G+C content conservation. Moreover, we show that time of replication is globally conserved despite numerous large-scale genome rearrangements. We systematically identify rearrangement fusion points and demonstrate that replication time can be locally diverged at these loci. Conversely, rearrangements are shown to be correlated with early replication and physical chromosomal proximity. These results suggest that large chromosomal domains of coordinated replication are shuffled by evolution while conserving the large-scale nuclear architecture of the genome.

During S-phase of the cell cycle, chromosomal DNA is replicated in a complex process involving the coordinated activity of thousands of replication forks, each of which duplicates a long stretch of DNA. Recent experiments revealed that the genome is replicating as a mosaic of large-scale early and late chromosomal domains and that this high-level domain organization is correlated with genomic properties like gene density and nucleotide composition. We compared genome-wide replication time maps of compatible human and mouse cells and revealed that their organization into replication domains is highly conserved despite the numerous large-scale genome rearrangements separating the two species. Analysis of recent chromosomal interaction data shows that regions with similar time of replication are more frequently interacting with each other than expected. The data also show that evolutionary rearrangements have predominantly occurred between regions that have similar time of replication and higher-than-expected chromosomal proximity. Our data suggests that the genome, while being continuously rearranged by evolution, maintains a conserved domain organization. Whether this conservation is driven by selection, or is a consequence of the rearrangement process itself, can be resolved by enhancing the comparative approach proposed here.

Differentiation and large scale spatial organization of the genome

The spatial organization of the genome plays an important role in the regulation of nuclear functions and undergoes large scale changes during differentiation. These changes in the nuclear distribution of chromatin are, in a complex way, related to transcriptional status and epigenetic modifications. Recent studies emphasize the roles that gene promoters and alterations in replication timing play in the spatial reorganization of chromatin during cell differentiation. Changes in the association of chromatin regions with the nuclear lamina also emerge as a significant factor of transcriptional regulation. New results suggest that the spatial organization of chromatin in embryonic stem cells may be important for maintenance of the pluripotent state, whereas the nuclear architecture of differentiated cells facilitates formation of transcriptionally active zones with shared transcription and splicing machinery.

In search of lost timing

Researchers learn more about how cells schedule DNA replication by determining when they tear up the timetable.

G2 phase chromatin lacks determinants of replication timing

Chromatin spatial organization helps establish the replication timing decision point at early G1. However, at G2, although retained, chromatin organization is no longer necessary or sufficient to maintain the replication timing program.

DNA replication in all eukaryotes follows a defined replication timing program, the molecular mechanism of which remains elusive. Using a Xenopus laevis egg extract replication system, we previously demonstrated that replication timing is established during early G1 phase of the cell cycle (timing decision point [TDP]), which is coincident with the repositioning and anchorage of chromatin in the newly formed nucleus. In this study, we use this same system to show that G2 phase chromatin lacks determinants of replication timing but maintains the overall spatial organization of chromatin domains, and we confirm this finding by genome-wide analysis of rereplication in vivo. In contrast, chromatin from quiescent cells retains replication timing but exhibits disrupted spatial organization. These data support a model in which events at the TDP, facilitated by chromatin spatial organization, establish determinants of replication timing that persist independent of spatial organization until the process of chromatin replication during S phase erases those determinants.