Positive correlation between DNA polymerase alpha-primase pausing and mutagenesis within polypyrimidine/polypurine microsatellite sequences

Pathogenic CANVAS (AAGGG)n repeats stall DNA replication due to the formation of alternative DNA structures

CANVAS is a recently characterized repeat expansion disease, most commonly caused by homozygous expansions of an intronic (A2G3)n repeat in the RFC1 gene. There are a multitude of repeat motifs found in the human population at this locus, some of which are pathogenic and others benign. In this study, we conducted structure-functional analyses of the pathogenic (A2G3)n and nonpathogenic (A4G)n repeats. We found that the pathogenic, but not the nonpathogenic, repeat presents a potent, orientation-dependent impediment to DNA polymerization in vitro. The pattern of the polymerization blockage is consistent with triplex or quadruplex formation in the presence of magnesium or potassium ions, respectively. Chemical probing of both repeats in vitro reveals triplex H-DNA formation by only the pathogenic repeat. Consistently, bioinformatic analysis of S1-END-seq data from human cell lines shows preferential H-DNA formation genome-wide by (A2G3)n motifs over (A4G)n motifs. Finally, the pathogenic, but not the nonpathogenic, repeat stalls replication fork progression in yeast and human cells. We hypothesize that the CANVAS-causing (A2G3)n repeat represents a challenge to genome stability by folding into alternative DNA structures that stall DNA replication.

Pathogenic CANVAS (AAGGG)n repeats stall DNA replication due to the formation of alternative DNA structures

CANVAS is a recently characterized repeat expansion disease, most commonly caused by homozygous expansions of an intronic (A2G3)n repeat in the RFC1 gene. There are a multitude of repeat motifs found in the human population at this locus, some of which are pathogenic and others benign. In this study, we conducted structure-functional analyses of the main pathogenic (A2G3)n and the main nonpathogenic (A4G)n repeats. We found that the pathogenic, but not the nonpathogenic, repeat presents a potent, orientation-dependent impediment to DNA polymerization in vitro. The pattern of the polymerization blockage is consistent with triplex or quadruplex formation in the presence of magnesium or potassium ions, respectively. Chemical probing of both repeats in supercoiled DNA reveals triplex H-DNA formation by the pathogenic repeat. Consistently, bioinformatic analysis of the S1-END-seq data from human cell lines shows preferential H-DNA formation genome-wide by (A2G3)n motifs over (A4G)n motifs in vivo. Finally, the pathogenic, but not the non-pathogenic, repeat stalls replication fork progression in yeast and human cells. We hypothesize that CANVAS-causing (A2G3)n repeat represents a challenge to genome stability by folding into alternative DNA structures that stall DNA replication.

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.

Dynamic alternative DNA structures in biology and disease

Repetitive elements in the human genome, once considered 'junk DNA', are now known to adopt more than a dozen alternative (that is, non-B) DNA structures, such as self-annealed hairpins, left-handed Z-DNA, three-stranded triplexes (H-DNA) or four-stranded guanine quadruplex structures (G4 DNA). These dynamic conformations can act as functional genomic elements involved in DNA replication and transcription, chromatin organization and genome stability. In addition, recent studies have revealed a role for these alternative structures in triggering error-generating DNA repair processes, thereby actively enabling genome plasticity. As a driving force for genetic variation, non-B DNA structures thus contribute to both disease aetiology and evolution.

Translesion polymerase eta both facilitates DNA replication and promotes increased human genetic variation at common fragile sites

Common fragile sites (CFSs) are difficult-to-replicate genomic regions that form gaps and breaks on metaphase chromosomes under replication stress. They are hotspots for chromosomal instability in cancer. Repetitive sequences located at CFS loci are inefficiently copied by replicative DNA polymerase (Pol) delta. However, translesion synthesis Pol eta has been shown to efficiently polymerize CFS-associated repetitive sequences in vitro and facilitate CFS stability by a mechanism that is not fully understood. Here, by locus-specific, single-molecule replication analysis, we identified a crucial role for Pol eta (encoded by the gene POLH) in the in vivo replication of CFSs, even without exogenous stress. We find that Pol eta deficiency induces replication pausing, increases initiation events, and alters the direction of replication-fork progression at CFS-FRA16D in both lymphoblasts and fibroblasts. Furthermore, certain replication pause sites at CFS-FRA16D were associated with the presence of non-B DNA-forming motifs, implying that non-B DNA structures could increase replication hindrance in the absence of Pol eta. Further, in Pol eta-deficient fibroblasts, there was an increase in fork pausing at fibroblast-specific CFSs. Importantly, while not all pause sites were associated with non-B DNA structures, they were embedded within regions of increased genetic variation in the healthy human population, with mutational spectra consistent with Pol eta activity. From these findings, we propose that Pol eta replicating through CFSs may result in genetic variations found in the human population at these sites.

PARTICLE - The RNA podium for genomic silencers

Radiation exposure can evoke cellular stress responses. Emerging recognition that long non-coding RNAs (lncRNAs) act as regulators of gene expression has broadened the spectra of molecules controlling the genomic landscape upon alterations in environmental conditions. Knowledge of the mechanisms responding to low dose irradiation (LDR) exposure is very limited yet most likely involve subtle ancillary molecular pathways other than those protecting the cell from direct cellular damage. The discovery that transcription of the lncRNA PARTICLE (promoter of MAT2A- antisense radiation-induced circulating lncRNA; PARTICL) becomes dramatically instigated within a day after LDR exposure introduced a new gene regulator onto the biological landscape. PARTICLE affords an RNA binding platform for genomic silencers such as DNA methyltransferase 1 and histone tri-methyltransferases to reign in the expression of tumor suppressors such as its neighboring MAT2A in cis as well as WWOX in trans. In silico evidence offers scope to speculate that PARTICLE exploits the abundance of Hoogsten bonds that exist throughout mammalian genomes for triplex formation, presumably a vital feature within this RNA silencer. PARTICLE may provide a buffering riboswitch platform for S-adenosylmethionine. The correlation of PARTICLE triplex formation sites within tumor suppressor genes and their abundance throughout the genome at cancer-related hotspots offers an insight into potential avenues worth exploring in future therapeutic endeavors.

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.

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.

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.

Effect of dC → d(m5C) substitutions on the folding of intramolecular triplexes with mixed TAT and C+GC base triplets

Oligonucleotide-directed triple helix formation has been recognized as a potential tool for targeting genes with high specificity. Cystosine methylation in the 5' position is both ubiquitous and a stable regulatory modification, which could potentially stabilize triple helix formation. In this work, we have used a combination of calorimetric and spectroscopic techniques to study the intramolecular unfolding of four triplexes and two duplexes. We used the following triplex control sequence, named Control Tri, d(AGAGAC5TCTCTC5TCTCT), where C5 are loops of five cytosines. From this sequence, we studied three other sequences with dC → d(m5C) substitutions on the Hoogsteen strand (2MeH), Crick strand (2MeC) and both strands (4MeHC). Calorimetric studies determined that methylation does increase the thermal and enthalpic stability, leading to an overall favorable free energy, and that this increased stability is cumulative, i.e. methylation on both the Hoogsteen and Crick strands yields the largest favorable free energy. The differential uptake of protons, counterions and water was determined. It was found that methylation increases cytosine protonation by shifting the apparent pKa value to a higher pH; this increase in proton uptake coincides with a release of counterions during folding of the triplex, likely due to repulsion from the increased positive charge from the protonated cytosines. The immobilization of water was not affected for triplexes with methylated cytosines on their Hoogsteen or Crick strands, but was seen for the triplex where both strands are methylated. This may be due to the alignment in the major groove of the methyl groups on the cytosines with the methyl groups on the thymines which causes an increase in structural water along the spine of the triplex.

Effect of Loop Length and Sequence on the Stability of DNA Pyrimidine Triplexes with TAT Base Triplets

We report the thermodynamic contributions of loop length and loop sequence to the overall stability of DNA intramolecular pyrimidine triplexes. Two sets of triplexes were designed: in the first set, the C₅ loop closing the triplex stem was replaced with ^5'-CT_nC loops (n = 1-5), whereas in the second set, both the duplex and triplex loops were replaced with a ^5'-GCAA or ^5'-AACG tetraloop. For the triplexes with a ^5'-CT_nC loop, the triplex with five bases in the loop has the highest stability relative to the control. A loop length lower than five compromises the strength of the base-pair stacks without decreasing the thermal stability, leading to a decreased enthalpy, whereas an increase in the loop length leads to a decreased enthalpy and a higher entropic penalty. The incorporation of the GCAA loop yielded more stable triplexes, whereas the incorporation of AACG in the triplex loop yielded a less stable triplex due to an unfavorable enthalpy term. Thus, addition of the GCAA tetraloop can cause an increase in the thermodynamics of the triplex without affecting the sequence or melting behavior and may result in an additional layer of genetic regulation.

Functional Mechanisms of Microsatellite DNA in Eukaryotic Genomes

Microsatellite repeat DNA is best known for its length mutability, which is implicated in several neurological diseases and cancers, and often exploited as a genetic marker. Less well-known is the body of work exploring the widespread and surprisingly diverse functional roles of microsatellites. Recently, emerging evidence includes the finding that normal microsatellite polymorphism contributes substantially to the heritability of human gene expression on a genome-wide scale, calling attention to the task of elucidating the mechanisms involved. At present, these are underexplored, but several themes have emerged. I review evidence demonstrating roles for microsatellites in modulation of transcription factor binding, spacing between promoter elements, enhancers, cytosine methylation, alternative splicing, mRNA stability, selection of transcription start and termination sites, unusual structural conformations, nucleosome positioning and modification, higher order chromatin structure, noncoding RNA, and meiotic recombination hot spots.

Role of specialized DNA polymerases in the limitation of replicative stress and DNA damage transmission

Replication stress is a strong and early driving force for genomic instability and tumor development. Beside replicative DNA polymerases, an emerging group of specialized DNA polymerases is involved in the technical assistance of the replication machinery in order to prevent replicative stress and its deleterious consequences. During S-phase, altered progression of the replication fork by endogenous or exogenous impediments induces replicative stress, causing cells to reach mitosis with genomic regions not fully duplicated. Recently, specific mechanisms to resolve replication intermediates during mitosis with the aim of limiting DNA damage transmission to daughter cells have been identified. In this review, we detail the two major actions of specialized DNA polymerases that limit DNA damage transmission: the prevention of replicative stress by non-B DNA replication and the recovery of stalled replication forks.

On avoided words, absent words, and their application to biological sequence analysis

Background

The deviation of the observed frequency of a word w from its expected frequency in a given sequence x is used to determine whether or not the word is avoided. This concept is particularly useful in DNA linguistic analysis. The value of the deviation of w, denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{dev}(w)$$\end{document}dev(w), effectively characterises the extent of a word by its edge contrast in the context in which it occurs. A word w of length \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k>2$$\end{document}k>2 is a \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}ρ-avoided word in x if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textit{dev}(w) \le \rho $$\end{document}dev(w)≤ρ, for a given threshold \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho < 0$$\end{document}ρ<0. Notice that such a word may be completely absent from x. Hence, computing all such words naïvely can be a very time-consuming procedure, in particular for large k.

Results

In this article, we propose an \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(n)$$\end{document}O(n)-time and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(n)$$\end{document}O(n)-space algorithm to compute all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}ρ-avoided words of length k in a given sequence of length n over a fixed-sized alphabet. We also present a time-optimal \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {O}(\sigma n)$$\end{document}O(σn)-time algorithm to compute all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}ρ-avoided words (of any length) in a sequence of length n over an integer alphabet of size \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document}σ. In addition, we provide a tight asymptotic upper bound for the number of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\rho $$\end{document}ρ-avoided words over an integer alphabet and the expected length of the longest one. We make available an implementation of our algorithm. Experimental results, using both real and synthetic data, show the efficiency and applicability of our implementation in biological sequence analysis.

Conclusions

The systematic search for avoided words is particularly useful for biological sequence analysis. We present a linear-time and linear-space algorithm for the computation of avoided words of length k in a given sequence x. We suggest a modification to this algorithm so that it computes all avoided words of x, irrespective of their length, within the same time complexity. We also present combinatorial results with regards to avoided words and absent words.

Regulation of DNA Replication through Natural Impediments in the Eukaryotic Genome

All living organisms need to duplicate their genetic information while protecting it from unwanted mutations, which can lead to genetic disorders and cancer development. Inaccuracies during DNA replication are the major cause of genomic instability, as replication forks are prone to stalling and collapse, resulting in DNA damage. The presence of exogenous DNA damaging agents as well as endogenous difficult-to-replicate DNA regions containing DNA–protein complexes, repetitive DNA, secondary DNA structures, or transcribing RNA polymerases, increases the risk of genomic instability and thus threatens cell survival. Therefore, understanding the cellular mechanisms required to preserve the genetic information during S phase is of paramount importance. In this review, we will discuss our current understanding of how cells cope with these natural impediments in order to prevent DNA damage and genomic instability during DNA replication.

Effects of Replication and Transcription on DNA Structure-Related Genetic Instability

Many repetitive sequences in the human genome can adopt conformations that differ from the canonical B-DNA double helix (i.e., non-B DNA), and can impact important biological processes such as DNA replication, transcription, recombination, telomere maintenance, viral integration, transposome activation, DNA damage and repair. Thus, non-B DNA-forming sequences have been implicated in genetic instability and disease development. In this article, we discuss the interactions of non-B DNA with the replication and/or transcription machinery, particularly in disease states (e.g., tumors) that can lead to an abnormal cellular environment, and how such interactions may alter DNA replication and transcription, leading to potential conflicts at non-B DNA regions, and eventually result in genetic stability and human disease.

Divalent ions attenuate DNA synthesis by human DNA polymerase α by changing the structure of the template/primer or by perturbing the polymerase reaction

DNA polymerases (pols) are sophisticated protein machines operating in the replication, repair and recombination of genetic material in the complex environment of the cell. DNA pol reactions require at least two divalent metal ions for the phosphodiester bond formation. We explore two understudied roles of metals in pol transactions with emphasis on polα, a crucial enzyme in the initiation of DNA synthesis. We present evidence that the combination of many factors, including the structure of the template/primer, the identity of the metal, the metal turnover in the pol active site, and the influence of the concentration of nucleoside triphosphates, affect DNA pol synthesis. On the poly-dT70 template, the increase of Mg(2+) concentration within the range typically used for pol reactions led to the severe loss of the ability of pol to extend DNA primers and led to a decline in DNA product sizes when extending RNA primers, simulating the effect of "counting" of the number of nucleotides in nascent primers by polα. We suggest that a high Mg(2+) concentration promotes the dynamic formation of unconventional DNA structure(s), thus limiting the apparent processivity of the enzyme. Next, we found that Zn(2+) supported robust polα reactions when the concentration of nucleotides was above the concentration of ions; however, there was only one nucleotide incorporation by the Klenow fragment of DNA pol I. Zn(2+) drastically inhibited polα, but had no effect on Klenow, when Mg(2+) was also present. It is possible that Zn(2+) perturbs metal-mediated transactions in pol active site, for example affecting the step of pyrophosphate removal at the end of each pol cycle necessary for continuation of polymerization.

New Perspectives on DNA and RNA Triplexes As Effectors of Biological Activity

Since the first description of the canonical B-form DNA double helix, it has been suggested that alternative DNA, DNA–RNA, and RNA structures exist and act as functional genomic elements. Indeed, over the past few years it has become clear that, in addition to serving as a repository for genetic information, genomic DNA elicits biological responses by adopting conformations that differ from the canonical right-handed double helix, and by interacting with RNA molecules to form complex secondary structures. This review focuses on recent advances on three-stranded (triplex) nucleic acids, with an emphasis on DNA–RNA and RNA–RNA interactions. Emerging work reveals that triplex interactions between noncoding RNAs and duplex DNA serve as platforms for delivering site-specific epigenetic marks critical for the regulation of gene expression. Additionally, an increasing body of genetic and structural studies demonstrates that triplex RNA–RNA interactions are essential for performing catalytic and regulatory functions in cellular nucleoprotein complexes, including spliceosomes and telomerases, and for enabling protein recoding during programmed ribosomal frameshifting. Thus, evidence is mounting that DNA and RNA triplex interactions are implemented to perform a range of diverse biological activities in the cell, some of which will be discussed in this review.

Sequence differences at orthologous microsatellites inflate estimates of human-chimpanzee differentiation

Background

Microsatellites---contiguous arrays of 2–6 base-pair motifs---have formed the cornerstone of population-genetic studies for over two decades. Their genotype data typically takes the form of PCR fragment lengths obtained using locus-specific primer pairs to amplify the genomic region encompassing the microsatellite. Recently, we reported a dataset of 5,795 human and 84 chimpanzee individuals with genotypes at 246 human-derived autosomal microsatellites as a resource to facilitate interspecies comparisons. A major assumption underlying this dataset is that PCR amplicons at orthologous microsatellites are commensurable between species.

Results

We find this assumption to be frequently incorrect owing to discordance in microsatellite organization and variability, as well as nontrivial length imbalances caused by small species-specific indels in microsatellite flanking sequences. Converting PCR fragment lengths into the repeat numbers they represent at 138 microsatellites whose organization and variability was found to be highly similar in both species, we show that interspecies incommensurability among PCR amplicons can inflate F_ST and D_PS estimates by up to 10.6%. Separate investigations of determinants of microsatellite variability in humans and chimpanzees uncover similar patterns with mean and maximum numbers of repeats, as well as numbers and ranges of distinct alleles, all important factors in predicting heterozygosity. In contrast, across microsatellites, numbers of repeats were significantly smaller in chimpanzees than in humans, while numbers and ranges of distinct alleles were instead larger.

Conclusions

Our findings have fundamental implications for interspecies comparisons using microsatellites and offer new opportunities for more accurate comparisons of patterns of human and chimpanzee genetic variation in numerous areas of application.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-990) contains supplementary material, which is available to authorized users.

Microsatellite interruptions stabilize primate genomes and exist as population-specific single nucleotide polymorphisms within individual human genomes

Interruptions of microsatellite sequences impact genome evolution and can alter disease manifestation. However, human polymorphism levels at interrupted microsatellites (iMSs) are not known at a genome-wide scale, and the pathways for gaining interruptions are poorly understood. Using the 1000 Genomes Phase-1 variant call set, we interrogated mono-, di-, tri-, and tetranucleotide repeats up to 10 units in length. We detected ∼26,000–40,000 iMSs within each of four human population groups (African, European, East Asian, and American). We identified population-specific iMSs within exonic regions, and discovered that known disease-associated iMSs contain alleles present at differing frequencies among the populations. By analyzing longer microsatellites in primate genomes, we demonstrate that single interruptions result in a genome-wide average two- to six-fold reduction in microsatellite mutability, as compared with perfect microsatellites. Centrally located interruptions lowered mutability dramatically, by two to three orders of magnitude. Using a biochemical approach, we tested directly whether the mutability of a specific iMS is lower because of decreased DNA polymerase strand slippage errors. Modeling the adenomatous polyposis coli tumor suppressor gene sequence, we observed that a single base substitution interruption reduced strand slippage error rates five- to 50-fold, relative to a perfect repeat, during synthesis by DNA polymerases α, β, or η. Computationally, we demonstrate that iMSs arise primarily by base substitution mutations within individual human genomes. Our biochemical survey of human DNA polymerase α, β, δ, κ, and η error rates within certain microsatellites suggests that interruptions are created most frequently by low fidelity polymerases. Our combined computational and biochemical results demonstrate that iMSs are abundant in human genomes and are sources of population-specific genetic variation that may affect genome stability. The genome-wide identification of iMSs in human populations presented here has important implications for current models describing the impact of microsatellite polymorphisms on gene expression.

Microsatellites are short tandem repeat DNA sequences located throughout the human genome that display a high degree of inter-individual variation. This characteristic makes microsatellites an attractive tool for population genetics and forensics research. Some microsatellites affect gene expression, and mutations within such microsatellites can cause disease. Interruption mutations disrupt the perfect repeated array and are frequently associated with altered disease risk, but they have not been thoroughly studied in human genomes. We identified interrupted mono-, di-, tri- and tetranucleotide MSs (iMS) within individual genomes from African, European, Asian and American population groups. We show that many iMSs, including some within disease-associated genes, are unique to a single population group. By measuring the conservation of microsatellites between human and chimpanzee genomes, we demonstrate that interruptions decrease the probability of microsatellite mutations throughout the genome. We demonstrate that iMSs arise in the human genome by single base changes within the DNA, and provide biochemical data suggesting that these stabilizing changes may be created by error-prone DNA polymerases. Our genome-wide study supports the model in which iMSs act to stabilize individual genomes, and suggests that population-specific differences in microsatellite architecture may be an avenue by which genetic ancestry impacts individual disease risk.

Impact of alternative DNA structures on DNA damage, DNA repair, and genetic instability

Repetitive genomic sequences can adopt a number of alternative DNA structures that differ from the canonical B-form duplex (i.e. non-B DNA). These non-B DNA-forming sequences have been shown to have many important biological functions related to DNA metabolic processes; for example, they may have regulatory roles in DNA transcription and replication. In addition to these regulatory functions, non-B DNA can stimulate genetic instability in the presence or absence of DNA damage, via replication-dependent and/or replication-independent pathways. This review focuses on the interactions of non-B DNA conformations with DNA repair proteins and how these interactions impact genetic instability.

The C-terminal domain of the DNA polymerase catalytic subunit regulates the primase and polymerase activities of the human DNA polymerase α-primase complex

The initiation of DNA synthesis during replication of the human genome is accomplished primarily by the DNA polymerase α-primase complex, which makes the RNA-DNA primers accessible to processive DNA pols. The structural information needed to understand the mechanism of regulation of this complex biochemical reaction is incomplete. The presence of two enzymes in one complex poses the question of how these two enzymes cooperate during priming of DNA synthesis. Yeast two-hybrid and direct pulldown assays revealed that the N-terminal domain of the large subunit of primase (p58N) directly interacts with the C-terminal domain of the catalytic subunit of polα (p180C). We found that a complex of the C-terminal domain of the catalytic subunit of polα with the second subunit (p180C-p70) stimulated primase activity, whereas the whole catalytically active heterodimer of polα (p180ΔN-p70) inhibited RNA synthesis by primase. Conversely, the polα catalytic domain without the C-terminal part (p180ΔN-core) possessed a much higher propensity to extend the RNA primer than the two-subunit polα (p180ΔN-p70), suggesting that p180C and/or p70 are involved in the negative regulation of DNA pol activity. We conclude that the interaction between p180C, p70, and p58 regulates the proper primase and polymerase function. The composition of the template DNA is another important factor determining the activity of the complex. We have found that polα activity strongly depends on the sequence of the template and that homopyrimidine runs create a strong barrier for DNA synthesis by polα.

A call to arms: on refining Plasmodium vivax microsatellite marker panels for comparing global diversity

BACKGROUND

Microsatellite (MS) markers have become an important tool for studying the population diversity, evolutionary history and multiplicity of infection (MOI) of malaria parasite infections. MS are typically selected on the basis of being highly polymorphic. However, it is known that the polymorphic potential (mutability) of each marker can vary as much as two orders of magnitude, which radically changes how diversity is represented in the genome from one marker to the next. Over the past decade, approximately 240 Plasmodium vivax MS have been published, comprising nine major panels of markers. Inconsistent usage of each panel has resulted in a surfeit of descriptive genetic diversity data that are largely incomparable between populations. The objective of this study was to statistically evaluate the quality of individual MS markers in order to validate a refined panel of markers that will provide a balanced picture of P. vivax population diversity.

METHODS

All previously published data, including genetic diversity indices, MS parameters, and population parameters, were assembled from 18 different global studies into a flat file to facilitate statistical analysis and modelling using JMP® Genomics 6.0 (SAS Institute Inc, Cary, NC, USA). Statistical modeling was employed to down-select markers with extreme variation among the mean number of alleles, expected heterozygosity, maximum repeat length and/or chromosomal location of the repeat. Individual MS were analysed by step-down whole model linear regression and standard least squares fit models, both stratified by annual parasite incidence to identify MS markers with values significantly different from the mean.

RESULTS

Of the 42 MS under evaluation in this study, 18 (nine high priority) were identified as ideal candidates for measuring population diversity between global regions, while five (two high priority) additional markers were identified as candidates for MOI studies.

CONCLUSIONS

MS diversity was found to be a function of endemicity and motif structure. Evaluation of individual MS permitted the assembly of a refined panel of markers that can be reliably utilized in the field to compare population structures between global regions.

Distinct mutational behaviors differentiate short tandem repeats from microsatellites in the human genome

A tandem repeat's (TR) propensity to mutate increases with repeat number, and can become very pronounced beyond a critical boundary, transforming it into a microsatellite (MS). However, a clear understanding of the mutational behavior of different TR classes and motifs and related mechanisms is lacking, as is a consensus on the existence of a boundary separating short TRs (STRs) from MSs. This hinders our understanding of MSs' mutational properties and their effective use as genetic markers. Using indel calls for 179 individuals from 1000 Genomes Pilot-1 Project, we determined polymorphism incidence for four major TR classes, and formalized its varying relationship with repeat number using segmented regression. We observed a biphasic regime with a transition from a faster to a slower exponential growth at 9, 5, 4, and 4 repeats for mono-, di-, tri-, and tetranucleotide TRs, respectively. We used an in vitro mutagenesis assay to evaluate the contribution of strand slippage errors to mutability. STRs and MSs differ in their absolute polymorphism levels, but more importantly in their rates of mutability growth. Although strand slippage is a major factor driving mononucleotide polymorphism incidence, dinucleotide polymorphism incidence is greater than that expected due to strand slippage alone, indicating that additional cellular factors might be driving dinucleotide mutability in the human genome. Leveraging on hundreds of human genomes, we present the first comprehensive, genome-wide analysis of TR mutational behavior, encompassing several motif sizes and compositions.

Mature microsatellites: mechanisms underlying dinucleotide microsatellite mutational biases in human cells

Dinucleotide microsatellites are dynamic DNA sequences that affect genome stability. Here, we focused on mature microsatellites, defined as pure repeats of lengths above the threshold and unlikely to mutate below it in a single mutational event. We investigated the prevalence and mutational behavior of these sequences by using human genome sequence data, human cells in culture, and purified DNA polymerases. Mature dinucleotides (≥10 units) are present within exonic sequences of >350 genes, resulting in vulnerability to cellular genetic integrity. Mature dinucleotide mutagenesis was examined experimentally using ex vivo and in vitro approaches. We observe an expansion bias for dinucleotide microsatellites up to 20 units in length in somatic human cells, in agreement with previous computational analyses of germ-line biases. Using purified DNA polymerases and human cell lines deficient for mismatch repair (MMR), we show that the expansion bias is caused by functional MMR and is not due to DNA polymerase error biases. Specifically, we observe that the MutSα and MutLα complexes protect against expansion mutations. Our data support a model wherein different MMR complexes shift the balance of mutations toward deletion or expansion. Finally, we show that replication fork progression is stalled within long dinucleotides, suggesting that mutational mechanisms within long repeats may be distinct from shorter lengths, depending on the biochemistry of fork resolution. Our work combines computational and experimental approaches to explain the complex mutational behavior of dinucleotide microsatellites in humans.

Tumor-specific microsatellite instability: do distinct mechanisms underlie the MSI-L and EMAST phenotypes?

Microsatellite DNA sequences display allele length alterations or microsatellite instability (MSI) in tumor tissues, and MSI is used diagnostically for tumor detection and classification. We discuss the known types of tumor-specific MSI patterns and the relevant mechanisms underlying each pattern. Mutation rates of individual microsatellites vary greatly, and the intrinsic DNA features of motif size, sequence, and length contribute to this variation. MSI is used for detecting mismatch repair (MMR)-deficient tumors, which display an MSI-high phenotype due to genome-wide microsatellite destabilization. Because several pathways maintain microsatellite stability, tumors that have undergone other events associated with moderate genome instability may display diagnostic MSI only at specific di- or tetranucleotide markers. We summarize evidence for such alternative MSI forms (A-MSI) in sporadic cancers, also referred to as MSI-low and EMAST. While the existence of A-MSI is not disputed, there is disagreement about the origin and pathologic significance of this phenomenon. Although ambiguities due to PCR methods may be a source, evidence exists for other mechanisms to explain tumor-specific A-MSI. Some portion of A-MSI tumors may result from random mutational events arising during neoplastic cell evolution. However, this mechanism fails to explain the specificity of A-MSI for di- and tetranucleotide instability. We present evidence supporting the alternative argument that some A-MSI tumors arise by a distinct genetic pathway, and give examples of DNA metabolic pathways that, when altered, may be responsible for instability at specific microsatellite motifs. Finally, we suggest that A-MSI in tumors could be molecular signatures of environmental influences and DNA damage. Importantly, A-MSI occurs in several pre-neoplastic inflammatory states, including inflammatory bowel diseases, consistent with a role of oxidative stress in A-MSI. Understanding the biochemical basis of A-MSI tumor phenotypes will advance the development of new diagnostic tools and positively impact the clinical management of individual cancers.

DNA polymerase kappa microsatellite synthesis: two distinct mechanisms of slippage-mediated errors

Microsatellite tandem repeats are frequent sites of strand slippage mutagenesis in the human genome. Microsatellite mutations often occur as insertion/deletion of a repeat motif (unit-based indels), and increase in frequency with increasing repeat length after a threshold is reached. We recently demonstrated that DNA polymerase κ (Pol κ) produces fewer unit-based indel errors within dinucleotide microsatellites than does polymerase δ. Here, we examined human Pol κ's error profile within microsatellite alleles of varying sequence composition and length, using an in vitro HSV-tk gap-filling assay. We observed that Pol κ displays relatively accurate synthesis for unit-based indels, using di- and tetranucleotide repeat templates longer than the threshold length. We observed an abrupt increase in the unit-based indel frequency when the total microsatellite length exceeds 28 nucleotides, suggesting that extended Pol κ protein-DNA interactions enhance fidelity of the enzyme when synthesizing these microsatellite alleles. In contrast, Pol κ is error-prone within the HSV-tk coding sequence, producing frequent single-base errors in a manner that is highly biased with regard to sequence context. Single-nucleotide errors are also created by Pol κ within di- and tetranucleotide repeats, independently of the microsatellite allele length and at a frequency per nucleotide similar to the frequency of single base errors within the coding sequence. These single-base errors represent the mutational signature of Pol κ, and we propose them a mechanism independent of homology-stabilized slippage. Pol κ's dual fidelity nature provides a unique research tool to explore the distinct mechanisms of slippage-mediated mutagenesis.

Mechanism of replicative DNA polymerase delta pausing and a potential role for DNA polymerase kappa in common fragile site replication

Common fragile sites (CFSs) are hot spots of chromosomal breakage, and CFS breakage models involve perturbations of DNA replication. Here, we analyzed the contribution of specific repetitive DNA sequence elements within CFSs to the inhibition of DNA synthesis by replicative and specialized DNA polymerases (Pols). The efficiency of in vitro DNA synthesis was quantitated using templates corresponding to regions within FRA16D and FRA3B harboring AT-rich microsatellite and quasi-palindrome (QP) sequences. QPs were predicted to form stems of ~75-100% self-homology, separated by 3-9 bases of intervening sequences. Analysis of DNA synthesis progression by human Pol δ demonstrated significant synthesis perturbation both at [A](n) and [TA](n) repeats in a length-dependent manner and at short (<40 base pairs) QP sequences. DNA synthesis by the Y-family polymerase κ was significantly more efficient than Pol δ through both types of repetitive elements. Using DNA trap experiments, we show that Pol δ pauses within CFS sequences are sites of enzyme dissociation, and dissociation was observed in the presence of RFC-loaded PCNA. We propose that enrichment of microsatellite and QP elements at CFS regions contributes to fragility by perturbing replication through multiple mechanisms, including replicative Pol pausing and dissociation. Our finding that Pol δ dissociates at specific CFS sequences is significant, since dissociation of the replication machinery and inability to efficiently recover the replication fork can lead to fork collapse and/or formation of double-strand breaks in vivo. Our biochemical studies also extend the potential involvement of Y-family polymerases in CFS maintenance to include polymerase κ.

Beyond translesion synthesis: polymerase κ fidelity as a potential determinant of microsatellite stability

Microsatellite DNA synthesis represents a significant component of human genome replication that must occur faithfully. However, yeast replicative DNA polymerases do not possess high fidelity for microsatellite synthesis. We hypothesized that the structural features of Y-family polymerases that facilitate accurate translesion synthesis may promote accurate microsatellite synthesis. We compared human polymerases κ (Pol κ) and η (Pol η) fidelities to that of replicative human polymerase δ holoenzyme (Pol δ4), using the in vitro HSV-tk assay. Relative polymerase accuracy for insertion/deletion (indel) errors within 2-3 unit repeats internal to the HSV-tk gene concurred with the literature: Pol δ4 >> Pol κ or Pol η. In contrast, relative polymerase accuracy for unit-based indel errors within [GT](10) and [TC](11) microsatellites was: Pol κ ≥ Pol δ4 > Pol η. The magnitude of difference was greatest between Pols κ and δ4 with the [GT] template. Biochemically, Pol κ displayed less synthesis termination within the [GT] allele than did Pol δ4. In dual polymerase reactions, Pol κ competed with either a stalled or moving Pol δ4, thereby reducing termination. Our results challenge the ideology that pol κ is error prone, and suggest that DNA polymerases with complementary biochemical properties can function cooperatively at repetitive sequences.

The in vitro fidelity of yeast DNA polymerase δ and polymerase ε holoenzymes during dinucleotide microsatellite DNA synthesis

Elucidating the sources of genetic variation within microsatellite alleles has important implications for understanding the etiology of human diseases. Mismatch repair is a well described pathway for the suppression of microsatellite instability. However, the cellular polymerases responsible for generating microsatellite errors have not been fully described. We address this gap in knowledge by measuring the fidelity of recombinant yeast polymerase δ (Pol δ) and ɛ (Pol ɛ) holoenzymes during synthesis of a [GT/CA] microsatellite. The in vitro HSV-tk forward assay was used to measure DNA polymerase errors generated during gap-filling of complementary GT(10) and CA(10)-containing substrates and ∼90 nucleotides of HSV-tk coding sequence surrounding the microsatellites. The observed mutant frequencies within the microsatellites were 4 to 30-fold higher than the observed mutant frequencies within the coding sequence. More specifically, the rate of Pol δ and Pol ɛ misalignment-based insertion/deletion errors within the microsatellites was ∼1000-fold higher than the rate of insertion/deletion errors within the HSV-tk gene. Although the most common microsatellite error was the deletion of a single repeat unit, ∼ 20% of errors were deletions of two or more units for both polymerases. The differences in fidelity for wild type enzymes and their exonuclease-deficient derivatives were ∼2-fold for unit-based microsatellite insertion/deletion errors. Interestingly, the exonucleases preferentially removed potentially stabilizing interruption errors within the microsatellites. Since Pol δ and Pol ɛ perform not only the bulk of DNA replication in eukaryotic cells but also are implicated in performing DNA synthesis associated with repair and recombination, these results indicate that microsatellite errors may be introduced into the genome during multiple DNA metabolic pathways.

Telomeric repeat mutagenicity in human somatic cells is modulated by repeat orientation and G-quadruplex stability

Telomeres consisting of tandem guanine-rich repeats can form secondary DNA structures called G-quadruplexes that represent potential targets for DNA repair enzymes. While G-quadruplexes interfere with DNA synthesis in vitro, the impact of G-quadruplex formation on telomeric repeat replication in human cells is not clear. We investigated the mutagenicity of telomeric repeats as a function of G-quadruplex folding opportunity and thermal stability using a shuttle vector mutagenesis assay. Since single-stranded DNA during lagging strand replication increases the opportunity for G-quadruplex folding, we tested vectors with G-rich sequences on the lagging versus the leading strand. Contrary to our prediction, vectors containing human [TTAGGG]₁₀ repeats with a G-rich lagging strand were significantly less mutagenic than vectors with a G-rich leading strand, after replication in normal human cells. We show by UV melting experiments that G-quadruplexes from ciliates [TTGGGG]₄ and [TTTTGGGG]₄ are thermally more stable compared to human [TTAGGG]₄. Consistent with this, replication of vectors with ciliate [TTGGGG]₁₀ repeats yielded a 3-fold higher mutant rate compared to the human [TTAGGG]₁₀ vectors. Furthermore, we observed significantly more mutagenic events in the ciliate repeats compared to the human repeats. Our data demonstrate that increased G-quadruplex opportunity (repeat orientation) in human telomeric repeats decreased mutagenicity, while increased thermal stability of telomeric G-quadruplexes was associated with increased mutagenicity.

What is a microsatellite: a computational and experimental definition based upon repeat mutational behavior at A/T and GT/AC repeats

Microsatellites are abundant in eukaryotic genomes and have high rates of strand slippage-induced repeat number alterations. They are popular genetic markers, and their mutations are associated with numerous neurological diseases. However, the minimal number of repeats required to constitute a microsatellite has been debated, and a definition of a microsatellite that considers its mutational behavior has been lacking. To define a microsatellite, we investigated slippage dynamics for a range of repeat sizes, utilizing two approaches. Computationally, we assessed length polymorphism at repeat loci in ten ENCODE regions resequenced in four human populations, assuming that the occurrence of polymorphism reflects strand slippage rates. Experimentally, we determined the in vitro DNA polymerase-mediated strand slippage error rates as a function of repeat number. In both approaches, we compared strand slippage rates at tandem repeats with the background slippage rates. We observed two distinct modes of mutational behavior. At small repeat numbers, slippage rates were low and indistinguishable from background measurements. A marked transition in mutability was observed as the repeat array lengthened, such that slippage rates at large repeat numbers were significantly higher than the background rates. For both mononucleotide and dinucleotide microsatellites studied, the transition length corresponded to a similar number of nucleotides (approximately 10). Thus, microsatellite threshold is determined not by the presence/absence of strand slippage at repeats but by an abrupt alteration in slippage rates relative to background. These findings have implications for understanding microsatellite mutagenesis, standardization of genome-wide microsatellite analyses, and predicting polymorphism levels of individual microsatellite loci.

DNA slippage occurs at microsatellite loci without minimal threshold length in humans: a comparative genomic approach

The dynamics of microsatellite, or short tandem repeats (STRs), is well documented for long, polymorphic loci, but much less is known for shorter ones. For example, the issue of a minimum threshold length for DNA slippage remains contentious. Model-fitting methods have generally concluded that slippage only occurs over a threshold length of about eight nucleotides, in contradiction with some direct observations of tandem duplications at shorter repeated sites. Using a comparative analysis of the human and chimpanzee genomes, we examined the mutation patterns at microsatellite loci with lengths as short as one period plus one nucleotide. We found that the rates of tandem insertions and deletions at microsatellite loci strongly deviated from background rates in other parts of the human genome and followed an exponential increase with STR size. More importantly, we detected no lower threshold length for slippage. The rate of tandem duplications at unrepeated sites was higher than expected from random insertions, providing evidence for genome-wide action of indel slippage (an alternative mechanism generating tandem repeats). The rate of point mutations adjacent to STRs did not differ from that estimated elsewhere in the genome, except around dinucleotide loci. Our results suggest that the emergence of STR depends on DNA slippage, indel slippage, and point mutations. We also found that the dynamics of tandem insertions and deletions differed in both rates and size at which these mutations take place. We discuss these results in both evolutionary and mechanistic terms.

Every microsatellite is different: Intrinsic DNA features dictate mutagenesis of common microsatellites present in the human genome

Microsatellite sequences are ubiquitous in the human genome and are important regulators of genome function. Here, we examine the mutational mechanisms governing the stability of highly abundant mono-, di-, and tetranucleotide microsatellites. Microsatellite mutation rate estimates from pedigree analyses and experimental models range from a low of approximately 10(-6) to a high of approximately 10(-2) mutations per locus per generation. The vast majority of observed mutational variation can be attributed to features intrinsic to the allele itself, including motif size, length, and sequence composition. A greater than linear relationship between motif length and mutagenesis has been observed in several model systems. Motif sequence differences contribute up to 10-fold to the variation observed in human cell mutation rates. The major mechanism of microsatellite mutagenesis is strand slippage during DNA synthesis. DNA polymerases produce errors within microsatellites at a frequency that is 10- to 100-fold higher than the frequency of frameshifts in coding sequences. Motif sequence significantly affects both polymerase error rate and specificity, resulting in strand biases within complementary microsatellites. Importantly, polymerase errors within microsatellites include base substitutions, deletions, and complex mutations, all of which produced interrupted alleles from pure microsatellites. Postreplication mismatch repair efficiency is affected by microsatellite motif size and sequence, also contributing to the observed variation in microsatellite mutagenesis. Inhibition of DNA synthesis within common microsatellites is highly sequence-dependent, and is positively correlated with the production of errors. DNA secondary structure within common microsatellites can account for some DNA polymerase pause sites, and may be an important factor influencing mutational specificity.

Role of TLS DNA polymerases eta and kappa in processing naturally occurring structured DNA in human cells

Accurate DNA replication during S-phase is fundamental to maintain genome integrity. During this critical process, replication forks frequently encounter obstacles that impede their progression. While the regulatory pathways which act in response to exogenous replication stress are beginning to emerge, the mechanisms by which fork integrity is maintained at naturally occurring endogenous replication-impeding sequences remains obscure. Notably, little is known about how cells replicate through special chromosomal regions containing structured non-B DNA, for example, G4 quartets, known to hamper fork progression or trigger chromosomal rearrangements. Here, we have investigated the role in this process of the human translesion synthesis (TLS) DNA polymerases of the Y-family (pol eta, pol iota, and pol kappa), specialized enzymes known to synthesize DNA through DNA damage. We show that depletion by RNA interference of expression of the genes for Pol eta or Pol kappa, but not Pol iota, sensitizes U2OS cells treated with the G4-tetraplex interactive compound telomestatin and triggers double-strand breaks in HeLa cells harboring multiple copies of a G-rich sequence from the promoter region of the human c-MYC gene, chromosomally integrated as a transgene. Moreover, we found that downregulation of Pol kappa only raises the level of DSB in HeLa cells containing either one of two breakage hotspot structured DNA sequences in the chromosome, the major break region (Mbr) of BCL-2 gene and the GA rich region from the far right-hand end of the genome of the Kaposi Sarcoma associated Herpesvirus. These data suggest that naturally occurring DNA structures are physiological substrates of both pol eta and pol kappa. We discuss these data in the light of their downregulation in human cancers.

Models for chromosomal replication-independent non-B DNA structure-induced genetic instability

Regions of genomic DNA containing repetitive nucleotide sequences can adopt a number of different structures in addition to the canonical B-DNA form: many of these non-B DNA structures are causative factors in genetic instability and human disease. Although chromosomal DNA replication through such repetitive sequences has been considered a major cause of non-B form DNA structure-induced genetic instability, it is also observed in non-proliferative tissues. In this review, we discuss putative mechanisms responsible for the mutagenesis induced by non-B DNA structures in the absence of chromosomal DNA replication.

DNA triple helices: biological consequences and therapeutic potential

DNA structure is a critical element in determining its function. The DNA molecule is capable of adopting a variety of non-canonical structures, including three-stranded (i.e. triplex) structures, which will be the focus of this review. The ability to selectively modulate the activity of genes is a long-standing goal in molecular medicine. DNA triplex structures, either intermolecular triplexes formed by binding of an exogenously applied oligonucleotide to a target duplex sequence, or naturally occurring intramolecular triplexes (H-DNA) formed at endogenous mirror repeat sequences, present exploitable features that permit site-specific alteration of the genome. These structures can induce transcriptional repression and site-specific mutagenesis or recombination. Triplex-forming oligonucleotides (TFOs) can bind to duplex DNA in a sequence-specific fashion with high affinity, and can be used to direct DNA-modifying agents to selected sequences. H-DNA plays important roles in vivo and is inherently mutagenic and recombinogenic, such that elements of the H-DNA structure may be pharmacologically exploitable. In this review we discuss the biological consequences and therapeutic potential of triple helical DNA structures. We anticipate that the information provided will stimulate further investigations aimed toward improving DNA triplex-related gene targeting strategies for biotechnological and potential clinical applications.

High frequency of microsatellites in S. cerevisiae meiotic recombination hotspots

Background

Microsatellites are highly abundant in eukaryotic genomes but their function and evolution are not yet well understood. Their elevated mutation rate makes them ideal markers of genetic difference, but high levels of unexplained heterogeneity in mutation rates among microsatellites at different genomic locations need to be elucidated in order to improve the power and accuracy of the many types of study that use them as genetic markers. Recombination could contribute to this heterogeneity, since while replication errors are thought to be the predominant mechanism for microsatellite mutation, meiotic recombination is involved in some mutation events. There is also evidence suggesting that microsatellites could function as recombination signals. The yeast S. cerevisiae is a useful model organism with which to further explore the link between microsatellites and recombination, since it is very amenable to genetic study, and meiotic recombination hotspots have been mapped throughout its entire genome.

Results

We examined in detail the relationship between microsatellites and hotspots of meiotic double-strand breaks, the precursors of meiotic recombination, throughout the S. cerevisiae genome. We included all tandem repeats with motif length (repeat period) between one and six base pairs. Long, short and two-copy arrays were considered separately. We found that long, mono-, di- and trinucleotide microsatellites are around twice as frequent in hot than non-hot intergenic regions. The associations are weak or absent for repeats with less than six copies, and also for microsatellites with 4–6 base pair motifs, but high-copy arrays with motif length greater than three are relatively very rare throughout the genome. We present evidence that the association between high-copy, short-motif microsatellites and recombination hotspots is not driven by effects on microsatellite distribution of other factors previously linked to both recombination and microsatellites, including transcription, GC-content and transposable elements.

Conclusion

Our findings suggest that a mutation bias relating to recombination hotspots causing repeats to form and grow, and/or regulation of a subset of hotspots by simple sequences, may be significant processes in yeast. Some previous evidence has cast doubt on both of these possibilities, and as a result they have not been explored on a large scale, but the strength of the association we report suggests that they deserve further experimental testing.

DNA polymerase kappa produces interrupted mutations and displays polar pausing within mononucleotide microsatellite sequences

Microsatellites are ubiquitously present in eukaryotic genomes and are implicated as positive factors in evolution. At the nucleotide level, microsatellites undergo slippage events that alter allele length and base changes that interrupt the repetitive tract. We examined DNA polymerase errors within a [T]₁₁ microsatellite using an in vitro assay that preferentially detects mutations other than unit changes. We observed that human DNA polymerase kappa (Pol κ) inserts dGMP and dCMP within the [T]₁₁ mononucleotide repeat, producing an interrupted 12-bp allele. Polymerase β produced such interruptions at a lower frequency. These data demonstrate that DNA polymerases are capable of directly producing base interruptions within microsatellites. At the molecular level, expanded microsatellites have been implicated in DNA replication fork stalling. Using an in vitro primer extension assay, we observed sequence-specific synthesis termination by DNA polymerases within mononucleotides. Quantitatively, intense, polar pausing was observed for both pol κ and polymerase α-primase within a [T]₁₁ allele. A mechanism is proposed in which pausing results from DNA bending within the duplex stem of the nascent DNA. Our data support the concept of a microsatellite life-cycle, and are consistent with the models in which DNA sequence or secondary structures contributes to non-uniform rates of replication fork progression.

Telomere ResQue and preservation--roles for the Werner syndrome protein and other RecQ helicases

Werner syndrome is an autosomal recessive disorder resulting from loss of function of the RecQ helicase, WRN protein. WS patients prematurely develop numerous clinical symptoms and diseases associated with aging early in life and are predisposed to cancer. WRN protein and many other RecQ helicases in general, seem to function during DNA replication in the processing of stalled replication forks. Genetic, cellular and biochemical evidence support roles for WRN in proper replication and repair of telomeric DNA, and indicate that telomere dysfunction contributes to the WS disease pathology.

Error-free bypass of 2-hydroxyadenine by human DNA polymerase lambda with Proliferating Cell Nuclear Antigen and Replication Protein A in different sequence contexts

1,2-dihydro-2-oxoadenine (2-OH-A), a common DNA lesion produced by reactive oxygen species, is a strong replicative block for several DNA polymerases (DNA pols). We have previously shown that various bases can be misincorporated opposite the 2-OH-A lesion and the type of mispairs varies with either the sequence context or the type of DNA pol tested. Here, we have analysed the ability of the human pol family X member DNA pol λ, to bypass the 2-OH-A lesion. DNA pol λ can perform error-free bypass of 2-OH-A when this lesion is located in a random sequence, whereas in a repeated sequence context, even though bypass was also largely error-free, misincorporation of dGMP could be observed. The fidelity of translesion synthesis of 2-OH-A in a repeated sequence by DNA pol λ was enhanced by the auxiliary proteins Proliferating Cell Nuclear Antigen (PCNA) and Replication Protein A (RP-A). We also found that the DNA pol λ active site residue tyrosine 505 determined the nucleotide selectivity opposite 2-OH-A. Our data show, for the first time, that the 2-OH-A lesion can be efficiently and faithfully bypassed by a human DNA pol λ in combination with PCNA and RP-A.

Replication fork regression in repetitive DNAs

Among several different types of repetitive sequences found in the human genome, this study has examined the telomeric repeat, necessary for the protection of chromosome termini, and the disease-associated triplet repeat (CTG)·(CAG)_n. Evidence suggests that replication of both types of repeats is problematic and that a contributing factor is the repetitive nature of the DNA itself. Here we have used electron microscopy to investigate DNA structures formed at replication forks on large model DNAs containing these repeat sequences, in an attempt to elucidate the contributory effect that these repetitive DNAs may have on their replication. Visualization of the DNA revealed that there is a high propensity for a paused replication fork to spontaneously regress when moving through repetitive DNAs, and that this results in a four-way chickenfoot intermediate that could present a significant block to replication in vivo, possibly leading to unwanted recombination events, amplifications or deletions.

Non-B DNA structure-induced genetic instability

Repetitive DNA sequences are abundant in eukaryotic genomes, and many of these sequences have the potential to adopt non-B DNA conformations. Genes harboring non-B DNA structure-forming sequences increase the risk of genetic instability and thus are associated with human diseases. In this review, we discuss putative mechanisms responsible for genetic instability events occurring at these non-B DNA structures, with a focus on hairpins, left-handed Z-DNA, and intramolecular triplexes or H-DNA. Slippage and misalignment are the most common events leading to DNA structure-induced mutagenesis. However, a number of other mechanisms of genetic instability have been proposed based on the finding that these structures not only induce expansions and deletions, but can also induce DNA strand breaks and rearrangements. The available data implicate a variety of proteins, such as mismatch repair proteins, nucleotide excision repair proteins, topoisomerases, and structure specific-nucleases in the processing of these mutagenic DNA structures. The potential mechanisms of genetic instability induced by these structures and their contribution to human diseases are discussed.

Evidence for interplay among yeast replicative DNA polymerases alpha, delta and epsilon from studies of exonuclease and polymerase active site mutations

BACKGROUND

DNA polymerase epsilon (Pol epsilon) is essential for S-phase replication, DNA damage repair and checkpoint control in yeast. A pol2-Y831A mutation leading to a tyrosine to alanine change in the Pol epsilon active site does not cause growth defects and confers a mutator phenotype that is normally subtle but strong in a mismatch repair-deficient strain. Here we investigate the mechanism responsible for the mutator effect.

RESULTS

Purified four-subunit Y831A Pol epsilon turns over more deoxynucleoside triphosphates to deoxynucleoside monophosphates than does wild-type Pol epsilon, suggesting altered coordination between the polymerase and exonuclease active sites. The pol2-Y831A mutation suppresses the mutator effect of the pol2-4 mutation in the exonuclease active site that abolishes proofreading by Pol epsilon, as measured in haploid strain with the pol2-Y831A,4 double mutation. Analysis of mutation rates in diploid strains reveals that the pol2-Y831A allele is recessive to pol2-4. In addition, the mutation rates of strains with the pol2-4 mutation in combination with active site mutator mutations in Pol delta and Pol alpha suggest that Pol epsilon may proofread certain errors made by Pol alpha and Pol delta during replication in vivo.

CONCLUSIONS

Our data suggest that Y831A replacement in Pol epsilon reduces replication fidelity and its participation in chromosomal replication, but without eliminating an additional function that is essential for viability. This suggests that other polymerases can substitute for certain functions of polymerase epsilon.