Sami Lapps PDF
Sami Lapps PDF
Sami Lapps PDF
74:661–682, 2004
The Saami are regarded as extreme genetic outliers among European populations. In this study, a high-resolution
phylogenetic analysis of Saami genetic heritage was undertaken in a comprehensive context, through use of ma-
ternally inherited mitochondrial DNA (mtDNA) and paternally inherited Y-chromosomal variation. DNA variants
present in the Saami were compared with those found in Europe and Siberia, through use of both new and previously
published data from 445 Saami and 17,096 western Eurasian and Siberian mtDNA samples, as well as 127 Saami
and 2,840 western Eurasian and Siberian Y-chromosome samples. It was shown that the “Saami motif” variant of
mtDNA haplogroup U5b is present in a large area outside Scandinavia. A detailed phylogeographic analysis of
one of the predominant Saami mtDNA haplogroups, U5b1b, which also includes the lineages of the “Saami motif,”
was undertaken in 31 populations. The results indicate that the origin of U5b1b, as for the other predominant
Saami haplogroup, V, is most likely in western, rather than eastern, Europe. Furthermore, an additional haplogroup
(H1) spread among the Saami was virtually absent in 781 Samoyed and Ob-Ugric Siberians but was present in
western and central European populations. The Y-chromosomal variety in the Saami is also consistent with their
European ancestry. It suggests that the large genetic separation of the Saami from other Europeans is best explained
by assuming that the Saami are descendants of a narrow, distinctive subset of Europeans. In particular, no evidence
of a significant directional gene flow from extant aboriginal Siberian populations into the haploid gene pools of
the Saami was found.
Estonians. At present, ∼60,000 Saami live in the north- and the variations in the mtDNA coding region of Eu-
ern regions of Norway, Sweden, and Finland, as well as ropean populations have been exploited to generate a
in the Russian Kola Peninsula (Haetta 1996). Soon after common phylogenetic nomenclature for mtDNA vari-
the beginning of the retreat of the ice sheets covering ants in Europe (Torroni et al. 1996; Macaulay et al.
the area, in the 8th–10th millennia before present (BP), 1999). Classification of the Saami mtDNA lineages re-
populations of hunters and fishermen, the producers of vealed that the absolute majority of these are clustered
the Mesolithic Komsa and Fosna-Hensbacka cultures, in a subset of the European mtDNA pool (Torroni et
inhabited the coastal region of Scandinavia, extending al. 1998; Villems et al. 1998), where two haplogroups—
well into Finland and to the Kola Peninsula (Kozlowski V and U—cover ∼80% of the variability of this subset
and Bandi 1984; Nygaard 1989; Sumkin 1990). The (Torroni et al. 1998; Tambets et al. 2000).
linguistic affiliation of these pioneer settlers of the north In contrast to the predominance of European mtDNA
is largely unknown, but it has been suggested that they haplogroups observed among the Saami, nearly half of
are the descendants of the Ahrensburgian population, their Y chromosomes share a TatC allele (haplogroup
which migrated toward the north from western Europe, N3, according to the nomenclature of the Y Chromosome
along the Atlantic coast of Norway. It has been proposed Consortium [YCC 2002]) with most Finno-Ugric and
that they might have been the ancestors of the present- Siberian populations. This variant is found at high fre-
day Saami (e.g., Sumkin 1990). Another presumably im- quencies among Siberian populations, such as the Yakuts
portant component in the postglacial recolonization of and the Buryats, but is virtually absent in western and
northern Fennoscandia came from the east, via Karelia Mediterranean Europe; even among the Norwegians and
and Finland. It has been associated with the movement the Swedes, populations that have historically lived in
of Mesolithic populations, carriers of post-Swiderian close proximity to the Saami, it is found at frequencies
cultures, to the north. Starting in the Neolithic period, of only 4%–8% (Zerjal et al. 2001; Passarino et al.
the northern population came into contact with tribes 2002). High frequencies of the TatC allele have also been
of territories lying to the south (e.g., Sumkin 1990). observed in Baltic (30%–40%) and Volga-Finnic–speak-
Thus, according to archeological data, the present-day ing populations (20%–50%) (Zerjal et al. 1997; Rootsi
Saami population might have been shaped in different et al. 2000; Rosser et al. 2000; Semino et al. 2000; Lai-
times both by the eastern and western influences.
tinen et al. 2002). These findings have been interpreted
Analyses of classic chromosomal marker variation
according to the classic view that a substantial element
have demonstrated that the genetic distances between
of the Saami (and other European Finno-Ugric–speaking
the Saami and other European populations are signifi-
populations) genetic lineages originated in a recent mi-
cantly larger than between any other pair of European
gration from Asia (Zerjal et al. 1997, 2001).
populations (Cavalli-Sforza et al. 1994). This distinctive
Another half of the Saami paternal lineages have pri-
pattern has suggested to some that the genetic com-
marily mutations M170, SRY-1532, or M173 (Semino
position of the Saami population arose from extensive
et al. 2000; Raitio et al. 2001; Wells et al. 2001). Ac-
admixture between Caucasoid and Mongoloid popu-
cordingly, they are identified as haplogroups I, R1a, and
lations. Estimates of the relative contribution of each
have varied from an equal input to ∼80% Caucasoid R1b (YCC 2002). Y chromosomes possessing those mu-
and 20% Mongoloid genes (Guglielmino et al. 1990; tations are widely spread in European populations (Ros-
Cavalli-Sforza et al. 1994, pp 272–273). However, not ser et al. 2000; Semino et al. 2000). Whereas SRY-1532
all classic genetic studies support an idea of an extensive and M173 are present at moderate frequencies in some
Caucasoid-Mongoloid admixture in the Saami genetic Siberian populations as well, M170 Y chromosomes are
background (e.g., Beckman et al. 1988, 1993). very rare there (Wells et al. 2001; Karafet et al. 2002).
Similarly, studies of mtDNA have identified large ge- Here, we first analyze mtDNA variation in the Swed-
netic distances between the Saami and other Europeans, ish Saami and improve the resolution level of their Y-
including the Finns (Sajantila and Pääbo 1995; Sajantila chromosomal haplotypes. Second, we examine the hy-
et al. 1995). Likewise, Lahermo et al. (1996) found no pothetical descent of the Saami gene pool from Siberian
overlap between Saami and the remaining European ancestors. We exploit extensive data on mtDNA and Y-
mtDNA patterns and concluded that the Saami and the chromosome variation in Eurasia, in particular among
Finns must have different genetic histories. One alter- Uralic-, Indo-European–, and Altaic-speaking popula-
native hypothesis to explain the presence of genetic dif- tions in eastern Europe and western Siberia. Finally, we
ferences and language similarities in the Finns and the refine the phylogenetic topology of a frequent Saami
Saami involves a language shift by the Finns from Indo- mtDNA haplogroup, U5b, to identify the founder hap-
European to Finno-Ugric (Sajantila and Pääbo 1995). lotypes. We draw conclusions about the possible origin
The combined sequence information from the first of mtDNA and NRY haplotypes present in the Saami
hypervariable segment (HVS-I) in the D-loop of mtDNA gene pool.
Tambets et al.: mtDNA and NRY Variation of the Saami 663
Subjects and Methods Karafet et al. 2002). The YAP and 12f2 polymorphisms
were identified following the procedures of Hammer and
Population Samples Horai (1995) and Casanova et al. (1985), respectively.
Furthermore, 35 Swedish Saami DNA samples that have
We analyzed mtDNA sequences from 17,541 Eurasian previously been analyzed for Y-chromosomal variation
individuals sampled from 65 different Eurasian popu- by Rosser et al. (2000) were also typed in the pres-
lations (table 1). These data consist of previously pub- ent study for the biallelic markers M89, M52, M130
lished sequences (n p 9,154) and new sequences an- (RSP4Y), M170, M173, M178, M201, and M269 (YCC
alyzed in this study (n p 8,387). Comparisons of dif- 2002).
ferent Saami subpopulations are presented separately in
table 2. Y-chromosomal markers were also analyzed Sequencing
from 32 different populations. Data from 1,598 of a
total of 2,967 individuals reported in table 3 were taken PCR-amplified products were purified using shrimp al-
from published sources, whereas the remaining 1,369 kaline phosphatase and exonuclease treatment follow-
DNA samples were extracted and analyzed in the present ing Kaessmann (1999). These were sequenced using the
study. Further details about the populations other than DYEnamic ET terminator cycle sequencing kit (Amer-
the Saami will be published elsewhere. Blood samples sham Pharmacia Biotech) on an ABI 377 DNA Sequencer.
were collected from healthy unrelated individuals after Sequences were aligned and analyzed using the Genetics
obtaining informed consent. DNA was extracted using Computer Group Wisconsin Package.
the phenol-chloroform method, as used by Sambrook
Data Analysis
(1989).
Principal-component (PC) analysis was undertaken
mtDNA Analysis using the POPSTR program, kindly provided by H. Har-
pending (see details described by Richards et al. [2002]).
HVS-I was sequenced between nucleotide positions Only those haplogroups that had a noticeable impact
(nps) 16024 and 16383 of the revised Cambridge Ref- on the scatterplot were used for the analysis (the sum
erence Sequence (Andrews et al. 1999) in 73 Swedish of the absolute values for both coordinates of each allele
Saami and 8,314 other mtDNA samples. RFLP analysis was 10.2 [M, V, U5, H, A, J, and U4 for mtDNA and
of diagnostic mtDNA positions was performed, and Q, N2, I, R1b, R1a, and N3 for the Y chromosome]).
mtDNA haplogroups were assigned to each sample by Bayesian 95% credible regions (CRs) for haplogroup
use of published criteria (Torroni et al. 1992, 1994, frequencies were calculated with the computer program
1996; Richards et al. 1998; Macaulay et al. 1999; Schurr SAMPLING, kindly provided by V. Macaulay. The di-
et al. 1999; Finnilä et al. 2000, 2001; Herrnstad et al. versity of mtDNA haplotypes was estimated as by Nei
2002; Kivisild et al. 2002). In addition, all U5b genomes (1987).
(n p 160) from 31 populations were typed by sequenc- Phylogenetic networks of mtDNA HVS-I haplotypes
ing for the np 5656 ArG mutation that differentiates were constructed using the program Network 3.1.1.1
between individuals belonging to haplogroups U5b* and (Fluxus Engineering Web site). Different weights were
U5b1. All identified U5b1 individuals (5656G) were fur- assigned to substitutions, as in the study by Richards et
ther sequenced to assess the variation at nps 7385 (A or al. (1998). Coalescence-age calculations and SDs were
G) and 10927 (T or C). A summary of these data is estimated following the formulae of Forster et al. (1996)
given in table 4. and Saillard et al. (2000).
Table 2
Frequencies (%) of mtDNA and Y-Chromosomal Haplogroups in Different Saami
Subpopulations
FREQUENCY (95% CR)a
(%)
Swedish Saamib Finnish Saamic Norwegian Saamid Kola Saamie
HAPLOGROUP (n p 98/35) (n p 69/69) (n p 278/ …) (n p … /23)
mtDNA:
H: 3.1 (1.1–8.6) 2.9 (.9–9.9) 4.7 (2.8–7.8) …
H1 .0 (.0–3.0) .0 (.0–4.2) 4.0 (2.2–6.9) …
W 1.0 (.2–5.5) .0 (.0–4.2) 1.4 (.6–3.6) …
T .0 (.0–3.0) .0 (.0–4.2) .4 (.1–2.0) …
U: 26.5 (18.8–36.1) 43.5 (32.4–55.3) 57.6 (51.7–63.2) …
U5a .0 (.0–3.0) 2.9 (.9–9.9) .7 (.2–2.6) …
U5b 26.5 (18.8–36.1) 40.6 (29.8–52.4) 56.8 (51.0–62.5) …
V 68.4 (58.6–76.7) 37.7 (27.2–49.5) 33.1 (27.8–38.8) …
M: 1.0 (.2–5.5) 15.9 (9.2–26.4) 2.9 (1.5–5.6) …
D5 .0 (.0–3.0) 8.7 (4.1–17.7) 2.9 (1.5–5.6) …
Z 1.0 (.2–5.5) 7.2 (3.2–15.9) .0 (.0–1.1) …
NRY:
E .0 (.0–8.0) .0 (.0–4.2) … 8.7 (2.7–27.0)
F* 5.7 (1.8–18.7) … .0 (.0–11.7)
I 31.4 (18.6–48.1) 40.6 (29.8–52.4)f … 17.4 (7.1–37.4)
J .0 (.0–.8) … 4.3 (1.0–21.1)
N3 37.1 (23.1–53.8) 55.1 (43.3–66.3) … 39.1 (22.1–59.4)
R1b 5.7 (1.8–18.7) 1.4 (.3–7.7) … 8.7 (2.7–27.0)
R1a 20.0 (10.1–36.0) 2.9 (.9–9.9) … 21.7 (9.8–42.2)
a
Sample size (n) is given separately for mtDNA/NRY data sets, respectively.
b
mtDNA data are from this study and inferred from Sajantila et al. (1995); NRY data from
Rosser et al. (2000) are further analyzed for markers M89, M52, M130 (RSP4Y), M170, M173,
M178, M201, and M269 (YCC 2002) in the present study.
c
mtDNA data are inferred from Sajantila et al. (1995); NRY data are inferred from Raitio et
al. (2001).
d
mtDNA data are inferred from Sajantila et al. (1995), Dupuy and Olaisen (1996), and Del-
ghandi et al. (1998).
e
NRY data are from Wells et al. (2001).
f
Frequency is that of haplogroups F, I, and J combined.
The phylogenetic network in figure 1 relates all the The frequencies of the most widely spread haplo-
HVS-I sequence haplotypes in the three geographic sub- groups among the Saami in other analyzed populations
populations of the Saami. It also indicates a limited num- are shown in table 1. Both haplogroups V and U5b are
ber of highly frequent haplotypes, which are, in most spread at moderate frequencies across Europe, from Ibe-
cases, shared by all analyzed geographic subpopulations ria to the Ural Mountains. In contrast, among 393 Ob-
(see table 2). Thus, 98% of the Saami mtDNA pool is Ugric speakers (including 98 Mansi, published elsewhere
encompassed within haplogroups V, U5b, H, Z, and D5. by Derbeneva et al. [2002b]) and 388 Samoyeds (in-
The two largest subsets, V and U5b, account for a lion’s cluding 58 Nenets, published elsewhere by Saillard et al.
share (89%) of the Saami mtDNA pool. The represen- [2000]), only one haplogroup V–carrying individual was
tative proportion of individual haplogroups varies sub- found in the Mansi sample, whereas not a single Saami
stantially among different Saami subpopulations (table variant of U5b (U5b1b1) was identified there (table 1).
2). Haplogroup V is by far the most frequent haplo- These two variants of maternal lineages are virtually
group in the Swedish Saami and is present at significantly absent in the other 1,199 native Siberians analyzed in
lower frequencies (P ! .0001) in Norwegian and Finnish this study.
subpopulations. We note that, in all three Saami sub- The third-most-frequent haplogroup among the Saami
populations, U5b was virtually the only subclade of the is H, which is present at a frequency of 4%, ∼10-fold
otherwise frequent and divergent western Eurasian hap- lower than that of other North European populations
logroup U5 (Tambets et al. 2003). It is the most frequent (Richards et al. 2000). The majority of the Saami hap-
haplogroup among Norwegian and Finnish subpopu- logroup H lineages (61%) contain transitions at nps 73
lations (table 2). and 16162. According to information obtained from
Table 3
Frequencies (%) of Selected Y-Chromosomal Haplogroups in the Saami and in Other Eurasian Populations
FREQUENCY OFHAPLOGROUP
(%)
REGION AND LANGUAGE SAMPLE
POPULATIONa FAMILYb SIZE Ic R1bd R1a N3 N2 Q C REFERENCE(S)
Scandinavia:
*Saamie U-FU 127 (58)f 25.9f 3.9 11.0 47.2 0 0 0 Raitio et al. 2001; Wells et al. 2001; present study
Swedes IE-G 141 48.2 22.0 18.4 2.8 0 0 0 Present study
Norwegians IE-G 72 40.3 27.8 23.6 6.9 0 0 0 Passarino et al. 2002
Danes IE-G 194 38.7 36.1 16.5 .5 0 0 0 Sanchez et al. 2003
Northern Europe:
*Finns U-FU 38 28.9 0 7.9 63.2 0 0 0 Zerjal et al. 2001
Estonians U-FU 209 18.2 9.1 33.5 30.6 0 0 0 Present study
Latvians IE-B 86 7.0 9.3 38.4 41.9 0 0 0 Present study
North Russians IE-S 77 20.8 0 40.3 28.6 0 0 0 Wells et al. 2001
Eastern Europe:
Russians IE-S 61 13.1 21.3 42.6 8.2 8.2 0 0 Karafet et al. 2002 (supplementary data obtained by request from the author)
Ukrainians IE-S 50 18.0 2.0 54.0 6.0 0 0 0 Semino et al. 2000
Polese IE-S 93 16.1 13.4 55.9 3.2 0 0 0 Present study
Volga-Ural region:
Maris U-FU 111 8.1 2.7 47.7 31.5 9.9 0 0 Semino et al. 2000; present study
Mordvin U-FU 83 19.3 13.3 26.5 16.9 2.4 0 0 Present study
Komis U-FU 94 5.3 16.0 33.0 22.3 12.8 0 0 Present study
Udmurts U-FU 87 1.1 2.3 10.3 56.3 28.7 0 0 Present study
Chuvashes A-T 79 11.4 3.8 31.6 17.7 10.1 0 1.3 Present study
Tatars A-T 126 4.0 8.7 34.1 18.3 4.8 .8 1.6 Present study
Southern Europe:
Croats IE-S 109 37.6 15.6 33.9 0 0 0 0 Barac et al. 2003
Western Europe:
French IE-R 61 16.4 59.0 0 0 0 0 0 Semino et al. 2000; present study
Central Europe:
Germans IE-G 16 37.5 0 0 0 0 0 0 Semino et al. 2000
Hungarians U-FU 113 28.3 20.4 20.4 .9 0 2.6 .9 Present study
Siberia:
Khants U-FU 47 0 19.1 4.3 38.3 38.3 0 0 Karafet et al. 2002 (supplementary data obtained by request from the author)
Nganasans U-SA 38 0 0 0 0 92.1 0 5.3 Karafet et al. 2002 (supplementary data obtained by request from the author)
Nenets U-SA 148 0 0 0 40.5 56.8 1.4 0 Karafet et al. 2002 (supplementary data obtained by request from the author)
Selkups U-SA 131 0 6.1 19.1 0 6.9 66.4 1.5 Karafet et al. 2002 (supplementary data obtained by request from the author)
Kets Isolate 48 0 0 0 0 0 93.7 6.2 Karafet et al. 2002 (supplementary data obtained by request from the author)
Dolgans A-T 67 1.5 1.5 16.4 22.4 11.9 0 37.3 Karafet et al. 2002 (supplementary data obtained by request from the author)
Yakuts A-T 155 1.3 1.9 1.9 88.4 0 0 3.2 Karafet et al. 2002 (supplementary data obtained by request from the author); present study
Buryats A-M 81 0 1.2 1.2 28.4 2.5 0 60.5 Karafet et al. 2002 (supplementary data obtained by request from the author)
Evenks A-TN 96 5.2 0 1.0 16.7 3.1 4.2 67.7 Karafet et al. 2002 (supplementary data obtained by request from the author)
Evens A-TN 31 3.2 0 6.5 12.9 0 0 74.2 Karafet et al. 2002 (supplementary data obtained by request from the author)
Altaians A-T 98 0 0 46.9 1.0 2.0 17.3 22.4 Karafet et al. 2002 (supplementary data obtained by request from the author)
Total 2,967
a
Populations for which data are deduced from published sources are indicated with an asterisk.
b
Language-family codes for each population are as follows: U p Uralic (FU p Finno-Ugric; SA p Samoyedic); IE p Indo-European (B p Baltic; G p Germanic; R p Romance; S p Slavic); A
p Altaic (M p Mongolic; T p Turkic, TN p Tungusic).
c
Haplogroup I is, in the case of data deduced from literature, defined by the characteristic STR pattern of haplogroup 2 Y chromosomes (Barac et al. 2003; present study).
d
Haplogroup R1b is, in the case of data deduced from literature, defined here by M173(xSRY-1532) as studied western Europeans (except three Danes, who belong to R1* [Sanchez et al. 2003])
from this clade have been shown to share the additional mutation M269 (Cruciani et al. 2002; present study).
e
35 Swedish Saami and 93 Poles from the study by Rosser et al. (2000) were further analyzed for markers M52, M89, M130, M170, M173, M178, M201, and M269.
f
In the estimation of haplogroup I proportions, only the data from the present study and from Wells et al. (2001) were used. The corresponding sample sizes are given in parentheses.
Table 4
Phylogeographic Distribution of Haplogroup U5b Subclades (Excluding U5b1b1)
NO. OF INDIVIDUALS IN POPULATIONc
a
WSE SC EE Other BASE AT NP
1
5 7 0
U5b 6 3 9
HVS-I SUB 5 8 2
SEQUENCEa HGb al bo br cr cz fr gr hu it si sk sl sw ba es ko la mo po ru ta ud uk mc ng go tu sh ku td nd 6 5 7
270 U5b* 1 2 3 A ND ND
051-168-189-270-362 U5b* 1 1 A ND ND
086-192-270-304 U5b* 1 1 A ND ND
093-172-189-270-362 U5b* 1 1 A ND ND
093-234-270-325 U5b* 1 1 A ND ND
668
093-258-270-292-362 U5b* 1 1 A ND ND
145-189-270 U5b* 1 1 A ND ND
169-172-189-270 U5b* 1 1 A ND ND
169CA-192-235-270-304 U5b* 1 1 A ND ND
189-192-201-270-398 U5b* 1 1 A ND ND
189-192-217-234-270-398 U5b* 2 2 A ND ND
189-192-234-256-270-311-362-398 U5b* 1 1 A ND ND
189-192-266-270-398 U5b* 1 1 A ND ND
189-192-270 U5b* 1 1 1 1 1 5 A ND ND
189-192-270-398 U5b* 1 1 1 3 A ND ND
189-239CA-270 U5b* 1 1 A ND ND
189-261-270 U5b* 1 1 A ND ND
189-270 U5b* 1 3 2 6 A ND ND
189-270-288-292 U5b* 1 1 A ND ND
189-270-311 U5b* 1 1 1 3 A ND ND
189-292-311 U5b* 1 1 A ND ND
192-212-270 U5b* 1 1 A ND ND
192-224-261 U5b* 1 1 A ND ND
192-235-270-304 U5b* 1 1 A ND ND
192-257-270-304 U5b* 1 1 A ND ND
192-270 U5b* 1 1 1 3 A ND ND
192-270-278-304 U5b* 1 1 A ND ND
192-270-304 U5b* 1 1 2 A ND ND
192-311 U5b* 1 1 1 4 2 9 A ND ND
222-224-270-311 U5b* 1 1 A ND ND
224-270-311 U5b* 1 1 2 A ND ND
234-270 U5b* 1 1 A ND ND
258-270-292-362 U5b* 1 1 A ND ND
261-270 U5b* 1 1 1 3 A ND ND
261-270-390 U5b* 1 1 A ND ND
051-189-270 U5b1* 2 2 G A T
086-189-192-270-311-336 U5b1* 3 3 G A T
092-189-270-294 U5b1* 1 1 G A T
093-182-189-192-270 U5b1* 1 1 G A T
093-189-270 U5b1* 1 1 2 G A T
093-189-270-362 U5b1* 2 2 G A T
129-189-223-270 U5b1* 1 1 G A T
129-189-270 U5b1* 1 1 G A T
129-192-270 U5b1* 1 1 G A T
147-189-270 U5b1* 1 3 1 5 G A T
174-189-223-270 U5b1* 1 1 G A T
174-189-270-311 U5b1* 1 1 2 G A T
182-189-192-270 U5b1* 1 1 G A T
186-189-270-293 U5b1* 1 1 G A T
189-192-218-270-320 U5b1* 1 1 G A T
189-192-270 U5b1* 6 3 1 10 G A T
189-192-270-292 U5b1* 1 1 G A T
189-192-270-311 U5b1* 1 1 G A T
189-221-223-270 U5b1* 1 1 G A T
189-221-270 U5b1* 1 1 2 G A T
189-223-270 U5b1* 1 1 G A T
189-270 U5b1* 8 1 1 1 5 2 1 5 2 26 G A T
189-270-311-336 U5b1* 1 1 G A T
192-224-261-270 U5b1* 1 1 G A T
192-270 U5b1* 1 1 1 1 1 5 G A T
192-270-296 U5b1* 1 1 G A T
192-270-319 U5b1* 1 1 G A T
074-172-189-270 U5b1b 1 1 G G C
074-189-270 U5b1b 1 1 G G C
093-129-189-270 U5b1b 1 2 1 4 G G C
093-189-270 U5b1b 1 1 2 G G C
669
093-189-270-301 U5b1b 1 1 G G C
093-189-270-301-311 U5b1b 1 1 G G C
173-189-270-235 U5b1b 1 1 G G C
189-192-270 U5b1b 1 1 2 G G C
189-270 U5b1b 1 1 1 1 1 5 G G C
189-270-325 U5b1b 1 1 G G C
Total 3 5 1 20 4 13 3 2 1 9 1 1 9 4 14 11 6 5 1 16 5 1 13 3 1 2 2 1 1 2 160
NOTE.—mtDNAs belonging to subclade U5b1b1 were screened selectively for the presence of the mutations in nps 5656, 7385, and 10927. All of those (18 Saami, 5 individuals from eastern Europe, and 5 individuals from
western and southern Europe) possessed aforementioned transitions.
a
Mutations in sequences are relative to the revised Cambridge Reference Sequence (Andrews et al. 1999). HVS-I sequences are given minus 16,000; only tranversions are further indicated. ND p not determined.
b
SubHG p subhaplogroup.
c
Geographic subdivision codes are as follows: WSE p western and southern Europe; SC p Scandinavian Europe; EE p eastern Europe. Population codes are as follows: al p Albanians; bo p Bosnians; br p Bretons; cr p
Croats; cz p Czechs; fr p French; gr p Greeks; hu p Hungarians; it p Italians; si p Sicilians; sk p Slovaks; sl p Slovenians; sw p Swedes; ba p Bashkirs; es p Estonians; ko p Komis; la p Latvians; mo p Mordvin; po
p Poles; ru p Russians; ta p Tatarians; ud p Udmurts; uk p Ukrainians; mc p Moroccans; ng p Nogays; go p Georgians; tu p Turks; sh p Shors; ku p Kumyks; td p Tadjik. In addition to population data from table 1,
data of Bretons, Moroccans, Shors, Tadjiks, and Turks are shown (authors’ unpublished data).
d
Number of subjects analyzed in haplotype.
Figure 1 Phylogenetic network of 445 Saami mtDNA HVS-I sequences. Only those coding-region markers that have been analyzed in this
study are shown (for finer resolution of haplogroup U subbranches, see fig. 3). Observed mutations are numbered according to the revised Cambridge
Reference Sequence (Anderson et al. 1981; Andrews et al. 1999). The tree is rooted in haplogroup L3. Variable nps are shown on the links; nucleotide
change is specified by suffixes only for transversions. Nodes indicate different haplotypes and have sizes that are proportional to their frequencies.
HVS-I sequences are from Sajantila et al. (1995), using the corrections given by Bandelt et al. (2001) (n p 114 ); from Dupuy and Olaisen (1996)
(n p 197); from Delghandi et al. (1998) (n p 61 ); and from the present study (n p 73 , underlined samples). The analysis of coding region variation
has been performed for only the 73 samples analyzed in the present study; haplogroup affiliations for data published elsewhere are inferred from
HVS-I sequences. sf p Saami from Finland, sn p Saami from Norway, ss p Saami from Sweden. Note that the presence of three transversions
reported in a single haplotype (sf) from haplogroup V (Sajantila et al. 1995) has not been reconfirmed by an independent study.
Tambets et al.: mtDNA and NRY Variation of the Saami 671
complete sequences, these lineages belong to subhaplo- I, the classification of the samples was more problematic.
group H1, defined by the coding region transition at np However, the STR pattern of those haplogroups in com-
3010 (Finnilä et al. 2001). We notice here that HVS-I bination with the typed biallelic markers is, in most cases,
haplotypes with a transition at np 16162 in the Volga- haplogroup-specific enough to make the misclassification
Ural area people (Bermisheva et al. 2002) and other of the samples highly unlikely. The probability of mis-
European populations studied by us belong exclusively classification is higher in the data set from Raitio et al.
to this branch of subhaplogroup H1 (authors’ unpub- (2001), in which only some informative biallelic markers
lished data). This variant of mtDNA can be found both and no STRs were scored. Therefore, these samples are
in western and eastern Europe, being as frequent in Ger- left at a lower resolution level in table 2. Detailed criteria
manic-speaking Scandinavians and the Germans as used in the process of deduction are given in the footnotes
among the Norwegian Saami but absent in the Finnish of table 3.
and Swedish Saami (table 2). Though present in several Three major haplogroups—N3, I, and R1a—comprise
Volga-Finnic populations, H1 haplotypes with mutation 180% of the Saami Y-chromosomal gene pool (table 2)
16162 are extremely rare or absent in almost 2,000 Si- and appear, likewise, as the major haplogroups in north-
berian Ugric-, Samoyedic-, and Altaic-speaking people eastern Europe (table 3). The nearly equal distribution
(table 1) as well as in central Asians (Comas et al. 1998; of the Y-chromosomal variants among different Saami
Metspalu et al. 1999)—that is to say that this clade has populations strongly supports their common paternal
an overwhelmingly European phylogeographic pattern. history. The comparison of the Y-chromosomal hap-
Its presence only in the Norwegian Saami sample may logroup frequencies in different subpopulations shows
have been generated by admixture with the Norwegian that, similar to mitochondrial haplogroup frequencies,
population. This inference is supported by the presence some Y-chromosomal haplogroups (J and R1a) exhibit
of one H1 haplotype, with HVS-I motif 16162–16189, notable variation between geographic subpopulations of
found among the Saami, in the Norwegian sample (Hel- the Saami (table 2).
gason et al. 2001). This haplotype has not been found Haplogroup N3, the most frequent haplogroup in the
in any other population analyzed in the present study. Saami population, is distributed in eastern European and
Only a minor portion of the Saami maternal lineages northern Asian populations but it is rare or absent in
(average ∼5%) that exhibit restricted diversity belong western Europe (table 3). All analyzed Swedish Saami
to haplogroups that are characteristic of Asian pop- N3 lineages fall into subcluster N3a, defined by M178
ulations—that is, D5 and Z (table 1). These eastern (YCC 2002). Although N3a is widespread in Siberia,
Eurasian haplogroups are significantly more frequent other haplogroups, characteristic of Samoyedic-speaking
(P ! .05) among the Finnish Saami compared with Nor- and other Siberian populations (such as C and Q), are
wegian and Swedish Saami samples (table 2). Such fluc- either almost absent in Baltic-Finnic populations, in-
tuations in haplogroup frequencies could be due to ge- cluding the Saami, or are only sporadic, as for haplo-
netic drift or just due to stochastic variation in relatively group N2, which is found only among Volga region
small samples. At the same time, the similar pattern of Finnic speakers (table 3).
mtDNA haplogroup distribution found in different sub- The second largest haplogroup, I, corresponds to
populations of the Saami provides evidence of their com- roughly one-third of the Saami Y-chromosomal lineages.
mon genetic background. It is widespread in Europe but virtually absent in Asian
populations (table 3). The defining mutation of haplo-
Saami Y-Chromosomal Haplogroup Pattern group I, M170, most likely arose in Europe (Semino et
and Frequencies al. 2000). It is well dispersed over the continent, as well
as among Volga River basin Finnic and Turkic speakers
The Y-chromosomal variation of 35 Swedish Saami, (table 3). At the current level of resolution, the phylo-
published by Rosser et al. (2000), was further analyzed genetic reconstructions do not identify the geographic
to improve the resolution of paternal lineage clusters ac- origin of haplogroup I in Europe. Nevertheless, its vir-
cording to the YCC (2002). For comparison, data from tual absence among Samoyeds as well as among Ugric-
different European and Siberian populations were used speaking Mansis and Khants suggests that these popu-
(table 3). To include the data published before the high- lations have not shared a recent ancestry with the Saami.
resolution nomenclature of biallelic markers (YCC 2002) Indeed, the high frequencies of haplogroup I in the Nor-
was available, a combination of the results of the analysis wegians, Swedes, Finns, and the Saami (table 3) suggest
of biallelic markers and polymorphisms of STRs was used, that haplogroup I represents the heritage from the very
when possible (Zerjal et al. 2001), to infer the haplogroup first settlers of Fennoscandia. Analyses of microsatellite
frequencies presented in table 3. In some cases (haplo- variation of the Saami haplogroup I Y chromosomes
groups N3 and R1a), the updating of the data was simple reveals that 9 of 10 Saami M170 chromosomes share
and unambiguous. In the case of haplogroups R1b and the same modal microsatellite haplotype (6/9) or its de-
672 Am. J. Hum. Genet. 74:661–682, 2004
rivatives with other Nordic haplogroup I chromosomes Figure 2B shows a similar analysis undertaken on the
(authors’ unpublished data). The latter corresponds to Y-chromosomal haplogroup frequencies. Here, the first
haplogroup 2 chromosomes with alleles 14-23-10-11- PC is determined, on the one hand, by the western Eu-
13, defined by DYS19-390-391-392-393, respectively, in rasian–specific haplogroups I, R1a, and R1b and, on the
the studies by Helgason et al. (2000), Zerjal et al. (2001), other hand, by Asian-specific haplogroup C and by the
Weale et al. (2002), and Capelli et al. (2003). These largely Siberian-specific N2, which both separate all
haplotypes form a subclade of haplogroup I that is dif- Samoyeds and Ob-Ugric-speaking populations (except
ferent from the variety of equally abundant M170 chro- Selkups) from the Saami. In contrast to that of other
mosomes found in the Adriatic coastal region, where the Samoyedic-speaking populations, the Selkup Y-chromo-
Dinaric Modal Haplotype (16-24-11-11-13) is common somal pool is dominated by eastern Eurasian haplogroup
(Barac et al. 2003). Q, which explains their solitary location on the second
Haplogroup R1a, encompassing 11% of the Saami Y- dimension of the PC graph. Two Saami subpopulations
chromosomal gene pool, is a frequent Y-chromosomal map close together and are placed in the immediate vi-
variety in eastern Europe, but it is also found in Siberia, cinity of the Finns, Volga-Finnic, and Volga-Turkic pop-
the Altai Mountains (table 3), and India (Kivisild et al. ulations, as well as the Russians. The latter two are
2003). Like haplogroup I, it is absent among Nenets, characterized by haplogroup N3 at moderate frequencies
the largest Samoyed population, living partly in north- (table 3). Furthermore, at this resolution, Germanic-
east Europe. Haplogroup R1a is present in the Swedish speaking populations (including Scandinavian Norwe-
and Kola Saami at frequencies that are comparable to gians and Swedes) and the French form an additional
those observed among other Scandinavian and north- cluster that is close to other European populations in-
eastern European populations, but it is relatively rare in cluded in the analysis.
Finland, both in the Saami and non-Saami populations
(table 3). Phylogeography of Haplogroup U5b1b
Haplogroup R1b is present in the Saami population and the “Saami-Specific” Motif (U5b1b1)
at a level of 4%. It is the most frequent haplogroup in
western Europe, especially among the French, and it is Haplogroup U5b is found at low frequencies all over
relatively frequent in Scandinavian Germanic-speaking Europe (Richards et al. 1998). A distinct phylogenetic
populations, but it has not been found among the Finns. subbranch of U5b carries the “Saami-specific” HVS-I
In the Volga-Ural region populations, R1b is present at motif defined by 16144-16189-16270 in the study by
frequencies similar to those of the Saami (table 3). Sajantila et al. (1995) and has been named “U5b1” by
Haplogroups J and E, found solely among Kola Saami Richards et al. (1998). All U5b sequences share two
(table 2), may have arisen as a recent contribution from diagnostic synonymous substitutions at nps 7768 and
the neighboring northern Russian population, since these 14182 (see fig. 3A). The major U5b subset found among
Y-chromosomal variants are present there (Wells et al. the Finns contains a derived G allele at noncoding np
2001). 5656 and a synonymous substitution at np 12618 (Fin-
nilä et al. 2001). Through use of complete mitochondrial
PC Analysis sequence data, it has previously been demonstrated that
the Finnish sequences containing the “Saami-specific”
The results of the PC analysis based on mtDNA hap- HVS-I motif share an additional synonymous mutation
logroup frequencies (fig. 2A) illustrate why the Saami at np 10927, which is often associated with another
are identified as mitochondrial genetic “outliers” in Eu- transition at np 7385 (Finnilä et al. 2001).
rope. Here, it is the relative haplogroup proportions and To study the phylogeography of the “Saami-specific”
not the phylogeographic affiliations that are distinctive. branch of U5b, an analysis of np 5656 variation among
The first PC (fig. 2A) reflects predominantly the differ- the relevant genomes was performed. The results reveal
ences between frequencies of haplogroups characteristic (see table 4) that 95 of 160 samples carried the derived
of the eastern and western parts of Eurasia (H, J, V, and 5656G allele and belong, thus, to a subhaplogroup re-
U5 vs. A and M). All Saami subpopulations cluster to- ferred to as “U5b1.” We note that this subhaplogroup
gether as a distinct group, but the distance between the is different from the U5b1 branch defined by Richards
Saami and native Siberians is much greater than between et al. (1998) using HVS-I data alone. Here, “U5b1”
the Saami and other Europeans. Different distributions refers to a deeper phylogenetic node of U5b (see fig. 3A).
of the populations along the axis of the second PC arise The transition at np 7385 was always found to be as-
primarily from the different proportions of haplogroups sociated with a transition at np 10927 and is referred
H, U5, and V in the populations. In addition, the lack to as “U5b1b.” This cluster also includes the “Saami-
of haplogroups J, U4, and A among the Saami plays a specific” subclade, which we name here as “U5b1b1”
relatively important role here. (fig. 3). Notice that the sequence motif 16189-16192-
Tambets et al.: mtDNA and NRY Variation of the Saami 673
Figure 2 PC analysis based on mtDNA (A) and Y-chromosomal (B) haplogroup frequencies in some European and Siberian populations.
Only the seven haplogroups (squares) having the main impact on the scatter plot have been used (see the “Subjects and Methods” section).
Numbers in parenthesis indicate the proportion of the total genetic information retained by a given PC. The Saami population is subdivided
according to the present location of the subpopulations. Population data (triangles) are the same as in tables 1 and 3 and are listed here in
alphabetical order as follows: at p Altaians, ch p Chuvashes, ev p Evenks, fi p Finns, fr p French, ge p Germans, kh p Khants, km p
Khants and Mansis, ko p Komis, ne p Nenets, nn p Nganasans, no p Norwegians, sf p Saami from Finland, sk p Saami from Kola
Peninsula, Russia, sn p Saami from Norway, ss p Saami from Sweden, se p Selkups, sw p Swedes, vt p Volga-Uralic Turkic-speakers
(Bashkirs and Chuvashes). Inclusion of all populations listed in tables 1 and 3 did not change the overall outcome (data not shown).
16270 in U5b1 appears to be associated with at least three the exploration of the phylogeography of haplogroup
independent subhaplogroups of U5b: U5b* (5656A), U5b1b ambiguous, and these sequences were not in-
U5b1*/U5b1a (5656G, 7385A, and 10927T), and U5b1b cluded in the analysis.
(5656G, 7385G, and 10927C) (fig. 3; table 4). Therefore, The 5656G allele, including all major subsets of U5b1,
reliance on published HVS-I sequences alone, without is broadly distributed both in western and eastern Eu-
regard to the relevant coding region information, makes rope (fig. 3A; table 4). The U5b1b subclade is found all
Figure 3 A, Phylogenetic network of U5b1b lineages based on HVS-I sequences and its position in the phylogeny of haplogroup U.
Sequence information from Herrnstadt et al. (2002) and Finnilä et al. (2001) has been used for the coding region and HVS-II (see also table
4). The nucleotide positions relative to the revised Cambridge Reference Sequence (Anderson et al. 1981; Andrews et al. 1999), at which two
nodes differ, are listed along links. Nucleotide changes are specified by suffixes only for transversions; “⫹” indicates an insertion. Note that
we have redefined subclade U5a of Finnilä et al. (2001) and Herrnstadt et al. (2002) as “U5b2,” on the basis of its position in the phylogenetic
tree. U5b1b1 haplotypes are shown as the square labeled as “Saami motif” and are further refined in panel B. B, Phylogenetic network of 330
U5b1b1 lineages based on HVS-I sequences. The star indicates the basal node (transitions in nps 16144, 16189, and 16270). Population sizes
and U5b1b1 frequencies are shown in table 1. Ar p Armenians, Ba p Bashkirs, Bo p Bosnians, Ch p Chuvashes, Cr p Croats, Cz p Czechs,
Es p Estonians, Fi p Finns, Fr p French, Hu p Hungarians, Ka p Karelians, Kb p Kabardians, Ko p Komis, La p Latvians, Li p
Lithuanians, Ma p Maris, Mc p Moroccans, Mo p Mordvin, Ng p Nogays, No p Norwegians, Po p Poles, Sa p Saami, Sl p Slovaks,
Ru p Russians, Si p Sicilians, Sw p Swedes, Ta p Tatars, Uk p Ukrainians. For further information, see the legend to figure 1.
Tambets et al.: mtDNA and NRY Variation of the Saami 675
Table 5
Coalescence Ages for Different Subsets of Subhaplogroup U5b1b1
Tc DT d
Set of Sequences Considereda n rb (years) (years)
U5b1b1 (without Sa) 118 25/118 p .212 4,300 1,400
U5b1b1 in Scandinavia (without Sa) 55 6/55 p .109 2,200 2,200
U5b1b1 in Scandinavia (with Sa,
without subfounder 16148) 187 11/187 p .059 1,200 700
U5b1b1 in eastern Europe (Es/Li/
Ma/Mo/Ko/Ru/Uk/Po/Ba/Ch) 47 13/47 p .277 5,600 1,700
U5b1b1 in southern and western
Europe (Bo/Cr/Sl/Hu/Fr) 11 3/11 p .273 5,500 3,200
a
Population codes are the same as those used in figure 3.
b
Average mutational distance to the founder haplotype of the cluster.
c
Coalescence time, calculated by means of r with a mutation rate of 1 transition
per 20,180 years in the HVS-I segment between nps 16090 and 16365 (Forster et
al. 1996).
d
SD for r, calculated as in the study by Saillard et al. (2000).
over Europe, but it occurs in western and central Europe all non-Saami U5b1b1 lineages in Europe (see fig. 3B) is
with notable sequence variation. The haplotype diversity 4,300 Ⳳ 1,400 years BP (table 5). This is approximately
(excluding subclade U5b1b1) is 0.96. We note that the same as that for the eastern or western-southern
U5b1b was also identified in a northwestern African European subsets but is older than that observed for
population, among the Moroccans (fig. 3A). In contrast, Scandinavians.
its diversity in eastern Europe is much lower. There, the
haplotype diversity (excluding subclade U5b1b1) is 0.79, Discussion
whereas most of the U5b1b sequences in eastern Europe
belong to the U5b1b1 branch (fig. 3A and 3B). In the
Samoyed/Siberian Heritage of the Saami Population?
Eurasian cohort, U5b1b (other than U5b1b1) is absent
from native Siberians and is notably absent from the Uralic-speaking Samoyeds, Khants, Mansis, and Altaic-
Ob-Ugric populations and the Samoyeds. speaking Siberians virtually lack the European mtDNA
Most importantly, the data indicate that U5b1b1, the haplogroups V and U5b1b1 that predominate in the
only subcluster of U5b in the Saami population, is present Saami mtDNA pool (fig. 1; table 1). Eastern Eurasian
outside the Scandinavian-Baltic and the Volga-Uralic mtDNA variants in the Saami are represented by a re-
regions—namely, in the French, Croatian, Bosnian, Slov- stricted set of lineages that belong to superhaplogroup M.
enian, Czech, Russian, Ukrainian, Polish, and Hungarian In this respect, the Saami do not differ markedly from
mtDNA pools—and, further, is present even in the Cau- Finnic-speaking Karelians, Maris, Komis, Udmurts, or
casus, among the Nogay, Kabardinian, and Armenian northern Russians, all of whom possess haplogroups of
mtDNAs (fig. 3B). Thus, this subhaplogroup is much eastern Eurasian origin at similar frequencies (table 1).
more widely distributed than was believed previously (Sa- This minor part of the Saami mtDNA pool consists of
jantila et al. 1995). Two principal haplotypes, differing two branches of the eastern Eurasian mtDNA tree—D5
by one mutational step, define the Saami U5b1b1 (fig. and Z1. According to published data, the frequency of
3B). One, with transitions at nps 16144, 16189, and haplogroup D5 is relatively high in China (Yao et al.
16270, can be identified as the founder haplotype because 2002). D5 is also present among Mongols and Siberians
of its presence in many other populations. The second, (Kolman et al. 1996; Derbeneva et al. 2002b). However,
containing an additional transition at np 16148, is so far the Saami haplogroup D5 lineages, with the HVS-I mo-
exclusive to the Saami population. This subfounder com- tif 16126-16136-16360 and its derivatives (defined as
prises 38% of their U5b1b1 mtDNAs and is present in “D5b” by Derenko et al. 2003), have been identified only
all studied subpopulations of the Saami. Furthermore, no in some northern and eastern European populations
other derivatives of the U5b1b1 founder node present in (among Karelians, Finns, Estonians, North-Russians, and
the Saami sample was identified in other Scandinavians, Komis) and in some Siberian populations but not in Sam-
including the Finns, irrespective of the previous obser- oyeds (table 1). This suggests, again, the lack of gene flow
vation that the northern-central Finnish population is rel- from Samoyeds to eastern Europe.
atively rich in U5b1b1 (Meinilä et al. 2001). Thus, the Haplogroup Z, a subcluster of the M8 clade within
“leakage” of the Saami U5b1b1 to neighboring popula- the haplogroup M family of mtDNA (Kivisild et al.
tions seems to be rather limited. The coalescence time of 2002), is found at highest frequencies in the northeastern
676 Am. J. Hum. Genet. 74:661–682, 2004
Asians: the Itelmens and Koryaks (Schurr et al. 1999). than Siberia, is a possible origin of the earliest expansion
It is also present in several Siberian populations, includ- of this haplogroup in northern Eurasia. Third, the lack
ing the Altaic people (table 1). Though not identified in of Y-chromosomal haplogroup C in Saami contrasts
a large data set of the Yakuts (Fedorova et al. 2003; with its high frequency among Tungusic-speaking native
Pakendorf et al. 2003), haplogroup Z has been observed Siberians (such as the Evenks and the Evens) as well as
among several Finnic- and Turkic-speaking populations among Mongolic-speaking Mongols, the Buryats, the
of the Volga-Ural region (Bermisheva et al. 2002). It is Kalmyks, and the Turkic-speaking Kazakhs and Uzbeks
curious that it is more frequent there in Finnic- than (Wells et al. 2001; Karafet et al. 2002; authors’ unpub-
in Turkic-speaking populations. The absence of haplo- lished data). Therefore, without introducing specific ad-
group Z from most of the Siberian Uralic-speaking pop- ditional ad hoc scenarios, these observations make it
ulations (Samoyedic-speaking Nenets and Selkups, as unlikely that there was recent Y-chromosomal flow from
well as Siberian Ob-Ugric-speaking Khants and Mansis) these Siberian populations into the gene pool of the
(table 1) is therefore striking. We note that all haplo- Saami.
group Z lineages that are found in eastern Europe belong
to a subhaplogroup Z1, characterized by transitions at Origins of the Saami mtDNA and Y Chromosomes
nps 151, 10325, and 16129 within the Z phylogeny
(Kong et al. 2003a, 2003b). A matching HVS-I founder Scenarios involving extremes of genetic drift, such as
haplotype has been observed in the Koryak and the Itel- that due to repeated bottlenecks, could explain how the
men populations (Schurr et al. 1999). The limited di- Saami mtDNA pool evolved as a narrow subset of that
versity of haplogroup Z in Europe suggests its relatively found in other European populations. Indeed, there are
recent spread west of the Urals. good reasons to believe that even a much larger Finnish
The sister clade of Z—haplogroup C—is far more di- population went through several “bottlenecks” in its
verse and frequent than haplogroup Z in eastern Eu- demographic history (Nevanlinna 1972; de la Chapelle
rasians (Derenko et al. 2003), as well as in the popu- and Wright 1998; Kittles et al. 1999; Peltonen et al.
lations of the Volga-Ural region (Bermisheva et al. 2002), 2000). Likewise, it is possible that, when the proto-
but it is absent among the Saami. Here, again, the lack Saami gene pool was in statu nascendi, it was restricted
of haplogroup C in the Saami is specifically indicative to only a few basic mtDNA haplotypes that were carried
of the absence of possible genetic links between the by the founding settlers—a plausible scenario during the
Saami and the Samoyeds. Indeed, the only Samoyedic- “Paleolithic isolation.”
speaking population that also inhabits northeastern Eu- Torroni et al. (1998, 2001) have suggested that the
rope and is thus geographically close to the Saami, the spread of haplogroup V in Scandinavia and in eastern
Nenets, has haplogroup C as the most frequent (∼30% Europe is due to its late Pleistocene/early Holocene
frequency) variant of mtDNA. expansion from a Franco-Cantabrian glacial refugium.
The predominant Saami Y-chromosomal haplogroup Haplogroup V shows high diversity and the presence of
N3 has a nearly uniform circumarctic distribution in pre-V lineages in Iberia and the islands of Croatia (Tolk
Eurasia (table 3). The closely related N2 lineages are et al. 2000; Torroni et al. 2001). However, the new data
frequent in Siberian and Volga-Uralic populations. Thus, on Volga River basin and Slavic-speaking populations
it is likely that haplogroup N variation represents a pre- (Bermisheva et al. 2002; Malyarchuk et al. 2002; present
historic link between the Siberian and eastern European/ study) show that haplogroup V is also well present in
proto-Finnic populations via their paternal heritage. The eastern Europeans. Furthermore, haplogroup V lineages
improved resolution of the Y-chromosomal phylogenetic with HVS-I transitions 16153 and 16298 (fig. 1) that
tree (Jobling and Tyler-Smith 2003) reveals an ancestral are frequent in the Saami population are much more
node shared by haplogroups N and O, with the latter widespread in eastern than in western Europe (Torroni
restricted largely to eastern Asia. This connection is in- et al. 2001; Bermisheva et al. 2002; authors’ unpublished
triguing, but it is still unclear when and where this com- data). This indicates that haplogroup V might have
mon ancestor first appeared. Nevertheless, one does not reached Fennoscandia via central/eastern Europe (see fig.
need to postulate a recent Siberian flow of Y chromo- 4A). Such a scenario is indirectly supported by the ab-
somes into the Saami gene pool to explain their high N3 sence, among the Saami, of the pre-V mtDNAs that are
frequency. First, such a flow from Samoyedic-speaking characteristic of southwestern Europeans and north-
aboriginal Siberians to the Saami Y-chromosomal pool western Africans and that are also present in Germans
would predict the presence there of haplogroup N2 and/ (Torroni et al. 2001).
or haplogroup Q, widely spread in Samoyeds (Karafet The phylogeography and the ancestry of the other
et al. 2002). Second, the much higher diversity of N3 components of the Saami mitochondrial profile have so
in eastern Europe than in Siberia (Villems et al. 1998; far not been well understood. Here, we showed that
Rootsi et al. 2000) suggests that eastern Europe, rather mtDNA haplogroup U5b1b1—the set of lineages with
Tambets et al.: mtDNA and NRY Variation of the Saami 677
the northernmost Uralic-speaking populations—Sibe- bull DM, Howell N (1999) Reanalysis and revision of the
rian Ugric and Samoyedic speakers—or among any other Cambridge reference sequence for human mitochondrial
aboriginal Siberians. The Samoyeds are the least genet- DNA. Nat Genet 23:147
ically close to the Saami among the people of the Uralic Baasner A, Schafer C, Junge A, Madea B (1998) Polymorphic
sites in human mitochondrial DNA control region se-
language family, whereas nearly all of the mtDNA and
quences: population data and maternal inheritance. Forensic
Y-chromosomal heritage of the Saami can be adequately Sci Int 98:169–178
explained within the European pools of the two haploid Bandelt HJ, Lahermo P, Richards M, Macaulay V (2001) De-
genetic systems. This genetics-based reconstruction (fig. tecting errors in mtDNA data by phylogenetic analysis. Int
4) is in agreement with the reconstruction of the spread J Legal Med 115:64–69
of Ahrensburgian and Swiderian Mesolithic technologies Barac L, Pericic M, Klaric IM, Rootsi S, Janicijevic B, Kivisild
in northern Europe, linking it with population expansion T, Parik J, Rudan I, Villems R, Rudan P (2003) Y chro-
that can be likely traced back to the post–Last Glacial mosomal heritage of Croatian population and its island iso-
Maximum recolonization of the European north (Tor- lates. Eur J Hum Genet 11:535–542
roni et al. 2001; Tambets et al. 2003). The results also Beckman G, Beckman L, Sikstrom C (1993) Serum comple-
stress that the grouping of populations according to lan- ment (C3, BF, C4) types in Swedish Saamis. Hum Hered 43:
362–365
guage families should be used exclusively only in a lin-
Beckman L, Beckman G, Nylander PO (1988) Gc subtypes in
guistic context. Finns, Swedes and Swedish Lapps. Hum Hered 38:18–21
Bermisheva M, Tambets K, Villems R, Khusnutdinova E (2002)
Acknowledgments [Diversity of mitochondrial DNA haplotypes in ethnic pop-
ulations of the Volga-Ural region of Russia.] Mol Biol
We thank Tatyana Karafet and Boris Malyarchuk, for useful (Mosk) 36:990–1001
information; Henry Harpending, for the program POPSTR; Bertranpetit J, Sala J, Calafell F, Underhill PA, Moral P, Comas
Vincent Macaulay, for the program SAMPLING; Ille Hilpus D (1995) Human mitochondrial DNA variation and the or-
and Jaan Lind, for technical assistance; and Charles Kurland igin of Basques. Ann Hum Genet 59:63–81
and Thomas Gilbert, for helpful discussion and comments. We Cali F, Le Roux MG, D’Anna R, Flugy A, De Leo G, Chiavetta
are grateful to two anonymous reviewers for their suggestions V, Ayala GF, Romano V (2001) mtDNA control region and
and advice. The research of R.V. was supported by Estonian RFLP data for Sicily and France. Int J Legal Med 114:229–
basic research grant 514 and European Commission Direc- 231
torate General Research grant ICA1CT20070006. The re- Capelli C, Redhead N, Abernethy JK, Gratrix F, Wilson JF,
search of T.K. was supported by Estonian basic research grant Moen T, Hervig T, Richards M, Stumpf MP, Underhill PA,
5574. The work of E.K. was supported by the Russian Foun- Bradshaw P, Shaha A, Thomas MG, Bradman N, Goldstein
dation for Basic Research (project number 01-04-48487a) and DB (2003) A Y chromosome census of the British isles. Curr
the Ministry of Sciences and Technology of Russia. M.G., S.Z., Biol 13:979–984
and L.O. received support from expedition grants from the Casanova M, Leroy P, Boucekkine C, Weissenbach J, Bishop
Siberian Branch of the Russian Academy of Sciences (1992– C, Fellous M, Purrello M, Fiori G, Siniscalco M (1985) A
1997) and the Russian Foundation of Basic Research (project human Y-linked DNA polymorphism and its potential for
number 02-06-80524-a), and the research of P.R. received sup- estimating genetic and evolutionary distance. Science 230:
port from project number 0196005 of the Ministry of Science 1403–1406
and Technology of the Republic of Croatia. Cavalli-Sforza LL, Menozzi P, Piazza A (1994) The history
and geography of human genes. Princeton University Press,
Electronic-Database Information Princeton
Cavalli-Sforza LL, Piazza A (1993) Human genomic diversity
The URL for data presented herein is as follows: in Europe: a summary of recent research and prospects for
the future. Eur J Hum Genet 1:3–18
Fluxus Engineering, http://www.fluxus-engineering.com/ (for Comas D, Calafell F, Mateu E, Perez-Lezaun A, Bosch E, Mar-
Network 3.1.1.1) tinez-Arias R, Clarimon J, Facchini F, Fiori G, Luiselli D,
Pettener D, Bertranpetit J (1998) Trading genes along the
References silk road: mtDNA sequences and the origin of Central Asian
populations. Am J Hum Genet 63:1824–1838
Akey JM, Sosnoski D, Parra E, Dios S, Hiester K, Su B, Bonilla Corte-Real HB, Macaulay VA, Richards MB, Hariti G, Issad
C, Jin L, Shriver MD (2001) Melting curve analysis of SNPs MS, Cambon-Thomsen A, Papiha S, Bertranpetit J, Sykes
(McSNP): a gel-free and inexpensive approach for SNP ge- BC (1996) Genetic diversity in the Iberian Peninsula deter-
notyping. Biotechniques 30:358–362, 364, 366–367 mined from mitochondrial sequence analysis. Ann Hum Ge-
Anderson S, Bankier AT, Barrell BG, de Bruijn MH, Coulson net 60:331–350
AR, Drouin J, Eperon IC, Nierlich DP, Roe BA, Sanger F, Crespillo M, Luque JA, Paredes M, Fernandez R, Ramirez E,
Schreier PH, Smith AJ, Staden R, Young IG (1981) Sequence Valverde JL (2000) Mitochondrial DNA sequences for 118
and organization of the human mitochondrial genome. Na- individuals from northeastern Spain. Int J Legal Med 114:
ture 290:457–465 130–132
Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turn- Cruciani F, Santolamazza P, Shen P, Macaulay V, Moral P,
Tambets et al.: mtDNA and NRY Variation of the Saami 679
Olckers A, Modiano D, Holmes S, Destro-Bisol G, Coia V, and evolution of Native American mtDNA variation: a re-
Wallace DC, Oefner PJ, Torroni A, Cavalli-Sforza LL, Scoz- appraisal. Am J Hum Genet 59:935–945
zari R, Underhill PA (2002) A back migration from Asia to Francalacci P, Bertranpetit J, Calafell F, Underhill PA (1996)
sub-Saharan Africa is supported by high-resolution analysis Sequence diversity of the control region of mitochondrial
of human Y-chromosome haplotypes. Am J Hum Genet 70: DNA in Tuscany and its implications for the peopling of
1197–1214 Europe. Am J Phys Anthropol 100:443–460
Danan C, Sternberg D, Van Steirteghem A, Cazeneuve C, Du- Guglielmino CR, Piazza A, Menozzi P, Cavalli-Sforza LL
quesnoy P, Besmond C, Goossens M, Lissens W, Amselem (1990) Uralic genes in Europe. Am J Phys Anthropol 83:
S (1999) Evaluation of parental mitochondrial inheritance 57–68
in neonates born after intracytoplasmic sperm injection. Am Haetta OM (1996) The Sami: an indigenous people of the
J Hum Genet 65:463–473 Arctic. Davvi Girji, Kárásjohka/Karasjoki, Vaasa
de la Chapelle A, Wright FA (1998) Linkage disequilibrium Hammer MF, Horai S (1995) Y chromosomal DNA variation
mapping in isolated populations: the example of Finland and the peopling of Japan. Am J Hum Genet 56:951–962
revisited. Proc Natl Acad Sci USA 95:12416–12423. Helgason A, Hickey E, Goodacre S, Bosnes V, Stefansson K,
Delghandi M, Utsi E, Krauss S (1998) Saami mitochondrial Ward R, Sykes B (2001) mtDNA and the islands of the North
DNA reveals deep maternal lineage clusters. Hum Hered 48: Atlantic: estimating the proportions of Norse and Gaelic
108–114 ancestry. Am J Hum Genet 68:723–737
Derbeneva OA, Starikovskaya EB, Volod’ko NV, Wallace DC, Helgason A, Sigurethardottir S, Nicholson J, Sykes B, Hill EW,
Sukernik RI (2002a) [Mitochondrial DNA variation in Kets Bradley DG, Bosnes V, Gulcher JR, Ward R, Stefansson K
and Nganasans and the early peoples of Northern Eurasia]. (2000) Estimating Scandinavian and Gaelic ancestry in the
Genetika 38:1554–1560 male settlers of Iceland. Am J Hum Genet 67:697–717
Derbeneva OA, Starikovskaya EB, Wallace DC, Sukernik RI Herrnstadt C, Elson JL, Fahy E, Preston G, Turnbull DM,
(2002b) Traces of early Eurasians in the Mansi of northwest Anderson C, Ghosh SS, Olefsky JM, Beal MF, Davis RE,
Siberia revealed by mitochondrial DNA analysis. Am J Hum Howell N (2002) Reduced-median-network analysis of com-
Genet 70:1009–1014 plete mitochondrial DNA coding-region sequences for the
Derenko MV, Grzybowski T, Malyarchuk BA, Dambueva IK, major African, Asian, and European haplogroups. Am J
Denisova GA, Czarny J, Dorzhu CM, Kakpakov VT, Mis- Hum Genet 70:1152–1171
cicka-Sliwka D, Wozniak M, Zakharov IA (2003) Diversity Hofmann S, Jaksch M, Bezold R, Mertens S, Aholt S, Paprotta
of mitochondrial DNA lineages in south Siberia. Ann Hum A, Gerbitz KD (1997) Population genetics and disease
Genet 67:391–411 susceptibility: characterization of central European haplo-
Dimo-Simonin N, Grange F, Taroni F, Brandt-Casadevall C, groups by mtDNA gene mutations, correlation with D loop
Mangin P (2000) Forensic evaluation of mtDNA in a pop- variants and association with disease. Hum Mol Genet 6:
ulation from south west Switzerland. Int J Legal Med 113: 1835–1846
89–97 Jobling MA, Tyler-Smith C (2003) The human Y chromosome:
Di Rienzo A, Wilson AC (1991) Branching pattern in the evo- an evolutionary marker comes of age. Nat Rev Genet 4:
lutionary tree for human mitochondrial DNA. Proc Natl 598–612
Acad Sci USA 88:1597–1601 Kaessmann H, Heissig F, von Haeseler A, Pääbo S (1999) DNA
Donner J (1995) The quaternary history of Scandinavia. Cam- sequence variation in a non-coding region of low recom-
bridge University Press, Cambridge bination on the human X chromosome. Nat Genet 22:78–
Dupuy BM, Olaisen B (1996) mtDNA sequences in the Nor- 81
wegian Saami and main population. In: Carracedo A, Brink- Karafet TM, Osipova LP, Gubina MA, Posukh OL, Zegura
mann B, Bär W (eds) Advances in forensic haemogenetics, SL, Hammer MF (2002) High levels of Y-chromosome dif-
vol 6. Springer-Verlag, Berlin, Heidelberg, New York, pp ferentiation among native Siberian populations and the ge-
23–25 netic signature of a boreal hunter-gatherer way of life. Hum
Fedorova SA, Bermisheva MA, Villems R, Maksimova NR, Biol 74:761–789
Khusnutdinova EK (2003) [Analysis of mitochondrial DNA Kittles RA, Bergen AW, Urbanek M, Virkkunen M, Linnoila
lineages in Yakuts.] Mol Biol (Mosk) 37:643–653 M, Goldman D, Long JC (1999) Autosomal, mitochondrial,
Finnilä S, Hassinen IE, Ala-Kokko L, Majamaa K (2000) Phy- and Y chromosome DNA variation in Finland: evidence for
logenetic network of the mtDNA haplogroup U in northern a male-specific bottleneck. Am J Phys Anthropol 108:381–
Finland based on sequence analysis of the complete coding 399
region by conformation-sensitive gel electrophoresis. Am J Kivisild T, Rootsi S, Metspalu M, Mastana S, Kaldma K, Parik
Hum Genet 66:1017–1026 J, Metspalu E, Adojaan M, Tolk HV, Stepanov V, Golge M,
Finnilä S, Lehtonen MS, Majamaa K (2001) Phylogenetic net- Usanga E, Papiha SS, Cinnioglu C, King R, Cavalli-Sforza
work for European mtDNA. Am J Hum Genet 68:1475– L, Underhill PA, Villems R (2003) The genetic heritage of
1484 the earliest settlers persists both in Indian tribal and caste
Forster P, Cali F, Rohl A, Metspalu E, D’Anna R, Mirisola M, populations. Am J Hum Genet 72:313–332
De Leo G, Flugy A, Salerno A, Ayala G, Kouvatsi A, Villems Kivisild T, Tolk H-V, Parik J, Wang Y, Papiha SS, Bandelt H-
R, Romano V (2002) Continental and subcontinental dis- J, Villems R (2002) The emerging limbs and twigs of the
tributions of mtDNA control region types. Int J Legal Med East Asian mtDNA tree. Mol Biol Evol 19:1737–1751
116:99–108 Kolman C, Sambuughin N, Bermingham E (1996) Mitochon-
Forster P, Harding R, Torroni A, Bandelt H-J (1996) Origin drial DNA analysis of Mongolian populations and impli-
680 Am. J. Hum. Genet. 74:661–682, 2004
cations for the origin of New World founders. Genetics 142: Nygaard S (1989) The stone age of Northern Scandinavia: a
1321–1334 review. J World Prehistory 3:71–116
Kong QP, Yao YG, Liu M, Shen SP, Chen C, Zhu CL, Palan- Opdal SH, Rognum TO, Vege A, Stave AK, Dupuy BM, Ege-
ichamy MG, Zhang YP (2003a) Mitochondrial DNA se- land T (1998) Increased number of substitutions in the D-
quence polymorphisms of five ethnic populations from loop of mitochondrial DNA in the sudden infant death syn-
northern China. Hum Genet 113:391–405 drome. Acta Paediatr 87:1039–1044
Kong QP, Yao YG, Sun C, Bandelt HJ, Zhu CL, Zhang YP Orekhov V, Poltoraus A, Zhivotovsky LA, Spitsyn V, Ivanov
(2003b) Phylogeny of East Asian mitochondrial DNA line- P, Yankovsky N (1999) Mitochondrial DNA sequence di-
ages inferred from complete sequences. Am J Hum Genet versity in Russians. FEBS Letters 445:197–201
73:671–676 Pakendorf B, Wiebe V, Tarskaia LA, Spitsyn VA, Soodyall H,
Kozlowski J, Bandi HG (1984) The paleohistory of circum- Rodewald A, Stoneking M (2003) Mitochondrial DNA ev-
polar arctic colonization. Arctic 37:359–372 idence for admixed origins of central Siberian populations.
Lahermo P, Sajantila A, Sistonen P, Lukka M, Aula P, Peltonen Am J Phys Anthropol 120:211–224
L, Savontaus ML (1996) The genetic relationship between Parson W, Parsons TJ, Scheithauer R, Holland MM (1998)
the Finns and the Finnish Saami (Lapps): analysis of nuclear Population data for 101 Austrian Caucasian mitochondrial
DNA and mtDNA. Am J Hum Genet 58:1309–1322 DNA d-loop sequences: application of mtDNA sequence
Laitinen V, Lahermo P, Sistonen P, Savontaus ML (2002) Y- analysis to a forensic case. Int J Legal Med 111:124–132
chromosomal diversity suggests that Baltic males share com- Passarino G, Cavalleri GL, Lin AA, Cavalli-Sforza LL, Bor-
mon Finno-Ugric-speaking forefathers. Hum Hered 53:68– resen-Dale AL, Underhill PA (2002) Different genetic com-
78 ponents in the Norwegian population revealed by the anal-
Lutz S, Weisser HJ, Heizmann J, Pollak S (1998) Location and ysis of mtDNA and Y chromosome polymorphisms. Eur J
frequency of polymorphic positions in the mtDNA control Hum Genet 10:521–529
region of individuals from Germany. Int J Legal Med 111: Peltonen L, Palotie A, Lange K (2000) Use of population iso-
67–77 lates for mapping complex traits. Nat Rev Genet 1:182–190
Macaulay VA, Richards MB, Hickey E, Vega E, Cruciani F, Pereira L, Prata MJ, Amorim A (2000) Diversity of mtDNA
Guida V, Scozzari R, Bonne-Tamir B, Sykes B, Torroni A lineages in Portugal: not a genetic edge of European vari-
(1999) The emerging tree of West Eurasian mtDNAs: a syn- ation. Ann Hum Genet 64:491–506
thesis of control-region sequences and RFLPs. Am J Hum Pfeiffer H, Brinkmann B, Huhne J, Rolf B, Morris AA, Steighner
Genet 64:232–249 R, Holland MM, Forster P (1999) Expanding the forensic
Malyarchuk BA, Derenko MV (2001) Mitochondrial DNA German mitochondrial DNA control region database: genetic
variability in Russians and Ukrainians: implications to the diversity as a function of sample size and microgeography.
origin of the Eastern Slavs. Ann Hum Genet 65:63–78 Int J Legal Med 112:291–298
Malyarchuk BA, Grzybowski T, Derenko MV, Czarny J, Pfeiffer H, Forster P, Ortmann C, Brinkmann B (2001) The
Drobnic K, Miscicka-Sliwka D (2003) Mitochondrial DNA results of an mtDNA study of 1,200 inhabitants of a German
variability in Bosnians and Slovenians. Ann Hum Genet 67: village in comparison to other Caucasian databases and its
412–425 relevance for forensic casework. Int J Legal Med 114:169–
Malyarchuk BA, Grzybowski T, Derenko MV, Czarny J, Woz- 172
niak M, Miscicka-Sliwka D (2002) Mitochondrial DNA var- Piercy R, Sullivan KM, Benson N, Gill P (1993) The appli-
iability in Poles and Russians. Ann Hum Genet 66:261–283 cation of mitochondrial DNA typing to the study of white
Mathias N, Bayes M, Tyler-Smith C (1994) Highly informative Caucasian genetic identification. Int J Legal Med 106:85–
compound haplotypes for the human Y chromosome. Hum 90
Mol Genet 3:115–123 Pult I, Sajantila A, Simanainen J, Georgiev O, Schaffner W,
Meinilä M, Finnilä S, Majamaa K (2001) Evidence for mtDNA Pääbo S (1994) Mitochondrial DNA sequences from Swit-
admixture between the Finns and the Saami. Hum Hered zerland reveal striking homogeneity of European popula-
52:160–170 tions. Biol Chem Hoppe Seyler 375:837–840
Metspalu E, Kivisild T, Kaldma K, Parik J, Reidla M, Tambets Puzyrev VP, Stepanov VA, Golubenko MV, Puzyrev KV, Max-
K, Villems R (1999) The trans-Caucasus and the expansion imova NR, Kharkov VN, Spiridonova MG, Nogovitsyna
of the Caucasoid-specific human mitochondrial DNA. In: Pa- AN (2003) mtDNA and Y-chromosome lineages in the Ya-
piha S, Deka R, Chakraborty R (eds) Genomic diversity: ap- kut population. Genetika 39:975–981
plication in human population genetics. Kluwer Academic/ Raitio M, Lindroos K, Laukkanen M, Pastinen T, Sistonen
Plenum Publishers, New York, pp 121–134 P, Sajantila A, Syvanen A (2001) Y-chromosomal SNPs
Mogentale-Profizi N, Chollet L, Stevanovitch A, Dubut V, in Finno-Ugric-speaking populations analyzed by minise-
Poggi C, Pradie MP, Spadoni JL, Gilles A, Beraud-Colomb quencing on microarrays. Genome Res 11:471–482
E (2001) Mitochondrial DNA sequence diversity in two Richards M, Corte-Real H, Forster P, Macaulay V, Wilkinson-
groups of Italian Veneto speakers from Veneto. Ann Hum Herbots H, Demaine A, Papiha S, Hedges R, Bandelt HJ,
Genet 65:153–166 Sykes B (1996) Paleolithic and Neolithic lineages in the Eu-
Nei M (1987) Molecular evolutionary genetics. Columbia Uni- ropean mitochondrial gene pool. Am J Hum Genet 59:185–
versity Press, New York, pp 145–163 203
Nevanlinna HR (1972) The Finnish population structure. A Richards M, Macaulay V, Hickey E, Vega E, Sykes B, Guida
genetic and genealogical study. Hereditas 71:195–236 V, Rengo C, et al (2000) Tracing European founder lineages
Tambets et al.: mtDNA and NRY Variation of the Saami 681
in the Near Eastern mtDNA pool. Am J Hum Genet 67: Shen P, Wells RS (2003) A novel Y-chromosome variant puts
1251–1276 an upper limit on the timing of first entry into the Americas.
Richards M, Macaulay V, Torroni A, Bandelt HJ (2002) In Am J Hum Genet 73:700–705
search of geographical patterns in European mitochondrial Semino O, Passarino G, Oefner PJ, Lin AA, Arbuzova S, Beck-
DNA. Am J Hum Genet 71:1168–1174 man LE, De Benedictis G, Francalacci P, Kouvatsi A, Lim-
Richards MB, Macaulay VA, Bandelt H-J, Sykes BC (1998) borska S, Marcikiae M, Mika A, Mika B, Primorac D, San-
Phylogeography of mitochondrial DNA in western Europe. tachiara-Benerecetti AS, Cavalli-Sforza LL, Underhill PA
Ann Hum Genet 62:241–260 (2000) The genetic legacy of Paleolithic Homo sapiens sap-
Rootsi S, Kivisild T, Tambets K, Adojaan M, Parik J, Reidla iens in extant Europeans: a Y chromosome perspective. Sci-
M, Metspalu E, Laos S, Tolk HV, Villems R (2000) On the ence 290:1155–1159
phylogeographic context of sex-specific genetic markers of Sumkin VJ (1990) On the ethnogenesis of the Saami: an ar-
Finno-Ugric populations. In: Künnap A (ed) The roots of chaeological view. Acta Borealia 7:3–20
peoples and languages of Northern Eurasia II and III. Uni- Tagliabracci A, Turchi C, Buscemi L, Sassaroli C (2001) Poly-
versity of Tartu, Division of Uralic Languages/Societas His- morphism of the mitochondrial DNA control region in Ital-
toriae Fenno-Ugricae, Tartu, pp 148–164 ians. Int J Legal Med 114:224–228
Rosser ZH, Zerjal T, Hurles ME, Adojaan M, Alavantic D, Tambets K, Kivisild T, Metspalu E, Parik J, Kaldma K, Laos
Amorim A, Amos W, et al (2000) Y-chromosomal diversity S, Tolk HV, Gölge M, Demirtas H, Geberhiwot T, Papiha
in Europe is clinal and influenced primarily by geography, SS, De Stefano GF, Villems R (2000) The topology of the
rather than by language. Am J Hum Genet 67:1526–1543 maternal lineages of the Anatolian and Trans-Caucasus pop-
Rousselet F, Mangin P (1998) Mitochondrial DNA polymor- ulations and the peopling of the Europe: some preliminary
phisms: a study of 50 French Caucasian individuals and considerations. In: Renfrew C, Boyle K (eds) Archaeoge-
application to forensic casework. Int J Legal Med 111:292– netics: DNA and the population prehistory of Europe.
298 McDonald Institute for Archaeological Research Mono-
Saillard J, Evseva I, Tranebjaerg L, Norby S (2000) Mito- graph Series, Cambridge University Press, Cambridge, pp
chondrial DNA diversity among Nenets. In: Renfrew C, 219–235
Boyle K (eds) Archaeogenetics: DNA and the population Tambets K, Tolk HV, Kivisild T, Metspalu E, Parik J, Reidla
prehistory of Europe. McDonald Institute for Archaeolog- M, Voevoda M, Damba L, Bermisheva M, Khusnutdinova
ical Research Monograph Series, Cambridge University E, Golubenko M, Stepanov V, Puzyrev V, Usanga E, Rudan
Press, Cambridge, pp 255–258 P, Beckman L, Villems R (2003) Complex signals for pop-
Saillard J, Forster P, Lynnerup N, Bandelt H-J, Nørby S (2000) ulation expansions in Europe and beyond. In: Bellwood P,
mtDNA variation among Greenland Eskimos: the edge of Renfrew C (eds) Examining the farming/language dispersal
the Beringian expansion. Am J Hum Genet 67:718–726 hypothesis. McDonald Institute for Archaeological Research
Sajantila A, Lahermo P, Anttinen T, Lukka M, Sistonen P, Monograph Series, Cambridge University Press, Cambridge,
Savontaus ML, Aula P, Beckman L, Tranebjaerg L, Gedde- pp 449–458
Dahl T, Issel-Tarver L, DiRienzo A, Pääbo S (1995) Genes Tolk HV, Pericic M, Barac L, Klaric IM, Janicijevic B, Rudan
and languages in Europe: an analysis of mitochondrial line- I, Parik J, Villems R, Rudan P (2000) mtDNA haplogroups
ages. Genome Res 5:42–52 in the populations of Croatian Adriatic Islands. Coll An-
Sajantila A, Pääbo S (1995) Language replacement in Scan- tropol 24:267–280
dinavia. Nat Genet 11:359–360 Torroni A, Bandelt H-J, D’Urbano L, Lahermo P, Moral P,
Sajantila A, Salem AH, Savolainen P, Bauer K, Gierig C, Pääbo Sellitto D, Rengo C, et al (1998) mtDNA analysis reveals a
S (1996) Paternal and maternal DNA lineages reveal a bot- major late Paleolithic population expansion from south-
tleneck in the founding of the Finnish population. Proc Natl western to northeastern Europe. Am J Hum Genet 62:1137–
Acad Sci USA 93:12035–12039 1152
Salas A, Comas D, Lareu MV, Bertranpetit J, Carracedo A Torroni A, Bandelt H-J, Macaulay V, Richards M, Cruciani F,
(1998) mtDNA analysis of the Galician population: a genetic Rengo C, Martinez-Cabrera V, et al (2001) A signal, from
edge of European variation. Eur J Hum Genet 6:365–375 human mtDNA, of postglacial recolonization in Europe. Am
Sambrook J, Fritsch EF, Maniatis T (1989) Molecular cloning: J Hum Genet 69:844–852
a laboratory manual. Cold Spring Harbor Laboratory Press, Torroni A, Huoponen K, Francalacci P, Petrozzi M, Morelli
Cold Spring Harbor, NY L, Scozzari R, Obinu D, Savontaus ML, Wallace DC (1996)
Sammallahti P (1998) The Saami languages: an introduction. Classification of European mtDNAs from an analysis of
Davvi Girji, Kárásjohka/Karasjoki, Vaasa three European populations. Genetics 144:1835–1850
Sanchez JJ, Børsting C, Hallenberg C, Buchard A, Hernandez Torroni A, Miller JA, Moore LG, Zamudio S, Zhuang J,
A, Morling N (2003) Multiplex PCR and minisequencing Droma T, Wallace DC (1994) Mitochondrial DNA analysis
of SNPs—a model with 35 Y chromosome SNPs. Forensic in Tibet: implications for the origin of the Tibetan popu-
Sci Int 137:74–84 lation and its adaptation to high altitude. Am J Phys An-
Schurr TG, Sukernik RI, Starikovskaya YB, Wallace DC (1999) thropol 93:189–199
Mitochondrial DNA variation in Koryaks and Itel’men: Torroni A, Schurr TG, Yang CC, Szathmary EJ, Williams RC,
population replacement in the Okhotsk Sea-Bering Sea re- Schanfield MS, Troup GA, Knowler WC, Lawrence DN,
gion during the Neolithic. Am J Phys Anthropol 108:1–39 Weiss KM, Wallace DC (1992) Native American mitochon-
Seielstad M, Yuldasheva N, Singh N, Underhill P, Oefner P, drial DNA analysis indicates that the Amerind and the Na-
682 Am. J. Hum. Genet. 74:661–682, 2004
dene populations were founded by two independent migra- Whitfield LS, Sulston JE, Goodfellow PN (1995) Sequence var-
tions. Genetics 130:153–162 iation of the human Y chromosome. Nature 378:379–380
Underhill PA, Passarino G, Lin AA, Shen P, Mirazon Lahr M, Yao YG, Kong QP, Bandelt HJ, Kivisild T, Zhang YP (2002)
Foley R, Oefner PJ, Cavalli-Sforza LL (2001) The phylo- Phylogeographic differentiation of mitochondrial DNA in
geography of Y chromosome binary haplotypes and the or- Han Chinese. Am J Hum Genet 70:635–651
igins of modern human populations. Ann Hum Genet 65: Y Chromosome Consortium (2002) A nomenclature system
43–62 for the tree of human Y-chromosomal binary haplogroups.
Villems R, Adojaan M, Kivisild T, Metspalu E, Parik J, Pielberg Genome Res 12:339–348
G, Rootsi S, Tambets K, Tolk HV (1998) Reconstruction of Zerjal T, Beckman L, Beckman G, Mikelsaar AV, Krumina A,
maternal lineages of Finno-Ugric speaking people and some
Kucinskas V, Hurles ME, Tyler-Smith C (2001) Geograph-
remarks on their paternal inheritance. In: Wiik K, Julku K
ical, linguistic, and cultural influences on genetic diversity:
(eds) The roots of peoples and languages of Northern Eur-
Y-chromosomal distribution in Northern European popu-
asia I. Societas Historiae Fenno-Ugricae, Turku, pp 180–
lations. Mol Biol Evol 18:1077–1087
200
Weale ME, Weiss DA, Jager RF, Bradman N, Thomas MG Zerjal T, Dashnyam B, Pandya A, Kayser M, Roewer L, Santos
(2002) Y chromosome evidence for Anglo-Saxon mass mi- FR, Schiefenhovel W, Fretwell N, Jobling MA, Harihara S,
gration. Mol Biol Evol 19:1008–1021 Shimizu K, Semjidmaa D, Sajantila A, Salo P, Crawford MH,
Wells RS, Yuldasheva N, Ruzibakiev R, Underhill PA, Evseeva Ginter EK, Evgrafov OV, Tyler-Smith C (1997) Genetic re-
I, Blue-Smith J, Jin L, et al (2001) The Eurasian heartland: lationships of Asians and Northern Europeans, revealed by
a continental perspective on Y-chromosome diversity. Proc Y-chromosomal DNA analysis. Am J Hum Genet 60:1174–
Natl Acad Sci USA 98:10244–10249 1183