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

Osteological profiling of femoral diaphysis and neck in aquatic, semiaquatic, and terrestrial carnivores and rodents: effects of body size and locomotor habits

The increased limb bone density documented previously for aquatic tetrapods has been proposed to be an adaptation to overcome buoyancy during swimming and diving. It can be achieved by increasing the amount of bone deposition or by reducing the amount of bone resorption, leading to cortical thickening, loss of medullary cavity, and compaction of trabecular bone. The present study examined the effects of locomotor habit, body size, and phylogeny on the densitometric, cross-sectional, and biomechanical traits of femoral diaphysis and neck in terrestrial, semiaquatic, and aquatic carnivores, and in terrestrial and semiaquatic rodents (12 species) by using peripheral quantitative computed tomography, three-point bending, and femoral neck loading tests. Groupwise differences were analyzed with the univariate generalized linear model and the multivariate linear discriminant analysis supplemented with hierarchical clustering. While none of the individual features could separate the lifestyles or species adequately, the combinations of multiple features produced very good or excellent classifications and clusterings. In the phocid seals, the aquatic niche allowed for lower femoral bone mineral densities than expected based on the body mass alone. The semiaquatic mammals mostly had high bone mineral densities compared to the terrestrial species, which could be considered an adaptation to overcome buoyancy during swimming and shallow diving. Generally, it seems that different osteological properties at the levels of mineral density and biomechanics could be compatible with the adaptation to aquatic, semiaquatic, or terrestrial niches.

Evidence supporting an evolutionary trade-off between material properties and architectural design in Anolis lizard long bones

In biology, "many-to-one mapping" occurs when multiple morphological forms can meet a particular functional demand. Knowledge of this mapping is crucial for understanding how selection on performance shapes the evolution of morphological diversity. Past research has focused primarily on the potential for geometrically alternative morphological designs to produce equivalent performance outcomes. Here, we ask whether the material properties of biological tissues hold similar potential. Through a phylogenetic comparative study of Anolis lizards, we show that the architectural design and mineral density of the femur trade off in a many-to-one functional system, yielding a morphospace featuring parallel isolines in size-relative bending strength. Anole femur evolution has largely tracked a narrow band of strength isolines over phylogenetic timescales, suggesting that geometry and mineral content shape the course of macroevolution through compensatory effects on performance. Despite this conserved evolutionary relationship, insular and continental species evolve strong bones differently, likely reflecting underlying ecological differences. Mainland anoles, which exhibit fast-paced life histories, typically have femora with lower mineralization and thinner walls than island species, which exhibit the opposite strategy. Together, our results reveal an overlooked dimension in the relationship between form and function, expanding our understanding of how many-to-one mapping can shape patterns of phenotypic diversity.

The Role of Locomotory Ancestry on Secondarily Aquatic Transitions

Land-to-sea evolutionary transitions are great transformations where terrestrial amniote clades returned to aquatic environments. These secondarily aquatic amniote clades include charismatic marine mammal and marine reptile groups, as well as countless semi-aquatic forms that modified their terrestrial locomotor anatomy to varying degrees to be suited for swimming via axial and/or appendicular propulsion. The terrestrial ancestors of secondarily aquatic groups would have started off swimming strikingly differently from one another given their evolutionary histories, as inferred by the way modern terrestrial amniotes swim. With such stark locomotor functional differences between reptiles and mammals, we ask if this impacted these transitions. Axial propulsion appears favored by aquatic descendants of terrestrially sprawling quadrupedal reptiles, with exceptions. Appendicular propulsion is more prevalent across the aquatic descendants of ancestrally parasagittal-postured mammals, particularly early transitioning forms. Ancestral terrestrial anatomical differences that precede secondarily aquatic invasions between mammals and reptiles, as well as the distribution of axial and appendicular swimming in secondarily aquatic clades, may indicate that ancestral terrestrial locomotor anatomy played a role, potentially in both constraint and facilitation, in certain aquatic locomotion styles. This perspective of the land-to-sea transition can lead to new avenues of functional, biomechanical, and developmental study of secondarily aquatic transitions.

A macroevolutionary common-garden experiment reveals differentially evolvable bone organization levels in slow arboreal mammals

Eco-morphological convergence, i.e., similar phenotypes evolved in ecologically convergent taxa, naturally reproduces a common-garden experiment since it allows researchers to keep ecological factors constant, studying intrinsic evolutionary drivers. The latter may result in differential evolvability that, among individual anatomical parts, causes mosaic evolution. Reconstructing the evolutionary morphology of the humerus and femur of slow arboreal mammals, we addressed mosaicism at different bone anatomical spatial scales. We compared convergence strength, using it as indicator of evolvability, between bone external shape and inner structure, with the former expected to be less evolvable and less involved in convergent evolution, due to anatomical constraints. We identify several convergent inner structural traits, while external shape only loosely follows this trend, and we find confirmation for our assumption in measures of convergence magnitude. We suggest that future macroevolutionary reconstructions based on bone morphology should include structural traits to better detect ecological effects on vertebrate diversification.

Bridging Performance and Adaptive Landscapes to Understand Long-Term Functional Evolution

AbstractUnderstanding functional adaptation demands an integrative framework that captures the complex interactions between form, function, ecology, and evolutionary processes. In this review, we discuss how to integrate the following two distinct approaches to better understand functional evolution: (1) the adaptive landscape approach (ALA), aimed at finding adaptive peaks for different ecologies, and (2) the performance landscape approach (PLA), aimed at finding performance peaks for different ecologies. We focus on the Ornstein-Uhlenbeck process as the evolutionary model for the ALA and on biomechanical modeling to estimate performance for the PLA. Whereas both the ALA and the PLA have each given insight into functional adaptation, separately they cannot address how much performance contributes to fitness or whether evolutionary constraints have played a role in form-function evolution. We show that merging these approaches leads to a deeper understanding of these issues. By comparing the locations of performance and adaptive peaks, we can infer how much performance contributes to fitness in species' current environments. By testing for the relevance of history on phenotypic variation, we can infer the influence of past selection and constraints on functional adaptation. We apply this merged framework in a case study of turtle shell evolution and explain how to interpret different possible outcomes. Even though such outcomes can be quite complex, they represent the multifaceted relations among function, fitness, and constraints.

Size And Locomotor Ecology Have Differing Effects on the External and Internal Morphologies of Squirrel (Rodentia: Sciuridae) Limb Bones

Mammals exhibit a diverse range of limb morphologies that are associated with different locomotor ecologies and structural mechanics. Much remains to be investigated, however, about the combined effects of locomotor modes and scaling on the external shape and structural properties of limb bones. Here, we used squirrels (Sciuridae) as a model clade to examine the effects of locomotor mode and scaling on the external shape and structure of the two major limb bones, the humerus and femur. We quantified humeral and femoral morphologies using 3D geometric morphometrics and bone structure analyses on a sample of 76 squirrel species across their four major ecotypes. We then used phylogenetic generalized linear models to test how locomotor ecology, size, and their interaction influenced morphological traits. We found that size and locomotor mode exhibit different relationships with the external shape and structure of the limb bones, and that these relationships differ between the humerus and femur. External shapes of the humerus and, to a lesser extent, the femur are best explained by locomotor ecology rather than by size, whereas structures of both bones are best explained by interactions between locomotor ecology and scaling. Interestingly, the statistical relationships between limb morphologies and ecotype were lost when accounting for phylogenetic relationships among species under Brownian motion. That assuming Brownian motion confounded these relationships is not surprising considering squirrel ecotypes are phylogenetically clustered; our results suggest that humeral and femoral variation partitioned early between clades and their ecomorphologies were maintained to the present. Overall, our results show how mechanical constraints, locomotor ecology, and evolutionary history may enact different pressures on the shape and structure of limb bones in mammals.

Fossorial mammals emphasise the forelimb muscle moment arms used for digging: New indices for reconstruction of the digging ability and behaviours in extinct taxa

Among fossorial mammals, forelimbs are major digging apparatuses for dwelling, sheltering and foraging underground. Forelimb-diggers have independently evolved in many lineages of mammals; thus, the method of digging with forelimbs varies by taxon. Therefore, the reconstruction of digging behaviours in extinct animals leads us to understand the evolutionary process of fossorial adaptation in each lineage. However, no morphological index was found to reconstruct if, or how, extinct taxa dug with forelimbs. In this study, we used the shoulder and elbow muscle moment arms in relation to the out-force lever on the manus as indices of the efficiency of motions. The mechanical advantage of two shoulder motions (medial rotation and retraction) and three elbow motions (extension, flexion and adduction) was measured in 381 extant mammal specimens representing 332 species, 279 genera, 103 families and 24 orders. Assuming that both forelimb-digging and -paddling in water require relatively high-output moment arm efficiency, the studied taxa were categorised into four groups based on the presence or absence of forelimb-digging and -paddling abilities. We found that the efficiencies of all five muscle moment arms in the forelimb-diggers and -paddlers were higher than those of the non-diggers and non-paddlers. Furthermore, among the forelimb-diggers, the taxa that dig compact substrates or frequently burrow tend to emphasise the muscle moment arms compared to the taxa that dig loose substrates or dig less frequently. The comparison among the 53 extant forelimb-diggers revealed that the efficiency marked among the five muscle moment arms reflects the difference in digging strategy: humeral rotation diggers emphasise the shoulder medial rotator and elbow adductor, hook-and-pull diggers emphasise the shoulder retractor and elbow flexor and scratch diggers emphasise the shoulder retractor and elbow extensor. We propose that these indices will be powerful tools for reconstructing the fossorial behaviours of extinct mammals. Applying these indices to extinct taxa, Ceratogaulus, Ernanodon, Metacheiromys and Prozaedyus are capable of more efficient forelimb-digging, and each may have adopted different digging strategies.

SegmentGeometry: a tool for measuring second moment of area in 3D slicer

Second moment of area is a measure of how well the cross-section of a beam will resist bending because of its shape. Many have used second moment of area to investigate the mechanical adaptations of biological structures from stingray jaws to animal limb bones. In this context it is important to acknowledge the assumptions of beam theory, in which second moment of area plays a key role, if reasonable results are desired. For example, to minimize shear the structure should be at least ten times longer than it is wide and deflection should be minimal. Analyzing the internal geometry of biological structures has never been easier or more accessible given the wide, and growing availability of micro-CT scans. Here, we offer a guide on the care that needs to be taken when interpreting second moment of area, and present open-access, open-source software that can process hundreds if not thousands of structures in a short time frame. SegmentGeometry, an extension for the open-source imaging platform 3D Slicer, iterates slice-by-slice through 3D structures calculating second moment of area and other cross-sectional properties. We analyzed two case studies to demonstrate the power of this tool and to highlight interpretations that can be gleaned from second moment of area. Second moment of area is just one part of the Euler-Bernoulli beam theory and considering the full equation would greatly increase the number and diversity of questions that can be answered.

Influence of captivity and selection on limb long bone cross-sectional morphology of reindeer

The emergence of pastoralism and animal husbandry has been a critical point in the history of human evolution. Beyond profound behavioural changes in domesticated animals compared to wild ones, characterising the morphological changes associated with domestication process remains challenging. Because reindeer (Rangifer tarandus) can be considered to still be in the early phases of the domestication process, the study of modern populations provides a unique opportunity to examine the impact of captivity and selective breeding on skeletal changes. In this work, we investigated the morphological changes in long limb bone cross-sections using 137 wild and domestic reindeer individuals bred in free-range, in captivity or used for racing and pulling. The shape and shaft cortical thickness of the six long limb bones (i.e., humerus, radioulna, metacarpal, femur, tibia and metatarsal) were measured using a 2D-geometric morphometrics approach taking into account subspecies, sex, body mass and lifestyle differences. These bones are important to understanding functional morphological changes because they can provide information on feeding and locomotor behaviours, as well as on body propulsion and weight bearing. Apart from the effects of taxonomy, etho-ecology and sex, we have found that captivity and selection induced important variations in the size and body mass of modern reindeer. Our results also showed that patterns of variation in cortical bone thickness of long limb bone cross-sections were strongly impacted by body mass and human-imposed restrictions in roaming. This demonstrates that bone cross-sections can provide information on changes in locomotor, reproductive and feeding behaviours induced by the domestication process. These results are valuable not only for (paleo) biologists studying the impact of captivity and selection in ungulates but also for archaeologists exploring the origins of domestication and early herding strategies.

Differing effects of size and lifestyle on bone structure in mammals

BACKGROUND

Mammals are a highly diverse group, with body mass ranging from 2 g to 170 t, and encompassing species with terrestrial, aquatic, aerial, and subterranean lifestyles. The skeleton is involved in most aspects of vertebrate life history, but while previous macroevolutionary analyses have shown that structural, phylogenetic, and functional factors influence the gross morphology of skeletal elements, their inner structure has received comparatively little attention. Here we analysed bone structure of the humerus and mid-lumbar vertebrae across mammals and their correlations with different lifestyles and body size.

RESULTS

We acquired bone structure parameters in appendicular and axial elements (humerus and mid-lumbar vertebra) from 190 species across therian mammals (placentals + marsupials). Our sample captures all transitions to aerial, fully aquatic, and subterranean lifestyles in extant therian clades. We found that mammalian bone structure is highly disparate and we show that the investigated vertebral structure parameters mostly correlate with body size, but not lifestyle, while the opposite is true for humeral parameters. The latter also show a high degree of convergence among the clades that have acquired specialised (non-terrestrial) lifestyles.

CONCLUSIONS

In light of phylogenetic, size, and functional factors, the distribution of each investigated structural parameter reveals patterns explaining the construction of appendicular and axial skeletal elements in mammalian species spanning most of the extant diversity of the clade in terms of body size and lifestyle. These patterns should be further investigated with analyses focused on specific lifestyle transitions that would ideally include key fossils.

Phylogeny, function and ecology in the deep evolutionary history of the mammalian forelimb

Mammals are the only living members of the larger clade Synapsida, which has a fossil record spanning 320 Ma. Despite the fact that much of the ecological diversity of mammals has been considered in the light of limb morphology, the ecological comparability of mammals to their fossil forerunners has not been critically assessed. Because of the wide use of limb morphology in testing ecomorphological hypothesis about extinct tetrapods, we sought: (i) to estimate when in synapsid history, modern mammals become analogues for predicting fossil ecologies; (ii) to document examples of ecomorphological convergence; and (iii) to compare the functional solutions of distinct synapsid radiations. We quantitatively compared the forelimb shapes of the multiple fossil synapsid radiations to a broad sample of extant Mammalia representing a variety of divergent locomotor ecologies. Our results indicate that each synapsid radiation explored different areas of morphospace and arrived at functional solutions that reflected their distinctive ancestral morphologies. This work counters the narrative of non-mammalian synapsid forelimb evolution as a linear progression towards more mammalian morphologies. Instead, a disparate array of early-evolving shapes subsequently contracted towards more mammal-like forms.

Differing limb functions and their potential influence upon the diversification of the mustelid hindlimb skeleton

Though form-function relationships of the mammalian locomotor system have been investigated for over a century, recent models of trait evolution have hitherto been seldom used to identify likely evolutionary processes underlying the locomotor system’s morphological diversity. Using mustelids, an ecologically diverse carnivoran lineage, I investigated whether variation in hindlimb skeletal morphology functionally coincides with climbing, digging, swimming and generalized locomotor habits by using 15 linear traits of the femur, tibia, fibula, calcaneum and metatarsal III across 44 species in a principal component analysis. I subsequently fit different models of Brownian motion and adaptive trait diversification individually to each trait. Climbing, digging and swimming mustelids occupy distinct regions of phenotypic space characterized by differences in bone robustness. Models of adaptive and neutral evolution are, respectively, the best fits for long bone lengths and muscle in-levers, suggesting that different kinds of traits may be associated with different evolutionary processes. However, simulations based upon models of best fit reveal low statistical power to rank the models. Though differences in mustelid hindlimb skeletal morphology appear to coincide with locomotor habits, further study, with sampling expanded beyond the Mustelidae, is necessary to better understand to what degree adaptive evolution shapes morphological diversity of the locomotor system.

Humeral diaphysis structure across mammals

Long bones comprise articular ends (epiphyses) joined by transitional metaphyses and a diaphysis (shaft). The structure of the latter is often viewed as regularly tubular across tetrapods (limbed vertebrates). However, assessments of the bone structure along the whole diaphysis are rare. Here, I assess whole-diaphysis profiles of global compactness (bone fraction) of 164 species of extant and extinct therian mammals (marsupials + placentals) in a phylogenetically informed context. Generally terrestrial, mammals have acquired multiple times the highly specialized aerial, fully aquatic, and subterranean lifestyles, allowing to potentially associate specific traits with these lifestyles. I show that there is a consistent increase in global compactness along the diaphysis in most mammals. This pattern is modified in a limited number of specialized species: all aerial clades (gliders and bats) have rather uniform and low values, while cetaceans' humeral diaphysis is marked by a slightly more compact mid-diaphyseal region. Among subterranean clades, structure alterations are most obvious in fossorial talpids (true moles) and their highly modified humerus. These results call for the investigation of bone structure in whole skeletal elements of key fossils in order to reconstruct the patterns of evolutionary modifications associated with lifestyle transitions.

The comparative method is not macroevolution: across-species evidence for within-species process

It is common for studies that employ the comparative method for the study of adaptation, i.e. documentation of potentially adaptive across-species patterns of trait-environment or trait-trait correlation, to be designated as "macroevolutionary." Authors are justified in using "macroevolution" in this way by appeal to definitions such as "evolution above the species level." I argue that regarding the comparative method as "macroevolutionary" is harmful because it hides in serious ways the true causal content of hypotheses tested with the comparative method. The comparative method is a means of testing hypotheses of adaptation and their alternatives. Adaptation is a population level phenomenon, involving heritable interindividual variation that is associated with fitness differences. For example, given heritable intrapopulational variation, more streamlined individuals in populations of fast-moving aquatic animals have higher locomotory efficiency and thus better survivorship and more resources directed to reproduction than less streamlined ones. Direct evidence consistent with this population-level scenario includes the observation that many unrelated species of fast-moving aquatic animals have similar streamlined shapes, an example of the comparative method. Crucial to note in this example is that although the data are observed across species, the comparative method for studying adaptation tests hypotheses regarding standard population-level natural selection with no content that can be construed as "macro." Even less "macro," individual-level developmental dynamics can limit or bias the range of variants available for selection. Calling any of these studies "macroevolutionary" implies that some additional process is at work, shrouding the need to test adaptation hypotheses and study the range of variants that can be produced in development.

Phylogeny and foraging behaviour shape modular morphological variation in bat humeri

Bats show a remarkable ecological diversity that is reflected both in dietary and foraging guilds (FGs). Cranial ecomorphological adaptations linked to diet have been widely studied in bats, using a variety of anatomical, computational and mathematical approaches. However, foraging-related ecomorphological adaptations and the concordance between cranial and postcranial morphological adaptations remain unexamined in bats and limited to the interpretation of traditional aerodynamic properties of the wing (e.g. wing loading [WL] and aspect ratio [AR]). For this reason, the postcranial ecomorphological diversity in bats and its drivers remain understudied. Using 3D virtual modelling and geometric morphometrics (GMM), we explored the phylogenetic, ecological and biological drivers of humeral morphology in bats, evaluating the presence and magnitude of modularity and integration. To explore decoupled patterns of variation across the bone, we analysed whole-bone shape, diaphyseal and epiphyseal shape. We also tested whether traditional aerodynamic wing traits correlate with humeral shape. By studying 37 species from 20 families (covering all FGs and 85% of dietary guilds), we found similar patterns of variation in whole-bone and diaphyseal shape and unique variation patterns in epiphyseal shape. Phylogeny, diet and FG significantly correlated with shape variation at all levels, whereas size only had a significant effect on epiphyseal morphology. We found a significant phylogenetic signal in all levels of humeral shape. Epiphyseal shape significantly correlated with wing AR. Statistical support for a diaphyseal-epiphyseal modular partition of the humerus suggests a functional partition of shape variability. Our study is the first to show within-structure modular morphological variation in the appendicular skeleton of any living tetrapod. Our results suggest that diaphyseal shape correlates more with phylogeny, whereas epiphyseal shape correlates with diet and FG.

Investigating the impact of captivity and domestication on limb bone cortical morphology: an experimental approach using a wild boar model

The lack of bone morphological markers associated with the human control of wild animals has prevented the documentation of incipient animal domestication in archaeology. Here, we assess whether direct environmental changes (i.e. mobility reduction) could immediately affect ontogenetic changes in long bone structure, providing a skeletal marker of early domestication. We relied on a wild boar experimental model, analysing 24 wild-born specimens raised in captivity from 6 months to 2 years old. The shaft cortical thickness of their humerus was measured using a 3D morphometric mapping approach and compared with 23 free-ranging wild boars and 22 pigs from different breeds, taking into account sex, mass and muscle force differences. In wild boars we found that captivity induced an increase in cortical bone volume and muscle force, and a topographic change of cortical thickness associated with muscular expression along a phenotypic trajectory that differed from the divergence induced by selective breeding. These results provide an experimental proof of concept that changes in locomotor behaviour and selective breeding might be inferred from long bones morphology in the fossil and archaeological record. These trends need to be explored in the archaeological record and further studies are required to explore the developmental changes behind these plastic responses.

New insights into the giant mustelids (Mammalia, Carnivora, Mustelidae) from Langebaanweg fossil site (West Coast Fossil Park, South Africa, early Pliocene)

Giant mustelids are a paraphyletic group of mustelids found in the Neogene of Eurasia, Africa and North America. Most are known largely from dental remains, with their postcranial skeleton mostly unknown. Here, we describe new craniodental and postcranial remains of the large lutrine Sivaonyx hendeyi and the leopard-size gulonine Plesiogulo aff. monspessulanus from the early Pliocene site Langebaanweg, South Africa. The new material of the endemic S. hendeyi, includes upper incisors and premolars, and fragmentary humerus, ulna and a complete astragalus. Its postcrania shares more traits with the living Aonyx capensis than the late Miocene Sivaonyx beyi from Chad. Sivaonyx hendeyi could therefore be tentatively interpreted as a relatively more aquatic taxon than the Chadian species, comparable to A. capensis. The new specimens of Plesiogulo comprise two edentulous maxillae, including one of a juvenile individual with incomplete decidual dentition, and a fragmentary forelimb of an adult individual. The new dental measurements point to this form being amongst the largest specimens of the genus. Both P3-4 differs from the very large species Plesiogulo botori from late Miocene of Kenya and Ethiopia. This confirms the existence of two distinct large species of Plesiogulo in Africa during the Mio/Pliocene, P. botori in the Late Miocene of Eastern Africa (6.1–5.5 Ma) and Plesiogulo aff. monspessulanus at the beginning of the Pliocene in southern Africa (5.2 Ma). Lastly, we report for the first time the presence of both Sivaonyx and Plesiogulo in MPPM and LQSM at Langebaanweg, suggesting that the differences observed from the locality may be produced by sedimentation or sampling biases instead of temporal replacement within the carnivoran guild.

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.

Trabecular bone architecture in the stylopod epiphyses of mustelids (Mammalia, Carnivora)

Mustelidae, a carnivoran clade that includes for instance weasels, badgers, otters and martens, has undergone several evolutionary transitions of lifestyle, resulting in specializations for fossorial, natatorial and scansorial locomotion, in addition to more generalized species. The family is therefore regarded as offering an adequate framework for morpho-functional analyses. However, the architecture of the epiphyseal trabecular bone, which is argued to be particularly responsive to the biomechanical environment, has never been studied. Here, we quantify trabecular bone parameters of the proximal and distal epiphyses of the humerus and femur in 29 species of mustelids and assess the differences of these parameters among groups defined a priori based on the aforementioned locomotor types. The parameters are assessed in a phylogenetic framework, taking into account the potential effect on an individual's body mass. The range of variation described by the acquired parameters is relatively restricted when compared to that of other clades. Generalists, however, are featuring a wider range of variation than the other types. While clear discrimination of locomotor types is difficult, some differences were highlighted by our analysis, such as a greater bone fraction associated with the natatorial taxa, which we discuss in a functional context.