Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Highly Accessed Research article

Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa

Erwan Pennarun1*, Toomas Kivisild12, Ene Metspalu1, Mait Metspalu1, Tuuli Reisberg1, Jean-Paul Moisan3, Doron M Behar13, Sacha C Jones4 and Richard Villems15

Author Affiliations

1 Estonian Biocentre and Department of Evolutionary Biology, University of Tartu, Tartu, Estonia

2 Division of Biological Anthropology, University of Cambridge, Cambridge, United Kingdom

3 Molecular Medicine Laboratory, Rambam Health Care Campus, Haifa, Israel

4 McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, United Kingdom

5 Estonian Academy of Sciences, Tallinn, Estonia

For all author emails, please log on.

BMC Evolutionary Biology 2012, 12:234  doi:10.1186/1471-2148-12-234

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/12/234


Received:17 February 2012
Accepted:20 October 2012
Published:3 December 2012

© 2012 Pennarun et al.; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License ( http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

A Southwest Asian origin and dispersal to North Africa in the Early Upper Palaeolithic era has been inferred in previous studies for mtDNA haplogroups M1 and U6. Both haplogroups have been proposed to show similar geographic patterns and shared demographic histories.

Results

We report here 24 M1 and 33 U6 new complete mtDNA sequences that allow us to refine the existing phylogeny of these haplogroups. The resulting phylogenetic information was used to genotype a further 131 M1 and 91 U6 samples to determine the geographic spread of their sub-clades. No southwest Asian specific clades for M1 or U6 were discovered. U6 and M1 frequencies in North Africa, the Middle East and Europe do not follow similar patterns, and their sub-clade divisions do not appear to be compatible with their shared history reaching back to the Early Upper Palaeolithic. The Bayesian Skyline Plots testify to non-overlapping phases of expansion, and the haplogroups’ phylogenies suggest that there are U6 sub-clades that expanded earlier than those in M1. Some M1 and U6 sub-clades could be linked with certain events. For example, U6a1 and M1b, with their coalescent ages of ~20,000–22,000 years ago and earliest inferred expansion in northwest Africa, could coincide with the flourishing of the Iberomaurusian industry, whilst U6b and M1b1 appeared at the time of the Capsian culture.

Conclusions

Our high-resolution phylogenetic dissection of both haplogroups and coalescent time assessments suggest that the extant main branching pattern of both haplogroups arose and diversified in the mid-later Upper Palaeolithic, with some sub-clades concomitantly with the expansion of the Iberomaurusian industry. Carriers of these maternal lineages have been later absorbed into and diversified further during the spread of Afro-Asiatic languages in North and East Africa.

Keywords:
mtDNA haplogroups M1 and U6; Afro-Asiatic languages; North Africa

Background

The North African mitochondrial DNA (mtDNA) genetic pool has been shown to reflect influence from different regions, with sizeable portions of lineages from Sub-Saharan Africa, the Middle East, and others that diversified perhaps first in Europe [1-10], a pattern also shown with autosomal data [11]. The geographic patterns of some of the haplogroups that constitute the North African mtDNA pool have been singled out as being more informative about early population histories than others; for example, the variation in haplogroup U6 [1,12], a haplogroup that has been termed “the main indigenous North African cluster” [13], and, to a lesser extent the variation in M1, which is more predominantly found in Eastern Africa/Ethiopia [14-16]. U6 and M1 both share the feature of being African-specific sub-clades of haplogroups otherwise spread only in non-African populations. Indeed, whilst most U clades are found in North Africa and in Eurasia, as far as the Ganges Basin, U6 is virtually restricted to North (West) Africa. For macro-haplogroup M, this African connection is even more puzzling, as haplogroups belonging to M are mostly found only in South, Central and East Asia, the Americas and Oceania, where no M1 has yet been reported.

The Palaeolithic archaeological record of North Africa is spatially and temporally diverse, revealing a variety of technological shifts during the later Pleistocene period. The Aterian, a regional variant of the Middle Palaeolithic (or Middle Stone Age), was previously thought to have existed ~40,000–20,000 years ago (KYA), and argued to mark the earliest modern humans in North Africa. These dates have been drastically reassessed and the upper bound is now closer to ~115 KYA [17] or even as old as ~145 KYA [18]. The transition from the Middle Palaeolithic to Upper Palaeolithic in North Africa is characterised by the appearance of the “Dabban”, an industry that is restricted to Cyrenaica in northeast Libya and represented at the caves of Hagfet ed Dabba and Haua Fteah [19]. Whilst a techno-typological shift occurred within the Dabban ~33 KYA [19], starker changes in the archaeological record occurred throughout North Africa and Southwest Asia ~23-20 KYA, represented by the widespread appearance of backed bladelet technologies. The appearance of these backed bladelet industries more or less coincides with the timing of the Last Glacial Maximum (LGM) (~23-18 KYA), including: ~21 KYA in Upper Egypt [20]; ~20 KYA at Haua Fteah with the Oranian [21]; the Iberomaurusian expansion in the Jebel Gharbi ~20 KYA [22]; and the first Iberomaurusian at Tamar Hat in Algeria ~20 KYA [23]. The earliest Iberomaurusian sites in Morocco appear to be only slightly younger ~18 KYA [24]. Whilst backed bladelet production is broadly shared across the different regions of North and East Africa, there was also a level of regional cultural diversity during this period, possibly mirroring a diversification of populations. The Sahara Desert expanded considerably during the LGM, perhaps concentrating human groups along the North African coastal belt and the Nile Valley. Climatic conditions improved in North Africa ~15 KYA, marking the beginning of a dramatic arid-to-humid transition [25]. This increase in humidity may have opened up ecological corridors, connecting North and Sub-Saharan Africa and allowing population dispersals between the two regions. An additional arid-humid transition occurred at 11.5–11 KYA [25]; this period coincides with a widespread change in the archaeological record that marks the beginning of Capsian lithic technologies. The Capsian is argued to have developed in situ in North Africa, marking a continuity from the Iberomaurusian and Oranian into the Capsian [21,24,26].

Given the geographical specificity of mtDNA haplogroups U6 and M1, some studies have investigated their potential implication in the peopling of North Africa [5,27-30], whilst some earlier studies assumed that M1 diverged from other M lineages prior to the early dispersals of Homo sapiens out of Africa ~60–70 KYA [14,15]. However, most research that has followed explains its presence in Africa by a back-migration from Asia [5,31]. Dating of the U6 and M1 variation in African and Middle Eastern populations has been at the centre of the debate on the timing of the back-migration to Africa and, in particular, whether these haplogroups co-dispersed with certain archaeological cultures or languages. A thorough study by Olivieri and co-authors [29] proposed that both M1 and U6 were involved in an early dispersal, 40–45KYA, from Southwest Asia to North Africa in association with the first arrival of anatomically modern humans in the Mediterranean region. Considering this time frame, it was suggested, furthermore, that the spread of Aurignacian culture in Europe and the Dabban industry in North Africa derived from the same source. This interpretation was questioned by Forster and Romano who, referring to the geographic correlates, questioned this evidence and proposed that, alternatively, the spread of these haplogroups could be potentially be explained by more recent events, perhaps contemporary to the dispersal of populations speaking Afro-Asiatic (AA) languages [32].

In this study, we re-evaluate the timeframe for M1 and U6 variations and their patterns of geographic spread at the resolution of complete mtDNA sequences using a range of phylogeographic and statistical methods. We try to assess to what extent the phylogeographies of U6 and M1 are correlated with each other and, indirectly, with the spread of AA languages. In order to address these questions, a survey of more than 5700 mtDNAs was undertaken, covering a broad geographic region encompassing North and East Africa, the Near and Middle East and the Caucasus. 24 M1 and 33 U6 complete mtDNA sequences were determined and, with the refined phylogenetic trees for M1 and U6 drawn, we use this information to genotype a further 131 M1 and 91 U6 samples of different geographic origin.

Results

Phylogeny, phylogeography and coalescent estimates of M1 and U6

Our genotyping of haplogroup U6 and M1 defining markers, analysed in combination with published data, confirmed earlier findings that these two haplogroups are present all over the Mediterranean Basin: both are particularly prevalent in the southern Mediterranean and M1 reaches as far away as East Africa (Figures  1a and b). Yet, some of their peak frequencies only partially overlap in Northwest Africa. In contrast to high frequencies of M1 sub-clades, haplogroup U6 is rare in East/Northeast Africa and the Middle East, and is virtually absent in the Caucasus (Table  1). Nevertheless, both haplogroups are by and large confined to the area where AA languages are spoken nowadays, being rare or absent in areas where other language families are dominant (Figure  1).

Table 1. Frequency of Haplogroups M1 and U6 in the geographic regions from this study

thumbnailFigure 1. Spatial distribution of haplogroup M1 and U6, with languages’ phyla. Frequency maps were obtained using Surfer 8 (Golden Software, Inc.). The Kriging procedure was used and the dataset was congregated with existing ones [29] and updated with the present study, as well as the data in [27,28]. Figure  1a: frequency map for haplogroup M1. Figure  1b: frequency map for haplogroup U6. Red dots indicate the populations geographic locations.

Concerning the estimated coalescent ages, Table  2 shows an excerpt of the Additional file 1, and contains only some coalescent ages relevant in a broader context, whilst Figure  2 shows a schematic tree of M1 and U6 phylogenies (See Additional file 2, Additional file 3, Additional file 4 and Additional file 5 for detailed phylogenies). The use of a different method (e.g. using only the synonymous mutations rather than all the mutations present in the mtDNA coding region; see [33]) for estimating molecular coalescent ages gives younger results than previously published [27-29] for both haplogroups with the coalescence of U6 at ~35 KYA and M1 at ~29 KYA. U6 is mostly prevalent in Northwest Africa (Additional file 4 and Additional file 5), a similar occurrence for M1b, which contrasts with M1a, the most diverse sub-clade of M1, for which most of its sub-clades are prevalent in East Africa. Both M1b and M1a have close coalescent ages around the LGM: ~20 and ~21 KYA respectively. M1a1 is the most diverse clade of M1a and is found in virtually all the populations where M1 has been sampled (except in Guinea-Bissau). Again, a variety of its sub-clades are more frequent in East Africa and, interestingly, a large subset of M1a1 samples could not be ascribed to any of its known sub-clades (Additional file 3). It is noteworthy to point out that all the Caucasian samples fall into just one sub-clade, M1a1b2, with no variation present at the intermediate level of resolution (Additional file 3), signature of a likely founder effect.

Additional file 1. – Coalescent age estimates for M1 and U6 and their most frequent sub-clades. Soaresa: These estimates include some sequences that are not complete, and are given just for indication, see the left panel for estimates based only on complete sequences.

Format: XLSX Size: 16KB Download fileOpen Data

Additional file 2. – Phylogenetic tree based on 105 M1 full sequences. All positions are scored against the RSRS [58] and are transitions, unless followed by a capital letter that marks the resulting transversion. Indels are scored with i or d, heteroplasmies follow the IUB code and reversal to ancestral state by an exclamation mark (!), double back mutations by two exclamation marks (!!). The positions are colour coded according to their status: purple – non-coding; blue – non-synonymous; and black – synonymous. Variations in the C tracts were mostly ignored (i.e., 16182C, 16193C, 309+2C, etc.) unless stated on the tree. The box containing the sample ID is colour coded according to the publications from which they were retrieved (See the main text for the full reference), and below it their geographic origin is colour coded (See Additional file 7 for the specifics). Sequences available only for the coding region, or for which some parts are missing, are flagged with a yellow mark under the geographic origin. The order for the root mutation(s) for M1a1g, M1a1h, M1a7 and M1b2c were deduced from additional partial sequencing (See Additional file 3).

Format: XLSX Size: 41KB Download fileOpen Data

Additional file 3. – Network based on 236 M1 samples. All positions are scored against the RSRS [58] and are transitions, unless followed by a capital letter that marks the resulting transversion. Indels are scored with i or d, heteroplasmies follow the IUB code and reversal to ancestral state by an exclamation mark (!), double back mutations by two exclamation marks (!!). The positions are colour coded according to their status: purple – non-coding; blue – non-synonymous; and black – synonymous. Variations in the C tracts were ignored (i.e., 16182C, 16193C, 309+2C, etc.). The box containing the sample ID is colour coded according to the publications which they are from (See the main text for the full reference), and below it their geographic origin is colour coded (See Additional file 7 for the specifics).

Format: XLSX Size: 49KB Download fileOpen Data

Additional file 4. – Phylogenetic tree based on 139 U6 full sequences. All positions are scored against the RSRS [58] and are transitions, unless followed by a capital letter that marks the resulting transversion. Indels are scored with i or d, heteroplasmies follow the IUB code and reversal to ancestral state by an exclamation mark (!). The positions are colour coded according to their status: purple – non-coding; blue – non-synonymous; and black – synonymous. Variations in the C tracts were mostly ignored (i.e., 16182C, 16193C, 309+2C, etc.) unless stated on the tree. The box containing the sample ID is colour coded according to the publications which they are from (See the main text for the full reference), and below it their geographic origin is colour coded (See Additional file 7 for the specifics). Sequences available only for the coding region are flagged with a yellow mark under the geographic origin. The potential reticulation created by position 150 between sub-clades U6a3a and U6a3c was resolved on the more frequent occurrence of 150 in various different haplogroups’ backgrounds (See [72]). We refined here the phylogeny of the Canary-specific branch formerly known as U6b1 [27,29]. There is an array of 2 common mutations before the branch splits into the so-called Canary-specific branch and one apparently specific to Northwest Africa. We propose therefore to rename U6b1a as U6b1a1 to comply with the revised phylogeny. The mutations order of some clades (U6a1a1b, U6a1a2, U6a2b, U6a2b1, U6a3d, U6a3d1, U6a3d1a, U6a6a, U6b1b1) was deduced for additional partial typing (see Additional file 5).

Format: XLSX Size: 64KB Download fileOpen Data

Additional file 5. – Network based on 230 U6 samples. All positions are scored against the RSRS [58] and are transitions, unless followed by a capital letter that marks the resulting transversion. Indels are scored with i or d, heteroplasmies follow the IUB code and reversal to ancestral state by an exclamation mark (!). The positions are colour coded according to their status: purple – non-coding; blue – non-synonymous; and black – synonymous. Variations in the C tracts were ignored (i.e., 16182C, 16193C, 309+2C, etc.). The box containing the sample ID is colour coded according to the publications which they are from (See the main text for the full reference), and below it their geographic origin is colour coded (See Additional file 7 for the specifics). The reticulation created by position 150 in U6a3’s clades is left unresolved.

Format: XLSX Size: 63KB Download fileOpen Data

Table 2. Coalescent age estimates for M1 and U6 and some of the most frequent sub-clades

thumbnailFigure 2. Schematic tree of Haplogroup M1 and U6. The tree, rooted in L3, shows the major sub-haplogroups of M1 and U6. The branching is phylogenetically correct, but the branches length is not accurate.

The most diversified sub-clade of U6 is U6a, largely due to the richness of its sub-clades in Northwest Africa. One of its sub-clades, U6a2, has been so far detected only in East African and Middle Eastern populations. Contrary to M1, various clades of U6 predate the LGM, including U6a, which is very close to the overall age of U6 (~33 KYA vs. ~36 KYA). Confirming some previous observations [27,29,30], U6b and U6c were confined in our samples to Northwest Africa.

Bayesian Skyline Plot analyses

We tested our panel of full sequences for expansion signal(s) using Bayesian Skyline Plots (BSP) that estimate past effective population size (Ne) dynamics on the basis of sequence data [35]. The method does not rely on any pre-specified parametric model of demographic history. However, its results should be taken with caution, as the curve representing Ne could also reflect changes in the sub-structure of the population rather than its true size variation [36,37], and that the reconstruction of Ne might also be biased by the purifying selection acting on the mtDNA genome [33,34,38,39]. Yet, as here both lineages have a similar ratio of non-synonymous to synonymous mutations (0,63), this effect is not likely to explain differences that we have found. Figure  3 displays the BSPs for M1 and U6. For each simulation, the median of the other haplogroup is overlaid for comparison. We also indicate the coalescent ages and the 95% CI of some sub-clades based on the full genome as in [34], hence the coalescent ages reported in this section may differ with the ones from the previous section. The rate by Soares et al.[34] is applied here as the entire mtDNA genomic sequence is used for the BSP analyses, whereas the rate by Loogväli et al.[33] applies only to the coding region. Nonetheless, the two different approaches offer similar estimates (see Table  2 and Additional file 1 for more detail).

thumbnailFigure 3. Bayesian Skyline Plot for Haplogroups M1 and U6. The BSPs show the variation of the Effective Population Size (Ne) through Time for M1 (Figure  3a) and U6 (Figure  3b) based on the full mitochondrial genomes. The axis scales are identical for both plots. For comparison, the median of the second haplogroup is shown in grey, but not the 95% HPD. Overlaid on the plots are the coalescent ages of some relevant sub-haplogroups, with the vertical bars indicating the calculated coalescent ages (using the calculator from [34]) and the horizontal ones their 95% confidence interval.

For U6, the initial expansion seems to more or less coincide with the ~26–27 KYA estimated coalescent age (based on full sequences) of U6a, the most diverse and prevalent sub-clade of U6. This expansion appears to have continued at a somewhat equal rate, gradually slowing down, until the curve even drops slightly, and eventually a new expansion phase takes place around ~6–7 KYA. For M1, the slope of the curve is steeper, with two clearly visible expansion phases. The first inflexion is ~22 KYA, slightly older than the estimated coalescent ages for M1a and M1b, with a strong increase until reaching a plateau at ~15 KYA. The second phase occurs at ~10–11 KYA, a time around which the estimated coalescent ages of various sub-clades of M1 fall (e.g., M1b1 and M1a1b). By directly comparing the median curves of U6 and M1, representing the past population dynamics extracted from the molecular data, it appears unlikely that the demographic histories of these haplogroups entirely overlap, both in terms of the timing of expansion phases, as well as the magnitude of these expansions.

Mantel correlation tests

To explore whether the frequencies of M1 and U6 across a geographic range of populations correlate with languages we used Mantel correlation tests. Notably, when M1 and U6 are grouped, or with U6 alone, no significant correlation is found, neither between genes, nor geography, nor language (Table  3). A correlation is found both between geography and language only for M1, being higher with geography than with language. To see which M1 clade contributes the most to this signal, the tests were done with M1a and M1b sorted separately. No correlation could be found between M1b and geography and/or language, whilst M1a was significantly correlated both with geography and language.

Table 3. Mantel test to assess the correlation between genes and geography and/or language

Discussion

Origins of M1 and U6, their implications in the colonisation of North Africa, and some of its archaeological landmarks

A Southwest Asian origin has been proposed for U6 and M1 [27-29]. Yet, this claim remains speculative unless some novel “earlier” Southwest Asian-specific clades, distinct from the known haplogroups, are found in which the described so far M1 and U6 lineages are nested. Claims for basal mutations shared with M1 have recently been made in the case of haplogroup M51 and M20 (both East Asian-specific clades [40,41]): They share a root mutation (C14110T) with M1. However, one should be cautious with phylogenetic inferences drawn from these findings because this mutation is not unique in the phylogeny of mtDNA: it also occurs in the background of non-M haplogroups and therefore identity by descent within haplogroup M remains uncertain. Unfortunately, the sampling of extant populations of Africa and West Asia may not solve the question of their origin.

Assuming that M1 and U6 were introduced to Africa by a dispersal event from Asia, it would be difficult to accept their involvement in the first demographic spread of anatomically modern humans around 40–45 KYA, as suggested by Olivieri et al. (2006), [29] who associated these two clades with the spread of Dabban industry in Africa. It has indeed been previously suggested that the colonisation of North Africa from the Levant took place during the early Upper Paleolithic, as marked by the “Dabban” industry in North Africa [42]. However, comparison of early Upper Palaeolithic artefacts from Haua Fteah and Ksar Akil does not support the notion that the early Dabban of Cyrenaica is an evidence of a population migration from the Levant into North Africa [43]. Marks [44] also noted differences between the two areas in terms of the methods of blade production, further arguing against a demographic connection between the regions. Likewise, the new coalescent date estimates for M1 obtained in this study are not compatible with the model implying the spread of M1 in Africa during the Early Upper Palaeolithic, 40–45 KYA.

Given the sequence data from 242 complete sequences and genotype data of 222 mtDNAs, we were unable to find conclusive evidence that any of the geographic regions of Africa or the Middle East would stand out as being uniquely or even significantly enriched with deep-rooted clades of U6 and M1 not found elsewhere. Whilst several U6 sub-clades seem to be confined to Northwest Africa, this pattern may be the result of drift and founder effects over many thousands of years and does not necessarily suggest that Northwest Africa was the geographic source of U6 dispersals in Africa. Similarly in the case of M1b1, the Northwest African frequency pattern is apparent, whilst its counterpart, M1a, is widely spread around the Mediterranean Basin, and its current diversity is highest in East Africa. The age estimates of M1b and U6a1 (~20 KYA) together with their Northwest African-spread patterns are more consistent with their appearance during or after the spread of the Iberomaurusian culture, rather than explainable by an earlier spread of the Dabban industry. Furthermore, there is no evidence that the Dabban industry spread to NW Africa, as indicated earlier [43,44]. When taking the most recent common ancestor estimates of mtDNA haplogroups at face value and comparing them with relevant archeological horizons, then the Capsian culture also appears to be a possible candidate for the co-spread of sub-clades U6b and M1b1.

Although mtDNA is a single locus, some parallels concerning the African expansion of M1 and U6 can be drawn from autosomal data. In a recent study, Behar and colleagues explored the genome-wide diversity of the Jewish Diaspora with regard to that of their host populations, as well as the Middle East [45]. In their supplemental figure four, results of analyses undertaken with the software ADMIXTURE are shown, and specifically at K=10, an ancestry component depicted in deep purple colour appears. Interestingly, its proportion is particularly high amongst Mozabite Berbers, who have very high frequencies of M1 and U6 [12]. This deep purple colour is also present at a fairly high frequency amongst Moroccans, and to a lesser extent amongst Ethiopians, both Jewish and non-Jewish, and Egyptians. Its proportion in the Near Eastern populations is by far smaller than in the African ones.

Mimicry of M1 and U6

A mimicry between U6 and M1 has been suggested [28,29]. Both are likely derived from a non-African ancestral clade at a similar time depth and both are largely confined to North and East Africa and the Middle East in their present-day geographic distribution. It seems, however, that the mimicry breaks down when analysing in further detail the coalescent times and frequency patterns of their sub-clades. Even at the general level, U6 is hardly found outside Northwest Africa, whilst M1 is ubiquitous throughout North Africa, East Africa and the Middle East, reaching also northern Caucasus. The coalescent age for U6a is almost 10 000 years older than that for either M1a or M1b, and most of its sub-clades coalesce before or around the LGM. In contrast, most of the estimates for M1a and M1b sub-clades are post-LGM. Also, the BSP analyses show that M1 and U6 have probably experienced different molecular histories. While the curves representing the median Ne for U6 and M1 overlap when taking the 95% HPD into account, the median curves themselves do differ. The earlier age of U6 is apparent, and though the U6 median follows a rather steady rate until declining, M1 bears testimony to two distinct expansion events. Although Hg U6 also experienced two expansion events, they do not superimpose on those of M1. It should be noted that the U6 curve should be taken with precaution as close to one half of the full U6 sequences used are from Europe. When taking into account the geography and running the BSP simulations by separate regions, it appears that the decline around 8–9 KYA is actually almost entirely driven by the European sequences (See Additional file 6). Unfortunately, it was not possible to ascertain if some of the signals present for M1 are also regional, because the number of regional sequences is too low. However, the proportion of “geographic outliers” in M1 is lower than in the case of U6.

Additional file 6. – BSP for U6 based on North African and European sequences separately. For the North African and European sequences, only a few independent runs were done to ascertain that convergence was reached. The 10 convergence runs for all sets of sequences are shown for comparison.

Format: PDF Size: 665KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

M1, U6 and the Afro-Asiatic language family

It has been proposed that M1 and U6, or some of their sub-clades, could be linked with the spread of AA languages [27,29,31]. Some of the main criteria for this are due to their geographical spread broadly overlapping with regions where AA languages are spoken today. There are currently two hypotheses about where AA languages originated. One places it in Northeast Africa, on the coast of the Red Sea [46,47], linking the reconstructed proto-Afro-Asiatic vocabulary to pre-Neolithic cultures in the Levant and their predecessors in southeast Egypt and northeastern Sudan, whilst the second places it in the Levant [48] , and emphasises the Neolithic component in the Afro-Asiatic cognates. Notably, even the earliest time frame (~10 KYA or more) considered by the linguists [47,49] for the earliest splits in the language family are more recent than the ages of U6 or M1 and their major sub-clades. However, if the sub-clades of M1 and U6 were to be involved in the dispersal event associated with the Afro-Asiatic languages they had to exist at the moment of the launch of this event, and therefore the fact that these sub-clades are older makes them plausible candidates for such dispersal. However, when considering M1 and U6 as a whole, or U6 alone, no correlation with language (and geography) was found with the current data, indicating for U6 that its expansion was not concomitant with that of the AA.

Concerning haplogroup M1 individually, a significant correlation with languages was observed. Furthermore, within M1, it appears that the correlation is mostly due to M1a. However, given the small sample size of M1b, any potential signal correlating with language might not be detectable. Interestingly, M1a has a likely East African origin, but its coalescent age of ~21 KYA still largely predates that of the proto-AA. Maybe a sub-clade of M1a would still give a similar correlation, but there are not sufficient samples to allow splitting M1a into its various sub-clades, and to test for a correlation. Although we found a correlation, limited sample sizes do not allow drawing unambiguous connection between genes and languages. Furthermore, it is also possible that this putative sub-clade of M1 does not testify for the expansion of AA speaking people, but was already present among the people who inhabited the area before the spread of the AA languages.

Conclusions

Our analyses do not support the model according to which mtDNA haplogroups M1 and U6 represent an early dispersal event of anatomically modern humans at around 40–45 KYA in association with the spread of Dabban industry in North Africa as proposed earlier [28,29]. A West Asian origin for these haplogroups still remains a viable hypothesis as sister clades of U (and ancestral to it, macro-hg N (including R)) and M are spread overwhelmingly outside Africa, notably in Eurasia, even though the phylogeographic data on extant populations do not present a clear support for it. Our estimates of coalescent times and demographic analyses of U6 and M1 variations suggest that their spread in North and East Africa is largely due to a number of demographic events, predominantly occurring at the end of the peak of as well as after the LGM, but largely before the Holocene. Hence, some of the topologically earliest sub-clades of U6 and M1 may have been involved in the origin and spread of the essentially North African Iberomaurusian culture, and the observed correlations with languages make it likely that the North and East African carriers of the two matrilineages have been absorbed into the expanding Afro-Asiatic languages-speaking people in the area, but in phylogeographically differential ways.

Methods

Samples

From over 5700 samples spanning Europe and countries around the Mediterranean Basin and beyond, 153 M1 and 121 U6 samples were identified based on their HVSI variation and then confirmed by RFLP (all unrelated individuals, who gave their informed consent). Samples from the literature/GenBank were retrieved, including: 77 M1 (2 from [50], 1 from [51], 3 from [52], 8 from [28], 1 from [53], 2 from [39], 1 from [54], 2 from [55], 2 from [56], 51 from [29], 1 from [57], 3 from [58] and 3 from [59]); and 93 U6 (1 from [60], 6 samples from Family Tree DNA deposited in GenBank, 1 from [61], 2 from [53], 1 from [62], 12 from [27], 30 from [29], 2 from [57], 39 from [30] and 7 from [58]). 9 samples were corrected (See Additional file 7 for the corrected positions) compared to their current GenBank entry at the time of this article’s submission, including 2 from [28], 1 from [55] and 5 from [27] (Dr. Vicente Cabrera’s personal communication). Also, 2 M1 and 3 U6 samples were kindly provided by Family Tree DNA (with some U6 samples having a potential match to sequences deposited in GenBank, see Additional file 7 and its legend for more details), bringing the total to 236 and 230 samples for M1 and U6 respectively (See Additional file 7 for detailed information). All the work complied with the Helsinki Declaration of Ethical Principles (59th WMA General Assembly, Seoul October 2008). The Estonian Basic Research project SF0182474 was approved by the Research Ethics Committee of the Estonian Biocentre.

Additional file 7. – List of the 466 samples. GeoBroad abbreviations are as follow: WA: West Africa; EA: East Africa; EUR: Europe; NE: Near/Middle East; NWA: North-West Africa; NEA: North-East Africa. See the main text for the full references. In the case of 4 samples originally provided by Familly Tree DNA, 3 samples have an identical sequence that matches an entry in GenBank, and as they cannot be differentiated, they have not been separately deposited into GenBank. For the last sample, there are two entries in GenBank with an identical sequence, and thus that sample as well has not been deposited into GenBank.

Format: XLSX Size: 38KB Download fileOpen Data

Sequencing, SNP typing

The 153 M1 samples from this study have been screened for approximately 2 kb of coding region in 4 separate fragments (between nps 700–1080, 6250–6990, 12590–13146, 14750–15580) chosen to cover some SNP-defining sub-clades of M1 based on previous knowledge [16,55]. 22 samples were fully sequenced following previously published protocol [63], and slightly modified. Based on the tree drawn from 105 full (or nearly full) sequences (Additional file 2), some SNPs have been typed in order to place precisely all the samples on the tree (See Additional file 3 and Additional file 8 for the full typing information).

Additional file 8. – Genotyping information for 153 M1 samples.

Format: XLSX Size: 69KB Download fileOpen Data

For the 121 U6 samples, several fragments have been amplified to type SNPs of interest based on the samples’ HVS I information (See Additional file 9 for full typing information), as well as from the tree based on 139 full sequences (See Additional file 4).

Additional file 9. – Genotyping information for 121 U6 samples.

Format: XLSX Size: 90KB Download fileOpen Data

Phylogenic tree, network

The trees and network were drawn by hand and checked with Network 4.5.1.0 ( http://www.fluxus-technology.com/ webcite[64]). If needed, a weighing scheme was used for highly recurrent polymorphisms.

Coalescent age estimates

For the coalescent age calculations, the rho (ρ) statistic and standard deviation were used as in [65,66] but see [67] for a critical assessment of it. Different rates were used: For the coding region, rate [33] is used, and for the full genome, estimates were calculated with the calculator provided in [34]. For all calculations, 2 M1 samples from [50] and the 3 M1 from [52] were discarded – the first ones were missing several portions of the coding region, and the second ones seemed to exhibit potential sequence errors (See [68] for details). For the full genome calculations, further samples were discarded, as their control region is not reported (M1: 1 from [53] and 2 from [39]; U6: 2 from [53], and 1 from [62]).

Mantel test

The haplotype of each sample was composed of all the polymorphisms detected in the coding region during the genotyping of each haplogroup, with the missing polymorphisms assumed to be similar to the RSRS [58], plus the control region (16024–16400). The HVS II was excluded of the haplotype, as it was only sequenced in some samples and, unlike for the coding region, it cannot be reasonably assumed that a specific polymorphism is absent in a different sample. For some populations the sample size was small, in which case they were grouped with a close geographic neighbour sharing the same language family. If this was not possible, the samples/populations were excluded (See Table  3 for details). The genetic distance matrices were based on Slatkin’s linearised FSTs. Because the language families present in the data are too divergent to rank and order them, we used a binary approach, with populations (or grouped populations) speaking a language from the same language family given a distance of 0, and a distance of 1 otherwise. Mantel tests were done with Arlequin 3.5.1.2 [69], with 100,000 permutations.

Bayesian Skyline Plot

The Bayesian Skyline Plot (BSP) [35] is a graphical depiction of the variation of the effective population size (Ne) through time. BSP analyses were performed with the software BEAST v1.5.4 [70]. The GTR substitution model was used with a gamma distribution, plus invariant sites. An uncorrelated lognormal relaxed clock was applied [71]. The whole mitochondrial genome was used, and runs were performed for 40,000,000 generations with 20 groups. In order to assure that convergence was reached, several independent runs were done (See Additional files 10 and Additional file 11 for M1 and U6, respectively). Also, the impact of the number of groups, which was user defined, was explored (See Additional file 12 and Additional file 13 for M1 and U6, respectively) by increments of 5, from 5 to 50. The axes were converted into their final units (effective population size vs. time) with a rate of 1,695 × 10-8[34] and a generation time of 25 years. But in order to take into account the purifying pressure acting on the whole molecule, ρ was deduced from the data, and then entered into the calculator provided in [34], resulting in a time scale which can be put in comparison with the coalescent ages calculated in the same way. Accordingly, the samples which were not available in full, with some missing parts, or which might suffer from errors (See the paragraph on coalescent age estimates) were not included in the analyses. Additional file 14 and Additional file 15 show the differences when the overall rate vs. the rate taking into account purifying selection [34] are used for plotting the results; the major impact being upon the time axis, whereas the impact on the effective population size or the overall shape of the curve are minimal.

Additional file 10. – 10 independent BSP runs for M1 with 20 groups. All runs were performed using the same parameters.

Format: PDF Size: 894KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 11. – 10 independent BSP run analyses for U6 with 20 groups. All runs were performed using the same parameters.

Format: PDF Size: 893KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 12. – BSP for M1 with groups varying from 5 to 50 groups, in increments of 5.

Format: PDF Size: 445KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 13. – BSP for U6 with groups varying from 5 to 50 groups, in increments of 5.

Format: PDF Size: 459KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 14. – BSP for M1 with the corrected rate versus uncorrected. The uncorrected rate use a rate of 1,695 x 10-8[34], and the corrected rate was deduced with the deduced rho values from the time, using the calculator from [34].

Format: PDF Size: 350KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Additional file 15. – BSP for U6 with the corrected rate versus uncorrected. The uncorrected rate use a rate of 1,695 × 10-8[34], and the corrected rate was deduced with the deduced rho values from the time, using the calculator from [34].

Format: PDF Size: 347KB Download file

This file can be viewed with: Adobe Acrobat ReaderOpen Data

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

RV, TK and JPM conceived the study, EP and DMB performed the full sequencing, EP, EM and TR genotyped the samples, EP performed the statistical analyses. EP, TK and RV interpreted the results and drafted the manuscript, MM designed some figures, SCJ revised the manuscript focusing on the archeological aspect. All authors read and approved the final manuscript.

Acknowledgements

In fond memory of the late Professor Jean-Paul Moisan, without whose contribution this study would not have been possible. We thank Ille Hilpus, Viljo Soo and Chrystelle Richard for technical assistance. We wish to thank Vicente M. Cabrera for providing the revised sequences from several of his studies, and Alexei Drummond and Simon Ho for guidance on BSP analyses and helpful discussions, as well as Phillip Endicott. We thank as well the EU European Regional Development Fund through the Centre of Excellence in Genomics to Estonian Biocentre and University of Tartu. The project was supported by Estonian Basic Research grant SF0182474; Tartu University grant (PBGMR06901) to TK; Estonian Science Foundation grants (7858) to EM and (8973) to MM; Estonian Ministry of Education and Research (0180142s08) and European Commission grant (ECOGENE205419) to MM and RV.

References

  1. Rando JC, Pinto F, Gonzalez AM, Hernandez M, et al.: Mitochondrial DNA analysis of northwest African populations reveals genetic exchanges with European, near-eastern, and sub-Saharan populations.

    Ann Hum Genet 1998, 62(Pt 6):531-550. PubMed Abstract | Publisher Full Text OpenURL

  2. Plaza S, Calafell F, Helal A, Bouzerna N, Lefranc G, et al.: Joining the pillars of Hercules: mtDNA sequences show multidirectional gene flow in the western Mediterranean.

    Ann Hum Genet 2003, 67(4):312-328. PubMed Abstract | Publisher Full Text OpenURL

  3. Coudray C, Olivieri A, Achilli A, Pala M, Melhaoui M, Cherkaoui M, El-Chennawi F, Kossmann M, Torroni A, Dugoujon JM: The complex and diversified mitochondrial gene pool of Berber populations.

    Ann Hum Genet 2009, 73:196-214. PubMed Abstract | Publisher Full Text OpenURL

  4. Brakez Z, Bosch E, Izaabel H, Akhayat O, Comas D, et al.: Human mitochondrial DNA sequence variation in the Moroccan population of the Souss area.

    Ann Hum Biol 2001, 28(3):295-307. PubMed Abstract | Publisher Full Text OpenURL

  5. Richards M, Rengo C, Cruciani F, Gratrix F, Wilson JF, et al.: Extensive female-mediated gene flow from sub-Saharan Africa into near eastern Arab populations.

    Am J Hum Genet 2003, 72(4):1058-1064. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  6. Cherni L, Loueslati BY, Pereira L, Ennafaa H, Amorim A, El Gaaied AB: Female gene pools of Berber and Arab neighboring communities in central Tunisia: microstructure of mtDNA variation in North Africa.

    Hum Biol 2005, 77(1):61-70. PubMed Abstract | Publisher Full Text OpenURL

  7. Fadhlaoui-Zid K, Plaza S, Calafell F, Ben Amor M, Comas D, AB Eg: Mitochondrial DNA heterogeneity in Tunisian Berbers.

    Ann Hum Genet 2004, 68:222-233. PubMed Abstract | Publisher Full Text OpenURL

  8. Atkinson QD, Gray RD, Drummond AJ: mtDNA variation predicts population size in humans and reveals a major Southern Asian chapter in human prehistory.

    Mol Biol Evol 2008, 25(2):468-474. PubMed Abstract | Publisher Full Text OpenURL

  9. Pereira L, Cerny V, Cerezo M, Silva NM, Hajek M, et al.: Linking the sub-Saharan and West Eurasian gene pools: maternal and paternal heritage of the Tuareg nomads from the African Sahel.

    Eur J Hum Genet 2010, 18(8):915-923. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  10. Fadhlaoui-Zid K, Rodríguez-Botigué L, Naoui N, Benammar-Elgaaied A, et al.: Mitochondrial DNA structure in North Africa reveals a genetic discontinuity in the Nile Valley.

    Am J Phys Anthropol 2011, 145(1):107-117. PubMed Abstract | Publisher Full Text OpenURL

  11. Henn BM, Botigue LR, Gravel S, Wang W, et al.: Genomic ancestry of North Africans supports back-to-Africa migrations.

    PLoS Genet 2012, 8(1):e1002397. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  12. Corte-Real HB, Macaulay VA, Richards MB, et al.: Genetic diversity in the Iberian Peninsula determined from mitochondrial sequence analysis.

    Ann Hum Genet 1996, 60(Pt 4):331-350. PubMed Abstract OpenURL

  13. Macaulay VA, Richards MB, Hickey E, et al.: The emerging tree of West Eurasian mtDNAs: a synthesis of control-region sequences and RFLPs.

    Am J Hum Genet 1999, 64(1):232-249. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  14. Quintana-Murci L, Semino O, Bandelt HJ, Passarino G, et al.: Genetic evidence of an early exit of Homo sapiens sapiens from Africa through eastern Africa.

    Nat Genet 1999, 23(4):437-441. PubMed Abstract | Publisher Full Text OpenURL

  15. Passarino G, Semino O, Quintana-Murci L, Excoffier L, Hammer M, et al.: Different genetic components in the Ethiopian population, identified by mtDNA and Y-chromosome polymorphisms.

    Am J Hum Genet 1998, 62(2):420-434. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  16. Kivisild T, Reidla M, Metspalu E, Rosa A, Brehm A, et al.: Ethiopian mitochondrial DNA heritage: tracking gene flow across and around the gate of tears.

    Am J Hum Genet 2004, 75(5):752-770. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  17. Barton RNE, Bouzouggar A, Collcutt SN, Schwenninger J-L, Clark-Balzan L: OSL dating of the Aterian levels at Dar es-Soltan I (Rabat, Morocco) and implications for the dispersal of modern Homo sapiens.

    Quaternary Sci Rev 2009, 28:1914-1931. Publisher Full Text OpenURL

  18. Richter D, Moser J, Nami M, Eiwanger J, Mikdad A: New chronometric data from Ifri n’Ammar (Morocco) and the chronostratigraphy of the Middle Paleolithic in the Western Maghreb.

    J Hum Evol 2010, 59:672-679. PubMed Abstract | Publisher Full Text OpenURL

  19. McBurney C: The Haua Fteah (Cyrenaica) and the stone age of the south-east Mediterranean. Cambridge: Cambridge University Press; 1967. OpenURL

  20. Close AE: Backed bladelets are a foreign country.

    Archeol Papers Am Anthropol Assoc 2002, 12(1):31-44. OpenURL

  21. Barker G, Antoniadou A, Armitage S, Brooks I, Candy I, Connell K, Douka K, Drake N, Farr L, Hill E, et al.: The Cyrenaican Prehistory Project 2010: the fourth season of investigations of the Haua Fteah cave and its landscape, and further results from the 2007–2009 fieldwork.

    Libyan Stud 2010, 41:63-88. OpenURL

  22. Barich B, Garcea E: Ecological patterns in the upper Pleistocene and Holocene in the Jebel Gharbi, Northern Libya: chronology, climate and human occupation.

    Afr Archaeol Rev 2008, 25(1):87-97. Publisher Full Text OpenURL

  23. Merzoug S, Sari L: Re-examination of the Zone I material from Tamar Hat (Algeria): zooarchaeological and technofunctional analyses.

    Afr Archaeol Rev 2008, 25(1):57-73. Publisher Full Text OpenURL

  24. Barton N, Bouzouggar A, Humphrey L, Berridge P, Collcutt S, Gale R, Parfitt S, Parker A, Rhodes E, Schwenninger J-L: Human burial evidence from Hattab II cave and the question of continuity in late Pleistocene–Holocene mortuary practices in Northwest Africa.

    Camb Archaeol J 2008, 18(02):195-214. OpenURL

  25. Gasse F: Hydrological changes in the African tropics since the last glacial maximum.

    Quaternary Sci Rev 2000, 19(1–5):189-211. OpenURL

  26. Garcea EAA, Giraudi C: Late quaternary human settlement patterning in the Jebel Gharbi.

    J Hum Evol 2006, 51(4):411-421. PubMed Abstract | Publisher Full Text OpenURL

  27. Maca-Meyer N, Gonzalez AM, Pestano J, Flores C, et al.: Mitochondrial DNA transit between West Asia and North Africa inferred from U6 phylogeography.

    BMC Genet 2003, 4(1):15. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. Gonzalez AM, Larruga JM, Abu-Amero KK, Shi Y, Pestano J, Cabrera VM: Mitochondrial lineage M1 traces an early human backflow to Africa.

    BMC Genomics 2007, 8:223. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  29. Olivieri A, Achilli A, Pala M, Battaglia V, Fornarino S, et al.: The mtDNA legacy of the Levantine early Upper Palaeolithic in Africa.

    Science 2006, 314(5806):1767-1770. PubMed Abstract | Publisher Full Text OpenURL

  30. Pereira L, Silva NM, Franco-Duarte R, Fernandes V, Pereira JB, Costa MD, Martins H, Soares P, Behar DM, Richards MB, et al.: Population expansion in the North African late Pleistocene signalled by mitochondrial DNA haplogroup U6.

    BMC Evol Biol 2010, 10:390. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  31. Forster P: Ice Ages and the mitochondrial DNA chronology of human dispersals: a review.

    Philos Trans R Soc Lond B Biol Sci 2004, 359(1442):255-264. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Forster P, Romano V: Timing of a back-migration into Africa.

    Science 2007, 316(5821):50-53. PubMed Abstract | Publisher Full Text OpenURL

  33. Loogvali EL, Kivisild T, Margus T, Villems R: Explaining the imperfection of the molecular clock of hominid mitochondria.

    PLoS One 2009, 4(12):e8260. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  34. Soares P, Ermini L, Thomson N, Mormina M, Rito T, et al.: Correcting for purifying selection: an improved human mitochondrial molecular clock.

    Am J Hum Genet 2009, 84(6):740-759. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  35. Drummond AJ, Rambaut A, Shapiro B, Pybus OG: Bayesian coalescent inference of past population dynamics from molecular sequences.

    Mol Biol Evol 2005, 22(5):1185-1192. PubMed Abstract | Publisher Full Text OpenURL

  36. Pannell JR, Charlesworth B: Effects of metapopulation processes on measures of genetic diversity.

    Philos Trans R Soc Lond B Biol Sci 2000, 355(1404):1851-1864. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  37. Wakeley J, Aliacar N: Gene genealogies in a metapopulation.

    Genetics 2001, 159(2):893-905. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  38. Subramanian S: Temporal trails of natural selection in human mitogenomes.

    Mol Biol Evol 2009, 26(4):715-717. PubMed Abstract | Publisher Full Text OpenURL

  39. Kivisild T, Shen P, Wall DP, Do B, Sung R, et al.: The role of selection in the evolution of human mitochondrial genomes.

    Genetics 2006, 172(1):373-387. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  40. Peng MS, Quang HH, Dang KP, Trieu AV, Wang HW, Yao YG, Kong QP, Zhang YP: Tracing the Austronesian Footprint in Mainland Southeast Asia: A Perspective from Mitochondrial DNA.

    Mol Biol Evol 2010, 27(10):2417-2430.

    Oct

    PubMed Abstract | Publisher Full Text OpenURL

  41. Kong QP, Sun C, Wang HW, Zhao M, Wang WZ, Zhong L, Hao XD, Pan H, Wang SY, Cheng YT, et al.: Large-scale mtDNA screening reveals a surprising matrilineal complexity in east Asia and its implications to the peopling of the region.

    Mol Biol Evol 2011, 28(1):513-522. PubMed Abstract | Publisher Full Text OpenURL

  42. Bar-Yosef O: The upper Paleolithic revolution.

    Annu Rev Anthropol 2002, 31:363-393. Publisher Full Text OpenURL

  43. Iovita RP: Reevaluating Connections between the Early Upper Paleolithic of Northeast Africa and the Levant: Technological Differences between the Dabban and the Emiran. In Transitions in Prehistory: Essays in Honor of Ofer Bar-Yosef. Edited by Shea JJaL DE. Oxford: Oxboz; 2009:125-142. OpenURL

  44. Marks AE: Problems in Prehistory: North Africa and the Levant. Edited by Wendorf F, Marks AE. Dallas: Southern Methodist University Press; 1975:439. OpenURL

  45. Behar DM, Yunusbayev B, Metspalu M, Metspalu E, Rosset S, Parik J, Rootsi S, Chaubey G, Kutuev I, Yudkovsky G, et al.: The genome-wide structure of the Jewish people.

    Nature 2010, 466(7303):238-242. PubMed Abstract | Publisher Full Text OpenURL

  46. Diakonoff IM: Afrasian languages. Moscow: Nauka Press; 1988. OpenURL

  47. Ehret C, Keita SO, Newman P: The origins of Afroasiatic.

    Science 2004, 306(5702):1680.

    author reply 1680

    PubMed Abstract | Publisher Full Text OpenURL

  48. Militarev A: The Prehistory of a Dispersal: the Proto-Afrasian (Afroasiatic) Farming Lexicon. In Examining the Farming/Language Dispersal Hypothesis. Edited by Bellwood PRC. Cambridge: McDonald Institute for Archaeological Research; 2003:135-150. OpenURL

  49. Militarev A: The Jewish Conundrum in World History. Boston: Academic Studies Press; 2010. OpenURL

  50. Abu-Amero KK, Larruga JM, Cabrera VM, Gonzalez AM: Mitochondrial DNA structure in the Arabian Peninsula.

    BMC Evol Biol 2008, 8:45. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  51. Behar DM, Metspalu E, Kivisild T, Rosset S, Tzur S, Hadid Y, Yudkovsky G, Rosengarten D, Pereira L, Amorim A, et al.: Counting the founders: the matrilineal genetic ancestry of the Jewish Diaspora.

    PLoS One 2008, 3(4):e2062. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  52. Gonder MK, Mortensen HM, Reed FA, de Sousa A, Tishkoff SA: Whole-mtDNA Genome Sequence Analysis of Ancient African Lineages.

    Mol Biol Evol 2007, 24(3):757-768. PubMed Abstract | Publisher Full Text OpenURL

  53. Herrnstadt C, Elson JL, Fahy E, Preston G, Turnbull DM, et al.: Reduced-median-network analysis of complete mitochondrial DNA coding-region sequences for the major African, Asian, and European haplogroups.

    Am J Hum Genet 2002, 70(5):1152-1171. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  54. Kujanova M, Pereira L, Fernandes V, Pereira JB, Cerny V: Near eastern neolithic genetic input in a small oasis of the Egyptian Western Desert.

    Am J Phys Anthropol 2009, 140(2):336-346. PubMed Abstract | Publisher Full Text OpenURL

  55. Maca-Meyer N, Gonzбlez AM, Larruga JM, Flores C, et al.: Major genomic mitochondrial lineages delineate early human expansions.

    BMC Genet 2001, 2(1):13. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  56. Malyarchuk BA, Derenko M, Perkova M, Grzybowski T, Vanecek T, Lazur J: Reconstructing the phylogeny of African mitochondrial DNA lineages in Slavs.

    Eur J Hum Genet 2008, 16(9):1091-1096. PubMed Abstract | Publisher Full Text OpenURL

  57. Pereira L, Goncalves J, Franco-Duarte R, Silva J, Rocha T, et al.: No Evidence for an mtDNA role in sperm motility: data from complete sequencing of asthenozoospermic males.

    Mol Biol Evol 2007, 24(3):868-874. PubMed Abstract | Publisher Full Text OpenURL

  58. Behar DM, van Oven M, Rosset S, Metspalu M, Loogvali EL, et al.: A "Copernican" reassessment of the human mitochondrial DNA tree from its root.

    Am J Hum Genet 2012, 90(4):675-684. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  59. Pereira L, Cerny V, Cerezo M, Silva NM, Hajek M, Vasikova A, Kujanova M, Brdicka R, Salas A: Linking the sub-Saharan and West Eurasian gene pools: maternal and paternal heritage of the Tuareg nomads from the African Sahel.

    Eur J Hum Genet 2010, 18(8):915-923. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  60. Fraumene C, Belle EMS, Castri L, Sanna S, Mancosu G, et al.: High Resolution Analysis and Phylogenetic Network Construction Using Complete mtDNA Sequences in Sardinian Genetic Isolates.

    Mol Biol Evol 2006, 23(11):2101-2111. PubMed Abstract | Publisher Full Text OpenURL

  61. Hartmann A, Thieme M, Nanduri LK, Stempfl T, Moehle C, Kivisild T, Oefner PJ: Validation of microarray-based resequencing of 93 worldwide mitochondrial genomes.

    Hum Mutat 2009, 30(1):115-122. PubMed Abstract | Publisher Full Text OpenURL

  62. Hinttala R, Smeets R, Moilanen JS, Ugalde C, Uusimaa J, Smeitink JA, Majamaa K: Analysis of mitochondrial DNA sequences in patients with isolated or combined oxidative phosphorylation system deficiency.

    J Med Genet 2006, 43(11):881-886. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  63. Rieder MJ, Taylor SL, Tobe VO, Nickerson DA, et al.: Automating the identification of DNA variations using quality-based fluorescence re-sequencing: analysis of the human mitochondrial genome.

    Nucleic Acids Res 1998, 26(4):967-973. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  64. Bandelt HJ, Forster P, Rцhl A: Median-joining networks for inferring intraspecific phylogenies.

    Mol Biol Evol 1999, 16(1):37-48. PubMed Abstract | Publisher Full Text OpenURL

  65. Forster P, Harding R, Torroni A, Bandelt HJ: Origin and evolution of Native American mtDNA variation: a reappraisal.

    Am J Hum Genet 1996, 59(4):935-945. PubMed Abstract | PubMed Central Full Text OpenURL

  66. Saillard J, Forster P, Lynnerup N, Bandelt HJ, Nшrby , et al.: mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion.

    Am J Hum Genet 2000, 67(3):718-726. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  67. Cox MP: Accuracy of molecular dating with the rho statistic: deviations from coalescent expectations under a range of demographic models.

    Hum Biol 2008, 80(4):335-357. PubMed Abstract | Publisher Full Text OpenURL

  68. Yao YG, Salas A, Logan I, Bandelt HJ: mtDNA data mining in GenBank needs surveying.

    Am J Hum Genet 2009, 85(6):929-933.

    author reply 933

    PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  69. Excoffier L, Lischer HEL: Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows.

    Mol Ecol Resour 2010, 10(3):564-567. PubMed Abstract | Publisher Full Text OpenURL

  70. Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees.

    BMC Evol Biol 2007, 7:214. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  71. Drummond AJ, Ho SY, Phillips MJ, Rambaut A, et al.: Relaxed phylogenetics and dating with confidence.

    PLoS Biol 2006, 4(5):e88. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  72. Pereira L, Freitas F, Fernandes V, Pereira JB, Costa MD, Costa S, Maximo V, Macaulay V, Rocha R, Samuels DC: The diversity present in 5140 human mitochondrial genomes.

    Am J Hum Genet 2009, 84(5):628-640. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL