High-Throughput Amplicon-Based Copy Number Detection of 11 Genes in Formalin-Fixed Paraffin-Embedded Ovarian Tumour Samples by MLPA-Seq

Whilst next generation sequencing can report point mutations in fixed tissue tumour samples reliably, the accurate determination of copy number is more challenging. The conventional Multiplex Ligation-dependent Probe Amplification (MLPA) assay is an effective tool for measurement of gene dosage, but is restricted to around 50 targets due to size resolution of the MLPA probes. By switching from a size-resolved format, to a sequence-resolved format we developed a scalable, high-throughput, quantitative assay. MLPA-seq is capable of detecting deletions, duplications, and amplifications in as little as 5ng of genomic DNA, including from formalin-fixed paraffin-embedded (FFPE) tumour samples. We show that this method can detect BRCA1, BRCA2, ERBB2 and CCNE1 copy number changes in DNA extracted from snap-frozen and FFPE tumour tissue, with 100% sensitivity and >99.5% specificity.     

Extended RAS and BRAF Mutation Analysis Using Next-Generation Sequencing

Somatic mutations in KRAS, NRAS, and BRAF genes are related to resistance to anti-EGFR antibodies in colorectal cancer. We have established an extended RAS and BRAF mutation assay using a next-generation sequencer to analyze these mutations. Multiplexed deep sequencing was performed to detect somatic mutations within KRAS, NRAS, and BRAF, including minor mutated components. We first validated the technical performance of the multiplexed deep sequencing using 10 normal DNA and 20 formalin-fixed, paraffin-embedded (FFPE) tumor samples. To demonstrate the potential clinical utility of our assay, we profiled 100 FFPE tumor samples and 15 plasma samples obtained from colorectal cancer patients. We used a variant calling approach based on a Poisson distribution. The distribution of the mutation-positive population was hypothesized to follow a Poisson distribution, and a mutation-positive status was defined as a value greater than the significance level of the error rate (α = 2 x 10^-5). The cut-off value was determined to be the average error rate plus 7 standard deviations. Mutation analysis of 100 clinical FFPE tumor specimens was performed without any invalid cases. Mutations were detected at a frequency of 59% (59/100). KRAS mutation concordance between this assay and Scorpion-ARMS was 92% (92/100). DNA obtained from 15 plasma samples was also analyzed. KRAS and BRAF mutations were identified in both the plasma and tissue samples of 6 patients. The genetic screening assay using next-generation sequencer was validated for the detection of clinically relevant RAS and BRAF mutations using FFPE and liquid samples.     

Multi-Purpose Utility of Circulating Plasma DNA Testing in Patients with Advanced Cancers

Tumor genomic instability and selective treatment pressures result in clonal disease evolution; molecular stratification for molecularly targeted drug administration requires repeated access to tumor DNA. We hypothesized that circulating plasma DNA (cpDNA) in advanced cancer patients is largely derived from tumor, has prognostic utility, and can be utilized for multiplex tumor mutation sequencing when repeat biopsy is not feasible. We utilized the Sequenom MassArray System and OncoCarta panel for somatic mutation profiling. Matched samples, acquired from the same patient but at different time points were evaluated; these comprised formalin-fixed paraffin-embedded (FFPE) archival tumor tissue (primary and/or metastatic) and cpDNA. The feasibility, sensitivity, and specificity of this high-throughput, multiplex mutation detection approach was tested utilizing specimens acquired from 105 patients with solid tumors referred for participation in Phase I trials of molecularly targeted drugs. The median cpDNA concentration was 17 ng/ml (range: 0.5–1600); this was 3-fold higher than in healthy volunteers. Moreover, higher cpDNA concentrations associated with worse overall survival; there was an overall survival (OS) hazard ratio of 2.4 (95% CI 1.4, 4.2) for each 10-fold increase in cpDNA concentration and in multivariate analyses, cpDNA concentration, albumin, and performance status remained independent predictors of OS. These data suggest that plasma DNA in these cancer patients is largely derived from tumor. We also observed high detection concordance for critical ‘hot-spot’ mutations (KRAS, BRAF, PIK3CA) in matched cpDNA and archival tumor tissue, and important differences between archival tumor and cpDNA. This multiplex sequencing assay can be utilized to detect somatic mutations from plasma in advanced cancer patients, when safe repeat tumor biopsy is not feasible and genomic analysis of archival tumor is deemed insufficient. Overall, circulating nucleic acid biomarker studies have clinically important multi-purpose utility in advanced cancer patients and further studies to pursue their incorporation into the standard of care are warranted.     

Single-Tubed Wild-Type Blocking Quantitative PCR Detection Assay for the Sensitive Detection of Codon 12 and 13 KRAS Mutations

The high degree of intra-tumor heterogeneity has meant that it is important to develop sensitive and selective assays to detect low-abundance KRAS mutations in metastatic colorectal carcinoma (mCRC) patients. As a major potential source of tumor DNA in the aforementioned genotyping assays, it was necessary to conduct an analysis on both the quality and quantity of DNA extracted from formalin-fixed paraffin-embedded (FFPE). Therefore, four commercial FFPE DNA extraction kits were initially compared with respect to their ability to facilitate extraction of amplifiable DNA. The results showed that TrimGen kits showed the greatest performance in relation to the quality and quantity of extracted FFPE DNA solutions. Using DNA extracted by TrimGen kits as a template for tumor genotyping, a real-time wild-type blocking PCR (WTB-PCR) assay was subsequently developed to detect the aforementioned KRAS mutations in mCRC patients. The results showed that WTB-PCR facilitated the detection of mutated alleles at a ratio of 1:10,000 (i.e. 0.01%) wild-type alleles. When the assay was subsequently used to test 49 mCRC patients, the results showed that the mutation detection levels of the WTB-PCR assay (61.8%; 30/49) were significantly higher than that of traditional PCR (38.8%; 19/49). Following the use of the real-time WTB-PCR assay, the ΔC _q method was used to quantitatively analyze the mutation levels associated with KRAS in each FFPE sample. The results showed that the mutant levels ranged from 53.74 to 0.12% in the patients analyzed. In conclusion, the current real-time WTB-PCR is a rapid, simple, and low-cost method that permits the detection of trace amounts of the mutated KRAS gene.     

Profiling Critical Cancer Gene Mutations in Clinical Tumor Samples

Background

Detection of critical cancer gene mutations in clinical tumor specimens may predict patient outcomes and inform treatment options; however, high-throughput mutation profiling remains underdeveloped as a diagnostic approach. We report the implementation of a genotyping and validation algorithm that enables robust tumor mutation profiling in the clinical setting.

Methodology

We developed and implemented an optimized mutation profiling platform (“OncoMap”) to interrogate ∼400 mutations in 33 known oncogenes and tumor suppressors, many of which are known to predict response or resistance to targeted therapies. The performance of OncoMap was analyzed using DNA derived from both frozen and FFPE clinical material in a diverse set of cancer types. A subsequent in-depth analysis was conducted on histologically and clinically annotated pediatric gliomas. The sensitivity and specificity of OncoMap were 93.8% and 100% in fresh frozen tissue; and 89.3% and 99.4% in FFPE-derived DNA. We detected known mutations at the expected frequencies in common cancers, as well as novel mutations in adult and pediatric cancers that are likely to predict heightened response or resistance to existing or developmental cancer therapies. OncoMap profiles also support a new molecular stratification of pediatric low-grade gliomas based on BRAF mutations that may have immediate clinical impact.

Conclusions

Our results demonstrate the clinical feasibility of high-throughput mutation profiling to query a large panel of “actionable” cancer gene mutations. In the future, this type of approach may be incorporated into both cancer epidemiologic studies and clinical decision making to specify the use of many targeted anticancer agents.     

Host choice and human blood index of Anopheles pseudopunctipennis in a village of the Andean valleys of Bolivia

Background

Knockdown resistance (kdr) is a well-characterized mechanism of resistance to pyrethroid insecticides in many insect species and is caused by point mutations of the pyrethroid target site the para -type sodium channel. The presence of kdr mutations in Anopheles gambiae, the most important malaria vector in Africa, has been monitored using a variety of molecular techniques. However, there are few reports comparing the performance of these different assays. In this study, two new high-throughput assays were developed and compared with four established techniques.

Methods

Fluorescence-based assays based on 1) TaqMan probes and 2) high resolution melt (HRM) analysis were developed to detect kdr alleles in An. gambiae. Four previously reported techniques for kdr detection, Allele Specific Polymerase Chain Reaction (AS-PCR), Heated Oligonucleotide Ligation Assay (HOLA), Sequence Specific Oligonucleotide Probe – Enzyme-Linked ImmunoSorbent Assay (SSOP-ELISA) and PCR-Dot Blot were also optimized. The sensitivity and specificity of all six assays was then compared in a blind genotyping trial of 96 single insect samples that included a variety of kdr genotypes and African Anopheline species. The relative merits of each assay was assessed based on the performance in the genotyping trial, the length/difficulty of each protocol, cost (both capital outlay and consumable cost), and safety (requirement for hazardous chemicals).

Results

The real-time TaqMan assay was both the most sensitive (with the lowest number of failed reactions) and the most specific (with the lowest number of incorrect scores). Adapting the TaqMan assay to use a PCR machine and endpoint measurement with a fluorimeter showed a slight reduction in sensitivity and specificity. HRM initially gave promising results but was more sensitive to both DNA quality and quantity and consequently showed a higher rate of failure and incorrect scores. The sensitivity and specificity of AS-PCR, SSOP-ELISA, PCR Dot Blot and HOLA was fairly similar with a small number of failures and incorrect scores.

Conclusion

The results of blind genotyping trials of each assay indicate that where maximum sensitivity and specificity are required the TaqMan real-time assay is the preferred method. However, the cost of this assay, particularly in terms of initial capital outlay, is higher than that of some of the other methods. TaqMan assays using a PCR machine and fluorimeter are nearly as sensitive as real-time assays and provide a cost saving in capital expenditure. If price is a primary factor in assay choice then the AS-PCR, SSOP-ELISA, and HOLA are all reasonable alternatives with the SSOP-ELISA approach having the highest throughput.     

Abbott RealTime HBV Assay Is More Sensitive in Detection of Low Viral Load and Little Impacted by Drug Resistant Mutation in Chronic Hepatitis B Patients under Nucleot(s)ide Analogues Therapy

Background

Selection of drug-resistant strains may lead to failure of HBV antiviral therapy. There is little information whether there is detection difference in drug resistant mutations between different viral load assays of HBV.

Objectives

This study is aimed to investigate whether there is drug-resistant strains related detection difference between Abbott RealTime HBV (RealTime) and CobasAmpliPrep/CobasTaqMan HBV assays 2.0 (TaqMan).

Study Design

One hundred and thirty-four CHB patients who received HBV anti-viral therapy were enrolled. HBV virological markers were tested 3 months apart regularly. Serum HBV DNA levels were determined using the TaqMan and RealTime. YMDD (rt180M and rt204V) mutation was checked in patients who experienced virologic breakthrough (VBT).

Results

The correlation of HBV DNA observed between the RealTime and TaqMan was good for all 571 samples (R² = 0.797; P<0.001). However, the correlation in the 434 samples with HBV DNA level <3 log₁₀ IU/ml was not as good as in all samples (R² = 0.457). Overall, 21.5% of samples had a detection difference of ≥1 log₁₀ IU/ml with 91.9% of these having HBV DNA level <3 log₁₀ IU/ml. Twenty-four patients experiencedVBT. Three of these patients had acquired the YMDD mutation and exhibited discordant viral load results between the two methods tested. In each case, persistent HBV DNA was detected by RealTime and undetectable with TaqMan. Of the patients who experienced a VBT and had acquired YMDD mutation, 4.7% had undetectable HBV DNA by TaqMan while all were detectable with RealTime.

Conclusions

RealTime assay is more sensitive and is little impacted by the development of drug resistant mutation.     

Detection of Ralstonia solanacearum Strains with a Quantitative, Multiplex, Real-Time, Fluorogenic PCR (TaqMan) Assay

A fluorogenic (TaqMan) PCR assay was developed to detect Ralstonia solanacearum strains. Two fluorogenic probes were utilized in a multiplex reaction; one broad-range probe (RS) detected all biovars of R. solanacearum, and a second more specific probe (B2) detected only biovar 2A. Amplification of the target was measured by the 5′ nuclease activity of Taq DNA polymerase on each probe, resulting in emission of fluorescence. TaqMan PCR was performed with DNA extracted from 42 R. solanacearum and genetically or serologically related strains to demonstrate the specificity of the assay. In pure cultures, detection of R. solanacearum to ≥10² cells ml⁻¹ was achieved. Sensitivity decreased when TaqMan PCR was performed with inoculated potato tissue extracts, prepared by currently recommended extraction procedures. A third fluorogenic probe (COX), designed with the potato cytochrome oxidase gene sequence, was also developed for use as an internal PCR control and was shown to detect potato DNA in an RS-COX multiplex TaqMan PCR with infected potato tissue. The specificity and sensitivity of the assay, combined with high speed, robustness, reliability, and the possibility of automating the technique, offer potential advantages in routine indexing of potato tubers and other plant material for the presence of R. solanacearum.     

High throughput interrogation of somatic mutations in high grade serous cancer of the ovary

Background

Epithelial ovarian cancer is the most lethal of all gynecologic malignancies, and high grade serous ovarian cancer (HGSC) is the most common subtype of ovarian cancer. The objective of this study was to determine the frequency and types of point somatic mutations in HGSC using a mutation detection protocol called OncoMap that employs mass spectrometric-based genotyping technology.

Methodology/Principal Findings

The Center for Cancer Genome Discovery (CCGD) Program at the Dana-Farber Cancer Institute (DFCI) has adapted a high-throughput genotyping platform to determine the mutation status of a large panel of known cancer genes. The mutation detection protocol, termed OncoMap has been expanded to detect more than 1000 mutations in 112 oncogenes in formalin-fixed paraffin-embedded (FFPE) tissue samples. We performed OncoMap on a set of 203 FFPE advanced staged HGSC specimens. We isolated genomic DNA from these samples, and after a battery of quality assurance tests, ran each of these samples on the OncoMap v3 platform. 56% (113/203) tumor samples harbored candidate mutations. Sixty-five samples had single mutations (32%) while the remaining samples had ≥2 mutations (24%). 196 candidate mutation calls were made in 50 genes. The most common somatic oncogene mutations were found in EGFR, KRAS, PDGRFα, KIT, and PIK3CA. Other mutations found in additional genes were found at lower frequencies (<3%).

Conclusions/Significance

Sequenom analysis using OncoMap on DNA extracted from FFPE ovarian cancer samples is feasible and leads to the detection of potentially druggable mutations. Screening HGSC for somatic mutations in oncogenes may lead to additional therapies for this patient population.     

Mrassf1a-Pap,a Novel Methylation-Based Assay for the Detection of Cell-Free Fetal DNA in Maternal Plasma

Objectives

RASSF1A has been described to be differentially methylated between fetal and maternal DNA and can therefore be used as a universal sex-independent marker to confirm the presence of fetal sequences in maternal plasma. However, this requires highly sensitive methods. We have previously shown that Pyrophosphorolysis-activated Polymerization (PAP) is a highly sensitive technique that can be used in noninvasive prenatal diagnosis. In this study, we have used PAP in combination with bisulfite conversion to develop a new universal methylation-based assay for the detection of fetal methylated RASSF1A sequences in maternal plasma.

Methods

Bisulfite sequencing was performed on maternal genomic (g)DNA and fetal gDNA from chorionic villi to determine differentially methylated regions in the RASSF1A gene using bisulfite specific PCR primers. Methylation specific primers for PAP were designed for the detection of fetal methylated RASSF1A sequences after bisulfite conversion and validated.

Results

Serial dilutions of fetal gDNA in a background of maternal gDNA show a relative percentage of ∼3% can be detected using this assay. Furthermore, fetal methylated RASSF1A sequences were detected both retrospectively as well as prospectively in all maternal plasma samples tested (n = 71). No methylated RASSF1A specific bands were observed in corresponding maternal gDNA. Specificity was further determined by testing anonymized plasma from non-pregnant females (n = 24) and males (n = 21). Also, no methylated RASSF1A sequences were detected here, showing this assay is very specific for methylated fetal DNA. Combining all samples and controls, we obtain an overall sensitivity and specificity of 100% (95% CI 98.4%–100%).

Conclusions

Our data demonstrate that using a combination of bisulfite conversion and PAP fetal methylated RASSF1A sequences can be detected with extreme sensitivity in a universal and sex-independent manner. Therefore, this assay could be of great value as an addition to current techniques used in noninvasive prenatal diagnostics.     

An Integrative Analysis of PIK3CA Mutation,PTEN, and INPP4B Expression in Terms of Trastuzumab Efficacy in HER2-Positive Breast Cancer

The phosphoinositide-3-kinase (PI3K) pathway is commonly deregulated in breast cancer through several mechanisms, including PIK3CA mutation and loss of phosphatase and tensin homolog (PTEN) and inositol polyphosphate 4-phosphatase-II (INPP4B). We aimed to evaluate the predictive relevance of these biomarkers to trastuzumab efficacy in HER2-positive disease. We evaluated the effect of trastuzumab in 43 breast cancer patients with HER2-overexpression who received neoadjuvant treatment. PIK3CA mutation was examined by direct sequencing and digital PCR assay, and PIK3CA copy number was assessed by digital PCR assay of pretreatment tissues. PTEN, pAkt, and INPP4B were assessed by immunohistochemistry. Direct sequencing detected mutant DNA in 21% of all patients, but the incidence increased to 49% using digital PCR. The pathological complete response (pCR) rate in patients with PIK3CA mutations was 29% compared with 67% for those without PIK3CA mutations (P = 0.093), when the mutation was defined as positive if the mutant proportion was more than 10% of total genetic content by digital PCR. Low PTEN expression was associated with less pCR compared to high expression (33% versus 72%, P = 0.034). There were no significant associations of PIK3CA copy number, pAKt, or INPP4B with trastuzumab efficacy. In multivariate analysis, activation of the PI3K pathway due to either PIK3CA mutation or low PTEN were related to poorer response to trastuzumab (OR of predictive pCR was 0.11, 95%CI; 0.03–0.48). In conclusion, activating the PI3K pathway is associated with low pCR to trastuzumab-based treatment in HER2-positive breast cancer. Combined analysis of PIK3CA mutation and PTEN expression may serve as critical indicators to identify patients unlikely to respond to trastuzumab.     

Oxidative Stress Is Not a Major Contributor to Somatic Mitochondrial DNA Mutations

The accumulation of somatic mitochondrial DNA (mtDNA) mutations is implicated in aging and common diseases of the elderly, including cancer and neurodegenerative disease. However, the mechanisms that influence the frequency of somatic mtDNA mutations are poorly understood. To develop a simple invertebrate model system to address this matter, we used the Random Mutation Capture (RMC) assay to characterize the age-dependent frequency and distribution of mtDNA mutations in the fruit fly Drosophila melanogaster. Because oxidative stress is a major suspect in the age-dependent accumulation of somatic mtDNA mutations, we also used the RMC assay to explore the influence of oxidative stress on the somatic mtDNA mutation frequency. We found that many of the features associated with mtDNA mutations in vertebrates are conserved in Drosophila, including a comparable somatic mtDNA mutation frequency (∼10⁻⁵), an increased frequency of mtDNA mutations with age, and a prevalence of transition mutations. Only a small fraction of the mtDNA mutations detected in young or old animals were G∶C to T∶A transversions, a signature of oxidative damage, and loss-of-function mutations in the mitochondrial superoxide dismutase, Sod2, had no detectable influence on the somatic mtDNA mutation frequency. Moreover, a loss-of-function mutation in Ogg1, which encodes a DNA repair enzyme that removes oxidatively damaged deoxyguanosine residues (8-hydroxy-2′-deoxyguanosine), did not significantly influence the somatic mtDNA mutation frequency of Sod2 mutants. Together, these findings indicate that oxidative stress is not a major cause of somatic mtDNA mutations. Our data instead suggests that somatic mtDNA mutations arise primarily from errors that occur during mtDNA replication. Further studies using Drosophila should aid in the identification of factors that influence the frequency of somatic mtDNA mutations.     

Deletion of mitochondrial DNA bypassing a chromosomal gene needed for maintenance of the killer plasmid of yeast

Strains of Saccharomyces cerevisiae carrying a 1.4 x 10⁶ dalton double-stranded (ds) RNA in virus-like particles (the killer plasmid or virus) secrete a toxin that is lethal to strains not carrying this plasmid (virus). The mak10 gene is one of 24 chromosomal genes (called pets, mak1, mak2,...) that are needed to maintain and replicate the killer plasmid. We report here isolation of spontaneous and induced mutants in which the killer plasmid is maintained and replicated in spite of a defect in the mak10 gene. The bypass (or suppressor) mutations in these strains are in the mitochondrial genome. Respiratory deficiency produced by various chromosomal pet mutations, by chloramphenicol, or by antimycin A, does not bypass the mak10–1 mutation. Several spontaneous mak10–1 killer strains have about 12-fold more of the killer plasmid ds RNA than do wild-type killers. Although the absence of mitochondrial DNA bypasses mak10–1, it does not bypass pets–1, mak1–1, mak3–1, mak4–1, mak5–1, mak6–1, mak7–1, or mak8–1.     

Yeast Cycloheximide-resistant crl Mutants Are Proteasome Mutants Defective in Protein Degradation

In 1988 McCusker and Haber generated a series of mutants which are resistant to the minimum inhibitory concentration of the protein synthesis inhibitor cycloheximide. These cycloheximide-resistant, temperature-sensitive (crl) mutants, in addition, exhibited other pleiotropic phenotypes, e.g., incorrect response to starvation, hypersensitivity against amino acid analogues, and other protein synthesis inhibitors. Temperature sensitivity of one of these mutants, crl3–2, had been found to be suppressed by a mutation, SCL1–1, which resided in an α-type subunit of the 20S proteasome. We cloned the CRL3 gene by complementation and found CRL3 to be identical to the SUG1/CIM3 gene coding for a subunit of the 19S cap complex of the 26S proteasome. Another mutation, crl21, revealed to be allelic with the 20S proteasomal gene PRE3. crl3–2 and crl21 mutant cells show significant defects in proteasome-dependent proteolysis, whereas the SCL1–1 suppressor mutation causes partial restoration of crl3–2-induced proteolytic defects. Notably, cycloheximide resistance was also detected for other proteolytically deficient proteasome mutants (pre1–1, pre2–1, pre3–1, pre4–1). Moreover, proteasomal genes were found within genomic sequences of 9 of 13 chromosomal loci to which crl mutations had been mapped. We therefore assume that most if not all crl mutations reside in the proteasome and that phenotypes found are a result of defective protein degradation.     

Mutation enrichment in human DNA samples via UV-mediated cross-linking

Detection of low-level DNA mutations can reveal recurrent, hotspot genetic changes of clinical relevance to cancer, prenatal diagnostics, organ transplantation or infectious diseases. However, the high excess of wild-type (WT) alleles, which are concurrently present, often hinders identification of salient genetic changes. Here, we introduce UV-mediated cross-linking minor allele enrichment (UVME), a novel approach that incorporates ultraviolet irradiation (∼365 nm UV) DNA cross-linking either before or during PCR amplification. Oligonucleotide probes matching the WT target sequence and incorporating a UV-sensitive 3-cyanovinylcarbazole nucleoside modification are employed for cross-linking WT DNA. Mismatches formed with mutated alleles reduce DNA binding and UV-mediated cross-linking and favor mutated DNA amplification. UV can be applied before PCR and/or at any stage during PCR to selectively block WT DNA amplification and enable identification of traces of mutated alleles. This enables a single-tube PCR reaction directly from genomic DNA combining optimal pre-amplification of mutated alleles, which then switches to UV-mediated mutation enrichment-based DNA target amplification. UVME cross-linking enables enrichment of mutated KRAS and p53 alleles, which can be screened directly via Sanger sequencing, high-resolution melting, TaqMan genotyping or digital PCR, resulting in the detection of mutation allelic frequencies of 0.001–0.1% depending on the endpoint detection method. UV-mediated mutation enrichment provides new potential for mutation enrichment in diverse clinical samples.     

Whole Genome Sequencing of Mutation Accumulation Lines Reveals a Low Mutation Rate in the Social Amoeba Dictyostelium discoideum

Spontaneous mutations play a central role in evolution. Despite their importance, mutation rates are some of the most elusive parameters to measure in evolutionary biology. The combination of mutation accumulation (MA) experiments and whole-genome sequencing now makes it possible to estimate mutation rates by directly observing new mutations at the molecular level across the whole genome. We performed an MA experiment with the social amoeba Dictyostelium discoideum and sequenced the genomes of three randomly chosen lines using high-throughput sequencing to estimate the spontaneous mutation rate in this model organism. The mitochondrial mutation rate of 6.76×10⁻⁹, with a Poisson confidence interval of 4.1×10⁻⁹ − 9.5×10⁻⁹, per nucleotide per generation is slightly lower than estimates for other taxa. The mutation rate estimate for the nuclear DNA of 2.9×10⁻¹¹, with a Poisson confidence interval ranging from 7.4×10⁻¹³ to 1.6×10⁻¹⁰, is the lowest reported for any eukaryote. These results are consistent with low microsatellite mutation rates previously observed in D. discoideum and low levels of genetic variation observed in wild D. discoideum populations. In addition, D. discoideum has been shown to be quite resistant to DNA damage, which suggests an efficient DNA-repair mechanism that could be an adaptation to life in soil and frequent exposure to intracellular and extracellular mutagenic compounds. The social aspect of the life cycle of D. discoideum and a large portion of the genome under relaxed selection during vegetative growth could also select for a low mutation rate. This hypothesis is supported by a significantly lower mutation rate per cell division in multicellular eukaryotes compared with unicellular eukaryotes.     

Mutagenic mechanisms of cancer-associated DNA polymerase ϵ alleles

A single amino acid residue change in the exonuclease domain of human DNA polymerase ϵ, P286R, is associated with the development of colorectal cancers, and has been shown to impart a mutator phenotype. The corresponding Pol ϵ allele in the yeast Saccharomyces cerevisiae (pol2-P301R), was found to drive greater mutagenesis than an entirely exonuclease-deficient Pol ϵ (pol2–4), an unexpected phenotype of ultra-mutagenesis. By studying the impact on mutation frequency, type, replication-strand bias, and sequence context, we show that ultra-mutagenesis is commonly observed in yeast cells carrying a range of cancer-associated Pol ϵ exonuclease domain alleles. Similarities between mutations generated by these alleles and those generated in pol2–4 cells indicate a shared mechanism of mutagenesis that yields a mutation pattern similar to cancer Signature 14. Comparison of POL2 ultra-mutator with pol2-M644G, a mutant in the polymerase domain decreasing Pol ϵ fidelity, revealed unexpected analogies in the sequence context and strand bias of mutations. Analysis of mutational patterns unique to exonuclease domain mutant cells suggests that backtracking of the polymerase, when the mismatched primer end cannot be accommodated in the proofreading domain, results in the observed insertions and T>A mutations in specific sequence contexts.