Skip to main content

The mitochondrial lineage U8a reveals a Paleolithic settlement in the Basque country

Abstract

Background

It is customary, in population genetics studies, to consider Basques as the direct descendants of the Paleolithic Europeans. However, until now there has been no irrefutable genetic proof to support this supposition. Even studies based on mitochondrial DNA (mtDNA), an ideal molecule for constructing datable maternal genealogies, have failed to achieve this. It could be that incoming gene flow has replaced the Basque ancient lineages but it could also be that these lineages have not been detected due to a lack of resolution of the Basque mtDNA genealogies. To assess this possibility we analyzed here the mtDNA of a large sample of autochthonous Basques using mtDNA genomic sequencing for those lineages that could not be unequivocally classified by diagnostic RFLP analysis and control region (HVSI and HVSII) sequencing.

Results

We show that Basques have the most ancestral phylogeny in Europe for the rare mitochondrial subhaplogroup U8a. Divergence times situate the Basque origin of this lineage in the Upper Palaeolithic. Most probably, their primitive founders came from West Asia. The lack of U8a lineages in Africa points to an European and not a North African route of entrance. Phylogeographic analysis suggest that U8a had two expansion periods in Europe, the first, from a south-western area including the Iberian peninsula and Mediterranean France before 30,000 years ago, and the second, from Central Europe around 15,000–10,000 years ago.

Conclusion

It has been demonstrated, for the first time, that Basques show the oldest lineages in Europe for subhaplogroup U8a. Coalescence times for these lineages suggest their presence in the Basque country since the Upper Paleolithic. The European U8 phylogeography is congruent with the supposition that Basques could have participated in demographic re-expansions to repopulate central Europe in the last interglacial periods.

Background

Considering Basques as the direct descendants of Paleolithic Europeans has become a multidisciplinary premise. However, there is no irrefutable evidence for this supposition. The Basque country has a well represented archaeological record in Paleolithic and Mesolithic periods [1] but Archaeology can seldom differentiate an "in situ" cultural evolution from successive waves of new incomers. Basques speak a non Indo-European language with no close affinities with any other extant language but, even if their roots could be found, they would not reach the Paleolithic deepness due to the fast rate of change of languages. Classical population genetic studies, showed the Basques as one of the major outliers in Europe [2]. Nevertheless, these results can be explained by genetic drift which implies isolation but not necessarily an old history for that population. Lack of recombination and the fast mutation rate made mtDNA the ideal molecule to construct maternal genealogies, which frame in time and space the evolution and dispersion of human populations. However, until now, mtDNA studies on the Basques have only confirmed its low genetic diversity in a common Western Europe background [35]. It has been proposed, on the basis of their geographic distributions, that several mitochondrial lineages as V [6, 7], and the H1 and H3 subgroups [4, 5, 8] are markers of a Paleolithic human dispersal from southwestern Europe, including the Basque country, to Northeast Europe. However, diversities for these lineages are not higher in Basques than in Central Europeans. It could be possible that this lack of distinctness in Basques is real, in fact, even small levels of gene flow during enough time might have replaced the majority of their ancient lineages [9], but it could also be possible that this uniformity is due to a lack of resolution of the Basque mtDNA genealogies [10].

To deal with this possibility, we analyzed a sample of 211 unrelated Basques using the hypervariable segment of the mtDNA control region (HVSI-II) and diagnostic RFLP analysis, and sequenced the complete mitochondrial DNA of the rare lineages.

Results and discussion

The analysis of the Basque sample showed three haplotypes (CRS, 16342, and 16278 16311) that by their mutated positions in the control region had uncertain subhaplogroup adscription, but that by diagnostic RFLPs (+12308 Hinf I) belonged to haplogroup U/K. Also, we found one individual (16146 16189 16342) that belongs to the scarce U8a subhaplogroup. Complete sequencing of the four rare U haplotypes and their inclusion in a phylogenetic tree [11] with all published U complete sequences (data not shown) allowed their correct subhaplogroup affiliation. A more schematic tree (Fig.1) shows that one of the lineages (Bq24) belongs to subhaplogroup K1a1 [12], being a back-mutation of the diagnostic position 16224 its main peculiarity. Its most related K1a1 complete sequences are the eleven found by Herrnstadt et al. [13] and Finn 153 in Finnilä et al. [14]. Applying a mutation-rate of 1.26 × 10-8 [15] to their average sequence divergence [16] a radiation upper bound of 12 ± 4 Ky is obtained for this group. The other three lineages clustered into the rare subhaplogroup U8a [17]. This subhaplogroup can be RFLP diagnosed as -7055 Alu I. What is outstanding of these sequences is their great genetic diversity that extends the range of all known U8a European sequences (Fig.1). This Basque diversity specially contrasts with the lack of variation in the Finn sequences. Furthermore, the phylogenetic radiation of their U8a lineages (Fig. 1 and 2) is characteristic of an old population without recent exponential growth [18]. In fact, the most ancestral sequence (Bq1820) indicates that U8a lineages could have been in the Basque country since 28 ± 9 Ky, and that the other Basque lineages, belonging to the U8a1 subgroup (RFLP diagnosed as -3737 Hph I and + 5235 MspA1 I), participated in a more recent European expansion around 13 ± 5 Ky, similar to that estimated for K1a, and congruent with a re-expansion from an Iberian refuge when glaciers retreated in Europe proposed for other mtDNA clades [4, 5, 7]. Although all the U8a complete sequences belong to Europeans, the ancestral radiation of haplogroup U most probably occurred in western Asia shortly after the out of Africa episode [19], with early branch expansions to India (U2), Europe (U5) and North Africa (U6). U8 may be considered another main branch with a broad geographic range. Its first split separated U8a from U8b/K around 57 ± 11 Ky. Relatively short in time a new subdivision gave the sister clades U8b and K [17]. Until now there was only one completely sequenced U8b subject [20]. The addition of our Jordan 767 sequence to the tree (Fig. 1) gives a branching point for U8b, defined by transitions 6546, 6599 and 12771. Two of them, 6546 and 12771, can be MnlI-RFLP-detected. Curiously, the similar number of substitutions to the coalescence point of the three U8 branches U8a (5), U8b(5), K(6), suggests that all of them radiated at a similar age, supporting the hypothesis that, most probably, global climatic changes favored human expansions simultaneously at a continental scale [21]. When only RFLPs and/or partial sequence data are available, U8a haplotypes in general can be identified by the -7055 Alu I RFLP or the 73, 282 HVSII motif, and the majority of U8a1 derivates in particular by the16146, 16342 HVSI motif. In turn, U8b can be classified by the 16189, 16234 HVSI motif. From a total of 20,563 sequences studied 10,677 could be, unequivocally, U8a analyzed (see Additional file 1) and referenced (see Additional file 2). For the remaining (9,886), in general, only their U8a1 assignation was possible. Analysis of 19,133 Eurasian, and 1,430 North African published and unpublished RFLP/HVSI/HVSII sequences showed a scattered but widespread U8a/1 distribution that is restricted to Europe. Its frequency (See Additional file 1 for the European distribution of subhaplogroup U8a, and Additional file 2 for references) ranges from 0 in the majority of samples to 8% (although with a 95% coefficient range around 2.9–21.4) in the region of Var in Southeast France [22]. In spite of its moderate sample size an important characteristic of this region is its high U8a polymorphism as all the three lineages detected are different. The U8a eastern boundaries seem to be in the Volga region near the Urals [23]. U8b is also a quantitatively minor clade that partially overlaps with U8a in Europe. However, its presence in the Caucasus [24], Iran [17], the Near East [25], and North Africa [26], where U8a has not been detected, attests a more southern geographic distribution. The third sister clade K, is the most widespread and abundant covering the U8a and U8b ranges [24] and even reaching India [27]. A network [28] built with the 48 U8a sequences found (Fig. 2), could be rooted and resolved attending to the phylogeny of the U8 complete sequences. Its most ancestral node is represented by Bq1820 and the only Anatolian lineage assignable to haplogroup U8a, both carrying the CRS motif in HVSI, and transitions 73, 282 in HVSII. This ancient connection might trace the hypothetic route followed by the U8a ancestor from West Asia to the Basque country. The absence of U8a in North-Africa and its extremely rare presence in the eastern Mediterranean area further reinforces this continental route of entrance against a southern alternative. It is deduced from the network (Fig. 2) that a first U8a radiation in Europe affected Iberia, Central Europe and reached the Baltic. A second, U8a1, broader expansion further enlarged its range to Russia and Scotland where the U8a diversities are lower than in the central area (Table 1). As the total rooted network does not achieve the star genealogy, to calculate coalescence ages [29], we estimated the average distance and coalescence age to the U8a1 founder haplotype (16146 16342) and to the ancestral U8a haplotype (CRS) independently, representing, in the last case, the U8a1 radiation by only one basic lineage. Time estimations for the younger U8a1 expansion was 14 ± 5 Ky and 23 ± 14 Ky for the U8a subset that, if added to first, would give a total time for the U8a coalescence of around 37 ± 14 Ky. Notice that both HVSI estimations are higher than, but not significantly different, from those calculated using complete sequences. To compare the U8a diversity (p ± σ) among regions, we grouped the European populations in different Paleolithic areas [24, 30]. The greatest diversity was found in the Iberian Peninsula when Basques are included, followed by the North Central area (Table 1). These data agree with the primary and secondary origins of expansion proposed on phylogenetic grounds, weakening the possibility that Basques would have obtain their total U8a diversity through recent immigrations from other European areas and reinforcing the hypothesis that the first U8a radiation in Europe happened in an area in which the Basque country was included.

Figure 1
figure 1

Phylogenetic tree based on complete U8 sequences. Numbers along links refer to nucleotide positions. Open boxes are nodes from which other (not shown) sequences branch. A, C, indicate transversions; "d" deletions and "i" insertions. Star has the following mutations with respect to rCRS: 73, 263, 311i, 750, 1438, 2706, 4769, 7028, 8860, 11719, 14766 and 15326 and the following ones respect to L3*: 8701, 9540, 10398, 10873, 12705, 15301, 16223, 16519. Subject origins are: Dutch (Dut86; [38]), Spanish (Spa) and Italian (Ita) [20], 3 Finns (Fin; [14]). British-Australian (BriT7; Obayashi T, Tanaka M, personal communication), Jordan (Jor767) and 4 Basques (Bq 1820, 30, 47, 24)

Figure 2
figure 2

Reduced median network relating U8a HVSI/II sequences. Star is CRS for HVSI and 73, 282 for HVSII. Numbers along links refer to nucleotide positions minus 16000; Underlined subjects were complete sequenced. The broken line is the less probable link in accordance with completed sequences (Fig.1). Codes are: AUS = Austrian; BQ = Basque; BRI = British-australian; EST = Estonian; FIN = Finn; FRA = French; GER = German; KAR = Karelian; LIT = Lithuanian; NCE = North-central European; NEE = Northeast European; NOR = Norwegian; NSP = Northeast Spanish; POL = Polish; POR = Portuguese; RUS = Russian; SCT = Scotch; SIC = Sicilian; SLO = Slovenian; SPA = Spanish; SSP = South Spanish; SWZ = Swiss; TUK = Turk; VOL = Volga-Ural.

Table 1 Gene diversity (p) and frequency of U8a/1 in different European macro-regions.

Conclusion

In summary, the analysis of U8a lineages supports the idea that Basques have lived in their country since the Paleolithic, and that they could have participated in demographic re-expansions to repopulate central Europe in the last interglacial periods. Furthermore, these primitive U8a founders most probably reached the Basque area from the East through Europe and not through North Africa. However, the fact that we can trace some Basque lineages back to the Paleolithic does not support the generalized supposition that the present day Basque population is the best representative of Paleolithic Europeans. First of all, U8a haplotypes only represent 1% of the present day Basque maternal pool, therefore, a complex set of different mtDNA lineages with possible different histories are left unstudied. In addition, there is empiric evidence that Basques have received recent male gene flow from adjacent areas [31], and even possible maternal North African influences predating the Muslim Iberian invasion [32]. Furthermore, ancient DNA studies on Basque historic and prehistoric samples [33] have detected important mtDNA haplogroup frequency fluctuations along different periods. Definitively, like other European populations, Basques have also suffered migration and genetic drift effects throughout its long history.

Methods

Samples

DNA isolated from bucal swabs or blood samples from 211 autochthonous, unrelated Basques from the Iberian provinces were analyzed. Appropriate informed consent to anonymously use their data was obtained from all the individuals sampled.

HVSI-II and RFLPs

Total DNA was PCR amplified as in Pinto et al [34], and directly sequenced for both complementary strands as detailed in Rando et al [35]. A sequence of 978 bp of the HVSI-II of the mtDNA control region, from position 15997 to 00408 [36] was determined and sorted into defined haplogroups [24]. To confirm this HVS-based haplogroup classification, all individuals assigned to a specific haplogroup were additionally tested by restriction analysis of the diagnostic coding region mutations proposed to unambiguously classify sequences into haplogroups [24].

Complete mtDNA sequences

Four Basques (three U8a, one K1) and one Jordan (U8b) rare lineages belonging to the U/K haplogroup were fully sequenced. The complete mitochondrial DNAs (mtDNA) were amplified by PCR using primer pairs already described [19]. Amplified products were sequenced for both complementary strands with the Big Dye Terminator Cycle sequencing kit (Applied Biosystems). Sequencing reactions were analyzed on an Applied Biosystems 3100 DNA analyzer.

Genetic analyses

Genetic diversity (p ± σ) was estimated as the average number of nucleotide differences between two sequences [37], using the HVSI region in the range from 16070 to 16365 nucleotide positions.

Phylogenetic relationships among complete mtDNA sequences, and among control region mtDNA sequences, were established using the reduced median network algorithm [11]. In addition to our five sequences, seven lineages were added: Dutch (Dut86; GenBank: DQ112821) [[38] and P. Shen, pers. comm], Spanish (Spa, GenBank: AY882392) and Italian (Ita, GenBank: AY882393) [20], 3 Finns (Fin; GenBank: AY339551, AY339552 and AY339553) [14] and British-Australian (BriT7; [T. Obayashi, M. Tanaka, pers. comm.]).

The presence of U8a was tested in published sequences by the 16146 16342 HVSI motif and/or the 73 282 HVSII motif, and U8b by the 16189 16234 motif. 19,133 Eurasian and 1,430 North-African published and unpublished HVSI/HVSII sequences were analyzed.

The average sequence divergence [16] for complete sequences was converted in time applying a mutation rate of 1.26 × 10-8 [15].

Time estimations, based on control region mtDNA, were calculated as the mean divergence ρ [39] from inferred ancestral sequence types and converted into time by assuming that one transition within np 16090–16365 corresponds to 20,180 years [29].

Accesion numbers

The five complete mitochondrial DNA sequences are registered in [GenBank: DQ200801, DQ200802, DQ200803, DQ200804, and DQ200805].

References

  1. Barandiaran I: El paleolítico y el mesolítico. Prehistoria de la Península Ibérica. Edited by: Barandiaran I, Marti B, Del Rincón MA, Maya JL. 1998, Barcelona, Spain: Ariel, 1-120.

    Google Scholar 

  2. Cavalli-Sforza LL, Menozzi P, Piazza A: The history and geography of human genes. 1994, New Jersey: Princeton University Press

    Google Scholar 

  3. Bertranpetit J, Sala J, Calafell F, Underhill PA, Moral P, Comas D: Human mitochondrial DNA variation and the origin of Basques. Ann Hum Genet. 1995, 59: 63-81.

    Article  PubMed  CAS  Google Scholar 

  4. Achilli A, Rengo C, Magri C, Battaglia V, Olivieri A, Scozzari R, Cruciani F, Zeviani M, Briem E, Carelli V, Moral P, Dugoujon JM, Roostalu U, Loogväli EL, Kivisild T, Bandelt HJ, Richards M, Villems R, Santachiara-Benerecetti AS, Semino O, Torroni A: The molecular dissection of mtDNA haplogroup H confirms that the Franco-Cantabrian glacial refuge was a major source for the European gene pool. Am J Hum Genet. 2004, 75: 910-918. 10.1086/425590.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  5. Pereira L, Richards M, Goios A, Alonso A, Albarrán C, García O, Behar DM, Gölge M, Hatina J, Al-Gazali L, Bradley D, Macaulay V, Amorim A: High-resolution mtDNA evidence for the late-glacial resettlement of Europe from an Iberian refugium. Genome Res. 2005, 15: 19-24. 10.1101/gr.3182305.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  6. Torroni A, Bandelt H-J, D'Urbano L, Lahermo P, Moral P, Sellitto D, Rengo C, Forster P, Savontaus M-L, Bonné-Tamir B, Scozzari R: mtDNA analysis reveals a major late paleolithic population expansion from southwestern to northeastern Europe. Am J Hum Genet. 1998, 62: 1137-1152. 10.1086/301822.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  7. Torroni A, Bandelt H-J, Macaulay V, Richards M, Cruciani F, Rengo C, Martinez-Cabrera V, Villems R, Kivisild T, Metspalu E, Parik J, Tolk HV, Tambets K, Forster P, Karger B, Francalacci P, Rudan P, Janicijevic B, Rickards O, Savontaus ML, Huoponen K, Laitinen V, Koivumaki S, Sykes B, Hickey E, Novelletto A, Moral P, Sellitto D, Coppa A, Al-Zaheri N, Santachiara-Benerecetti AS, Semino O, Scozzari R: A signal, from human mtDNA, of postglacial recolonization in Europe. Am J Hum Genet. 2001, 69: 844-852. 10.1086/323485.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  8. Loogväli E-L, Roostalu U, Malyarchuk BA, Derenko MV, Kivisild T, Metspalu E, Tambets K, Reidla M, Tolk H-V, Parik J, Pennarun E, Laos S, Lunkina A, Golubenko M, Barac L, Pericic M, Balanovsky PO, Gusar V, Khusnutdinova EK, Stepanov V, Puzyrev V, Rudan P, Balanovska EV, Grechanina E, Richard C, Moisan J-P, Chaventré A, Anagnou NP, Pappa KI, Michalodimitrakis EN, Claustres M, Gölge M, Mikerezi I, Usanga E, Villems R: Disuniting uniformity: A pied cladistic canvas of mtDNA haplogroup H in Eurasia. Mol Biol Evol. 2004, 21: 2012-2021. 10.1093/molbev/msh209.

    Article  PubMed  Google Scholar 

  9. Harpending HC, Eller E: Human diversity and its history. The Biology of biodiversity. Edited by: Kato M, Takahata N. 1999, Tokyo: Springer-Verlag, 301-314.

    Google Scholar 

  10. Richards M, Macaulay V: The mitochondrial gene tree comes of age. Am J Hum Genet. 2001, 68: 1315-1320. 10.1086/320615.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  11. Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.

    Article  PubMed  CAS  Google Scholar 

  12. Palanichamy MG, 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-978. 10.1086/425871.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  13. 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.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  14. Finnilä S, Lehtonen MS, Majamaa K: Phylogenetic network for European mtDNA. Am J Hum Genet. 2001, 68: 1475-1484. 10.1086/320591.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 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. Proc Natl Acad Sci USA. 2003, 100: 171-176. 10.1073/pnas.0136972100.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  16. Saillard J, Forster P, Lynnerup N, Bandelt H-J, Nurby S: MtDNA variation among Greenland Eskimos: the edge of the Beringian expansion. Am J Hum Genet. 2000, 67: 718-726. 10.1086/303038.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  17. 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.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  18. Harpending HC, Rogers AR: Genetic perspectives on human origins and differentiation. Ann Rev Genomics Hum Genet. 2000, 1: 361-385. 10.1146/annurev.genom.1.1.361.

    Article  CAS  Google Scholar 

  19. 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.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  20. Achilli A, Rengo C, Battaglia V, Pala M, Olivieri A, Fornarino S, Magri C, Scozzari R, Babudri N, Santachiara-Benerecetti AS, Bandelt HJ, Semino O, Torroni A: Saami and Berbers–An unexpected mitochondrial DNA link. Am J Hum Genet. 2005, 76: 883-886. 10.1086/430073.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  21. Tanaka M, Cabrera VM, González AM, Larruga JM, Takeyasu T, Fuku N, Guo L-J, Hirose R, Fujita Y, Kurata M, Shinoda K-i, 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.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  22. Dubut V, Chollet L, Murail P, Cartault F, Beraud-Colomb E, Serre M, Mogentale-Profizi N: MtDNA polymorphisms in five French groups: importance of regional sampling. Eur J Hum Genet. 2003, 12: 293-300. 10.1038/sj.ejhg.5201145.

    Article  Google Scholar 

  23. Bermisheva MA, Tambets K, Villems R, Khusnutdinova EK: Diversity of mitochondrial DNA haplogroups in ethnic populations of the Volga-Ural region. Mol Biol. 2002, 36: 802-812. 10.1023/A:1021677708482.

    Article  CAS  Google Scholar 

  24. Richards M, Macaulay V, Hickey E, Vega E, Sykes B, Guida V, Rengo C, Rengo C, Sellito 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, Novelleto A, Oppenheim A, Norby S, Al-Zaheri N, Santachiara-Benerecetti S, Scozzari R, Torroni A, Bandelt H-J: Tracing European founder lineages in the Near Eastern mtDNA pool. Am J Hum Genet. 2000, 67: 1251-1276.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  25. Thomas MG, Weale ME, Jones AL, Richards M, Skorecki K, Torroni A, Scozzari R, Gratrix F, Tarekegn A, Wilson JF, Capelli C, Bradman N, Goldstein DB: Founding mothers of Jewish communities: geographically separated Jewish groups were independently founded by very few female ancestors. Am J Hum Genet. 2002, 70: 1411-1420. 10.1086/340609.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  26. Cherni L, Loueslati BY, Pereira L, Ennafaâ H, Amorim A, El Gaaied ABA: Female gene pools of Berber and Arab neighboring communities in Central Tunisia: Microstructure of mtDNA variation in North Africa. Hum Biol. 2005, 77: 61-70. 10.1353/hub.2005.0028.

    Article  PubMed  Google Scholar 

  27. Metspalu M, Kivisild T, Metspalu E, Parik J, Hudjashov G, Kaldma K, Serk P, Karmín M, Behar DM, Gilbert MT, 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 human. BMC Genet. 2004, 5: 26-10.1186/1471-2156-5-26.

    Article  PubMed  PubMed Central  Google Scholar 

  28. Bandelt H-J, Forster P, Sykes BC, Richards MB: Mitochondrial portraits of human populations using median networks. Genetics. 1995, 141: 743-753.

    PubMed  CAS  PubMed Central  Google Scholar 

  29. Forster P, Harding R, Torroni A, Bandelt H-J: Origin and evolution of Native American mtDNA variation: a reappraisal. Am J Hum Genet. 1996, 59: 935-945.

    PubMed  CAS  PubMed Central  Google Scholar 

  30. Gamble C: The Paleolithic societies of Europe. 1999, Cambridge: Cambridge University Press

    Google Scholar 

  31. Hurles ME, Veitia R, Arroyo E, Armenteros M, Bertranpetit J, Pérez-Lezaun A, Bosch E, Shlumukova M, Cambon-Thamsen A, McElreavey K, López de Munain A, Rohl A, Wilson LJ, Singh L, Pandya A, Santos FR, Tyler-Smith C, Jobling MA: Recent male-mediated gene flow over a linguistic barrier in Iberia, suggested by analysis of a Y-chromosomal DNA polymorphism. Am J Hum Genet. 1999, 65: 1437-1448. 10.1086/302617.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  32. Alzualde A, Izaguirre N, Alonso S, Alonso A, Albarrán C, Azkarate A, De la Rúa C: Insights into the "isolation" of the Basques: mtDNA lineages from the historical site of Aldaieta (6)th–7th c. AD). Am J Phys Anthropol in press. 2005

    Google Scholar 

  33. Alzualde A, Izagirre N, Alonso S, Alonso A, de la Rúa C: Temporal mitochondrial DNA variation in the Basque country: Influence of post-Neolithic events. Ann Hum Genet. 2005, 69 (Pt 6): 665-79. 10.1046/j.1529-8817.2005.00170.x.

    Article  PubMed  CAS  Google Scholar 

  34. Pinto F, González AM, Hernández M, Larruga JM, Cabrera VM: Genetic relationship between the Canary Islanders and their African and Spanish ancestors inferred from mitochondrial DNA sequences. Ann Hum Genet. 1996, 60: 321-330.

    Article  PubMed  CAS  Google Scholar 

  35. Rando JC, Pinto F, González AM, Hernández M, Larruga JM, Cabrera VM, 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-550. 10.1046/j.1469-1809.1998.6260531.x.

    Article  PubMed  CAS  Google Scholar 

  36. 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.

    Article  PubMed  CAS  Google Scholar 

  37. Nei M: Molecular evolutionary genetics. 1987, New York: Columbia University Press

    Google Scholar 

  38. Kivisild T, Shen P, Wall DP, Do B, Sung R, Davis KK, Passarino G, Underhill PA, Scharfe C, Torroni A, Scozzari R, Modiano D, Coppa A, de Knjiff P, Feldman MW, Cavalli-Sforza LL, Oefner PJ: The role of selection in the evolution of human mitochondrial. Genetics. 2006, 172: 373-387. 10.1534/genetics.105.043901.

    Article  PubMed  CAS  PubMed Central  Google Scholar 

  39. Morral N, Bertranpetit J, Estivill X, Nunes V, Casals T, Giménez J, Reis A, Varon-Mateeva R, Macek M, Kalaydjieva L, Angelicheva D, Dancheva R, Romeo G, Russo MP, Garnerone S, Restagno G, Ferrari M, Magnani C, Claustres M, Desgeoges M, Schwartz M, Dallapiccola B, Novelli G, Ferec C, de Arce M, Nemeti M, Kere J, Anvret M, Dahl N, Kadasi L: The origin of the major cystic fibrosis mutation (Δ F508) in European populations. Nat Genet. 1994, 7: 169-175. 10.1038/ng0694-169.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgements

We thank Dr. P. Shen and Drs. T. Obayashi and M. Tanaka, for providing us information about the Dutch (Dut86) and British-Australian (BriT7) sequences, respectively.

This study was supported by grants BMC2001-3511 from Ministerio de Ciencia y Tecnología and COF2002-015 from Gobierno de Canarias to V.M.C

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Ana M González.

Additional information

Authors' contributions

We have distinguished the following principal steps in the work:

1) Conception and design

2) Collection of samples

3) Sequencing of the HVSI-II regions and RFLPs analysis

4) Complete mtDNA sequencing of rare lineages

5) Collection of data to assess the distribution of U8a/b haplogroups in Europe, Near-East, Asia and Africa

6) Phylogenetic reconstruction for HVSI, and complete mtDNA sequences

7) Discussion of the results obtained

8) To draft the manuscript (text, table, figures and additional files)

AMG participated in 1, 4, 5, 6, 7 and 8 steps. OG participated in 1,2, 3, 5, and 7 steps. JML participated in 1, 3, 4, 5, 7 and 8 steps. VMC participated in 1, 2, 4, 6, 7 and 8 steps. All authors read and approved the final manuscript.

Electronic supplementary material

12864_2005_507_MOESM1_ESM.xls

Additional File 1: European distribution of subhaplogroup U8a/1. File González_add1.xls have data about sample and percentage of U8a/1 haplogroup in different populations from Europe and North-Africa. (XLS 55 KB)

12864_2005_507_MOESM2_ESM.doc

Additional File 2: References cited in additional file 1. File González_add2.doc shows a list of references cited in additional file 1. (DOC 62 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Authors’ original file for figure 3

Rights and permissions

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

Reprints and permissions

About this article

Cite this article

González, A.M., García, O., Larruga, J.M. et al. The mitochondrial lineage U8a reveals a Paleolithic settlement in the Basque country. BMC Genomics 7, 124 (2006). https://doi.org/10.1186/1471-2164-7-124

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2164-7-124

Keywords