Using 3D MRI Bone Shape to Predict Pre-Osteoarthritis of the Knee 2 Years After Anterior Cruciate Ligament Reconstruction

BACKGROUND

Anterior cruciate ligament (ACL) injury increases risks for osteoarthritis (OA), a poorly modifiable and disabling condition. Joint changes of potentially reversible pre-OA have been described just 2 years after ACL reconstruction (ACLR) when early bone shape changes have also been reported.

PURPOSE

This study evaluates relationships between interlimb differences in tibiofemoral bone shape derived from statistical shape modeling (SSM) of magnetic resonance imaging (MRI) and participant factors on patient-reported outcomes 2 years after unilateral ACLR.

STUDY DESIGN

Cross-sectional study; Level of evidence, 3.

METHODS

SSM-derived tibiofemoral bone shape and subchondral bone area were assessed from bilateral knee MRI scans of 72 participants with unilateral ACLR (mean age, 34 ± 11 years; 32 women) and compared with a reference cohort of 398 older individuals without OA (mean age, 50 ± 3 years; 213 women). Multivariable logistic regression models examined relationships between participant and surgical factors with interlimb differences in bone shapes or subchondral bone areas. Relationships between patient-reported outcomes and the interlimb differences in bone shape and subchondral area were examined using similar models.

RESULTS

Bone shape scores and subchondral bone areas were greater (more OA-like) in ACLR knees than uninjured contralateral knees in every bone metric tested (P≤ .001). Interlimb differences in femur shape scores of participants with ACLR were 65% greater (P < .001) than those of the significantly older reference cohort. Taller height, medial meniscal tears, and decreasing age were associated with larger interlimb differences in shape scores and subchondral areas (P < .05). Bone-patellar tendon-bone (BPTB) autograft recipients demonstrated greater interlimb subchondral area differences compared with allograft recipients (P < .05). Interlimb differences for hamstring autograft recipients did not differ from those with BPTB or allograft. Greater interlimb differences in medial femur subchondral areas were associated with worse patient-reported Knee injury and Osteoarthritis Outcome Score Symptoms (R = 0.27; P = .040).

CONCLUSION

Even in the absence of radiographic OA, just 2 years after unilateral ACLR patients showed greater bone shape scores and subchondral areas consistent with pre-OA in their ACLR knees. Furthermore, greater medial femur bone areas were weakly associated with worse symptoms. Patients who are younger, are taller, have meniscal tears, or have BPTB grafts may be at increased risk for bony asymmetries 2 years after ACLR.

Understanding the form and function in Chinese bound foot from last-generation cases

Purpose: Foot adaptation in the typically developed foot is well explored. In this study, we aimed to explore the form and function of an atypical foot, the Chinese bound foot, which had a history of over a thousand years but is not practised anymore. Methods: We evaluated the foot shape and posture via a statistical shape modelling analysis, gait plantar loading distribution via gait analysis, and bone density adaptation via implementing finite element simulation and bone remodelling prediction. Results: The atypical foot with binding practice led to increased foot arch and vertically oriented calcaneus with larger size at the articulation, apart from smaller metatarsals compared with a typically developed foot. This shape change causes the tibia, which typically acts as a load transfer beam and shock absorber, to extend its function all the way through the talus to the calcaneus. This is evident in the bound foot by i) the reduced center of pressure trajectory in the medial-lateral direction, suggesting a reduced supination-pronation; ii) the increased density and stress in the talus-calcaneus articulation; and iii) the increased bone growth in the bound foot at articulation joints in the tibia, talus, and calcaneus. Conclusion: Knowledge from the last-generation bound foot cases may provide insights into the understanding of bone resorption and adaptation in response to different loading profiles.

Determining the variation in premaxillary and dentary bone morphology that may underlie beak shape between two pure layer lines

Beak treatment is an effective method of reducing the damage inflicted by severe feather pecking (SFP) but there is significant pressure to eliminate these treatments and rely solely on alternative strategies. Substantial variation in beak shape exists within non-beak treated layer flocks and beak shape appears to be heritable. There is the potential to use this pre-existing variation and genetically select for hens whose beak shapes are less apt to cause damage during SFP. To do this, we must first understand the range of phenotypes that exist for both the external beak shape and the bones that provide its structure. The objective of this study was to determine the variation in premaxillary (within the top beak) and dentary (within the bottom beak) bone morphology that exists in 2 non-beak treated pure White Leghorn layer lines using geometric morphometrics to analyze radiographs. Lateral head radiographs were taken of 825 hens and the premaxillary and dentary bones were landmarked. Landmark coordinates were standardized by Procrustes superimposition and the covariation was analyzed by principal components analysis and multivariate regression using Geomorph (an R package). Three principal components (PCs) explained 85% of total premaxillary bone shape variation and showed that the shape ranged from long and narrow with pointed bone tips to short and wide with more curved tips. Two PCs explained 81% of total dentary bone shape variation. PC1 described the dentary bone length and width and PC2 explained the angle between the bone tip and its articular process. For both bones, shape was significantly associated with bone size and differed significantly between the two lines. Bone size accounted for 42% of the total shape variation for both bones. Together, the results showed a range of phenotypic variation in premaxillary and dentary bone shape, which in turn may influence beak shape. These bone phenotypes will guide further quantitative genetic and behavioral analyses that will help identify which beaks shapes cause the least damage when birds engage in SFP.

Multiparametric MRI characterization of knee articular cartilage and subchondral bone shape in collegiate basketball players

Magnetic resonance imaging (MRI) is commonly used to evaluate the morphology of the knee in athletes with high-knee impact; however, complex repeated loading of the joint can lead to biochemical and structural degeneration that occurs before visible morphological changes. In this study, we utilized multiparametric quantitative MRI to compare morphology and composition of articular cartilage and subchondral bone shape between young athletes with high-knee impact (basketball players; n = 40) and non-knee impact (swimmers; n = 25). We implemented voxel-based relaxometry to register all cases to a single reference space and performed a localized compositional analysis of T _1ρ - and T ₂ -relaxation times on a voxel-by-voxel basis. Additionally, statistical shape modeling was employed to extract differences in subchondral bone shape between the two groups. Evaluation of cartilage composition demonstrated a significant prolongation of relaxation times in the medial femoral and tibial compartments and in the posterolateral femur of basketball players in comparison to relaxation times in the same cartilage compartments of swimmers. The compositional analysis also showed depth-dependent differences with prolongation of the superficial layer in basketball players. For subchondral bone shape, three total modes were found to be significantly different between groups and related to the relative sizes of the tibial plateaus, intercondylar eminences, and the curvature and concavity of the patellar lateral facet. In summary, this study identified several characteristics associated with a high-knee impact which may expand our understanding of local degenerative patterns in this population.

Deletion of Fibroblast growth factor 9 globally and in skeletal muscle results in enlarged tuberosities at sites of deltoid tendon attachments

BACKGROUND

The growth of most bony tuberosities, like the deltoid tuberosity (DT), rely on the transmission of muscle forces at the tendon-bone attachment during skeletal growth. Tuberosities distribute muscle forces and provide mechanical leverage at attachment sites for joint stability and mobility. The genetic factors that regulate tuberosity growth remain largely unknown. In mouse embryos with global deletion of fibroblast growth factor 9 (Fgf9), the DT size is notably enlarged. In this study, we explored the tissue-specific regulation of DT size using both global and targeted deletion of Fgf9.

RESULTS

We showed that cell hypertrophy and mineralization dynamics of the DT, as well as transcriptional signatures from skeletal muscle but not bone, were influenced by the global loss of Fgf9. Loss of Fgf9 during embryonic growth led to increased chondrocyte hypertrophy and reduced cell proliferation at the DT attachment site. This endured hypertrophy and limited proliferation may explain the abnormal mineralization patterns and locally dysregulated expression of markers of endochondral development in Fgf9^null attachments. We then showed that targeted deletion of Fgf9 in skeletal muscle leads to postnatal enlargement of the DT.

CONCLUSION

Taken together, we discovered that Fgf9 may play an influential role in muscle-bone cross-talk during embryonic and postnatal development.

Longitudinal analysis of the contribution of 3D patella and trochlear bone shape on patellofemoral joint osteoarthritic features

To explore bone shape features that are associated with patellofemoral joint (PFJ) osteoarthritic features. Thirty subjects with PFJ degeneration (six males, 53.2 ± 9.8 years) and 23 controls (12 males, 48.1 ± 10.6 years) were included. Magnetic resonance (MR) assessment was performed to provide bone segmentation, morpholgocial grading, and cartilage relaxation times. In addition, subject self-reported symptoms were reported. Logistic regressions were used to identify the shape features that were associated with the presence and worsening of PFJ morphological lesions over 3 years, and worsening of self-reported symptoms. Statistical parametric mapping was used to evaluate the associations between shape features and cartilage relaxation times at 3 years. Results indicated that subjects with PFJ degeneration exhibited a trochlea with longer lateral condyle and shallower trochlear groove (adjusted odds ratio [OR] = 0.30; 95% confidence interval [CI]: 0.10, 0.86; P = .025). Subjects with worsening of PFJ degeneration exhibited a patella with equally distributed facets (adjusted OR = 3.14; 95% CI: 1.05, 9.37; P = .040) and lateral bump (adjusted OR = 0.14; 95% CI: 0.02, 0.83; P = .030). No shape features were associated with worsening of self-reported symptoms. Elevated T_1ρ and T₂ times at 3 years were associated with a patella with a lateral hook, equally distributed facets, round and thick as well as a trochlea larger in size (R = 0.38~0.46, P = .015~.025). The study demonstrated the ability of 3D statistical shape modeling to quantify patella and trochlear bone shape features that are associated with the presence and progression of PFJ osteoarthritic features.

Differences in Tibia Shape in Organically Reared Chicken Lines Measured by Means of Geometric Morphometrics

Simple Summary

Organic poultry production should use only genetic lines and animals resistant to disease and well adapted to live outdoor, according to principles, rules, and requirements of organic farming systems. When broiler’s walking performance is reduced animals are not suitable for outdoor rearing. There is a straightforward relationship between bone health and growth rate in broilers. Body and breast weight play an important role in leg disorders. During the last decades, genetic selection has led to high producing broilers over the time. Unfortunately, fast growth may negatively influence correct leg development, reducing walking performance, and raising welfare issues. Leg abnormalities could represent a criterion for the choice of genetic lines suitable for organic production. A method for their early detection was developed in this study by means of Geometric Morphometrics (GM) that represents a tool for bone shape analysis and its correlation with walking capability. A valuable information emerged from the present study in relation to broiler intrinsic adaptability to organic production.

Abstract

In the present study, the conformation of the tibia of seven genetic lines of broilers was analyzed by Geometric Morphometrics and correlated to carcass weight and walking ability. The used chicken genetic lines were classified as fast, medium, or slow growing and ranked for their walking ability. Six chicken types were reared in an organic farm and slaughtered at 81 days of age while one slow-growing and highly walking line (Naked Neck) was reared in a commercial farm and used as external reference for moving activity and growth speed. A mixed landmarks and semi-landmarks model was applied to the study of tibia shape. Results of this study showed that: (i) body weight gain was positively correlated to the curvature of the antero-posterior axis of the tibia; (ii) the shape of the tibia and the active walking behavior were significantly correlated; (iii) walking and not-walking genetic lines could be discriminated in relation to the overall shape of the tibia; (iv) a prevalence of static behavior was correlated to a more pronounced curvature of the antero-posterior axis of the tibia. Results of this study revealed that the walking genetic types have a more functional and natural tibia conformation. This easy morphologic method for evaluating tibia shape could help to characterize the adaptability of genotypes to organic and outdoor rearing.

Learning osteoarthritis imaging biomarkers from bone surface spherical encoding

PURPOSE

To learn bone shape features from spherical bone map of knee MRI images using established convolutional neural networks (CNN) and use these features to diagnose and predict osteoarthritis (OA).

METHODS

A bone segmentation model was trained on 25 manually annotated 3D MRI volumes to segment the femur, tibia, and patella from 47 078 3D MRI volumes. Each bone segmentation was converted to a 3D point cloud and transformed into spherical coordinates. Different fusion strategies were performed to merge spherical maps obtained by each bone. A total of 41 822 merged spherical maps with corresponding Kellgren-Lawrence grades for radiographic OA were used to train a CNN classifier model to diagnose OA using bone shape learned features. Several OA Diagnosis models were tested and the weights for each trained model were transferred to the OA Incidence models. The OA incidence task consisted of predicting OA from a healthy scan within a range of eight time points, from 1 y to 8 y. The validation performance was compared and the test set performance was reported.

RESULTS

The OA Diagnosis model had an area-under-the-curve (AUC) of 0.905 on the test set with a sensitivity and specificity of 0.815 and 0.839. The OA Incidence models had an AUC ranging from 0.841 to 0.646 on the test set for the range from 1 y to 8 y.

CONCLUSION

Bone shape was successfully used as a predictive imaging biomarker for OA. This approach is novel in the field of deep learning applications for musculoskeletal imaging and can be expanded to other OA biomarkers.

[18 F]-Sodium Fluoride PET/MR Imaging for Bone-Cartilage Interactions in Hip Osteoarthritis: A Feasibility Study

This study characterized the distribution of [¹⁸ F]-sodium fluoride (NaF) uptake and blood flow in the femur and acetabulum in hip osteoarthritis (OA) patients to find associations between bone remodeling and cartilage composition in the presence of morphological abnormalities using simultaneous positron emission tomography and magnetic resonance imaging (PET/MR), quantitative magnetic resonance imaging (MRI) and femur shape modeling. Ten patients underwent a [¹⁸ F]-NaF PET/MR dynamic scan of the hip simultaneously with: (i) fast spin-echo CUBE for morphology grading and (ii) T_1ρ /T₂ magnetization-prepared angle-modulated partitioned k-space spoiled gradient echo snapshots for cartilage, bone segmentation, bone shape modeling, and T_1ρ /T₂ quantification. The standardized uptake values (SUVs) and Patlak kinetic parameter (K_pat ) were calculated for each patient as PET outcomes, using an automated post-processing pipeline. Shape modeling was performed to extract the variations in bone shapes in the patients. Pearson's correlation coefficients were used to study the associations between bone shapes, PET outcomes, and patient reported pain. Direct associations between quantitative MR and PET evidence of bone remodeling were established in the acetabulum and femur. Associations of shaft thickness with SUV in the femur (p = 0.07) and K_pat in the acetabulum (p = 0.02), cam deformity with acetabular score (p = 0.09), osteophytic growth on the femur head with K_pat (p = 0.01) were observed. Pain had increased correlations with SUV in the acetabulum (p = 0.14) and femur (p = 0.09) when shaft thickness was accounted for. This study demonstrated the ability of [¹⁸ F]-NaF PET-MRI, 3D shape modeling, and quantitative MRI to investigate cartilage-bone interactions and bone shape features in hip OA, providing potential investigative tools to diagnose OA. © 2019 The Authors. Journal of Orthopaedic Research^® published by Wiley Periodicals, Inc. on behalf of Orthopaedic Research Society J Orthop Res 37:2671-2680, 2019.

MR study of longitudinal variations in proximal femur 3D morphological shape and associations with cartilage health in hip osteoarthritis

The goal of this study was to use quantitative MRI analysis to longitudinally observe the relationship between 3D proximal femur shape and hip joint degenerative changes. Forty-six subjects underwent unilateral hip MR imaging at three time points (baseline, 18 and 36 months). 3D shape analysis, hip cartilage T_1ρ /T₂ relaxation time quantification, and SHOMRI MRI grading were performed at each time point. Subjects were grouped based on KL, SHOMRI, and HOOS pain scores. Associations between these score groupings, time, and longitudinal variation in shape, were analyzed using a generalized estimating equation. One-way ANCOVA was conducted to evaluate change in shape as a predictor of the worsening of degenerative changes at 36 months. Our results demonstrated that subjects displayed an increase in the volume of the femoral head and neck (Mode 3) over time. This shape mode was significantly more prevalent in patients that reported pain. Longitudinal changes in this shape mode also served as borderline predictors of elevated T_1ρ values (p = 0.055) and of cartilage lesions (p = 0.068). Subjects showed a change in the Femoral Neck Anteversion angle (FNA) over time (Mode 6). This shape mode showed a significant interaction with the presence of cartilage lesions. The results of this study suggest that specific variations in bone shape quantified through 3D-MRI based Statistical Shape modeling show an observable relationship with hip joint compositional and morphological changes. The shapes observed lead to early degenerative changes, which may lead into OA, thus confirming the important role of bone shape changes in the pathogenesis of OA. © 2018 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res.

Study of the interactions between proximal femur 3d bone shape, cartilage health, and biomechanics in patients with hip Osteoarthritis

In this study quantitative MRI and gait analysis were used to investigate the relationships between proximal femur 3D bone shape, cartilage morphology, cartilage biochemical composition, and joint biomechanics in subject with hip Osteoarthritis (OA). Eighty subjects underwent unilateral hip MR-imaging: T1ρ and T2 relaxation times were extracted through voxel based relaxometry and bone shape was assessed with 3D MRI-based statistical shape modeling. In addition, 3D gait analysis was performed in seventy-six of the studied subjects. Associations between shape, cartilage lesion presence, severity, and cartilage T1ρ and T2 were analyzed with linear regression and statistical parametric mapping. An ad hoc analysis was performed to investigate biomechanics and shape associations. Our results showed that subjects with a higher neck shaft angle in the coronal plane (higher mode 1, coxa valga), thicker femoral neck and a less spherical femoral head (higher mode 5, pistol grip) exhibited more severe acetabular and femoral cartilage abnormalities, showing different interactions with demographics factors. Subjects with coxa valga also demonstrated a prolongation of T1ρ and T2. Subjects with pistol grip deformity exhibited reduced hip internal rotation angles and subjects with coxa valga exhibited higher peak hip adduction moment and moment impulse. The results of this study establish a clear relationship between 3D proximal femur shape variations and markers of hip joint degeneration-morphological, compositional, well as insight on the possible interactions with demographics and biomechanics, suggesting that 3D MRI-based bone shape maybe a promising biomarker of early hip joint degeneration. © 2017 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 36:330-341, 2018.

Analysis of the articular cartilage T1ρ and T2 relaxation times changes after ACL reconstruction in injured and contralateral knees and relationships with bone shape

The objectives of this study were twofold: (1) to evaluate the longitudinal change in cartilage T_1ρ and T₂ 6- and 12-months after ACL reconstruction (ACLR) in both reconstructed and intact contralateral knees with the aim of validating the role of the contralateral knee as an internal control in longitudinal studies; (2) to explore relationships between bone shape at the time of injury and the progression of T_1ρ and T₂ over 12-months after ACLR. T_1ρ and T₂ cartilage relaxation times and 3D MRI-based statistical shape modeling (SSM) of tibia and femur were computed for both knees of forty ACL-injured patients and 15 healthy controls. ACL subjects were scanned 8.4 ± 6.4 weeks after injury (2.4 ± 3.7 weeks prior to ACLR), 6- and 12-months after ACLR. Longitudinal changes in T_1ρ and T₂ values were assessed using linear mixed model, and partial correlation coefficients were calculated between bone shape and longitudinal changes in T_1ρ and T₂ values. Significant longitudinal increases in T_1ρ and T₂ values were observed in reconstructed and contralateral knees 6-months after ACLR. Tibial bone shape features, associated with the medial plateau height and width, were observed to be correlated with cartilage T_1ρ and T₂ progression in reconstructed knees. Our results suggest that caution should be used in considering contralateral knee as internal controls in longitudinal ACL studies and 3D MRI-based-SSM might serve as an imaging biomarker for the early stratification of patients at risk for developing post-traumatic accelerated cartilage degeneration and potentially osteoarthritis after ACL tear. © 2016 Orthopaedic Research Society. Published by Wiley Periodicals, Inc. J Orthop Res 35:707-717, 2017.

Gap Junction in the Teleost Fish Lineage: Duplicated Connexins May Contribute to Skin Pattern Formation and Body Shape Determination

Gap junctions are intercellular channels that allow passage of ions and small molecules between adjacent cells. Gap junctions in vertebrates are composed of connexons, which are an assembly of six proteins, connexins. Docking of two connexons on the opposite cell surfaces forms a gap junction between the cytoplasm of two neighboring cells. Connexins compose a family of structurally related four-pass transmembrane proteins. In mammals, there are ~20 connexins, each of which contributes to unique permeability of gap junctions, and mutations of some connexin-encoding genes are associated with human diseases. Zebrafish has been predicted to contain 39 connexin-encoding genes; the high number can be attributed to gene duplication during fish evolution, which resulted in diversified functions of gap junctions in teleosts. The determination of body shapes and skin patterns in animal species is an intriguing question. Mathematical models suggest principle mechanisms explaining the diversification of animal morphology. Recent studies have revealed the involvement of gap junctions in fish morphological diversity, including skin pattern formation and body shape determination. This review focuses on connexins in teleosts, which are integrated in the mathematical models explaining morphological diversity of animal skin patterns and body shapes.

Evidence for enhanced characterization of cortical bone using novel pQCT shape software

Bone shape, mass, structural geometry, and material properties determine bone strength. This study describes novel software that uses peripheral quantitative computed tomography (pQCT) images to quantify cortical bone shape and investigates whether the combination of shape-sensitive and manufacturer's software enhances the characterization of tibiae from contrasting populations. Existing tibial pQCT scans (4% and 50% sites) from Gambian (n=38) and British (n=38) women were used. Bone mass, cross-sectional area (CSA), and geometry were determined using manufacturer's software; cross-sectional shape was quantified using shape-sensitive software. At 4% site, Gambian women had lower total bone mineral content (BMC: -15.4%), CSA (-13.4%), and trabecular bone mineral density (BMD: -19%), but higher cortical subcortical BMD (6.1%). At 50% site, Gambian women had lower cortical BMC (-7.6%), cortical CSA (-12.6%), and mean cortical thickness (-15.0%), but higher cortical BMD (4.9%) and endosteal circumference (8.0%). Shape-sensitive software supported the finding that Gambian women had larger tibial endosteal circumference (9.8%), thinner mean cortical thickness (-26.5%) but smaller periosteal circumference (-5.6%). Shape-sensitive software revealed that Gambian women had tibiae with shorter maximum width (-7.6%) and thinner cortices (-22% to -41.2%) and more closely resembled a circle or ellipse. Significant differences remained after adjusting for age, height, and weight. In conclusion, shape-sensitive software enhanced the characterization of tibiae in 2 contrasting groups of women.