Phylogeographic analysis of mitochondrial DNA haplogroup F2 in China reveals T12338C in the initiation codon of the ND5 gene not to be pathogenic

Toward a mtDNA locus-specific mutation database using the LOVD platform

The Human Variome Project (HVP) is a global effort to collect and curate all human genetic variation affecting health. Mutations of mitochondrial DNA (mtDNA) are an important cause of neurogenetic disease in humans; however, identification of the pathogenic mutations responsible can be problematic. In this article, we provide explanations as to why and suggest how such difficulties might be overcome. We put forward a case in support of a new Locus Specific Mutation Database (LSDB) implemented using the Leiden Open-source Variation Database (LOVD) system that will not only list primary mutations, but also present the evidence supporting their role in disease. Critically, we feel that this new database should have the capacity to store information on the observed phenotypes alongside the genetic variation, thereby facilitating our understanding of the complex and variable presentation of mtDNA disease. LOVD supports fast queries of both seen and hidden data and allows storage of sequence variants from high-throughput sequence analysis. The LOVD platform will allow construction of a secure mtDNA database; one that can fully utilize currently available data, as well as that being generated by high-throughput sequencing, to link genotype with phenotype enhancing our understanding of mitochondrial disease, with a view to providing better prognostic information.

Median network analysis of defectively sequenced entire mitochondrial genomes from early and contemporary disease studies

Sequence analysis of the mitochondrial genome has become a routine method in the study of mitochondrial diseases. Quite often, the sequencing efforts in the search of pathogenic or disease-associated mutations are affected by technical and interpretive problems, caused by sample mix-up, contamination, biochemical problems, incomplete sequencing, misdocumentation and insufficient reference to previously published data. To assess data quality in case studies of mitochondrial diseases, it is recommended to compare any mtDNA sequence under consideration to their phylogenetically closest lineages available in the Web. The median network method has proven useful for visualizing potential problems with the data. We contrast some early reports of complete mtDNA sequences to more recent total mtDNA sequencing efforts in studies of various mitochondrial diseases. We conclude that the quality of complete mtDNA sequences generated in the medical field in the past few years is somewhat unsatisfactory and may even fall behind that of pioneer manual sequencing in the early nineties. Our study provides a paradigm for an a posteriori evaluation of sequence quality and for detection of potential problems with inferring a pathogenic status of a particular mutation.

Distilling artificial recombinants from large sets of complete mtDNA genomes

Background

Large-scale genome sequencing poses enormous problems to the logistics of laboratory work and data handling. When numerous fragments of different genomes are PCR amplified and sequenced in a laboratory, there is a high immanent risk of sample confusion. For genetic markers, such as mitochondrial DNA (mtDNA), which are free of natural recombination, single instances of sample mix-up involving different branches of the mtDNA phylogeny would give rise to reticulate patterns and should therefore be detectable.

Methodology/Principal Findings

We have developed a strategy for comparing new complete mtDNA genomes, one by one, to a current skeleton of the worldwide mtDNA phylogeny. The mutations distinguishing the reference sequence from a putative recombinant sequence can then be allocated to two or more different branches of this phylogenetic skeleton. Thus, one would search for two (or three) near-matches in the total mtDNA database that together best explain the variation seen in the recombinants. The evolutionary pathway from the mtDNA tree connecting this pair together with the recombinant then generate a grid-like median network, from which one can read off the exchanged segments.

Conclusions

We have applied this procedure to a large collection of complete human mtDNA sequences, where several recombinants could be distilled by our method. All these recombinant sequences were subsequently corrected by de novo experiments – fully concordant with the predictions from our data-analytical approach.

Coding region mitochondrial DNA SNPs: targeting East Asian and Native American haplogroups

We have developed a single PCR multiplex SNaPshot reaction that consists of 32 coding region SNPs that allows (i) increasing the discrimination power of the mitochondrial DNA (mtDNA) typing in forensic casework, and (ii) haplogroup assignments of mtDNA profiles in both human population studies (e.g. anthropological) and medical research. The selected SNPs target the East Asian phylogeny, including its Native American derived branches. We have validated this multiplex assay by genotyping a sample of East Asians (Taiwanese) and Native Americans (Argentineans). In addition to the coding SNP typing, we have sequenced the complete control region for the same samples. The genotyping results (control region plus SNaPshot profiles) are in good agreement with previous human population genetic studies (based on e.g. complete sequencing) and the known mtDNA phylogeny. We observe that the SNaPshot method is reliable, rapid, and cost effective in comparison with other techniques of multiplex SNP genotyping. We discuss the advantages of our SNP genotyping selection with respect to previous attempts, and we highlight the importance of using the known mtDNA phylogeny as a framework for SNP profile interpretation and as a tool to minimize genotyping errors.

What is a 'novel' mtDNA mutation--and does 'novelty' really matter?

The hunt for pathogenic mitochondrial DNA (mtDNA) mutations is often fueled by the seeming novelty of mutations that are either nonsynonymous or affect the protein synthesis machinery in patients. In order to determine the novelty of a detected mutation, the working geneticist nearly always consults MITOMAP--often exclusively. By reanalyzing some case studies of refractory anemia with ring sideroblasts, prostate cancer, and hearing impairment, we demonstrate that the practice of solely relying on MITOMAP can be most misleading. A notorious example is the T1243C mutation, which was assessed to be novel and deemed to be associated with some (rare) disease simply because researchers did not realize that T1243C defines a deep branch in the Eurasian mtDNA phylogeny. The majority of 'novel' mutations suspected of being pathogenic are in actual fact known (and presumably neutral) polymorphisms (although unknown to MITOMAP), and this becomes glaringly evident when proper database searches and straightforward Internet queries are carried out.

mtDNA mutation C1494T, haplogroup A, and hearing loss in Chinese

Mutation C1494T in mitochondrial 12S rRNA gene was recently reported in two large Chinese families with aminoglycoside-induced and nonsyndromic hearing loss (AINHL) and was claimed to be pathogenic. This mutation, however, was first reported in a sample from central China in our previous study that was aimed to reconstruct East Asian mtDNA phylogeny. All these three mtDNAs formed a subclade defined by mutation C1494T in mtDNA haplogroup A. It thus seems that mutation C1494T is a haplogroup A-associated mutation and this matrilineal background may contribute a high risk for the penetrance of mutation C1494T in Chinese with AINHL. To test this hypothesis, we first genotyped mutation C1494T in 553 unrelated individuals from three regional Chinese populations and performed an extensive search for published complete or near-complete mtDNA data sets (>3000 mtDNAs), we then screened the C1494T mutation in 111 mtDNAs with haplogroup A status that were identified from 1823 subjects across China. The search for published mtDNA data sets revealed no other mtDNA besides the above-mentioned three carrying mutation C1494T. None of the 553 randomly selected individuals and the 111 haplogroup A mtDNAs was found to bear this mutation. Therefore, our results suggest that C1494T is a very rare event. The mtDNA haplogroup A background in general is unlikely to play an active role in the penetrance of mutation C1494T in AINHL.

Updating the East Asian mtDNA phylogeny: a prerequisite for the identification of pathogenic mutations

Knowledge about the world phylogeny of human mitochondrial DNA (mtDNA) is essential not only for evaluating the pathogenic role of specific mtDNA mutations but also for performing reliable association studies between mtDNA haplogroups and complex disorders. In the past few years, the main features of the East Asian portion of the mtDNA phylogeny have been determined on the basis of complete sequencing efforts, but representatives of several basal lineages were still lacking. Moreover, some recently published complete mtDNA sequences did apparently not fit into the known phylogenetic tree and conflicted with the established nomenclature. To refine the East Asian mtDNA tree and resolve data conflicts, we first completely sequenced 20 carefully selected mtDNAs--likely representatives of novel sub-haplogroups--and then, in order to distinguish diagnostic mutations of novel haplogroups from private variants, we applied a 'motif-search' procedure to a large sample collection. The novel information was incorporated into an updated East Asian mtDNA tree encompassing more than 1000 (near-) complete mtDNA sequences. A reassessment of the mtDNA data from a series of disease studies testified to the usefulness of such a refined mtDNA tree in evaluating the pathogenicity of mtDNA mutations. In particular, the claimed pathogenic role of mutations G3316A, T3394C, A4833G and G15497A appears to be most questionable as those initial claims were derived from anecdotal findings rather than e.g. appropriate association studies. Following a guideline based on the phylogenetic knowledge as proposed here could help avoiding similar problems in the future.

A reappraisal of complete mtDNA variation in East Asian families with hearing impairment

In a number of recent studies, we summarized the obvious errors and shortcomings that can be spotted in many (if not most) mitochondrial DNA (mtDNA) data sets published in medical genetics. We have reanalyzed here the complete mtDNA genome data published in various recent reports of East Asian families with hearing impairment, using a phylogenetic approach, in order to demonstrate the persistence of lab-specific mistakes in mtDNA genome sequencing in cases where those caveats were (deliberately) neglected. A phylogenetic reappraisal of complete mtDNAs with mutation A1555G (or G11778A) indeed supports the suggested lack of association between haplogroup background and phenotypic presentation of these mutations in East Asians. In contrast, the claimed pathogenicity of mutation T1095C in Chinese families with hearing impairment seems unsupported, basically because this mutation is rather basal in the mtDNA phylogeny, being specific to haplogroup M11 in East Asia. The roles of other haplogroup specific or associated variants, such as A827G, T961C, T1005C, in East Asian subjects with aminoglycoside-induced and non-syndromic hearing loss are also unclear in view of the known mtDNA phylogeny.

A critical reassessment of the role of mitochondria in tumorigenesis

BACKGROUND

Mitochondrial DNA (mtDNA) is being analyzed by an increasing number of laboratories in order to investigate its potential role as an active marker of tumorigenesis in various types of cancer. Here we question the conclusions drawn in most of these investigations, especially those published in high-rank cancer research journals, under the evidence that a significant number of these medical mtDNA studies are based on obviously flawed sequencing results.

METHODS AND FINDINGS

In our analyses, we take a phylogenetic approach and employ thorough database searches, which together have proven successful for detecting erroneous sequences in the fields of human population genetics and forensics. Apart from conceptual problems concerning the interpretation of mtDNA variation in tumorigenesis, in most cases, blocks of seemingly somatic mutations clearly point to contamination or sample mix-up and, therefore, have nothing to do with tumorigenesis.

CONCLUSION

The role of mitochondria in tumorigenesis remains unclarified. Our findings of laboratory errors in many contributions would represent only the tip of the iceberg since most published studies do not provide the raw sequence data for inspection, thus hindering a posteriori evaluation of the results. There is no precedent for such a concatenation of errors and misconceptions affecting a whole subfield of medical research.

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.

Traces of archaic mitochondrial lineages persist in Austronesian-speaking Formosan populations

Genetic affinities between aboriginal Taiwanese and populations from Oceania and Southeast Asia have previously been explored through analyses of mitochondrial DNA (mtDNA), Y chromosomal DNA, and human leukocyte antigen loci. Recent genetic studies have supported the "slow boat" and "entangled bank" models according to which the Polynesian migration can be seen as an expansion from Melanesia without any major direct genetic thread leading back to its initiation from Taiwan. We assessed mtDNA variation in 640 individuals from nine tribes of the central mountain ranges and east coast regions of Taiwan. In contrast to the Han populations, the tribes showed a low frequency of haplogroups D4 and G, and an absence of haplogroups A, C, Z, M9, and M10. Also, more than 85% of the maternal lineages were nested within haplogroups B4, B5a, F1a, F3b, E, and M7. Although indicating a common origin of the populations of insular Southeast Asia and Oceania, most mtDNA lineages in Taiwanese aboriginal populations are grouped separately from those found in China and the Taiwan general (Han) population, suggesting a prevalence in the Taiwanese aboriginal gene pool of its initial late Pleistocene settlers. Interestingly, from complete mtDNA sequencing information, most B4a lineages were associated with three coding region substitutions, defining a new subclade, B4a1a, that endorses the origin of Polynesian migration from Taiwan. Coalescence times of B4a1a were 13.2 +/- 3.8 thousand years (or 9.3 +/- 2.5 thousand years in Papuans and Polynesians). Considering the lack of a common specific Y chromosomal element shared by the Taiwanese aboriginals and Polynesians, the mtDNA evidence provided here is also consistent with the suggestion that the proto-Oceanic societies would have been mainly matrilocal.

A single nucleotide substitution that abolishes the initiator methionine codon of the GLDC gene is prevalent among patients with glycine encephalopathy in Jerusalem

Glycine encephalopathy (GE) (non-ketotic hyperglycinemia) is an autosomal recessive neurometabolic disease caused by defective activity of the glycine cleavage system. Clinically, patients present usually in the neonatal period with hypotonia, encephalopathy, hiccups and breath arrests with or without overt seizures. GE is considered rare, but its incidence is relatively high in several geographical areas around the world. We report a novel mutation causing GE in six extended Arab families, all from a small suburban village (population 5,000). A methionine to threonine change in the initiation codon of the glycine decarboxylase gene led to markedly reduced glycine decarboxylase mRNA levels and abolished glycine cleavage system activity.

Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia

To resolve the phylogeny of the autochthonous mitochondrial DNA (mtDNA) haplogroups of India and determine the relationship between the Indian and western Eurasian mtDNA pools more precisely, a diverse subset of 75 macrohaplogroup N lineages was chosen for complete sequencing from a collection of >800 control-region sequences sampled across India. We identified five new autochthonous haplogroups (R7, R8, R30, R31, and N5) and fully characterized the autochthonous haplogroups (R5, R6, N1d, U2a, U2b, and U2c) that were previously described only by first hypervariable segment (HVS-I) sequencing and coding-region restriction-fragment-length polymorphism analysis. Our findings demonstrate that the Indian mtDNA pool, even when restricted to macrohaplogroup N, harbors at least as many deepest-branching lineages as the western Eurasian mtDNA pool. Moreover, the distribution of the earliest branches within haplogroups M, N, and R across Eurasia and Oceania provides additional evidence for a three-founder-mtDNA scenario and a single migration route out of Africa.