Extraordinary claims require extraordinary evidence in asserted mtDNA biparental inheritance.

A breakthrough article published in PNAS by Luo et al. challenges a central dogma in biology which states that the mitochondrial DNA (mtDNA) in humans is inherited exclusively from the mother. We re-analyzed original FASTQ files and results reported by Luo et al. to investigate methodological issues (e.g. nuclear mitochondrial DNA or NUMTs, DNA rearrangements) that could lead to biological misinterpretations. A comprehensive analysis of their data reveals several methodological and analytical issues that must be carefully addressed before challenging the current paradigm. We first show that the probability of the findings described by the authors is extremely small (most likely below 10-37). The sequencing replicates from the same donors show aberrations in the variants detected that need further investigation to exclude contributions from other sources or methodological artifacts. Applying the principle of reductio ad absurdum, we demonstrate that the nuclear factor invoked by the authors to explain the phenomenon would need to be extraordinarily complex and precise to preclude linear accumulation of mtDNA lineages across generations, which would make the appearance of mixed haplotypes a much more frequent event in the population. We discuss alternate scenarios that explain findings of the same nature as reported by Luo et al., in the context of in-vitro fertilization and therapeutic mtDNA replacement ooplasmic transplantation.

HaploGrep 2: mitochondrial haplogroup classification in the era of high-throughput sequencing.

Mitochondrial DNA (mtDNA) profiles can be classified into phylogenetic clusters (haplogroups), which is of great relevance for evolutionary, forensic and medical genetics. With the extensive growth of the underlying phylogenetic tree summarizing the published mtDNA sequences, the manual process of haplogroup classification would be too time-consuming. The previously published classification tool HaploGrep provided an automatic way to address this issue. Here, we present the completely updated version HaploGrep 2 offering several advanced features, including a generic rule-based system for immediate quality control (QC). This allows detecting artificial recombinants and missing variants as well as annotating rare and phantom mutations. Furthermore, the handling of high-throughput data in form of VCF files is now directly supported. For data output, several graphical reports are generated in real time, such as a multiple sequence alignment format, a VCF format and extended haplogroup QC reports, all viewable directly within the application. In addition, HaploGrep 2 generates a publication-ready phylogenetic tree of all input samples encoded relative to the revised Cambridge Reference Sequence. Finally, new distance measures and optimizations of the algorithm increase accuracy and speed-up the application. HaploGrep 2 can be accessed freely and without any registration at http://haplogrep.uibk.ac.at.

The saga of the many studies wrongly associating mitochondrial DNA with breast cancer.

BACKGROUND: A large body of genetic research has focused on the potential role that mitochondrial DNA (mtDNA) variants might play on the predisposition to common and complex (multi-factorial) diseases. It has been argued however that many of these studies could be inconclusive due to artifacts related to genotyping errors or inadequate design. METHODS: Analyses of the data published in case-control breast cancer association studies have been performed using a phylogenetic-based approach. Variation observed in these studies has been interpreted in the light of data available on public resources, which now include over >27,000 complete mitochondrial sequences and the worldwide phylogeny determined by these mitogenomes. Complementary analyses were carried out using public datasets of partial mtDNA sequences, mainly corresponding to control-region segments. RESULTS: By way of example, we show here another kind of fallacy in these medical studies, namely, the phenomenon of SNP-SNP interaction wrongly applied to haploid data in a breast cancer study. We also reassessed the mutually conflicting studies suggesting some functional role of the non-synonymous polymorphism m.10398A>G (ND3 subunit of mitochondrial complex I) in breast cancer. In some studies, control groups were employed that showed an extremely odd haplogroup frequency spectrum compared to comparable information from much larger databases. Moreover, the use of inappropriate statistics signaled spurious "significance" in several instances. CONCLUSIONS: Every case-control study should come under scrutiny in regard to the plausibility of the control-group data presented and appropriateness of the statistical methods employed; and this is best done before potential publication.

The case for the continuing use of the revised Cambridge Reference Sequence (rCRS) and the standardization of notation in human mitochondrial DNA studies.

Since the determination in 1981 of the sequence of the human mitochondrial DNA (mtDNA) genome, the Cambridge Reference Sequence (CRS), has been used as the reference sequence to annotate mtDNA in molecular anthropology, forensic science and medical genetics. The CRS was eventually upgraded to the revised version (rCRS) in 1999. This reference sequence is a convenient device for recording mtDNA variation, although it has often been misunderstood as a wild-type (WT) or consensus sequence by medical geneticists. Recently, there has been a proposal to replace the rCRS with the so-called Reconstructed Sapiens Reference Sequence (RSRS). Even if it had been estimated accurately, the RSRS would be a cumbersome substitute for the rCRS, as the new proposal fuses--and thus confuses--the two distinct concepts of ancestral lineage and reference point for human mtDNA. Instead, we prefer to maintain the rCRS and to report mtDNA profiles by employing the hitherto predominant circumfix style. Tree diagrams could display mutations by using either the profile notation (in conventional short forms where appropriate) or in a root-upwards way with two suffixes indicating ancestral and derived nucleotides. This would guard against misunderstandings about reporting mtDNA variation. It is therefore neither necessary nor sensible to change the present reference sequence, the rCRS, in any way. The proposed switch to RSRS would inevitably lead to notational chaos, mistakes and misinterpretations.

Haplogrouping mitochondrial DNA sequences in Legal Medicine/Forensic Genetics.

Haplogrouping refers to the classification of (partial) mitochondrial DNA (mtDNA) sequences into haplogroups using the current knowledge of the worldwide mtDNA phylogeny. Haplogroup assignment of mtDNA control-region sequences assists in the focused comparison with closely related complete mtDNA sequences and thus serves two main goals in forensic genetics: first is the a posteriori quality analysis of sequencing results and second is the prediction of relevant coding-region sites for confirmation or further refinement of haplogroup status. The latter may be important in forensic casework where discrimination power needs to be as high as possible. However, most articles published in forensic genetics perform haplogrouping only in a rudimentary or incorrect way. The present study features PhyloTree as the key tool for assigning control-region sequences to haplogroups and elaborates on additional Web-based searches for finding near-matches with complete mtDNA genomes in the databases. In contrast, none of the automated haplogrouping tools available can yet compete with manual haplogrouping using PhyloTree plus additional Web-based searches, especially when confronted with artificial recombinants still present in forensic mtDNA datasets. We review and classify the various attempts at haplogrouping by using a multiplex approach or relying on automated haplogrouping. Furthermore, we re-examine a few articles in forensic journals providing mtDNA population data where appropriate haplogrouping following PhyloTree immediately highlights several kinds of sequence errors.

Mitochondrial genomes from modern horses reveal the major haplogroups that underwent domestication.

Archaeological and genetic evidence concerning the time and mode of wild horse (Equus ferus) domestication is still debated. High levels of genetic diversity in horse mtDNA have been detected when analyzing the control region; recurrent mutations, however, tend to blur the structure of the phylogenetic tree. Here, we brought the horse mtDNA phylogeny to the highest level of molecular resolution by analyzing 83 mitochondrial genomes from modern horses across Asia, Europe, the Middle East, and the Americas. Our data reveal 18 major haplogroups (A-R) with radiation times that are mostly confined to the Neolithic and later periods and place the root of the phylogeny corresponding to the Ancestral Mare Mitogenome at ~130-160 thousand years ago. All haplogroups were detected in modern horses from Asia, but F was only found in E. przewalskii--the only remaining wild horse. Therefore, a wide range of matrilineal lineages from the extinct E. ferus underwent domestication in the Eurasian steppes during the Eneolithic period and were transmitted to modern E. caballus breeds. Importantly, now that the major horse haplogroups have been defined, each with diagnostic mutational motifs (in both the coding and control regions), these haplotypes could be easily used to (i) classify well-preserved ancient remains, (ii) (re)assess the haplogroup variation of modern breeds, including Thoroughbreds, and (iii) evaluate the possible role of mtDNA backgrounds in racehorse performance.

Rapid coastal spread of First Americans: novel insights from South America's Southern Cone mitochondrial genomes.

It is now widely agreed that the Native American founders originated from a Beringian source population ~15-18 thousand years ago (kya) and rapidly populated all of the New World, probably mainly following the Pacific coastal route. However, details about the migration into the Americas and the routes pursued on the continent still remain unresolved, despite numerous genetic, archaeological, and linguistic investigations. To examine the pioneering peopling phase of the South American continent, we screened literature and mtDNA databases and identified two novel mitochondrial DNA (mtDNA) clades, here named D1g and D1j, within the pan-American haplogroup D1. They both show overall rare occurrences but local high frequencies, and are essentially restricted to populations from the Southern Cone of South America (Chile and Argentina). We selected and completely sequenced 43 D1g and D1j mtDNA genomes applying highest quality standards. Molecular and phylogeographic analyses revealed extensive variation within each of the two clades and possibly distinct dispersal patterns. Their age estimates agree with the dating of the earliest archaeological sites in South America and indicate that the Paleo-Indian spread along the entire longitude of the American double continent might have taken even <2000 yr. This study confirms that major sampling and sequencing efforts are mandatory for uncovering all of the most basal variation in the Native American mtDNA haplogroups and for clarification of Paleo-Indian migrations, by targeting, if possible, both the general mixed population of national states and autochthonous Native American groups, especially in South America.

Current next generation sequencing technology may not meet forensic standards.

In a Nature paper of 2010, the concern was raised that intra-individual mtDNA variation may be more pronounced than previously believed, in that heteroplasmies are common and vary markedly from tissue to tissue. This claim taken at face value would have considerable impact on forensic casework. It turns out however that the employed technology detected the germ-line variation relative to the reference sequence only incompletely: on average at least five mutations were missed per sample, as an in silico reassessment of the data reveals. Before one can really set out to access to entire mtDNA genome data with relative ease for forensic purposes, one needs careful calibration studies under strict forensic conditions-or might have to wait for another generation.

Is mitochondrial tRNA(phe) variant m.593T>C a synergistically pathogenic mutation in Chinese LHON families with m.11778G>A?

Mitochondrial transfer RNA (mt-tRNA) mutations have been reported to be associated with a variety of diseases. In a previous paper that studied the mtDNA background effect on clinical expression of Leber's hereditary optic neuropathy (LHON) in 182 Chinese families with m.11778G>A, we found a strikingly high frequency (7/182) of m.593T>C in the mitochondrially encoded tRNA phenylalanine (MT-TF) gene in unrelated LHON patients. To determine the potential role of m.593T>C in LHON, we compared the frequency of this variant in 479 LHON patients with m.11778G>A, 843 patients with clinical features of LHON but without the three known primary mutations, and 2374 Han Chinese from the general populations. The frequency of m.593T>C was higher in LHON patients (14/479) than in suspected LHON subjects (12/843) or in general controls (49/2374), but the difference was not statistically significant. The overall penetrance of LHON in families with both m.11778G>A and m.593T>C (44.6%) was also substantially higher than that of families with only m.11778G>A (32.9%) (P = 0.083). Secondary structure prediction of the MT-TF gene with the wild type or m.593T>C showed that this nucleotide change decreases the free energy. Electrophoretic mobility of the MT-TF genes with the wild type or m.593T>C transcribed in vitro further confirmed the change of secondary structure in the presence of this variant. Although our results could suggest a modest synergistic effect of variant m.593T>C on the LHON causing mutation m.11778G>A, the lack of statistical significance probably due to the relatively small sample size analyzed, makes necessary more studies to confirm this effect.

The archaeogenetics of Europe.

A new timescale has recently been established for human mitochondrial DNA (mtDNA) lineages, making mtDNA at present the most informative genetic marker system for studying European prehistory. Here, we review the new chronology and compare mtDNA with Y-chromosome patterns, in order to summarize what we have learnt from archaeogenetics concerning five episodes over the past 50,000 years which significantly contributed to the settlement history of Europe: the pioneer colonisation of the Upper Palaeolithic, the Late Glacial re-colonisation of the continent from southern refugia after the Last Glacial Maximum, the postglacial re-colonization of deserted areas after the Younger Dryas cold snap, the arrival of Near Easterners with an incipient Neolithic package, and the small-scale migrations along continent-wide economic exchange networks beginning with the Copper Age. The available data from uniparental genetic systems have already transformed our view of the prehistory of Europe, but our knowledge of these processes remains limited. Nevertheless, their legacy remains as sedimentary layers in the gene pool of modern Europeans, and our understanding of them will improve substantially when more mtDNAs are completely sequenced, the Y chromosome more thoroughly analysed, and haplotype blocks of the autosomal genome become amenable to phylogeographic studies.

Amerindian mitochondrial DNA haplogroups predominate in the population of Argentina: towards a first nationwide forensic mitochondrial DNA sequence database.

The study presents South American mitochondrial DNA (mtDNA) data from selected north (N = 98), central (N = 193) and south (N = 47) Argentinean populations. Sequence analysis of the complete mtDNA control region (CR, 16024-576) resulted in 288 unique haplotypes ignoring C-insertions around positions 16193, 309, and 573; the additional analysis of coding region single nucleotide polymorphisms enabled a fine classification of the described lineages. The Amerindian haplogroups were most frequent in the north and south representing more than 60% of the sequences. A slightly different situation was observed in central Argentina where the Amerindian haplogroups represented less than 50%, and the European contribution was more relevant. Particular clades of the Amerindian subhaplogroups turned out to be nearly region-specific. A minor contribution of African lineages was observed throughout the country. This comprehensive admixture of worldwide mtDNA lineages and the regional specificity of certain clades in the Argentinean population underscore the necessity of carefully selecting regional samples in order to develop a nationwide mtDNA database for forensic and anthropological purposes. The mtDNA sequencing and analysis were performed under EMPOP guidelines in order to attain high quality for the mtDNA database.

Mitochondrial genome evidence reveals successful Late Paleolithic settlement on the Tibetan Plateau.

Due to its numerous environmental extremes, the Tibetan Plateau--the world's highest plateau--is one of the most challenging areas of modern human settlement. Archaeological evidence dates the earliest settlement on the plateau to the Late Paleolithic, while previous genetic studies have traced the colonization event(s) to no earlier than the Neolithic. To explore whether the genetic continuity on the plateau has an exclusively Neolithic time depth, we studied mitochondrial DNA (mtDNA) genome variation within 6 regional Tibetan populations sampled from Tibet and neighboring areas. Our results confirm that the vast majority of Tibetan matrilineal components can trace their ancestry to Epipaleolithic and Neolithic immigrants from northern China during the mid-Holocene. Significantly, we also identified an infrequent novel haplogroup, M16, that branched off directly from the Eurasian M founder type. Its nearly exclusive distribution in Tibetan populations and ancient age (>21 kya) suggest that M16 may represent the genetic relics of the Late Paleolithic inhabitants on the plateau. This partial genetic continuity between the Paleolithic inhabitants and the contemporary Tibetan populations bridges the results and inferences from archaeology, history, and genetics.

High mitochondrial DNA stability in B-cell chronic lymphocytic leukemia.

BACKGROUND: Chronic Lymphocytic Leukemia (CLL) leads to progressive accumulation of lymphocytes in the blood, bone marrow, and lymphatic tissues. Previous findings have suggested that the mtDNA could play an important role in CLL. METHODOLOGY/PRINCIPAL FINDINGS: The mitochondrial DNA (mtDNA) control-region was analyzed in lymphocyte cell DNA extracts and compared with their granulocyte counterpart extract of 146 patients suffering from B-Cell CLL; B-CLL (all recruited from the Basque country). Major efforts were undertaken to rule out methodological artefacts that would render a high false positive rate for mtDNA instabilities and thus lead to erroneous interpretation of sequence instabilities. Only twenty instabilities were finally confirmed, most of them affecting the homopolymeric stretch located in the second hypervariable segment (HVS-II) around position 310, which is well known to constitute an extreme mutational hotspot of length polymorphism, as these mutations are frequently observed in the general human population. A critical revision of the findings in previous studies indicates a lack of proper methodological standards, which eventually led to an overinterpretation of the role of the mtDNA in CLL tumorigenesis. CONCLUSIONS/SIGNIFICANCE: Our results suggest that mtDNA instability is not the primary causal factor in B-CLL. A secondary role of mtDNA mutations cannot be fully ruled out under the hypothesis that the progressive accumulation of mtDNA instabilities could finally contribute to the tumoral process. Recommendations are given that would help to minimize erroneous interpretation of sequencing results in mtDNA studies in tumorigenesis.

Mitochondrial haplogroup U5b3: a distant echo of the epipaleolithic in Italy and the legacy of the early Sardinians.

There are extensive data indicating that some glacial refuge zones of southern Europe (Franco-Cantabria, Balkans, and Ukraine) were major genetic sources for the human recolonization of the continent at the beginning of the Holocene. Intriguingly, there is no genetic evidence that the refuge area located in the Italian Peninsula contributed to this process. Here we show, through phylogeographic analyses of mitochondrial DNA (mtDNA) variation performed at the highest level of molecular resolution (52 entire mitochondrial genomes), that the most likely homeland for U5b3-a haplogroup present at a very low frequency across Europe-was the Italian Peninsula. In contrast to mtDNA haplogroups that expanded from other refugia, the Holocene expansion of haplogroup U5b3 toward the North was restricted by the Alps and occurred only along the Mediterranean coasts, mainly toward nearby Provence (southern France). From there, approximately 7,000-9,000 years ago, a subclade of this haplogroup moved to Sardinia, possibly as a result of the obsidian trade that linked the two regions, leaving a distinctive signature in the modern people of the island. This scenario strikingly matches the age, distribution, and postulated geographic source of a Sardinian Y chromosome haplogroup (I2a2-M26), a paradigmatic case in the European context of a founder event marking both female and male lineages.

The multifaceted origin of taurine cattle reflected by the mitochondrial genome.

A Neolithic domestication of taurine cattle in the Fertile Crescent from local aurochsen (Bos primigenius) is generally accepted, but a genetic contribution from European aurochsen has been proposed. Here we performed a survey of a large number of taurine cattle mitochondrial DNA (mtDNA) control regions from numerous European breeds confirming the overall clustering within haplogroups (T1, T2 and T3) of Near Eastern ancestry, but also identifying eight mtDNAs (1.3%) that did not fit in haplogroup T. Sequencing of the entire mitochondrial genome showed that four mtDNAs formed a novel branch (haplogroup R) which, after the deep bifurcation that gave rise to the taurine and zebuine lineages, constitutes the earliest known split in the mtDNA phylogeny of B. primigenius. The remaining four mtDNAs were members of the recently discovered haplogroup Q. Phylogeographic data indicate that R mtDNAs were derived from female European aurochsen, possibly in the Italian Peninsula, and sporadically included in domestic herds. In contrast, the available data suggest that Q mtDNAs and T subclades were involved in the same Neolithic event of domestication in the Near East. Thus, the existence of novel (and rare) taurine haplogroups highlights a multifaceted genetic legacy from distinct B. primigenius populations. Taking into account that the maternally transmitted mtDNA tends to underestimate the extent of gene flow from European aurochsen, the detection of the R mtDNAs in autochthonous breeds, some of which are endangered, identifies an unexpected reservoir of genetic variation that should be carefully preserved.

Contamination and sample mix-up can best explain some patterns of mtDNA instabilities in buccal cells and oral squamous cell carcinoma.

The study of somatic DNA instabilities constitutes a debatable topic because different causes can lead to seeming DNA alteration patterns between different cells or tissues from the same individual. Carcinogenesis or the action of a particular toxic could generate such patterns, and this is in fact the leitmotif of a number of studies on mitochondrial DNA (mtDNA) instability. Patterns of seeming instabilities could also arise from technical errors at any stage of the analysis (DNA extraction, amplification, mutation screening/sequencing, and documentation). Specifically, inadvertent DNA contamination or sample mixing would yield mosaic variation that could be erroneously interpreted as real mutation differences (instabilities) between tissues from the same individual. From the very beginning, mtDNA studies comparing cancerous to non-cancerous tissues have suffered from such mosaic results. We demonstrate here that the phylogenetic linkage of whole arrays of mtDNA mutations provides strong evidence of artificial recombination in previous studies on buccal cells and oral squamous cell carcinoma.

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.

Distinctive Paleo-Indian migration routes from Beringia marked by two rare mtDNA haplogroups.

BACKGROUND: It is widely accepted that the ancestors of Native Americans arrived in the New World via Beringia approximately 10 to 30 thousand years ago (kya). However, the arrival time(s), number of expansion events, and migration routes into the Western Hemisphere remain controversial because linguistic, archaeological, and genetic evidence have not yet provided coherent answers. Notably, most of the genetic evidence has been acquired from the analysis of the common pan-American mitochondrial DNA (mtDNA) haplogroups. In this study, we have instead identified and analyzed mtDNAs belonging to two rare Native American haplogroups named D4h3 and X2a. RESULTS: Phylogeographic analyses at the highest level of molecular resolution (69 entire mitochondrial genomes) reveal that two almost concomitant paths of migration from Beringia led to the Paleo-Indian dispersal approximately 15-17 kya. Haplogroup D4h3 spread into the Americas along the Pacific coast, whereas X2a entered through the ice-free corridor between the Laurentide and Cordilleran ice sheets. The examination of an additional 276 entire mtDNA sequences provides similar entry times for all common Native American haplogroups, thus indicating at least a dual origin for Paleo- Indians. CONCLUSIONS: A dual origin for the first Americans is a striking novelty from the genetic point of view, and it makes plausible a scenario positing that within a rather short period of time, there may have been several entries into the Americas from a dynamically changing Beringian source. Moreover, this implies that most probably more than one language family was carried along with the Paleo-Indians.