Association of Bone Marrow Lesion Localisation with Weight Bearing Pain in People with Knee Osteoarthritis: Data From The Osteoarthritis Initiative

Background: Subchondral bone marrow lesions (BMLs) detected on MRI in knee osteoarthritis (OA) are associated with knee pain. The prevalence and progression of BMLs are increased by mechanical knee load. However, associations of BML location with weight-bearing knee pain are currently unknown. In this study, we aim to demonstrate associations of BML location with weight-bearing knee pain in knee OA. Methods: We analyzed 1412 and 582 varus knees from cross-sectional and longitudinal Osteoarthritis Initiative datasets, respectively. BML scores were analysed for 5 anatomical regions (median and lateral femorotibial, medial and lateral patellofemoral, and subspinous). Weight-bearing and non-weight-bearing pain scores were derived from WOMAC pain items. Correlation and negative binomial regression models were used for analysis of associations between the BML scores and pain at baseline, and changes in the BML scores and changes in pain after 24-month follow up. Results: Greater BML scores at medial femorotibial and lateral patellofemoral compartments were associated with total and weight-bearing pain scores, and statistical signicance was retained after adjusting for BML scores at the other 4 joint compartments and other OA features, as well as for non-weight-bearing pain, age, sex and Body Mass Index (BMI) (femorotibial; B=0.08 (95%CI, 0.01-0.09) p=0.02, patellofemoral; B=0.13 (95%CI, 0.03-0.23) p=0.01). BML scores were not signicantly associated with non-weight-bearing pain after adjustment for weightbearing pain and other demographic and OA features. Change over 24 months in BML score in the medial femorotibial compartment was signicantly associated with change in weight-bearing pain after adjusting for non-weight-bearing pain, age, sex, baseline weight-bearing pain, BMI, and BML at the other 4 joint compartments (B=0.10 (95%CI, 0.02-0.18) p=0.01). Conclusions: BML size at the medial

in total knee BML volume were associated with increased knee pain severity [3]. Greater BML scores were associated with weight-bearing pain and less so with non-weight-bearing pain [2].
The exact pathophysiology of BMLs is still under debate. Histologically, BMLs are characterized by bone marrow necrosis, trabecular abnormalities, bone marrow broses, edema, cellular in ltration and vascular proliferation [4,5]. Microarray analysis of subchondral BMLs in OA demonstrated upregulation of genes implicated in neurogenesis, osteochondral turnover and in ammation [5]. BMLs might contribute to OA pain by this generation of chemical factors, sensitising nerves within the subchondral bone. Sensitised subchondral nerves might be activated by biomechanical forces, for example while standing, walking or climbing stairs. If so, then BMLs in weight-bearing rather than non-weight-bearing regions would be expected to speci cally associate with weight-bearing pain. However, associations of BML location with weight-bearing knee pain are currently unknown.
OA structural change particularly affects weight-bearing components of the joint and mechanical knee load increases the prevalence and progression of BMLs [6]. Knee OA is most commonly associated with varus angulation, with increased force and OA structural change affecting the medial femorotibial compartment. BMLs present predominantly medially in varus and laterally in valgus lower limbs [7].
BMLs typically occur in regions of the joint with the most severe structural change, for example, fullthickness cartilage lesions, and associations between BML scores and cartilage volume were reported [8].
A previous systematic review [6] suggested that a greater medial or lateral load was related to compartment-speci c BMLs. Associations between BMLs, cartilage damage and varus angulation, and between weight-bearing and non-weight-bearing pain, necessitate careful phenotyping of large numbers of individuals in order to explore speci c associations between subchondral BMLs in weight-bearing components of the knee and weight-bearing pain. Longitudinal cohorts can help indicate whether speci cally localized BMLs might mediate weight-bearing pain, and support the development of BML scores for the strati cation of treatments aiming to improve weight-bearing pain.
The purpose of this study was to assess the association between BML location and weight-bearing knee pain using the Osteoarthritis Initiative (OAI) dataset. We hypothesized that BML size at the medial femorotibial joint compartment in people with varus angulation is associated with severity of weightbearing pain, independent of non-weight-bearing pain and the other OA related MRI features, and that increasing BML size over time is associated with increased weight-bearing pain severity.

Study design and subjects
Subjects were selected from the OAI, a publicly available multi-center, longitudinal, prospective observational study of knee OA. The OAI dataset includes data of MRI and radiographic images, and Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) questionnaires [9]. Figure 1 shows a owchart of participant selection. Of the 4796 participants enrolled in the OAI study, varus knee OA with BML score information of the MRI Osteoarthritis Knee Score (MOAKS) [10], WOMAC questionnaires, femoro-tibial angle (FTA) and body mass index (BMI) were selected. For subjects with bilateral symptomatic knee OA, one knee was randomly selected using random number generator. We analyzed data from 1412 and 582 knees in the cross-sectional and longitudinal OAI datasets, respectively.

Knee pain assessment
Knee pain was assessed at baseline and 24-month follow up using the WOMAC questionnaires. A recent factor analysis of the ve WOMAC pain questions suggested that there are two constructs within the ve questions [11]. Based on the result, we de ned weight-bearing pain as the sum of three WOMAC questions (pain on climbing stairs, on walking, and on standing) and non-weight-bearing pain as the sum of two WOMAC questions (pain in bed and when sitting or lying). The time points of interest for crosssectional and longitudinal study were the baseline visit, and baseline and 24-month follow up visits.

Statistical analysis
Statistical analyses were performed with SPSS version 26 software. Spearman's correlation coe cients were used for correlations between BML score and weight-bearing pain. Pearson's correlation coe cients were used correlations for between changes in BML scores and changes in weight-bearing pain. Negative binomial regression was performed for analysis of associations between weight-bearing or no-weightbearing pain or total WOMAC score and scores of BML, articular cartilage, osteophytes, Hoffa's synovitis, effusion and medial meniscus extrusion. P<0.05 indicated statistical signi cance.

Results
The participant characteristics in this cross-sectional and longitudinal study are shown in Table 1.  Greater BML, cartilage, osteophyte, Hoffa's synovitis or effusion-synovitis scores were each associated with total pain scores adjusted for age, sex and BMI (Table 2). Associations between BML, osteophyte and effusion-synovitis scores and total pain scores remained signi cant after also adjusting for the other OA related MRI features and anterior cruciate ligament tears (Table 2). Greater BML, cartilage, osteophyte, Hoffa's synovitis or effusion-synovitis scores were each associated with weight bearing pain adjusted for non-weight bearing pain, age, sex, BMI (Table 3). Only association between BML and weight-bearing pain remained signi cant after adjusting for the other OA related MRI features and anterior cruciate ligament tears (Table 3). Separate models explored association of BML localization with weight-bearing pain. Greater BML scores at medial femorotibial and lateral patellofemoral compartments were signi cantly associated with greater weight-bearing pain, and these associations retained signi cance after adjustment for BML scores at the other 4 joint compartments (Table 3). Greater BML, cartilage, osteophyte, Hoffa's synovitis or effusion-synovitis scores were each also associated with non-weight-bearing pain adjusted for weight-bearing pain, age, sex, BMI (Table 4). However, only medial meniscus extrusion score remained signi cantly associated with non-weightbearing pain after adjusting for the other OA-related MRI features and anterior cruciate ligament tears (Table 4). Changes from baseline to 24-onth follow up of variables in the longitudinal study are shown in Table 5.
Changes in total and site-speci c BML scores were associated with change in weight-bearing pain. Speci cally, change in BML score at medial femorotibial compartments was signi cantly associated with change in weight-bearing pain after adjusting for sex, baseline weight-bearing pain, age, BMI, BML at the other four joint compartments and change in non-weight-bearing pain ( Table 6). No signi cant association between change in medial meniscus extrusion score and change in non-weight-bearing pain was found after adjusting for sex, baseline non-weight-bearing pain, age, BMI and medial meniscus extrusion score, and change in weight-bearing pain (B = 0.15 (95%CI, -0.15-0.44) p = 0.33).

Discussion
In cross-sectional analyses, we rstly demonstrated that BML size at the medial femorotibial joint compartment in varus knee OA participants was associated with weight-bearing pain severity, and that this association was not explained by non-weight-bearing pain, other OA related MRI features, age, sex or BMI. In longitudinal analyses, increasing BML size at the medial femorotibial joint compartment was associated with increased weight-bearing pain severity. Our ndings support the hypothesis that BMLs increase OA knee pain due to biomechanical factors acting through the affected subchondral bone.
Greater BML scores at medial femorotibial compartments were signi cantly associated with greater total and weight-bearing pain, even after adjusting for non-weight-bearing pain and the other OA related MRI features. Change in BML score at the medial femorotibial compartment was signi cantly associated with change in weight-bearing pain. In contrast to these ndings, BML scores were not signi cantly associated with non-weight-bearing pain in our fully adjusted models. Hence, we conclude that BMLs at weightbearing components are speci cally associated with weight-bearing pain. We previously demonstrated that subchondral bone histopathology characteristic of BMLs, occurring in middle third of medial tibial plateau (an important weight-bearing area), was associated with knee OA pain, not dependent on chondropathy and synovitis [12].
BMLs at the lateral patellofemoral joint compartment were also associated with weight-bearing pain. Patellar BMLs have been associated with any knee pain, and with pain when going up or down stairs [13].
The patellofemoral joint is one of the most commonly affected compartments in knee OA, and varus knee deformity has been associated with worsening of patellofemoral OA, especially in the lateral facet [14]. We extend these ndings to identify that BMLs at the lateral, and not at medial patellofemoral joint compartments were signi cantly associated with the severity of weight-bearing pain in people with varus angulation. This again suggests a biomechanical explanation linking BMLs to OA pain. However, patellofemoral BML score change was not associated with pain change, suggesting that BML changes at the patellofemoral joint has a smaller impact on changes in weight-bearing pain than does BML change at the femorotibial joint.
BMLs represent regions of subchondral bone histologically characterized as displaying increased bone turnover and expression of factors that can increase nerve sensitization [5]. We previously showed that nerve growth factor expression within osteochondral channels and subchondral osteoclast density each is associated with knee OA pain [12], and calcitonin gene-related peptide immunoreactive sensory nerves within osteochondral channels are associated with pain in human and rat knee OA [15]. Activation of sensory nerves in subchondral bone might contribute to weight-bearing pain in knee OA. BMLs might be an imaging biomarker for pathology which sensitises subchondral nerves, and therefore increases weight-bearing pain in knee OA. Further investigation is needed to clarify the cellular and molecular factors which might mediate the observed association between BMLs and weight-bearing pain.
Our results suggested that BMLs might mediate mechanically-induced pain such as during weightbearing. Biomechanical factors may reciprocally contribute to the pathogenesis of BMLs. Previous studies [6] reported that increased mechanical load due to malalignment of the knee joint is a risk factor for incident or enlarging BMLs in the femorotibial joint. Meniscus and cartilage can be shock absorbers that protect subchondral bone from overloading. Associations between meniscal pathology [16,17] or cartilage loss [16,18] and increasing BML size might be mediated by altered biomechanical forces through the subchondral bone. Patellofemoral bracing, which may reduce contact stress across patellofemoral joint, reduced BML volume at patellofemoral joint compartments in patients with painful patellofemoral OA [19]. BMLs therefore might result from increased mechanical load imposed on subchondral bone. However, lateral wedge insoles, which might reduce the load in the medial femorotibial joint compartment during walking [20], did not signi cantly changes in BMLs and pain in medial knee OA [21], and the extent of biomechanical unloading that is required to reverse BML pathology remains uncertain. Treatments with large effects on mechanical load such as high tibial osteotomy and valgus bracing might be needed to reduce BMLs at medial femorotibial joint compartments.
Regarding pharmacotherapy to BMLs, a previous study reported that the bisphosphonate zoledronic acid reduced BML size and knee pain in people with OA [22]. However, recent meta-analysis [23] and randomized controlled trial [24] did not support analgesic effects of bisphosphonates in knee OA. These studies support that biomechanical unloading is more effective than current pharmacotherapy for BML associated pain. Knee OA pain has multiple sources, and targeting BMLs would be expected to have greatest bene t in those cases where BMLs are the predominant driver of weight-bearing pain.
In this study, medial meniscus extrusion score was signi cantly associated with non-weight-bearing pain severity after adjusting for weight-bearing pain and the other OA related MRI. A recent study reported that patients with medial submeniscal ap tear complained of pains during sleep, but not during daytime activities [25]. Pain associated with meniscal tears might be caused by increased mechanical load on meniscus not only during loaded knee but also during unloaded knee exion-extension motion. However, change in meniscus extrusion score was not associated with change in non-weight bearing pain, which suggests that meniscus extrusion may have a relative small impact on non-weight-bearing pain.
Not only BMLs, but also synovitis [26], effusion-synovitis [27], cartilage defects [28] and osteophytes [29] have previously been associated with knee OA pain. We also demonstrated that osteophyte and effusion scores were independently associated with total pain score, even after adjusting for the other OA related MRI features. Signi cant associations of cartilage and Hoffa's synovitis scores did not retain statistical signi cance after adjusting for the other OA related MRI features including BML score, suggesting that some of these observed associations may be explained by other closely associated OA features.
Our study has several potential limitations. We investigated only 2 time points (baseline and 24 months). A previous study using OAI dataset demonstrated that changes in total knee BML volume after 24 months were positively associated with changes in knee pain severity [3]. However, multiple assessments with shorter intervals are needed to evaluate how uctuations of BML size relate to changes in weightbearing pain. BML size was assessed semi-quantitatively and, although our study had a large sample size, quantitative BML measurements might have provided additional information. Direct measurement of mechanical loading through joint compartments was not possible in the current study, and interventional studies would be required to con rm our conclusions that BMLs are biomarkers for pathology which causes biomechanically-induced nociceptive pain.
In conclusion, we demonstrated that BML size at medial femorotibial joint compartments in varus knee OA was associated with weight-bearing pain, that this association was speci c for weight-bearing rather than non-weight bearing pain, and was over and above any effects of other OA-related MRI features, age, sex and BMI. Our ndings suggest that speci c association of weight-bearing pain with BMLs in weightbearing compartments of the knee indicates that BMLs aggravate biomechanical factors leading to OA pain.

Declarations
Funding This study was supported by a grant from Versus Arthritis (Centre initiative grant number 20777).
Competing interests D.A. Walsh: Grants from Versus Arthritis, while the study was being conducted; grants from P zer Ltd, other from P zer Ltd, personal fees from GlaxoSmithKline, outside the submitted work.
F. McWilliams: grants from P zer Ltd.
The remaining authors have no con icts of interest to declare.
Author Contributions Figure 1 Subject selection from OAI database n = number of participants included in the analysis.

Supplementary Files
This is a list of supplementary les associated with this preprint. Click to download. Supplementaltext.docx