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.

Architecture of the femoral and tibial diaphyses in relation to body mass and composition: Research from whole-body CT scans of adult humans.

OBJECTIVES: Recent investigations have evaluated the influence of body composition on long bones in order to overcome the limits of body mass (BM) estimation methods and eventually lead to studying nutrition in past populations. Knowing how fat mass (FM) and fat-free mass (FFM) impact the skeleton would also enhance the understanding of mobility, activity, and locomotion derived from bone architecture. We investigated the relationship between BM and composition, and the architecture of the entire tibial and femoral diaphyses in an adult sample representative of a wide range of variation in age, BM, and composition. MATERIALS AND METHODS: Body composition was measured directly from 78 whole-body CT scans for which the age, sex, BM, and stature were recorded. The entire diaphyseal thickness, volume, curvature, and cross-sectional geometry parameters of both the femur and tibia were numerically extracted. RESULTS: FM correlates with large portions of the femoral thickness in females only. FFM correlates with the femoral diaphysis in males but not in females. FFM correlates with the tibia architecture in both sexes, while FM is correlated in males exclusively. DISCUSSION: BM and body components influence the architecture of the diaphysis of lower limb long bones in sex-specific patterns that are mostly reflected in their thickness and can be recorded, in some cases, for their strength, rigidity, and volume. Our results suggest that (1) long bone diaphyses should be thoroughly studied, as a whole, when possible; and (2) BM and body components should be accounted for when deriving activity, mobility, or locomotion patterns from cortical bone.

Sex-specific functional adaptation of the femoral diaphysis to body composition.

OBJECTIVES: The human femoral diaphysis is often used to reconstruct loading histories (mobility, activity, body mass). The proximal femur is known to be differentially affected by changes in total fat-mass (FM), fat-free mass (FFM), and body fat percentage (BF%), but the adaptation of the entire diaphysis to body composition has not been thoroughly characterized to date. Understanding how the femoral diaphysis adapts to body components would benefit biomechanical interpretations of the femoral variation and nutrition-related studies. METHODS: Combining various methods from clinical nutrition, biological anthropology, and geometric morphometrics, we evaluated the correlation of measures taken on the entire femoral diaphysis with estimated FM, FFM, and BF% from 61 CT scans (17 females, 44 males). The sample was predominantly composed of people with obesity. RESULTS: Cortical area of the cross-sections and local cortical thickness showed high correlation with BF% in particular, in females only. The curvature significantly decreased with FM and BF% in both sexes. The lowest correlations are found with FFM. CONCLUSIONS: The observed sexual dimorphism is consistent with differing aging processes; cortical bone decreases in females through endosteal resorption while it remains almost constant in males who compensate for endosteal resorption by periosteal apposition on the diaphyseal surface. The functional adaptation to compressive forces indicates a systemic endosteal apposition of bone material with increased BF% and FM in females only. FM and BF% are linked to a straighter femur in both sexes, suggesting an optimization of the resistance to compressive loads by distributing them more linearly along the entire diaphysis.

Semiautomatic extraction of cortical thickness and diaphyseal curvature from CT scans.

The understanding of locomotor patterns, activity schemes, and biological variations has been enhanced by the study of the geometrical properties and cortical bone thickness of the long bones measured using CT scan cross-sections. With the development of scanning procedures, the internal architecture of the long bones can be explored along the entire diaphysis. Recently, several methods that map cortical thickness along the whole femoral diaphysis have been developed. Precise homology is vital for statistical examination of the data; however, the repeatability of these methods is unknown and some do not account for the curvature of the bones. We have designed a semiautomatic workflow that improves the morphometric analysis of cortical thickness, including robust data acquisition with minimal user interaction and considering the bone curvature. The proposed algorithm also performs automatic landmark refinement and rigid registration on the extracted morphometric maps of the cortical thickness. Because our algorithm automatically reslices the diaphysis into 100 cross-sections along the medial axis and uses an adaptive thresholding method, it is usable on CT scans that contain soft tissues as well as on bones that have not been oriented specifically prior to scanning. Our approach exhibits considerable robustness to error in user-supplied landmarks, suppresses distortion caused by the curvature of the bones, and calculates the curvature of the medial axis.

Detecting menarcheal status through dental mineralization stages?

Menarche is an indicator frequently used to study variation in growth, development, and related health conditions among members of living populations. As a life event, menarche is often associated with changes in an individual's social identity. The reproductive lifespan, which for females starts with menarche, is a paramount feature of palaeodemographic studies. Determination of menarche status from the skeletal remains of individuals of past populations can be obtained by assessing the developmental status of the iliac crest, as well as the hand and wrist bones, which are, unlike teeth, often poorly recovered in bioarchaeological contexts. The present study seeks to evaluate the link between dental mineralization and menarche in a population of known menarche status. The relationship between permanent teeth mineralization and menarche status was investigated by using data of developing permanent teeth (167 radiographs) rated in accordance with the well-known standards of Demirjian et al. and Moorrees et al. collected among 73 living French females of known menarcheal status. Using correlation ratios, GLMM and CART algorithm, menarcheal status is correlated with mineralization of the premolars. Menarcheal status is predicted correctly for 92 and 77% of radiographs of the learning and validation samples, respectively. Although promising, the results require caution prior to generalization to other populations. The age of menarche in this particular sample may simply coincide with the development of the premolars in this particular sample. Therefore, further investigation applied to populations with various mean ages of menarche is required in order to provide new evidence of variation in human growth and development from the correspondence between the mineralization of the permanent teeth and menarche.

Body composition estimation from selected slices: equations computed from a new semi-automatic thresholding method developed on whole-body CT scans.

BACKGROUND: Estimating volumes and masses of total body components is important for the study and treatment monitoring of nutrition and nutrition-related disorders, cancer, joint replacement, energy-expenditure and exercise physiology. While several equations have been offered for estimating total body components from MRI slices, no reliable and tested method exists for CT scans. For the first time, body composition data was derived from 41 high-resolution whole-body CT scans. From these data, we defined equations for estimating volumes and masses of total body AT and LT from corresponding tissue areas measured in selected CT scan slices. METHODS: We present a new semi-automatic approach to defining the density cutoff between adipose tissue (AT) and lean tissue (LT) in such material. An intra-class correlation coefficient (ICC) was used to validate the method. The equations for estimating the whole-body composition volume and mass from areas measured in selected slices were modeled with ordinary least squares (OLS) linear regressions and support vector machine regression (SVMR). RESULTS AND DISCUSSION: The best predictive equation for total body AT volume was based on the AT area of a single slice located between the 4th and 5th lumbar vertebrae (L4-L5) and produced lower prediction errors (|PE| = 1.86 liters, %PE = 8.77) than previous equations also based on CT scans. The LT area of the mid-thigh provided the lowest prediction errors (|PE| = 2.52 liters, %PE = 7.08) for estimating whole-body LT volume. We also present equations to predict total body AT and LT masses from a slice located at L4-L5 that resulted in reduced error compared with the previously published equations based on CT scans. The multislice SVMR predictor gave the theoretical upper limit for prediction precision of volumes and cross-validated the results.

Body mass estimation from the skeleton: An evaluation of 11 methods.

Estimating an individual body mass (BM) from the skeleton is a challenge for forensic anthropology. However, identifying someone's BMI (Body Mass Index) category, i.e. underweight, normal, overweight or obese, could contribute to identification. Individual BM is also known to influence the age-at-death estimation from the skeleton. Several methods are regularly used by both archaeologists and forensic practitioners to estimate individual BM. The most commonly used methods are based on femoral head breadth, or stature and bi-iliac breadth. However, those methods have been created from mean population BMs and are therefore meant to estimate the average BM of a population. Being that they are based on individual BM data and estimated femoral cortical areas, the newest published methods are supposed to be more accurate. We evaluated the accuracy and reliability of the most commonly used and most recent BM estimation methods (n=11) on a sample of 64 individuals. Both sexes and all BMI categories are represented, as well as a wide range of BM. Ages in this sample range from 20 to 87 years of age. Absolute and real differences between actual BM and estimated BM were assessed; they determined the accuracy for individual BM estimation and for average BM estimation of a population, respectively. The proportion of the sample whose estimated BM falls within ±10% and ±20% of their actual BM determines the reliability of the methods in our sample for, respectively, individual BM estimation and average BM of a population. The tested methods result in an absolute difference of 11kg-26kg±10kg with regards to prediction of individuals actual BM. The real differences are very variable from method to method, ranging from -14kg to 25kg. None of the tested methods is able to estimate BM of half of the sample within ±10% of their actual BM but most of them can estimate BM of more than half of the sample within ±20% of their actual BM. The errors increase with increasing BM, demonstrating a bias in all the methods. No bone variable tested correlated with BM. BMI categories were correctly predicted for less than 50% of the sample in most cases. In conclusion, our study demonstrates that the 11 methods tested are not suited for estimating individual BM or for predicting BMI categories. However, they are accurate and reliable enough for estimating the average BM of a population.