Haohao Zhang‡
a,
Junhao Liang‡a,
Yaxin Hana,
Jiajing Tiana,
Yahui Tua,
Rui Fanab,
Wenli Zhuab,
Zhaofeng Zhang*abc and
Haifeng Zhao*de
aDepartment of Nutrition and Food Hygiene, School of Public Health, Peking University, Haidian District, Beijing 100191, People's Republic of China. E-mail: zhangzhaofeng@bjmu.edu.cn; Fax: +86-10-82801575
bBeijing's Key Laboratory of Food Safety Toxicology Research and Evaluation, Beijing 100191, People's Republic of China
cInstitute of Medical Technology, Peking University Health Science Center, Beijing 100191, People's Republic of China
dNutritional and Food Sciences Research Institute, Department of Nutrition and Food Hygiene, School of Public Health, Shanxi Medical University, Taiyuan 030001, People's Republic of China. E-mail: haifengzao75@163.com; Tel: +86-351-3985907
eMOE Key Laboratory of Coal Environmental Pathogenicity and Prevention, School of Public Health, Shanxi Medical University, Taiyuan 030001, People's Republic of China
First published on 24th July 2025
Background. This study was aimed at investigating the association between a plant-based diet and the risk of arthritis and at identifying a strategy that achieves the ambitious goal of healthy aging. Methods. The nationwide cohort study included 10059 adults aged 65 and older from the 2008–2018 waves of the Chinese longitudinal healthy longevity survey (CLHLS). Dietary intake was collected using a simplified food frequency questionnaire and used to calculate the plant-based diet index (PDI). Arthritis was defined as participants self-reporting suffering from arthritis. Time-dependent Cox regression model was used to calculate hazard ratios (HRs) and 95% confidence intervals (CIs) for the risk of arthritis. Interaction analysis was used to explore the interaction between PDI and exercise status. Stratified analyses were used to examine factors that may modify the association. Results. During a median follow-up period of 4.2 years, 1482 participants who were free of arthritis at baseline reported arthritis. The highest quality of plant-based diet was associated with a 16.0% decrease in the risk of suffering from arthritis (HR: 0.840, 95% CI: 0.757, 0.932). Interaction analysis showed that participants with a high PDI and exercise had a significantly lower risk of arthritis. Stratified analysis showed that the association between PDI and arthritis was significant among participants living in rural areas. Conclusions. Greater adherence to a plant-based diet may help delay the onset of arthritis symptoms. Promoting plant-based dietary patterns may be a strategy to reduce arthritis incidence and improve healthy life expectancy.
The most prevalent forms of arthritis include rheumatoid arthritis (RA) (ICD-11: FA20), osteoarthritis (OA) (ICD-11: FA0Z), gouty arthritis (ICD-11: FA25), and polymyalgia rheumatica (ICD-11: FA22). According to the global burden of disease study 2021, as of 2020, the global population comprises approximately 595 million individuals with OA, 17.6 million with RA, and 55.8 million with gout. Furthermore, the incidence and prevalence of arthritis are progressively increasing over time. Globally, the total number of OA cases has increased by 132.2%, the prevalence of RA has increased by 14.1%, and the prevalence of gout has increased by 22.5% compared to those in 1990. It is projected that by 2050, 31.7 million people will have RA globally, the total number of prevalent cases of gout will reach 95.8 million, and cases of OA in the knee will increase by 74.9%, the hand by 48.6%, the hip by 78.6%, and other types of OA by 95.1%.4–6 Arthritis has emerged as a prominent contributor to disability, particularly among the elderly, and has assumed a significant role in global public health.
Despite the identification of potential signaling pathways, the current lack of a definitive pharmaceutical intervention for arthritis can be attributed to its multifaceted etiology. Extensive research has revealed a correlation between prevalent forms of arthritis and the release of specific inflammatory mediators.7 Notably, myeloid cells, particularly macrophages, assume a pivotal role in the progression of arthritis from a preclinical to a clinical state. Upon activation, macrophages secrete cytokines such as interleukin-1 (IL-1), tumor necrosis factor-α (TNF-α), and granulocyte-macrophage colony-stimulating factor. Substantial therapeutic benefits can be achieved through the inhibition of these cytokines.8 Currently, the management of arthritis predominantly relies on non-steroidal anti-inflammatory drugs, glucocorticoids, and artificial joint replacement. Nevertheless, individuals subjected to prolonged administration of these medications or surgical intervention not only endure the associated adverse effects but also shoulder a significant financial strain.
The World Health Organization is aligning its efforts with the Global Strategy and Action Plan on Ageing and Health 2016–2020, as well as the United Nations Decade of Healthy Ageing (2021–2030), in order to promote healthy ageing worldwide. The government of China introduced the 14th Five-Year Plan for Healthy Ageing. However, the advancement of healthy ageing policies is impeded by the absence of effective coping strategies for arthritis. Additionally, the limited availability of medical resources has shifted the current focus towards dietary interventions. The modulation of systemic inflammation through dietary patterns has been observed to yield positive health outcomes in various chronic diseases.9–11 A study conducted by Jun Li et al., involving a follow-up period of 5291
518 person-years, suggests that reducing the inflammatory potential of dietary patterns could be an effective approach to managing cardiovascular disease.12 Mazzucca et al. conducted a prospective analysis to examine the effect of different food intakes on RA onset through the use of the UK Biobank.13 However, it is important to acknowledge that the existing body of research on arthritis and dietary patterns has certain limitations.14–16 Most studies have relied on cross-sectional designs, which may limit the ability to establish causality. Additionally, the small sample sizes and the inclusion of populations from diverse regions with varying genetic backgrounds may restrict the generalizability of the findings to a global context.
The plant-based diet encompasses a vegetarian dietary approach that incorporates a wide range of foods such as fruits, vegetables, legumes, grains, nuts, and beans, among others.17 This dietary pattern has been extensively investigated for its numerous anti-inflammatory components, which have shown promise in various disease contexts.18 Notably, research conducted in China has demonstrated the potential of plant-based diets in mitigating the adverse effects of PM2.5 on cognitive decline and reducing the risk of all-cause mortality.19,20 Building upon these findings, this study aims to investigate whether a plant-based dietary pattern, utilizing a nationally representative database from China, can serve as a preventive measure against arthritis. Our objective was to discover a straightforward and practical approach to implementing a method that would enhance healthy life expectancy, thereby attaining the ambitious objective of promoting healthy aging.
Our study used the 2008 waves of the CLHLS and their follow-up waves in 2011, 2014 and 2018. We excluded participants who were younger than 65 years old, missed values for arthritis, missed dietary patterns and covariates, and had arthritis at baseline. In total, 10059 participants were included in the analysis. More details on participant inclusion and exclusion can be found in ESI Fig. S1.†
We also constructed a healthful plant-based diet index (hPDI) to explore whether a healthful plant-based diet had a greater influence on arthritis. Based on the plant-based diet index, salt-preserved vegetables and sugar were scored 5, 4, 3, 2 or 1. Refined grain was scored 1 and whole grain was scored 5. The total score ranged from 16 to 80. The population was also divided into two groups. More details on the construction and scoring of the PDI and hPDI can be found in ESI Table S1.†
Follow-up ended on the date of suffering from arthritis, death, loss to follow-up, or end of study period. We used the Kaplan–Meier cumulative incidence curve and the Schoenfeld residuals method to test the assumption for the proportional hazards model. Cox proportional hazard models were used to assess the association between plant-based diet patterns and the risk of developing arthritis during the follow-up. The participants were stratified by PDI scores and individual food groups to explore the association between the plant-based diet and arthritis and the association between individual foods and arthritis. Using the low PDI group as a reference, we calculated the hazard ratios (HRs) and 95% confidence intervals (CIs) of arthritis across different PDI categories. For all outcomes, models without any adjusted covariates (model-not adjusted), models adjusted only for age and sex (model 1-adjusting age and sex), models adjusted for age, sex, resident type, educational attainment, financial status (model 2-further adjusting sociological factors), models adjusted for age, sex, exercise status, smoking status, drinking status, BMI, comorbidity (model 3-further adjusting demographic variables) and models adjusted for age, sex, resident type, educational attainment, financial status, exercise status, smoking status, drinking status, BMI, comorbidity (model 4-multivariate adjustment) were constructed. We also used restricted cubic splines, which were fitted with 4 knots, to explore the nonlinearity between PDI and arthritis. The Wald test was used to assess whether the observed relationships were linear or nonlinear.
For hPDI, we also constructed models to calculate the HRs and 95% CIs of arthritis across categories of hPDI, which aimed to determine whether a healthful plant-based diet was more effective than a plant-based diet in preventing arthritis.
To ensure the stability of the results, a sensitivity analysis was performed. First, considering death as a competing risk factor for arthritis outcome, the multi-factor competitive risk model was applied as a sensitivity analysis. Therefore, multi-factor competing risks regression models were calculated. Second, our study aimed to examine the correlation between PDI and arthritis, utilizing various definitions of arthritis, such as hospital diagnosis and a combination of patient-reported symptoms and hospital diagnosis. Third, we conducted a stratified analysis based on demographic factors, including sex, type of residence, educational attainment, financial status, exercise habits, smoking habits, and alcohol consumption.
To explore the interaction between PDI and exercise status, interaction analysis was conducted. Participants were divided into four subgroups based on PDI and exercise status (do not exercise and exercise). Multi-adjusted time-dependent Cox regression was used to examine the interaction. For individual foods, multi-adjusted time-dependent Cox regression was used to calculate HRs and 95% CIs to explore the association between the frequency of each food intake and arthritis. Statistical analysis was performed using R software. Statistical significance was considered at P < 0.05 in two-sided tests.
Characteristicsb | Low PDI (N = 4474) | High PDI (N = 5585) | Total (N = 10![]() |
P-Value |
---|---|---|---|---|
a Continuous variables are presented as mean (SD) and median [min, max], and the categorical variables are presented as n (%).b Low PDI indicated participants in low PDI group, and High PDI indicates participants in high PDI group. Participants were divided into low PDI group and high PDI group based on the median PDI score of 49. PDI – plant-based diet index, BMI – body mass index. | ||||
Age (year) | <0.001 | |||
Mean (SD) | 89.4 (11.0) | 85.4 (11.3) | 87.2 (11.3) | |
Median [min, max] | 91.0 [65.0, 114] | 86.0 [65.0,116] | 89.0 [65.0, 116] | |
Sex | <0.001 | |||
Male | 1866 (41.7%) | 2652 (47.5%) | 4518 (44.9%) | |
Female | 2608 (58.3%) | 2933 (52.5%) | 5541 (55.1%) | |
Resident type | 0.013 | |||
City and town | 1594 (35.6%) | 2149 (38.5%) | 3743 (37.2%) | |
Rural | 2880 (64.4%) | 3436 (61.5%) | 6316 (62.8%) | |
BMI | <0.001 | |||
Mean (SD) | 19.9 (3.45) | 20.6 (3.58) | 20.3 (3.54) | |
Median [min, max] | 19.5 [10.8,40.8] | 20.2 [11.1,55.5] | 20.0 [10.8,55.5] | |
BMI type | <0.001 | |||
Abnormal | 2178 (48.7%) | 2460 (44.0%) | 4638 (46.1%) | |
Normal | 2296 (51.3%) | 3125 (56.0%) | 5421 (53.9%) | |
Educational attainment (year) | <0.001 | |||
<1 | 2913 (65.1%) | 3399 (60.9%) | 6312 (62.7%) | |
≥1 | 1561 (34.9%) | 2186 (39.1%) | 3747 (37.3%) | |
Finance status (Yuan) | 0.928 | |||
<10![]() |
1907 (42.6%) | 2402 (43.0%) | 4309 (42.8%) | |
≥10![]() |
2567 (57.4%) | 3183 (57.0%) | 5750 (57.2%) | |
Exercise status | 0.002 | |||
Exercise | 1649 (36.9%) | 2252 (40.3%) | 3901 (38.8%) | |
Don't exercise | 2825 (63.1%) | 3333 (59.7%) | 6158 (61.2%) | |
Smoke status | <0.001 | |||
Smoke | 1391 (31.1%) | 2119 (37.9%) | 3510 (34.9%) | |
Don't smoke | 3083 (68.9%) | 3466 (62.1%) | 6549 (65.1%) | |
Drink status | <0.001 | |||
Drink | 1298 (29.0%) | 1901 (34.0%) | 3199 (31.8%) | |
Don't drink | 3176 (71.0%) | 3684 (66.0%) | 6860 (68.2%) | |
Comorbidity | 0.995 | |||
No comorbidity | 2459 (55.0%) | 3081 (55.2%) | 5540 (55.1%) | |
1 comorbidity | 1307 (29.2%) | 1638 (29.3%) | 2945 (29.3%) | |
2 comorbidities | 501 (11.2%) | 599 (10.7%) | 1100 (10.9%) | |
≥3 comorbidities | 207 (4.6%) | 267 (4.8%) | 474 (4.7%) |
PDI | ||||
---|---|---|---|---|
Low PDI | High PDI | P-Value | P for trendf | |
a Model was not adjusted.b Model 1 was adjusted for age and sex.c Model 2 was adjusted for age, sex, resident type, educational attainment and financial status.d Model 3 was adjusted for age, sex, exercise status, smoking status, drinking status, BMI and comorbidity.e Model 4 was adjusted for age, sex, resident type, educational attainment, financial status, exercise status, smoking status, drinking status, BMI and comorbidity.f Test for trend based on a variable containing the median value for each quintile. | ||||
Range | 26–48 | 49–68 | ||
Case | 662 | 820 | ||
Person-years | 17![]() |
25![]() |
||
Case/1000 person-years | 37.7 | 32.1 | ||
Modela | 1 (ref) | 0.821 (0.741, 0.909) | <0.01 | <0.01 |
Model 1b | 1 (ref) | 0.839 (0.756, 0.931) | <0.01 | <0.01 |
Model 2c | 1 (ref) | 0.838 (0.755, 0.929) | <0.01 | <0.01 |
Model 3d | 1 (ref) | 0.841 (0.758, 0.933) | <0.01 | <0.01 |
Model 4e | 1 (ref) | 0.840 (0.757, 0.932) | <0.01 | <0.01 |
Based on the unadjusted model, participants with the highest quality of plant-based diet (PDI range: 49–68) were associated with a 17.9% decrease in the risk of arthritis (HR: 0.821, 95% CI: 0.741, 0.909), compared to those with the lowest quality (PDI range: 26–48). In both age- and sex-adjusted analyses and multivariable-adjusted analysis, the association between PDI and the risk of arthritis showed a similar trend but with a slightly diminished magnitude (HR: 0.839, 95% CI: 0.756, 0.931 in model 1; HR: 0.838, 95% CI: 0.755, 0.929 in model 2; HR: 0.841, 95% CI: 0.758, 0.933 in model 3; HR: 0.840, 95% CI: 0.757, 0.932 in model 4). Compared to a plant-based diet, a healthy plant-based diet did not show superiority (ESI Table S3†). The HRs (95% CIs) of arthritis according to continuous PDI are shown in ESI Table S2.† The Kaplan–Meier cumulative incidence curve is shown in ESI Fig. S2.†
Overall, there was no nonlinear relation between PDI and the risk of arthritis (P for non-linearity = 0.93). The risk of arthritis was decreased in participants with higher PDI (Fig. 1). When PDI was below 49, the diet pattern acted as a risk factor for arthritis. In contrast, when PDI was 49 or higher, the diet pattern acted as a protective factor for arthritis.
PDI | Exercise status | Case | Person-years | Case/1000 person-years | HR (95% CI)a | P-Value |
---|---|---|---|---|---|---|
a Models were adjusted for age, sex, resident type, educational attainment, financial status, exercise status, smoking status, drinking status, BMI and comorbidity. | ||||||
Low | Don't exercise | 394 | 10![]() |
37.3 | 1 (ref) | |
Low | Exercise | 268 | 6982 | 38.4 | 0.987 (0.841, 1.16) | 0.875 |
High | Don't exercise | 476 | 14![]() |
32.3 | 0.856 (0.748, 0.980) | 0.024 |
High | Exercise | 344 | 10![]() |
31.7 | 0.808 (0.694, 0.941) | <0.01 |
![]() | ||
Fig. 3 Individual food analyses. Models were adjusted for age, sex, residence, educational attainment, financial status, exercise status, smoking status, drinking status, BMI and comorbidity. |
Our findings suggest that a plant-based diet can delay the onset of arthritis symptoms, mainly due to its high content of polyphenols, flavonoids, antioxidant vitamins, and dietary fiber, while its low saturated fat content also plays a role.27,28 Notably, polyphenols possess the ability to counteract free radicals and inhibit the activation of NF-κB, thereby suppressing inflammatory reactions.29–31 Additionally, vegetable oil is a valuable source of monounsaturated fatty acids, such as oleic acid. Experimental studies have demonstrated that oral oleic acid can be converted into oleoylethanolamide, which may serve as a promising therapeutic agent for many inflammatory disorders in the small intestine.32 Consequently, this inhibitory effect can impede the transition of arthritis from a state of heightened susceptibility to the manifestation of clinical symptoms.8 Conversely, saturated fatty acids derived from animal sources have been observed to elevate the levels of pro-inflammatory cytokines, including TNF-α and IL-6, potentially accelerating the onset or exacerbating the severity of arthritis symptoms.33
Several researchers reported that exercise has a beneficial impact on RA.34,35 To investigate the relationship between plant-based diet and exercise, we conducted an interaction analysis. The findings indicated that enhancing the quality of a plant-based diet could reduce the risk of arthritis and enhance the synergistic effects of the plant-based diet and regular exercise through multiple mechanisms. Our results provide a scientific basis for the effectiveness of combining diet and exercise in relieving and preventing arthritis symptoms. High PDI diet and exercise are feasible and economic lifestyle interventions that have important implications for individuals with lower economic status or limited access to good healthcare. We suggest that the elderly adopt a high PDI diet combined with exercise to prevent arthritis.
In our study, we conducted a stratified analysis to identify potential factors that may modify the effect of a plant-based diet. Firstly, our findings indicate that individuals residing in rural areas with a high PDI diet experienced a significant reduction in the risk of arthritis; a pattern that was also observed among those with lower levels of education. These disparities may be attributed to the income gap between urban and rural areas, as well as to the unequal distribution of medical resources in rural regions.36 This inequality may lead to the failure in timely diagnosis and treatment of early symptoms of arthritis, thereby exacerbating the disease progression. Additionally, compared to urban regions, rural food supplies are relatively constrained. Against this backdrop, a plant-based diet, as a simple, feasible, and cost-effective intervention, plays a more pronounced role in alleviating arthritis symptoms and delaying disease progression. Similarly, participants with lower levels of education may have a simpler understanding of a healthy diet and place a greater emphasis on the accessibility and affordability of food. Therefore, such participants may be more likely to adopt a plant-based diet as a healthy lifestyle recommendation, thus demonstrating a stronger effect in reducing the risk of arthritis. Secondly, our findings indicate a significant reduction in arthritis risk among individuals with a normal BMI who adhered to a plant-based diet. However, no discernible impact was observed among those classified as underweight, overweight, or obese. This discrepancy in effects across different BMI groups may be attributed to elevated levels of inflammatory factors, which potentially counteract the beneficial effects of a plant-based diet in individuals with obesity.37 Consequently, these results underscore the importance of maintaining a healthy body weight in the elderly population. This notion aligns with the China Healthy Lifestyle for All initiative, which emphasizes the importance of maintaining a healthy body weight. Thirdly, our research indicates that adopting a plant-based diet has a greater efficacy in mitigating the susceptibility to arthritis among individuals who engage in smoking and excessive alcohol consumption. The deleterious effects of smoking and high-dose alcohol intake on inflammatory cytokine levels and inflammation exacerbation are widely recognized.38 Consequently, it is imperative for elderly individuals who partake in smoking and alcohol consumption to embrace a plant-based diet as a preventive measure against the onset of arthritis. Fourth, our study reveals a significant correlation between the adoption of a plant-based diet and decreased arthritis susceptibility among individuals with a limited number of comorbidities. This finding may potentially be attributed to the observed link between comorbidity and inflammation.39 To enhance the preventive measures against arthritis, we recommend that older adults engage in physical fitness activities and incorporate a plant-based diet into their daily regimen.
Our study possessed several notable strengths. (1) It was the first investigation to explore the correlation between PDI and arthritis risk, utilizing a nationally representative sample of older adults with a substantial sample size. (2) We meticulously accounted for various covariates, performed multiple stratified analyses, and conducted sensitivity analyses to guarantee the robustness and reliability of our findings. (3) We presented a locally feasible, economically viable, and socially acceptable approach to dietary interventions for arthritis. This pathway holds potential for effectively promoting health, poverty alleviation and rural revitalization strategies in China, thereby contributing to the goals of the Healthy China 2030 Plan.
Prudence was imperative in the interpretation and practical implementation of our findings. Firstly, a limitation was the relatively advanced age of the respondents, which could introduce potential biases such as recall bias, survival bias, and competing risk of mortality. We conducted a sensitivity analysis utilizing a competing-risks model, which yielded consistent results. Secondly, due to the lack of precise details regarding food consumption, we employed a basic scale. It was elucidated that utilizing intake frequency instead of intake itself was a dependable and efficacious approach. Thirdly, it should be noted that the identification of participants afflicted with arthritis relied on self-reporting rather than hospital diagnosis, thus presenting a limitation. Nonetheless, to address this concern, we conducted a comprehensive analysis incorporating multi-adjusted time-dependent Cox regression, employing various definitions of arthritis. Encouragingly, the outcome of this analysis demonstrated the effectiveness of the plant-based diet. Fourthly, due to the observational nature of this study, we controlled for some potential confounders but cannot rule out residual or unmeasured confounding. Nevertheless, during multiple model adjustments, the higher-quality plant-based diet consistently showed potential in reducing arthritis risk. Finally, our study did not differentiate between distinct subtypes of arthritis. Nonetheless, our findings support the notion that a plant-based diet holds advantages for arthritis overall, making it a viable option for dissemination in expansive rural regions.
CLHLS | Chinese Longitudinal Healthy Longevity Survey |
PDI | Plant-based diet index |
hPDI | Healthful plant-based diet index |
HRs | Hazard ratios |
CIs | Confidence intervals |
RA | Rheumatoid arthritis |
OA | Osteoarthritis |
IL-1 | Interleukin-1 |
TNF-α | Tumor necrosis factor-α |
SD | Standard deviation |
NF-κB | Nuclear factor kappa-light-chain-enhancer of activated B cells |
IL-6 | Interleukin-6 |
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d5fo02192h |
‡ These authors contributed equally to this paper and should be considered as co-first authors. |
This journal is © The Royal Society of Chemistry 2025 |