Cranial deformation and genetic diversity in three adolescent male individuals from the Great Migration Period from Osijek, eastern Croatia

Daniel Fernandes; Kendra Sirak; Olivia Cheronet; Rachel Howcroft; Mislav Čavka; Dženi Los; Josip Burmaz; Ron Pinhasi; Mario Novak

doi:10.1371/journal.pone.0216366

Abstract

Three individuals dating to the Great Migration Period (5^th century CE) were discovered in a pit at the Hermanov vinograd site in Osijek, Croatia. We were inspired to study these individuals based on their unusual burial context as well as the identification of two different types of artificial cranial deformation in two of the individuals. We combine bioarchaeological analysis with radiographic imaging, stable isotopes analysis, and ancient DNA to analyze their dietary patterns, molecular sex, and genetic affinities in the context of the archaeological data and their bioarchaeological attributes. While all three individuals were adolescent males with skeletal evidence of severe malnutrition and similar diets, the most striking observation is that they had major differences in their genetic ancestry. Results of the genetic analyses of the nuclear ancient DNA data for these individuals indicate that the individual without artificial cranial deformation shows broadly West Eurasian associated-ancestry, the individual with tabular oblique-type has East Asian ancestry and the third individual with circular erect-type has Near Eastern associated-ancestry. Based on these results, we speculate that artificial cranial deformation type may have been a visual indicator membership in a specific cultural group, and that these groups were interacting intimately on the Pannonian Plain during the Migration Period.

Citation: Fernandes D, Sirak K, Cheronet O, Howcroft R, Čavka M, Los D, et al. (2019) Cranial deformation and genetic diversity in three adolescent male individuals from the Great Migration Period from Osijek, eastern Croatia. PLoS ONE 14(8): e0216366. https://doi.org/10.1371/journal.pone.0216366

Editor: Peter F. Biehl, University at Buffalo - The State University of New York, UNITED STATES

Received: January 11, 2019; Accepted: July 16, 2019; Published: August 21, 2019

Copyright: © 2019 Fernandes 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: All relevant data are within the manuscript and its Supporting Information files. The raw data has been deposited for public access in the European Nucleotide Database with the accession number PRJEB33257.

Funding: DF was supported by an Irish Research Council Post-Graduate grant GOIPG/2013/36. RP was supported by the European Research Council Starting Grant ERC-2010-StG 263441. MN was supported by the Irish Research Council Postdoctoral Fellowship grant GOIPD/2013/1, and by the Croatian Science Fund grant HRZZ IP-2016-06-1450. http://research.ie/; https://erc.europa.eu/; http://www.hrzz.hr/. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Two of the co-authors (Dženi Los and Josip Burmaz) are employed by a commercial company Kaducej Ltd. The funders provided support in the form of salaries for authors DF and MN, but did not have any additional role in the 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. Two of the co-authors (Dženi Los and Josip Burmaz) are employed by a commercial company Kaducej Ltd. This commercial affiliation does not alter our adherence to PLOS ONE policies on sharing data and materials.

Introduction

Humans exhibit similarities in cultural practices that span across space and throughout time. One such practice is Artificial Cranial Deformation (ACD), a widespread cultural phenomenon performed to denote group and/or individual identity; that is, to distinguish particular people from others or provide visible evidence of status, nobility, or affiliation to a certain class or group [1–4]. This practice of permanently shaping the cranial vault has been documented by anthropologists to have occurred on every inhabited continent, possibly dating back to the Late Paleolithic and practiced at least up to the 20^th century [5–10].

ACD is an irreversible and deliberate act performed by adults on infants, reflecting an investment of time and an effort to achieve an ascribed identity [1, 11]. Unlike other types of body modification, ACD is not associated with any other known rite of passage, making it an indelible symbol of social identity [11]. As such, it serves to foster intra-group solidarity and acts as a visual manifestation of inter-group cultural differences [11]. Taking advantage of the plasticity of the infant skull, ACD is carried out by application of constant external pressure on the head from the first days of life, often using boards, pads, bags of earth or clay, or specially-made headdresses, to achieve a desired shape [5, 12]. Because the cranium is a functional matrix comprising numerous interacting units, artificially-inflicted restriction of cranial growth along a particular axis will result in compensatory growth along a less restricted or unrestricted axis, ultimately causing an alteration of cranial shape [13–15]. Growth restriction of the cranial vault may also affect growth in adjacent areas, such as the cranial base or face [13, 14, 16–18]. While cranial morphology is affected, the brain enclosed inside the artificially deformed cranium achieves a volume comparable to that enclosed inside an unmodified cranium [15].

Historically, anthropologists have employed various classification systems for artificially-induced cranial deformations based on external vault morphology, though the definitions of deformation types are not consistent across the literature [14]. Some studies establish unique characteristic categories based on specific observations of the available samples, employing classification types such as “posterior flattening,” “bilobed,” and “circumferential” [19], or “occipital,” “coronal-occipital,” “occipital-parietal,” and “circular” [20]. Others emphasize more general morphological types, broadly comparing the anteroposteriorly-deformed type compared to the circumferentially-deformed type [5]. Both of these broad groups contain subgroups of “erecta” and “obliqua” types referring to the orientation of the deforming pressure, and therefore to the position of the occipital plane in relation to the facial plane. Specifically, the erect type deformation causes the occipital bone to take a direction perpendicular to the Frankfort plane, while the oblique type deformation causes the occipital to take a direction oblique in relation to the Frankfort plane [14, 21–23].

As important as classifying ACD by deformation type is understanding its origin and meaning in an archaeological context. The practice of ACD has ancient roots in Eurasia, appearing to arise independently at multiple points in time [1]. Archaeological evidence suggests that ACD on a broader scale was introduced to Central Europe with the migration of the Huns [2, 24]. During the Great Migration Period (5^th/6^th century CE), this custom became widespread in the region of the Pannonian Plain (modern-day Austria, Croatia, Hungary, Romania, Serbia, Slovakia, and Slovenia), and was particularly popular among the Sarmatian, Alan, Gothic, Longobard, Gepidic, and Hun populations [4, 25]. The ACD detected in these populations has been categorized within the “Danube Basin Group” of Eurasian ACD, part of the six-group system based on archaeological, chronological, and geographic contexts by Nemeskéri [26]. Within this broad group, mainly represented by Hun and Germanic crania dated to the 5^th/6^th centuries CE, the extent of deformation varies from slightly to heavily deformed and includes diverse types of modifications (i.e., different cranial shapes) [4]. Though detected across Pannonian Plain populations, the meaning of this cultural practice to these groups has not yet been clearly elucidated. While some earlier studies suggested that the deformed skull shape was initially a visible indication of nobility, and later became somewhat of a “trend” [27, 28], Molnár and colleagues more recently proposed that the custom of ACD in this region might be an indication of the social status [4].

Rescue archaeological excavations recently conducted at the well-known Hermanov vinograd site in Osijek, eastern Croatia, revealed a single pit, without any traces of the contemporary settlement, which contained a large number of commingled and fragmented human and animal skeletal remains together with the pottery fragments dated to the 5^th century CE, an era known as the Great Migration Period. Three partial human skeletons were discovered in this pit. Although the human remains belonging to the three individuals were haphazardly positioned without any consistencies in positioning and were only partially preserved, they all had complete crania. Analysis of the crania revealed that two of the individuals exhibited evidence of ACD, the first evidence of ACD detected in Croatia before the Avar Period. The discovery of this unusual archaeological feature and the unexpected characteristics of the remains found within it raised a number of questions about these individuals, primarily regarding their geographic origins and the potential meaning of ACD in this particular context.

Historic sources indicate that during the 5^th and 6^th centuries CE, the region around modern-day city of Osijek (the Roman province Pannonia Secunda) was under the Hunnish rule (until 454 CE), then under the Ostrogothic rule (till 473 CE), and finally under the Gepid rule, with some shorter interruptions, until 567 CE when they were replaced by the Avars [29]. The available archaeological data also support this hypothesis [30]. While archaeological information from Osijek does not clearly assign the individuals discussed in this study as belonging to one population or another, it is possible that these individuals were Huns, Ostrogoths or Gepids.

Recent studies of ACD have offered unprecedented insight into the origins of ACD, the biological and social correlates of this phenomenon, the long-term effects of this procedure on skeletal development and morphology, and details of ancient individuals who exhibit ACD, including their molecular sex and population genetic affinity. For example, a study by Cottin et al. aimed to reinvestigate the presumed effect of ACD on basicranial and masticatory elements by applying a 3D geometric morphometric approach to CT scans [31] while a study by Veeramah and colleagues [32] generated genomic data from 41 individuals with and without ACD that provided clear evidence of female-biased long-distance migration during the Migration Period in Europe. The mentioned studies led us to focus our investigation on the three individuals from Osijek by taking an approach which combines traditional bioarchaeological analysis with radiographic (CT) imaging, stable isotope analysis, and ancient DNA analysis to explore patterns associated with ACD at Osijek.

In this context, our present analysis of three archaeological individuals from the Hermanov vinograd site, two of whom exhibit ACD, has three aims: (i) determine the type of ACD in each individual; (ii) assess similarities and differences between the analyzed individuals in terms of age at death, sex, and overall health; and (iii) explore genetic affinities and components of ancestry in the studied individuals by means of ancient DNA analysis.

Results

Though the skeletons were only partially preserved, the degree of preservation of the elements that were recovered during excavation was excellent. In the case of SU 259, an almost complete cranium with mandible, both scapulae, and most long bones (minus the left radius and ulna, and both fibulae) together with most of the ribs were present; SU 261 was represented by a complete cranium without the mandible, and coupled by the first ribs; and SU 750 was represented by a complete cranium with mandible, right scapula, five left and five right ribs, one thoracic and two lumbar vertebrae, the sacrum, the right innominate, and both femora. All three individuals were adolescents: estimated age at death for SU 259 was 14–16 years, for SU 261 was 12–16 years, and for SU 750 was 12–14 years. The skeletons exhibit several pathological changes: SU 259 had healed sub-periosteal new bone formation on both tibiae, SU 261 had healed porotic hyperostosis on both parietal bones, and SU 750 had healed porotic hyperostosis on both parietal bones and linear enamel hypoplasias on all anterior teeth. No traces of injuries (ante- or peri-mortem) or post-mortem cut-marks were recorded on any of the analyzed skeletons. The skeletal remains were relatively robust, and cranial morphology suggested that all three individuals were males. However, due to widely-known issues in sexing non-adult skeletons from both forensic and archaeological contexts based on metric and non-metric traits, we looked to confirm our sex assessment using results obtained by ancient DNA (aDNA) analysis.

The cranium belonging to the individual SU 261 did not exhibit any traces of ACD, while morphological alternations of the crania consistent with ACD were recorded for SU 259 and SU 750. The cranial length of SU 259 is considerably increased and the general shape axis is dislocated posteriorly above the Frankfort horizontal plane (Fig 1A–1C). The cranium exhibits a depressed and strongly inclined frontal bone with a significant thickening of the posterior part visible on the anterior-posterior sagittal CT scan (Fig 1D); no traces of flattening are present on the occipital bone. The supraorbital ridges are almost non-existent, and all cranial sutures are completely fused as a result of the artificial modification. These observations suggest that SU 259 exhibits the tabular oblique type of ACD. It is widely accepted that the modifications resulting in tabular oblique deformation are carried out by the application of anterior-posterior compression using hard and rigid tools [3, 4]. In the case from Osijek, the modification was likely performed using a rigid tool in the frontal region only, as no traces of flattening are present on the occipital bone.

Download:

Fig 1.

A-C) CT reconstruction showing the artificially deformed cranium belonging to individual SU 259, frontal and lateral views. The shape axis is dislocated posteriorly above the Frankfort horizontal plane while the cranium exhibits a depressed and strongly inclined frontal bone strongly indicating tabular oblique type of deformation. D) X-ray of the same cranium (lateral view) showing a significant thickening of the posterior part.

https://doi.org/10.1371/journal.pone.0216366.g001

In contrast, the cranium belonging to the individual SU 750 exhibits a pronounced flattening of the frontal bone resulting in a remarkable growth of cranium height. Approximately 30 mm anterior of bregma, the frontal bone is inclined at a 45-degree angle and is completely flattened up to about 60 mm superior of the supraorbital ridges (Fig 2A–2D). However, no similar changes were recorded on the occipital bone. These observations suggest that SU 750 exhibits the circular erect type of ACD. To achieve the circular erect deformation observed in SU 750, more flexible tools such as bandages, bands, tapes, and/or headdress were likely used [3, 4]. While it is possible that bandages were used to shape the cranium of SU 750, there are no macroscopically visible grooves caused by bandaging.

Download:

Fig 2.

A-C) CT reconstruction showing the artificially deformed cranium belonging to individual SU 751, frontal and lateral views. A pronounced flattening of the frontal bone resulting in a remarkable growth of cranium height suggests the circular erect type of ACD. D) The flattening and significant thickening of the frontal bone anterior of bregma is visible on the X-ray image (lateral view).

https://doi.org/10.1371/journal.pone.0216366.g002

Analysis of carbon and nitrogen stable isotopes was performed on rib fragments belonging to all three individuals. While the bone sample from the individual SU 259 did not provide enough collagen to be analyzed, the samples from SU 750 and SU 261 did produce enough collagen and were analyzed as having C:N ratios of 3.2 (SU 750) and 3.3 (SU 261), indicating good protein preservation. Both had high δ¹³C values (SU 261 −17.0‰ and SU 750 −15.5‰) and very similar δ¹⁵N values (SU 261 9.5‰ and SU 750 9.9‰), suggesting a mixed terrestrial C3/C4 diet with a heavy reliance on resources such as millet, and a relatively low animal protein intake.

Authentic aDNA was recovered for all individuals (details presented in Table 1). Because only low-coverage shotgun sequencing was performed, approximately 500,000 human reads were identified for each individual, corresponding to endogenous DNA contents between 66–68%. The DNA sequences were examined for authenticity using deamination frequencies at the end of the DNA molecule and average sequence length. Deamination frequencies ranged between 15–17% on the 5’-end and 10–13% on the 3’-end of the DNA molecules, while sequence lengths averaged between 50–53bp with a standard deviation of 12–13bp (Table 1). Both of these metrics match the expected patterns for ancient DNA, which is predicted to exhibit a high frequency of DNA deamination at the terminal ends of the molecules and short sequence lengths due to chemical damage to DNA over time. Negligible human reads were detected in any of the blanks suggesting extremely low levels of cross-contamination. Between 4,086 and 5,838 SNPs were recovered for the three individuals (Table 1). The raw data has been deposited for public access in the European Nucleotide Database with the accession number PRJEB33257. Analysis of genetic sex was performed using a method that determines sex by considering the ratio of sequences aligning to the X and Y chromosomes [33]. Results were consistent with the assessment made using skeletal morphology, finding that all three individuals were males (Table 1).

Download:

Table 1. Sequencing results and sample details.

https://doi.org/10.1371/journal.pone.0216366.t001

The low number of SNPs used allowed us to associate individuals only with broad geographical areas when we performed PCA; however, PCA results clearly suggest starkly different population affinities among the individuals. SU 259 projects near to modern East Asian populations (the closest being Cambodian and Tu), while SU 261 overlaps with Near/Middle Eastern populations (represented, for example, by Palestinian, Syrian, and Lebanese), and SU 750 projects between populations from the Caucasus, Europe, and the Near/Middle East (Fig 3A). Though it lacks the power to elucidate detailed relationships between the individuals under study and with present-day populations, PCA provides a qualitative indication that, at this resolution, is able to distinguish populations in sub-continental groupings, demonstrating that these three individuals do not share close genetic affinities with each other.

Download:

Fig 3.

A) Principal component analysis of ancient samples (wheels and diamonds) projected onto variation of modern Eurasian and North-African populations. B) Model-based ancestry clustering with ADMIXTURE using K = 5 as the number of ancestral populations. Most informative modern populations where individual components were maximized are plotted.

https://doi.org/10.1371/journal.pone.0216366.g003

With the intent of confirming that these results were not an artefact of low-resolution analysis, a secondary PCA was generated by reducing the pool of SNPs of the modern populations and the Early Medieval ancient individuals from Veeramah and colleagues [32] to the combined SNPs of our three individuals (totaling 15,698), and again projecting all ancient samples on to the PCs calculated using this dataset of reduced modern populations (S1 Fig). The distribution of the modern populations and the projected position of the three Osijek individuals were virtually unchanged, while the cluster of previously published Bavarian individuals expanded mostly on the PC1 axis, but while maintaining their West Eurasian affinities. This demonstrated that reducing the number of SNPs did not affect the sub-continental population groupings, and therefore suggests no significant analytical bias for these categorizations. Nevertheless, a conservative approach to this analysis should be taken, which is to consider both SU 261 and SU 750 as part of the West Eurasian cluster, and SU 259 part of the East Eurasian cluster.

The results from ADMIXTURE model-based clustering analysis are consistent with the differences in genetic affinities first detected using PCA (Fig 3A and 3B). Due to the low number of SNPs used in this analysis (ca. 5,000 SNPs per individual), only major ancestry components can be considered, as the small amounts of ancestry seen could be the result of statistical noise. Since cross-validation errors plateaued from K = 5 to K = 7, the replicate chosen was the one that maximized the main components seen on our three individuals, in this case at K = 5. ADMIXTURE results show that SU 259 primarily exhibits components of ancestry related to modern populations from East and Northeast Asia (including Ami, Mongola and Nganasan), consistent with the PCA results that suggested this individual shared the most affinity with populations from these regions (Fig 3B). The results of outgroup-f₃ tests of the form f₃(Yoruba, SU 259; East Asians/West Eurasians) confirm these affinities by showing that SU 259 shares the most alleles with East Asian populations such as Hezhen or Ulchi (Fig 4). Even when considering 2 standard deviations, the f₃ estimates for West Eurasian populations do not overlap those for East Asia, providing further support for the Asian origin of SU 259. In contrast, individuals SU 261 and SU 750 do not exhibit significant amounts of East Asian ancestry. The main difference between these latter two individuals is the presence of a large component seen maximized in populations from Southeastern Europe, the Near/Middle East, and North Africa (including Algerian, Druze or Assyrian) in SU261, consistent with the PCA results that suggest this individual clusters closer to Near/Middle Eastern groups than to European groups (Fig 3A). In contrast, the most prominent component of SU 750 is one that is related to a component maximized in present-day continental European populations (including Lithuanian, Basque, or Hungarian). These latter two individuals show ancestries similar to the ones seen in the Early Medieval Bavarians from Veeramah and colleagues [32] (Fig 3B), which also included the presence of individuals of Southern European/Near Eastern ancestries.

Download:

Fig 4. Outgroup-f3 results for individual SU 259 with thick and thin bars representing 1 and 2 standard deviations, respectively.

The results show this individual shares higher affinities with East Asian than with West Eurasian, even when the standard deviations are considered.

https://doi.org/10.1371/journal.pone.0216366.g004

Discussion

Although the purpose of the pit excavated at Osijek—Hermanov vinograd remains unclear, we were able to successfully retrieve valuable biological information from the individuals whose remains were recovered from this pit.

As previously mentioned, all individuals shared notable similarities in terms of sex, age and general health: (i) all were males, (ii) all were adolescents between 12 and 16 years old at the time of death, and (iii) all displayed evidence of skeletal pathologies, including porotic hyperostosis, linear enamel hypoplasia and new sub-periosteal bone formation, though all of the lesions seemed to be inactive at the time of death.

The fact that all three studied individuals from Osijek were all adolescent males of almost the same age (between 12 and 16 years) may be suggestive of a cultural, regional, or temporal style of interment, though more data will be required to test this hypothesis. Though they appear to have been placed in the pit without clear and consistent positioning, it is possible that they were deposited with a certain intent rather than randomly disposed, as their remains were accompanied by faunal remains and pottery fragments. This style of interment and the characteristics of the remains recovered have also been seen in another Migration Period context; specifically, the sex and age of the Osijek individuals resemble the sex and age at death of the individual recovered from the Migration Period pit at the Ernei—Carieră site in Romania [34]. This feature in certain aspects resembles the case from Osijek (fragmented skeleton of a young male without any cut-marks in a pit accompanied by animal bones and pottery dated to the 4^th/5^th c. CE); however, it also displays some differences such as the presence of numerous metal objects (weapons, tools and jewelry) and traces of burning [34]. Crișan and Lăzărescu characterized the feature from Ernei—Carieră as a ritual site “most probably involving beliefs regarding agricultural activities” [34]. Based upon the present state of research, it would be far-fetched to assume that the young males from Osijek and Ernei were human sacrifices; however, one has to bear in mind that certain cultures such as the ancient Moche of Peru preferred young adult males as their victims [35], and therefore, that this is a valid hypothesis that will benefit from future research.

The other striking similarity between the young males from Osijek is the fact that all three experienced ill-health during their childhood that was severe and lasting enough to be observed on their skeletal remains. Importantly, they all managed to survive and recover from physiological insults occurring during the early stages of their lives based on the healed state of these pathologies. Individuals SU 261 (no ACD) and SU 750 (ACD) display healed porotic hyperostosis on the cranial vault, a condition recognized by the occurrence of porosity on the outer table of the cranial vault (most often parietal bones and the occipital bone) and associated with acquired or genetic anemia caused by inadequate nutrition, metabolic or blood disorders, infectious disease, parasitism and weanling diarrhea [36–38]. SU 259 (ACD) displays healed sub-periosteal new bone formation on both tibiae, a woven bone formation that is macroscopically recognized as osseous plaques with demarcated margins or irregular elevations of bone surfaces [39]. The occurrence of this pathological change is usually associated with ‘non-specific infections’, but also with the conditions like birth trauma, metabolic disorders, hypervitaminosis A, leukemia, and infantile cortical hyperostosis [40]. Finally, SU 750 shows evidence of linear enamel hypoplasia which is an indicator of subadult stress usually recognized as one or more horizontal lines on the labial tooth surface. Due to the fact that these defects cannot be remodeled after the formation of enamel makes them excellent indicators of subadult stress in the first seven years of life [41, 42]. The most common explanations for the occurrence of linear enamel hypoplasia in the past are childhood disease and malnutrition [39, 43], poor sanitary conditions [44], and the effects of weaning [44, 45].

Based on the pathology data presented above, it is clear that the three individuals in question suffered from a number of long-term conditions, mostly associated with nutritional deficiency, during their childhood. According to the contemporaneous historic sources, such as Procopius and Enodius, frequent episodes of hunger and starvation were a common occurrence during the Migration Period in southern Pannonia [46]. Extended and severe periods of hunger were likely to have had a significant negative impact on the general health of the majority of population, as indicated by the results of conventional bioarchaeological studies conducted on skeletal remains recovered from several Late Roman/Migration Period sites in the region [47–49]. More specifically, these studies showed that at some sites, such as Zmajevac [47] and Štrbinci [48], approximately half of the analyzed skeletons exhibit evidence of linear enamel hypoplasia, cribra orbitalia/porotic hyperostosis and sub-periosteal new bone formation. It is therefore unlikely that these three individuals from Osijek were intentionally starved; instead, it is probable that other members of the community would have exhibited these same skeletal insults that suggest the embodiment of long-term stress.

Beside similarities in skeletal evidence of pathologies, similarities in diet were also identified; specifically, carbon and nitrogen stable isotopes analysis showed that at least two of the individuals (SU 261 and SU 750) had a mixed terrestrial C3/C4 diet with a heavy reliance on resources such as millet. This is almost identical to the results obtained by Hakenbeck and colleagues [50] in their study of dietary habits of the historically documented Huns and other nomadic groups from five Hungarian sites dated to the 5^th century CE. Similar results were obtained by the authors who studied diet of the somewhat later Avar populations from the region; for example, Noche-Dowdy [51], and Vidal Ronchas and colleagues [52] recorded almost identical values in several Avar sites from Hungary and Croatia dated to the 6^th-9^th c. CE indicating a mixed diet containing C3 and C4 resources with a relatively low proportion of animal protein. Taken together, these studies suggest that C4 resources such as millet played the important role in everyday diet of various populations during the Migration Period in the region. This assumption is additionally strengthened by the presence of broomcorn millet in the nearby Bosnia and Herzegovina during this period [53], but also at the Avar site of Nuštar in eastern Croatia [54].

However, while they appear similar in terms of pathology and diet, two notable and potentially interrelated differences between these three individuals deserve further discussion and attempts at interpretation. First, two individuals from the same burial context at Osijek exhibit two different types of ACD. Second, distinct biogeographic ancestries are observed for the three individuals recovered from the pit.

The practice of ACD was particularly popular in Central Europe during the Migration Period, and so it is not unexpected that it would be identified at Osijek. However, while hundreds of cases of ACD from this region have been published, most of which have been associated with the Huns or various east Germanic tribes, including Gepids, Goths, and Alans [4, 25, 32, 55–57], only one case from Croatia, assigned to the Avar Period, has been published so far [58]. This makes our report of ACD in these individuals from Osijek very important in further elucidating the connection between Migration Period populations in Croatia and those in the rest of Central Europe. However, even more notable than the identification of ACD in the individuals from Osijek is that among the three individuals excavated from the same feature, two exhibited unique types of ACD. Specifically, SU 259 exhibited a type of ACD known as tabular oblique deformation, while the morphology of the cranium of SU 750 strongly suggested circular erect deformation [4].

As previously mentioned, both of these types of ACD are part of the “Danube Basin Group”, based on contexts by Nemeskéri [26]. The center of this group is estimated to be in present-day Hungary, but regional sub-centers corresponding to differing types and extent of deformation have been distinguished in the surrounding area [4]. In addition to Osijek, tabular oblique-type cranial deformation has been identified elsewhere on the Pannonian Plain during the Migration Period. Specifically, Molnár et al. [4] identify this type of ACD in the northeastern part of the Great Hungarian Plain in present-day Hungary at the site of Nyíregyháza-Rozsrétszőlő, a site possibly associated with Gepids or Huns. Circular deformation has been attributed to the Huns, who practiced this deformation type at a very high rate and spread its practice from the Ural Mountains to the Danube River following the Hun invasion of Eurasia [1]. Torres-Rouff and Yablonsky [1] report that Eurasian peoples, irrespective of their genetic origin and local customs, tried to emulate their prestigious Hun conquerors both culturally and physically, and that the circular-type deformation that arose as an element of Hun culture quickly lost its specific ethnic content, making it difficult to associate this type of deformation with a specific group. In addition to Osijek, erect cranial deformation has been identified at other sites on the Pannonian Plain during the Migration Period. For example, Molnár et al. [4] also identify this type of ACD at the site of Ároktő Csík-gát, where all individuals with ACD were of this type.

It remains unknown why two individuals who exhibit similarities in terms of sex and age at death, buried in the same context, would exhibit such different types of ACD. Upon identifying two distinct types of ACD in the Osijek individuals, we wondered if the distinct manifestations of ACD might have been a way of identifying the individuals as members of different cultural groups. We undertook ancient DNA analysis to further explore this idea, specifically using low-coverage shotgun sequencing data obtained using Next-Generation Sequencing to investigate the population affinities and biogeographic ancestry components of each individual. One of the individuals exhibiting ACD, SU 259, shows the greatest affinity to present-day Eastern Asians. This is the first case that we know of in which an individual from Migration Period Europe is shown to have major components of ancestry associated with present-day populations from East Asia, and expands the findings of Veeramah and colleagues [32], who identified possible components of East Asian ancestry in three Early Medieval Bavarian individuals, and Neparáczki et al. [59] who found mitogenomes of Asian origin in the conquering Hungarians of the 10^th century CE.

On the other hand, the individuals SU 261 and SU 750 exhibit components of ancestry close to modern Europeans and Near Easterners, respectively, a pattern of individual variation very similar to the individuals from Veeramah et al. [32]. The individual SU 261, who was estimated as having possible Southern European or Near Eastern ancestry is a male and did not show signs of ACD, clearly contrasting with the ACD female individuals found in Bavaria who had the same population affinity [32]. Further analysis of these individuals after deeper sequencing will enable higher precision population genomic analyses, such as admixture modelling, and will likely draw a clearer picture of these possible interactions.

Conclusion

The presented study reports an unusual archaeological feature with human remains belonging to three individuals from eastern Croatia dated to the Migration Period. Based on the archaeological context and a direct radiocarbon date, the feature was dated to the 5^th/6^th century CE and is likely associated with the Huns or Germanic tribes, a hypothesis supported by historical records. In addition to deriving from the same context, the three disarticulated and commingled skeletons found in the feature exhibit numerous similarities in biological terms: they belong to three adolescent males (between 12 and 16 years old at the time of death), they show similar pathologies suggesting a sustained and severe experience of stress, and finally, at least two individuals had a very similar diet. However, there is one significant morphological difference that distinguishes these individuals: two individuals show two different types of ACD (tabular oblique and circular erect deformations), while the cranium belonging to the third individual does not show any signs of artificial deformation. The results of genetic analyses point to an additional difference that cannot be discerned with the naked eye: these three individuals did not share the same biogeographic ancestry. More precisely, SU 259 (an individual with ACD) shows mainly East Asian ancestry, and is, to our knowledge, the first individual from the Migration Period with a majority of his ancestry originating in East Asia to be found in Europe. On the other hand, SU 750 (another individual with ACD) and SU 261 show West Eurasian ancestries similar to the patterns of population ancestry variation seen in Early Medieval Bavaria [32], though we note that our interpretations are made using low-resolution genomic data and are therefore conservative.

Although there are several hypotheses about the meaning of ACD in Pannonian Plain during the Migration Period, including an indication of social status or an indelible fashion statement, the cases from Osijek could suggest a different explanation. Based on the presented data, it seems that different types of ACD (or the lack thereof) might be associated with affiliation with a particular cultural group (at least in the case of Osijek), leading us to consider that populations with different ancestries and potentially different cultures were interacting intimately on the Pannonian Plain during the Great Migration Period. This leads us to ask: is this a random peculiarity or part of a larger-scale pattern of association between ACD type and group membership? We believe that future multi- and inter-disciplinary studies combining archaeology, bioarchaeology, history, stable isotopes analysis, and ancient DNA, conducted on a larger skeletal sample from a wider region will aid us in answering this question.

Materials and methods

Archaeological context

The human osteological remains presented in this paper were excavated in summer 2013 by Kaducej Ltd. at the Hermanov vinograd archaeological site in the city of Osijek in eastern Croatia (the exact geographic coordinates are 45° 32' 37'' N, 18° 40' 13'' E). While the 2013 excavations were a rescue campaign conducted as part of the construction of the Osijek southern bypass road [60], this famous archaeological site has been known since the end of the 19^th century when first archaeological excavations took place [61].

All necessary permits were obtained for the study, which complied with all relevant regulations. The permissions for the archaeological fieldwork and the scientific analyses of the human skeletal remains were issued by the Conservation Department of the Ministry of Culture of the Republic of Croatia in Osijek.

Hermanov vinograd is located in southwestern part of Osijek, and it is a predominantly Neolithic site associated with a Sopot culture settlement. The presence of several vertically-stratified construction phases suggest a period of active use that lasted for about 1,000 years [61]. More recent excavations at the Hermanov vinograd site revealed some archaeological features associated with later periods, mostly late Roman and early medieval times (5^th/6^th century CE). One of these features was a circular pit (designated as burial 2), approximately 1.2 m in diameter and about 1.2 m deep. Pottery fragments from the feature suggested that it dated to the 5^th century CE (the Migration Period). No traces of a settlement associated with this pit were found, making it most likely that this pit was used by a group of nomadic people living in the area during this period.

In addition to large quantities of animal bones and pottery fragments, mostly-disarticulated and partially-preserved human remains belonging to three individuals (SU 259, SU 261, and SU 750) were identified in the feature. Based on their lack of intentional positioning, the human remains appeared to be haphazardly deposited in the pit (Fig 5A and 5B). Radiocarbon dating of a human bone sample (Beta-435512) confirmed the use of this archaeological feature between 415 and 560 CE (Cal BP 1535 to 1390, 95% probability), consistent with the estimation made using pottery fragments, and confirming that the pit was in use during the Migration Period.

Download:

Fig 5. The pit from Osijek—Hermanov vinograd site during the excavation.

A) The upper layer containing commingled skeletal remains, mostly human. B) The lower layer with commingled human and animal remains.

https://doi.org/10.1371/journal.pone.0216366.g005

Bioarchaeological analysis, radiographic imaging, and stable isotopes analysis

Conventional bioarchaeological analysis was conducted at the bioarchaeological laboratory of the Institute for Anthropological Research in Zagreb, Croatia. The skeletal remains are permanently stored at the Institute for Anthropological Research under reference number O-HV-1. Due to the disarticulated and partially-preserved nature of the remains, re-individualization of the skeletal remains was conducted prior to any other analysis. The re-individualization process was carried out using two types of data. The first is archaeological context, i.e. the position of the various skeletal elements in different layers of the pit. The second type of data used was taphonomic characteristics such as the state of preservation of the bone cortex as well as the coloring of the bone. Age at death was assessed for each individual using a frequently-utilized combination of skeletal features. This included the observation of changes occurring during the development and formation of deciduous and permanent teeth [62, 63], the degree of bone ossification [64], and the length of the diaphysis of long bones [65]. Also contributing to the estimation of age at death was an examination of the ectocranial suture fusion [66].

The biological sex of the three individuals was established based on macroscopic examination focusing on the differences in pelvic and cranial morphology between adult males and females [67–69]. All pathological changes were recorded according to criteria described by Aufderheide & Rodríguez-Martín [15], and Ortner [37]. For the purpose of this study, we used an ACD classification system previously proposed by Molnár et al. [4] in their analysis of Hun-Germanic (5^th/6^th century CE) individuals from the Great Hungarian Plain. This classification system includes four different types of ACD: (i) tabular oblique, where cranial length and width both increase considerably and the general shape axis is dislocated posteriorly, (ii) tabular erect, where cranial height and width are increased as a result of anteroposterior compression and the general shape axis is approximately orthogonal in relation to the Frankfort horizontal plane, (iii) circular oblique, where a slight/pronounced flattening of the frontal and occipital bones and backward inclination in relation to the Frankfort plane is observed, and (iv) circular erect, where there is a slight or pronounced flattening of the frontal bone with a remarkable growth of the skull in height [4].

Radiographic imaging (CT scanning) of the two artificially deformed crania was conducted at the Department of Diagnostic and Interventional Radiology, University Hospital Dubrava in Zagreb, Croatia. The imaging was done utilizing a Multidetector computerized tomography (MDCT) unit (Emotion 16, Siemens AG Medical Solutions, Erlangen, Germany). Scanning parameters were 16 × 1.2 mm collimation and 1.0 reconstruction increment (RI) with 130 kV and 190 mAs respectively. Three-dimensional (3D) Volume Rendering Technique (VRT), Maximum Intensity Projection (MIP) and Multiplanar Reconstructions (MPR) were done with OsiriX MD Imaging software, v 7.0.4. (Pixmeo, Geneva, Switzerland).

Bone samples from the ribs of all three individuals were selected for the purpose of stable isotope analysis (nitrogen and carbon). Collagen for stable isotope analysis was extracted at University College Dublin, Ireland. Rib samples were demineralized in 0.5M HCl at 4C and then gelatinized in pH3 water at 70C for approximately 48 hours. The resulting solutions were filtered using Ezee filters prior to being freeze-dried. Carbon and nitrogen stable isotope analysis was conducted at University of Bradford Stable Light Isotope Laboratory, using a DeltaPLUS XL continuous flow isotope ratio mass spectrometer coupled via a ConFlo-III interface to a Thermo Flash EA 1112 elemental analyzer.

Ancient DNA analysis

Bone powders from the cochlear part of the petrous of the three individuals were collected based on research showing that this the part of the skeleton preserves more endogenous DNA than other skeletal elements [70]. All bone processing, extraction, library preparation and quality control took place in dedicated ancient DNA facilities at University College Dublin, Ireland in adherence with strict anti-contamination protocols that included the physical separation of pre- and post-PCR spaces, a one-way rule for moving between ancient and modern laboratories, and the use of UV-irradiation in all working spaces [71–73].

DNA was extracted from the bone powder [74] and sequencing libraries were built for whole-genome shotgun sequencing [75]. Negative blanks were included at each step. The libraries from the three individuals were pooled into equimolar concentrations and sequenced together with the blanks on a single-end Illumina MiSeq (50-cycle v2 kit) run with a sequencing length of 65 base pairs (bp). The software cutadapt v1.1276 [76] was used to remove sequencing adapters allowing for a minimum overlap of 1 bp between the read and the adapter and minimum read length to 17bp, and bwa v.0.7.5a-r405 [77] was used to align the sequences to the human reference genome hg19, using a disabled seed and a minimum quality score of 30. The alignment files were then sorted and duplicate sequences were removed with the software samtools v0.1.19-96b5f2294a [78].

The genetic sex of each individual was assessed following the chromosome ratio method of Skoglund and colleagues [33], and single nucleotide polymorphisms (SNPs) were called using the Pileup tool of the Genome Analyzer Tool Kit’s (GATK) v.3.3-0-g37228af [79] based on the 616,938 autosomal genomic positions of the Human Origins dataset (retrieved from Lazaridis et al. [80]). For authenticity assessment, we calculated deamination frequencies of the terminal bases of the DNA sequences using the mapDamage2 software [81] as well as the average length of all sequences, recognizing that ancient DNA will be more damaged and fragmented into shorter molecules than modern DNA. The number of human reads in the sequenced negative blanks was also considered in our evaluation of authenticity and contamination. We tried to estimate contamination based on heterozygosity of the X-chromosome in male individuals using the software ANGSD but not enough sites were covered for such tests (median of 0).

For analysis of genetic affinity between our ancient individuals and other present-day and ancient populations of interest, we used the smartpca tool from the Eigensoft package [82], which projects ancient individuals onto a principal component analysis (PCA) calculated from modern Eurasian and North African populations. The modern populations consisted of the following, from [83, 84], and retrieved from [80]: Abkhasian, Adygei, Albanian, Algerian, Altaian, Ami, Armenian, Assyrian, Balkar, Balochi, Basque, Belarusian, Brahui, Bulgarian, Burusho, Cambodian, Canary_Islander, Chechen, Chuvash, Croatian, Cypriot, Czech, Dai, Daur, Druze, Egyptian, English, Estonian, Finnish, French, Georgian, Greek, Han, Hazara, Hezhen, Hungarian, Icelandic, Iranian, Iranian_Bandari, Italian_North, Italian_South, Japanese, Jordanian, Kalash, Kalmyk, Korean, Koryak, Kumyk, Kyrgyz, Lahu, Lebanese, Lezgin, Lithuanian, Makrani, Maltese, Miao, Mongola, Mordovian, Naxi, Nganasan, Nogai, North_Ossetian, Norwegian, Orcadian, Oroqen, Palestinian, Punjabi, Romanian, Russian, Saami_WGA, Sardinian, Saudi, Scottish, Selkup, She, Sicilian, Sindhi, Spanish, Spanish_North, Syrian, Tajik, Tu, Tubalar, Tujia, Tunisian, Turkish, Turkmen, Tuvinian, Ukrainian, Ulchi, Uzbek, Xibo, Yakut, Yi. As ancient individuals we used the contemporary ones from [32], who also exhibited cranial deformation. Although Amorim et al. [85] published ancient data from a similar period as the one studied in this manuscript we believe our current analysis benefits more from the data from Veeramah et al. [32] for the following reasons: a) no ACD has been described (or at least mentioned) in the two studied sites; b) the ancestry profiles of the individuals analyzed did not show any signs of East Asian ancestry (the one individual that possibly did, CL31, was found to have very high levels of contamination). We therefore believe Veeramah et al. [32] offers a more significant dataset considering the questions we tried to answer in our study. After pruning the dataset for linkage disequilibrium in the plink software with the option—indep-pairwise 200 25 0.4, unsupervised population structure was assessed with ADMIXTURE [86], for which we ran 5 replicates of the analysis for each K value (i.e., number of ancestral populations) between 2 and 8, and then examined the replicates with the lowest cross-validation error. Outgroup-f₃ tests were used to estimate allele sharing between individuals/populations with the software qp3Pop from the AdmixTools package [83].

Supporting information

S1 Fig. Principal component analysis replicating Fig 4A on a reduced dataset composed of the combined 15,698 SNPs of the three studied individuals.

https://doi.org/10.1371/journal.pone.0216366.s001

(TIF)

Acknowledgments

Open access funding provided by University of Vienna.

References

1. Torres-Rouff C, Yablonsky LT. Cranial vault modification as a cultural artifact: A comparison of the Eurasian steppes and the Andes. Homo. 2005;56: 1–16. pmid:15901115
- View Article
- PubMed/NCBI
- Google Scholar
2. Tubbs RS, Salter EG, Oakes WJ. Artificial deformation of the human skull: A review. Clin. Anat. 2006;19: 372–377. pmid:16092127
- View Article
- PubMed/NCBI
- Google Scholar
3. Tiesler V. The Bioarchaeology of artificial cranial modifications. New York: Springer; 2014.
4. Molnár M, János I, Szucs L, Szathmáry L. Artificially deformed crania from the Hun-Germanic Period (5^th-6^th century AD) in northeastern Hungary: Historical and morphological analysis. Neurosurg. Focus. 2014;36: 1–9.
- View Article
- Google Scholar
5. Dingwall EJ. Artificial cranial deformation: A contribution to the study of ethnic mutilation. London: John Bale, Sons & Danielsson; 1931.
6. Weidenreich F. On the earliest representatives of modern mankind recovered on the soil of East Asia. Pek. Nat. Hist. Bull. 1939;13: 161–174.
- View Article
- Google Scholar
7. Hasluck M. Head deformation in the Near East. Man. 1947;47: 130–131.
- View Article
- Google Scholar
8. Blackwood B, Danby PM. A study of artificial cranial deformation in New Britain. J. R. Anthropol. Inst. 1955;85: 173–191.
- View Article
- Google Scholar
9. Meiklejohn C, Agelarakis A, Akkermans PA, Smith PEL, Solecki R. Artificial cranial deformation in the Proto-neolithic and Neolithic Near East and its possible origin: Evidence from four sites. Paléorient. 1992;18: 83–97.
- View Article
- Google Scholar
10. Durband A. Artificial cranial deformation in Kow Swamp 1 and 5: A response to Curnoe (2007). Homo. 2008;59: 261–269. pmid:18674757
- View Article
- PubMed/NCBI
- Google Scholar
11. Torres-Rouff C. A bioarchaeological analysis of crania from Pachacamac, Peru. Am. J. Phys. Anthropol. 2003;Suppl. 36: 210.
- View Article
- Google Scholar
12. Allison M, Gerzten E, Munizaga J, Santoro C, Focacci G. La practica de la deformacion craneana entre los pueblos Andinos Precolombinos. Chungara. 1981;7: 238–260.
- View Article
- Google Scholar
13. Moss ML. The pathogenesis of artificial cranial deformation. Am. J. Phys. Anthropol. 1958;16: 269–286. pmid:13649897
- View Article
- PubMed/NCBI
- Google Scholar
14. Anton SC. Intentional cranial vault deformation and induced changes of the cranial base and face. Am. J. Phys. Anthropol. 1989;79: 253–267. pmid:2662783
- View Article
- PubMed/NCBI
- Google Scholar
15. Aufderheide AC, Rodríguez-Martín C. The Cambridge encyclopedia of human paleopathology. Cambridge: Cambridge University Press; 1998.
16. McNeill WR, Newton GN. Cranial base morphology in association with intentional cranial vault deformation. Am. J. Phys. Anthropol. 1965;23: 241–254. pmid:5862964
- View Article
- PubMed/NCBI
- Google Scholar
17. Brown P. Artificial cranial deformation: A component in the variation in Pleistocene Australian Aboriginal crania. Archaeol. Ocean. 1981;16: 156–167.
- View Article
- Google Scholar
18. Rhode M, Arriaza BT. Influence of cranial deformation on facial morphology among prehistoric South Central Andean populations. Am. J. Phys. Anthropol. 2006;130: 462–470. pmid:16444729
- View Article
- PubMed/NCBI
- Google Scholar
19. Pomeroy E, Stock J, Zakrzewski S, Lahr M. A metric study of three types of artificial cranial modification from north-Central Peru. Int. J. Osteoarchaeol. 2010;20: 317–334.
- View Article
- Google Scholar
20. Zhirov EV. Ob iskusstvennoi deformatsii golovi. Kratkie soobshenia Instituta Istorii materialnoi kulturi. 1940;8: 81–88.
- View Article
- Google Scholar
21. Dembo A, Imbelloni J. Deformaciones intencionales del cuerpo humano de carácter étnico. Buenos Aires: José Anesi; 1938.
22. Romero J. Dental mutilation, trephination and cranial deformation. In: Stewart TD, editor. Handbook of Middle American Indians, Vol. 9. Austin: University of Texas Press; 1970. pp. 55–67.
23. Munizaga JR. Intentional cranial deformation in the precolumbian populations of Ecuador. Am. J. Phys. Anthropol. 1976;45: 687–694.
- View Article
- Google Scholar
24. Schijman E. Artificial cranial deformation in newborns in the pre-Columbian Andes. Child. Nerv. Syst. 2005;21: 945–950.
- View Article
- Google Scholar
25. Bereczki Z, Marcsik A. Artificial cranial deformation in Hungary. In: Mednikova M, editor. Artificial deformation of human head in Eurasian past. Moscow: Russian Academy of Sciences; 2006. pp. 96–114.
26. Nemeskéri J. Über den Künstlich Deformierten Schädel von Schöningen, Kr. Helmstedt (5.-6. Jahrhundert). Nachrichten aus Niedersachsens Urgeschichte. 1976;45: 129–154.
- View Article
- Google Scholar
27. Kiszely I. The origins of artificial cranial formation in Eurasia from the sixth millennium B.C. to the seventh century A.D. Oxford: British Archaeological Reports; 1978.
28. Hankó KI, Kiszely I. Mesterségesen deformált kora népvándorlás kori férfi koponya Tamási-Adorjánpusztáról. A Szekszárdi Béri Balogh Ádám Múzeum Évkönyve. 1972;2/3: 67–84.
- View Article
- Google Scholar
29. Gračanin H. Gepidi, Heruli, Langobardi i južna Panonija. Scrinia Slavonica. 2007;7: 7–64.
- View Article
- Google Scholar
30. Gračanin H. Južna Panonija u kasnoj antici i ranom srednjovjekovlju (od konca 4. do konca 11. stoljeća). Zagreb: Plejada; 2011.
31. Cottin M, Hossein Khonsari R, Friess M. Assessing cranial plasticity in humans: The impact of artificial deformation on masticatory and basicranial structures. C. R Palevol. 2017;16: 545–556.
- View Article
- Google Scholar
32. Veeramah KR, Rott A, Groß M, van Dorp L, López S, Kirsanow K, et al. Population genomic analysis of elongated skulls reveals extensive female-biased immigration in Early Medieval Bavaria. Proc. Natl. Acad. Sci. U. S. A. 2018;115: 3494–3499. pmid:29531040
- View Article
- PubMed/NCBI
- Google Scholar
33. Skoglund P, Storå J, Götherström A, Jakobsson M. Accurate sex identification of ancient human remains using DNA shotgun sequencing. J. Archaeol. Sci. 2013;40: 4477–4482.
- View Article
- Google Scholar
34. Crișan C, Lăzărescu VA. Forgotten ritual connected to agriculture in the Early Migration Period—feature G27 from Ernei “Carieră” (Mureș County). Ephemeris Napocensis. 2010;20: 221–256.
- View Article
- Google Scholar
35. Verano J. Paleopathology and human sacrifice in the Moche culture of ancient Peru. New Orleans: Tulane University; 2001.
36. Stuart-Macadam P. Porotic hyperostosis: representative of a childhood condition. Am. J. Phys. Anthropol. 1985;66: 391–398. pmid:3887936
- View Article
- PubMed/NCBI
- Google Scholar
37. Ortner DJ. Identification of pathological conditions in human skeletal remains. San Diego: Academic Press; 2003.
38. Walker PL, Bathurst RR, Richman R, Gjerdrum T, Andrushko VA. The causes of porotic hyperostosis and cribra orbitalia: a reappraisal of the iron-deficiency anemia hypothesis. Am. J. Phys. Anthropol. 2009;139: 109–125. pmid:19280675
- View Article
- PubMed/NCBI
- Google Scholar
39. Larsen CS. Bioarchaeology. Interpreting behavior from the human skeleton. Cambridge: Cambridge University Press; 1997.
40. Wheeler SM. Nutritional and disease stress of juveniles from the Dakhleh Oasis, Egypt. Int. J. Osteoarchaeol. 2012;22: 219–234.
- View Article
- Google Scholar
41. Goodman AH, Armelagos GJ, Rose JC. Enamel hypoplasias as indicators of stress in three prehistoric populations from Illinois. Hum. Biol. 1980;52: 515–528. pmid:7005071
- View Article
- PubMed/NCBI
- Google Scholar
42. Hillson S. Tooth development in human evolution and bioarchaeology. Cambridge: Cambridge University Press; 2014.
43. Roberts CA, Manchester K. Archaeology of disease. Ithaca: Cornell University Press; 2005.
44. Łukasik S, Krenz-Niedbala M. Age of linear enamel hypoplasia formation based on Massler and colleagues’ and Reid and Dean’s standards in a Polish sample dated to 13th–18th century CE. Homo. 2014;65: 296–310. pmid:24767821
- View Article
- PubMed/NCBI
- Google Scholar
45. Blakey ML, Leslie TE, Reidy JP. Frequency and chronological distribution of dental enamel hypoplasia in enslaved Americans: A test of the weaning hypothesis. Am. J. Phys. Anthropol. 194;95: 371–383. pmid:7864059
- View Article
- PubMed/NCBI
- Google Scholar
46. Andrić S. Southern Pannonia during the age of the Great Migrations. Scrinia Slavonica. 2002;2: 117–167.
- View Article
- Google Scholar
47. Šlaus M, Pećina Šlaus N, Brkić H. Life stress on the Roman limes in continental Croatia. Homo. 2004;54: 240–263. pmid:15216669
- View Article
- PubMed/NCBI
- Google Scholar
48. Novak M, Premužić Z, Vyroubal V. Šlaus M. Bioarchaeology of the late Antique population from Štrbinci. Arh. rad. raspr. 2009;16: 265–326.
- View Article
- Google Scholar
49. Miladinović-Radmilović N, Vulović D, Đukić K. Health status of children in ancient Sirmium. Starinar. 2016;66: 65–80.
- View Article
- Google Scholar
50. Hakenbeck SE, Evans J, Chapman H, Fóthi E. Practising pastoralism in an agricultural environment: an isotopic analysis of the impact of the Hunnic incursions on Pannonian populations. PLOS ONE. 2017;12: e0173079. pmid:28328951
- View Article
- PubMed/NCBI
- Google Scholar
51. Noche-Dowdy L. Multi-isotope analysis to reconstruct dietary and migration patterns of an Avar population from Sajópetri, Hungary, AD 568–895. M.Sc. Thesis. University of South Florida. 2015.
52. Vidal-Ronchas R, Rajić Šikanjić P, Premužić Z, Rapan Papeša A, Lightfoot E. Diet, sex, and social status in the Late Avar period: stable isotope investigations at Nuštar cemetery, Croatia. Archaeol. Anthropol. Sci. 2018;
- View Article
- Google Scholar
53. Motuzaite-Matuzeiciute G, Staff R, Hunt HV, Xinyi L, Jones K. The early chronology of broomcorn millet (Panicum miliaceum) in Europe. Antiquity. 2013;87: 1073–1085.
- View Article
- Google Scholar
54. Rapan Papeša A, Kenéz A, Pető Á. The archaeobotanical assessment of grave samples from the Avar age cemetery of Nuštar (Eastern Croatia). Pril. Inst. arheol. Zagrebu. 2015;32: 261–288.
- View Article
- Google Scholar
55. Mikić Ž. Deux nécropoles de la Grande migration des peuples à Viminacium. Balcanica. 2007;38: 45–55.
- View Article
- Google Scholar
56. Miriţoiu N. Cranii deformate din necropola de la Piatra Frecăţei. Peuce (Serie Nouă)–Studii şi cercetari de istorie şi arheologie. 2011;9: 539–574.
- View Article
- Google Scholar
57. Schmölzer A. The long-heads—strangers of the east? Artificial cranial deformation in Austria. In: Franković F, Popović A, Fileš M, Kramberger J, editors. Proceedings of the international graduate and doctoral student conference, October 14–15 2016, Faculty of Humanities and Social Sciences, Zagreb, Croatia. Zagreb: Faculty of Humanities and Social Sciences, University of Zagreb; 2018. pp. 58–69.
58. Šlaus M. The bioarchaeology of continental Croatia. Oxford: Archaeopress; 2002.
59. Neparáczki E, Maróti Z, Kalmár T, Kocsy K, Maár K, Bihari P, et al. Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya origin in the conquering Hungarians. PLOS ONE; 2018;13: e0205920. pmid:30335830
- View Article
- PubMed/NCBI
- Google Scholar
60. Los Dž. Osijek—Hermanov vinograd (AN 6 i AN 7). In: Wiewegh Z, editor. Hrvatski arheološki godišnjak 10/2013. Zagreb: Ministarstvo kulture Republike Hrvatske; 2013. pp. 27–34.
61. Šimić J. Hermanov vinograd—arheološko nalazište mlađeg kamenog doba u Osijeku. Osijek: Muzej Slavonije; 2008.
62. Moorrees CFA, Fanning EA, Hunt EE. Age variation of formation stages for ten permanent teeth. J. Dent. Res. 1963;42: 1490–1502. pmid:14081973
- View Article
- PubMed/NCBI
- Google Scholar
63. Gustafson G, Koch G. Age estimation up to 16 years of age based on dental development. Odontol. Revy. 1974:25; 297–306. pmid:4530955
- View Article
- PubMed/NCBI
- Google Scholar
64. Scheuer JL, Black S. Developmental juvenile osteology. San Diego: Academic Press; 2000.
65. Maresh MM. Measurements from roentgenograms. In: McCammon RW, editor. Human growth and development. Springfield: C.C. Thomas; 1970. pp. 157–200.
66. Meindl RS, Lovejoy CO. Ectocranial suture closure: a revised method for the determination of skeletal age at death based on the lateral-anterior sutures. Am. J. Phys. Anthropol. 1985;68: 57–66. pmid:4061602
- View Article
- PubMed/NCBI
- Google Scholar
67. Krogman WM, Işcan MY. The human skeleton in forensic medicine. Springfield: C.C. Thomas; 1986.
68. Buikstra JE, Ubelaker DH. Standards for data collection from human skeletal remains. Fayetteville: Arkansas Archaeological Survey; 1994.
69. Bass WM. Human osteology. A laboratory and field manual of the human skeleton. Columbia: Missouri Archaeological Society; 1995.
70. Pinhasi R, Fernandes D, Sirak K, Novak M, Connell S, Alpaslan-Roodenberg S, et al. Optimal Ancient DNA yields from the inner ear part of the human petrous bone. PLOS ONE. 2015;10: e0129102. pmid:26086078
- View Article
- PubMed/NCBI
- Google Scholar
71. Cooper A, Poinar HN. Ancient DNA: Do it right or not at all. Science. 2000;289: 1139. pmid:10970224
- View Article
- PubMed/NCBI
- Google Scholar
72. Champlot S, Berthelot C, Pruvost M, Bennett EA, Grange T, Geigl EM. An efficient multistrategy DNA decontamination procedure of PCR reagents for hypersensitive PCR applications. PLOS ONE. 2010;5: e13042. pmid:20927390
- View Article
- PubMed/NCBI
- Google Scholar
73. Llamas B, Valverde G, Fehren-Schmitz L, Weyrich LS, Cooper A, Haak W. From the field to the laboratory: Controlling DNA contamination in human ancient DNA research in the high-throughput sequencing era. Sci. Technol. Archaeol. Res. 2017;3: 1–14.
- View Article
- Google Scholar
74. Dabney J, Knapp M, Glocke I, Gansauge MT, Weihmann A, Nickel B, et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. U. S. A. 2013;110: 15758–15763. pmid:24019490
- View Article
- PubMed/NCBI
- Google Scholar
75. Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 2010;6: db.prot5448.
- View Article
- Google Scholar
76. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011;17: 10–12.
- View Article
- Google Scholar
77. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25: 1754–1760. pmid:19451168
- View Article
- PubMed/NCBI
- Google Scholar
78. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25: 2078–2079. pmid:19505943
- View Article
- PubMed/NCBI
- Google Scholar
79. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20: 1297–1303. pmid:20644199
- View Article
- PubMed/NCBI
- Google Scholar
80. Lazaridis I, Nadel D, Rollefson G, Merrett DC, Rohland N, Mallick S, et al. Genomic insights into the origin of farming in the ancient Near East. Nature. 2016;536: 419–424. pmid:27459054
- View Article
- PubMed/NCBI
- Google Scholar
81. Jónsson H, Ginolhac A, Schubert M, Johnson PLF, Orlando L. mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics. 2013;29: 1682–1684. pmid:23613487
- View Article
- PubMed/NCBI
- Google Scholar
82. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLOS Genet. 2006;12: e190.
- View Article
- Google Scholar
83. Patterson N, Moorjani P, Luo Y, Mallick S, Rohland N, Zhan Y, et al. Ancient admixture in human history. Genetics. 2012;192: 1065–1093. pmid:22960212
- View Article
- PubMed/NCBI
- Google Scholar
84. Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014;513: 409–413. pmid:25230663
- View Article
- PubMed/NCBI
- Google Scholar
85. Amorim CEG, Vai S, Posth C, Modi A, Koncz I, Hakenbeck S, et al. Understanding 6th-century barbarian social organization and migration through paleogenomics. Nat Commun. 2018;9: 3547. pmid:30206220
- View Article
- PubMed/NCBI
- Google Scholar
86. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19: 1655–1664. pmid:19648217
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 1. Torres-Rouff C, Yablonsky LT. Cranial vault modification as a cultural artifact: A comparison of the Eurasian steppes and the Andes. Homo. 2005;56: 1–16. pmid:15901115
View Article
PubMed/NCBI
Google Scholar

[2] View Article

[3] PubMed/NCBI

[4] Google Scholar

[ref2] 2. Tubbs RS, Salter EG, Oakes WJ. Artificial deformation of the human skull: A review. Clin. Anat. 2006;19: 372–377. pmid:16092127
View Article
PubMed/NCBI
Google Scholar

[6] View Article

[7] PubMed/NCBI

[8] Google Scholar

[ref3] 3. Tiesler V. The Bioarchaeology of artificial cranial modifications. New York: Springer; 2014.

[ref4] 4. Molnár M, János I, Szucs L, Szathmáry L. Artificially deformed crania from the Hun-Germanic Period (5^th-6^th century AD) in northeastern Hungary: Historical and morphological analysis. Neurosurg. Focus. 2014;36: 1–9.
View Article
Google Scholar

[11] View Article

[12] Google Scholar

[ref5] 5. Dingwall EJ. Artificial cranial deformation: A contribution to the study of ethnic mutilation. London: John Bale, Sons & Danielsson; 1931.

[ref6] 6. Weidenreich F. On the earliest representatives of modern mankind recovered on the soil of East Asia. Pek. Nat. Hist. Bull. 1939;13: 161–174.
View Article
Google Scholar

[15] View Article

[16] Google Scholar

[ref7] 7. Hasluck M. Head deformation in the Near East. Man. 1947;47: 130–131.
View Article
Google Scholar

[18] View Article

[19] Google Scholar

[ref8] 8. Blackwood B, Danby PM. A study of artificial cranial deformation in New Britain. J. R. Anthropol. Inst. 1955;85: 173–191.
View Article
Google Scholar

[21] View Article

[22] Google Scholar

[ref9] 9. Meiklejohn C, Agelarakis A, Akkermans PA, Smith PEL, Solecki R. Artificial cranial deformation in the Proto-neolithic and Neolithic Near East and its possible origin: Evidence from four sites. Paléorient. 1992;18: 83–97.
View Article
Google Scholar

[24] View Article

[25] Google Scholar

[ref10] 10. Durband A. Artificial cranial deformation in Kow Swamp 1 and 5: A response to Curnoe (2007). Homo. 2008;59: 261–269. pmid:18674757
View Article
PubMed/NCBI
Google Scholar

[27] View Article

[28] PubMed/NCBI

[29] Google Scholar

[ref11] 11. Torres-Rouff C. A bioarchaeological analysis of crania from Pachacamac, Peru. Am. J. Phys. Anthropol. 2003;Suppl. 36: 210.
View Article
Google Scholar

[31] View Article

[32] Google Scholar

[ref12] 12. Allison M, Gerzten E, Munizaga J, Santoro C, Focacci G. La practica de la deformacion craneana entre los pueblos Andinos Precolombinos. Chungara. 1981;7: 238–260.
View Article
Google Scholar

[34] View Article

[35] Google Scholar

[ref13] 13. Moss ML. The pathogenesis of artificial cranial deformation. Am. J. Phys. Anthropol. 1958;16: 269–286. pmid:13649897
View Article
PubMed/NCBI
Google Scholar

[37] View Article

[38] PubMed/NCBI

[39] Google Scholar

[ref14] 14. Anton SC. Intentional cranial vault deformation and induced changes of the cranial base and face. Am. J. Phys. Anthropol. 1989;79: 253–267. pmid:2662783
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

[ref15] 15. Aufderheide AC, Rodríguez-Martín C. The Cambridge encyclopedia of human paleopathology. Cambridge: Cambridge University Press; 1998.

[ref16] 16. McNeill WR, Newton GN. Cranial base morphology in association with intentional cranial vault deformation. Am. J. Phys. Anthropol. 1965;23: 241–254. pmid:5862964
View Article
PubMed/NCBI
Google Scholar

[46] View Article

[47] PubMed/NCBI

[48] Google Scholar

[ref17] 17. Brown P. Artificial cranial deformation: A component in the variation in Pleistocene Australian Aboriginal crania. Archaeol. Ocean. 1981;16: 156–167.
View Article
Google Scholar

[50] View Article

[51] Google Scholar

[ref18] 18. Rhode M, Arriaza BT. Influence of cranial deformation on facial morphology among prehistoric South Central Andean populations. Am. J. Phys. Anthropol. 2006;130: 462–470. pmid:16444729
View Article
PubMed/NCBI
Google Scholar

[53] View Article

[54] PubMed/NCBI

[55] Google Scholar

[ref19] 19. Pomeroy E, Stock J, Zakrzewski S, Lahr M. A metric study of three types of artificial cranial modification from north-Central Peru. Int. J. Osteoarchaeol. 2010;20: 317–334.
View Article
Google Scholar

[57] View Article

[58] Google Scholar

[ref20] 20. Zhirov EV. Ob iskusstvennoi deformatsii golovi. Kratkie soobshenia Instituta Istorii materialnoi kulturi. 1940;8: 81–88.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref21] 21. Dembo A, Imbelloni J. Deformaciones intencionales del cuerpo humano de carácter étnico. Buenos Aires: José Anesi; 1938.

[ref22] 22. Romero J. Dental mutilation, trephination and cranial deformation. In: Stewart TD, editor. Handbook of Middle American Indians, Vol. 9. Austin: University of Texas Press; 1970. pp. 55–67.

[ref23] 23. Munizaga JR. Intentional cranial deformation in the precolumbian populations of Ecuador. Am. J. Phys. Anthropol. 1976;45: 687–694.
View Article
Google Scholar

[65] View Article

[66] Google Scholar

[ref24] 24. Schijman E. Artificial cranial deformation in newborns in the pre-Columbian Andes. Child. Nerv. Syst. 2005;21: 945–950.
View Article
Google Scholar

[68] View Article

[69] Google Scholar

[ref25] 25. Bereczki Z, Marcsik A. Artificial cranial deformation in Hungary. In: Mednikova M, editor. Artificial deformation of human head in Eurasian past. Moscow: Russian Academy of Sciences; 2006. pp. 96–114.

[ref26] 26. Nemeskéri J. Über den Künstlich Deformierten Schädel von Schöningen, Kr. Helmstedt (5.-6. Jahrhundert). Nachrichten aus Niedersachsens Urgeschichte. 1976;45: 129–154.
View Article
Google Scholar

[72] View Article

[73] Google Scholar

[ref27] 27. Kiszely I. The origins of artificial cranial formation in Eurasia from the sixth millennium B.C. to the seventh century A.D. Oxford: British Archaeological Reports; 1978.

[ref28] 28. Hankó KI, Kiszely I. Mesterségesen deformált kora népvándorlás kori férfi koponya Tamási-Adorjánpusztáról. A Szekszárdi Béri Balogh Ádám Múzeum Évkönyve. 1972;2/3: 67–84.
View Article
Google Scholar

[76] View Article

[77] Google Scholar

[ref29] 29. Gračanin H. Gepidi, Heruli, Langobardi i južna Panonija. Scrinia Slavonica. 2007;7: 7–64.
View Article
Google Scholar

[79] View Article

[80] Google Scholar

[ref30] 30. Gračanin H. Južna Panonija u kasnoj antici i ranom srednjovjekovlju (od konca 4. do konca 11. stoljeća). Zagreb: Plejada; 2011.

[ref31] 31. Cottin M, Hossein Khonsari R, Friess M. Assessing cranial plasticity in humans: The impact of artificial deformation on masticatory and basicranial structures. C. R Palevol. 2017;16: 545–556.
View Article
Google Scholar

[83] View Article

[84] Google Scholar

[ref32] 32. Veeramah KR, Rott A, Groß M, van Dorp L, López S, Kirsanow K, et al. Population genomic analysis of elongated skulls reveals extensive female-biased immigration in Early Medieval Bavaria. Proc. Natl. Acad. Sci. U. S. A. 2018;115: 3494–3499. pmid:29531040
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref33] 33. Skoglund P, Storå J, Götherström A, Jakobsson M. Accurate sex identification of ancient human remains using DNA shotgun sequencing. J. Archaeol. Sci. 2013;40: 4477–4482.
View Article
Google Scholar

[90] View Article

[91] Google Scholar

[ref34] 34. Crișan C, Lăzărescu VA. Forgotten ritual connected to agriculture in the Early Migration Period—feature G27 from Ernei “Carieră” (Mureș County). Ephemeris Napocensis. 2010;20: 221–256.
View Article
Google Scholar

[93] View Article

[94] Google Scholar

[ref35] 35. Verano J. Paleopathology and human sacrifice in the Moche culture of ancient Peru. New Orleans: Tulane University; 2001.

[ref36] 36. Stuart-Macadam P. Porotic hyperostosis: representative of a childhood condition. Am. J. Phys. Anthropol. 1985;66: 391–398. pmid:3887936
View Article
PubMed/NCBI
Google Scholar

[97] View Article

[98] PubMed/NCBI

[99] Google Scholar

[ref37] 37. Ortner DJ. Identification of pathological conditions in human skeletal remains. San Diego: Academic Press; 2003.

[ref38] 38. Walker PL, Bathurst RR, Richman R, Gjerdrum T, Andrushko VA. The causes of porotic hyperostosis and cribra orbitalia: a reappraisal of the iron-deficiency anemia hypothesis. Am. J. Phys. Anthropol. 2009;139: 109–125. pmid:19280675
View Article
PubMed/NCBI
Google Scholar

[102] View Article

[103] PubMed/NCBI

[104] Google Scholar

[ref39] 39. Larsen CS. Bioarchaeology. Interpreting behavior from the human skeleton. Cambridge: Cambridge University Press; 1997.

[ref40] 40. Wheeler SM. Nutritional and disease stress of juveniles from the Dakhleh Oasis, Egypt. Int. J. Osteoarchaeol. 2012;22: 219–234.
View Article
Google Scholar

[107] View Article

[108] Google Scholar

[ref41] 41. Goodman AH, Armelagos GJ, Rose JC. Enamel hypoplasias as indicators of stress in three prehistoric populations from Illinois. Hum. Biol. 1980;52: 515–528. pmid:7005071
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref42] 42. Hillson S. Tooth development in human evolution and bioarchaeology. Cambridge: Cambridge University Press; 2014.

[ref43] 43. Roberts CA, Manchester K. Archaeology of disease. Ithaca: Cornell University Press; 2005.

[ref44] 44. Łukasik S, Krenz-Niedbala M. Age of linear enamel hypoplasia formation based on Massler and colleagues’ and Reid and Dean’s standards in a Polish sample dated to 13th–18th century CE. Homo. 2014;65: 296–310. pmid:24767821
View Article
PubMed/NCBI
Google Scholar

[116] View Article

[117] PubMed/NCBI

[118] Google Scholar

[ref45] 45. Blakey ML, Leslie TE, Reidy JP. Frequency and chronological distribution of dental enamel hypoplasia in enslaved Americans: A test of the weaning hypothesis. Am. J. Phys. Anthropol. 194;95: 371–383. pmid:7864059
View Article
PubMed/NCBI
Google Scholar

[120] View Article

[121] PubMed/NCBI

[122] Google Scholar

[ref46] 46. Andrić S. Southern Pannonia during the age of the Great Migrations. Scrinia Slavonica. 2002;2: 117–167.
View Article
Google Scholar

[124] View Article

[125] Google Scholar

[ref47] 47. Šlaus M, Pećina Šlaus N, Brkić H. Life stress on the Roman limes in continental Croatia. Homo. 2004;54: 240–263. pmid:15216669
View Article
PubMed/NCBI
Google Scholar

[127] View Article

[128] PubMed/NCBI

[129] Google Scholar

[ref48] 48. Novak M, Premužić Z, Vyroubal V. Šlaus M. Bioarchaeology of the late Antique population from Štrbinci. Arh. rad. raspr. 2009;16: 265–326.
View Article
Google Scholar

[131] View Article

[132] Google Scholar

[ref49] 49. Miladinović-Radmilović N, Vulović D, Đukić K. Health status of children in ancient Sirmium. Starinar. 2016;66: 65–80.
View Article
Google Scholar

[134] View Article

[135] Google Scholar

[ref50] 50. Hakenbeck SE, Evans J, Chapman H, Fóthi E. Practising pastoralism in an agricultural environment: an isotopic analysis of the impact of the Hunnic incursions on Pannonian populations. PLOS ONE. 2017;12: e0173079. pmid:28328951
View Article
PubMed/NCBI
Google Scholar

[137] View Article

[138] PubMed/NCBI

[139] Google Scholar

[ref51] 51. Noche-Dowdy L. Multi-isotope analysis to reconstruct dietary and migration patterns of an Avar population from Sajópetri, Hungary, AD 568–895. M.Sc. Thesis. University of South Florida. 2015.

[ref52] 52. Vidal-Ronchas R, Rajić Šikanjić P, Premužić Z, Rapan Papeša A, Lightfoot E. Diet, sex, and social status in the Late Avar period: stable isotope investigations at Nuštar cemetery, Croatia. Archaeol. Anthropol. Sci. 2018;
View Article
Google Scholar

[142] View Article

[143] Google Scholar

[ref53] 53. Motuzaite-Matuzeiciute G, Staff R, Hunt HV, Xinyi L, Jones K. The early chronology of broomcorn millet (Panicum miliaceum) in Europe. Antiquity. 2013;87: 1073–1085.
View Article
Google Scholar

[145] View Article

[146] Google Scholar

[ref54] 54. Rapan Papeša A, Kenéz A, Pető Á. The archaeobotanical assessment of grave samples from the Avar age cemetery of Nuštar (Eastern Croatia). Pril. Inst. arheol. Zagrebu. 2015;32: 261–288.
View Article
Google Scholar

[148] View Article

[149] Google Scholar

[ref55] 55. Mikić Ž. Deux nécropoles de la Grande migration des peuples à Viminacium. Balcanica. 2007;38: 45–55.
View Article
Google Scholar

[151] View Article

[152] Google Scholar

[ref56] 56. Miriţoiu N. Cranii deformate din necropola de la Piatra Frecăţei. Peuce (Serie Nouă)–Studii şi cercetari de istorie şi arheologie. 2011;9: 539–574.
View Article
Google Scholar

[154] View Article

[155] Google Scholar

[ref57] 57. Schmölzer A. The long-heads—strangers of the east? Artificial cranial deformation in Austria. In: Franković F, Popović A, Fileš M, Kramberger J, editors. Proceedings of the international graduate and doctoral student conference, October 14–15 2016, Faculty of Humanities and Social Sciences, Zagreb, Croatia. Zagreb: Faculty of Humanities and Social Sciences, University of Zagreb; 2018. pp. 58–69.

[ref58] 58. Šlaus M. The bioarchaeology of continental Croatia. Oxford: Archaeopress; 2002.

[ref59] 59. Neparáczki E, Maróti Z, Kalmár T, Kocsy K, Maár K, Bihari P, et al. Mitogenomic data indicate admixture components of Central-Inner Asian and Srubnaya origin in the conquering Hungarians. PLOS ONE; 2018;13: e0205920. pmid:30335830
View Article
PubMed/NCBI
Google Scholar

[159] View Article

[160] PubMed/NCBI

[161] Google Scholar

[ref60] 60. Los Dž. Osijek—Hermanov vinograd (AN 6 i AN 7). In: Wiewegh Z, editor. Hrvatski arheološki godišnjak 10/2013. Zagreb: Ministarstvo kulture Republike Hrvatske; 2013. pp. 27–34.

[ref61] 61. Šimić J. Hermanov vinograd—arheološko nalazište mlađeg kamenog doba u Osijeku. Osijek: Muzej Slavonije; 2008.

[ref62] 62. Moorrees CFA, Fanning EA, Hunt EE. Age variation of formation stages for ten permanent teeth. J. Dent. Res. 1963;42: 1490–1502. pmid:14081973
View Article
PubMed/NCBI
Google Scholar

[165] View Article

[166] PubMed/NCBI

[167] Google Scholar

[ref63] 63. Gustafson G, Koch G. Age estimation up to 16 years of age based on dental development. Odontol. Revy. 1974:25; 297–306. pmid:4530955
View Article
PubMed/NCBI
Google Scholar

[169] View Article

[170] PubMed/NCBI

[171] Google Scholar

[ref64] 64. Scheuer JL, Black S. Developmental juvenile osteology. San Diego: Academic Press; 2000.

[ref65] 65. Maresh MM. Measurements from roentgenograms. In: McCammon RW, editor. Human growth and development. Springfield: C.C. Thomas; 1970. pp. 157–200.

[ref66] 66. Meindl RS, Lovejoy CO. Ectocranial suture closure: a revised method for the determination of skeletal age at death based on the lateral-anterior sutures. Am. J. Phys. Anthropol. 1985;68: 57–66. pmid:4061602
View Article
PubMed/NCBI
Google Scholar

[175] View Article

[176] PubMed/NCBI

[177] Google Scholar

[ref67] 67. Krogman WM, Işcan MY. The human skeleton in forensic medicine. Springfield: C.C. Thomas; 1986.

[ref68] 68. Buikstra JE, Ubelaker DH. Standards for data collection from human skeletal remains. Fayetteville: Arkansas Archaeological Survey; 1994.

[ref69] 69. Bass WM. Human osteology. A laboratory and field manual of the human skeleton. Columbia: Missouri Archaeological Society; 1995.

[ref70] 70. Pinhasi R, Fernandes D, Sirak K, Novak M, Connell S, Alpaslan-Roodenberg S, et al. Optimal Ancient DNA yields from the inner ear part of the human petrous bone. PLOS ONE. 2015;10: e0129102. pmid:26086078
View Article
PubMed/NCBI
Google Scholar

[182] View Article

[183] PubMed/NCBI

[184] Google Scholar

[ref71] 71. Cooper A, Poinar HN. Ancient DNA: Do it right or not at all. Science. 2000;289: 1139. pmid:10970224
View Article
PubMed/NCBI
Google Scholar

[186] View Article

[187] PubMed/NCBI

[188] Google Scholar

[ref72] 72. Champlot S, Berthelot C, Pruvost M, Bennett EA, Grange T, Geigl EM. An efficient multistrategy DNA decontamination procedure of PCR reagents for hypersensitive PCR applications. PLOS ONE. 2010;5: e13042. pmid:20927390
View Article
PubMed/NCBI
Google Scholar

[190] View Article

[191] PubMed/NCBI

[192] Google Scholar

[ref73] 73. Llamas B, Valverde G, Fehren-Schmitz L, Weyrich LS, Cooper A, Haak W. From the field to the laboratory: Controlling DNA contamination in human ancient DNA research in the high-throughput sequencing era. Sci. Technol. Archaeol. Res. 2017;3: 1–14.
View Article
Google Scholar

[194] View Article

[195] Google Scholar

[ref74] 74. Dabney J, Knapp M, Glocke I, Gansauge MT, Weihmann A, Nickel B, et al. Complete mitochondrial genome sequence of a Middle Pleistocene cave bear reconstructed from ultrashort DNA fragments. Proc. Natl. Acad. Sci. U. S. A. 2013;110: 15758–15763. pmid:24019490
View Article
PubMed/NCBI
Google Scholar

[197] View Article

[198] PubMed/NCBI

[199] Google Scholar

[ref75] 75. Meyer M, Kircher M. Illumina sequencing library preparation for highly multiplexed target capture and sequencing. Cold Spring Harb. Protoc. 2010;6: db.prot5448.
View Article
Google Scholar

[201] View Article

[202] Google Scholar

[ref76] 76. Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet. J. 2011;17: 10–12.
View Article
Google Scholar

[204] View Article

[205] Google Scholar

[ref77] 77. Li H, Durbin R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009;25: 1754–1760. pmid:19451168
View Article
PubMed/NCBI
Google Scholar

[207] View Article

[208] PubMed/NCBI

[209] Google Scholar

[ref78] 78. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25: 2078–2079. pmid:19505943
View Article
PubMed/NCBI
Google Scholar

[211] View Article

[212] PubMed/NCBI

[213] Google Scholar

[ref79] 79. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20: 1297–1303. pmid:20644199
View Article
PubMed/NCBI
Google Scholar

[215] View Article

[216] PubMed/NCBI

[217] Google Scholar

[ref80] 80. Lazaridis I, Nadel D, Rollefson G, Merrett DC, Rohland N, Mallick S, et al. Genomic insights into the origin of farming in the ancient Near East. Nature. 2016;536: 419–424. pmid:27459054
View Article
PubMed/NCBI
Google Scholar

[219] View Article

[220] PubMed/NCBI

[221] Google Scholar

[ref81] 81. Jónsson H, Ginolhac A, Schubert M, Johnson PLF, Orlando L. mapDamage2.0: fast approximate Bayesian estimates of ancient DNA damage parameters. Bioinformatics. 2013;29: 1682–1684. pmid:23613487
View Article
PubMed/NCBI
Google Scholar

[223] View Article

[224] PubMed/NCBI

[225] Google Scholar

[ref82] 82. Patterson N, Price AL, Reich D. Population structure and eigenanalysis. PLOS Genet. 2006;12: e190.
View Article
Google Scholar

[227] View Article

[228] Google Scholar

[ref83] 83. Patterson N, Moorjani P, Luo Y, Mallick S, Rohland N, Zhan Y, et al. Ancient admixture in human history. Genetics. 2012;192: 1065–1093. pmid:22960212
View Article
PubMed/NCBI
Google Scholar

[230] View Article

[231] PubMed/NCBI

[232] Google Scholar

[ref84] 84. Lazaridis I, Patterson N, Mittnik A, Renaud G, Mallick S, Kirsanow K, et al. Ancient human genomes suggest three ancestral populations for present-day Europeans. Nature. 2014;513: 409–413. pmid:25230663
View Article
PubMed/NCBI
Google Scholar

[234] View Article

[235] PubMed/NCBI

[236] Google Scholar

[ref85] 85. Amorim CEG, Vai S, Posth C, Modi A, Koncz I, Hakenbeck S, et al. Understanding 6th-century barbarian social organization and migration through paleogenomics. Nat Commun. 2018;9: 3547. pmid:30206220
View Article
PubMed/NCBI
Google Scholar

[238] View Article

[239] PubMed/NCBI

[240] Google Scholar

[ref86] 86. Alexander DH, Novembre J, Lange K. Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 2009;19: 1655–1664. pmid:19648217
View Article
PubMed/NCBI
Google Scholar

[242] View Article

[243] PubMed/NCBI

[244] Google Scholar

Figures

Abstract

Introduction

Results

Discussion

Conclusion

Materials and methods

Archaeological context

Bioarchaeological analysis, radiographic imaging, and stable isotopes analysis

Ancient DNA analysis

Supporting information

S1 Fig. Principal component analysis replicating Fig 4A on a reduced dataset composed of the combined 15,698 SNPs of the three studied individuals.

Acknowledgments

References