Correlation between Bone Mineral Density and Progression of Hip Osteoarthritis in Adult Men and Women in Bulgaria—Results from a 7-Year Study

Changes in clinical presentation, radiographic progression (RP), bone mineral density (BMD), bone turnover (BT), and cartilage turnover (CT) markers were compared in two groups of patients with hip osteoarthritis (HOA) over a period of 7 years. Each group consisted of 150 patients, including a control group on standard-of-care therapy (SC) with simple analgesics and physical exercises, and a study group (SG) on standard-of-care therapy supplemented by vitamin D3 and intravenous administration of zoledronic acid (5 mg) yearly for 3 consecutive years. Patient groups were homogenized regarding the following: (1) radiographic grade (RG), including 75 patients with hip OA RG II according to the Kellgren–Lawrence grading system (K/L), and 75 with RG III on K/L; (2) radiographic model (RM), as each of the K/L grades was subdivided into three subgroups consisting of 25 patients of different RMs: atrophic (‘A’), intermediate (‘I’), and hypertrophic (‘H’); (3) gender-equal ratio of men and women in each subgroup (Female/Male = 15/10). The following parameters were assessed: (1) clinical parameters (CP), pain at walking (WP-VAS 100 mm), functional ability (WOMAC-C), and time to total hip replacement (tTHR); (2) radiographic indicators(RI)—joint space width (JSW) and speed of joint space narrowing (JSN), changes in BMD (DXA), including proximal femur (PF-BMD), lumbar spine (LS-BMD), and total body (TB-BMD); (3) laboratory parameters (LP)—vitamin D3 levels and levels of BT/CT markers. RV were assessed every 12 months, whereas CV/LV were assessed every 6 months. Results: Cross-sectional analysis (CsA) at baseline showed statistically significant differences (SSD) at p < 0.05 in CP (WP, WOMAC-C); BMD of all sites and levels of CT/BT markers between the ‘A’ and ‘H’ RM groups in all patients. Longitudinal analysis (LtA) showed SSD (p < 0.05) between CG and SG in all CP (WP, WOMAC-C, tTHR) parameters of RP (mJSW, JSN), BMD of all sites, and levels of CT/BT markers for all ‘A’ models and in 30% of ‘I’-RMs (those with elevated markers for BT/CT at baseline and during the observation period). Conclusion: The presence of SSD at baseline (‘A’ vs. ‘H’) supported the thesis that at least two different subgroups of HOA exist: one associated with ‘A’ and the other with ‘H’ models. D3 supplementation and the intravenous administration of bisphosphonate were the treatment strategies that slowed down RP and postponed tTHR by over 12 months in the ‘A’ and ‘I’ RM with elevated BT/CT markers.

Patients with the following characteristics were excluded from the study: (1) secondary HOA; (2) severe deviations in the weight-bearing axis; (3) presence of synovitis or effusion; (4) signs of rapidly progressing HOA in 'A'-RM; (5) any intraarticular treatment or treatment with sulfate sugars, biocollagen, hyaluronic acid, diacerein, or avocado and soybean unsaponifiables within the 6 months prior to the baseline visit; (6) age above 60 and below 70; (7) body mass index (BMI) less than 21 kg/m 2 or more than 28 kg/m 2 ; (8) poorly controlled internal diseases, including hypertension, diabetes, cardiovascular, and cerebrovascular disease.
The age restrictions were introduced in order to exclude mild hip dysplasia (clinical presentation before 55 years of age) [28], the effect of hormonal changes in the early menopause (50-55 years) over SB [34], and the high rate of HOA in people over 60 years of age [18,22], as well as due to data regarding life duration (decreasing the chance of successfully finishing the 7-year follow-up period). The restriction in BMI and severe deviations in the weight-bearing axis were introduced in order to eliminate the effect of these variables on the RP of HOA [26,27]. The need for a fixed gender distribution in the patient groups was discussed previously [29][30][31][32][33][34][35][36][37].

Study Design
The sample size of each patient group and subgroup was calculated using the methodology described by J. Wittes [44], considering a 2% possible loss of patients in the treatment groups, as compared to 1% in the control group. Restricted block randomization [45] was applied during the screening and selection process, resulting in the formation of the following six blocks: K/L-II'A'; K/L-II'I'; K/L-II'H'; K/L-III'A'; K/L-III'I'; and K/L-III'H'. From each block of 50 patients, subgroups with a fixed female to male ratio (15/10) were randomly generated using a computer program (Table 1).

Blinding
The rheumatologists performing the clinical follow-up were blinded to the group types and numbers, as well as the laboratory and radiographic (DXA and X-ray) follow-up results. The medication was administered by the hospital pharmacist, who knew the randomization group and number but was blinded to the clinical and radiological data. The radiologists were blinded to the design, randomization group and number, and clinical data of the patients. The study coordinators (external, non-medical personnel provided by the Bulgarian Drug Agency) were responsible for contact with the patients (calling for the visits) and knew which randomization number corresponded to which patient, but were blinded to the clinical and radiological data of the patients. Only the principal Life 2023, 13, 421 4 of 21 monitor (external medical personnel provided by the Bulgarian Drug Agency) had access to the randomization group numbers and all data from CP and RI, but they did not know the patients. A summary of the research design, including the patient groups, methods, follow-up, and analyses, is presented in Table 1.   tTHR-time to conversion to total hip replacement; D3-25-hydroxy vitamin D; s-CTX-I-serum-beta-isomerized carboxy-terminal cross-linking telopeptide of type I collagen; u-CTX-II-urine -C-terminal crosslinking telopeptides of collagen type II; mJSW-mean joint space width; JSN-joint space narrowing; BMD-bone mineral density; PF/LS/TB-proximal femur, lumbar spine, total body.

Physical Examination
Patient physical examinations and recording of the clinical data were performed by a board-certified rheumatologist. At each patient visit, the following factors were evaluated: vital signs, height, weight (BMI), walking pain (VAS-100 mm) [41], functional ability (WOMAC-C) [42], the presence of adverse events, and quality of life (SF-36 and PtGA) [46]. Treatment responses by OMERACT-OARSI set of responder criteria [47] and minimal clinically important improvement (MCII) [48] were also assessed.

Radiological Examinations
The radiographic images were taken in an upright weight-bearing position, anterior-posterior projection, with a slight (15-20 degrees) inner rotation of the feet, which was ensured by using a 'V '-shaped pad placed 100 cm away from the source, with a perpendicular ray, focused on 4 cm above the symphysis. The following parameters were assessed: RG on K/L grade [37]; RM-'H', 'I', 'A'; JSW-the measurement of each joint was performed manually using software for measuring distances in digital radiographic images at three points: superolateral, apical, and superomedial ( Figure 1). For the statistical analyses, the mean value of the three distances was used (mean joint space width-mJSW) posterior projection, with a slight (15-20 degrees) inner rotation of the feet, which wa ensured by using a 'V '-shaped pad placed 100 cm away from the source, with a perpen dicular ray, focused on 4 cm above the symphysis. The following parameters were as sessed: RG on К/L grade [37]; RM-'Н', 'I', 'А'; JSW-the measurement of each joint wa performed manually using software for measuring distances in digital radiographic im ages at three points: superolateral, apical, and superomedial ( Figure 1). For the statistica analyses, the mean value of the three distances was used (mean joint space width-mJSW and the annual speed of JSN-mm/year (JSNM12 = mJSW-M0 -mJSWM-12), according to the recommendations of the 2004 Barcelona Consensus Group [29].  [29], and in a patient from our study (right image).
The measurements and interpretation of the results from the radiographic and DXA investigations were conducted by two separate certified by ISCD radiologists, who wer blinded to the design and clinical and laboratory data, and with very good inter-reade reliability (intraclass correlation coefficient ICC of 0.918, 95% CI: 0.846-0.960) and PABAK (prevalence-adjusted and bias-adjusted kappa) values for X-ray/DXA reading of 0.860 and 0.880, respectively.  [29], and in a patient from our study (right image).
DXA measurements were conducted using Lunar Prodigy Primo-en CORE, version 17, according to the methodology of ISCD, including recommendations for calibration, measurement, and interpretation of results from 2015 [30]. The following parameters were assessed: PF-BMD; LS-BMD; TB-BMD (proximal femur, lumbar spine, total body); and parameters for bone geometry for HSA: HAL (hip axis length); NSA (neck shaft angle); CSA (cross-sectional area); CSMI (cross-sectional moment of inertia); MNW (minimal neck width); FN-CT (cortical thickness of the femoral neck); FS-CT (cortical thickness of femoral shaft); SM (section modulus); BR (buckling ratio) (Figures 2-4).   The measurements and interpretation of the results from the radiographic and DXA investigations were conducted by two separate certified by ISCD radiologists, who were blinded to the design and clinical and laboratory data, and with very good inter-reader reliability (intraclass correlation coefficient ICC of 0.918, 95% CI: 0.846-0.960) and PABAK  DXA measurements: dual femur-this scan was used to obtain the BMD parameters: oral neck (FN-BMD); femoral head (FH-BMD); total hip BMD (TH-BMD); parameters of bon ometry (HAL-hip axis length; NSA-neck shaft angle; MNW-minimal neck width; CSA-c sectional area; CSMI-cross-sectional moment of inertia; FN-CT-cortical thickness of femoral FS-CT-cortical thickness of femoral shaft; SM-section modules; BR-buckling ratio) and als the comparisons between the mentioned parameters-left/right (target/non-target) joint.    DXA measurements: dual femur-this scan was used to obtain the BMD parameters: oral neck (FN-BMD); femoral head (FH-BMD); total hip BMD (TH-BMD); parameters of bon ometry (HAL-hip axis length; NSA-neck shaft angle; MNW-minimal neck width; CSA-c sectional area; CSMI-cross-sectional moment of inertia; FN-CT-cortical thickness of femoral FS-CT-cortical thickness of femoral shaft; SM-section modules; BR-buckling ratio) and als the comparisons between the mentioned parameters-left/right (target/non-target) joint.   . DXA measurements: total body-this scan was used to obtain the BMD parameters from the different regions (head-BMD; arms-BMD; legs-BMD; trunk-BMD; ribs-BMD; spine-BMD; pelvis-BMD) and a total body score (TB-BMD).

Follow-Up
The following parameters were assessed for the patients in both groups (CG and SG): RV (JSW/JSN; BMD-DXA); CV-(WP; F -WOMAC-C; PtGA; tTHR); LV-levels of vitamin D3 and levels of BT/CT markers. RV were assessed every 12 months, whereas CV and LV were assessed every 6 months.

Statistical Analyses
The data were analyzed using Statistical Package for the Social Sciences (SPSS) version 21. The results were presented as numbers and percentages (%) for the qualitative variables and as medians with interquartile ranges (IQRs) for the continuous variables. For the detection of statistically significant differences (SSDs), Mann-Whitney and Kruskal-Wallis non-parametric tests were used. If SSD was present, post hoc multiple comparisons with Dunn-Bonferroni were performed. Pearson correlation analysis and multiple logistic regression were used to analyze the effect of BMD and bone geometry markers on radiographic progression. All statistical tests were two-tailed and performed at a level of significance (α) of 0.05 and a power of 80%. Exact p-values were used to interpret the results, with a p < 0.05 meaning a statistically significant change.

Results from the Cross-Sectional Study
There were no significant differences in the values of the clinical, laboratory, and radiographic parameters at baseline between the SG and CG groups in all three models ( Table 2).
All 'A' models from both RGs of the two groups had decreased BMD at all measurement sites (PF-BMD; LS-BMD; TB-BMD) and 70% of them met the ISCD criteria [30] for osteoporosis. On the other hand, all 'I' and 'H' models were with normal ('I') or slightly increased ('H') local BMD (p < 0.05), as compared to the controls without HOA ( Figure 5). , scull (Head), and total body (Total) for each RM. Z scores were used to compare to controls without HOA after standardization for gender, age, and BMI. *-a statistically significant difference (p<0.05) when comparing to controls without HOA after standardization for gender, age, and BMI.
All 'А' models had increased levels of CT or BT markers as compared to the 'I' model. The latter showed increased CT/BT markers as compared to the 'Н' model, p < 0.05 (Table 2 and Figure 6).
Biomarkers of cartilage degeneration (u-CTX-II) normotrophic and atrophic). The results are shown as Z scores of femoral neck (FN), lumbar spine (LS), scull (Head), and total body (Total) for each RM. Z scores were used to compare to controls without HOA after standardization for gender, age, and BMI. *-a statistically significant difference (p < 0.05) when comparing to controls without HOA after standardization for gender, age, and BMI.      All 'A' models had increased levels of CT or BT markers as compared to the 'I' model. The latter showed increased CT/BT markers as compared to the 'H' model, p < 0.05 (Table 2 and Figure 6). , and total body (Total) for each RM. Z scores were used to compare to controls without HOA after standardization for gender, age, and BMI. *-a statistically significant difference (p<0.05) when comparing to controls without HOA after standardization for gender, age, and BMI.
All 'А' models had increased levels of CT or BT markers as compared to the 'I' model. The latter showed increased CT/BT markers as compared to the 'Н' model, p < 0.05 (Table 2 and Figure 6). The comparison between 'А' and 'Н' models showed statistically significant differences in CP (p < 0.01), the levels of ВТ/СТ markers, and BMD and bone geometry markers (a wider femoral neck and increased FSI, p < 0.05), and no significant differences in mJSW. The comparisons of 'H' vs. 'I' and of 'I' vs. 'A' showed some differences, but they were not statistically significant (p > 0.05), even at the interim analyses (Table 3). The comparison between 'A' and 'H' models showed statistically significant differences in CP (p < 0.01), the levels of BT/CT markers, and BMD and bone geometry markers (a wider femoral neck and increased FSI, p < 0.05), and no significant differences in mJSW. The comparisons of 'H' vs. 'I' and of 'I' vs. 'A' showed some differences, but they were not statistically significant (p > 0.05), even at the interim analyses (Table 3).

Results from the Longitudinal Study
Serious adverse effects were not observed in either group, and there was no loss of patients during the 7-year follow-up period. Transient myalgia (2-6 h with spontaneous resolution) after the intravenous administration of ZA was observed in 46 patients (31%) from the SG group.
The loss of whole subgroups was observed at the following times: M36-CG-III'A'; M42-CG-III'I'/SG-III'A'; M54-CG-III'H' + II'A' /SG-II'H' + III'I'; M66-CG-II'I'/SG-II'A'; and M78-CG-II'H'/SG-II'H'+II'I' (Figure 7). This was an expected result of HOA's natural and therapy-modified evolution, rather than the loss of subjects. The final clinical and radiographic investigations were completed, with a clear path toward conversion to total hip replacement (THR) and the possibility of follow-up.  In the CG group, the natural evolution of HOA did not demonstrate unexpected results. While the CT and BT markers were constantly increasing, CP (WP and F) worsened over time, along with advances in morphological changes (mJSW) and decreases in BMD at all measurement visits ( Figure 8 and Table 4). In the SG group, a delay in the decrease in BMD followed by a transient increase in BMD was observed in all RMs from both RGs, albeit at different rates: at month 12 (M12) in the 'А' and 'I' models and at М36 in the 'Н' model. This increase was followed by a rapid decrease in BMD, starting from M48 (12 months after the last administration of ZA), at the same rate of change as in the SC group (Table 4). BMD changes were accompanied  In the CG group, the natural evolution of HOA did not demonstrate unexpected results. While the CT and BT markers were constantly increasing, CP (WP and F) worsened over time, along with advances in morphological changes (mJSW) and decreases in BMD at all measurement visits ( Figure 8 and Table 4).  In the CG group, the natural evolution of HOA did not demonstrate unexpected results. While the CT and BT markers were constantly increasing, CP (WP and F) worsened over time, along with advances in morphological changes (mJSW) and decreases in BMD at all measurement visits ( Figure 8 and Table 4). In the SG group, a delay in the decrease in BMD followed by a transient increase in BMD was observed in all RMs from both RGs, albeit at different rates: at month 12 (M12) in the 'А' and 'I' models and at М36 in the 'Н' model. This increase was followed by a rapid decrease in BMD, starting from M48 (12 months after the last administration of ZA),  Each group consisted of 25 subjects; CG-control group; SG-study group; WP-pain at walking; VAS-visual analogue scale -100 mm.; WOMAC-C-WOMAC-function scale; s-CTX-I-serum-beta-isomerized carboxyterminal cross-linking telopeptide of type I collagen in nanogram per milliliter; u-CTX-II-urine-C-terminal crosslinking telopeptides of collagen type II, presented as corrected concentrations of uCTX-II for urinary creatinine concentration in ng/mmol Cr. mJSW-mean joint space width. Values associated with the OARSI responses are marked with * and these with MCII responses, with # .

TH-BMD changes over 36 months (CG vs. SG) in the different radiological models
In the SG group, a delay in the decrease in BMD followed by a transient increase in BMD was observed in all RMs from both RGs, albeit at different rates: at month 12 (M12) in the 'A' and 'I' models and at M36 in the 'H' model. This increase was followed by a rapid decrease in BMD, starting from M48 (12 months after the last administration of ZA), at the same rate of change as in the SC group (Table 4). BMD changes were accompanied by similar changes in the levels of BT/CT markers, namely a decrease followed by an increase after M48 at the same rate as in the CG group. The described changes in CT/BT and BMD were accompanied by transient but statistically significant differences in the values of CP and RI (SG vs. CG).

Within Group Comparisons
At baseline (M0) in both patient groups (SC and SG), a significant difference existed between the endmost RMs ('H' vs. 'A') in CP (WP/F). No significant differences were found in mJSW, including the end RMs. During the follow-up, at points M12 (CP's/ JSN) and M36 (mJSW), significant differences were observed between all RMs of the SC group (Tables 5 and 6). and at M36 (K/L-III) ( Table 5). Regarding the radiological indicators, the comparisons in the SG showed no significant differences ('H' vs. 'I' vs. 'A') in the slow progression stage (K/L-II) until the end of the follow-up. In the rapidly progressing stage (K/L-III), the intragroup comparisons of 'H' vs. 'I' vs. 'A' showed no significant differences in JSN, despite the presence of differences in mJSW (M12/M36), suggesting a slowing down of the speed of JSN (Table 6). Therapeutic responses, according to the OMERACT-OARSI set of responder criteria [49] and Tubach F. (minimal clinically important improvement -MCII) [50] were registered in 'A' and 30% of 'I' RMs from both radiological grades, K/L-II and III (CG vs. SG) in CP (WP; WOMAC-C; PtGA). These responses were first detected at month six after the first application of ZA and persisted until the sixth (OMERACT-OARSI set) or twelfth month (MCII) after the last ZA administration (Table 4).

Between Group (SC vs. SG) Comparisons
The first occurrence of significant differences in JSN between the two groups (SG and CG) was detected at month 12 after the first application of ZA in the 'A' models of both RGs and in 30% of 'I' models, namely those with elevated markers for BT and CT at baseline and during the observation period. These differences were accompanied by analogous ones in mJSW with a later onset at M36 and persisted during all follow-up visits, until the end of the study (Table 7). Statistically significant differences between the study and control groups were observed in the 'A' and 'I' models from both RGs regarding the length of time before THR had to be performed. Due to the longer natural evolution of HOA in the K/L-II grade (K/L-II'A' vs. K/L-III'A') these differences were highly significant (p < 0.001) ( Table 8). By using multiple logistic regression, the observed changes in BMD (CG vs. SG) were identified as the main reason (Table 9) for the observed changes in CP and RP. The decrease in their values by one standard deviation (SD) was associated with an accelerated progression (OR = 6.561, p < 0.001; OR = 6.495, p < 0.001, respectively).

Discussion
The methodology of computer-simulated models, which is widely used today [49,50], had been not developed in detail at the time of the current study's planning (2011)(2012)(2013), resulting in the planning and execution of an experimental study instead of a simulation study.
The gender-related differences in the levels of CB and CT markers and radiological indicators (mJSW/BMD) in men and women with the same values of CP (WP/F), radiological grade, and model, are well-known and have been discussed in the literature [29][30][31][32][33][34]. However, these gender differences, as well as the differences in CP, levels of the CT-BT markers, and BMD between 'A' and 'H' models of the same RG, suggest that patient groups should be homogenized not only by RG, but also by RM, with a fixed ratio in the number of men and women in each group. Furthermore, the significant differences in CP, CT/BT markers, and BMD in the absence of corresponding changes in mJSW between 'A' and 'H' RMs of the same RG suggest that the value of the CP (WP/F) is dependent not only on mJSW (loss of AC), but also on the quality (changes) in the SB.
Differences in bone geometry, BMD, levels of BT/CT markers, and bone histomorphometry between the different RMs of HOA were reported in a number of studies [10,[24][25][26]. 'H' models had a wider femoral neck, increased BMD at all measurement points, and increased FSI compared to 'A' models, which had a narrower femoral neck, decreased BMD at all measurement points, and decreased FSI. The levels of BT/CT markers showed delayed or accelerated bone resorption in the 'H' as compared to 'A' models [24,36], whereas histomorphometric studies of specimens from THR surgeries showed smaller bone volumes and thinner trabeculae in 'A' models [25].
One of the objectives of our study was to examine the relation between RM and BMD in the Bulgarian population and the treatment options available to delay RP in patients with different types of RM of HOA through vitamin D3 supplementation and intravenous administration of ZA (which causes a change in BMD).
The results from the comparison between bone geometry parameters, BMD, and the levels of serum CTX-I and urine CTH-II at baseline in our study were similar to those reported in other similar studies [10,24]. The abovementioned results (ours and from other studies) suggest the presence of at least two different types of HOA, associated with 'A' and 'H' models.
At baseline, within-group comparisons in both groups (SC and SG) showed significant differences between the endmost RMs ('H' vs. 'A') in CP (WP/F) and an absence of significant differences in mJSW, even between the endmost models.
During the follow-up, the natural evolution of the SC group led to the appearance of significant differences at M12 in the clinical parameters and JSN and at M36 in mJSW among all RMs. The absence of such differences ('H' vs. 'I') and their later onset ('I' vs. 'A' at M60 for K/L-II and at M36 for K/L-III) in the SG group could be attributed to the therapy-modified (D3 + ZA) evolution of the HOA. These factors suggest that the changes in BMD due to the anti-resorptive therapy equalized the rate of the radiographic progression of the models with reduced BMD ('A') to those with a normal ('I') or slightly increased ('H') BMD.
At baseline, all 'A' models in both RGs in the two groups (SC and SG) had decreased BMD at all measurement sites (PF-BMD; LS-BMD; TB-BMD), all 'I' models had normal BMD, and all 'H' models had normal or slightly increased local BMD.
The analysis of the BMD change over the follow-up period demonstrated that the BMD of the control group decreased, reaching the normal, or lower than the normal, BMD of patients without HOA (at all measurement sites) as early as the second ('I') or the third year ('H'), which was accompanied by a decrease in vitamin D3 levels and an increase in serum CTX-I and urinary CTX-II. A possible explanation for the changes in vitamin D3, s-CTX-I, and u-CTX-II levels could be the decreased physical activity and lifestyle changes related to HOA progression. Nevertheless, the observed significant differences at baseline (CsA) in TB-BMD and head BMD (which are not related to the biomechanics in HOA), cannot be explained in this way, but rather by the presence of different types of HOA, as was discussed earlier.
A possible explanation for the changes in vitamin D3, s-CTX-I, and u-CTX-II levels is the decreased physical activity and lifestyle changes related to HOA progression. Nevertheless, the observed significant differences at baseline (CsA) in TB-BMD and head BMD, which are not related to the biomechanics in HOA, cannot be explained in this way, but rather by the presence of different types of HOA, as was discussed earlier.
The longitudinal significant differences in RP (JSN-M12/mJSW-M36) between the two groups (SC vs. SG), as well as in the tTHR between the 'A' and 'I' models of both groups, showed a clear relationship with the detected changes in BMD (FN-BMD; TH-BMD. This relation correlates by grade with the well-known impact of these parameters on the fracture risk. In our study, a decrease in FN-BMD and TH-BMD by 1 SD was associated with accelerated radiographic progression and increased fracture risk.
To our knowledge, no other studies have examined the effect of intravenous ZA on the progression of HOA [51]. The only study evaluating the effect of a bisphosphonate administration on the progression of HOA used oral alendronate [52]. In the present study, we observed a significant improvement in VAS and WOMAC scores and BMD, with a parallel decrease in the values of sCTX-I and uCTX-II, but no significant delay in the progression of HOA, defined as a decrease in JSW >0.30 mm or conversion to total hip arthroplasty [52].
At the same time, a number of studies have analyzed the efficacy of ZA administration in reducing postoperative complications after THR [53]. There have also been a large number of studies similar to ours but on the efficacy of ZA administration in knee OA (KOA) [54]. Those KOA studies showed contradictory results, which, in addition to the failure of alendronate to slow the radiographic progression of HOA [52], are most likely due to the inconsistency of the study designs used; the discussed influence of the radiographic model and gender on the clinical and radiographic variables.
Considering the rates of the different RMs of HOA in the general population (60% 'I'-RM; 30% 'H'-RM; 10% 'A'-RM) [23,24], and this study's results, which showed that the treatment was effective in all 'A' models and in 30% of 'I' models (those with persistently elevated BP/CT), it can be extrapolated that 30% of patients with HOA would benefit from a combined treatment of D3 with ZA. At first glance, this percentage does not seem so high, but considering the prevalence of HOA among the population over 60 years of age, 0.1-0.2% of people in the general population would benefit from the treatment strategy used in our study.
With regard to when an anti-resorptive therapy in patients with 'A' models of HOA should be started, our recommendation is as early as possible. At baseline, all 'A' models of HOA had DXA scores of osteopenia or osteoporosis (70%) at all measurement sites and increased levels of BT and CT markers.
All 'I' models with increased BT and CT markers had slightly reduced BMD at all measurement sites, as compared to the controls without HOA, without fulfilling the ISCD criteria for osteopenia [30]. However, within the two-year follow-up, these values changed and fulfilled the criteria for osteopenia in 100% of patients and for osteoporosis in 56% of them.
Our findings suggest that patients with 'I' models of HOA who have elevated BP and CT markers should be followed up on an annual basis. If a significant decrease in BMD is detected on two consecutive measurements, antiresorptive therapy should be started, even if the DXA scores do not meet the ISCD criteria for osteoporosis.
However, it is important to remember that all 'H' and 'I' models of HOA with normal BT and CT markers, which accounts for 70% of HOA patients, do not benefit significantly from anti-resorptive therapy.
Our study has several limitations that need to be addressed. One of them is the small number of bone (serum CTX-I) and cartilage (urinary -II) turnover markers, which could not adequately demonstrate the composite mechanism of the processes occurring at the osteochondral junction. Another limitation stems from the dual nature of the DXA investigation used for HSA, in the absence of parallel data from a three-dimensional imaging modality (QTC). Another area for improvement is the calculation of mJSW and JSN based on manual measurements of JSW at three points rather than MRI-based assessments.
In conclusion, the significant differences in CP, BMD, bone geometry parameters, and levels of BT and CT markers observed at baseline between the 'A' and 'H' models support the hypothesis of the existence of two distinct subgroups of HOA. The one associated with 'A' models and the other with 'H' models exhibit specific characteristics of ACbreakdown, SB-remodeling, and clinical and radiological progression, and a different response to anti-resorptive drugs. In 'A' and 'I' models with increased BT and CT markers, D3 supplementation and intravenous administration of ZA constitute a treatment strategy that slows RP and delays tTHR for more than a year.