Neural network and kinetic modelling of human genome replication reveal replication origin locations and strengths.

In human and other metazoans, the determinants of replication origin location and strength are still elusive. Origins are licensed in G1 phase and fired in S phase of the cell cycle, respectively. It is debated which of these two temporally separate steps determines origin efficiency. Experiments can independently profile mean replication timing (MRT) and replication fork directionality (RFD) genome-wide. Such profiles contain information on multiple origins' properties and on fork speed. Due to possible origin inactivation by passive replication, however, observed and intrinsic origin efficiencies can markedly differ. Thus, there is a need for methods to infer intrinsic from observed origin efficiency, which is context-dependent. Here, we show that MRT and RFD data are highly consistent with each other but contain information at different spatial scales. Using neural networks, we infer an origin licensing landscape that, when inserted in an appropriate simulation framework, jointly predicts MRT and RFD data with unprecedented precision and underlies the importance of dispersive origin firing. We furthermore uncover an analytical formula that predicts intrinsic from observed origin efficiency combined with MRT data. Comparison of inferred intrinsic origin efficiencies with experimental profiles of licensed origins (ORC, MCM) and actual initiation events (Bubble-seq, SNS-seq, OK-seq, ORM) show that intrinsic origin efficiency is not solely determined by licensing efficiency. Thus, human replication origin efficiency is set at both the origin licensing and firing steps.

Chromatin stiffening underlies enhanced locus mobility after DNA damage in budding yeast.

DNA double-strand breaks (DSBs) induce a cellular response that involves histone modifications and chromatin remodeling at the damaged site and increases chromosome dynamics both locally at the damaged site and globally in the nucleus. In parallel, it has become clear that the spatial organization and dynamics of chromosomes can be largely explained by the statistical properties of tethered, but randomly moving, polymer chains, characterized mainly by their rigidity and compaction. How these properties of chromatin are affected during DNA damage remains, however, unclear. Here, we use live cell microscopy to track chromatin loci and measure distances between loci on yeast chromosome IV in thousands of cells, in the presence or absence of genotoxic stress. We confirm that DSBs result in enhanced chromatin subdiffusion and show that intrachromosomal distances increase with DNA damage all along the chromosome. Our data can be explained by an increase in chromatin rigidity, but not by chromatin decondensation or centromeric untethering only. We provide evidence that chromatin stiffening is mediated in part by histone H2A phosphorylation. Our results support a genome-wide stiffening of the chromatin fiber as a consequence of DNA damage and as a novel mechanism underlying increased chromatin mobility.

Cooperativity in the annealing of DNA origamis.

DNA based nanostructures built on a long single stranded DNA scaffold, known as DNA origamis, offer the possibility to organize various molecules at the nanometer scale in one pot experiments. The folding of the scaffold is guaranteed by the presence of short, single stranded DNA sequences (staples), that hold together separate regions of the scaffold. In this paper, we modelize the annealing-melting properties of these DNA constructions. The model captures important features such as the hysteresis between melting and annealing, as well as the dependence upon the topology of the scaffold. We show that cooperativity between staples is critical to quantitatively explain the folding process of DNA origamis.

Direct visualization of transient thermal response of a DNA origami.

The DNA origami approach enables the construction of complex objects from DNA strands. A fundamental understanding of the kinetics and thermodynamics of DNA origami assembly is extremely important for building large DNA structures with multifunctionality. Here both experimental and theoretical studies of DNA origami melting were carried out in order to reveal the reversible association/disassociation process. Furthermore, by careful control of the temperature cycling via in situ thermally controlled atomic force microscopy, the self-assembly process of a rectangular DNA origami tile was directly visualized, unveiling key mechanisms underlying their structural and thermodynamic features.

Folding of small origamis.

A model that preserves the known thermodynamic properties of double stranded DNA is introduced to study the formation of more complex DNA constructions, such as small origamis or Holliday junctions. We show that the thermodynamic behaviour of these complex DNA constructions is not only given by their sequence but also by their topology.

FORK-seq: Single-Molecule Profiling of DNA Replication.

Most genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200 nucleotide resolution using a nanopore current interpretation tool allowing the quantification of BrdU incorporation. Along pulse-chased replication intermediates from Saccharomyces cerevisiae, we can orient replication tracks and reproduce population-based replication directionality profiles. Additionally, we can map individual initiation and termination events. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

Genome-wide mapping of individual replication fork velocities using nanopore sequencing.

Little is known about replication fork velocity variations along eukaryotic genomes, since reference techniques to determine fork speed either provide no sequence information or suffer from low throughput. Here we present NanoForkSpeed, a nanopore sequencing-based method to map and extract the velocity of individual forks detected as tracks of the thymidine analogue bromodeoxyuridine incorporated during a brief pulse-labelling of asynchronously growing cells. NanoForkSpeed retrieves previous Saccharomyces cerevisiae mean fork speed estimates (≈2 kb/min) in the BT1 strain exhibiting highly efficient bromodeoxyuridine incorporation and wild-type growth, and precisely quantifies speed changes in cells with altered replisome progression or exposed to hydroxyurea. The positioning of >125,000 fork velocities provides a genome-wide map of fork progression based on individual fork rates, showing a uniform fork speed across yeast chromosomes except for a marked slowdown at known pausing sites.

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.

FORK-seq: replication landscape of the Saccharomyces cerevisiae genome by nanopore sequencing.

Genome replication mapping methods profile cell populations, masking cell-to-cell heterogeneity. Here, we describe FORK-seq, a nanopore sequencing method to map replication of single DNA molecules at 200-nucleotide resolution. By quantifying BrdU incorporation along pulse-chased replication intermediates from Saccharomyces cerevisiae, we orient 58,651 replication tracks reproducing population-based replication directionality profiles and map 4964 and 4485 individual initiation and termination events, respectively. Although most events cluster at known origins and fork merging zones, 9% and 18% of initiation and termination events, respectively, occur at many locations previously missed. Thus, FORK-seq reveals the full extent of cell-to-cell heterogeneity in DNA replication.

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

The time-dependent rate [Formula: see text] 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 [Formula: see text] 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 [Formula: see text] 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.

Inferring the physical properties of yeast chromatin through Bayesian analysis of whole nucleus simulations.

BACKGROUND: The structure and mechanical properties of chromatin impact DNA functions and nuclear architecture but remain poorly understood. In budding yeast, a simple polymer model with minimal sequence-specific constraints and a small number of structural parameters can explain diverse experimental data on nuclear architecture. However, how assumed chromatin properties affect model predictions was not previously systematically investigated. RESULTS: We used hundreds of dynamic chromosome simulations and Bayesian inference to determine chromatin properties consistent with an extensive dataset that includes hundreds of measurements from imaging in fixed and live cells and two Hi-C studies. We place new constraints on average chromatin fiber properties, narrowing down the chromatin compaction to ~53-65 bp/nm and persistence length to ~52-85 nm. These constraints argue against a 20-30 nm fiber as the exclusive chromatin structure in the genome. Our best model provides a much better match to experimental measurements of nuclear architecture and also recapitulates chromatin dynamics measured on multiple loci over long timescales. CONCLUSION: This work substantially improves our understanding of yeast chromatin mechanics and chromosome architecture and provides a new analytic framework to infer chromosome properties in other organisms.

Connecting localized DNA strand displacement reactions.

Logic circuits based on DNA strand displacement reactions have been shown to be versatile enough to compute the square root of four-bit numbers. The implementation of these circuits as a set of bulk reactions faces difficulties which include leaky reactions and intrinsically slow, diffusion-limited reaction rates. In this paper, we consider simple examples of these circuits when they are attached to platforms (DNA origamis). As expected, constraining distances between DNA strands leads to faster reaction rates. However, it also induces side-effects that are not detectable in the solution-phase version of this circuitry. Appropriate design of the system, including protection and asymmetry between input and fuel strands, leads to a reproducible behaviour, at least one order of magnitude faster than the one observed under bulk conditions.

How to build a yeast nucleus.

Biological functions including gene expression and DNA repair are affected by the 3D architecture of the genome, but the underlying mechanisms are still unknown. Notably, it remains unclear to what extent nuclear architecture is driven by generic physical properties of polymers or by specific factors such as proteins binding particular DNA sequences. The budding yeast nucleus has been intensely studied by imaging and biochemical techniques, resulting in a large quantitative data set on locus positions and DNA contact frequencies. We recently described a quantitative model of the interphase yeast nucleus in which chromosomes are represented as passively moving polymer chains. This model ignores the DNA sequence information except for specific constraints at the centromeres, telomeres, and the ribosomal DNA (rDNA). Despite its simplicity, the model accounts for a large majority of experimental data, including absolute and relative locus positions and contact frequency patterns at chromosomal and subchromosomal scales. Here, we also illustrate the model's ability to reproduce observed features of chromatin movements. Our results strongly suggest that the dynamic large-scale architecture of the yeast nucleus is dominated by statistical properties of randomly moving polymers with a few sequence-specific constraints, rather than by a large number of DNA-specific factors or epigenetic modifications. In addition, we show that our model accounts for recently measured variations in homologous recombination efficiency, illustrating its potential for quantitatively understanding functional consequences of nuclear architecture.

Guided assembly of tetramolecular G-quadruplexes.

Nucleic acids are finding applications in nanotechnology as nanomaterials, mechanical devices, templates, and biosensors. G-quadruplex DNA, formed by π-π stacking of guanine (G) quartets, is an attractive alternative to regular B-DNA because of the kinetic and thermodynamic stability of quadruplexes. However, they suffer from a fatal flaw: the rules of recognition, i.e., the formation of a G-quartet in which four identical bases are paired, prevent the controlled assembly between different strands, leading to complex mixtures. In this report, we present the solution to this recognition problem. The proposed design combines two DNA elements: parallel-stranded duplexes and a quadruplex core. Parallel-stranded duplexes direct controlled assembly of the quadruplex core, and their strands present convenient points of attachments for potential modifiers. The exceptional stability of the quadruplex core provides integrity to the entire structure, which could be used as a building block for nucleic acid-based nanomaterials. As a proof of principle for the design's versatility, we assembled quadruplex-based 1D structures and visualized them using atomic force and transmission electron microscopy. Our findings pave the way to broader utilization of G-quadruplex DNA in structural DNA nanomaterials.

Modeling the mechanical properties of DNA nanostructures.

We discuss generalizations of a previously published coarse-grained description [Mergell et al., Phys. Rev. E 68, 021911 (2003)] of double stranded DNA (dsDNA). The model is defined at the base-pair level and includes the electrostatic repulsion between neighbor helices. We show that the model reproduces mechanical and elastic properties of several DNA nanostructures (DNA origamis). We also show that electrostatic interactions are necessary to reproduce atomic force microscopy measurements on planar DNA origamis.