Objectives We aimed to evaluate the association between fibre intake and its food sources, and the risk of total knee replacement (TKR) due to severe knee osteoarthritis (KOA).
Methods We used data from the Singapore Chinese Health Study, a prospective cohort study that recruited 63 257 participants aged 45–74 years from 1993 to 1998. At baseline, we assessed diet using a validated 165-item semiquantitative food frequency questionnaire, together with body mass index (BMI) and lifestyle factors. Incident TKR cases were identified via record linkage with nationwide hospital discharge database through 2017.
Results There were 2816 cases of incident TKR due to severe KOA. The total fibre intake at baseline was not associated with the risk of TKR after adjustment for confounders. Among the food sources of fibre, higher intake of legumes was associated with a lower risk of TKR in a dose-dependent manner; compared with those having the lowest quartile intake, HR (95% CI) was 0.86 (0.76, 0.96) for those having the highest quartile intake (p for trend=0.004). This association was consistent after including BMI in the model and homogeneous across BMI categories. The consumption of other fibre sources, namely grain products, nuts and seeds, soy food, fruits and vegetables, was not associated with the risk of TKR.
Conclusion Intake of legumes, but not total fibre, was associated with a reduced risk of TKR. Further research is needed to replicate our findings and to evaluate possible biological mechanisms that could explain the effect of dietary legumes on pathogenesis or progression of KOA.
Data availability statement
Data are available upon reasonable request. All data are available upon reasonable request with email to corresponding author.
This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.
Statistics from Altmetric.com
What is already known about this subject?
The association between diet and development of knee osteoarthritis (KOA) may be confounded or mediated by body mass index, the latter being a strong risk factor.
Recent finding that dietary fibre may reduce the risk of KOA in Western populations needs to be replicated in leaner Asian populations.
What does this study add?
In a relatively lean Asian population, total fibre intake was not associated with the risk of total knee replacement (TKR) for severe KOA.
Among the food sources of fibre, higher intake of legumes was associated with a lower risk of TKR in a dose-dependent manner.
How might this impact on clinical practice or further developments?
Further studies are needed to replicate this finding and also to study biological plausibility underlying the effect of dietary legumes on development of osteoarthritis.
Knee osteoarthritis (KOA) is the most prevalent arthritis affecting 10% of the general population and up to 30% of elderly.1 2 KOA causes knee pain, swelling and stiffness, leading to activity limitation, participation restriction, emotional disturbance and eventual loss of independence. It is one of the leading causes of disability worldwide.3 4 Despite the huge impact on affected individuals and the society, there is a dearth of effective and safe pharmacological treatment,5 and total knee replacement (TKR) is often the only option for end-stage severe KOA6 Obesity is the major risk factor for KOA, associated with a 2-fold to over 10-fold increase in the risk of KOA.7 8 Apart from the increase in biomechanical loading, obesity may also increase systemic inflammation, which is thought to play an important role in the pathoetiologies of KOA.9 Dietary fibre is the edible parts of plants that are resistant to digestion in the small intestine and partially fermented in the colon. Dietary patterns rich in food sources of fibre, such as vegetables, legumes and fruits, have been associated with a lower risk of KOA and higher cartilage thickness.10 11 Clinical trials have demonstrated the beneficial effect of dietary fibre on weight loss and reducing systemic inflammation.12 13 Thus, dietary fibre may have biological effects on the preservation of cartilage and reduction in the risk of KOA.
Recently, data from the Osteoarthritis Initiative (OAI) of 4796 participants and from the Framingham Offspring Osteoarthritis Study of 1268 persons in the USA showed that dietary fibre was associated with a lower risk of incident symptomatic KOA after adjustment for body mass index (BMI).14 Specifically, in the OAI study, the mean BMI of those in the lowest quartile intake of fibre (29.4 kg/m2) was 1.8 kg/m2 higher than those in the highest quartile (27.6 kg/m2).14 BMI is a well-established major risk factor of KOA,7 8 15 and while overweight or obesity obviously increases biomechanical loading onto weight-bearing joints, chronic inflammation associated with obesity has also been hypothesised to be an additional mechanism for the biological plausibility linking increased BMI to the development of KOA.16 As those with higher fibre intake had lower body weight in the OAI study, the authors later reported that after controlling for BMI, this inverse association between fibre and incident symptomatic KOA, though still reaching statistical significance, was attenuated. Hence, they concluded that this inverse association was partially mediated by BMI.
As diet is integrally related to BMI, the association between diet and KOA may indeed be confounded or mediated by BMI in observational studies conducted in populations with a wide BMI range. As BMI is such a strong risk factor of KOA, it is difficult to eliminate residual confounding effect attributed by adjusting for BMI with statistical modelling alone, and this is particularly problematic when there is a wide disparity in BMI among the ordinal categories of dietary intake.
In this study, we evaluated the association between fibre intake and the risk of TKR due to severe KOA in a relatively lean population nested in the Singapore Chinese Health Study. In a leaner population with a lower variability in BMI, differences in the mean BMI by dietary intake will be comparatively smaller, thereby reducing the confounding effect of BMI on the diet–KOA risk association. In addition, we evaluated the individual association between specific food sources of fibre in this population, namely grain products, legumes, nuts and seeds, soy food, fruits and vegetables, and the risk of TKR.
Research design and methods
The Singapore Chinese Health Study is a prospective cohort of 63 257 Singaporean Chinese, aged 45–74 years during recruitment between 1993 and 1998 in Singapore. Participants were recruited from public housing estates, where 86% of Singapore’s population lived during the time of recruitment. Study participants were restricted to two major dialect groups in Singapore, the Hokkien and Cantonese, who originated from the Fujian and Guangdong Provinces in Southern China, respectively.
Baseline exposure assessment
Information at baseline was collected by face-to-face interviews using structured questionnaires.17 Information collected included age at recruitment, sex, education level, alcohol consumption, cigarette smoking, habitual physical activity, hours of sitting, sleeping hours and comorbidities. Habitual dietary intake of the past year was collected using a validated 165-item semiquantitative Food Frequency Questionnaire (FFQ) that covered food in the following categories: rice, noodles, bread and cereals, meats (red meat comprising pork, beef and mutton; poultry comprising chicken and duck; and seafood comprising fish and shellfish), soy foods, legumes, nuts and seeds, fruits, vegetables, dairy products, beverages, condiments and preserved foods. For major food items, participants were asked to report their habitual intake with eight food-frequency categories (from never or hardly ever to two or more times a day) and three portion sizes that were illustrated using actual eating utensils and accompanying photographs. Daily nutrient and energy intakes were computed from the Singapore Food Composition Table that was developed specifically for this cohort.17 The FFQ was subsequently validated in a subset of 810 participants by performing a repeat administration of the FFQ and comparing results from two 24-hour recalls with their results from the FFQ. The validation study using these two methods showed similar distributions, with the majority of the mean values of energy and nutrients obtained from the FFQ and from the 24-hour recalls falling within 10% of each other. For energy and nutrients, the correlation coefficients between these two methods ranged from 0.24 to 0.79,17 which were comparable to those from previous validation studies in diverse populations.18
In this Chinese cohort in an Asian country, the main contributors of fibre were grain products (30.6%), fruits and related juices (27.6%), vegetables and related juices (21.0%), soy (3.4%), nuts and seeds (1.2%) and legumes (0.8%), which accounted for 84.6% of the total dietary fibre. The grain products included noodles and pasta, rice, bread and pancakes, breakfast cereals, biscuits and crackers. Compared with the percentage contribution from food sources of dietary fibre in three US-based studies,14 19 the contributions from grains, nuts and legumes were considerably lower in this Singapore Chinese cohort. Instead, contributions from fruits, vegetables and other miscellaneous food sources were higher in our cohort. In the validation study, the mean daily intakes of total fibre assessed by 24-hour recalls and the FFQ were 12.2 g and 12.5 g, respectively. The correlation coefficients between FFQ-based and 24-hour recall-based intakes for dietary fibre (kcal-adjusted) were 0.72, 0.67, 0.66 and 0.65 among Cantonese men, Cantonese women, Hokkien men and Hokkien women, respectively. Furthermore, the intake of fibre in this cohort was comparable to the intake level of the participants in the Third National Health and Nutrition Examination Survey data in the USA.17
BMI was calculated from weight and height, which were self-reported at baseline interview. A total of 10 349 cohort subjects (16%) did not report either weight or height. For participants with either missing weight or height, missing values were imputed using this equation: weight=y-intercept+gradient×height. The values for y-intercept and gradient were obtained from sex-specific weight–height regression lines drawn from all participants in the cohort with known weights and heights. If both weight and height were missing, the missing height was assigned as the sex-specific median and the missing weight was calculated from the equation.20
Identification of incident cases of severe KOA
We first identified TKR cases in this cohort and only included those that were first-occurrence as incident cases. We then used the diagnosis codes to include primary KOA and exclude cases with diagnosis codes that indicated or suggested that the TKR was done for other diseases or injuries of the knee joint. We identified participants who underwent TKR due to severe KOA via record linkage with the nationwide MediClaim System hospital discharge database.21 22 Although it is used for subsequent financial claims, and hence its name, this database, regardless of the claim status, has captured surgical procedures and up to three diagnoses per patient for hospital admissions in all public and private hospitals in Singapore since 1990. The Singapore Ministry of Health conducts regular check on the MediClaims Database to ensure its completeness.23 We identified cases of TKR in this cohort using the operation codes SB010K, SB012K and SB013K (unilateral right, unilateral left and bilateral TKR) for linkage up to December 2011 and operation codes SB716K and SB809K (TKR unilateral and TKR bilateral) for linkage from January 2012 to December 2017. We only included information for the first TKR in cases with more than one episode. Among the TKR cases, we further included only cases with a diagnosis of KOA using International Classification of Diseases (ICD)-9 code 715 and all subcategories for linkage up to December 2011, and ICD-10 code M17 and all subcategories for linkage from January 2012 to December 2017. We excluded all cases with any mention of diagnoses that included septic arthritis, osteomyelitis, villonodular synovitis, rheumatoid arthritis, psoriatic arthritis, ankylosing spondylitis and other inflammatory arthritis, or a secondary cause of KOA, such as avascular or aseptic necrosis of joint, meniscus or ligament injuries, old tear or rupture or injuries, delay in development, other acquired deformities of knee and neuritis. Although we did not ask for the history of knee injury in our study, we excluded TKR cases with osteoarthritis secondary to other underlying conditions, including severe injury, if this was mentioned in the medical records, so that we could focus on TKR done for primary KOA. After excluding 128 participants who had undergone TKR prior to recruitment, we identified 2816 cases of incident TKR due to severe KOA in this cohort. Deaths in this cohort were ascertained through record linkage analysis with the nationwide Singapore Registry of Births and Deaths. As of 31 December 2011, only 47 (0.07%) subjects from this cohort were known to be lost to follow-up due to migration out of Singapore or for other reasons ascertained from our follow-up interviews in this cohort.
We calculated person-years for each participant from date of recruitment until TKR, time of death, lost-to-follow-up or 31 December 2017, whichever came first. We used the residual method to derive energy-adjusted intake of food and fibre for each participant.24 We used Cox proportional hazards regression to estimate the HRs and 95% CIs for the association of quartiles of total fibre intake, as well as food sources of fibre, with the risk of TKR. Proportional hazards assumption was tested using Schoenfeld residuals test, and no violation was observed. The linear trend was tested by assigning participants the median value of the quartile they belonged to and treating this as a continuous variable in the model.24 25 As total energy intake is often related to disease risk because of their associations with body size and physical activity, we further adjusted for total energy intake as a covariate in a multivariate nutrient density model to address the potential problem of confounding by total energy intake in the nutrient density approach.24 In model 1, we adjusted for sex, age at interview, years of interview (1993–1995, 1996–1998), total energy intake (kcal/day), dialect (Hokkien, Cantonese) and educational level (none, primary, secondary and more). In model 2, we further adjusted for BMI (kg/m2), total physical activity duration (hours/week), sleep duration (hours/day) and sitting duration (hours/day), smoking status (never, past, current) and a history of hypertension, diabetes, heart attack and stroke (yes/no). These covariates were chosen because they either are known risk factors of KOA in the literature15 or have been reported to affect the risk of TKR in our cohort.8 21 26 27 Finally, to investigate the association between individual food sources of fibre and the risk of TKR, we added a model 3 that included mutual adjustment for all the other food sources.
We tested for interaction between dietary variables and BMI or sex by adding their cross-product terms into regression models. For BMI, we used the cut-point of 23 kg/m2 to divide participants into lean and heavy categories. This cut-point was recommended for the definition of cardiovascular risk for Asians by the WHO.28 Finally, sensitivity analysis was conducted by excluding participants with imputed BMI (10 349 participants), those with an extreme energy intake of <600 or >3000 kcal/day for women (n=584) and <700 or>3700 kcal/day for men (n=476), and those with a follow-up time of less than 5 years (333 TKR cases and 3495 non-cases). We also conducted an analysis that only included cases with TKR done 5–10 years after recruitment and the rest of the cohort participants with 5–10 years of follow-up. To further examine for the confounding effect of BMI and comorbidities separately, we conducted sensitivity analyses to compare the risk estimates from models with and without the inclusion of BMI or comorbidities. We also used the method described by Lange et al29 in mediation analysis to estimate the direct effects in the association between dietary fibre/legumes and the risk of TKR, as well as the indirect effects mediated by BMI. We conducted all statistical analyses using SAS V.9.4 (SAS Institute). All reported p values were two-sided, and p<0.05 was considered statistically significant in the models for total fibre. To adjust for multiple testing, the Bonferroni-adjusted p<0.008 was considered statistically significant in models that studied the association between six fibre food sources and the risk of KOA.
After a median follow-up of 19.6 years, 2816 cases of incident TKR due to severe KOA were identified among 63 129 participants included in the analysis (figure 1). The mean (SD) age at TKR was 55.7 (6.9) years, and women accounted for 82% of TKR cases (table 1). Participants who had TKR had higher BMI (25.0 (3.0) kg/m2) at baseline compared with the rest of the cohort (23.0 (3.2) kg/m2). Participants with TKR had higher prevalence of hypertension and smoking, but lower prevalence of diabetes than those without TKR. Participants who had TKR participated in less physical activity, sitting for longer duration, but sleeping about the same amount of time compared with those without TKR.
Table 2 shows the comparison of the demographic, lifestyle and dietary factors by quartiles of total fibre intake. Compared with those with the lowest quartile, participants with the highest quartile of total fibre intake were more likely to be women, have higher education and abstain from smoking. The mean BMI was very similar among quartiles of total fibre intake, being 23.0 kg/m2 for those in the lowest quartile and 23.2 kg/m2 for those in the highest quartile. The characteristics of participants by various food sources of dietary fibre are summarised in online supplemental tables 1 and 2.
In model 1, which included adjustment for sociodemographic variables and total energy intake, the total fibre intake was associated with a stepwise increase in the risk of TKR (p for trend=0.003); compared with those in the lowest quartile, the HR (95% CI) for the risk of TKR for participants in the highest quartile of total fibre intake was 1.18 (1.05 to 1.32) (table 3). In model 2, after further adjustment for BMI and other lifestyle factors, the risk estimates in the higher quartiles of total fibre intake were attenuated, and the associations were no longer statistically significant in any of the higher quartiles or in the overall trend. The HRs (95% CIs) for second, third and fourth quartiles of total fibre intake were 1.06 (0.95 to 1.19), 1.09 (0.97 to 1.22) and 1.11 (0.99 to 1.24) compared with the lowest quartile.
In the evaluation of various food sources of fibre, we noted a stepwise reduction in the risk of TKR with increasing intake of legumes; compared with the lowest quartile of legume intake (referent), HR (95% CI) was 0.98 (0.86 to 1.09) for the second quartile, 0.95 (0.84 to 1.06) for the third quartile and 0.86 (0.76 to 0.96) for the highest quartile (p for trend=0.004) in model 3 (table 3). There was initially a stepwise reduction in the risk of TKR with increasing intake of nuts and seeds in model 1, but the association was no longer statistically significant after adjustment with other confounders and food sources of fibre. The other sources of fibre, namely grain products, soy food, fruits and vegetables, were not significantly associated with the risk of TKR.
In the analysis stratified by BMI or sex, the intake of total fibre was not associated with the risk of TKR in any stratum, and the risk estimates were not dissimilar by sex or BMI categories (p for interaction ≥0.32). Conversely, the inverse association between the intake of legumes and the risk of TKR remained statistically significant in women and also in those with a BMI ≥23 kg/m2. Although the inverse associations in men and in those with BMI <23 kg/m2 were not statistically significant, likely due to small case numbers, the differences in risk estimates between men and women, or between the two BMI categories, were not statistically significant (p for interaction ≥0.53) (online supplemental table 3).
In sensitivity analyses that excluded participants with imputed BMI (online supplemental table 4), those with extreme energy intake (online supplemental table 5) and those with less than 5 years of follow-up (online supplemental table 6), we also did an analysis including only participants with 5–10 years of follow-up (online supplemental table 7); the results remained materially unchanged compared with the findings based on the whole cohort.
In the sensitivity analysis comparing models with and without adjustment for BMI, the risk estimates in these two models were essentially unchanged, indicating that BMI did not have a significant mediation effect on the fibre–TKR association. Specifically, comparing extreme quartiles, the HR (95% CI) for the association between total fibre and risk of TKR was 1.11 (0.99 to 1.24) in both models, whereas the HR (95% CI) for the association between legumes and risk of TKR was 0.82 (0.73 to 0.92) in the model without BMI and 0.86 (0.76 to 0.96) in the model that adjusted for BMI. Predictably, in the mediation analysis, there was no significant indirect effect mediated by BMI for the association between intake of fibre or legumes and the risk of TKR (data not shown). Finally, we repeated the analysis using models with and without comorbidities (hypertension, diabetes, heart attack and stroke), and the risk estimates from these two models were materially the same (online supplemental table 8).
In this prospective cohort of a relatively lean population, we did not observe any association between total dietary fibre and risk of TKR due to severe KOA after adjusting for BMI and other confounders. However, among the food source of fibre, higher intake of legumes was associated with a lower risk of TKR in a dose-dependent manner. These results remained robust to different sensitivity analyses.
The strength of our study is that the BMI levels across quartiles of fibre intake were quite homogeneous among the relatively leaner participants in this prospective cohort study in Singapore. Hence, expectantly, the associations between total fibre or food intake and risk of TKR were essentially unchanged in the models with and without BMI adjustment, thus indicating that any mediation or confounding effect from BMI was minimal. Another strength of our study is its prospective population-based cohort design with a large sample size over a long follow-up time. TKR cases were identified via linkage with a comprehensive nationwide database that had essentially captured a high number of incident TKR cases for severe KOA in this cohort.
To our best knowledge, this is the first study to specifically report an inverse relationship between legume intake and risk of TKR due to severe KOA. The association remained the same after adjusting for BMI and other food sources of fibre. Although Dai et al also reported a dose-dependent inverse association between nut and legumes and incident symptomatic and incident radiographic KOA from the Framingham Offspring cohort,14 nuts and legumes were studied as a single item. In this study, we studied nuts and legumes as separate food items, and although we observed an inverse association of nuts and seeds with the risk of TKR in the preliminary multivariable model, this association was not statistically significant in the fully adjusted model. On the contrary, the inverse association between legumes and the risk of TKR was robust and remained statistically significant in the fully adjusted model.
Legumes are edible seeds or pods of the plant family Fabaceae (or Leguminosae). Examples of common legumes are beans, peas, chickpeas, lentils and lupin beans. Legumes are a rich source of saponin, phenolic compounds, phytochemicals, vitamins and minerals that have antioxidant properties30–32 and anti-inflammatory effects that could possibly reduce low-grade systemic inflammation33 and activate mitochondrial oxidation.34 The generation of reactive oxygen species (ROS) and lipid peroxidation products has been an implication in the pathogenesis of the arthritic joint.35 Hence, it is possible that phytochemicals in legumes could reduce ROS production and salvage mitochondria dysfunction36 37 to slow down progression of osteoarthritis. However, further research is dreadfully needed to evaluate these hypotheses.
We acknowledge a few limitations in the study. First, we studied TKR as a surrogate for severe KOA, and although it represented a relevant clinical state of KOA when surgical intervention is indicated due to symptomatic, functional and structural severity in addition to pain,38 we could have missed severe cases of KOA that did not undergo TKR. Second, dietary intakes, BMI and comorbidities were self-reported, and measurement errors may likely lead to non-differential misclassification and an underestimation of the observed association. Nonetheless, self-reported BMI has generally been reported to be valid for epidemiological studies across different populations.39 Third, diet and BMI were only assessed once at baseline; therefore, we lacked information about potential changes in these exposures that could have occurred during the follow-up period. Although FFQs tend to overestimate and underestimate the absolute intake of various nutrients and foods, for most nutrient and foods, the ability of the FFQ to rank subjects according to their dietary intake has been shown to be acceptable to good, and this is the basis of all epidemiological studies that have used FFQ to study diet–disease relationships.40 There was a 64-fold difference in the median intake of legumes between the lowest and the highest quartile in the Singapore Chinese Health Study cohort that was sufficient to demonstrate the inverse legume–TKR association in a dose-dependent manner. Finally, residual confounding cannot be completely ruled out in our study due to the limitation of the observational design.
In conclusion, in this prospective cohort study, we provided novel epidemiology evidence that although total fibre intake was not associated with severe KOA, increased intake of legumes could possibly reduce the risk of severe KOA in a dose-dependent manner. Further research is needed to replicate our findings, and the biological mechanisms underlying the effect of dietary legumes on pathogenesis or progression of osteoarthritis warrant further investigations.
Data availability statement
Data are available upon reasonable request. All data are available upon reasonable request with email to corresponding author.
The Institutional Review Board of the National University of Singapore read and approved the study protocol. All participants signed informed consents prior to recruitment.
Presented at Manuscripts based on work previously presented at a conference and published as a conference abstract in EULAR 2021.
Contributors YYL and WPK designed and conducted the research; WKP, JMY, LWA and KBT collected the data; AJ and WPK analysed data; YYL and WPK wrote the paper; and WPK had primary responsibility for final content. All authors critically reviewed and approved the final manuscript.
Funding This study was supported by the National Institutes of Health, USA (NIH R01 CA144034 and UM1 CA182876). YYL and WKP were supported by the National Medical Research Council, Singapore (NMRC/CSA-Inv/0022/2017 and MOH-CSASI19nov-0001, respectively).
Competing interests None declared.
Provenance and peer review Not commissioned; externally peer reviewed.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.