Skip to main content

Trabecular architecture in the sciuromorph femoral head: allometry and functional adaptation

Abstract

Background

Sciuromorpha (squirrels and close relatives) are diverse in terms of body size and locomotor behavior. Individual species are specialized to perform climbing, gliding or digging behavior, the latter being the result of multiple independent evolutionary acquisitions. Each lifestyle involves characteristic loading patterns acting on the bones of sciuromorphs. Trabecular bone, as part of the bone inner structure, adapts to such loading patterns. This network of thin bony struts is subject to bone modeling, and therefore reflects habitual loading throughout lifetime. The present study investigates the effect of body size and lifestyle on trabecular structure in Sciuromorpha.

Methods

Based upon high-resolution computed tomography scans, the femoral head 3D inner microstructure of 69 sciuromorph species was analyzed. Species were assigned to one of the following lifestyle categories: arboreal, aerial, fossorial and semifossorial. A cubic volume of interest was selected in the center of each femoral head and analyzed by extraction of various parameters that characterize trabecular architecture (degree of anisotropy, bone volume fraction, connectivity density, trabecular thickness, trabecular separation, bone surface density and main trabecular orientation). Our analysis included evaluation of the allometric signals and lifestyle-related adaptation in the trabecular parameters.

Results

We show that bone surface density, bone volume fraction, and connectivity density are subject to positive allometry, and degree of anisotropy, trabecular thickness, and trabecular separation to negative allometry. The parameters connectivity density, bone surface density, trabecular thickness, and trabecular separation show functional signals which are related to locomotor behavior. Aerial species are distinguished from fossorial ones by a higher trabecular thickness, lower connectivity density and lower bone surface density. Arboreal species are distinguished from semifossorial ones by a higher trabecular separation.

Conclusion

This study on sciuromorph trabeculae supplements the few non-primate studies on lifestyle-related functional adaptation of trabecular bone. We show that the architecture of the femoral head trabeculae in Sciuromorpha correlates with body mass and locomotor habits. Our findings provide a new basis for experimental research focused on functional significance of bone inner microstructure.

Background

Sciuromorpha represent an expedient group for studying morphological adaptations to functional constraints, as they display diverse locomotor habits. They comprise more than 300 species with a wide range in body size [1] and which have adopted different locomotor behaviors (termed ’lifestyles’ in this study), which can be classified into arboreal, fossorial, aerial, and semifossorial. Whereas arboreal species live and nest in trees and spend little time on the ground, fossorial ones are adapted to a ground-dwelling lifestyle, digging burrows and nesting in the ground [2]. The fossorial lifestyle most probably evolved multiple times independently in the sciuromorph clade [2, 3]. Aerial (gliding) squirrels acquired a patagium enabling the animals to glide over long distances [4]. Members of the genus Tamias (chipmunks) adopted an intermediate lifestyle between arboreal and fossorial. They dig subterranean burrows for nesting, but climb trees when escaping predators [5]. They are thus called ’semifossorial’ in this study.

These lifestyle categories are expected to differ in directional variability, frequency and amount of loading. Arboreal habitats are diverse with respect to substrate slope, branch- or trunk diameter and texture. They require both climbing and jumping abilities, and impose diverse loading patterns on the limb bones (e.g. [6]). Aerial sciuromorphs, being mostly arboreal themselves, are found in habitats of similar diversity. However, gliding requires increased bone lightness and endurance in spanning the patagium. Additionally, they need to sporadically generate particularly large take-off forces and withstand particularly large landing forces (e.g. [7]). Fossorial species are expected to experience a more uniform loading than arboreal ones, as terrestrial locomotion and digging behavior are most demanding in developing high speeds and large unidirectional forces, respectively (e.g. [6]).

Different locomotor habits demand appropriate morphological adaptations to the characteristic loading patterns. Trabecular, or cancellous, bone enables such adaptations to loads acting on the skeleton of vertebrates [813]. The thin bony struts, called trabeculae, build an anisotropic microstructure that offers both stability and a light-weight construction [14]. The principal trabecular orientation often corresponds to the main stress trajectories within the bone (e.g. [15]) and therefore provides resistance against loads acting on it. The initial structure of trabecular bone in individuals is genetically regulated, and thus inherited [16]. However, trabecular bone adjusts to changes in loading during life with high sensitivity, which becomes visible in both the orientation and the dimensions of the trabeculae [12, 17]. Trabecular bone parameters furthermore scale allometrically [1820]. In humans and other primates, the trabecular architecture in long bones is related to locomotor patterns [13, 21, 22]. The trabecular bone structure has been used as an indicator for reconstruction of locomotor habits in extinct taxa [23, 24]. However, it remains controversial whether trabecular bone reflects lifestyle in primates, and if so, how. In contrast to the above mentioned findings, other studies in anthropoids reported the femoral trabecular microstructure to be independent of locomotor behavior (e.g. [25]). A broader sampling of vertebrate bone morphology is required to address the question of how trabecular architecture is related to locomotor behavior.

Since research on lifestyle-related functional adaptation of bone microstructure mainly focuses on primates, little has been done in non-primate mammals so far. Recent studies in Xenarthra [26] and Odontoceti [27] expanded the sampling of mammalian 3D trabecular architecture in context of lifestyle. In rodents, however, the functional role of trabeculae is mostly studied experimentally. In rats and mice, the effect of physical exercise on trabecular bone has been investigated (e.g. [8, 9]). Experiments in rats have shown that unloading of hind limbs causes bone loss and degradation towards a lighter trabecular structure within three weeks [28]. Trabecular bone loss has also been shown to occur in hibernating ground squirrels [29]. Yet research on trabecular bone in rodents mainly focuses on the hormonal influence on bone microstructure, rather than on the biomechanical constraints [3033].

This study investigates the effect of body mass and locomotor habits on femoral head trabecular bone structure in Sciuromorpha. Our aim is to elucidate the relationship of trabecular architecture to body mass and lifestyle-related loading conditions. In particular, we seek to identify whether the different loading regimes are reflected by trabecular architecture, testing the hypothesis that the more multidirectional loading in arboreal locomotion affects trabecular structure differently than the more unidirectional loading in fossorial lifestyle (cf. [6]). We furthermore aim at determining whether semifossorial species resemble rather the arboreal or fossorial sciuromorphs. Three-dimensional data analysis of high resolution computed tomography (CT-) scans was conducted in order to characterize trabecular bone in the femoral head of 69 sciuromorph species from different lifestyle categories. The body mass of these species ranged from approximately 16.5 g to 3000 g (Myosciurus pumilio [34] and Marmota marmota [35], respectively). A cubic volume of interest was selected in each femoral head for calculation of degree of anisotropy DA, main direction of trabeculae MDT (cf. [26]), connectivity density ConnD, mean trabecular thickness TbTh, mean trabecular separation TbSp, bone volume fraction BV/TV and bone surface density BS/BV. These trabecular parameters were used for analyses of the allometric and functional signal. We show that trabecular architecture relates to body mass. The data reveal functional signals related to lifestyle in some of the trabecular parameters. This study provides new findings on the significance of external loading patterns on trabecular bone structure in Sciuromorpha that will be valuable for experimental research aimed at elucidating the significance of bone microstructure in functional contexts.

Methods

Specimens

Sciuromorph femora have been chosen from the following museum collections for high resolution CT-scans: AMNH, American Museum of Natural History, New York City, USA; CeNaK, Centrum für Naturkunde, Universität Hamburg, Germany; MNHN, Muséum national d’Histoire naturelle, Paris, France; PMJ, Phyletisches Museum Jena, Germany; SMF, Naturmuseum Senckenberg, Frankfurt am Main, Germany; SMNS, Staatliches Museum für Naturkunde, Stuttgart, Germany; SNMNH, Smithsonian National Museum of Natural History, Washington, D.C., USA; UMMZ, University of Michigan Museum of Zoology, Ann Arbor, USA; ZMB, Museum für Naturkunde Berlin, Germany; ZSM, Zoologische Staatssammlung München, Germany.

See Additional file 1 for a detailed list of all specimens. CT-scans were acquired with resolutions of 6.4 - 42.5 μm voxel size and saved as 8-bit or 16-bit grayscale image stacks in tif format. Juvenile specimens (defined as those with an epiphyseal plate at the femoral head visible on the scans) were excluded from the analyses. After applying quality thresholds (described below), the final dataset covered all three sciuromorph families (Sciuridae, Gliridae, Aplodontiidae), including 33 genera and 69 species. Lifestyle categorization based on previous literature [3638] yielded the following composition (see Fig. 1): 27 arboreal species (representing Sciurillinae, Protoxerini, Sciurini, Ratufinae, Callosciurinae, Gliridae), 19 fossorial species (representing Xerini, Aplodontia rufa, Marmotini excl. Tamias), 15 semifossorial species (representing Tamias) and 8 aerial species (representing Pteromyini).

Fig. 1
figure 1

Phylogenetic tree of the sciuromorph species analyzed in this study. Colors indicate the lifestyle classification

3D data analysis procedure

Analysis of CT-scan data was done using the Fiji contribution of ImageJ vers. 1.51k [39, 40] and the plugin BoneJ vers. 1.4.2 [41]. An analysis Java script (Additional file 2) was composed and applied to all specimens. Prior to the analysis, each femur was oriented within the stack coordinate system using the ’Rotate’ function and the ’Reslice’ function in Fiji. The proximodistal axis of the femoral shaft was thereby aligned along the z-axis of the coordinate system (defined then as the proximodistal axis also for the femoral head). The mediolateral and anteroposterior axes were aligned along the x- and y-axis, respectively, using the femoral condyles as reference. A cubic volume of interest (VOI) was then selected in the center of the femoral head (Fig. 2). The central slice was defined as the plane that divides the femoral head into a distal and a proximal half (Fig. 2b). The VOI originated from this slice and was expanded as a cube in all three dimensions as much as possible without including cortical bone. The center of the VOI on the mid-slice was thereby positioned such that it allows for a maximum spread over the eight corners of the VOI. Scans with an insufficient resolution were excluded. We used a relative resolution as the number of pixels covering the mean trabecular thickness TbTh (relative resolution = TbTh /pixel size) [42, 43] and discarded all scans which fell below a relative resolution of five (cf. [42]). To allow for application of phylogenetically informed methods, we included only those species in the analysis for which the exact position within the phylogenetic tree is known (shown in Fig. 1).

Fig. 2
figure 2

Workflow for the selection of the cubic volume of interest (VOI) in the femoral head. Selection of the VOI is shown with a specimen of Tamias rufus (SNMNH, 564127). a Micro-CT coronal slice, 16-bit grayscale image. b Selection of the cubic VOI in the center of the femoral head. The horizontal orange line indicates the central slice which divides the femoral head and the VOI into a proximal and a distal half. The VOI was expanded in three dimensions as much as possible without including cortical bone. c 3D view of the final selection after cropping the image stack. Analysis of trabecular architecture was executed on this substack. Note that these pictures serve as a demonstration for the VOI selection but do not conform with the proper analysis done on this specimen. Scale bars: 1 mm

All further analyses were done using the BoneJ plugin. For preparation the VOI was binarized (using the ’Optimise Threshold’ function) and purified (’Purify’ function). The calculated parameters for characterization of the trabecular architecture were the following: degree of anisotropy DA with vector of main direction of trabeculae MDT (’Anisotropy’ function), connectivity density ConnD (’Connectivity’ function), mean trabecular thickness TbTh and trabecular separation TbSp (’Thickness’ function) as well as bone volume fraction BV/TV and bone surface density BS/BV (’Volume Fraction’ function and ’Isosurface’ function, respectively). DA is calculated in BoneJ by using the mean intercept length method [44]. The parameter indicates to what extent the trabeculae are oriented in a preferential direction. The measure as used here ranges from zero (fully isotropic, no preferred direction) to one (fully anisotropic). The main orientation of trabeculae MDT is given by the eigenvector corresponding to the lowest eigenvalue, which reflects the major axis of the ellipsoid defining the trabecular anisotropy [44]. All vectors were projected onto the proximal hemisphere of the femoral head and those of right femora were mirrored to the left side. MDT is a three-dimensional parameter (x,y,z) with x describing the mediolateral, y the anteroposterior and z the proximodistal contribution to the main orientation. Note that MDT can be compared among the specimens because the femora are all oriented in the same reference coordinate system given by the image axes. ConnD is calculated in BoneJ using the Euler characteristic Δχ with ConnD=(1−Δχ)/stack volume. The parameter represents a measure for the number of trabeculae per volume [45]. Accurate analysis of trabeculae implies a continuum assumption for the trabecular bone, which requires a minimal number of trabeculae within the VOI [46]. Thus all samples with an absolute connectivity (Conn) of less than 50 were discarded. TbTh and TbSp are calculated as the mean thickness of the fore- and background in the 3D VOI stack. Thickness at a point is thereby defined as the diameter of the greatest sphere fitting within fore- or background structure and containing this point. BV/TV is measured using a voxel-based algorithm that calculates the ratio of foreground voxels (representing bone tissue) and total number of voxels in the sample. The measure gives the proportion of volume that is occupied by bone. BS/BV was calculated by measuring the bone surface area and dividing it by the bone volume of the sample.

Phylogenetic tree

The sciuromorph phylogenetic tree used herein (Fig. 1) is based on a previous study [47] and complemented with species from the TimeTree database [48] (Gliridae, Tamiops swinhoei, Tamiops rodolphii, Aplodontia rufa, Heliosciurus rufobrachium) using Mesquite version 3.04 [49]. Given the apparent correlation of the lifestyle categories with the phylogeny (see Fig. 1), we conducted subsequent analyses with phylogenetically informed methods.

Allometry analysis

For all trabecular parameters the relation with VOI edge length vl was tested. Since the actual individual body masses were unknown for the specimens, vl served as a proxy for body mass. A test with available body masses taken from literature for a subset of 64 species (93% of the species covered in this study) revealed significant correlation of vl with body mass (r2=0.92, p<10−44 with r being the Pearson correlation coefficient). Most importantly, vl scales isometrically with body mass m. We found mvl2.95 with mvl3 being the expected scaling under isometry. The observed scaling (exponent 2.95) did not differ significantly from the isometric (exponent 3) case (p=0.56). Thus, we consider vl a suitable proxy for body mass, a predominant factor in loading.

The relation of any trabecular parameter tp with vl was quantified by calculating the scaling exponent a such that tpvla and accordingly log(tp)a·log(vl). All raw data were natural log-transformed beforehand. Linear regression analysis was performed in R vers. 3.4.2 [50] using the gls()-function in the ’nlme’ package vers. 3.1.131 [51]. Regressions were done phylogenetically informed (corPagel()-function in the ’ape’ package vers. 5.0 [52]). The linear correlation between the log-transformed data and log-transformed vl was quantified by calculating the slope a and the coefficient of determination r2. To identify the type of scaling (negative or positive allometry or isometry), we tested for significant deviation of the observed scaling exponent (a obs ) from the expected scaling exponent for the case of isometry (a iso ). To do so, we set the isometric regression line as baseline (such that the isometric slope becomes zero for all parameters) by subtracting the value expected under isometry (vl·a iso ) from each data point. We then performed the regression as described above, testing the null hypothesis a obs =0. Slopes significantly higher or lower than zero indicate positive or negative allometry, respectively. Slopes not significantly deviating from zero indicate isometry or near-isometry. For BV/TV and DA, a iso =0 (Table 1), as these measures are dimensionless ratios, which do not change with increasing size under assumption of isometry. On the other hand TbTh and TbSp (unit: length1) are absolute linear measures which would increase under isometry with a iso =1. For BS/BV (unit: length−1) applies a iso =−1 (the volume BV increases more rapidly than BS). The isometric scaling exponent for ConnD (unit: length−3) is a iso =−3. Under isometric growing of the trabeculae (model of scaling with constant trabecular geometry [53]), fewer trabeculae fit into one unit volume with increasing size. For each parameter, a p-value was calculated for testing the null hypothesis a obs =a iso as described above. Thus, isometry was assumed when p>0.05. Otherwise we assumed positive allometry for a obs >a iso and negative allometry for a obs <a iso .

Table 1 Allometry in the sciuromorph femoral head trabeculae - relation between trabecular parameters and VOI edge length

Lifestyle analysis

In order to test for differences in trabecular architecture among the four lifestyle categories (arboreal, fossorial, aerial and semifossorial), we performed a pairwise statistical analysis on the individual trabecular parameters.

The significance level for all statistical tests was fixed at α=0.05. Each pair of lifestyles was compared for each tp by regression in R vers. 3.4.2 [50] with the gls()-function in the ’nlme’ package vers. 3.1.131 [51] while taking correlation with phylogeny into account (corPagel()-function in the ’ape’ package vers. 5.0 [52]). Since all parameters were correlated with vl, the body mass proxy was included as covariate in the function to separate the influence of lifestyle and vl on the parameter (R formula: tpvl+lifestyle). Because this procedure of multiple statistical tests increases the rate of type I errors, we corrected the p-values with a Benjamini-Hochberg method [54] (p.adjust()-function in the ’stats’ package in R). Only the corrected p-values are reported herein. Since this analysis procedure takes correlation with phylogeny and body mass into account, it is ensured that the confounding effect of these factors is reduced. Hence, we are able to analyze and interpret the data with regard to an effect of lifestyle on trabecular properties.

Results

Main orientations of trabeculae are distributed around anteroposterior axis

In order to investigate how the trabeculae are oriented within the femoral head, their main direction MDT was displayed as vectors in spherical coordinates (Fig. 3). These vectors, originating in the center of the VOI, were projected onto the proximal hemisphere of the femoral head. We observed no apparent differences in MDT among the lifestyle categories (Fig. 3). The trabeculae of most specimens were oriented in the proximo-antero-medial direction. The main variation in the whole dataset was observed along the anteroposterior axis.

Fig. 3
figure 3

Trabecular orientation. The main direction of trabecular anisotropy (MDT) is represented by vectors originating in the center of the femoral head and projected onto its proximal hemisphere. A femur, seen in proximal view, is depicted at lower right to illustrate the orientation (g.tr./l.tr., greater/lesser trochanter; f.h., femoral head). Data points close to the center are oriented more proximodistally. Colors indicate the lifestyle categories, circle size the degree of anisotropy (DA), ranging from 0.451 (smallest circle) to 0.883 (largest circle). The data show the results for left femora (right femora were mirrored)

Femoral head trabecular bone parameters scale allometrically

To test whether the trabecular parameters depend on body mass, the correlation of each of them with the edge length vl of the cubic VOI was analyzed by calculating the coefficient of determination r2 and the scaling exponent a. To identify allometric or isometric scaling, we compared the observed scaling exponents with the expected scaling exponents for the case of isometry (a iso , see Table 1).

We observed positive allometry for BV/TV, ConnD and BS/BV and negative allometry for DA, TbTh and TbSp (Table 1, Fig. 4). We conclude from the negative allometry in absolute TbTh and TbSp that relative TbTh and TbSp were lower with higher vl. Thus, trabeculae tend to be relatively thinner and more densely packed in larger animals (higher ConnD). Comparison of the scaling of TbSp and TbTh revealed that TbTh increased more rapidly with vl than TbSp did, which is consistent with the positive allometry in BV/TV. Taken together, these results reveal a significant allometric signal in all six trabecular parameters.

Fig. 4
figure 4

Allometric scaling in trabecular parameters. Each plot shows the relation between the trabecular parameter and edge length of the cubic VOI (vl), the latter serving as proxy for body mass (see text). Raw data were plotted on log-log scale. r2 - coefficient of determination, p - p-value of the regression. a-c Degree of anisotropy DA, trabecular thickness TbTh and trabecular separation TbSp scale with negative allometry. d-f Bone volume fraction BV/TV, connectivity density ConnD and bone surface density BS/BV scale with positive allometry

Femoral head trabecular bone reflects lifestyle

To analyze how the individual trabecular parameters relate to lifestyle, a phylogenetically informed pairwise statistical analysis of the log-transformed data was performed by comparing each two of the four lifestyle groups while accounting for the effect of mass. We identified functional signals in TbTh, ConnD, BS/BV, and TbSp (Fig. 5a-d): Aerial species are distinguished from fossorial ones by a higher TbTh (p(ae,fo)=2×10−7), a lower ConnD (p(ae,fo)=0.0007) and a lower BS/BV (p(ae,fo)=2.2×10−5). Arboreal species have a higher TbSp than semifossorial species have (p(ar,sf)=4.1×10−5). DA and BV/TV showed no significant differences among lifestyles (Fig. 5e-f). However, although not being significant, DA tends to be higher in arboreal species. When considering a ’lifestyle gradient’ as the sequence aerial > arboreal > semifossorial > fossorial, we observe the general trend of decreasing TbTh and TbSp and increasing ConnD and BS/BV towards the fossorial lifestyle (Fig. 5). When comparing semifossorial species with arboreal and fossorial ones, we observed that they generally resemble the fossorial sciuromorphs more closely than they do the arboreal ones (mostly apparent in ConnD, TbSp and DA, see Fig. 5). Taken together, these results reveal lifestyle-related differences in four of six analyzed trabecular parameters (TbTh, ConnD, BS/BV, and TbSp).

Fig. 5
figure 5

Lifestyle-related adaptations of trabecular parameters. The boxplots illustrate distributions of the size-corrected (s.c., residuals after allometry regression) data of the trabecular parameters in the lifestyle groups aerial ae, arboreal ar, semifossorial sf and fossorial fo. Significant differences between groups, as revealed by phylogenetically informed regression, are indicated by the horizontal lines (p-values were corrected for multiple testing, for detailed explanation see text). a-c Aerial species are distinguished from fossorial ones by a higher trabecular thickness TbTh, lower connectivity density ConnD and lower bone surface density BS/BV. d Trabecular separation TbSp is higher in arboreal species than in semifossorial ones. e-f Degree of anisotropy DA and bone volume fraction BV/TV do not show significant functional signals

Discussion

The aim of this study was to characterize trabecular microstructure of the sciuromorph femoral head and to investigate its relation to body mass and locomotor behavior. The sciuromorph trabeculae are preferentially oriented in proximo-antero-medial direction within the femoral head and are distributed around its anteroposterior axis. The tested trabecular parameters show different scaling effects: whereas DA, TbTh and TbSp scale with negative allometry, a positive allometry was revealed for BV/TV, ConnD and BS/BV. Phylogenetically informed analysis identified functional signals in TbTh, TbSp, ConnD and BS/BV. It distinguished aerial Sciuromorpha from fossorial ones (higher TbTh, lower ConnD and lower BS/BV in aerial species) and arboreal Sciuromorpha from semifossorial ones (higher TbSp in arboreal species). Semifossorial sciuromorphs resembled the fossorial species rather than the arboreal ones.

The observed proximomedial orientation of the trabeculae is on average roughly parallel to the femoral neck. The anteroposterior distribution of main trabecular orientation is in agreement with observations in primate femora [10, 22]. We suggest that this pattern of trabecular orientation may indicate that the predominant pro-/retraction angle at high loading phases during locomotion is a determining factor of this trabecular orientation, both in primates and sciuromorphs.

Allometric scaling has been observed in trabecular architecture of different bones in a large variety of taxa [1820]. The allometric scaling patterns of trabecular parameters observed for the sciuromorph femoral head trabeculae conform with observations in primates [19]. Other studies reported BV/TV and DA to be independent of body size [18, 20]. Our result for the scaling of ConnD (a=−2.138) qualitatively conforms with the observation reported for the primate femur (a=−1.22, see [19]). However, in that previous study, ConnD has been considered as a shape variable which does not change with body size and scales with a=0 in the isometric case. Hence, the scaling has been reported as negative allometry. This would be the case if trabecular size (e.g. TbTh) remains constant with increasing body size and solely the trabecular number increases (model of ’constant trabecular size’ as against ’constant trabecular geometry’ [53]). If ConnD is interpreted as number of trabeculae per volume, as herein, and trabecular size scales isometrically, the isometric scaling exponent should be a=−3 and the scaling be positively allometric. Positive allometry in BV/TV, as observed herein, is also observed in other rodents, mainly due to variation in trabecular number rather than variation in TbTh, the latter being the main reason in humans [20].

Previous studies drew different conclusions as to whether trabecular bone structure reflects locomotor habits. The femoral neck trabecular architecture in nonhuman anthropoid primates has been reported by one study to be independent of locomotor mode, phylogenetic background and even body mass [25]. In contrast, other studies, including the present one, observed a relation between trabecular bone and locomotion [10, 13, 26, 55, 56]. We identified TbTh, TbSp, ConnD, and BS/BV in the sciuromorph femoral head microstructure to be correlated with locomotor habits. It has also been shown that the proximal femur of human foraging populations has higher BV/TV and TbTh and a lower BS/BV compared to agriculturalist groups, revealing a correlation between high mobility and relative bone stiffness [13]. Other comparative studies observed functional adaptations related to leaping behavior in BV/TV, DA or MDT in different primate species [10, 55, 56]. Furthermore, it has been shown that trabecular bone reflects basic primate locomotor categories, such as brachiation, quadrupedal/bipedal walking, climbing, and jumping [21, 22]. Trabecular bone has already been used in reconstructing locomotor modes of fossil primates [23, 56].

However, the approach of using bone microstructure as an indicator for locomotor habits has limitations. In addition to body mass and lifestyle, ontogenetic stage of an animal might also influence trabecular architecture (e.g., [57]). In collection-based research age cannot be considered, as the precise age of death of the specimens is usually not known. However, it has been shown in mice that trabecular parameters change mostly during the first six months of life and then pass into a stable phase [58]. By excluding specimens whose epiphyseal plate was not fused, we avoided including too young specimens that might not have reached that state. For studying allometry, an estimation of body mass is desirable. Equations for that purpose, based upon bone articular and diaphyseal structure, exist but should be used cautiously [59]. Thus, directly available measures, such as femoral head dimensions, are commonly used as body size proxy (e.g., [18, 19]). In this study, we decided to use the length of the VOI (vl) as a proxy for body size, as we could show that this parameter scales isometrically with body mass, a predominant factor in loading. The architecture of trabecular bone is highly variable not only among species but even among individuals of the same species [24, 53]. Our dataset, composed of solely one or two individuals per species, does not allow an assessment of this intraspecies variance. However, since the dataset comprises a broad sampling on the level of genera, lifestyle and body mass, it still allows us to draw conclusions about trabecular structure in a general functional context.

This study conducted research on lifestyle-related functional adaptation of trabecular architecture in Sciuromorpha and thereby supplements the few non-primate studies of that kind. The different body sizes and locomotor habits within the Sciuromorpha account for naturally varying characteristic loading patterns acting on the bones of the animals. These loading patterns conferred structural differences in the femoral head trabecular microstructure, as revealed by phylogenetically informed GLS regression.

The observed trends support our hypothesis that the loading patterns related to climbing and digging behavior are differently reflected by trabecular parameters. The apparent bone density is related to bone stability [60] and biomechanical loading [13]. This may be due to the fact that trabecular bone, a material which needs to absorb energy during impact loading, can increase stability and absorbing capacity through increasing the apparent bone density [61]. This can be achieved through (1) an increase in trabecular number, as typical for rodents, or (2) an increase in trabecular thickness, as typical for humans [13, 20]. When considering a ’lifestyle gradient’ (without relation to the evolutionary history) as the sequence aerial > arboreal > semifossorial > fossorial, we observed the mechanism (1) (increase in trabecular number, thus in ConnD) as a trend toward fossorial sciuromorphs and mechanism (2) (increase in TbTh) as a trend towards aerial sciuromorphs. Which functional constraints constitute these different mechanisms of structural adaptation? On the one hand, in aerial sciuromorphs, bones (while needing to retain light-weight structure) need to resist high peak loads when the animals land after gliding [7]. Thus, the bone is loaded with a low frequency, but with a high magnitude, a loading pattern that particularly triggers adaptive response in bones [62, 63]. Here, the need to avoid fractures within the trabecular network may dominate and account for the mechanism of increased trabecular thickness in aerial sciuromorphs (cf. [22]). On the other hand, fossorial species are expected to experience lower peak forces, but at a higher frequency. The hindlimbs provide propulsion during running and ensure stability during digging activity of the forelimbs [64]. Here, the need to retain more sustained material stability against moderate loading may account for the mechanism of increased trabecular number (cf. [11, 21]). Thus, trabeculae in fossorial sciuromorphs are thinner and more closely packed (as observed as a trend in semifossorial species as well), a pattern that we found also for large-sized sciuromorphs when analyzing allometric scaling of trabeculae. This comes along with an increased BS/BV, which may facilitate transfer of minerals and calcium mobilization at the surfaces of trabeculae [53] and might thus be an adaptive mechanism observed both in large sized animals and in animals with particularly sustained physical activity.

Conclusions

This study characterizes trabecular architecture and its relation to body mass and locomotor behavior in the sciuromorph femoral head. The results reveal allometric scaling in the set of trabecular parameters analyzed herein. The trabecular microstructure reflects lifestyle-related functional signals. Similar results have been documented in a few mammalian clades (primates [13, 21, 22], Odontoceti [27], Xenarthra [26]) and thus suggest to provide important insight into the functional morphology of mammals in general.

Subsequent research ought to focus on identifying the different mechanisms of functional adaptation in trabecular microstructure in animals of different body sizes and lifestyles. Our findings could be verified with experimental studies including setups that imitate different habitats and therefore different loading patterns applied on living Sciuromorpha. Thus, one could investigate experimentally whether trabecular parameters are indeed differently affected by loading patterns related to climbing and digging, as suggested in this study. This would further elucidate the functional significance of bone internal structure in Sciuromorpha.

Abbreviations

BS/BV:

Bone surface density

BV/TV:

Bone volume fraction

ConnD:

Connectivity density

CT:

computed tomography

DA:

Degree of anisotropy

MDT:

Main direction of trabeculae

TbSp:

Trabecular separation

TbTh:

Trabecular thickness

tp:

Trabecular parameter

vl:

VOI edge length

VOI:

Volume of interest

References

  1. Wilson DE, Reeder DM. Mammal species of the world: a taxonomic and geographic reference. vol. 1 3rd ed.Baltimore: JHU Press; 2005.

    Google Scholar 

  2. Steppan SJ, Storz BL, Hoffmann RS. Nuclear DNA phylogeny of the squirrels (Mammalia: Rodentia) and the evolution of arboreality from c-myc and RAG1. Mol Phylogenet Evol. 2004; 30(3):703–19.

    Article  PubMed  CAS  Google Scholar 

  3. Rocha RG, Leite YLR, Costa LP, Rojas D. Independent reversals to terrestriality in squirrels (Rodentia: Sciuridae) support ecologically mediated modes of adaptation. J Evol Biol. 2016; 29(12):2471–9.

    Article  PubMed  CAS  Google Scholar 

  4. Thorington RW, Heaney LR. Body proportions and gliding adaptations of flying squirrels (Petauristinae). J Mammal. 1981; 62(1):101–14.

    Article  Google Scholar 

  5. Clarke MF, da Silva KB, Lair H, Pocklington R, Kramer DL, McLaughlin RL. Site Familiarity Affects Escape Behaviour of the Eastern Chipmunk, Tamias striatus. Oikos. 1993; 66(3):533–7.

    Article  Google Scholar 

  6. Burr DB, Ruff CB Johnson. Structural adaptations of the femur and humerus to arboreal and terrestrial environments in three species of macaque. Am J Phys Anthropol. 1989; 79(3):357–67.

    Article  PubMed  CAS  Google Scholar 

  7. Paskins KE, Bowyer A, Megill WM, Scheibe JS. Take-off and landing forces and the evolution of controlled gliding in northern flying squirrels Glaucomys sabrinus. J Exp Biol. 2007; 210(8):1413–23.

    Article  PubMed  Google Scholar 

  8. Søgaard CH, Danielsen CC, Thorling EB, Mosekilde L. Long-term exercise of young and adult female rats: Effect on femoral neck biomechanical competence and bone structure. J Bone Miner Res. 1994; 9(3):409–16.

    Article  PubMed  Google Scholar 

  9. Joo YI, Sone T, Fukunaga M, Lim SG, Onodera S. Effects of endurance exercise on three-dimensional trabecular bone microarchitecture in young growing rats. Bone. 2003; 33(4):485–93.

    Article  PubMed  Google Scholar 

  10. Ryan TM, Ketcham RA. Angular orientation of trabecular bone in the femoral head and its relationship to hip joint loads in leaping primates. J Morphol. 2005; 265(3):249–63.

    Article  PubMed  Google Scholar 

  11. Plochocki JH, Riscigno CJ, Garcia M. Functional adaptation of the femoral head to voluntary exercise. Anat Rec A: Discov Mol Cell Evol Biol. 2006; 288(7):776–81.

    Article  Google Scholar 

  12. Barak MM, Lieberman DE, Hublin JJ. A Wolff in sheep’s clothing: Trabecular bone adaptation in response to changes in joint loading orientation. Bone. 2011; 49(6):1141–51.

    Article  PubMed  Google Scholar 

  13. Ryan TM, Shaw CN. Gracility of the modern Homo sapiens skeleton is the result of decreased biomechanical loading. Proc Nat Acad Sci. 2014; 112(2):372–7.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  14. Galante J, Rostoker W, Ray R. Physical properties of trabecular bone. Calcif Tissue Res. 1970; 5(1):236–46.

    Article  PubMed  CAS  Google Scholar 

  15. Venieratos D, Papadopoulos N, Anastassiou J, Katritsis E. A quantitative estimation of the divergence between the trabecular system and the stress trajectories in the upper end of the human femoral bone. Anat Anz. 1987; 163(4):301–10.

    PubMed  CAS  Google Scholar 

  16. Turner CH, Hsieh YF, Müller R, Bouxsein ML, Baylink DJ, Rosen CJ, et al. Genetic regulation of cortical and trabecular bone strength and microstructure in inbred strains of mice. J Bone Miner Res. 2000; 15(6):1126–31.

    Article  PubMed  CAS  Google Scholar 

  17. Pontzer H, Lieberman D, Momin E, Devlin MJ, Polk J, Hallgrimsson B, et al. Trabecular bone in the bird knee responds with high sensitivity to changes in load orientation. J Exp Biol. 2006; 209(1):57–65.

    Article  PubMed  CAS  Google Scholar 

  18. Doube M, Kłosowski MM, Wiktorowicz-Conroy AM, Hutchinson JR, Shefelbine SJ. Trabecular bone scales allometrically in mammals and birds. Proc R Soc B Biol Sci. 2011; 278(1721):3067–73.

    Article  Google Scholar 

  19. Ryan TM, Shaw CN. Trabecular bone microstructure scales allometrically in the primate humerus and femur. Proc R Soc B Biol Sci. 2013; 280(1758):20130172.

    Article  CAS  Google Scholar 

  20. Barak MM, Lieberman DE, Hublin JJ. Of mice, rats and men: Trabecular bone architecture in mammals scales to body mass with negative allometry. J Struct Biol. 2013; 183(2):123–31.

    Article  PubMed  Google Scholar 

  21. Ryan TM, Shaw CN. Unique suites of trabecular bone features characterize locomotor behavior in human and non-human anthropoid primates. PLoS ONE. 2012; 7(7):e41037.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  22. Saparin P, Scherf H, Hublin JJ, Fratzl P, Weinkamer R. Structural adaptation of trabecular bone revealed by position resolved analysis of proximal femora of different primates. Anat Rec. 2010; 294(1):55–67.

    Article  Google Scholar 

  23. Rook L, Bondioli L, Köhler M, Moyà-Solà S, Macchiarelli R. Oreopithecus was a bipedal ape after all: evidence from the iliac cancellous architecture. Proc Nat Acad Sci. 1999; 96(15):8795–9.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  24. Saers JP, Cazorla-Bak Y, Shaw CN, Stock JT, Ryan TM. Trabecular bone structural variation throughout the human lower limb. J Hum Evol. 2016; 97:97–108.

    Article  PubMed  Google Scholar 

  25. Fajardo RJ, Müller R, Ketcham RA, Colbert M. Nonhuman anthropoid primate femoral neck trabecular architecture and its relationship to locomotor mode. Anat Rec. 2007; 290(4):422–36.

    Article  Google Scholar 

  26. Amson E, Arnold P, van Heteren AH, Canoville A, Nyakatura JA. Trabecular architecture in the forelimb epiphyses of extant xenarthrans (Mammalia). Front Zool. 2017; 14(1):52.

    Article  PubMed  PubMed Central  Google Scholar 

  27. Rolvien T, Hahn M, Siebert U, Püschel K, Wilke HJ, Busse B, et al. Vertebral bone microarchitecture and osteocyte characteristics of three toothed whale species with varying diving behaviour. Sci Rep. 2017; 7:1604.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  28. David V, Laroche N, Boudignon B, Lafage-Proust MH, Alexandre C, Ruegsegger P, et al. Noninvasive in vivo monitoring of bone architecture alterations in hindlimb-unloaded female rats using novel three-dimensional microcomputed tomography. J Bone Miner Res. 2003; 18(9):1622–31.

    Article  PubMed  Google Scholar 

  29. McGee-Lawrence ME, Stoll DM, Mantila ER, Fahrner BK, Carey HV, Donahue SW. Thirteen-lined ground squirrels (Ictidomys tridecemlineatus) show microstructural bone loss during hibernation but preserve bone macrostructural geometry and strength. J Exp Biol. 2011; 214(8):1240–7.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Miller SC, Wronski T. Long-term osteopenic changes in cancellous bone structure in ovariectomized rats. Anat Rec. 1993; 236(3):433–41.

    Article  PubMed  CAS  Google Scholar 

  31. Dempster D, Birchman R, Xu R, Lindsay R, Shen V. Temporal changes in cancellous bone structure of rats immediately after ovariectomy. Bone. 1995; 16(1):157–61.

    Article  PubMed  CAS  Google Scholar 

  32. Bagi C, Ammann P, Rizzoli R, Miller S. Effect of estrogen deficiency on cancellous and cortical bone structure and strength of the femoral neck in rats. Calcif Tissue Int. 1997; 61(4):336–44.

    Article  PubMed  CAS  Google Scholar 

  33. Waarsing JH, Day JS, Verhaar JA, Ederveen AG, Weinans H. Bone loss dynamics result in trabecular alignment in aging and ovariectomized rats. J Orthop Res. 2006; 24(5):926–35.

    Article  PubMed  Google Scholar 

  34. Emmons LH. A Note on the Forefoot of Myosciurus pumilio. J Mammal. 1979; 60(2):431.

    Article  Google Scholar 

  35. Arnold W, Heldmaier G, Ortmann S, Pohl H, Ruf T, Steinlechner S. Ambient temperatures in hibernacula and their energetic consequences for alpine marmots Marmota marmota. J Therm Biol. 1991; 16(4):223–6.

    Article  Google Scholar 

  36. Thorington Jr RW, Koprowski JL, Steele MA, Whatton JF. Squirrels of the world. Baltimore: JHU Press; 2012.

    Google Scholar 

  37. Carraway LN, Verts BJ. Aplodontia rufa. Mamm Species. 1993; 431:1–10.

    Google Scholar 

  38. Nowak RM. Walker’s Mammals of the World, vol. 1. 1999.

  39. Schneider CA, Rasband WS, Eliceiri KW. NIH Image to ImageJ: 25 years of image analysis. Nat Methods. 2012; 9(7):671.

    Article  PubMed  PubMed Central  CAS  Google Scholar 

  40. Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, et al. Fiji: an open-source platform for biological-image analysis. Nat Methods. 2012; 9(7):676–82.

    Article  PubMed  CAS  Google Scholar 

  41. Doube M, Kłosowski MM, Arganda-Carreras I, Cordelières FP, Dougherty RP, Jackson JS, et al. BoneJ: Free and extensible bone image analysis in ImageJ. Bone. 2010; 47(6):1076–9.

    Article  PubMed  PubMed Central  Google Scholar 

  42. Sode M, Burghardt AJ, Nissenson RA, Majumdar S. Resolution dependence of the non-metric trabecular structure indices. Bone. 2008; 42(4):728–36.

    Article  PubMed  Google Scholar 

  43. Kivell TL, Skinner MM, Lazenby R, Hublin JJ. Methodological considerations for analyzing trabecular architecture: an example from the primate hand. J Anat. 2011; 218(2):209–25.

    Article  PubMed  Google Scholar 

  44. Harrigan T, Mann R. Characterization of microstructural anisotropy in orthotropic materials using a second rank tensor. J Mater Sci. 1984; 19(3):761–7.

    Article  CAS  Google Scholar 

  45. Odgaard A, Gundersen HJ. Quantification of connectivity in cancellous bone, with special emphasis on 3-D reconstructions. Bone. 1993; 14:173–82.

    Article  PubMed  CAS  Google Scholar 

  46. Harrigan TP, Jasty M, Mann RW, Harris WH. Limitations of the continuum assumption in cancellous bone. J Biomech. 1988; 21(4):269–75.

    Article  PubMed  CAS  Google Scholar 

  47. Zelditch ML, Li J, Tran LAP, Swiderski DL. Relationships of diversity, disparity, and their evolutionary rates in squirrels (Sciuridae). Evolution. 2015; 69(5):1284–300.

    Article  PubMed  Google Scholar 

  48. Kumar S, Stecher G, Suleski M, Hedges SB. TimeTree: a resource for timelines, timetrees, and divergence times. Mol Biol Evol. 2017; 34(7):1812–1819.

    Article  PubMed  Google Scholar 

  49. Maddison WP, Maddison DR. Mesquite: a modular system for evolutionary analysis. 2017. http://mesquiteproject.org. Accessed 6 Dec 2017.

  50. R Core Team. R: A Language and Environment for Statistical Computing. Vienna; 2017. https://www.R-project.org/. Accessed 6 Dec 2017.

  51. Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team. nlme: Linear and Nonlinear Mixed Effects Models. 2017. R package version 3.1-131. https://CRAN.R-project.org/package=nlme. Accessed 6 Dec 2017.

  52. Paradis E, Blomberg S, Bolker B, Brown J, Claude J, Cuong HS, et al. ape: Analyses of Phylogenetics and Evolution; 2017. R package version 5.0. https://CRAN.R-project.org/package=ape. Accessed 6 Dec 2017.

  53. Swartz SM, Parker A, Huo C. Theoretical and empirical scaling patterns and topological homology in bone trabeculae. J Exp Biol. 1998; 201(4):573–90.

    PubMed  CAS  Google Scholar 

  54. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995; 57(1):289–300.

    Google Scholar 

  55. MacLatchy L, Müller R. A comparison of the femoral head and neck trabecular architecture of Galago and Perodicticus, using micro-computed tomography (μCT). J Hum Evol. 2002; 43(1):89–105.

    Article  PubMed  Google Scholar 

  56. Ryan TM, Ketcham RA. Femoral head trabecular bone structure in two omomyid primates. J Hum Evol. 2002; 43(2):241–63.

    Article  PubMed  Google Scholar 

  57. Acquaah F, Brown KAR, Ahmed F, Jeffery N, Abel RL. Early trabecular development in human vertebrae: overproduction, constructive regression, and refinement. Front Endocrinol. 2015; 6:67.

    Article  Google Scholar 

  58. Glatt V, Canalis E, Stadmeyer L, Bouxsein ML. Age-related changes in trabecular architecture differ in female and male C57BL/6J mice. J Bone Miner Res. 2007; 22(8):1197–207.

    Article  PubMed  Google Scholar 

  59. Ruff CB. Long bone articular and diaphyseal structure in Old World monkeys and apes. II: estimation of body mass. Am J Phys Anthropol. 2003; 120(1):16–37.

    Article  PubMed  Google Scholar 

  60. van Lenthe GH, Müller R. Prediction of failure load using micro-finite element analysis models: toward in vivo strength assessment. Drug Discov Today Technol. 2006; 3(2):221–9.

    Article  PubMed  Google Scholar 

  61. Currey JD. Chpt. 4.7. In: The mechanical adaptations of bones. Princeton: Princeton University Press: 1984. p. 133–4.

    Chapter  Google Scholar 

  62. Brodt MD, Silva MJ. Aged mice have enhanced endocortical response and normal periosteal response compared with young-adult mice following 1 week of axial tibial compression. J Bone Miner Res. 2010; 25(9):2006–15.

    Article  PubMed  PubMed Central  Google Scholar 

  63. Sugiyama T, Meakin LB, Browne WJ, Galea GL, Price JS, Lanyon LE. Bones’ adaptive response to mechanical loading is essentially linear between the low strains associated with disuse and the high strains associated with the lamellar/woven bone transition. J Bone Miner Res. 2012; 27(8):1784–93.

    Article  PubMed  PubMed Central  Google Scholar 

  64. Lehmann WH. The forelimb architecture of some fossorial rodents. J Morphol. 1963; 113(1):59–76.

    Article  PubMed  CAS  Google Scholar 

Download references

Acknowledgments

We are grateful to Miriam Zelditch for providing the phylogeny. We thank Falk Mielke for debugging and contributing parts of the analysis scripts and for valuable feedback on the manuscript. We are thankful to Karl Jepsen and Ingrid Hilger who gave us access to their CT-scan facilities and Lauran M. Smith and Brigitte Grobis for their help with the scanning procedure. We want to thank the following curators, collection managers and research assistants who gave us access to their collections: Violaine Nicolas (MNHN); Martin S. Fischer and Mathias Krüger (PMJ); Frieder Mayer and Christiane Funk (ZMB); Thomas M. Kaiser (CeNaK); Stefan Merker (SMNS); Irina Ruf and Katrin Krohmann (SMF); Darrin P. Lunde, Nicole R. Edmison and Esther M. Langan (SNMNH); Neil Duncan and Eileen Westwig (AMNH); Cody Thompson (UMMZ); Laura Abraczinskas (MSU) and Lawrence Heaney (FM). We thank two anonymous reviewers for their detailed comments and suggestions which have greatly improved this manuscript.

Funding

The project was funded by the DFG (EXC 1027 Bild Wissen Gestaltung. Ein interdisziplinäres Labor). EA was funded by the Alexander von Humboldt Foundation. PA was funded by the Max Planck Society.

Availability of data and materials

The dataset generated during the current study is provided in Additional file 1.

Author information

Authors and Affiliations

Authors

Contributions

MM conducted the 3D data analysis, subsequent statistical analysis and wrote the manuscript. EA and JAN conceived of this study. EA, JAN and JW gave crucial input regarding methods and analysis. EA, PA, AHvH and JW acquired the CT scans. EA and JW prepared the phylogenetic tree. All authors gave valuable feedback on the manuscript and approved the final version for submission.

Corresponding author

Correspondence to Maja Mielke.

Ethics declarations

Ethics approval and consent to participate

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

List of all analyzed specimens and their trabecular parameters. Complete list of the specimens analyzed in this study (.csv-file), including information about their origin (museum and collection number), their relative resolution, the CT-scanner, their sex, their assigned lifestyle category and the raw data characterizing the VOIs (MDT, DA, ConnD, TbTh, TbSp, BV/TV, BS/BV). (CSV 13 kb)

Additional file 2

3D data analysis Java script for BoneJ in Fiji. Java script, implemented for the extraction of trabecular parameters of an image stack with the BoneJ plugin in Fiji. With BoneJ being installed it is ready to use after opening it in Fiji (File/Open...). (0.763 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License(http://creativecommons.org/licenses/by/4.0/), 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(http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Mielke, M., Wölfer, J., Arnold, P. et al. Trabecular architecture in the sciuromorph femoral head: allometry and functional adaptation. Zoological Lett 4, 10 (2018). https://doi.org/10.1186/s40851-018-0093-z

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s40851-018-0093-z

Keywords