Abstract The genetic composition of the medieval populations of Central Europe has been poorly investigated to date. In particular, the region of modern-day Slovakia is a blank spot in archaeogenetic research. This paper reports the study of mitochondrial DNA (mtDNA) in ancient samples from the 9th–12th centuries originating from the cemeteries discovered in Nitra-Šindolka and Čakajovce, located in western Slovakia (Central Europe). This geographical region is interesting to study because its medieval multi-ethnic population lived in the so-called contact zone of the territory of the Great Moravian and later Hungarian state formations. We described 16 different mtDNA haplotypes in 19 individuals, which belong to the most widespread European mtDNA haplogroups: H, J, T, U and R0. Using comparative statistical and population genetic analyses, we showed the differentiation of the European gene pool in the medieval period. We also demonstrated the heterogeneous genetic characteristics of the investigated population and its affinity to the populations of modern Europe.

Citation: Csákyová V, Szécsényi-Nagy A, Csősz A, Nagy M, Fusek G, Langó P, et al. (2016) Maternal Genetic Composition of a Medieval Population from a Hungarian-Slavic Contact Zone in Central Europe. PLoS ONE 11(3): e0151206. https://doi.org/10.1371/journal.pone.0151206 Editor: Michael Hofreiter, University of York, UNITED KINGDOM Received: October 8, 2015; Accepted: February 23, 2016; Published: March 10, 2016 Copyright: © 2016 Csákyová 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: Sequences were submitted to NCBI GenBank under the accession numbers KU739137 – KU739155. Funding: This research was realized in cooperation with the Constantine the Philosopher University in Nitra, the Institute of Archaeology, Research Centre for the Humanities, Hungarian Academy of Sciences, the Institute of Archaeology of the Slovak Academy of Sciences, and the University of Janos Selye in Komárno, in realizing the national grant of the Ministry of Education, Science, Research and Sport of the Slovak Republic and Slovak Academy of Sciences (Vedecká grantová agentúra Ministerstva školstva, vedy, výskumu a športu Slovenskej republiky a Slovenskej akadémie vied) VEGA 1/0897/12 „Archaeogenetic research of contact zone from the 10th century in Slovakia“ (MaB) and later supported also by the project of University Grant Agency of the Constantine the Philosopher University in Nitra UGA VIII/5/2015 (VC). The grants funded the research but had no role in study desing data collection and analysis. URLs: https://www.minedu.sk/rozpis-dotacii-na-nove-a-pokracujuce-projekty-vega-na-rok-2012/; https://www.ukf.sk/projekty/univerzitna-grantova-agentura. Competing interests: The authors have declared that no competing interests exist.

Introduction The territory of Central Europe–including the Carpathian Basin–was a place of great migration events of populations in the past. Hungarians arrived into Central Europe from across the Carpathian Mountains and settled in the Danubian Basin in 895–896 AD [1,2]. This region had been settled before the Hungarians’ arrival by Dacians, Celts, Romans, Sarmatians, Huns, German tribes (Goths, Gepids, Lombards), Avars and others, but the majority of the indigenous population was Slavic. The arrival of the Slavs and their settlement on the outskirts of the Carpathian Basin had already begun in the late migration period (6th century) during which they came into contact with Lombards and Gepids [3]. At present, vast territories of East-Central and South-Eastern Europe are inhabited by Slavic populations [4]. Three groups of present-day Slavs are identified based on 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) [5]. The modern Hungarian population, with a different cultural, linguistic and geographic origin, possesses a central location between these groups of Slavs. Archaeologists commonly associate certain grave goods, settlement patterns and burial customs of the cemeteries investigated at Nitra-Šindolka and Čakajovce as typically medieval plebeian mixed populations (so-called “Bijelo Brdo culture” [6,7]) of medieval Slavs and Hungarians (Magyars). The Nitra-Šindolka cemetery, dated to the 10th–11th centuries and located in western Slovakia (Fig 1), was excavated in 1985–1986 during the construction of the highway around the town of Nitra and consists of two parts [8]. The medieval Čakajovce cemetery originating in the 9th–12th centuries was excavated in 1974 and 1976–1986 [9]. PPT PowerPoint slide

PowerPoint slide PNG larger image

larger image TIFF original image Download: Fig 1. Location of medieval cemeteries from the contact zone of Central Europe. (1) Čakajovce (Slovakia). (2) Nitra-Šindolka (Slovakia). (3) Lébény (Hungary). (Illustration). https://doi.org/10.1371/journal.pone.0151206.g001 The observed archaeological material culture is open to widely different interpretations, and in particular it is unclear whether and to what extent the attributes (e.g. grave goods and burial traditions) are indicators of ethnic and social identity, and which of these materials are linked to population movements [10]. Ancient DNA provides direct genetic evidence for past demographic events. Mitochondrial DNA (mtDNA) is a convenient tool for studying the migration of populations in the past (e.g. [11–16]) since due to its exclusive maternal inheritance it does not recombine [17], so we can monitor the mutations which have accumulated over the centuries. Therefore, human mtDNA contains a molecular recording of the genealogical history and of the migrations of women who transmitted mitochondria through the generations [15]. This information obtained from parallel biological and genetic analysis, together with archaeological and historical information and interpretations, helps us to better understand the history of human populations, their origins, migration patterns and family relationships. The European medieval period is poorly investigated genetically; only a few hundred data have been published for the 6th–16th centuries from Hungary, Poland, Italy, Spain and Scandinavia [10,18–25]. Furthermore, this dataset encompasses a thousand years, and is characterized by migrations and admixture of populations, by the depopulating effects of the Black Death, as well as by formations of new medieval states and ethnicities. Most of the medieval populations were heterogeneous and shared common West Eurasian mtDNA haplogroups, independent of their geographical locations and origin [10,18–25]. Consequently, more genetic data is needed to understand migration and to reconstruct the genetic composition of populations living in Europe a thousand years ago. The aim of our study was to determine the genetic composition of maternal lineages in ancient populations from the medieval period of present-day Slovakia, and to characterize relationships to other medieval and modern human populations from Europe and Asia, using population genetics analysis. Genetic exploration of medieval populations in Central Europe is of paramount importance, because people from mixed populations of different ethnic origins lived contemporaneously in this area, making it difficult to distinguish them through archaeological records.

Results We obtained mtDNA sequences from 19 of 28 analyzed samples originating from the medieval Nitra-Šindolka and Čakajovce cemeteries (western Slovakia). We could not extract DNA from four samples, and we excluded a further three samples from Nitra-Šindolka and one sample from Čakajovce from the statistical analysis because of bad DNA preservation or ambiguous haplotype results. The HVR I sequences spanned the range of minimum np 16059–16400 and maximum np 16040–16421. The DNA of 19 individuals was extracted and amplified at least twice per individual from different skeletal elements, and the HVR I fragments and other coding region positions were reproduced in several PCR products per extracts. The results of these replicates along with the negative control results suggest the haplotypes are authentic. The medieval sequences encompass almost the entire range of Western Eurasian macro-haplogroups: H, J, T, R0 and U (S2 Table). However, we observed lower frequencies of haplogroup H and higher frequencies of haplogroups J and T in the medieval population compared with modern-day Europeans, e.g. present-day Slovaks. The haplogroup frequencies in the medieval population of Slovakia are the following: H 25%, J 20%, T 30% (T1 15% and T2 15%), R0 5% and U 20% (U4 5%, U5a 5% and U8 10%). Pairwise genetic distances were calculated between 11 medieval populations and between the same medieval populations and a set of modern-day populations from Eurasia. Of the medieval populations, the medieval Poles alone showed non-significant difference from the medieval population of present-day Slovakia (SVK_med) with F ST = 0.02702 and p = 0.08474 ± 0.0025 (S6 Table and Fig 2). Interestingly, pairwise F ST values of modern-day populations indicated a non-significant difference (p > 0.05) only between the medieval Slovakian population and the population from Iraq (F ST = 0.01372, p = 0.08105 ± 0.0081). The present-day Slovaks have a slightly larger genetic distance from medieval Slovaks than other Eastern and Central European modern populations: F ST = 0.05492, p = 0.000 ± 0.000 (S5 Table). The MDS plot (stress value = 0.1130, S1 Fig) based on linearized Slatkin F ST values contains aggregation of most of the medieval and modern European populations along coordinates 1 and 2; only the Asian populations of Pakistan, Uzbekistan, Kazakhstan and Kyrgyzstan are separated, forming an eastern cluster along coordinate 1. The ancient Xiongnu population from Mongolia is the farthest population along this coordinate. The medieval population from Slovakia is situated along coordinates 1 and 3 among modern European populations, but is differentiated from them along coordinate 2. In contrast, the other medieval European populations are close to modern populations along coordinates 1 and 2, and differentiated along coordinate 3, which implies their stronger affinity to modern European populations compared with the Slovakian dataset. PCA was calculated based on mtDNA haplogroup frequencies of medieval and modern populations. PCA of medieval populations shows a clustering of populations from Slovakia, Lombards (Hungary) and Poland along the PC1 and PC2 components (39.6% of the variance is displayed) (Fig 3A). The medieval populations and the Vikings of North Europe are clustered together, as well as the medieval populations from South Europe. The farthest populations along the PC1 component are the medieval Mongolians and Cumanians from Hungary, whereas the medieval Hungarians are the utmost population along PC2. Hierarchical clustering shows a very similar diagram to the PCA, except for one phenomenon, namely that the medieval population from Hungary clusters with Lombards from Hungary and with medieval populations from Poland and Slovakia, which we can explain with the interaction of all PCs in this analysis (Fig 3B). The PCA plots of modern Eurasian populations indicate alignment of the European populations (except Saami) including our investigated medieval population on the PC1 and PC2 components (variance: PC1 = 24.45%, PC2 = 16.60%). However, the different frequencies of haplogroups J, T1, T2, U4, U5a and U8 result in differentiation from all present-day populations along the third component (variance PC3 = 10.71%) (Figs 4 and S2). The SHA shows that the medieval population from Slovakia shares the majority of haplotypes with medieval populations from Spain and Italy, and that the relatively short genetic distance from the medieval Poles does not mean significant lineage sharing at the HVR I level (S3 Fig).

Discussion In this study we investigated a population living in the Hungarian-Slavic contact zone of the Carpathian Basin in the 9th–12th centuries AD. Using well established protocols for aDNA analysis, we have obtained the first mtDNA dataset (to the best of our knowledge) concerning the genetic variation of the HVR I region and six nucleotide polymorphisms in the coding region of mtDNA. These markers enabled us to assign the mtDNA haplogroups of the samples. The representative haplogroups of the investigated medieval population belong to common modern European (West Eurasian) haplogroups. However, their frequencies are quite different from present-day populations. Contrarily, the detected frequency of haplogroup H in the medieval population from Slovakia was lower (25%) than it is in European modern-day populations (including the Slovaks), where this most frequent haplogroup has a frequency of 45.4% [14,15,40]. The most ancient European haplogroup U [41,42] occurs with a similar frequency in the medieval population (20%) and in modern Slovaks (15%) [40]. Haplogroups J and T occurred in the 9–12th centuries two or three times more frequently than in present-day Europe [32,40,43]. The same phenomenon is perceptible in the Lombard population from the western Hungarian site Szólád with regard to the haplogroups H and J, which have either lower (in case of H) or higher (in case of J) frequency compared to other modern European populations [19]. This similarity results in the clustering of the medieval populations of Slovakia and Lombards from Hungary in the PCA (Fig 3). The Lombard people were of diverse origin and their cemetery was in use only for a short period of time in the 6th century. The research team of Alt et al. [19] suggested that people moved as a group whilst also integrating unaffiliated individuals. Based on the observed diversity of the mtDNA haplogroups, they came to the conclusion that this population includes numerous characteristic Central European lineages as well as a few rare haplogroups. Moreover, our PCA, MDS and SHA presented two populations sharing the same cultural label “Lombards” from different locations (Italy and Hungary). Interestingly, they are relatively far from each other in genetic context, in spite of the historical records of Lombard migration across Europe. This probably follows from the characteristic Central European mtDNA lineages of the Lombards from Hungary [19] in contrast to the Lombards from Italy, who had correlation with modern-day populations inhabiting the same Italian geographical area [10]. Furthermore, the PCA showed that the medieval population of Poland clusters with the medieval population of Slovakia and the Lombards from Hungary (Fig 3), which we explain with similar mtDNA haplogroup composition of medieval Poles to the modern Slavic population spread over a major area of Central and Eastern Europe [20]. The medieval Slovakian samples show closer ties with Polish samples of Slavic origin than with 10th century Hungarian samples from Hungary [18]. However, the ongoing analysis of further 10th century samples from the region will probably result in a more accurate estimation of the affinity of the medieval populations. The medieval population under study shares mtDNA haplogroups with common European populations, which is represented in the PCA of modern Eurasian populations. The connection is especially strong to modern Central and South Europeans (Fig 4). The MDS plot showed clustering of the medieval population from Slovakia and other European populations especially along the third coordinate, as well as population genetic differentiation from Asia along the first and second coordinates, although some F ST values were not statistically significant (p > 0.05) (S5 Table and S1 Fig). From the F ST calculation a different result was obtained by the SHA, where a high frequency of shared haplotypes was shown between medieval populations from Slovakia, Spain and Italy. Summarizing, we consider the studied medieval population of Slovakia to be a mixture of groups of different genetic origins. However, more genetic data from Central Europe would be needed to investigate this epoch in detail.

Conclusion The comparative analysis of the medieval sample sets suggests that the present-day European gene pool could have been formed in the medieval era. This is supported by the presence of a majority of haplogroups and haplotypes in the medieval period which still exist in modern-day Europe. Genetic maternal lineages of the medieval population of present-day Slovakia were diverse and rather similar to medieval Lombards from Hungary and to medieval Slavs from Poland, as well as to modern Europeans, which supports the historical and archaeological standpoint of mixed populations of medieval Slavs and Magyars in the cemeteries investigated. However, the results could have been biased due to the small sample size, and genetic drift could also unnoticeably influence our results. Furthermore, we assessed population interactions and migrations only along the maternal lineages; the paternal side of the medieval population history remained unexamined in this study. The collection of archaeological relics from genetically examined Slovakian cemeteries overlaps the image drawn by population genetics analysis. The finds of the 10th century cemeteries originated in the material culture of the preceding era, implying that Hungarian groups arriving in the territory blended in with the local population, genetically as well as culturally; this interaction produced the medieval population of the area. Undoubtedly, further research, including the analysis of the whole mitochondrial genome and nuclear DNA markers, is necessary in the investigation of the genetic connections between medieval and modern populations and in the clarification of the genetic origin of Slavs and Hungarians. Nevertheless, the study presented here of the first successful description of medieval mtDNA variability from Slovakia contributes an important dataset to help reveal the genetic diversity of medieval Europe.

Acknowledgments The authors are grateful to Guido Brandt from the Max Planck Institute for the Science of Human History, Jena, Germany for collecting the modern reference data and to Milan Hanuliak from the Institute of Archaeology of the Slovak Academy of Sciences, Nitra, Slovakia for the mediation of materials from the Čakajovce cemetery. They furthermore thank Alena Šefčáková from the Slovak National Museum, Bratislava, Slovakia for helping and lending the ancient bone material from the depository of the museum.

Author Contributions Conceived and designed the experiments: M. Bauerová MN. Performed the experiments: VC AC. Analyzed the data: VC ASN AC. Contributed reagents/materials/analysis tools: M. Bauer PM MN. Wrote the paper: VC ASN. Selected and collected the samples: MN BGM GF PL VC. Collected the reference data: ASN VC. Performed the population genetic statistical analyses: ASN VC. Discussed the results and commented on the manuscript: M. Bauerová MN BGM GF PL VC M. Bauer PM ASN AC.