Mutation rate and specificity analysis of tetranucleotide microsatellite DNA alleles in somatic human cells

Mismatch repair is a double-edged sword in the battle against microsatellite instability

Roughly 3% of the human genome consists of microsatellites or tracts of short tandem repeats (STRs). These STRs are often unstable, undergoing high-frequency expansions (increases) or contractions (decreases) in the number of repeat units. Some microsatellite instability (MSI) is seen at multiple STRs within a single cell and is associated with certain types of cancer. A second form of MSI is characterised by expansion of a single gene-specific STR and such expansions are responsible for a group of 40+ human genetic disorders known as the repeat expansion diseases (REDs). While the mismatch repair (MMR) pathway prevents genome-wide MSI, emerging evidence suggests that some MMR factors are directly involved in generating expansions in the REDs. Thus, MMR suppresses some forms of expansion while some MMR factors promote expansion in other contexts. This review will cover what is known about the paradoxical effect of MMR on microsatellite expansion in mammalian cells.

Genome-wide characterization of microsatellite DNA in fishes: survey and analysis of their abundance and frequency in genome-specific regions

BACKGROUND

Microsatellite repeats are ubiquitous in organism genomes and play an important role in the chromatin organization, regulation of gene activity, recombination and DNA replication. Although microsatellite distribution patterns have been studied in most phylogenetic lineages, they are unclear in fish species.

RESULTS

Here, we present the first systematic examination of microsatellite distribution in coding and non-coding regions of 14 fish genomes. Our study showed that the number and type of microsatellites displayed nonrandom distribution for both intragenic and intergenic regions, suggesting that they have potential roles in transcriptional or translational regulation and DNA replication slippage theories alone were insufficient to explain the distribution patterns. Our results showed that microsatellites are dominant in non-coding regions. The total number of microsatellites ranged from 78,378 to 1,012,084, and the relative density varied from 4925.76 bp/Mb to 25,401.97 bp/Mb. Overall, (A + T)-rich repeats were dominant. The dependence of repeat abundance on the length of the repeated unit (1-6 nt) showed a great similarity decrease, whereas more tri-nucleotide repeats were found in exonic regions than tetra-nucleotide repeats of most species. Moreover, the incidence of different repeated types appeared species- and genomic-specific. These results highlight potential mechanisms for maintaining microsatellite distribution, such as selective forces and mismatch repair systems.

CONCLUSIONS

Our data could be beneficial for the studies of genome evolution and microsatellite DNA evolutionary dynamics, and facilitate the exploration of microsatellites structural, function, composition mode and molecular markers development in these species.

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.

A time- and cost-effective strategy to sequence mammalian Y Chromosomes: an application to the de novo assembly of gorilla Y

The mammalian Y Chromosome sequence, critical for studying male fertility and dispersal, is enriched in repeats and palindromes, and thus, is the most difficult component of the genome to assemble. Previously, expensive and labor-intensive BAC-based techniques were used to sequence the Y for a handful of mammalian species. Here, we present a much faster and more affordable strategy for sequencing and assembling mammalian Y Chromosomes of sufficient quality for most comparative genomics analyses and for conservation genetics applications. The strategy combines flow sorting, short- and long-read genome and transcriptome sequencing, and droplet digital PCR with novel and existing computational methods. It can be used to reconstruct sex chromosomes in a heterogametic sex of any species. We applied our strategy to produce a draft of the gorilla Y sequence. The resulting assembly allowed us to refine gene content, evaluate copy number of ampliconic gene families, locate species-specific palindromes, examine the repetitive element content, and produce sequence alignments with human and chimpanzee Y Chromosomes. Our results inform the evolution of the hominine (human, chimpanzee, and gorilla) Y Chromosomes. Surprisingly, we found the gorilla Y Chromosome to be similar to the human Y Chromosome, but not to the chimpanzee Y Chromosome. Moreover, we have utilized the assembled gorilla Y Chromosome sequence to design genetic markers for studying the male-specific dispersal of this endangered species.

Recent insertion/deletion (reINDEL) mutations: increasing awareness to boost molecular-based research in ecology and evolution

Today, the comparative analysis of DNA molecules mainly uses information inferred from nucleotide substitutions. Insertion/deletion (INDEL) mutations, in contrast, are largely considered uninformative and discarded, due to our lacking knowledge on their evolution. However, including rather than discarding INDELs would be relevant to any research area in ecology and evolution that uses molecular data. As a practical approach to better understanding INDEL evolution in general, we propose the study of recent INDEL (reINDEL) mutations – mutations where both ancestral and derived state are seen in the sample. The precondition for reINDEL identification is knowledge about the pedigree of the individuals sampled. Sound reINDEL knowledge will allow the improved modeling needed for including INDELs in the downstream analysis of molecular data. Both microsatellites, currently still the predominant marker system in the analysis of populations, and sequences generated by next-generation sequencing, a promising and rapidly developing range of technologies, offer the opportunity for reINDEL identification. However, a 2013 sample of animal microsatellite studies contained unexpectedly few reINDELs identified. As most likely explanation, we hypothesize that reINDELs are underreported rather than absent and that this underreporting stems from common reINDEL unawareness. If our hypothesis applies, increased reINDEL awareness should allow gathering data rapidly. We recommend the routine reporting of either the absence or presence of reINDELs together with standardized key information on the nature of mutations when they are detected and the use of the keyword “reINDEL” to increase visibility in both instances of successful and unsuccessful search.

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.

Clinical and biochemical function of polymorphic NR0B1 GGAA-microsatellites in Ewing sarcoma: a report from the Children's Oncology Group

BACKGROUND

The genetics involved in Ewing sarcoma susceptibility and prognosis are poorly understood. EWS/FLI and related EWS/ETS chimeras upregulate numerous gene targets via promoter-based GGAA-microsatellite response elements. These microsatellites are highly polymorphic in humans, and preliminary evidence suggests EWS/FLI-mediated gene expression is highly dependent on the number of GGAA motifs within the microsatellite.

OBJECTIVES

Here we sought to examine the polymorphic spectrum of a GGAA-microsatellite within the NR0B1 promoter (a critical EWS/FLI target) in primary Ewing sarcoma tumors, and characterize how this polymorphism influences gene expression and clinical outcomes.

RESULTS

A complex, bimodal pattern of EWS/FLI-mediated gene expression was observed across a wide range of GGAA motifs, with maximal expression observed in constructs containing 20-26 GGAA motifs. Relative to white European and African controls, the NR0B1 GGAA-microsatellite in tumor cells demonstrated a strong bias for haplotypes containing 21-25 GGAA motifs suggesting a relationship between microsatellite function and disease susceptibility. This selection bias was not a product of microsatellite instability in tumor samples, nor was there a correlation between NR0B1 GGAA-microsatellite polymorphisms and survival outcomes.

CONCLUSIONS

These data suggest that GGAA-microsatellite polymorphisms observed in human populations modulate EWS/FLI-mediated gene expression and may influence disease susceptibility in Ewing sarcoma.

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.

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.

Inhibition of DNA synthesis facilitates expansion of low-complexity repeats: is strand slippage stimulated by transient local depletion of specific dNTPs?

Simple DNA repeats (trinucleotide repeats, micro- and minisatellites) are prone to expansion/contraction via formation of secondary structures during DNA synthesis. Such structures both inhibit replication forks and create opportunities for template-primer slippage, making these repeats unstable. Certain aspects of simple repeat instability, however, suggest additional mechanisms of replication inhibition dependent on the primary DNA sequence, rather than on secondary structure formation. I argue that expanded simple repeats, due to their lower DNA complexity, should transiently inhibit DNA synthesis by locally depleting specific DNA precursors. Such transient inhibition would promote formation of secondary structures and would stabilize these structures, facilitating strand slippage. Thus, replication problems at simple repeats could be explained by potentiated toxicity, where the secondary structure-driven repeat instability is enhanced by DNA polymerase stalling at the low complexity template DNA.

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.

Factors affecting germline mutations in a hypervariable microsatellite: a comparative analysis of six species of swallows (Aves: Hirundinidae)

Microsatellites mutate frequently by replication slippage. Empirical evidence shows that the probability of such slippage mutations may increase with the length of the repeat region as well as exposure to environmental mutagens, but the mutation rate can also differ between the male and female germline. It has been hypothesized that more intense sexual selection or sperm competition can also lead to elevated mutation rates, but the empirical evidence is inconclusive. Here, we analyzed the occurrence of germline slippage mutations in the hypervariable pentanucleotide microsatellite locus HrU10 across six species of swallow (Aves: Hirundinidae). These species exhibit marked differences in the length range of the microsatellite, as well as differences in the intensity of sperm competition. We found a strong effect of microsatellite length on the probability of mutation, but no residual effect of species or their level of sperm competition when the length effect was accounted for. Neither could we detect any difference in mutation rate between tree swallows (Tachycineta bicolor) breeding in Hamilton Harbour, Ontario, an industrial site with previous documentation of elevated mutation rates for minisatellite DNA, and a rural reference population. However, our cross-species analysis revealed two significant patterns of sex differences in HrU10 germline mutations: (1) mutations in longer alleles occurred typically in the male germline, those in shorter alleles in the female germline, and (2) male germline mutations were more often expansions than contractions, whereas no directional bias was evident in the female germline. These results indicate some fundamental differences in male and female gametogenesis affecting the probability of slippage mutations. Our study also reflects the value of a comparative, multi-species approach for locus-specific mutation analyses, through which a wider range of influential factors can be assessed than in single-species studies.

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.

Defective mismatch repair, microsatellite mutation bias, and variability in clinical cancer phenotypes

Microsatellite instability is associated with 10% to 15% of colorectal, endometrial, ovarian, and gastric cancers, and has long been used as a diagnostic tool for hereditary nonpolyposis colorectal carcinoma-related cancers. Tumor-specific length alterations within microsatellites are generally accepted to be a consequence of strand slippage events during DNA replication, which are uncorrected due to a defective postreplication mismatch repair (MMR) system. Mutations arising within microsatellites associated with critical target genes are believed to play a causative role in the evolution of MMR-defective tumors. In this review, we summarize current evidence of mutational biases within microsatellites arising as a consequence of intrinsic DNA sequence effects as well as variation in MMR efficiency. Microsatellite mutational biases are generally not considered during clinical testing; however, we suggest that such biases may be clinically significant as a factor contributing to phenotypic variation among microsatellite instability-positive tumors.

Sequence determinants of human microsatellite variability

BACKGROUND

Microsatellite loci are frequently used in genomic studies of DNA sequence repeats and in population studies of genetic variability. To investigate the effect of sequence properties of microsatellites on their level of variability we have analyzed genotypes at 627 microsatellite loci in 1,048 worldwide individuals from the HGDP-CEPH cell line panel together with the DNA sequences of these microsatellites in the human RefSeq database.

RESULTS

Calibrating PCR fragment lengths in individual genotypes by using the RefSeq sequence enabled us to infer repeat number in the HGDP-CEPH dataset and to calculate the mean number of repeats (as opposed to the mean PCR fragment length), under the assumption that differences in PCR fragment length reflect differences in the numbers of repeats in the embedded repeat sequences. We find the mean and maximum numbers of repeats across individuals to be positively correlated with heterozygosity. The size and composition of the repeat unit of a microsatellite are also important factors in predicting heterozygosity, with tetra-nucleotide repeat units high in G/C content leading to higher heterozygosity. Finally, we find that microsatellites containing more separate sets of repeated motifs generally have higher heterozygosity.

CONCLUSIONS

These results suggest that sequence properties of microsatellites have a significant impact in determining the features of human microsatellite variability.

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.

Human postmeiotic segregation 2 exhibits biased repair at tetranucleotide microsatellite sequences

The mismatch repair (MMR) system plays a major role in removing DNA polymerization errors, and loss of this pathway results in hereditary cancers characterized by microsatellite instability. We investigated microsatellite stability during DNA replication within human postmeiotic segregation 2 (hPMS2)-deficient and proficient human lymphoblastoid cell lines. Using a shuttle vector assay, we measured mutation rates at reporter cassettes containing defined mononucleotide, dinucleotide, and tetranucleotide microsatellite sequences. A mutator phenotype was observed in the hPMS2-deficient cell line. The mutation rate of vectors containing [G/C](10) or [GT/CA](10) alleles was elevated 20-fold to 40-fold in hPMS2-deficient cells, relative to an hPMS2-expressing cell line. We observed a 6-fold and 12-fold relative increase in mutation rate of [TTTC/AAAG](9) and [TTCC/AAGG](9) sequences, respectively, in hPMS2-deficient cells. Mutational specificity analyses suggested that repair by hPMS2 is biased. In the absence of hPMS2, a greater number of microsatellite expansion versus deletion mutations was observed, and expansion rates of the tetranucleotide alleles were similar. In the presence of hPMS2, we observed a 29-fold decrease in the [TTCC/AAGG](9) expansion rate but only a 6-fold decrease for the [TTTC/AAAG](9) allele. Our data indicate that hPMS2 is more protective of tetranucleotide expansions than deletions and that hPMS2 displays a sequence bias, wherein [TTCC/AAGG] sequences are stabilized to a greater extent than [TTTC/AAAG]. Our results allow for greater accuracy during identification of MMR defects by providing a mutational signature characteristic of hPMS2 defect. This study also provides clues to possible mechanisms of repair by hPMS2 in the context of the MMR system.

RecQ and RecG helicases have distinct roles in maintaining the stability of polypurine.polypyrimidine sequences

DNA triplex structures can block the replication fork and result in double-stranded DNA breaks (DSBs). RecQ and RecG helicases may be important for replication of such sequences as RecQ resolves synthetic triplex DNA structures and RecG mediates replication restart by fork regression. Primer extension on an 88 bp triplex-forming polypurine.polypyrimidine (Pu.Py) tract from the PKD1 gene demonstrated that RecQ, but not RecG, facilitated primer extension by T7 DNA polymerase. A high-throughput, dual plasmid screening system using isogenic bacterial lines deficient in RecG, RecQ, or both, revealed that RecQ deficiency increased mutation to sequence flanking this 88 bp tract by eight to ten-fold. Although RecG facilitated small deletions in an 88 bp mirror repeat-containing sequence, it was absolutely required to maintain a 2.5 kb Pu.Py tract containing multiple mirror repeats. These results support a two-tiered model where RecQ facilitates fork progression through triplex-forming tracts and, failing processivity, RecG is critical for replication fork restart.

Microsatellite evolution: Mutations, sequence variation, and homoplasy in the hypervariable avian microsatellite locus HrU10

Background

Microsatellites are frequently used genetic markers in a wide range of applications, primarily due to their high length polymorphism levels that can easily be genotyped by fragment length analysis. However, the mode of microsatellite evolution is yet not fully understood, and the role of interrupting motifs for the stability of microsatellites remains to be explored in more detail. Here we present a sequence analysis of mutation events and a description of the structure of repeated regions in the hypervariable, pentanucleotide microsatellite locus HrU10 in barn swallows (Hirundo rustica) and tree swallows (Tachycineta bicolor).

Results

In a large-scale parentage analysis in barn swallows and tree swallows, broods were screened for mutations at the HrU10 locus. In 41 cases in the barn swallows and 15 cases in the tree swallows, mutations corresponding to the loss or gain of one or two repeat units were detected. The parent and mutant offspring alleles were sequenced for 33 of these instances (26 in barn swallows and 7 in tree swallows). Replication slippage was considered the most likely mutational process. We tested the hypothesis that HrU10, a microsatellite with a wide allele size range, has an increased probability of introductions of interruptive motifs (IMs) with increasing length of the repeated region. Indeed, the number and length of the IMs was strongly positively correlated with the total length of the microsatellite. However, there was no significant correlation with the length of the longest stretch of perfectly repeated units, indicating a threshold level for the maximum length of perfectly repeated pentanucleotide motifs in stable HrU10 alleles. The combination of sequence and pedigree data revealed that 15 barn swallow mutations (58%) produced alleles that were size homoplasic to other alleles in the data set.

Conclusion

Our results give further insights into the mode of microsatellite evolution, and support the assumption of increased slippage rate with increased microsatellite length and a stabilizing effect of interrupting motifs for microsatellite regions consisting of perfect repeats. In addition, the observed extent of size homoplasy may impose a general caution against using hypervariable microsatellites in genetic diversity measures when alleles are identified by fragment length analysis only.

The genome-wide determinants of human and chimpanzee microsatellite evolution

Mutation rates of microsatellites vary greatly among loci. The causes of this heterogeneity remain largely enigmatic yet are crucial for understanding numerous human neurological diseases and genetic instability in cancer. In this first genome-wide study, the relative contributions of intrinsic features and regional genomic factors to the variation in mutability among orthologous human-chimpanzee microsatellites are investigated with resampling and regression techniques. As a result, we uncover the intricacies of microsatellite mutagenesis as follows. First, intrinsic features (repeat number, length, and motif size), which all influence the probability and rate of slippage, are the strongest predictors of mutability. Second, mutability increases nonuniformly with length, suggesting that processes additional to slippage, such as faulty repair, contribute to mutations. Third, mutability varies among microsatellites with different motif composition likely due to dissimilarities in secondary DNA structure formed by their slippage intermediates. Fourth, mutability of mononucleotide microsatellites is impacted by their location on sex chromosomes vs. autosomes and inside vs. outside of Alu repeats, the former confirming the importance of replication and the latter suggesting a role for gene conversion. Fifth, transcription status and location in a particular isochore do not influence microsatellite mutability. Sixth, compared with intrinsic features, regional genomic factors have only minor effects. Finally, our regression models explain approximately 90% of variation in microsatellite mutability and can generate useful predictions for the studies of human diseases, forensics, and conservation genetics.

Escherichia coli DNA polymerase IV contributes to spontaneous mutagenesis at coding sequences but not microsatellite alleles

Slipped strand mispairing during DNA synthesis is one proposed mechanism for microsatellite or short tandem repeat (STR) mutation. However, the DNA polymerase(s) responsible for STR mutagenesis have not been determined. In this study, we investigated the effect of the Escherichia colidinB gene product (Pol IV) on mononucleotide and dinucleotide repeat stability, using an HSV-tk gene episomal reporter system for microsatellite mutations. For the control vector (HSV-tk gene only) we observed a statistically significant 3.5-fold lower median mutation frequency in dinB(-) than dinB(+) cells (p<0.001, Wilcoxon Mann Whitney Test). For vectors containing an in-frame mononucleotide allele ([G/C](10)) or either of two dinucleotide alleles ([GT/CA](10) and [TC/AG](11)) we observed no statistically significant difference in the overall HSV-tk mutation frequency observed between dinB(+) and dinB(-) strains. To determine if a mutational bias exists for mutations made by Pol IV, mutational spectra were generated for each STR vector and strain. No statistically significant differences between strains were observed for either the proportion of mutational events at the STR or STR specificity among the three vectors. However, the specificity of mutational events at the STR alleles in each strain varied in a statistically significant manner as a consequence of microsatellite sequence. Our results indicate that while Pol IV contributes to spontaneous mutations within the HSV-tk coding sequence, Pol IV does not play a significant role in spontaneous mutagenesis at [G/C](10), [GT/CA](10), or [TC/AG](11) microsatellite alleles. Our data demonstrate that in a wild type genetic background, the major factor influencing microsatellite mutagenesis is the allelic sequence composition.

Microdeletions and microinsertions causing human genetic disease: common mechanisms of mutagenesis and the role of local DNA sequence complexity

In the Human Gene Mutation Database (www.hgmd.org), microdeletions and microinsertions causing inherited disease (both defined as involving < or = 20 bp of DNA) account for 8,399 (17%) and 3,345 (7%) logged mutations, in 940 and 668 genes, respectively. A positive correlation was noted between the microdeletion and microinsertion frequencies for 564 genes for which both microdeletions and microinsertions are reported in HGMD, consistent with the view that the propensity of a given gene/sequence to undergo microdeletion is related to its propensity to undergo microinsertion. While microdeletions and microinsertions of 1 bp constitute respectively 48% and 66% of the corresponding totals, the relative frequency of the remaining lesions correlates negatively with the length of the DNA sequence deleted or inserted. Many of the microdeletions and microinsertions of more than 1 bp are potentially explicable in terms of slippage mutagenesis, involving the addition or removal of one copy of a mono-, di-, or trinucleotide tandem repeat. The frequency of in-frame 3-bp and 6-bp microinsertions and microdeletions was, however, found to be significantly lower than that of mutations of other lengths, suggesting that some of these in-frame lesions may not have come to clinical attention. Various sequence motifs were found to be over-represented in the vicinity of both microinsertions and microdeletions, including the heptanucleotide CCCCCTG that shares homology with the complement of the 8-bp human minisatellite conserved sequence/chi-like element (GCWGGWGG). The previously reported indel hotspot GTAAGT and its complement ACTTAC were also found to be overrepresented in the vicinity of both microinsertions and microdeletions, thereby providing a first example of a mutational hotspot that is common to different types of gene lesion. Other motifs overrepresented in the vicinity of microdeletions and microinsertions included DNA polymerase pause sites and topoisomerase cleavage sites. Several novel microdeletion/microinsertion hotspots were noted and some of these exhibited sufficient similarity to one another to justify terming them "super-hotspot" motifs. Analysis of sequence complexity also demonstrated that a combination of slipped mispairing mediated by direct repeats, and secondary structure formation promoted by symmetric elements, can account for the majority of microdeletions and microinsertions. Thus, microinsertions and microdeletions exhibit strong similarities in terms of the characteristics of their flanking DNA sequences, implying that they are generated by very similar underlying mechanisms.

Evolution of hypervariable microsatellites in apomictic polyploid lineages of Ranunculus carpaticola: directional bias at dinucleotide loci

Microsatellites are widely used in genetic and evolutionary analyses, but their own evolution is far from simple. The mechanisms maintaining the mutational patterns of simple repeats and the typical stable allele-frequency distributions are still poorly understood. Asexual lineages may provide particularly informative models for the indirect study of microsatellite evolution, because their genomes act as complete linkage groups, with mutations being the only source of genetic variation. Here, we study the direction of accumulated dinucleotide microsatellite mutations in wild asexual lineages of hexaploid Ranunculus carpaticola. Whereas the overall number of contractions is not significantly different from that of expansions, the within-locus frequency of contractions, but not of expansions, significantly increases with allele length. Moreover, within-locus polymorphism is positively correlated with allele length, but this relationship is due solely to the influence of contraction mutations. Such asymmetries may explain length constraints generally observed with microsatellites and are consistent with stable, bell-shaped allele-frequency distributions. Although apomictic and allohexaploid, the R. carpaticola lineages show mutational patterns resembling the trends observed in a broad range of organisms, including sexuals and diploids, suggesting that, even if not of germline origin, the mutations in these apomicts may be the consequence of similar mechanisms.

Activating mutations in the tyrosine kinase domain of the epidermal growth factor receptor are associated with improved survival in gefitinib-treated chemorefractory lung adenocarcinomas

PURPOSE

Activating mutations in the tyrosine kinase domain of the epidermal growth factor receptor (EGFR) confer a strong sensitivity to gefitinib, a selective tyrosine kinase inhibitor of EGFR.

EXPERIMENTAL DESIGN

We examined EGFR mutations at exons 18, 19, and 21 in tumor tissue from 68 gefitinib-treated, chemorefractory, advanced non-small cell lung cancer patients from the United States, Europe, and Asia and in a highly gefitinib-sensitive non-small cell lung cancer cell line and correlated their presence with response and survival. In addition, in a subgroup of 28 patients for whom the remaining tumor tissue was available, we examined the relationship among EGFR mutations, CA repeats in intron 1 of EGFR, EGFR and caveolin-1 mRNA levels, and increased EGFR gene copy numbers.

RESULTS

Seventeen patients had EGFR mutations, all of which were in lung adenocarcinomas. Radiographic response was observed in 16 of 17 (94.1%) patients harboring EGFR mutations, in contrast with 6 of 51 (12.6%) with wild-type EGFR (P < 0.0001). Probability of response increased significantly in never smokers, patients receiving a greater number of prior chemotherapy regimens, Asians, and younger patients. Median survival was not reached for patients with EGFR mutations and was 9.9 months for those with wild-type EGFR (P = 0.001). EGFR mutations tended to be associated with increased numbers of CA repeats and increased EGFR gene copy numbers but not with EGFR and caveolin-1 mRNA overexpression (P = not significant).

CONCLUSIONS

The presence of EGFR mutations is a major determinant of gefitinib response, and targeting EGFR should be considered in preference to chemotherapy as first-line treatment in lung adenocarcinomas that have demonstrable EGFR mutations.

Hypermutation Rate Normalized by Chronological Time

It is generally believed that in cells undergoing Ig somatic hypermutation, more cell divisions result in more mutations. This is because DNA synthesis and replication is thought to play roles in the known mechanisms-cytidine deamination and subsequent conversion to thymidine, uracil-DNA glycosylase-mediated repair, mismatch repair, and DNA synthesis by error-prone polymerases. In this study, we manipulated the number of cell generations by varying the rate at which cultures of a mouse cell line were replenished with fresh medium. We found that the frequency of mutants does not necessarily increase with the number of cell generations. On the contrary, a greater number of divisions can lead to a lower frequency of mutants, indicating that cell division is not a rate-limiting step in the hypermutation process. Thus, when comparing mutation rates, we suggest that rates are more appropriately expressed as mutations per day than per cell generation.

Highly polymorphic repeat region in the CETP promoter induces unusual DNA structure

Genetic variation in the human cholesteryl ester transfer protein (CETP) promoter is associated with HDL cholesterol levels and cardiovascular disease with much of the genetic variation in CETP attributed to the promoter region. In this region, there are several single nucleotide polymorphisms as well as a variable length tandem repeat located 1946 base pairs upstream of the CETP transcription start that is highly polymorphic with respect to both length and sequence. There are more than 10 different long alleles and these vary in their repeat structure. We find that the short allele of this repeat is associated with high HDL cholesterol levels in vivo (P<0.0001). In males, this association is independent of the nearby -629 polymorphism. In addition, the variable length GAAA repeat can stimulate an adjacent GGGGA repeat to form a structure that hinders DNA amplification and sequencing. This structure also has an effect in vivo as shown by orientation effects and cloning efficiency in Escherichia coli.

Modification of breast cancer risk in young women by a polymorphic sequence in the egfr gene

The regulation of the epidermal growth factor receptor (egfr) gene in human cancer is not yet fully understood. Recent data on a polymorphic CA repeat located at the 5'-regulatory sequence in intron 1 of the egfr gene [egfr CA simple sequence repeat (SSR) I] point to a possible inheritance of cancer risk associated with the egfr gene. Furthermore, we have detected frequent allelic imbalances restricted to the egfr CA SSR I in breast cancer tissue and nontumorous breast tissue adjacent to invasive and in situ breast cancer representing amplifications. Therefore, we conducted a population-based case-control study to assess the relationship between the egfr polymorphism and breast cancer risk. Cases with a first primary breast cancer by age 50 years and age-matched population controls provided information on known and suspected risk factors. The allelic length of the egfr CA SSR was determined in 616 cases and 1072 population-sampled controls. Genotypes were categorized for analysis by allele length. Multivariate logistic regression was used to compare genotype distributions, accounting for other risk factors, and to investigate gene-environment interactions. We found a modifying effect, albeit no main effect, of the allelic length of the egfr polymorphism on breast cancer risk. The presence of two long alleles (>/==" BORDER="0">19 CA) was associated with a significantly elevated odds ratio (OR) of 10.4 [95% confidence interval (CI), 1.85-58.70] among women with a first-degree family history of breast cancer (P = 0.015 for interaction). The risk increase associated with high red meat consumption (OR, 10.68; 95% CI, 1.57-72.58) and the protective effect of high vegetable intake (OR, 0.07; 95% CI, 0.004-1.07) was also most pronounced among carriers of two long alleles (>/==" BORDER="0">19 CA). The length of the egfr CA SSR may increase the risk for familial breast cancers, and its effect could be modulated by dietary factors.

Y-chromosomal microsatellite mutation rates: differences in mutation rate between and within loci

Precise estimates of mutation rates at Y-chromosomal microsatellite STR (short tandem repeat) loci make an important basis for paternity diagnostics and dating of Y chromosome lineage origins. There are indications of considerable locus mutation rate variability between (inter-) and within (intra-) loci. We have studied nine Y-STR loci-DYS19, DYS389I/II, DYS390, DYS391, DYS392, DYS393, DYS385, and DYS388-in 1,766 father-son pairs of confirmed paternity (a total of 15,894 meioses). Five biallelic markers were also analyzed in the fathers-Tat, YAP, 12f2, SRY1532, and 92R7-defining haplogroups 1, 2, 3, 4, 9, and 16, respectively. A total of 36 fragment length mutations were observed: 24 gains (22 single-step, two double-step) and 12 single-step losses. Thus, there was a significant surplus of gains (p=0.045). Overall, the mutation rate was positively correlated to STR repeat length and there was a significant relative excess of losses in long alleles and gains in short alleles (p=0.043). In contrast to the situation in autosomal STR loci and in MSY-1, no noteworthy correlation between mutation rate and the father's age at the child's birth was observed. We observed significant interlocus differences in Y-STR mutation rates (p<0.01). The number of observed mutations ranged from zero in DYS392 to eight in DYS391 and DYS390. We have also demonstrated obvious differences in mutation rates between the haplogroups studied (p=0.024), a phenomenon that is a reflection of the dependence of mutation rate on allele size. Our study has thus demonstrated the necessity of not only locus-specific, but even allele-specific, mutation rate estimates for forensic and population genetic purposes, and provides a considerable basis for such estimates.

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

Microsatellite DNA sequences are ubiquitous in the human genome, and mutation rates of these repetitive sequences vary with respect to DNA sequence as well as length. We have analyzed polymerase-DNA interactions as a function of microsatellite sequence, using polypyrimidine/polypurine di- and tetranucleotide alleles representative of those found in the human genome. Using an in vitro primer extension assay and the mammalian DNA polymerase alpha-primase complex, we have observed a polymerase termination profile for each microsatellite that is unique to that allele. Interestingly, a periodic termination profile with an interval size (9-11 nucleotides) unrelated to microsatellite unit length was observed for the [TC](20) and [TTCC](9) templates. In contrast, a unit-punctuated polymerase termination profile was found for the longer polypurine templates. We detected strong polymerase pauses within the [TC](20) allele at low reaction pH which were eliminated by the addition of deaza-dGTP, consistent with these specific pauses being a consequence of triplex DNA formation during DNA synthesis. Quantitatively, a strand bias was observed in the primer extension assay, in that polymerase synthesis termination is more intense when the polypurine sequence serves as the template, relative to its complementary polypyrimidine sequence. The HSV-tk forward mutation assay was utilized to determine the corresponding polymerase alpha-primase error frequencies and specificities at the microsatellite alleles. A higher microsatellite polymerase error frequency (50x10(-4) to 60x10(-4)) was measured when polypurine sequences serve as templates for DNA synthesis, relative to the polypyrimidine template (18x10(-4)). Thus, a positive correlation exists between polymerase alpha-primase pausing and mutagenesis within microsatellite DNA alleles.