- Research article
- Open Access
Proximal tibial trabecular bone mineral density is related to pain in patients with osteoarthritis
Arthritis Research & Therapyvolume 19, Article number: 200 (2017)
Our objective was to examine the relationships between proximal tibial trabecular (epiphyseal and metaphyseal) bone mineral density (BMD) and osteoarthritis (OA)-related pain in patients with severe knee OA.
The knee was scanned preoperatively using quantitative computed tomography (QCT) in 42 patients undergoing knee arthroplasty. OA severity was classified using radiographic Kellgren-Lawrence scoring and pain was measured using the pain subsection of the Western Ontario and McMaster Universities Arthritis Index (WOMAC). We used three-dimensional image processing techniques to assess tibial epiphyseal trabecular BMD between the epiphyseal line and 7.5 mm from the subchondral surface and tibial metaphyseal trabecular BMD 10 mm distal from the epiphyseal line. Regional analysis included the total epiphyseal and metaphyseal region, and the medial and lateral epiphyseal compartments. The association between total WOMAC pain scores and BMD measurements was assessed using hierarchical multiple regression with age, sex, and body mass index (BMI) as covariates. Statistical significance was set at p < 0.05.
Total WOMAC pain was associated with total epiphyseal BMD adjusted for age, sex, and BMI (p = 0.013) and total metaphyseal BMD (p = 0.017). Regionally, total WOMAC pain was associated with medial epiphyseal BMD adjusted for age, sex, and BMI (p = 0.006).
These findings suggest that low proximal tibial trabecular BMD may have a role in OA-related pain pathogenesis.
Knee osteoarthritis (OA) is a debilitating and painful disease characterized by changes in cartilage and subchondral bone. Pain is a complex combination of social, psychological and biological factors , and is often the primary sign that a patient may be afflicted with OA . Unfortunately, the local biological pain pathogenesis within the knee joint is poorly understood  as it could be related to many structural factors (e.g., altered joint alignment , bone marrow lesions (BMLs) , or cysts ). Knee OA is commonly characterized by altered subchondral properties, including altered subchondral bone thickness , bone volume fraction , and volumetric bone mineral density (BMD) . Importantly, altered BMD may disrupt local innervation  and/or the local mechanical behavior of bone , and thus may be a factor in OA-related knee pain.
To date, research investigating association between OA-related knee pain and bone has focused primarily on bone near the subchondral surface (e.g., subchondral cortical and subchondral trabecular bone) [12, 13]. Adjacent trabecular bone (e.g., epiphyseal bone, metaphyseal bone) is also affected by OA , with observations of thinner trabeculae, lower bone volume fraction, and lower density with progressing OA severity [14,15,16]. To date, there are no studies reporting relationships between epiphyseal or metaphyseal trabecular BMD and pain. A recent finite element (FE) study conducted by Amini et al.  has suggested that low epiphyseal trabecular bone density in OA [14,15,16], which is directly linked to the elastic modulus of epiphyseal bone , may explain OA proximal tibiae being less stiff than normal . Importantly, a less stiff proximal tibia would result in higher bone deformation potentially explaining (at least to some degree) OA-related knee pain.
A clear understanding of pain pathogenesis is crucial for rational therapeutic targeting . Further, as pain is the reason patients seek medical care, rational treatment targeting requires specific understanding of which structures contribute to pain . With the aim of furthering our understanding of potential factors that may influence knee pain, the objective of this study was to investigate relationships between proximal tibial epiphyseal and metaphyseal trabecular BMD and OA-related knee pain.
In total 42 participants with OA were recruited prior to total knee replacement (TKR) (17 male, 25 female; mean age 64, SD ± 10.1 years; mean body mass index (BMI) 28.7 ± 3.7; 18 left, 24 right) . Study exclusion criteria included pregnant women, patients having a revision replacement instead of primary knee replacement, and patients with a prior history of bone pathologic change at the knee joint. The Institutional Research Board of the New England Baptist Hospital approved the study. Informed consent was obtained from all study participants.
OA severity was classified using Kellgren-Lawrence (KL) scoring ; participants had severity scores of 2–4. Pain severity was measured at the affected knee joint using the pain subsection of the Western Ontario McMasters Osteoarthritis Index (WOMAC) . Participants were asked to assess the level of pain in the affected knee joint within the past 24 hours while walking on a flat surface, going up or down stairs, nocturnal pain at night in bed, sitting or lying down, and standing upright using a 5-point Likert scale (0–4). Individual element pain scores were then summed for a possible WOMAC pain score of 20. Summed pain scores ranged from 4 to 16. We also used the Self-Administered Comorbidity Questionnaire  to assess participants for any potential confounding comorbidities (e.g., diabetes mellitus or heart disease).
Computed tomography (CT) scan acquisition
We used a single-energy clinical CT scanner (Lightspeed 4-slice, General Electric, Milwaukee, WI, USA) for bone imaging. A solid quantitative CT (QCT) reference phantom of known bone mineral densities (Model 3 T, Mindways Software Inc, Austin, TX, USA) was placed under the participants and included in all CT scans. Participants were oriented supine within the CT gantry and both legs were simultaneously scanned. Scans included the distal femur, patella, proximal tibia, and the 66% tibial shaft site proximal to the distal tibial endplate . Only the proximal tibia and the 66% tibial shaft site were used in the current analysis.
CT scanning parameters included: 120 kVp tube voltage; 150 mAs tube current-time product; axial scanning plane; 0.625-mm isotropic voxel size (0.625 slice thickness, 0.625 mm × 0.625 mm in-plane pixel size); ~ 250 slices; and ~ 60s scan time. A standard bone kernel (BONE) was used for CT image post-processing. The effective radiation dose was ~ 0.073 mSv per scan, estimated using shareware software (CT-DOSE, National Board of Health, Herley, Denmark). For comparison, the average effective radiation dose during a transatlantic flight from Europe to North America is about 0.05 mSv .
CT image analysis
We used a custom algorithm developed specifically for this study (Matlab 2010b; MathWorks, Natick, MA, USA), combined with manual segmentation to determine epiphyseal and metaphyseal trabecular BMD. We considered the epiphyseal region (subarticular region) as the proximal tibial volume between the subchondral surface and the epiphyseal line . A single user (WDB) performed all segmentations and analyses. As this algorithm was developed specifically for this study, we assessed repeatability in a precision study performed on an independent sample of healthy participants and participants with OA  using recommended methods . In summary, 14 participants were scanned three times with repositioning between each scan (42 scans, 28 degrees of freedom (DOF)). The repeatability expressed as precision error, of each BMD measurement was assessed using root mean square coefficients of variation (CV%) and ranged from 0.7% to 3.6%.
To derive BMD, grayscale Hounsfield units (HU) were converted to equivalent volumetric BMD (mg/cm3 K2HPO4) using subject-specific linear regression equations developed from known densities ranging from − 50 to 375 mg/cm3 K2HPO4 within the QCT phantom included in each individual axial image (r 2 > 0.99)  and interpolation to determine equivalent volumetric BMD values. Higher density values were linearly extrapolated (Fig. 1a). Subject-specific half maximum height thresholds  were then determined to define the proximal tibial subchondral and cortical surfaces. Two 3D image volumes were built, one including the entire proximal tibia as previously described [13, 28] and another by segmenting individual serial images using semi-automatic region growing and manual correction at the epiphyseal line (Fig. 1b). Both sets of imaged volumes were segmented using commercial software (Analyze10.0; Mayo Foundation, Rochester, MN, USA) and an interactive touch-screen tablet (Cintiq 21UX; Wacom, Krefeld, Germany). Imaged volumes were reoriented to a neutral position where medial and lateral plateaus were approximately parallel. We then divided the imaged volumes into medial and lateral compartments, measured by using 40% of the maximum medial-lateral axis of each respective side  (Fig. 1c).
To ensure that trabecular BMD measurements did not include cysts (which would lead to arbitrarily low measures of BMD) [13, 30] or peripheral high-density cortical bone, the most proximal 7.5-mm region (relative to the subchondral surface) was removed from the segmentations (Fig. 1d), as was 2.5 mm of peripheral cortical bone (Fig. 1d). The 7.5-mm depth was based upon observed cyst locations from our earlier work [13, 30] and work by Chiba et al. , which limited depth analyses to 5 mm from the subchondral surface. In extreme cases, large cysts extended from the subchondral cortical region (0 − 2.5 mm) through the subchondral trabecular region (2.5–5 mm) and occasionally into depths greater than 5 mm from the subchondral surface. By using a conservative 7.5-mm depth from the subchondral surface, we ensured the exclusion of large cysts from our analysis. Following material removal, we measured epiphyseal trabecular BMD from the 7.5-mm depth to the epiphyseal line (Fig. 1e), which was located approximately 15 mm from the subchondral surface. Metaphyseal trabecular BMD was measured 10 mm distal to the epiphyseal line (Fig. 1e).
We included cortical BMD of the tibial shaft (66% of the tibial length, proximal from the distal tibial plateau)  to assess whether associations with pain were systemic or joint-specific. More specifically, if similar associations between pain and BMD were observed at the proximal tibia and tibial shaft, this would indicate systemic effects with low BMD being a plausible secondary effect of other factors, such as mechanical loading, nutrition or medication . Tibial shaft cortical BMD was segmented using subject-specific half-maximum-height thresholds, and measured using commercial software (Analyze10.0; Mayo Foundation, Rochester, MN, USA).
We first checked all underlying assumptions for multiple linear regression (assumptions of linear relationships, homoscedasticity, independency and normality of residuals) using standardized residual scatter plots, P-P plots, and histograms . We identified any outliers using the modified Thompson tau (τ) test .
We report univariate correlation coefficients (Pearson) between pain, BMD, age, sex, and BMI and illustrate associations between pain and BMD with scatter plots and coefficients of determination (R 2) from linear regression. We used hierarchical multiple linear regression analyses to explain the variance in total WOMAC pain. We selected age, sex, and BMI as covariates for our base model based on observed correlation in univariate analysis (age and WOMAC pain) and literature (age, sex, and BMI) evaluating relationships between BMD and pain [12, 35]. All BMD measurements (total and regional epiphyseal BMD, total metaphyseal BMD, and tibial shaft cortical BMD) were individually added to our base model. We assessed multicollinearity between all independent variables in each model using variance inflation factor (VIF), setting the maximum tolerance value as 10. We report adjusted R 2, change in R 2 from the base model (Δ), standardized beta (β)-coefficients, and p values. Statistical significance was defined as p < 0.05, and analyses were performed using SPSS 21.0 (IBM, Armonk, NY, USA).
Characteristics of all study participants, including age, sex, BMI, KL grades, joint space narrowing (JSN) score, non-weight-bearing alignment scores, and BMD measurements are shown in Table 1. As per the modified Thompson τ test , we identified a single outlier based on the total WOMAC pain score with a τ value outside of the sample’s rejection zone (τ > 5.56), and removed it from the analysis. All underlying assumptions for linear regression were appropriately met. There was no evidence of multicollinearity between independent variables in any of our models. Unadjusted relationships between total WOMAC pain and total or regional epiphyseal or metaphyseal BMD measurements are presented in Fig. 2. Pearson correlation analyses in all participants, and in male and female patients are presented in Additional files 1, 2 and 3: Tables S1–S3.
Regression models predicting variance in pain are presented in Table 2. After adding total epiphyseal BMD to the base model, (of age, sex, and BMI) the coefficient of determination (R 2) for total pain improved (ΔR 2 = 0.12). Our models improved when medial epiphyseal BMD (ΔR 2 = 0.15) and metaphyseal BMD (ΔR 2 = 0.12) were independently added to our base model. There was no association between cortical BMD at the 66% tibial site and pain.
Our regression models suggested that tibial epiphyseal and metaphyseal BMD independently explained variance in total pain in patients with OA prior to TKR, whereby patients with lower BMD tended to have higher levels of pain. Regionally, our models indicated that medial epiphyseal BMD was a significant predictor of total OA-related pain, again whereby lower BMD was associated with higher levels of pain. These findings suggest that there may be potentially overlooked characteristics in proximal tibial BMD, such as trabecular BMD, which may have a role in the pathogenesis of OA-related pain.
The study findings support our previous research (using the same cohort), which investigated links between OA-related nocturnal pain and subchondral cortical and subchondral trabecular bone near the subchondral surface (0–10 mm from the surface). This previous study found a (nonsignificant) trend toward low medial BMD  in patients with severe nocturnal pain, which is in agreement with the study findings of low medial epiphyseal and total epiphyseal and metaphyseal BMD in patients with high levels of pain. We , and others , however, have questioned whether our previously observed trend toward low medial BMD was due to the presence of cysts or diminished bone architecture and/or mineralization. Subsequent follow-up analyses indicated that both cysts and BMD were independently associated with pain . The novelty of this study was that we focused our analyses in epiphyseal and metaphyseal trabecular regions largely void of cysts to determine any potential independent associations between BMD and pain.
Of note, the study findings both support and contrast with the previous study which also identified high lateral focal BMD in the subchondral trabecular region (2.5–10 mm below the surface) in patients with severe nocturnal pain . High lateral focal BMD may be explained by the presence of BMLs, chondro-protection, or altered loading. First, prior research in this cohort identified a positive association between nocturnal pain and BMLs . Given that BMLs have higher local BMD than surrounding bone tissue , a positive association between nocturnal pain and BMD is foreseeable. Future research needs to evaluate whether high focal BMD measurements exactly overlay the BML locations. Second, high lateral focal BMD may be a consequence of chondro-protection developed via low trabecular bone density. To explain, recent finite element (FE) simulations indicated that reduced proximal tibial trabecular bone density results in lower subchondral bone stiffness  and lower cartilage stresses , the latter presumably due to improved congruence between articulations . As many of the study participants had evidence of medial OA, low trabecular BMD may be a physiologic response to lessen medial cartilage stress. At the same time, this chondro-protective process would also naturally transfer more load to the lateral compartment since the two compartments function in parallel. This altered loading should result in loading-induced adaptation; specifically higher lateral BMD near the subchondral surface to meet the mechanical demands of higher load transmission. Third, many of the study participants with evidence of medial OA may be self-altering their knee kinematics and stance to off-load the medial compartment, with the aim of alleviating joint pain. This altered loading could result in loading-induced adaptation with higher lateral BMD and lower medial BMD . Fourth, as higher BMD appears to be focused in subchondral regions (<10 mm from the tibial surface) , joint load may be primarily transferred through the subchondral cortical endplate and subchondral trabecular bone to the peripheral cortex, off-loading epiphyseal and metaphyseal trabecular bone, thus explaining lower BMD in these regions. However, this explanation warrants further research given that we did not find association between pain and alignment . Studies using subject-specific FE modeling are needed to investigate load transmission and subchondral bone stiffness at different stages of pain severity and disease progression.
In this study we report a significant association between age and pain assessed by WOMAC, whereby older participants reported lower pain. Specifically, younger male patients reported higher WOMAC pain scores (Additional file 2: Table S2). We recommend further analysis in larger longitudinal studies to evaluate if this finding is unique to this sample or if this is more widespread within patients with OA. It is also worthwhile noting that we report no associations between age and BMD (Table 2, Additional files 1, 2 and 3: Tables S1–S3). This is in agreement with previous OA research reporting no association between age and BMD  or age and bone volume fraction [42, 43]. Although there is consensus that bone loss is associated with normal aging , this association appears not to pertain to bone tissue within the joint in OA. In support of this, there was no collinearity concern between BMD and age in our models predicting variance in pain. To further explore these associations, we ran the analysis with BMD as the dependent variable, pain as the independent variable, and age, sex, and BMI as covariates. These models suggested pain to be an independent predictor of BMD (Additional file 4: Table S4).
This study has certain limitations. First, pain severity and assessment was based on the entire knee joint, including all joint surfaces (tibiofemoral and patellofemoral) and tissues (e.g., bone, menisci, and synovium), and it is uncertain if pain originated within the proximal tibial bony structure, other tissues, or a combination of tissues. Second, although OA severity was homogeneous across study participants, all were in late stages of OA and it may not be possible to apply our findings to patients in the early stages of OA. Third, our study sample size was small (n = 41). Further analysis with larger samples including healthy participants and participants with various stages of OA severity and pain, are needed to verify these preliminary study findings. Of note, our sample comprised participants with severe OA (primarily with KL scores of 3–4). This limited range constrained our ability to include it in the statistical model. Also, with a basic rule of a minimum of 10 events (or samples) per predictor , we were limited to four predictors (independent variables) in each model: one independent variable (BMD) and three covariates (age, sex, and BMI), and thus other known predictors of pain were not assessed or investigated (e.g., smoking/alcohol history , activity level , mental health status , and specific medications). Of note, we attempted to account for possible differences in physical activity (mechanical loading/unloading) through use of cortical BMD measures at the 66% tibial shaft site. Previous work has identified differences in tibial shaft cortical BMD between highly active individuals (e.g., sprinters, endurance runners, triple-jumpers, high-jumpers, and hurdlers) and less active controls . However, in this study, we did not note any associations between pain and tibial shaft cortical BMD, potentially indicating, at least to some degree, similar levels of activity and mechanical loading amongst study participants. Fourth, our 0.625-mm isotropic voxel size prevented assessment of trabecular microarchitecture and limited us to measurements of volumetric BMD. Accordingly, it is unclear if low BMD is due to trabecular thinning or wide trabecular spacing. For future research, it would be advantageous to investigate links between pain and trabecular microarchitecture with advanced texture analysis and smaller voxel sizes .
In this study we present statistically significant relationships as opposed to clinically significant relationships. As a statistically significant relationship does not measure the clinical effect of a result , it is important to consider the clinical effect that changes in epiphyseal or metaphyseal BMD may have on OA-related knee pain. According to Angst et al. , the minimal clinically important difference for OA-related pain is a change in WOMAC score greater than 6% of its maximum value (which is 20 for WOMAC). In other words, a change in pain will not be perceived unless the WOMAC score changes by 1.2 points. With this in mind, we can identify the BMD change that will correspond with a 1.2-point change in pain. Based on our model, a 44 g/cm3 reduction in epiphyseal or metaphyseal BMD will be marked by a perceived change in pain status. Assuming an average BMD of 100 g/cm3 for epiphyseal and metaphyseal bone, this would equate with ~ 50% change in density. Accordingly, a rational therapeutic approach would be to monitor bone while simultaneously striving to maintain bone and limit bone loss. Density changes in these regions could be monitored using QCT, dual-energy x-ray absorptiometry (DXA) or radiography. With regards to maintaining bone, potentially, this could be achieved through exercise interventions or pharmacological therapies. Our preliminary findings may also be clinically important for TKR preparation and planning. Patients with low preoperative BMD have been shown to be at higher risk of implant failure by loosening or migration , higher risk of revision surgery , and risk of failure following revision procedures . Current tibial implant design components typically include a single central post, which is inserted through the tibial epiphysis and extends into the tibial shaft. Based on our findings, there may be low quantities of bone stock in individuals with higher levels of OA-related pain, potentially placing them at risk of inadequate osseo-integration and implant fixation  and possibly implant loosening . As there is an expected normal decrease in tibial BMD during healing , reduced amounts of tibial epiphyseal bony support structure prior to TKR could compromise implant fixation and success in the early stages, potentially compromising long-term implant success. It may be beneficial to use imaging and complementary image-processing techniques to evaluate preoperative bone density, especially in the commonly overlooked tibial epiphyseal and metaphyseal regions, to compliment customized surgical approaches in patients with higher levels of pain.
In our study, low tibial epiphyseal and metaphyseal BMD, and low medial epiphyseal BMD, was associated with OA-related pain in patients with severe OA prior to TKR. This study suggests that there may be overlooked characteristics within trabecular bone that may be related to the pathogenesis of OA-related pain in patients with severe OA. These preliminary findings from current and previous studies  may be valuable in guiding outcome selection in OA studies addressing subchondral bone and pain, particularly in determining regions of interest of the proximal tibia for potential epidemiological studies.
Bone mineral density
Bone marrow lesion
Coefficients of variation
Dual-energy x-ray absorptiometry
Joint space narrowing
Quantitative computed tomography
Total knee replacement
Variance inflation factor
Western Ontario and McMaster Universities Arthritis Index
Hawker GA, Stewart L, French MR, Cibere J, Jordan JM, March L, et al. Understanding the pain experience in hip and knee osteoarthritis - an OARSI/OMERACT initiative. Osteoarthr Cartil. 2008;16:415–22.
Hunter DJ, Felson DT. Osteoarthritis. Br Med J. 2006;332:639–42.
Dieppe P, Lohmander L. Pathogenesis and management of pain in osteoarthritis. Lancet. 2005;365:965–73.
Lo GH, Harvey WF, McAlindon TE. Associations of varus thrust and alignment with pain in knee osteoarthritis. Arthritis Rheum. 2012;64:2252–9.
Lo GH, McAlindon TE, Niu J, Zhang Y, Beals C, Dabrowski C, et al. Bone marrow lesions and joint effusion are strongly and independently associated with weight-bearing pain in knee osteoarthritis: data from the osteoarthritis initiative. Osteoarthr Cartil. 2009;17:1562–9.
Kornaat PR, Bloem JL, Ceulemans RYT, Riyazi N, Rosendaal FR, Nelissen RG, et al. Osteoarthritis of the knee: association between clinical findings and MR imaging findings. Radiology. 2006;239:811–7.
Buckland-Wright C. Subchondral bone changes in hand and knee osteoarthritis detected by radiography. Osteoarthritis Cartil. 2004;12:10.
Bobinac D, Spanjol J, Zoricic S, Maric I. Changes in articular cartilage and subchondral bone histomorphometry in osteoarthritic knee joints in humans. Bone. 2003;32:284–90.
Zysset PK, Sonny M, Hayes WC. Morphology-mechanical property relations in trabecular bone of the osteoarthritic proximal tibia. J Arthroplasty. 1994;9:203–16.
Mach DB, Rogers SD, Sabino MC, Luger NM, Schwei MJ, Pomonis JD, et al. Origins of skeletal pain: sensory and sympathetic innervation of the mouse femur. Neuroscience. 2002;113:155–66.
Radin EL, Paul IL, Rose RM. Role of mechanical factors in pathogenesis of primary osteoarthritis. Lancet. 1972;1:519–22.
Akamatsu Y, Mitsugi N, Taki N, Ashi HK, Saito T. Medial versus lateral condoyle bone mineral density ratios in a cross-sectional study: a potential marker for medial knee osteoarthritis severity. Arthritis Care Res. 2012;64:1036–45.
Burnett WD, Kontulainen S, McLennan C, Hazel D, Talmo C, Hunter D, et al. Knee osteoarthritis patients with severe nocturnal pain have altered subchondral tibial bone mineral density. Osteoarthr Cartil. 2015;23:1483–90.
Ding M, Danielsen CC, Hvid I. Bone density does not reflect mechanical properties in early-stage arthrosis. Acta Orthop Scand. 2001;72:181–5.
Day JS, Ding M, van der Linden JC, Hvid I, Sumner DR, Weinans H. A decreased subchondral trabecular bone tissue elastic modulus is associated with pre-arthritic cartilage damage. J Orthop Res. 2001;19:914–8.
Bennell KL, Creaby MW, Wrigley TV, Hunter DJ. Tibial subchondral trabecular volumetric bone density in medial knee joint osteoarthritis using peripheral quantitative computed tomography technology. Arthritis Rheum. 2008;58:2776–85.
Amini M, Nazemi SM, Lanovaz J, Kontulainen S, Masri BA, Wilson DR, et al. Individual and combined effects of OA-related subchondral bone alterations on proximal tibial surface stiffness: a parametric finite element modeling study. Med Eng Phys. 2015;37:783–91.
Carter DR, Hayes WC. The compressive behavior of bone as a two-phase porous structure. J Bone Joint Surg Am. 1977;59:954–62.
Neogi T. Clinical significance of bone changes in osteoarthritis. Ther Adv Musculoskelet Dis. 2012;4:259–67.
Kellgren JH, Lawrence JS. Radiological assessment of osteo-arthrosis. Ann Rheum Dis. 1957;16:494–502.
Bellamy N, Buchanan WW, Goldsmith CH, Campbell J, Stitt LW. Validation study of WOMAC: a health status instrument for measuring clinically important patient relevant outcomes to antirheumatic drug therapy in patients with osteoarthritis of the hip or knee. J Rheumatol. 1988;15:1833–40.
Sangha O, Stucki G, Liang MH, Fossel AH, Katz JN. The Self-Admininstered Comorbidity Questionnaire: a new method to assess comorbidity for clinical and health services research. Arthritis Rheum. 2003;49:156–63.
Kontulainen SA, Johnston JD, Liu D, Leung C, Oxland TR, McKay HA. Strength indices from pQCT imaging predict up to 85% of variance in bone failure properties at tibial epiphysis and diaphysis. J Musculoskelet Nueronal Interact. 2008;8:401–9.
UNSCEAR 2000 Report to the General Assembly - Annex B: Exposures from natural radiation sources. United Nations Scientific Committee on the Effects of Atomic Radiation.
Ljunggren AE. Variations in the relationship between the diaphysis and the epiphyses of the tibia. Acta Morphol Neerl Scand. 1976;14:101–37.
Johnston JD, McLennan CE, Hunter DJ, Wilson DR. In vivo precision of a depth-specific topographic mapping technique in the CT analysis of osteoarthritic and normal proximal tibial subchondral bone density. Skelet Radiol. 2010;40:1057–64.
Glüer CC, Blake G, Lu Y, Blunt BA, Jergas M, Genant HK. Accurate assessment of precision errors: how to measure the reproducibility of bone densitometry techniques. Osteoporos Int. 1995;5:262–70.
Johnston JD, Masri BA, Wilson DR. Computed tomography topographic mapping of subchondral density (CT-TOMASD) in osteoarthritic and normal knees: methodological development and preliminary findings. Osteoarthr Cartil. 2009;17:1319–26.
Spoor CF, Zonneveld FW, Macho GA. Linear measurements of cortical bone and dental enamel by computed tomography: applications and problems. Am J Phys Anthropol. 1993;91:469–84.
Chen Y, Huang Y-C, Lu WW. Is subchondral bone mineral density associated with nocturnal pain in knee osteoarthritis patients? Osteoarthritis Cartil. 2015;23:2297–8.
Chiba K, Burghardt AJ, Osaki M, Majumdar S. Three-dimensional analysis of subchondral cysts in hip osteoarthritis: an ex vivo HR-pQCT study. Bone. 2014;66:140–5.
Frank AW, Labas MC, Johnston JD, Kontulainen SA. Site-specific variance in radius and tibia bone strength as determined by muscle size and body mass. Physiother Can. 2012;64:292–301.
Schneider A, Hommel G, Blettner M. Linear regression analysis. Dtsch Arztebl Int. 2010;107:776–82.
Thompson R. A note on restricted maximum likelihood estimation with an alternative outlier model. J R Stat Soc B Methodol. 1985;47:53–5.
Dore D, Ding C, Jones G. A pilot study of the reproducibility and validity of measuring knee subchondral bone density in the tibia. Osteoarthritis Cartil. 2008;16:1539–44.
Burnett W, Kontulainen S, McLennan C, Hazel D, Talmo C, Hunter D, et al. Response to Letter to the Editor: ‘Is subchondral bone mineral density associated with nocturnal pain in knee osteoarthritis patients?’. Osteoarthritis Cartil. 2015;23:2299–301.
Seah S, Wheaton D, Li L, Dyke JP, Talmo C, Harvey WF, et al. The relationship of tibial bone perfusion in knee osteoarthritis. Osteoarthr Cartil. 2012;20:1527–33.
Lowitz T, Museyko O, Bousson V, Laouisset L, Kalendar WA, Laredo JD, et al. Bone marrow lesions identified by MRI in knee osteoarthritis are associated with locally increased bone mineral density measured by QCT. Osteoarthr Cartil. 2013;21:957–64.
Venäläinen MS, Mononen ME, Jurvelin JS, Töyräs J, Virén T, Korhonen RK. Importance of material properties and porosity of bone on mechanical response of articular cartilage in human knee joint: a two-dimensional finite element study. J Biomech Eng. 2014;136:121005.
Hertz VHH. Uber die Berührung fester elastischer Körper. J Reine Angew Math. 1882;92:156–71. German.
Wada M, Maezawa Y, Baba H, Shimada S, Sasaki S, Nose Y. Relationships among bone mineral densities, static alignment and dynamic load in patients with medial compartment knee osteoarthritis. Rheumatology (Oxford). 2001;40:499–505.
Li G, Zheng Q, Landao-Bassonga E, Cheng TS, Pavlos NJ, Ma Y, et al. Influence of age and gender on microarchitecture and bone remodeling in subchondral bone of the osteoarthritic femoral head. Bone. 2015;77:91–7.
Perilli E, Baleani M, Ohman C, Baruffaldi F, Viceconti M. Structural parameters and mechanical strength of cancellous bone in the femoral head in osteoarthritis do not depend on age. Bone. 2007;41:760–8.
Mazess RB. On aging bone loss. Clin Orthop Relat Res. 1982;165:239–52.
Peduzzi P, Concato J, Kemper E, Holford TR, Feinstein AR. A simulation study of the number of events per variable in logistic regression analysis. J Clin Epidemiol. 1996;49:1373–9.
Miranda H, Viikari-Juntra E, Martikainen R, Riihimaki H. A prospective study on knee pain and its risk factors. Osteoarthritis Cartil. 2002;10:623–30.
Wesseling J, Welsing PMJ, Bierma-Zeinstra SMA, Dekker J, Gorter KJ, Kloppenburg M, et al. Impact of self-reported comorbidity on physical and mental health status in early sympotmatic osteoarthritis: the CHECK (Cohort Hip and Cohort Knee) study. Rheumatology. 2013;52:180–8.
Rantalainen T, Nikander R, Daly RM, Heinonen A, Sievanen H. Exercise loading and cortical bone distribution at the tibial shaft. Bone. 2011;48:786–91.
Wasserstein R, Lazar N. The ASA’s statements on p-values: context, process, and purpose. Am Stat. 2016;70:129–33.
Angst F, Aeschlimann A, Stucki G. Smallest detectable and minimal clinically important differences of rehabilitation intervention with their implications for required sample sizes using WOMAC and SF-36 Quality of Life measurement instruments in patients with osteoarthritis of the lower extremities. Arthritis Care Res (Hoboken). 2001;45:384–91.
Petersen MM, Nielsen PT, Lebech A, Toksvig-Larsen S, Lund B. Preoperative bone mineral density of the proximal tibia and migration of the tibial compartment after uncemented total knee arthroplasty. J Arthroplasty. 1999;14:77–81.
Levitz CL, Lotke PA, Karp JS. Long-term changes in bone mineral density following total knee replacement. Clin Orthop Relat Res. 1995;321:68–72.
Mavrogenis AF, Dimitriou R, Parvizi J, Babis GC. Biology of implant osseointegration. J Musculoskelet Nueronal Interact. 2009;9:61–71.
Sharkey PF, Hozack WJ, Rothman RH, Shastri S, Jacoby SM. Why are total knee arthroplasties failing today? Clin Orthop Relat Res. 2002;404:7–13.
Ritter MA, Davis KE, Small SR, Merchun JG, Farris A. Trabecular bone density of the proximal tibia as it relates to failure of a total knee replacement. Bone Joint J. 2014;96-B:1503–9.
This project was funded through support from the Canadian Arthritis Network (CAN) and New England Baptist Hospital Research Funding Awards.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
This study was approved by the Institutional Research Board of the New England Baptist Hospital. Informed consent was obtained from all study participants.
Consent for publication
The authors declare that they have no competing interests. The corresponding author had full access to all the data in the study and had final responsibility for the decision to submit for publication.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Coefficients (r) with 95% confidence intervals for correlation between all model variables for all included participants (n = 41). Significant associations are in bold. (DOCX 12 kb)
Coefficients (r) with 95% confidence intervals for correlation between all model variables for included male participants (n = 17). Significant associations are in bold. (DOCX 12 kb)
Coefficients (r) with 95% confidence intervals for correlation between all model variables for included female participants (n = 24). Significant associations are in bold. (DOCX 12 kb)
Adjusted coefficients of determination (R 2) and standardized beta coefficients (β) of the base model I (age, sex, and BMI) and base model II (age, sex, BMI and WOMAC pain) to predict variance in bone mineral density (BMD) at the total and regional proximal tibia. Significant R 2 and β values are in bold; p values are in parentheses. (DOCX 13 kb)