The Association of Blood Biomarkers and Body Mass Index in Knee Osteoarthritis: A Cross-Sectional Study.

OBJECTIVE: Despite massive efforts, there are no diagnostic blood biomarkers for knee osteoarthritis (KOA). This study investigated several candidate diagnostic biomarkers and the metabolic phenotype in end-stage KOA in the context of obesity.
DESIGN: In this cross-sectional study, adult patients undergoing knee arthroplasty were enrolled and KOA severity was assessed using the Lequesne index. Blood biomarkers with an important role in obesity, the metabolic syndrome, or KOA (oxidized form of low-density lipoprotein [oxLDL], advanced glycation end product [AGE], soluble AGE receptor [sRAGE], fatty acid binding protein 4 [FABP4], phospholipase A2 group IIA [PLA2G2A], fibroblast growth factor 23 [FGF-23], ghrelin, leptin, and resistin) were measured using enzyme-linked immunosorbent assay (ELISA; n = 70) or Luminex technique (subgroup of n = 35). H1-NMR spectroscopy was used for the quantification of metabolite levels (subgroup of n = 31). The hip-knee-ankle angle was assessed. Multivariable and multivariate regression analysis was used to examine the relationship of biomarkers with body mass index (BMI) and KOA severity in complete case and multiple imputation analysis.
RESULTS: While most of the investigated biomarkers were not associated with KOA severity, FABP4 and leptin were found to correlate with BMI and gender. Resistin was associated with Lequesne index in complete case analysis. Using a targeted metabolomics approach, BMI-dependent changes in the metabolome were hardly visible.
CONCLUSIONS: Our findings confirm studies on FABP4, leptin, and resistin with regard to obesity and the metabolic syndrome. There was no association of the investigated biomarkers with KOA severity, most likely due to the patient selection (end-stage KOA patients). Based on this absence of BMI-dependent changes in the metabolome, we might assume that BMI is not correlated with KOA severity in this specific patient group.

Introduction

Background

Knee osteoarthritis (KOA) is a major cause of disability and pain worldwide. Despite massive efforts over the past decades, no blood serum marker has emerged as diagnostic marker so far.

Based on recent research results, inflammation, obesity, and the metabolic syndrome have been described as risk factors for KOA manifestation and progression of OA.

Obesity, defined as a body mass index (BMI) of greater than 30, is a silent epidemic associated with several chronic conditions and affecting more than 2.1 billion individuals worldwide.[3,4] It is also a major contributing factor to the pathogenesis and progression of KOA due to adipocyte-derived systemic inflammation. Obesity-related osteoarthritis has been described as a complex, multifactorial condition that can cause significant impact on quality of life of affected patients. There is a close correlation of BMI and clinical and functional consequences of KOA. This observation is relevant when considering KOA treatment strategies based on obesity. Investigating candidate biomarkers to assess a possible correlation to BMI and/or KOA severity seems to be an important step in finding diagnostic blood biomarkers.

In addition, obesity represents a key finding for the metabolic syndrome (metS). The metS is defined as a combination of several disorders associated with cardiovascular risk. It causes local and systemic changes that induce a pro-inflammatory state leading to oxidative stress and eventually cartilage degradation.

Although adipocytokines and the metabolome are considered to be important pathogenic factors in obesity-related OA, there has been little research of these factors as diagnostic biomarkers.

Leptin, resistin, and grehlin play a major role in obesity and metS. These markers contribute to the imbalance of joint homeostasis and are involved in the pathogenesis of osteoarthritis.[8,9]

The oxidized form of low-density lipoprotein (oxLDL) is not only a marker of metabolic disease, but plays a major role in cartilage degeneration.

Fatty acid binding protein 4 (FABP4) is a adipokine that is closely associated with obesity and metabolic disease. In knock-out models and after pharmaceutical inhibition of FABP4, OA induced by high-fat diet in mice was alleviated.

The receptor/ligand system of advanced glycation end products (AGEs) and their cell receptor (RAGE) are similarly involved in the pathophysiology of metabolic diseases. The AGE-RAGE interaction results in increased generation of oxygen radicals and pro-inflammatory cytokines. Circulating soluble AGE receptors (sRAGE) interact with AGEs to counterbalance the negative effects of AGEs-RAGE interaction.

The finding that PLA2G2A is highly abundant in biological fluids of patients suffering from inflammatory diseases, such as arthritis, sepsis, and myocardial infarctions, emphasizes the important role of PLA2G2A in inflammation. Elevated levels of PLA2G2A have been described as biomarker for cardiovascular disease and the metS.

Fibroblast growth factor 23 (FGF-23) is a regulator of cartilage differentiation and might also play a role in the pathogenesis of osteoarthritis.

Previous metabolomics studies have identified changes in branched-amino acid (BCAA) and arginine levels in serum and synovial fluid of KOA patients. Another study, using an untargeted metabolomics approach, found several differences in obese KOA patients compared with nonobese KOA patients.

To the best of our knowledge, these peptides and proteins with important roles in cartilage matrix turnover (PLA2G2A), glucose metabolism (AGE, soluble RAGE), and chondrocyte differentiation (FGF-23) have not been studied extensively in the context of obesity and end-stage KOA before.

In this article, we examined abovementioned biomarkers as markers of burden of disease in end-stage KOA, measured using the Lequesne index. This study was conducted to investigate whether any of these biomarkers are associated with obesity, the metS, and KOA and to test their ability in characterizing KOA severity.

Methods

This study was conducted as a cross-sectional study at the Department of Orthopedics and Trauma of the Medical University of Graz, Austria, from January 2019 to April 2020. The procedures described were in accordance with the ethical standards of the ethics committee of the Medical University of Vienna (IRB#2029/2016) and the Medical University of Graz (IRB#31-133ex18/19) and with the Helsinki Declaration.

Study Population

Adult patients undergoing knee arthroplasty were enrolled in the “Better Life in Osteoarthritis Registry” (BLOAR) after informed consent. The BLOAR is an Austrian multicenter registry for OA patients with the goal to further research in the prevention and treatment of osteoarthritis. A total of 70 patients were grouped based on BMI (group 1 = “underweight”: BMI < 20 [n = 19], group 2 = “normal weight”: 20 to 30 [n = 32], group 3 = “obese”: BMI > 30 [n = 19]) and chosen randomly (enrollment based on date of surgery and available resources on the day of enrollment) from the registry for analysis in this study.

Blood samples were taken preoperatively and centrifuged within 1 hour. Serum was stored at −20°C until further processing. To achieve a more distinct separation, and due to limited resources, only patients with the lowest and highest BMI in the underweight and obese group, respectively, were further analyzed with metabolic phenotyping and Luminex assays.

Patient Characteristics

Height, weight, and abdominal circumference and blood pressure measurements were taken. For the purpose of this study, obesity was defined as a BMI of greater than 30, normal weight as a BMI of 20 to 30, and underweight as a BMI less than 20. Previous nonoperative and operative treatment (including physiotherapy, physical therapy, weight-reduction program, intra-articular injections, topical and systemic analgesics) was recorded. In women, hormonal status was recorded (previous gynecologic surgeries, oral contraceptive pills, hormone-replacement treatment). Comorbidities and medication were recorded based on history taking and medication was grouped based on mechanism of action. Based on the National Cholesterol Education Program (NCEP) definition of the metabolic syndrome, patients were diagnosed with metabolic syndrome when 3 or more of the following conditions were met: (1) on antihypertensive medication, known arterial hypertension, or systolic office measurement greater than 130 mmHg; (2) known dyslipidemia or on lipid-lowering drugs; (3) known diabetes mellitus or on antidiabetic drugs; and (4) abdominal obesity (waist circumference > 102 cm for men or > 88 cm for women). The Lequesne index, a patient-reported outcome index, is a questionnaire with 3 sections on pain and discomfort, walking distance, and activities of daily living. A higher score is associated with a greater disability. This index was used to assess KOA severity. The ASA physical status classification system was used to assess the overall patient health status. The ASA physical status classification is a system that can be used to assess the health status and evaluate perioperative risks.

X-Ray Measurements

Preoperative full-length limb standing x-ray images were taken and analyzed for the hip-knee-ankle angle. Measurements were performed using imageJ (version 1.52j) by P.S. and by M.F. in a subset of 16 patients. Interrater agreement was assessed using intraclass correlation coefficient 3,A (ICC3,A). Interrater agreement was excellent (ICC3,A = 0.99, P = 0).

Enzyme-Linked Immunosorbent Assay (ELISA)

The following ready-to-use sandwich ELISA kits were used according to the manufacturer’s instruction: oxidized low-density lipoprotein cholesterol (oxLDL; Mercodia, Uppsala, Sweden), advanced glycation end products (AGE; Cell Biolabs, San Diego, CA), soluble receptor for advanced glycation end products (sRAGE; BioVendor, Brno, Czech Republic), fatty acid binding protein 4 (FABP4; BioVendor), and membrane-bound phospholipase A2 group IIA (PLA2G2A; RayBiotech Life, Peachtree Corners, GA). The AGE-sRAGE ratio was calculated based on these measurements. All measurements were performed in duplicates at 450 nm with a microplate reader (Infinite F50 from Tecan, Austria).

XMAP Human Bone Metabolism Magnetic Bead Panel

Using the Luminex xMAP platform in a magnetic bead format, we simultaneously analyzed the following targets from human serum samples in a subgroup of 35 patients who were selected based on BMI (see above; BMI group 1: n = 16, BMI group 3: n = 19): fibroblast growth factor 23 (FGF-23), ghrelin, leptin, and resistin. For detection, we used the commercially available ProcartaPlex (Thermo Fisher, Waltham, MA) on a Bio-Plex 200 system (Bio-Rad Laboratories, Hercules, CA) according to the manufacturer’s instructions. There was no cross-reactivity between the antibodies for an analyte and any of the other analytes in this panel. Measurement of median fluorescence intensity (MFI) was performed using the Bio-Plex Manager software, version 4.1 (Bio-Rad Laboratories) and used for further analysis.

Metabolic Phenotyping

In another, partly overlapping subgroup of 31 patients (BMI group 1: n = 15, BMI group 3: n = 16), changes in metabolic phenotypes were explored using nuclear magnetic resonance (NMR) spectroscopy (Bruker BioSpin GmbH, Rheinstetten, Germany). Human serum samples from KOA patients were lyophilized and 500 µl of NMR buffer (5.56 g Na2HPO4, 0.4 g trimethylsilylpropanoic acid [TSP], 0.2 g NaN3, in 400 ml of D2O; pH 7.4) were added. All NMR experiments were performed at a temperature of 310 K on a Bruker AVANCE Neo Ultrashield 600 MHz spectrometer equipped with a triple resonance probe head and processed as described previously. Bruker TopSpin version 3.1 (Bruker BioSpin GmbH) was used for NMR data acquisition and automatic processing (exponential line broadening of 0.3 Hz, phasing, and referencing to TSP at 0.0 ppm). Regions around water, TSP, and remaining MeOH signals were excluded. The spectra for all samples were further analyzed, deconvoluted, and quantified automatically, using the R package ASICS. Probabilistic quotient normalization was performed to correct for sample metabolite dilution. Concentrations shown are normalized concentrations.

Statistical Analysis

All analyses were performed using R version 4.0.3 (2020-10-10) on Manjaro Linux 5.6.15-1. The following R packages were used: mice, lmerTest, tidyverse, visreg, irr, ASICS, ClustOfVars, car.[23,25 -31] For Luminex data analysis, a collection of R scripts provided by Breen et al. was used.

In the ELISA data set, hierarchical cluster analysis, based on principal component analysis (PCA) for a mixture of continuous and categorical variables (PCAMIX), was performed on variables to gain a better understanding of the relationship among variables collected.30 For the purpose of this study, we cut the cluster dendogram based on aggregation levels and the Rand index into a total of 11 clusters.

Missing data and sensitivity analysis

There were missing data in the ELISA biomarkers and the Lequesne index (all less than 16%, see Suppl. Tables S1 and S2). It was assumed that data were missing completely at random (MCAR). Details on assumed mechanisms of missingness are outlined in the Supplement. Based on the MCAR assumption, missing values were imputed using multiple imputation chained equation (MICE) with the R package mice. All analyses were performed on complete cases and the multiply imputed data sets. To explore violations of the MCAR assumption, sensitivity analysis was performed. The δ adjustment technique was used for this purpose:[28,33] A δ of ±15% and ±30% of the complete case mean was added to the imputed values to understand how deviations from the MCAR assumptions influence analysis.

Inferential analysis

Normality was assessed using Shapiro-Wilk tests, Q-Q, and density plots. Means were compared using a Student t test or Wilcoxon rank sum test where appropriate. Correlations were assessed using Spearman rank correlation.

To answer our research questions, the following regression analyses were performed: (1) Multivariable linear regression was performed on each ELISA biomarker as dependent variable. All models were adjusted to the known risk factors such as age, gender, BMI, and hip-knee-ankle angle. To understand the relationship with the Lequesne index and the metabolic syndrome, these variables were also added to all models. To detect multicollinearity in fixed-effects models, the variance inflation factor (VIF) was calculated. All VIF values were smaller than 2. (2) Logarithmic transformation was applied to Luminex fluorescence values and these were used in linear mixed models as outcome variable.[21,32] Patient identity was used as cluster variable in a model with random intercepts. Models were adjusted for age, gender, BMI, and HKAA. To understand the association with the Lequesne index and the metabolic syndrome, these variables were also added to the model. (3) For metabolic phenotyping, data were analyzed using PCA. The number of principal components was determined using the cumulative variance (exceeding 50% of cumulative variance). Partial least squares regression and discriminant Analysis (PLSR, PLS-DA) using the R package ASICS was used to test the relationship of metabolite concentrations and BMI group and test scores. PLS and PLS-DA models were assessed using R2 to estimate goodness of fit and Q2 to estimate predictive performance. The number of components added to the model was determined based on R2 and Q2 scores using 7-fold cross-validation. The predictive performance of the full model was assessed based on the cumulative Q2 score using permutation testing (1000 permutations). All models were assessed using diagnostic plots.

All P values were adjusted using the Benjamini-Hochberg correction as implemented in R and rounded to 3 digits, with a value of <.05 indicating significance.

Results

A total of 70 patients were enrolled. Patient characteristics are shown in . The mean age was 70.09 ± 8.92 and mean BMI was 30.42 ± 6.27. A total of 50 women (71.43%) were included. Supplementary Tables S3-A and S3-B show the subgroup patient characteristics.

Table 1.

Study Population by BMI Group.

	Underweight(n = 19)	Normal(n = 32)	Obese(n = 19)	Total(N = 70)
Age	75.16 (6.38)	70.28 (9.41)	64.68 (7.33)	70.09 (8.92)
BMI	22.88 (1.63)	30.18 (2.62)	38.35 (3.07)	30.42 (6.27)
Female	15 (78.9%)	20 (62.5%)	15 (78.9%)	50 (71.4%)
HKAA	7.40 (3.94)	7.63 (4.36)	8.52 (4.02)	7.81 (4.12)
ASA
1	1 (5.3%)	3 (9.4%)	0 (0.0%)	4 (5.7%)
2	9 (47.4%)	17 (53.1%)	8 (42.1%)	34 (48.6%)
3	9 (47.4%)	10 (31.2%)	11 (57.9%)	30 (42.9%)
4	0 (0.0%)	2 (6.2%)	0 (0.0%)	2 (2.9%)
Lequesne index	14.30 (3.73)	12.79 (3.82)	14.83 (2.92)	13.80 (3.61)
TUGT	14.07 (3.84)	12.89 (6.38)	14.71 (5.92)	13.71 (5.67)
Diabetes mellitus	4 (21.1%)	5 (15.6%)	3 (15.8%)	12 (17.1%)
Hyperlipidemia	6 (37.5%)	10 (38.5%)	4 (21.1%)	20 (32.8%)
Abdominal circumference	89.08 (7.30)	106.96 (11.22)	122.76 (9.90)	107.63 (15.95)
CAD	4 (21.1%)	5 (15.6%)	3 (15.8%)	12 (17.1%)
Hysterectomy	0 (0.0%)	5 (17.2%)	7 (36.8%)	12 (19.0%)

BMI = body mass index; HKAA = hip-knee-ankle angle; ASA = ASA physical status classification system; CAD = coronary artery disease; TUGT = timed up and go test.

Hierarchical Cluster Analysis

Data on 66 variables were collected. Hierarchical cluster analysis groups highly correlated variables in clusters to give a better understanding of the meaning of each variable. shows the cluster dendogram. The clusters and variables are listed in Supplementary Table S4.

Figure 1.

Dendogram of hierarchical cluster analysis of variables. A synthetic variable is computed on the basis of principal component analysis for a mixture of continuous and categorical variables (PCAMIX). Variables are grouped based on the correlation with this variable. The cluster dendogram was cut (red line) into 11 clusters based on the agglomeration level and the Rand index. The dendogram height indicates the agglomeration level and thus lower height indicates more similarity. OCP = oral contraceptive pill; GI disorder = gastrointestinal disorder; EMF = dielectric heating; TENS = transcutaneous electrical nerve stimulation; CAD = coronary artery disease; ASA = ASA physical status classification system; HKAA = hip-knee-ankle angle; MS = metabolic syndrome; numbers indicate cluster number.

Hierarchical cluster analysis showed that FABP4 level was closely related to gender and height, whereas sRAGE and oxLDL were closely connected among each other and with hysterectomy. Furthermore, AGE and PLA2G2A revealed a close correlation with each other as well as with osteoporosis. The Lequesne index was in the same cluster as the HKAA. Conservative treatment options (physiotherapy and physical and occupational therapy) formed an additional cluster. These findings were as expected and showed the close connection of variables.

Association of ELISA/Luminex Measurements With Clinical Parameters

FABP4 level was associated with BMI and gender, but not Lequesne index, in KOA.

In univariable analysis, we found a significant positive correlation of BMI and FABP4 levels (ρ = 0.4851, P = 0). When adjusting for age, gender, HKAA, and metS in complete cases, FABP4 was still associated with BMI and gender (see Suppl. Table S5). The same result could be observed in the imputed data sets, where BMI and gender were associated with FABP4 after adjustment (see , ). Thus, higher FABP4 levels were found in women and obese patients. No association with the Lequesne index was observed.

Table 2.

Results of Linear Regression With FABP4 as the Dependent Variable in Multiple Imputation Models.

Variable	Crude Est.	Adj. Est.	2.5%	97.5%	Adj. P Value
Age	−0.35	0.169	–0.355	0.693	.521
BMI	1.64	1.610	0.851	2.370	.001***
Gender (male)	–19.34	–17.201	–26.052	–8.349	.001***
HKAA	0.26	–0.372	–1.411	0.666	.521
Lequesne index	1.25	0.858	–0.416	2.131	.425
No metabolic syndrome	–8.49	–3.871	–12.659	4.916	.521

The model was adjusted for the age, gender, BMI, HKAA, and metabolic syndrome. P values are adjusted. There is a significant positive association of FABP4 levels and BMI and female gender. Crude est. = crude estimate; Adj. est. = adjusted estimate; 2.5% and 97.5% = 2.5% and 97.5% confidence intervals; no metabolic syndrome = patient did not meet metS criteria as outlined in the text; FABP4 = fatty acid binding protein 4; BMI = body mass index; HKAA = hip-knee-ankle angle.

Figure 2.

Top: Association of gender and FABP4 as predicted by linear regression. FABP4 is higher in women. Bottom: Correlation of BMI and FABP4 as predicted by regression. FABP4 increases with BMI. FABP4 = fatty acid binding protein 4; BMI = body mass index.

In the sensitivity analysis, results were stable (see Suppl. Tables S6-A and S6-B). Thus, even if the imputed values deviated by up to 30% from the mean of the observed values, we would still observe a similar effect.

All other biomarkers explored by using the ELISA technique (oxLDL, AGE, sRAGE, PLA2G2A) did not show a significant association with BMI, Lequesne index, or metS (Suppl. Tables S7-A to S7-E).

Leptin was associated with gender and BMI group

In complete cases, leptin was associated with BMI group and gender (Suppl. Table S8). The same results were observed after multiple imputation ( ). This association was further examined using sensitivity analysis and remained stable (Suppl. Tables S9-A and S9-B). Leptin levels were higher in obesity compared with underweight and in women compared with men. shows the relationship of BMI and gender as predicted by the linear mixed model.

Table 3.

Results of Linear Mixed Models on Imputed Data Showing Interactions Only.

Coefficient	Unadj. Est.	Adj. Est.	2.5%	97.5%	Adj. P Value
Ghrelin: Age	0.002	–0.004	–0.058	0.049	.972
Ghrelin: BMI group 3	–0.135	–0.222	–1.173	0.729	.972
Ghrelin: Gender (male)	0.052	0.175	–0.723	1.073	.972
Ghrelin: HKAA	–0.013	–0.008	–0.114	0.098	.972
Ghrelin: Lequesne Index	–0.012	–0.019	–0.141	0.102	.972
Ghrelin: No metabolic syndrome	–0.125	–0.251	–1.083	0.581	.972
Leptin: Age	–0.031	0.028	–0.026	0.081	.787
Leptin: BMI group 3	1.014	1.642	0.688	2.596	.011*
Leptin: Gender (male)	–1.039	–1.418	–2.323	–0.513	.021*
Leptin: HKAA	0.033	0.040	–0.069	0.149	.972
Leptin: Lequesne Index	0.005	–0.031	–0.182	0.121	.972
Leptin: No metabolic syndrome	–0.565	–0.535	–1.383	0.313	.669
Resistin: Age	–0.003	–0.028	–0.082	0.026	.787
Resistin: BMI group 3	–0.321	–0.785	–1.752	0.183	.391
Resistin: Gender (male)	–0.279	–0.234	–1.137	0.668	.972
Resistin: HKAA	0.030	0.127	0.020	0.234	.113
Resistin: Lequesne Index	–0.097	–0.166	–0.301	–0.030	.113
Resistin: No metabolic syndrome	0.169	–0.028	–0.884	0.827	.972

The model was adjusted for age, gender, BMI group, HKAA, and metabolic syndrome. There is a significant positive association of leptin levels and obesity (BMI group 3) and female gender. Unadj. est. = unadjusted estimate; Adj. est. = adjusted estimate; 2.5% and 97.5% = 2.5% and 97.5% confidence intervals; ‘:’ = interaction term; BMI = body mass index; HKAA = hip-knee-ankle angle.

P < .05.

Figure 3.

The relationship of BMI group (top) and gender (bottom) with fluorescence intensity per cytokine is shown as predicted by linear mixed modeling. Leptin levels are higher in women in obese patients. BMI = body mass index; BMI group 1 = underweight; BMI group 3 = obesity; 0 = female gender.

Resistin was associated with Lequesne index in complete case analysis

In complete cases, resistin showed an association with Lequesne index (Suppl. Table S8). After applying multiple imputation, this association was not observed, however ( ). This was confirmed in the sensitivity analysis (Suppl. Tables S9-A and S9-B) demonstrating that only under the most were able to explain the data. extreme imputation scenario (+30%), resistin was associated with Lequesne index. Thus, this indicated that resistin levels were probably not correlated with Lequesne index.

FGF-23 and ghrelin were not associated with BMI or Lequesne index

Neither FGF-23 nor ghrelin was associated with BMI, metS, or Lequesne index in complete case analysis and multiple imputation (Suppl. Table S8, , Suppl. Tables S9-A and S9-B).

Serum probes from the 2 BMI groups were analyzed for their metabolite composition and summarized in a heat map ( ). While, on the heatmap, there is clearly a pattern, we were not able to detect an association with BMI or metS.

Figure 4.

Heat map of metabolite concentrations based on BMI group (dark green = underweight, light green = obese). While there is obviously a pattern in metabolite concentrations, we could not detect an association with BMI group or metabolic syndrome in our study. Numbers at the bottom = sample number; BMI = body mass index; IMG = inositol monophosphate; GMP = guanosine monophosphate; GTP = guanosine-5-triphosphate; NADP = nicotinamide adenine dinucleotide phosphate.

To understand the relationship of metabolites with clinical parameters, PCA, group comparison with t tests, and PLS or PLS-DA were performed.

PCA results did not indicate a possible difference in metabolites of BMI groups 1 (underweight) and 3 (obese; see ). In univariate analysis, we did not detect any differences. Accordingly, in PLS-DA, we were not able to predict BMI group membership (R2Y = 0.96, Q2Y = 0.27, PR2Y = .27, PQ2 = .043).

Figure 5.

PCA score plot of the first 2 principal components. Samples were colored according to their BMI group. There is an overlap of BMI groups 1 and 3. “t1, t2” = first 2 principal components (explained variance) and deviation from the respective component (axis); s1-s31 = patient samples; PCA = principal component analysis; BMI = body mass index.

Using PLS and PLS-DA, we were also not able to predict scores of the Lequesne index or metabolic syndrome (Suppl. Tables S10-S12).

Discussion

In this study, we investigated several candidate biomarkers, as well as the metabolome, in the context of obesity, in end-stage KOA. A cross-sectional clinical/lab study was conducted in end-stage KOA patients using serum biomarkers as measured by ELISA/Luminex technique or NMR spectroscopy to conduct a multivariable and a multivariate analysis for possible associations with clinical parameters such as BMI and the Lequesne index. Our hypothesis was that biomarker levels change with respect to BMI and that biomarkers were associated with KOA severity.

We found that higher FABP4 and leptin were significantly associated with obesity and female gender. In agreement with previous studies, we did not detect an association of leptin and the Lequesne index, but BMI group and gender. Resistin did not show an association with Lequesne index in end-stage KOA patients. Finally, we were not able to detect BMI-dependent changes in the metabolome of KOA patients.

FABP4 is one of several members of a family of lipid chaperons. It is mainly expressed in adipocytes and is involved in intracellular fatty acid transportation. There is a strong correlation between FABP4 and obesity and FABP4 expression is higher in females than males. In mice models of induced obesity, FABP4 knock-out and pharmaceutical inhibition alleviated cartilage degeneration. In KOA patients, FABP4 concentrations are higher compared with non-KOA controls and are associated with obesity and KOA severity on the Kellgren-Lawrence scale. In agreement with all these findings, our results confirm higher FABP4 levels in obese females in KOA, independent of age, HKAA, and metabolic syndrome. We did not detect an association with the Lequesne index, however. Only a homogeneous group of end-stage KOA patients on the day before knee arthroplasty was included in our study and this might explain why no such association could be observed.

PLA2G2A has been implicated in the pathogenesis of OA as a cartilage degrading enzyme. To our knowledge, no study has previously assessed PLA2G2A as a blood biomarker for KOA. Our study did not demonstrate an association of PLA2G2A with BMI and/or KOA severity.

Hyperglycemia leads to an accumulation of advanced glycation end products (AGE), proteins modified by nonenzymatic glycosylation that cause local inflammation by binding to membrane-bound receptors (mRAGE). The interaction of AGEs with soluble receptors (sRAGE) is associated with anti-inflammatory responses by blocking interaction with mRAGE. Our results did not indicate an association with BMI or KOA severity. A study on fasting in KOA also found no association of changes in pain and function scores and AGE/sRAGE.

oxLDL, the oxidized form of LDL-cholesterol, is a significant regulator of cartilage degeneration : In KOA, oxLDL levels are higher compared with controls. Furthermore, in another study, oxLDL was found to be positively correlated with KOA severity. This study did not adjust for confounders, however. Our study did not demonstrate such a relationship, neither with BMI nor with KOA severity, after adjusting for several important confounders. Furthermore, in contrast to the aforementioned study, only end-stage KOA patients were included for this investigation and differences in KOA severity might be minimal.

Leptin is mostly produced by adipocytes and leptin levels are positively associated with obesity and female gender.[39,40] Leptin levels are increased in KOA compared with controls and there is an association with BMI and gender.[41 -43] Our results confirm that obese female patients have higher serum leptin levels. In agreement with a previous study, we could not detect an association with KOA severity.[42,43] This contrasts with another study including patients at different stages of KOA and might be explained by a more rigorous adjustment for confounding covariates in our study and the fact that we included end-stage KOA patients only.

In KOA patients with joint effusion, resistin is associated with severity even after adjusting for several confounders in serum and synovial fluid.[43,44] Another case-control study in obese, female patients found higher resistin levels in KOA patients and resistin levels were positively correlated with WOMAC scores. Another study in different KOA stages did not find an association of pain and resistin levels in serum, but synovial fluid only. In contrast to aforementioned study in female patients with joint effusion, in our study resistin levels were weakly, inversely associated with the Lequesne index in complete case analysis. This negative correlation was observed in complete case analysis and, in the most extreme scenario, in multiple imputation only. Our results are supported by a study showing that serum resistin levels are not associated with preoperative visual-analog scale (VAS) pain levels. Although not significant, in that study the association was also a negative one. Similarly, a previous report found that serum resistin levels were not associated with cartilage damage in end-stage KOA. It is unclear why there was a negative association in complete case analysis, if any, when compared with studies of KOA at different stages showing a positive association. It is possible that missing data introduced a bias in our data that was mitigated through multiple imputation. To summarize, our results show that resistin is not associated with Lequesne index.

In KOA, ghrelin was found to be associated with Western Ontario and McMaster Universities Osteoarthritic Index (WOMAC), even after adjusting for BMI, age, and gender. We could not confirm this finding, possibly due to more rigorous adjustment applied in our study and due to the inclusion of end-stage KOA patients only.

FGF-23 plays an important role in chondrocyte differentiation. FGF-23 is upregulated in osteoarthritic cartilage compared with normal controls. FGF-23 levels were found to be higher in FGF-23 KOA patients compared with normal patients in serum and there was a correlation of KOA severity and FGF-23 level. In contrast, in our study we could not detect an association of serum FGF-23 level with KOA severity nor BMI. The aforementioned study did not control for confounding factors and this might explain the different findings.

Previous studies, using a targeted or untargeted metabolomics approaches, showed significant changes, particularly in branched-chain amino acids (BCAA) and arginine pathways, in KOA compared with control in agreement with epidemiologic studies indicating that changes in BCAA serum concentrations are associated with diabetes, cardiovascular disease, and metabolic syndrome.[15,50] In another study, obese KOA patients showed higher oxidative stress levels compared with nonobese KOA patients. In contrast, in our study we did not detect any BMI-dependent changes in end-stage KOA patients. Again, the reason might be that only end-stage KOA patients were included, there was no control group, and differences in metabolites may be minimal within this group. Furthermore, while the abovementioned study showed BMI-dependent changes in KOA patients, that study used an untargeted approach with mass spectroscopy and these differences in the methodology might explain the discrepancies compared with our findings.

Limitations

This study has the following limitations: Only a small sample size was enrolled in this study. Due to resource limitations, there were missing data. This could be mitigated by multiple imputation and sensitivity analysis, as demonstrated, however. Furthermore, blood samples were taken at random and patient fasting state was unknown. This might influence ELISA, Luminex, and NMR spectroscopy results. However, based on the agreement with previous studies, we are confident that our results are valid. Finally, our data are applicable to patients undergoing knee arthroplasty for end-stage KOA only.

Conclusion

This study investigated the relationship of several previously not extensively examined biomarkers with BMI and KOA severity. We could confirm the association of BMI and gender with FABP4 and leptin. None of the explored biomarkers was associated with KOA severity and no BMI-dependent changes in the metabolome of end-stage KOA patients were observed.

Click here for additional data file.

Supplemental material, sj-docx-1-car-10.1177_19476035211069251 for The Association of Blood Biomarkers and Body Mass Index in Knee Osteoarthritis: A Cross-Sectional Study by Paul Schadler, Birgit Lohberger, Bettina Thauerer, Martin Faschingbauer, Werner Kullich, Martin Helmut Stradner, Andreas Leithner, Valentin Ritschl, Maisa Omara and Bibiane Steinecker-Frohnwieser in CARTILAGE