The dynamics of genome replication using deep sequencing

The budding yeast Fkh1 Forkhead associated (FHA) domain promotes a G1-chromatin state and the activity of chromosomal DNA replication origins

In Saccharomyces cerevisiae, the forkhead (Fkh) transcription factor Fkh1 (forkhead homolog) enhances the activity of many DNA replication origins that act in early S-phase (early origins). Current models posit that Fkh1 acts directly to promote these origins' activity by binding to origin-adjacent Fkh1 binding sites (FKH sites). However, the post-DNA binding functions that Fkh1 uses to promote early origin activity are poorly understood. Fkh1 contains a conserved FHA (forkhead associated) domain, a protein-binding module with specificity for phosphothreonine (pT)-containing partner proteins. At a small subset of yeast origins, the Fkh1-FHA domain enhances the ORC (origin recognition complex)-origin binding step, the G1-phase event that initiates the origin cycle. However, the importance of the Fkh1-FHA domain to either chromosomal replication or ORC-origin interactions at genome scale is unclear. Here, S-phase SortSeq experiments were used to compare genome replication in proliferating FKH1 and fkh1-R80A mutant cells. The Fkh1-FHA domain promoted the activity of ≈ 100 origins that act in early to mid- S-phase, including the majority of centromere-associated origins, while simultaneously inhibiting ≈ 100 late origins. Thus, in the absence of a functional Fkh1-FHA domain, the temporal landscape of the yeast genome was flattened. Origins are associated with a positioned nucleosome array that frames a nucleosome depleted region (NDR) over the origin, and ORC-origin binding is necessary but not sufficient for this chromatin organization. To ask whether the Fkh1-FHA domain had an impact on this chromatin architecture at origins, ORC ChIPSeq data generated from proliferating cells and MNaseSeq data generated from G1-arrested and proliferating cell populations were assessed. Origin groups that were differentially regulated by the Fkh1-FHA domain were characterized by distinct effects of this domain on ORC-origin binding and G1-phase chromatin. Thus, the Fkh1-FHA domain controlled the distinct chromatin architecture at early origins in G1-phase and regulated origin activity in S-phase.

The budding yeast Fkh1 Forkhead associated (FHA) domain promoted a G1-chromatin state and the activity of chromosomal DNA replication origins

In Saccharomyces cerevisiae, the forkhead (Fkh) transcription factor Fkh1 (forkhead homolog) enhances the activity of many DNA replication origins that act in early S-phase (early origins). Current models posit that Fkh1 acts directly to promote these origins' activity by binding to origin-adjacent Fkh1 binding sites (FKH sites). However, the post-DNA binding functions that Fkh1 uses to promote early origin activity are poorly understood. Fkh1 contains a conserved FHA (forkhead associated) domain, a protein-binding module with specificity for phosphothreonine (pT)-containing partner proteins. At a small subset of yeast origins, the Fkh1-FHA domain enhances the ORC (origin recognition complex)-origin binding step, the G1-phase event that initiates the origin cycle. However, the importance of the Fkh1-FHA domain to either chromosomal replication or ORC-origin interactions at genome scale is unclear. Here, S-phase SortSeq experiments were used to compare genome replication in proliferating FKH1 and fkh1-R80A mutant cells. The Fkh1-FHA domain promoted the activity of 100 origins that act in early to mid- S-phase, including the majority of centromere-associated origins, while simultaneously inhibiting 100 late origins. Thus, in the absence of a functional Fkh1-FHA domain, the temporal landscape of the yeast genome was flattened. Origins are associated with a positioned nucleosome array that frames a nucleosome depleted region (NDR) over the origin, and ORC-origin binding is necessary but not sufficient for this chromatin organization. To ask whether the Fkh1-FHA domain had an impact on this chromatin architecture at origins, ORC ChIPSeq data generated from proliferating cells and MNaseSeq data generated from G1-arrested and proliferating cell populations were assessed. Origin groups that were differentially regulated by the Fkh1-FHA domain were characterized by distinct effects of this domain on ORC-origin binding and G1-phase chromatin. Thus, the Fkh1-FHA domain controlled the distinct chromatin architecture at early origins in G1-phase and regulated origin activity in S-phase.

Rationally designed chromosome fusion does not prevent rapid growth of Vibrio natriegens

DNA replication is essential for the proliferation of all cells. Bacterial chromosomes are replicated bidirectionally from a single origin of replication, with replication proceeding at about 1000 bp per second. For the model organism, Escherichia coli, this translates into a replication time of about 40 min for its 4.6 Mb chromosome. Nevertheless, E. coli can propagate by overlapping replication cycles with a maximum short doubling time of 20 min. The fastest growing bacterium known, Vibrio natriegens, is able to replicate with a generation time of less than 10 min. It has a bipartite genome with chromosome sizes of 3.2 and 1.9 Mb. Is simultaneous replication from two origins a prerequisite for its rapid growth? We fused the two chromosomes of V. natriegens to create a strain carrying one chromosome with a single origin of replication. Compared to the parental, this strain showed no significant deviation in growth rate. This suggests that the split genome is not a prerequisite for rapid growth.

Integrating high-throughput analysis to create an atlas of replication origins in Trypanosoma cruzi in the context of genome structure and variability

Trypanosoma cruzi is the etiologic agent of the most prevalent human parasitic disease in Latin America, Chagas disease. Its genome is rich in multigenic families that code for virulent antigens and are present in the rapidly evolving genomic compartment named Disruptive. DNA replication is a meticulous biological process in which flaws can generate mutations and changes in chromosomal and gene copy numbers. Here, integrating high-throughput and single-molecule analyses, we were able to identify Predominant, Flexible, and Dormant Orc1Cdc6-dependent origins as well as Orc1Cdc6-independent origins. Orc1Cdc6-dependent origins were found in multigenic family loci, while independent origins were found in the Core compartment that contains conserved and hypothetical protein-coding genes, in addition to multigenic families. In addition, we found that Orc1Cdc6 density is related to the firing of origins and that Orc1Cdc6-binding sites within fired origins are depleted of a specific class of nucleosomes that we previously categorized as dynamic. Together, these data suggest that Orc1Cdc6-dependent origins may contribute to the rapid evolution of the Disruptive compartment and, therefore, to the success of T. cruzi infection and that the local epigenome landscape is also involved in this process.IMPORTANCETrypanosoma cruzi, responsible for Chagas disease, affects millions globally, particularly in Latin America. Lack of vaccine or treatment underscores the need for research. Parasite's genome, with virulent antigen-coding multigenic families, resides in the rapidly evolving Disruptive compartment. Study sheds light on the parasite's dynamic DNA replication, discussing the evolution of the Disruptive compartment. Therefore, the findings represent a significant stride in comprehending T. cruzi's biology and the molecular bases that contribute to the success of infection caused by this parasite.

Sir2 and Fun30 regulate ribosomal DNA replication timing via Mcm helicase positioning and nucleosome occupancy

The association between late replication timing and low transcription rates in eukaryotic heterochromatin is well-known, yet the specific mechanisms underlying this link remain uncertain. In Saccharomyces cerevisiae, the histone deacetylase Sir2 is required for both transcriptional silencing and late replication at the repetitive ribosomal DNA arrays (rDNA). We have previously reported that in the absence of SIR2, a derepressed RNA PolII repositions MCM replicative helicases from their loading site at the ribosomal origin, where they abut well-positioned, high-occupancy nucleosomes, to an adjacent region with lower nucleosome occupancy. By developing a method that can distinguish activation of closely spaced MCM complexes, here we show that the displaced MCMs at rDNA origins have increased firing propensity compared to the non-displaced MCMs. Furthermore, we found that both, activation of the repositioned MCMs and low occupancy of the adjacent nucleosomes critically depend on the chromatin remodeling activity of FUN30. Our study elucidates the mechanism by which Sir2 delays replication timing, and it demonstrates, for the first time, that activation of a specific replication origin in vivo relies on the nucleosome context shaped by a single chromatin remodeler.

Genome replication in asynchronously growing microbial populations

Biological cells replicate their genomes in a well-planned manner. The DNA replication program of an organism determines the timing at which different genomic regions are replicated, with fundamental consequences for cell homeostasis and genome stability. In a growing cell culture, genomic regions that are replicated early should be more abundant than regions that are replicated late. This abundance pattern can be experimentally measured using deep sequencing. However, a general quantitative theory linking this pattern to the replication program is still lacking. In this paper, we predict the abundance of DNA fragments in asynchronously growing cultures from any given stochastic model of the DNA replication program. As key examples, we present stochastic models of the DNA replication programs in budding yeast and Escherichia coli. In both cases, our model results are in excellent agreement with experimental data and permit to infer key information about the replication program. In particular, our method is able to infer the locations of known replication origins in budding yeast with high accuracy. These examples demonstrate that our method can provide insight into a broad range of organisms, from bacteria to eukaryotes.

RNase H2 degrades toxic RNA:DNA hybrids behind stalled forks to promote replication restart

R-loops represent a major source of replication stress, but the mechanism by which these structures impede fork progression remains unclear. To address this question, we monitored fork progression, arrest, and restart in Saccharomyces cerevisiae cells lacking RNase H1 and H2, two enzymes responsible for degrading RNA:DNA hybrids. We found that while RNase H-deficient cells could replicate their chromosomes normally under unchallenged growth conditions, their replication was impaired when exposed to hydroxyurea (HU) or methyl methanesulfonate (MMS). Treated cells exhibited increased levels of RNA:DNA hybrids at stalled forks and were unable to generate RPA-coated single-stranded (ssDNA), an important postreplicative intermediate in resuming replication. Similar impairments in nascent DNA resection and ssDNA formation at HU-arrested forks were observed in human cells lacking RNase H2. However, fork resection was fully restored by addition of triptolide, an inhibitor of transcription that induces RNA polymerase degradation. Taken together, these data indicate that RNA:DNA hybrids not only act as barriers to replication forks, but also interfere with postreplicative fork repair mechanisms if not promptly degraded by RNase H.

The in vivo measurement of replication fork velocity and pausing by lag-time analysis

An important step towards understanding the mechanistic basis of the central dogma is the quantitative characterization of the dynamics of nucleic-acid-bound molecular motors in the context of the living cell. To capture these dynamics, we develop lag-time analysis, a method for measuring in vivo dynamics. Using this approach, we provide quantitative locus-specific measurements of fork velocity, in units of kilobases per second, as well as replisome pause durations, some with the precision of seconds. The measured fork velocity is observed to be both locus and time dependent, even in wild-type cells. In this work, we quantitatively characterize known phenomena, detect brief, locus-specific pauses at ribosomal DNA loci in wild-type cells, and observe temporal fork velocity oscillations in three highly-divergent bacterial species.

DNA read count calibration for single-molecule, long-read sequencing

There are many applications in which quantitative information about DNA mixtures with different molecular lengths is important. Gene therapy vectors are much longer than can be sequenced individually via short-read NGS. However, vector preparations may contain smaller DNAs that behave differently during sequencing. We have used two library preparations each for Pacific Biosystems (PacBio) and Oxford Nanopore Technologies NGS to determine their suitability for quantitative assessment of varying sized DNAs. Equimolar length standards were generated from E. coli genomic DNA. Both PacBio library preparations provided a consistent length dependence though with a complex pattern. This method is sufficiently sensitive that differences in genomic copy number between DNA from E. coli grown in exponential and stationary phase conditions could be detected. The transposase-based Oxford Nanopore library preparation provided a predictable length dependence, but the random sequence starts caused the loss of original length information. The ligation-based approach retained length information but read frequency was more variable. Modeling of E. coli versus lambda read frequency via cubic spline smoothing showed that the shorter genome could be used as a suitable internal spike-in for DNAs in the 200 bp to 10 kb range, allowing meaningful QC to be carried out with AAV preparations.

Systematic Analysis of Copy Number Variations in the Pathogenic Yeast Candida parapsilosis Identifies a Gene Amplification in RTA3 That is Associated with Drug Resistance

We analyzed the genomes of 170 C. parapsilosis isolates and identified multiple copy number variations (CNVs). We identified two genes, RTA3 (CPAR2_104610) and ARR3 (CPAR2_601050), each of which was the target of multiple independent amplification events. Phylogenetic analysis shows that most of these amplifications originated only once. For ARR3, which encodes a putative arsenate transporter, 8 distinct CNVs were identified, ranging in size from 2.3 kb to 10.5 kb with 3 to 23 copies. For RTA3, 16 distinct CNVs were identified, ranging in size from 0.3 kb to 4.5 kb with 2 to ~50 copies. One unusual amplification resulted in a DUP-TRP/INV-DUP structure similar to some human CNVs. RTA3 encodes a putative phosphatidylcholine (PC) floppase which is known to regulate the inward translocation of PC in Candida albicans. We found that an increased copy number of RTA3 correlated with resistance to miltefosine, an alkylphosphocholine drug that affects PC metabolism. Additionally, we conducted an adaptive laboratory evolution experiment in which two C. parapsilosis isolates were cultured in increasing concentrations of miltefosine. Two genes, CPAR2_303950 and CPAR2_102700, coding for putative PC flippases homologous to S. cerevisiae DNF1 gained homozygous protein-disrupting mutations in the evolved strains. Overall, our results show that C. parapsilosis can gain resistance to miltefosine, a drug that has recently been granted orphan drug designation approval by the United States Food and Drug Administration for the treatment of invasive candidiasis, through both CNVs or loss-of-function alleles in one of the flippase genes. IMPORTANCE Copy number variations (CNVs) are an important source of genomic diversity that have been associated with drug resistance. We identify two unusual CNVs in the human fungal pathogen Candida parapsilosis. Both target a single gene (RTA3 or ARR3), and they have occurred multiple times in multiple isolates. The copy number of RTA3, a putative floppase that controls the inward translocation of lipids in the cell membrane, correlates with resistance to miltefosine, a derivative of phosphatidylcholine (PC) that was originally developed as an anticancer drug. In 2021, miltefosine was designated an orphan drug by the United States Food and Drug Administration for the treatment of invasive candidiasis. Importantly, we find that resistance to miltefosine is also caused by mutations in flippases, which control the outward movement of lipids, and that many C. parapsilosis isolates are prone to easily acquiring an increased resistance to miltefosine.

Karyotype engineering reveals spatio-temporal control of replication firing and gene contacts

Eukaryotic genomes vary in terms of size, chromosome number, and genetic complexity. Their temporal organization is complex, reflecting coordination between DNA folding and function. Here, we used fused karyotypes of budding yeast to characterize the effects of chromosome length on nuclear architecture. We found that size-matched megachromosomes expand to occupy a larger fraction of the enlarged nucleus. Hi-C maps reveal changes in the three-dimensional structure corresponding to inactivated centromeres and telomeres. De-clustering of inactive centromeres results in their loss of early replication, highlighting a functional correlation between genome organization and replication timing. Repositioning of former telomere-proximal regions on chromosome arms exposed a subset of contacts between flocculin genes. Chromatin reorganization of megachromosomes during cell division remained unperturbed, and it revealed that centromere-rDNA contacts in anaphase, extending over 0.3 Mb on wild-type chromosome, cannot exceed ∼1.7 Mb. Our results highlight the relevance of engineered karyotypes to unveiling relationships between genome organization and function.

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.

Theulot et al. introduce NanoForkSpeed, a nanopore sequencing-based method to map individual replication fork velocities on entire genomes. NFS shows that fork speed is uniform across yeast chromosomes except for a marked slowdown at pausing sites.

Differential condensation of sister chromatids acts with Cdc6 to ensure asynchronous S-phase entry in Drosophila male germline stem cell lineage

During Drosophila melanogaster male germline stem cell (GSC) asymmetric division, preexisting old versus newly synthesized histones H3 and H4 are asymmetrically inherited. However, the biological outcomes of this phenomenon have remained unclear. Here, we tracked old and new histones throughout the GSC cell cycle through the use of high spatial and temporal resolution microscopy. We found unique features that differ between old and new histone-enriched sister chromatids, including differences in nucleosome density, chromosomal condensation, and H3 Ser10 phosphorylation. These distinct chromosomal features lead to their differential association with Cdc6, a pre-replication complex component, and subsequent asynchronous DNA replication initiation in the resulting daughter cells. Disruption of asymmetric histone inheritance abolishes differential Cdc6 association and asynchronous S-phase entry, demonstrating that histone asymmetry acts upstream of these critical cell-cycle progression events. Furthermore, disruption of these GSC-specific chromatin features leads to GSC defects, indicating a connection between histone inheritance, cell-cycle progression, and cell fate determination.

Genome-wide sequencing analysis of Sgs1, Exo1, Rad51, and Srs2 in DNA repair by homologous recombination

Homologous recombination is essential to maintain genome stability in response to DNA damage. Here, we have used genome-wide sequencing to quantitatively analyze at nucleotide resolution the dynamics of DNA end resection, re-synthesis, and gene conversion at a double-strand break. Resection initiates asymmetrically in an MRX-independent manner before proceeding steadily in both directions. Sgs1, Exo1, Rad51, and Srs2 differently regulate the rate and symmetry of early and late resection. Exo1 also ensures the coexistence of resection and re-synthesis, while Srs2 guarantees a constant and symmetrical DNA re-polymerization. Gene conversion is MMR independent, spans only a minor fraction of the resected region, and its unidirectionality depends on Srs2. Finally, these repair factors prevent the development of alterations remote from the DNA lesion, such as subtelomeric instability, duplication of genomic regions, and over-replication of Ty elements. Altogether, this approach allows a quantitative analysis and a direct genome-wide visualization of DNA repair by homologous recombination.

Cell-Cycle-Dependent Chromatin Dynamics at Replication Origins

Origins of DNA replication are specified by the ordered recruitment of replication factors in a cell-cycle–dependent manner. The assembly of the pre-replicative complex in G1 and the pre-initiation complex prior to activation in S phase are well characterized; however, the interplay between the assembly of these complexes and the local chromatin environment is less well understood. To investigate the dynamic changes in chromatin organization at and surrounding replication origins, we used micrococcal nuclease (MNase) to generate genome-wide chromatin occupancy profiles of nucleosomes, transcription factors, and replication proteins through consecutive cell cycles in Saccharomyces cerevisiae. During each G1 phase of two consecutive cell cycles, we observed the downstream repositioning of the origin-proximal +1 nucleosome and an increase in protected DNA fragments spanning the ARS consensus sequence (ACS) indicative of pre-RC assembly. We also found that the strongest correlation between chromatin occupancy at the ACS and origin efficiency occurred in early S phase, consistent with the rate-limiting formation of the Cdc45–Mcm2-7–GINS (CMG) complex being a determinant of origin activity. Finally, we observed nucleosome disruption and disorganization emanating from replication origins and traveling with the elongating replication forks across the genome in S phase, likely reflecting the disassembly and assembly of chromatin ahead of and behind the replication fork, respectively. These results provide insights into cell-cycle–regulated chromatin dynamics and how they relate to the regulation of origin activity.

Checkpoint-mediated DNA polymerase ε exonuclease activity curbing counteracts resection-driven fork collapse

DNA polymerase ε (Polε) carries out high-fidelity leading strand synthesis owing to its exonuclease activity. Polε polymerase and exonuclease activities are balanced, because of partitioning of nascent DNA strands between catalytic sites, so that net resection occurs when synthesis is impaired. In vivo, DNA synthesis stalling activates replication checkpoint kinases, which act to preserve the functional integrity of replication forks. We show that stalled Polε drives nascent strand resection causing fork functional collapse, averted via checkpoint-dependent phosphorylation. Polε catalytic subunit Pol2 is phosphorylated on serine 430, influencing partitioning between polymerase and exonuclease active sites. A phosphormimetic S430D change reduces exonucleolysis in vitro and counteracts fork collapse. Conversely, non-phosphorylatable pol2-S430A expression causes resection-driven stressed fork defects. Our findings reveal that checkpoint kinases switch Polε to an exonuclease-safe mode preventing nascent strand resection and stabilizing stalled replication forks. Elective partitioning suppression has implications for the diverse Polε roles in genome integrity maintenance.

The capacity of origins to load MCM establishes replication timing patterns

Loading of the MCM replicative helicase at origins of replication is a highly regulated process that precedes DNA replication in all eukaryotes. The stoichiometry of MCM loaded at origins has been proposed to be a key determinant of when those origins initiate replication during S phase. Nevertheless, the genome-wide regulation of MCM loading stoichiometry and its direct effect on replication timing remain unclear. In order to investigate why some origins load more MCM than others, we perturbed MCM levels in budding yeast cells and, for the first time, directly measured MCM levels and replication timing in the same experiment. Reduction of MCM levels through degradation of Mcm4, one of the six obligate components of the MCM complex, slowed progression through S phase and increased sensitivity to replication stress. Reduction of MCM levels also led to differential loading at origins during G1, revealing origins that are sensitive to reductions in MCM and others that are not. Sensitive origins loaded less MCM under normal conditions and correlated with a weak ability to recruit the origin recognition complex (ORC). Moreover, reduction of MCM loading at specific origins of replication led to a delay in their replication during S phase. In contrast, overexpression of MCM had no effects on cell cycle progression, relative MCM levels at origins, or replication timing, suggesting that, under optimal growth conditions, cellular MCM levels are not limiting for MCM loading. Our results support a model in which the loading capacity of origins is the primary determinant of MCM stoichiometry in wild-type cells, but that stoichiometry is controlled by origins' ability to recruit ORC and compete for MCM when MCM becomes limiting.

Genome-wide analysis of DNA replication and DNA double-strand breaks using TrAEL-seq

Faithful replication of the entire genome requires replication forks to copy large contiguous tracts of DNA, and sites of persistent replication fork stalling present a major threat to genome stability. Understanding the distribution of sites at which replication forks stall, and the ensuing fork processing events, requires genome-wide methods that profile replication fork position and the formation of recombinogenic DNA ends. Here, we describe Transferase-Activated End Ligation sequencing (TrAEL-seq), a method that captures single-stranded DNA 3' ends genome-wide and with base pair resolution. TrAEL-seq labels both DNA breaks and replication forks, providing genome-wide maps of replication fork progression and fork stalling sites in yeast and mammalian cells. Replication maps are similar to those obtained by Okazaki fragment sequencing; however, TrAEL-seq is performed on asynchronous populations of wild-type cells without incorporation of labels, cell sorting, or biochemical purification of replication intermediates, rendering TrAEL-seq far simpler and more widely applicable than existing replication fork direction profiling methods. The specificity of TrAEL-seq for DNA 3' ends also allows accurate detection of double-strand break sites after the initiation of DNA end resection, which we demonstrate by genome-wide mapping of meiotic double-strand break hotspots in a dmc1Δ mutant that is competent for end resection but not strand invasion. Overall, TrAEL-seq provides a flexible and robust methodology with high sensitivity and resolution for studying DNA replication and repair, which will be of significant use in determining mechanisms of genome instability.

Rtt109 slows replication speed by histone N-terminal acetylation

The wrapping of DNA around histone octamers challenges processes that use DNA as their template. In vitro, DNA replication through chromatin depends on histone modifiers, raising the possibility that cells modify histones to optimize fork progression. Rtt109 is an acetyl transferase that acetylates histone H3 before its DNA incorporation on the K56 and N-terminal residues. We previously reported that, in budding yeast, a wave of histone H3 K9 acetylation progresses ∼3–5 kb ahead of the replication fork. Whether this wave contributes to replication dynamics remained unknown. Here, we show that the replication fork velocity increases following deletion of RTT109, the gene encoding the enzyme required for the prereplication H3 acetylation wave. By using histone H3 mutants, we find that Rtt109-dependent N-terminal acetylation regulates fork velocity, whereas K56 acetylation contributes to replication dynamics only when N-terminal acetylation is compromised. We propose that acetylation of newly synthesized histones slows replication by promoting replacement of nucleosomes evicted by the incoming fork, thereby protecting genome integrity.

A practical guide to cancer subclonal reconstruction from DNA sequencing

Subclonal reconstruction from bulk tumor DNA sequencing has become a pillar of cancer evolution studies, providing insight into the clonality and relative ordering of mutations and mutational processes. We provide an outline of the complex computational approaches used for subclonal reconstruction from single and multiple tumor samples. We identify the underlying assumptions and uncertainties in each step and suggest best practices for analysis and quality assessment. This guide provides a pragmatic resource for the growing user community of subclonal reconstruction methods.

Genome duplication in Leishmania major relies on persistent subtelomeric DNA replication

DNA replication is needed to duplicate a cell’s genome in S phase and segregate it during cell division. Previous work in Leishmania detected DNA replication initiation at just a single region in each chromosome, an organisation predicted to be insufficient for complete genome duplication within S phase. Here, we show that acetylated histone H3 (AcH3), base J and a kinetochore factor co-localise in each chromosome at only a single locus, which corresponds with previously mapped DNA replication initiation regions and is demarcated by localised G/T skew and G4 patterns. In addition, we describe previously undetected subtelomeric DNA replication in G2/M and G1-phase-enriched cells. Finally, we show that subtelomeric DNA replication, unlike chromosome-internal DNA replication, is sensitive to hydroxyurea and dependent on 9-1-1 activity. These findings indicate that Leishmania’s genome duplication programme employs subtelomeric DNA replication initiation, possibly extending beyond S phase, to support predominantly chromosome-internal DNA replication initiation within S phase.

Isolation and analysis of rereplicated DNA by Rerep-Seq

Changes in gene copy number contribute to genomic instability, the onset and progression of cancer, developmental abnormalities and adaptive potential. The origins of gene amplifications have remained elusive; however, DNA rereplication has been implicated as a source of gene amplifications. The inability to determine which sequences are rereplicated and under what conditions have made it difficult to determine the validity of the proposed models. Here we present Rerep-Seq, a technique that selectively enriches for rereplicated DNA in preparation for analysis by DNA sequencing that can be applied to any species. We validated Rerep-Seq by simulating DNA rereplication in yeast and human cells. Using Rerep-Seq, we demonstrate that rereplication induced in Saccharomyces cerevisiae by deregulated origin licensing is non-random and defined by broad domains that span multiple replication origins and topological boundaries.

The role of Rif1 in telomere length regulation is separable from its role in origin firing

To examine the established link between DNA replication and telomere length, we tested whether firing of telomeric origins would cause telomere lengthening. We found that RIF1 mutants that block Protein Phosphatase 1 (PP1) binding activated telomeric origins but did not elongate telomeres. In a second approach, we found overexpression of ∆N-Dbf4 and Cdc7 increased DDK activity and activated telomeric origins, yet telomere length was unchanged. We tested a third mechanism to activate origins using the sld3-A mcm5-bob1 mutant that de-regulates the pre-replication complex, and again saw no change in telomere length. Finally, we tested whether mutations in RIF1 that cause telomere elongation would affect origin firing. We found that neither rif1-∆1322 nor rif1_HOOK affected firing of telomeric origins. We conclude that telomeric origin firing does not cause telomere elongation, and the role of Rif1 in regulating origin firing is separable from its role in regulating telomere length.

Sir2 mitigates an intrinsic imbalance in origin licensing efficiency between early- and late-replicating euchromatin

A eukaryotic chromosome relies on the function of multiple spatially distributed DNA replication origins for its stable inheritance. The spatial location of an origin is determined by the chromosomal position of an MCM complex, the inactive form of the DNA replicative helicase that is assembled onto DNA in G1-phase (also known as origin licensing). While the biochemistry of origin licensing is understood, the mechanisms that promote an adequate spatial distribution of MCM complexes across chromosomes are not. We have elucidated a role for the Sir2 histone deacetylase in establishing the normal distribution of MCM complexes across Saccharomyces cerevisiae chromosomes. In the absence of Sir2, MCM complexes accumulated within both early-replicating euchromatin and telomeric heterochromatin, and replication activity within these regions was enhanced. Concomitantly, the duplication of several regions of late-replicating euchromatin were delayed. Thus, Sir2-mediated attenuation of origin licensing within both euchromatin and telomeric heterochromatin established the normal spatial distribution of origins across yeast chromosomes important for normal genome duplication.

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.

DNA polymerase ε relies on a unique domain for efficient replisome assembly and strand synthesis

DNA polymerase epsilon (Pol ε) is required for genome duplication and tumor suppression. It supports both replisome assembly and leading strand synthesis; however, the underlying mechanisms remain to be elucidated. Here we report that a conserved domain within the Pol ε catalytic core influences both of these replication steps in budding yeast. Modeling cancer-associated mutations in this domain reveals its unexpected effect on incorporating Pol ε into the four-member pre-loading complex during replisome assembly. In addition, genetic and biochemical data suggest that the examined domain supports Pol ε catalytic activity and symmetric movement of replication forks. Contrary to previously characterized Pol ε cancer variants, the examined mutants cause genome hyper-rearrangement rather than hyper-mutation. Our work thus suggests a role of the Pol ε catalytic core in replisome formation, a reliance of Pol ε strand synthesis on a unique domain, and a potential tumor-suppressive effect of Pol ε in curbing genome re-arrangements.

Comparative Analysis of the Minimum Number of Replication Origins in Trypanosomatids and Yeasts

Single-celled eukaryote genomes predominantly replicate through multiple origins. Although origin usage during the S-phase has been elucidated in some of these organisms, few studies have comparatively approached this dynamic. Here, we developed a user-friendly website able to calculate the length of the cell cycle phases for any organism. Next, using a formula developed by our group, we showed a comparative analysis among the minimum number of replication origins (MO) required to duplicate an entire chromosome within the S-phase duration in trypanosomatids (Trypanosoma cruzi, Leishmania major, and Trypanosoma brucei) and yeasts (Saccharomyces cerevisiae and Schizosaccharomyces pombe). Using the data obtained by our analysis, it was possible to predict the MO required in a situation of replication stress. Also, our findings allow establishing a threshold for the number of origins, which serves as a parameter for genome approaches that map origins. Moreover, our data suggest that when compared to yeasts, trypanosomatids use much more origins than the minimum needed. This is the first time a comparative analysis of the minimum number of origins has been successfully applied. These data may provide new insight into the understanding of the replication mechanism and a new methodological framework for studying single-celled eukaryote genomes.

Budding yeast complete DNA synthesis after chromosome segregation begins

To faithfully transmit genetic information, cells must replicate their entire genome before division. This is thought to be ensured by the temporal separation of replication and chromosome segregation. Here we show that in 20–40% of unperturbed yeast cells, DNA synthesis continues during anaphase, late in mitosis. High cyclin-Cdk activity inhibits DNA synthesis in metaphase, and the decrease in cyclin-Cdk activity during mitotic exit allows DNA synthesis to finish at subtelomeric and some difficult-to-replicate regions. DNA synthesis during late mitosis correlates with elevated mutation rates at subtelomeric regions, including copy number variation. Thus, yeast cells temporally overlap DNA synthesis and chromosome segregation during normal growth, possibly allowing cells to maximize population-level growth rate while simultaneously exploring greater genetic space.

In the S phase of the cell cycle, the full genome needs to be replicated before cell division occurs. Here, authors show that in budding yeast DNA synthesis is completed after chromosome segregation begins.

Too Much of a Good Thing: How Ectopic DNA Replication Affects Bacterial Replication Dynamics

Each cell division requires the complete and accurate duplication of the entire genome. In bacteria, the duplication process of the often-circular chromosomes is initiated at a single origin per chromosome, resulting in two replication forks that traverse the chromosome in opposite directions. DNA synthesis is completed once the two forks fuse in a region diametrically opposite the origin. In some bacteria, such as Escherichia coli, the region where forks fuse forms a specialized termination area. Polar replication fork pause sites flanking this area can pause the progression of replication forks, thereby allowing forks to enter but not to leave. Transcription of all required genes has to take place simultaneously with genome duplication. As both of these genome trafficking processes share the same template, conflicts are unavoidable. In this review, we focus on recent attempts to add additional origins into various ectopic chromosomal locations of the E. coli chromosome. As ectopic origins disturb the native replichore arrangements, the problems resulting from such perturbations can give important insights into how genome trafficking processes are coordinated and the problems that arise if this coordination is disturbed. The data from these studies highlight that head-on replication–transcription conflicts are indeed highly problematic and multiple repair pathways are required to restart replication forks arrested at obstacles. In addition, the existing data also demonstrate that the replication fork trap in E. coli imposes significant constraints to genome duplication if ectopic origins are active. We describe the current models of how replication fork fusion events can cause serious problems for genome duplication, as well as models of how such problems might be alleviated both by a number of repair pathways as well as the replication fork trap system. Considering the problems associated both with head-on replication-transcription conflicts as well as head-on replication fork fusion events might provide clues of how these genome trafficking issues have contributed to shape the distinct architecture of bacterial chromosomes.

Mec1 Is Activated at the Onset of Normal S Phase by Low-dNTP Pools Impeding DNA Replication

The Mec1 and Rad53 kinases play a central role during acute replication stress in budding yeast. They are also essential for viability in normal growth conditions, but the signal that activates the Mec1-Rad53 pathway in the absence of exogenous insults is currently unknown. Here, we show that this pathway is active at the onset of normal S phase because deoxyribonucleotide triphosphate (dNTP) levels present in G₁ phase may not be sufficient to support processive DNA synthesis and impede DNA replication. This activation can be suppressed experimentally by increasing dNTP levels in G₁ phase. Moreover, we show that unchallenged cells entering S phase in the absence of Rad53 undergo irreversible fork collapse and mitotic catastrophe. Together, these data indicate that cells use suboptimal dNTP pools to detect the onset of DNA replication and activate the Mec1-Rad53 pathway, which in turn maintains functional forks and triggers dNTP synthesis, allowing the completion of DNA replication.

The Beacon Calculus: A formal method for the flexible and concise modelling of biological systems

Biological systems are made up of components that change their actions (and interactions) over time and coordinate with other components nearby. Together with a large state space, the complexity of this behaviour can make it difficult to create concise mathematical models that can be easily extended or modified. This paper introduces the Beacon Calculus, a process algebra designed to simplify the task of modelling interacting biological components. Its breadth is demonstrated by creating models of DNA replication dynamics, the gene expression dynamics in response to DNA methylation damage, and a multisite phosphorylation switch. The flexibility of these models is shown by adapting the DNA replication model to further include two topics of interest from the literature: cooperative origin firing and replication fork barriers. The Beacon Calculus is supported with the open-source simulator bcs (https://github.com/MBoemo/bcs.git) to allow users to develop and simulate their own models.

Simulating a model of a biological system can suggest ideas for future experiments and help ensure that conclusions about a mechanism are consistent with data. The Beacon Calculus is a new language that makes modelling simple by allowing users to simulate a biological system in only a few lines of code. This simplicity is critical as it allows users the freedom to come up with new ideas and rapidly test them. Models written in the Beacon Calculus are also easy to modify and extend, allowing users to add new features to the model or incorporate it into a larger biological system. We demonstrate the breadth of applications in this paper by applying the Beacon Calculus to DNA replication and DNA damage repair, both of which have implications for genome stability and cancer. We also apply it to multisite phosphorylation, which is important for cellular signalling. To enable users to create their own models, we created the open-source Beacon Calculus simulator bcs (https://github.com/MBoemo/bcs.git) which is easy to install and is well-supported by documentation and examples.

Genomic methods for measuring DNA replication dynamics

Genomic DNA replicates according to a defined temporal program in which early-replicating loci are associated with open chromatin, higher gene density, and increased gene expression levels, while late-replicating loci tend to be heterochromatic and show higher rates of genomic instability. The ability to measure DNA replication dynamics at genome scale has proven crucial for understanding the mechanisms and cellular consequences of DNA replication timing. Several methods, such as quantification of nucleotide analog incorporation and DNA copy number analyses, can accurately reconstruct the genomic replication timing profiles of various species and cell types. More recent developments have expanded the DNA replication genomic toolkit to assays that directly measure the activity of replication origins, while single-cell replication timing assays are beginning to reveal a new level of replication timing regulation. The combination of these methods, applied on a genomic scale and in multiple biological systems, promises to resolve many open questions and lead to a holistic understanding of how eukaryotic cells replicate their genomes accurately and efficiently.

DNA copy-number measurement of genome replication dynamics by high-throughput sequencing: the sort-seq, sync-seq and MFA-seq family

Genome replication follows a defined temporal programme that can change during cellular differentiation and disease onset. DNA replication results in an increase in DNA copy number that can be measured by high-throughput sequencing. Here we present a protocol to determine genome replication dynamics using DNA copy-number measurements. Cell populations can be obtained in three variants of the method. First, sort-seq reveals the average replication dynamics across S phase in an unperturbed cell population; FACS is used to isolate replicating and non-replicating subpopulations from asynchronous cells. Second, sync-seq measures absolute replication time at specific points during S phase using a synchronized cell population. Third, marker frequency analysis can be used to reveal the average replication dynamics using copy-number analysis in any proliferating asynchronous cell culture. These approaches have been used to reveal genome replication dynamics in prokaryotes, archaea and a wide range of eukaryotes, including yeasts and mammalian cells. We have found this approach straightforward to apply to other organisms and highlight example studies from across the three domains of life. Here we present a Saccharomyces cerevisiae version of the protocol that can be performed in 7-10 d. It requires basic molecular and cellular biology skills, as well as a basic understanding of Unix and R.

The evolutionary plasticity of chromosome metabolism allows adaptation to constitutive DNA replication stress

Many biological features are conserved and thus considered to be resistant to evolutionary change. While rapid genetic adaptation following the removal of conserved genes has been observed, we often lack a mechanistic understanding of how adaptation happens. We used the budding yeast, Saccharomyces cerevisiae, to investigate the evolutionary plasticity of chromosome metabolism, a network of evolutionary conserved modules. We experimentally evolved cells constitutively experiencing DNA replication stress caused by the absence of Ctf4, a protein that coordinates the enzymatic activities at replication forks. Parallel populations adapted to replication stress, over 1000 generations, by acquiring multiple, concerted mutations. These mutations altered conserved features of two chromosome metabolism modules, DNA replication and sister chromatid cohesion, and inactivated a third, the DNA damage checkpoint. The selected mutations define a functionally reproducible evolutionary trajectory. We suggest that the evolutionary plasticity of chromosome metabolism has implications for genome evolution in natural populations and cancer.

Esc2 promotes telomere stability in response to DNA replication stress

Telomeric regions of the genome are inherently difficult-to-replicate due to their propensity to generate DNA secondary structures and form nucleoprotein complexes that can impede DNA replication fork progression. Precisely how cells respond to DNA replication stalling within a telomere remains poorly characterized, largely due to the methodological difficulties in analysing defined stalling events in molecular detail. Here, we utilized a site-specific DNA replication barrier mediated by the ‘Tus/Ter’ system to define the consequences of DNA replication perturbation within a single telomeric locus. Through molecular genetic analysis of this defined fork-stalling event, coupled with the use of a genome-wide genetic screen, we identified an important role for the SUMO-like domain protein, Esc2, in limiting genome rearrangements at a telomere. Moreover, we showed that these rearrangements are driven by the combined action of the Mph1 helicase and the homologous recombination machinery. Our findings demonstrate that chromosomal context influences cellular responses to a stalled replication fork and reveal protective factors that are required at telomeric loci to limit DNA replication stress-induced chromosomal instability.

Checkpoint inhibition of origin firing prevents DNA topological stress

A universal feature of DNA damage and replication stress in eukaryotes is the activation of a checkpoint-kinase response. In S-phase, the checkpoint inhibits replication initiation, yet the function of this global block to origin firing remains unknown. To establish the physiological roles of this arm of the checkpoint, we analyzed separation of function mutants in the budding yeast Saccharomyces cerevisiae that allow global origin firing upon replication stress, despite an otherwise normal checkpoint response. Using genetic screens, we show that lack of the checkpoint-block to origin firing results in a dependence on pathways required for the resolution of topological problems. Failure to inhibit replication initiation indeed causes increased DNA catenation, resulting in DNA damage and chromosome loss. We further show that such topological stress is not only a consequence of a failed checkpoint response but also occurs in an unperturbed S-phase when too many origins fire simultaneously. Together we reveal that the role of limiting the number of replication initiation events is to prevent DNA topological problems, which may be relevant for the treatment of cancer with both topoisomerase and checkpoint inhibitors.

The effects of manipulating levels of replication initiation factors on origin firing efficiency in yeast

Chromosome replication in Saccharomyces cerevisiae is initiated from ~300 origins that are regulated by DNA sequence and by the limited abundance of six trans-acting initiation proteins (Sld2, Sld3, Dpb11, Dbf4, Sld7 and Cdc45). We set out to determine how the levels of individual factors contribute to time of origin activation and/or origin efficiency using induced depletion of single factors and overexpression of sets of multiple factors. Depletion of Sld2 or Sld3 slows growth and S phase progression, decreases origin efficiency across the genome and impairs viability as a result of incomplete replication of the rDNA. We find that the most efficient early origins are relatively unaffected by depletion of either Sld2 or Sld3. However, Sld3 levels, and to a lesser extent Sld2 levels, are critical for firing of the less efficient early origins. Overexpression of Sld3 simultaneously with Sld2, Dpb11 and Dbf4 preserves the relative efficiency of origins. Only when Cdc45 and Sld7 are also overexpressed is origin efficiency equalized between early- and late-firing origins. Our data support a model in which Sld3 together with Cdc45 (and/or Sld7) is responsible for the differential efficiencies of origins across the yeast genome.

Eukaryotic chromosome duplication begins at sites called origins of replication along the chromosomal DNA. A conserved property of eukaryotic origins is that they vary in efficiency—the proportion of cells in a population in which they “fire”—and in the average time of activation within S phase, but the molecular details underlying this variation are not well understood. Previous work has shown that limiting concentrations of a set of conserved replication initiation proteins referred to as “SSDDCS” (Sld2, Sld3, Dbf4, Dpb11, Cdc45, and Sld7) are rate limiting for origin activation in budding yeast and other eukaryotes; combined overexpression of these proteins increases and/or advances origin firing. However, it remained possible that different factors affect different aspects of origin activation (e.g., timing vs. efficiency). Here, by depleting individual factors or by overexpressing sets of factors in budding yeast, we demonstrate that it is levels of Sld3, Cdc45 and/or Sld7 levels are primarily responsible for modulating the differences in relative origin efficiency and timing. This work gives further insights into what shapes the landscape of genome duplication.

Rif1 Binding and Control of Chromosome-Internal DNA Replication Origins Is Limited by Telomere Sequestration

The Saccharomyces cerevisiae telomere-binding protein Rif1 plays an evolutionarily conserved role in control of DNA replication timing by promoting PP1-dependent dephosphorylation of replication initiation factors. However, ScRif1 binding outside of telomeres has never been detected, and it has thus been unclear whether Rif1 acts directly on the replication origins that it controls. Here, we show that, in unperturbed yeast cells, Rif1 primarily regulates late-replicating origins within 100 kb of a telomere. Using the chromatin endogenous cleavage ChEC-seq technique, we robustly detect Rif1 at late-replicating origins that we show are targets of its inhibitory action. Interestingly, abrogation of Rif1 telomere association by mutation of its Rap1-binding module increases Rif1 binding and origin inhibition elsewhere in the genome. Our results indicate that Rif1 inhibits replication initiation by interacting directly with origins and suggest that Rap1-dependent sequestration of Rif1 increases its effective concentration near telomeres, while limiting its action at chromosome-internal sites.

Comprehensive Analysis of Replication Origins in Saccharomyces cerevisiae Genomes

DNA replication initiates from multiple replication origins (ORIs) in eukaryotes. Discovery and characterization of replication origins are essential for a better understanding of the molecular mechanism of DNA replication. In this study, the features of autonomously replicating sequences (ARSs) in Saccharomyces cerevisiae have been comprehensively analyzed as follows. Firstly, we carried out the analysis of the ARSs available in S. cerevisiae S288C. By evaluating the sequence similarity of experimentally established ARSs, we found that 94.32% of ARSs are unique across the whole genome of S. cerevisiae S288C and those with high sequence similarity are prone to locate in subtelomeres. Subsequently, we built a non-redundant dataset with a total of 520 ARSs, which are based on ARSs annotation of S. cerevisiae S288C from SGD and then supplemented with those from OriDB and DeOri databases. We conducted a large-scale comparison of ORIs among the diverse budding yeast strains from a population genomics perspective. We found that 82.7% of ARSs are not only conserved in genomic sequence but also relatively conserved in chromosomal position. The non-conserved ARSs tend to distribute in the subtelomeric regions. We also conducted a pan-genome analysis of ARSs among the S. cerevisiae strains, and a total of 183 core ARSs existing in all yeast strains were determined. We extracted the genes adjacent to replication origins among the 104 yeast strains to examine whether there are differences in their gene functions. The result showed that the genes involved in the initiation of DNA replication, such as orc3, mcm2, mcm4, mcm6, and cdc45, are conservatively located adjacent to the replication origins. Furthermore, we found the genes adjacent to conserved ARSs are significantly enriched in DNA binding, enzyme activity, transportation, and energy, whereas for the genes adjacent to non-conserved ARSs are significantly enriched in response to environmental stress, metabolites biosynthetic process and biosynthesis of antibiotics. In general, we characterized the replication origins from the genome-wide and population genomics perspectives, which would provide new insights into the replication mechanism of S. cerevisiae and facilitate the design of algorithms to identify genome-wide replication origins in yeast.

Roles for DNA polymerase δ in initiating and terminating leading strand DNA replication

Most current evidence indicates that DNA polymerases ε and δ, respectively, perform the bulk of leading and lagging strand replication of the eukaryotic nuclear genome. Given that ribonucleotide and mismatch incorporation rates by these replicases influence somatic and germline patterns of variation, it is important to understand the details and exceptions to this overall division of labor. Using an improved method to map where these replicases incorporate ribonucleotides during replication, here we present evidence that DNA polymerase δ universally participates in initiating leading strand synthesis and that nascent leading strand synthesis switches from Pol ε to Pol δ during replication termination. Ribonucleotide maps from both the budding and fission yeast reveal conservation of these processes. These observations of replisome dynamics provide important insight into the mechanisms of eukaryotic replication and genome maintenance.

Condensin-Mediated Chromosome Folding and Internal Telomeres Drive Dicentric Severing by Cytokinesis

In Saccharomyces cerevisiae, dicentric chromosomes stemming from telomere fusions preferentially break at the fusion. This process restores a normal karyotype and protects chromosomes from the detrimental consequences of accidental fusions. Here, we address the molecular basis of this rescue pathway. We observe that tandem arrays tightly bound by the telomere factor Rap1 or a heterologous high-affinity DNA binding factor are sufficient to establish breakage hotspots, mimicking telomere fusions within dicentrics. We also show that condensins generate forces sufficient to rapidly refold dicentrics prior to breakage by cytokinesis and are essential to the preferential breakage at telomere fusions. Thus, the rescue of fused telomeres results from a condensin- and Rap1-driven chromosome folding that favors fusion entrapment where abscission takes place. Because a close spacing between the DNA-bound Rap1 molecules is essential to this process, Rap1 may act by stalling condensins.

Capturing the dynamics of genome replication on individual ultra-long nanopore sequence reads

Replication of eukaryotic genomes is highly stochastic, making it difficult to determine the replication dynamics of individual molecules with existing methods. We report a sequencing method for the measurement of replication fork movement on single molecules by detecting nucleotide analog signal currents on extremely long nanopore traces (D-NAscent). Using this method, we detect 5-bromodeoxyuridine (BrdU) incorporated by Saccharomyces cerevisiae to reveal, at a genomic scale and on single molecules, the DNA sequences replicated during a pulse-labeling period. Under conditions of limiting BrdU concentration, D-NAscent detects the differences in BrdU incorporation frequency across individual molecules to reveal the location of active replication origins, fork direction, termination sites, and fork pausing/stalling events. We used sequencing reads of 20-160 kilobases to generate a whole-genome single-molecule map of DNA replication dynamics and discover a class of low-frequency stochastic origins in budding yeast. The D-NAscent software is available at https://github.com/MBoemo/DNAscent.git .

Bayesian inference of origin firing time distributions, origin interference and licencing probabilities from Next Generation Sequencing data

DNA replication is a stochastic process with replication forks emanating from multiple replication origins. The origins must be licenced in G1, and the replisome activated at licenced origins in order to generate bi-directional replication forks in S-phase. Differential firing times lead to origin interference, where a replication fork from an origin can replicate through and inactivate neighbouring origins (origin obscuring). We developed a Bayesian algorithm to characterize origin firing statistics from Okazaki fragment (OF) sequencing data. Our algorithm infers the distributions of firing times and the licencing probabilities for three consecutive origins. We demonstrate that our algorithm can distinguish partial origin licencing and origin obscuring in OF sequencing data from Saccharomyces cerevisiae and human cell types. We used our method to analyse the decreased origin efficiency under loss of Rat1 activity in S. cerevisiae, demonstrating that both reduced licencing and increased obscuring contribute. Moreover, we show that robust analysis is possible using only local data (across three neighbouring origins), and analysis of the whole chromosome is not required. Our algorithm utilizes an approximate likelihood and a reversible jump sampling technique, a methodology that can be extended to analysis of other mechanistic processes measurable through Next Generation Sequencing data.

DNA Replication Timing Enters the Single-Cell Era

In mammalian cells, DNA replication timing is controlled at the level of megabase (Mb)-sized chromosomal domains and correlates well with transcription, chromatin structure, and three-dimensional (3D) genome organization. Because of these properties, DNA replication timing is an excellent entry point to explore genome regulation at various levels and a variety of studies have been carried out over the years. However, DNA replication timing studies traditionally required at least tens of thousands of cells, and it was unclear whether the replication domains detected by cell population analyses were preserved at the single-cell level. Recently, single-cell DNA replication profiling methods became available, which revealed that the Mb-sized replication domains detected by cell population analyses were actually well preserved in individual cells. In this article, we provide a brief overview of our current knowledge on DNA replication timing regulation in mammals based on cell population studies, outline the findings from single-cell DNA replication profiling, and discuss future directions and challenges.

A transcriptome analysis of the ameliorate effect of Cyclocarya paliurus triterpenoids on ethanol stress in Saccharomyces cerevisiae

Saccharomyces cerevisiae (S. cerevisiae) plays a critical role in ethanol fermentation. However, during the fermentation, yeast cells are exposed to the accumulation of ethanol, which significantly affect the cell growth and the target product yield. In the present work, we employed RNA-sequence (RNA-seq) to investigate the ameliorate effect of Cyclocarya paliurus (C. paliurus) triterpenoids on S. cerevisiae under the ethanol stress. After C. paliurus triterpenoids intervention (0.3% v/v), 84 differentially expressed genes (DEGs) were identified, including 39 up-regulated and 45 down-regulated genes. The addition of triterpenoids decreased the filamentous and invasive growth of cells, and benefit to the redox balance and glycolysis. This study offers a global view through transcriptome analysis to understand the molecular response to ethanol in Sc131 by the treatment of C. paliurus triterpenoids, which may be helpful to enhance ethanol tolerance of S. cerevisiae in the fermentation of Chinese fruit wine.

ChECing out Rif1 action in freely cycling cells

In buddying yeast, like all eukaryotes examined so far, DNA replication is under temporal control, such that some origins fire early and some late during S phase. This replication timing program is established in G1 phase, where chromatin states are thought to prevent binding of key-limiting initiation factors at late-firing origins. Although many factors are involved in replication initiation, a new player, Rif1, has recently entered the scene, with a spate of papers revealing a global role for the protein in the control of replication initiation timing from yeasts to humans. Since budding yeast Rif1 was known to bind only to telomeric and silent mating loci regions, it remained controversial whether Rif1 acts directly at replication origins or instead influences origin activity indirectly. In this perspective, we discuss our recent finding that Rif1 binds directly to the replication origins that it controls. In this study, we also found that Rif1's regulatory activity at origins is best revealed by an assay (sort-seq) that measures replication in unperturbed, freely cycling cultures, as opposed to commonly used protocols in which cells are first blocked in the G1 phase of the cell cycle by mating pheromone, then released into a synchronous S phase. Finally, we discuss how the sequestration of Rif1 at telomeres, through an interaction with the arrays of Rap1 molecules bound there, plays an important role in limiting Rif1's action primarily to telomere-proximal replication origins.

Mrc1 and Rad9 cooperate to regulate initiation and elongation of DNA replication in response to DNA damage

The S-phase checkpoint maintains the integrity of the genome in response to DNA replication stress. In budding yeast, this pathway is initiated by Mec1 and is amplified through the activation of Rad53 by two checkpoint mediators: Mrc1 promotes Rad53 activation at stalled forks, and Rad9 is a general mediator of the DNA damage response. Here, we have investigated the interplay between Mrc1 and Rad9 in response to DNA damage and found that they control DNA replication through two distinct but complementary mechanisms. Mrc1 rapidly activates Rad53 at stalled forks and represses late-firing origins but is unable to maintain this repression over time. Rad9 takes over Mrc1 to maintain a continuous checkpoint signaling. Importantly, the Rad9-mediated activation of Rad53 slows down fork progression, supporting the view that the S-phase checkpoint controls both the initiation and the elongation of DNA replication in response to DNA damage. Together, these data indicate that Mrc1 and Rad9 play distinct functions that are important to ensure an optimal completion of S phase under replication stress conditions.

Rapid high-resolution measurement of DNA replication timing by droplet digital PCR

Genomes are replicated in a reproducible temporal pattern. Current methods for assaying allele replication timing are time consuming and/or expensive. These include high-throughput sequencing which can be used to measure DNA copy number as a proxy for allele replication timing. Here, we use droplet digital PCR to study DNA replication timing at multiple loci in budding yeast and human cells. We establish that the method has temporal and spatial resolutions comparable to the high-throughput sequencing approaches, while being faster than alternative locus-specific methods. Furthermore, the approach is capable of allele discrimination. We apply this method to determine relative replication timing across timing transition zones in cultured human cells. Finally, multiple samples can be analysed in parallel, allowing us to rapidly screen kinetochore mutants for perturbation to centromere replication timing. Therefore, this approach is well suited to the study of locus-specific replication and the screening of cis- and trans-acting mutants to identify mechanisms that regulate local genome replication timing.

Genomic Analysis of DNA Double-Strand Break Repair in Escherichia coli

Counting DNA whole genome sequencing reads is providing new insight into DNA double-strand break repair (DSBR) in the model organism Escherichia coli. We describe the application of RecA chromatin immunoprecipitation coupled to genomic DNA sequencing (RecA-ChIP-seq) and marker frequency analysis (MFA) to analyze the genomic consequences of DSBR. We provide detailed procedures for the preparation of DNA and the analysis of data. We compare different ways of visualizing ChIP data and show that alternative protocols for the preparation of DNA for MFA differentially affect the recovery of branched DNA molecules containing Holliday junctions.