Structural organization of human replication timing domains

Genome architecture plasticity underlies DNA replication timing dynamics in cell differentiation

During the S-phase of eukaryotic cell cycle, DNA is replicated in a dedicatedly regulated temporal order, with regions containing active and inactive genes replicated early and late, respectively. Recent advances in sequencing technology allow us to explore the connection between replication timing (RT), histone modifications, and three-dimensional (3D) chromatin structure in diverse cell types. To characterize the dynamics during cell differentiation, corresponding sequencing data for human embryonic stem cells and four differentiated cell types were collected. By comparing RT and its extent of conservation before and after germ layer specification, the human genome was partitioned into distinct categories. Each category is then subject to comparisons on genomic, epigenetic, and chromatin 3D structural features. As expected, while constitutive early and late replication regions showed active and inactive features, respectively, dynamic regions with switched RT showed intermediate features. Surprisingly, although early-to-late replication and late-to-early replication regions showed similar histone modification patterns in hESCs, their structural preferences were opposite. Specifically, in hESCs, early-to-late replication regions tended to appear in the B compartment and large topologically associated domains, while late-to-early replication regions showed the opposite. Our results uncover the coordinated regulation of RT and 3D genome structure that underlies the loss of pluripotency and lineage commitment and indicate the importance and potential roles of genome architecture in biological processes.

High-throughput analysis of single human cells reveals the complex nature of DNA replication timing control

DNA replication initiates from replication origins firing throughout S phase. Debate remains about whether origins are a fixed set of loci, or a loose agglomeration of potential sites used stochastically in individual cells, and about how consistent their firing time is. We develop an approach to profile DNA replication from whole-genome sequencing of thousands of single cells, which includes in silico flow cytometry, a method for discriminating replicating and non-replicating cells. Using two microfluidic platforms, we analyze up to 2437 replicating cells from a single sample. The resolution and scale of the data allow focused analysis of replication initiation sites, demonstrating that most occur in confined genomic regions. While initiation order is remarkably similar across cells, we unexpectedly identify several subtypes of initiation regions in late-replicating regions. Taken together, high throughput, high resolution sequencing of individual cells reveals previously underappreciated variability in replication initiation and progression.

Chromatin Network Analyses: Towards Structure-Function Relationships in Epigenomics

Recent technological advances have allowed us to map chromatin conformation and uncover the genome’s spatial organization of the genome inside the nucleus. These experiments have revealed the complexities of genome folding, characterized by the presence of loops and domains at different scales, which can change across development and in different cell types. There is strong evidence for a relationship between the topological properties of chromatin contacts and cellular phenotype. Chromatin can be represented as a network, in which genomic fragments are the nodes and connections represent experimentally observed spatial proximity of two genomically distant regions in a specific cell type or biological condition. With this approach we can consider a variety of chromatin features in association with the 3D structure, investigating how nuclear chromatin organization can be related to gene regulation, replication, malignancy, phenotypic variability and plasticity. We briefly review the results obtained on genome architecture through network theoretic approaches. As previously observed in protein-protein interaction networks and many types of non-biological networks, external conditions could shape network topology through a yet unidentified structure-function relationship. Similar to scientists studying the brain, we are confronted with a duality between a spatially embedded network of physical contacts, a related network of correlation in the dynamics of network nodes and, finally, an abstract definition of function of this network, related to phenotype. We summarise major developments in the study of networks in other fields, which we think can suggest a path towards better understanding how 3D genome configuration can impact biological function and adaptation to the environment.

Organization of DNA Replication Origin Firing in Xenopus Egg Extracts: The Role of Intra-S Checkpoint

During cell division, the duplication of the genome starts at multiple positions called replication origins. Origin firing requires the interaction of rate-limiting factors with potential origins during the S(ynthesis)-phase of the cell cycle. Origins fire as synchronous clusters which is proposed to be regulated by the intra-S checkpoint. By modelling the unchallenged, the checkpoint-inhibited and the checkpoint protein Chk1 over-expressed replication pattern of single DNA molecules from Xenopus sperm chromatin replicated in egg extracts, we demonstrate that the quantitative modelling of data requires: (1) a segmentation of the genome into regions of low and high probability of origin firing; (2) that regions with high probability of origin firing escape intra-S checkpoint regulation and (3) the variability of the rate of DNA synthesis close to replication forks is a necessary ingredient that should be taken in to account in order to describe the dynamic of replication origin firing. This model implies that the observed origin clustering emerges from the apparent synchrony of origin firing in regions with high probability of origin firing and challenge the assumption that the intra-S checkpoint is the main regulator of origin clustering.

Human ORC/MCM density is low in active genes and correlates with replication time but does not delimit initiation zones

Eukaryotic DNA replication initiates during S phase from origins that have been licensed in the preceding G1 phase. Here, we compare ChIP-seq profiles of the licensing factors Orc2, Orc3, Mcm3, and Mcm7 with gene expression, replication timing, and fork directionality profiles obtained by RNA-seq, Repli-seq, and OK-seq. Both, the origin recognition complex (ORC) and the minichromosome maintenance complex (MCM) are significantly and homogeneously depleted from transcribed genes, enriched at gene promoters, and more abundant in early- than in late-replicating domains. Surprisingly, after controlling these variables, no difference in ORC/MCM density is detected between initiation zones, termination zones, unidirectionally replicating regions, and randomly replicating regions. Therefore, ORC/MCM density correlates with replication timing but does not solely regulate the probability of replication initiation. Interestingly, H4K20me3, a histone modification proposed to facilitate late origin licensing, was enriched in late-replicating initiation zones and gene deserts of stochastic replication fork direction. We discuss potential mechanisms specifying when and where replication initiates in human cells.

Resolution of Complex Issues in Genome Regulation and Cancer Requires Non-Linear and Network-Based Thermodynamics

The apparent lack of success in curing cancer that was evidenced in the last four decades of molecular medicine indicates the need for a global re-thinking both its nature and the biological approaches that we are taking in its solution. The reductionist, one gene/one protein method that has served us well until now, and that still dominates in biomedicine, requires complementation with a more systemic/holistic approach, to address the huge problem of cross-talk between more than 20,000 protein-coding genes, about 100,000 protein types, and the multiple layers of biological organization. In this perspective, the relationship between the chromatin network organization and gene expression regulation plays a fundamental role. The elucidation of such a relationship requires a non-linear thermodynamics approach to these biological systems. This change of perspective is a necessary step for developing successful 'tumour-reversion' therapeutic strategies.

Replication timing in Drosophila and its peculiarities in polytene chromosomes

Drosophila melanogaster is one of the popular model organisms in DNA replication studies. Since the 1960s, DNA replication of polytene chromosomes has been extensively studied by cytological methods. In the recent two decades, the progress in our understanding of DNA replication was associated with new techniques. Use of fluorescent dyes increased the resolution of cytological methods significantly. High-throughput methods allowed analysis of DNA replication on a genome scale, as well as its correlation with chromatin structure and gene activi ty. Precise mapping of the cytological structures of polytene chromosomes to the genome assembly allowed comparison of replication between polytene chromosomes and chromosomes of diploid cells. New features of replication characteristic for D. melanogaster were described for both diploid and polytene chromosomes. Comparison of genomic replication profiles revealed a significant similarity between Drosophila and other well-studi ed eukaryotic species, such as human. Early replication is often confined to intensely transcribed gene-dense regions characterized by multiple replication initiation sites. Features of DNA replication in Drosophila might be explained by a compact genome. The organization of replication in polytene chromosomes has much in common with the organization of replication in chromosomes in diploid cells. The most important feature of replication in polytene chromosomes is its low rate and the dependence of S-phase duration on many factors: external and internal, local and global. The speed of replication forks in D. melanogaster polytene chromosomes is affected by SUUR and Rif1 proteins. It is not known yet how universal the mechanisms associated with these factors are, but their study is very promising.

The Protective Role of Dormant Origins in Response to Replicative Stress

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

The eukaryotic bell-shaped temporal rate of DNA replication origin firing emanates from a balance between origin activation and passivation

The time-dependent rate I(t) of origin firing per length of unreplicated DNA presents a universal bell shape in eukaryotes that has been interpreted as the result of a complex time-evolving interaction between origins and limiting firing factors. Here, we show that a normal diffusion of replication fork components towards localized potential replication origins (p-oris) can more simply account for the I(t) universal bell shape, as a consequence of a competition between the origin firing time and the time needed to replicate DNA separating two neighboring p-oris. We predict the I(t) maximal value to be the product of the replication fork speed with the squared p-ori density. We show that this relation is robustly observed in simulations and in experimental data for several eukaryotes. Our work underlines that fork-component recycling and potential origins localization are sufficient spatial ingredients to explain the universality of DNA replication kinetics.

Before a cell can divide, it must duplicate its DNA. In eukaryotes – organisms such as animals and fungi, which store their DNA in the cell’s nucleus – DNA replication starts at specific sites in the genome called replication origins. At each origin sits a protein complex that will activate when it randomly captures an activating protein that diffuses within the nucleus. Once a replication origin activates or “fires”, the complex then splits into two new complexes that move away from each other as they duplicate the DNA. If an active complex collides with an inactive one at another origin, the latter is inactivated – a phenomenon known as origin passivation. When two active complexes meet, they release the activating proteins, which diffuse away and eventually activate other origins in unreplicated DNA.

The number of origins that activate each minute divided by the length of unreplicated DNA is referred to as the “rate of origin firing”. In all eukaryotes, this rate – also known as I(t) – follows the same pattern. First, it increases until more than half of the DNA is duplicated. Then it decreases until everything is duplicated. This means that, if plotted out, the graph of origin firing rate would always be a bell-shaped curve, even for organisms with genomes of different sizes that have different numbers of origins. The reason for this universal shape remained unclear.

Scientists had tried to create numerical simulations that model the rate of origin firing. However, for these simulations to reproduce the bell-shape curve, a number of untested assumptions had to be made about how DNA replication takes place. In addition, these models ignored the fact that it takes time to replicate the DNA between origins.

To take this time into account, Arbona et al. instead decided to model the replication origins as discrete and distinct entities. This way of building the mathematical model succeeded in reproducing the universal bell curve shape without additional assumptions. With this simulation, the balance between origin activation and passivation is enough to achieve the observed pattern.

The new model also predicts that the maximum rate of origin firing is determined by the speed of DNA replication and the density of origins in the genome. Arbona et al. verified this prediction in yeast, fly, frog and human cells – organisms with different sized genomes that take between 20 minutes and 8 hours to replicate their DNA. Lastly, the prediction also held true in yeast treated with hydroxyurea, an anticancer drug that slows DNA replication.

A better understanding of DNA replication can help scientists to understand how this process is perturbed in cancers and how drugs that target DNA replication can treat these diseases. Future work will explore how the 3D organization of the genome affects the diffusion of activating proteins within the cell nucleus.

Similarity in replication timing between polytene and diploid cells is associated with the organization of the Drosophila genome

Morphologically, polytene chromosomes of Drosophila melanogaster consist of compact “black” bands alternating with less compact “grey” bands and interbands. We developed a comprehensive approach that combines cytological mapping data of FlyBase-annotated genes and novel tools for predicting cytogenetic features of chromosomes on the basis of their protein composition and determined the genomic coordinates for all black bands of polytene chromosome 2R. By a PCNA immunostaining assay, we obtained the replication timetable for all the bands mapped. The results allowed us to compare replication timing between polytene chromosomes in salivary glands and chromosomes from cultured diploid cell lines and to observe a substantial similarity in the global replication patterns at the band resolution level. In both kinds of chromosomes, the intervals between black bands correspond to early replication initiation zones. Black bands are depleted of replication initiation events and are characterized by a gradient of replication timing; therefore, the time of replication completion correlates with the band length. The bands are characterized by low gene density, contain predominantly tissue-specific genes, and are represented by silent chromatin types in various tissues. The borders of black bands correspond well to the borders of topological domains as well as to the borders of the zones showing H3K27me3, SUUR, and LAMIN enrichment. In conclusion, the characteristic pattern of polytene chromosomes reflects partitioning of the Drosophila genome into two global types of domains with contrasting properties. This partitioning is conserved in different tissues and determines replication timing in Drosophila.

Genome-wide analysis of replication timing by next-generation sequencing with E/L Repli-seq

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.

Optical computed tomography for spatially isotropic four-dimensional imaging of live single cells

Quantitative three-dimensional (3D) computed tomography (CT) imaging of living single cells enables orientation-independent morphometric analysis of the intricacies of cellular physiology. Since its invention, x-ray CT has become indispensable in the clinic for diagnostic and prognostic purposes due to its quantitative absorption-based imaging in true 3D that allows objects of interest to be viewed and measured from any orientation. However, x-ray CT has not been useful at the level of single cells because there is insufficient contrast to form an image. Recently, optical CT has been developed successfully for fixed cells, but this technology called Cell-CT is incompatible with live-cell imaging due to the use of stains, such as hematoxylin, that are not compatible with cell viability. We present a novel development of optical CT for quantitative, multispectral functional 4D (three spatial + one spectral dimension) imaging of living single cells. The method applied to immune system cells offers truly isotropic 3D spatial resolution and enables time-resolved imaging studies of cells suspended in aqueous medium. Using live-cell optical CT, we found a heterogeneous response to mitochondrial fission inhibition in mouse macrophages and differential basal remodeling of small (0.1 to 1 fl) and large (1 to 20 fl) nuclear and mitochondrial structures on a 20- to 30-s time scale in human myelogenous leukemia cells. Because of its robust 3D measurement capabilities, live-cell optical CT represents a powerful new tool in the biomedical research field.

Epigenetic Transitions and Knotted Solitons in Stretched Chromatin

The spreading and regulation of epigenetic marks on chromosomes is crucial to establish and maintain cellular identity. Nonetheless, the dynamic mechanism leading to the establishment and maintenance of tissue-specific, epigenetic pattern is still poorly understood. In this work we propose, and investigate in silico, a possible experimental strategy to illuminate the interplay between 3D chromatin structure and epigenetic dynamics. We consider a set-up where a reconstituted chromatin fibre is stretched at its two ends (e.g., by laser tweezers), while epigenetic enzymes (writers) and chromatin-binding proteins (readers) are flooded into the system. We show that, by tuning the stretching force and the binding affinity of the readers for chromatin, the fibre undergoes a sharp transition between a stretched, epigenetically disordered, state and a crumpled, epigenetically coherent, one. We further investigate the case in which a knot is tied along the chromatin fibre, and find that the knotted segment enhances local epigenetic order, giving rise to "epigenetic solitons" which travel and diffuse along chromatin. Our results point to an intriguing coupling between 3D chromatin topology and epigenetic dynamics, which may be investigated via single molecule experiments.

Cell-to-cell variability and robustness in S-phase duration from genome replication kinetics

Genome replication, a key process for a cell, relies on stochastic initiation by replication origins, causing a variability of replication timing from cell to cell. While stochastic models of eukaryotic replication are widely available, the link between the key parameters and overall replication timing has not been addressed systematically. We use a combined analytical and computational approach to calculate how positions and strength of many origins lead to a given cell-to-cell variability of total duration of the replication of a large region, a chromosome or the entire genome. Specifically, the total replication timing can be framed as an extreme-value problem, since it is due to the last region that replicates in each cell. Our calculations identify two regimes based on the spread between characteristic completion times of all inter-origin regions of a genome. For widely different completion times, timing is set by the single specific region that is typically the last to replicate in all cells. Conversely, when the completion time of all regions are comparable, an extreme-value estimate shows that the cell-to-cell variability of genome replication timing has universal properties. Comparison with available data shows that the replication program of three yeast species falls in this extreme-value regime.

Replicating Large Genomes: Divide and Conquer

Complete duplication of large metazoan chromosomes requires thousands of potential initiation sites, only a small fraction of which are selected in each cell cycle. Assembly of the replication machinery is highly conserved and tightly regulated during the cell cycle, but the sites of initiation are highly flexible, and their temporal order of firing is regulated at the level of large-scale multi-replicon domains. Importantly, the number of replication forks must be quickly adjusted in response to replication stress to prevent genome instability. Here we argue that large genomes are divided into domains for exactly this reason. Once established, domain structure abrogates the need for precise initiation sites and creates a scaffold for the evolution of other chromosome functions.

Initial high-resolution microscopic mapping of active and inactive regulatory sequences proves non-random 3D arrangements in chromatin domain clusters

BACKGROUND

The association of active transcription regulatory elements (TREs) with DNAse I hypersensitivity (DHS[+]) and an 'open' local chromatin configuration has long been known. However, the 3D topography of TREs within the nuclear landscape of individual cells in relation to their active or inactive status has remained elusive. Here, we explored the 3D nuclear topography of active and inactive TREs in the context of a recently proposed model for a functionally defined nuclear architecture, where an active and an inactive nuclear compartment (ANC-INC) form two spatially co-aligned and functionally interacting networks.

RESULTS

Using 3D structured illumination microscopy, we performed 3D FISH with differently labeled DNA probe sets targeting either sites with DHS[+], apparently active TREs, or DHS[-] sites harboring inactive TREs. Using an in-house image analysis tool, DNA targets were quantitatively mapped on chromatin compaction shaped 3D nuclear landscapes. Our analyses present evidence for a radial 3D organization of chromatin domain clusters (CDCs) with layers of increasing chromatin compaction from the periphery to the CDC core. Segments harboring active TREs are significantly enriched at the decondensed periphery of CDCs with loops penetrating into interchromatin compartment channels, constituting the ANC. In contrast, segments lacking active TREs (DHS[-]) are enriched toward the compacted interior of CDCs (INC).

CONCLUSIONS

Our results add further evidence in support of the ANC-INC network model. The different 3D topographies of DHS[+] and DHS[-] sites suggest positional changes of TREs between the ANC and INC depending on their functional state, which might provide additional protection against an inappropriate activation. Our finding of a structural organization of CDCs based on radially arranged layers of different chromatin compaction levels indicates a complex higher-order chromatin organization beyond a dichotomic classification of chromatin into an 'open,' active and 'closed,' inactive state.

Multi-scale structural community organisation of the human genome

BACKGROUND

Structural interaction frequency matrices between all genome loci are now experimentally achievable thanks to high-throughput chromosome conformation capture technologies. This ensues a new methodological challenge for computational biology which consists in objectively extracting from these data the structural motifs characteristic of genome organisation.

RESULTS

We deployed the fast multi-scale community mining algorithm based on spectral graph wavelets to characterise the networks of intra-chromosomal interactions in human cell lines. We observed that there exist structural domains of all sizes up to chromosome length and demonstrated that the set of structural communities forms a hierarchy of chromosome segments. Hence, at all scales, chromosome folding predominantly involves interactions between neighbouring sites rather than the formation of links between distant loci.

CONCLUSIONS

Multi-scale structural decomposition of human chromosomes provides an original framework to question structural organisation and its relationship to functional regulation across the scales. By construction the proposed methodology is independent of the precise assembly of the reference genome and is thus directly applicable to genomes whose assembly is not fully determined.

Predicting chromatin architecture from models of polymer physics

We review the picture of chromatin large-scale 3D organization emerging from the analysis of Hi-C data and polymer modeling. In higher mammals, Hi-C contact maps reveal a complex higher-order organization, extending from the sub-Mb to chromosomal scales, hierarchically folded in a structure of domains-within-domains (metaTADs). The domain folding hierarchy is partially conserved throughout differentiation, and deeply correlated to epigenomic features. Rearrangements in the metaTAD topology relate to gene expression modifications: in particular, in neuronal differentiation models, topologically associated domains (TADs) tend to have coherent expression changes within architecturally conserved metaTAD niches. To identify the nature of architectural domains and their molecular determinants within a principled approach, we discuss models based on polymer physics. We show that basic concepts of interacting polymer physics explain chromatin spatial organization across chromosomal scales and cell types. The 3D structure of genomic loci can be derived with high accuracy and its molecular determinants identified by crossing information with epigenomic databases. In particular, we illustrate the case of the Sox9 locus, linked to human congenital disorders. The model in-silico predictions on the effects of genomic rearrangements are confirmed by available 5C data. That can help establishing new diagnostic tools for diseases linked to chromatin mis-folding, such as congenital disorders and cancer.

Deciphering DNA replication dynamics in eukaryotic cell populations in relation with their averaged chromatin conformations

We propose a non-local model of DNA replication that takes into account the observed uncertainty on the position and time of replication initiation in eukaryote cell populations. By picturing replication initiation as a two-state system and considering all possible transition configurations, and by taking into account the chromatin’s fractal dimension, we derive an analytical expression for the rate of replication initiation. This model predicts with no free parameter the temporal profiles of initiation rate, replication fork density and fraction of replicated DNA, in quantitative agreement with corresponding experimental data from both S. cerevisiae and human cells and provides a quantitative estimate of initiation site redundancy. This study shows that, to a large extent, the program that regulates the dynamics of eukaryotic DNA replication is a collective phenomenon that emerges from the stochastic nature of replication origins initiation.