Skip to main content
Download PDF

Continuous care intervention with carbohydrate restriction improves physical function of the knees among patients with type 2 diabetes: a non-randomized study

BMC Musculoskeletal Disorders volume 23, Article number: 297 (2022) Cite this article

Abstract

Background

In a previous study, we assessed a novel, remotely monitored carbohydrate restricted diet regimen including nutritional ketosis in patients with type 2 diabetes and reported significant improvements in weight, glycemic control, abdominal fat and inflammation from baseline to 2 years. Knee outcome measures were collected as a secondary outcome in the trial. This study aims to assess the effect of this intervention on knee functional scores and to identify if changes in weight, central abdominal fat (CAF), glycemic status and high sensitivity C-reactive protein (hsCRP) were associated with its improvement.

Methods

This prospective analysis included continuous care intervention (CCI, n = 173) and usual care (UC, n = 69) trial participants with type 2 diabetes that reported knee pain at baseline. Knee outcome measures included the Knee injury and Osteoarthritis Outcome Score (KOOS) pain, symptoms, activities of daily living (ADL), sports and recreation function, and knee-related quality of life subscales, and total KOOS score were assessed from baseline to 2 years. Missing data at each time point were replaced with multiple imputation under the assumption of missing at random. To assess if the primary analysis of the knee scores changed under plausible missing not at random assumptions, sensitivity analysis was also performed using pattern mixture models. In CCI, we also assessed factors associated with the improvement of knee scores.

Results

In the primary analysis, CCI participants demonstrated a statistically significant improvement in total KOOS and all KOOS individual subscale scores at 1 year and maintained through 2 years as opposed to UC patients who showed no significant changes from baseline to 2 years. The significant improvement in total KOOS and its individual subscale scores from baseline to 2 years remained relatively stable in CCI in the sensitivity analysis under different missing not at random scenarios confirming the robustness of the findings from the primary analysis. Approximately 46% of the CCI participants met the 10 points minimal clinically important change at 2 years. A reduction in CAF was associated with improvement in total KOOS and KOOS ADL, while a decrease in hsCRP was associated with improvement in KOOS symptoms scores.

Conclusion

A very low carbohydrate intervention including nutritional ketosis resulted in significant improvements in knee pain and function among patients with T2D. The improvements in knee function were likely secondary to a reduction in central adiposity and inflammation. Future research on the applicability of this intervention in radiographically confirmed OA patients is important.

Trial registration

Clinical trial registration: NCT02519309 (10/08/2015).

Peer Review reports

Background

Osteoarthritis (OA) is a painful and debilitating disease process affecting roughly 54.4 million adults (22.7%) in the United States annually [1]. Knee OA is the most common form and is particularly morbid in its effects, leading to decreased mobility, difficulty performing self-care and performing usual activities of daily living [2]. Among adults with diabetes, the prevalence of OA increases to an estimated 48.1% [3]. Although the exact mechanistic link between diabetes and OA is poorly understood, diabetes is recognized as a substantial predisposing risk factor for OA and increased OA severity [4,5,6]. Type 2 diabetes is a chronic metabolic disorder mainly characterized by hyperglycemia and insulin resistance. It is suggested that the persistent high level of glucose in those with diabetes alter connective tissue extracellular matrix, induces cell injury and damages synovium leading to musculoskeletal disorder [5, 7]. Likewise, low-grade systemic inflammation is observed in both diabetes and OA, and it is possibly linked to the pathophysiology of these two conditions [5]. In addition to its predisposing effects on OA, diabetes is also an independent amplifier of knee osteoarthritic pain [6, 8], and its effect on pain is independent of weight [9, 10]. However, the impact of glycemic control on knee pain is debatable with some reporting an association between elevated HbA1c and knee pain [11] while others reporting no evidence of association [12].

Among nonsurgical treatment options for knee OA, weight loss is routinely recommended and has been shown to significantly improve knee pain and function [13,14,15]. In fact, weight loss is one of the few nonpharmacologic treatments recommended by the American Academy of Orthopaedic Surgeons in their most recent guidelines for treatment of the osteoarthritic knee [13]. The Intensive Diet and Exercise for Arthritis (IDEA) trial demonstrated a 40% reduction in symptoms related to knee OA with a > 10% decrease in overall bodyweight [16]. It is likely that the effect of weight loss seen in these studies is due to a decrease in the knee joint loading. Specifically, thigh and abdominal fat content, in addition to overall body mass have been linked to significantly increased knee loading force [17]. Additional studies utilizing magnetic resonance imaging have also demonstrated that weight loss reduces actual overall chondral degeneration over time in patients with knee OA [18]. Likewise, the pivotal Look AHEAD trial found that weight loss leads to decreased knee pain and improved patient reported knee function scores even in overweight patients without a confirmed diagnosis of OA [19].

While the relationship between bodyweight and knee osteoarthritis has been well established, more recent studies have explored diet as an independent modulator of knee OA. Preliminary studies suggest that diet may be an independent variable influencing knee pain in osteoarthritic mice [20]. Furthermore, ketogenic diet modification has been proposed as a means of lowering overall body inflammation [21], which has been shown to be beneficial in reducing many types of inflammatory pain [22]. Among proposed mechanisms for the anti-inflammatory effects of ketosis is the inhibitory effect beta-hydroxybutyrate plays in the NLRP3 inflammasome-mediated inflammatory response [23]. Carbohydrate restriction and ketogenic diets are also widely used in the clinical management of obesity and diabetes [24,25,26]. Carbohydrate restriction to less than 30 to 50 g per day is often necessary to induce nutritional ketosis, defined as blood beta-hydroxybutyrate in the range of 0.5 to 3.0 mM [27,28,29,30]. Physiologically during nutritional ketosis, insulin secretion is reduced to a level that promotes fat oxidation and utilization of ketones as the main energy supply [30].

While ketogenic diets have been shown to be successful in treating diabetes and improving some forms of pain [22, 24,25,26], and diabetes may be an independent amplifier of knee pain [6, 8], there are currently no studies evaluating the potential independent effects of a carbohydrate restricted diet on knee pain and function in diabetes. Our recent study in patients with type 2 diabetes following a remotely monitored carbohydrate restricted diet regimen inducing nutritional ketosis reported significant improvements in weight and glycemic control from baseline to 2 years while deprescribing antihyperglycemic medication [27,28,29]. We also reported a significant reduction in inflammation at 2 years [27, 29]. Subjective data on patient reported knee pain and function scores were collected from these patients as a secondary outcome at each assessment time point. The knee assessment was performed with the Knee Injury and Osteoarthritis Outcome Score (KOOS) questionnaire. The KOOS is a validated patient-reported metric designed to assess the short-term and long-term symptoms and function of patients with knee injury and/or osteoarthritis [31,32,33]. It is a multidimensional questionnaire measuring five knee-related domains including pain, symptoms, activities of daily living, sports and recreation function and knee-related quality of life.

We postulate that a carbohydrate restriction intervention would improve knee pain and function scores and it is possibly driven by the improvement in weight, glycemic control, and inflammation. Therefore, in this study, we aim to evaluate the effect of this carbohydrate restriction intervention on subjective knee pain and function scores that were secondarily collected as part of the trial. In addition, we also assessed and identified factors that were associated with the knee improvements in the intervention. For the association analysis, we included outcome measures that were previously reported and commonly associated with knee pain and OA [16,17,18,19].

Methods

Study design

The full detail of the study design, inclusion and exclusion criteria and the results from primary outcomes have been previously published [32, 33]. This study was approved by the Franciscan Health Lafayette Institutional Review Board and registered under Clinical trials.gov identifier NCT02519309. Eligible participants were aged between 21 and 65 years with a diagnosis of T2D and a BMI > 25 kg/m2 and all participants provided written informed consent. The trial was an open-label, non-randomized, controlled, longitudinal study with patients divided into two groups: remotely monitored continuous care intervention (CCI) and usual care (UC). In the CCI regimen, participants self-selected to receive educational content either through on-site group meetings (CCI-onsite) or virtual delivery (CCI-virtual). The medical treatment and educational content were the same between the CCI-onsite and CCI-virtual, except the CCI-onsite group had additional clinic group meetings that ranged from weekly for 12 weeks to quarterly in the second year (26–28). All participants in the CCI arm were given access to a mobile health application (app) that allowed them to communicate with a team of healthcare professionals consisting of a personal health coach and health care provider (physician or nurse practitioner). These patients were educated and guided remotely to achieve nutritional ketosis, which typically included a daily dietary carbohydrate restriction to less than 30 g. Patients were also encouraged to measure body weight daily and input fingerstick blood glucose and beta hydroxybutyrate concentrations in the app allowing for continuous monitoring of patients’ diet adherence, progress and to guide their healthcare team to make appropriate adjustments. Since there were no significant differences in the knee outcomes including total KOOS and its individual subscale scores between the two education groups, the CCI-onsite and CCI-virtual groups were combined for the analyses in this study.

The usual care group was recruited separately. This group consisted of T2D patients who began participating in a local diabetes education program with treatment by their primary care physician or endocrinologist and diet counseling by registered dietitians. No modification was made to their care for the study. This group served as a reference of diabetes disease progression under routine care within the same geographical area and health system.

Study participants in knee assessment

The recruitment of patients into the study were mainly based on the inclusion and exclusion criteria listed for the primary outcomes. There were 262 CCI and 87 UC participants with T2D enrolled at baseline for the main study. The study was not designed for assessing knee outcome; the knee subjective data were collected as a secondary outcome in the trial along with other questionnaires assessing sleep and depression symptoms. The KOOS questionnaires were administered to participants who completed their laboratory assessments by the study coordinator during their baseline, 1- and 2-years follow-up visit in the clinic. The patients completed paper-based questionnaires after their clinic visit. At baseline, there were some patients who failed to complete the questionnaire or sections of the questionnaire. At 1 and 2 years, in addition to patients not completing the questionnaire, missingness in the knee questionnaire data was also related to patients withdrawing from the study, and patients failed to attend their clinic follow-up visit after completion of their laboratory assessment.

For the knee assessment, participants who reported no knee pain at baseline (based on their response to question P1 “How often do you experience knee pain” of the KOOS questionnaire) were excluded from analysis and those who reported knee pain at baseline or missing any baseline data were retained for further analysis (Supplementary Fig. 1, 173 CCI; 141 with data and 32 missing, 69 UC; 52 with data and 17 missing). Missing data at each assessment time point were multiple imputed for the final analysis.

Outcome measures

Demographic and clinical variables

Both CCI and UC participants underwent comprehensive laboratory analysis and clinical anthropometric assessment at baseline, 1 year, and 2 years. Details of the methods were previously published [27, 29]. Blood analytes were processed at a laboratory certified by the Clinical Laboratory Improvement Amendments (CLIA). CCI participants were further assessed for body composition changes at baseline, 1 and 2 years using dual X-ray absorptiometry (DXA) (Lunar GE Prodigy, Madison, WI) analyzed using GE Encore software (v11.10, Madison, WI). The data assessing safety, effectiveness, and sustainability of the two intervention arms (CCI and UC) on glycemic status, diabetes medication use, weight and other metabolic parameters were previously published. Some of these clinical variables that significantly improved in the CCI arm from baseline to 2 years [27, 29] including body weight, body composition (central abdominal fat), glycemic status (HbA1c and fasting glucose) and inflammation (high sensitivity C-reactive protein and white blood cells) were included in this study for a post-hoc association analysis in the CCI group.

Knee injury and Osteoarthritis Outcome Score (KOOS)

Clinical assessment of knee symptoms was performed with the Knee injury and Osteoarthritis Outcome Score (KOOS) survey. The KOOS is a multi-dimensional questionnaire measuring five knee-related domains including pain, symptoms, activities of daily living, sports and recreation function, and knee-related quality of life. The KOOS and its individual sub-scores were used as a subjective measure to gauge patient knee pain and function with or without a diagnosis of knee osteoarthritis. Each section in the KOOS survey is scored on a normalized 0–100 scale (with 100 categorized as the best with no knee problems and 0 as the worse knee outcome). The total KOOS score was calculated as the averaged individual KOOS subscale scores (0–100). The Cronbach’s α coefficients between the KOOS individual sub-scores and total KOOS scores were excellent in all three time points (> 0.95).

Statistical analyses

All data analyses were conducted using IBM SPSS Statistics for Windows, version 26.0 (IBM Corp, Armonk, N.Y., USA). We assessed the normality and linearity of the knee KOOS data including KOOS Symptoms, KOOS Pain, KOOS Activities of Daily Living (KOOS ADL), KOOS Sports and Recreation (KOOS Sport), KOOS Quality of Life (KOOS QoL) subscores and total KOOS score from all three time points using skewness and kurtosis [34].

Descriptive statistics of the baseline characteristics were performed, and independent sample student t-tests were used to analyze differences in baseline characteristics between CCI and UC, and completers and non-completers (including dropouts). The difference in the proportion of missing data at each time point between CCI versus UC was assessed using Chi-Square test.

Due to high degree of missingness in the knee outcome data. The primary assessment of KOOS total and subscores was performed under missing at random (MAR) assumption using multiply imputed dataset. A standard multiple imputation (MI) was performed to generate 50 imputed datasets. For within-group longitudinal changes from baseline to 2 years and between-group differences (CCI vs UC) at 1 and 2 years in the KOOS total and sub-scores, we conducted linear mixed-effects model (LMMs) with time as a repeated measure. The models included time, group effect (CCI vs UC) and time by group interaction as fixed effects and several baseline and demographic covariates (as listed in the Table 1 footnote). A maximum likelihood-approach was used for parameter estimation and the correlations between repeated measures were accounted for using an unstructured covariance. The LMMs were carried out for each multiple imputed datasets as an intent-to-treat analysis and the results were combined using Rubin’s rules. A second completers only analysis was also performed on participants with available data at baseline and 2 years.

Table 1 Adjusted mean changes over time in participants with knee pain at baseline [CCI T2D versus UC T2D]
Full size table

To further investigate the possibility of the data missing not at random (MNAR), several sensitivity analyses were conducted using pattern-mixture models (PMMs) applying different scenarios of plausible departures from missing at random (MAR) assumption. The imputed values in the 50 MI datasets were multiplied by a sensitivity parameter delta (δ) representing the average difference between the observed values collected from the completers and those missing any values. For the total KOOS and subscores data, we reported a significant increase of at least 10% from baseline in CCI in the primary analysis using MI data with MAR assumption. Based on the primary analysis, we estimated a difference of 5 to 10% is considered as a reliable δ to assess the different MNAR scenarios. A comprehensive step-by-step detailed description of the sensitivity analysis using PMMs is included in the Supplementary Method Section.

As an additional exploratory analysis only in the CCI group, we then assessed predictors associated with KOOS improvement. This analysis was performed among those who completed the study with available knee data at 2 years. We first performed unadjusted Pearsons’ correlations between changes in total KOOS and its subscale scores with changes in glycemic, weight, abdominal fat and inflammation at 1 and 2 years. Variables with significant correlations were then included in multivariate linear mixed-effects regression models with random intercepts and random time slopes. Six different models were generated to test the association between the total KOOS score and its individual subscale scores with these potential predictors. In these models, the time trend was controlled, and the predictors were modeled as a time-varying covariate. The models also included predictors by time interaction term to assess concurrent association over time between the predictor variables and KOOS. The models also included several baseline and demographic covariates (as listed in the Table 2 footnote).

Table 2 Regression coefficients (β), 95% confidence intervals and p-values for each predictor from linear mixed effects regression models associated with total and individual KOOS subscale scores
Full size table

Additionally, to further explore clinical improvement, we assessed the efficacy of the interventions using minimal clinically important change (MIC) at 2 years among those who completed the study with available knee data at baseline and 2 years. We utilized a cut-off of 10 (31,33) for MIC; those achieving a > + 10 increase in total KOOS score from baseline to 2 years were classified as “improvement” and those having a lower score or < 10 increase in score were classified as “no improvement”. Proportions of participants with improvement and no improvement were calculated in both CCI and UC and between-group differences were analyzed using logistic regression adjusted for age, gender, race, and diabetes duration. Factors associated with likelihood of achieving “improvement” versus “no improvement” in CCI were assessed with logistic regression.

For the primary aim of the paper assessing longitudinal changes in the total and individual subscale scores, a Bonferroni adjusted p-value =0.05/6 variables = 0.008 was used to determine statistical significance. For all other exploratory analyses, a p-value = 0.05 was used to determine statistical significance.

Results

Baseline participants characteristics

Patient baseline characteristics for knee scores and other demographic, metabolic and inflammatory variables among CCI and UC participants reporting knee pain at baseline are provided in supplementary Table 1. There was a total of 173 CCI and 69 UC participants reporting knee pain at baseline. Participant groups were consistent with each other in their baseline demographic, biochemical variables, and total KOOS and subscale scores, except for body weight and white blood cell count (Supplementary Table 1). However, using standardized differences there were 8 variables showing imbalance between CCI and UC participants including age, gender, diabetes duration, weight, white blood cell count, KOOS pain, KOOS activities of daily living and KOOS quality of life. There were no significant differences between completers and non-completers in the baseline total KOOS and subscale scores in CCI and UC, except that KOOS QoL was significantly greater in non-completers versus completers in UC.

Missing data

Missing data were encountered in all assessment time points even at baseline in both the intervention arms (Supplementary Fig. 1). The greatest level of missingness was observed at 2 years where only 59.5% of CCI and 53.6% of UC participants had complete KOOS data. Missing data at 2 years were mainly due to patients who had withdrawn from the study (28.9% in CCI and 23.2% in UC) and patients missed the follow-up clinic appointment or patients failed to complete the questionnaire (12.1% in CCI and 23.2% in UC). The distribution of missing data was consistent between CCI and UC at baseline (18.5% in CCI vs 27.5% in UC) and 2 years (40.5% in CCI and 46.4% in UC) but at 1-year, missing data were more common in CCI than UC (35.8% in CCI vs 17.4% in UC, p = 0.003).

Primary post-intervention analysis in total KOOS and subscale scores

Among participants with knee pain at baseline (CCI, n = 173; UC, n = 69), CCI demonstrated a statistically significant improvement in total KOOS and all KOOS individual subscale scores at 1 year compared to UC patients [Table 1; total KOOS increased at 1 year = 9.9 (95% CI: 6.9, 12.5) in CCI, decreased in UC = − 0.5 (95% CI: − 4.1, 5.0)]. These improvements remained significant from baseline at 2 years [Table 1; CCI = 10.0 (95% CI: 7.4, 13.5) increase, UC = − 3.7 (95% CI: − 8.7, 0.9) decrease) with no significant differences between 1 and 2 years. There were approximately 10% increase in total KOOS and its subscale scores at 1 and 2 years (Fig. 1) in CCI. Similar improvements were not observed in UC (Fig. 1). The overall mean improvement in total KOOS and its subscale scores in CCI at 1 and 2 years also reached the minimal clinically important change of 8–10 points, except for KOOS Symptoms (7.0 at 1 year and 7.5 at 2 years).

Fig. 1
figure 1

Percentage changes in the total and its individual subscale Knee Injury and Osteoarthritis Scores (KOOS) from baseline to 2 years in Continuous Care Intervention and Usual Care

Full size image

The CCI group also demonstrated significantly higher scores compared to their UC counterparts at 1 year (p < 0.01) and 2 years (p < 0.01) in total KOOS as well as in all KOOS subscales (Table 1). Results of completers only were consistent with the intent-to-treat analysis.

MNAR sensitivity analysis

The longitudinal changes in KOOS total and its subscores under different MNAR modified missing data assumption are listed in Supplementary Table 2. The PMM sensitivity analysis demonstrated that the longitudinal changes in KOOS total and its subscores remain relatively stable in both CCI and UC under different MNAR assumption (MAR departures). In all the MAR departure scenarios, the CCI showed significant improvements in total KOOS (Fig. 2A) and its subscores from baseline to 2 years (Supplementary Table 2) but not in UC (Fig. 2B, Supplementary Table 2) as originally demonstrated in the primary MI dataset under MAR assumption. For example, even when the missing values were estimated to be 10% lower in CCI and 10% higher in UC, we observed similar trend of overtime change in total KOOS and its subscores in both CCI and UC. This confirms that the missing pattern between the two treatment groups CCI and UC are similar where both groups equally experience MNAR bias that cancels out when we performed the sensitivity analysis. This is further corroborated by the equal distribution of missing data or attrition between CCI and UC at baseline and 2 years.

Fig. 2
figure 2

Longitudinal change in total Knee Injury and Osteoarthritis Scores (KOOS) from baseline to 2 years under different MNAR scenarios. A. Continuous Care Intervention. 1: Primary imputed values under missing at random assumption. 0.95: Primary imputed values estimated 5% lower. 0.90: Primary imputed values estimated 10% lower. B. Usual Care. 1: Primary imputed values under missing at random assumption. 0.95: Primary imputed values estimated 5% lower. 0.90: Primary imputed values estimated 10% lower. 1.10: Primary imputed values estimated 10% higher

Full size image

Identifying factors associated with knee score improvement in CCI

Unadjusted Pearson correlation revealed significant inverse correlations between change in weight, central abdominal fat (CAF) and hs C-reactive protein (hsCRP) with change in total KOOS and several of the KOOS subscale scores. Changes in HbA1c, fasting glucose, and white blood cell were not correlated with changes in total KOOS and its individual subscale scores.

We then performed multivariate linear mixed-effects regression models that included all three potential predictors: weight, CAF and hsCRP with random intercepts and random time slopes. CAF was the only factor significantly and negatively associated with improvement of total KOOS (β= − 5.0, p = 0.04) and KOOS ADL (β= − 4.8, p = 0.04) scores, where a decrease in CAF was associated with higher total KOOS and KOOS ADL scores (Table 2). There were no significant associations between weight, CAF, or hsCRP with change in KOOS pain, sport and quality of life. However, for KOOS symptoms, hsCRP was the only factor significantly associated with its improvement with β= − 0.9 (p-0.03), where a decrease or lower hsCRP was associated with higher KOOS Symptom score (Table 2).

KOOS Minimal Clinically Important Change (MIC)

At 2 years, 45.5% of the CCI participants and 14.8% of UC participants were classified as total KOOS score ‘improvement”. As an exploratory analysis, the odds of achieving the MIC in total KOOS was 5.4 times (95% CI: 1,7, 17.9, p = 0.01) greater in CCI than UC and greater reduction in both CAF and hsCRP at 2 years were associated with greater likelihood of achieving the MIC (CAF, OR = 0.11, p = 0.02 and hsCRP, OR = 0.79, p = 0.03) in CCI.

Discussion

The purpose of this secondary analysis was to evaluate the effect of a novel, remotely monitored carbohydrate restricted dietary intervention including nutritional ketosis on subjective knee pain and function scores in patients with type 2 diabetes (T2D), and to identify factors associated with improvement. In this paper, we also conducted sensitivity analysis using pattern mixture-models (PMMs) to assess the issue of missingness in our data and the likelihood of the data missing not at random (MNAR) [35, 36]. Missingness is a common issue encountered in patient reported outcome measures especially when the data is collected as a secondary endpoint [37] and if paper-based data collection is used [38, 39].

Our primary analysis results suggest that a medically supervised carbohydrate restricted diet and continuous care intervention statistically and clinically improved total KOOS and all its individual subscores at 1 year and that these improvements are largely maintained out to 2 years among patients with knee pain. In contrast, knee function and pain in the usual care arm did not change significantly from baseline showing a gradually worsening trend commonly observed in patients with T2D [4, 5, 40]. In accordance with the primary analysis, our MNAR sensitivity analysis using PMMs also confirmed that the continuous care intervention showed a statistically significant longitudinal improvement in total KOOS and all its individual subscores from baseline to 2 years even when we assumed a differential attrition pattern between CCI (those missing data had worse outcomes) and UC (those missing data had better outcomes) [36, 41]. The sensitivity analysis confirms that the primary post-intervention analysis for total KOOS and its subscores can be interpreted as robust to MNAR assumption and the overall interpretation of the primary analysis did not change [36, 41]. Our findings were consistent with a recent study that reported a significant decrease in pain scores for patients with knee OA via a low carbohydrate diet compared to a low fat and a routine diet [42].

Compared to prior publications, while our analysis demonstrates that weight reduction correlated with improvement in overall KOOS by Pearson correlation, changes in weight alone were not predictive of improvement in overall KOOS by multivariate linear mixed-effects regression. Instead, improvements in CAF and hsCRP were each independent predictors associated with significant improvements in overall KOOS at 2 years. Other studies have demonstrated a strong correlation between weight loss and improved patient reported knee pain and function scores in overweight and T2D patients [19]. The Look AHEAD trial demonstrated an intensive lifestyle intervention (ILI) involving diet/exercise education and supervision led to significant weight loss compared to usual care in a T2D patient population [43]. The study showed statistically significant improvement in Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) knee scores in participants as well as a decrease in overall inflammation as measured by hsCRP [21, 44]. The IDEA trial showed similar results in WOMAC scores and reported reduction in inflammation measured by IL-6 levels in an overweight and obese population without diabetes [16]. The current study confirms and builds upon these findings in that it demonstrates similar improvement in patient reported knee pain and function (KOOS) at 1 and 2 years as well as reaffirms a potential role of reduced inflammation, in this case measured by hsCRP. Our study also further demonstrates that CCI participants were not only able to achieve statistically significant improvement in knee pain and function but also approximately 46% of the CCI participants had a clinically meaningful improvement of over 10 points in total KOOS score [31]. Additionally, CCI participants were more likely to achieve this improvement than their UC counterparts. Achieving clinically meaningful improvements suggests that CCI participants may be more likely to increase physical activities [31], and hence achieve additional benefits from fitness associated with alleviating knee pain and improving knee function [19, 45].

While many studies have demonstrated improved knee pain with reduced body weight [15, 16, 19], this study is unique because, in addition to demonstrating a correlation with weight loss and improved overall KOOS score, it identifies central abdominal fat as the only anthropometric measure to predict improvement in total KOOS scores. This coincides with recent findings by Pan et al. noting a correlation between greater central obesity and increased knee pain among individuals with metabolic syndrome [46]. Furthermore, our results were in accordance with other studies that also assessed body composition indices [17, 47]. These studies revealed that changes in body fat or percentage body fat explained a greater proportion of variance in change in knee pain and function than weight loss alone [17, 47]. Given the fact that weight and adiposity are strongly correlated with each other, it is possible that studies reporting weight loss as a powerful predictor of knee improvement [15, 16, 19] may have also encountered reduction in body fat and central adiposity. Although a precise mechanistic link between central adipose tissue and knee pain has yet to be delineated, a number of mechanisms have been proposed for the relationship between obesity and pain [48]. Adipose tissue is also known to play an important role in producing proinflammatory cytokines. This presents a possible mechanism for increased knee pain and a potential explanation for the predictive nature of decreased abdominal fat and reduced serum hsCRP in improved knee pain and function [49, 50]. Additional studies are needed to further characterize the associations between adipose tissue, inflammatory markers, and knee pain.

The strengths of this study are prospective analysis, long term follow-up in a relatively large sample size, and extensive laboratory and anthropometric analysis at multiple time points. Among its potential weaknesses, this is a nonrandomized study performed in a single healthcare system which both represent potential areas for unintentional bias. Therefore, the interpretation of the primary findings was limited to the longitudinal change in total KOOS and its subscores in CCI and UC. The knee score assessment was also limited by the level of missingness in the data. There are potential unmeasured confounders that are not accounted for in the sensitivity analysis and the true missing data trajectory cannot be confirmed from the data. Additional exploratory aims were included to understand the impact of the intervention on knee pain and scores mainly for future study planning. Furthermore, there is an absence of information regarding radiographically confirmed diagnoses of OA in the study subjects. This challenges the direct applicability of our findings to patients with OA as we do not know which of the participants have radiographic OA and which do not. However, given that the causes of knee pain in this population with a high initial BMI are varied and often present in the absence of radiographically significant OA, this study does provide clinical application in that it demonstrates improved generic knee pain and function via a ketogenic diet without regard to an underlying diagnosis.

Conclusions

A continuous care intervention including nutritional ketosis was associated with significant improvements in knee pain and function in patients with T2D. While these improvements were observed concurrently with weight loss, this study suggested that the change in knee function were possibly mediated by reduction in central adiposity and inflammation. The findings from this study prompt future research on the applicability of this intervention in radiographically confirmed OA patients.

Availability of data and materials

The datasets analyzed in the current study are available from the corresponding authors on reasonable request.

Abbreviations

T2D:

Type 2 diabetes

CCI:

Continuous care intervention

UC:

Usual care

OA:

Osteoarthritis

KOOS:

Knee Injury and Osteoarthritis Outcome Score

CAF:

Central abdominal fat

hsCRP:

High sensitivity C-reactive protein

IDEA:

Intensive diet and exercise for arthritis

BMI:

Body mass index

CLIA:

Clinical laboratory improvement amendments

DXA:

Dual X-ray absorptiometry

MIC:

Minimal clinically important change

ILI:

Intensive lifestyle intervention

WOMAC:

Western Ontario and McMaster Universities Osteoarthritis Index

References

  1. Barbour KE, Helmick CG, Boring M, Brady TJ. Vital signs: prevalence of doctor-diagnosed arthritis and arthritis-attributable activity limitation-United States, 2013-2015. MMWR Morb Mortal Wkly Rep. 2017;66:246–53.

    Article  Google Scholar 

  2. Clynes MA, Jameson KA, Edwards MH, Cooper C, Dennison EM. Impact of osteoarthritis on activities of daily living: does joint site matter? Aging Clin Exp Res. 2019;31:1049–56. https://doi.org/10.1007/s40520-019-01163-0.

    Article  PubMed  PubMed Central  Google Scholar 

  3. Cheng YJ, Imperatore G, Caspersen CJ, Gregg EW, Albrright AL, Helmick CG. Prevalence of diagnosed arthritis and arthritis-attributable activity limitation among adults with and without diagnosed diabetes: United States 2008-2010. Diabetes Care. 2012;35:1686–91.

    Article  Google Scholar 

  4. Eymard F, Parsons C, Edwards MH, Petit-Dop F, Reginster JY, Bruyere O, et al. Diabetes is a risk factor for knee osteoarthritis progression. Osteoarthr Cartil. 2015;23:851–9.

    Article  CAS  Google Scholar 

  5. King KB, Rosenthal AK. The adverse effects of diabetes on osteoarthritis update on clinical evidence and molecular mechanisms. Osteoarthr Cartil. 2015;23:841–50.

    Article  CAS  Google Scholar 

  6. Schett G, Kleyer A, Perricone C, Sahinbegovic E, Iagnocco A, Zwerina J, et al. Diabetes is an independent predictor for severe osteoarthritis: results from a longitudinal cohort study. Diabetes Care. 2013;36:403–9.

    Article  Google Scholar 

  7. Laiguillon MC, Courties A, Houard X, Auclair M, Sautet A, Capeau J, et al. Characterization of diabetic osteoarthritic cartilage and role of high glucose environment on chondrocyte activation: toward pathophysiological delineation of diabetes mellitus-related osteoarthritis. Osteoarthr Cartil. 2015;23:1513–22.

    Article  Google Scholar 

  8. Eitner A, Pester J, Vogel F, Marintschev I, Lehmann T, Hofmann GO, et al. Pain sensation in human osteoarthritic knee joints is strongly enhanced by diabetes mellitus. Pain. 2017;158:1743–53.

    Article  Google Scholar 

  9. Eitner A, Culvernor AG, Wirth W, Schaible HG, Eckstein F. Impact of diabetes mellitus on knee osteoarthritis pain and physical and mental status: data from the osteoarthritis initiative. Arthritis Care Res (Hoboken). 2021;73:540–8.

    Article  Google Scholar 

  10. Alenazi AM, Alshehri MM, Alothman S, Alqahtani BA, Rucker J, Sharma N, et al. The association of diabetes with knee pain severity and distribution in people with knee osteoarthritis using data from the Osteoarthristis initiative. Sci Rep. 2020;10:3985. https://doi.org/10.1038/s41598-020-60989-1.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Charen DA, Solomon D, Zubizaretta N, Poeran J, Colvin AC. Examining the association of knee pain with modifiable cardiometabolic risk factors. Arthritis Care Res. 2020. https://doi.org/10.1002/acr.24423.

  12. Garessus ED, de Mutsert R, Visser AW, Rosandaal FR, Kloppenburg M. No association between impaired glucose metabolism and osteoarthritis. Osteoarthr Cartil. 2016;24:1541–7.

    Article  CAS  Google Scholar 

  13. American Academy of Orthopaedic Surgeons: Treatment of Osteoarthritis of the Knee-3rd Edition: Evidence-Based Guideline. https://www.aaos.org/globalassets/quality-and-practice-resources/osteoarthritis-of-the-knee/oak3cpg.pdf. Published Aug 31, 2021.

  14. Taylor N. Nonsurgical Management of Osteoarthritis Knee Pain in the older adult. Clin Geriatr Med. 2017;33:41–51.

    Article  Google Scholar 

  15. Alrushud AS, Rushton AB, Kanavaki AM, Greig CA. Effect of physical activity and dietary restriction interventions on weight loss and the musculoskeletal function of overweight and obese older adults with knee osteoarthritis: a systematic review and mixed method data synthesis. BMJ Open. 2017;7:e014537. https://doi.org/10.1136/bmjopen-2016-014537.

    Article  PubMed  PubMed Central  Google Scholar 

  16. Messier SP, Mihalko SL, Legault C, Miller GD, Nicklas BJ, DeVita P, et al. Effects of intensive diet and exercise on knee joint loads, inflammation, and clinical outcomes among overweight and obese adults with knee osteoarthritis: the IDEA randomized clinical trial. JAMA. 2013;310:1263–73.

    Article  CAS  Google Scholar 

  17. Messier SP, Beavers DP, Loeser RF, Carr JJ, Khajanchi S, Legault C, et al. Knee joint loading in knee osteoarthritis: influence of abdominal and thigh fat. Med Sci Sports Exerc. 2014;46(9):1677–83.

    Article  Google Scholar 

  18. Gersing AS, Schwaiger BJ, Nevitt MC, Zarnowski J, Jospeh GB, Feuerriegel G, et al. Weight loss regimen in obese and overweight individuals is associated with reduced cartilage degeneration: 96-month data from the osteoarthritis initiative. Osteoarthr Cartil. 2019;27:863–70.

    Article  CAS  Google Scholar 

  19. Foy CG, Lewis CE, Hairston KG, Miller GD, Lang W, Jakicic JM, et al. Intensive lifestyle intervention improves physical function among obese adults with knee pain: findings from the look AHEAD trial. Obesity (Silver Spring). 2011;19:83–93.

    Article  Google Scholar 

  20. Loredo-Perez AA, Montalvo-Blanco CE, Hernandez-Gonzalez LI, Anaya-Reyes M, Fernandez Del Valle-Laisequilla C, Reyes-Garcia JG, et al. High-fat diet exacerbates pain-like behaviors and periarticular bone loss in mice with CFA-induced knee arthritis. Obesity (Silver Spring). 2016;24:1106–15.

    Article  CAS  Google Scholar 

  21. Forsythe CE, Phinney SD, Fernandez ML, Quann EE, Wood RJ, Bibus DM, et al. Comparison of low fat and low carbohydrate diets on circulating fatty acid composition and markers of inflammation. Lipids. 2008;43:65–77. https://doi.org/10.1007/s11745-007-3132-7.

    Article  CAS  PubMed  Google Scholar 

  22. Masino SA, Ruskin DN. Ketogenic diets and pain. J Child Neurol. 2013;28:993–1001.

    Article  Google Scholar 

  23. Youm YH, Nguyen KY, Grant RW, Goldberg EL, Bodogai M, Kim D, et al. The ketone metabolite β-hydroxybutyrate blocks NLRP3 inflammasome-mediated inflammatory disease. Nat Med. 2015;2:263–9.

    Article  Google Scholar 

  24. Yamada Y, Uchida J, Izumi H, Tsukamoto Y, Inoue G, Watanabe Y, et al. A non-calorie-restricted low-carbohydrate diet is effective as an alternative therapy for patients with type 2 diabetes. Intern Med. 2014;53:13–9.

    Article  CAS  Google Scholar 

  25. Saslow LR, Daubenmier JJ, Moskowitz JT, Kim S, Murphy EJ, Phinney SD, et al. Twelve-month outcomes of a randomized trial of a moderate-carbohydrate versus very low-carbohydrate diet in overweight adults with type 2 diabetes mellitus or prediabetes. Nutr Diabetes. 2017;7:304. https://doi.org/10.1038/s41387-017-0006-9.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Volek JS, Phinney SD, Forsythe CE, Quann EE, Wood RJ, Puglisi MJ, et al. Carbohydrate restriction has a more favorable impact on the metabolic syndrome than a low fat diet. Lipids. 2009;44:297–309.

    Article  CAS  Google Scholar 

  27. Hallberg SJ, McKenzie AL, Williams PT, Bhanpuri NH, Peters AL, Campbell WW, et al. Effectiveness and safety of a novel care model for the management of type 2 diabetes at 1 year: an open-label, non-randomized, controlled study. Diabetes Ther. 2018;9:583–612.

    Article  Google Scholar 

  28. Bhanpuri NH, Hallberg SJ, Williams PT, McKenzie AL, Ballard KD, Campbell WW, et al. Cardiovascular disease risk factor responses to a type 2 diabetes care model including nutritional ketosis induced by sustained carbohydrate restriction at 1 year: an open label, non-randomized, controlled study. Cardiovasc Diabetol. 2018;17:56. https://doi.org/10.1186/s12933-018-0698-8.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  29. Athinarayanan SJ, Adams RN, Hallberg SJ, McKenzie AL, Bhanpuri NH, Campbell WW, et al. Long-term effects of a novel continuous remote care intervention including nutritional ketosis for the management of type 2 diabetes: a 2-year non-randomized clinical trial. Front Endocrinol (Laussanne). 2019;10:348. https://doi.org/10.3389/fendo.2019.00348.

    Article  Google Scholar 

  30. Poff A, Koutnik A, Egan B. Nutritional ketosis with ketogenic diets or exogenous ketones: features, convergence and divergence. Curr Sports Med Rep. 2020;19(7):251–9.

    Article  Google Scholar 

  31. Roos EM, Lohmander LS. The knee injury and osteoarthritis outcome score (KOOS): from joint injury to osteoarthritis. Health Qual Life Outcomes. 2003;1:64. https://doi.org/10.1186/1477-7525-1-64.

    Article  PubMed  PubMed Central  Google Scholar 

  32. Rodriguez-Merchan EC. Knee instruments and rating scales designed to measure outcomes. J Orthop Traumatol. 2012;13:1–6.

    Article  Google Scholar 

  33. Collins NJ, Misra D, Felson DT, Crossley KM, Roos EM. Measures of knee function: international knee documentation committee (IKDC) subjective knee evaluation form, knee injury and osteoarthritis outcome score (KOOS), knee injury and osteoarthritis outcome score physical function short form (KOOS-PS), knee outcome survey activities of daily living scale (KOS-ADL), Lysholm knee scoring scale, Oxford knee score (OKS), Western Ontario and McMaster universities osteoarthritis index (WOMAC), activity rating scale (ARS), and Tegner activity score (TAS). Arthritis Care Res. 2011;63(Suppl 11):S208–28.

    Article  Google Scholar 

  34. Kline RB. Principles and practice of structural equation modeling. 3rd ed. New York: The Guilford Press; 2011.

    Google Scholar 

  35. Permutt T. Sensitivity analysis for missing data in regulatory submissions. Stat Med. 2016;35:2876–9. https://doi.org/10.1002/sim.6753.

    Article  PubMed  Google Scholar 

  36. Ratitch B, O'Kelly M, Tosiello R. Missing data in clinical trials: from clinical assumptions to statistical analysis using pattern mixture models. Pharm Stat. 2013;12(6):337–47. https://doi.org/10.1002/pst.1549.

    Article  PubMed  Google Scholar 

  37. Gomes M, Gutacker N, Bojke C, Street A. Addressing missing data in patient-reported outcome measures (PROMS): implications for the use of PROMS for comparing provider performance. Health Econ. 2016;25(5):515–28. https://doi.org/10.1002/hec.3173.

    Article  PubMed  Google Scholar 

  38. Hagan TL, Belcher SM, Donovan HS. Mind the mode: differences in paper vs. web-based survey modes among women with Cancer. J Pain Symptom Manag. 2017;54(3):368–75. https://doi.org/10.1016/j.jpainsymman.2017.07.005.

    Article  Google Scholar 

  39. Ebert JF, Huibers L, Christensen B, Christensen MB. Paper- or Web-Based Questionnaire Invitations as a Method for Data Collection: Cross-Sectional Comparative Study of Differences in Response Rate, Completeness of Data, and Financial Cost. J Med Internet Res. 2018;20(1):e24. https://doi.org/10.2196/jmir.8353. Published 2018 Jan 23.

    Article  PubMed  PubMed Central  Google Scholar 

  40. Neumann J, Hofmann FC, Heilmeier U, Ashmeik W, Tang K, Gersing AS, et al. Type 2 diabetes patients have accelerated cartilage matrix degeneration compared to diabetes free controls: data from the osteoarthritis initiative. Osteoarthr Cartil. 2018;26:751–61.

    Article  CAS  Google Scholar 

  41. Cro S, Morris T, Kenward M, Carpenter J. Sensitivity analysis for clinical trials with missing continuous outcome data using controlled multiple imputation: a practical guide. Stat Med. 2020;39. https://doi.org/10.1002/sim.8569.

  42. Strath LJ, Jones CD, George AP, Lukens SL, Morrison SA, Soleymani T, et al. The effect of low-carbohydrate and low-fat diets on pain in individuals with knee osteoarthritis. Pain Med. 2020;21(1):150–60.

    Article  Google Scholar 

  43. Look AHEAD Research Group, Wadden TA, West DS, Delahanty L, Jakicic J, Rejeski J, et al. The look AHEAD study: a description of the lifestyle intervention and the evidence supporting it. Obesity (Silver Spring). 2006;14:737–52.

    Article  Google Scholar 

  44. Belalcazar LM, Reboussin DM, Haffner SM, Hoogeveen RC, Kriska AM, Schwenke DC, et al. Look AHEAD research group. A 1-year lifestyle intervention for weight loss in individuals with type 2 diabetes reduces high C-reactive protein levels and identifies metabolic predictors of change: from the look AHEAD (action for health in diabetes). Diabetes Care. 2010;33:2297–303.

    Article  Google Scholar 

  45. Kraus VB, Sprow K, Powell KE, Buchner D, Bloodgood B, Piercy K, et al. Effects of physical activity in knee and hip osteoarthritis: a systematic umbrella review. Med Sci Sports Exerc. 2019;51:1324–39.

    Article  Google Scholar 

  46. Pan F, Tien K, Cicuttini F, Jones G. Metabolic syndrome and trajectory of knee pain in older adults. Osteoarthr Cartil. 2020;28:45–52.

    Article  CAS  Google Scholar 

  47. Toda Y, Toda T, Takemura S, Wada T, Morimoto T, Ogawa R. Change in body fat, but not body weight or metabolic correlates of obesity, is related to symptomatic relief of obese patients with knee osteoarthritis after a weight control program. J Rheumatol. 1998;25:2181–6.

    CAS  PubMed  Google Scholar 

  48. Narouze S, Souzdalnitski D. Obesity and chronic pain: systematic review of prevalence and implications for pain practice. Reg Anesth Pain Med. 2015;40:91–111.

    Article  Google Scholar 

  49. Urban H, Little CB. The role of fat and inflammation in the pathogenesis and management of osteoarthritis. Rheumatology (Oxford). 2018;57:iv10–21.

    Article  CAS  Google Scholar 

  50. Azamar-Llamas D, Hernandez-Molina G, Ramos-Avalos B, Furuzawa-Carballeda J. Adipokine contribution to the pathogenesis of osteoarthritis. Mediat Inflamm. 2017;2017:5468023. https://doi.org/10.1155/2017/5468023.

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank the participants for their involvement in the study. We also appreciate the effort of the clinical care team, including health coaches and health care providers in providing the required care and treatment for the participants to achieve their dietary and health goals.

Funding

Virta Health Corp.

Author information

Authors and Affiliations

  1. Department of Orthopedics, Harbor-UCLA, 1000 W Carson St, Torrance, CA, 90502, USA

    Kade S. Lyman & John P. Andrawis

  2. Virta Health, 501 Folsom Street, San Francisco, CA, 94105, USA

    Shaminie J. Athinarayanan, Amy L. McKenzie, Camy L. Pearson, Rebecca N. Adams, Sarah J. Hallberg, Jeff S. Volek & Stephen D. Phinney

  3. Indiana University Health Arnett, Lafayette, IN, USA

    Sarah J. Hallberg

  4. School of Medicine, Indiana University, Indianapolis, IN, USA

    Sarah J. Hallberg

  5. Abbott Diabetes Care, Alameda, CA, 94502, USA

    James P. McCarter

  6. Department of Genetics, Washington University School of Medicine, St. Louis, MO, USA

    James P. McCarter

  7. Department of Human Sciences, The Ohio State University, Columbus, OH, USA

    Jeff S. Volek

  8. Torrance Memorial Physician Network, 23560 Crenshaw Blvd Unit 102, Torrance, CA, 90505, USA

    John P. Andrawis

Authors
  1. Kade S. Lyman
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Shaminie J. Athinarayanan
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Amy L. McKenzie
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Camy L. Pearson
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Rebecca N. Adams
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. Sarah J. Hallberg
    View author publications

    You can also search for this author in PubMed Google Scholar

  7. James P. McCarter
    View author publications

    You can also search for this author in PubMed Google Scholar

  8. Jeff S. Volek
    View author publications

    You can also search for this author in PubMed Google Scholar

  9. Stephen D. Phinney
    View author publications

    You can also search for this author in PubMed Google Scholar

  10. John P. Andrawis
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

KLM interpreted the data and drafted the manuscript. SJA contributed to the study analysis plan, collected data, analyzed, interpreted the data and drafted the manuscript, ALM contributed to the study design and collected data, RNA collected data, SJH contributed to the study design and collected data, JPM, JSV and SDP contributed to the study design. JPA contributed to the study analysis planning, interpreted the data, and edited the draft. All authors revised and provided final approval of the manuscript.

Corresponding author

Correspondence to John P. Andrawis.

Ethics declarations

Ethics approval and consent to participate

All procedures in this study involving human participants were performed in accordance with the ethical standards of the Institutional Board and the ethical guidelines of the 1964 Helsinki declaration and its updates. This study was approved by the Franciscan Health Lafayette Institutional Review Board, Lafayette, IN, USA and a written informed consent was obtained from all study participant and/or their legal guardian(s) to participate in the study.

Consent for publication

Not applicable.

Competing interests

Athinarayanan, McKenzie, Adams, Pearson, Hallberg, and Phinney have been employed by Virta Health Corp and were offered stock options. McCarter holds stocks from Virta Health Corp. Phinney and Volek are founders of Virta Health Corp. Andrawis and Lyman declare no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Information

Additional file 1: Supplementary Method Section. Supplementary Table 1.

Baseline characteristics of participants with knee pain in CCI and UC and comparison of completers and dropouts or missing data. Supplementary Table 2. Pattern-mixture model sensitivity analysis based on multiple imputation and different MNAR assumption for missing KOOS total and subscores data. Supplementary Figure 1. Patient flowchart from baseline to 2 years with knee score data.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Lyman, K.S., Athinarayanan, S.J., McKenzie, A.L. et al. Continuous care intervention with carbohydrate restriction improves physical function of the knees among patients with type 2 diabetes: a non-randomized study. BMC Musculoskelet Disord 23, 297 (2022). https://doi.org/10.1186/s12891-022-05258-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12891-022-05258-0

Keywords

BMC Musculoskeletal Disorders

ISSN: 1471-2474

Contact us