Mitochondrial Fission Promotes the Continued Clearance of Apoptotic Cells by Macrophages.

Clearance of apoptotic cells (ACs) by phagocytes (efferocytosis) prevents post-apoptotic necrosis and dampens inflammation. Defective efferocytosis drives important diseases, including atherosclerosis. For efficient efferocytosis, phagocytes must be able to internalize multiple ACs. We show here that uptake of multiple ACs by macrophages requires dynamin-related protein 1 (Drp1)-mediated mitochondrial fission, which is triggered by AC uptake. When mitochondrial fission is disabled, AC-induced increase in cytosolic calcium is blunted owing to mitochondrial calcium sequestration, and calcium-dependent phagosome formation around secondarily encountered ACs is impaired. These defects can be corrected by silencing the mitochondrial calcium uniporter (MCU). Mice lacking myeloid Drp1 showed defective efferocytosis and its pathologic consequences in the thymus after dexamethasone treatment and in advanced atherosclerotic lesions in fat-fed Ldlr-/- mice. Thus, mitochondrial fission in response to AC uptake is a critical process that enables macrophages to clear multiple ACs and to avoid the pathologic consequences of defective efferocytosis in vivo.

Activity-driven local ATP synthesis is required for synaptic function.

Cognitive function is tightly related to metabolic state, but the locus of this control is not well understood. Synapses are thought to present large ATP demands; however, it is unclear how fuel availability and electrical activity impact synaptic ATP levels and how ATP availability controls synaptic function. We developed a quantitative genetically encoded optical reporter of presynaptic ATP, Syn-ATP, and find that electrical activity imposes large metabolic demands that are met via activity-driven control of both glycolysis and mitochondrial function. We discovered that the primary source of activity-driven metabolic demand is the synaptic vesicle cycle. In metabolically intact synapses, activity-driven ATP synthesis is well matched to the energetic needs of synaptic function, which, at steady state, results in ∼10(6) free ATPs per nerve terminal. Despite this large reservoir of ATP, we find that several key aspects of presynaptic function are severely impaired following even brief interruptions in activity-stimulated ATP synthesis.

iATPSnFR2: A high-dynamic-range fluorescent sensor for monitoring intracellular ATP.

We developed a significantly improved genetically encoded quantitative adenosine triphosphate (ATP) sensor to provide real-time dynamics of ATP levels in subcellular compartments. iATPSnFR2 is a variant of iATPSnFR1, a previously developed sensor that has circularly permuted superfolder green fluorescent protein (GFP) inserted between the ATP-binding helices of the ε-subunit of a bacterial F0-F1 ATPase. Optimizing the linkers joining the two domains resulted in a ~fivefold to sixfold improvement in the dynamic range compared to the previous-generation sensor, with excellent discrimination against other analytes, and affinity variants varying from 4 µM to 500 µM. A chimeric version of this sensor fused to either the HaloTag protein or a suitable spectrally separated fluorescent protein provides an optional ratiometric readout allowing comparisons of ATP across cellular regions. Subcellular targeting the sensor to nerve terminals reveals previously uncharacterized single-synapse metabolic signatures, while targeting to the mitochondrial matrix allowed direct quantitative probing of oxidative phosphorylation dynamics.

A 3.8-million-year-old hominin cranium from Woranso-Mille, Ethiopia.

The cranial morphology of the earliest known hominins in the genus Australopithecus remains unclear. The oldest species in this genus (Australopithecus anamensis, specimens of which have been dated to 4.2-3.9 million years ago) is known primarily from jaws and teeth, whereas younger species (dated to 3.5-2.0 million years ago) are typically represented by multiple skulls. Here we describe a nearly complete hominin cranium from Woranso-Mille (Ethiopia) that we date to 3.8 million years ago. We assign this cranium to A. anamensis on the basis of the taxonomically and phylogenetically informative morphology of the canine, maxilla and temporal bone. This specimen thus provides the first glimpse of the entire craniofacial morphology of the earliest known members of the genus Australopithecus. We further demonstrate that A. anamensis and Australopithecus afarensis differ more than previously recognized and that these two species overlapped for at least 100,000 years-contradicting the widely accepted hypothesis of anagenesis.

Trabecular bone ontogeny tracks neural development and life history among humans and non-human primates.

Trabecular bone-the spongy bone inside marrow cavities-adapts to its mechanical environment during growth and development. Trabecular structure can therefore be interpreted as a functional record of locomotor behavior in extinct vertebrates. In this paper, we expand upon traditional links between form and function by situating ontogenetic trajectories of trabecular bone in four primate species into the broader developmental context of neural development, locomotor control, and ultimately life history. Our aim is to show that trabecular bone structure provides insights into ontogenetic variation in locomotor loading conditions as the product of interactions between increases in body mass and neuromuscular maturation. Our results demonstrate that age-related changes in trabecular bone volume fraction (BV/TV) are strongly and linearly associated with ontogenetic changes in locomotor kinetics. Age-related variation in locomotor kinetics and BV/TV is in turn strongly associated with brain and body size growth in all species. These results imply that age-related variation in BV/TV is a strong proxy for both locomotor kinetics and neuromuscular maturation. Finally, we show that distinct changes in the slope of age-related variation in bone volume fraction correspond to the age of the onset of locomotion and the age of locomotor maturity. Our findings compliment previous studies linking bone development to locomotor mechanics by providing a fundamental link to brain development and life history. This implies that trabecular structure of fossil subadults can be a proxy for the rate of neuromuscular maturation and major life history events like locomotor onset and the achievement of adult-like locomotor repertoires.

Brain energy metabolism: A roadmap for future research.

Although we have learned much about how the brain fuels its functions over the last decades, there remains much still to discover in an organ that is so complex. This article lays out major gaps in our knowledge of interrelationships between brain metabolism and brain function, including biochemical, cellular, and subcellular aspects of functional metabolism and its imaging in adult brain, as well as during development, aging, and disease. The focus is on unknowns in metabolism of major brain substrates and associated transporters, the roles of insulin and of lipid droplets, the emerging role of metabolism in microglia, mysteries about the major brain cofactor and signaling molecule NAD+, as well as unsolved problems underlying brain metabolism in pathologies such as traumatic brain injury, epilepsy, and metabolic downregulation during hibernation. It describes our current level of understanding of these facets of brain energy metabolism as well as a roadmap for future research.

The control of release probability at nerve terminals.

Exocytosis is a fundamental membrane fusion process by which the soluble or membrane-associated cargoes of a secretory vesicle are delivered to the extracellular milieu or the cell surface. While essential for all organs, the brain relies on a specialized form of exocytosis to mediate information flow throughout its vast circuitry. Neurotransmitter-laden synaptic vesicles fuse with the plasma membrane on cue with astonishing speed in a probabilistic process that is both tightly regulated and capable of a fascinating array of plasticities. Here, we examine progress in the molecular understanding of synaptic vesicle fusion and its control.

Financial impact of donation after circulatory death heart transplantation: A single-center analysis.

INTRODUCTION: Clinical success of donation after circulatory death (DCD) heart transplantation is leading to growing adoption of this technique. In comparison to procurement from a brain-dead donor, DCD requires additional resources. The economic impact of DCD heart transplantation from the hospital perspective is not well known. METHODS: We compared the financial data of patients who received DCD allografts to those who received a DBD organ at our institution from January 1, 2021 to December 31, 2022. We also compared the cost of ex-situ machine perfusion to in-situ organ perfusion employed during DCD recovery. RESULTS: We performed 58 DBD and 22 DCD heart-alone transplantations during the study period. Out of 22 DCD grafts, 16 were recovered with thoracoabdominal normothermic regional perfusion (TA-NRP) and six with direct procurement followed by normothermic machine perfusion (DP-NMP). The contribution margin per case for DBD versus DCD was $234,362 and $235,440 (P = .72). The direct costs did not significantly differ between the two groups ($171,949 and 186,250; P = .49). In comparing the two methods of procuring hearts from DCD donors, the direct cost of TA-NRP was $155,955 in comparison to $223,399 for DP-NMP (P = .21). This difference translated into a clinically meaningful but not statistically significant greater contribution margin for TA-NRP ($242, 657 vs. $175,768; P = .34). CONCLUSIONS: Our data showed that the adoption of DCD procurement did not have a negative financial impact on the contribution margin in our institution. Programs considering starting DCD heart transplantation, and those who are currently performing DCD procurement should evaluate their own financial situation.

Elite rhetoric can undermine democratic norms.

Democratic stability depends on citizens on the losing side accepting election outcomes. Can rhetoric by political leaders undermine this norm? Using a panel survey experiment, we evaluate the effects of exposure to multiple statements from former president Donald Trump attacking the legitimacy of the 2020 US presidential election. Although exposure to these statements does not measurably affect general support for political violence or belief in democracy, it erodes trust and confidence in elections and increases belief that the election is rigged among people who approve of Trump's job performance. These results suggest that rhetoric from political elites can undermine respect for critical democratic norms among their supporters.

Molecular circuitry of endocytosis at nerve terminals.

Presynaptic terminals are specialized compartments of neurons responsible for converting electrical signals into secreted chemicals. This self-renewing process of chemical synaptic transmission is accomplished by the calcium-triggered fusion of neurotransmitter-containing vesicles with the plasma membrane and subsequent retrieval and recycling of vesicle components. Whereas the release of neurotransmitters has been studied for over 50 years, the process of synaptic vesicle endocytosis has remained much more elusive. The advent of imaging techniques suited to monitor membrane retrieval at presynaptic terminals and the discovery of the molecules that orchestrate endocytosis have revolutionized our understanding of this critical trafficking event.

Growth and development of trabecular structure in the calcaneus of Japanese macaques (Macaca fuscata) reflects locomotor behavior, life history, and neuromuscular development.

Bone structure dynamically adapts to its mechanical environment throughout ontogeny by altering the structure of trabecular bone, the three-dimensional mesh-like structure found underneath joint surfaces. Trabecular structure, then, can provide a record of variation in loading directions and magnitude; and in ontogenetic samples, it can potentially be used to track developmental shifts in limb posture. We aim to broaden the analysis of trabecular bone ontogeny by incorporating interactions between ontogenetic variation in locomotor repertoire, neuromuscular maturation, and life history. We examine the associations between these variables and age-related variation in trabecular structure in the calcaneus of Japanese macaques (Macaca fuscata). We used high-resolution micro-computed tomography scanning to image the calcaneus in a cross-sectional sample of 34 juvenile M. fuscata aged between 0 and 7 years old at the Primate Research Institute, Japan. We calculated whole bone averages of standard trabecular properties and generated whole-bone morphometric maps of bone volume fraction and Young's modulus. Trabecular structure becomes increasingly heterogeneous in older individuals. Bone volume fraction (BV/total volume [TV]) decreases during the first month of life and increases afterward, coinciding with the onset of independent locomotion in M. fuscata. At birth, primary Young's modulus is oriented orthogonal to the ossification center, but after locomotor onset bone structure becomes stiffest in the direction of joint surfaces and muscle attachments. Age-related variation in bone volume fraction is best predicted by an interaction between the estimated percentage of adult brain size, body mass, and locomotor onset. To explain our findings, we propose a model where interactions between age-related increases in body weight and maturation of the neuromuscular system alter the loading environment of the calcaneus, to which the internal trabecular structure dynamically adapts. This model cannot be directly tested based on our cross-sectional data. However, confirmation of the model by longitudinal experiments and in multiple species would show that trabecular structure can be used both to infer behavior from fossil morphology and serve as a valuable proxy for neuromuscular maturation and life history events like locomotor onset and the achievement of an adult-like gait. This approach could significantly expand our knowledge of the biology and behavior of fossil species.

Impact of Initial Warfarin Dosing on Time in Therapeutic Range for Postoperative Left Ventricular Assist Device Patients.

ABSTRACT: Initial warfarin dosing and time in therapeutic range (TTR) are poorly characterized for early post-operative left ventricular assist device (LVAD) patients. This study evaluated TTR after LVAD implantation compared between patients receiving low-dose (<3 mg) and high-dose (≥3 mg) warfarin. This single-center, retrospective analysis included 234 LVAD patients who received warfarin within 5 days of implantation. The primary outcome was TTR during the 5 days following first international normalized ratio (INR) ≥2 compared between low-dose and high-dose groups. Secondary outcomes were hospital and intensive care unit length of stay, time to first INR ≥2, TTR after first INR ≥2, and reinitiation of parenteral anticoagulation. No difference in TTR was detected between warfarin groups (57.2% vs. 62.7%, P = 0.13). Multivariable analysis did not detect any factors predictive of TTR during the primary outcome timeframe, but age and body mass index were associated with the warfarin dose. The low-dose group received a mean warfarin dose of 1.9 mg (±0.64 mg), and the high dose group received 4.34 mg (±1.38 mg). Cohort TTR during the primary outcome timeframe was 60.5% and 56.5% for hospitalization. The low-dose group had longer intensive care unit length of stay, shorter time to therapeutic INR, and more frequently reinitiated parenteral anticoagulation. Patients with recent LVAD implantation are complex and have diverse warfarin sensitivity factors, which did not allow for optimal warfarin dose detection, although half of all patients received doses between 2.04 mg and 4.33 mg. Individualized dosing should be used, adjusting for patient-specific factors such as age, body mass index, and drug interactions.

Protic Ionic Liquid Cation Alkyl Chain Length Effect on Lysozyme Structure.

Solvents that stabilize protein structures can improve and expand their biochemical applications, particularly with the growing interest in biocatalytic-based processes. Aiming to select novel solvents for protein stabilization, we explored the effect of alkylammonium nitrate protic ionic liquids (PILs)-water mixtures with increasing cation alkyl chain length on lysozyme conformational stability. Four PILs were studied, that is, ethylammonium nitrate (EAN), butylammonium nitrate (BAN), hexylammonium nitrate (HAN), and octylammonium nitrate (OAN). The surface tension, viscosity, and density of PIL-water mixtures at low to high concentrations were firstly determined, which showed that an increasing cation alkyl chain length caused a decrease in the surface tension and density as well as an increase in viscosity for all PIL solutions. Small-angle X-ray scattering (SAXS) was used to investigate the liquid nanostructure of the PIL solutions, as well as the overall size, conformational flexibility and changes to lysozyme structure. The concentrated PILs with longer alkyl chain lengths, i.e., over 10 mol% butyl-, 5 mol% hexyl- and 1 mol% octylammonium cations, possessed liquid nanostructures. This detrimentally interfered with solvent subtraction, and the more structured PIL solutions prevented quantitative SAXS analysis of lysozyme structure. The radius of gyration (Rg) of lysozyme in the less structured aqueous PIL solutions showed little change with up to 10 mol% of PIL. Kratky plots, SREFLEX models, and FTIR data showed that the protein conformation was maintained at a low PIL concentration of 1 mol% and lower when compared with the buffer solution. However, 50 mol% EAN and 5 mol% HAN significantly increased the Rg of lysozyme, indicating unfolding and aggregation of lysozyme. The hydrophobic interaction and liquid nanostructure resulting from the increased cation alkyl chain length in HAN likely becomes critical. The impact of HAN and OAN, particularly at high concentrations, on lysozyme structure was further revealed by FTIR. This work highlights the negative effect of a long alkyl chain length and high concentration of PILs on lysozyme structural stability.

New Pliocene hominin remains from the Leado Dido'a area of Woranso-Mille, Ethiopia.

Fossiliferous deposits at Woranso-Mille span the period when Australopithecus anamensis gave rise to Australopithecus afarensis (3.8-3.6 Ma) and encompass the core of the A. afarensis range (ca. 3.5-3.2 Ma). Within the latter period, fossils described to date include the intriguing but taxonomically unattributed Burtele foot, dentognathic fossils attributed to Australopithecus deyiremeda, and one specimen securely attributed to A. afarensis (the Nefuraytu mandible). These fossils suggest that at least one additional hominin lineage lived alongside A. afarensis in the Afar Depression. Here we describe a collection of hominin fossils from a new locality in the Leado Dido'a area of Woranso-Mille (LDD-VP-1). The strata in this area are correlated to the same chron as those in the Burtele area (C2An.3n; 3.59-3.33 Ma), and similar in age to the Maka Sands and the Basal through lower Sidi Hakoma Members of the Hadar Formation. We attribute all but one of the LDD hominin specimens to A. afarensis, based on diagnostic morphology of the mandible, maxilla, canines, and premolars. The LDD specimens generally fall within the range of variation previously documented for A. afarensis but increase the frequency of some rare morphological variants. However, one isolated M3 is extremely small, and its taxonomic affinity is currently unknown. The new observations support previous work on temporal trends in A. afarensis and demonstrate that the large range of variation accepted for this species is present even within a limited spatiotemporal range. The value added with this sample lies in its contribution to controlling for spatiotemporal differences among site samples in the A. afarensis hypodigm and its contemporaneity with non-A. afarensis specimens at Woranso-Mille.

Unique foot posture in Neanderthals reflects their body mass and high mechanical stress.

Neanderthal foot bone proportions and morphology are mostly indistinguishable from those of Homo sapiens, with the exception of several distinct Neanderthal features in the talus. The biomechanical implications of these distinct talar features remain contentious, fueling debate around the adaptive meaning of this distinctiveness. With the aim of clarifying this controversy, we test phylogenetic and behavioral factors as possible contributors, comparing tali of 10 Neanderthals and 81 H. sapiens (Upper Paleolithic and Holocene hunter-gatherers, agriculturalists, and postindustrial group) along with the Clark Howell talus (Omo, Ethiopia). Variation in external talar structures was assessed through geometric morphometric methods, while bone volume fraction and degree of anisotropy were quantified in a subsample (n = 45). Finally, covariation between point clouds of site-specific trabecular variables and surface landmark coordinates was assessed. Our results show that although Neanderthal talar external and internal morphologies were distinct from those of H. sapiens groups, shape did not significantly covary with either bone volume fraction or degree of anisotropy, suggesting limited covariation between external and internal talar structures. Neanderthal external talar morphology reflects ancestral retentions, along with various adaptations to high levels of mobility correlated to their presumably unshod hunter-gatherer lifestyle. This pairs with their high site-specific trabecular bone volume fraction and anisotropy, suggesting intense and consistently oriented locomotor loading, respectively. Relative to H.sapiens, Neanderthals exhibit differences in the talocrural joint that are potentially attributable to cultural and locomotor behavior dissimilarity, a talonavicular joint that mixes ancestral and functional traits, and a derived subtalar joint that suggests a predisposition for a pronated foot during stance phase. Overall, Neanderthal talar variation is attributable to mobility strategy and phylogenesis, while H. sapiens talar variation results from the same factors plus footwear. Our results suggest that greater Neanderthal body mass and/or higher mechanical stress uniquely led to their habitually pronated foot posture.

Microemulsion-Assisted Templating of Metal-Stabilized Poly(ethylene glycol) Nanoparticles.

Poly(ethylene glycol) (PEG) is well known to endow nanoparticles (NPs) with low-fouling and stealth-like properties that can reduce immune system clearance in vivo, making PEG-based NPs (particularly sub-100 nm) of interest for diverse biomedical applications. However, the preparation of sub-100 nm PEG NPs with controllable size and morphology is challenging. Herein, we report a strategy based on the noncovalent coordination between PEG-polyphenolic ligands (PEG-gallol) and transition metal ions using a water-in-oil microemulsion phase to synthesize sub-100 nm PEG NPs with tunable size and morphology. The metal-phenolic coordination drives the self-assembly of the PEG-gallol/metal NPs: complexation between MnII and PEG-gallol within the microemulsions yields a series of metal-stabilized PEG NPs, including 30-50 nm solid and hollow NPs, depending on the MnII/gallol feed ratio. Variations in size and morphology are attributed to the changes in hydrophobicity of the PEG-gallol/MnII complexes at varying MnII/gallol ratios based on contact angle measurements. Small-angle X-ray scattering analysis, which is used to monitor the particle size and intermolecular interactions during NP evolution, reveals that ionic interactions are the dominant driving force in the formation of the PEG-gallol/MnII NPs. pH and cytotoxicity studies, and the low-fouling properties of the PEG-gallol/MnII NPs confirm their high biocompatibility and functionality, suggesting that PEG polyphenol-metal NPs are promising systems for biomedical applications.

New species from Ethiopia further expands Middle Pliocene hominin diversity.

Middle Pliocene hominin species diversity has been a subject of debate over the past two decades, particularly after the naming of Australopithecus bahrelghazali and Kenyanthropus platyops in addition to the well-known species Australopithecus afarensis. Further analyses continue to support the proposal that several hominin species co-existed during this time period. Here we recognize a new hominin species (Australopithecus deyiremeda sp. nov.) from 3.3-3.5-million-year-old deposits in the Woranso-Mille study area, central Afar, Ethiopia. The new species from Woranso-Mille shows that there were at least two contemporaneous hominin species living in the Afar region of Ethiopia between 3.3 and 3.5 million years ago, and further confirms early hominin taxonomic diversity in eastern Africa during the Middle Pliocene epoch. The morphology of Au. deyiremeda also reinforces concerns related to dentognathic (that is, jaws and teeth) homoplasy in Plio-Pleistocene hominins, and shows that some dentognathic features traditionally associated with Paranthropus and Homo appeared in the fossil record earlier than previously thought.

Combinations of trabecular and cortical bone properties distinguish various loading modalities between athletes and controls.

OBJECTIVES: Variation in trabecular and cortical bone properties is often used to infer habitual behavior in the past. However, the structures of both types of bone are rarely considered together and may even contradict each other in functional interpretations. We examine trabecular and cortical bone properties in various athletes and sedentary controls to clarify the associations between combinations of cortical and trabecular bone properties and various loading modalities. MATERIALS AND METHODS: We compare trabecular and cortical bone properties using peripheral quantitative computed tomography scans of the tibia between groups of 83 male athletes (running, hockey, swimming, cricket) and sedentary controls using Bayesian multilevel models. We quantify midshaft cortical bone rigidity and area (J, CA), midshaft shape index (Imax/Imin), and mean trabecular bone mineral density (BMD) in the distal tibia. RESULTS: All groups show unique combinations of biomechanical properties. Cortical bone rigidity is high in sports that involve impact loading (cricket, running, hockey) and low in nonimpact loaded swimmers and controls. Runners have more anteroposteriorly elliptical midshafts compared to other groups. Interestingly, all athletes have greater trabecular BMD compared to controls, but do not differ credibly among each other. DISCUSSION: Results suggest that cortical midshaft hypertrophy is associated with impact loading while trabecular BMD is positively associated with both impact and nonimpact loading. Midshaft shape is associated with directionality of loading. Individuals from the different categories overlap substantially, but group means differ credibly, suggesting that nuanced group-level inferences of habitual behavior are possible when combinations of trabecular and cortical bone are analyzed.

Automated resolution independent method for comparing in vivo and dry trabecular bone.

OBJECTIVES: Variation in human trabecular bone morphology can be linked to habitual behavior, but it is difficult to investigate in vivo due to the radiation required at high resolution. Consequently, functional interpretations of trabecular morphology remain inferential. Here we introduce a method to link low- and high-resolution CT data from dry and fresh bone, enabling bone functional adaptation to be studied in vivo and results compared to the fossil and archaeological record. MATERIALS AND METHODS: We examine 51 human dry bone distal tibiae from Nile Valley and UK and two pig tibiae containing soft tissues. We compare low-resolution peripheral quantitative computed tomography (pQCT) parameters and high-resolution micro CT (µCT) in homologous single slices at 4% bone length and compare results to our novel Bone Ratio Predictor (BRP) method. RESULTS: Regression slopes between linear attenuation coefficients of low-resolution pQCT images and bone area/total area (BA/TA) of high-resolution µCT scans differ substantially between geographical subsamples, presumably due to diagenesis. BRP accurately predicts BA/TA (R2 = .97) and eliminates the geographic clustering. BRP accurately estimates BA/TA in pigs containing soft tissues (R2 = 0.98) without requiring knowledge of true density or phantom calibration of the scans. DISCUSSION: BRP allows automated comparison of image data from different image modalities (pQCT, µCT) using different energy settings, in archeological bone and wet specimens. The method enables low-resolution data generated in vivo to be compared with the fossil and archaeological record. Such experimental approaches would substantially improve behavioral inferences based on trabecular bone microstructure.

Using point clouds to investigate the relationship between trabecular bone phenotype and behavior: An example utilizing the human calcaneus.

OBJECTIVES: The objective of this study is to demonstrate a new method for analyzing trabecular bone volume fraction and degree of anisotropy in three dimensions. METHODS: We use a combination of automatic mesh registration, point-cloud correspondence registration, and P-value corrected univariate statistical tests to compare bone volume fraction and degree of anisotropy on a point by point basis across the entire calcaneus of two human groups with different subsistence strategies. RESULTS: We found that the patterns of high and low bone volume fraction and degree of anisotropy distribution between the Black Earth (hunter-gatherers) and Norris Farms (mixed-strategy agriculturalists) are very similar, but differ in magnitude. The hunter-gatherers exhibit higher levels of bone volume fraction and less anisotropic trabecular bone organization. Additionally, patterns of bone volume fraction and degree of anisotropy in the calcaneus correspond well with biomechanical expectations of relative forces experienced during walking and running. CONCLUSIONS: We conclude that comparing site-specific, localized differences in trabecular bone variables such as bone volume fraction and degree of anisotropy in three-dimensions is a powerful analytical tool. This method makes it possible to determine where similarities and differences between groups are located within the whole skeletal element of interest. The visualization of multiple variables also provides a way for researchers to see how the trabecular bone variables interact within the morphology, and allows for a more nuanced understanding of how they relate to one another and the broader mechanical environment.