- Research article
- Open Access
Carriers of mitochondrial DNA macrohaplogroup L3 basal lineages migrated back to Africa from Asia around 70,000 years ago
© The Author(s). 2018
- Received: 10 August 2017
- Accepted: 5 June 2018
- Published: 19 June 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.
- Human evolution
- Mitochondrial DNA
- Haplogroup L3
- Haplogroup E
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–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–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–23] and East Asia [24–26] but also from the regions along the hypothetical southern route, such as India [27–30], mainland southeastern Asia [31–33], island southeast Asia [34–39], New Guinea, North Island, Melanesia and Australia [40–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–47]. Others postulate ancient local population differentiation in each region without any evidence of the shared ancestry expected by the southern dispersal model [48–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–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, v22.214.171.124, 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.
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
MtDNA and NRY mean values for MRCA, Out-of-Africa and haplogroup coalescences
184 ± 61.0
71.0 ± 12.0a
L3'4'6: 95.8 ± 14.0
71.0 ± 8.0b
L3'4: 84.1 ± 8.6
171.5 ± 13.7
93.9 ± 25.3
DE: 69.0 ± 6.7
E: 65.5 ± 8.5
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 phylogeography of the L3 and E lineages inside Africa
MtDNA and Y-Chromosome mean haplogroup frequencies in the major regions of the African continent
AMOVA and k-mean clustering results
Variance % within populations
Frequency values for k-means centroids in 1 to 5 classes
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–110], eastern African forested areas around the Great Lakes [111–114], and the semi-desert open spaces of South Africa [115–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–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–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–55].
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 conceived and designed the study, analyzed the data and wrote the manuscript. PM edited and submitted mtDNA sequences, designed figures, and contributed to the data analysis independently confirming the analysis results. KKAA carried out the sequencing of the Arabian and eastern African samples and made corrections on the manuscript. JML carried out the sequencing of the Canary Islands and western African samples, and contributed to the collection of published data and their analysis. All the authors read and approved the final manuscript.
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.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Cann RL, Stoneking M, Wilson AC. Mitochondrial DNA and human evolution. Nature. 325:31–6.Google Scholar
- Pääbo S. The diverse origins of the human gene pool. Nat Rev Genet. 2015;16:313.PubMedView ArticleGoogle Scholar
- Ingman M, Kaessmann H, Pääbo S, Gyllensten U. Mitochondrial genome variation and the origin of modern humans. Nature. 2000;408:708–13.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Scally A, Durbin R. Revising the human mutation rate: implications for understanding human evolution. Nat Rev Genet. 2012;13:745–53.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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;Google Scholar
- 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.PubMedView ArticlePubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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;5Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Ingman M, Gyllensten U. Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines. Genome Res. 2003;13:1600–6.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Holst Pellekaan S. Genetic evidence for the colonization of Australia. Quat Int. 2013;285:44–56.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.Google Scholar
- 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;Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Fregel R, Delgado S. HaploSearch: a tool for haplotype-sequence two-way transformation. Mitochondrion. 2011;11:366–7.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticlePubMed CentralGoogle Scholar
- 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.Google Scholar
- Van Oven M, Kayser M. Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat. 2009;30:E386–94.PubMedView ArticleGoogle Scholar
- Bandelt H-J, Forster P, Rӧhl A. Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999;16:37–48.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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;Google Scholar
- Abu-Amero KK, Larruga JM, Cabrera VM, González AM. Mitochondrial DNA structure in the Arabian Peninsula. BMC Evol Biol. 2008;8:1.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Underhill PA, Kivisild T. Use of y chromosome and mitochondrial DNA population structure in tracing human migrations. Annu Rev Genet. 2007;41:539–64.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- Otte M. Arguments for population movement of anatomically modern humans from Central Asia to Europe. Rethinking Human Revol. 2007:359–66.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.Google Scholar
- 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;8Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Flores C, Maca-Meyer N, Perez JA, Hernandez M, Cabrera VM. Y-chromosome differentiation in Northwest Africa. Hum Biol. 2001;73:513–24.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Soodyall H, Makkan H, Haycock P, Naidoo T. The genetic prehistory of the Khoe and San. South Afr Humanit. 2008;20:37–48.Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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;5Google Scholar
- Watson E, Forster P, Richards M, Bandelt HJ. Mitochondrial footprints of human expansions in Africa. Am J Hum Genet. 1997;61:691–704.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- Bräuer G. The origin of modern anatomy: by speciation or intraspecific evolution? Evol Anthropol Issues News Reviews. 2008;17:22–37.View ArticleGoogle Scholar
- Garrigan D, Hammer MF. Reconstructing human origins in the genomic era. Nat Rev Genet. 2006;7:669–80.PubMedView ArticleGoogle Scholar
- MISSING:richter2017age: MISSING:richter2017age. 2018.Google Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedPubMed CentralView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- 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.PubMedView ArticleGoogle Scholar
- Vishnyatsky LB. The Paleolithic of Central Asia. J World Prehist. 1999;13:69–122.View ArticleGoogle Scholar
- Glantz MM. The history of hominin occupation of Central Asia in review. In: Asian Paleoanthropology: Springer; 2011. p. 101–12.Google Scholar