3D quantitative comparative analysis of long bone diaphysis variations in microanatomy and cross-sectional geometry

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.

Locomotor and postural diversity among reptiles viewed through the prism of femoral microanatomy: Palaeobiological implications for some Permian and Mesozoic taxa

The water-to-land transition by the first tetrapod vertebrates represents a key stage in their evolution. Selection pressures exerted by this new environment on animals led to the emergence of new locomotor and postural strategies that favoured access to different ecological niches and contributed to their evolutionary success. Today, amniotes show great locomotor and postural diversity, particularly among Reptilia, whose extant representatives include parasagittally locomoting erect and crouched bipeds (birds), sub-parasagittal 'semi-erect' quadrupeds (crocodylians) and sprawling quadrupeds (squamates and turtles). But the different steps leading to such diversity remain enigmatic and the type of locomotion adopted by many extinct species raises questions. This is notably the case of certain Triassic taxa such as Euparkeria and Marasuchus. The exploration of the bone microanatomy in reptiles could help to overcome these uncertainties. Indeed, this locomotor and postural diversity is accompanied by great microanatomical disparity. On land, the bones of the appendicular skeleton support the weight of the body and are subject to multiple constraints that partly shape their external and internal morphology. Here we show how microanatomical parameters measured in cross-section, such as bone compactness or the position of the medullocortical transition, can be related to locomotion. We hypothesised that this could be due to variations in cortical thickness. Using statistical methods that take phylogeny into account (phylogenetic flexible discriminant analyses), we develop different models of locomotion from a sample of femur cross-sections from 51 reptile species. We use these models to infer locomotion and posture in 7 extinct reptile taxa for which they remain debated or not fully clear. Our models produced reliable inferences for taxa that preceded and followed the quadruped/biped and sprawling/erect transitions, notably within the Captorhinidae and Dinosauria. For taxa contemporary with these transitions, such as Terrestrisuchus and Marasuchus, the inferences are more questionable. We use linear models to investigate the effect of body mass and functional ecology on our inference models. We show that body mass seems to significantly impact our model predictions in most cases, unlike the functional ecology. Finally, we illustrate how taphonomic processes can impact certain microanatomical parameters, especially the eccentricity of the section, while addressing some other potential limitations of our methods. Our study provides insight into the evolution of enigmatic locomotion in various early reptiles. Our models and methods could be used by palaeontologists to infer the locomotion and posture in other extinct reptile taxa, especially when considered in combination with other lines of evidence.

How does bone microanatomy and musculature covary? An investigation in the forelimb of two species of martens (Martes foina, Martes martes)

The long bones and associated musculature play a prominent role in the support and movement of the body and are expected to reflect the associated mechanical demands. But in addition to the functional response to adaptive changes, the conjoined effects of phylogenetic, structural and developmental constraints also shape the animal's body. In order to minimise the effect of the aforementioned constraints and to reveal the biomechanical adaptations in the musculoskeletal system to locomotor mode, we here study the forelimb of two closely related martens: the arboreal pine marten (Martes martes) and the more terrestrial stone marten (Martes foina), focusing on their forelimb muscle anatomy and long bone microanatomy; and, especially, on their covariation. To do so, we quantified muscle data and bone microanatomical parameters and created 3D and 2D maps of the cortical thickness distribution for the three long bones of the forelimb. We then analysed the covariation of muscle and bone data, both qualitatively and quantitatively. Our results reveal that species-specific muscular adaptations are not clearly reflected in the microanatomy of the bones. Yet, we observe a global thickening of the bone cortex in the radius and ulna of the more arboreal pine marten, as well a stronger flexor muscle inserting on its elbow. We attribute these differences to variation in their locomotor modes. Analyses of our 2D maps revealed a shift of cortical thickness distribution pattern linked to ontogeny, rather than species-specific patterns. We found that although intraspecific variation is not negligible, species distinction was possible when taking muscular and bone microanatomical data into consideration. Results of our covariation analyses suggest that the muscle-bone correlation is linked to ontogeny rather than to muscular strength at zones of insertion. Indeed, if we find a correlation between cortical thickness distribution and the strength of some muscles in the humerus, that is not the case for the others and in the radius and ulna. Cortical thickness distribution appears rather linked to bone contact zones and ligament insertions in the radius and ulna, and to some extent in the humerus. We conclude that inference on muscle from bone microanatomy is possible only for certain muscles in the humerus.

Sesamoid bones also show functional adaptation in their microanatomy-The example of the patella in Perissodactyla

The patella is the largest sesamoid bone of the skeleton. It is strongly involved in the knee, improving output force and velocity of the knee extensors, and thus plays a major role in locomotion and limb stability. However, the relationships between its structure and functional constraints, that would enable a better understanding of limb bone functional adaptations, are poorly known. This contribution proposes a comparative analysis, both qualitative and quantitative, of the microanatomy of the whole patella in perissodactyls, which show a wide range of morphologies, masses, and locomotor abilities, in order to investigate how the microanatomy of the patella adapts to evolutionary constraints. The inner structure of the patella consists of a spongiosa surrounded by a compact cortex. Contrary to our expectations, there is no increase in compactness with bone size, and thus body size and weight, but only an increase in the tightness of the spongiosa. No particular thickening of the cortex associated with muscle insertions is noticed but a strong thickening is observed anteriorly at about mid-length, where the strong intermediate patellar ligament inserts. The trabeculae are mainly oriented perpendicularly to the posterior articular surface, which highlights that the main stress is anteroposteriorly directed, maintaining the patella against the femoral trochlea. Conversely, anteriorly, trabeculae are rather circumferentially oriented, following the insertion of the patellar ligament and, possibly also, of the quadriceps tendon. A strong variation is observed among perissodactyl families but also intraspecifically, which is in accordance with previous studies suggesting a higher variability in sesamoid bones. Clear trends are nevertheless observed between the three families. Equids have a much thinner cortex than ceratomorphs. Rhinos and equids, both characterized by a development of the medial border, show an increase in trabecular density laterally suggesting stronger stresses laterally. The inner structure in tapirs is more homogeneous despite the absence of medial development of the medial border with no "compensation" of the inner structure, which suggests different stresses on their knees associated with a different morphology of their patellofemoral joint.

Functional Morphology and Morphological Diversification of Hind Limb Cross-Sectional Traits in Mustelid Mammals

Locomotor habits in mammals are strongly tied to limb bones’ lengths, diameters, and proportions. By comparison, fewer studies have examined how limb bone cross-sectional traits relate to locomotor habit. Here, we tested whether climbing, digging, and swimming locomotor habits reflect biomechanically meaningful differences in three cross-sectional traits rendered dimensionless— cross-sectional area (CSA), second moments of area (SMA), and section modulus (MOD)—using femora, tibiae, and fibulae of 28 species of mustelid. CSA and SMA represent resistance to axial compression and bending, respectively, whereas MOD represents structural strength. Given the need to counteract buoyancy in aquatic environments and soil’s high density, we predicted that natatorial and fossorial mustelids have higher values of cross-sectional traits. For all three traits, we found that natatorial mustelids have the highest values, followed by fossorial mustelids, with both of these groups significantly differing from scansorial mustelids. However, phylogenetic relatedness strongly influences diversity in cross-sectional morphology, as locomotor habit strongly correlates with phylogeny. Testing whether hind limb bone cross-sectional traits have evolved adaptively, we fit Ornstein–Uhlenbeck (OU) and Brownian motion (BM) models of trait diversification to cross-sectional traits. The cross-sectional traits of the femur, tibia, and fibula appear to have, respectively, diversified under a multi-rate BM model, a single rate BM model, and a multi-optima OU model. In light of recent studies on mustelid body size and elongation, our findings suggest that the mustelid body plan—and perhaps that of other mammals—is likely the sum of a suite of traits evolving under different models of trait diversification.

What about limb long bone nutrient canal(s)? - a 3D investigation in mammals

The nutrient arteries, located in the long bone diaphysis, are the major blood supply to long bones, especially during the early phases of growth and ossification. Their intersection with the central axis of the medullary area corresponds to the ossification center, and their opening on the outer bone surface to the nutrient foramen. Nutrient arteries/foramen have essentially been analyzed in humans, and only to a much lesser extent in a few mammals. Some studies have taken measurements of the nutrient foramen; others have investigated the shape and orientation of the nutrient canals, although only partially. No studies have analyzed the nutrient canal in three dimensions inside the bone and the relationships between nutrient foramen, nutrient canal, growth, and physiology require further investigation. The current study proposes to investigate in three dimensions the shape of the nutrient canal in stylopod bones of various mammals. Qualitative and quantitative parameters are defined to discuss the diversity in, for example, morphology, orientation, and diameter encountered, resorting to two different datasets to maximize differences within mammals and then analyze variation within morphologically and phylogenetically closer taxa. This study highlights a strong intraspecific variation for various parameters, with limited biological signal, but also shows trends. It notably provides evidence that canals are generally more numerous and relatively thinner in less elongated bones. Moreover, it shows that the growth center is located distally in the humerus and proximally in the femur, and that the canals are essentially oriented towards the faster growing end, so that the nutrient foramen does not indicate the location of the growth center. This result seems general in mammals but cannot be generalized outside of Mammalia. Further analyses of the features of nutrient arteries in reptiles are required to make comparisons with the trends observed in mammals.

Overall Bone Structure as Assessed by Slice-by-Slice Profile

Quantifying the inner structure of bones is central to various analyses dealing with the phenotypic evolution of animals with an ossified skeleton. Computed tomography allows to assess the repartition of bone tissue within an entire skeletal element. Two parameters of importance for such analyses are the global compactness (Cg) and total cross-sectional area (Tt.Ar). However, no open-source, time-efficient methods are available to acquire these parameters for whole bones. A methodology to assess the variation of these parameters along a profile following one of the studied bone’s anatomical axes is also wanting. Here I present an ImageJ macro and associated R script to automatically acquire Cg and Tt.Ar along an axis of the skeletal element of interest using a slice-by-slice approach. No manual segmentation is required and several bones can be present on the analysed scan, as long as the bone of interest is isolated and the largest element on each slice. While some bias might be involved by the automatic acquisition, semi-automatic slice exclusion and correction procedures can be used to efficiently account for it. As a test case, µCT-data was gathered for the mid-lumbar vertebra of over 70 mammals. The two evaluated correction procedures proved to perform equally well, with a slight advantage for the one relying on the exclusion of local outliers. The presented macro allows to efficiently build a dataset concerned with the quantification of bone inner structure. The code being readily available, further improvement of the methodology and adjustment to particular needs can be easily performed.

Morphological diversification of biomechanical traits: mustelid locomotor specializations and the macroevolution of long bone cross-sectional morphology

BACKGROUND

Morphological diversity of limb bone lengths, diameters, and proportions in mammals is known to vary strongly with locomotor habit. It remains less well known how different locomotor habits are correlated with cross-sectional traits of the limb skeleton, such as cross-sectional area (CSA), second moments of area (SMA), and section modulus (MOD) and whether these traits have evolved adaptively. CSA and SMA represent the bone's resistance to axial compression and bending, respectively, whereas MOD represents bone structural strength related to shape. Sampling 28 species of mustelids, a carnivoran lineage with diverse locomotor habits, we tested for differences in humeral, radial, and ulnar cross-sectional traits among specialists for climbing, digging, and swimming, in addition to generalists. Given that the limbs of digging specialists function in the dense substance of soil, and that swimming specialists need to counteract buoyancy, we predicted that these mustelids with these specializations should have the greatest values of cross-sectional traits.

RESULTS

We analyzed cross-sectional traits (calculated via μCT scanning and rendered dimensionless) in 5% increments along a bone's length and found significant differences among locomotor habits, though differences in ulnar cross-sectional traits were fewer than those for the humerus and radius. Swimming specialists had the greatest values of cross-sectional traits, followed by digging specialists. Climbing specialists had the lowest values of cross-sectional traits. However, phylogenetic affinity underlies these results. Fitting models of trait evolution to CSA and SMA revealed that a multi-rate Brownian motion model and a multi-optima Ornstein-Uhlenbeck model are the best-fitting models of evolution for these traits. However, inspection of α-values uncovered that many of the OU models did not differ from a Brownian motion model.

CONCLUSIONS

Within Mustelidae, differences in limb function and locomotor habit influence cross-sectional traits in ways that produce patterns that may diverge from adaptive patterns exhibited by external traits (e.g., bone lengths) of the mammalian limb skeleton. These results suggest that not all the traits of a single organ evolve under a single evolutionary process and that models of trait evolution should be fit to a range of traits for a better understanding of the evolution of the mammalian locomotor system.