- Research article
Carriers of mitochondrial DNA macrohaplogroup L3 basal lineages migrated back to Africa from Asia around 70,000 years ago
BMC Evolutionary Biologyvolume 18, Article number: 98 (2018)
The main unequivocal conclusion after three decades of phylogeographic mtDNA studies is the African origin of all extant modern humans. In addition, a southern coastal route has been argued for to explain the Eurasian colonization of these African pioneers. Based on the age of macrohaplogroup L3, from which all maternal Eurasian and the majority of African lineages originated, the out-of-Africa event has been dated around 60-70 kya. On the opposite side, we have proposed a northern route through Central Asia across the Levant for that expansion and, consistent with the fossil record, we have dated it around 125 kya. To help bridge differences between the molecular and fossil record ages, in this article we assess the possibility that mtDNA macrohaplogroup L3 matured in Eurasia and returned to Africa as basal L3 lineages around 70 kya.
The coalescence ages of all Eurasian (M,N) and African (L3 ) lineages, both around 71 kya, are not significantly different. The oldest M and N Eurasian clades are found in southeastern Asia instead near of Africa as expected by the southern route hypothesis. The split of the Y-chromosome composite DE haplogroup is very similar to the age of mtDNA L3. An Eurasian origin and back migration to Africa has been proposed for the African Y-chromosome haplogroup E. Inside Africa, frequency distributions of maternal L3 and paternal E lineages are positively correlated. This correlation is not fully explained by geographic or ethnic affinities. This correlation rather seems to be the result of a joint and global replacement of the old autochthonous male and female African lineages by the new Eurasian incomers.
These results are congruent with a model proposing an out-of-Africa migration into Asia, following a northern route, of early anatomically modern humans carrying pre-L3 mtDNA lineages around 125 kya, subsequent diversification of pre-L3 into the basal lineages of L3, a return to Africa of Eurasian fully modern humans around 70 kya carrying the basal L3 lineages and the subsequent diversification of Eurasian-remaining L3 lineages into the M and N lineages in the outside-of-Africa context, and a second Eurasian global expansion by 60 kya, most probably, out of southeast Asia. Climatic conditions and the presence of Neanderthals and other hominins might have played significant roles in these human movements. Moreover, recent studies based on ancient DNA and whole-genome sequencing are also compatible with this hypothesis.
Based on molecular genetics analyses, the hypothesis of a recent African origin of modern humans that occurred approximately 200 thousand years age (kya) was formulated three decades ago . Today this hypothesis is widely accepted. There is also multidisciplinary agreement that the out-of-Africa expansion of modern humans promoted the extinction of other hominins in Eurasia with only a minor introgression of their genomes into modern human DNA . However, despite the enormous quantity of data accumulated to date, mainly from the analysis of mtDNA and Y-chromosome haploid markers, there is a lack of consensus regarding the time(s) of modern human dispersal out of Africa and route(s) followed. All of the indigenous mtDNA diversity outside Africa is comprised into clades M and N, which are derivative branches of the African haplogroup L3 [3,4,5]. This fact places the genetic time frame for the out of Africa dispersal at approximately 55-70 kya, which is the coalescence age of haplogroup L3 . Similarly, recent Y-chromosome sequence analysis detected a cluster of major non-African founder haplogroups originating within a short time interval at 47-52 kya . However, estimates based on a molecular clock depend on the mutation rate employed [8, 9]. These temporal windows for the exit of modern humans from Africa conflict with fossil, archaeological and ancient DNA data on the timing of the initial migration of anatomically modern humans (AMHs) out of Africa. Skeletal remains unearthed in the Skhul and Qafzeh caves demonstrated that early modern humans were present in the Levant between 125 and 80 kya . The discovery of modern human teeth in southern China dated to 120-80 kya , also supports the presence of AMHs in eastern Asia during this period. Several archaeological studies uncovered Middle Stone Age (MSA) lithic assemblages, dated at approximately 125-75 kya, in different regions of the Arabian Peninsula, presenting affinities with northeastern African assemblages of the same period [12,13,14]. These findings suggest that African AMHS may have extended their geographic range to eastern and northern Arabia, as well as south Asia, as early as 125 kya, long before the time frame of the migration that shaped the modern global genetic pool as suggested by molecular data.
Ancient DNA (aDNA) analysis is an important tool in the reconstruction of past human history. Based on aDNA analyses, the Neanderthal introgression into modern humans in Europe has been dated to within 35-65 kya , which is well within the molecular clock-based time frame established for the African exit of modern humans. However, an ancient gene flow from early modern humans into the ancestors of eastern Neanderthals more than 100 kya was recently reported . These data evidenced that early modern humans and ancestors of Neanderthals from the Siberian Altai region interbred much earlier than previously thought. Furthermore, whole-genome based studies place the split of Eurasian from African populations at 88-112 kya , and the presence of AMHS out of Africa has been documented before 75 kya . A way to reconcile these contradictory pieces of evidence is to state that all these ancient movements out of Africa, prior to 70 kya, did not significantly contribute genetically to present-day human populations . However, major efforts should be dedicated to resolving the conflicting evidence.
Concerning the potential routes followed by modern humans out of Africa, there are two, non-mutually exclusive main alternatives: a northern dispersal along the Nile-Sinai corridor, and a southern dispersal from the Horn of Africa across the Bab al Mandeb Strait. In early mtDNA phylogeographic studies, the virtual absence of mtDNA haplogroup M in the Levant, and its presence in Ethiopia, southern Arabia, the Indian subcontinent and East Asia, rendered M the first genetic indicator of a southern-route exit from eastern Africa . Shortly after these studies, based on the rarity of mtDNA haplogroup N(xR) in India, and its continuous presence above the Himalayas, we proposed an additional northern route through the Levant . Since that time, intensive and extensive research on mtDNA has been carried out, on populations not only from Central [20,21,22,23] and East Asia [24,25,26] but also from the regions along the hypothetical southern route, such as India [27,28,29,30], mainland southeastern Asia [31,32,33], island southeast Asia [34,35,36,37,38,39], New Guinea, North Island, Melanesia and Australia [40,41,42,43,44]. The most unexpected results were that some N haplogroups in southern China (N10, N11) were older than the oldest N western Asia lineages (N1, N2), and that some M haplogroups in Melanesia (M27, Q) were older than the oldest Indian M lineages (M2, M33). Different researchers have provided conflicting interpretations of these results. Some perceived them as confirming a rapid, southern coastal spread of modern humans from Africa [30, 45,46,47]. Others postulate ancient local population differentiation in each region without any evidence of the shared ancestry expected by the southern dispersal model [48,49,50]. The existence of a northern route, deduced from the phylogeography of macrohaplogroup N , has received additional support from the fossil record , whole-genome studies comparing Egyptian and Ethiopian populations , and the fact that all non-African populations present a signal of Neanderthal introgression . However, we realized that what actually macrohaplogroup N suggests is a human movement from southeastern Asia to western Asia . We observed the same tendency for macrohaplogroup M, in this case, expanding westwards to India . A similar trend was observed for macrohaplogroup R, the main sister branch of N . Thus, we confirmed that macrohaplogroup M and N indicated, major southern and northern expansions, respectively, of modern humans but, in the opposite sense we had predicted previously . Studies based on Y-chromosome sequences also suggested southeastern Asia as an early center of human expansions [56,57,58,59]. If one accepts that basal L3 lineages (M, N) evolved independently in southeastern Asia and not in Africa or near the borders of the African continent where the remaining L3 lineages expanded, one is confronted with the question of where the basal trunk of L3 evolved. A gravitating midpoint between eastern Africa and southeastern Asia would situate the origin of L3 in inner Asia, with possible opposite-direction expansions back to Africa and forward to eastern Asia. This possibility has been modeled, along with other options, by others obtaining the highest likelihood value among competing models ; however, it has not received the attention it deserves. The parallelism of this early back-to-Africa scenario of mtDNA haplogroup L3 with that proposed for the Y-chromosome haplogroup E  is striking.
In this work, we have improved the phylogeny of some African mtDNA L3 subclades. In addition, we have compared the frequency distribution of the younger African mtDNA haplogroup L3 and Y-chromosome haplogroup E in the main regions of the African continent. With these data at hand, we assess the possibility of the following scenario: L3 exited from Africa as a pre-L3 lineage that evolved as basal L3 in inner Asia. From there, it expanded, returning to Africa as well as expanding to southeastern Asia, giving rise to the African L3 branches in eastern Africa and the M and N L3 Eurasian branches in southeastern Asia, respectively. This model, which implies an earlier exit of modern humans out of Africa, has been tested against independent results from other disciplines.
A total of 69 complete mtDNA genomes were sequenced in this study (Additional file 1: Table S1). They comprise the main African L haplogroups, excepting L6. To remedy this lack, 12 previously published and complete L6 sequences were included in our phylogenetic tree (Additional file 2: Figure S1). The different branches of haplogroup L3 are represented by 45 of these sequences. To establish the relative frequency of mtDNA macrohaplogroup L3 in the main African regions, a total of 25,203 partial and total publicly available mtDNA sequences were screened. Of them, 1,138 represent our unpublished data (Additional file 1: Table S2). To establish the relative frequency of Y-chromosome macrohaplogroup E in the main African regions, a total of 21,286 publicly available Y-chromosome African samples were screened. Of them, 737 represent our unpublished data (Additional file 1: Table S2). All of our samples were collected in the Canary Islands or Saudi Arabia from academic and health-care centers. The procedure of human population sampling adhered to the tenets of the Declaration of Helsinki, and written consent was obtained from all participants before they participated in the study. The study underwent a formal review and was approved by the College of Medicine Ethical Committee of the King Saud University (proposal N° 09-659) and by the Ethics Committee for Human Research at the University of La Laguna (proposal NR157).
Total DNA was isolated from buccal or blood samples using the POREGENE DNA isolation kit from Gentra Systems (Minneapolis, USA). The PCR conditions and procedures for sequencing the mtDNA genome were as previously published . Successfully amplified products were sequenced for both complementary strands using the DYEnamic™ETDye terminator kit (Amersham Biosciences). Samples were run on MegaBACE™ 1000 sequencer (Amersham Biosciences) according to the manufacturer’s protocol. The 69 new complete mtDNA sequences have been deposited in GenBank under the accession numbers MF621062 to MF621130 (Additional file 1: Table S1).
Compilation of Previously published data
Sequences belonging to specific mtDNA L haplogroups were obtained from public databases such as NCBI, MITOMAP, the-1000 Genomes Project and the literature. We searched for mtDNA lineages directly by using diagnostic SNPs (http://www.mitomap.org/foswiki/bin/view/MITOMAP/WebHome), or by submitting short fragments including the diagnostic SNPs to a BLAST search (http://blast.st-va.ncbi.nlm.nih.gov/Blast.cgi). Haplotypes extracted from the literature were transformed into sequences using the HaploSearch program (http://www.haplosite.com/haplosearch/process/) . Sequences were manually aligned and compared to the rCRS  with the BioEdit Sequence Alignment program . Haplogroup assignment was performed by hand, with screening for diagnostic positions or motifs at both hyper-variable and coding regions performed whenever possible. Sequence alignment and haplogroup assignment were carried out twice by two independent researchers and any discrepancy was resolved according to the PhyloTree database (Build 17; http://www.phylotree.org/) . For the screening of the Y-chromosome haplogroup E, we considered samples as belonging to this haplogroup if they were found positive for, at least, the diagnostic DE-YAP or E-M40, E-M96 markers.
The phylogenetic tree was constructed by means of the Network program, v188.8.131.52, using, the Reduced Median algorithm, Median Joining algorithm and Steiner (MP) algorithm in sequential order . Remaining reticulations were manually resolved. Haplogroup branches were named following the nomenclature proposed by the PhyloTree database . Our coalescence ages were estimated by using statistics rho  and Sigma  and the calibration rate proposed by Soares et al. .
To calculate the total mean age of each haplogroup, we compiled all of the different estimation ages from the literature without taking into account the mtDNA sequence segment analyzed, the mutation rate considered, or the observed partial overlapping of the samples used. In cases where the same sample set was used to calculate its age by using multiple methods, we selected the age calculated based on the rho statistic as the most generalized method (Additional file 1: Table S3). To calculate haplogroup mean coalescence ages for the non-recombining region of the Y-chromosome (NRY), we compiled estimations preferably based on single nucleotide polymorphisms (SNPs) obtained by sequencing. When different mutation rates were used in the same study, we chose the age calculated based on the slowest mutation rate (Additional file 1: Table S4).
In this study, we focus on the earliest periods of the out-of-Africa spread of modern humans and the likely return to Africa of the carriers of primary mtDNA L3 and Y-chromosome E lineages.
As the phylogeography of the different branches of these lineages has been extensively studied by other authors to reveal more recent human movements on the continent, we focus here on the continental distributions in the major African regions. For phylogeographic purposes, we divided the African continent into the following eight major regions: 1. Northwest Africa (including Morocco, West Sahara, Algeria and Tunisia), 2. Northeast Africa (including Libya and Egypt), 3. West Sahel (including Mauritania, Mali, and Niger), 4. East Sahel (including Chad, Sudan, Ethiopia, Somalia, and Eritrea), 5. West Guinea (including Senegambia, Guinea-Bissau, Guinea-Conakry, Sierra-Leona, Liberia, Ivory-Coast, Burkina-Faso, Ghana, Togo, Benin, and Nigeria), 6. Central Africa (including Cameroon, Central African Republic (CAR), Congo Democratic Republic (CDR), Congo-Brazzaville, Gabon and Equatorial Guinea), 7. East Guinea (including Uganda, Rwanda, Kenya, and Tanzania), and 8. Southern Africa (including Angola, Zambia, Malawi, Mozambique, Zimbabwe, Botswana, Namibia and South African Republic (SAR)).
To evaluate the level of geographic structure in the mtDNA macrohaplogroup L3 and the Y-chromosome macrohaplogroup E in Africa, we performed AMOVA and K-means clustering analyses. We used the GenAlEx6.5 software to implement AMOVA and XLSTAT statistical software to perform the K-means clustering analysis. The possible associations between the frequencies of mtDNA macrohaplogroup L3 and those of the Y-chromosome macrohaplogroup E, for both the whole African continent and each of its principal geographic subdivisions were tested by Pearson correlation analyses using XLSTAT software. As an extensive overlap exists among the expansion ages of the L3 branches with those of the widespread African mtDNA macrohaplogroup L2, the global frequencies of L2 were included in the majority of the phylogeographic analyses performed.
Results and discussion
Phylogeny and affinities of our African complete sequences
Overall, our 69 mtDNA complete sequences (Additional file 2: Figure S1) could be allocated into previously defined clades in the PhyloTree database. Their closest affinities were with other sequences of the same haplogroups. Thus, our Kenyan L3a1a (Kn028) sequence shares tip mutations 514, 3796 and 4733 with a Tanzanian sequence (EF184630) but only mutation 514 with a Somalian sequence (JN655813) of the same clade. The Sudanese L3b1a (Su238) sequence shares the very conservative transition at 12557 with an L3b sequence (KF055324) from an African-American glaucoma patient . Our L3b1a2 (Su002) sequence has matches at 195, 12490 and 16311 with several African sequences (EU092669, EU092744, EU092795, EU092825, EU9355449) with which it composes a new branch, L3b1a2a, defined by these three transitions. Similarly, the L3f2a1 (Su004) sequence has matches at mutated positions 6182, 8676, 9731, 12280, 12354 and 13105 with other published Senegalese sequences (JN655832, JN655841) with which composes a new derived branch. We expected L sequences detected in the Canary Islands to have their closest relatives among sequences from the African continent. This was observed in some cases; for example, the L3d1b3 (Go764) sequence from La Gomera island shares tip transitions 14040 and 16256 with an Ovimbundu isolate (KJ185837) from Angola . However, unexpectedly, the Canarian sequence TF0005, allocated to the L3f1b subclade, has its closest relatives in the Iberian Peninsula, sharing the 8994 transition with two Asturian L3f1b sequences (KJ959229, KJ959230) . Furthermore, the L3x2 (TF116) sequence from Tenerife shares all of its terminal variants(650, 7933, 8158, 15519, 16261) with sequences from Galicia (HQ675033, JN214446) and Andalusia (KT819228), not with African sequences. Saudi Arabia has been identified as an important receptor of mtDNA Eurasian lineages, as well as those of African origin. Arab sequences belonging to the L3i1a (AR429) and L3x1a1 (AR260) haplogroups have their closest relatives with sequences JN655780 and DQ341067, respectively, from nearby Ethiopia, and the L3h1b1 (AR381) sequence is identical to a previously published Yemeni isolate (KM986547). However, the L3h1b2 (AR221) sequence is most related to the JQ044990 lineage from Burkina Faso , with which it shares particular transitions at positions 7424, 13194, 16192 and 16218. The affinities of the Arab L1c2b1a'b (AR1252) with other sequences are the most unexpected. This sequence, particularly characterized by the presence of an insertion of 11 nucleotides at the 16029 position in the control region, has an exact match with an L1 isolate from the Dominican Republic (DQ341059). Its closest relatives in Africa, while lacking the above-mentioned insertion, are found in Angola (KJ185814) and Zambia (KJ185662) among Bantu-speakers . The control region of this AR1252 isolate was previously published (KP960821). Concerning the less frequent L4, L5, and L6 clades, our L4b1a (Iv136) sequence from the Ivory Coast shares tip mutations 789, 7166 and 14935 with geographically nearby sequences (JQ044848, JQ045081) from Burkina Faso . Similarly, the Arab L4a2 (AR1116) sequence is closely related to other African L4a2 sequences (EU092799, EU092800), and the L4b2a1 (AR197) isolate is identical to a sequence (KM986608) from Yemen . From the analysis of partial sequences [53, 74], we can be certain that representatives of branches L4a1, L4a2 and L4b2 exist in Saudi Arabia. However, we have not yet detected sequences belonging to the large Sudanese L4b1b clade (Additional file 2: Figure S1). Published [53, 74] and unpublished data allow us to confirm that L5a is represented in Saudi Arabia by at least one lineage that has the following haplotype in/nearby the coding region: 15884, 16093, 16129, 16148, 16166, 16187, 16189, 16223, 16265C, 16311, 16355, 16362/ 73, 152, 182, 195, 198, 247, 263, 315iC, 455i2T, 459iC, 513, 522dCA, 709, 750, 769, 825A, 851, and 930. It represents 0.58% of the Saudi mtDNA gene pool. There is also a L5b lineage characterized by mutations at 15927, 16111, 16129, 16148, 16166, 16187, 16189, 16223, 16233, 16254, 16265C, 16278, 16360, 16519/ 73, 195, 247, 249d, 263, 315iC, 459iC, 501, 535, 750, 769, and 825A, with minor presence (0.09%) in Saudi Arabia. Furthermore, although we lack complete L6 sequences, we can confirm based on partial sequencing and specific SNP analysis, that L6a1 (0.13%) and L6b (0.09%) lineages are also present in the Saudi population [53, 74]. With 12 complete sequences, L3h is the best-represented haplogroup in our phylogeny (Additional file 2: Figure S1). In it, we have provisionally defined some new branches as follows: The retromutation at position 16223 defines a Sudanese L3h1a1a branch. Two clades, L3h1a2a1a and L3h1a2a1b are characterized by transitions 3892, 7705, and 15346 and transitions 5108 and 16165, respectively. An additional subclade, defined by transitions 7310, 13153, 14407 and transversion 9824A is provisionally named L3h1a2b1. In addition, after introducing the AR221 sequence, the old branch L3h1b2 is characterized only by transitions 294, 8842, 9758, 12882, 13437, 16129, and 16362.
Our only discrepancy with the PhyloTree phylogeny, regards the rare and old L5 clade. We have identified a new branch, provisionally named L5c (Additional file 2: Figure S1). In light of the information provided by this lineage, the PhyloTree L5b node, that joints haplogroups L5b1 and L5b2 by sharing retromutations at positions 182, 13105 and 16311 and a transition at position 16254 appears to lack phylogenetic robustness. The PhyloTree L5b2 clade is more appropriatelly considered a sister branch of our new Sudanese L5c sequence (Su412), which share a retromutation at the 195 position and transitions at the 6527 and 11809 positions.
Despite the small sample sizes employed, our coalescence age results fall within the standard deviations calculated for the different haplogroups (Additional file 1: Table S3). However, the mean age value for the L3 macrohaplogroup in Africa (71 ± 12 kya), which theoretically marks the upper-bound time for the out-of-Africa movement of modern humans, falls short compared to those estimated from the fossil record in the Levant .
The Eurasian origin of mtDNA macrohaplogroup L3
The southern route hypothesis proposes that the Eurasian branches (M and N) of macrohaplogroup L3 differentiated in or near the African continent and rapidly spread across the Asian peninsulas to reach Australia and Melanesia . Under this hypothesis, it is expected that, the coalescence ages of haplogroups should generally decrease from Africa to Australia. However, we have demonstrated that this is not the case [53,54,55]. Just on the contrary, the oldest M and N haplogroups are from southern China and Australasia, not India, and the associations between longitudinal geographic distances and relative ages of the M and N haplogroups run, contrary to expectation, from east to west [53, 54]. This presents us with a dilemma: It appears that two gravity centers of L3 expansion exist, one in Africa and one in southeastern Asia. A geographic equidistant midpoint would situate the primary radiation of L3 in India if a southern route were followed by the African colonizers and above the Himalayas, between Tibet and Pamir, if the northern route was followed. Furthermore, as the coalescence age of the African L3 branches and that of the Eurasian L3 (MN) branches are very similar (Table 1), at approximately 71 kya, the temporal and spatial midpoints might also coincide. As the group of modern humans that hypothetically returned to Africa is expected to include both females and males, searching for Y-chromosome phylogenetic and phylogeographic information might provide additional information.
The Eurasian origin of Y-chromosome haplogroup E
An origin in Asia and return to Africa was proposed, long ago, for the Y-chromosome African haplogroup E . This hypothesis was based on the derived state of its African YAP+ haplotypes 4 and 5 (haplogroup E) with respect to the ancestral Asian YAP+ haplotype 3 (haplogroup D). The later discovery of new markers evidenced that D and E were sister branches of the YAP+ node. Haplogroup D showed the derived status for M174 and the ancestral status for M40 in Asia, whereas haplogroup E was characterized by the derived status for M40 and the M174 ancestral status throughout Africa, thus the migratory sense between continents of both haplogroups could not be assured . A few YAP+ individuals, ancestral for both markers, were detected in West Africa  and in Tibet . Although assigned to the para-haplogroup DE*, its actual ancestral state could not be confirmed. Furthermore, a new mutation (P143) united the two other Eurasian haplogroups C and F as brothers and, in turn, DE and CF were joined to an ancestral node defined by mutations M168 and M294 . The initial scenario proposed to explain this complicated situation was that two independent migrations out of Africa occurred: one marked by D, and the other marked by the CF pair of lineages . However, a new interpretation arose following the discovery of more than 60,000 single nucleotide variants by next generation sequencing techniques. The most parsimonious interpretation of the Y-chromosome phylogeny constructed with these variants is that the predominant African haplogroup E arose outside the continent and back-migrated to Africa . The DE split as a lower bound (69.0 ± 14.7 kya) and the radiation of E into Africa as an upper bound (65.5 ± 8.5 kya) yields a time frame highly coincidental with the dates estimated for the mtDNA haplogroup L3 expansions (Table 1). Furthermore, the spatial distribution of the residual Y-chromosome haplogroup D in Asia is a good indicator of the geographical location of the putative DE split. The highest frequency and diversity of D are in the Tibet region. Although D is also present at low incidence throughout southeast Asia, two other centers with notable frequency are in Japan and the Andaman Islands [78, 81, 82], possibly representing edge relic areas of a wider distribution long ago. There are no native D lineages in India, decreasing the possibility that this subcontinent was the center of the DE partition, which argues against a southern route. It is likely that the divergence of the Y-chromosome D and E haplogroups occurred up to the Himalayas and in or westward of the Tibet which coincides with the hypothetical bifurcation center proposed for the mtDNA L3 macrohaplogroup.
Coetaneous east and westward expansions of modern humans from Central Asia
The above described coincidental divergences of female and male lineages occurred during a glaciation period (70 - 100 Kya), it is possible that cold climatic conditions forced humans southwards and that, upon being confronted with the Himalayas, humans dispersed across southeastern and southwestern Asia. It is likely that this climatic change also led the Neanderthals to extend their southern range and, consequently increase their geographic overlap with humans and, possibly, Denisovans, outcompeting them in the search for resources (Figure 1). This southward retreat was more extensive in the west, as evidenced by the total occupation of the Levant by the Neanderthals approximately 70 kya  and the forced return of modern humans carrying the mtDNA L3 and Y-chromosome E basal lineages to Africa. However, 20 ky later the situation changed. Modern humans advanced westwards from inner Asia, displacing Neanderthals along the way and colonizing East Asia, South Asia, and Central Asia, from where they reached the Levant approximately 50 kya  and Europe shortly thereafter [85,86,87]. This westward colonization of modern human was also proposed based on an archaeological perspective [87, 88]. Under the proposed model, early modern humans would have left Africa much earlier than the time frame proposed by geneticists under mitochondrial molecular clock restrictions . In an attempt to resolve this discrepancy, one could fix this period at the L3'4 or L3'4'6 mtDNA coalescence node, i.e., at approximately 80 to 100 kya (Table 1), such that mutations 769, 1018 and 16311, which define the basal L3* lineage, would have arisen after the movement out of Africa. Accordingly, the exit of the companion men could be dated at the split of branch CDEF-M168 from B-M181 approximately 86-120 kya [59, 90]. However, given the inaccuracies of the molecular clock, we prefer to rely on fossil and climatic records to establish the timing of the out-of-Africa movement of early modern humans across the Levant, which yields approximately 125 kya as the most likely period.
The phylogeography of the L3 and E lineages inside Africa
The continental mean frequencies of the mtDNA and NRY haplogroups in the six main regions of the African continent are presented in Table 2. The mtDNA haplogroup L3 is more frequent in sub-Saharan Africa than in North Africa or the Sahel. In contrast, the Eurasian mtDNA haplogroups, including M1 and U6, are more frequent in northern Africa and the Sahel than in sub-Saharan Africa. In general, the Y-chromosome haplogroup E is more frequent in western than eastern Africa, while the Eurasian Y haplogroups show the opposite trend. This geographic distribution confirms that the history of Africa is marked by multiple Eurasian migratory waves that brought the first carriers of female haplogroup L3* and male haplogroup E* basal lineages into Africa. It has been suggested that the few sub-Saharan haplogroups present in northern Africa are the result of recent historical incorporations . Studies of ancient DNA in the area appear to confirm this scenario both in the northwest  and the northeast . The fact that L3k, the only autochthonous L3 lineage in northern Africa, has only a residual presence in the area also supports this scenario . Under this supposition, male E lineages, currently present in North Africa, would have reached the area as a secondary wave accompanying Eurasian female lineages as M1 and U6. The main indigenous E clades in the region, E-M81 in the northwest and E-M78 in the northeast, are derived from haplogroup E-M35, which also includes the European Mediterranean and western Asian branches E-V13 and E-V22 . In conflict with earlier studies that considered a Paleolithic implantation of E-M81 in the Maghreb [96, 97], later work suggested that the low microsatellite diversity of this clade in northwest Africa is better explained as the result of Neolithic or post-Neolithic gene flow episodes from the Near East . However, the subsequent discovery of a new sister branch of E-M81, named E-V257 , without Near Eastern roots but present in the European western and central Mediterranean shores and in Cameroon and Kenyan populations [99, 100], weakened support for the suggested Levantine connection. Furthermore, E-M81 and E-M78 precursors are very old lineages that accumulated 23 and 16 mutations, respectively, in their basal branches. It has been reported that E-M78 radiated in eastern Africa between 20 and 15 kya while E-M81 did not, most likely because it was already in the Maghreb at that time. This timing coincides with the expansion age of the mtDNA U6a haplogroup in the area . Thus, a recent re-expansion after a large bottleneck is the best explanation for the low variance of E-M81 in present day . The persistence of an even older male demographic substrate in this area was evidenced by the detection of representatives of the deepest Y phylogenetic clades A0 and A1a in the region . There is general consensus in attributing an eastern African origin to the initial expansion of the NRY haplogroup E in the continent . Curiously, ancestral E* lineages have been detected in the Arabian Peninsula  and the Levant . Regardless of its origin, haplogroup E shows lower frequencies in northeastern Africa and the eastern Sahel than in the west. The opposite pattern is observed in the frequencies of Y Eurasian haplogroups in the same areas (Table 2), which suggests later stronger Eurasian male gene flow throughout the northeastern side.
Detailed frequencies for mtDNA haplogroups L2 and L3 and Y-chromosome haplogroup E, all around the African continent, are listed in Additional file 1: Table S2. We have included L2 in the analysis because it is the sister clade of the composite eastern African node L3'4'6 that, through consecutive range expansions, promoted the exit of the L3 precursor out of Africa. In addition, within Africa, several L2-derived spreads coincide in time with the later expansion of L3 branches in the continent. Furthermore, there is a suggestive positive association between the mean frequency estimates for L2 and the Y-chromosome haplogroup E across the major African regions (Table 2). In fact, there are significant positive correlations between E haplogroup frequency and both L3 (r = 0.400; p < 0.0001) and L2 (r = 0.347; p < 0.0001) mtDNA haplogroup frequency across Africa. An even stronger correlation is observed between E haplogroup frequency and the combined frequency of the two L2 and L3 haplogroup (r = 0.477; p < 0.0001). However, the strength of the association varies among the different regions considered, with no significant correlation in northern Africa, as can be expected from the discussion above. The correlation is barely significant in the Sahel region (r = 0.246; p = 0.045), but is highly significant in the remaining regions, particularly in southern Africa (r= 0.615; p < 0.0001). Nevertheless, these correlations are only slightly associated with geography, as indicated by the AMOVA results showing that only 4% of the variance is explained by differences among regions (Table 3). It appears that the E and L2/L3 expansions were strongest in western Sahel and western Guinea, where they replaced the majority of the oldest mtDNA (L0 and L1) and Y-chromosome (A and B) lineages (Table 2). We applied the k-means clustering algorithm to our L3 and E frequency data (Additional file 1: Table S2). The consecutive partitioning of the samples into clusters functions to minimize the variance within groups while augmenting the variance among them (Table 3). At k = 5, less than 20% of the variance is due to differences within clusters. At this value of k, the five classes obtained have centroid means for L3 and E that minimize the mean-square distance of the samples grouped to this center (Table 4 and Additional file 1: Tables S5, S6 and S7). Class I, characterized by relatively low frequencies for both the L3 and E haplogroups, encompasses the majority of the Khoesan-speaking groups from South Africa, Namibia, and Angola and the Hadza from Tanzania, as well as several pygmy groups from Cameroon, Gabon, CAR and Congo, such as the Baka and the Babinga. In addition to having different geographic locations, these groups are differentiated by the frequencies of other haplogroups. Thus, samples from Khoesan-speaking groups harbor high frequencies of mtDNA L0d and L0k haplogroups and Y-chromosome A lineages, while the Central African pygmies are characterized by the highest frequencies of mtDNA L1c and Y-chromosome B-M60 lineages. The Hadza share with pygmy groups high percentages of B-M60 chromosomes and the highest frequency of mtDNA haplogroup L4 among the entire African continent. Other groups that belong to this class are the majority of Nilotes from Sudan and Uganda, characterized by their high frequencies of mtDNA haplogroup L2a1, and samples from several Afro-Asiatic-speaking groups from Egypt and Sudan which also show high levels of L2a1 lineages and have many mtDNA L0a1 and Y-chromosome B representatives. The remaining of the samples from Khoesan-speaking groups fall into Class II having a low frequency of L3 and an intermediate frequency of E (Table 4), suggesting male-biased gene flow from western sub-Saharan Africans. This class also contains the remaining central Africa pygmies, including Sanga, Mbenzele, Biaka, and Mbuti, which all harbor high frequencies of Y B-M60 and L1c lineages. This class includes mainly northeast African and Ethiopian samples of individuals speaking Afro-Asiatic languages; some Nilo-Saharan-speaking groups such as the Fur from Sudan, the Anuak from Ethiopia and the Maasai from Kenya; and Bantu-speakers such as the Shona from Zimbabwe and Botswana and the Tswana from Botswana and South Africa. However, most of the northwestern African Afro-Asiatic speaking groups fall into Class III, defined by low frequencies of L3 and high frequencies of E (Table 4). A general low frequency of mtDNA haplogroup L2 is also a characteristic of this class. The majority of Berber and Tuareg samples belong to this class, including the Gossi, Tamashek, and Douentza from Mali and the Gorom from Burkina Faso. Interestingly, the click-speaking Sandawe and the Nilo-Saharan speaking Datog from Tanzania also belong to Class III. On the maternal side, these Tanzanian samples are characterized by relatively high frequencies of haplogroups L0a and L4. However, the most abundant components of this class are samples from Niger-Congo speakers, including the majority of the Senegalese samples as well as southern African-specific-Bantu-speakers such as the Zulu and Xhosa. Curiously, the click-speaking Xeg Khoesan and Khwe belong to this cluster, indicating substantial gene flow from Bantu-speaking immigrants. Such flow was reported long ago for the Khwe, who were found to be more closely linked to non-Khoesan-speaking Bantu populations . Another class dominated by Niger-Congo speakers is Class IV, the largest class. It includes samples that possess intermediate frequencies of L3 and high frequencies of E (Table 4). Linia and Kanembou from Chad, Rimaibe from Burkina-Faso, Songhai from Nigeria, and Masalit from Sudan are the only Nilo-Saharan speakers in this class. All the western African countries are represented by different Niger-Congo speaking groups, including the Bateke pygmies from Congo. There are also instances of eastern or southern African Bantu representatives. Finally, Class V includes those samples that present high frequencies for both the L3 and E haplogroups (Table 4). Again, Niger-Congo speakers make up the majority of this class, although Nilo-Saharan from western countries such as Menaka of Mali, Diffa of Niger and Kanuri from Nigeria and those of eastern countries such as Bongor from Chad and Luo from Kenya are included in this class. The two best-represented countries are the western sub-Saharan Africa Nigeria and Cameroon, which provided most of the samples of Niger-Congo speaking individuals. Nevertheless, there are Bantu-specific speakers from Kenya and southern Africa in this class. This group includes the click-speakers Damara from Namibia and South Africa, who genetically have been associated with Bantu-speaking groups rather than to other Khoesan-speaking groups .
The above results show that the positive correlation found between Y-chromosome haplogroup E and mtDNA haplogroup L3 (and L2) lineages is not strongly associated with either geography nor language. It is better explained as the result of a gradual substitution of the most basal mtDNA (L0, L1, L5) and Y-chromosome (A, B) lineages by the phylogenetically younger clades L2/L3, and E, respectively, throughout Africa. The data also suggest important sex-biased dispersal between populations. These evident gene replacements in Africa have been mainly attributed to recent geographic range expansions of pastoralist and agriculturalist populations from eastern and western Africa at the expense of the hunter-gatherers inhabitants of the Central Africa rainforest [108,109,110], eastern African forested areas around the Great Lakes [111,112,113,114], and the semi-desert open spaces of South Africa [115,116,117,118]. Under our hypothesis of an early return to Africa from Eurasia of basal mtDNA L3 and Y-chromosome E lineages and their expansion approximately 70 kya, into East and subsequently West Africa, these lineage replacements must have begun very early. It appears that in this first expansion, mtDNA haplogroup L2 was incorporated by female assimilation, whereas their hypothetical Y-chromosome haplogroup B counterparts were outcompeted by the incoming E chromosomes. An ancient expansion from a Central African source into eastern Africa at 70-50 kya has been associated with haplogroup L2 . Similarly, an early expansion within Africa 60-80 kya involving L3 and possibly L2 was detected long ago . This latter expansion was considered the crucial event in the exit of modern humans from Africa into Eurasia. However, our proposition is that it signaled a backflow from Eurasia and subsequent expansion into Africa.
Moreover, our hypothesis of an early return and subsequent expansion inside Africa of carriers of L3 and E haplogroups might help explain, the Neanderthal introgression detected in the western African Yoruba [121, 122] and in northern African Tunisian Berbers .
A new mtDNA model of the origin and dispersal of Homo sapiens
At mtDNA level, data accumulated during the last thirty years, including those contributed by studies of ancient DNA, allow us to propose a more detailed model of the origin and worldwide spread of modern humans than the ones proposed three decades ago. There are three fossil series in northwest, northeast, and southern Africa that chronologically and morphologically recapitulate the evolution of Homo sapiens from early archaic humans approximately 600 kya to early modern humans by 200 kya . This ancestral geographic structure is compatible with the higher diversity values observed in African populations compared to the in out-of-Africa populations . The recent dating of Middle Stone Age tools (315 ± 34 kya) and of early modern human fossils (286 ± 32 kya) from Jebel Irhoud in Morocco places the emergence of our species and of the Middle Stone Age, near to each other in time and long before the age of approximately 200 kya previously suggested for the common origin of all humans in eastern Africa . These data coincide in time with the existence of an old Y-chromosome lineage (A00) detected in samples of western-central African ascendance and dated 338 kya (95% CI: 237-581 kya). This age is much older than common estimates based on the Y-chromosome and mtDNA TMRCAs . The fact that the following more divergent Y-chromosome A lineages (A0, A1a) also have a western-central African location, strongly supports this region as the origin of an ancestral human population from which the ancestors of early modern humans emerged [90, 103]. The most ancient splits and expansions of mtDNA lineages are also situated at the hypothetical origin of all extant maternal lineages around this area. Although the earliest L0 clade diverged approximately 145 kya (Additional file 1: Table S3) and had its first expansions in southern Africa (L0d, L0k), the subsequent splits gave rise to L1 and L5 approximately 131 kya and 123 kya, respectively, spreading to western and eastern Africa respectively. These long- range African dispersals place the putative origin of L somewhere in Central Africa (Figure 1a). The same "center-of-gravity" argument has been used by other authors to suggest a Central African origin . While ancestral southern African Khoesan-speaking populations maintain high frequencies of primary L0d and k lineages [94, 106, 128, 129] and while the L1c haplogroup is dominant in the hunter-gatherer populations of central-western Africa [108, 109], L5 in eastern Africa has today only a marginal presence [114, 130]. This marginal presence is likely due to its displacement following incoming waves of better-adapted incomers. The presence of L5 in southern Africa and in eastern Mbuti pygmies [70, 109, 118, 128] is the result of later migrations. It is likely that the next split, approximately 100 kya, also occurred in Central Africa, resulting in sister clusters L2 and L3'4'6 that produced initial westward and eastward expansions, respectively (Fig. 1a). Although the oldest L2 lineages have been sampled in western Africa , today, as result of successive expansions within the continent, this clade has a pan-African range . In eastern Africa, the cluster L3'4'6 was the progenitor of the entire Eurasian maternal diversity. Its first lineage was haplogroup L6, which presently is a rare eastern lineage with a deep founder age (about 100 kya) but a rather recent expansion (about 25 kya). It has been found at frequencies below 1% in Egyptians , Somalis , Kenyans , and eastern Nilotes from Uganda . A higher frequency is observed in Ethiopia (3.15 ± 1.15 %) and the maximum frequency (15.8%) is observed in Ongota, an extinguishing linguistic isolate of uncertain adscription . Outside Africa, L6 has not been detected in the Levant . It is present in the Arabian Peninsula at frequencies below 1% in Saudi samples and as high as 12% in some Yemeni samples . Based on the L6 phylogeny (Additional file 2: Figure S1), it appears that not all of the Yemeni lineages are a subset of the eastern African lineages as there is at least one for which its common node coincides with the expansion of the whole haplogroup. Based on its peculiar phylogeography, the possibility that L6 originated from the same out-of-Africa southern migration that colonized Eurasia has been suggested . If this were the case, this early L6 expansion would provide genetic support to the reported presence of modern humans in the Arabian Peninsula, approximately 125 kya based on archaeological evidence [12,13,14]. This suggestion also receives climatic support as this period coincides with humid environmental conditions in Arabia . However, it appears that this potential human expansion did not extend beyond the peninsula, as L6-derived lineages have not yet been detected across Eurasia. It is likely that the return to arid conditions caused the decline of populations carrying the L6 lineage who had to retreat to refuge areas such as the highlands of Yemen and Ethiopia until more favorable conditions allowed their subsequent recovery in eastern Africa and Yemen. The long mutational stem that precedes the expansion of L6 (Additional file 2: Figure S1), is consistent with such a strong, long-term bottleneck. Subsequent lineage bifurcation produced the ancestors of the L3 and L4 haplogroups (Additional file 2: Figure S1). Presently, the highest frequencies and diversities of L4 are found in eastern Africa, but the haplogroup has spread over the entire continent (Table 3). Furthermore, it has been detected at frequencies below 1% in the Levant  and Arabian Peninsula [74, 139]. Some populations show outstanding frequencies of L4 likely due to drift effects. In western Africa, Samoya (28.6%) and Kassena (21.2%) speakers of the Gur linguistic family, show high L4 frequencies . In Ethiopia, samples of the Omotic-speaking Hamer (18.2%), the Cushitic-speaking Daasanach (22.2%), and the Nilotic-speaking Gumuz (24.0%) and Nyangatom (21.6%) also show high frequencies [130, 140].The Tanzanian click-speaking Hadza (58%) and Sandawe (43%) show the highest frequencies of L4 in Africa [111,112,113]. This latter finding, together with the elevated frequencies of the Y-chromosome haplogroup B-M112 among the Hadza (50%) and Sandawe (15%) , suggests that the human expansions from the North were those that most strongly influenced the gene pool of these groups. Considering the age of bifurcation from L3 (around 95 kya), these L4 expansions could have occurred before our proposed return of basal L3 lineages to Africa. However, as the main expansions of its descendant clusters L4a (54.8 kya) and L4b (48.9 kya) [94, 139] occurred around the same time window as those of the majority of the L3 and L2 branches in Africa, the most probable explanation is that improved climatic conditions after 60 kya motivated demographic growth on the African continent. Noticed that the evidence of a first expansion of L3 into East Africa  is likewise in support of the out-of-Africa scenario than of an Eurasian back-flow as proposed here. We hypothetically situated the L3'4 node in northeast Africa or the Near East (Fig. 1a) to allow an out-of-Africa expansion of the pre-L3 clade. The Y-chromosome CDEF ancestor must be this clade's male counterpart. Other female and male lineages might have moved with them; if so, they presumably, became extinct without contributing to either the maternal or paternal gene pools of the living human populations of Eurasia.
Under the scenario proposed here, early anatomically modern humans traveled out of Africa approximately 125 kya with a simple Middle Stone Age technology that was not superior to that manufactured by the Neanderthals. Favored by mild climatic conditions, these African pioneers progressed through West Asia and further followed the southern boundary of Central Asia, potentially overlapping en route the southern geographic range occupied by the now extinct hominins. Furthermore, as Neanderthals and early modern humans are associated with the same lithic industries in the Levant, it seems plausible that a new vision of the Early- and Mid-Middle Paleolithic fossil and archaeological records of those regions [87, 142, 143] might uncover the path followed by those early African colonizers. Under favorable conditions for both hominin groups, we might predict limited exchange of skills, lithic technology, and sex. However, when glacial environments became dominant after 75 kya, Neanderthals retreated southwards, possibly displacing humans as they went. Confronted with the northern foothills of the Himalayas, humans moved in two directions: westwards, eventually to return to Africa; and eastwards, eventually reaching southeastern Asia via China (Fig. 1b). The second part of this model has been outlined in preceding articles [53,54,55].
Cann RL, Stoneking M, Wilson AC. Mitochondrial DNA and human evolution. Nature. 325:31–6.
Pääbo S. The diverse origins of the human gene pool. Nat Rev Genet. 2015;16:313.
Ingman M, Kaessmann H, Pääbo S, Gyllensten U. Mitochondrial genome variation and the origin of modern humans. Nature. 2000;408:708–13.
Maca-Meyer N, González AM, Larruga JM, Flores C, Cabrera VM. Major genomic mitochondrial lineages delineate early human expansions. BMC Genet. 2001;2:1.
Herrnstadt C, Elson JL, Fahy E, Preston G, Turnbull DM, Anderson C, Ghosh SS, Olefsky JM, Beal MF, Davis RE, Howell N. 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:1152–71.
Soares P, Ermini L, Thomson N, Mormina M, Rito T, Röhl A, Salas A, Oppenheimer S, Macaulay V, Richards MB. Correcting for purifying selection: an improved human mitochondrial molecular clock. Am J Hum Genet. 2009;84:740–59.
Karmin M, Saag L, Vicente M, Sayres MAW, Järve M, Talas UG, Rootsi S, Ilumäe A-M, Mägi R, Mitt M, et al. A recent bottleneck of Y chromosome diversity coincides with a global change in culture. Genome Res. 2015;25:459–66.
Scally A, Durbin R. Revising the human mutation rate: implications for understanding human evolution. Nat Rev Genet. 2012;13:745–53.
Fu Q, Mittnik A, Johnson PL, Bos K, Lari M, Bollongino R, Sun C, Giemsch L, Schmitz R, Burger J, et al. A revised timescale for human evolution based on ancient mitochondrial genomes. Curr Biol 2013, 23:553–559.
Shea JJ, Bar-Yosef O. Who Were The Skhul/Qafzeh People? An Archaeological Perspective on Eurasia’s Oldest Modern Humans. J Isr Prehistoric Soc. 2005;35:451-468.
Liu W, Martinón-Torres M, Cai Y, Xing S, Tong H, Pei S, Sier MJ, Wu X, Edwards RL, Cheng H, et al. The earliest unequivocally modern humans in southern China. Nature. 2015;
Armitage SJ, Jasim SA, Marks AE, Parker AG, Usik VI, Uerpmann H-P. The southern route “out of Africa”: evidence for an early expansion of modern humans into Arabia. Science. 2011;331:453–6.
Rose JI, Usik VI, Marks AE, Hilbert YH, Galletti CS, Parton A, Geiling JM, Cerný V, Morley MW, Roberts RG. The Nubian complex of Dhofar, Oman: an African middle stone age industry in southern Arabia. PLoS One. 2011;6:e28239.
Petraglia MD, Alsharekh A, Breeze P, Clarkson C, Crassard R, Drake NA, Groucutt HS, Jennings R, Parker AG, Parton A, et al. Hominin dispersal into the Nefud desert and Middle Palaeolithic settlement along the Jubbah palaeolake, northern Arabia. PLoS One. 2012;7:e49840.
Sankararaman S, Patterson N, Li H, Pääbo S, Reich D. The date of interbreeding between Neandertals and modern humans. PLoS Genet. 2012;8:e1002947.
Kuhlwilm M, Gronau I, Hubisz MJ, de Filippo C, Prado-Martinez J, Kircher M, Fu Q, Burbano HA, Lalueza-Fox C, de La Rasilla M, et al. Ancient gene flow from early modern humans into Eastern Neanderthals. Nature. 2016;530:429–33.
Xing J, Watkins WS, Hu Y, Huff CD, Sabo A, Muzny DM, Bamshad MJ, Gibbs RA, Jorde LB, Yu F. Genetic diversity in India and the inference of Eurasian population expansion. Genome Biol. 2010;11:R113.
Pagani L, Lawson DJ, Jagoda E, Mӧrseburg A, Eriksson A, Mitt M, Clemente F, Hudjashov G, DeGiorgio M, Saag L, et al. Genomic analyses inform on migration events during the peopling of Eurasia. Nature. 2016;538:238–42.
Quintana-Murci L, Semino O, Bandelt H-J, Passarino G, McElreavey K, Santachiara-Benerecetti AS. Genetic evidence of an early exit of Homo sapiens sapiens from Africa through eastern Africa. Nat Genet. 1999;23:437–41.
Quintana-Murci L, Chaix R, Wells RS, Behar DM, Sayar H, Scozzari R, Rengo C, Al-Zahery N, Semino O, Santachiara-Benerecetti AS, et al. Where west meets east: the complex mtDNA landscape of the southwest and Central Asian corridor. Am J Hum Genet. 2004;74:827–45.
Qin Z, Yang Y, Kang L, Yan S, Cho K, Cai X, Lu Y, Zheng H, Zhu D, Fei D, et al. A mitochondrial revelation of early human migrations to the Tibetan Plateau before and after the last glacial maximum. Am J Phys Anthropol. 2010;143:555–69.
Kang L, Zheng H-X, Chen F, Yan S, Liu K, Qin Z, Liu L, Zhao Z, Li L, Wang X, et al. mtDNA lineage expansions in Sherpa population suggest adaptive evolution in Tibetan highlands. Mol Biol Evol. 2013:mst147.
Wang C-C, Wang L-X, Shrestha R, Zhang M, Huang X-Y, Hu K, Jin L, Li H: Genetic structure of Qiangic populations residing in the western Sichuan corridor. PLoS One 2014, 9:e103772.
Kong Q-P, Yao Y-G, Sun C, Bandelt H-J, Zhu C-L, Zhang Y-P. Phylogeny of east Asian mitochondrial DNA lineages inferred from complete sequences. Am J Hum Genet. 2003;73:671–6.
Tanaka M, Cabrera VM, González AM, Larruga JM, Takeyasu T, Fuku N, Guo L-J, Hirose R, Fujita Y, Kurata M, Shinoda K, Umetsu K, Yamada Y, Oshida Y, Sato Y, Hattori N, Mizuno Y, Arai Y, Hirose N, Ohta S, Ogawa O, Tanaka Y, Kawamori R, Shamoto-Nagai M, Maruyama W, Shimokata H, Suzuki R, Shimodaira H. Mitochondrial genome variation in eastern Asia and the peopling of Japan. Genome Res. 2004;14:1832–50.
Derenko M, Malyarchuk B, Denisova G, Perkova M, Rogalla U, Grzybowski T, Khusnutdinova E, Dambueva I, Zakharov I. Complete mitochondrial DNA analysis of eastern Eurasian haplogroups rarely found in populations of northern Asia and eastern Europe. PLoS One. 2012;7:e32179.
Metspalu M, Kivisild T, Metspalu E, Parik J, Hudjashov G, Kaldma K, Serk P, Karmin M, Behar DM, Gilbert MTP, et al. Most of the extant mtDNA boundaries in south and southwest Asia were likely shaped during the initial settlement of Eurasia by anatomically modern humans. BMC Genet. 2004;5:1.
Gounder Palanichamy M, Sun C, Agrawal S, Bandelt H-J, Kong Q-P, Khan F, Wang C-Y, Chaudhuri TK, Palla V, Zhang Y-P. Phylogeny of mitochondrial DNA macrohaplogroup N in India, based on complete sequencing: implications for the peopling of South Asia. Am J Hum Genet. 2004;75:966–78.
Cordaux R, Saha N, Bentley GR, Aunger R, Sirajuddin S, Stoneking M. Mitochondrial DNA analysis reveals diverse histories of tribal populations from India. Eur J Hum Genet. 2003;11:253–64.
Chandrasekar A, Kumar S, Sreenath J, Sarkar BN, Urade BP, Mallick S, Bandopadhyay SS, Barua P, Barik SS, Basu D. others: Updating phylogeny of mitochondrial DNA macrohaplogroup m in India: dispersal of modern human in South Asian corridor. PLoS One. 2009;4:e7447.
Hill C, Soares P, Mormina M, Macaulay V, Meehan W, Blackburn J, Clarke D, Raja JM, Ismail P, Bulbeck D, et al. Phylogeography and ethnogenesis of aboriginal Southeast Asians. Mol Biol Evol. 2006;23:2480–91.
Summerer M, Horst J, Erhart G, Weißensteiner H, Schӧnherr S, Pacher D, Forer L, Horst D, Manhart A, Horst B, et al. Large-scale mitochondrial DNA analysis in Southeast Asia reveals evolutionary effects of cultural isolation in the multi-ethnic population of Myanmar. BMC Evol Biol. 2014;14:1.
Li Y-C, Wang H-W, Tian J-Y, Liu L-N, Yang L-Q, Zhu C-L, Wu S-F, Kong Q-P, Zhang Y-P. Ancient inland human dispersals from Myanmar into interior East Asia since the Late Pleistocene. Sci Rep. 2015;5
Hill C, Soares P, Mormina M, Macaulay V, Clarke D, Blumbach PB, Vizuete-Forster M, Forster P, Bulbeck D, Oppenheimer S, et al. A mitochondrial stratigraphy for island southeast Asia. Am J Hum Genet. 2007;80:29–43.
Mona S, Grunz KE, Brauer S, Pakendorf B, Castrì L, Sudoyo H, Marzuki S, Barnes RH, Schmidtke J, Stoneking M, et al. Genetic admixture history of Eastern Indonesia as revealed by Y-chromosome and mitochondrial DNA analysis. Mol Biol Evol. 2009;26:1865–77.
Tabbada KA, Trejaut J, Loo J-H, Chen Y-M, Lin M, Mirazón-Lahr M, Kivisild T, De Ungria MCA. Philippine mitochondrial DNA diversity: a populated viaduct between Taiwan and Indonesia? Mol Biol Evol. 2010;27:21–31.
Gunnarsdóttir ED, Li M, Bauchet M, Finstermeier K, Stoneking M. High-throughput sequencing of complete human mtDNA genomes from the Philippines. Genome Res. 2011;21:1–11.
Tumonggor MK, Karafet TM, Hallmark B, Lansing JS, Sudoyo H, Hammer MF, Cox MP. The Indonesian archipelago: an ancient genetic highway linking Asia and the Pacific. J Hum Genet. 2013;58:165–73.
Gomes SM, Bodner M, Souto L, Zimmermann B, Huber G, Strobl C, Rӧck AW, Achilli A, Olivieri A, Torroni A, et al. Human settlement history between Sunda and Sahul: a focus on East Timor (Timor-Leste) and the Pleistocenic mtDNA diversity. BMC Genomics. 2015;16:1.
Ingman M, Gyllensten U. Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines. Genome Res. 2003;13:1600–6.
Van Holst Pellekaan SM, Ingman M, Roberts-Thomson J, Harding RM. Mitochondrial genomics identifies major haplogroups in Aboriginal Australians. Am J Phys Anthropol. 2006;131:282–94.
Hudjashov G, Kivisild T, Underhill PA, Endicott P, Sanchez JJ, Lin AA, Shen P, Oefner P, Renfrew C, Villems R, et al. Revealing the prehistoric settlement of Australia by Y chromosome and mtDNA analysis, vol. 104: Proceedings of the National Academy of Sciences; 2007. p. 8726–30.
Friedlaender JS, Friedlaender FR, Hodgson JA, Stoltz M, Koki G, Horvat G, Zhadanov S, Schurr TG, Merriwether DA. Melanesian mtDNA complexity. PLoS One. 2007;2:e248.
Duggan AT, Evans B, Friedlaender FR, Friedlaender JS, Koki G, Merriwether DA, Kayser M, Stoneking M. Maternal history of Oceania from complete mtDNA genomes: contrasting ancient diversity with recent homogenization due to the Austronesian expansion. Am J Hum Genet. 2014;94:721–33.
Macaulay V, Hill C, Achilli A, Rengo C, Clarke D, Meehan W, Blackburn J, Semino O, Scozzari R, Cruciani F, et al. Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science. 2005;308:1034–6.
Kong Q-P, Sun C, Wang H-W, Zhao M, Wang W-Z, Zhong L, Hao X-D, Pan H, Wang S-Y, Cheng Y-T, 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:513–22.
Fernandes V, Alshamali F, Alves M, Costa MD, Pereira JB, Silva NM, Cherni L, Harich N, Cerny V, Soares P, et al. The Arabian cradle: mitochondrial relicts of the first steps along the southern route out of Africa. Am J Hum Genet. 2012;90:347–55.
Cordaux R, Weiss G, Saha N, Stoneking M. The northeast Indian passageway: a barrier or corridor for human migrations? Mol Biol Evol. 2004;21:1525–33.
Merriwether DA, Hodgson JA, Friedlaender FR, Allaby R, Cerchio S, Koki G, Friedlaender JS. Ancient mitochondrial M haplogroups identified in the Southwest Pacific. Proc Natl Acad Sci U S A. 2005;102:13034–9.
Van Holst Pellekaan S. Genetic evidence for the colonization of Australia. Quat Int. 2013;285:44–56.
Pagani L, Schiffels S, Gurdasani D, Danecek P, Scally A, Chen Y, Xue Y, Haber M, Ekong R, Oljira T, et al. Tracing the route of modern humans out of Africa by using 225 human genome sequences from Ethiopians and Egyptians. Am J Hum Genet. 2015;96:986–91.
Sankararaman S, Mallick S, Dannemann M, Prüfer K, Kelso J, Pääbo S, Patterson N, Reich D. The genomic landscape of Neanderthal ancestry in present-day humans. Nature. 2014;507:354–7.
Fregel R, Cabrera V, Larruga JM, Abu-Amero KK, González AM. Carriers of Mitochondrial DNA Macrohaplogroup N Lineages Reached Australia around 50,000 Years Ago following a Northern Asian Route. PLoS One. 2015;10:e0129839.
Marrero P, Abu-Amero KK, Larruga JM, Cabrera VM. Carriers of human mitochondrial DNA macrohaplogroup M colonized India from southeastern Asia. BMC Evol Biol. 2016;16:246.
Larruga JM, Marrero P, Abu-Amero KK, Golubenko MV, Cabrera VM. Carriers of mitochondrial DNA macrohaplogroup R colonized Eurasia and Australasia from a southeast Asia core area. BMC Evol Biol. 2017;17:115.
Magoon GR, Banks RH, Rottensteiner C, Schrack BE, Tilroe VO, Grierson AJ. Generation of high-resolution a priori Y-chromosome phylogenies using “next-generation” sequencing data. In: bioRxiv; 2013. p. 000802.
Karafet TM, Mendez FL, Sudoyo H, Lansing JS, Hammer MF. Improved phylogenetic resolution and rapid diversification of Y-chromosome haplogroup K-M526 in Southeast Asia. Eur J Hum Genet. 2014;
Hallast P, Batini C, Zadik D, Delser PM, Wetton JH, Arroyo-Pardo E, Cavalleri GL, De Knijff P, Bisol GD, Dupuy BM, et al. The Y-chromosome tree bursts into leaf: 13,000 high-confidence SNPs covering the majority of known clades. Mol Biol Evol. 2015;32:661–73.
Poznik GD, Xue Y, Mendez FL, Willems TF, Massaia A, Sayres MAW, Ayub Q, McCarthy SA, Narechania A, Kashin S, et al. Punctuated bursts in human male demography inferred from 1,244 worldwide Y-chromosome sequences. Nat Genet. 2016;48:593–9.
Groucutt HS, Petraglia MD, Bailey G, Scerri EM, Parton A, Clark-Balzan L, Jennings RP, Lewis L, Blinkhorn J, Drake NA, et al. Rethinking the dispersal of Homo sapiens out of Africa. Evol Anthropol Issues News Reviews. 2015;24:149–64.
Hammer MF, Karafet T, Rasanayagam A, Wood ET, Altheide TK, Jenkins T, Griffiths RC, Templeton AR, Zegura SL. Out of Africa and back again: nested cladistic analysis of human Y chromosome variation. Mol Biol Evol. 1998;15:427–41.
Fregel R, Delgado S. HaploSearch: a tool for haplotype-sequence two-way transformation. Mitochondrion. 2011;11:366–7.
Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, Howell N. Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet. 1999;23:147.
Hall TA. BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser. 1999;41:95–8.
Van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009;30:E386–94.
Bandelt H-J, Forster P, Rӧhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.
Forster P, Harding R, Torroni A, Bandelt HJ. Origin and evolution of Native American mtDNA variation: a reappraisal. Am J Hum Genet. 1996;59:935–45.
Saillard J, Forster P, Lynnerup N, Bandelt HJ, Nørby S. mtDNA variation among Greenland Eskimos: the edge of the Beringian expansion. Am J Hum Genet. 2000;67:718–26.
Collins DW, Gudiseva HV, Trachtman BT, Jerrehian M, Gorry T, Merritt WT, Rhodes AL, Sankar PS, Regina M, Miller-Ellis E, O’Brien JM. Mitochondrial sequence variation in African-American primary open-angle glaucoma patients. PLoS One. 2013;8:e76627.
Barbieri C, Vicente M, Oliveira S, Bostoen K, Rocha J, Stoneking M, Pakendorf B. Migration and interaction in a contact zone: mtDNA variation among Bantu-speakers in southern Africa. PLoS One. 2014;9:e99117.
Pardiñas AF, Martínez JL, Roca A, García-Vazquez E, López B. Over the sands and far away: interpreting an Iberian mitochondrial lineage with ancient Western African origins. Am J Hum Biol. 2014;26:777–83.
Barbieri C, Whitten M, Beyer K, Schreiber H, Li M, Pakendorf B. Contrasting maternal and paternal histories in the linguistic context of Burkina Faso. Mol Biol Evol. 2011;29:1213–23.
Vyas DN, Kitchen A, Miró-Herrans AT, Pearson LN, Al-Meeri A, Mulligan CJ. Bayesian analyses of Yemeni mitochondrial genomes suggest multiple migration events with Africa and Western Eurasia. Am J Phys Anthropol. 2015;
Abu-Amero KK, Larruga JM, Cabrera VM, González AM. Mitochondrial DNA structure in the Arabian Peninsula. BMC Evol Biol. 2008;8:1.
Mellars P, Gori KC, Carr M, Soares PA, Richards MB. Genetic and archaeological perspectives on the initial modern human colonization of southern Asia. Proc Natl Acad Sci U S A. 2013;110:10699–704.
Underhill PA, Passarino G, Lin AA, Shen P, Lahr MM, Foley RA, Oefner PJ, Cavalli-Sforza LL. The phylogeography of Y chromosome binary haplotypes and the origins of modern human populations. Ann Hum Genet. 2001;65:43–62.
Weale ME, Shah T, Jones AL, Greenhalgh J, Wilson JF, Nymadawa P, Zeitlin D, Connell BA, Bradman N, Thomas MG. Rare deep-rooting Y chromosome lineages in humans: lessons for phylogeography. Genetics. 2003;165:229–34.
Shi H, Zhong H, Peng Y, Dong Y-L, Qi X-B, Zhang F, Liu L-F, Tan S-J, Ma RZ, Xiao C-J, et al. Y chromosome evidence of earliest modern human settlement in East Asia and multiple origins of Tibetan and Japanese populations. BMC Biol. 2008;6:1.
Karafet TM, Mendez FL, Meilerman MB, Underhill PA, Zegura SL, Hammer MF. New binary polymorphisms reshape and increase resolution of the human Y chromosomal haplogroup tree. Genome Res. 2008;18:830–8.
Underhill PA, Kivisild T. Use of y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu Rev Genet. 2007;41:539–64.
Su B, Xiao C, Deka R, Seielstad MT, Kangwanpong D, Xiao J, Lu D, Underhill P, Cavalli-Sforza L, Chakraborty R, et al. Y chromosome haplotypes reveal prehistorical migrations to the Himalayas. Hum Genet. 2000;107:582–90.
Thangaraj K, Singh L, Reddy AG, Rao VR, Sehgal SC, Underhill PA, Pierson M, Frame IG, Hagelberg E. Genetic affinities of the Andaman Islanders, a vanishing human population. Curr Biol. 2003;13:86–93.
Shea JJ. Transitions or turnovers? Climatically-forced extinctions of Homo sapiens and Neanderthals in the east Mediterranean Levant. Quat Sci Rev. 2008;27:2253–70.
Hershkovitz I, Marder O, Ayalon A, Bar-Matthews M, Yasur G, Boaretto E, Caracuta V, Alex B, Frumkin A, Goder-Goldberger M, et al. Levantine cranium from Manot Cave (Israel) foreshadows the first European modern humans. Nature. 2015;520:216–219.
Rougier H, Milota S, Rodrigo R, Gherase M, Sarcina L, Moldovan O, Zilhão J, Constantin S, Franciscus RG, Zollikofer CPE, Ponce de León M, Trinkaus E. Peştera cu Oase 2 and the cranial morphology of early modern Europeans. Proc Natl Acad Sci U S A. 2007;104:1165–70.
Benazzi S, Slon V, Talamo S, Negrino F, Peresani M, Bailey SE, Sawyer S, Panetta D, Vicino G, Starnini E, et al. The makers of the Protoaurignacian and implications for Neandertal extinction. Science. 2015;348:793–6.
Marks AE. Comments after four decades of research on the Middle to Upper Paleolithic transition. Mitteilungen der Gesellschaft für Urgeschichte. 2005;14:81–6.
Otte M. Arguments for population movement of anatomically modern humans from Central Asia to Europe. Rethinking Human Revol. 2007:359–66.
Soares P, Alshamali F, Pereira JB, Fernandes V, Silva NM, Afonso C, Costa MD, Musilová E, Macaulay V, Richards MB, Cerny V, Pereira L. The Expansion of mtDNA Haplogroup L3 within and out of Africa. Mol Biol Evol. 2012;29:915–27.
Scozzari R, Massaia A, Trombetta B, Bellusci G, Myres NM, Novelletto A, Cruciani F. An unbiased resource of novel SNP markers provides a new chronology for the human Y chromosome and reveals a deep phylogenetic structure in Africa. Genome Res. 2014;24:535–44.
Harich N, Costa MD, Fernandes V, Kandil M, Pereira JB, Silva NM, Pereira L. The trans-Saharan slave trade-clues from interpolation analyses and high-resolution characterization of mitochondrial DNA lineages. BMC Evol Biol. 2010;10:138.
Kefi R, Hechmi M, Naouali C, Jmel H, Hsouna S, Bouzaid E, Abdelhak S, Beraud-Colomb E, Stevanovitch A. On the origin of Iberomaurusians: new data based on ancient mitochondrial DNA and phylogenetic analysis of Afalou and Taforalt populations. Mitochondrial DNA Part A. 2016:1–11.
Schuenemann VJ, Peltzer A, Welte B, van Pelt WP, Molak M, Wang C-C, Furtwängler A, Urban C, Reiter E, Nieselt K, et al. Ancient Egyptian mummy genomes suggest an increase of Sub-Saharan African ancestry in post-Roman periods. Nat Commun. 2017;8
Behar DM, Villems R, Soodyall H, Blue-Smith J, Pereira L, Metspalu E, Scozzari R, Makkan H, Tzur S, Comas D, et al. The dawn of human matrilineal diversity. Am J Hum Genet. 2008;82:1130–40.
Cruciani F, La Fratta R, Trombetta B, Santolamazza P, Sellitto D, Colomb EB, Dugoujon J-M, Crivellaro F, Benincasa T, Pascone R, et al. Tracing past human male movements in northern/eastern Africa and western Eurasia: new clues from Y-chromosomal haplogroups E-M78 and J-M12. Mol Biol Evol. 2007;24:1300–11.
Bosch E, Calafell F, Comas D, Oefner PJ, Underhill PA, Bertranpetit J. High-resolution analysis of human Y-chromosome variation shows a sharp discontinuity and limited gene flow between northwestern Africa and the Iberian Peninsula. Am J Hum Genet. 2001;68:1019–29.
Flores C, Maca-Meyer N, Perez JA, Hernandez M, Cabrera VM. Y-chromosome differentiation in Northwest Africa. Hum Biol. 2001;73:513–24.
Arredi B, Poloni ES, Paracchini S, Zerjal T, Fathallah DM, Makrelouf M, Pascali VL, Novelletto A, Tyler-Smith C. A predominantly neolithic origin for Y-chromosomal DNA variation in North Africa. Am J Hum Genet. 2004;75:338–45.
Trombetta B, Cruciani F, Sellitto D, Scozzari R. A new topology of the human Y chromosome haplogroup E1b1 (E-P2) revealed through the use of newly characterized binary polymorphisms. PLoS One. 2011;6:e16073.
Trombetta B, D’Atanasio E, Massaia A, Ippoliti M, Coppa A, Candilio F, Coia V, Russo G, Dugoujon J-M, Moral P, et al. Phylogeographic refinement and large scale genotyping of human Y chromosome haplogroup E provide new insights into the dispersal of early pastoralists in the African continent. Genome Biol Evol. 2015;7:1940–50.
Secher B, Fregel R, Larruga JM, Cabrera VM, Endicott P, Pestano JJ, González AM. The history of the North African mitochondrial DNA haplogroup U6 gene flow into the African, Eurasian and American continents. BMC Evol Biol. 2014;14:109.
Semino O, Magri C, Benuzzi G, Lin AA, Al-Zahery N, Battaglia V, Maccioni L, Triantaphyllidis C, Shen P, Oefner PJ, et al. Origin, diffusion, and differentiation of Y-chromosome haplogroups E and J: inferences on the neolithization of Europe and later migratory events in the Mediterranean area. Am J Hum Genet. 2004;74:1023–34.
Cruciani F, Trombetta B, Massaia A, Destro-Bisol G, Sellitto D, Scozzari R. A revised root for the human Y chromosomal phylogenetic tree: the origin of patrilineal diversity in Africa. Am J Hum Genet. 2011;88:814–8.
Abu-Amero KK, Hellani A, González AM, Larruga JM, Cabrera VM, Underhill PA. Saudi Arabian Y-Chromosome diversity and its relationship with nearby regions. BMC Genet. 2009;10:1.
Zalloua PA, Xue Y, Khalife J, Makhoul N, Debiane L, Platt DE, Royyuru AK, Herrera RJ, Hernanz DFS, Blue-Smith J, et al. Y-chromosomal diversity in Lebanon is structured by recent historical events. Am J Hum Genet. 2008;82:873–82.
Chen Y-S, Olckers A, Schurr TG, Kogelnik AM, Huoponen K, Wallace DC. mtDNA variation in the South African Kung and Khwe—and their genetic relationships to other African populations. Am J Hum Genet. 2000;66:1362–83.
Soodyall H, Makkan H, Haycock P, Naidoo T. The genetic prehistory of the Khoe and San. South Afr Humanit. 2008;20:37–48.
Quintana-Murci L, Quach H, Harmant C, Luca F, Massonnet B, Patin E, Sica L, Mouguiama-Daouda P, Comas D, Tzur S. Maternal traces of deep common ancestry and asymmetric gene flow between Pygmy hunter-gatherers and Bantu-speaking farmers. Proc Natl Acad Sci. 2008;105:1596–601.
Batini C, Lopes J, Behar DM, Calafell F, Jorde LB, Van der Veen L, Quintana-Murci L, Spedini G, Destro-Bisol G, Comas D. Insights into the demographic history of African Pygmies from complete mitochondrial genomes. Mol Biol Evol. 2010;28:1099–110.
Anagnostou P, Battaggia C, Capocasa M, Boschi I, Brisighelli F, Batini C, Spedini G, Destro-Bisol G. Reevaluating a Model of Gender-Biased Gene Flow among Sub-Saharan Hunter-Gatherers and Farmers. Hum Biol. 2013;85:597–606.
Knight A, Underhill PA, Mortensen HM, Zhivotovsky LA, Lin AA, Henn BM, Louis D, Ruhlen M, Mountain JL. African Y chromosome and mtDNA divergence provides insight into the history of click languages. Curr Biol. 2003;13:464–73.
Henn BM, Gignoux C, Lin AA, Oefner PJ, Shen P, Scozzari R, Cruciani F, Tishkoff SA, Mountain JL, Underhill PA. Y-chromosomal evidence of a pastoralist migration through Tanzania to southern Africa. Proceedings of the, vol. 105: National Academy of Sciences; 2008. p. 10693–8.
Tishkoff SA, Gonder MK, Henn BM, Mortensen H, Knight A, Gignoux C, Fernandopulle N, Lema G, Nyambo TB, Ramakrishnan U, et al. History of click-speaking populations of Africa inferred from mtDNA and Y chromosome genetic variation. Mol Biol Evol. 2007;24:2180–95.
Gomes V, Pala M, Salas A, Álvarez-Iglesias V, Amorim A, Gómez-Carballa A, Carracedo Á, Clarke DJ, Hill C, Mormina M, et al. Mosaic maternal ancestry in the Great Lakes region of East Africa. Hum Genet 2015, 134:1013.
Salas A, Richards M, De la Fe T, Lareu M-V, Sobrino B, Sánchez-Diz P, Macaulay V, Carracedo Á. The making of the African mtDNA landscape. Am J Hum Genet. 2002;71:1082–111.
Coelho M, Sequeira F, Luiselli D, Beleza S, Rocha J. On the edge of Bantu expansions: mtDNA, Y chromosome and lactase persistence genetic variation in southwestern Angola. BMC Evol Biol. 2009;9:80.
Pour NA, Plaster CA, Bradman N. Evidence from Y-chromosome analysis for a late exclusively eastern expansion of the Bantu-speaking people. Eur J Hum Genet. 2013;21:423–9.
Marks SJ, Montinaro F, Levy H, Brisighelli F, Ferri G, Bertoncini S, Batini C, Busby GB, Arthur C, Mitchell P, et al. Static and moving frontiers: the genetic landscape of Southern African Bantu-speaking populations. Mol Biol Evol. 2014;32:29–43.
Silva M, Alshamali F, Silva P, Carrilho C, Mandlate F, Trovoada MJ, Černý V, Pereira L, Soares P. 60,000 years of interactions between Central and Eastern Africa documented by major African mitochondrial haplogroup L2. Sci Rep. 2015;5
Watson E, Forster P, Richards M, Bandelt HJ. Mitochondrial footprints of human expansions in Africa. Am J Hum Genet. 1997;61:691–704.
Prüfer K, Racimo F, Patterson N, Jay F, Sankararaman S, Sawyer S, Heinze A, Renaud G, Sudmant PH, De Filippo C, et al. The complete genome sequence of a Neanderthal from the Altai Mountains. Nature. 2014;505:43–9.
Sánchez-Quinto F, Botigué LR, Civit S, Arenas C, Ávila-Arcos MC, Bustamante CD, Comas D, Lalueza-Fox C. North African populations carry the signature of admixture with Neandertals. PLoS One. 2012;7:e47765.
Bräuer G. The origin of modern anatomy: by speciation or intraspecific evolution? Evol Anthropol Issues News Reviews. 2008;17:22–37.
Garrigan D, Hammer MF. Reconstructing human origins in the genomic era. Nat Rev Genet. 2006;7:669–80.
MISSING:richter2017age: MISSING:richter2017age. 2018.
Mendez FL, Krahn T, Schrack B, Krahn A-M, Veeramah KR, Woerner AE, Fomine FLM, Bradman N, Thomas MG, Karafet TM, et al. An African American paternal lineage adds an extremely ancient root to the human Y chromosome phylogenetic tree. Am J Hum Genet. 2013;92:454–9.
Rito T, Richards MB, Fernandes V, Alshamali F, Cerny V, Pereira L, Soares P. The first modern human dispersals across Africa. PLoS One. 2013;8:e80031.
Barbieri C, Güldemann T, Naumann C, Gerlach L, Berthold F, Nakagawa H, Mpoloka SW, Stoneking M, Pakendorf B. Unraveling the complex maternal history of Southern African Khoisan populations. Am J Phys Anthropol. 2014;153:435–48.
Chan EK, Hardie R-A, Petersen DC, Beeson K, Bornman RM, Smith AB, Hayes VM. Revised timeline and distribution of the earliest diverged human maternal lineages in southern Africa. PLoS One. 2015;10:e0121223.
Boattini A, Castrì L, Sarno S, Useli A, Cioffi M, Sazzini M, Garagnani P, De Fanti S, Pettener D, Luiselli D. mtDNA variation in East Africa unravels the history of afro-asiatic groups. Am J Phys Anthropol. 2013;150:375–85.
Rando J, Pinto F, Gonzalez A, Hernandez M, Larruga J, Cabrera V, Bandelt H-J. Mitochondrial DNA analysis of Northwest African populations reveals genetic exchanges with European, Near-Eastern, and sub-Saharan populations. Ann Hum Genet. 1998;62:531–50.
Elmadawy MA, Nagai A, Gomaa GM, Hegazy HM, Shaaban FE, Bunai Y. Investigation of mtDNA control region sequences in an Egyptian population sample. Legal Med. 2013;15:338–41.
Mikkelsen M, Fendt L, Rӧck AW, Zimmermann B, Rockenbauer E, Hansen AJ, Parson W, Morling N. Forensic and phylogeographic characterisation of mtDNA lineages from Somalia. Int J Legal Med. 2012;126:573.
Batai K, Babrowski KB, Arroyo JP, Kusimba CM, Williams SR. Mitochondrial DNA diversity in two ethnic groups in Southeastern Kenya: perspectives from the northeastern periphery of the Bantu expansion. Am J Phys Anthropol. 2013;150:482–91.
Badro DA, Douaihy B, Haber M, Youhanna SC, Salloum A, Ghassibe-Sabbagh M, Johnsrud B, Khazen G, Matisoo-Smith E, Soria-Hernanz DF, et al. Y-chromosome and mtDNA genetics reveal significant contrasts in affinities of modern Middle Eastern populations with European and African populations. PLoS One. 2013;8:e54616.
Kivisild T, Reidla M, Metspalu E, Rosa A, Brehm A, Pennarun E, Parik J, Geberhiwot T, Usanga E, Villems R. Ethiopian mitochondrial DNA heritage: tracking gene flow across and around the gate of tears. Am J Hum Genet. 2004;75:752–70.
Rosenberg T, Preusser F, Fleitmann D, Schwalb A, Penkman K, Schmid T, Al-Shanti M, Kadi K, Matter A. Humid periods in southern Arabia: windows of opportunity for modern human dispersal. Geology. 2011;39:1115–8.
Bekada A, Fregel R, Cabrera VM, Larruga JM, Pestano J, Benhamamouch S, González AM. Introducing the Algerian mitochondrial DNA and Y-chromosome profiles into the North African landscape. PLoS One. 2013;8:e56775.
Fernandes V, Triska P, Pereira JB, Alshamali F, Rito T, Machado A, Fajkošová Z, Cavadas B, Černý V, Soares P, et al. Genetic stratigraphy of key demographic events in Arabia. PLoS One. 2015;10:e0118625.
Poloni ES, Naciri Y, Bucho R, Niba R, Kervaire B, Excoffier L, Langaney A, Sanchez-Mazas A. Genetic evidence for complexity in ethnic differentiation and history in East Africa. Ann Hum Genet. 2009;73:582–600.
Wood ET, Stover DA, Ehret C, Destro-Bisol G, Spedini G, McLeod H, Louie L, Bamshad M, Strassmann BI, Soodyall H, et al. Contrasting patterns of Y chromosome and mtDNA variation in Africa: evidence for sex-biased demographic processes. Eur J Hum Genet. 2005;13:867–76.
Vishnyatsky LB. The Paleolithic of Central Asia. J World Prehist. 1999;13:69–122.
Glantz MM. The history of hominin occupation of Central Asia in review. In: Asian Paleoanthropology: Springer; 2011. p. 101–12.
We are grateful to Dra. Ana M. González for her experimental contribution and bright ideas brought to this work. This research was supported by Grant n° CGL2010-16195 from the Spanish Ministerio de Ciencia e Innovación to JML.
Availability of data and materials
The sequence set supporting the results of this article are available in the GenBank repository (MF621062 to MF621130), Additional file 1: Table S1, and Additional file 2: Figure S1. References for the haplogroup frequencies used in this study are listed in Additional file 1: Table S2. All results obtained from our statistical analysis are presented in tables and figures of this article and in the additional files.
VMC is actually retired.
Ethics approval and consent to participate
The procedure of human population sampling adhered to the tenets of the Declaration of Helsinki. Written consent was recorded from all participants prior to taking part in the study. The study underwent formal review and was approved by the College of Medicine Ethical Committee of the King Saud University (proposal N°09-659), and by the Ethics Committee for Human Research at the University of La Laguna (proposal NR157).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Table S1. Complete mtDNA macrohaplogroup L sequences. Table S2. Frequencies of mtDNA haplogroups L2 and L3 and Y-chromosome haplogroup E lineages across Africa. Table S3. Coalescence ages in thousand years (kya) with 95% coefficient intervals (CI), or standard deviations, for the main mitochondrial DNA African haplogroups. Table S4. Coalescence ages in thousand years (kya) with 95% coefficient intervals (CI), or standard deviations, for Y-chromosome most recent common ancestor (MRCA), the out-of-Africa event, and the splits of haplogroup DE and E. Table S5. Population clustering into five classes. Table S6. k-means cluster results using African populations characterized by mtDNA L3 and Y-chromosome E haplogroup frequencies. Table S7. k-means cluster results using African populations characterized by mtDNA L2 and L3 and Y-chromosome E haplogroup frequencies. (XLSX 403 kb)
Figure S1. Phylogenetic tree of mtDNA macrohaplogroup L complete African sequences produced in this study. (XLSX 71 kb)