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.

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.

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.

Mitochondrial portraits of human populations using median networks.

Analysis of variation in the hypervariable region of mitochondrial DNA (mtDNA) has emerged as an important tool for studying human evolution and migration. However, attempts to reconstruct optimal intraspecific mtDNA phylogenies frequently fail because parallel mutation events partly obscure the true evolutionary pathways. This makes it inadvisable to present a single phylogenetic tree at the expense of neglecting equally acceptable ones. As an alternative, we propose a novel network approach for portraying mtDNA relationships. For small sample sizes (< approximately 50), an unmodified median network contains all most parsimonious trees, displays graphically the full information content of the sequence data, and can easily be generated by hand. For larger sample sizes, we reduce the complexity of the network by identifying parallelisms. This reduction procedure is guided by a compatibility argument and an additional source of phylogenetic information: the frequencies of the mitochondrial haplotypes. As a spin-off, our approach can also assist in identifying sequencing errors, which manifest themselves in implausible network substructures. We illustrate the advantages of our approach with several examples from existing data sets.

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.

Weak hierarchies associated with similarity measures--an additive clustering technique.

A new and apparently rather useful and natural concept in cluster analysis is studied: given a similarity measure on a set of objects, a sub-set is regarded as a cluster if any two objects a, b inside this sub-set have greater similarity than any third object outside has to at least one of a, b. These clusters then form a closure system which can be described as a hypergraph without triangles. Conversely, given such a system, one may attach some weight to each cluster and then compose a similarity measure additively, by letting the similarity of a pair be the sum of weights of the clusters containing that particular pair. The original clusters can be reconstructed from the obtained similarity measure. This clustering model is thus located between the general additive clustering model of Shepard and Arabie (1979) and the standard hierarchical model. Potential applications include fitting dendrograms with few additional nonnested clusters and simultaneous representation of some families of multiple dendrograms (in particular, two-dendrogram solutions), as well as assisting the search for phylogenetic relationships by proposing a somewhat larger system of possibly relevant "family groups", from which an appropriate choice (based on additional insight or individual preferences) remains to be made.

Split decomposition: a new and useful approach to phylogenetic analysis of distance data.

In order to analyze the structure inherent to a matrix of dissimilarities (such as evolutionary distances) we propose to use a new technique called split decomposition. This method accurately dissects the given dissimilarity measure as a sum of elementary "split" metrics plus a (small) residue. The split summands identify related groups which are susceptible to further interpretation when casted against the available biological information. Reanalysis of previously published ribosomal RNA data sets using split decomposition illustrate the potential of this approach.

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.

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.

Origin and evolution of Native American mtDNA variation: a reappraisal.

The timing and number of prehistoric migrations involved in the settlement of the American continent is subject to intense debate. Here, we reanalyze Native American control region mtDNA data and demonstrate that only an appropriate phylogenetic analysis accompanied by an appreciation of demographic factors allows us to discern different migrations and to estimate their ages. Reappraising 574 mtDNA control region sequences from aboriginal Siberians and Native Americans, we confirm in agreement with linguistic, archaeological and climatic evidence that (i) the major wave of migration brought one population, ancestral to the Amerinds, from northeastern Siberia to America 20,000-25,000 years ago and (ii) a rapid expansion of a Beringian source population took place at the end of the Younger Dryas glacial phase approximately 11,300 years ago, ancestral to present Eskimo and Na-Dene populations.

Mitochondrial footprints of human expansions in Africa.

mtDNA studies support an African origin for modern Eurasians, but expansion events within Africa have not previously been investigated. We have therefore analyzed 407 mtDNA control-region sequences from 13 African ethnic groups. A number of sequences (13%) were highly divergent and coalesced on the "mitochondrial Eve" in Africans. The remaining sequences also ultimately coalesced on this sequence but fell into four major clusters whose starlike phylogenies testify to demographic expansions. The oldest of these African expansions dates to approximately 60,000-80,000 years ago. Eurasian sequences are derived from essentially one sequence within this ancient cluster, even though a diverse mitochondrial pool was present in Africa at the time.

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.

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.

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.

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.

Phylogeography of mitochondrial DNA in western Europe.

For most of the past century, prehistorians have had to rely on the fossil and archaeological records in order to reconstruct the past. In the last few decades, this evidence has been substantially supplemented from classical human genetics. More recently, phylogenetic analyses of DNA sequences that incorporate geographical information have provided a high-resolution tool for the investigation of prehistoric demographic events, such as founder effects and population expansions. These events can be dated using a molecular clock when the mutation rate and founder haplotypes are known. We have previously applied such methods to sequence data from the mitochondrial DNA control region, to suggest that most extant mitochondrial sequences in western Europe have a local ancestry in the Early Upper Palaeolithic, with a smaller proportion arriving from the Near East in the Neolithic. Here, we describe a cladistic notation for mitochondrial variation and expand upon our earlier analysis to present a more detailed portrait of the European mitochondrial record.

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.