Mitochondrial DNA structure in the Arabian Peninsula
© Abu-Amero et al; licensee BioMed Central Ltd. 2008
Received: 17 September 2007
Accepted: 12 February 2008
Published: 12 February 2008
Two potential migratory routes followed by modern humans to colonize Eurasia from Africa have been proposed. These are the two natural passageways that connect both continents: the northern route through the Sinai Peninsula and the southern route across the Bab al Mandab strait. Recent archaeological and genetic evidence have favored a unique southern coastal route. Under this scenario, the study of the population genetic structure of the Arabian Peninsula, the first step out of Africa, to search for primary genetic links between Africa and Eurasia, is crucial. The haploid and maternally inherited mitochondrial DNA (mtDNA) molecule has been the most used genetic marker to identify and to relate lineages with clear geographic origins, as the African Ls and the Eurasian M and N that have a common root with the Africans L3.
To assess the role of the Arabian Peninsula in the southern route, we genetically analyzed 553 Saudi Arabs using partial (546) and complete mtDNA (7) sequencing, and compared the lineages obtained with those present in Africa, the Near East, central, east and southeast Asia and Australasia. The results showed that the Arabian Peninsula has received substantial gene flow from Africa (20%), detected by the presence of L, M1 and U6 lineages; that an 18% of the Arabian Peninsula lineages have a clear eastern provenance, mainly represented by U lineages; but also by Indian M lineages and rare M links with Central Asia, Indonesia and even Australia. However, the bulk (62%) of the Arabian lineages has a Northern source.
Although there is evidence of Neolithic and more recent expansions in the Arabian Peninsula, mainly detected by (preHV)1 and J1b lineages, the lack of primitive autochthonous M and N sequences, suggests that this area has been more a receptor of human migrations, including historic ones, from Africa, India, Indonesia and even Australia, than a demographic expansion center along the proposed southern coastal route.
The hypothesis that modern humans originated in Africa and later migrated out to Eurasia replacing there archaic humans [1, 2] has continued to gain support from genetic contributions [3–6]. Anthropologically, the most ancient presence of modern humans out of Africa has been documented in the Levant about 95–125 kya [7, 8], and in Australia about 50–70 kya . Based on archaeological  and classic genetic studies [11, 12] two dispersals from Africa were proposed: A northern route that reached western and central Asia through the Near East, and a Southern route that, coasting Asia, reached Australia. However, ages for these dispersals were very tentative. The first phylogeographic analysis using complete mtDNA genomic sequences dated the out of Africa migrations around of 55–70 kya, when two branches, named M and N, of the African macrohaplogroup L3 radiation supposedly began the Eurasian colonization [5, 6]. A more recent analysis, based on a greater number of sequences, pushed back the lower bound of the out-of-Africa migration, signed by the L3 radiation, to around 85 kya . This date is no so far from the above commented presence of modern humans in the Levant about 100–125 kya. Interestingly, this migration is also in frame with the putative presence of modern humans in Eritrean coasts , and corresponds with an interglacial period (OIS 5), when African faunas expanded to the Levant . After that, it seems that, at least in the Levant, there was a long period of population bottleneck, as there is no modern human evidence in the area until 50 kyr later, again in a relatively warm period (OIS 3). This contraction phase might be reflected in the basal roots of M and N lineages by the accumulation of 4 and 5 mutations before their next radiation around 60 kya .
A total of 365 different mtDNA haplotypes were observed in 553 Saudi Arab sequences. 299 of them (82%) could have been detected using only the HVSI sequence information and 66 (18%) when the HVSII information was also taken into account. Additional analysis of diagnostic positions allowed the unequivocal assortment of the majority (96%) of the haplotypes into subhaplogroups [see Additional file 1]. However, 11 haplotypes were classified at the HV/R level, 3 assigned to macrohaplogroups L3*, M* and N* respectively, and only one was left unclassified [see Additional file 1]. The most probable origin of these Saudi haplotypes deserves a more detailed analysis.
Macrohaplogroup L lineages
Sub-Saharan Africa L lineages in Saudi Arabia account for 10% of the total. χ2 analyses showed that there is not significant regional differentiation in this Country. However, there is significant heterogeneity (p < 0.001) when all the Arabian Peninsula countries are compared. This is mainly due to the comparatively high frequency of sub-Saharan lineages in Yemen (38%) compared to Oman-Qatar (16%) and to Saudi Arabia-UAE (10%). Most probably, the higher frequencies shown in southern countries reflect their greater proximity to Africa, separated only by the Bab al Mandab strait. However, when attending to the relative contribution of the different L haplogroups, Qatar, Saudi Arabia and Yemen are highly similar for their L3 (34%), L2 (36%) and L0 (21%) frequencies whereas in Oman and UAE the bulk of L lineages belongs to L3 (72%). In this enlarged sample of Saudi Arabs, representatives of all the recently defined East African haplogroups L4 , L5 , L6  and L7 , have been found. The only L4 Saudi haplotype belongs to the L4a1 subclade defined by 16207T/C transversion. Although it has no exact matches its most related types are found in Ethiopia . Four L5 lineages have been found in Saudi Arabia but all have the same haplotype that belongs to the L5a1 subclade defined in the HVSI region by the 16355–16362 motif . It has matches in Egypt and Ethiopia. L6 was found the most abundant clade in Yemen . It has been now detected in Saudi Arabia but only once. This haplotype (16048-16223-16224-16243-16278-16311) differs from all the previous L6 lineages by the presence of mutation 16243. In addition it lacks the 16362 transition that is carried by all L6 lineages from Yemen but has the ancestral 16048 mutation only absent in one Yemeni lineage . This Saudi type adds L6 variability to Arabia, because until now L6 was only represented by a very abundant and a rare haplotype in Yemen. Attending to the most probable geographic origin of the sub-Saharan Africa lineages in Saudi Arabia, 33 (61%) have matches with East Africa, 7 (13%) with Central or West Africa whereas the rest 14 (26%) have not yet been found in Africa. Nevertheless, half of them belong to haplogroups with Western Africa origin and the other half to haplogroups with eastern Africa adscription [35, 30]. It is supposed that the bulk of these African lineages reached the area as consequence of slave trade, but more ancient historic contacts with northeast Africa are also well documented [36, 30, 31].
Macrohaplogroup M lineages
M lineages in Saudi Arabia account for 7% of the total. Half of them belong to the M1 African clade. There is no significant heterogeneity within Saudi Arabia regions nor among Arabian Peninsula countries for the total M frequency. However, when we compared the frequency of the African clade M1 against that of the other M clades of Asiatic provenance, it was significantly greater in western Arabian Peninsula regions than in the East (χ2 = 12.53 d.f. = 4 p < 0'05).
Inclusion of rare Saudi and other published African M1 sequences into the M1 genomic phylogenetic tree
Inclusion of rare Saudi Asiatic M sequences into the macrohaplogroup M tree
The Saudi sequence 201 deserves special mention (Figure 3). It was previously tentatively related to the Indian M34 clade because both share the 3010 transition. However, it was stated that due to the high recurrence of 3010 most probably the 201 sequence would belong to a yet undefined clade . The recent study of new Australian lineages  has allowed us to find out an interesting link between their Australian M14 lineage and our Saudi 201 sequence (Figure 3). The authors related M14 to the Melanesian clade M28  because both share the 1719–16148 motif . We think that the alternative motif shared with the Saudi lineage, 234-4216-6962, (Figure 3) is stronger, as 1719 and 16148 transitions are more recurrent than 234, 4216 and 6962 . Therefore, we think that the last three mutations defined the true root of the Australian M14 clade and relate it to a Saudi Arab sequence.
Macrohaplogroup N lineages
All the main non R West Eurasian haplogroups that directly spread from the root of macrohaplogroup N (N1a, N1b, N1c, I, W, X) were present in Saudi Arabia albeit in low frequencies [see Additional file 2]. Summing up, their total frequency only represents 12% of the Saudi pool. There is no heterogeneity among Saudi regions for the joint distribution of these haplogrups, which is in contrast with the high differentiation observed (χ2 = 13.9; d.f. 3; p < 0.01) when all the Arabian Peninsula countries were compared. However, this is mainly due to the high frequency in UAE of unclassified N* lineages . The three N1 subgroups have similar frequencies (2.4%) in Saudi Arabia [see Additional file 2], but haplotypic diversity (h) reach different levels being comparatively high for N1b (0.89) or N1a (0.83) and lower for N1c (0.63). N1a is also very diverse (0.89) in Yemen  but there are no haplotypic matches between both Arabian countries, and only one (16147G-16172-16223-16248-16355) of the 7 Saudi haplotypes has an exact match in Ethiopia. The high diversity of N1a in the Arabian Peninsula, Ethiopia and Egypt raises the possibility that this area was a secondary center of expansion for this haplogroup. However, the highest diversity for N1b and N1c are in Turkey, and Kurds and Iranians, respectively [see Additional file 3]. Only haplogroup N1c shows significant (χ2 = 12.8; d.f. 3; p < 0.01) regional distribution in Saudi being the highest frequencies in the Northern region [see Additional file 2]. Haplogroup I has been only detected in Central and Southeastern regions in Saudi Arabia and in all Arabian Peninsula countries with the exception of UAE . Haplogroup W has not been found in the western Saudi region nor in the Southern Arabian Yemen and Oman countries [30, 32]. Finally, haplogroup X is present in all the Saudi regions and in all Arabian countries excepting Oman [30, 32]. Two geographically well differentiated X branches can be distinguished by HVSII positions . The North African specific subclade X1 was defined by the 146 transition and the Eurasian specific subclade X2 by the 195 transition. As it was already found in Yemen , all the X haplotypes found in Saudi Arabia have the 195 transition, falling within the Eurasian branch, which discards East African introductions. Curiously, only the basic haplotype (16189-16223-16278) has matches with other regions.
Macrohaplogroup R lineages
Macrohaplogroup R is the main branch of N and their major subclades (H, J-T, K-U) embraced the majority of the West Eurasian mtDNA lineages. In Saudi Arabia R derivates represent a (70.5 ± 2.4) % of the total having a homogeneous regional distribution. This contrasts with the significant heterogeneity (χ2 = 46.1; d.f. 4; p < 0.001) found for the whole Arabian Peninsula, although it is due to the low frequency of R in Yemen (46.2%). The Western Asia haplogroup H is the most abundant haplogroup in Europe and the Near East . However in the Arabian Peninsula its mean frequency (9.4 ± 1.1) is moderate, reaching its highest value in Oman (13.3%) and its lowest, again, in Yemen (7.2%). Within Saudi Arabia, the highest frequencies are in the Northern and Central regions (9.3%) and the lowest in the West region (2.8%), being CRS, H2a1 and H6 the most abundant subclades which confirm other authors results . Haplogroup T shows regional heterogeneity in Saudi Arabia (χ2 = 10.4; d.f. 3; p < 0.05) and has significantly lower frequencies in Southern Yemen and Oman countries (χ2 = 5.8; p < 0.05). Furthermore, whereas subclade T3 is the most abundant in the Saudi Central region, subclades T1 and T5 are so in Northern and Western regions. Haplogroup U comprises numerous branches (U1 to U9 and K) that have different geographic distributions [47, 16, 17, 49]. In Saudi Arabia all of them have representatives albeit in minor frequencies, K (4%) and U3 (2.3%) being the most abundant clades. There is no geographical heterogeneity for the total U distribution in Saudi Arabia. Nevertheless, it is significantly different among the Arabian Peninsula countries (χ2 = 12.0; d.f. 4; p < 0.05), with Southern countries showing higher frequencies than the others [see Additional file 3]. Analysis of specific haplogroups reveals some interesting geographic distributions in the Arabian Peninsula. The European clade U2e and the rare clade U9 (χ2 = 10.5; d.f. 4; p < 0.05) could have reached the Arabian Peninsula from northern areas. On the contrary, the Indian clade U2 (χ2 = 34.5; d.f. 4; p < 0.001), U3, U4, U7 (χ2 = 11.7; d.f. 4; p < 0.05) and K (χ2 = 10.5; d.f. 4; p < 0.05), most probably came from the East. Finally, the North African clade U6, had a Western provenance. Haplogroup (preHV)1 was phylogenetically and phylogeographically studied in detail previously . However, this enlarged Saudi sample has allowed a regional analysis in Saudi Arabia. An heterogeneity test showed that (preHV)1 is significantly (χ2 = 8.5; d.f. 3; p < 0.05) more abundant in Northern (18.6%) and Central (21.8%) regions, and has particularly low frequency in the West region (8.3%). Extending the analysis to the whole Arabian Peninsula the heterogeneity grows considerably (χ2 = 30.9; d.f. 4; p < 0.001). This is mainly due to the comparatively lower frequencies in Yemen, Qatar and UAE [see Additional file 3]. Most probably, this haplogroup reached Arabia from the North , extending its geographic range to the South using mainly internal instead of coastal routes. As a whole, haplogroup J reaches its highest frequency in Saudi Arabia , where its regional distribution is also significantly heterogeneous (χ2 = 15.0; d.f. 3; p < 0.01), but opposite to that found for (preHV)1. For the J, the West (37.5%) and Southeast (25.7%) regions have higher frequencies than the Central (17.6%) and North (16.3%) regions. Heterogeneity in the whole Peninsula is also significant (χ2 = 16.5; d.f. 4; p < 0.01) being Saudi Arabia (21%) and Qatar (17.8%) the two countries with the highest J frequencies. However, the subclade distribution is different in each country. Subclade J1b is the main contributor (9.4%) in Saudi Arabia while other J subclades account for 14.5% in Qatar [see Additional file 3]. With the Qatar exception, J1b is the most frequent subclade in the Arabian Peninsula [see Additional file 3]. So for, it deserves more detailed phylogenetic and phylogeographic studies.
Phylogeny of haplogroup J1b
Phylogeography of haplogroup J1b
K haplotype diversity (in %) and π (× 1000) diversity for the total J1b, J1b1a1 and J1b1a1a subhaplogroups.
3.59 ± 2.65
7.91 ± 4.98
8.67 ± 5.23
6.32 ± 4.02
2.50 ± 2.33
3.78 ± 3.02
3.97 ± 2.84
1.59 ± 1.56
Population based comparisons
Eurasian and African influences in the Arabian Peninsula
Although until recent times the majority of the Saudi Arabia population was nomad, a moderate level of mitochondrial genetic structure has been found amongst its different regions. This heterogeneity grew considerably when all the Arabian Peninsula countries were included in the AMOVA analysis. It seems that the main cause of this diversity is the unequal influence that the different areas received from their geographically closest neighbors. This fact is graphically reflected in the MDS plot (Figure 6) where all the Arabian Peninsula samples are compared with samples from East Africa, the Near East and the Caucasus areas. The clustering of all the Saudi regions clearly shows that, in comparison to other geographically more distant populations, they form a rather homogenous entity, as was previously suggested from analysis based on classical markers . The more distant positions of Yemen, grouped with African samples, and the UAE, and in a lesser degree the Qatar and Oman, proximity to Near East countries, reflect their different frequencies of African and Eurasian lineages in their respective mitochondrial pools. Roughly, the African contribution to whole Arabian Peninsula accounts for 20% of its lineages if, in addition to all the L haplogroups, the North African M1 and U6 clades are added. However, the western and southern areas have received significantly stronger influences than the rest. Particularly, Yemen has the largest contribution of L lineages . So, most probably, this area was the entrance gate of a portion of these lineages in prehistoric times, which participated in the building of the primitive Arabian population. Later, received gene flows from North Africa and the Near East, and suffered expansions and retractions in humid or arid climatic periods. These fluctuations are also reflected in the frequent loss of diversity for several African clades as the L6 in Yemen  or the L5 in Saudi Arabia. However, the presence of western Africa L lineages and the different composition of L subclades in the African pool of different countries might reflect unequal participation of the primitive and the recent slave trade substrates in their respective African components.
An important group of the Arabian Peninsula lineages (18%), comprising representatives of the majority of the U clades, R2, and Central Asian, Indian, and Indonesian M lineages, seem to have their origins in the East, reaching the Arabian Peninsula through Iran where, in contrast to the Near East, the U clades (29%) have the highest frequency instead of the H (17%) group . Congruently, this Eastern gene flow had a significantly stronger impact in the Eastern and Southern areas of the Arabian Peninsula. However, the bulk of the Arabian N and R lineages (62%) had a Northern source. Haplogroups (preHV)1 and J1b were the main contributors of this gene flow. Nevertheless, its present day geographic distributions in the Arabian Peninsula are different. Whereas (preHV)1 presents significant higher frequencies in the North and Central Saudi regions and in Oman, J1b shows its highest frequencies in the more peripheral West and Southeast Saudi regions. It seems that at least haplogroups H, N1c and subclade T3 could have followed the (preHV)1 internal way of dispersion, while the T1 and T5 branches of haplogroup T and other branches of haplogroup J followed the peripheral route of clade J1b. Attending to the radiation ages of (preHV)1 and J1b clades and their derivate branches, striking similarities but also differences can be observed. The first expansion of both clades in the Near East had similar Paleolithic ages around 20,000 years ago. However whereas the ancestral HVSI motif of the (preHV)1 expansion was barely present in Saudi Arabia , the ancestral HVSI motif of the J1b radiation had an important incidence in that area (Figure 5) suggesting an active role in Arabia of the first J1b spread but not for that of (preHV)1. The succeeding most important radiations of both clades, (preHV)1a1 and J1b1a1 had, again, similar ages around 10,000 years that place them in Neolithic times. Now, in both cases, there is a shortage or absence of the ancestral motif in Arabia discarding this area as a radiation center. However, it participated in the (preHV)1a1 spread  but not in the J1b1a1 one (Figure 5). Finally, the third more abundant subclades, (preHV)1b rooted by 16304  and J1b rooted by 16136 (Figure 5) had the Arabian Peninsula as the most probable source of expansion. Nevertheless, whereas the J1b branch TMRCA (11,099 ± 8,381 years ago) was contemporary to that of the northern J1b1a1, the recalculated age of the (preHV)1b branch (by adding all the new HVSI sequences found in the present survey to the ones previously used ), was of only 4,036 ± 2,211 years ago which situates this expansion in the Bronze Age. These results could be satisfactorily explained if we admit an older Paleolithic implantation in Saudi Arabia of the J1b clade that, perhaps, with some other N and L clades would form the primitive population. Posterior (preHV)1 subclade radiations, accompanied by other clades, penetrated from the North using internal routes and even had secondary spreads in central Arabia diluting the J1b frequencies in these areas and causing its peripheral distribution.
Genomic dissection of rare M lineages
By genomic sequencing of seven M lineages (Accession numbers: EU370391–97), it has been demonstrated that the majority of the rare M lineages detected in Saudi Arabia (Figure 3) have Indian roots. However, the link found between the M Saudi 201 sequence and an M14 Australian sequence is puzzling. Although at first sight it could be taken as a signal of the connection between the two utmost ends of the southern route, it seems not to be the case. First, both lineages share three basal positions and this hypothetical link would considerably delay the arrival age of M in comparison to that of East Asia. It would be improbable that similar Australian links with other M lineages mainly from India were not found. Third, if the Arab lineage had such an old implantation in the Arabian Peninsula some detectable autochthonous radiation should be expected. Most probably, the M42 sequence belongs to an Australian clade and its related lineage found in Saudi Arabia is also of Australian origin. Historical links as those invoked to explain the presence of Indian and Indonesian sequences in the Arabian Peninsula pool should also be valid for this case. In our opinion, the camel trade between Saudi Arabia and Australia  could be a probable historic cause of this link. Future detection in Aboriginal Australians of other M42 lineages will confirm the Australian origin of this clade and its radiation age in that Continent. However, the link between the East Asia M10 clade  and the Australian M42 clade, if not due to convergence, seems to be more interesting as it would confirm, once more, the rapid expansion of macrohaplogroup M all along the Asian coasts [6, 13]. The lack of autochthonous M and N lineages in the present day Arabian Peninsula populations confirms that this area was not a place of demographic expansion in the dispersal out of Africa .
Although there is evidence of Neolithic and more recent expansions in the Arabian Peninsula, mainly detected by (preHV)1 and J1b lineages, the lack of primitive autochthonous M and N sequences, suggests that this area has been more a receptor of human migrations, including historic ones, from Africa, India, Indonesia and even Australia, than a demographic expansion center along the proposed southern coastal route.
Buccal swabs or peripheral blood were obtained from 553 (120 of them previously published in Abu-Amero et al. ) maternally unrelated Saudi Arabs all whose known ancestors were of Saudi Arabian origin. The main Saudi Arabian geographic regions were sampled (Figure 1 and Additional file 1). Sequence analysis was performed of mtDNA regulatory region hypervariable segment I (HVSI) and hypervariable segment II (HVSII) and of haplogroup diagnostic mutations using RFLPs or partial sequencing [see Additional file 1]. In addition, genomic mtDNA sequencing was carried out in 7 individuals of uncertain or interesting haplogroup adscription. For population and phylogeographic comparison, we used 21,808 published or unpublished partial sequences from Europe (11,174), South Asia (2,746), Caucasus (1,638), North Africa (1,009), East Africa (888), Near East (2,001), Arabian Peninsula (1,129) and Jews (1,223), as detailed in Additional files 4 and 5. Informed consent was obtained from all individuals.
Total DNA was isolated from buccal and blood samples using the PUREGENE DNA isolation kit from Gentra Systems (Minneapolis, USA). HVSI and HVSII segments were PCR amplified using primers pairs L15840/H16401 and L16340/H408, respectively, as previously described . Genomic mtDNA sequences and segments including diagnostic positions were amplified using a set of 32 separate PCRs and cycling conditions as detailed elsewhere . Successfully amplified products were sequenced for both complementary strands using the DYEnamic™ ET dye terminator kit (Amersham Biosciences), and samples were run on MegaBACE 1000 (Amersham Biosciences) according to the manufacturer protocol.
Classification into sub-haplogroups was performed using nomenclature previously described for African [35, 30, 34] and for Eurasian [47, 49, 40, 24, 20, 37, 31, 26] sequences. Published sequences employed for comparative genetic analysis were re-classified into sub-clades using the same criteria in order to permit comparisons.
Haplotype diversity was calculated as h  and as K (haplotype number/sample size quotient). Only HVSI positions from 16069 to 16385 were used for genetic comparisons of partial sequences with other published data. Genetic variation was apportioned within and among geographic regions using AMOVA by means of ARLEQUIN2 . Four regions (North, Central, West and South-East) were considered to assess the Saudi Arabia genetic structure (Figure 1 and Additional file 2). For more extended geographic comparisons the following areas were taken into account: Arabian Peninsula (including Saudi Arabia, Qatar, UAE, Oman, Yemen and Bedouin Arabs), North-East Africa (including samples from Egypt, Nubian, Sudan, Ethiopia, and Kenya), and Near East (containing samples from Druze, Iran, Iraq, Jordan, Kurds, Palestine, Syria and Turkey), as detailed in Additional file 3. Pairwise FST distances between populations were calculated from haplogroup and haplotype frequencies, and their significance assessed by a nonparametric permutation test (ARLEQUIN2). Multidimensional scaling (MDS) plots were obtained with SPSS version 13.0 (SPSS Inc., Chicago, Illinois). Phylogenetic relationships among HVSI and genomic mtDNA sequences were established using the reduced median network algorithm . In addition to our 7 genomic mtDNA sequences, 7, 12 and 23 published complete or nearly complete mtDNA sequences were used to establish the M1 (Figure 2), M (Figure 3) and J1b (Figure 4) phylogenies, respectively.
Only substitutions in the coding region were taken into account for complete sequences, excluding insertions and deletions. The mean number of substitutions per site compared to the most common ancestor (ρ) of each clade and its standard error were calculated following Morral et al.  and Saillard et al.  respectively, and converted into time using previously published substitution rates [5, 52]. For HVSI, the age of clusters or expansions was calculated as the mean divergence (ρ) from inferred ancestral sequence types [59, 60] and converted into time by assuming that one transition within np 16090–16365 corresponds to 20,180 years .
The seven new complete mitochondrial DNA sequences are registered under GenBank accession numbers: EU370391–EU370397
This research was supported by grant BFU2006-04490 from Ministerio de Ciencia y Tecnología to J.M.L. We would like to thank Rebecca S. Just and coauthors for allowing us the use of their unpublished and valuable data.
- Bräuer G: A craniological approach to the origin of anatomically modern Homo sapiens in Africa and implications for the appearance of modern Europeans. The origins of Modern Humans: A World Survey of the Fossil Evidence. Edited by: Stringer CB, Mellars. 1984, New York: P Alan R. Liss, 327-410.Google Scholar
- Stringer CB, Andrews T: Genetic and fossil evidence for the origin of modern humans. Science. 1988, 239: 1263-1268. 10.1126/science.3125610.View ArticlePubMedGoogle Scholar
- Cann RL, Stoneking M, Wilson AC: Mitochondrial DNA and human evolution. Nature. 1987, 325: 31-36. 10.1038/325031a0.View ArticlePubMedGoogle Scholar
- Vigilant L, Stoneking M, Harpending H, Hawkes K, Wilson AC: African populations and the evolution of human mitochondrial DNA. Science. 1991, 253: 1503-1507. 10.1126/science.1840702.View ArticlePubMedGoogle Scholar
- Ingman M, Kaessmann H, Paabo S, Gyllensten U: Mitochondrial genome variation and the origin of modern humans. Nature. 2000, 408: 708-713. 10.1038/35047064. Erratum in: Nature 2001, 410: 611View ArticlePubMedGoogle 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: 13-10.1186/1471-2156-2-13.PubMed CentralView ArticlePubMedGoogle Scholar
- Valladas H, Reyss J, Joron J, Valladas G, Bar-Yosef O, Vandermeersch B: Thermolumininescence datin of Mousterian Proto-Cro-Magnon remains from Israel and the origin of modern man. Nature. 1988, 331: 614-616. 10.1038/331614a0.View ArticleGoogle Scholar
- Mercier N, Valladas H, Bar-Yosef O, Vandermeersch B, Stringer C, Joron J-L: Thermoluminescence date for the Mousterian burial site of Es Skhul, Mt. Carmel. Journal of Archaeological Science. 1993, 20: 169-174. 10.1006/jasc.1993.1012.View ArticleGoogle Scholar
- Thorne A, Gruën R, Mortimer G, Spooner NA, Simpson JJ, McCulloch M, Taylor L, Curnoe D: Australia's oldest human remains : age of the Lake Mungo 3 skeleton. J Hum Evol. 1999, 36: 591-612. 10.1006/jhev.1999.0305.View ArticlePubMedGoogle Scholar
- Foley RA, Lahr MM: Mode 3 technologies and the evolution of modern humans. Cambridge Archaeological Journal. 1997, 19: 3-36.View ArticleGoogle Scholar
- Nei M, Roychoudhury AK: Evolutionary relationships of human populations on a global scale. Mol Biol Evol. 1993, 10: 927-943.PubMedGoogle Scholar
- Cavalli-Sforza LL, Menozzi P, Piazza A: The history and geography of human genes. 1994, Princeton: University PressGoogle Scholar
- Macaulay V, Hill C, Achilli A, Rengo C, Clarke D, Meehan W, Blackburn J, Semino O, Scozzari R, Cruciani F, Taha A, Shaari NK, Raja JM, Ismail P, Zainuddin Z, Goodwin W, Bulbeck D, Bandelt HJ, Oppenheimer S, Torroni A, Richards M: Single, rapid coastal settlement of Asia revealed by analysis of complete mitochondrial genomes. Science. 2005, 308: 1034-1036. 10.1126/science.1109792.View ArticlePubMedGoogle Scholar
- Walter RC, Buffler RT, Bruggemann JH, Guilaume MMM, Berhe SM, Negassi B, Libsekal Y, Cheng H, Edwards RL, von Cosel R, Néraudeau D, Gagnon M: Early human occupation of the Red Sea coast of Eritrea during the last interglacial. Nature. 2000, 405: 65-69. 10.1038/35011048.View ArticlePubMedGoogle Scholar
- Tchernov E: Biochronology, paleoecology, and dispersal events of hominids in the southern Levant. The evolution and dispersal of modern humans in Asia. Edited by: Akazawa T, Aoiki K, Kimura T. 1992, Tokyo: Hokusen-sha, 149-188.Google Scholar
- Kivisild T, Rootsi S, Metspalu M, Mastana S, Kaldma K, Parik J, Metspalu E, Adojaan M, Tolk HV, Stepanov V, Golge M, Usanga E, Papiha SS, Cinnioglu C, King R, Cavalli-Sforza L, Underhill PA, Villems R: The genetic heritage of the earliest settlers persists both in Indian tribal and caste populations. Am J Hum Genet. 2003, 72: 313-332. 10.1086/346068.PubMed CentralView ArticlePubMedGoogle Scholar
- Metspalu M, Kivisild T, Metspalu E, Parik J, Hudjashov G, Kaldma K, Serk P, Carmin M, Behar DM, Gilbert MTP, Endicott P, Mastana S, Papiha SS, Skorecki K, Torroni A, Villems R: 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: 26-10.1186/1471-2156-5-26.PubMed CentralView ArticlePubMedGoogle Scholar
- Palanichamy MG, Sun C, Agrawal S, Bandelt HJ, Kong QP, Khan F, Wang CY, Chaudhuri TK, Palla V, Zhang YP: 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-978. 10.1086/425871.PubMed CentralView ArticlePubMedGoogle Scholar
- Rajkumar R, Banerjee J, Gunturi HB, Trivedi R, Kashyap VK: Phylogeny and antiquity of M macrohaplogroup inferred from complete mtDNA sequence of Indian specific lineages. BMC Evol Biol. 2005, 5: 26-10.1186/1471-2148-5-26.PubMed CentralView ArticlePubMedGoogle Scholar
- Sun C, Kong QP, Palanichamy MG, Agrawal S, Bandelt HJ, Yao YG, Khan F, Zhu CL, Chaudhuri TK, Zhang YP: The dazzling array of basal branches in the mtDNA macrohaplogroup M from India as inferred from complete genomes. Mol Biol Evol. 2006, 23: 683-690. 10.1093/molbev/msj078.View ArticlePubMedGoogle Scholar
- Thangaraj K, Chanbey G, Singh VK, Canniarajan A, Thanseem I, Reddy AG, Singh L: In situ origin of deep rooting lineages on Indian mitochondrial macrohaplogroup M in India. BMC Genomics. 2006, 7: 151-10.1186/1471-2164-7-151.PubMed CentralView ArticlePubMedGoogle Scholar
- Ingman M, Gyllensten U: Mitochondrial genome variation and evolutionary history of Australian and New Guinean aborigines. Genome Res. 2003, 13: 1600-1606. 10.1101/gr.686603.PubMed CentralView ArticlePubMedGoogle Scholar
- Friedlaender J, Schurr T, Gentz F, Koki G, Friedlaender F, Horvat G, Babb P, Cerchio S, Kaestle F, Schanfield M, Deka R, Yanagihara R, Merriwether DA: Expanding Southwest Pacific mitochondrial haplogroups P and Q. Mol Biol Evol. 2005, 22: 1506-17. 10.1093/molbev/msi142. Erratum in: Mol Biol Evol 2005, 22: 2313View ArticlePubMedGoogle 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 USA. 2005, 102: 13034-13039. 10.1073/pnas.0506195102. Erratum in: Proc Natl Acad Sci USA 2005, 102: 16904PubMed CentralView ArticlePubMedGoogle 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-294. 10.1002/ajpa.20426.View ArticlePubMedGoogle Scholar
- Hudjashov G, Kivisild T, Underhill PA, Endicott P, Sanchez JJ, Lin AA, Shen P, Oefner P, Renfrew C, Villems R, Forster P: Revealing the prehistoric settlement of Australia by Y chromosome and mtDNA analysis. PNAS. 2007, 104: 8726-8730. 10.1073/pnas.0702928104.PubMed CentralView ArticlePubMedGoogle Scholar
- Forster P, Torroni A, Renfrew C, Röhl A: Phylogenetic star contraction applied to Asian and Papuan mtDNA evolution. Mol Biol Evol. 2001, 18: 1864-1881.View ArticlePubMedGoogle Scholar
- Forster P: Ice Ages and the mitochondrial DNA chronology of human dispersals: a review. Phil Trans R Soc Lond B. 2004, 359: 255-264. 10.1098/rstb.2003.1394.View ArticleGoogle Scholar
- Oppenheimer S: Out of Eden: the peopling of the world. Constable, London. Republished as: The real Eve: modern man's journey out of Africa. 2003, New York: Carroll and GrafGoogle 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-770. 10.1086/425161.PubMed CentralView ArticlePubMedGoogle Scholar
- Abu-Amero KK, González AM, Larruga JM, Bosley TM, Cabrera VM: Eurasian and African mitochondrial DNA influences in the Saudi Arabian population. BMC Evolutionary Biology. 2007, 7: 32-10.1186/1471-2148-7-32.PubMed CentralView ArticlePubMedGoogle Scholar
- Rowold DJ, Luis JR, Terreros MC, Herrera RJ: Mitochondrial DNA geneflow indicates preferred usage of the Levant Corridor over the Horn of Africa passageway. J Hum Genet. 2007, 52: 436-447. 10.1007/s10038-007-0132-7.View ArticlePubMedGoogle Scholar
- Shen P, Lavi T, Kivisild T, Chou V, Sengun D, Gefel D, Shpirer I, Woolf E, Hillel J, Feldman MW, Oefner PJ: Reconstruction of patrilineages and matrilineages of Samaritans and other Israeli populations from Y-chromosome and mitochondrial DNA sequence variation. Hum Mutation. 2004, 24: 248-260. 10.1002/humu.20077.View ArticleGoogle Scholar
- Torroni A, Achilli A, Macaulay V, Richards M, Bandelt H-J: Harvesting the fruit of the human mtDNA tree. Trends Genet. 2006, 22: 339-345. 10.1016/j.tig.2006.04.001.View ArticlePubMedGoogle Scholar
- Salas A, Richards M, De la Fe T, Lareu MV, Sobrino B, Sanchez-Diz P, Macaulay V, Carracedo A: The making of the African mtDNA landscape. Am J Hum Genet. 2002, 71: 1082-1111. 10.1086/344348.PubMed CentralView ArticlePubMedGoogle Scholar
- Richards M, Rengo C, Cruciani F, Gratrix F, Wilson Jf, Scozzari R, Macaulay V, Torroni A: Extensive female-mediated gene flow from sub-Saharan Africa into near eastern Arab populations. Am J Hum Genet. 2003, 72: 1058-1064. 10.1086/374384.PubMed CentralView ArticlePubMedGoogle Scholar
- Olivieri A, Achilli A, Pala M, Battaglia V, Fornarino S Al-Zahery N, Scozzari R, Cruciani F, Behar DM, Dugoujon JM, Coudray C, Santachiara-Benerecetti AS, Semino O, Bandelt HJ, Torroni A: The mtDNA legacy of the Levantine early Upper Palaeolithic in Africa. Science. 2006, 314: 1767-1770. 10.1126/science.1135566.View ArticlePubMedGoogle Scholar
- González AM, Larruga JM, Abu-Amero KK, Shi Y, Pestano J, Cabrera VM: Mitochondrial lineage M1 traces an early human backflow to Africa. BMC Genomics. 2007, 8: 223-10.1186/1471-2164-8-223.PubMed CentralView ArticlePubMedGoogle Scholar
- Gonder MK, Mortensen HM, Reed FA, de Sousa A, Tishkoff SA: Whole-mtDNA Genome Sequence Analysis of Ancient African Lineages. Mol Biol Evol. 2007, 24: 757-768. 10.1093/molbev/msl209.View ArticlePubMedGoogle Scholar
- Tanaka M, Cabrera VM, González AM, Larruga JM, Takeyasu T, Fuku N, Guo LJ, 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-1850. 10.1101/gr.2286304.PubMed CentralView ArticlePubMedGoogle Scholar
- Comas D, Plaza S, Wells RS, Yuldaseva N, Lao O, Calafell F, Bertranpetit J: Admixture, migrations, and dispersals in Central Asia: evidence from maternal DNA lineages. Eur J Hum Genet. 2004, 12: 495-504. 10.1038/sj.ejhg.5201160.View ArticlePubMedGoogle Scholar
- Tommaseo-Ponzetta M, Attimonelli M, De Robertis M, Tanzariello F, Saccone C: Mitochondrial DNA variability of West New Guinea populations. Am J Phys Anthropol. 2002, 117: 49-67. 10.1002/ajpa.10010.View ArticlePubMedGoogle Scholar
- Husson L: Indonesians in Saudi Arabia for worship and work. Rev Eur Migr Int. 1997, 13: 125-147.View ArticlePubMedGoogle Scholar
- Kivisild T, Shen P, Wall DP, Do B, Sung R, Davis K, Passarino G, Underhill PA, Scharfe C, Torroni A, Scozzari R, Modiano D, Coppa A, de Knijff P, Feldman M, Cavalli-Sforza LL, Oefner PJ: The role of selection in the evolution of human mitochondrial genomes. Genetics. 2006, 172: 373-387. 10.1534/genetics.105.043901.PubMed CentralView ArticlePubMedGoogle Scholar
- MITOMAP. [http://www.mitomap.org]
- Reidla M, Kivisild T, Metspalu E, Kaldma K, Tambets K, Tolk HV, Parik J, Loogvali EL, Derenko M, Malyarchuk B, Bermisheva M, Zhadanov S, Pennarun E, Gubina M, Golubenko M, Damba L, Fedorova S, Gusar V, Grechanina E, Mikerezi I, Moisan JP, Chaventre A, Khusnutdinova E, Osipova L, Stepanov V, Voevoda M, Achilli A, Rengo C, Rickards O, De Stefano GF, Papiha S, Beckman L, Janicijevic B, Rudan P, Anagnou N, Michalodimitrakis E, Koziel S, Usanga E, Geberhiwot T, Herrnstadt C, Howell N, Torroni A, Villems R: Origin and diffusion of mtDNA haplogroup X. Am J Hum Genet. 2003, 73: 1178-1190. 10.1086/379380.PubMed CentralView ArticlePubMedGoogle Scholar
- Richards M, Macaulay V, Hickey E, Vega E, Sykes B, Guida V, Rengo C, Sellitto D, Cruciani F, Kivisild T, Villems R, Thomas M, Rychkov S, Rychkov O, Rychkov Y, Gölge M, Dimitrov D, Hill E, Bradley D, Romano V, Cali F, Vona G, Demaine A, Papiha S, Triantaphyllidis C, Stefanescu G, Hatina J, Belledi M, Di Rienzo A, Novelletto A, Oppenheim A, Norby S, Al-Zaheri N, Santachiara-Benerecetti S, Scozari R, Torroni A, Bandelt HJ: Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet. 2000, 67: 1251-1276. 10.1086/321197.PubMed CentralView ArticlePubMedGoogle Scholar
- Roostalu U, Kutuev I, Loogväli E-L, Metspalu E, Tambets K, Reidla M, Khusnutdinova EK, Usanga E, Kivisild T, Villems R: Origin and expansion of haplogroup H, the dominant human mitochondrial DNA lineage in West Eurasia: the Near Eastern and Caucasian perspective. Mol Biol Evol. 2006, 24: 436-448. 10.1093/molbev/msl173.View ArticlePubMedGoogle Scholar
- Quintana-Murci L, Chaix R, Wells RS, Behar DM, Sayar H, Scozzari R, Rengo C, Al-Zahery N, Semino O, Santachiara-Benerecetti AS, Coppa A, Ayub Q, Mohyuddin A, Tyler-Smith C, Qasim Mehdi S, Torroni A, McElreavey K: Where west meets east: the complex mtDNA landscape of the southwest and Central Asian corridor. Am J Hum Genet. 2004, 74: 827-845. 10.1086/383236.PubMed CentralView ArticlePubMedGoogle Scholar
- Finnilä S, Lehtonen MS, Majamaa K: Phylogenetic network for European mtDNA. Am J Hum Genet. 2001, 68: 1475-1484. 10.1086/320591.PubMed CentralView ArticlePubMedGoogle Scholar
- Forster P, Romano V, Calì F, Röhl A, Hurles M: mtDNA markers for Celtic and Germanic Language Areas in the British Isles. Traces of ancestry. Edited by: Martin Jones McDonald Institute Monographs. 2004, Published by McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, UK, 68: 99-114.Google Scholar
- Mishmar D, Ruiz-Pesini E, Golik P, Macaulay V, Clark AG, Hosseini S, Brandon M, Easley K, Chen E, Brown MD, Sukernik RI, Olckers A, Wallace DC: Natural selection shaped regional mtDNA variation in humans. P Natl Acad Sci USA. 2003, 100: 171-176. 10.1073/pnas.0136972100.View ArticleGoogle Scholar
- Di Rienzo A, Wilson AC: Branching pattern in the evolutionary tree for human mitochondrial DNA. Proc Nat Acad Sci USA. 1991, 88: 1597-1601. 10.1073/pnas.88.5.1597.PubMed CentralView ArticlePubMedGoogle Scholar
- Clark A: Camels Down Under. 1988, Saudi Aramco World, 16-23.Google Scholar
- Field JS, Lahr MM: Assessment of the Southern Dispersal: GIS-based analyses of potential routes at oxygen isotopic stage 4. J World Prehistory. 2006, 19: 1-45. 10.1007/s10963-005-9000-6.View ArticleGoogle Scholar
- Nei M: Molecular evolutionary genetics. 1987, New York: Columbia University PressGoogle Scholar
- Schneider S, Roessli D, Excoffier L: Arlequin ver. 2.000: A software for population genetics data analysis. Genetics and Biometry Laboratory. Switzerland: University of Geneva
- Bandelt H-J, Forster P, Rohl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.View ArticlePubMedGoogle Scholar
- Morral N, Bertranpetit J, Estivill X, Nunes V, Casals T, Giménez J, Reis A, Varon-Mateeva R, Macek M, Kalaydjieva L: The origin of the major cystic fibrosis mutation (ΔF508) in European populations. Nat Genet. 1994, 7: 169-175. 10.1038/ng0694-169.View ArticlePubMedGoogle Scholar
- Saillard J, Forster P, Lynnerup N, Bandelt HJ Norby S: MtDNA variation among Greenland Eskimos: the edge of the Beringian expansion. Am J Hum Genet. 2000, 67: 718-726. 10.1086/303038.PubMed CentralView ArticlePubMedGoogle 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-945.PubMed CentralPubMedGoogle Scholar
- Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, Schreier PH, Smith AJ, Staden R, Young IG: Sequence and organisation of the human mitochondrial genome. Nature. 1981, 290: 457-465. 10.1038/290457a0.View ArticlePubMedGoogle 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-10.1038/13779.View ArticlePubMedGoogle 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-1171. 10.1086/339933. Erratum in: Am J Hum Genet 2002, 71: 448–449PubMed CentralView ArticlePubMedGoogle Scholar
- Coble MD, Just RS, O'Callaghan JE, Letmanyi IH, Peterson CT, Irwin JA, Parsons TJ: Single nucleotide polymorphisms over the entire mtDNA genome that increase the power of forensic testing in Caucasians. Int J Legal Med. 2004, 118: 137-146. 10.1007/s00414-004-0427-6.View ArticlePubMedGoogle Scholar
- Rose G, Passarino G, Carrieri G, Altomare K, Greco V, Bertolini S, Bonafe M, Franceschi C, De Benedictis G: Paradoxes in longevity: sequence analysis of mtDNA haplogroup J in centenarians. Eur J Hum Genet. 2001, 9: 701-707. 10.1038/sj.ejhg.5200703.View ArticlePubMedGoogle Scholar
- Esteitie N, Hinttala R, Wibom R, Nilsson H, Hance N, Naess K, Tear-Fahnehjelm K, von Dobeln U, Majamaa K, Larsson NG: Secondary metabolic effects in complex I deficiency. Ann Neurol. 2005, 58: 544-552. 10.1002/ana.20570.View ArticlePubMedGoogle Scholar
- Just RS, Diegoli TM, Saunier JL, Irwin JA, Parsons TJ: Complete mitochondrial genome sequences for 265 African American and U.S. "Hispanic" individuals. Forensic Sci Intl: Genet.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.