Microstructural differences in the osteochondral unit of terrestrial and aquatic mammals
Abstract
During evolution, animals have returned from land to water, adapting with morphological modifications to life in an aquatic environment. We compared the osteochondral units of the humeral head of marine and terrestrial mammals across species spanning a wide range of body weights, focusing on microstructural organization and biomechanical performance. Aquatic mammals feature cartilage with essentially random collagen fiber configuration, lacking the depth-dependent, arcade-like organization characteristic of terrestrial mammalian species. They have a less stiff articular cartilage at equilibrium with a significantly lower peak modulus, and at the osteochondral interface do not have a calcified cartilage layer, displaying only a thin, highly porous subchondral bone plate. This totally different constitution of the osteochondral unit in aquatic mammals reflects that accommodation of loading is the primordial function of the osteochondral unit. Recognizing the crucial importance of the microarchitecture-function relationship is pivotal for understanding articular biology and, hence, for the development of durable functional regenerative approaches for treatment of joint damage, which are thus far lacking.
Editor's evaluation
It is important to determine microstructure-function relationships among different animal species, such as terrestrial and aquatic mammals, since it will help us understand articular biology and inform disease treatment. The authors compared the microstructure and biomechanical property of the osteochondral bone of the humeral head and found the cartilage of aquatic animals have a less stiff cartilage, a more random alignment of collagen fibers, and a lack of a calcified cartilage layer at the cartilage-bone interface. The specific composition of the osteochondral bone in aquatic mammals also reflects the changes in loading.
https://doi.org/10.7554/eLife.80936.sa0Introduction
In the course of evolution, mammals have gone from land to water at several occasions, resulting in separate lineages ranging from Cetacea, which include whales and dolphins, to diverging lines, such as the polar bear (Ursus maritimus) and sea otter (Enhydra lutris) (Uhen, 2007). Polar bears and sea otters have maintained their terrestrial form despite spending a significant portion of their life in water, while other species, in their evolutionary process, have become permanent residents of the aquatic environment, undergoing considerable morphological adaptations to marine life (Uhen, 2007). This long adaptation process that started in the Eocene, is particularly evident in Cetacea (Uhen, 2007; Reidenberg, 2007), where thoracic limbs became flippers, hind limbs regressed and disappeared, and vertebral bones adopted a more uniform morphology (Cozzi et al., 2010; Thewissen et al., 2006).
The main drive for these changes was obviously meeting the physical demands of life in an aquatic environment, where locomotion is largely disentangled from weight-bearing (Turnbull and Cowan, 1999), and thus radically different than as experienced on land. In fact, in water, buoyancy changes the influence of gravity dramatically and water is denser and more viscous than air, influencing resistance (Dejours, 1987). Therefore, locomotion in water happens when the mechanism that propels the animal forward is powerful enough to counter the effects of drag. Different mammalian species specialized in the form of diverse morphologies of their body and limbs (or fins), however, always with the common denominator of the development of propulsive surfaces that oscillate to generate thrust. Swimming of aquatic mammals differs significantly in mechanical efficiency and stress distribution from locomotion in terrestrial mammals relying on limb movement. In particular, pinnipeds, cetaceans and sirenians descend from mammals that evolved to move limbs in a vertical plane, and consequently use a hydrodynamic lift-based momentum exchange to move through water (Farnum, 2007).
In all mammals, regardless whether they live on land or in water, locomotion is effectuated by the musculoskeletal system. Its framework is composed by the skeleton, whose bony components articulate via joints. The function of synovial joints in particular, is to ensure the correct articulation of adjacent bones. In addition, it minimizes friction between the bony components by providing a smooth articulating surface and lubrication to accommodate the biomechanical forces generated by locomotion, transmitting and dampening them (Sophia Fox et al., 2009). The fundamental element that enables physiological function of the synovial joint is the osteochondral unit, composed of articular cartilage and bone. In the front extremities, synovial joints are the principal joint type in terrestrial mammals; in aquatic mammals they have been preserved during evolution at the articulating ends of the humerus, radius, and ulna (Sanchez and Berta, 2010), but have been replaced by fibrous joints in the distal joints of flippers, where bony elements are connected by dense connective tissue (Thewissen, 2009).
Earlier comparative studies on the musculoskeletal system of terrestrial and aquatic mammals have largely focused on bone differences. This is probably because bone is a tissue known to respond extensively and quickly to changes in loading, as described by Wolff’s law and Frost’s mechanostat theory (Wolff, 1893; Frost, 2001). Studies on the macroscopic structures have shown that the bones of the cetacean arm and forearm adopted an hourglass-like form, and that the medullary cavity disappeared to strengthen flipper resistance for steering during swimming (de Buffrénil and Schoevaert, 1988; Felts and Spurrell, 1965). At tissue level, evolutionary adaptations of the bone have been described for both terrestrial (Doube et al., 2010; Mancini et al., 2019) and aquatic mammals (Gray et al., 2007), with a particular focus on how aquatic mammals manage buoyancy by variations in their structural bone density (Gray et al., 2007; Taylor, 2000). However, besides a few studies on some joint diseases (Nganvongpanit et al., 2017; Turnbull and Cowan, 1999), there is very limited knowledge of specific adaptations of the architecture of the cartilage extracellular matrix in aquatic mammals. Moreover, subchondral bone microstructure and related mechanical characteristics, which provide insight in essential aspects of the functionality of the osteochondral unit in joints, have not yet been studied in detail.
The osteochondral unit is the pivotal element of the joint, with a fundamental clinical relevance in highly prevalent joint disorders as osteoarthritis (OA), for which no effective, durable treatment exists yet (Lepage et al., 2019). In terrestrial mammals, comparative studies on the osteochondral unit have shown that the biochemical components of the tissues are strongly preserved across a wide range of species (Malda et al., 2013). A negative allometric relationship exists between cartilage thickness and body mass (Malda et al., 2013) and it is the trabecular bone structure that adapts to increasing body mass (Mancini et al., 2019). At present, comparative studies of the osteochondral unit of terrestrial and aquatic mammalian species do not exist, despite the high clinical relevance of the structure and the potential fundamental insights such studies would provide.
The aim of this study was, therefore, to investigate if and how the structure of both the cartilage and the bone component of the osteochondral unit differs between mammals living on land or in water. Given that the biological function of the osteochondral unit is largely structural in nature and determined by the composition and the architecture of the extracellular matrix (ECM) of both the cartilage and the bone part, the present work focused on analyzing the microstructural composition and architectural features of the osteochondral units from the humeral head of six aquatic and nine terrestrial mammalian species.
Results
Histological and polarized light microscopy analysis of cartilage
Articular cartilage thickness varied widely between species, ranging from 209 μm in the rat to 7660 μm in the sperm whale (Physeter macrocephalus). Total cartilage thickness correlated with body mass (BM) with a negative allometric relationship (R2=0.91, a=0.28, Figure 1A). This relationship was maintained when assessing separately the terrestrial (R2=0.80, a=0.26, green, Figure 1B) and aquatic mammals (R2=0.96, a=0.26, blue, Figure 1B). The relative thickness of the superficial, middle and deep layers of the cartilage (expressed as percentage of total thickness) in terrestrial and aquatic mammals was compared. The analysis showed significant differences in relative layer thickness, with the deep layer being relatively thicker in terrestrial than in aquatic mammals, and the middle layer being relatively thicker in aquatic than in terrestrial mammals. No significant differences were found in superficial layer relative thickness (Figure 1C).
Polarized Light Microscopy (PLM) enabled the visualization of the different collagen fiber orientations (Figure 2A–C), revealing that in terrestrial mammals, the fibers are predominantly oriented perpendicularly to the surface in the deep layer and parallel to the surface in the superficial layer, whereas for the aquatic mammals this arcade orientation of fibers could not be distinguished (Figure 2D–E). However, some alignment in the superficial layer was also observed (Figure 2B), albeit not to the same extend as in the samples of terrestrial mammals.
Biomechanical analysis
Characterization of biomechanical behavior showed a marked difference in the stress relaxation curve of aquatic and terrestrial cartilage samples (Figure 3A). Analysis of time dependent mechanical properties showed that peak modulus was significantly higher in terrestrial mammals than in aquatic mammals (Figure 3C, p<0.05), while at equilibrium, cartilage appeared to be only slightly stiffer in terrestrial mammals (Figure 3D). Initial stress relaxation appeared slower for aquatic mammals and steeper for terrestrial mammals (Figure 3A, Figure 3—figure supplement 1). No significant differences were observed in the phase of slow relaxation (Figure 3—figure supplement 1). In addition, comparison of the relationship of peak and equilibrium modulus and body mass showed, with this sample size, no significant differences (Figure 3E and F).
Histological and micro-computed tomography analysis of the interface
Imaging through histology and micro-CT showed that in aquatic mammals (Figure 4A and B) the interface between the hyaline articular cartilage and the subchondral bone was lacking a layer of calcified cartilage and a dense subchondral plate, as is typically observed at the cartilage-bone interface in terrestrial mammals (Figure 4C; data of additional species provided in Figure 4—figure supplement 1, and additional images cartilage bone interface of the sperm whale provided in Figure 4—figure supplement 2). In fact, the cartilage-bone transition in aquatic mammals is a direct transition from hyaline cartilage to trabecular bone, without intermediate forms of both articular cartilage and subchondral bone, as is the case in terrestrial animals (Figure 4D and E). Bone volume over total volume (BV/TV) in the subchondral area (immediately underneath the cartilage-bone interface) averaged at 80.4 ± 16.8% for aquatic mammals, which was significantly lower than in terrestrial mammals (98.5 ± 0.8%, p<0.01, Figure 4F).
The trabecular bone below the subchondral plate under the central loading area of the humeri was denser, with thicker trabeculae in terrestrial mammals (34.5 ± 11.1% for terrestrial and 29.4 ± 10.3% for aquatic mammals) (Figure 5A, B, C, D, E and F). Trabecular thickness averaged at 211.7±141.1 μm in terrestrial mammals and at 120.0±27.9 μm in aquatic mammals. Trabecular BV/TV increased with size in terrestrial mammals (R2=0.44, Figure 5G, green) and decreased with size in aquatic mammals (R2=0.88, Figure 5G, blue). Thickness of the bone trabeculae (Tb.Th) increased with size in terrestrial mammals (Figure 5H, R2=0.65, a=0.12), but was independent of size in aquatic mammals.
Discussion
In the present study, it was investigated how osteochondral tissue has evolved differently in terrestrial and aquatic mammals to accommodate the different demands of life and locomotion on land or in water. Here, the characteristics of articular cartilage of aquatic mammals are described for the first time and they appeared to be widely different from all known terrestrial species in which this has been investigated. In contrast to the arcade-like (Benninghoff, 1925) collagen arrangement in the articular cartilage of terrestrial mammals, articular cartilage collagen fibers did not display this alignment in aquatic mammals. Moreover, the cartilage-bone interface of aquatic mammals entirely lacks the typical calcified layer present in all terrestrial species in which the osteochondral unit has been investigated. Further, it also lacks a dense subchondral plate, featuring an abrupt transition from hyaline cartilage to subchondral trabecular bone. Aquatic mammals lack the intermediate transitional forms of either cartilage (the calcified layer) or bone (the subchondral plate) that are typical of terrestrial species. Interestingly, for semi-aquatic species, such as the harbor seal (Phoca vitulina), there appears still to be some densification of the subchondral bone region, but this needs further investigation. Nevertheless, the trabeculae of the subchondral bone are thinner in aquatics than in terrestrials and they do not increase in size with increasing body mass, in contrast with what has been observed in terrestrial mammals.
Previously, it was shown that cartilage biochemistry is remarkably preserved among species, and that cartilage tissue displays microstructural features that allow functionality in mice as much as in elephants (Malda et al., 2013). One of the important features in the accommodation of different forces generated by body mass is the thickness of the tissue. Our group previously reported that articular cartilage thickness scaled with negative allometry in relationship to body mass (a=0.28), based on the analysis of nearly 60 terrestrial mammalian species (Malda et al., 2013). Interestingly, in the analysis presented here, it was observed that this correlation was similar in aquatic mammals, suggesting that this relationship might be dependent of factors other than loading, such as possibly metabolic limitations (Holliday et al., 2010). In hyaline cartilage, an avascular tissue in which nutrient supply is dependent on diffusion and fluid flow across a poroviscoelastic matrix (Lu and Mow, 2008), metabolite limitation is known to be related to tissue thickness (Malda et al., 2013; White and Seymour, 2003). In the sperm whale, an almost direct transition from hyaline cartilage to trabecular subchondral bone was observed, making diffusion from the rich vascular network of the trabecular bone to the hyaline cartilage probably as easy as diffusion from the synovial fluid at the articular surface of that hyaline cartilage. The increased porosity, and thus the greater potential for transfer of solutes, may permit the great cartilage thickness seen in some aquatic mammals (over 7 mm in the largest whale [Balaenoptera physalus]). As an avascular tissue, cartilage is dependent on diffusion, so critical thickness will be larger if diffusion of nutrients can be provided more efficiently from two sides. The ratio of nutrient exchange is likely to be different in aquatic mammals, as at the few sites where the hyaline cartilage is in direct contact with the subchondral bone in terrestrial animals, solute exchange has been shown to be fivefold higher than through the calcified cartilage (Arkill and Winlove, 2008).
Simon et al. did not observe a consistent correlation between cartilage thickness and estimated compressive stress on the joint in a study conducted on five species of terrestrial mammals that varied widely in body mass (mouse, rat, dog, sheep, and cow; Simon, 1970). When studying the thicknesses of the different zones that make up articular cartilage, significant differences were observed between aquatic and terrestrial mammals in the deep and middle layers. However, this was not the case for the superficial layer. This supports the notion that the primary role of the superficial layer of cartilage is to protect the deeper layers from shear forces and to distribute impact forces and loads among the tissue (Julkunen et al., 2007; Bevill et al., 2010). Importantly, the load distribution of forces is essential in both aquatic and terrestrial species. Unlike shear forces, compressive loads will, however, be substantially lower in aquatic species, which would be a likely explanation for the relatively thinner deep layer.
A significant difference in the predominance of fiber orientation was found between aquatic and terrestrial mammals in all cartilage layers: collagen fiber orientation in terrestrial species followed the classic Benninghoff arcade model (Benninghoff, 1925). This is known to be conserved in all terrestrial species in which the configuration has been investigated (Kääb et al., 1998) and in which orientation is principally parallel to the surface in the superficial layer, and is changing in the middle zone as the fibers are bending towards an orientation perpendicular to the subchondral bone in the deep layer. The deep layer is known to increase in relative thickness with increasing body mass and therefore increasing compressive loading (Felts and Spurrell, 1965). Conversely, in aquatic mammals, orientation of the collagen fibers throughout the whole thickness of the cartilage was more random. In an aquatic mammal, the total biomechanical loading of the joints will be less because of the absence of gravity and in a relative sense, shear forces, which perhaps are less than in case of terrestrial locomotion, but are still generated by swimming (Dejours, 1987; Williams and Worthy, 2009; Williams, 1999; Williams, 2001), will become more important. This could potentially explain the presence of only a thin superficial layer. Moreover, it is likely that the heterogeneity in the alignment of the collagen fibers in aquatic mammals represents the strongly reduced need for a specifically organized orientation as seen in terrestrial mammals.
The different collagen fiber orientation observed in the aquatic mammals compared to terrestrial mammals is also reflected in the biomechanical behavior of cartilage. To analyze this effect, both peak and equilibrium modulus of cartilage were determined. These parameters are commonly used to describe the mechanical functionality of native and engineered cartilage because they provide a better understanding of how the constituents of cartilage contribute to its overall mechanical function (Mow et al., 1984; Klika et al., 2016). The mechanical responses observed in the aquatic mammals are indicative of an overall less stiff cartilage tissue, compared to that in terrestrial mammals. Furthermore, the absence of a defined collagen network structure in the aquatic mammals is likely to provide less resistance to osmotic swelling, as in terrestrial mammal’s cartilage up to two thirds of the elastic modulus in compression arises from the electrostatic contributions and osmotic swelling that is restricted inside the organized collagen matrix (Lai et al., 1991; Canal Guterl et al., 2010).
The lack of a calcified cartilage layer and a compact subchondral plate, which are typical features of the osteochondral unit in terrestrial animals (Sophia Fox et al., 2009), was highly remarkable. The main function of the calcified layer in the osteochondral unit is to provide a gradual transition in stiffness from hyaline cartilage to subchondral bone and to serve as anchor point for the collagen fibers in the deep layer of the hyaline cartilage, that is to serve as a foundation for the Benninghoff arcades (Sophia Fox et al., 2009; Benninghoff, 1925; Madry et al., 2010). The subchondral plate has a role in distributing the compressive forces acting on the joint. As discussed earlier, in aquatic mammals compressive loading is hardly present and it may well be that firm anchoring of the collagen fibrils is much less important for this reason. Detailed studies using high resolution electron microscopy might shed some light on this issue.
Bone density deeper under the cartilage tissue (trabecular BV/TV) was relatively higher in terrestrial mammals and increased with body mass, while the opposite was observed in aquatic mammals. For the terrestrial animals the relationship, this is in accordance with the earlier observation in these animals trabecular bone increases with increasing in size (Barak et al., 2013). The second observation, can be attributed to the lack of loading and the fact that in aquatic mammals living in deep waters, in which the absence of gravity obviates the need for such adaptation to increasing body mass, it has been shown that low bone density enables dynamic buoyancy control (Gray et al., 2007), explaining the opposite tendency in terrestrials. According to Dumont et al., bone microanatomy carries an ecological signal that holds information about habitat and locomotion patterns, which is particularly evident in tetrapods that have developed advanced secondary adaptations to life in water (Dumont et al., 2013). In the matter of buoyancy control of aquatic mammals, species recently adapted to aquatic life or living in shallow waters, such as sea cows (sirenians), have bones that are denser than in terrestrial forms to allow static control (ballast), as opposed to a lower bone density typically associated with the dynamic buoyancy control required by animals that exhibit deep diving behavior for longer periods (Gray et al., 2007; Wall, 1983).
Moreover, bone trabecular thickness behaved in a similar manner. This parameter also increased with size in terrestrial mammals, with a trend similar to previous reports in literature (Mancini et al., 2019; Doube et al., 2011), while in aquatic mammals trabecular thickness appeared to scale independently. This independence can possibly be explained by the lack of compressive forces necessary to elicit a proportionate response in the bone for trabecular thickness as suggested by Wolff’s law (Wolff, 1893; Frost, 2001). In fact, Rolvien et al. investigated characteristics of the vertebral bodies of three toothed whale species sperm whale, orca (Orcinus orca) and harbor porpoise (Phocoena phocoena) and found that bone volume fraction (BV/TV) did not scale with body mass, but that trabeculae were thicker and fewer in number in larger whale species (Rolvien et al., 2017). Although research has shown the relationship between bone microanatomical organization, habitats and locomotion patterns in different taxa (Dumont et al., 2013), this was less clear when comparing different species of whales, particularly in view of the variability of for example diving-depth behavior of the different animals (Rolvien et al., 2017).
Interestingly, the histological sections of the aquatic mammals, displayed small islands of cartilaginous tissue embedded within the subchondral bone, and in the larger aquatic species, such as the sperm whale, large islands of cartilage tissue were found deep in the trabecular bone. These appeared as large lacunae upon first inspection with micro-CT, and an initial explanation was that these were a possible areas of osteocyte death due to barotraumas, which have been anecdotally described in deep-diving mammals (Moore and Early, 2004). However, upon histological evaluation, these islands appeared to be composed of healthy cartilage-like tissue naturally embedded in the bone, devoid of evident signs of bone remodeling, not suggesting any form of pathology. The areas may, therefore, represent incomplete osteochondral ossification, well-known in the horse and other species as osteochondrosis (Thomas and Barnes, 2015; McCoy et al., 2013) similarly to what was found and described by Rolvien et al. in the central areas of the vertebrae of sperm whales (Rolvien et al., 2017).
The number of individuals and species from which samples could be obtained was based on species availability and practical (such as the possibility of tissue harvest shortly after death from aquatic mammals) and does not compare to the numbers of earlier multispecies investigations in terrestrial animals (Mancini et al., 2019; Malda et al., 2013), which could be considered as a major limitation of this study. However, the included aquatic species covered a wide range of body masses (from 54 to 48,000 kg) and can, given the high consistency of the findings and the clear and unambiguous landslide differences with terrestrial species, be considered representative, especially in view of the very limited and subtle differences observed between a large number of different terrestrial species in the earlier studies (Mancini et al., 2019; Malda et al., 2013). Taken together, the findings of the current study add significantly to our knowledge of the basic aspects of the biology of the osteochondral unit and clearly reflect the crucial influence of loading. We know from evolution that life started to develop in the water and that later species came on land to develop further (Colbert and Morales, 2001). We also know that all mammals have initially evolved on land and that some of them returned to the water to become what we now call aquatic mammals (Colbert and Morales, 2001). Where we do not have histological data from the joints of the land-borne species from which our current aquatic mammalian species have evolved, it is reasonable to assume that their osteochondral units were similar to the common format of all currently existing terrestrial mammalian species. The huge structural differences between the osteochondral units of the currently living aquatic and terrestrial mammals make hence clear that accommodating loading is the most important evolutionary pressure on the osteochondral unit of terrestrial mammals. Further, it highlights that this challenge is dealt with by the architecture and structural composition of the osteochondral unit and not by molecular composition or cellular diversity. The last two elements of the osteochondral unit are similar in aquatics and terrestrials. It is the specific tissue architecture and structural organization of especially the cartilage-bone interface that got lost in aquatic mammals, once the loading-based evolutionary pressure disappeared. This increased insight in the fundamental biology of the structural components of the osteochondral unit, obtained by observing Nature, is of great importance from a tissue repair perspective, as it is a natural guide for which considerations should be prioritized in the quest for recapitulation of the native template for regenerative purposes. In this sense, the study presents strong support for a recent call to rethink the current paradigms in articular cartilage regeneration (Malda et al., 2019), through focusing on solutions with a proper architecture that allow restoration of (biomechanical) function, rather than on mere cellular activity in terms of production of extracellular matrix components.
Methods
Collection of materials and tissues
To investigate the morphology of cartilage tissue and of the interface between cartilage and bone, histological analysis was performed on osteochondral samples taken from fresh humeral heads (Figure 6) of a variety of terrestrial and aquatic mammals (Table 1). Osteochondral tissue samples were harvested postmortem from the weight bearing central area of the humeral head of adult animals sent for necropsy to the division of Pathology, Faculty of Veterinary Medicine, Utrecht University, in the Netherlands. For aquatic large species, samples were harvested on site directly from animals that were found dead on the coasts of the Netherlands.
Animal species, age and body mass were recorded (or estimated in the case of whales based on species and body size), and macroscopic pictures of the joints were taken. Joints demonstrating macroscopic or microscopic signs of cartilage degeneration were excluded; animals displaying signs of incomplete endochondral ossification were identified as skeletally immature and excluded as well.
In total 84 tissue samples (34 for histology and polarized light microscopy, 30 for micro-CT and 20 for biomechanical testing) were harvested from the humeral heads of mammals belonging to 15 different species, 9 terrestrial and 6 aquatic (Table 2); samples for histology were fixed in formalin 4%, while samples for micro-CT analysis were stored in 70% ethanol; samples for biomechanical testing were frozen immediately after harvest in Tissue-tek (Sakura Finitek, USA) until use.
Histological preparation and analysis
Samples were fixed using 4% formalin, decalcified with Formical-2000 (EDTA/formic acid; Decal Chemical Corporation, Tallman, NY), dehydrated, cleared in xylene, embedded in paraffin and cut with a microtome to yield 5 μm sections. Sections were stained with fast green and Safranin-O for measurements of cartilage thickness (distance from the surface to the interface with the subchondral bone). Average thickness of cartilage tissue for each sample was determined by averaging 3 measurements per image taken from different locations of the section. Digital images were analyzed using ImageJ software (Schindelin et al., 2012).
Polarized light microscopy
To evaluate the orientation of the cartilaginous collagen network, Polarized Light Microscopy (PLM) was chosen for its capacity to visualize the orientation of anisotropic materials (Changoor et al., 2011; Yarker et al., 1983; Speer and Dahners, 1979). Histological sections were stained with picrosirius red to stain the collagen fibers of the cartilage extracellular matrix (Junqueira et al., 1979). All microscopic images were acquired through an Olympus DP73 digital camera using Cell^F software (Matrix Optics, Malaysia). A combination of U-ANT and U-POT filters mounted on a light microscope (Olympus BX51, Olympus) was used for the PLM measurements. The two cross polarizers were used so that highly ordered collagen fibers that were perpendicular or tangential to the articular surface appeared bright or red, while fibers with other orientations (non-birefringent) appeared darkest. Collagen fiber orientation visualized via PLM was used for the classification of the superficial, middle and deep layer of cartilage (Changoor et al., 2011; Rieppo et al., 2008). The area with the minimum birefringence value was identified as the border between the superficial and the middle zones, whereas the deep zone was considered to begin when the orientation angle values reached a plateau (typically close to 90 degrees with respect to cartilage surface) (Julkunen et al., 2007; Arokoski et al., 1996). Micrographs acquired via PLM were converted to 8-bit images and the angle of orientation of the collagen structures was calculated within each of the three zones via automated image analysis and Fourier spectrum analysis, using the Directionality plug-in from the FIJI software (Schindelin et al., 2012). Data were plotted as histograms of frequency for each representative angle, and as radial plots (0° represent fibers parallel to the articulating surface) (Figure 7). For each of the histograms, the peaks centered at 0° and ±90° and were integrated to calculate their area using the OriginPro 8 software package (OriginLab, USA). Finally, cartilage relative layer thicknesses (%) was calculated based on the data obtained with PLM.
Biomechanical testing and data analysis
The mechanical properties of the native cartilage of the osteochondral plugs of the different species were assessed by uniaxial unconfined compression using a universal testing machine (Zwick Z010, Germany) equipped with a 20 N load cell. Stress relaxation tests were performed by first applying a pre-load of 0.01 N to test samples and then strained to 15% at a rate of 10 µm/s, followed by a relaxation period of 900 s. Tests were conducted on cylindrical osteochondral plugs (Table 2). Prior to testing, samples were fixed to the bottom of a custom-made polycarbonate container using a cyanoacrylate-based adhesive and all tests were performed in PBS to approximate physiological conditions.
Peak and equilibrium stresses were calculated from each engineered stress-strain relaxation curve at peak or equilibrium, respectively. Here, stress is defined as the force divided by the specimen’s unloaded cross-sectional area, and strain as the ratio between the original sample thickness and the displacement of compression platen. To quantify the relaxation response, a piecewise exponential function was fitted to the obtained stress-time curves, following the method by Castilho et al., 2019. Briefly, the fitted curves consist of three segments corresponding to an initial loading phase, and a fast and a slow relaxation phase. The formulae for the exponential curves for the fast and slow relaxation phases used are as follows,
where the coefficients A and B are statistical parameters quantifying the shape of the relaxation curves and τ is a time constant that determines the rate of stress relaxation. The fitting procedure used minimizes the root mean square error between these curves and the measured stress, subject to the constraints that the fitted curves match the peak and equilibrium stress and that the two curves match at t=100 s.
Micro-Computed Tomography
For the detailed analysis of the microstructural features of the subchondral and trabecular bone, Micro-Computed Tomography (micro-CT) was selected for the accurate measurement of micron-sized structures that constitute the bony tissue (Fajardo and Müller, 2001; Holdsworth and Thornton, 2002). Micro-CT images of the osteochondral cores were obtained with a micro-CT scanner (Quantum FX, Perkin Elmer, USA, voxel size = 20 μm3). The automatically reconstructed micro-CT images were subsequently converted to series of 2D TIFF images and were binarized using local thresholding (Bernsen technique). BoneJ software (Doube et al., 2010) was used to determine the trabecular thickness (Tb.Th) and the bone volume fraction (Tb BV/TV) of the bone immediately underneath the cartilage (subchondral bone; determined as previously described Mancini et al., 2019), and at the center of the bone of the osteochondral cores (trabecular bone).
Statistical analysis
Comparisons between multiple groups were performed using a one-way ANOVA combined with post-hoc t-tests with Bonferroni correction. A one-way ANOVA with Tukey’s post hoc test was used for comparison in biomechanical analysis. A regression analysis using a power curve fit was performed for correlations between body mass and the parameters of cartilage thickness, peak modulus, trabecular thickness, and trabecular BV/TV. When applied, normality and homogeneity were checked with Shapiro-Wilks and Levene’s tests. The level of statistical significance was set at p<0.05. Statistical analysis was performed using GraphPad Prism 8 (GraphPad Software, Dotmatics, USA).
Data availability
Data generated or analysed during this study are included in the manuscript and supporting file; Source Data files have been provided for Figure 1.
References
-
Solute transport in the deep and calcified zones of articular cartilageOsteoarthritis and Cartilage 16:708–714.https://doi.org/10.1016/j.joca.2007.10.001
-
Of mice, rats and men: trabecular bone architecture in mammals scales to body mass with negative allometryJournal of Structural Biology 183:123–131.https://doi.org/10.1016/j.jsb.2013.04.009
-
Form und bau der gelenkknorpel in ihren beziehungen zur funktionZeitschrift Für Anatomie Und Entwicklungsgeschichte 76:43–63.https://doi.org/10.1007/BF02134417
-
A polarized light microscopy method for accurate and reliable grading of collagen organization in cartilage repairOsteoarthritis and Cartilage 19:126–135.https://doi.org/10.1016/j.joca.2010.10.010
-
BookColbert’s Evolution of the Vertebrates: A History of the Backboned Animals through TimeJohn Wiley & Sons.
-
Diving adaptations of the cetacean skeletonOpen Zoology Journal 2:24–32.https://doi.org/10.2174/1874336600902010024
-
ConferenceWater and air physical characteristics and their physiological consequencesComparative physiology: Life in water and on land 9.
-
Trabecular bone scales allometrically in mammals and birdsProceedings of the Royal Society B 278:3067–3073.https://doi.org/10.1098/rspb.2011.0069
-
Three-dimensional analysis of nonhuman primate trabecular architecture using micro-computed tomographyAmerican Journal of Physical Anthropology 115:327–336.https://doi.org/10.1002/ajpa.1089
-
BookPostnatal Growth of Fins and Limbs through Endochondral OssificationFins into limbs: evolution, development, and transformation.https://doi.org/10.7208/9780226313405-011
-
Structural orientation and density in cetacean humeriThe American Journal of Anatomy 116:171–203.https://doi.org/10.1002/aja.1001160109
-
From Wolff’s law to the Utah paradigm: insights about bone physiology and its clinical applicationsThe Anatomical Record 262:398–419.https://doi.org/10.1002/ar.1049
-
Sink or swim? Bone density as a mechanism for buoyancy control in early cetaceansThe Anatomical Record 290:638–653.https://doi.org/10.1002/ar.20533
-
Micro-CT in small animal and specimen imagingTrends in Biotechnology 20:S34–S39.https://doi.org/10.1016/S0167-7799(02)02004-8
-
Picrosirius staining plus polarization microscopy, a specific method for collagen detection in tissue sectionsThe Histochemical Journal 11:447–455.https://doi.org/10.1007/BF01002772
-
An overview of multiphase cartilage mechanical modelling and its role in understanding function and pathologyJournal of the Mechanical Behavior of Biomedical Materials 62:139–157.https://doi.org/10.1016/j.jmbbm.2016.04.032
-
A triphasic theory for the swelling and deformation behaviors of articular cartilageJournal of Biomechanical Engineering 113:245–258.https://doi.org/10.1115/1.2894880
-
Biomechanics of articular cartilage and determination of material propertiesMedicine and Science in Sports and Exercise 40:193–199.https://doi.org/10.1249/mss.0b013e31815cb1fc
-
The basic science of the subchondral boneKnee Surgery, Sports Traumatology, Arthroscopy 18:419–433.https://doi.org/10.1007/s00167-010-1054-z
-
Rethinking articular cartilage regeneration based on a 250-year-old statementNature Reviews. Rheumatology 15:571–572.https://doi.org/10.1038/s41584-019-0278-7
-
Articular osteochondrosis: a comparison of naturally-occurring human and animal diseaseOsteoarthritis and Cartilage 21:1638–1647.https://doi.org/10.1016/j.joca.2013.08.011
-
Fluid transport and mechanical properties of articular cartilage: a reviewJournal of Biomechanics 17:377–394.https://doi.org/10.1016/0021-9290(84)90031-9
-
Anatomical adaptations of aquatic mammals: the anatomical record: advances in integrative anatomy and evolutionary biologyAdvances in Integrative Anatomy and Evolutionary Biology 290:507–513.https://doi.org/10.1002/ar.20541
-
Practical considerations in the use of polarized light microscopy in the analysis of the collagen network in articular cartilageMicroscopy Research and Technique 71:279–287.https://doi.org/10.1002/jemt.20551
-
Comparative anatomy and evolution of the odontocete forelimbMarine Mammal Science 26:140–160.https://doi.org/10.1111/j.1748-7692.2009.00311.x
-
Fiji: an open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
Scale effects in animal joints i articular cartilage thickness and compressive stressArthritis and Rheumatism 13:244–256.https://doi.org/10.1002/art.1780130305
-
The collagenous architecture of articular cartilage: correlation of scanning electron microscopy and polarized light microscopy observationsClinical Orthopaedics and Related Research pp. 267–275.
-
The bone joint osteochondrosis in extant and fossil marine mammalsContributions in Science 523:1–35.https://doi.org/10.5962/p.226778
-
Synovial joint disease in wild cetaceansJournal of Wildlife Diseases 35:511–518.https://doi.org/10.7589/0090-3558-35.3.511
-
Evolution of marine mammals: back to the sea after 300 million yearsAnatomical Record 290:514–522.https://doi.org/10.1002/ar.20545
-
The correlation between high limb-bone density and aquatic habits in recent mammalsJournal of Paleontology 57:197–207.
-
The evolution of cost efficient swimming in marine mammals: limits to energetic optimizationPhilosophical Transactions of the Royal Society of London. Series B 354:193–201.https://doi.org/10.1098/rstb.1999.0371
-
Marine Mammal Biology: An Evolutionary ApproachAnatomy and physiology: the challenge of aquatic living, Marine Mammal Biology: An Evolutionary Approach, John Wiley & Sons.
-
Das gesetz der transformation der knochenDMW - Deutsche Medizinische Wochenschrift 19:1222–1224.https://doi.org/10.1055/s-0028-1144106
-
Birefringence of articular cartilage and the distribution on collagen fibril orientationsConnective Tissue Research 11:207–213.https://doi.org/10.3109/03008208309004857
Article and author information
Author details
Funding
European Commission (3099622 (FP7))
- Irina AD Mancini
- P René van Weeren
- Jos Malda
Dutch Arthritis Society (LLP12)
- Riccardo Levato
- Miguel Dias Castilho
- P René van Weeren
- Jos Malda
Dutch Arthritis Society (LLP22)
- Riccardo Levato
- Miguel Dias Castilho
- P René van Weeren
- Jos Malda
Dutch Ministry of Nature, Agriculture and Food Quality (140000353)
- Lonneke L IJsseldijk
Dutch Ministry of Nature, Agriculture and Food Quality (WOT04-009-045)
- Lonneke L IJsseldijk
Dutch Ministry of Nature, Agriculture and Food Quality (HD3611/BO11018.02 065)
- Lonneke L IJsseldijk
European Community’s Seventh Framework Programme (FP7/2007– 2013 (HydroZONES))
- P René van Weeren
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors acknowledge the staff of the division of Pathology, from the Department of Biomolecular Health Sciences, of the faculty of Veterinary Medicine, in particular R Wagensveld-van den Dikkenberg, for their valuable help with the sample collection. The authors would also like to thank P Kamminga and the Naturalis Biodiversity Center in Leiden, the Netherlands, for access to their collection of mammalian skeletons. Post-mortem investigations of cetaceans in the Netherlands are commissioned by the Dutch Ministry of Nature, Agriculture and Food Quality (project grant numbers: 140000353; WOT04-009-045, HD3611/BO11018.02 065). The research leading to these results has received funding from the European Community’s Seventh Framework Programme (FP7/2007– 2013) under grant agreement 309962 (HydroZONES), and the Dutch Arthritis Association (ReumaNL, LLP-12 and LLP-22).
Copyright
© 2023, Mancini et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 645
- views
-
- 154
- downloads
-
- 0
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Developmental Biology
The morphogen FGF8 establishes graded positional cues imparting regional cellular responses via modulation of early target genes. The roles of FGF signaling and its effector genes remain poorly characterized in human experimental models mimicking early fetal telencephalic development. We used hiPSC-derived cerebral organoids as an in vitro platform to investigate the effect of FGF8 signaling on neural identity and differentiation. We found that FGF8 treatment increases cellular heterogeneity, leading to distinct telencephalic and mesencephalic-like domains that co-develop in multi-regional organoids. Within telencephalic regions, FGF8 affects the anteroposterior and dorsoventral identity of neural progenitors and the balance between GABAergic and glutamatergic neurons, thus impacting spontaneous neuronal network activity. Moreover, FGF8 efficiently modulates key regulators responsible for several human neurodevelopmental disorders. Overall, our results show that FGF8 signaling is directly involved in both regional patterning and cellular diversity in human cerebral organoids and in modulating genes associated with normal and pathological neural development.
-
- Developmental Biology
Wnt signaling plays crucial roles in embryonic patterning including the regulation of convergent extension (CE) during gastrulation, the establishment of the dorsal axis, and later, craniofacial morphogenesis. Further, Wnt signaling is a crucial regulator of craniofacial morphogenesis. The adapter proteins Dact1 and Dact2 modulate the Wnt signaling pathway through binding to Disheveled. However, the distinct relative functions of Dact1 and Dact2 during embryogenesis remain unclear. We found that dact1 and dact2 genes have dynamic spatiotemporal expression domains that are reciprocal to one another suggesting distinct functions during zebrafish embryogenesis. Both dact1 and dact2 contribute to axis extension, with compound mutants exhibiting a similar CE defect and craniofacial phenotype to the wnt11f2 mutant. Utilizing single-cell RNAseq and an established noncanonical Wnt pathway mutant with a shortened axis (gpc4), we identified dact1/2-specific roles during early development. Comparative whole transcriptome analysis between wildtype and gpc4 and wildtype and dact1/2 compound mutants revealed a novel role for dact1/2 in regulating the mRNA expression of the classical calpain capn8. Overexpression of capn8 phenocopies dact1/2 craniofacial dysmorphology. These results identify a previously unappreciated role of capn8 and calcium-dependent proteolysis during embryogenesis. Taken together, our findings highlight the distinct and overlapping roles of dact1 and dact2 in embryonic craniofacial development, providing new insights into the multifaceted regulation of Wnt signaling.