Quantitative genetics of skeletal nonmetric traits in the rhesus macaques on Cayo Santiago. II. Phenotypic, genetic, and environmental correlations between traits

Exploring Kinship Within a Late Meroitic to Medieval Cemetery in Sudan Using k-Modes Clustering

OBJECTIVES

Intracemetery patterns at the Qinifab School site, a late Meroitic to Christian period cemetery (c. 250-1400 CE) between the fourth and fifth cataracts of the Nile River in Sudan, are modeled to explore the role of kinship in cemetery organization during significant sociopolitical and religious shifts.

MATERIALS AND METHODS

Sixty-eight cranial and 36 dental nonmetric traits were examined among 67 adults. A machine-learning clustering method (k-modes clustering algorithm), new to biological anthropology, was used to detect patterning among biological, demographic, and temporal data and validated with inter-individual Mahalanobis distances and mixed data principal components analysis using only biological data.

RESULTS

When three validated clusters were examined, a pattern emerged that aligned with the archaeological context of the cemetery when time period is considered. One cluster concentrated at the oldest end of the site appeared to be a founding group, with the other two groups equally comprising later periods. Some mortuary patterning by sex was visually identified; a few late to Post-Meroitic males grouped together in one cluster.

DISCUSSION

Individuals from Clusters 1 and 3 were buried near each other infrequently, which suggests two different kin groups may have intermarried with the third, or one group held a status that allowed identification with the other clusters. Kinship was not a determinant for males with archery items or individuals with incisor avulsion, as they were found in all three clusters. The Qinifab School cemetery reflects temporal organization, a multilocal residence pattern, and a lack of sex stratification for the medieval period.

Heritability in the Rhesus Macaque (Macaca mulatta) Vertebral Column

The vertebral column plays a central role in primate locomotion and positional behavior. Understanding its evolution, therefore, has the potential to clarify evolutionary processes that have occurred in the primate lineage as well as the specific behaviors of extinct primates. However, to understand primate vertebral anatomy, it is important to determine how much of this anatomy is heritable and how much develops as a response to environmental factors during life. We estimated heritability for vertebral counts as well as typical cervical, thoracic, and lumbar elements from 210 individuals from the pedigreed Cayo Santiago Macaca mulatta skeletal collection. We found moderate heritability of vertebral counts (h2 = 0.216-0.326), but with strong heritability of the type of variation (e.g., a tendency to meristic or homeotic change) in the vertebral count (h2 = 0.599), suggesting a possible explanation for high variability in vertebral numbers among the hominoids in particular. The moderate heritability of vertebral count also suggests that vertebral count is an unsuitable metric for estimating the ancestral state for some taxa. We found strong heritability in the morphology of cervical and upper lumbar zygapophyseal facets (h2 = 0.548-0.550) and the thoracic spinous processes (h2 = 0.609-0.761), including high heritability of the spinous process angle in the upper thoracic and upper lumbar elements (h2 = 0.649-0.752). We suggest these are related to maintaining stability in the cervical and lumbar regions, and reducing motion in the thoracic region, respectively. We propose that spinous processes may contain greater phylogenetic information, whereas transverse processes may contain greater information of function 'in life'. We also found important size effects, suggesting that size is the most heritable component of overall form and largely responsible for intertrait differences. This suggests that it is inappropriate to indiscriminately remove size effects from morphological comparisons.

Coming to the Caribbean: Eighty-five years of rhesus macaques (Macaca mulatta) at Cayo Santiago-A rare nonhuman primate model for the studies of adaptation, diseases, genetics, natural disasters, and resilience

The Cayo Santiago rhesus macaque colony represents one of the most important nonhuman primate resources since their introduction to the Caribbean area in 1938. The 85 years of continuing existence along with the comprehensive database of the rhesus colony and the derived skeletal collections have provided and will continue to provide a powerful tool to test hypotheses about adaptive and evolutionary mechanisms in both biology and medicine.

History of Health at Cayo Santiago-An Investigation of Environmental and Genetic Influences on the Skeletal Remains of the Introduced Rhesus Macaque (Macaca mulatta) Colony

The Cayo Santiago rhesus macaque colony is a renowned primate population that has experienced significant natural and anthropogenic ecological variation in their 85-year history. Demographic and familial information is also tracked and collated for the majority of monkeys. Thus, the health history of rhesus macaques at Cayo Santiago should reflect the impacts of both environmental and genetic factors. In this study, we utilized a sample of skeletal remains comprised of 2787 individuals (1571 females, 1091 males), born between 1938 and 2017 from the derived skeletal collection of the primate colony to assess variation in survivorship, pathology, bone mineral density (BMD), and dental eruption status, in the context of hurricane impacts, nutritional fluctuations, and matriline genealogy. Results demonstrated that rhesus macaques at Cayo Santiago exhibit a range of skeletal pathologies that encompass biomedical and archaeological significance, multiple etiologies, severities, locations, and types, in addition to a secular trend of declining BMD that is hypothesized to reflect decreasing physical activity levels under increasing population densities. Specifically, hurricane impacts were found to increase the rate of systemic disease, decrease BMD in young adults, and delay eruption of the primary dentition. Certain matrilines exhibited heightened levels of systemic disease at early ages while others exhibited greater rates of congenital disease. Early-life adversity, through the experience of major hurricanes, may enhance inflammatory pathways, heightening the risk of disease and accelerating the aging process leading to reduced BMD. Such impacts may underly greater levels of observed infection post-hurricane through intensification of pathogen transmission and disease rates brought on by hurricane-adaptive social strategies that favor closer proximity. Familial susceptibility to disease indicates heritable host genetic factors are likely influencing disease patterning in the population. A cluster of congenital diseases may most convincingly illustrate this, or alternatively reflects low levels of genetic diversity in the population.

Human population dynamics in Upper Paleolithic Europe inferred from fossil dental phenotypes

Despite extensive archaeological research, our knowledge of the human population history of Upper Paleolithic Europe remains limited, primarily due to the scarce availability and poor molecular preservation of fossil remains. As teeth dominate the fossil record and preserve genetic signatures in their morphology, we compiled a large dataset of 450 dentitions dating between ~47 and 7 thousand years ago (ka), outnumbering existing skeletal and paleogenetic datasets. We tested a range of competing demographic scenarios using a coalescent-based machine learning Approximate Bayesian Computation (ABC) framework that we modified for use with phenotypic data. Mostly in agreement with but also challenging some of the hitherto available evidence, we identified a population turnover in western Europe at ~28 ka, isolates in western and eastern refugia between ~28 and 14.7 ka, and bottlenecks during the Last Glacial Maximum. Methodologically, this study marks the pioneering application of ABC to skeletal phenotypes, paving the way for exciting future research avenues.

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.

Biological distances between Late Holocene hunter‐gatherers from the lower Paraná wetland (Argentina)

INTRODUCTION

This study assesses the biological affinity between Late Holocene hunter-gatherers that inhabited the region of the lower Paraná wetland (Argentina) at the end of the Late Holocene (1800-700 years BP).

METHODS

Cranial nonmetric traits were registered in 159 adult individuals of both sexes. Biodistances were estimated between environmental and archaeological units through Smith's Mean Measure of Divergence, Mahalanobis distance, multidimensional scaling, F_st index, and R-matrix. Biological distances concerning temporal and geographic distances between the archaeological sites were also evaluated through Mantel and partial Mantel tests.

RESULTS

The environmental units showed statistically significant biodistances according to a distance spatial pattern but low genetic differentiation between groups. Additionally, some archaeological sites with median radiocarbon dates around 700 years BP also presented statistically significant biodistances.

CONCLUSIONS

Gene flow and genetic drift could be responsible for the observed biodistance results. These processes gave rise to different subpopulations of hunter-gatherers with a shared material culture, whose territorial limits could have been reinforced around 700 years BP or earlier.

Multi-dimensional machine learning approaches for fruit shape phenotyping in strawberry

BACKGROUND

Shape is a critical element of the visual appeal of strawberry fruit and is influenced by both genetic and non-genetic determinants. Current fruit phenotyping approaches for external characteristics in strawberry often rely on the human eye to make categorical assessments. However, fruit shape is an inherently multi-dimensional, continuously variable trait and not adequately described by a single categorical or quantitative feature. Morphometric approaches enable the study of complex, multi-dimensional forms but are often abstract and difficult to interpret. In this study, we developed a mathematical approach for transforming fruit shape classifications from digital images onto an ordinal scale called the Principal Progression of k Clusters (PPKC). We use these human-recognizable shape categories to select quantitative features extracted from multiple morphometric analyses that are best fit for genetic dissection and analysis.

RESULTS

We transformed images of strawberry fruit into human-recognizable categories using unsupervised machine learning, discovered 4 principal shape categories, and inferred progression using PPKC. We extracted 68 quantitative features from digital images of strawberries using a suite of morphometric analyses and multivariate statistical approaches. These analyses defined informative feature sets that effectively captured quantitative differences between shape classes. Classification accuracy ranged from 68% to 99% for the newly created phenotypic variables for describing a shape.

CONCLUSIONS

Our results demonstrated that strawberry fruit shapes could be robustly quantified, accurately classified, and empirically ordered using image analyses, machine learning, and PPKC. We generated a dictionary of quantitative traits for studying and predicting shape classes and identifying genetic factors underlying phenotypic variability for fruit shape in strawberry. The methods and approaches that we applied in strawberry should apply to other fruits, vegetables, and specialty crops.

A Parametric Bootstrap for the Mean Measure of Divergence

For more than 50 years the Mean Measure of Divergence (MMD) has been one of the most prominent tools used in anthropology for the study of non-metric traits. However, one of the problems, in anthropology including palaeoanthropology (more often there), is the lack of big enough samples or the existence of samples without sufficiently measured traits. Since 1969, with the advent of bootstrapping techniques, this issue has been tackled successfully in many different ways. Here, we present a parametric bootstrap technique based on the fact that the transformed θ, obtained from the Anscombe transformation to stabilize the variance, nearly follows a normal distribution with standard deviation $\sigma = 1 / \sqrt{N + 1/2}$σ=1/N+1/2, where N is the size of the measured trait. When the probabilistic distribution is known, parametric procedures offer more powerful results than non-parametric ones. We profit from knowing the probabilistic distribution of θ to develop a parametric bootstrapping method. We explain it carefully with mathematical support. We give examples, both with artificial data and with real ones. Our results show that this parametric bootstrap procedure is a powerful tool to study samples with scarcity of data.

The prevalence of under-nutrition among the tribal children in India: a systematic review

Tribal population of India constitutes 8.6% of the total population. They are more susceptible to undernutrition which is recognized as a prevalent health problem mainly because of their uncertainty of food supply, which has serious long term consequences for the child and adversely influences the development of the nation. The assessment of nutritional status is paramount importance because it impels to identify malnutrition which is an aggravation of morbidity and mortality. This review was aimed to find out the overall nutritional status of tribal preschool children. It also seek the major socio-cultural causes which influence the nutritional status from bio-cultural perspectives among the tribal children in India in order to make foundation of new research. In the present review, 41 papers on nutritional status of tribal preschool children published from January 1st 2000 till date, have been identified and included for analysis. To analyze the data, meta-analysis was done using MedCalc software.

The meta-analysis revealed the average rate of prevalence of underweight, stunting and wasting among the preschool tribal children of India was 42.96%; 44.82% and 23.69%. Among the 41 different studies these rates vary among the different tribal groups of different states associated with their socio-economic status (10%), their cultures of food consumption (10%), maternal education (15%), child feeding practices (20%), dietary deficit during pregnancy (25%) and poor nutrition of the child (52%).

Prevalence rate of under-nutrition and stunting is relatively high in comparison to WHO in India whereas, rate of wasting was lower in comparison to national level which reflect that immediate nutritional status was poor but chronic deficiency of nutrition was less. Nutritional education and short term appropriately planned nutritional intervention programs may also be useful for enhancing their nutritional status.

Non-metric traits, physiological stress indicators and paleopathological lesion on human skeletal remains from an early modern cemetery in Wyszyński Street, Wrocław, Poland (15th–18th centuries AD)

The anthropological analysis and assessment of the living conditions of historical populations should be comprehensive. Due to the scarcity of the well-preserved skeletal remains of a population buried at the cemetery in Wyszyński Street, Wrocław, each piece of information is particularly important in providing the complementary information about living and health conditions of late medieval and early modern inhabitants of Wrocław. This study aims to assess the living conditions of the late medieval and early modern Wrocław inhabitants. This aim was reached by analyzing the frequency of non-metrical cranial and postcranial traits, physiological stress indicators, and pathological lesions. The importance of these traits for population studies is also discussed. For 98 skeletons (22.4% male and 37.8% female) from a cemetery located in Wyszyński Street, (15th –18th centuries AD), the frequencies (p) of the following traits were analyzed: 38 non-metric cranial and 9 post-cranial traits, the morphological indicators of physiological stress (i.e. cribra orbitalia and enamel hypoplasia), caries and pathological lesions on bones. Among the non-metric cranial traits, only nine (23.7%) were not observed. The most frequent (from p= 0.30 to p=0.65) were foramen parietale, incisura frontalis accesoria seu foramen frontale accesorius, foramen supraorbitale, M3 mandibulare, ossicula suturae lambdoideae, and foramen mastoideum extrasuturale. Only three non-metric post-cranial traits were observed, but with low frequency (p=0.01 and 0.02). Cribra orbitalia was present in 37.0% of the skeletons while enamel hypoplasia was present in only 2.6%. Susceptibility to caries occurred in 68% of the cases while carries lesions occurred in 31.2%. Among pathological changes, the most frequent were degenerative changes (21.6%) and injuries (19.6%). Owing to the poorly preserved remains, the complete assessment of the living conditions of the population was difficult. The number of the observed traits would have probably increased, had the number of complete and well-preserved skeletons been higher. The results of the study suggest that the living conditions of the studied population were poor. It is hypothesized that the cemetery was a burial place for prisoners and victims of various epidemics that affected the city. But since the historical sources of this cemetery are scarce, it is hard to unequivocally state its purpose with certainty.

Additive genetic variation in the craniofacial skeleton of baboons (genus Papio) and its relationship to body and cranial size

OBJECTIVES

Determining the genetic architecture of quantitative traits and genetic correlations among them is important for understanding morphological evolution patterns. We address two questions regarding papionin evolution: (1) what effect do body and cranial size, age, and sex have on phenotypic (V_P ) and additive genetic (V_A ) variation in baboon crania, and (2) how might additive genetic correlations between craniofacial traits and body mass affect morphological evolution?

MATERIALS AND METHODS

We use a large captive pedigreed baboon sample to estimate quantitative genetic parameters for craniofacial dimensions (EIDs). Our models include nested combinations of the covariates listed above. We also simulate the correlated response of a given EID due to selection on body mass alone.

RESULTS

Covariates account for 1.2-91% of craniofacial V_P . EID V_A decreases across models as more covariates are included. The median genetic correlation estimate between each EID and body mass is 0.33. Analysis of the multivariate response to selection reveals that observed patterns of craniofacial variation in extant baboons cannot be attributed solely to correlated response to selection on body mass, particularly in males.

DISCUSSION

Because a relatively large proportion of EID V_A is shared with body mass variation, different methods of correcting for allometry by statistically controlling for size can alter residual V_P patterns. This may conflate direct selection effects on craniofacial variation with those resulting from a correlated response to body mass selection. This shared genetic variation may partially explain how selection for increased body mass in two different papionin lineages produced remarkably similar craniofacial phenotypes.

A COMPARISON OF GENETIC AND PHENOTYPIC CORRELATIONS

Genetic variances and correlations lie at the center of quantitative evolutionary theory. They are often difficult to estimate, however, due to the large samples of related individuals that are required. I investigated the relationship of genetic- and phenotypic-correlation magnitudes and patterns in 41 pairs of matrices drawn from the literature in order to determine their degree of similarity and whether phenotypic parameters could be used in place of their genetic counterparts in situations where genetic variances and correlations cannot be precisely estimated. The analysis indicates that squared genetic correlations were on average much higher than squared phenotypic correlations and that genetic and phenotypic correlations had only broadly similar patterns. These results could be due either to biological causes or to imprecision of genetic-correlation estimates due to sampling error. When only those studies based on the largest sample sizes (effective sample size of 40 or more) were included, squared genetic-correlation estimates were only slightly greater than their phenotypic counterparts and the patterns of correlation were strikingly similar. Thus, much of the dissimilarity between phenotypic- and genetic-correlation estimates seems to be due to imprecise estimates of genetic correlations. Phenotypic correlations are likely to be fair estimates of their genetic counterparts in many situations. These further results also indicate that genetic and environmental causes of phenotypic variation tend to act on growth and development in a similar manner.

Discovering phenotypic causal structure from nonexperimental data

The evolutionary potential of organisms depends on how their parts are structured into a cohesive whole. A major obstacle for empirical studies of phenotypic organization is that observed associations among characters usually confound different causal pathways such as pleiotropic modules, interphenotypic causal relationships and environmental effects. The present article proposes causal search algorithms as a new tool to distinguish these different modes of phenotypic integration. Without assuming an a priori structure, the algorithms seek a class of causal hypotheses consistent with independence relationships holding in observational data. The technique can be applied to discover causal relationships among a set of measured traits and to distinguish genuine selection from spurious correlations. The former application is illustrated with a biological data set of rat morphological measurements previously analysed by Cheverud et al. (Evolution 1983, 37, 895).

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.

Free-ranging Cayo Santiago rhesus monkeys (Macaca mulatta): III. Dental eruption patterns and dental pathology

This article describes the dental eruption patterns, dentition, and dental wear, including tooth loss and breakage, of the free-ranging population of rhesus monkeys (Macaca mulatta) on Cayo Santiago (CS), Puerto Rico, ranging from 24 hr to 25 years old. Of the 694 monkeys on the island in the year 1985, 688 (99.1%; 366 males, 322 females) were captured and the dentition of 685 subjects (98.7% of the total population; 366 males, 319 females) was examined. Animals ranged in age from less than 24 hr to 331 months (27.58 years), encompassing the entire life span of the CS macaques. Results demonstrated that the first deciduous teeth appeared as early as the third day of life and that the sequence of dental eruption was comparable to the pattern observed in laboratory rhesus. However, there were slight differences in the age of eruption of individual teeth. For example, the canines and third molars erupted about a year later in the CS macaques compared to some laboratory rhesus. Overall, CS rhesus had good oral health and dental condition although tooth wear, loss, and breakage were common in aged animals, especially in males. This report, combined with earlier studies on morphological characteristics and skeletal remains of the CS macaques, provides the basis for further studies on the biology, genetics, life history, and effects of the environment on rhesus monkeys.

Sex and age related non-metric variation of the human sternum in a Northwest Indian postmortem sample: a pilot study

Nonmetric traits are quasi-continuous skeletal traits that cannot be measured, but are observed as ranks or scores. In the present study, 8 nonmetric traits of sternum were examined from 343 sternums collected from the autopsy cases. More males, particularly older subjects, had a mesosternal foramen. The frequency of mesosternal foramen, arch-shaped prominence on the manubrium and radial strips on mesosternum significantly increased but that of lateral projection of manubrium decreased in the 30+ year age-group (older individuals). The denticulate costal incisures were found mostly in females and younger subjects; however, the crescent incisures and arch shaped prominence of manubrium were noticed mainly in males and older subjects. Majority of subjects had 'deep and narrow' jugular notch and xiphi-sternum was 'bifid' or 'oblong' in males and 'small sized' in females and younger subjects. Xiphi-sternum was generally found absent in females and younger subjects. The different shapes of costal incisures, jugular notch and xiphoid process were found to have significant sex and age differences. The number of cases with fused elements significantly increased with advancing age in both sexes. The sternal elements' fusion showed a variable pattern and hence was not found to be a reliable criterion. Chi-square analysis revealed significant differences between two age-groups of individual sex in the distribution of 3 of the 8 nonmetric traits. Significant sex and age differences were noticed in the frequency of occurrence of mesosternal foramina between the younger and older age groups in both the sexes. The logistic regression analysis of the scores classified 73.8% sternums to their sex category and 70.0% sternums to their age-group. Though present findings may not be enough for absolute personal identification, these traits can still serve as important identification tools if antemortem radiographs/CT or MRI scans are available for comparisons, or help classify an unknown sternum to its sex or age category based on the morphological features discussed in present study.

Skeletal growth and the changing genetic landscape during childhood and adulthood

Growth, development, and decline of the human skeleton are of central importance to physical anthropology. All processes of skeletal growth (longitudinal growth as well as gains and losses of bone mass) are subjected to environmental and genetic influences. These influences, and their relative contributions to the phenotype, can be asserted at any stage of life. We present here the gross phenotypic and genetic landscapes of four skeletal traits, and show how they vary across the life span. Phenotypic sex differences are found in bone diameter and cortical index (a ratio of cortical thickness over bone diameter) at a very early age and continue throughout most of life. Sexual dimorphism in summed cortical thickness and bone length, however, is not evident until shortly after the pubertal growth spurt. Genetic contributions (heritability) to these skeletal phenotypes are generally moderate to high. Bone length and bone diameter (which both scale with body size) tend to have the highest heritability, with heritability of bone length fairly stable across ages (with a notable dip in early childhood) and that of bone diameter peaking in early childhood. The bone traits summed cortical thickness and cortical index that may better reflect bone mass, a more plastic phenomenon, have slightly lower genetic influences, on average. Results from our phenotypic and genetic landscapes serve three key purposes: 1) demonstration of the integrated nature of the genetic and environmental underpinnings of skeletal form, 2) identification of periods of bone's relative sensitivity to genetic and environmental influences, 3) and stimulation of hypotheses predicting the effects of exposure to environmental variables on the skeleton, given variation in the underlying genetic architecture.

Zygomaticomaxillary suture shape analyzed with digital morphometrics: reassessing patterns of variation in American Indian and European populations

Typological classification of human zygomaticomaxillary suture (ZMS) shape is often used in forensic assessment of ancestry, following earlier studies reporting higher frequencies of "angled" sutures among American Indians and higher frequencies of "curved" sutures among Caucasians. In this paper we present a new method of digital morphometrics to quantify and compare ZMS shape in 60 American Indian and 60 European crania. Suture outlines were recorded as three-dimensional (3D) contours on digital models of adult male and female crania created with a portable 3D laser scanner. Each contour was represented by about four hundred point coordinates, which were transformed via Fourier analysis into amplitude coefficients suitable for use in linear discriminant analysis. Discriminant functions were created that accurately predicted group membership for 83% of the crania in the sample, after leave-one-out cross-validation. The results were compared with traditional typological classifications based on visual evaluation of ZMS shape, and the contour-based method was found to be more effective than the typological approach. However, the distribution of ZMS types within the two sample groups did not conform to previously reported patterns. This discrepancy indicates that ZMS shape may reflect not only genetic factors, but also environmental factors such as diet and stress. In addition, some evidence for sexual dimorphism in the zygomaticomaxillary complex was observed. Based on these findings, we recommend caution when using ZMS shape analysis in forensic ancestry determination.

Relationship of cranial robusticity to cranial form, geography and climate in Homo sapiens

Variation in cranial robusticity among modern human populations is widely acknowledged but not well-understood. While the use of "robust" cranial traits in hominin systematics and phylogeny suggests that these characters are strongly heritable, this hypothesis has not been tested. Alternatively, cranial robusticity may be a response to differences in diet/mastication or it may be an adaptation to cold, harsh environments. This study quantifies the distribution of cranial robusticity in 14 geographically widespread human populations, and correlates this variation with climatic variables, neutral genetic distances, cranial size, and cranial shape. With the exception of the occipital torus region, all traits were positively correlated with each other, suggesting that they should not be treated as individual characters. While males are more robust than females within each of the populations, among the independent variables (cranial shape, size, climate, and neutral genetic distances), only shape is significantly correlated with inter-population differences in robusticity. Two-block partial least-squares analysis was used to explore the relationship between cranial shape (captured by three-dimensional landmark data) and robusticity across individuals. Weak support was found for the hypothesis that robusticity was related to mastication as the shape associated with greater robusticity was similar to that described for groups that ate harder-to-process diets. Specifically, crania with more prognathic faces, expanded glabellar and occipital regions, and (slightly) longer skulls were more robust than those with rounder vaults and more orthognathic faces. However, groups with more mechanically demanding diets (hunter-gatherers) were not always more robust than groups practicing some form of agriculture.

An examination of Nubian and Egyptian biological distances: support for biological diffusion or in situ development?

Many authors have speculated on Nubian biological evolution. Because of the contact Nubians had with other peoples, migration and/or invasion (biological diffusion) were originally thought to be the biological mechanism for skeletal changes in Nubians. Later, a new hypothesis was put forth, the in situ hypothesis. The new hypothesis postulated that Nubians evolved in situ, without much genetic influence from foreign populations. This study examined 12 Egyptian and Nubian groups in an effort to explore the relationship between the two populations and to test the in situ hypothesis. Data from nine cranial nonmetric traits were assessed for an estimate of biological distance, using Mahalanobis D(2) with a tetrachoric matrix. The distance scores were then input into principal coordinates analysis (PCO) to depict the relationships between the two populations. PCO detected 60% of the variation in the first two principal coordinates. A plot of the distance scores revealed only one cluster; the Nubian and Egyptian groups clustered together. The grouping of the Nubians and Egyptians indicates there may have been some sort of gene flow between these groups of Nubians and Egyptians. However, common adaptation to similar environments may also be responsible for this pattern. Although the predominant results in this study appear to support the biological diffusion hypothesis, the in situ hypothesis was not completely negated.

A genomewide linkage scan for quantitative trait loci influencing the craniofacial complex in baboons (Papio hamadryas spp.)

Numerous studies have detected significant contributions of genes to variation in development, size, and shape of craniofacial traits in a number of vertebrate taxa. This study examines 43 quantitative traits derived from lateral cephalographs of 830 baboons (Papio hamadryas) from the pedigreed population housed at the Southwest National Primate Research Center. Quantitative genetic analyses were conducted using the SOLAR analytic platform, a maximum-likelihood variance components method that incorporates all familial information for parameter estimation. Heritability estimates were significant and of moderate to high magnitude for all craniofacial traits. Additionally, 14 significant quantitative trait loci (QTL) were identified for 12 traits from the three developmental components (basicranium, splanchnocranium, and neurocranium) of the craniofacial complex. These QTL were found on baboon chromosomes (and human orthologs) PHA1 (HSA1), PHA 2 (HSA3), PHA4 (HSA6), PHA11 (HSA12), PHA13 (HSA2), PHA16 (HSA17), and PHA17 (HSA13) (PHA, P. hamadryas; HSA, Homo sapiens). This study of the genetic architecture of the craniofacial complex in baboons provides the groundwork needed to establish the baboon as an animal model for the study of genetic and nongenetic influences on craniofacial variation.

Why are rare traits unilaterally expressed?: Trait frequency and unilateral expression for cranial nonmetric traits in humans

Based on an analysis of nonmetric trait databases from several large skeletal series in Northern Europe and South America, representing 27 bilateral traits, we report a predictable relationship between the frequency of nonmetric traits and the probability that they are expressed bilaterally. In a wider sampling of traits and populations, this study thus confirms the findings of an earlier study by Ossenberg ([1981] Am. J. Phys. Anthropol. 54:471-479), which reported the same relationship for two mandibular traits. This trend was previously explained by extending the multifactorial threshold model for discontinuous traits to incorporate either separate thresholds for unilateral or bilateral expression, or by a fuzzy threshold in which the probability of bilateral expression increases away from the median threshold value. We show that the trend is produced under the standard multifactorial threshold model for discontinuous traits simply if the within-individual or developmental instability variance remains relatively constant across the range of liability. Under this assumption, the number of individuals in which one side but not the other is pushed over the threshold for trait formation will be a larger proportion of the number of individuals expressing the trait when the trait frequency is low. As trait frequency increases, the significance of within-individual variance as a determinant of trait formation decreases relative to the genetic and among-individual environmental variance. These results have implications for interpreting nonmetric trait data as well as for understanding the prevalence of unilateral vs. bilateral expression of a wide variety of discontinuous traits, including dysmorphologies in humans.

Nonmetric cranial trait variation and prehistoric biocultural change in the Azapa Valley, Chile

Historically, interpretations of both biological and cultural change within the prehistoric Azapa Valley, northern Chile, have cited large-scale population movements, with replacement from complex societies located in the adjacent highlands to the east. Biological estimates of this change have traditionally relied upon biodistance estimates, using craniofacial measures of both deformed and nondeformed skulls. In order to evaluate whether large-scale prehistoric migrations occurred in the Azapa Valley, we examine biodistance results from nonmetric cranial traits for eight mortuary samples that represent all time periods of prehistoric occupation of the valley. None of the mean measures of divergences (MMD) among mortuary samples examined by this study were significant. These results suggest biological continuity in the Azapa Valley during 5,000 years of prehistory, with nonsignificant gene flow during the late Middle Horizon (AD 750-1100) and Late Intermediate (AD 1100-1476) periods. Biodistance results also suggest endogamy within coastal and inland populations. The broader implications of these results are discussed.

Concordance of cranial and dental morphological traits and evidence for endogamy in ancient Egypt

A biological affinities study based on frequencies of cranial nonmetric traits in skeletal samples from three cemeteries at predynastic Naqada, Egypt, confirms the results of a recent nonmetric dental morphological analysis. Both cranial and dental traits analyses indicate that the individuals buried in a cemetery characterized archaeologically as high status are significantly different from individuals buried in two other, apparently nonelite cemeteries and that the nonelite samples are not significantly different from each other. A comparison with neighbouring Nile Valley skeletal samples suggests that the high status cemetery represents an endogamous ruling or elite segment of the local population at Naqada, which is more closely related to populations in northern Nubia than to neighbouring populations in southern Egypt.

Size and shape of human cranial sutures--a new scoring method

A method for the differentiation of sutural patterns of the human cranial vault is introduced. Three criteria of differentiation are considered, one for size and two for shape: 1) maximal shape extension; 2) basic configuration; 3) secondary protrusion. The method is illustrated here for the coronal and lambdoid sutures of 70 recent Italian skulls (35 adult males and 35 adult females). Differences between coronal and lambdoid sutural size and shape can be detected analytically; for example, the coronal suture commonly shows lesser degrees of shape extension, a simpler basic configuration, and an absence of secondary protrusion. Heterogeneity within each suture, as well as a relationship among corresponding sections and between the three criteria adopted, have been also observed; symmetry predominates for both the sutures, and sexual differences are slight.

Genetic correlations between sides and heritability of asymmetry for nonmetric traits in rhesus macaques on Cayo Santiago

The use of nonmetric traits for estimation of biological distance is a long-standing practice in biological anthropology. Nonmetric traits can be scored using either the individual or the side of the individual as the unit of measure. If sides of the individual are genetically correlated the use of sides would produce redundant genetic information. For this reason, Korey (Am. J. Phys. Anthropol. 53:19-23, 1980) argues for the use of individuals as the unit of measure for nonmetric traits. Ossenberg (Am. J. Phys, Anthropol. 54:471-479, 1981), however, argues that bilateral occurrence of nonmetric traits indicates greater genetic liability for the trait and that therefore the sides are the more biologically correct unit of measure. Genetic correlations for 13 cranial nonmetric traits are estimated for a sample of rhesus macaque skeletons from Cayo Santiago. In addition, heritability of asymmetry is estimated for these 13 traits as a test of Ossenberg's contention that asymmetry is genetically influenced. Significant genetic correlations between sides support Korey's contention that nonmetric traits should be scored by individual. Only two asymmetry heritabilities were significantly different from zero, providing no significant support for Ossenberg's contention that asymmetry is genetically determined. Our results support the theory that asymmetry represents a measure of the ability of an organism to buffer stresses. Therefore, a measure of the heritability of asymmetry is a measure of the heritability of the ability to buffer stresses. This ability does not appear to be heritable in this sample.