Craniometric analysis of European Upper Palaeolithic and Mesolithic samples supports discontinuity at the Last Glacial Maximum

Bioanthropological analysis of human remains from the archaic and classic period discovered in Puyil cave, Mexico

OBJECTIVES

Determine the geographic place of origin and maternal lineage of prehistoric human skeletal remains discovered in Puyil Cave, Tabasco State, Mexico, located in a region currently populated by Olmec, Zoque and Maya populations.

MATERIALS AND METHODS

All specimens were radiocarbon (14 C) dated (beta analytic), had dental modifications classified, and had an analysis of 13 homologous reference points conducted to evaluate artificial cranial deformation (ACD). Following DNA purification, hypervariable region I (HVR-1) of the mitogenome was amplified and Sanger sequenced. Finally, Next Generation Sequencing (NGS) was performed for total DNA. Mitochondrial DNA (mtDNA) variants and haplogroups were determined using BioEdit 7.2 and IGV software and confirmed with MITOMASTER and WebHome softwares.

RESULTS

Radiocarbon dating (14 C) demonstrated that the inhabitants of Puyil Cave lived during the Archaic and Classic Periods and displayed tabular oblique and tabular mimetic ACD. These pre-Hispanic remains exhibited five mtDNA lineages: A, A2, C1, C1c and D4. Network analysis revealed a close genetic affinity between pre-Hispanic Puyil Cave inhabitants and contemporary Maya subpopulations from Mexico and Guatemala, as well as individuals from Bolivia, Brazil, the Dominican Republic, and China.

CONCLUSIONS

Our results elucidate the dispersal of pre-Hispanic Olmec and Maya ancestors and suggest that ACD practices are closely related to Olmec and Maya practices. Additionally, we conclude that ACD has likely been practiced in the region since the Middle-Archaic Period.

Inferring human neutral genetic variation from craniodental phenotypes

There is a growing consensus that global patterns of modern human cranial and dental variation are shaped largely by neutral evolutionary processes, suggesting that craniodental features can be used as reliable proxies for inferring population structure and history in bioarchaeological, forensic, and paleoanthropological contexts. However, there is disagreement on whether certain types of data preserve a neutral signature to a greater degree than others. Here, we address this unresolved question and systematically test the relative neutrality of four standard metric and nonmetric craniodental data types employing an extensive computational genotype-phenotype comparison across modern populations from around the world. Our computation draws on the largest existing data sets currently available, while accounting for geographically structured environmental variation, population sampling uncertainty, disparate numbers of phenotypic variables, and stochastic variation inherent to a neutral model of evolution. Our results reveal that the four data types differentially capture neutral genomic variation, with highest signals preserved in dental nonmetric and cranial metric data, followed by cranial nonmetric and dental metric data. Importantly, we demonstrate that combining all four data types together maximizes the neutral genetic signal compared with using them separately, even with a limited number of phenotypic variables. We hypothesize that this reflects a lower level of genetic integration through pleiotropy between, compared to within, the four data types, effectively forming four different modules associated with relatively independent sets of loci. Therefore, we recommend that future craniodental investigations adopt holistic combined data approaches, allowing for more robust inferences about underlying neutral genetic variation.

Midfacial Morphology and Neandertal-Modern Human Interbreeding

Ancient DNA from, Neandertal and modern human fossils, and comparative morphological analyses of them, reveal a complex history of interbreeding between these lineages and the introgression of Neandertal genes into modern human genomes. Despite substantial increases in our knowledge of these events, the timing and geographic location of hybridization events remain unclear. Six measures of facial size and shape, from regional samples of Neandertals and early modern humans, were used in a multivariate exploratory analysis to try to identify regions in which early modern human facial morphology was more similar to that of Neandertals, which might thus represent regions of greater introgression of Neandertal genes. The results of canonical variates analysis and hierarchical cluster analysis suggest important affinities in facial morphology between both Middle and Upper Paleolithic early modern humans of the Near East with Neandertals, highlighting the importance of this region for interbreeding between the two lineages.

Early Alpine occupation backdates westward human migration in Late Glacial Europe

Before the end of the Last Glacial Maximum (LGM, ∼16.5 ka ago)¹ set in motion major shifts in human culture and population structure,² a consistent change in lithic technology, material culture, settlement pattern, and adaptive strategies is recorded in Southern Europe at ∼18-17 ka ago. In this time frame, the landscape of Northeastern Italy changed considerably, and the retreat of glaciers allowed hunter-gatherers to gradually recolonize the Alps.^3-6 Change within this renewed cultural frame (i.e., during the Late Epigravettian phase) is currently associated with migrations favored by warmer climate linked to the Bølling-Allerød onset (14.7 ka ago),^7-11 which replaced earlier genetic lineages with ancestry found in an individual who lived ∼14 ka ago at Riparo Villabruna, Italy, and shared among different contexts (Villabruna Cluster).⁹ Nevertheless, these dynamics and their chronology are still far from being disentangled due to fragmentary evidence for long-distance interactions across Europe.¹² Here, we generate new genomic data from a human mandible uncovered at Riparo Tagliente (Veneto, Italy), which we directly dated to 16,980-16,510 cal BP (2σ). This individual, affected by focal osseous dysplasia, is genetically affine to the Villabruna Cluster. Our results therefore backdate by at least 3 ka the diffusion in Southern Europe of a genetic component linked to Balkan/Anatolian refugia, previously believed to have spread during the later Bølling/Allerød event. In light of the new genetic evidence, this population replacement chronologically coincides with the very emergence of major cultural transitions in Southern and Western Europe.

Gravettian cranial morphology and human group affinities during the European Upper Palaeolithic

Archaeologically defined Upper Palaeolithic (UP, 45,000-10,000 years ago) "cultures" are often used as proxies to designate fossil populations. While recent genomic studies have partly clarified the complex relationship between European UP "cultures" and past population dynamics, they leave open numerous questions regarding the biological characterization of these human groups, especially regarding the Mid-UP period (MUP, 33,000-24,000 years ago), which encompasses a pan-European cultural mosaic (Gravettian) with several regional facies. Here, we analyse a large database of well-dated and well-preserved UP crania, including MUP specimens from South-West France (SWF) and Moravia, using 3D geometric morphometrics to test for human group affinities. Our results show that the Gravettian makers from these two regions form a remarkably phenetically homogeneous sample which is different from, and more homogeneous than, the Late UP sample. Those results are congruent with genomic studies indicating a genetic continuity within the Gravettian manufacturers and a discontinuity marked by the Last Glacial Maximum (LGM). Moreover, our study expands the geographical range of the MUP phenetic continuity to SWF, for which aDNA data are scarce, and clarifies the post-LGM European population structure in SWF, with a possible dual ancestry stemming from different LGM refugia.

Prevalence of cranial trauma in Eurasian Upper Paleolithic humans

OBJECTIVES

This study characterizes patterns of cranial trauma prevalence in a large sample of Upper Paleolithic (UP) fossil specimens (40,000-10,000 BP).

MATERIALS AND METHODS

Our sample comprised 234 individual crania (specimens), representing 1,285 cranial bones (skeletal elements), from 101 Eurasian UP sites. We used generalized linear mixed models (GLMMs) to assess trauma prevalence in relation to age-at-death, sex, anatomical distribution, and between pre- and post-Last Glacial Maximum (LGM) samples, while accounting for skeletal preservation.

RESULTS

Models predicted a mean cranial trauma prevalence of 0.07 (95% CI 0.003-0.19) at the level of skeletal elements, and of 0.26 (95% CI 0.08-0.48) at the level of specimens, each when 76-100% complete. Trauma prevalence increased with skeletal preservation. Across specimen and skeletal element datasets, trauma prevalence tended to be higher for males, and was consistently higher in the old age group. We found no time-specific trauma prevalence patterns for the two sexes or age cohorts when comparing samples from before and after the LGM. Samples showed higher trauma prevalence in the vault than in the face, with vault remains being affected predominantly in males.

DISCUSSION

Cranial trauma prevalence in UP humans falls within the variation described for Mesolithic and Neolithic samples. According to our current dataset, UP males and females were exposed to slightly different injury risks and trauma distributions, potentially due to different activities or behaviors, yet both sexes exhibit more trauma among the old. Environmental stressors associated with climatic changes of the LGM are not reflected in cranial trauma prevalence. To analyze trauma in incomplete skeletal remains we propose GLMMs as an informative alternative to crude frequency calculations.

Morphological evidence for early dog domestication in the European Pleistocene: New evidence from a randomization approach to group differences

The antiquity of the wolf/dog domestication has been recently pushed back in time from the Late Upper Paleolithic (~14,000 years ago) to the Early Upper Paleolithic (EUP; ~36,000 years ago). Some authors questioned this early dog domestication claiming that the putative (EUP) Paleolithic dogs fall within the morphological range of recent wolves. In this study, we reanalyzed a data set of large canid skulls using unbalanced- and balanced-randomized discriminant analyses to assess whether the putative Paleolithic dogs are morphologically unique or whether they represent a subsample of the wolf morpho-population. We evaluated morphological differences between 96 specimens of the 4 a priori reference groups (8 putative Paleolithic dogs, 41 recent northern dogs, 7 Pleistocene wolves, and 40 recent northern wolves) using discriminant analysis based on 5 ln-transformed raw and allometrically size-adjusted cranial measurements. Putative Paleolithic dogs are classified with high accuracies (87.5 and 100.0%, cross-validated) and randomization experiment suggests that these classification rates cannot be exclusively explained by the small and uneven sample sizes of reference groups. It indicates that putative Upper Paleolithic dogs may represent a discrete canid group with morphological signs of domestication (a relatively shorter skull and wider palate and braincase) that distinguish them from sympatric Pleistocene wolves. The present results add evidence to the view that these specimens could represent incipient Paleolithic dogs that were involved in daily activities of European Upper Paleolithic forager groups.

Virtual reconstruction of the Upper Palaeolithic skull from Zlatý Kůň, Czech Republic: Sex assessment and morphological affinity

The incomplete cranium discovered at the Zlatý kůň site in the Bohemian Karst is a rare piece of skeletal evidence of human presence in Central Europe during the Late Glacial period. The relative position of cranial fragments was restored and missing parts of the cranium were virtually reconstructed using mirroring and the Thin-plate splines algorithm. The reconstruction allowed us to collect principal cranial measurements, revise a previous unfounded sex assignment and explore the specimen's morphological affinity. Visual assessment could not reliably provide a sexual diagnosis, as such methods have been developed on modern populations. Using a population-specific approach developed on cranial measurements collected from the literature on reliably sexed European Upper Palaeolithic specimens, linear discriminant analysis confirmed previous assignment to the female sex. However, caution is necessary with regard to the fact that it was assessed from the skull. The Zlatý kůň specimen clearly falls within the range of Upper Palaeolithic craniometric variation. Despite the shift in cranial variation that accompanied the Last Glacial Maximum (LGM), the Zlatý kůň skull exhibits a morphological affinity with the pre-LGM population. Several interpretations are proposed with regard to the complex population processes that occurred after the LGM in Europe.

Dental phenotypic shape variation supports a multiple dispersal model for anatomically modern humans in Southeast Asia

The population history of anatomically modern humans (AMH) in Southeast Asia (SEA) is a highly debated topic. The impact of sea level variations related to the Last Glacial Maximum (LGM) and the Neolithic diffusion on past population dispersals are two key issues. We have investigated competing AMH dispersal hypotheses in SEA through the analysis of dental phenotype shape variation on the basis of very large archaeological samples employing two complementary approaches. We first explored the structure of between- and within-group shape variation of permanent human molar crowns. Second, we undertook a direct test of competing hypotheses through a modeling approach. Our results identify a significant LGM-mediated AMH expansion and a strong biological impact of the spread of Neolithic farmers into SEA during the Holocene. The present work thus favors a "multiple AMH dispersal" hypothesis for the population history of SEA, reconciling phenotypic and recent genomic data.

Morphological change in cranial shape following the transition to agriculture across western Eurasia

The Neolithic transition brought about fundamental social, dietary and behavioural changes in human populations, which, in turn, impacted skeletal morphology. Crania are shaped through diverse genetic, ontogenetic and environmental factors, reflecting various elements of an individual’s life. To determine the transition’s effect on cranial morphology, we investigated its potential impact on the face and vault, two elements potentially responding to different influences. Three datasets from geographically distant regions (Ukraine, Iberia, and the Levant plus Anatolia) were analysed. Craniometric measurements were used to compare the morphology of pre-transition populations with that of agricultural populations. The Neolithic transition corresponds to a statistically significant increase only in cranial breadth of the Ukrainian vaults, while facial morphology shows no consistent transformations, despite expected changes related to the modification of masticatory behaviour. The broadening of Ukrainian vaults may be attributable to dietary and/or social changes. However, the lack of change observed in the other geographical regions and the lack of consistent change in facial morphology are surprising. Although the transition from foraging to farming is a process that took place repeatedly across the globe, different characteristics of transitions seem responsible for idiosyncratic responses in cranial morphology.

Gradual decline in mobility with the adoption of food production in Europe

Increased sedentism during the Holocene has been proposed as a major cause of decreased skeletal robusticity (bone strength relative to body size) in modern humans. When and why declining mobility occurred has profound implications for reconstructing past population history and health, but it has proven difficult to characterize archaeologically. In this study we evaluate temporal trends in relative strength of the upper and lower limb bones in a sample of 1,842 individuals from across Europe extending from the Upper Paleolithic [11,000-33,000 calibrated years (Cal y) B.P.] through the 20th century. A large decline in anteroposterior bending strength of the femur and tibia occurs beginning in the Neolithic (∼ 4,000-7,000 Cal y B.P.) and continues through the Iron/Roman period (∼ 2,000 Cal y B.P.), with no subsequent directional change. Declines in mediolateral bending strength of the lower limb bones and strength of the humerus are much smaller and less consistent. Together these results strongly implicate declining mobility as the specific behavioral factor underlying these changes. Mobility levels first declined at the onset of food production, but the transition to a more sedentary lifestyle was gradual, extending through later agricultural intensification. This finding only partially supports models that tie increased sedentism to a relatively abrupt Neolithic Demographic Transition in Europe. The lack of subsequent change in relative bone strength indicates that increasing mechanization and urbanization had only relatively small effects on skeletal robusticity, suggesting that moderate changes in activity level are not sufficient stimuli for bone deposition or resorption.