Estimation of the Spontaneous Mutation Rate per Nucleotide Site in a Drosophila melanogaster Full-Sib Family

We employed deep genome sequencing of two parents and 12 of their offspring to estimate the mutation rate per site per generation in a full-sib family of Drosophila melanogaster recently sampled from a natural population. Sites that were homozygous for the same allele in the parents and heterozygous in one or more offspring were categorized as candidate mutations and subjected to detailed analysis. In 1.23 × 10⁹ callable sites from 12 individuals, we confirmed six single nucleotide mutations. We estimated the false negative rate in the experiment by generating synthetic mutations using the empirical distributions of numbers of nonreference bases at heterozygous sites in the offspring. The proportion of synthetic mutations at callable sites that we failed to detect was <1%, implying that the false negative rate was extremely low. Our estimate of the point mutation rate is 2.8 × 10⁻⁹ (95% confidence interval = 1.0 × 10⁻⁹ − 6.1 × 10⁻⁹) per site per generation, which is at the low end of the range of previous estimates, and suggests an effective population size for the species of ∼1.4 × 10⁶. At one site, point mutations were present in two individuals, indicating that there had been a premeiotic mutation cluster, although surprisingly one individual had a G→A transition and the other a G→T transversion, possibly associated with error-prone mismatch repair. We also detected three short deletion mutations and no insertions, giving a deletion mutation rate of 1.2 × 10⁻⁹ (95% confidence interval = 0.7 × 10⁻⁹ − 11 × 10⁻⁹).     

Mitochondrial DNA Haplogroups M7b1′2 and M8a Affect Clinical Expression of Leber Hereditary Optic Neuropathy in Chinese Families with the m.11778G→A Mutation

Leber hereditary optic neuropathy (LHON) is the most extensively studied mitochondrial disease, with the majority of the cases being caused by one of three primary mitochondrial DNA (mtDNA) mutations. Incomplete disease penetrance and gender bias are two features of LHON and indicate involvement of additional genetic or environmental factors in the pathogenesis of the disorder. Haplogroups J, K, and H have been shown to influence the clinical expression of LHON in subjects harboring primary mutations in European families. However, whether mtDNA haplogroups would affect the penetrance of LHON in East Asian families has not been evaluated yet. By studying the penetrance of LHON in 1859 individuals from 182 Chinese families (including one from Cambodia) with the m.11778G→A mutation, we found that haplogroup M7b1′2 significantly increases the risk of visual loss, whereas M8a has a protective effect. Analyses of the complete mtDNA sequences from LHON families with m.11778G→A narrow the association of disease expression to m.12811T→C (Y159H) in the NADH dehydrogenase 5 gene (MT-ND5) in haplogroup M7b1′2 and suggest that the specific combination of amino acid changes (A20T-T53I) in the ATP synthase 6 protein (MT-ATP6) caused by m.8584G→A and m.8684C→T might account for the beneficial background effect of M8a. Protein secondary-structure prediction for the MT-ATP6 with the two M8a-specific amino acid changes further supported our inferences. These findings will assist in further understanding the pathogenesis of LHON and guide future genetic counseling in East Asian patients with m.11778G→A.     

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.     

Ultra-Sensitive Sequencing Reveals an Age-Related Increase in Somatic Mitochondrial Mutations That Are Inconsistent with Oxidative Damage

Mitochondrial DNA (mtDNA) is believed to be highly vulnerable to age-associated damage and mutagenesis by reactive oxygen species (ROS). However, somatic mtDNA mutations have historically been difficult to study because of technical limitations in accurately quantifying rare mtDNA mutations. We have applied the highly sensitive Duplex Sequencing methodology, which can detect a single mutation among >10⁷ wild type molecules, to sequence mtDNA purified from human brain tissue from both young and old individuals with unprecedented accuracy. We find that the frequency of point mutations increases ∼5-fold over the course of 80 years of life. Overall, the mutation spectra of both groups are comprised predominantly of transition mutations, consistent with misincorporation by DNA polymerase γ or deamination of cytidine and adenosine as the primary mutagenic events in mtDNA. Surprisingly, G→T mutations, considered the hallmark of oxidative damage to DNA, do not significantly increase with age. We observe a non-uniform, age-independent distribution of mutations in mtDNA, with the D-loop exhibiting a significantly higher mutation frequency than the rest of the genome. The coding regions, but not the D-loop, exhibit a pronounced asymmetric accumulation of mutations between the two strands, with G→A and T→C mutations occurring more often on the light strand than the heavy strand. The patterns and biases we observe in our data closely mirror the mutational spectrum which has been reported in studies of human populations and closely related species. Overall our results argue against oxidative damage being a major driver of aging and suggest that replication errors by DNA polymerase γ and/or spontaneous base hydrolysis are responsible for the bulk of accumulating point mutations in mtDNA.     

Pathogenic mitochondrial DNA mutations are common in the general population

Mitochondrial DNA (mtDNA) mutations are a major cause of genetic disease, but their prevalence in the general population is not known. We determined the frequency of ten mitochondrial point mutations in 3168 neonatal-cord-blood samples from sequential live births, analyzing matched maternal-blood samples to estimate the de novo mutation rate. mtDNA mutations were detected in 15 offspring (0.54%, 95% CI = 0.30–0.89%). Of these live births, 0.00107% (95% CI = 0.00087–0.0127) harbored a mutation not detected in the mother's blood, providing an estimate of the de novo mutation rate. The most common mutation was m.3243A→G. m.14484T→C was only found on sub-branches of mtDNA haplogroup J. In conclusion, at least one in 200 healthy humans harbors a pathogenic mtDNA mutation that potentially causes disease in the offspring of female carriers. The exclusive detection of m.14484T→C on haplogroup J implicates the background mtDNA haplotype in mutagenesis. These findings emphasize the importance of developing new approaches to prevent transmission.     

Rates and Genomic Consequences of Spontaneous Mutational Events in Drosophila melanogaster

Because spontaneous mutation is the source of all genetic diversity, measuring mutation rates can reveal how natural selection drives patterns of variation within and between species. We sequenced eight genomes produced by a mutation-accumulation experiment in Drosophila melanogaster. Our analysis reveals that point mutation and small indel rates vary significantly between the two different genetic backgrounds examined. We also find evidence that ∼2% of mutational events affect multiple closely spaced nucleotides. Unlike previous similar experiments, we were able to estimate genome-wide rates of large deletions and tandem duplications. These results suggest that, at least in inbred lines like those examined here, mutational pressures may result in net growth rather than contraction of the Drosophila genome. By comparing our mutation rate estimates to polymorphism data, we are able to estimate the fraction of new mutations that are eliminated by purifying selection. These results suggest that ∼99% of duplications and deletions are deleterious—making them 10 times more likely to be removed by selection than nonsynonymous mutations. Our results illuminate not only the rates of new small- and large-scale mutations, but also the selective forces that they encounter once they arise.     

Mutagenic effects of 2-hydroxy-dATP on replication in a HeLa extract: induction of substitution and deletion mutations

The mutagenicity of an oxidized form of dATP, 2-hydroxydeoxyadenosine 5′-triphosphate (2-OH-dATP), was examined using an SV40 origin-dependent in vitro replication system with a HeLa extract. 2-OH-dATP induced mutations in a dose-dependent manner and elicited substitution and deletion mutations. Of the substitutions, a G·C→A·T transition including a tandem (CC→TT) mutation was mainly observed. This result agrees with our previous observation that mammalian DNA polymerase α misincorporates the oxidized nucleotide opposite C, but is in contrast to the finding that 2-OH-dATP elicits G·C→T·A transversions in Escherichia coli. This type of mutation was also elicited, but to a lesser extent. Interestingly, the mutagenicity of 2-OH-dATP was enhanced in the presence of 2-hydroxydeoxyadenosine 5′-diphosphate, an inhibitor of the MTH1 protein, suggesting that this protein functions in the hydrolysis of 2-OH-dATP in the replication reaction mixture, and probably in living cells. These results indicate that 2-OH-dATP is mutagenic and that its mutagenicity is suppressed by the MTH1 protein in mammalian cells.     

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.     

A Highlights from MBoC Selection: Genetic mosaic analysis of a deleterious mitochondrial DNA mutation in Drosophila reveals novel aspects of mitochondrial regulation and function

Various human diseases are associated with mitochondrial DNA (mtDNA) mutations, but heteroplasmy—the coexistence of mutant and wild-type mtDNA—complicates their study. We previously isolated a temperature-lethal mtDNA mutation in Drosophila, mt:CoI^T300I, which affects the cytochrome c oxidase subunit I (CoI) locus. In the present study, we found that the decrease in cytochrome c oxidase (COX) activity was ascribable to a temperature-dependent destabilization of cytochrome a heme. Consistently, the viability of homoplasmic flies at 29°C was fully restored by expressing an alternative oxidase, which specifically bypasses the cytochrome chains. Heteroplasmic flies are fully viable and were used to explore the age-related and tissue-specific phenotypes of mt:CoI^T300I. The proportion of mt:CoI^T300I genome remained constant in somatic tissues along the aging process, suggesting a lack of quality control mechanism to remove defective mitochondria containing a deleterious mtDNA mutation. Using a genetic scheme that expresses a mitochondrially targeted restriction enzyme to induce tissue-specific homoplasmy in heteroplasmic flies, we found that mt:CoI^T300I homoplasmy in the eye caused severe neurodegeneration at 29°C. Degeneration was suppressed by improving mitochondrial Ca²⁺ uptake, suggesting that Ca²⁺ mishandling contributed to mt:CoI^T300I pathogenesis. Our results demonstrate a novel approach for Drosophila mtDNA genetics and its application in modeling mtDNA diseases.     

Mitogenomic analysis of a 50-generation chicken pedigree reveals a rapid rate of mitochondrial evolution and evidence for paternal mtDNA inheritance

Mitochondrial genomes represent a valuable source of data for evolutionary research, but studies of their short-term evolution have typically been limited to invertebrates, humans and laboratory organisms. Here we present a detailed study of 12 mitochondrial genomes that span a total of 385 transmissions in a well-documented 50-generation pedigree in which two lineages of chickens were selected for low and high juvenile body weight. These data allowed us to test the hypothesis of time-dependent evolutionary rates and the assumption of strict maternal mitochondrial transmission, and to investigate the role of mitochondrial mutations in determining phenotype. The identification of a non-synonymous mutation in ND4L and a synonymous mutation in CYTB, both novel mutations in Gallus, allowed us to estimate a molecular rate of 3.13 × 10⁻⁷ mutations/site/year (95% confidence interval 3.75 × 10⁻⁸–1.12 × 10⁻⁶). This is substantially higher than avian rate estimates based upon fossil calibrations. Ascertaining which of the two novel mutations was present in an additional 49 individuals also revealed an instance of paternal inheritance of mtDNA. Lastly, an association analysis demonstrated that neither of the point mutations was strongly associated with the phenotypic differences between the two selection lines. Together, these observations reveal the highly dynamic nature of mitochondrial evolution over short time periods.     

The role of selection in the evolution of human mitochondrial genomes

High mutation rate in mammalian mitochondrial DNA generates a highly divergent pool of alleles even within species that have dispersed and expanded in size recently. Phylogenetic analysis of 277 human mitochondrial genomes revealed a significant (P < 0.01) excess of rRNA and nonsynonymous base substitutions among hotspots of recurrent mutation. Most hotspots involved transitions from guanine to adenine that, with thymine-to-cytosine transitions, illustrate the asymmetric bias in codon usage at synonymous sites on the heavy-strand DNA. The mitochondrion-encoded tRNAThr varied significantly more than any other tRNA gene. Threonine and valine codons were involved in 259 of the 414 amino acid replacements observed. The ratio of nonsynonymous changes from and to threonine and valine differed significantly (P = 0.003) between populations with neutral (22/58) and populations with significantly negative Tajima's D values (70/76), independent of their geographic location. In contrast to a recent suggestion that the excess of nonsilent mutations is characteristic of Arctic populations, implying their role in cold adaptation, we demonstrate that the surplus of nonsynonymous mutations is a general feature of the young branches of the phylogenetic tree, affecting also those that are found only in Africa. We introduce a new calibration method of the mutation rate of synonymous transitions to estimate the coalescent times of mtDNA haplogroups.     

Calibrating the Human Mutation Rate via Ancestral Recombination Density in Diploid Genomes

The human mutation rate is an essential parameter for studying the evolution of our species, interpreting present-day genetic variation, and understanding the incidence of genetic disease. Nevertheless, our current estimates of the rate are uncertain. Most notably, recent approaches based on counting de novo mutations in family pedigrees have yielded significantly smaller values than classical methods based on sequence divergence. Here, we propose a new method that uses the fine-scale human recombination map to calibrate the rate of accumulation of mutations. By comparing local heterozygosity levels in diploid genomes to the genetic distance scale over which these levels change, we are able to estimate a long-term mutation rate averaged over hundreds or thousands of generations. We infer a rate of 1.61 ± 0.13 × 10⁻⁸ mutations per base per generation, which falls in between phylogenetic and pedigree-based estimates, and we suggest possible mechanisms to reconcile our estimate with previous studies. Our results support intermediate-age divergences among human populations and between humans and other great apes.     

Rare missense and synonymous variants in UBE1 are associated with X-linked infantile spinal muscular atrophy

X-linked infantile spinal muscular atrophy (XL-SMA) is an X-linked disorder presenting with the clinical features hypotonia, areflexia, and multiple congenital contractures (arthrogryposis) associated with loss of anterior horn cells and infantile death. To identify the XL-SMA disease gene, we performed large-scale mutation analysis in genes located between markers DXS8080 and DXS7132 (Xp11.3–Xq11.1). This resulted in detection of three rare novel variants in exon 15 of UBE1 that segregate with disease: two missense mutations (c.1617 G→T, p.Met539Ile; c.1639 A→G, p.Ser547Gly) present each in one XL-SMA family, and one synonymous C→T substitution (c.1731 C→T, p.Asn577Asn) identified in another three unrelated families. Absence of the missense mutations was demonstrated for 3550 and absence of the synonymous mutation was shown in 7914 control X chromosomes; therefore, these results yielded statistical significant evidence for the association of the synonymous substitution and the two missense mutations with XL-SMA (p = 2.416 × 10⁻¹⁰, p = 0.001815). We also demonstrated that the synonymous C→T substitution leads to significant reduction of UBE1 expression and alters the methylation pattern of exon 15, implying a plausible role of this DNA element in developmental UBE1 expression in humans. Our observations indicate first that XL-SMA is part of a growing list of neurodegenerative disorders associated with defects in the ubiquitin-proteasome pathway and second that synonymous C→T transitions might have the potential to affect gene expression.     

A Heterozygous Truncating Mutation in RRM2B Causes Autosomal-Dominant Progressive External Ophthalmoplegia with Multiple mtDNA Deletions

Autosomal-dominant progressive external ophthalmoplegia (adPEO) is a mitochondrial disorder that is characterized by accumulation of multiple mitochondrial DNA (mtDNA) deletions in postmitotic tissues. The disorder is heterogeneous, with five known nuclear disease genes that encode the proteins ANT1, Twinkle, POLG, POLG2, and OPA1. Defects in these proteins affect mtDNA maintenance, probably leading to stalled replication forks, consequent mtDNA deletion formation, and progressive respiratory chain deficiency. Here we present a large adPEO family with multiple mtDNA deletions, whose disease was not explained by mutations in any of the known adPEO loci. We mapped the disease locus in this family to chromosome 8q22.1-q23.3. The critical linkage region contained the RRM2B gene, which encodes the small subunit of the ribonucleotide reductase p53R2, which has previously been shown to be essential for the maintenance of mtDNA copy number. Mutation screening of RRM2B revealed a heterozygous nonsense mutation in exon 9 (c.979C→T [p.R327X]) in all affected individuals that was absent in 380 control chromosomes. The same mutation was found to segregate in another adPEO family. The mutant mRNA escaped nonsense-mediated decay and resulted in a protein with truncation of 25 highly conserved C-terminal amino acids essential for the interaction with the ribonucleotide reductase subunit R1. We conclude that dominant-negative or gain-of-function mutations in RRM2B are a cause of multiple mtDNA deletions and adPEO.     

Detection of Ultra-Rare Mitochondrial Mutations in Breast Stem Cells by Duplex Sequencing

Long-lived adult stem cells could accumulate non-repaired DNA damage or mutations that increase the risk of tumor formation. To date, studies on mutations in stem cells have concentrated on clonal (homoplasmic) mutations and have not focused on rarely occurring stochastic mutations that may accumulate during stem cell dormancy. A major challenge in investigating these rare mutations is that conventional next generation sequencing (NGS) methods have high error rates. We have established a new method termed Duplex Sequencing (DS), which detects mutations with unprecedented accuracy. We present a comprehensive analysis of mitochondrial DNA mutations in human breast normal stem cells and non-stem cells using DS. The vast majority of mutations occur at low frequency and are not detectable by NGS. The most prevalent point mutation types are the C>T/G>A and A>G/T>C transitions. The mutations exhibit a strand bias with higher prevalence of G>A, T>C, and A>C mutations on the light strand of the mitochondrial genome. The overall rare mutation frequency is significantly lower in stem cells than in the corresponding non-stem cells. We have identified common and unique non-homoplasmic mutations between non-stem and stem cells that include new mutations which have not been reported previously. Four mutations found within the MT-ND5 gene (m.12684G>A, m.12705C>T, m.13095T>C, m.13105A>G) are present in all groups of stem and non-stem cells. Two mutations (m.8567T>C, m.10547C>G) are found only in non-stem cells. This first genome-wide analysis of mitochondrial DNA mutations may aid in characterizing human breast normal epithelial cells and serve as a reference for cancer stem cell mutation profiles.     

A novel mutator of Escherichia coli carrying a defect in the dgt gene, encoding a dGTP triphosphohydrolase

A novel mutator locus in Escherichia coli was identified from a collection of random transposon insertion mutants. Several mutators in this collection were found to have an insertion in the dgt gene, encoding a previously characterized dGTP triphosphohydrolase. The mutator activity of the dgt mutants displays an unusual specificity. Among the six possible base pair substitutions in a lacZ reversion system, the G·C→C·G transversion and A·T→G·C transition are strongly enhanced (10- to 50-fold), while a modest effect (two- to threefold) is also observed for the G·C→A·T transition. Interestingly, a two- to threefold reduction in mutant frequency (antimutator effect) is observed for the G·C→T·A transversion. In the absence of DNA mismatch repair (mutL) some of these effects are reduced or abolished, while other effects remain unchanged. Analysis of these effects, combined with the DNA sequence contexts in which the reversions take place, suggests that alterations of the dGTP pools as well as alterations in the level of some modified dNTP derivatives could affect the fidelity of in vivo DNA replication and, hence, account for the overall mutator effects.     

High mutation rates in the mitochondrial genomes of Daphnia pulex

Despite the great utility of mitochondrial DNA (mtDNA) sequence data in population genetics and phylogenetics, key parameters describing the process of mitochondrial mutation (e.g., the rate and spectrum of mutational change) are based on few direct estimates. Furthermore, the variation in the mtDNA mutation process within species or between lineages with contrasting reproductive strategies remains poorly understood. In this study, we directly estimate the mtDNA mutation rate and spectrum using Daphnia pulex mutation-accumulation (MA) lines derived from sexual (cyclically parthenogenetic) and asexual (obligately parthenogenetic) lineages. The nearly complete mitochondrial genome sequences of 82 sexual and 47 asexual MA lines reveal high mtDNA mutation rate of 1.37 × 10(-7) and 1.73 × 10(-7) per nucleotide per generation, respectively. The Daphnia mtDNA mutation rate is among the highest in eukaryotes, and its spectrum is dominated by insertions and deletions (70%), largely due to the presence of mutational hotspots at homopolymeric nucleotide stretches. Maximum likelihood estimates of the Daphnia mitochondrial effective population size reveal that between five and ten copies of mitochondrial genomes are transmitted per female per generation. Comparison between sexual and asexual lineages reveals no statistically different mutation rates and highly similar mutation spectra.     

Private Mitochondrial DNA Variants in Danish Patients with Hypertrophic Cardiomyopathy

Hypertrophic cardiomyopathy (HCM) is a genetic cardiac disease primarily caused by mutations in genes coding for sarcomeric proteins. A molecular-genetic etiology can be established in ~60% of cases. Evolutionarily conserved mitochondrial DNA (mtDNA) haplogroups are susceptibility factors for HCM. Several polymorphic mtDNA variants are associated with a variety of late-onset degenerative diseases and affect mitochondrial function. We examined the role of private, non-haplogroup associated, mitochondrial variants in the etiology of HCM. In 87 Danish HCM patients, full mtDNA sequencing revealed 446 variants. After elimination of 312 (69.9%) non-coding and synonymous variants, a further 109 (24.4%) with a global prevalence > 0.1%, three (0.7%) haplogroup associated and 19 (2.0%) variants with a low predicted in silico likelihood of pathogenicity, three variants: MT-TC: m.5772G>A, MT-TF: m.644A>G, and MT-CYB: m.15024G>A, p.C93Y remained. A detailed analysis of these variants indicated that none of them are likely to cause HCM. In conclusion, private mtDNA mutations are frequent, but they are rarely, if ever, associated with HCM.