- Research article
- Open Access
Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa
BMC Evolutionary Biology volume 12, Article number: 234 (2012)
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.
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.
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.
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 . 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” , 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  or even as old as ~145 KYA . 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 . Whilst a techno-typological shift occurred within the Dabban ~33 KYA , 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 KYA at Haua Fteah with the Oranian ; the Iberomaurusian expansion in the Jebel Gharbi ~20 KYA ; and the first Iberomaurusian at Tamar Hat in Algeria ~20 KYA . The earliest Iberomaurusian sites in Morocco appear to be only slightly younger ~18 KYA . 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 . 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 ; 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  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 .
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.
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).
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 ) 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.
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 . 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 , hence the coalescent ages reported in this section may differ with the ones from the previous section. The rate by Soares et al. is applied here as the entire mtDNA genomic sequence is used for the BSP analyses, whereas the rate by Loogväli et al. applies only to the coding region. Nonetheless, the two different approaches offer similar estimates (see Table 2 and Additional file 1 for more detail).
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.
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),  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 . 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 . Marks  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 . 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 . 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.
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  , 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.
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.
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 , 1 from , 3 from , 8 from , 1 from , 2 from , 1 from , 2 from , 2 from , 51 from , 1 from , 3 from  and 3 from ); and 93 U6 (1 from , 6 samples from Family Tree DNA deposited in GenBank, 1 from , 2 from , 1 from , 12 from , 30 from , 2 from , 39 from  and 7 from ). 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 , 1 from  and 5 from  (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.
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 , 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).
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).
Phylogenic tree, network
Coalescent age estimates
For the coalescent age calculations, the rho (ρ) statistic and standard deviation were used as in [65, 66] but see  for a critical assessment of it. Different rates were used: For the coding region, rate  is used, and for the full genome, estimates were calculated with the calculator provided in . For all calculations, 2 M1 samples from  and the 3 M1 from  were discarded – the first ones were missing several portions of the coding region, and the second ones seemed to exhibit potential sequence errors (See  for details). For the full genome calculations, further samples were discarded, as their control region is not reported (M1: 1 from  and 2 from ; U6: 2 from , and 1 from ).
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 , 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 188.8.131.52 , with 100,000 permutations.
Bayesian Skyline Plot
The Bayesian Skyline Plot (BSP)  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 . The GTR substitution model was used with a gamma distribution, plus invariant sites. An uncorrelated lognormal relaxed clock was applied . 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 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 , 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  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.
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.
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. 10.1046/j.1469-1809.2003.00039.x.
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. 10.1111/j.1469-1809.2008.00493.x.
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. 10.1080/030144601300119106.
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. 10.1086/374384.
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. 10.1353/hub.2005.0028.
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. 10.1046/j.1529-8817.2004.00096.x.
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. 10.1093/molbev/msm277.
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. 10.1038/ejhg.2010.21.
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. 10.1002/ajpa.21472.
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-10.1371/journal.pgen.1002397.
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.
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. 10.1086/302204.
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. 10.1038/70550.
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. 10.1086/301702.
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. 10.1086/425161.
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. 10.1016/j.quascirev.2009.03.010.
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. 10.1016/j.jhevol.2010.07.024.
McBurney C: The Haua Fteah (Cyrenaica) and the stone age of the south-east Mediterranean. 1967, Cambridge: Cambridge University Press
Close AE: Backed bladelets are a foreign country. Archeol Papers Am Anthropol Assoc. 2002, 12 (1): 31-44.
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.
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. 10.1007/s10437-008-9020-6.
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. 10.1007/s10437-008-9028-y.
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.
Gasse F: Hydrological changes in the African tropics since the last glacial maximum. Quaternary Sci Rev. 2000, 19 (1–5): 189-211.
Garcea EAA, Giraudi C: Late quaternary human settlement patterning in the Jebel Gharbi. J Hum Evol. 2006, 51 (4): 411-421. 10.1016/j.jhevol.2006.05.002.
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-10.1186/1471-2156-4-15.
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-10.1186/1471-2164-8-223.
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. 10.1126/science.1135566.
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-10.1186/1471-2148-10-390.
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. 10.1098/rstb.2003.1394.
Forster P, Romano V: Timing of a back-migration into Africa. Science. 2007, 316 (5821): 50-53. 10.1126/science.316.5821.50.
Loogvali EL, Kivisild T, Margus T, Villems R: Explaining the imperfection of the molecular clock of hominid mitochondria. PLoS One. 2009, 4 (12): e8260-10.1371/journal.pone.0008260.
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. 10.1016/j.ajhg.2009.05.001.
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. 10.1093/molbev/msi103.
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. 10.1098/rstb.2000.0740.
Wakeley J, Aliacar N: Gene genealogies in a metapopulation. Genetics. 2001, 159 (2): 893-905.
Subramanian S: Temporal trails of natural selection in human mitogenomes. Mol Biol Evol. 2009, 26 (4): 715-717. 10.1093/molbev/msp005.
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.
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. 10.1093/molbev/msq131. Oct
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. 10.1093/molbev/msq219.
Bar-Yosef O: The upper Paleolithic revolution. Annu Rev Anthropol. 2002, 31: 363-393. 10.1146/annurev.anthro.31.040402.085416.
Iovita RP: Reevaluating Connections between the Early Upper Paleolithic of Northeast Africa and the Levant: Technological Differences between the Dabban and the Emiran. Transitions in Prehistory: Essays in Honor of Ofer Bar-Yosef. Edited by: Shea JJaL DE. 2009, Oxford: Oxboz, 125-142.
Marks AE: Problems in Prehistory: North Africa and the Levant. Edited by: Wendorf F, Marks AE. 1975, Dallas: Southern Methodist University Press, 439
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. 10.1038/nature09103.
Diakonoff IM: Afrasian languages. 1988, Moscow: Nauka Press
Ehret C, Keita SO, Newman P: The origins of Afroasiatic. Science. 2004, 306 (5702): 1680-author reply 1680
Militarev A: The Prehistory of a Dispersal: the Proto-Afrasian (Afroasiatic) Farming Lexicon. Examining the Farming/Language Dispersal Hypothesis. Edited by: Bellwood PRC. 2003, Cambridge: McDonald Institute for Archaeological Research, 135-150.
Militarev A: The Jewish Conundrum in World History. 2010, Boston: Academic Studies Press
Abu-Amero KK, Larruga JM, Cabrera VM, Gonzalez AM: Mitochondrial DNA structure in the Arabian Peninsula. BMC Evol Biol. 2008, 8: 45-10.1186/1471-2148-8-45.
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-10.1371/journal.pone.0002062.
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. 10.1093/molbev/msl209.
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. 10.1086/339933.
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. 10.1002/ajpa.21078.
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-10.1186/1471-2156-2-13.
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. 10.1038/ejhg.2008.70.
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. 10.1093/molbev/msm004.
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. 10.1016/j.ajhg.2012.03.002.
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. 10.1038/ejhg.2010.21.
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. 10.1093/molbev/msl084.
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. 10.1002/humu.20816.
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. 10.1136/jmg.2006.042168.
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. 10.1093/nar/26.4.967.
Bandelt HJ, Forster P, Rцhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16 (1): 37-48. 10.1093/oxfordjournals.molbev.a026036.
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.
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. 10.1086/303038.
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. 10.3378/1534-6617-80.4.335.
Yao YG, Salas A, Logan I, Bandelt HJ: mtDNA data mining in GenBank needs surveying. Am J Hum Genet. 2009, 85 (6): 929-933. 10.1016/j.ajhg.2009.10.023. author reply 933
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. 10.1111/j.1755-0998.2010.02847.x.
Drummond AJ, Rambaut A: BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007, 7: 214-10.1186/1471-2148-7-214.
Drummond AJ, Ho SY, Phillips MJ, Rambaut A, et al: Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006, 4 (5): e88-10.1371/journal.pbio.0040088.
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. 10.1016/j.ajhg.2009.04.013.
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.
The authors declare that they have no competing interests.
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.
Electronic supplementary material
Additional file 2: – Phylogenetic tree based on 105 M1 full sequences. All positions are scored against the RSRS  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). (XLSX 41 KB)
Additional file 3: – Network based on 236 M1 samples. All positions are scored against the RSRS  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). (XLSX 49 KB)
Additional file 4: – Phylogenetic tree based on 139 U6 full sequences. All positions are scored against the RSRS  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 ). 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). (XLSX 64 KB)
Additional file 5: – Network based on 230 U6 samples. All positions are scored against the RSRS  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. (XLSX 63 KB)
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. (PDF 665 KB)
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. (XLSX 38 KB)
Additional file 14: – BSP for M1 with the corrected rate versus uncorrected. The uncorrected rate use a rate of 1,695 x 10-8, and the corrected rate was deduced with the deduced rho values from the time, using the calculator from . (PDF 350 KB)
Additional file 15: – BSP for U6 with the corrected rate versus uncorrected. The uncorrected rate use a rate of 1,695 × 10-8, and the corrected rate was deduced with the deduced rho values from the time, using the calculator from . (PDF 347 KB)
About this article
Cite this article
Pennarun, E., Kivisild, T., Metspalu, E. et al. Divorcing the Late Upper Palaeolithic demographic histories of mtDNA haplogroups M1 and U6 in Africa. BMC Evol Biol 12, 234 (2012). https://doi.org/10.1186/1471-2148-12-234
- mtDNA haplogroups M1 and U6
- Afro-Asiatic languages
- North Africa