Pseudomitochondrial genome haunts disease studies.

The accidental amplification of nuclear mitochondrial pseudogenes (NUMTs) can pose a serious problem for mitochondrial disease studies. This report shows that the mutation spectrum left by spurious amplification of a NUMT can be detected because it usually differs considerably from the authentic natural spectrum. This study examined the problem introduced by an ND5 gene NUMT that was recorded in a proband with hearing loss and reviews other disease studies erroneously reporting NUMT variation as genuine mutations in their patients. NUMTs can emerge in population genetic studies, as exemplified here by cases in this study and from published sources. Appropriate database searches and a phylogenetic approach can prevent hasty claims for novelty of mitochondrial DNA (mtDNA) variants inadvertently derived from NUMTs and help to direct investigators to the real source.

Consistent treatment of length variants in the human mtDNA control region: a reappraisal.

In forensic science, as well as in molecular anthropology and medical genetics, human mitochondrial DNA (mtDNA) variation is being recorded by aligning mtDNA sequences to the revised Cambridge reference sequence (rCRS). This task is straightforward for the vast majority of nucleotide positions but appears to be difficult for some short sequence stretches, namely, in regions displaying length variation. Earlier guidelines for imposing a unique alignment relied on binary alignment to a standard sequence (the rCRS) and used additional priority rules for resolving ambiguities. It turns out, however, that these rules have not been applied rigorously and led to inconsistent nomenclature. There is no way to adapt the priority rules in a reasonable way because binary alignment to a standard sequence is bound to produce artificial alignments that may place sequences separated by a single mutation at mismatch distance larger than 1. To remedy the situation, we propose a phylogenetic approach for multiple alignment and resulting notation.

Phylogeographic investigations: the role of trees in forensic genetics.

The human mitochondrial DNA (mtDNA) genome is commonly analyzed in various disciplines, such as population, medical, and forensic genetics, but conceptual and scientific exchange between them is still limited. Here we review several aspects of the mtDNA phylogeny that are particularly--but not exclusively--of interest to the forensic community. Among the issues that arise, we emphasize the importance of integrating evolutionary concepts into the forensic routine. We also discuss topics such as mtDNA mutation-rate heterogeneity and the weight of evidence, ethnic affiliations of mtDNA profiles, and the abuse of reference databases. Finally, we show the usefulness of coding-region variation in a forensic context.

Quality assessment of DNA sequence data: autopsy of a mis-sequenced mtDNA population sample.

Published DNA data sets constitute a body of sequencing results resting in silico that are supposed to reflect the variation of (once) living cells. In cases where the DNA variation reported is suspected to be fraught with artefacts, an autopsy of the full body of data is needed to clarify the amount and causes of mis-sequencing. In this paper we elaborate on strategies that allow a clear-cut identification of the problems in severely flawed mtDNA data. This approach is applied, by way of example, to a data set of HVS-I sequences from the Caucasus, published by Nasidze & Stoneking in 2001. These data bear numerous ambiguous nucleotide positions and suffer from an even higher number of phantom mutations, indicating that severe biochemical problems adversely influenced those sequencing results at the time. Furthermore, systematic omission of sequences with a long C-stretch (incurred by a transition at position 16189) must have severely biased the data set. Since no complete correction of these data has appeared to date, this example of mis-sequencing necessitates circumstantial evidence that is bullet-proof.

More evidence for non-maternal inheritance of mitochondrial DNA?

BACKGROUND: A single case of paternal co-transmission of mitochondrial DNA (mtDNA) in humans has been reported so far. OBJECTIVE: To find potential instances of non-maternal inheritance of mtDNA. METHODS: Published medical case studies (of single patients) were searched for irregular mtDNA patterns by comparing the given haplotype information for different clones or tissues with the worldwide mtDNA database as known to date-a method that has proved robust and reliable for the detection of flawed mtDNA sequence data. RESULTS: More than 20 studies were found reporting clear cut instances with mtDNAs of different ancestries in single individuals. As examples, cases are reviewed from recent published reports which, at face value, may be taken as evidence for paternal inheritance of mtDNA or recombination. CONCLUSIONS: Multiple types (or recombinant types) of quite dissimilar mitochondrial DNA from different parts of the known mtDNA phylogeny are often reported in single individuals. From re-analyses and corrigenda of forensic mtDNA data, it is apparent that the phenomenon of mixed or mosaic mtDNA can be ascribed solely to contamination and sample mix up.

Artificial recombination in forensic mtDNA population databases.

Artificial recombination of two or more mitochondrial DNA fragments from different samples would constitute a serious cause of error in forensic DNA typing, and yet one can demonstrate that such events have happened in the preparation of several published mtDNA databases. Focussed database searches, phylogenetic analysis, and network representations can highlight mosaic patterns and thus pinpoint sample mix-up. Therefore, we suggest that this approach should be applied to data prior to publication in order to uncover such errors in time.

Identification of Native American founder mtDNAs through the analysis of complete mtDNA sequences: some caveats.

In this study, a detailed analysis of both previously published and new data was performed to determine whether complete, or almost complete, mtDNA sequences can resolve the long-debated issue of which Asian mtDNAs were founder sequences for the Native American mtDNA pool. Unfortunately, we now know that coding region data and their analysis are not without problems. To obtain and report reasonably correct sequences does not seem to be a trivial task, and to discriminate between Asian and Native American mtDNA ancestries may be more complex than previously believed. It is essential to take into account the effects of mutational hot spots in both the control and coding regions, so that the number of apparent Native American mtDNA founder sequences is not erroneously inflated. As we report here, a careful analysis of all available data indicates that there is very little evidence that more than five founder mtDNA sequences entered Beringia before the Last Glacial Maximum and left their traces in the current Native American mtDNA pool.

Mitochondrial portrait of the Cabo Verde archipelago: the Senegambian outpost of Atlantic slave trade.

In order to study the matrilineal genetic composition in Cabo Verde (Republic of Cape Verde), an archipelago that used to serve as a Portuguese entrepôt of the Atlantic slave trade, we have analysed a total of 292 mtDNAs sampled from the seven inhabited islands for the hypervariable segment I (HVS-I) and some characteristic RFLPs of the coding regions. The different settlement history of the northwestern group of the islands is well reflected in the mtDNA pool. The total Cabo Verde sample clearly displays the characteristic mitochondrial features of the Atlantic fringe of western Africa and testifies to almost no mitochondrial input from the Portuguese colonizers.

Detecting errors in mtDNA data by phylogenetic analysis.

Sequencing and documenting a sample of homologous DNA stretches is prone to copying errors in a way rather analogous to the biological replication process. Previous attempts at obtaining representative mtDNA sequences, typically of the control region, for evolutionary studies or forensic purposes have yielded rather unsatisfactory results in many cases. The key ingredient in pinpointing problems with given data is the phylogenetic analysis of closely related mtDNAs within the framework of an established worldwide phylogeny that is supported by coding region information. We develop some general rules by which likely errors in data tables can readily be detected without rereading whole sequences repeatedly. Following these guidelines, one can expect to lower the error rate by at least an order of magnitude, although it will still be hard to beat the mitochondrial gamma polymerase in precision.

A signal, from human mtDNA, of postglacial recolonization in Europe.

Mitochondrial HVS-I sequences from 10,365 subjects belonging to 56 populations/geographical regions of western Eurasia and northern Africa were first surveyed for the presence of the T-->C transition at nucleotide position 16298, a mutation which has previously been shown to characterize haplogroup V mtDNAs. All mtDNAs with this mutation were then screened for a number of diagnostic RFLP sites, revealing two major subsets of mtDNAs. One is haplogroup V proper, and the other has been termed "pre*V," since it predates V phylogenetically. The rather uncommon pre*V tends to be scattered throughout Europe (and northwestern Africa), whereas V attains two peaks of frequency: one situated in southwestern Europe and one in the Saami of northern Scandinavia. Geographical distributions and ages support the scenario that pre*V originated in Europe before the Last Glacial Maximum (LGM), whereas the more recently derived haplogroup V arose in a southwestern European refugium soon after the LGM. The arrival of V in eastern/central Europe, however, occurred much later, possibly with (post-)Neolithic contacts. The distribution of haplogroup V mtDNAs in modern European populations would thus, at least in part, reflect the pattern of postglacial human recolonization from that refugium, affecting even the Saami. Overall, the present study shows that the dissection of mtDNA variation into small and well-defined evolutionary units is an essential step in the identification of spatial frequency patterns. Mass screening of a few markers identified using complete mtDNA sequences promises to be an efficient strategy for inferring features of human prehistory.

Phylogeography of the human mitochondrial haplogroup L3e: a snapshot of African prehistory and Atlantic slave trade.

The mtDNA haplogroup L3e, which is identified by the restriction site +2349 MboI within the Afro-Eurasian superhaplogroup L3 (-3592 HpaI), is omnipresent in Africa but virtually absent in Eurasia (except for neighbouring areas with limited genetic exchange). L3e was hitherto poorly characterised in terms of HVS-I motifs, as the ancestral HVS-I type of L3e cannot be distinguished from the putative HVS-I ancestor of the entire L3 (differing from the CRS by a transition at np 16223). An MboI screening at np 2349 of a large number of Brazilian and Caribbean mtDNAs (encompassing numerous mtDNAs of African ancestry), now reveals that L3e is subdivided into four principal clades, each characterised by a single mutation in HVS-I, with additional support coming from HVS-II and partial RFLP analysis. The apparently oldest of these clades (transition at np 16327) occurs mainly in central Africa and was probably carried to southern Africa with the Bantu expansion(s). The most frequent clade (transition at np 16320) testifies to a pronounced expansion event in the mid-Holocene and seems to be prominent in many Bantu groups from all of Africa. In contrast, one clade (transition at np 16264) is essentially restricted to Atlantic western Africa (including Cabo Verde). We propose a tentative L3e phylogeny that is based on 197 HVS-I sequences. We conclude that haplogroup L3e originated in central or eastern Africa about 46,000 (+/-14,000) years ago, and was a hitchhiker of much later dispersal and local expansion events, with the rise of food production and iron smelting. Enforced migration of African slaves to the Americas translocated L3e mitochondria, the descendants of which in Brazil and the Caribbean still reflect their different regional African ancestries.

Tracing European founder lineages in the Near Eastern mtDNA pool.

Founder analysis is a method for analysis of nonrecombining DNA sequence data, with the aim of identification and dating of migrations into new territory. The method picks out founder sequence types in potential source populations and dates lineage clusters deriving from them in the settlement zone of interest. Here, using mtDNA, we apply the approach to the colonization of Europe, to estimate the proportion of modern lineages whose ancestors arrived during each major phase of settlement. To estimate the Palaeolithic and Neolithic contributions to European mtDNA diversity more accurately than was previously achievable, we have now extended the Near Eastern, European, and northern-Caucasus databases to 1,234, 2, 804, and 208 samples, respectively. Both back-migration into the source population and recurrent mutation in the source and derived populations represent major obstacles to this approach. We have developed phylogenetic criteria to take account of both these factors, and we suggest a way to account for multiple dispersals of common sequence types. We conclude that (i) there has been substantial back-migration into the Near East, (ii) the majority of extant mtDNA lineages entered Europe in several waves during the Upper Palaeolithic, (iii) there was a founder effect or bottleneck associated with the Last Glacial Maximum, 20,000 years ago, from which derives the largest fraction of surviving lineages, and (iv) the immigrant Neolithic component is likely to comprise less than one-quarter of the mtDNA pool of modern Europeans.

mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion.

The Eskimo-Aleut language phylum is distributed from coastal Siberia across Alaska and Canada to Greenland and is well distinguished from the neighboring Na Dene languages. Genetically, however, the distinction between Na Dene and Eskimo-Aleut speakers is less clear. In order to improve the genetic characterization of Eskimos in general and Greenlanders in particular, we have sequenced hypervariable segment I (HVS-I) of the mitochondrial DNA (mtDNA) control region and typed relevant RFLP sites in the mtDNA of 82 Eskimos from Greenland. A comparison of our data with published sequences demonstrates major mtDNA types shared between Na Dene and Eskimo, indicating a common Beringian history within the Holocene. We further confirm the presence of an Eskimo-specific mtDNA subgroup characterized by nucleotide position 16265G within mtDNA group A2. This subgroup is found in all Eskimo groups analyzed so far and is estimated to have originated <3,000 years ago. A founder analysis of all Eskimo and Chukchi A2 types indicates that the Siberian and Greenland ancestral mtDNA pools separated around the time when the Neo-Eskimo culture emerged. The Greenland mtDNA types are a subset of the Alaskan mtDNA variation: they lack the groups D2 and D3 found in Siberia and Alaska and are exclusively A2 but at the same time lack the A2 root type. The data are in agreement with the view that the present Greenland Eskimos essentially descend from Alaskan Neo-Eskimos. European mtDNA types are absent in our Eskimo sample.

Median networks: speedy construction and greedy reduction, one simulation, and two case studies from human mtDNA.

Molecular data sets characterized by few phylogenetically informative characters with a broad spectrum of mutation rates, such as intraspecific control-region sequence variation of human mitochondrial DNA (mtDNA), can be usefully visualized in the form of median networks. Here we provide a step-by-step guide to the construction of such networks by hand. We improve upon a previously implemented algorithm by outlining an efficient parametrized strategy amenable to large data sets, greedy reduction, which makes it possible to reconstruct some of the confounding recurrent mutations. This entails some postprocessing as well, which assists in capturing more parsimonious solutions. To simplify the creation of the resulting network by hand, we describe a speedy approach to network construction, based on a careful planning of the processing order. A coalescent simulation tailored to human mtDNA variation in Eurasia testifies to the usefulness of reduced median networks, while highlighting notorious problems faced by all phylogenetic methods in this context. Finally, we discuss two case studies involving the comparison of characters in the two hypervariable segments of the human mtDNA control region in the light of the worldwide control-region sequence database, as well as additional restriction fragment length polymorphism information. We conclude that only a minority of the mutations that hit the second segment occur at sites that would have a mutation rate comparable to those at most sites in the first segment. Discarding the known "noisy" sites of the second segment enhances the analysis.

The ancestry of Brazilian mtDNA lineages.

We have analyzed 247 Brazilian mtDNAs for hypervariable segment (HVS)-I and selected restriction fragment-length-polymorphism sites, to assess their ancestry in different continents. The total sample showed nearly equal amounts of Native American, African, and European matrilineal genetic contribution but with regional differences within Brazil. The mtDNA pool of present-day Brazilians clearly reflects the imprints of the early Portuguese colonization process (involving directional mating), as well as the recent immigrant waves (from Europe) of the last century. The subset of 99 mtDNAs from the southeastern region encompasses nearly all mtDNA haplogroups observed in the total Brazilian sample; for this regional subset, HVS-II was analyzed, providing, in particular, some novel details of the African mtDNA phylogeny.

Genetic evidence of an early exit of Homo sapiens sapiens from Africa through eastern Africa.

The out-of-Africa scenario has hitherto provided little evidence for the precise route by which modern humans left Africa. Two major routes of dispersal have been hypothesized: one through North Africa into the Levant, documented by fossil remains, and one through Ethiopia along South Asia, for which little, if any, evidence exists. Mitochondrial DNA (mtDNA) can be used to trace maternal ancestry. The geographic distribution and variation of mtDNAs can be highly informative in defining potential range expansions and migration routes in the distant past. The mitochondrial haplogroup M, first regarded as an ancient marker of East-Asian origin, has been found at high frequency in India and Ethiopia, raising the question of its origin. (A haplogroup is a group of haplotypes that share some sequence variations.) Its variation and geographical distribution suggest that Asian haplogroup M separated from eastern-African haplogroup M more than 50,000 years ago. Two other variants (489C and 10873C) also support a single origin of haplogroup M in Africa. These findings, together with the virtual absence of haplogroup M in the Levant and its high frequency in the South-Arabian peninsula, render M the first genetic indicator for the hypothesized exit route from Africa through eastern Africa/western India. This was possibly the only successful early dispersal event of modern humans out of Africa.

Median-joining networks for inferring intraspecific phylogenies.

Reconstructing phylogenies from intraspecific data (such as human mitochondrial DNA variation) is often a challenging task because of large sample sizes and small genetic distances between individuals. The resulting multitude of plausible trees is best expressed by a network which displays alternative potential evolutionary paths in the form of cycles. We present a method ("median joining" [MJ]) for constructing networks from recombination-free population data that combines features of Kruskal's algorithm for finding minimum spanning trees by favoring short connections, and Farris's maximum-parsimony (MP) heuristic algorithm, which sequentially adds new vertices called "median vectors", except that our MJ method does not resolve ties. The MJ method is hence closely related to the earlier approach of Foulds, Hendy, and Penny for estimating MP trees but can be adjusted to the level of homoplasy by setting a parameter epsilon. Unlike our earlier reduced median (RM) network method, MJ is applicable to multistate characters (e.g., amino acid sequences). An additional feature is the speed of the implemented algorithm: a sample of 800 worldwide mtDNA hypervariable segment I sequences requires less than 3 h on a Pentium 120 PC. The MJ method is demonstrated on a Tibetan mitochondrial DNA RFLP data set.

Phylogeographic patterns of mtDNA reflecting the colonization of the Canary Islands.

Although the Canary Islands were settled by humans, possibly of Berber origin, as late as 2500 years ago, the precise course and numbers of early migrations to the archipelago remain controversial. We have therefore analysed mtDNA variation (HVS-I as well as selected RFLP sites) in 300 individuals from the seven Canary Islands. The distribution and variation across the islands in a specific mtDNA clade of Northwest African ancestry suggest that there was one dominant initial settlement process that affected all the islands, from east to west. This indicates that a certain genetic affinity of present-day Canary Islanders to Northwest African Berbers mainly stems from the autochthonous population rather than slaves captured on the neighbouring African coast. The slave trade after the European conquest left measurable, though minor, traces in the mtDNA pool of the Canary Islands, which in its majority testifies to the European immigration.