Mechanical and morphometric approaches to body mass estimation in rhesus macaques: A test of skeletal variables.

OBJECTIVES: Estimation of body mass from skeletal metrics can reveal important insights into the paleobiology of archeological or fossil remains. The standard approach constructs predictive equations from postcrania, but studies have questioned the reliability of traditional measures. Here, we examine several skeletal features to assess their accuracy in predicting body mass. MATERIALS AND METHODS: Antemortem mass measurements were compared with common skeletal dimensions from the same animals postmortem, using 115 rhesus macaques (male: n = 43; female: n = 72). Individuals were divided into training (n = 58) and test samples (n = 57) to build and assess Ordinary Least Squares or multivariate regressions by residual sum of squares (RSS) and AIC weights. A leave-one-out approach was implemented to formulate the best fit multivariate models, which were compared against a univariate and a previously published catarrhine body-mass estimation model. RESULTS: Femur circumference represented the best univariate model. The best model overall was composed of four variables (femur, tibia and fibula circumference and humerus length). By RSS and AICw, models built from rhesus macaque data (RSS = 26.91, AIC = -20.66) better predicted body mass than did the catarrhine model (RSS = 65.47, AIC = 20.24). CONCLUSION: Body mass in rhesus macaques is best predicted by a 4-variable equation composed of humerus length and hind limb midshaft circumferences. Comparison of models built from the macaque versus the catarrhine data highlight the importance of taxonomic specificity in predicting body mass. This paper provides a valuable dataset of combined somatic and skeletal data in a primate, which can be used to build body mass equations for fragmentary fossil evidence.

Conceptual issues in hominin taxonomy: Homo heidelbergensis and an ethnobiological reframing of species.

Efforts to name and classify Middle Pleistocene Homo, often referred to as "Homo heidelbergensis" are hampered by confusing patterns of morphology but also by conflicting paleoanthropological ideologies that are embedded in approaches to hominin taxonomy, nomenclature, and the species concept. We deconstruct these issues to show how the field's search for a "real" species relies on strict adherence to pre-Darwinian essentialist naming rules in a post-typological world. We then examine Middle Pleistocene Homo through the framework of ethnobiology, which examines on how Indigenous societies perceive, classify, and name biological organisms. This research reminds us that across human societies, taxonomies function to (1) identify and classify organisms based on consensus pattern recognition and (2) construct a stable nomenclature for effective storage, retrieval and communication of information. Naming Middle Pleistocene Homo as a "real" species cannot be verified with the current data; and separating regional groups into distinct evolutionary lineages creates taxa that are not defined by readily perceptible or universally salient differences. Based on ethnobiological studies of this kind of patterning, referring to these hominins above the level of the species according to their generic category with modifiers (e.g., "European Middle Pleistocene Homo") is consistent with observed human capabilities for cognitive differentiation, is both necessary and sufficient given the current data, and will allow for the most clear communication across ideologies going forward.

Lack of biological mortality bias in the timing of dental formation in contemporary children: Implications for the study of past populations.

OBJECTIVES: Biological mortality bias is the idea that individuals who perish (non-survivors) are biologically distinct from those who survive (survivors). If biological mortality bias is large enough, bioarchaeological studies of nonsurvivors (skeletal samples) cannot accurately represent the experiences of the survivors of that population. This effect is particularly problematic for the study of juvenile individuals, as growth is particularly sensitive to environmental insults. In this study, we test whether biological mortality bias exists in one dimension of growth, namely dental development. MATERIALS AND METHODS: Postmortem computed tomography scans of 206 children aged 12 years and younger at death were collected from two institutions in the United States and Australia. The sample was separated into children dying from natural causes as proxies for non-survivors and from accidental causes as proxies for survivors. Differences in the timing of dental development were assessed using sequential logistic regressions between dental formation stages and residual analysis of dental minus chronological age. RESULTS: No consistent delay in age of attainment of dental stages was documented between survivors and non-survivors. Delays between survivors and non-survivors in dental relative to chronological age were greatest for infants, and were greater for females than for males. DISCUSSION: Lack of biological mortality bias in dental development reinforces confidence in juvenile age estimates and therefore in skeletal growth profiles and growth studies. As dental development is known to be less environmentally sensitive than skeletal growth and development, further studies should examine biological mortality bias in long bone length.

A comparison of subadult skeletal and dental development based on living and deceased samples.

OBJECTIVES: A fundamental assumption in biological anthropology is that living individuals will present with different growth than non-survivors of the same population. The aim is to address the question of whether growth and development data of non-survivors are reflective of the biological consequences of selective mortality and/or stress. MATERIALS AND METHODS: The study compares dental development and skeletal growth collected from radiographic images of contemporary samples of living and deceased individuals from the United States (birth to 20 years) and South Africa (birth to 12 years). Further evaluation of deceased individuals is used to explore differential patterns among manners of death (MOD). RESULTS: Results do not show any significant differences in skeletal growth or dental development between living and deceased individuals. However, in the South African deceased sample the youngest individuals exhibited substantially smaller diaphyseal lengths than the living sample, but by 2 years of age the differences were negligible. In the US sample, neither significant nor substantial differences were found in dental development or diaphyseal length according to MOD and age (>2 years of age), though some long bones in individuals <2 years of age did show significant differences. No significant differences were noted in diaphyseal length according to MOD and age in the SA sample. DISCUSSION: The current findings refute the idea that contemporary deceased and living individuals would present with differential growth and development patterns through all of ontogeny as well as the assumptions linking short stature, poor environments, and MOD.

Paleodemographic analysis of age at death for a population of Black Sea Scythians: An exploration by using Bayesian methods.

OBJECTIVES: Studies of the demography of past populations involving deterministic life tables can be criticized for ignoring the errors of estimation. Bayesian methods offer an alternative, by focusing on the uncertainty of the estimates, although their results are often sensitive to the choice of prior distributions. The aim of this study is to explore a range of Bayesian methods for estimating age at death for a population of nomadic warriors-Scythians from the Black Sea region. MATERIALS AND METHODS: In total, skeletons of 312 individuals (93 children and 219 adults) from Glinoe (Moldova), dated to the 5th-2nd century BCE, were examined. We unified the age categories corresponding to different aging methods, allowing an application of a probabilistic assessment of the age categorization. A hierarchical Bayesian multinomial-Dirichlet-Dirichlet model was applied, with a hypothetical, subjective reference population, a real reference population, and no reference. RESULTS: Stationary-population life expectancy was estimated as 27.7 years (95% CI: 25.1-30.3) for a newborn (e0 ), and 16.4 years (14.0-19.0) for 20-year-olds (e20 ), although with high uncertainty, and sensitive to the model specification. Slight differences in longevity between different social strata and between the Classical and Late chronological periods were found, although with high estimation errors. A more robust finding, confirming earlier studies, was a high probability of death in young adulthood, which could depend on Scythian lifestyle (conflicts, wars). DISCUSSION: Our study shows a way to overcome some limitations of broad age categorization by using the Bayesian approach with alternative model specifications, allowing to assess the impact of reference populations.

Ingestive behaviors in bearded capuchins (Sapajus libidinosus).

The biomechanical and adaptive significance of variation in craniodental and mandibular morphology in fossil hominins is not always clear, at least in part because of a poor understanding of how different feeding behaviors impact feeding system design (form-function relationships). While laboratory studies suggest that ingestive behaviors produce variable loading, stress, and strain regimes in the cranium and mandible, understanding the relative importance of these behaviors for feeding system design requires data on their use in wild populations. Here we assess the frequencies and durations of manual, ingestive, and masticatory behaviors from more than 1400 observations of feeding behaviors video-recorded in a wild population of bearded capuchins (Sapajus libidinosus) at Fazenda Boa Vista in Piauí, Brazil. Our results suggest that ingestive behaviors in wild Sapajus libidinosus were used for a range of food material properties and typically performed using the anterior dentition. Coupled with previous laboratory work indicating that ingestive behaviors are associated with higher mandibular strain magnitudes than mastication, these results suggest that ingestive behaviors may play an important role in craniodental and mandibular design in capuchins and may be reflected in robust adaptations in fossil hominins.

Covariation of fetal skull and maternal pelvis during the perinatal period in rhesus macaques and evolution of childbirth in primates.

A large brain combined with an upright posture in humans has resulted in a high cephalopelvic proportion and frequently obstructed labor. Fischer and Mitteroecker [B. Fischer, P. Mitteroecker, Proc. Natl. Acad. Sci. U.S.A. 112, 5655-5660 (2015)] proposed that the morphological covariations between the skull and pelvis could have evolved to ameliorate obstructed labor in humans. The availability of quantitative data of such covariation, especially of the fetal skull and maternal pelvis, however, is still scarce. Here, we present direct evidence of morphological covariations between the skull and pelvis using actual mother-fetus dyads during the perinatal period of Macaca mulatta, a species that exhibits cephalopelvic proportions comparable to modern humans. We analyzed the covariation of the three-dimensional morphology of the fetal skull and maternal pelvis using computed tomography-based models. The covariation was mostly observed at the pelvic locations related to the birth canal, and the forms of the birth canal and fetal skull covary in such a way that reduces obstetric difficulties. Therefore, cephalopelvic covariation could have evolved not only in humans, but also in other primate taxa in parallel, or it could have evolved already in the early catarrhines.

Historical human remains identification through maternal and paternal genetic signatures in a founder population with extensive genealogical record.

OBJECTIVES: We describe a method to identify human remains excavated from unmarked graves in historical Québec cemeteries by combining parental-lineage genetic markers with the whole-population genealogy of Québec contained in the BALSAC database. MATERIALS AND METHODS: The remains of six men were exhumed from four historical cemeteries in the province of Québec, Canada. DNA was extracted from the remains and genotyped to reveal their mitochondrial and Y-chromosome haplotypes, which were compared to a collection of haplotypes of genealogically-anchored modern volunteers. Maternal and paternal genealogies were searched in the BALSAC genealogical record for parental couples matching the mitochondrial and the Y-chromosome haplotypic signatures, to identify candidate sons from whom the remains could have originated. RESULTS: Analysis of the matching genealogies identified the parents of one man inhumed in the cemetery of the investigated parish during its operating time. The candidate individual died in 1833 at the age of 58, a plausible age at death in light of osteological analysis of the remains. DISCUSSION: This study demonstrates the promising potential of coupling genetic information from living individuals to genealogical data in BALSAC to identify historical human remains. If genetic coverage is increased, the genealogical information in BALSAC could enable the identification of 87% of the men (n = 178,435) married in Québec before 1850, with high discriminatory power in most cases since >75% of the parental couples have unique biparental signatures in most regions. Genotyping and identifying Québec's historical human remains are a key to reconstructing the genomes of the founders of Québec and reinhuming archeological remains with a marked grave.

Particle size matters: The effect of particle size on carbon and oxygen isotope composition of bone hydroxyapatite.

OBJECTIVES: Stable isotope studies often focus on hydroxyapatite (bioapatite) to answer questions of paleodiet, paleomobility, and paeloenvironment. This study seeks to determine the effect that sample particle size (in particular SA:V, or surface area to volume ratios) has on measured carbon and oxygen stable isotope values (δ13 C and δ18 O) in bone hydroxyapatite. MATERIALS AND METHODS: Previously ground Homo sapiens sapiens cortical bone samples were subdivided using geological screens to obtain three separate sub-samples, differing only in their particle size. These aliquots (n = 60) were then treated using established protocols to remove any exogenous organic material (2.5% NaOH) and adsorbed carbonates (0.1 M CH3 COOH), and analyzed for δ13 C and δ18 O using a Kiel-IV Carbonate Device coupled to a Thermo-Finnigan DeltaPlus IRMS. RESULTS: Data obtained indicate that decreased particle size leads to increases in both δ13 C and δ18 O, with oxygen isotope values being more dramatically affected. Specifically, it is possible to produce isotopic shifts of as much as 1.0‰ and 4.0‰ for δ13 C and δ18 O, respectively, solely by analyzing different sized particles from the same individual, bone, and sample. DISCUSSION: Based upon the variability seen between different size fractions from the same sample, it is clear that particle size has a meaningful impact on carbon and oxygen isotope composition. We attribute these shifts to the differential adsorption or precipitation of environmental carbon and oxygen during pretreatment. We recommend that particle size be added to the list of potential variables affecting isotope composition, alongside other factors including diagenesis, reagent concentration, and treatment time. We would also note that while most individuals exhibit consistent changes, some do not, and thus further investigation into these phenomena is warranted.

Improving anatomical stature estimation method. The relationship between living stature and intervertebral disc thickness.

Anatomical stature estimation methods reconstruct stature for skeletal specimens by adding up the heights of skeletal elements contributing to stature. In addition, these estimations factor in a certain amount of soft tissue known as "soft tissue correction". Our study focuses on the relationship between living stature and one of the major soft tissue contributors to stature: the intervertebral disc thickness/height. The purpose of this study was to clarify whether intervertebral disc thickness is greater in tall individuals and whether there is a linear correlation between stature and intervertebral disc height. To conduct this study, we utilized a subsample of the Northern Finland Birth Cohort of 1966 (n = 12,058) with known stature. We measured vertebral heights and intervertebral disc heights from low back MRI examination performed at the age of 46 years (n = 200). All subjects were considered healthy with no spinal injuries or pathologies. Our results clearly indicate that stature and intervertebral disc height have positive, statistically significant association. According to our results it is advisable to take into account the individual's skeletal height when soft tissue corrections for anatomical stature estimations are performed. Further studies utilizing full body MRI are needed to produce more accurate soft tissue corrections.

Determining primates bite force from histological tooth sections.

OBJECTIVES: The ability to accurately estimate bite force (BF) in extant and fossil primates is valuable to biological anthropologists. BF is generally evaluated using complex jaw musculature and lever arm analyses employing numerous assumptions and requiring complete cranial morphology. Here, a simple method to determine BF from data measured on histological sections of fossil teeth is proposed. METHODS: Published sections of molar teeth encompassing 27 different extinct and extant primates dating back to as early as 17 million years ago were examined. Focusing on the cusp region, the extracted data include characteristic enamel thickness dc and dentin horn angle φ. The occlusal force needed to fracture a cusp, PF , was determined from these variables with the aid of a finite element stress analysis similarly to a previous study on postcanine human teeth. The bite force was obtained by linking BF to PF using a universal constant. RESULTS: The measured variables dc and φ are conclusively linked. This link produces a virtually constant fracture force PF and in turn bite force BF for all cusps in the molar row. An explicit formula tying BF to dc and φ was derived. For nonhominin taxa the bite force, molar crown area, and body mass are found to be intimately related. The case of hominins is more involved. The so determined BF is gender-averaged, with the bite force of males estimated to be ≈12% greater than that of females. CONCLUSIONS: The use of "fracture mechanics" concepts from mechanics of materials facilitates determination of critical bite force in primates based on characteristic enamel thickness dc and dentin horn angle φ as extracted from histological sections of molar teeth. This novel approach enables quantitative insight into the role played by crown area, body mass and bite force on evolutionary trends. The conclusive link between cuspal enamel thickness and dentin horn angle facilitates optimal food processing without hindering cusp resilience. The proposed approach may be extended to mammals having asymmetric cusp structures.

Height inequality and socioeconomic implications in Korea: analysis of individuals born between 1890 and 1919.

This study examined height inequality as an indicator of income inequality during the colonial period (1910-1945) in Korea. Data were taken from a sample of 1796 male prisoners from a wide range of geographical locations and with varied socioeconomic backgrounds. Height inequality was measured using the coefficient of variation of height (CV) for each birth decade. The results indicated that height inequality, as measured by the CV, increased slightly from 3.32 to 3.35 for the birth decades 1890-99 and 1900-09, then jumped to 3.50 for the birth decade 1910-19. Considering the Kuznets curve, the presented results have socioeconomic implications for Japan's impact in Korea, at least during the early colonial period.

Comparison of extraction methods for recovering ancient microbial DNA from paleofeces.

OBJECTIVES: Paleofeces are valuable to archeologists and evolutionary biologists for their potential to yield health, dietary, and host information. As a rich source of preserved biomolecules from host-associated microorganisms, they can also provide insights into the recent evolution and changing ecology of the gut microbiome. However, there is currently no standard method for DNA extraction from paleofeces, which combine the dual challenges of complex biological composition and degraded DNA. Due to the scarcity and relatively poor preservation of paleofeces when compared with other archeological remains, it is important to use efficient methods that maximize ancient DNA (aDNA) recovery while also minimizing downstream taxonomic biases. METHODS: In this study, we use shotgun metagenomics to systematically compare the performance of five DNA extraction methods on a set of well-preserved human and dog paleofeces from Mexico (~1,300 BP). RESULTS: Our results show that all tested DNA extraction methods yield a consistent microbial taxonomic profile, but that methods optimized for ancient samples recover significantly more DNA. CONCLUSIONS: These results show promise for future studies that seek to explore the evolution of the human gut microbiome by comparing aDNA data with those generated in modern studies.

Reliability of Physiological Responses Induced by Basic Emotions: A Pilot Study.

BACKGROUND: Although emotion-specific autonomic responses based on the discrete theory of emotion have been widely studied, studies on the reliability of physiological responses to emotional stimuli are limited. In this study, we aimed to assess the reliability of physiological changes induced by the six basic emotions (happiness, sadness, anger, fear, disgust, and surprise) that were measured during 10 weekly repeated experiments. METHODS: Twelve college students participated, and in each experiment, physiological signals were collected before and while participants were watching emotion-provoking film clips. Additionally, the participants self-evaluated the emotions that they experienced during the film presentation at the end of each emotional stimulus. To avoid adaptation of participants to identical stimuli during repeated measurements, we used 10 different film clips for each emotion, and thus a total of 60 film clips over 10 weeks were used. Physiological features, such as skin conductance level (SCL), fingertip temperature (FT), heart rate (HR), and blood volume pulse (BVP), were extracted from the physiological signals. Two reliability indices, Cronbach's alpha and intraclass correlation coefficient, were calculated from the physiological features to assess internal consistency and interrater reliability, respectively. RESULTS: We found that SCL, HR, and BVP measured during the emotion-provoking phase over the 10 weekly sessions were more reliable than those assessed at baseline. Furthermore, SCL, HR, and BVP from the emotion-provoking phase exhibited excellent internal consistency and interrater reliability. CONCLUSIONS: Our findings suggest that these features can be used as reliable physiological indices in emotion studies. The results also support the significance of physiological signals as meaningful indicators for emotion recognition in HCI (human computer interface) area.

Identifying Archaeological Bone via Non-Destructive ZooMS and the Materiality of Symbolic Expression: Examples from Iroquoian Bone Points.

Today, practical, functional and symbolic choices inform the selection of raw materials for worked objects. In cases where we can discern the origin of worked bone, tooth, ivory and antler objects in the past, we assume that similar choices are being made. However, morphological species identification of worked objects is often impossible due to the loss of identifying characteristics during manufacture. Here, we describe a novel non-destructive ZooMS (Zooarchaeology by Mass Spectrometry) method which was applied to bone points from Pre-Contact St. Lawrence Iroquoian village sites in southern Quebec, Canada. The traditional ZooMS technique requires destructive analysis of a sample, which can be problematic when dealing with artefacts. Here we instead extracted proteins from the plastic bags in which the points had been stored. ZooMS analysis revealed hitherto unexpected species, notably black bear (Ursus americanus) and human (Homo sapiens sapiens), used in point manufacture. These surprising results (confirmed through genomic sequencing) highlight the importance of advancing biomolecular research in artefact studies. Furthermore, they unexpectedly and exceptionally allow us to identify and explore the tangible, material traces of the symbolic relationship between bears and humans, central to past and present Iroquoian cosmology and mythology.

Living with observational data in biological anthropology.

The establishment of cause and effect relationships is a fundamental objective of scientific research. Many lines of evidence can be used to make cause-effect inferences. When statistical data are involved, alternative explanations for the statistical relationship need to be ruled out. These include chance (apparent patterns due to random factors), confounding effects (a relationship between two variables because they are each associated with an unmeasured third variable), and sampling bias (effects due to preexisting properties of compared groups). The gold standard for managing these issues is a controlled randomized experiment. In disciplines such as biological anthropology, where controlled experiments are not possible for many research questions, causal inferences are made from observational data. Methods that statisticians recommend for this difficult objective have not been widely adopted in the biological anthropology literature. Issues involved in using statistics to make valid causal inferences from observational data are discussed.

Signatures of stress: Pilot study of accentuated laminations in porcine enamel.

OBJECTIVES: To examine the correspondence between histology of defective enamel and documented developmental stressors using a pig proxy; with the ultimate goal of differentiating, within enamel, stressor types including birth, weaning, surgery, and social disruption. MATERIALS AND METHODS: Lower first molars were removed from five female and four male domestic pigs, thin-sectioned and examined with phase-contrast microscopy for daily and accentuated laminations. Known stressor events include birth; processing (nonanesthetized clipping of needle teeth, tail docking, and castration) within 5 days of birth; weaning in nursery; penned vaccination 2-4 days later; transfer to grower barn; additional vaccination. Timing of accentuated laminations was derived from counts of daily laminations from birth and between accentuated laminations; or, more commonly, from measurements between accentuated laminations divided by average daily lamination width. RESULTS: Acknowledging our small sample, we confirm daily periodicity of laminations in pig enamel. Lamination width varies among sexes (males wider), cusps, crown decile, and enamel depth (wider toward surface). Accentuated laminations occurred at reconstructed median ages of 3.7 days after birth (cf. up to 5 days for "processing"), 19.7 days (cf. 18-25 days for weaning), 4.5 days (cf. 2-4 days for vaccination) and ~39.5 days (cf. 63 days for transfer). DISCUSSION: Encouragingly, the timing of known stressors (birth, surgical processing, weaning, and vaccination) can be determined with high precision, in good thin sections, from accentuated laminations in pig enamel. Timing of transfer was poorly estimated, likely reflecting lesser severity and the occurrence of undocumented stress events in the nursery.

Diaphysator: An online application for the exhaustive cartography and user-friendly statistical analysis of long bone diaphyses.

The cross-sectional geometry (CSG) of long bone diaphyses is used in bioanthropology to evaluate their resistance to biomechanical constraints and to infer life-history-related patterns such as mobility, activity specialization or intensity, sexual dimorphism, body mass and proportions. First limited by technical analytical constraints to the analysis of one or two cross sections per bone, it has evolved into the analysis of cross sections of the full length of the diaphyseal part of long bones. More recently, researchers have developed analytical tools to map the cortical thickness of entire diaphyses to evaluate locomotor signatures. However, none of these analytical tools are easy to use for scientists who are not familiar with computer programming, and some statistical procedures-such as mapping the correlation coefficients of the diaphyseal thickness with various parameters have yet to be made available. Therefore, we developed an automated and open-source application that renders those analyses (both CSG and cortical thickness) in a semiautomated and user friendly manner. This application, called "Diaphysator", is associated with another free software ("Extractor", presented in Dupej et al. (2017). American Journal of Physical Anthropology, 164, 868-876). Diaphysator can be used as an online application (https://diaphysator.shinyapps.io/maps) or as a package for R statistical software. Along with the mean maps of cortical thickness and mean CSG parameter graphs, the users can evaluate the correlations and partial correlations of both CSG parameters at every cross section along the diaphyseal length, and cortical thickness data points of the entire diaphysis, with any factor such as age, sex, stature, and body mass.

A New Tool for Digital Alignment in Virtual Anthropology.

The study of the fossil record is fundamental to understand the evolution of traits. Because fossil remains are often fragmented and/or deformed by taphonomic processes, a preliminary realignment of their constituent parts is often necessary to properly interpret their shapes. In virtual anthropology, these procedures are carried out using digital models of the remains. We present a new semi-automatic alignment R software, Digital Tool for Alignment (DTA), which uses the shape information contained in a reference sample to find the best alignment solution for the disarticulated regions. We tested DTA on three different case-studies: (1) a sample of 14 primate species including both male and female individuals, (2) a simulated, disarticulated skull of Homo sapiens, and (3) a real disarticulated human fossil specimen, Amud 1 (Homo neanderthalensis). In the first case study, we simulated disarticulation directly on digital models of the primate skulls and tested alignment quality as a function of phylogenetic proximity, sex, and body size. In the second, we compared DTA to manual alignments conducted for the same digital models. Finally, we performed DTA on a real-world case study. We found that phylogenetic proximity provides is the most important factor for alignment efficiency. However, sex and allometric effects might also be important and should therefore be taken into account at selecting reference models for alignments. DTA performs at least as well as manual alignments. Yet, as compared to manual procedures, it is faster, requires no prior anatomical knowledge and expertise and allows indefinite manipulation of the fossil items. Anat Rec, 302:1104-1115, 2019. © 2019 Wiley Periodicals, Inc.

Quantification of anatomical variation at the atlanto-occipital articulation: morphometric resolution of commingled human remains within the repatriation documentation process.

Within many institutional collections are skeletal and mummified human remains representing a part of our species' adaptation and evolution to various biocultural environments. Archaeologically recovered individuals come from deep into our past, and possess information that provides insight into population history, genetics, diet, health and other questions relevant to all living peoples. Academic concerns have been raised regarding the reinterment of these collections due to the rise of the international repatriation movement, the passage of various laws and implementation of institutional policies. While all potential research questions cannot be anticipated, the proactive documentation of collections is one way to ensure primary data are maintained for future study. This paper explores developments in digitization technology that allow the archive of virtual copies of human remains, and an example of how anatomical and archaeological collections can be digitized towards pragmatic research goals. The anatomical variability of the human atlanto-occipital (AO) articular surfaces was studied using non-metric categorical shape, 2D measurement and 3D morphometric analyses to provide reference standards for the reassociation of individuals from commingled skeletal remains, such as found in some archaeological sites or forensic investigations including mass grave or mass disaster recovery scenes. Results suggest that qualitative shape observations and caliper-derived measurements of the articulating AO condyles tend to display significant sexual dimorphism and biological ancestry-related size and shape differences. Variables derived from a scanned 3D mesh, such as condylar angle and articular surface curvature, quantify biomechanical variation and display a stronger congruency within individuals. It is recommended that a two-stage approach involving initial screening and identification of possible reassociation candidates is accomplished with a linear osteometric approach, followed by 3D laser scanning of the candidate joint surfaces for morphometric analyses to confirm reassociations when destructive DNA typing is not allowed or otherwise impractical due to cost or other resource restrictions.