Replication of [AT/TA]25 Microsatellite Sequences by Human DNA Polymerase δ Holoenzymes Is Dependent on dNTP and RPA Levels.

Fragile sites are unstable genomic regions that are prone to breakage during stressed DNA replication. Several common fragile sites (CFS) contain A+T-rich regions including perfect [AT/TA] microsatellite repeats that may collapse into hairpins when in single-stranded DNA (ssDNA) form and coincide with chromosomal hotspots for breakage and rearrangements. While many factors contribute to CFS instability, evidence exists for replication stalling within [AT/TA] microsatellite repeats. Currently, it is unknown how stress causes replication stalling within [AT/TA] microsatellite repeats. To investigate this, we utilized FRET to characterize the structures of [AT/TA]25 sequences and also reconstituted lagging strand replication to characterize the progression of pol δ holoenzymes through A+T-rich sequences. The results indicate that [AT/TA]25 sequences adopt hairpins that are unwound by the major ssDNA-binding complex, RPA, and the progression of pol δ holoenzymes through A+T-rich sequences saturated with RPA is dependent on the template sequence and dNTP concentration. Importantly, the effects of RPA on the replication of [AT/TA]25 sequences are dependent on dNTP concentration, whereas the effects of RPA on the replication of A+T-rich, nonstructure-forming sequences are independent of dNTP concentration. Collectively, these results reveal complexities in lagging strand replication and provide novel insights into how [AT/TA] microsatellite repeats contribute to genome instability.

Replication of [AT/TA]25 microsatellite sequences by human DNA polymerase δ holoenzymes is dependent on dNTP and RPA levels.

Difficult-to-Replicate Sequences (DiToRS) are natural impediments in the human genome that inhibit DNA replication under endogenous replication. Some of the most widely-studied DiToRS are A+T-rich, high "flexibility regions," including long stretches of perfect [AT/TA] microsatellite repeats that have the potential to collapse into hairpin structures when in single-stranded DNA (ssDNA) form and are sites of recurrent structural variation and double-stranded DNA (dsDNA) breaks. Currently, it is unclear how these flexibility regions impact DNA replication, greatly limiting our fundamental understanding of human genome stability. To investigate replication through flexibility regions, we utilized FRET to characterize the effects of the major ssDNA-binding complex, RPA, on the structure of perfect [AT/TA]25 microsatellite repeats and also re-constituted human lagging strand replication to quantitatively characterize initial encounters of pol δ holoenzymes with A+T-rich DNA template sequences. The results indicate that [AT/TA]25 sequences adopt hairpin structures that are unwound by RPA and pol δ holoenzymes support dNTP incorporation through the [AT/TA]25 sequences as well as an A+T-rich, non-structure forming sequence. Furthermore, the extent of dNTP incorporation is dependent on the sequence of the DNA template and the concentration of dNTPs. Importantly, the effects of RPA on the replication of [AT/TA]25 sequences are dependent on the concentration of dNTPs, whereas the effects of RPA on the replication of an A+T-rich, non-structure forming sequence are independent of dNTP concentration. Collectively, these results reveal complexities in lagging strand replication and provide novel insights into how flexibility regions contribute to genome instability.

Accurate sequencing of DNA motifs able to form alternative (non-B) structures.

Approximately 13% of the human genome at certain motifs have the potential to form noncanonical (non-B) DNA structures (e.g., G-quadruplexes, cruciforms, and Z-DNA), which regulate many cellular processes but also affect the activity of polymerases and helicases. Because sequencing technologies use these enzymes, they might possess increased errors at non-B structures. To evaluate this, we analyzed error rates, read depth, and base quality of Illumina, Pacific Biosciences (PacBio) HiFi, and Oxford Nanopore Technologies (ONT) sequencing at non-B motifs. All technologies showed altered sequencing success for most non-B motif types, although this could be owing to several factors, including structure formation, biased GC content, and the presence of homopolymers. Single-nucleotide mismatch errors had low biases in HiFi and ONT for all non-B motif types but were increased for G-quadruplexes and Z-DNA in all three technologies. Deletion errors were increased for all non-B types but Z-DNA in Illumina and HiFi, as well as only for G-quadruplexes in ONT. Insertion errors for non-B motifs were highly, moderately, and slightly elevated in Illumina, HiFi, and ONT, respectively. Additionally, we developed a probabilistic approach to determine the number of false positives at non-B motifs depending on sample size and variant frequency, and applied it to publicly available data sets (1000 Genomes, Simons Genome Diversity Project, and gnomAD). We conclude that elevated sequencing errors at non-B DNA motifs should be considered in low-read-depth studies (single-cell, ancient DNA, and pooled-sample population sequencing) and in scoring rare variants. Combining technologies should maximize sequencing accuracy in future studies of non-B DNA.

Special Issue "DNA Replication/Repair, and the DNA Damage Response in Human Disease".

Mutations of numerous genes involved in DNA replication, DNA repair, and DNA damage response (DDR) pathways lead to a variety of human diseases, including aging and cancer [...].

Nontraditional Roles of DNA Polymerase Eta Support Genome Duplication and Stability.

DNA polymerase eta (Pol η) is a Y-family polymerase and the product of the POLH gene. Autosomal recessive inheritance of POLH mutations is the cause of the xeroderma pigmentosum variant, a cancer predisposition syndrome. This review summarizes mounting evidence for expanded Pol η cellular functions in addition to DNA lesion bypass that are critical for maintaining genome stability. In vitro, Pol η displays efficient DNA synthesis through difficult-to-replicate sequences, catalyzes D-loop extensions, and utilizes RNA-DNA hybrid templates. Human Pol η is constitutively present at the replication fork. In response to replication stress, Pol η is upregulated at the transcriptional and protein levels, and post-translational modifications regulate its localization to chromatin. Numerous studies show that Pol η is required for efficient common fragile site replication and stability. Additionally, Pol η can be recruited to stalled replication forks through protein-protein interactions, suggesting a broader role in replication fork recovery. During somatic hypermutations, Pol η is recruited by mismatch repair proteins and is essential for VH gene A:T basepair mutagenesis. Within the global context of repeat-dense genomes, the recruitment of Pol η to perform specialized functions during replication could promote genome stability by interrupting pure repeat arrays with base substitutions. Alternatively, not engaging Pol η in genome duplication is costly, as the absence of Pol η leads to incomplete replication and increased chromosomal instability.

Variation in G-quadruplex sequence and topology differentially impacts human DNA polymerase fidelity.

G-quadruplexes (G4s), a type of non-B DNA, play important roles in a wide range of molecular processes, including replication, transcription, and translation. Genome integrity relies on efficient and accurate DNA synthesis, and is compromised by various stressors, to which non-B DNA structures such as G4s can be particularly vulnerable. However, the impact of G4 structures on DNA polymerase fidelity is largely unknown. Using an in vitro forward mutation assay, we investigated the fidelity of human DNA polymerases delta (δ4, four-subunit), eta (η), and kappa (κ) during synthesis of G4 motifs representing those in the human genome. The motifs differ in sequence, topology, and stability, features that may affect DNA polymerase errors. Polymerase error rate hierarchy (δ4 < κ < Î·) is largely maintained during G4 synthesis. Importantly, we observed unique polymerase error signatures during synthesis of VEGF G4 motifs, stable G4s which form parallel topologies. These statistically significant errors occurred within, immediately flanking, and encompassing the G4 motif. For pol δ4, the errors were deletions, insertions and complex errors within the G4 or encompassing the G4 motif and surrounding sequence. For pol η, the errors occurred in 3' sequences flanking the G4 motif. For pol κ, the errors were frameshift mutations within G-tracts of the G4. Because these error signatures were not observed during synthesis of an antiparallel G4 and, to a lesser extent, a hybrid G4, we suggest that G4 topology and/or stability could influence polymerase fidelity. Using in silico analyses, we show that most polymerase errors are predicted to have minimal effects on predicted G4 stability. Our results provide a unique view of G4s not previously elucidated, showing that G4 motif heterogeneity differentially influences polymerase fidelity within the motif and flanking sequences. Thus, our study advances the understanding of how DNA polymerase errors contribute to G4 mutagenesis.

Impact of G-Quadruplexes and Chronic Inflammation on Genome Instability: Additive Effects during Carcinogenesis.

Genome instability is an enabling characteristic of cancer, essential for cancer cell evolution. Hotspots of genome instability, from small-scale point mutations to large-scale structural variants, are associated with sequences that potentially form non-B DNA structures. G-quadruplex (G4) forming motifs are enriched at structural variant endpoints in cancer genomes. Chronic inflammation is a physiological state underlying cancer development, and oxidative DNA damage is commonly invoked to explain how inflammation promotes genome instability. We summarize where G4s and oxidative stress overlap, with a focus on DNA replication. Guanine has low ionization potential, making G4s vulnerable to oxidative damage. Impacts to G4 structure are dependent upon lesion type, location, and G4 conformation. Occasionally, G4s pose a challenge to replicative DNA polymerases, requiring specialized DNA polymerases to maintain genome stability. Therefore, chronic inflammation creates a dual challenge for DNA polymerases to maintain genome stability: faithful G4 synthesis and bypassing unrepaired oxidative lesions. Inflammation is also accompanied by global transcriptome changes that may impact mutagenesis. Several studies suggest a regulatory role for G4s within cancer- and inflammatory-related gene promoters. We discuss the extent to which inflammation could influence gene regulation by G4s, thereby impacting genome instability, and highlight key areas for new investigation.

Selection and thermostability suggest G-quadruplexes are novel functional elements of the human genome.

Approximately 1% of the human genome has the ability to fold into G-quadruplexes (G4s)-noncanonical strand-specific DNA structures forming at G-rich motifs. G4s regulate several key cellular processes (e.g., transcription) and have been hypothesized to participate in others (e.g., firing of replication origins). Moreover, G4s differ in their thermostability, and this may affect their function. Yet, G4s may also hinder replication, transcription, and translation and may increase genome instability and mutation rates. Therefore, depending on their genomic location, thermostability, and functionality, G4 loci might evolve under different selective pressures, which has never been investigated. Here we conducted the first genome-wide analysis of G4 distribution, thermostability, and selection. We found an overrepresentation, high thermostability, and purifying selection for G4s within genic components in which they are expected to be functional-promoters, CpG islands, and 5' and 3' UTRs. A similar pattern was observed for G4s within replication origins, enhancers, eQTLs, and TAD boundary regions, strongly suggesting their functionality. In contrast, G4s on the nontranscribed strand of exons were underrepresented, were unstable, and evolved neutrally. In general, G4s on the nontranscribed strand of genic components had lower density and were less stable than those on the transcribed strand, suggesting that the former are avoided at the RNA level. Across the genome, purifying selection was stronger at stable G4s. Our results suggest that purifying selection preserves the sequences of functional G4s, whereas nonfunctional G4s are too costly to be tolerated in the genome. Thus, G4s are emerging as fundamental, functional genomic elements.

Overexpression of oncogenic H-Ras in hTERT-immortalized and SV40-transformed human cells targets replicative and specialized DNA polymerases for depletion.

DNA polymerases play essential functions in replication fork progression and genome maintenance. DNA lesions and drug-induced replication stress result in up-regulation and re-localization of specialized DNA polymerases η and κ. Although oncogene activation significantly alters DNA replication dynamics, causing replication stress and genome instability, little is known about DNA polymerase expression and regulation in response to oncogene activation. Here, we investigated the consequences of mutant H-RAS G12V overexpression on the regulation of DNA polymerases in h-TERT immortalized and SV40-transformed human cells. Focusing on DNA polymerases associated with the replication fork, we demonstrate that DNA polymerases are depleted in a temporal manner in response to H-RAS G12V overexpression. The polymerases targeted for depletion, as cells display markers of senescence, include the Pol α catalytic subunit (POLA1), Pol δ catalytic and p68 subunits (POLD1 and POLD3), Pol η, and Pol κ. Both transcriptional and post-transcriptional mechanisms mediate this response. Pol η (POLH) depletion is sufficient to induce a senescence-like growth arrest in human foreskin fibroblast BJ5a cells, and is associated with decreased Pol α expression. Using an SV-40 transformed cell model, we observed cell cycle checkpoint signaling differences in cells with H-RasG12V-induced polymerase depletion, as compared to Pol η-deficient cells. Our findings contribute to our understanding of cellular events following oncogene activation and cellular transformation.

Non-B DNA: a major contributor to small- and large-scale variation in nucleotide substitution frequencies across the genome.

Approximately 13% of the human genome can fold into non-canonical (non-B) DNA structures (e.g. G-quadruplexes, Z-DNA, etc.), which have been implicated in vital cellular processes. Non-B DNA also hinders replication, increasing errors and facilitating mutagenesis, yet its contribution to genome-wide variation in mutation rates remains unexplored. Here, we conducted a comprehensive analysis of nucleotide substitution frequencies at non-B DNA loci within noncoding, non-repetitive genome regions, their ±2 kb flanking regions, and 1-Megabase windows, using human-orangutan divergence and human single-nucleotide polymorphisms. Functional data analysis at single-base resolution demonstrated that substitution frequencies are usually elevated at non-B DNA, with patterns specific to each non-B DNA type. Mirror, direct and inverted repeats have higher substitution frequencies in spacers than in repeat arms, whereas G-quadruplexes, particularly stable ones, have higher substitution frequencies in loops than in stems. Several non-B DNA types also affect substitution frequencies in their flanking regions. Finally, non-B DNA explains more variation than any other predictor in multiple regression models for diversity or divergence at 1-Megabase scale. Thus, non-B DNA substantially contributes to variation in substitution frequencies at small and large scales. Our results highlight the role of non-B DNA in germline mutagenesis with implications to evolution and genetic diseases.

Sequence and Nuclease Requirements for Breakage and Healing of a Structure-Forming (AT)n Sequence within Fragile Site FRA16D.

Common fragile sites (CFSs) are genomic regions that display gaps and breaks in human metaphase chromosomes under replication stress and are often deleted in cancer cells. We studied an ∼300-bp subregion (Flex1) of human CFS FRA16D in yeast and found that it recapitulates characteristics of CFS fragility in human cells. Flex1 fragility is dependent on the ability of a variable-length AT repeat to form a cruciform structure that stalls replication. Fragility at Flex1 is initiated by structure-specific endonuclease Mus81-Mms4 acting together with the Slx1-4/Rad1-10 complex, whereas Yen1 protects Flex1 against breakage. Sae2 is required for healing of Flex1 after breakage. Our study shows that breakage within a CFS can be initiated by nuclease cleavage at forks stalled at DNA structures. Furthermore, our results suggest that CFSs are not just prone to breakage but also are impaired in their ability to heal, and this deleterious combination accounts for their fragility.

Long-read sequencing technology indicates genome-wide effects of non-B DNA on polymerization speed and error rate.

DNA conformation may deviate from the classical B-form in ∼13% of the human genome. Non-B DNA regulates many cellular processes; however, its effects on DNA polymerization speed and accuracy have not been investigated genome-wide. Such an inquiry is critical for understanding neurological diseases and cancer genome instability. Here, we present the first simultaneous examination of DNA polymerization kinetics and errors in the human genome sequenced with Single-Molecule Real-Time (SMRT) technology. We show that polymerization speed differs between non-B and B-DNA: It decelerates at G-quadruplexes and fluctuates periodically at disease-causing tandem repeats. Analyzing polymerization kinetics profiles, we predict and validate experimentally non-B DNA formation for a novel motif. We demonstrate that several non-B motifs affect sequencing errors (e.g., G-quadruplexes increase error rates), and that sequencing errors are positively associated with polymerase slowdown. Finally, we show that highly divergent G4 motifs have pronounced polymerization slowdown and high sequencing error rates, suggesting similar mechanisms for sequencing errors and germline mutations.

Genome-wide Identification of Structure-Forming Repeats as Principal Sites of Fork Collapse upon ATR Inhibition.

DNA polymerase stalling activates the ATR checkpoint kinase, which in turn suppresses fork collapse and breakage. Herein, we describe use of ATR inhibition (ATRi) as a means to identify genomic sites of problematic DNA replication in murine and human cells. Over 500 high-resolution ATR-dependent sites were ascertained using two distinct methods: replication protein A (RPA)-chromatin immunoprecipitation (ChIP) and breaks identified by TdT labeling (BrITL). The genomic feature most strongly associated with ATR dependence was repetitive DNA that exhibited high structure-forming potential. Repeats most reliant on ATR for stability included structure-forming microsatellites, inverted retroelement repeats, and quasi-palindromic AT-rich repeats. Notably, these distinct categories of repeats differed in the structures they formed and their ability to stimulate RPA accumulation and breakage, implying that the causes and character of replication fork collapse under ATR inhibition can vary in a DNA-structure-specific manner. Collectively, these studies identify key sources of endogenous replication stress that rely on ATR for stability.

DNA Polymerase Eta Prevents Tumor Cell-Cycle Arrest and Cell Death during Recovery from Replication Stress.

: Neoplastic transformation and genome instability are enhanced by replication stress, conditions that slow or stall DNA replication forks. Consequently, cancer cells require multiple enzymes and checkpoint signaling pathways to mitigate replication stress for their viability and proliferation. Targeting proteins that enhance cancer cell survival during replication stress is a recent approach in clinical strategies, especially when targets produce synthetic lethality. DNA polymerase eta (Pol η) has many key functions in genome stability, particularly for translesion synthesis. Here we demonstrate that endogenous Pol η displays significant protein induction and forms intense foci throughout the nucleus in response to replication stress induced by drugs that do not directly form DNA adducts. During replication stress, Pol η-deficient cells displayed hyperactivation of the ATR replication checkpoint and arrested late in the cell cycle. During recovery from replication stress, Pol η-deficient cells continue to display aberrant phenotypes, including delayed cell-cycle progression, apoptosis, and cell survival. Depletion or inhibition of ATR was synthetically lethal with Pol η deficiency, particularly when tumor cells were treated with replication stress-inducing drugs. Together our data expand knowledge of the cellular environments that increase endogenous Pol η expression beyond DNA damaging agents and demonstrate that Pol η regulation is central to the replication stress response. Because Pol η is aberrantly expressed in several tumor types, our results are critical for developing more effective chemotherapy approaches and identify coinhibition of Pol η and ATR as a potential therapeutic strategy. SIGNIFICANCE: This study demonstrates that replication stress upregulates Pol η (POLH) in tumor cells and reveals a role for Pol η in tumor cell recovery following replication stress.

Detours to Replication: Functions of Specialized DNA Polymerases during Oncogene-induced Replication Stress.

Incomplete and low-fidelity genome duplication contribute to genomic instability and cancer development. Difficult-to-Replicate Sequences, or DiToRS, are natural impediments in the genome that require specialized DNA polymerases and repair pathways to complete and maintain faithful DNA synthesis. DiToRS include non B-DNA secondary structures formed by repetitive sequences, for example within chromosomal fragile sites and telomeres, which inhibit DNA replication under endogenous stress conditions. Oncogene activation alters DNA replication dynamics and creates oncogenic replication stress, resulting in persistent activation of the DNA damage and replication stress responses, cell cycle arrest, and cell death. The response to oncogenic replication stress is highly complex and must be tightly regulated to prevent mutations and tumorigenesis. In this review, we summarize types of known DiToRS and the experimental evidence supporting replication inhibition, with a focus on the specialized DNA polymerases utilized to cope with these obstacles. In addition, we discuss different causes of oncogenic replication stress and its impact on DiToRS stability. We highlight recent findings regarding the regulation of DNA polymerases during oncogenic replication stress and the implications for cancer development.

DNA polymerases eta and kappa exchange with the polymerase delta holoenzyme to complete common fragile site synthesis.

Common fragile sites (CFSs) are inherently unstable genomic loci that are recurrently altered in human tumor cells. Despite their instability, CFS are ubiquitous throughout the human genome and associated with large tumor suppressor genes or oncogenes. CFSs are enriched with repetitive DNA sequences, one feature postulated to explain why these loci are inherently difficult to replicate, and sensitive to replication stress. We have shown that specialized DNA polymerases (Pols) η and κ replicate CFS-derived sequences more efficiently than the replicative Pol δ. However, we lacked an understanding of how these enzymes cooperate to ensure efficient CFS replication. Here, we designed a model of lagging strand replication with RFC loaded PCNA that allows for maximal activity of the four-subunit human Pol δ holoenzyme, Pol η, and Pol κ in polymerase mixing assays. We discovered that Pol η and κ are both able to exchange with Pol δ stalled at repetitive CFS sequences, enhancing Normalized Replication Efficiency. We used this model to test the impact of PCNA mono-ubiquitination on polymerase exchange, and found no change in polymerase cooperativity in CFS replication compared with unmodified PCNA. Finally, we modeled replication stress in vitro using aphidicolin and found that Pol δ holoenzyme synthesis was significantly inhibited in a dose-dependent manner, preventing any replication past the CFS. Importantly, Pol η and κ were still proficient in rescuing this stalled Pol δ synthesis, which may explain, in part, the CFS instability phenotype of aphidicolin-treated Pol η and Pol κ-deficient cells. In total, our data support a model wherein Pol δ stalling at CFSs allows for free exchange with a specialized polymerase that is not driven by PCNA.

Maintenance of Genome Integrity: How Mammalian Cells Orchestrate Genome Duplication by Coordinating Replicative and Specialized DNA Polymerases.

Precise duplication of the human genome is challenging due to both its size and sequence complexity. DNA polymerase errors made during replication, repair or recombination are central to creating mutations that drive cancer and aging. Here, we address the regulation of human DNA polymerases, specifically how human cells orchestrate DNA polymerases in the face of stress to complete replication and maintain genome stability. DNA polymerases of the B-family are uniquely adept at accurate genome replication, but there are numerous situations in which one or more additional DNA polymerases are required to complete genome replication. Polymerases of the Y-family have been extensively studied in the bypass of DNA lesions; however, recent research has revealed that these polymerases play important roles in normal human physiology. Replication stress is widely cited as contributing to genome instability, and is caused by conditions leading to slowed or stalled DNA replication. Common Fragile Sites epitomize "difficult to replicate" genome regions that are particularly vulnerable to replication stress, and are associated with DNA breakage and structural variation. In this review, we summarize the roles of both the replicative and Y-family polymerases in human cells, and focus on how these activities are regulated during normal and perturbed genome replication.

Reverse Transcription Errors and RNA-DNA Differences at Short Tandem Repeats.

Transcript variation has important implications for organismal function in health and disease. Most transcriptome studies focus on assessing variation in gene expression levels and isoform representation. Variation at the level of transcript sequence is caused by RNA editing and transcription errors, and leads to nongenetically encoded transcript variants, or RNA-DNA differences (RDDs). Such variation has been understudied, in part because its detection is obscured by reverse transcription (RT) and sequencing errors. It has only been evaluated for intertranscript base substitution differences. Here, we investigated transcript sequence variation for short tandem repeats (STRs). We developed the first maximum-likelihood estimator (MLE) to infer RT error and RDD rates, taking next generation sequencing error rates into account. Using the MLE, we empirically evaluated RT error and RDD rates for STRs in a large-scale DNA and RNA replicated sequencing experiment conducted in a primate species. The RT error rates increased exponentially with STR length and were biased toward expansions. The RDD rates were approximately 1 order of magnitude lower than the RT error rates. The RT error rates estimated with the MLE from a primate data set were concordant with those estimated with an independent method, barcoded RNA sequencing, from a Caenorhabditis elegans data set. Our results have important implications for medical genomics, as STR allelic variation is associated with >40 diseases. STR nonallelic transcript variation can also contribute to disease phenotype. The MLE and empirical rates presented here can be used to evaluate the probability of disease-associated transcripts arising due to RDD.