Replication Initiation Sites Research Articles

In vivo DNA replication dynamics unveil aging-dependent replication stress

The genome duplication program is affected by multiple factors invivo, including developmental cues, genotoxic stress, and aging. Here, we monitored DNA replication initiation dynamics in regenerating livers of young and old mice after partial hepatectomy to investigate the impact of aging. In young mice, the origin firing sites were well defined; the majority were located 10-50 kb upstream or downstream of expressed genes, and their position on the genome was conserved in human cells. Old mice displayed the same replication initiation sites, but origin firing was inefficient and accompanied by a replication stress response. Inhibitors of the ATR checkpoint kinase fully restored origin firing efficiency in the old mice but at the expense of an inflammatory response and without significantly enhancing the fraction of hepatocytes entering the cell cycle. These findings unveil aging-dependent replication stress and a crucial role of ATR in mitigating the stress-associated inflammation, a hallmark of aging.

Replication initiation sites and zones in the mammalian genome: Where are they located and how are they defined?

Eukaryotic DNA replication is a tightly controlled process that occurs in two main steps, i.e., licensing and firing, which take place in the G1 and S phases of the cell cycle, respectively. In Saccharomyces cerevisiae, the budding yeast, replication origins contain consensus sequences that are recognized and bound by the licensing factor Orc1–6, which then recruits the replicative Mcm2–7 helicase. By contrast, mammalian initiation sites lack such consensus sequences, and the mammalian ORC does not exhibit sequence specificity. Studies performed over the past decades have identified replication initiation sites in the mammalian genome using sequencing-based assays, raising the question of whether replication initiation occurs at confined sites or in broad zones across the genome. Although recent reports have shown that the licensed MCMs in mammalian cells are broadly distributed, suggesting that ORC-dependent licensing may not determine the initiation sites/zones, they are predominantly located upstream of actively transcribed genes. This review compares the mechanism of replication initiation in yeast and mammalian cells, summarizes the sequencing-based technologies used for the identification of initiation sites/zones, and proposes a possible mechanism of initiation-site/zone selection in mammalian cells. Future directions and challenges in this field are also discussed.

Genome-wide binding sites of Plasmodium falciparum mini chromosome maintenance protein MCM6 show new insights into parasite DNA replication

Multiple rounds of DNA replication take place in various stages of the life cycle in the human malaria parasite Plasmodium falciparum. Previous bioinformatics analysis has shown the presence of putative Autonomously Replicating Sequence (ARS) like sequences in the Plasmodium genome. However, the actual sites and frequency of replication origins in the P. falciparum genome based on experimental data still remain elusive. Minichromosome maintenance (MCM) proteins are recruited by the Origin recognition complex (ORC) to the origins of replication in eukaryotes including P. falciparum. We used PfMCM6 for chromatin immunoprecipitation followed by sequencing (ChIP-seq) in the quest for identification of putative replication origins in the parasite. PfMCM6 DNA binding sites annotation revealed high enrichment at exon regions. This is contrary to higher eukaryotes that show an inclination of origin sites towards transcriptional start sites. ChIP-seq results were further validated by ChIP-qPCR results as well as nascent strand abundance assay at the selected PfMCM6 enriched sites that also showed preferential binding of PfORC1 suggesting potential of these sites as origin sites. Further, PfMCM6 ChIP-seq data showed a positive correlation with previously published histone H4K8Ac genome-wide binding sites but not with H3K9Ac sites suggesting epigenetic control of replication initiation sites in the parasites. Overall, our data show the genome-wide distribution of PfMCM6 binding sites with their potential as replication origins in this deadly human pathogen that not only broadens our knowledge of parasite DNA replication and its unique biology, it may help to find new avenues for intervention processes.

Dynamic Patterns of Mammalian Mitochondrial DNA Replication Uncovered Using SSiNGLe-5'ES.

The proper replication of mitochondrial DNA is key to the maintenance of this crucial organelle. Multiple studies aimed at understanding the mechanisms of replication of the mitochondrial genome have been conducted in the past several decades; however, while highly informative, they were conducted using relatively low-sensitivity techniques. Here, we established a high-throughput approach based on next-generation sequencing to identify replication start sites with nucleotide-level resolution and applied it to the genome of mitochondria from different human and mouse cell types. We found complex and highly reproducible patterns of mitochondrial initiation sites, both previously annotated and newly discovered in this work, that showed differences among different cell types and species. These results suggest that the patterns of the replication initiation sites are dynamic and might reflect, in some yet unknown ways, the complexities of mitochondrial and cellular physiology. Overall, this work suggests that much remains unknown about the details of mitochondrial DNA replication in different biological states, and the method established here opens up a new avenue in the study of the replication of mitochondrial and potentially other genomes.

Prevalent and dynamic binding of the cell cycle checkpoint kinase Rad53 to gene promoters.

Replication of the genome must be coordinated with gene transcription and cellular metabolism, especially following replication stress in the presence of limiting deoxyribonucleotides. The Saccharomyces cerevisiae Rad53 (CHEK2 in mammals) checkpoint kinase plays a major role in cellular responses to DNA replication stress. Cell cycle regulated, genome-wide binding of Rad53 to chromatin was examined. Under replication stress, the kinase bound to sites of active DNA replication initiation and fork progression, but unexpectedly to the promoters of about 20% of genes encoding proteins involved in multiple cellular functions. Rad53 promoter binding correlated with changes in expression of a subset of genes. Rad53 promoter binding to certain genes was influenced by sequence-specific transcription factors and less by checkpoint signaling. However, in checkpoint mutants, untimely activation of late-replicating origins reduces the transcription of nearby genes, with concomitant localization of Rad53 to their gene bodies. We suggest that the Rad53 checkpoint kinase coordinates genome-wide replication and transcription under replication stress conditions.

Temporal regulation of head-on transcription at replication initiation sites

Head-on (HO) collisions between the DNA replication machinery and RNA polymerase over R-loop forming sequences (RLFS) are genotoxic, leading to replication fork blockage and DNA breaks. Current models suggest that HO collisions are avoided through replication initiation site (RIS) positioning upstream of active genes, ensuring co-orientation of replication fork movement and genic transcription. However, this model does not account for pervasive transcription, or intragenic RIS. Moreover, pervasive transcription initiation and CG-rich DNA is a feature of RIS, suggesting that HO transcription units (HO TUs) capable of forming R-loops might occur. Through mining phased GRO-seq data, and developing an informatics strategy to stringently identify RIS, we demonstrate that HO TUs containing RLFS occur at RIS in MCF-7 cells, and are downregulated at the G1/S phase boundary. Our analysis reveals a novel spatiotemporal relationship between transcription and replication, and supports the idea that HO collisions are avoided through transcriptional regulatory mechanisms.

Modeling the Dynamics of Eukaryotic DNA Synthesis in Remembrance of Tunde Ogunnaike.

Two things Tunde loved were dynamics and probability. The work described herein combined them both, which explains why Tunde invariably asked me each time we talked how this work was proceeding. However, as I've come to appreciate as reminiscent of a surprisingly large amount of work in almost any researcher's career, I did not complete a peer-reviewed article on the matter while he could see it. We were broadly motivated by analysis of data for total DNA content in single cells, across thousands of cells. From such data one can estimate the proportions of cells in different phases of the cell cycle by fitting a mixture model for subpopulations of G0/G1 phase cells (1 relative copy of the genome), S phase cells (between 1 and 2 relative copies of the genome), and G2/M phase cells (2 relative copies of the genome). Given an asynchronously cycling population, Gaussian models are reasonable for the G0/G1 and G2/M subpopulations, but an appropriate functional form for the S-phase subpopulation was unclear. Since the probability of observing an S-phase cell is intimately related to the dynamics of DNA replication, we worked to derive a model for DNA replication dynamics from first principles, resulting in a closed-form, analytic expression for the dynamics of DNA synthesis. While quite arguably a somewhat superfluous effort, there is a certain satisfaction and academic beauty to modeling systems from a first-principles approach, and it can sometimes lead to unexpected scientific insights. Yet, while mathematically elegant, there was a fundamental issue with a key assumption that the so-called inter-origin distance distribution (distances between DNA replication initiation sites) was time-invariant. First, I present the model as developed previously. Then, to address the time-invariant inter-origin distance distribution issue, I provide a treatment of time-varying inter-origin distance distributions that, while mathematically simple, provides (i) mechanistic predictions for how all the DNA in a fertilized frog egg can be replicated "on time" despite some inter-origin distances initially exceeding the corresponding amount of allowable time and (ii) evidence that, based only on data from DNA content versus time and average inter-origin distances, somatic cell DNA is parsed into distinct regions whose replication is temporally separate.

Genome scale metabolic model combined with single molecule real-time sequencing to analyze Actinomycete chromosomal heterogeneity

The rapid development of biotechnology has provided new perspectives to observe and helped to gradually understand the significance of genetic instability in Actinobacteria. High frequency deletions of extremities and abnormal methylation of chromosomes suggest there might be relevant between the two phenomena. With this suspicion, we used single molecule real-time (SMRT) sequencing to map the genome-level methylation of one branch of actinomycetes, Saccharopolyspora erythraea, which have ring-shaped chromosomes. S. erythraea used for analysis in this study shares the same highly unstable phenotypic traits, as evidenced by diverse spore morphology and fluctuating erythromycin production. Multiple amplification of genomic islands closes to the replication initiation site and 6-methyladenine (m6A) deletion in genomic islands suggest that the interaction between the restriction modification (R-M) system and transposable elements provides an explanation for the division of labor by genomic heterogeneity in actinomycetes.

Escherichia coli cell factories with altered chromosomal replication scenarios exhibit accelerated growth and rapid biomass production

BackgroundGenerally, bacteria have a circular genome with a single replication origin for each replicon, whereas archaea and eukaryotes can have multiple replication origins in a single chromosome. In Escherichia coli, bidirectional DNA replication is initiated at the origin of replication (oriC) and arrested by the 10 termination sites (terA–J).ResultsWe constructed E. coli derivatives with additional or ectopic replication origins, which demonstrate the relationship between DNA replication and cell physiology. The cultures of E. coli derivatives with multiple replication origins contained an increased fraction of replicating chromosomes and the cells varied in size. Without the original oriC, E. coli derivatives with double ectopic replication origins manifested impaired growth irrespective of growth conditions and enhanced cell size, and exhibited excessive and asynchronous replication initiation. The generation time of an E. coli strain with three replication origins decreased in a minimal medium supplemented with glucose as the sole carbon source. As well as cell growth, the introduction of additional replication origins promoted increased biomass production.ConclusionsBalanced cell growth and physiological stability of E. coli under rapid growth condition are affected by changes in the position and number of replication origins. Additionally, we show that, for the first time to our knowledge, the introduction of replication initiation sites to the chromosome promotes cell growth and increases protein production.

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.

Comprehensive analysis of DNA replication timing across 184 cell lines suggests a role for MCM10 in replication timing regulation.

Cellular proliferation depends on the accurate and timely replication of the genome. Several genetic diseases are caused by mutations in key DNA replication genes; however, it remains unclear whether these genes influence the normal program of DNA replication timing. Similarly, the factors that regulate DNA replication dynamics are poorly understood. To systematically identify trans-acting modulators of replication timing, we profiled replication in 184 cell lines from three cell types, encompassing 60 different gene knockouts or genetic diseases. Through a rigorous approach that considers the background variability of replication timing, we concluded that most samples displayed normal replication timing. However, mutations in two genes showed consistently abnormal replication timing. The first gene was RIF1, a known modulator of replication timing. The second was MCM10, a highly conserved member of the pre-replication complex. Cells from a single patient carrying MCM10 mutations demonstrated replication timing variability comprising 46% of the genome and at different locations than RIF1 knockouts. Replication timing alterations in the mutated MCM10 cells were predominantly comprised of replication delays and initiation site gains and losses. Taken together, this study demonstrates the remarkable robustness of the human replication timing program and reveals MCM10 as a novel candidate modulator of DNA replication timing.

In vitro amplification of whole large plasmids via transposon-mediated oriC insertion.

DNA amplification is a fundamental technique in molecular biology. The replication cycle reactionis a new method for amplification of large circular DNA having oriC sequences, which is a replication initiation site of the Escherichia coli chromosome. We here developed areplication cycle reaction-based method useful for amplification of various circular DNAs lacking oriC, even in the absence of any sequence information, via transposon-mediated oriC insertion to the circular DNA template. A 15-kb non-oriC plasmid was amplified from a very small amount of starting DNA (50fg, 1fM). The method was also applicable toGC-rich plasmid (69%) or large F-plasmid (230kb). This method thus provides a powerful tool to amplify various environmental circular DNAs.

Increased expression of Polδ does not alter the canonical replication program in vivo.

Background: In vitro experiments utilising the reconstituted Saccharomyces cerevisiae eukaryotic replisome indicated that the efficiency of the leading strand replication is impaired by a moderate increase in Polδ concentration. It was hypothesised that the slower rate of the leading strand synthesis characteristic for reactions containing two-fold and four-fold increased concentration of Polδ represented a consequence of a relatively rare event, during which Polδ stochastically outcompeted Polε and, in an inefficient manner, temporarily facilitated extension of the leading strand. Inspired by this observation, we aimed to determine whether similarly increased Polδ levels influence replication dynamics in vivo using the fission yeast Schizosaccharomyces pombe as a model system. Methods: To generate S. pombe strains over-expressing Polδ, we utilised Cre-Lox mediated cassette exchange and integrated one or three extra genomic copies of all four Polδ genes. To estimate expression of respective Polδ genes in Polδ-overexpressing mutants, we measured relative transcript levels of cdc1 + , cdc6 + (or cdc6 L591G ), cdc27 + and cdm1 + by reverse transcription followed by quantitative PCR (RT-qPCR). To assess the impact of Polδ over-expression on cell physiology and replication dynamics, we used standard cell biology techniques and polymerase usage sequencing. Results: We provide an evidence that two-fold and four-fold over-production of Polδ does not significantly alter growth rate, cellular morphology and S-phase duration. Polymerase usage sequencing analysis further indicates that increased Polδ expression does not change activities of Polδ, Polε and Polα at replication initiation sites and across replication termination zones. Additionally, we show that mutants over-expressing Polδ preserve WT-like distribution of replication origin efficiencies. Conclusions: Our experiments do not disprove the existence of opportunistic polymerase switches; however, the data indicate that, if stochastic replacement of Polε for Polδ does occur i n vivo, it represents a rare phenomenon that does not significantly influence canonical replication program.

Increased expression of Polδ does not alter the canonical replication program in vivo

Background: In vitro experiments utilising the reconstituted Saccharomyces cerevisiae eukaryotic replisome indicated that the efficiency of the leading strand replication is impaired by a moderate increase in Polδ concentration. It was hypothesised that the slower rate of the leading strand synthesis characteristic for reactions containing two-fold and four-fold increased concentration of Polδ represented a consequence of a relatively rare event, during which Polδ stochastically outcompeted Polε and, in an inefficient manner, temporarily facilitated extension of the leading strand. Inspired by this observation, we aimed to determine whether similarly increased Polδ levels influence replication dynamics in vivo using the fission yeast Schizosaccharomyces pombe as a model system. Methods: To generate S. pombe strains over-expressing Polδ, we utilised Cre-Lox mediated cassette exchange and integrated one or three extra genomic copies of all four Polδ genes. To estimate expression of respective Polδ genes in Polδ-overexpressing mutants, we measured relative transcript levels of cdc1 +, cdc6 + (or cdc6 L591G), cdc27 + and cdm1 + by reverse transcription followed by quantitative PCR (RT-qPCR). To assess the impact of Polδ over-expression on cell physiology and replication dynamics, we used standard cell biology techniques and polymerase usage sequencing. Results: We provide an evidence that two-fold and four-fold over-production of Polδ does not significantly alter growth rate, cellular morphology and S-phase duration. Polymerase usage sequencing analysis further indicates that increased Polδ expression does not change activities of Polδ, Polε and Polα at replication initiation sites and across replication termination zones. Additionally, we show that mutants over-expressing Polδ preserve WT-like distribution of replication origin efficiencies. Conclusions: Our experiments do not disprove the existence of opportunistic polymerase switches; however, the data indicate that, if stochastic replacement of Polε for Polδ does occur i n vivo, it represents a rare phenomenon that does not significantly influence canonical replication program.

Cation competition and recruitment around the c-kit1 G-quadruplex using polarizable simulations

Nucleic acid-ion interactions are fundamentally important to the physical, energetic, and conformational properties of DNA and RNA. These interactions help fold and stabilize highly ordered secondary and tertiary structures, such as G-quadruplexes (GQs), which are functionally relevant in telomeres, replication initiation sites, and promoter sequences. The c-kit proto-oncogene encodes for a receptor tyrosine kinase and is linked to gastrointestinal stromal tumors, mast cell disease, and leukemia. This gene contains three unique GQ-forming sequences that have proposed antagonistic effects on gene expression. The dominant GQ, denoted c-kit1, has been shown to decrease expression of c-kit transcripts, making the c-kit1 GQ a promising drug target. Toward disease intervention, more information is needed regarding its conformational dynamics and ion binding properties. Therefore, we performed molecular dynamics simulations of the c-kit1 GQ with K+, Na+, Li+, and mixed salt solutions using the Drude-2017 polarizable force field. We evaluated GQ structure, ion sampling, core energetics, ion dehydration and binding, and ion competition and found that each analysis supported the known GQ-ion specificity trend (K+ > Na+ > Li+). We also found that K+ ions coordinate in the tetrad core antiprismatically, whereas Na+ and Li+ align coplanar to guanine tetrads, partially because of their attraction to surrounding water. Further, we showed that K+ occupancy is higher around the c-kit1 GQ and its nucleobases than Na+ and Li+, which tend to interact with backbone and sugar moieties. Finally, we showed that K+ binding to the c-kit1 GQ is faster and more frequent than Na+ and Li+. Such descriptions of GQ-ion dynamics suggest the rate of dehydration as the dominant factor for preference of K+ by DNA GQs and provide insight into noncanonical nucleic acids for which little experimental data exist.

Increased expression of Polδ does not alter the canonical replication program in vivo.

Background: In vitro experiments utilising the reconstituted Saccharomyces cerevisiae eukaryotic replisome indicated that the efficiency of the leading strand replication is impaired by a moderate increase in Polδ concentration. It was hypothesised that the slower rate of the leading strand synthesis characteristic for reactions containing two-fold and four-fold increased concentration of Polδ represented a consequence of a relatively rare event, during which Polδ stochastically outcompeted Polε and, in an inefficient manner, temporarily facilitated extension of the leading strand. Inspired by this observation, we aimed to determine whether similarly increased Polδ levels influence replication dynamics in vivo using the fission yeast Schizosaccharomyces pombe as a model system. Methods: To generate S. pombe strains over-expressing Polδ, we utilised Cre-Lox mediated cassette exchange and integrated one or three extra genomic copies of all four Polδ genes. To estimate expression of respective Polδ genes in Polδ-overexpressing mutants, we measured relative transcript levels of cdc1 + , cdc6 + (or cdc6 L591G ), cdc27 + and cdm1 + by reverse transcription followed by quantitative PCR (RT-qPCR). To assess the impact of Polδ over-expression on cell physiology and replication dynamics, we used standard cell biology techniques and polymerase usage sequencing. Results: We provide an evidence that two-fold and four-fold over-production of Polδ does not significantly alter growth rate, cellular morphology and S-phase duration. Polymerase usage sequencing analysis further indicates that increased Polδ expression does not change activities of Polδ, Polε and Polα at replication initiation sites and across replication termination zones. Additionally, we show that mutants over-expressing Polδ preserve WT-like distribution of replication origin efficiencies. Conclusions: Our experiments do not disprove the existence of opportunistic polymerase switches; however, the data indicate that, if stochastic replacement of Polε for Polδ does occur i n vivo, it represents a rare phenomenon that does not significantly influence canonical replication program.

INO80C Remodeler Maintains Genomic Stability by Preventing Promiscuous Transcription at Replication Origins

SUMMARYThe proper coordination of transcription with DNA replication and repair is central for genomic stability. We investigate how the INO80C chromatin remodeling enzyme might coordinate these genomic processes. We find that INO80C co-localizes with the origin recognition complex (ORC) at yeast replication origins and is bound to replication initiation sites in mouse embryonic stem cells (mESCs). In yeast· INO80C recruitment requires origin sequences but does not require ORC· suggesting that recruitment is independent of pre-replication complex assembly. In both yeast and ESCs· INO80C co-localizes at origins with Mot1 and NC2 transcription factors· and genetic studies suggest that they function together to promote genome stability. Interestingly· nascent transcript sequencing demonstrates that INO80C and Mot1 prevent pervasive transcription through origin sequences· and absence of these factors leads to formation of new DNA double-strand breaks. We propose that INO80C and Mot1/NC2 function through distinct pathways to limit origin transcription· maintaining genomic stability.

Bioinformatical dissection of fission yeast DNA replication origins.

Replication origins in eukaryotes form a base for assembly of the pre-replication complex (pre-RC), thereby serving as an initiation site of DNA replication. Characteristics of replication origin vary among species. In fission yeast Schizosaccharomyces pombe, DNA of high AT content is a distinct feature of replication origins; however, it remains to be understood what the general molecular architecture of fission yeast origin is. Here, we performed ChIP-seq mapping of Orc4 and Mcm2, two representative components of the pre-RC, and described the characteristics of their binding sites. The analysis revealed that fission yeast efficient origins are associated with two similar but independent features: a ≥15 bp-long motif with stretches of As and an AT-rich region of a few hundred bp. The A-rich motif was correlated with chromosomal binding of Orc, a DNA-binding component in the pre-RC, whereas the AT-rich region was associated with efficient binding of the DNA replicative helicase Mcm. These two features, in combination with the third feature, a transcription-poor region of approximately 1 kb, enabled to distinguish efficient replication origins from the rest of chromosome arms with high accuracy. This study, hence, provides a model that describes how multiple functional elements specify DNA replication origins in fission yeast genome.

Conditional knockout of RAD51-related genes in Leishmania major reveals a critical role for homologous recombination during genome replication.

Homologous recombination (HR) has an intimate relationship with genome replication, both during repair of DNA lesions that might prevent DNA synthesis and in tackling stalls to the replication fork. Recent studies led us to ask if HR might have a more central role in replicating the genome of Leishmania, a eukaryotic parasite. Conflicting evidence has emerged regarding whether or not HR genes are essential, and genome-wide mapping has provided evidence for an unorthodox organisation of DNA replication initiation sites, termed origins. To answer this question, we have employed a combined CRISPR/Cas9 and DiCre approach to rapidly generate and assess the effect of conditional ablation of RAD51 and three RAD51-related proteins in Leishmania major. Using this approach, we demonstrate that loss of any of these HR factors is not immediately lethal but in each case growth slows with time and leads to DNA damage and accumulation of cells with aberrant DNA content. Despite these similarities, we show that only loss of RAD51 or RAD51-3 impairs DNA synthesis and causes elevated levels of genome-wide mutation. Furthermore, we show that these two HR factors act in distinct ways, since ablation of RAD51, but not RAD51-3, has a profound effect on DNA replication, causing loss of initiation at the major origins and increased DNA synthesis at subtelomeres. Our work clarifies questions regarding the importance of HR to survival of Leishmania and reveals an unanticipated, central role for RAD51 in the programme of genome replication in a microbial eukaryote.

Replication origin location might contribute to genetic variability in Trypanosoma cruzi

BackgroundDNA replication in trypanosomatids operates in a uniquely challenging environment, since most of their genomes are constitutively transcribed. Trypanosoma cruzi, the etiological agent of Chagas disease, presents high variability in both chromosomes size and copy number among strains, though the underlying mechanisms are unknown.ResultsHere we have mapped sites of DNA replication initiation across the T. cruzi genome using Marker Frequency Analysis, which has previously only been deployed in two related trypanosomatids. The putative origins identified in T. cruzi show a notable enrichment of GC content, a preferential position at subtelomeric regions, coinciding with genes transcribed towards the telomeres, and a pronounced enrichment within coding DNA sequences, most notably in genes from the Dispersed Gene Family 1 (DGF-1).ConclusionsThese findings suggest a scenario where collisions between DNA replication and transcription are frequent, leading to increased genetic variability, as seen by the increase SNP levels at chromosome subtelomeres and in DGF-1 genes containing putative origins.