Skip to main content

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



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.


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.


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.


Specializations in mammalian limb morphology are well documented and have been the subject of study for more than a century [21, 33, 66, 75, 91], with these studies primarily focusing on bone lengths, diameters, and the in-levers of muscles [9, 10, 12, 23, 27, 28, 39, 40, 53, 69, 78, 84, 87, 88, 97]. Notably, differing locomotor habits within mammals are associated with distinct limb morphologies. Scansorial/climbing mammals are characterized by limb skeletons with relatively elongate and gracile elements both proximally and distally, including elongate digits, whereas cursorial/running mammals are characterized by limb skeletons with gracile elements and elongate distal limb elements. In contrast, fossorial/digging mammals are characterized by more robust and relatively shorter long bones (for shorter muscle out-levers) and longer muscle in-levers (e.g., olecranon process), some of which exhibit large tuberosities (e.g., the deltoid ridge/tuberosity). Natatorial/swimming mammals are also characterized by robust long bones, longer olecranon processes, and, in species relying upon the hindlimb for swimming, elongate metatarsals and phalanges. Remarkably, locomotor adaptations similar to these traits are also known to occur in early mammaliaform taxa with regards to digging (e.g., [67]), swimming (e.g., [49]), and climbing (e.g., [70]).

In spite of the extensive work on adaptations of the external traits of the limb skeleton, particularly long bones, there has been noticeably less investigation into adaptations of traits associated with the internal geometry of long bones: cross-sectional area (CSA), second moment of area (SMA), and section modulus (MOD) (Fig. 1). CSA represents resistance to axial compression, whereas SMA represents resistance to bending about a specified axis. MOD represents structural strength due to cross-sectional shape, with higher values representing a bone being able to incur greater maximum bending moments relative to maximum mechanical stress. Adaptive evolution should occur among these traits, as differing locomotor habits most likely entail differing mechanical demands acting upon the limb skeleton [13, 39, 40, 78, 87, 88, 101]. For instance, the functioning of fossorial limbs in soil should be associated with higher values of CSA, SMA, and MOD due to the high density of soil (1.83–2.58 g/cm3; [85]). Notably, limbs specialized for digging have longer in-levers for stronger joint extension and often long bones with more robust external dimensions [22, 23, 53, 78, 84]. Similarly, the limbs of natatorial taxa are also associated with a greater robustness and more prominent in-levers [9, 78], and greater values of bone cross-sectional traits could also be beneficial to mitigating larger mechanical loads due to functioning in water, which also has a high density. However, a likely more critical biomechanical factor for natatorial tetrapods is the need to counteract buoyancy in aquatic environments. Among natatorial tetrapods, increased long bone cortical thickness and/or compactness has been concluded to counteract buoyancy during diving and aquatic foraging [5, 41,42,43], and greater cortical thickness and compactness should be reflected by greater CSA values, if not SMA and MOD values as well.

Fig. 1

Cross-sectional traits and example loading regimes. In A, a humerus showing the orientation of anatomical axes about which SMA was measured: cranio-caudal (Cr-Ca) and medio-lateral (M-L). The total bone tissue (dark grey) within the cross-section determines the bone’s cross-sectional area (CSA), whereas the distribution of bone tissue about specified axes determines the bone’s second moment of area (SMA). CSA and SMA withstand specific forms of mechanical loads (B). CSA determines resistance to forces causing axial compression along the length of the bone, whereas SMA determines resistance to bending moments that cause a bone to flex about a given axis. SMAML is associated with bending about the M-L axis (i.e., bending in the Cr-Ca plane), whereas SMACC is associated with bending about the Cr-Ca axis (i.e., bending in the M-L plane). Notably the farther bone tissue is from an axis of bending, the greater the resistance to bending. Thus in A, SMAML is greater than SMACC. Note that the example of bending in B is exaggerated for illustrative purposes, and all cross-sectional traits were measured with BoneJ 1.4.2 [17]

Testing for the potential for adaptations in bone cross-sectional traits in mammals requires a lineage that exhibits a range of locomotor behaviors. Mustelids are a species-rich lineage of carnivoran mammals that include scansorial, fossorial, and natatorial specialists, in addition to a more generalized locomotor habit. Within mustelids, the forelimbs play a role in each of their exhibited locomotor habits, including climbing [24, 27, 38], digging [71, 84], and swimming [9, 30], suggesting that the long bones of the forelimb are ideal for determining if divergent locomotor habits are associated with differences in bone cross-sectional traits. Here we test whether CSA, SMA, and MOD of the humerus, radius, and ulna differ (for a given body size) among fossorial, natatorial, scansorial, and generalized mustelids to understand biomechanical differences among mustelids of differing locomotor habits. More specifically, we predict that fossorial and natatorial mustelids will have the highest values of cross-sectional traits. Regarding fossorial mustelids, this prediction is due to their limbs functioning in soil, which is highly dense (1.83–2.58 g/cm3; [85]) and thus likely subjects a limb engaged in digging to high mechanical loads. Regarding natatorial mustelids, this prediction stems from the moderately high density of water (1.0 g/cm3), which could subject a limb engaged in swimming to higher mechanical loads, and the additional need for natatorial taxa to counteract buoyancy in aquatic environments. We also predict that forelimb cross-sectional traits have evolved adaptively within mustelids under selective regimes relating to limb function, and here fit competing models of trait diversification to test this.

Materials and methods

Cross-sectional traits were measured from x-ray computed tomography (CT) scans of the humerus, radius, and ulna at 5% increments along the length of each bone (Fig. 1) in 28 species of mustelid (Table 1; Fig. 2). To minimize the influence of differences in sample size among mustelid locomotor habits, we sampled seven taxa each from generalized, fossorial, natatorial, and scansorial mustelids. All sampled individuals were adults, as determined by the epiphyses being fully fused to the metaphyses. The only exception to this is the specimen representing Melogale orientalis, the Javan ferret badger (ZMB MAM 8949). Although it was a subadult, we chose to sample this specimen, as postcranial material for Melogale sp. is rare in museum collections.

Table 1 Scanned mustelid species alongside their locomotor habit. N = number of specimens for that species. Scan Location indicates the facility where CT scans were made: 1) Museum für Naturkunde Berlin, Berlin, Germany, 2) Royal Veterinary College, Hertfordshire, UK and 3) University of Chicago, Chicago, USA
Fig. 2

Phylogeny of the sampled Mustelidae. The color of terminal branches and taxon names indicate the locomotor habit for each mustelid species. To obtain this phylogeny, the phylogeny of Law et al. [63] was pruned to our sampled taxa

CT scans were generated at the Museum für Naturkunde Berlin, the University of Chicago, and the Royal Veterinary College (Table 1). µCT scans were made on a Phoenix|x-ray Nanotom (GE Sensing and Inspection Technologies GmbH, Wunstorf, Germany) in Berlin and a Phoenix|x-ray Nanotom and v|tome|x combination in Chicago, respectively. At the Royal Veterinary College, lower resolution medical CT scans were made on a GE LightSpeed 16 scanner (GE Medical Systems, Pollards Wood, UK). X-ray slices were reconstructed using the software datos|x-reconstruction version (GE Sensing and Inspection Technologies GmbH, Phoenix|x-ray) and resulting reconstructed bones were oriented in VG Studio Max 2.0 and 2.1 (Volume Graphics, Heidelberg, Germany).

For each μCT scan, DICOM image stacks of > 1000 slices were generated sequentially along a bone’s proximo-distal axis, whereas for medical quality scans, image stacks of 250–270 slices were generated. For each image stack, the number of the slices containing the proximal-most and distal-most extremes of an individual bone were identified. To obtain slices at 5% increments along the length of the bone, the difference was calculated between the slice numbers containing the proximal and distal extremes of a bone, and this value was divided by 20 to downsample the data to 21 slices. Slices analyzed for cross-sectional traits were segmented in ImageJ 1.51n [90]. For segmentation, the area of the image external to the perimeter of the bone was replaced with black pixels (greyscale value = 0) by using the freehand selection tool in combination with the Fill and Clear Outside commands. After segmentation, the 16-bit images produced by VG Studio were converted into 8-bit prior to data collection with BoneJ 1.4.2 [17]. Cross-sectional properties were then collected using BoneJ’s slice geometry function with an orientation of the cranio-caudal axis set to 270°. The minimum greyscale value associated with bone tissue was manually determined for each segmented slice, whereas the maximum greyscale value associated with bone tissues was 255, as all possible background noise and image artifacts had been previously removed with segmentation.

Using BoneJ, the following traits were measured: CSA, SMA, and MOD. In addition to their biomechanical relevance, CSA represents the total amount of bone tissue in a cross-section, whereas SMA and MOD represent the distribution of the total amount of bone tissue about specified axes. We chose two anatomical axes with regards to measuring SMA and MOD: cranio-caudal (CC) and medio-lateral (ML) (Fig. 1). With bone tissue being identified by maximum and minimum greyscale values, and pixels calibrated to mm within the DICOM image format, BoneJ calculated a bone’s CSA from pixels falling in the appropriate range of greyscale values. SMA is calculated as SMA = Σ(Ad2), where A is an unit of area and d is the distance of that unit from the cross-section’s neutral axis, and MOD is calculated as MOD = SMA/c, where c is the distance from the neutral axis to the farthest unit of bone tissue [7]. BoneJ calculates SMA and MOD using the greyscale values denoting bone tissue and pixel coordinates [17]. Given a difference in size of  two orders of magnitude exists in our sample (e.g., Mustela kathiah [0.21 kg] vs. Enhydra lutra [29.50 kg]; [63]), we rendered the trait values dimensionless as follows to facilitate comparison. We first reduced each trait to a single linear dimension, as CSA, SMA, and MOD have units of mm2, mm4, and mm3, respectively. We took the second root of CSA (i.e., CSA1/2), the fourth root of SMA (i.e., SMA1/4), and the third root of MOD (i.e., MOD1/3) to transform each of these traits into units of mm1. After reducing each trait to a single linear dimension, the trait value was further divided by the bone’s proximodistal (inter-articular) length to render it dimensionless. In addition to these traits, we also calculated a dimensionless ratio characterizing a bone’s relative resistance to bending vs. compression:

$$ R=\frac{SMA}{CSA} $$

This metric uses dimensionless SMA and dimensionless CSA and represents the ratio of the resistance to bending (i.e., SMA) to resistance to axial compression (i.e., CSA). This ratio was calculated separately for SMAML (i.e., RML) and SMACC (i.e., RCC). Comparison of this metric across locomotor habits can reveal whether particular locomotor habits are associated with long bones more predisposed to withstanding bending vs. compression.

Statistical analyses

The study of the evolution of biomechanical traits analytically requires a dual approach: an ahistorical approach focused on how the diversity of observed trait values relates to differences in biomechanical capability among species and a historical approach focused on how phylogeny influences biomechanical trait diversity and the underlying processes governing the evolution of these traits. An ahistorical approach is necessary as the mechanics imposed by a species’ ecological niche can only act upon the trait values possessed by the species and not upon a ‘phylogenetically corrected’ trait value. For instance, a bone’s internal stress due to compression (σcompression) would stem from a compressive force (F) acting on that bone’s CSA ((σcompression = F/CSA). Likewise, a bone’s internal stress due to bending (σbending) would stem from a bending moment (M), the bone’s SMA, and an additional term, y, which denotes the furthest distance of bone tissue from the bone’s neutral axis (σbending = My/SMA). These two calculations of bone stress apply regardless of the specific evolutionary processes underlying trait values. Species must be capable of meeting the biomechanical demands of their ecological niche, so an ahistorical analysis would also be informative for understanding how the biomechanical demands of species’ occupied niches are reflected in morphology.

However, as raised by Felsenstein [29], ahistorical analyses fail to address the influence of the shared ancestry of species upon observed trait values. Likewise, while an ahistorical analysis may contribute to the development of hypotheses regarding specific evolutionary processes (e.g., natural selection), they can in no way directly test for the past action of specific evolutionary processes upon trait diversification. To understand the role of phylogenetic non-independence in trait diversity, as well as to discern likely processes responsible for trait diversification, phylogenetic comparative methods must be employed to put biomechanical traits in a macroevolutionary context.

To address whether the differing biomechanical demands of fossorial, scansorial, natatorial, and generalist locomotor habits may be reflected in cross-sectional morphology, we tested for differences in cross-sectional traits by performing standard one-way ANOVAs using Tukey’s posthoc test (Psignificance ≤ 0.05) for each separate trait, with locomotor habit being the independent factor. Additionally, we tested for differences in cross-sectional traits at 5% increments along the bone’s length, resulting in 19 increments (5–95%). This allowed us to determine if the association between locomotor habit and cross-sectional traits varied along a bone’s length. As we were comparing cross-sectional traits for 19 increments along a bone’s length, we also performed Bonferroni corrections (PBonferroni = 0.05/19 = 0.0026). ANOVAs were performed in R vers. 3.3.1 [81]. To address the role of phylogenetic relatedness, we also performed phylogenetic ANOVAs following the methodology of Adams & Collyer [2] by using the R package geomorph [3]. We additionally assess the clustering of locomotor habit within mustelid phylogeny using two-block partial least squares to test for a correlation between mustelid phylogeny and locomotor habit [2].

To test whether differences in bone cross-sectional traits are likely due to adaptive evolution, we fitted three Ornstein-Uhlenbeck (OU) models of trait evolution to the cross-sectional traits for each 5% increment along a bone’s length: a single optimum (OU1), a three optima (OU3), and a four optima (OU4) model (Fig. 3). OU1 represents a single phenotypic optimum (which can roughly be thought of as a selective pressure) acting across all branches of the phylogeny, OU3 represents separate phenotypic optima for scansorial, natatorial, and remaining mustelids, and OU4 represents separate phenotypic optima for scansorial, natatorial, fossorial, and generalized mustelids. Additionally, for each OU model, we tested a counterpart Brownian motion (BM) model with distinct rates of BM evolution replacing OU phenotypic optima: a single rate (BM1), a three-rate (BM3), and a four-rate model (BM4). BM models were used because differing rates of evolution acting across the phylogeny can result in increased morphological disparity [74]. Lastly, we also fit an Early Burst (EB) model [36], a model in which the net rate of evolution decreases over time (or increases, depending on how the rate parameter is bounded). The purpose of fitting models (OU, BM, and EB) to each 5% increment measured was to determine if the evolutionary dynamics likely governing the evolution of cross-sectional morphology vary along a bone’s length. BM and OU models were fit in R using package OUwie [6], and the EB model was fit using the R package geiger [37]. The recently published phylogeny of Law et al. [63] (Fig. 2) was used for the fitting of trait evolution models, as well as to run our phylogenetic ANOVAs.

Fig. 3

Three hypothetical models of the evolution of cross-sectional morphology in Mustelidae. The phylogeny in A represents a single rate (BM1) or single phenotypic optimum (OU1) process acting across all branches of the phylogeny, both internal and external. The phylogeny in B distinguishes three distinct rates (BM3) or optima (OU3) acting at the terminal branches of the phylogeny, with natatorial, scansorial, and remaining mustelids (i.e., fossorial and generalist taxa) each having their rate/optimum. The phylogeny in C distinguishes four distinct rates (BM4) or optima (OU4) acting at the terminal branches of the phylogeny, one each for the four locomotor habits within Mustelidae. For the BM3/OU3 and BM4/OU4 models, rates and optima acting along internal branches were estimated using stochastic character mapping [8, 45], which reflects uncertainty in character states of internal branches (see main text)

While locomotor habits are known for the terminal taxa in the phylogeny (Table 1), the locomotor habits along the internal branches of the phylogeny are not. To incorporate uncertainty in ancestral states of locomotor habits into our analyses, we performed stochastic character mapping [8, 45] with the R package phytools [82]. In this method, the discrete character in question is randomly mapped onto the internal branches of the phylogeny and the models are then fitted to the tree, with model parameters and criteria being calculated. This procedure of random mapping of the discrete character and model fitting is then repeated over a number of iterations (for the current study: 500), and the model parameters are then averaged over all iterations. Best-fitting models were determined for CSA and SMA at each 5% increment of a bone's length by calculating Akaike’s Information Criterion for small samples (AICc: [47]) and taking the mean of this parameter across all iterations. Using the mean AICc values, Akaike weights [100] were then calculated to determine the relative fit of each model. Akaike weights were calculated with the R package qpcR [94]. As both SMA and MOD both reflect the distribution of total bone tissue about an axis of bending, we chose to model the evolution of only one of these traits.


For the humerus, radius, and ulna, long bone cross-sectional traits significantly differed among mustelid locomotor habits. Natatorial mustelids (e.g., otters) tended to have humeri with the greatest values of CSA, SMA (Fig. 4), and MOD (Additional file 1: Figure S1), followed by fossorial mustelids (e.g., badgers). Scansorial mustelids (e.g., martens) tended to have the lowest values. Notably across nearly the entire length of the humerus, natatorial and fossorial mustelids had significantly greater values of cross-sectional traits than scansorial mustelids. Additionally, natatorial mustelids had significantly greater values of cross-sectional traits than generalized mustelids (e.g., weasels) along largely the entire length of the humerus. Fossorial mustelids also had significantly greater SMACC and MODCC than generalist mustelids across almost the entire length of the humerus; significant differences in CSA, SMAML, and MODML between these two groups tended to be more localized along the length of the humerus. Primarily along the diaphysis, natatorial mustelids were found to have significantly greater CSA, SMAML and MODML than fossorial mustelids. In contrast, significant differences in SMACC and MODCC were largely nonexistent between fossorial and natatorial mustelids. Notably, dimensionless trait values for the subadult Melogale orientalis did not appear to be outliers and clustered closely with other fossorial mustelids, in particular Melogale moschata.

Fig. 4

Differences in humeral, radial, and ulnar CSA and SMA among mustelid locomotor habits. Differences in CSA, SMAML, and SMACC were tested at 5% increments along each bone’s length. CSA and SMA were rendered dimensionless by respectively taking the square and fourth root and dividing by bone length. Colored symbols indicate a significant difference (Adjusted P < 0.05) for the indicated pairwise comparison at a given increment. A lack of significant pairwise differences for a given increment indicates an overall ANOVA result of P > 0.0026 for that increment (the Bonferroni-corrected P-value). For the mechanical significance of these traits, see Fig. 1

Differences in cross-sectional traits of the radius were similar to those for the humerus. The most prevalent pairwise differences among locomotor groups were that natatorial mustelids significantly differed (i.e., greater values for traits) from both scansorial and generalized mustelids (Fig. 4 and Additional file 1: Figure S1). Additionally, natatorial mustelids significantly differed from fossorial mustelids with regards to CSA largely across the length of the radius; however, significant differences between these two locomotor habits were more regionalized to the diaphysis with regards to the remaining traits. Scansorial mustelids also exhibited significantly lower values of radial CSA than generalist taxa across the entire length of the radius.

Across the ulna’s length, significant differences in cross-sectional traits were recovered (Fig. 4 and Additional file 1: Figure S1). As was the case with the humerus and radius, the most ubiquitous pairwise differences consisted of natatorial mustelids having significantly greater values of cross-sectional traits than both scansorial and generalist mustelids. Natatorial mustelids also possessed greater values of cross-sectional traits than fossorial mustelids; however, significant differences between these groups were only pervasive across the ulna’s length with regards to CSA. Scansorial and fossorial mustelids also tended to differ from one another for ulnar traits; however, with the exception of CSA, significant differences between these groups did not extend across the length of the ulna.

Testing for significant differences among the ratio R calculated from dimensionless SMA and CSA revealed that only humeral RCC showed widespread differences among differing mustelid locomotor habits (Fig. 5). These differences consisted largely of fossorial mustelids being distinct from other mustelids towards more proximal (25–45% humeral length) and distal (70–85% humeral length) regions of the humerus and natatorial mustelids contrasting to other mustelids at roughly the mid-length of the humerus (45–65%). For humeral RML, there were no significant differences among groups. Likewise for the ratios of the radius and ulna, there were few, if any, significant differences along their lengths, indicating that the ratio of the resistance to bending to the resistance to compression did not significantly differ among mustelid locomotor habits. R values were largely below 1.0, indicating that mustelid long bones have a greater relative resistance to axial compression than bending.

Fig. 5

Differences in resistance ratio R for the humerus, radius, and ulna among mustelid locomotor habits. Dimensionless values of CSA and SMA were used to calculate R. Colored symbols indicate a significant difference (Adjusted P < 0.05) for the indicated pairwise comparison at a given increment. A lack of significant pairwise differences for a given increment indicates an overall ANOVA result of P > 0.0026 for that increment (the Bonferroni-corrected P-value). Ratios were analysed at 5% increments of bone length

Trait evolution

Phylogenetic ANOVA

Mustelid locomotor habit was significantly correlated with phylogeny (r = 0.912; P = 0.0001). Phylogenetic ANOVAs recovered did not recover any significant differences in humeral and ulnar cross-sectional properties, with the exception of humeral CSA (P = 0.0210) at 85% humeral length and ulnar SMACC (P = 0.0410) and MODML (P = 0.0490) at 20 and 30%, respectively. For the radius, there were significant differences in all cross-sectional traits at 5 and 10% of radial length (P ≤ 0.0250). Mustelid locomotor groups also significantly differed at 25% of radius length for CSA, SMAML, and MODML (P < 0.0300) and at 95% of radius length for CSA, SMACC, and MODCC (P < 0.0500).

Trait evolution

The results of model fitting indicated that forelimb cross-sectional traits most likely evolved under either a multi-rate Brownian motion or a multi-optima OU model with distinct rates or optima for at least natatorial and climbing mustelids (Fig. 6; Additional file 2: Tables S1-S3). With regards to humeral CSA, of the 19 increments for which traits were sampled, the most prevalent model of best fit was a four-optima OU model (OU4), as the best fitting model for 12 increments (Additional file 2: Table S1). The model with next highest number of increments was the three-optima OU model (OU3) with four increments, followed by the single and four-rate Brownian motion models (BM1 and BM4, respectively) with two and one increments, respectively. With regards to humeral SMAML, OU4 and BM4 were the best fitting models for eight and six increments, respectively, followed by a three-rate Brownian, which best fit four increments. For humeral SMACC, the best-fitting model was overwhelmingly OU4. In the case of radial CSA, the best-fitting model was a mix largely of BM3, BM4, and OU3 (four, five, and eight increments, respectively) (Additional file 2: Table S2), whereas for radial SMAML, the best-fitting models was largely OU3. In the case of radial SMACC. the best-fitting model was a mix consisting primarily of BM3, OU3, and BM4 (eight, five, and four increments, respectively). Ulnar CSA and SMACC were best fit by OU3 extensively across the length of the bone (Additional file 2: Table S3). OU3 was the best-fitting model for 11 increments of ulnar SMAML, followed by OU4, which was the best fit for four increments.

Fig. 6

Best-fitting models of trait evolution for humeral, radial, and ulnar cross-sectional traits. For each model, the total number of increments best fit by the model are presented as black bars. In the case of the Ornstein-Uhlenbeck models, OU3 and OU4, the superimposed red bars indicate the number of OU models that also had a significant value of α, representing the strength of selection

In the cases where an OU process was the best-fitting model, inspection of α-values revealed many instances where α could not be distinguished from a value of 0.0 (Fig. 6, red vs. black superimposed bars; Additional file 3: Tables S4-S6). In these instances, the fitted OU model did not substantially differ from a BM model of trait evolution. Only for humeral CSA and SMACC, did a reasonably large fraction of best-fitting OU models possess significant α-values (Additional file 3: Table S4). For radial and ulnar CSA and SMAML, only a couple of α-values proved significant, whereas for radial and ulnar SMAML and ulnar SMACC none of the α-values for best-fitting OU models were significant (Fig. 6; Additional file 3: Tables S5 and S6). Additionally, the confidence limits for α with regards to traits for all three bones – most markedly those of those of the ulna – suffered from implausibly wide confidence limits, indicating high difficulty in fitting OU models to ulnar CSA, SMAML, and SMACC.


Cross-sectional morphology and biomechanical advantage

Our results demonstrate that the cross-sectional morphology of long bones can differ among specialized locomotor habits in Mustelidae, a functionally diverse and speciose lineage within Carnivora. These findings fit well with broader patterns throughout the vertebrate skeleton, indicating linkages of form, function and behavior or performance; in other words, that bone geometry grossly reflects loading patterns. For example, differences in long bone cross-sectional traits have been reported in birds of differing locomotor modes [34, 92] and primates differing in slow climbing, suspensory, and leaping locomotor habits [16, 44, 86]. While many prior studies have focused on external bone dimensions and their relationships with higher-level biological factors such as locomotion, we have contributed a new, focused analysis of Mustelidae using the valuable perspective that analysis of internal (i.e. cross-sectional) bone dimensions can bring. Our findings have some general correspondences with similar analyses by Doube et al. [18,19,20]; cf. shapes of curves in our Figs. 4 and 5) and others, but the Mustelidae-specific insights are important and novel.

In line with our prediction, natatorial and fossorial mustelids tended to have greater values of cross-sectional traits than remaining mustelids (Fig. 4 and Additional file 1: Figure S1), with natatorial and scansorial mustelids possessing respectively the highest and lowest values of cross-sectional traits. Natatorial mustelids significantly differed from scansorial (red circles) and generalized mustelids (bronze triangles) in all cross-sectional traits of the humerus and radius, with these differences being rather extensive along these bones’ lengths. Natatorial mustelids also significantly differed from scansorial and generalized mustelids in ulnar cross-sectional traits, though to a much lesser extent, apart from ulnar CSA. In contrast, fossorial mustelids tended to significantly differ primarily from scansorial mustelids (navy circles) in humeral and radial cross-sectional traits; however, differences in ulnar SMA and MOD were not prevalent across the entire length of the ulna. We found that significant differences between fossorial and generalized mustelids (red triangles) only occurred in humeral SMA and MOD at localized regions along this bone’s length. Thus in mustelids, any possible locomotion-distinct phenotypes associated with cross-sectional morphology do not necessarily encompass all forelimb bones or all cross-sectional traits. Moreover, the four locomotor habits within Mustelidae sampled here likely are not each characterized by a distinct cross-sectional morphology, because generalized mustelids only rarely differed in cross-sectional traits from both scansorial and fossorial species (bronze circles and red triangles, respectively).

The low values of cross-sectional traits of the forelimb skeleton in scansorial mustelids correspond to the greater gracility of their forelimb skeleton [27, 53] and the relatively elongate and lightweight limbs of scansorial mammals in general [51, 52]. The gracile and elongate forelimb skeleton of martens, though not as extreme as in other scansorial carnivorans [64], likely confers advantages in bridging discontinuities in supports (e.g., tree branches) while climbing [13]. Moreover, Cartmill [13] argued that larger body sizes may hamper climbing ability; therefore it is also seems plausible that overly robust or more massive limbs may also be disadvantageous for a climbing lifestyle.

Greater values of cross-sectional traits strongly distinguishing natatorial mustelids from scansorial, generalized, and, to a lesser extent, fossorial mustelids (Fig. 4 and Additional file 1: Figure S1: red circles, bronze and navy triangles) indicate that otters have humeri, radii, and, to a lesser degree, ulnae with greater relative resistance to compression (i.e., CSA) and bending (i.e., SMA) and greater structural strength (i.e., MOD) than mustelids of other locomotor habits. The greater values of cross-sectional traits for natatorial mustelids would be advantageous for swimming by drag-based propulsion, though the degree to which forelimbs function in swimming varies among otter species. Notably, the forelimbs of sea otters (Enhydra lutris) do not play a role in swimming but are extensively involved in tool use and prey manipulation, such as hammering open or prying prey loose [50, 68]. It could be that the forces generated in this behavior could require a forelimb skeleton structurally stronger than other mustelids; however, there currently appears to be no published data on the mechanics of tool use in this species.

However, greater values of cross-sectional traits – and consequently the increased load resistance they offer – are likely not critical to swimming in mustelids. Recent work comparing bone loading in turtles, both during walking and swimming, found significantly lower bone strains during swimming than walking, likely due to buoyant forces removing the need for the limbs to support body weight despite their roles in providing thrust during locomotion [103, 104].

An alternative and more likely, though not mutually exclusive, explanation would be the need for thicker bones to help counteract buoyancy during subsurface swimming [42, 43]. Given that natatorial mustelids have the highest values of cross-sectional traits, it strongly suggests that the need to counteract buoyancy may have a stronger influence upon cross-sectional morphology than any increased resistance to the musculoskeletal loads imposed by specialized limb functions occurring in mustelids. An exception to the general trend among otters is the small-clawed otter (Amblonyx cinereus), which lies comfortably in the range of scansorial mustelids. Notably, this species possesses rather gracile long bones more comparable to scansorial mustelids [27, 53] than to other otters, with its humerus further lacking the strong anterior bowing characteristic of other otters ([9]; pers. obs.). Moreover, this species forages somewhat more terrestrially where it occurs sympatrically with Eurasian otters (Lutra lutra) and smooth-coated otters (Lutrogale perspicillata) [55], and the webbing is incomplete between its digits [60], so it could be considered less aquatic than other otter species.

Fossorial mustelids have high values of cross-sectional traits compared to scansorial and generalized mustelids (Fig. 4 and Additional file 1: Figure S1: brown curves), likely due to the limbs having to function in soil, which has a high density. Although the degree of fossoriality may vary among taxa [84], most badgers and other fossorial mustelids (e.g., zorilla, Ictonyx striatus) dig as a means of foraging and may dig their own burrows [56, 59, 73, 77, 93, 98]. However, some badgers display exceptional digging ability, including rapid digging [56], digging extensive burrow systems [83], digging a new den every night [65], and burying food items several times larger than themselves as a cache [31]. Interestingly, significant differences in SMA and MOD between fossorial and other mustelids were noticeably not as widespread in the ulna as in the humerus and radius. This is surprising given the insertion of the triceps muscle group, which is highly specialized with an angular head in mustelids [25], onto the olecranon process, and the triceps’ highly integral role in exerting force during scratch digging [39, 40, 71]. However, these distinct results for the ulna may be due to the trochlea of the humerus and the trochlear notch of the ulna restricting its movement to flexion and extension relative to the humerus regardless of specializations in limb function. Thus, the ulna cannot exhibit long axis rotation unlike the radius, and thus may experience a lower diversity of loading regimes than the latter bone. Moreover, given that the distal articulation of the radius has much broader contact with the carpus than that of the ulna, it could be possible that the radius receives more of the mechanical loads transmitted proximally along the forelimb by the manus than the ulna, and, if so, this may be reflected in the differences in radial cross-sectional morphology among mustelid locomotor habits. This discrepancy in results among the humerus, radius, and ulna suggests that the loading of limb bones during digging may be complex, with differing bones operating at different loads and safety factors (e.g. perhaps fitting the “mixed chain” hypothesis; [4]).

In addition to function, size may be another factor influencing forelimb morphology in mustelids. In particular, greater values of cross-sectional traits are generally associated with larger body sizes in many mammals [18, 19] and birds [20]. Plotting dimensionless values of mustelid cross-sectional traits against body mass reveals a complicated relationship with body size (Additional file 4: Figure S2). Otters, which include the most massive mustelids, appear to have an allometric trajectory distinct from other mustelids’. However, for a given body mass where multiple locomotor habits coincide, scansorial mustelids have smaller values of cross-sectional traits than either fossorial and generalist mustelids do. This differentiation of locomotor habits for a given body mass suggests that our results are not solely due to the influence of size (i.e. scaling). Rather our results appear subject to the mixed influences of locomotor habit and size.

Resistance to bending vs. compression

The ratio R revealed that, by and large, differences in locomotor habit are not associated with a trade-off in resistance to bending vs. compression (Fig. 5). Humeral RCC was an exception to this, with significant differences occurring among mustelids locomotor habits between 25 and 90% of humeral length. Notably, in contrast to our separate tests of individual cross-sectional traits, RCC distinguished fossorial and natatorial mustelids, with badgers having significantly greater values of RCC than otters. This result suggests that, in the case of badger humeri, possible selection with regards to fossoriality in mustelids may pertain more to the ratio of resistance to particular loading regimes than the absolute resistance to a single loading regime. Compared to otters, badgers exhibit humeri with relatively greater resistance to bending about the cranio-caudal axis (Fig. 1) relative to the total amount of bone tissue comprising their humeral cross-section. In other words, badgers have a wider distribution of bone tissue in their humeral cross-section than otters in spite of having a lower amount of overall bone tissue within their cross-section. This result concurs with our earlier finding of badgers having more robust forelimb long bones (in terms of external dimensions; [53]) and our current finding that badgers have lower values of CSA than otters.

Humeral RCC also distinguished natatorial mustelids from scansorial mustelids, with martens having greater SMACC relative to CSA, and fossorial mustelids from generalized mustelids, further suggesting that the humerus’ relative resistance to different loading regimes may distinguish mustelid locomotor habits. The low values of humeral RCC displayed by natatorial mustelids likely reflect the medio-laterally compressed humeral diaphysis of otters, with such compression being common for aquatic tetrapods [104]. These differences in RCC suggest there may be differences in incurred loading regime as forelimbs conduct different functions in mustelids. While this is an exciting topic of investigation, it unfortunately is beyond the scope of our study.

Apart from humeral RCC, there is striking uniformity among other ratios of R (Fig. 5) in Mustelidae. This uniformity suggests that the relative resistance to different loading regimes is not fundamental to functional specializations of the limb and that a single ‘design’ of relative loading resistance allows for disparate limb functions. Furthermore, the uniformity in R values suggests that distribution of bone tissue (i.e., SMA) relative to the total amount of bone tissue (i.e., CSA) of a cross-section may possibly be phylogenetically conserved, or biomechanically or developmentally constrained, at least for the radius and ulna. A conserved internal morphology of the ulna is particularly surprising when considering mammals more broadly, given the variability of the ulna’s external dimensions in terms of its reduction, relative olecranon length, and robustness in relation to specialized limb functions [87, 88], though admittedly mustelids in and of themselves do not display such wide extremes in ulnar morphology. It remains unclear if our findings would, however, relate to the mesopodium or autopodium (carpus/manus).

Evolution of cross-sectional morphology

Within Mustelidae, locomotor habit is intimately linked with phylogeny. Notably, natatorial species evolved from a single ancestor within our sampled mustelids, as is the case for scansorial species (Fig. 1). Among our sample, there is one instance of convergence in fossorial limb function (Ictonyx striatus), though until recently there was thought to be more convergence in fossoriality in Mustelidae [89]. The preponderance of non-significant results for phylogenetic ANOVAs further underscores that phylogeny is a strong component of the observed morphological variation in Mustelidae. However, our lack of significant findings with phylogenetic ANOVAs goes against the known biomechanical relevance of cross-sectional bone dimensions for many mustelid species, particularly otters (e.g., [42]). While phylogenetic ANOVAs are vital to address the influence of shared ancestry upon trait variation, such analyses by themselves could lead to faulty interpretations of how morphology relates to biomechanical function. In turn, while standard ANOVAs are able to discern morphological differences relevant to biomechanics, they obviously fail to address the role of phylogeny in trait variation. Thus, the pairing of both ahistorical and historical analyses is required for a more comprehensive view of the evolution of biomechanical traits.

Fitting models of trait evolution uncovered that the most likely pattern of evolution with regards to the cross-sectional traits of the humerus, radius, and ulna was either a multi-rate Brownian motion model (BM3/BM4) or a multi-optima Ornstein-Uhlenbeck model (OU3/OU4) (Fig. 6). These models distinguish either distinct rates of evolution (Brownian motion models) or evolution towards distinct phenotypic optima (Ornstein-Uhlenbeck models) for the differing locomotor habits within Mustelidae. Both of these models propose that natatorial and scansorial mustelids morphologically diverged from one another and remaining mustelids, either by evolving under differing rates of Brownian motion or towards distinct adaptive optima. Moreover, finding BM4 and OU4 as the best-fitting model indicates that each locomotor habit within Mustelidae is tied to a divergence in forelimb cross-sectional traits. This result is in line with these two locomotor habits being the extremes of cross-sectional morphology in mustelid long bones.

The prevalence of OU models as the best-fitting models would suggest that the locomotor diversity among mustelids is the result of evolution towards distinct phenotypes ‘optimal’ for the biomechanical demands of a given locomotor habit. However, inspection of α, commonly interpreted as the strength of selection in OU models [35], is crucial prior to accepting an OU model as the most plausible mode of evolution for a given trait [15]. When α does not significantly differ from 0.0, then the OU model is equivalent to a Brownian motion model [11]. Inspection of α-values in instances where OU models were the best-fitting models revealed numerous instances where α could not be distinguished from 0.0 (Fig. 6 and Additional file 3: Tables S4-S6). The outperformance by the OU model vs. the BM model in these instances was due to the additional parameters of the OU model affording the best description of the data’s variance outside of the model’s biological relevance [15].

It thus appears that the cross-sectional morphology of the mustelid humerus, radius, and ulna has evolved predominantly due to a multi-rate Brownian motion process. Under such a mode of evolution, the distinct cross-sectional morphologies of mustelid locomotor habits are associated with a distinct rate of phenotypic evolution, and it is possible that these differences in rate are associated with different constraints upon the evolution of cross-sectional morphology in mustelids (see [74, 96]). Such constraints regarding mustelid limbs could be the biomechanical benefits of thinner, and presumably more lightweight, bones associated with a climbing lifestyle or more robust bones associated with an aquatic lifestyle (see above). Alternatively, cross-sectional morphology may have been under selection at one point during mustelid evolution, with resulting changes in morphology being conserved among later divergences of mustelids (i.e., phylogenetic inertia in the trait). This would be in contrast to a continuous selective regime acting across the branches associated with those later divergences (as is the case in an OU model).

This overall result contrasts with the likely mode of evolution for the external dimensions of the forelimb skeleton (e.g., lengths, diameters, and muscle in-lever lengths). The external dimensions of the forelimb skeleton likely evolved adaptively, with adaptive peaks distinguishing scansorial from remaining mustelids in terms of the length of muscle in-levers and long bone gracility [53]. Then again, the contrasting results of the current study and those of Kilbourne [53] might be due to sample size. In the current study, we restricted our sample size to seven taxa per locomotor habit for a total of 28 taxa, whereas Kilbourne [53] sampled as many mustelid species as possible for a total of 41 taxa. However, another possible explanation may be that the different, functionally relevant traits within a single functioning organ may evolve by different processes in mustelids. These results also raise the question of how do differing traits, with different biomechanical functions (e.g., the mechanical advantage offered by muscle in-levers vs. the bending resistance offered by SMA), contribute to the overall adaptations occurring in a limb? This question merits future focus in trait evolution studies combining different kinds of traits, though current methods may be ill equipped to address it [1].


The cross-sectional morphologies of the humerus, radius, and ulna exhibit differences among mustelid locomotor specializations. In particular, scansorial and natatorial mustelids are strongly associated with more gracile and more robust cross-sections, respectively. Comparing the ratio of bending to compression resistance in mustelids finds a largely uniform ratio across locomotor habits, excluding humeral RCC, suggesting a possible constraint upon forelimb morphology in mustelids (e.g. conserved distribution of bone tissue relative to the anatomical axes for a given CSA). However, the phylogenetic relationships among mustelids have had considerable influence upon the disparity of mustelid cross-sectional traits. Cross-sectional traits appear to most likely have evolved according to a multi-rate Brownian motion process, with distinct rates of Brownian motion evolution for scansorial and natatorial mustelids. This is in contrast to previous findings when fitting trait diversification models to the external traits (e.g., lengths, diameters) of the forelimb skeleton, which appear to have diversified adaptively. This difference in results may be due to differences in sample size or the possibility that biomechanical traits with differing roles undergo substantially different modes of evolution within a given organ.


  1. 1.

    Adams DC, Collyer ML. Multivariate phylogenetic comparative methods: evaluations, comparisons, and recommendations. Syst Biol. 2018a;67:14–31.

    PubMed  Article  PubMed Central  Google Scholar 

  2. 2.

    Adams DC, Collyer ML. Phylogenetic ANOVA: group-clade aggregation, biological challenges, and a refined permutation procedure. Evolution. 2018b;72:1204–15.

    PubMed  Article  PubMed Central  Google Scholar 

  3. 3.

    Adams, D.C., Collyer, M.L., Kaliontzopoulou, A.. Geomorph: software for geometric morphometric analyses. R package version 3.0.6. 2018.

    Google Scholar 

  4. 4.

    Alexander RM. A theory of mixed chains applied to safety factors in biological systems. J Theor Biol. 1997;184:247–52.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  5. 5.

    Amson E, de Muizon C, Laurin M, Argot C, de Buffrénil V. Gradual adaptation of bone structure to aquatic lifestyle in extinct sloths from Peru. Proc R Soc B. 2014;281:20140192.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Beaulieu JM, Jhwueng D-C, Boettiger C, O’Meara BC. Modeling stabilizing selection: expanding the Ornstein-Uhlenbeck model of adaptive evolution. Evolution. 2012;66:2369–83.

    PubMed  PubMed Central  Article  Google Scholar 

  7. 7.

    Beer FP, Johnston RE, DeWolf JT. Mechanics of materials. 3rd ed. USA: McGraw-Hill; 2001.

    Google Scholar 

  8. 8.

    Bollback JP. SIMMAP: stochastic character mapping of discrete traits on phylogenies. BMC Bioinform. 2006;7:88.

    Article  CAS  Google Scholar 

  9. 9.

    Botton-Divet L, Cornette R, Fabre A-C, Herrel A, Houssaye A. Morphological analysis of long bones in semi-aquatic mustelids and terrestrial relatives. Integr Comp Biol. 2016;56:1298–309.

    PubMed  Article  PubMed Central  Google Scholar 

  10. 10.

    Botton-Divet L, Cornette R, Houssaye A, Fabre A-C, Herrel A. Swimming and running: a study of the convergence in long bone morphology among semi-aquatic mustelids (Carnivora: Mustelidae). Biol J Linn Soc. 2017;121:38–49.

    Article  Google Scholar 

  11. 11.

    Butler MA, King AA. Phylogenetic comparative analysis: a modeling approach for adaptive evolution. Am Nat. 2004;164:683–95.

    PubMed  Article  PubMed Central  Google Scholar 

  12. 12.

    Carrano MT. What, if anything, is a cursor? Categories versus continua for determining locomotor habit in mammals and dinosaurs. J Zool. 1999;247:29–42.

    Article  Google Scholar 

  13. 13.

    Cartmill M. Climbing. In: Hildebrand M, Bramble DM, Liem KF, Wake DB, editors. Functional vertebrate morphology. Cambridge: Belknap Press; 1985. p. 73–88.

    Google Scholar 

  14. 14.

    Clark TW, Anderson E, Douglas C, Strickland M. Martes Americana. Mamm Species. 1987;289:1–8.

    Article  Google Scholar 

  15. 15.

    Cooper N, Thomas GH, Venditti C, Meade A, Freckleton RP. A cautionary note on the use of Ornstein-Uhlenbeck models in macroevolutionary analyses. Biol J Linn Soc. 2016;118:64–77.

    Article  Google Scholar 

  16. 16.

    Demes B, Jungers WL, Selpien K. Body size, locomotion, and long bone-cross-sectional geometry in indriid primates. Amer J Phys Anthro. 1991;86:537–47.

    CAS  Article  Google Scholar 

  17. 17.

    Doube M, Klosowski MM, Arganda-Carreras I, Cordelières FP, Doughtery RP, Jackson JS, Schmid B, Hutchinson JR, Shefelbine SJ. BoneJ: free and extensible bone image analysis in ImageJ. Bone. 2010;47:1076–9.

    PubMed  PubMed Central  Article  Google Scholar 

  18. 18.

    Doube M, Felder AF, Chua MY, Lodhia K, Kłosowski M, Hutchinson JR, Shefelbine SJ. Limb bone scaling in hopping macropods and quadrupedal artiodactyls. R Soc Open Sci. 2018;5:180152.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Doube M, Wiktorowicz Conroy A, Christiansen P, Hutchinson JR, Shefelbine S. Three-dimensional geometric analysis of felid limb bone allometry. PLoS One. 2009;4:e4742.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Doube M, Yen SCW, Klosowski MM, Farke AA, Hutchinson JR, Shefelbine SJ. Whole bone scaling of the avian pelvic limb. J Anat. 2012;221:21–9.

    PubMed  PubMed Central  Article  Google Scholar 

  21. 21.

    Dublin LI. Adaptations to aquatic, arboreal, fossorial and cursorial habits in mammals. II Arboreal adaptations. Am Nat. 1903;37:731–6.

    Article  Google Scholar 

  22. 22.

    Elissamburu A, de Santis L. Forelimb proportions and fossorial adaptations in the scratch-digging rodent Ctenomys (Caviomorpha). J Mammal. 2011;92:683–9.

    Article  Google Scholar 

  23. 23.

    Elissamburu A, Vizcaíno SF. Limb proportions and adaptations in caviomorph rodents (Rodentia: Caviomorpha). J Zool. 2004;262:145–59.

    Article  Google Scholar 

  24. 24.

    Ercoli MD, Youlatos D. Integrating locomotion, postures and morphology: the case of the tayra, Eira Barbara (Carnivora, Mustelidae). Mamm Biol. 2016;81:464–76.

    Article  Google Scholar 

  25. 25.

    Ercoli MD, Álvarez A, Stefanini MI, Busker F, Morales MM. Muscular anatomy of the forelimbs of the lesser grison (Galictis cuja), and a functional and phylogenetic overview of Mustelidae and other Caniformia. J Mammal Evol. 2015;22:57–91.

  26. 26.

    Estes JA. Enhydra lutris. Mamm Species. 1980;133:1–8.

    Article  Google Scholar 

  27. 27.

    Fabre A-C, Cornette R, Slater G, Argot C, Peigné S, Goswami A, Pouydebat E. Getting a grip on the evolution of grasping in musteloidcarnivorans: a three-dimensional analysis of forelimb shape. J Evol Biol. 2013;26:1521–35.

    PubMed  Article  PubMed Central  Google Scholar 

  28. 28.

    Fabre A.-C, Cornette R, Goswami A, Peigne S. Do constraints associated with the locomotor habitat drive the evolution of forelimb shape? A case study in musteloid carnivorans. J. Anat. 2015;226:596–610.

    PubMed  PubMed Central  Article  Google Scholar 

  29. 29.

    Felsenstein J. Phylogenies and the Comparative Method. Amer Nat. 1985;125:1–15.

    Article  Google Scholar 

  30. 30.

    Fish FE. Association of propulsive swimming mode with behavior in river otters (Lutra canadensis). J Mammal. 1994;75:989–97.

    Article  Google Scholar 

  31. 31.

    Frehner EH, Buechley ER, Christensen T, Şekercioğlu ÇH. Subterranean caching of domestic cow (Bos taurus) carcasses by American badgers (Taxidea taxus) in the Great Basin Desert, Utah. W N Amer Nat. 2017;77:124–9.

    Article  Google Scholar 

  32. 32.

    Gorsuch WA, Larivière S. Vormela peregusna. Mamm Species. 2005;779:1–5.

    Article  Google Scholar 

  33. 33.

    Gregory WK. Notes on the principles of quadrupedal locomotion and on the mechanism in the limbs of hoofed animals. Ann N Y Acad Sci. 1912;22:287–94.

    Article  Google Scholar 

  34. 34.

    Habib MB, Ruff CB. The effects of locomotion on the structural characteristics of avian limb bones. Zool. J Linn Soc. 2008;153:601–624.

    Article  Google Scholar 

  35. 35.

    Hansen TF. Stabilizing selection and the comparative analysis of adaptation. Evolution. 1997;51:1341–135.

    PubMed  PubMed Central  Article  Google Scholar 

  36. 36.

    Harmon LJ, Losos JB, Davies TJ, Gillespie RG, Gittleman JL, Jennings WB, et al. Early bursts of body size and shape evolution are rare in comparative data. Evolution. 2010;64:2385–96.

    PubMed  PubMed Central  Google Scholar 

  37. 37.

    Harmon LJ, Weir JT, Brock CD, Glor RE, Challenger W. GEIGER: investigating evolutionary radiations. Bioinformatics. 2008;24:129–31.

    CAS  Article  Google Scholar 

  38. 38.

    Heinrich RE, Biknevicius AR. Skeletal allometry and interlimb scaling patterns in mustelid carnivorans. J Morphol. 1998;235:121–34.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  39. 39.

    Hildebrand M. Walking and running. In: Hildebrand M, Bramble DM, Liem KF, Wake DB, editors. Functional vertebrate morphology. Cambridge: Belknap Press; 1985a. p. 38–57.

    Chapter  Google Scholar 

  40. 40.

    Hildebrand M. Digging in quadrupeds. In: Hildebrand M, Bramble DM, Liem KF, Wake DB, editors. Functional vertebrate morphology. Cambridge: Belknap Press; 1985b. p. 89–109.

    Chapter  Google Scholar 

  41. 41.

    Houssaye A. “Pachyostosis” in aquatic amniotes: a review. Integr Zool. 2009;4:325–40.

    PubMed  Article  PubMed Central  Google Scholar 

  42. 42.

    Houssave A, Botton-Divet L. From land to water: evolutionary changes in long bone microanatomy of otters (Mammalia: Mustelidae). Biol J Linn Soc. 2018;125:240–9.

    Article  Google Scholar 

  43. 43.

    Houssaye A, Sander PM, Klein N. Adaptive patterns in aquatic amniote bone microanatomy – more complex than previously thought. Integr Comp Biol. 2016;56:1349–69.

    PubMed  Article  PubMed Central  Google Scholar 

  44. 44.

    Houssaye A, Taverne M, Cornette R. 3D quantitative comparative analysis of long bone diaphysis variations in microanatomy and crosssectional geometry. J Anat. 2018;232:836–849.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  45. 45.

    Huelsenbeck JP, Neilsen R, Bollback JP. Stochastic mapping of morphological characters. Syst Biol. 2003;52:131–58.

    PubMed  Article  PubMed Central  Google Scholar 

  46. 46.

    Hung N, Law CJ. Lutra lutra (Carnivora: Mustelidae). Mamm Species. 2016;48:109–22.

    Article  Google Scholar 

  47. 47.

    Hurvich CM, Tsai CL. Regression and time series model selection in small samples. Biometrika. 1989;76:297–307.

    Article  Google Scholar 

  48. 48.

    Hwang YT, Larivière S. Lutrogale perspicillata. Mamm Species. 2005;786:1–4.

    Article  Google Scholar 

  49. 49.

    Ji Q, Luo ZX, Yuan C-X, Tabrum AR. A swimming mammaliform from the middle Jurassic and ecomorphological diversification of early mammals. Science. 2006;311:1123–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  50. 50.

    Kenyon KW. The sea otter in the eastern Pacific Ocean. N Amer Fauna. 1969;68:1–352.

    Article  Google Scholar 

  51. 51.

    Kilbourne BM, Hoffman LC. Scale effects between body size and limb design in quadrupedal mammals. PLoS One. 2013;8:e78392.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  52. 52.

    Kilbourne BM, Hoffman LC. Energetic benefits and adaptations in mammalian limbs: scale effects and selective pressures. Evolution. 2015;69:1546–59.

    PubMed  Article  PubMed Central  Google Scholar 

  53. 53.

    Kilbourne BM. Selective regimes and functional anatomy in the mustelid forelimb: diversification towards specializations for climbing, digging, and swimming. Ecol Evol. 2017;7:8852–63.

    PubMed  PubMed Central  Article  Google Scholar 

  54. 54.

    King CM. Mustela erminea. Mamm Species. 1983;195:1–8.

  55. 55.

    Kruuk H. Otters: ecology, behavior, and conservation. New York: Oxford University Press; 2006.

    Book  Google Scholar 

  56. 56.

    Kruuk H, Mills MGL. Notes on food and foraging of the honey badger Mellivora capensis in the Kalahari Gemsbok National Park. Koedoe. 1983;26:153–7.

    Article  Google Scholar 

  57. 57.

    Larivière S. Lontra felina. Mamm Species. 1998;575:1–5.

    Google Scholar 

  58. 58.

    Larivière S. Lontra longicaudis. Mamm Species. 1999;609:1–5.

    Google Scholar 

  59. 59.

    Larivière S. Ictonyx striatus. Mamm Species. 2002;698:1–5.

    Article  Google Scholar 

  60. 60.

    Larivière S. Amblonyx cinereus. Mamm Species. 2003;720:1–5.

    Article  Google Scholar 

  61. 61.

    Larivière S, Jennings AP. Family Mustelidae. In: Wilson DE, MIttermeier RA, editors. Handbook of the Mammals of the World. Vol. 1. Carnivores. Barcelona: Lynx Edicions; 2009. p. 564–658.

    Google Scholar 

  62. 62.

    Law CJ. Mustela sibirica (Carnivora: Mustelidae). Mamm Species. 2018;50:109–18.

    Article  Google Scholar 

  63. 63.

    Law CJ, Slater GJ, Mehta RS. Lineage diversity and size disparity in Musteloidea: testing patterns of adaptive radiation using molecular and fossil-based methods. Syst Biol. 2018;67:127–44.

    PubMed  Article  PubMed Central  Google Scholar 

  64. 64.

    Leach D. The descriptive and comparative postcranial osteology of marten (Martes americana, Turton) and fisher (Martes pennanti, Erxleben): the appendicular skeleton. Can J Zool. 1977;55:199–214.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  65. 65.

    Long CA. Taxidea taxus. Mamm Species. 1973;26:1–4.

    Article  Google Scholar 

  66. 66.

    Lull RS. Adaptations to aquatic, arboreal, fossorial and cursorial habits in mammals. IV. Cursorial adaptations. Am Nat. 1904;38:1–11.

    Article  Google Scholar 

  67. 67.

    Luo ZX, Wible JR. A late Jurassic digging mammal and early mammalian diversification. Science. 2005;308:103–7.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  68. 68.

    Mann J, Patterson EM. Tool use by aquatic animals. Philos Trans R Soc B. 2013;368:20120424.

    Article  Google Scholar 

  69. 69.

    Maynard Smith J, Savage RJG. Some locomotory adaptations in mammals. Zool J Linnean Soc. 1956;42:603–22.

    Article  Google Scholar 

  70. 70.

    Meng Q-J, Ji Q, Zhang Y-G, Liu D, Grossnickle DM, Luo Z-X. An arboreal docodont from the Jurassic and mammaliaform ecological diversification. Science. 2015;347:764–748.

    CAS  PubMed  Article  PubMed Central  Google Scholar 

  71. 71.

    Moore AL, Budny JE, Russell AP, Butcher MT. Architectural specialization of the intrinsic thoracic limb musculature of the American badger (Taxidea taxus). J Morphol. 2013;274:35–48.

    PubMed  Article  PubMed Central  Google Scholar 

  72. 72.

    Noonan P, Prout S, Hayssen V. Pteronura brasiliensis (Carnivora: Mustelidae). Mamm Species. 2017;49:97–108.

    Article  Google Scholar 

  73. 73.

    Nowak RM. Walker’s carnivores of the world. Baltimore: Johns Hopkins University Press; 2005.

    Google Scholar 

  74. 74.

    O’Meara BC, Ané C, Sanderson MJ, Wainwright PC. Testing for different rates of continuous trait evolution using likelihood. Evolution. 2006;60:922–33.

    PubMed  Article  PubMed Central  Google Scholar 

  75. 75.

    Osburn RC. Adaptation to aquatic, arboreal, fossorial and cursorial habits in mammals. I. Aquatic adaptations. Am Nat. 1903;37:651–65.

    Article  Google Scholar 

  76. 76.

    Pasitschniak-Arts M, Larivière S. Gulo gulo. Mamm Species. 1995;499:1–10.

    Google Scholar 

  77. 77.

    Pigozzi G. Digging behaviour while foraging by the European badger, Meles meles, in a Mediterranean habitat. Ethology. 1989;83:121–8.

    Article  Google Scholar 

  78. 78.

    Polly PD. Limbs in mammalian evolution. In: Hall BK, editor. Fins into Limbs: Evolution, Development, and Transformations. Chicago: University of Chicago Press; 2007. p. 245–68.

    Google Scholar 

  79. 79.

    Powell RA. Martes pennanti. Mamm Species. 1981;156:1–6.

    Google Scholar 

  80. 80.

    Presley SJ. Eira Barbara. Mamm Species. 2000;636:1–6.

    Article  Google Scholar 

  81. 81.

    R Development Core Team. R: a language and environment for statistical computing (Vers. 3.3.1). Vienna: R Foundation for Statistical Computing; 2016. ISBN 3-900051-07-0. Available at

    Google Scholar 

  82. 82.

    Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Methods Ecol Evol. 2012;3:217–23.

    Article  Google Scholar 

  83. 83.

    Roper TJ. Badger Meles meles setts – architecture, internal environment and function. Mammal Rev. 1992;22:433–53.

    Article  Google Scholar 

  84. 84.

    Rose J, Moore A, Russell A, Butcher M. Functional osteology of the forelimb digging apparatus of badgers. J Mammal. 2014;95:543–58.

    Article  Google Scholar 

  85. 85.

    Rowell DL. Soil science: Methods & Applications. Abingdon: Routledge Press; 2014.

    Book  Google Scholar 

  86. 86.

    Ruff CB. Long bone articular and diaphyseal structure in Old World monkeys and apes. I: locomotor effects. Amer J Phys Anthro. 2002;119:305–42.

    Article  Google Scholar 

  87. 87.

    Samuels JX, Meachen JA, Sakai SA. Postcranial morphology and the locomotor habits of living and extinct carnivorans. J Morphol. 2013;274:121–46.

    PubMed  Article  PubMed Central  Google Scholar 

  88. 88.

    Samuels JX, Van Valkenburgh V. Skeletal indicators of locomotor adaptations in living and extinct rodents. J Morphol. 2008;269:1387–411.

    PubMed  Article  PubMed Central  Google Scholar 

  89. 89.

    Sato JJ, Wolsan M, Prevosti FJ, D’Elia G, Begg C, Begg K, et al. Evolutionary and biogeographic history of weasel-like carnivorans (Musteloidea). Mol Phylogenet Evol. 2012;63:745–57.

    PubMed  Article  PubMed Central  Google Scholar 

  90. 90.

    Schneider CA, Rasband WS, Eliceiri KW. NIH image to ImageJ: 25 years of image analysis. Nat Methods. 2012;9:671–5.

    CAS  PubMed  PubMed Central  Article  Google Scholar 

  91. 91.

    Shimer HW. Adaptations to aquatic, arboreal, fossorial and cursorial habits in mammals. III. Fossorial adaptations. Am Nat. 1903;37:819–25.

    Article  Google Scholar 

  92. 92.

    Simons EL, Hieronymous TL, O’Connor PM. Cross sectional geometry of the forelimb skeleton and flight mode in Pelecaniform birds. J Morphol. 2011;272:958–71.

    PubMed  Article  PubMed Central  Google Scholar 

  93. 93.

    Snead E, Hendrickson GO. Food habits of the badger in Iowa. J Mammal. 1942;23:380–91.

    Article  Google Scholar 

  94. 94.

    Speiss, A.-J.. qpcR: modelling and analysis of real-time PCR data. R package vers. 1.4-0. 2014.

    Google Scholar 

  95. 95.

    Storz JF, Wozencraft WC. Melogale moschata. Mamm Species. 1999;631:1–4.

    Google Scholar 

  96. 96.

    Thomas GH, Freckleton RP, Székely T. Comparative analyses of the influence of developmental mode on phenotypic diversification rates in shorebirds. Proc R Soc B. 2006;273:1619–24.

    PubMed  Article  PubMed Central  Google Scholar 

  97. 97.

    Van Valkenburgh B. Locomotor behaviors in living and extinct carnivores. J Vertebr Paleontol. 1987;7:162–82.

    Article  Google Scholar 

  98. 98.

    Vanderhaar JM, Ten Hwang Y. Mellivora capensis. Mamm Species. 2003;721:1–8.

    Article  Google Scholar 

  99. 99.

    Warton DI, Duursma Remko A, Falster DS, Taskinen S. smatr 3 - an R package for estimation and inference about allometric lines. Met. Ecol. Evol. 2012;3:257–259.

    Google Scholar 

  100. 100.

    Wagenmakers E-J, Farrell S. AIC model selection using Akaike weights. Psychon Bull Rev. 2004;11:192–6.

    PubMed  Article  PubMed Central  Google Scholar 

  101. 101.

    Webb PW Blake RW. Swimming. In Hildebrand M, Bramble DM, Liem KF & Wake DB (Eds.), Functional Vertebrate Morphology. Cambridge: Belknap Press. 1985;110-128.

  102. 102.

    Yensen E, Tarifa T. Galictis vittata. Mamm Species. 2003;727:1–8.

    Article  Google Scholar 

  103. 103.

    Young VKH, Blob RW. Limb bone loading in swimming turtles: changes in loading facilitate transitions from tubular to flipper-shaped limbs during aquatic invasions. Biol Lett. 2015;11:20150110.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  104. 104.

    Young VKH, Wienands CE, Wilburn P, Blob RW. Humeral loads during swimming and walking in turtles: implications for morphological change during aquatic reinvasions. J Exp Biol. 2017;220:3873–7.

    PubMed  Article  PubMed Central  Google Scholar 

Download references


We thank the following for specimen access: Frieder Mayer, Christiane Funk, Steffen Bock, and Detlef Willborn (Museum für Naturkunde Berlin); Larry Heaney, Bruce Patterson, Adam Ferguson, Rebecca Banasiak, and John Phelps (Field Museum of Natural History); Eline Lorenzen and Daniel Klingberg Johansson (Natural History Museum of Denmark); and Petr Benda (Národní Muzeum). For help with CT scanning and, in some cases, conducting scanning itself, we additionally thank Kristin Mahlow and Johannes Müller (Museum für Naturkunde Berlin) and Zhe-Xi Luo and April Neander (University of Chicago). We also thank Michael Doube for advice on working with BoneJ, Falko Glöckler for computational support, and Will Heflin for help with figures. This study was funded by the Deutsche Forschungsgemeinschaft (grant number Ki 1843/3-1), and the publication of this article was funded by the Open Access Fund of the Leibniz Association.


This study was funded in its entirety by the Deutsche Forschungsgemeinschaft (grant number Ki 1843/3-1).

Availability of data and materials

The data and code generated for this study are available at the online repositories Dryad ( and Morphosource (

Author information




BMK selected specimens, generated CT scans, collected and analyzed data, and wrote the manuscript. JRH generated CT scans, provided training in methodology, and wrote the manuscript. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Brandon M. Kilbourne.

Ethics declarations

Ethics approval and consent to participate

Not applicable

Consent for publication

Not applicable

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Figure S1. Differences in humeral radial, and ulnar section modulus among mustelid locomotor habits. Differences in MODML and MODCC were tested at 5% increments along each bone’s length. MOD was rendered dimensionless by respectively taking the third root and dividing by bone length. Colored symbols indicate a significant difference (Adjusted P < 0.05) for the indicated pairwise comparison at a given increment. A lack of significant pairwise differences for a given increment indicates an overall ANOVA result of P > 0.0026 for that increment (the Bonferroni-corrected P-value). (EPS 7340 kb)

Additional file 2:

Table S1. Akaike weights for trait diversification models fitted to humeral cross-sectional traits. Table S2. Akaike weights for trait diversification models fitted to radial cross-sectional traits. Table S3. Akaike weights for trait diversification models fitted to ulnar cross-sectional traits. (DOCX 58 kb)

Additional file 3:

Table S4. Significance and magnitude of alpha values for OU models that were determined to be the best fitting model for humeral cross-sectional traits. Table S5. Significance and magnitude of alpha values for OU models that were determined to be the best fitting model for radial cross-sectional traits. Table S6. Significance and magnitude of alpha values for OU models that were determined to be the best fitting model for ulnar cross-sectional traits. (DOCX 51 kb)

Additional file 4:

Figure S2. Allometric patterns present within mustelid humeral cross-sectional traits, here represented by dimensionless CSA and SMA plotted against body mass on a logarithmic scale. The rows of blue, red and grey numbers on each plot show the scaling exponents (and 95% confidence intervals) for each regression (for behavioral sub-groups or all of Mustelidae). All regressions shown have exponents > 0, and hence slight positive allometry is evident for all parameters quantified. Regressions were performed with the R package smatr [99]. (EPS 3570 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Kilbourne, B.M., Hutchinson, J.R. Morphological diversification of biomechanical traits: mustelid locomotor specializations and the macroevolution of long bone cross-sectional morphology. BMC Evol Biol 19, 37 (2019).

Download citation


  • Mustelidae
  • Locomotion
  • Trait evolution
  • Adaptation
  • Morphological diversification
  • Cross-sectional properties