Bone

PurposeOsteoporosis and associated hip fractures are a major public health concern. Dualenergy X-ray Absorptiometry (DXA) remains the diagnostic gold standard, but its areal (a) bone mineral density (BMD) measurements have limited sensitivity, as many fractures occur at T-scores above -2.5. Three-dimensional (3D) DXA provides compartment-specific volumetric parameters of the hip, potentially improving osteoporosis management. This study aimed to establish reference data for 3D-DXA parameters to improve osteoporosis management and investigate potential compartmental imbalances at the hip. MethodThis multicenter, cross-sectional, population-based observational study (SEIOMM-3D-DXA project), supported by the Spanish Society for Bone and Mineral Metabolism (SEIOMM), analyzed hip DXA scans from 1366 Spanish men and women across six centers. 3D-DXA analyses were conducted using the 3D-Shaper software (3D-Shaper Medical, Barcelona, Spain), producing estimates of trabecular volumetric (v) BMD and cortical surface (s) BMD. Age- and sex-specific reference curves were generated using the LMS method, and thresholds were established based on regression with T-score values. Moreover, trabecular vBMD and cortical sBMD Z-scores were calculated to evaluate potential compartmental imbalances. ResultsThe derived aBMD curves closely aligned with the NHANES III Caucasian reference. Sex-specific thresholds for trabecular vBMD and cortical sBMD were established for patient stratification. Z-score comparisons revealed significant discrepancies between trabecular and cortical compartments in 52.0% of females and 48.7% of males, underlining the importance of compartment-specific bone assessment. ConclusionsThis study establishes reference curves for clinical interpretation of 3D-DXA parameters and demonstrates the potential of 3D-DXA to capture compartmental imbalances at the hip. Mini AbstractIn this study, hip scans from over a thousand men and women in Spain were analyzed to create normative reference values for 3D-DXA parameters. These results can help doctors better stratify people based on the status of each part of the bone and improve the management of osteoporosis.

IntroductionTrabecular bone exhibits brittle behaviour governed by microscale deformation and damage processes, yet quantitative characterisation of crack progression remains challenging because classical fracture mechanics approaches do not apply to architecturally discontinuous porous tissues. This study evaluates whether synchrotron X-ray computed tomography (XCT) combined with digital volume correlation (DVC) can provide a practical experimental approach for quantifying crack opening behaviour in human trabecular bone. MethodSemicylindrical specimens harvested from femoral heads of hip-fracture donors (n = 5) and non-fracture controls (n = 5) underwent stepwise three-point-bending during XCT imaging. Full-field displacement maps enabled direct measurement of crack mouth opening displacement (CMOD), crack length (a), and their ratio, CMOD/a, used here as a geometry-normalised comparative descriptor of brittle response. Automated crack segmentation using phase-congruency crack detection (PCCD) was compared against manual measurements. ResultsXCT-DVC successfully resolved three-dimensional displacement discontinuities during crack initiation and propagation in all specimens. Hip-fracture donors exhibited significantly lower critical crack-opening ratios (CMOD/a)* than Controls (0.31 vs 0.47; p = 0.008) and reached mechanical instability at lower applied loads, consistent with a more brittle structural response under this test configuration. Despite these differences, total crack extension ({Delta}a*) was similar between groups. Automated crack tracking using phase-congruency-based segmentation showed excellent agreement with manual measurements (r{superscript 2} = 0.98), confirming reliable extraction of crack geometry from DVC displacement fields. DiscussionThese results indicate that XCT-DVC can provide a practical approach for quantifying crack-opening behaviour in trabecular bone when classical fracture-mechanics parameters are not applicable in anatomically constrained specimens. The reduced critical crack-opening ratios and earlier instability observed in Hip-fracture donors are consistent with a more brittle comparative mechanical response that is not captured by crack extension alone. The strong agreement between automated and manual crack measurements further supports displacement-based descriptors as reliable comparative indicators of brittle behaviour in porous, architecturally discontinuous tissues. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=76 SRC="FIGDIR/small/714043v1_ufig1.gif" ALT="Figure 1"> View larger version (28K): org.highwire.dtl.DTLVardef@31c5d7org.highwire.dtl.DTLVardef@1b3d9a4org.highwire.dtl.DTLVardef@95df7borg.highwire.dtl.DTLVardef@1834216_HPS_FORMAT_FIGEXP M_FIG O_FLOATNOGraphical abstractC_FLOATNO C_FIG

Background: A previous meta-analysis by Singh-Ospina et al. (2017) suggested that Gender affirming hormone treatment (GAHT) does not change transgender mens bone mineral density (BMD) at any clinically relevant site; emerging studies and advances in synthesis methods necessitate an updated evaluation. The primary aim was to update the bone measures of Singh-Ospina et al. (2017), with the secondary aim to expand measures to how GAHT affects musculoskeletal health. Methods: A systematic review with meta-analysis was conducted using studies published in English up to 31 July 2024, identified through three electronic databases (PubMed, Embase, SportDiscus), and final cross-referencing in summer 2025. Primary outcomes were longitudinal changes in femoral neck (FN), lumbar spine (LS), and total hip (TH) bone mineral density (BMD). Secondary outcomes included body composition and muscle strength. Standardised effect sizes (Hedges g) were pooled using the inverse heterogeneity (IVhet) model. Results: GAHT (4 years) was not associated with significant longitudinal changes in FN, LS, or TH BMD. In contrast, substantial anabolic effects were observed, including increases in BMI (g = 0.13), body mass (g = 0.18), fat-free mass (g = 0.59), and muscle strength (g = 0.86). Heterogeneity was high for muscle strength, FN and TH BMD, limiting confidence in pooled estimates. Conversely, changes in LS BMD, BMI, body mass and fat-free mass demonstrated low heterogeneity and greater consistency across studies. Conclusion: Masculinising GAHT does not negatively affect clinically relevant BMD sites while reliably increasing lean mass and muscle strength; however, the evidence base remains methodologically weak and highly variable, particularly for FN and TH. The need for continued clinical monitoring of bone health and muscle function, alongside high-quality longitudinal research incorporating advanced imaging modalities such as HR pQCT is emphasised. Strengthening the evidence base will be essential for clarifying long-term skeletal trajectories as transgender men age. PROSPERO registration: CRD42024573102

RNF43 and ZNRF3 are transmembrane E3 ubiquitin ligases that negatively regulate Wnt signaling by promoting ubiquitination and degradation of Frizzled receptors. Loss of either gene enhances Wnt/{beta}-catenin signaling and has been linked to tumorigenesis. Wnt signaling is a key regulator of skeletal development and bone homeostasis, and pharmacologic activation of this pathway is an established therapy for osteoporosis. In Xenopus laevis, simultaneous disruption of rnf43 and znrf3 results in supernumerary limb formation; however, their roles in mammalian limb development and skeletal maintenance remain unclear. We demonstrate that mice homozygous for null alleles of both Rnf43 and Znrf3 do not develop supernumerary limbs. Because activation of Wnt/{beta}-catenin signaling in osteoblasts increases bone mass, we hypothesized that osteoblast-specific deletion of Rnf43 and/or Znrf3 would produce a high-bone-mass phenotype. Instead, osteoblast-specific loss of Znrf3 resulted in age-and sex-dependent reductions in trabecular bone mass, characterized by decreased bone mineral density and bone volume fraction, reduced trabecular number, and increased trabecular separation. Cortical bone exhibited increased cross-sectional size with reduced cortical area fraction and altered structural properties, while tissue mineral density was unchanged. In contrast, deletion of Rnf43 had minimal skeletal effects, and combined deletion of both genes did not exacerbate the phenotype observed with loss of Znrf3 alone. These findings identify Znrf3 as the dominant functional paralog regulating bone architecture in mature osteoblasts and underscore the importance of evaluating skeletal geometry when modulating upstream Wnt regulators.

Osteoblasts generate bone by secreting collagen and mineralizing it in response to various signaling cues. We have previously shown that the majority of ATP generated by differentiated osteoblasts in response to glucose is through glycolysis in contrast to undifferentiated cells that are more dependent on oxidative phosphorylation. To confirm our previous findings, metabolomics was performed for unlabeled polar metabolites, revealing elevated glycolytic metabolites at the later stages of differentiation. Krebs cycle (TCA cycle) metabolites were also changed confirming metabolic rerouting with differentiation. We hypothesized that an increase in mitophagy shifts ATP generation towards glycolysis resulting in the observed bioenergetic and metabolic changes. Utilizing calvarial osteoblasts isolated from a mitophagy reporter mouse model (MitoQC), an increase in mitophagy and the mitophagy receptor, Bnip3, was observed with osteoblast differentiation. KD of Bnip3 in osteoblasts inhibited differentiation and mineralization arising from impaired mitochondrial function. In vivo, male Bnip3 null mice exhibited a significant decrease in osteoblast numbers resulting in lower bone mass. Mechanistically we identified decreased fusion and increased fission factors, impaired stress signaling and increased proapoptotic factors in the absence of Bnip3. These data demonstrate for the first time that BNIP3 expression and mitophagy during osteoblast differentiation are necessary for relieving mitochondrial stress to maintain optimal bone mass.

Purpose To determine whether clinically significant weight loss (>5% of body weight) is associated with slower 2-year knee cartilage degeneration in individuals with and without radiographic osteoarthritis. This study used a cartilage structural assessment score derived from the spatial distribution of cartilage thickness, referred to as the Cartilage Thickness Score (CTh-Score). It is based on cartilage thickness patterns and scores the cartilage between 0 and 100, with higher scores indicating greater severity. Methods We conducted a retrospective matched cohort study within the Osteoarthritis Initiative. High-resolution cartilage thickness maps (CTh-Maps), along with their corresponding CTh-Score, were extracted from a public repository. Participants with complete radiographic and MRI data at baseline and 24 months were stratified by baseline Kellgren-Lawrence (KL) grade into non-radiographic OA (non-ROA; KL<2) and radiographic OA (ROA; KL>=2). Within strata, cases (>5% 2-year weight loss) were propensity score-matched 1:2 to weight-stable controls on age, sex, height, weight, KL grade, joint space width (JSW), KOOS Pain, baseline CTh-Score, and mean cartilage thickness in the medial and lateral femoral and tibial compartments. The primary outcome was 2-year change (delta) in CTh-Score, where higher values indicate worsening. Secondary outcomes were delta JSW, delta regional mean cartilage thickness, and delta KOOS Pain. Non-parametric tests were used. Results We included 164 cases and 328 controls in non-ROA, and 266 cases and 532 controls in ROA. Median (interquartile range) weight loss was -6.10 kg (-8.90, -4.70) versus +0.30 kg (-1.30, 2.20) in non-ROA and -6.80 kg (-9.10, -5.02) versus +0.40 kg (-1.40, 2.82) in ROA (both p<0.001). Weight loss was associated with significantly smaller 2-year increases in CTh-Score: in non-ROA, median 1.58 (0.61, 6.53) vs 3.14 (0.44, 7.12) (p=0.005); in ROA, median 1.69 (0.97, 6.71) vs 2.90 (0.19, 7.38) (p=0.004). No between-group differences were detected for delta JSW or delta regional mean cartilage thickness in any of the 4 ROIs. A trend toward greater KOOS Pain improvement with weight loss was observed in ROA: 2.75 (-3.35, 13.40) vs 0.00 (-5.60, 8.40) (p=0.06). Conclusions Achieving >5% weight loss over 2 years is associated with approximately 50% lower progression in median cartilage degeneration, as assessed by CTh-Score, in both non-ROA and ROA. No change was observed with conventional structural metrics. These findings support weight management as a structural disease-modifying strategy and highlight CTh-Score as a sensitive endpoint.

Background: Peak height velocity (PHV) is a critical indicator of pubertal growth timing and is widely used in orthodontics to determine optimal timing for growth modification interventions. Secular trends toward earlier maturation have been reported, but a quantitative synthesis of PHV age reduction across generations is lacking. Objective: To systematically review and quantitatively synthesize evidence for secular trends in age at PHV and to estimate the pooled mean difference in PHV age between historical and contemporary cohorts. Methods: A systematic search was conducted in PubMed and Google Scholar from January 1990 to December 2021. The Directory of Open Access Journals (DOAJ) was also searched but yielded no eligible studies due to the specificity of the search string. Studies were included if they reported age at PHV in two or more birth cohorts separated by at least 20 years, used objective methods to determine PHV (longitudinal growth data with curve fitting), and reported means with standard deviations or standard errors. Risk of bias was assessed using the Newcastle-Ottawa Scale. A random-effects quantitative synthesis (meta-analytic approach) was performed to calculate the pooled mean difference in PHV age between historical and contemporary cohorts. Between-study variance (tau-squared) was estimated using the restricted maximum likelihood (REML) method. Heterogeneity was assessed using I-squared statistics. Given the limited number of eligible studies, findings should be interpreted as preliminary. Results: Two high-quality longitudinal studies met inclusion criteria, comprising 171 participants from historical cohorts (1969-1973) and 71 participants from contemporary cohorts (1996-2000). The pooled mean difference in PHV age was -0.48 years (95% CI: -0.72 to -0.24, P < 0.001), indicating that contemporary children reach PHV approximately 0.5 years earlier than their historical counterparts. PHV velocity showed a pooled increase of 0.71 cm/year (95% CI: 0.48 to 0.94, P < 0.001). Heterogeneity was low (I-squared = 0% for both analyses). Both studies were rated as low risk of bias. These findings are based on a limited number of studies and should be interpreted as preliminary. Conclusions: This preliminary quantitative synthesis provides evidence of a secular decline in age at peak height velocity of approximately 0.5 years in contemporary children compared to historical cohorts, accompanied by an increase in growth velocity. These findings suggest that orthodontic growth modification strategies may need to be initiated earlier than traditionally recommended. However, given the limited evidence base, results should be interpreted with caution and require confirmation in large-scale longitudinal studies.

Study design A retrospective case control study Objective To predict proximal junctional kyphosis (PJK) risk by normalizing individual vertebral bone strength using the ratio of vertebral Hounsfield unit (HU) values around the upper instrumented vertebrae (UIV). Summary of background data PJK poses a significant challenge in treating patients after adult spinal deformity (ASD) surgery. While the vertebral body HU value is associated with PJK risk, the optimal threshold remains unclear, and a relative assessment of HU values within individuals has not been conducted. Methods Data on patients who underwent corrective fusion of the middle to lower thoracic region of the pelvis for ASD were assessed. The 126 patients were categorized into PJK and non-PJK groups. We compared the patients' backgrounds, vertebral body HU, and junctional HU ratio, defined as the HU value of UIV+1 divided by the HU value of UIV (HUUIV+1/HUUIV). The UIV+2/UIV+1 HU ratio was calculated similarly. Results The PJK and non-PJK groups included 30 and 96 patients, respectively. After propensity score matching, 28 patients from each group were analyzed. HU values at UIV+2 and UIV+1 (117.0 {+/-} 46.6 vs 145.1 {+/-} 45.9, p=0.018, and 105.5 {+/-} 36.2 vs 147.3 {+/-} 44.9, p<0.001, respectively) were lower in the PJK group. Junctional HU ratio was significantly lower in the PJK group (0.88 {+/-} 0.18 vs 1.13 {+/-} 0.25, p<0.001), and receiver operating characteristic analysis showed that the junctional HU ratio had the highest discriminative ability (area under the curve 0.812). At the optimal cutoff value (HU ratio of 0.905), the sensitivity and specificity for PJK were 64.3% and 89.3%, respectively. Conclusions A low junctional HU ratio was strongly associated with PJK after ASD surgery. This parameter reflects the bone strength mismatch at the proximal junction and may help improve preoperative risk assessment and UIV selection.

Objective To evaluate the Cartilage Thickness Score (CTh-Score) as a quantitative measure of cartilage damage severity by assessing its association with three osteoarthritis (OA) milestones and comparing its performance with conventional morphometric measures (radiographic minimum joint space width (JSW) and regional average cartilage thickness). Methods Data were obtained from the Osteoarthritis Initiative (OAI) and the publicly available OAI CTh-Maps and CTh-Score dataset. Three matched case-control designs were used to represent major OA milestones: (i) incident radiographic OA onset, (ii) combined pain and structural progression, and (iii) knee replacement (KR) in the coming 2 years. Progression subjects were extracted from the FNIH Biomarkers Consortium cohort. Cases and controls were compared at 4 years (T-4Y), 2 years (T-2Y), and 0 years (T0) before the milestone. MRI-based CTh-Score and regional average cartilage thickness, as well as JSW, were analyzed cross-sectionally and longitudinally. Associations with case status were assessed using adjusted logistic regression models, and responsiveness was evaluated using longitudinal change and standardized response means. Results The onset cohort included 307 matched case-control pairs, the progression cohort 164 cases and 369 controls, and the KR cohort 81 cases and 324 controls. Across all three study designs, the CTh-Score significantly differentiated cases from controls at all timepoints. In the onset cohort, the CTh-Score was higher in future cases than controls at T-4Y (16.2 vs 12.6, p=0.007), T-2Y (23.5 vs 16.7, p<0.001), and T0 (39.8 vs 18.6, p<0.001), whereas JSW and regional thickness measures showed limited or later discrimination. Similar findings were observed for progression (43.2 vs 33.0 at T-4Y; p<0.001) and KR (55.4 vs 46.1 at T-4Y; p=0.02) cohorts. Longitudinally, CTh-Score changes differentiated cases from controls earlier and more consistently than JSW or regional average thickness, and its responsiveness was consistently the highest across OA milestones and time intervals. In adjusted models, the CTh-Score was independently associated with all outcomes at T-4Y and T-2Y, with odds ratios per standard deviation increase ranging from 1.3 to 2.2. Conclusion The CTh-Score captures high-resolution cartilage thickness patterns associated with OA onset, progression, and future knee replacement, outperforming conventional morphometric measures in early discrimination, responsiveness, and predictive association. These findings support CTh-Score as a sensitive quantitative marker of cartilage damage severity across the OA continuum.

Objective: To evaluate the feasibility of nurse-led ultrasound hip screening for newborns and infants during home visits, focusing on whether trained public health nurses (PHNs) can obtain interpretable images for orthopedic pediatric surgeons' diagnosis, imaging error patterns, immediate operational challenges, and follow-up results of infants with suspected developmental dysplasia of the hip (DDH). Design: Pilot prospective cohort study. Sample: Forty-two infants were screened. PHNs conducted ultrasound hip screenings during home visits. Measurements: Diagnostically interpretable images, as determined by two pediatric orthopedic surgeons. Results: Diagnostically interpretable images of 75/84 (89.3%) hips were obtained. Surgeons identified three error patterns: incomplete visualization of the ilium (n = 2), joint capsule (n = 1), or bony roof (n = 2). Infant crying was an operational challenge (n = 1). Thirty-three (78.6%) hips were normal, four (9.5%) had abnormal findings requiring abduction exercises, three (7.1%) were referred to a hospital, and two (4.8%) failed imaging. One hip was diagnosed with subluxation, which went undetected by physical or risk screening. Conclusion: Nurse-led ultrasound hip screening for newborns and infants during home visits is feasible and may aid in early DDH detection. Further studies should assess diagnostic accuracy, cost-effectiveness, and long-term outcomes.

BackgroundPatient radiation exposure in diagnostic radiology is an important concern for radiation protection and patient safety. Monitoring radiation dose levels during radiographic examinations is essential to ensure compliance with diagnostic reference levels (DRLs) and to optimize radiological practices. ObjectiveThe aim of this study was to evaluate patient radiation dose during conventional lumbar spine radiography and compare the obtained values with diagnostic reference levels. MethodsA descriptive cross-sectional multicenter study was conducted in four hospitals in the Sous Massa region, Morocco, between April and June 2017. Data were collected from 142 patients undergoing lumbar spine radiography examinations and from 20 radiology technicians. Exposure parameters including tube voltage, tube current, exposure time, focus-to-film distance, and field size were recorded. Entrance surface dose (ESD) was estimated using MICADO software, and dose area product (DAP) values were subsequently calculated. The 75th percentile values were determined and compared with diagnostic reference levels. ResultsThe regional 75th percentile ESD values were 5.33 mGy for the anteroposterior projection and 7.38 mGy for the lateral projection. Corresponding DAP values were 1840.9 mGy.cm2 and 2783.65 mGy.cm2, respectively. All obtained values were below the diagnostic reference levels used for comparison. However, variations between hospitals were observed, likely due to differences in imaging protocols and equipment. ConclusionRadiation doses associated with lumbar spine radiography in the evaluated hospitals were within acceptable limits according to diagnostic reference levels. Continuous monitoring of patient radiation exposure and optimization of radiographic techniques remain essential to ensure effective radiation protection.

The bone matrix is precisely maintained and optimized to resist fractures. However, aging and disease deteriorate the bone matrix and increase fragility. Individuals with type 2 diabetes (T2D) have an elevated risk of bone fracture despite apparently normal bone mass. The chronic hyperglycemia in T2D promotes the formation of advanced glycation end-products (AGEs) in the bone tissue and modify the matrix mechanics. AGEs also bind to its receptor, RAGE, to activate inflammation and alter homeostasis. Using a leptin-receptor deficient mouse model of diabetes, we used a combination of high-resolution methods across multiple scales to evaluate the microarchitectural-, material- and cellular- level changes affected by the modulation of RAGE. To demonstrate the relevance of RAGE, we genetically ablated RAGE (RAGE-null) before the onset of diabetes; and to demonstrate the potency of RAGE as a disease modifying therapy, a RAGE antagonist (FPS-ZM1) was administered after prolonged diabetes. Diabetes impaired bone microstructure, the homeostatic actions of bone cells, the bone matrix nanomechanics, and whole- bone strength. The constitutive ablation of RAGE in diabetic animals prevented AGEs accumulation and the decline of trabecular connectivity; protected against the loss of osteocyte lacunae density and morphology; and maintained the matrix nanomechanics and bone strength. The inhibition of RAGE after the onset of diabetes reversed AGE accumulation and loss of bone volume; rescued osteocyte lacunae density and osteoclast activity; and restored matrix nanomechanics and bone strength. These results suggest that RAGE is a viable therapeutic target for diabetes-mediated impairments of bone quality. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=129 SRC="FIGDIR/small/716153v1_ufig1.gif" ALT="Figure 1000"> View larger version (26K): org.highwire.dtl.DTLVardef@6d8a5borg.highwire.dtl.DTLVardef@1967b5borg.highwire.dtl.DTLVardef@7cb1edorg.highwire.dtl.DTLVardef@85491d_HPS_FORMAT_FIGEXP M_FIG C_FIG

Domestic dogs (Canis familiaris) display more morphological variation than any other mammal. Cranial morphology has been extensively studied, as have the relationships with function, development, genetics, veterinary medicine, and breed welfare. Postcrania remain comparatively understudied, despite well-documented breed-specific predispositions to musculoskeletal disease. Here, we apply three-dimensional landmark-free morphometrics to quantify the shape of 743 elements from 213 dogs, including the scapula, humerus, radius, ulna, pelvic girdle, femur, tibia, and fibula. We assess integration among limb elements and investigate drivers of shape variation within and between breeds. Across most breeds, limb bone shape is strikingly similar. Dachshunds, however, exhibit distinct morphology across all elements and one to two orders of magnitude greater variation than any other breed. Despite this disparity, integration remains high between all element pairs. Remarkably, we find no significant relationship between bone shape and body mass, age, or pathology, but comparison with historic specimens reveals marked changes in dachshund long bone shape over the past [~]150 years. These extreme differences are not shared by other sampled chondrodysplastic breeds, underscoring the need to understand morphological diversity beyond simple categorisation. These findings provide a quantitative framework for linking postcranial morphology with function, disease risk, and evidence-based improvements to canine welfare.

Background: Diagnosing the over 700 known rare bone diseases (RBDs) is inherently challenging and often requires extensive time and multiple clinical visits. Effective treatment, particularly for RBDs with approved therapies, depends on early and precise identification of the specific RBD type. Image recognition artificial intelligence (AI) has the potential to significantly enhance diagnostic processes and improve patient outcomes. Many of these disorders cause characteristic skeletal changes, especially in the hands, and are associated with growth abnormalities. Consequently, affected children routinely undergo hand radiographs for bone age assessment, making these images a widely available yet underutilized diagnostic resource. Materials and Methods: We retrospectively compiled 5,623 multi-institutional hand radiographs from 2,471 patients with 45 different RBDs and 1,382 unaffected controls. We trained two deep learning models: a binary classifier to differentiate between RBD and non-RBD hand radiographs, and a multi-class classifier covering ten RBDs (or RBD groups), using 5-fold cross-validation. Preprocessing included masking, normalization, and data augmentation. Additionally, we applied occlusion sensitivity mapping to visualize class-specific features and evaluated the learned representations through cosine-based retrieval and UMAP projections of the feature space. Results: The affected versus unaffected classifier achieved a balanced accuracy of 85.5% on the test dataset. The ten-class classifier reached a balanced (top-1) accuracy of 76.6%, with top-3 accuracy exceeding 90%. Disorders with highly distinctive phenotypes, such as achondroplasia, achieved accuracies above 95%, whereas phenotypically overlapping disorders, such as ACAN- and SHOX-related short stature, were more frequently confused. Feature space analysis showed that validation samples clustered closely with their respective training distributions, supporting the consistency and generalizability of the learned embeddings. Conclusion: This manuscript presents a proof of principle for the development of Bone2Gene, a next-generation phenotyping (NGP) tool for the detection and differential diagnosis of RBDs, currently based on hand radiographs. Ongoing efforts focus on expanding the dataset to include additional RBDs or RBD groups in the current multi-class classifier for differential diagnosis and to further evaluate its generalizability. The Bone2Gene study is open to collaboration.

Objectives: Knee osteoarthritis (KOA) is a leading cause of disability, yet which patients will experience structural decline remains unclear. Body mass index (BMI) and lower limb alignment are established risk factors for KOA, but their independent and interactive effects on compartment-specific cartilage loss and total knee replacement (TKR) have not been characterized at scale. Methods: We analyzed 5,832 limbs from 3,016 participants in the Osteoarthritis Initiative followed over 7 years. Cartilage thickness in the weight-bearing medial and lateral femur and tibia was quantified, and lower limb alignment was measured using hip-knee-ankle (HKA) angle obtained from full-limb radiographs. Linear mixed-effects models estimated the independent and interactive effects of BMI and lower limb alignment on longitudinal cartilage thinning, and mixed-effects logistic regression modeled TKR risk. Results: In the medial compartment, BMI and varus alignment interacted multiplicatively, with their combined effect exceeding the sum of independent contributions (femur: p = 0.011; tibia: p < 0.001). At +10 kg/m2 BMI and +10 degrees varus, the rate of medial femur cartilage thinning was 243.5% faster than the reference rate. In the lateral compartment, BMI and valgus alignment were independently associated with faster cartilage thinning, with no significant interaction. TKR risk increased exponentially with HKA deviation (odds ratio [OR] = 1.38 per 1 degree; ~five-fold at 5 degrees malalignment) but was not associated with BMI. Conclusion: BMI and lower limb alignment influence structural KOA progression through compartment-specific pathways. The multiplicative interaction in the medial compartment identifies high BMI combined with varus malalignment as a discrete high-risk phenotype, with implications for clinical risk stratification and disease-modifying intervention design.

Keratin-positive giant cell-rich tumor (KPGCT) is a newly described bone and soft tissue tumor. The tumor is characterized by scattered keratin-positive cells and the presence of HMGA2::NCOR2 fusions. It is not known if the HMGA2::NCOR2 fusion is located in the keratin-positive cells, and little is known about how KPGCT develops. KPGCT shares some histologic features with tenosynovial giant cell tumor (TGCT), a soft tissue tumor with CSF1 rearrangements. Single-nuclei RNA sequencing (snRNA-seq) and Xenium spatial transcriptomics were used to elucidate the mechanisms driving KPGCT and compare KPGCT to TGCT. We show that the neoplastic cells in KPGCT constitute only a minority of cells in the tumor, and that they co-express keratin, HMGA2 and CSF1. The neoplastic cells in KPGCT express no synovial markers, confirming KPGCT as a distinct entity, separate from TGCT. The bulk of the tumor consists of CSF1R-expressing macrophages and osteoclast-like giant cells, suggesting an important role for CSF1-CSF1R signaling. In addition, we find that the cells with the HMGA2 translocation show activation of the hippo signaling pathway, which is known to regulate CSF1 expression. We show that the CSF1-CSF1R axis, possibly regulated through the hippo signaling pathway, plays an important role in KPGCT. This axis likely stimulates the migration and proliferation of macrophages, which form the majority of cells in the tumor, as well as their differentiation into osteoclasts-like giant cells. These results provide a rationale for the use of CSF1R inhibitors, which have already shown efficacy in TGCT, as a therapy for KPGCT. SignificanceKeratin-positive giant cell-rich tumor (KPGCT) is a rare, newly described soft tissue and bone tumor. By examining this tumour on a single-cell level, we confirm the identity of the neoplastic cells on a molecular level, showing these form a minority of cells in the tumor. We show that activation of the hippo pathway in the neoplastic cells is a likely driver of tumorigenesis. Additionally, we show the neoplastic cells produce large amounts of CSF1, attracting the macrophages that form the majority of cells in the tumor. This finding gives supporting evidence for anecdotal reports of response to CSF1 inhibitor therapy. Finally, we identify key differences between KPGCT and tenosynovial giant cell tumor, a tumor that shares histological features with KPGCT.

The osteochondral junction is a specialized region ensuring the biomechanical and biological integration of the unmineralized articular cartilage with the subchondral bone through an intermediate layer of mineralized cartilage. This location is of clinical relevance, being the target of osteoarthritis. While aging is considered a risk factor for osteoarthritis, the interplay between microstructural and material changes during aging and predisposing to joint degeneration is not fully clear. This is especially true for mineralized cartilage, which remains understudied despite its critical role in load transfer from unmineralized articular cartilage to bone. We investigate age-related alterations of mineralized cartilage and subchondral bone in rat tibiae of adult and aged animals using a multimodal, high-resolution, correlative analysis. Our approach includes micro-computed tomography to measure microstructural features, second harmonic generation imaging to visualize collagen organization, quantitative backscattered electron imaging to map local mineral content, and nanoindentation to obtain mechanical properties. Mineralized cartilage and subchondral bone exhibited distinct age-related modifications. At the architectural level, the subchondral plate thickened and the trabecular network became coarser, those changes being different from those observed in the metaphysis. At the tissue level, mineralized cartilage was less mineralized than bone but exhibits a greater relative increase of mineral content with age, underlying differences in mineralization. A central observation is that aging led to an abrupt transition in mineral content and mechanical properties across the interface between unmineralized and mineralized cartilage, with a conceivable impact on stress localization. Overall, these changes may alter load transfer and contribute to age-related joint degeneration.

Objectives. The association between skeletal muscle gene expression and knee osteoarthritis (OA) was examined among older adult participants of the Study of Muscle, Mobility and Aging (SOMMA). Methods. Inclusion criteria included knee radiographs and bulk RNA sequencing (RNAseq) in vastus lateralis muscle, resulting in 523 participants (56% female). Radiographic knee OA was determined by Kellgren-Lawrence (KL) grades. Differential gene expression was analyzed using a control group (KL [&le;] 1, n = 326) and two nested case groups: (a) KL [&ge;] 2 (n = 197), (b) KL [&ge;] 3 (n = 112). Results. Compared with controls, there were 27 and 41 genes associated (FDR [&le;] 0.05) with KL [&ge;] 2 and KL [&ge;] 3, respectively, and 16 genes significantly associated in both contrasts. For 15 of the 16 genes, the association magnitude was larger with more severe OA (KL [&ge;] 3). Genes associated in both contrasts included brain-derived neurotrophic factor (BDNF) and interferon regulatory factor-2 (IRF2). Gene sets enriched in KL [&ge;] 2 and KL [&ge;] 3 contrasts included DNA repair and branched chain amino acid (BCAA) catabolism. Conclusions. Our results in older adult SOMMA participants indicate that knee OA is associated with genes and pathways expressed in skeletal muscle that are involved in pain sensitization, BCAA catabolism, muscle function preservation, calcium transport and storage, inflammation, and extracellular matrix remodeling. Additional longitudinal studies will be needed to determine how these genes could affect the progression of knee OA.

ObjectiveTo investigate the effects of different gait patterns on knee joint biomechanics and dynamic stability during stair ascent. MethodsFourteen healthy males were recruited to ascend stairs using two distinct gait patterns: the "single-step" (leading with the same leg) and "cross-step" (alternating legs) strategies. Kinematic and kinetic data were collected synchronously using a Qualisys infrared motion capture system and a Kistler 3D force plate. Dynamic stability was quantified using the Margin of Stability (MOS), and knee joint biomechanics were evaluated using Patellofemoral Joint Stress (PFJS) and other relevant metrics. ResultsThroughout the gait cycle, there was no significant difference in the Medio-Lateral (ML) MOS between the single-step and cross-step patterns (P=0.318). However, in the Anterior-Posterior (AP) direction, the MOS for both patterns remained negative and decreased over time, with the cross-step pattern exhibiting significantly lower AP MOS values than the single-step pattern (P=0.002). At the moment of left foot-off, significant differences were observed in the right knee joint angle, right knee joint moment, net joint moment, effective quadriceps muscle lever arm, Quadriceps Force (QF), the angle between the quadriceps tendon and patellar ligament, Patellofemoral Joint Force (PFJF), patellofemoral joint stress, and patellofemoral contact area (all P<0.001). ConclusionsDuring stair ascent, the cross-step pattern reduces body stability, thereby increasing the risk of backward falls. Furthermore, this pattern increases patellofemoral joint stress, subjecting the knee to greater loading. Therefore, it is recommended to enhance lower limb muscle strength through targeted training to reduce fall risk. Additionally, adopting a more cautious gait strategy (such as the single-step pattern) can help minimize patellofemoral joint loading and mitigate the risk of patellofemoral pain.

The anatomical localization of lymphatic vessels in bone remains controversial and has led to conflicting interpretations of skeletal lymphatic function. Here we assessed lymphatic identity and localization in bone using mouse genetic labeling, tissue clearance, and three-dimensional imaging. We analyzed long bones after extensive periosteum removal and identified Vegfr3+ blood vessels lacking Lyve1 expression within bone marrow, whereas Vegfr3+Lyve1+ lymphatic vessels were confined to residual periosteal regions. Genetic lineage tracing using Prox1-Cre/ER;mScarlet further confirmed that lymphatic vessels are absent from long bone marrow and restricted to periosteal compartments, particularly in fibrous but not cambial layers. Extending these analyses to the mandible, we observed Vegfr3+Lyve1+ lymphatic vessels localized to periarticular soft tissues surrounding the temporomandibular joint (TMJ), while mandibular bone marrow contained only Vegfr3+Lyve1- blood vessels and lacked Prox1 lineage-traced lymphatic vessels. Together, these findings establish that lymphatic vessels in bone are confined to periosteal and periarticular compartments and absent from bone marrow, providing a framework for interpreting lymphatic contributions to skeletal physiology and disease.