While numerous ancient human DNA datasets from across Europe have been published till date, modern-day Poland in particular, remains uninvestigated. Besides application in the reconstruction of continent-wide human history, data from this region would also contribute towards our understanding of the history of the Slavs, whose origin is hypothesized to be in East or Central Europe. Here, we present the first population-scale ancient human DNA study from the region of modern-day Poland by establishing mitochondrial DNA profiles for 23 samples dated to 200 BC – 500 AD (Roman Iron Age) and for 20 samples dated to 1000–1400 AD (Medieval Age). Our results show that mitochondrial DNA sequences from both periods belong to haplogroups that are characteristic of contemporary West Eurasia. Haplotype sharing analysis indicates that majority of the ancient haplotypes are widespread in some modern Europeans, including Poles. Notably, the Roman Iron Age samples share more rare haplotypes with Central and Northeast Europeans, whereas the Medieval Age samples share more rare haplotypes with East-Central and South-East Europeans, primarily Slavic populations. Our data demonstrates genetic continuity of certain matrilineages (H5a1 and N1a1a2) in the area of present-day Poland from at least the Roman Iron Age until present. As such, the maternal gene pool of present-day Poles, Czechs and Slovaks, categorized as Western Slavs, is likely to have descended from inhabitants of East-Central Europe during the Roman Iron Age.
Citation: Juras A, Dabert M, Kushniarevich A, Malmström H, Raghavan M, Kosicki JZ, et al. (2014) Ancient DNA Reveals Matrilineal Continuity in Present-Day Poland over the Last Two Millennia. PLoS ONE 9(10): e110839. https://doi.org/10.1371/journal.pone.0110839
Editor: David Caramelli, University of Florence, Italy
Received: May 20, 2014; Accepted: September 20, 2014; Published: October 22, 2014
Copyright: © 2014 Juras et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: The authors confirm that all data underlying the findings are fully available without restriction. All relevant data are within the paper and its Supporting Information files.
Funding: Funding provided by Polish Ministry of Science and Higher Education; grant N N303 406836, Estonian Research Council IUT24-1 and the European Regional Development Fund (European Union) through the Centre of Excellence in Genomics to Estonian Biocentre and University of Tartu. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Continuity of human occupation in the territory of Central Europe, modern-day Poland in particular, and its relation to the origins of the Slavs have been widely discussed in the archaeological, linguistic and historical literature; however, these questions still remain contentious –. At present, vast territories of East-Central and South-East Europe are inhabited by Slavic populations . Three groups of present-day Slavs are identified on the basis of their linguistic affinities: Western Slavs (Poles, Czechs and Slovaks), Eastern Slavs (Ukrainians, Belarusians and Russians) and Southern Slavs (Croatians, Bulgarians, Slovenians, Bosnians, Macedonians, Montenegrins and Serbians) . It is supposed that all Slavs, besides their linguistic affinity, also share a common place of origin, although the latter is still inconclusive .
Several hypotheses have been advanced regarding the origin and early migrations of Slavs, of which two – autochthonous and allochthonous – have predominated. According to the autochthonous hypothesis, territories around Oder and Vistula rivers (in present-day Poland) were continuously inhabited by ancestors of Slavs from the Roman Iron Age (0–400 AD), or perhaps even further back in time from the Bronze Age (3200–600 BC)  until the Medieval Age (500–1500 AD) . In contrast, the allochthonous theory suggests the discontinuity of settlements between Roman Iron Age and Medieval Age in the territory of present-day Poland. Allochthonists hypothesize that the Slavs originated in the Pripet and Middle Dnieper River basins in modern-day Ukraine, from where they migrated to the west and south of Europe in the beginning of 5th century AD and inhabited the lands of present-day Poland, which was previously occupied by Germanic tribes during the Roman Iron Age . However, morphological analyses of skeletal materials from present-day Poland have suggested a continuity between Roman Iron Age (represented by Przeworsk and Wielbark cultures) and Medieval Age populations ,  thus providing less support to the allochthonous model.
Genetic studies on present-day Slavic-speaking populations have also addressed the complex genetic history of the Slavs –. The comparison of the complete mitochondrial genome sequences revealed a number of lineages that seem specific for Central and Eastern Europe. Moreover, based on age estimations, the authors suggest a genetic continuity of some Slavic mitochondrial lineages from at least the Bronze Age .
Ancient DNA (aDNA) provides direct genetic evidence for past demographic events. Mitochondrial DNA (mtDNA) from skeletal remains has been particularly successful in reconstructing the evolutionary history of European populations (e.g. –). However, no large-scale aDNA study on putative ancestral populations of modern-day Slavs have been reported thus far. Ancient DNA datasets from regions geographically adjacent to present-day Poland are limited to Iron and Middle Age samples from Denmark –, Neolithic samples e.g. from the Linearbandkeramik culture from Germany  and Bronze Age samples from Ukraine, Bulgaria and Moldova .
Therefore, to provide fresh perspectives on the debate of genetic continuity in Central Europe during the last two millennia, and to contribute to the resolution of the complex origin of the Slavs, we present the first population-level ancient mtDNA analysis on samples originating from six archaeological sites in Poland. The studied samples date to the Roman Iron Age, represented by Wielbark and Przeworsk cultures and to the Medieval Age. Our aim is to determine matrilineal genetic structure of ancient populations in the area comprising contemporary Poland, their relationships to one other and to other ancient and modern human populations from Europe, and to investigate potential genetic continuity between populations spanning two millennia.
Materials and Methods
Archaeological sites and samples
The skeletal material studied here originated from burial sites located in present-day Poland dating to the Roman Iron Age and the Medieval Age (n = 72) (Figure 1). The Roman Iron Age samples (RoIA) comprised 38 human remains from cemeteries in Kowalewko (n = 11) and Rogowo (n = 13) assigned to the Wielbark culture, and two burial sites in Karczyn (n = 12) and Gąski (n = 2) assigned to the Przeworsk culture. The Wielbark culture extended in the north-eastern territories of contemporary Poland during the 1st to the 4th century AD. The Przeworsk culture was present in the Western, Central and Southern Poland from the 3rd century BC to the 5th century AD. Both Wielbark and Przeworsk cultures were dated to the Roman Iron Age based on the archaeological context. Medieval Age (ME) samples comprised 34 human remains recovered from cemeteries in Cedynia (n = 18) and Ostrów Lednicki (n = 16). Detailed information about each sample, the archaeological context of the burial sites and their geographic origins is presented in Table 1 and Table S1. The handing history of samples is not well recognized; however seems to be minimal due to obtained results since aDNA was obtained only from intact teeth that are thought to be less prone to modern human DNA contaminations than other skeletal parts . The permission for collecting samples for aDNA studies from all mentioned above archaeological sites was provided by the supervisor of the skeletal materials, Head of the Institute of Anthropology and Department of Human Evolutionary Biology at Adam Mickiewicz University in Poznan.
Figure 1. Locations of the Roman Iron Age (RoIA) and Medieval Age (ME) burial sites in the territory of present-day Poland.
RoIA (stars): Kowalewko (K), Karczyn (KA), Gąski (G), Rogowo (R). ME (triangles): Cedynia (C), Ostrów Lednicki (OL).
Table 1. Geographical origins and archaeological contexts of ancient samples analyzed in the present study.
Extraction of ancient DNA
Sampling was performed using disposable gloves, facemasks and body suits to minimize the risk of contamination from modern humans. Two teeth were collected from each individual. All pre-PCR work was conducted in laboratories dedicated exclusively to the analysis of low copy number DNA (Centre for GeoGenetics at the University of Copenhagen, Denmark and Ancient DNA Laboratory at the Adam Mickiewicz University in Poznan, Poland). The outer surface of the teeth was decontaminated using 0.1 M HCl, followed by drilling of the teeth and treatment of the powder with 0.5% NaOCl . DNA was extracted using a silica-column based method ,  which was modified by the addition of urea to the extraction buffer . DNA extractions from the two teeth from each individual were performed at separate times. Preparation of reagents and solutions was conducted under sterile conditions and appropriate precautions (e.g. UV irradiation) were taken to avoid modern DNA contamination. Negative controls were set up during extractions (one control for every four samples) and amplifications (one control for every eight PCR reactions). Three faunal samples, contemporary with the human skeletal samples, were retrieved from one of the locations (Kowalewko) and were used as controls to screen for contamination from modern human sources both during and post-excavation.
Mitochondrial DNA analysis
Seven sets of overlapping primer pairs (Table S2) were used to amplify 360 base pairs (bp) of the first hypervariable region (HVRI) of the mtDNA control region between nucleotide positions (nps) 16043-16403, according to the revised Cambridge Reference Sequence (rCRS) (NC_012920.1) . Haplogroup-diagnostic nps in HVRII and mtDNA coding region were amplified using eighteen primer pairs, respectively, and HVRI regions between nps 16043-16132 and 16307-16403 were amplified using M13-tailed primers (Table S2). PCR reactions were set up as follows: 2 µl DNA extract, 0.5 U Platinum Taq DNA Polymerase High Fidelity (Invitrogen), 1X High-Fidelity PCR Buffer, 2 mM MgSO4, 0.8 mg/ml RSA (Calbiochem), 200 µM each of dNTPs (Invitrogen), 500 nM of each primer and ddH2O up to 25 µl. The thermocycling conditions were as follows: initial denaturation at 94°C for 4 minutes; 42 cycles of 94°C for 30 seconds, 52°C–60°C depending on the primer pair, (Table S2) for 20 seconds, 68°C for 20 seconds; and final extension at 72°C for 10 minutes. PCR products were visualized through electrophoresis on 2% agarose gel.
Almost all analyzed fragments, including HVRI, HVR-II and coding region fragments containing diagnostic SNPs were cloned and sequenced. Only fragments comprising nps 16048-16132 and 16307-16403 were amplified with M13-tailed primers and directly sequenced without cloning. Amplicons were cloned into pCR2.1-TOPO vector and transformed into competent E. coli cells (One Shot E. coli) using TOPO TA cloning kit (Invitrogen), following the supplier’s instructions. More than twelve bacterial colonies from each cloning experiment were screened for inserts using M13 universal primers. At least four positive clones were sequenced for each case. PCR products were purified with exonuclease I and Fast alkaline phosphatase (Fermentas) and sequenced using BigDye Terminator v3.1 kit and ABI Prism 3130xl Genetic Analyzer (Applied Biosystems), following manufacturer’s instructions. Amplicons generated with M13-tailed primers were sequenced in both directions using M13 universal primers.
In order to determine consensus sequence and detect post-mortem damages or/and possible contaminations, alignment of mtDNA sequences was performed using BioEdit v.18.104.22.168 (http://www.mbio.ncsu.edu/BioEdit/bioedit.html). Polymorphic positions in HVRI, HVRII and coding region sequences were scored against the rCRS sequence. Haplogroups were assigned following the established hierarchy of mtDNA phylogeny .
Populations used in comparative analyses
In order to compare mtDNA profiles of RoIA and ME samples to modern populations, we compiled a database of published mtDNA diversity in contemporary Poles, Czechs, Slovaks, Belarusians, Ukrainians, Russians (from the European part), Slovenians, Bosnians, Croatians, Bulgarians, Serbs, Germans, Finns, Estonians, Latvians, Lithuanians, Macedonians and Swedes (Table S3). To determine genetic distances between ancient populations, we also compiled a dataset of published mtDNA diversity in populations from Iron Age, Early Christian period and Middle Age from present-day Denmark (Germanic tribes), and Neolithic population from Linear Pottery Culture (LBK) from Germany (Table S3).
For haplotype sharing analysis we considered sequences bounded by np 16043 to 16403 of the HVRI for both aDNA samples in this study and those from modern populations of our database. All matches between haplotypes were scored taking into account their haplogroup affiliation. Because different studies have used different phylogenetic resolution of mtDNA, for haplotype sharing analysis and for MDS we considered the phylogenetic resolution depth reached in our aDNA samples.
Haplotype sharing analysis
Haplotype sharing analysis was conducted in order to detect mtDNA haplotypes shared between RoIA and ME samples and modern Europeans. To this end, eighteen contemporary populations representing East-Central, Southeast and North Europe were used in the analysis (Table S3). Several populations (Slovaks/Czechs, Macedonians/Serbs, Bosnians/Slovenians/Croatians, Finns/Estonians and Lithuanians/Latvians) were pooled together so that the total sample size (N) for each population amounted to approximately 300 individuals (between 277 and 317). Pooling was performed according to the linguistic affinities and geographic location of populations. Likewise, the sample sizes for populations with N>300 samples were decreased by sampling 300 randomly selected individuals. Consequently, a total of twelve populations or groups of pooled populations were used in the analysis.
The presence or absence of a particular haplotype in a given contemporary population was marked as „1” or „0”, respectively, and the total number as well as frequency of ancient haplotypes found in each modern population was calculated. Furthermore, all haplotypes present in RoIA and ME samples were divided into three classes based on their incidence in the comparative modern populations: informative haplotypes, which were found in less than half of the comparative populations; non-informative haplotypes, which occurred in more than half of the comparative populations; and unique haplotypes, which did not have exact matches in contemporary populations. In cases where we reduced sample sizes, we have cross-checked for informative/unique haplotypes if they present in whole samples size.
We used the two-tailed z-test to assess statistical significance of shared informative haplotypes between ancient and modern populations . Nonparametric bootstrapping of 1000 replicates for each population was used to generate the confidence intervals for the percentage of all matches, informative matches, and non-informative matches. The analysis was performed in R (R Development Core Team 2013)  using boot library .
Population pairwise FST
MtDNA haplogroup frequencies were used to determine genetic distances between ancient and modern populations. MtDNA haplogroups of comparative modern Slavic groups as well as Iron Age and Middle Age populations from present-day Denmark – and Neolithic population (LBK) from present-day Germany  were used in this analysis. Slatkin’s linearized pairwise FST were calculated using Arlequin v.3.5 . Multidimensional Scaling (MDS) plot was built based on population pairwise FST values using Statistica v.10 StatSoft (2011).
Authenticity of ancient DNA results
Reproducible mtDNA sequences were obtained from 23 out of the 38 specimens from Roman Iron Age burial sites (Table 2), and from 20 out of the 34 samples from the Medieval Age sites (Table 3). The consensus nucleotide sequences were supported by an alignment of at least 10 clones deriving from partially overlapping amplicons obtained from at least two independent DNA extractions (nps 16050-16130, 16119-16196, 16181-16226 and 16209-16356) (Figure S1). Therefore, if the same polymorphic positions were detected in cloned sequences retrieved from two teeth of each individual, the results were considered to be authentic. In four individuals we did not obtain a fragment 16209-16356 from one of the two DNA samples extracted from the same individual. Thus a shorter fragment of HVRI (np 16249-16317) was amplified, cloned and sequenced. For eleven individuals a fragment comprising nps 16050-16130 was obtained. For samples which had been confirmed by cloning as containing aDNA free from contaminations of modern human DNA a direct sequencing of M13-tailed amplicons comprising nps 16048-16132 and 16307-16403 was performed. Twenty nine out of 72 samples failed during the PCR amplifications or showed inconsistent alignment of cloned sequences and were discarded from the analyses.
Table 2. MtDNA haplogroups (hg) identified in Roman Iron Age populations.
Table 3. MtDNA haplogroups identified in Medieval populations.
All cloned mtDNA inserts showed characteristic post-mortem ancient DNA damages, of which 99% wereby C>T and G>A transitions –. Sequences obtained with the same procedure in the two independent laboratories were consistent. No evidence of modern human DNA contamination was present in the faunal material used as negative controls.
MtDNA haplogroup composition in ancient populations
All haplotypes identified in RoIA and ME populations belonged to typical West Eurasian mtDNA haplogroups (hg): H, K, U, T, J, W, HV, X, HV0, R0, and N (Table 2 and 3). There were 18 different mtDNA haplotypes found among 23 RoIA individuals and 20 distinct mtDNA haplotypes identified among 20 ME samples.
Haplogroup (hg) H was the most abundant hg in the RoIA populations (60.9% of all genotyped RoIA individuals). Samples from Rogowo comprised the highest number of individuals (10 out of 12) assigned to hg H, of which three displayed HVRI sequence identical to the rCRS. Three haplotypes, originating from Kowalewko and Gąski burial sites, were assigned to three sub-branches of hg U (U5a, U5b and U3) (Table 2). Haplogroup W was only identified in three individuals from the Kowalewko site. Remaining hgs (T, J2a and N1a) were identified in one individual each.
Haplogroup (hg) H was also the most frequent hg among the ME samples (40% of all analyzed ME individuals) (Table 3). Other frequent hgs in the ME population were K, J and X, each occurring in two specimens. Haplogroups T, W, N1b, HV, HVO, and R0a were only observed in one individual each.
MtDNA haplotypes shared between ancient and modern populations
A set of 3595 modern mitochondrial haplotypes, including majority of Slavic populations from regions geographically adjacent to Poland, was used for comparative analysis. Haplotypes shared between ancient and modern populations are presented in Table S4.
Among haplotypes identified in the RoIA populations, seven were found to be frequent in most of the modern-day populations in our database and thus classified as non-informative. These widespread haplotypes belong to common West Eurasian hgs H, T and W. The 16223-16292 (hg W) haplotype represents a basal haplotype within the phylogeny of hg W. Eight RoIA haplotypes were infrequent in contemporary populations and were classified as informative (Table S4). The relative frequencies of the shared informative and non-informative RoIA haplotypes were calculated in each of twelve present-day populations. Three modern populations or groups of populations (Lithuanians and Latvians, Poles, and Czechs and Slovaks) were found to contain significantly higher percentages (p<0.05) of shared informative haplotypes with the RoIA samples compared to other present-day populations (Figure 2, Table S4). Notably, modern Poles shared the highest number (nine) of informative mtDNA haplotypes with the RoIA individuals. The remaining three haplotypes had no match in the screened modern populations and were classified as unique (Table S3). These unique haplotypes belonged to mtDNA hgs N1a (16147A-16223-16248-16320-16355), with (likely) back mutation (C>T) at the position 16172, U5a (16256-16263-16270) and W (16192-16223-16292-16399).
Figure 2. Frequencies (%) of haplotypes shared between Roman Iron Age individuals and sampled modern Europeans.
Grey bar- all shared haplotypes; green – informative haplotypes, which not differ significantly from the average frequency (p<0.05); graphite – informative haplotypes, which differ significantly from the average frequency (p<0.05). Abbreviations for populations: Poles (PL), Belarusians (BY), Russians (European part) (RU), Czechs/Slovaks (CZ, SK), Bosnians/Slovenians/Croatians (BH, SLO, HR), Bulgarians (BG), Macedonians/Serbs (MA, SI), Ukrainians (UA), Swedes (SE), Germans (DE), Lithuanians/Latvians (LT, LV), Estonians (EE), Finns (FI).
Among distinct mtDNA haplotypes in the ME samples, eight belonging to hgs H, T, K, and J were present in high frequencies in modern-day populations and were hence classified as non-informative (Table S4). The informative group consisted of nine haplotypes occurring in low frequencies in six modern populations (Table S4); however, only three out of the 12 tested modern populations (Bulgarians, Poles and Belarusians) were found to share significantly higher percentage (p<0.05) of informative haplotypes with the ME samples (Figure 3). Three haplotypes, characteristic of hg X4 (16183-16189-16223-16258-16266-16274-16278), hg HV0 (16129-16145-16298) and hg HV (16189-16271) did not have matches to the studied modern populations and were thus considered as unique to the ME samples.
Figure 3. Frequencies (%) of haplotypes (matches) shared between Medieval Age individuals and sampled modern Europeans.
Grey bar – all shared haplotypes; green bar- informative haplotypes, which not differ significantly from the average frequency; graphite bar – informative haplotypes, which differ significantly from the average frequency (p<0.05). For population abbreviations, see Figure 2.
MtDNA haplotypes shared between RoIA and ME populations
Comparison of mtDNA compositions of the RoIA and ME populations revealed two shared haplotypes. The first one was assigned to hg H with a mutation at position 16362. The second one belonged to hg H5 with a mutation at the diagnostic position 16304. Of the analyzed individuals assigned to hg H5, three RoIA haplotypes: 16153-16304 (sample R3), 16304 (sample K1), 16129-16304 (sample G1) (Table 2), and one ME haplotype (16304 (sample OL1)) (Table 3), further belonged to subhaplogroup H5a1 with defining mutation at coding region position 15833.
Pairwise distances and MDS
Pairwise genetic distances were calculated in order to reconstruct the genetic relationship between the ancient and modern populations. Pairwise FST values showed non-significant differences between the RoIA and the ME samples (p>0.01) (Table S5). The RoIA samples differed significantly when compared to Neolithic individuals (LBK, Germany), Ukrainians, Belarusians, Latvians and Finns (p<0.01). The ME individuals showed no significant genetic differences to other populations used in the analysis (p>0.01), with the exception of the Neolithic (LBK) population and Finns (p<0.01) (Table S5). Correspondingly, in resulted MDS plot both ancient samples from this study (RoIA, ME) are mostly differentiated from northeast modern Europeans and Neolithic LBK sample, while being within the conglomerate of the remaining populations (Figure 4).
Figure 4. Multidimensional Scaling plot based on FST values calculated from mitochondrial haplogroup frequencies in sampled European populations.
Modern Slavic populations and other populations adjacent to Poles (black): for populations abbreviations, see Figure 2; ancient comparative populations (green): Danish (Iron Age) (DIA), Danish (Medieval) (DM), Neolithic (LBK, Germany) (NEO); present study populations (red): Roman Iron Age (RoIA), Medieval Age (ME).
This is the first large-scale study presenting mtDNA profiles for 43 individuals recovered from Roman Iron Age and Medieval Age burial sites from modern-day Poland. All mtDNA hgs identified in RoIA and ME populations are observed in most modern Slavic populations –. However, distinct haplotypes detected in the ancient samples and corresponding frequencies in modern populations have the potential to address the question of continuity of mitochondrial lineages over time in the area of present-day Poland.
Haplotype sharing analysis indicates that the RoIA individuals shared the highest number of informative mtDNA haplotypes with present day Poles (Table S4). Of particular interest are three RoIA samples assigned to subhaplogroup H5a1, which were recovered from the Kowalewko (sample K1), the Gąski, and the Rogowo (samples G1 and R3) burial sites (see Figure 1). Recent studies on mtDNA hg H5 have revealed that phylogenetically older subbranches, H5a3, H5a4 and H5e, are observed primarily in modern populations from southern Europe, while the younger ones, including H5a1 that was found among RoIA individuals in our study, date to around 4.000 years ago (kya) and are found predominantly among Slavic populations of Central and East Europe, including contemporary Poles . Notably, we also found one ME sample belonging to subhaplogroup H5a1 (sample OL1 in Table 3). The presence of subclusters of H5a1 in four ancient samples belonging to both the RoIA and the ME periods, and in contemporary Poles, indicates the genetic continuity of this maternal lineage in the territory of modern-day Poland from at least Roman Iron Age i.e. ∼ 2 kya. Age estimates for other subhaplogroups of mtDNA hg H5 (Ha5a2, H5e1a and H5u1) as well as for U5a2b1, U5a2a and U4a are similar as for subhaplogroup H5a1 (∼4 kya) in Central and Eastern Europe , .
The evolutionary age of H5 sub-branches (∼4 kya)  also approximates the age of N1a1a2 subclade found in the RoIA population (sample KA2) (Table 2). The coalescence age of N1a1a2 is around 3.4–4 kya, making this haplotype one of the youngest sub-branches within hg N . The N1a1a2 haplotype found in one RoIA individual was classified as unique because no exact match was found among the twelve comparative populations or groups of populations used in the haplotype sharing test. Notably, a similar N1a1a2 haplotype carrying an additional transition at position 16172 was found in a modern-day Polish individual . Taken together, the presence of mtDNA subhaplogroups N1a1a2 and H5a1 in both the ancient populations as well as in studied modern Poles suggests a genetic continuity of certain matrilineages in the territory of present-day Poland, at least from the Roman Iron Age.
Haplotype sharing test revealed that RoIA populations share significantly higher number of haplotypes not only with present-day Poles, and Czechs and Slovaks (Slavic populations), but also with Lithuanians and Latvians (Table S4). Similarity in genetic compositions among ancient populations across Central and north-eastern Europeans may reflect their shared deep matrilineal history being consistent with demographic history of Europe in general . Additionally, more recent events such as migration and admixture might have also contributed to the observed genetic similarity. For instance, archaeological evidence and historical records show extensive cultural connections between populations over wide areas of Europe during the Iron Age . In particular, the Roman Empire extended throughout the latitudinal breadth of Europe in the Iron Age. During this period, trade routes crossed the territories of modern Poland and eastern regions, inhabited by Baltic tribes  thus enabling contacts between RoIA populations (in particular those belonging to Wielbark culture) and neighboring Baltic groups, resulting not only in cultural exchange but, possibly, also in gene flow.
The ME populations share significantly higher percentages of haplotypes with modern Poles, Belarusians and Bulgarians (p<0.05) (Table S4). Shared haplotypes over time might reflect genetic continuity between populations from RoIA, ME and modern Poles as well as a matrilineal gene flow towards east and south of Europe at the onset of the Middle Ages, coinciding with movements of Slavic groups .
Inclusion of ancient Danish samples for comparative analysis was prompted by the fact that these groups as well as those from the other parts of Scandinavia are thought to have descended from Germanic tribes that occupied the territory of modern-day Poland during the Roman Iron Age . The MDS plot showed large genetic differentiation between Iron Age individuals from Denmark and RoIA individuals from Poland (Figure 4), although the FST values were not statistically significant (p>0.01). Both populations shared haplotypes belonging to subhaplogroups U3a and T2. Notably, we did not identify mtDNA hg I, occurring at high frequency (13%) in ancient Iron and Middle Age Danish populations , among any of the 43 successfully amplified ancient individuals from Poland.
We note that due to the uniparental inheritance of the mitochondrial genomes and relatively small sample sizes in our study, our observation of the various mtDNA hgs in the ancient populations, both the presence/absence and frequencies, are likely to be prone to stochastic events. Hence, future research including nuclear DNA markers are necessary to further explore the genetic connections between these ancient and modern populations, including the question of a common origin and the extent of admixture, if any, between them.
Results of our study indicate genetic continuity of mitochondrial lineages between ancient and modern populations in the territory of contemporary Poland. In particular, presence of sub-clades of hg H5a1 among both RoIA and ME ancient samples and present-day Poles, and the identification of N1a1a2 haplotype in RoIA and contemporary Poles is consistent with the idea of continuity of maternal lineages from at least Roman Iron Age in the region. Our data demonstrates that present-day Western Slavs, among analyzed Europeans, exhibit a mtDNA profile that is more similar to one found among ancient inhabitants of Central Europe. This observation appears to be in concordance with the autochthonous hypothesis. Studies on the genetic profiles of other ancient Slavic populations, especially employing nuclear markers, are necessary for further resolution of the complex origin of the Slavs.
Figure_S1.pdffigshare1 / 6 download
Alignments of cloned aDNA sequences analyzed in this study. The first lines report the revised Cambridge Reference Sequence (rCRS) with the numbering of the nucleotide positions.
Alignments of cloned aDNA sequences analyzed in this study. The first lines report the revised Cambridge Reference Sequence (rCRS) with the numbering of the nucleotide positions.
Detailed information about all samples used in the study involving their repository names, geographical location, context and references.
PCR primers used in the present study.
Modern and ancient populations used in the comparative analysis.
Haplotype sharing analysis.
FST distances by Slatkin’s based on haplogroup frequencies (P<0.01).
The authors are grateful to Anna Wrzesińska from Museum of the First Piasts at Lednica for her help in collecting samples.
Conceived and designed the experiments: AJ JP EW HM MD. Performed the experiments: AJ MR. Analyzed the data: AJ EM AK JZK. Contributed reagents/materials/analysis tools: AJ MD EW. Wrote the paper: AJ AK MR HM JZK.
- 1.Pohl W (2006) Beginnings of Slavs, a few historical observations In: Urbańczyk P, editor. Not-Slavs about beginnings of Slavs. Warsaw: Mała Biblioteka – Poznańskie Towarzystwo Przyjaciół Nauki. pp.11–25. (In Polish).
- 2.Buko A (2006) Polish archaeology of early Middle Ages: findings, hypotheses, interpretations. Warsaw: Trio (In Polish).
- 3.Alinei M (2003) Interdisciplinary and linguistic evidence for Paleolithic continuity of Indo-European, Uralic and Altaic populations in Eurasia with an excursus on Slavic ethnogenesis. Quadermi di Semantica 24: 187–216.
- 4.Curta F (2001) The Making of the Slavs. History and archaeology of the lower Danube Region, c. 500–700. Cambridge studies in medieval life and thought Fourth series, 52. Cambridge; New York: Cambridge University Press: xxv, 463 p.
- 5.Barford PM (2001) The early Slavs: culture and society in Early medieval Eastern Europe. London: British Museum.
- 6.Sussex R, Cubberley P (2006) The Slavic languages. Cambridge: Cambridge University Press.
- 7.Kostrzewski J, Chmielewski W, Jazdzewski K (1965) The prehistory of Poland. Wroclaw: The Ossolineum (In Polish).
- 8.Makiewicz T (2001) Problem of the cultural continuation between antiquity and the early Middle Ages in the light of new ceramic materials from Greater Poland. In: Mańczak W, editor. Homeland of the Slavs. Cracow (In Polish).
- 9.Godłowski K (2000) The first settlements of the Slavs. Cracow: Institute of Archaeology, Jagiellonian University (In Polish).
- 10.Dąbrowski R (2007) Human populations from the territories of Oder and Vistula rivers in Roman Iron Age and early Medieval. Series of Anthropology no 23, Adam Mickiewicz University in Poznan (In Polish).
- 11.Piontek J, Iwanek B, Segeda S (2008) Anthropology and the origin of Slavs. Poznan: Institute of Anthropology, Adam Mickiewicz University in Poznan (In Polish).
- 12.Malyarchuk BA, Grzybowski T, Derenko MV, Czarny J, Drobnič K, et al. (2003) Mitochondrial DNA variability in Bosnians and Slovenians. Ann Hum Genet 67: 412–425.
- 13.Grzybowski T, Malyarchuk BA, Derenko MV, Perkova MA, Bednarek J, et al. (2007) Complex interactions of the Eastern and Western Slavic populations with other European groups as revealed by mitochondrial DNA analysis. Forensic Sci Int: Genetic 1: 141–147.
- 14.Malyarchuk BA, Grzybowski T, Derenko M, Perkova M, Vanecek T, et al. (2008) Mitochondrial DNA phylogeny in Eastern and Western Slavs. Mol Biol Evol 25: 1651–1658.
- 15.Mielnik-Sikorska M, Daca P, Malyarchuk B, Derenko M, Skonieczna K, et al. (2013) The history of Slavs inferred form complete mitochondrial genome sequences. PLoS ONE 8: e54360.
- 16.Malmström H, Gilbert MTP, Thomas MG, Brandström M, Storå J, et al. (2009) Ancient DNA Reveals Lack of Continuity between Neolithic Hunter-Gatheres and Contemporary Scandinavians. Curr Biol 19: 1758–1762.
- 17.Deguilloux MF, Soler L, Pemonge MH, Scarre C, Joussaume R, et al. (2011) News from the west: Ancient DNA from a french megalithic burial chamber. Am J Phys Anthropol 144: 108–118.
- 18.Hervella M, Neskutus I, Alonso S, Fregel R, Alonso A, et al. (2012) Ancient DNA from Hunter-Gatherer and Farmer Groups from Northern Spain Supports a Random Dispersion Model for the Neolithic Expansion into Europe. PLoS ONE 7: e34417.
- 19.Brotherton P, Haak W, Templeton J, Brandt G, Soubrier J, et al. (2012) Neolithic mitochondrial haplogroup H genomes and the genetic origins of Europeans. Nat Commun 4: 1764.
- 20.Brandt G, Haak W, Adler CJ, Roth C, Szécsényi-Nagy A, et al. (2013) Ancient DNA Reveals Key Stages in the Formation of Central European Mitochondrial Genetic Diversity. Science 342: 257–261.
- 21.Bollongino R, Nehlich O, Richards MP, Orschiedt J, Thomas MG, et al. (2013) 2000 years of parallel societies in Stone Age Central Europe. Science 342: 479–81.
- 22.Rudbeck L, Gilbert MT, Willerslev E, Hansen AJ, Lynnerup N, et al. (2005) mtDNA analysis of human remains from an early Danish Christian cemetery. Am J Phys Anthropol 128: 424–429.
- 23.Melchior L, Gilbert MTP, Kivisild T, Lynnerup N, Dissing J (2008) Rare mtDNA Haplogroups and Genetic Differences in Rich and poor Danish Iron-Age villages. Am J Phys Anthropol 135: 206–215.
- 24.Melchior L, Kivisilid T, Lynnerup N, Dissing J (2008) Evidence of Authentic DNA from Danish Viking Age Skeletons Untouched by Humans for 1,000 Years. PLoS ONE 3: 1–8.
- 25.Melchior L, Lynnerup N, Siegismund HR, Kivisild T, Dissing J (2010) Genetic diversity among ancient Nordic populations. PLoS ONE 5: e11898.
- 26.Haak W, Balanovsky O, Sanchez JJ, Koshel S, Zaporozhchenko V, et al. (2010) Ancient DNA from European early Neolithic farmers reveals their Eastern Affinities. PLoS Biology 8: e1000536.
- 27.Wilde S, Timpson A, Kirsanow K, Kaiser E, Kayser M, et al. (2014) Direct evidence for positive selection of skin, hair, and eye pigmentation in Europeans during the last 5,000 y. PNAS111: 4832–4837.
- 28.Pilli E, Modi A, Serpico C, Achilli A, Lancioni H, et al. (2013) Monitoring DNA Contamination in Handled vs. Directly Excavated Ancient Human Skeletal Remains. PLoS ONE 8: e52524.
- 29.Yang DY, Eng D, Waye JS, Dudar JC, Saunders SR (1998) Technical note: Improved DNA extraction from ancient bones using silica based spin columns. Am J Phys Anthropol 105: 539–54.
- 30.Malmström H, Svensson EM, Gilbert MTP, Willerslev E, Gotherstrom A, et al. (2007) More on Contamination: The Use of Asymmetric Molecular Behavior to Identify Authentic Ancient Human DNA. Mol Biol Evol 24: 998–1004.
- 31.Svensson EM, Anderung C, Baubliene J, Persson P, Malmström H, et al. (2007) Tracing genetic change over time using nuclear SNPs in ancient and modern cattle. Anim Genet 38: 378–383.
- 32.Andrews RM, Kubacka I, Chinnery PF, Lightowlers RN, Turnbull DM, et al. (1999) Reanalysis and revision of the Cambridge reference sequence for human mitochondrial DNA. Nat Genet 23: 147.
- 33.Van Oven M, Kayser M (2009) Updated comprehensive phylogenetic tree of global human mitochondrial DNA variation. Hum Mutat 30: E386–E394.
- 34.R Development Core Team (2013) R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Available: http://www.R-project.org/. Accessed 5 March 2014.
- 35.Canty A (2013) Boot: Bootstrap Functions. R package version 1.9.13. Available: Boot: http://statwww.epfl.ch/davison/BMA/library.html. Accessed 5 March 2014.
- 36.Excoffier L, Laval LG, Schneider S (2005) Arlequin ver. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online 1: 47–50.
- 37.Stiller M, Green RE, Ronan M, Simons JF, Du L, et al. (2006) Patterns of nucleotide misincorporations during enzymatic amplification and direct large-scale sequencing of ancient DNA. Proc Natl Acad Sci USA 103: 13578–13584.
- 38.Gilbert MT, Willerslev E, Hansen AJ, Barnes I, Rudbeck L, et al. (2003) Distribution patterns of postmortem damage in human mitochondrial DNA. Am J Hum Genet 72: 32–47.
- 39.Hansen AJ, Willerslev E, Wiuf C, Mourier T, Arctander P (2001) Statistical evidence for miscoding lesions in ancient DNA templates. Mol Biol Evol 18: 262–265.
- 40.Kushniarevich A, Sivitskaya L, Danilenko N, Novogrodskii T, Tsybovsky I, et al. (2013) Uniparental Genetic Heritage of Belarusians: Encounter of Rare Middle Matrilineages with a Central European Mitochondrial DNA Pool. PLoS ONE 8: e66499.
- 41.Karachanak S, Carossa V, Nesheva D, Olivieri A, Pala M, et al. (2011) Bulgarians vs the other European populations: a mitochondrial DNA perspective. Int J Legal Med 126: 497–503.
- 42.Malyarchuk B, Derenko M, Grzybowski T, Perkova M, Rogalla U, et al. (2010) The peopling of Europe from the Mitochondrial Haplogroup U5 perspective. PLoS ONE 5: e10285.
- 43.Zgonjanin D, Veselinović I, Kubat M, Furac I, Antov M, et al. (2010) Sequence polymorphism of the mitochondrial DNA control region in the population of Vojvodina Povince, Serbia. Leg Med (Tokyo) 12: 104–107.
- 44.Malyarchuk BA, Perkova MA, Derenko MV, Vanecek T, Lazur J, et al. (2008) Mitochondrial DNA variability in Slovaks, with application to the Roma origin. Ann Hum Genet 72: 228–40.
- 45.Zimmermann B, Brandstätter A, Duftner N, Niederwieser D, Spiroski M, et al. (2007) Mitochondrial DNA control region population data from Macedonia. Forensic Sci Int: Genetic 1: e4–9.
- 46.Malyarchuk BA, Vanecek T, Perkova MA, Derenko MV, Sip M (2006) Mitochondrial DNA variability in the Czech population, with application to the ethnic history of Slavs. Hum Biol 78: 681–696.
- 47.Branicki W, Kalista K, Kupiec T, Wolańska-Nowak P, Zołedziewska M, et al. (2005) Distribution of mtDNA haplogroups in a population sample from Poland. J Forensic Sci 50: 732–733.
- 48.Vanecek T, Vorel F, Sip M (2004) Mitochondrial D-loop hypervariable regions: Czech population data. Int J of Legal Med 118: 14–18.
- 49.Malyarchuck, B A, Grzybowski T, Derenko MV, Czarny J, Wozniak, et al (2003) Mitochondrial DNA variability in Poles and Russians. Ann Hum Genet 66: 261–83.
- 50.Belyaeva O, Bermisheva M, Khrunin A, Slominsky P, Bebyakova N, et al. (2003) Mitochondrial DNA variations in Russian and Belorussian populations. Hum Biol 75: 647–660.
- 51.Malyarchuk BA, Derenko MV (2001) Mitochondrial DNA variability in Russians and Ukrainians: implication to the origin of the Eastern Slavs. Ann Hum Genet 65: 63–78.
- 52.Palanichamy MG, Zhang CL, Mitra B, Malyarchuk B, Derenko M, et al. (2010) Mitochondrial haplogroup N1a phylogeography with implication to the origin of European farmers. BMC Evol Biol 10: 304.
- 53.Grzybowski T, Malyarchuk BA, Derenko MV, Perkova MA, Bednarek J, et al. (2007) Complex interactions of the Eastern and Western Slavic populations with other European groups as revealed by mitochondrial DNA analysis. Forensic Sci Int Genet 1: 141–147.
- 54.Soares P, Achilli A, Semino O, Davies W, Macaulay V, et al. (2010) The archeogenetics of Europe. Curr Biol 20: R174–R183.
- 55.Gimbutas M (1963) The Balts. London: Thames and Hudson Press. 286 p.
- 56.Sedov VV (2002) Slaviane: Istoriko-arheologicheskoe issledovanie. Moskva, Yazyki slavianskoi kultury (in Russian).