Bone microstructure in finite element modeling: the functional role of trabeculae in the femoral head of Sciurus vulgaris

How compactness affects long bone resistance to compression-An investigation into the rhinoceros humerus

The functional signal of bone internal structure has been widely studied. Isolated form-function relationships have often been assumed from the observation of presumed morphofunctional relationships, but have never been truly tested. Indeed, distinct bone microanatomical feature co-evolve in response to various constraints that are difficult to detangle. This study tested for the first time the impact of various microanatomical parameters taken one by one, plus some in pairs, on bone strength under compression using biomechanical modelling. We carried out finite element analyses on humerus models, obtained from a white rhinoceros, with different heterogeneous internal structures, and analysed the magnitude and distribution of von Mises stresses. These tests validated earlier hypotheses of form-function relationships about the greater resistance to compression provided by the thickening of the cortex and the filling of the medullary area by trabecular bone and highlighted the stronger impact of increasing trabecular bone compactness than of avoiding an open medullary cavity. By making it possible to estimate the relative impact of each parameter and of combinations of microanatomical features, they also showed the more limited impact of the trabecular bone compactness in the epiphyses to resist compression, and the fact that microanatomical changes of opposite but of similar amplitude impact can compensate each other, but that the impact of the sum of two negative microanatomical changes far exceeds the sum of the impacts of each of the two changes taken separately. These results contribute to a better understanding of bone adaptation and form-function relationships so that they later can be used with confidence for palaeobiological inferences on fossil specimens, contributing to a better understanding of skeletal evolution during the evolutionary history of vertebrates. They also highlight the potential of taking internal structure into account in the bone biomechanical analyses. In addition, they can be used in bioinspiration to design resistant structures subjected to compression.

A new biomechanical model of the mammal jaw based on load path analysis

The primary function of the tetrapod jaw is to transmit jaw muscle forces to bite points. The routes of force transfer in the jaw have never been studied but can be quantified using load paths - the shortest, stiffest routes from regions of force application to support constraints. Here, we use load path analysis to map force transfer from muscle attachments to bite point and jaw joint, and to evaluate how different configurations of trabecular and cortical bone affect load paths. We created three models of the mandible of the Virginia opossum, Didelphis virginiana, each with a cortical bone shell, but with different material properties for the internal spaces: (1) a cortical-trabecular model, in which the interior space is modeled with bulk properties of trabecular bone; (2) a cortical-hollow model, in which trabeculae and mandibular canal are modeled as hollow; and (3) a solid-cortical model, in which the interior is modeled as cortical bone. The models were compared with published in vivo bite force and bone strain data, and the load paths calculated for each model. The trabecular model, which is preferred because it most closely approximates the actual morphology, was best validated by in vivo data. In all three models, the load path was confined to cortical bone, although its route within the cortex varied depending on the material properties of the inner model. Our analysis shows that most of the force is transferred through the cortical, rather than trabecular bone, and highlights the potential of load path analysis for understanding form-function relationships in the skeleton.

Small skeletons show size-specific scaling: an exploration of allometry in the mammalian lumbar spine

Studies of vertebrate bone biomechanics often focus on skeletal adaptations at upper extremes of body mass, disregarding the importance of skeletal adaptations at lower extremes. Yet mammals are ancestrally small and most modern species have masses under 5 kg, so the evolution of morphology and function at small size should be prioritized for understanding how mammals subsist. We examined allometric scaling of lumbar vertebrae in the small-bodied Philippine endemic rodents known as cloud rats, which vary in mass across two orders of magnitude (15.5 g-2700 g). External vertebral dimensions scale with isometry or positive allometry, likely relating to body size and nuances in quadrupedal posture. In contrast to most mammalian trabecular bone studies, bone volume fraction and trabecular thickness scale with positive allometry and isometry, respectively. It is physiologically impossible for these trends to continue to the upper extremes of mammalian body size, and we demonstrate a fundamental difference in trabecular bone allometry between large- and small-bodied mammals. These findings have important implications for the biomechanical capabilities of mammalian bone at small body size; for the selective pressures that govern skeletal evolution in small mammals; and for the way we define 'small' and 'large' in the context of vertebrate skeletons.

Diversity and function of the fused anuran radioulna

In tetrapods, fusion between elements of the appendicular skeleton is thought to facilitate rapid movements during running, flying, and jumping. Although such fusion is widespread, frogs stand out because adults of all living species exhibit fusion of the zeugopod elements (radius and ulna, tibia and fibula), regardless of jumping ability or locomotor mode. To better understand what drives the maintenance of limb bone fusion in frogs, we use finite element modeling methods to assess the functional consequences of fusion in the anuran radioulna, the forearm bone of frogs that is important to both locomotion and mating behavior (amplexus). Using CT scans of museum specimens, measurement tools, and mesh‐editing software, we evaluated how different degrees of fusion between the radius and ulna affect the von Mises stress and bending resistance of the radioulna in three loading scenarios: landing, amplexus, and long‐axis loading conditions. We find that the semi‐fused state observed in the radioulna exhibits less von Mises stress and more resistance to bending than unfused or completely fused models in all three scenarios. Our results suggest that radioulna morphology is optimized to minimize von Mises stress across different loading regimes while also minimizing volume. We contextualize our findings in an evaluation of the diversity of anuran radioulnae, which reveals unique, permanent pronation of the radioulna in frogs and substantial variation in wall thickness. This work provides new insight into the functional consequences of limb bone fusion in anuran evolution.

Accuracy of bone segmentation and surface generation strategies analyzed by using synthetic CT volumes

Different kinds of bone measurements are commonly derived from computed-tomography (CT) volumes to answer a multitude of questions in biology and related fields. The underlying steps of bone segmentation and, optionally, polygon surface generation are crucial to keep the measurement error small. In this study, the performance of different, easily accessible segmentation techniques (global thresholding, automatic local thresholding, weighted random walk, neural network, and watershed) and surface generation approaches (different algorithms combined with varying degrees of simplification) was analyzed and recommendations for minimizing inaccuracies were derived. The different approaches were applied to synthetic CT volumes for which the correct segmentation and surface geometry were known. The most accurate segmentations of the synthetic volumes were achieved by setting a case-specific window to the gray value histogram and subsequently applying automatic local thresholding with appropriately chosen thresholding method and radius. Surfaces generated by the Amira® module Generate Lego Surface in combination with careful surface simplification were the most accurate. Surfaces with sub-voxel accuracy were obtained even for synthetic CT volumes with low contrast-to-noise ratios. Segmentation trials with real CT volumes supported the findings. Very accurate segmentations and surfaces can be derived from CT volumes by using readily accessible software packages. The presented results and derived recommendations will help to reduce the measurement error in future studies. Furthermore, the demonstrated strategies for assessing segmentation and surface qualities can be adopted to quantify the performance of new segmentation approaches in future studies.

Ecomorphology of the pectoral girdle in anurans (Amphibia, Anura): Shape diversity and biomechanical considerations

Frogs and toads (Lissamphibia: Anura) show a diversity of locomotor modes that allow them to inhabit a wide range of habitats. The different locomotor modes are likely to be linked to anatomical specializations of the skeleton within the typical frog Bauplan. While such anatomical adaptations of the hind limbs and the pelvic girdle are comparably well understood, the pectoral girdle received much less attention in the past. We tested for locomotor-mode-related shape differences in the pectoral girdle bones of 64 anuran species by means of micro-computed-tomography-based geometric morphometrics. The pectoral girdles of selected species were analyzed with regard to the effects of shape differences on muscle moment arms across the shoulder joint and stress dissipation within the coracoid. Phylogenetic relationships, size, and locomotor behavior have an effect on the shape of the pectoral girdle in anurans, but there are differences in the relative impact of these factors between the bones of this skeletal unit. Remarkable shape diversity has been observed within locomotor groups indicating many-to-one mapping of form onto function. Significant shape differences have mainly been related to the overall pectoral girdle geometry and the shape of the coracoid. Most prominent shape differences have been found between burrowing and nonburrowing species with headfirst and backward burrowing species significantly differing from one another and from the other locomotor groups. The pectoral girdle shapes of burrowing species have generally larger moment arms for (simulated) humerus retractor muscles across the shoulder joint, which might be an adaptation to the burrowing behavior. The mechanisms of how the moment arms were enlarged differed between species and were associated with differences in the reaction of the coracoid to simulated loading by physiologically relevant forces.

Deciphering an extreme morphology: bone microarchitecture of the hero shrew backbone (Soricidae: Scutisorex)

Biological structures with extreme morphologies are puzzling because they often lack obvious functions and stymie comparisons to homologous or analogous features with more typical shapes. An example of such an extreme morphotype is the uniquely modified vertebral column of the hero shrew Scutisorex, which features numerous accessory intervertebral articulations and massively expanded transverse processes. The function of these vertebral structures is unknown, and it is difficult to meaningfully compare them to vertebrae from animals with known behavioural patterns and spinal adaptations. Here, we use trabecular bone architecture of vertebral centra and quantitative external vertebral morphology to elucidate the forces that may act on the spine of Scutisorex and that of another large shrew with unmodified vertebrae (Crocidura goliath). X-ray micro-computed tomography (µCT) scans of thoracolumbar columns show that Scutisorex thori is structurally intermediate between C. goliath and S. somereni internally and externally, and both Scutisorex species exhibit trabecular bone characteristics indicative of higher in vivo axial compressive loads than C. goliath. Under compressive load, Scutisorex vertebral morphology is adapted to largely restrict bending to the sagittal plane (flexion). Although these findings do not solve the mystery of how Scutisorex uses its byzantine spine in vivo, our work suggests potentially fruitful new avenues of investigation for learning more about the function of this perplexing structure.