Human cranial anatomy and the differential preservation of population history and climate signatures

The relative correspondence of cranial and genetic distances in papionin taxa and the impact of allometric adjustments

The reconstruction of phylogenetic relationships in the primate fossil record is dependent upon a thorough understanding of the phylogenetic utility of craniodental characters. Here, we test three previously proposed hypotheses for the propensity of primate craniomandibular data to exhibit homoplasy, using a study design based on the relative congruence between cranial distance matrices and a consensus genetic distance matrix ("genetic congruence") for papionin taxa: 1) matrices based on cranial regions subjected to less masticatory strain are more genetically congruent than high-strain cranial matrices; 2) matrices based on cranial regions developing earlier in ontogeny are more genetically congruent than matrices based on regions that develop later; and 3) matrices based on cranial regions with greater anatomical/functional complexity are more genetically congruent than matrices based on anatomically simpler regions. Morphological distance matrices based on the shape of 15 different cranial regions, delineated on the basis of previous catarrhine studies, were statistically compared to a matrix of known genetic distances in papionins. Since sexual dimorphism and allometry are known to characterize this clade, several analytical iterations were conducted: 1) mixed-sex, male-only, and female-only analyses and 2) with and without an allometric scaling adjustment. Across all datasets, the chondrocranium matrix was the most consistently correlated with genetic distances, which is also consistent with previous studies of cercopithecoid taxa; however, there was no support for the internal predictions of the three hypotheses tested. Allometric scaling corrections had the largest impact on the genetic congruence of facial shape matrices, a result consistent with previous studies that have described facial homoplasy in papionin taxa. These findings differ from patterns described for hominoid taxa, suggesting that no single predictive criterion can explain phylogenetic utility of cranial datasets across catarrhine primate taxa. Many of the differences in morphological-genetic matrix correlations could result from different levels of phenotypic integration and evolvability in cercopithecoids and hominoids, suggesting that further study of these phenomena in extant primates is warranted.

Examining biological continuity across the late holocene occupation of the Aleutian Islands using cranial morphometrics and quantitative genetic permutation

OBJECTIVES

The number of distinct human migrations into the Aleutian Islands during the Holocene has been a recurrent debate in the anthropological literature. Stemming from Hrdlička's sorting of the prehistoric remains into two distinct populations based on archaeological context and cranial measurements, the human occupation of the Aleutian Islands has long been thought to be the consequence of two distinct human migrations, a Paleo-Aleut migration that provided the initial settlement of the islands, and a Neo-Aleut migration that replaced the original settlers around 1000 BP. This study examines the relationship of the Aleut cranial assemblages in the context of greater Alaskan population variability to assess the evidence for a substantial migration into the Aleutian Islands during the late Holocene.

MATERIALS AND METHODS

A battery of 29 cranial measurements that quantify global cranial shape were analyzed using Euclidean morphometric methods and quantitative genetic permutation methods to examine the plausibility for two distinct Aleut occupations ("Paleo-Aleut" and "Neo-Aleut"), the latter of which is held to share closer phenetic affinities to mainland Alaskan populations than the former. The Aleut skeletal assemblages were arranged according to temporal association, geographic location, and cranial typology, and analyzed within a comparative framework of mainland Alaskan samples using principal coordinates, biological distance and random skewers permutation methods.

RESULTS

Regardless of how the Aleut assemblages are divided, they show greater similarity to each other than to any of the mainland Alaskan assemblages. These findings are consistent across the methodological approaches.

DISCUSSION

The results obtained in this study provide no support for a cranial morphology-based subdivision of the Aleuts into two distinct samples, Hence, there is no evidence for a substantial population migration of so-called Neo-Aleuts, nor for a population replacement event of an extant Paleo-Aleut population by a mainland-affiliated Neo-Aleuts population at or after 1000 BP.

Morphological Covariation between the Maxillary Sinus and Midfacial Skeleton among Sub-Saharan and Circumpolar Modern Humans

OBJECTIVES

Maxillary sinus volume tracks ecogeographic differences in nasal form and may serve as a zone of accommodation for ontogenetic and evolutionary changes in nasal cavity breadth. However, little is known regarding how sinus volume is distributed within the midface. This study investigates morphological covariation between midfacial and sinus shape to better understand structural and functional relationships between the sinus, midface, and nasal cavity.

METHODS

Cranial and sinus models were rendered from CT scans of modern human samples from two disparate climates: sub-Saharan (South Africans [n = 15], West Africans [n = 17]), and circumpolar (Siberian Buriats [n = 18], Alaskan Inuit [n = 20]). Twenty-five 3D coordinate landmarks were placed on the models and subjected to generalized Procrustes analysis. Two-block partial least squares (2B-PLS) analysis was employed to identify patterns of covariation.

RESULTS

The 2B-PLS analysis indicates PLS1 (58.6% total covariation) relates to height and breadth relationships between the midface, nasal cavity, and maxillary sinus. Significant regional differences in PLS1 scores are evident: circumpolar samples possess taller/narrower noses with taller/wider sinuses compared to sub-Saharan samples. Importantly, PLS1 indicates that sinus breadth is not exclusively related to nasal cavity breadth; variation in lateral sinus expansion toward the zygoma represents an important contributing factor. PLS2 (16%) relates to supero-inferior positioning of the sinus within the midface. Allometric trends, while statistically significant, explain only a small portion of these covariation patterns.

CONCLUSIONS

These results suggest that the maxillary sinus serves as a zone of accommodation at the confluence of multiple facial components, potentially minimizing effects of morphological alterations to certain components on adjacent structures. Am J Phys Anthropol 160:483-497, 2016. © 2016 Wiley Periodicals, Inc.

A three-dimensional geometric morphometrics view of the cranial shape variation and population history in the New World

OBJECTIVES

Craniofacial variation in past and present Amerindians has been attributed to the effect of multiple founder events, or to one major migration followed by in situ differentiation and possibly recurrent contacts among Circum-Arctic groups. Our study aims to: (i) detect morphological differences that may indicate several migrations; (ii) test for the presence of genetic isolation; and (iii) test the correlation between shape data and competing settlement hypotheses by taking into account geography, chronology, climate effects, the presence of genetic isolation and recurrent gene flow.

METHODS

We analyzed a large sample of three-dimensional (3D) cranial surface scans (803 specimens) including past and modern groups from America and Australasia. Shape variation was investigated using geometric morphometrics. Differential external gene flow was evaluated by applying genetic concepts to morphometric data (Relethford-Blangero approach). Settlement hypotheses were tested using a matrix correlation approach (Mantel tests).

RESULTS

Our results highlight the strong dichotomy between Circum-Arctic and continental Amerindians as well as the impact of climate adaptation, and possibly recurrent gene flow in the Circum-Arctic area. There is also evidence for the impact of genetic isolation on phenetic variation in Baja California. Several settlement hypotheses are correlated with our data.

CONCLUSIONS

The three approaches used in this study highlight the importance of local processes especially in Baja California, and caution against the use of overly simplistic models when searching for the number of migration events. The results stress the complexity of the settlement of the Americas as well as the mosaic nature of the processes involved in this process. Am. J. Hum. Biol. 28:646-661, 2016. © 2016 Wiley Periodicals, Inc.

Virtual ancestor reconstruction: Revealing the ancestor of modern humans and Neandertals

The timing and geographic origin of the common ancestor of modern humans and Neandertals remain controversial. A poor Pleistocene hominin fossil record and the evolutionary complexities introduced by dispersals and regionalisation of lineages have fuelled taxonomic uncertainty, while new ancient genomic data have raised completely new questions. Here, we use maximum likelihood and 3D geometric morphometric methods to predict possible morphologies of the last common ancestor of modern humans and Neandertals from a simplified, fully resolved phylogeny. We describe the fully rendered 3D shapes of the predicted ancestors of humans and Neandertals, and assess their similarity to individual fossils or populations of fossils of Pleistocene age. Our results support models of an Afro-European ancestral population in the Middle Pleistocene (Homo heidelbergensis sensu lato) and further predict an African origin for this ancestral population.

Early South Americans Cranial Morphological Variation and the Origin of American Biological Diversity

Recent South Americans have been described as presenting high regional cranial morphological diversity when compared to other regions of the world. This high diversity is in accordance with linguistic and some of the molecular data currently available for the continent, but the origin of this diversity has not been satisfactorily explained yet. Here we explore if this high morphological variation was already present among early groups in South America, in order to refine our knowledge about the timing and origins of the modern morphological diversity. Between-group (Fst estimates) and within-group variances (trace of within-group covariance matrix) of the only two early American population samples available to date (Lagoa Santa and Sabana de Bogotá) were estimated based on linear craniometric measurements and compared to modern human cranial series representing six regions of the world, including the Americas. The results show that early Americans present moderate within-group diversity, falling well within the range of modern human groups, despite representing almost three thousand years of human occupation. The between-group variance apportionment is very low between early Americans, but is high among recent South American groups, who show values similar to the ones observed on a global scale. Although limited to only two early South American series, these results suggest that the high morphological diversity of native South Americans was not present among the first human groups arriving in the continent and must have originated during the Middle Holocene, possibly due to the arrival of new morphological diversity coming from Asia during the Holocene.

Developmental Changes in Morphology of the Middle and Posterior External Cranial Base in Modern Homo sapiens

The basicranium has been described as phylogenetically informative, developmentally stable, and minimally affected by external factors and consequently plays an important role in cranial size and shape in subadult humans. Here basicranial variation of subadults from several modern human populations was investigated and the impact of genetic relatedness on basicranial morphological similarities was investigated. Three-dimensional landmark data were digitized from subadult basicrania from seven populations. Published molecular data on short tandem repeats were statistically compared to morphological data from three ontogenetic stages. Basicranial and temporal bone morphology both reflect genetic distances in childhood and adolescence (5–18 years), but not in infancy (<5 years). The occipital bone reflects genetic distances only in adolescence (13–18 years). The sphenoid bone does not reflect genetic distances at any ontogenetic stage but was the most diagnostic region evaluated, resulting in high rates of correct classification among populations. These results suggest that the ontogenetic processes driving basicranial development are complex and cannot be succinctly summarized across populations or basicranial regions. However, the fact that certain regions reflect genetic distances suggests that the morphology of these regions may be useful in reconstructing population history in specimens for which direct DNA evidence is unavailable, such as archaeological sites.

Differential Scaling Patterns in Maxillary Sinus Volume and Nasal Cavity Breadth Among Modern Humans

Among modern humans, nasal cavity size and shape reflect its vital role in air conditioning processes. The ability for the nasal cavity to augment its shape, particularly in inferior breadth, likely relates to the surrounding maxillary sinuses acting as zones of accommodation. However, much is still unknown regarding how nasal and sinus morphology relate to each other and to overall craniofacial form, particularly across diverse populations with varying respiratory demands. As such, this study uses computed tomographic (CT) scans of modern human crania (N = 171) from nine different localities to investigate ecogeographic differences in (1) the interaction between maxillary sinus volume (MSV) and nasal cavity breadth (NCB) and (2) scaling patterns of MSV and NCB in relation to craniofacial size. Reduced major axis (RMA) regression reveals that all samples exhibit an inverse relationship between MSV and NCB, but statistical significance and the strength of that relationship is sample dependent. Individuals from cold-dry climates have larger MSVs with narrower NCBs, while smaller MSVs are associated with wider NCBs in hot-humid climates. MSV and NCB each scale with positive allometry relative to overall craniofacial size. However, sample differences are evident in the both the interaction between MSV and NCB, as well as their correlation with craniofacial size. While these results provide further support that the maxillary sinus and nasal cavity are integrated among populations from opposite ends of the climatic spectrum, additional epigenetic factors are needed to explain variation of these structures among populations from more intermediate climates.

Heritability of the Human Craniofacial Complex

Quantifying normal variation and the genetic underpinnings of anatomical structures is one of the main goals of modern morphological studies. However, the extent of genetic contributions to normal variation in craniofacial morphology in humans is still unclear. The current study addresses this gap by investigating the genetic underpinnings of normal craniofacial morphology. The sample under investigation consists of 75 linear and angular measurements spanning the entire craniofacial complex, recorded from lateral cephalographs of 1,379 participants in the Fels Longitudinal Study. Heritabilities for each trait were estimated using SOLAR, a maximum-likelihood variance components approach utilizing all pedigree information for parameter estimation. Trait means and mean effects of the covariates age, sex, age(2) , sex × age, and sex × age(2) were simultaneously estimated in the analytic models. All traits of the craniofacial complex were significantly heritable. Heritability estimates ranged from 0.10 to 0.60, with the majority being moderate. It is important to note that we found similar ranges of heritability occurring across the different functional/developmental components of the craniofacial complex, the splanchnocranium, the basicranium, and the neurocranium. This suggests that traits from different regions of the craniofacial complex are of comparable utility for the purposes of population history and phylogeny reconstruction. At the same time, this genetic influence on craniofacial morphology signals a caution to researchers of nongenetic studies to consider the implications of this finding when selecting samples for study given their project design and goals.

The cranial morphology of the Botocudo Indians, Brazil

The Botocudo Indians were hunter-gatherer groups that occupied the East-Central regions of Brazil decimated during the colonial period in the country. During the 19th century, craniometric studies suggested that the Botocudo resembled more the Paleoamerican population of Lagoa Santa than typical Native Americans groups. These results suggest that the Botocudo Indians might represent a population that retained the biological characteristics of early groups of the continent, remaining largely isolated from groups that gave origin to the modern Native South American variation. Moreover, recently, some of the Botocudo remains have been shown to have mitochondrial and autosomal DNA lineages currently found in Polynesian populations. Here, we explore the morphological affinities of Botocudo skulls within a worldwide context. Distinct multivariate analyses based on 32 craniometric variables show that 1) the two individuals with Polynesian DNA sequences have morphological characteristics that fall within the Polynesian and Botocudo variation, making their assignation as Native American specimens problematic, and 2) there are high morphological affinities between Botocudo, Early Americans, and the Polynesian series of Easter Island, which support the early observations that the Botocudo can be seen as retaining the Paleoamerican morphology, particularly when the neurocranium is considered. Although these results do not elucidate the origin of the Polynesian DNA lineages among the Botocudo, they support the hypothesis that the Botocudo represent a case of late survival of ancient Paleoamerican populations, retaining the morphological characteristics of ancestral Late Pleistocene populations from Asia.

Analysis of the human osseous nasal shape--population differences and sexual dimorphism

OBJECTIVES

In this study, the shape of the outer osseous nose in a German and a Chinese sample is analyzed using a dense set of semi-landmarks. Shape differences related to population and sex as well as directional and fluctuating asymmetry were statistically evaluated and also visualized.

MATERIALS AND METHODS

Shape differences in the bony nose were investigated between a large sample of CT scans of German (140 ♀, 127 ♂) and Chinese (135 ♀, 132 ♂) crania. We used semi-automatic methods to represent the shape of this region as a dense point-cloud, consisting of 370 three-dimensional bilateral coordinates. Both the symmetric and asymmetric modes of shape variation were addressed.

RESULTS

Strong differences in nasal shape were found between the two populations, while sex was found to play a minor role in explaining the observed shape variation. The expression of sexual dimorphism was similar in both populations. Differences attributed to population affinity and to sexual dimorphism were both found to affect the shape of the ossa nasalia and the projection of the spina nasalis. The correlation with population/sex was weak for directional asymmetry, but strong for fluctuating asymmetry. The nasal region is more asymmetric in Germans than in Chinese, with males displaying more asymmetry than females in both populations.

DISCUSSION

While the bony nose is well suited for predicting population affinity, regarding the populations under investigation, its value for sexing unknown individuals is rather moderate. The similar expression of sexual dimorphism in those otherwise very dissimilar populations indicates common factors responsible for these differences.

Late Pleistocene age and archaeological context for the hominin calvaria from GvJm-22 (Lukenya Hill, Kenya)

Kenya National Museums Lukenya Hill Hominid 1 (KNM-LH 1) is a Homo sapiens partial calvaria from site GvJm-22 at Lukenya Hill, Kenya, associated with Later Stone Age (LSA) archaeological deposits. KNM-LH 1 is securely dated to the Late Pleistocene, and samples a time and region important for understanding the origins of modern human diversity. A revised chronology based on 26 accelerator mass spectrometry radiocarbon dates on ostrich eggshells indicates an age range of 23,576-22,887 y B.P. for KNM-LH 1, confirming prior attribution to the Last Glacial Maximum. Additional dates extend the maximum age for archaeological deposits at GvJm-22 to >46,000 y B.P. (>46 kya). These dates are consistent with new analyses identifying both Middle Stone Age and LSA lithic technologies at the site, making GvJm-22 a rare eastern African record of major human behavioral shifts during the Late Pleistocene. Comparative morphometric analyses of the KNM-LH 1 cranium document the temporal and spatial complexity of early modern human morphological variability. Features of cranial shape distinguish KNM-LH 1 and other Middle and Late Pleistocene African fossils from crania of recent Africans and samples from Holocene LSA and European Upper Paleolithic sites.

Geometric morphometric study of the skull shape diversification in Sciuridae (Mammalia, Rodentia)

It is generally accepted that the high phenotypic diversity of mammals is a combined result of developmental constraint and ecological adaptation, although the influence of these endogenous and exogenous factors varies in different mammal groups. The rodent family Sciuridae represents an ideal candidate for examining phenotypic diversity in relation to phylogeny and ecological adaptations. In the present study, we investigate the effects of phylogeny and lifestyle on the skull shape in different species of Sciuridae by applying geometric morphometric methods. In addition, we investigate the importance of allometry on sciurid skull shape, because results from geometric morphometrics sometimes dispute those of traditional morphometry. Here, we identify significant associations between patristic distances obtained from molecular phylogeny and shape distances in all 3 views of the cranium and the lateral view of the mandible. Multivariate regression demonstrates that shape differences among lifestyle categories are substantial, especially in the dorsal and ventral structures after the influence of phylogeny is taken into account. Allometry plays an important role in the shape variation, although its importance on different skull structures varies. Our results indicate that complex structures of this highly diverse mammal group, which occupies different niches, are affected by ecological factors and developmental constraint.

The contribution of subsistence to global human cranial variation

Diet-related cranial variation in modern humans is well documented on a regional scale, with ample examples of cranial changes related to the agricultural transition. However, the influence of subsistence strategy on global cranial variation is less clear, having been confirmed only for the mandible, and dietary effects beyond agriculture are often neglected. Here we identify global patterns of subsistence-related human cranial shape variation. We analysed a worldwide sample of 15 populations (n = 255) with known subsistence strategies using 3-D landmark datasets designed to capture the shape of different units of the cranium. Results show significant correlations between global cranial shape and diet, especially for temporalis muscle shape and general cranial shape. Importantly, the differences between populations with either a plant- or an animal-based diet are more pronounced than those between agriculturalists and hunter-gatherers, suggesting that the influence of diet as driver of cranial variation is not limited to Holocene transitions to agricultural subsistence. Dental arch shape did not correlate with subsistence pattern, possibly indicating the high plasticity of this region of the face in relation to age, disease and individual use of the dentition. Our results highlight the importance of subsistence strategy as one of the factors underlying the evolution of human geographic cranial variation.

Environmental temperature impact on bone and cartilage growth

Environmental temperature can have a surprising impact on extremity growth in homeotherms, but the underlying mechanisms have remained elusive for over a century. Limbs of animals raised at warm ambient temperature are significantly and permanently longer than those of littermates housed at cooler temperature. These remarkably consistent lab results closely resemble the ecogeographical tenet described by Allen's "extremity size rule," that appendage length correlates with temperature and latitude. This phenotypic growth plasticity could have adaptive significance for thermal physiology. Shortened extremities help retain body heat in cold environments by decreasing surface area for potential heat loss. Homeotherms have evolved complex mechanisms to maintain tightly regulated internal temperatures in challenging environments, including "facultative extremity heterothermy" in which limb temperatures can parallel ambient. Environmental modulation of tissue temperature can have direct and immediate consequences on cell proliferation, metabolism, matrix production, and mineralization in cartilage. Temperature can also indirectly influence cartilage growth by modulating circulating levels and delivery routes of essential hormones and paracrine regulators. Using an integrated approach, this article synthesizes classic studies with new data that shed light on the basis and significance of this enigmatic growth phenomenon and its relevance for treating human bone elongation disorders. Discussion centers on the vasculature as a gateway to understanding the complex interconnection between direct (local) and indirect (systemic) mechanisms of temperature-enhanced bone lengthening. Recent advances in imaging modalities that enable the dynamic study of cartilage growth plates in vivo will be key to elucidating fundamental physiological mechanisms of long bone growth regulation.

Ecogeography, genetics, and the evolution of human body form

Genetic resemblances among groups are non-randomly distributed in humans. This population structure may influence the correlations between traits and environmental drivers of natural selection thus complicating the interpretation of the fossil record when modern human variation is used as a referential model. In this paper, we examine the effects of population structure and natural selection on postcranial traits that reflect body size and shape with application to the more general issue of how climate - using latitude as a proxy - has influenced hominin morphological variation. We compare models that include terms reflecting population structure, ascertained from globally distributed microsatellite data, and latitude on postcranial phenotypes derived from skeletal dimensions taken from a large global sample of modern humans. We find that models with a population structure term fit better than a model of natural selection along a latitudinal cline in all cases. A model including both latitude and population structure terms is a good fit to distal limb element lengths and bi-iliac breadth, indicating that multiple evolutionary forces shaped these morphologies. In contrast, a model that included only a population structure term best explained femoral head diameter and the crural index. The results demonstrate that population structure is an important part of human postcranial variation, and that clinally distributed natural selection is not sufficient to explain among-group differentiation. The distribution of human body form is strongly influenced by the contingencies of modern human origins, which calls for new ways to approach problems in the evolution of human variation, past and present.

Ancestry estimation in South Africa using craniometrics and geometric morphometrics

Population history and positive assortative mating directs gene flow in such a way that biological differences are recognized among groups. In turn, forensic anthropologists quantify biological differences to estimate ancestry. Some anthropologists argue that highly admixed population groups, such as South African coloureds, cannot achieve acceptable accuracies because within group variance is too large. Whereas ancestry estimation in South Africa has been limited to craniometric data from South African blacks and whites, the current study integrates craniometric and geometric morphometric data from the three largest South African groups. Crania from 377 South African individuals (black=158, white=112, and coloured=107) comprised the sample. Standard measurements were collected and the coordinate data were subjected to Generalized Procrustes Analysis (GPA), which resulted in size-free shape variables (ProCoords). A principal component analysis was used to combine the shape variation captured in the ProCoords (ProCoords PC). Linear discriminant analysis (LDA), using equal priors, stepwise variable selection and leave-one-out cross-validation, was conducted on the ProCoords, the ProCoords PCs, and the traditional craniometric data. The LDA using 18 stepwise selected ProCoords resulted in the highest cross-validated accuracy (89%). Utilization of geometric morphometric data emphasized that the relative location of cranial landmarks was more discriminating than simple linear distances. Regardless of high levels of genetic admixture, South African coloureds are a homogeneous group and morphologically distinct from other contemporaneous South African populations. Furthermore, the present study demonstrated a correspondence between peer-reported race and morphological differences in the crania of black, white, and coloured South Africans.

Admixture in Latin America: geographic structure, phenotypic diversity and self-perception of ancestry based on 7,342 individuals

The current genetic makeup of Latin America has been shaped by a history of extensive admixture between Africans, Europeans and Native Americans, a process taking place within the context of extensive geographic and social stratification. We estimated individual ancestry proportions in a sample of 7,342 subjects ascertained in five countries (Brazil, Chile, Colombia, México and Perú). These individuals were also characterized for a range of physical appearance traits and for self-perception of ancestry. The geographic distribution of admixture proportions in this sample reveals extensive population structure, illustrating the continuing impact of demographic history on the genetic diversity of Latin America. Significant ancestry effects were detected for most phenotypes studied. However, ancestry generally explains only a modest proportion of total phenotypic variation. Genetically estimated and self-perceived ancestry correlate significantly, but certain physical attributes have a strong impact on self-perception and bias self-perception of ancestry relative to genetically estimated ancestry.

Latin America has a history of extensive mixing between Native Americans and people arriving from Europe and Africa. As a result, individuals in the region have a highly heterogeneous genetic background and show great variation in physical appearance. Latin America offers an excellent opportunity to examine the genetic basis of the differentiation in physical appearance between Africans, Europeans and Native Americans. The region is also an advantageous setting in which to examine the interplay of genetic, physical and social factors in relation to ethnic/racial self-perception. Here we present the most extensive analysis of genetic ancestry, physical diversity and self-perception of ancestry yet conducted in Latin America. We find significant geographic variation in ancestry across the region, this variation being consistent with demographic history and census information. We show that genetic ancestry impacts many aspects of physical appearance. We observe that self-perception is highly influenced by physical appearance, and that variation in physical appearance biases self-perception of ancestry relative to genetically estimated ancestry.

Variation and signatures of selection on the human face

There has been much debate about why humans throughout the world differ in facial form. Previous studies of human skull morphology found levels of among-population differentiation that were comparable to those of neutral genetic markers, suggesting that genetic drift (neutral processes) played an important role in influencing facial differentiation. However, variation in soft-tissue morphology has not been studied in detail. In this study, we analyzed high-resolution 3D images of soft-tissue facial form in four Eurasian populations: Han Chinese, Tibetans, Uyghur and Europeans. A novel method was used to establish a high-density alignment across all of the faces, allowing facial diversity to be examined at an unprecedented resolution. These data exhibit signatures of population structure and history. However, among-population differentiation was higher for soft-tissue facial form than for genome-wide genetic loci, and high-resolution analyses reveal that the nose, brow area and cheekbones exhibit particularly strong signals of differentiation (Qst estimates: 0.3-0.8) between Europeans and Han Chinese. Our results suggest that local adaptation and/or sexual selection have been important in shaping human soft-tissue facial morphology.

Iwo Eleru's place among Late Pleistocene and Early Holocene populations of North and East Africa

The Iwo Eleru site in Nigeria preserves the only terminal Pleistocene fossil from tropical West Africa. The peoples of this region contributed to significant population movements throughout the continent during the Holocene. As such, characterizing the phenotype of Late Pleistocene West African populations is critical for disentangling the evolutionary signatures of a highly complex African population history and structure. Previous research approached the calvaria's morphology from a paleoanthropological perspective, noting its mosaic of archaic and modern neurocranial features and distinctiveness from Pleistocene fossil taxa and contemporary modern human samples. In this paper, I compare Iwo Eleru with contemporary Late Pleistocene Africans and also consider the specimen's affinities with Holocene populations of the central and western Sahara, Nile Valley, and East Africa. Craniometric data were recorded for 22 neurocranial dimensions and subjected to principal components analysis and Mahalanobis distance estimation. Multidimensional scaling of distances indicated that Iwo Eleru fell outside the observed range of variation of other terminal Pleistocene supra-equatorial African populations, confirming previous results that documented its divergence from Neanderthals, Upper Paleolithic Europeans, and modern Africans. The calvaria was also distinct from Holocene Saharan, Nile Valley, and East African populations, which suggests limited West African input into the Sahara during the African Humid Period. Results presented here bolster previous research that suggested Iwo Eleru's anatomy reflected either admixture with archaic humans or the long-term survival of populations with more archaic neurocranial anatomy until the end of the Pleistocene.

Phenotypic variation in infants, not adults, reflects genotypic variation among chimpanzees and bonobos

Studies comparing phenotypic variation with neutral genetic variation in modern humans have shown that genetic drift is a main factor of evolutionary diversification among populations. The genetic population history of our closest living relatives, the chimpanzees and bonobos, is now equally well documented, but phenotypic variation among these taxa remains relatively unexplored, and phenotype-genotype correlations are not yet documented. Also, while the adult phenotype is typically used as a reference, it remains to be investigated how phenotype-genotye correlations change during development. Here we address these questions by analyzing phenotypic evolutionary and developmental diversification in the species and subspecies of the genus Pan. Our analyses focus on the morphology of the femoral diaphysis, which represents a functionally constrained element of the locomotor system. Results show that during infancy phenotypic distances between taxa are largely congruent with non-coding (neutral) genotypic distances. Later during ontogeny, however, phenotypic distances deviate from genotypic distances, mainly as an effect of heterochronic shifts between taxon-specific developmental programs. Early phenotypic differences between Pan taxa are thus likely brought about by genetic drift while late differences reflect taxon-specific adaptations.

The interaction of neutral evolutionary processes with climatically-driven adaptive changes in the 3D shape of the human os coxae

Differences in the breadth of the pelvis among modern human populations and among extinct hominin species have often been interpreted in the light of thermoregulatory adaptation, whereby a larger pelvic girdle would help preserve body temperature in cold environments while a narrower pelvis would help dissipate heat in tropical climates. There is, however, a theoretical problem in interpreting a pattern of variation as evidence of selection without first accounting for the effects of neutral evolutionary processes (i.e., mutation, genetic drift and migration). Here, we analyse 3D configurations of 27 landmarks on the os coxae of 1494 modern human individuals representing 30 male and 23 female populations from five continents and a range of climatic conditions. We test for the effects of climate on the size and shape of the pelvic bone, while explicitly accounting for population history (i.e., geographically-mediated gene flow and genetic drift). We find that neutral processes account for a substantial proportion of shape variance in the human os coxae in both sexes. Beyond the neutral pattern due to population history, temperature is a significant predictor of shape and size variation in the os coxae, at least in males. The effect of climate on the shape of the pelvic bone, however, is comparatively limited, explaining only a small percentage of shape variation in males and females. In accordance with previous hypotheses, the size of the os coxae tends to increase with decreasing temperature, although the significance of the association is reduced when population history is taken into account. In conclusion, the shape and size of the human os coxae reflect both neutral evolutionary processes and climatically-driven adaptive changes. Neutral processes have a substantial effect on pelvic variation, suggesting such factors will need to be taken into account in future studies of human and fossil hominin coxal variation.

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

The Last Glacial Maximum (LGM) represents the most significant climatic event since the emergence of anatomically modern humans (AMH). In Europe, the LGM may have played a role in changing morphological features as a result of adaptive and stochastic processes. We use craniometric data to examine morphological diversity in pre- and post-LGM specimens. Craniometric variation is assessed across four periods--pre-LGM, late glacial, Early Holocene and Middle Holocene--using a large, well-dated, data set. Our results show significant differences across the four periods, using a MANOVA on size-adjusted cranial measurements. A discriminant function analysis shows separation between pre-LGM and later groups. Analyses repeated on a subsample, controlled for time and location, yield similar results. The results are largely influenced by facial measurements and are most consistent with neutral demographic processes. These findings suggest that the LGM had a major impact on AMH populations in Europe prior to the Neolithic.

Genomic and cranial phenotype data support multiple modern human dispersals from Africa and a southern route into Asia

Despite broad consensus on Africa as the main place of origin for anatomically modern humans, their dispersal pattern out of the continent continues to be intensely debated. In extant human populations, the observation of decreasing genetic and phenotypic diversity at increasing distances from sub-Saharan Africa has been interpreted as evidence for a single dispersal, accompanied by a series of founder effects. In such a scenario, modern human genetic and phenotypic variation was primarily generated through successive population bottlenecks and drift during a rapid worldwide expansion out of Africa in the Late Pleistocene. However, recent genetic studies, as well as accumulating archaeological and paleoanthropological evidence, challenge this parsimonious model. They suggest instead a "southern route" dispersal into Asia as early as the late Middle Pleistocene, followed by a separate dispersal into northern Eurasia. Here we test these competing out-of-Africa scenarios by modeling hypothetical geographical migration routes and assessing their correlation with neutral population differentiation, as measured by genetic polymorphisms and cranial shape variables of modern human populations from Africa and Asia. We show that both lines of evidence support a multiple-dispersals model in which Australo-Melanesian populations are relatively isolated descendants of an early dispersal, whereas other Asian populations are descended from, or highly admixed with, members of a subsequent migration event.

Detecting genetic association of common human facial morphological variation using high density 3D image registration

Human facial morphology is a combination of many complex traits. Little is known about the genetic basis of common facial morphological variation. Existing association studies have largely used simple landmark-distances as surrogates for the complex morphological phenotypes of the face. However, this can result in decreased statistical power and unclear inference of shape changes. In this study, we applied a new image registration approach that automatically identified the salient landmarks and aligned the sample faces using high density pixel points. Based on this high density registration, three different phenotype data schemes were used to test the association between the common facial morphological variation and 10 candidate SNPs, and their performances were compared. The first scheme used traditional landmark-distances; the second relied on the geometric analysis of 15 landmarks and the third used geometric analysis of a dense registration of ∼30,000 3D points. We found that the two geometric approaches were highly consistent in their detection of morphological changes. The geometric method using dense registration further demonstrated superiority in the fine inference of shape changes and 3D face modeling. Several candidate SNPs showed potential associations with different facial features. In particular, one SNP, a known risk factor of non-syndromic cleft lips/palates, rs642961 in the IRF6 gene, was validated to strongly predict normal lip shape variation in female Han Chinese. This study further demonstrated that dense face registration may substantially improve the detection and characterization of genetic association in common facial variation.

Heritability of human facial appearance is an intriguing question to the general public and researchers. Although it is known that some facial features are highly heritable, the exact genetic basis is unknown. Previous studies used simple linear measurements such as landmark distances, to evaluate the facial shape variation. Such approaches, although easy to carry out, may lack statistical power and miss complex morphological changes. In this study, we utilized a new 3D face registration method that enables subtle differences to be detected at high resolution 3D images. Based on this, we tried to test and characterize the associations of 10 candidate genetic variants to common facial morphological variations. Different types of phenotype data were extracted and compared in the association tests. Our results show that geometry based data performed better than simple distance based data. Furthermore, high density geometric data outstood the others in capturing small shape changes and modeling the 3D face visualization. Interestingly, a genetic variant from IRF6 gene, which is also a well-known risk factor of non-syndrome cleft lip, was found to strongly predispose the mouth shape in Han Chinese females.

Functional and morphological correlates of mandibular symphyseal form in a living human sample

Variation in recent human mandibular form is often thought to reflect differences in masticatory behavior associated with variation in food preparation and subsistence strategies. Nevertheless, while mandibular variation in some human comparisons appear to reflect differences in functional loading, other comparisons indicate that this relationship is not universal. This suggests that morphological variation in the mandible is influenced by other factors that may obscure the effects of loading on mandibular form. It is likely that highly strained mandibular regions, including the corpus, are influenced by well-established patterns of lower facial skeletal integration. As such, it is unclear to what degree mandibular form reflects localized stresses incurred during mastication vs. a larger set of correlated features that may influence bone distribution patterns. In this study, we examine the relationship between mandibular symphyseal bone distribution (i.e., second moments of area, cortical bone area) and masticatory force production (i.e., in vivo maximal bite force magnitude and estimated symphyseal bending forces) along with lower facial shape variation in a sample of n = 20 living human male subjects. Our results indicate that while some aspects of symphyseal form (e.g., wishboning resistance) are significantly correlated with estimates of symphyseal bending force magnitude, others (i.e., vertical bending resistance) are more closely tied to variation in lower facial shape. This suggests that while the symphysis reflects variation in some variables related to functional loading, the complex and multifactorial influences on symphyseal form underscores the importance of exercising caution when inferring function from the mandible especially in narrow taxonomic comparisons.

Cranial airways and the integration between the inner and outer facial skeleton in humans

The cranial airways are in the center of the human face. Therefore variation in the size and shape of these central craniofacial structures could have important consequences for the surrounding midfacial morphology during development and evolution. Yet such interactions are unclear because one school of thought, based on experimental and developmental evidence, suggests a relative independence (modularity) of these two facial compartments, whereas another one assumes tight morphological integration. This study uses geometric morphometrics of modern humans (N = 263) and 40 three-dimensional-landmarks of the skeletal nasopharynx and nasal cavity and outer midfacial skeleton to analyze these questions in terms of modularity. The sizes of all facial compartments were all strongly correlated. Shape integration was high between the cranial airways and the outer midfacial skeleton and between the latter and the anterior airway openings (skeletal regions close to and including piriform aperture). However, no shape integration was detected between outer midface and posterior airway openings (nasopharynx and choanae). Similarly, no integration was detected between posterior and anterior airway openings. This may reflect functional modularization of nasal cavity compartments related to respiratory physiology and differential developmental interactions with the face. Airway size likely relates to the energetics of the organism, whereas airways shape might be more indicative of respiratory physiology and climate. Although this hypothesis should be addressed in future steps, here we suggest that selection on morphofunctional characteristics of the cranial airways could have cascading effects for the variation, development, and evolution of the human face.

Evaluating developmental shape changes in Homo antecessor subadult facial morphology

The fossil ATD6-69 from Atapuerca, Spain, dated to ca. 900 ka (thousands of years ago) has been suggested to mark the earliest appearance of modern human facial features. However, this specimen is a subadult and the interpretation of its morphology remains controversial, because it is unclear how developmental shape changes would affect the features that link ATD6-69 to modern humans. Here we analyze ATD6-69 in an evolutionary and developmental context. Our modern human sample comprises cross-sectional growth series from four populations. The fossil sample covers human specimens from the Pleistocene to the Upper Paleolithic, and includes several subadult Early Pleistocene humans and Neanderthals. We digitized landmarks and semilandmarks on surface and CT scans and analyzed the Procrustes shape coordinates using multivariate statistics. Ontogenetic allometric trajectories and developmental simulations were employed in order to identify growth patterns and to visualize potential adult shapes of ATD6-69. We show that facial differences between modern and archaic humans are not exclusively allometric. We find that while postnatal growth further accentuates the differences in facial features between Neanderthals and modern humans, those features that have been suggested to link ATD6-69's morphology to modern humans would not have been significantly altered in the course of subsequent development. In particular, the infraorbital depression on this specimen would have persisted into adulthood. However, many of the facial features that ATD6-69 shares with modern humans can be considered to be part of a generalized pattern of facial architecture. Our results present a complex picture regarding the polarity of facial features and demonstrate that some modern human-like facial morphology is intermittently present in Middle Pleistocene humans. We suggest that some of the facial features that characterize recent modern humans may have developed multiple times in human evolution.

Skull and limb morphology differentially track population history and environmental factors in the transition to agriculture in Europe

The Neolithic transition in Europe was a complex mosaic spatio-temporal process, involving both demic diffusion from the Near East and the cultural adoption of farming practices by indigenous hunter-gatherers. Previous analyses of Mesolithic hunter-gatherers and Early Neolithic farmers suggest that cranial shape variation preserves the population history signature of the Neolithic transition. However, the extent to which these same demographic processes are discernible in the postcranium is poorly understood. Here, for the first time, crania and postcranial elements from the same 11 prehistoric populations are analysed together in an internally consistent theoretical and methodological framework. Results show that while cranial shape reflects the population history differences between Mesolithic and Neolithic lineages, relative limb dimensions exhibit significant congruence with environmental variables such as latitude and temperature, even after controlling for geography and time. Also, overall limb size is found to be consistently larger in hunter-gatherers than farmers, suggesting a reduction in size related to factors other than thermoregulatory adaptation. Therefore, our results suggest that relative limb dimensions are not tracking the same demographic population history as the cranium, and point to the strong influence of climatic, dietary and behavioural factors in determining limb morphology, irrespective of underlying neutral demographic processes.