Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Scientific Reports volume  12, Article number: 11707 (2022 ) Cite this article

The functional relationship between bone and cartilage is modulated by mechanical factors. Scarce data exist on the relationship between bone shape and the spatial distribution of cartilage thickness. The aim of the study was to characterise the coupled variation in knee bone morphology and cartilage thickness distributions in knees with healthy cartilage and investigate this relationship as a function of sex, height, body mass, and age. MR images of 51 knees from young adults (28.4 ± 4.1 years) were obtained from a previous study and used to train a statistical shape model of the femur, tibia, and patella and their cartilages. Five multiple linear regression models were fitted to characterise morphology as a function of sex, height, body mass, and age. A logistic regression classifier was fitted to characterise morphological differences between males and females, and tenfold cross-validation was performed to evaluate the models’ performance. Our results showed that cartilage thickness and its distribution were coupled to bone morphology. The first five shape modes captured over 90% of the variance and described coupled changes to the bone and spatial distribution of cartilage thickness. Mode 1 (size) was correlated to sex (p < 0.001) and height (p < 0.0001). Mode 2 (aspect ratio) was also correlated to sex (p = 0.006) and height (p = 0.017). Mode 4 (condylar depth) was correlated to sex only (p = 0.024). A logistic regression model trained on modes 1, 2, and 4 could classify sex with an accuracy of 92.2% (95% CI [81.1%, 97.8%]). No other modes were influenced by sex, height, body mass, or age. This study demonstrated the coupled relationship between bone and cartilage, showing that cartilage is thicker with increased bone size, diaphysis size, and decreased femoral skew. Our results show that sex and height influence bone shape and the spatial distribution of cartilage thickness in a healthy young adult population, but body mass and age do not.

The functional relationship between bone and cartilage is complex and is modulated by mechanical factors introduced by loading and motion1. Mature cartilage has a location-dependent histomorphology that is developed in response to its specific loading history1. The mechanical stresses and strains experienced by the cartilage over time influence the morphology and material properties of the tissue. Conversely, the morphology of bone and cartilage influences the instantaneous stresses and strains experienced in the joint2. Bone size and shape can affect the contact area and the lines of action of the muscles and supporting ligaments crossing the joint3, while cartilage thickness influences the mechanical properties of cartilage4 and the stresses and strains experienced in the tissue2,5,6. Characterising the relationship between bone and cartilage morphology in vivo is important to understand the knee joint’s functional anatomy and pathology.

Knee bone morphology has been characterised as a function of sex7,8, osteoarthritis (OA)9,10,11,12, and knee ligament injury13. Studies have also quantified cartilage thickness across individuals5,14,15,16 and found that cartilage is generally thicker where cartilage stresses are high6,17. Males appear to have greater volume of knee cartilage compared to females18,19, and cartilage morphology is in part mediated by size and body mass of the individual15,17,18. However, these studies have analysed cartilage and bone morphology in isolation or have simplified cartilage morphology to a singular value (such as cartilage volume or mean thickness), and thus, do not capture the complexity and nuance of the coupled variation in cartilage thickness distribution and bone morphology in three dimensions. This information may have important implications for understanding the onset of pathology such as OA, the development of surrogate and predictive models, and the design of implants.

Statistical shape modelling (SSM) characterises complex 3D morphologies by decomposing shape features into a set of statistically significant principal components that describe the main variations in the population20,21. The shape features can be extended to include scalar field values (such as cartilage thickness, and bone mineral density) to explore shape-field relationships. The principal components, or modes, can then be used to visualize variations and correlate to parameters such as patient demographics to quantitatively examine differences between cohorts10,22,23,24. To the best of our knowledge, statistical shape models have not been used to characterise the coupled morphological variations of adult knee bones and their cartilage thickness distribution in relation to subject demographics.

The purpose of this study was to develop a statistical shape and cartilage thickness field model of the adult knee, consisting of the femur, tibia, patella and their cartilages, to characterise the coupled variation in bone morphology and cartilage thickness as a function of sex, height, body mass, and age, in a cohort of 51 knees with healthy cartilage.

MR images (n = 51) of adult knees (28.4 ± 4.1 years, Table 1) were obtained from an ongoing study on patellofemoral pain16,25,26, where all participants were screened by a radiologist to ensure they had no cartilage damage or degenerative changes. All participants were advised on all aspects of the study and analysis of data before informed consent was obtained. Permission was obtained to analyse the data in this study. Ethics approval for this study (Reference #3346) was granted by the Stanford University Institutional Review Board, and the procedures followed were in accordance with the ethical standards of the responsible committee on human experimentation (institutional and national) and with the Helsinki Declaration of 1975, as revised in 2000. All methods were carried out in accordance with relevant guidelines and regulations. The bones, including the femur, patella, tibia, and their corresponding cartilages, including both the subchondral and articular surfaces, were manually segmented in Stradwin27 to produce triangulated surface meshes and point clouds (Fig. 1A). To correct for the varying field of view (FOV) in the imaging data, MAPClient28 was used to fit an existing statistical shape model (SSM)8, trained on a CT dataset (n = 204) from the Victoria Institute of Forensic Medicine, to the femurs and tibias (Fig. 1A). This fit consisted of an SSM fit followed with a local nodal fit to a root-mean-squared (RMS) error < 1 mm. The femurs were then cropped at a height equal to the epicondylar width and the tibias were cropped at a height equal to the width of the tibial plateau.

Schematic of workflow to generate cartilage thickness field and bone shape model. MR images (n = 51) of the knee were segmented, processed in MAPClient, and resampled (A) to produce bone and cartilage point clouds. Parametric or correspondent bone point clouds were obtained via an iterative fitting process (B). Subchondral bone node numbers on these correspondent point clouds were found by combining node numbers obtained from all subjects using a closest-point algorithm (C). Cartilage thickness maps were calculated (D) by computing the magnitude of the projection of the closest articular cartilage point to the normal vector of each subchondral bone node. Principal component analysis (E) was performed on features consisting of the corresponding nodal coordinates of the bone and the cartilage thickness per subchondral node to produce a statistical model of the cartilage thickness and bone shape.

A published method10,22,24 was adapted to maximise nodal correspondence in the bone point clouds and create a training set of parametric point clouds (Fig. 1B). One set of bone point clouds were selected and resampled in Meshlab29 to generate template point clouds with mean nearest neighbour distance of 0.5 mm and mesh density of ~ 2.6 points per mm2 for the femur, tibia, and patella (Fig. 1A). The template point clouds were iteratively fitted to the training set data via a series of coarse to fine fits using adaptive radial basis function fitting to generate a training set of maximally correspondent points30. Rigid alignment of the training set data to the templates was achieved by minimising the least squared distances of corresponding nodes. Principal component analysis (PCA) was performed on the correspondent nodes to generate an intermediate statistical shape model of the bones31. This process was repeated with the shape model as the template until the correspondent RMS error between the point cloud in the current and previous iteration was less than 0.01 mm.

The cartilage point clouds were resampled (Fig. 1A) to a density of ~ 16 points per mm2, resulting in ~ 200,000 points for the femoral cartilage, ~ 50,000 points for the patellar cartilage, and ~ 50,000 points each for the lateral and medial tibial cartilages. The cartilage point clouds for each subject in the training set were overlaid with their corresponding parametric bone point clouds. A closest point algorithm was used to identify bone nodes that were closest to the subchondral cartilage layer (Fig. 1C). These bone nodes, of all subjects, were combined to produce a list of subchondral bone nodes shared by the training set.

The cartilage thickness was then calculated for each subchondral bone node (Fig. 1D)32. The closest articular cartilage point to the normal of each subchondral bone node was projected onto the normal vector. The magnitude from the subchondral bone node to this projection was calculated to obtain the thickness at that subchondral bone node. This calculation was performed on the femur, tibia, and patella for all subjects in the training set (n = 51) to obtain correspondent maps of cartilage thickness.

Correspondent nodes and cartilage thickness at each node were then analysed using PCA. A 2D matrix was constructed with n-rows (n = 51 independent observations) and m-columns of features, consisting of mean centred nodal coordinates and thickness values. The matrix was normalised by dividing all columns by its standard deviation. PCA was performed on the resulting matrix to produce a statistical shape and field model (SSFM) of bone shape and cartilage thickness (Fig. 1E).

The SSFM scores of the first n modes, such that the cumulative variation captured was over 90%, were used for statistical analysis. Multiple linear regression was used to investigate the influence of sex, height, body mass, and age on the SSFM scores.

A multivariable logistic regression analysis33 was performed on the scores of the SSFM to characterise the coupled morphological differences in bone shape and cartilage thickness distribution between males and females. Principal components were added in a stepwise fashion until no more statistically significant improvement of the fit of the model was observed. Tenfold stratified cross-validation was performed to evaluate the models performance.

The first five principal components (or modes) of the SSFM of bone shape and cartilage thickness captured over 90%, of the variation in the morphology in the training set (Fig. 2).

Variation and cumulative variation represented by principal component number for the statistical bone shape and cartilage thickness model.

Multiple linear regression models were fitted for the first five modes of the SSFM (Table 2). The first mode(size) was strongly influenced by sex (p < 0.001) and height (p < 0.0001). The second mode was also influenced by sex (p = 0.006) and height (p = 0.017). The fourth mode was influenced by sex (p = 0.024). These modes are visualised in Figs. 3 and 4. No other modes were influenced by sex, height, body mass, or age.

Coupled variations in bone shape and cartilage thickness of the SSFM showing principal components that were influenced by sex and height (PC1 and PC2), or sex (PC4).

Coupled variations in bone shape of the SSFM showing principal components that were influenced by sex and height (PC1 and PC2) or sex (PC4).

The regions with the thickest cartilage in the femur were located in the trochlear groove, the femoral condyles, and along the medial condyle (Fig. 3). In the tibia, the thickest cartilage was found in the lateral tibial plateau. In the patella, the thickest cartilage was generally located centrally with some variations where the distribution bifurcated across the patella’s vertical ridge.

The first mode qualitatively explained differences in the overall size of the knee bones as well as the overall thickness of the cartilage (Fig. 3). From this model, we observed that bone size was coupled to cartilage thickness, with larger bones having thicker cartilage. Male sex and tall height were correlated with larger bones and thicker cartilage (more negative PC1) (Table 2).

The second mode described coupled changes associated with aspect ratio or width along the mediolateral (ML) direction of the bones and cartilage. In the bones (Fig. 4), we observed coupled changes in the size of the diaphysis relative to the ML width of the epiphyses in both the femur and tibia. In the tibia, we also saw a widening of the tibial plateau with decreased diaphysis size. This widening was accompanied with a wider distal femur and patella. Additionally, this mode captured changes in the anterior–posterior thickness of the patella, which was lower with decreased diaphysis size. In the cartilage, we observed a number of coupled changes across the bones (Fig. 3), including thinning of cartilage in both the femur and tibia with decreased diaphysis size, and a proximal–distal shift in the location of the thickest cartilage in both the patella and trochlea cartilage. Furthermore, thinning of the trochlear cartilage on the femur was coupled with thicker patellar cartilage on the lateral side of the vertical ridge and thinner cartilage on the medial side. Male sex and shorter height were correlated with high aspect ratio (more negative PC2) (Table 2).

In the bone, the fourth mode describes the depth of the femoral condyle in the anteroposterior (AP) direction. It also describes changes to the angle of the femoral shaft axis to the epicondylar axis (femoral shaft angle) and the mediolateral tilt of the tibial plateau. A deeper distal femur was accompanied with a larger femoral shaft angle and a flatter tibial plateau (Fig. 4). In addition, a deeper distal femur was accompanied with a less pronounced ridge located on the medial side of the patella. A more pronounced ridge was coupled with a bifurcation in cartilage thickness across the patellar ridge. This increased condylar depth was accompanied with thicker cartilage throughout and a more superiorly located cartilage boundary in the trochlea (Fig. 3). Male sex was correlated with greater condylar depth (more negative PC4) (Table 2).

The logistic regression model trained on principal components 1, 2 and 4 (Fig. 5) of the SSFM scores characterised the morphology of the bone and cartilage between males and females to an accuracy of 94.1% within the training set. Cross-validation of this model yielded a classification accuracy of 92.2% (95% CI [81.1%, 97.8%]) and an area under the curve (AUC) of 94.1%. The decision boundary (Fig. 5), shows that females had generally smaller bones (positive PC1) accompanied with lower aspect ratios (positive PC2), steeper mediolateral tibial tilt, more pronounced patellar ridge, and thinner distal femur in the AP direction (positive PC4).

3D scatter of principal components 1, 2, and 4, of the training set by sex, showing the decision boundary plane (grey) of the logistic regression model, and, a vector (green) that passes through the average male and female knee.

The purpose of this study was to develop a statistical shape and field model of the femur, tibia, patella, and their cartilages to characterise the coupled variation in bone morphology and cartilage thickness in a cohort of healthy knees. The bone shape and cartilage field model captured over 90% of the variance with the first five modes. This model was interrogated to explore the effects of sex, height, body mass, and age, on the morphology of the bones and cartilage thickness distribution, and finally used to characterise the morphological differences between males and females. We found that sex and height strongly influenced knee joint morphology, while body mass and age did not influence morphology in this young adult cohort.

Our results showed that cartilage thickness and its distribution was coupled to bone morphology, including both size (mode 1) and shape variations (mode 2 onwards).

The greatest influence on the overall cartilage thickness was bone size, where larger bones displayed thicker cartilage (Fig. 3, mode 1). Although this seems intuitive, studies across species have shown that cartilage thickness does not necessarily scale with size. For example, the equine cartilage is thinner than human cartilage in smaller equivalent joints34,35, despite equine cartilage experiencing much higher stresses than humans36,37,38,39. Within the lower limb joints of humans, Shepherd and Seedhom40 hypothesised that cartilage thickness is a function of joint congruency, due to the hip and ankle having thinner cartilage than the knee in their sample of 11 cadavers. This is believed to be related to how much stress the cartilage needs to bear6,17. We expect an increase in body mass with bone size to result in increased stress that requires thicker cartilage to bear. Studies have also reported a correlation between anthropomorphic measurements such as height and mass to mean cartilage thickness15,17,40, implying a relationship between bone size, stress, and cartilage thickness. Our results confirm a part of this finding in a larger sample (n = 51), where height correlated to increased overall cartilage thickness (mode 1). Previous studies indicated that mass only described some of the variability in cartilage thickness15,17,40. However, our results show that mass was not correlated with any of the modes of variation. This is unexpected, since mass is typically a predictor of muscle volume and thus, muscle forces, which dominate joint contact forces and are related to the load history during growth and adolescence when the bone and cartilage is undergoing a modelling process41. The literature has shown that incorporating information beyond mass, such as height, can enable better predictions of joint moments during locomotion and muscle volumes in the lower limb42. However, we found that none of the modes investigated were correlated with mass, mass × height, or BMI when height was included in the model.

Our model showed that cartilage thickness was coupled to bone shape variations. Knees with higher aspect ratio and smaller diaphyses (mode 2) possessed thinner cartilage overall and the location of the thickest cartilage in the patella and trochlea shifted proximally, whereas knees with lower aspect ratio and larger diaphyses possessed thicker cartilage. Increased diaphysis size increases the bone section modulus which typically results from increased load experienced during growth and development43. This increased load bearing may explain why this feature is accompanied with thicker cartilage.

The mediolateral slope of the tibial plateau and the angle of the femoral shaft axis were also coupled (mode 4) where a larger angle was accompanied with thicker distal femur in the AP direction. Recent studies have reported that a high Q-angle (similar to the femoral shaft axis angle) is negatively correlated with cartilage thickness44, especially in the lateral compartment45. Kusiak and Kawczyński45 hypothesised that this is due to increased compression imposed by the patellar cartilage on this condyle. Our results do not appear to support this hypothesis, but instead suggest that the patellar and trochlear cartilages get thicker with increased femoral shaft axis angle. Without knowledge of the muscle force distribution on the patella and the contact mechanics of the patellofemoral articulation, it is difficult to draw any conclusions regarding patellofemoral joint mechanobiology. Further modelling work might uncover these form-function relationships and sex differences.

Sexual dimorphism was present in three modes of the model (Fig. 5). Size (mode 1), aspect ratio (mode 2) and condylar depth (mode 4) of both the bone and cartilage were important for classifying sex. On average, this meant that females presented a smaller joint with thinner cartilage, with a lower aspect ratio and shallower distal femur compared to males. This is somewhat expected and conforms with previous studies that have demonstrated sexual dimorphism in bone morphology10,46. Kim et al.46 showed that four key morphometric measures including condylar width, width of the medial condyle, and the depths of the lateral condyle and intercondylar notch were important for sex classification. These measures were captured in the modes of our model, for example, size and aspect ratio capture all four key measures identified, while condylar depth captures the depth of the lateral condyle. Prior studies have also reported differences in cartilage volume between males and females, where males had larger cartilage volumes than females16,17,18,19,47. Our results add to these previous findings by giving further clarity as to how the spatial distribution of the cartilage thickness varies in the population along with variations to bone morphology. We expect further modelling work to reveal how cartilage stress is maintained across different bone and cartilage morphologies.

There are many implications of this work. It is important to consider both the cartilage morphology as well as the bone morphology as they both influence the mechanics. This is particularly important when considering joint contact models as the cartilage thickness influences the mechanical response of the cartilage and the internal stresses experienced by the tissue. In implant design, the articulating geometry needs to be considered, not just the subchondral bone. This method can be extended to include an OA cohort to compare the coupled changes in bone and cartilage morphology as a function of OA progression.

There are several limitations in this work that should be acknowledged. Firstly, the conclusions drawn in this paper are limited to knees with healthy cartilage, as determined by a musculoskeletal radiologist as part of the inclusion criteria from the original experimental study16. Some of these participants had patellofemoral pain, which may influence the results. Across the age-range of our study cohort we did not observe any morphological differences due to age. It is possible that widening the age-range of the study cohort would lead to detection of age-related differences. We expect to observe age-related morphological differences in the bone and cartilage in an adolescent population where growth and development occurs, and in an older population where OA degeneration may be occurring. Our methodology normalised for bone orientation, a factor that may influence the interaction between bones and hence the cartilage thickness. As the scope of this study was to characterise the coupled variation in bone morphology and cartilage thickness, we did not investigate how these relationships are affected when normalised for bone size or maximum cartilage thickness. Lastly, as with any investigation of bone and cartilage morphology using medical imaging data, our results are subject to manual segmentation errors which are typically < 0.5 mm for knee joint structures32.

In conclusion, we have characterised and described the coupled variation in bone morphology and cartilage thickness in a cohort of 51 young adult knees. We observed that cartilage is thicker with increased bone size, diaphysis size, and condylar depth. Sexual dimorphism was present in these three modes which may be used to understand the dimorphism that exists in OA, knee ligament injuries, and patellofemoral pain. Lastly, we found that bone morphology and cartilage thickness were strongly correlated with height but not body mass or age.

The data that support the findings of this study are openly available in SimTK.org Anatomical Knee at https://simtk.org/projects/anatomicalknee48.

Carter, D. R. et al. The mechanobiology of articular cartilage development and degeneration. Clin. Orthop. Relat. Res. 427, S69–S77. https://doi.org/10.1097/01.blo.0000144970.05107.7e (2004).

Anderson, A. E., Ellis, B. J., Maas, S. A. & Weiss, J. A. Effects of idealized joint geometry on finite element predictions of cartilage contact stresses in the hip. J. Biomech. 43, 1351–1357. https://doi.org/10.1016/j.jbiomech.2010.01.010 (2010).

Article  PubMed  PubMed Central  Google Scholar 

Arnold, A. S. & Delp, S. L. Rotational moment arms of the medial hamstrings and adductors vary with femoral geometry and limb position: Implications for the treatment of internally rotated gait. J. Biomech. 34, 437–447. https://doi.org/10.1016/S0021-9290(00)00232-3 (2001).

CAS  Article  PubMed  Google Scholar 

Thambyah, A., Nather, A. & Goh, J. Mechanical properties of articular cartilage covered by the meniscus. Osteoarthritis Cartilage 14, 580–588. https://doi.org/10.1016/j.joca.2006.01.015 (2006).

CAS  Article  PubMed  Google Scholar 

Ateshian, G. A., Soslowsky, L. J. & Mow, V. C. Quantitation of articular surface topography and cartilage thickness in knee joints using stereophotogrammetry. J. Biomech. 24, 761–776. https://doi.org/10.1016/0021-9290(91)90340-S (1991).

CAS  Article  PubMed  Google Scholar 

Li, G. et al. The cartilage thickness distribution in the tibiofemoral joint and its correlation with cartilage-to-cartilage contact. Clin. Biomech. 20, 736–744. https://doi.org/10.1016/j.clinbiomech.2005.04.001 (2005).

Zhang, J., Fernandez, J., Hislop-Jambrich, J. & Besier, T. F. Lower limb estimation from sparse landmarks using an articulated shape model. J. Biomech. 49, 3875–3881. https://doi.org/10.1016/j.jbiomech.2016.10.021 (2016).

Zhang, J. & Besier, T. F. Accuracy of femur reconstruction from sparse geometric data using a statistical shape model. Comput. Methods Biomech. Biomed. Eng. 20, 566–576. https://doi.org/10.1080/10255842.2016.1263301 (2017).

Bredbenner, T. L. et al. Statistical shape modeling describes variation in tibia and femur surface geometry between Control and Incidence groups from the Osteoarthritis Initiative database. J. Biomech. 43, 1780–1786. https://doi.org/10.1016/j.jbiomech.2010.02.015 (2010).

Article  PubMed  PubMed Central  Google Scholar 

Lynch, J. T. et al. Statistical shape modelling reveals large and distinct subchondral bony differences in osteoarthritic knees. J. Biomech. 93, 177–184. https://doi.org/10.1016/j.jbiomech.2019.07.003 (2019).

Felson, D. T. et al. The incidence and natural history of knee osteoarthritis in the elderly, the framingham osteoarthritis study. Arthritis Rheum. 38, 1500–1505. https://doi.org/10.1002/art.1780381017 (1995).

CAS  Article  PubMed  Google Scholar 

Wise, B. L. et al. Bone shape mediates the relationship between sex and incident knee osteoarthritis. BMC Musculoskelet. Disord. 19, 331. https://doi.org/10.1186/s12891-018-2251-z (2018).

Article  PubMed  PubMed Central  Google Scholar 

Pedoia, V. et al. Three-dimensional MRI-based statistical shape model and application to a cohort of knees with acute ACL injury. Osteoarthr. Cartil. 23, 1695–1703. https://doi.org/10.1016/j.joca.2015.05.027 (2015).

Cohen, Z. A. et al. Knee cartilage topography, thickness, and contact areas from MRI: In-vitro calibration and in-vivo measurements. Osteoarthr. Cartil. 7, 95–109. https://doi.org/10.1053/joca.1998.0165 (1999).

Eckstein, F., Winzheimer, M., Hohe, J., Englmeier, K. H. & Reiser, M. Interindividual variability and correlation among morphological parameters of knee joint cartilage plates: Analysis with three-dimensional MR imaging. Osteoarthr. Cartil. 9, 101–111. https://doi.org/10.1053/joca.2000.0365 (2001).

Draper, C. et al. Is cartilage thickness different in young subjects with and without patellofemoral pain?. Osteoarthr. Cartil. 14, 931–937. https://doi.org/10.1016/j.joca.2006.03.006 (2006).

Connolly, A., FitzPatrick, D., Moulton, J., Lee, J. & Lerner, A. Tibiofemoral cartilage thickness distribution and its correlation with anthropometric variables. Proc. Inst. Mech. Eng. H 222, 29–39. https://doi.org/10.1243/09544119JEIM306 (2008).

CAS  Article  PubMed  Google Scholar 

Ding, C., Cicuttini, F., Scott, F., Glisson, M. & Jones, G. Sex differences in knee cartilage volume in adults: Role of body and bone size, age and physical activity. Rheumatology 42, 1317–1323. https://doi.org/10.1093/rheumatology/keg374 (2003).

CAS  Article  PubMed  Google Scholar 

Faber, S. C. et al. Gender differences in knee joint cartilage thickness, volume and articular surface areas: Assessment with quantitative three-dimensional MR imaging. Skeletal Radiol. 30, 144–150. https://doi.org/10.1007/s002560000320 (2001).

CAS  Article  PubMed  Google Scholar 

Cootes, T. F., Cooper, D. H., Taylor, C. J. & Graham, J. Trainable method of parametric shape description. Image Vis. Comput. 10, 289–294. https://doi.org/10.1016/0262-8856(92)90044-4 (1992).

Dryden, I. L. & Mardia, K. V. Statistical shape analysis: Wiley series in probability and statistics (Wiley, 1998).

Schneider, M. et al. Men and women have similarly shaped carpometacarpal joint bones. J. Biomech. 48, 3420–3426. https://doi.org/10.1016/j.jbiomech.2015.05.031 (2015).

CAS  Article  PubMed  PubMed Central  Google Scholar 

Schneider, M. T. et al. Trapeziometacarpal joint contact varies between men and women during three isometric functional tasks. Med. Eng. Phys. 50, 43–49. https://doi.org/10.1016/j.medengphy.2017.09.002 (2017).

Article  PubMed  PubMed Central  Google Scholar 

Schneider, M. et al. Early morphologic changes in trapeziometacarpal joint bones with osteoarthritis. Osteoarthr. Cartil. 26, 1338–1344. https://doi.org/10.1016/j.joca.2018.06.008 (2018).

Pal, S. et al. Patellofemoral cartilage stresses are most sensitive to variations in vastus medialis muscle forces. Comput. Methods Biomech. Biomed. Eng. 22, 206–216. https://doi.org/10.1080/10255842.2018.1544629 (2019).

Besier, T. F. et al. The role of cartilage stress in patellofemoral pain. Med. Sci. Sports Exerc. 47, 2416. https://doi.org/10.1249/MSS.0000000000000685 (2015).

CAS  Article  PubMed  PubMed Central  Google Scholar 

Treece, G., Prager, R. & Gee, A. Stradwin 6.0. http://mi.eng.cam.ac.uk/~rwp/stradwin.

Zhang, J. et al. The MAP Client: User-Friendly Musculoskeletal Modelling Workflows. In: Bello, F., Cotin, S. (eds) Biomedical Simulation. ISBS. Lecture Notes in Computer Science 8789, 182–192. https://doi.org/10.1007/978-3-319-12057-7_21 (2014).

Corsini, M., Cignoni, P. & Scopigno, R. Efficient and flexible sampling with blue noise properties of triangular meshes. IEEE Trans. Visual Comput. Graphics 18, 914–924. https://doi.org/10.1109/TVCG.2012.34 (2012).

Zhang, J., Ackland, D. & Fernandez, J. Point-cloud registration using adaptive radial basis functions. Comput. Methods Biomech. Biomed. Eng. 21, 498–502. https://doi.org/10.1080/10255842.2018.1484914 (2018).

Zhang, J., Hislop-Jambrich, J. & Besier, T. F. Predictive statistical models of baseline variations in 3-D femoral cortex morphology. Med. Eng. Phys. 38, 450–457. https://doi.org/10.1016/j.medengphy.2016.02.003 (2016).

Rooks, N. et al. A method to compare heterogeneous types of bone and cartilage meshes. J. Biomech. Eng. 20, 20. https://doi.org/10.1115/1.4051281 (2021).

Demšar, J. et al. Orange: Data mining toolbox in Python. J. Mach. Learn. Res. 14, 2349–2353 (2013).

Frisbie, D. D., Cross, M. W. & McIlwraith, C. W. A comparative study of articular cartilage thickness in the stifle of animal species used in human pre-clinical studies compared to articular cartilage thickness in the human knee. Vet. Comp. Orthop. Traumatol. 19, 142–146. https://doi.org/10.1055/s-0038-1632990 (2006).

CAS  Article  PubMed  Google Scholar 

Firth, E. C. The response of bone, articular cartilage and tendon to exercise in the horse. J. Anat. 208, 513–526. https://doi.org/10.1111/j.1469-7580.2006.00547.x (2006).

Article  PubMed  PubMed Central  Google Scholar 

Frazer, L. L., Santschi, E. M. & Fischer, K. J. Impact of a void in the equine medial femoral condyle on bone stresses and peak contact pressures in a finite element model. Vet. Surg. 48, 237–246. https://doi.org/10.1111/vsu.13139 (2019).

Rao, C. et al. A statistical finite element model of the knee accounting for shape and alignment variability. Med. Eng. Phys. 35, 1450–1456. https://doi.org/10.1016/j.medengphy.2013.03.021 (2013).

Farrokhi, S., Keyak, J. H. & Powers, C. M. Individuals with patellofemoral pain exhibit greater patellofemoral joint stress: A finite element analysis study. Osteoarthr. Cartil. 19, 287–294. https://doi.org/10.1016/j.joca.2010.12.001 (2011).

Harrison, S. M., Whitton, R. C., Kawcak, C. E., Stover, S. M. & Pandy, M. G. Evaluation of a subject-specific finite-element model of the equine metacarpophalangeal joint under physiological load. J. Biomech. 47, 65–73. https://doi.org/10.1016/j.jbiomech.2013.10.001 (2014).

Shepherd, D. E. T. & Seedhom, B. B. Thickness of human articular cartilage in joints of the lower limb. Ann. Rheum. Dis. 58, 27–34. https://doi.org/10.1136/ard.58.1.27 (1999).

CAS  Article  PubMed  PubMed Central  Google Scholar 

Jones, G. et al. Knee articular cartilage development in children: A longitudinal study of the effect of sex, growth, body composition, and physical activity. Pediatr. Res. 54, 230–236. https://doi.org/10.1203/01.PDR.0000072781.93856.E6 (2003).

ADS  Article  PubMed  Google Scholar 

Handsfield, G. G., Meyer, C. H., Hart, J. M., Abel, M. F. & Blemker, S. S. Relationships of 35 lower limb muscles to height and body mass quantified using MRI. J. Biomech. 47, 631–638. https://doi.org/10.1016/j.jbiomech.2013.12.002 (2014).

Ruff, C. Growth in bone strength, body size, and muscle size in a juvenile longitudinal sample. Bone 33, 317–329. https://doi.org/10.1016/S8756-3282(03)00161-3 (2003).

Ayşe Aydemir, E., Hamarat, H. & Musmul, A. Relationship between Q-angle and articular cartilage in female patients with symptomatic knee osteoarthritis: Ultrasonographic and radiologic evaluation. Arch. Rheumatol. 32, 347. https://doi.org/10.5606/ArchRheumatol.2017.6145 (2017).

Kusiak, M. & Kawczyński, A. Ultrasonographic assessment of articular cartilage of the femoral condyle in patients with an increased Q-angle. J. Ultrasonogr. 18, 181. https://doi.org/10.15557/jou.2018.0027 (2018).

Kim, D.-I., Kwak, D.-S. & Han, S.-H. Sex determination using discriminant analysis of the medial and lateral condyles of the femur in Koreans. Forensic Sci. Int. 233, 121–125. https://doi.org/10.1016/j.forsciint.2013.08.028 (2013).

Cicuttini, F. et al. Gender differences in knee cartilage volume as measured by magnetic resonance imaging. Osteoarthr. Cartil. 7, 265–271. https://doi.org/10.1053/joca.1998.0200 (1999).

Schneider MT, Rooks N, Besier TF. (2021) Anatomical Knee. https://simtk.org/projects/anatomicalknee

This work was supported by the Auckland Bioengineering Institute and by the National Institutes of Health (Award numbers 1R01EB024573-01 and 1R01EB005790-01). We thank Christine Draper and Garry Gold for help with data collection. The content is solely the responsibility of the authors and does not necessarily represent the official views of the Auckland Bioengineering Institute or the National Institutes of Health.

Auckland Bioengineering Institute, The University of Auckland, Level 6, 70 Symonds Street, Auckland, 1010, New Zealand

Marco Tien-Yueh Schneider, Nynke Rooks & Thor Besier

Department of Engineering Science, The University of Auckland, Auckland, New Zealand

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

You can also search for this author in PubMed  Google Scholar

All authors have made significant contribution to the conception and design of the study, acquisition of data, and interpretation of data. Development of the computational models and, drafting of the paper were performed by M.T.S. T.B. and N.R. contributed to critical revision of the paper and approved the final version of the manuscript to be published.

Correspondence to Marco Tien-Yueh Schneider.

The authors declare no competing interests.

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

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Schneider, M.TY., Rooks, N. & Besier, T. Cartilage thickness and bone shape variations as a function of sex, height, body mass, and age in young adult knees. Sci Rep 12, 11707 (2022). https://doi.org/10.1038/s41598-022-15585-w

DOI: https://doi.org/10.1038/s41598-022-15585-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.

Scientific Reports (Sci Rep) ISSN 2045-2322 (online)

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.