Introduction

Osteoarthritis (OA) refers to one of the most common joint diseases worldwide, influencing around 10% of men and 18% of women over 60 years of age1. it can occur in both large and small joints throughout the body, with the primary manifestations including morning stiffness, pain, deformity and joint activity limitations, and in severe cases, it can lead to joint instability or limb disability, seriously impacting quality of life2.

Epidemiological studies have indicated that physical activity is closely related to OA3. Based on different activity types, physical activity types can be categorized into occupational activity, recreational activity, and transportation activity. Additionally, different physical activity levels can be calculated based on various factors such as activity intensity, frequency, duration, and type. Most previous studies consistently revealed that physical activity can lower the risk of developing OA in joints. However, during physical activity, the load on the joints of the human body increases, and different types of physical activity impose varying loads on the knees 4,5. The study found that when examining overall physical activity, the risk of OA did not increase6. One cohort study suggested that all kinds of physical activity are protective against knee OA7,8,9. However, in a recent study from a university, it was found that low-to-moderate physical intensity of activity is protective for OA10, while high-intensity physical intensity of activity provides no protection against OA11.

While many studies indicate that an active lifestyle can reduce the risk of OA, several studies have found that leisure physical activity and occupational physical activity have different effects on OA 12,13,14. Additionally, there is no definitive evidence supporting recommendations regarding physical activity types and levels. As a result, according to the National Health and Nutrition Examination Survey (NHANES) data in the USA, this study aims to explore the relationships between physical activity types (vigorous work activity, vigorous recreational activity, moderate work activity, moderate recreational activity, walking or bicycling) and levels (non, low, medium, high) and the risk of OA. Clarifying the association between specific physical activity types and levels and OA risk will provide insights and thus guide more _targeted physical activity recommendations for OA prevention and management. The findings are conducive to informing and optimizing public health strategies to reduce the personal and socioeconomic burden of OA.

Materials and methods

Study population

National Health and Nutrition Examination Survey (NHANE), a nationally representative cross-sectional survey conducted by the Centers for Disease Control and Prevention in the United States, employs a stratified, multi-stage random sampling method. The data analyzed in this study span the years 2007–2008, 2009–2010, 2011–2012, 2013–2014, 2015–2016, 2017–2018, and the ongoing period of 2017–2020.Among the 56,171 participants in NHANES from 2007 to 2020, rigorous assessments were conducted to ensure the integrity, consistency, and logical coherence of physical activity data. In this present study, a meticulous examination of the dataset led to the exclusion of individuals with missing, refused, or unknown information regarding physical activity, sedentary time, and durations exceeding 24 h (n = 24,376). Additionally, participants lacking, refusing, or with unknown covariate information (n = 8894). Ultimately, 22,901 participants remained with no missing or confounding information regarding key outcomes, exposures, or variables (Fig. 1).

Fig. 1
figure 1

Patient Enrollment Process Flowchart.

Diagnosis of OA

OA diagnosis data was from the “Medical Conditions” questionnaire section of the NHANES. First of all, participants were asked if doctor ever said they had arthritis. If they answered “yes,” they would be further asked to identify “which type of arthritis was it” (The arthritis was classified as OA, rheumatoid arthritis, psoriatic arthritis, and others based on NHANES questionnaire data).

Assessment of types of physical activity and levels of physical activity

The Global Physical Activity Questionnaire was used to collect respondent-level weekly physical activity information15. Physical activity levels were examined for three categories of physical activity participation: vigorous work activity/vigorous recreational activities, moderate work activity/moderate recreational activities, and walk/bicycle information. Vigorous work activity was defined as typical activity that induced significant increases in breathing or heart rate for at least 10 min continuously, such as hauling or lifting heavy weights, excavating, or building work. Vigorous recreational activities were characterized as high-intensity sports, fitness activities, or leisure activities that generate significant increases in respiration or heart rate, such as jogging or basketball. Moderate work activity was defined as any activity that induces minor increases in breathing or heart rate, such as brisk walking or carrying light burdens for at least 10 min continuously. Moderate leisure activities were defined as those that generate a slight rise in breathing or heart rate for at least 10 min continuously, such as brisk walking, biking, swimming, or volleyball. Walking or use of a bicycle was defined as a transportation to school/work or for shopping for at least 10 min continuously to get to and from places.

According to the World Health Organization (WHO) Analytical Guidelines, physical activity (PA) is converted into metabolic equivalent minutes (MET-min) per week for moderate to vigorous PA16. MET values vary by the type of exercise, and NHANES provides recommended MET values for each PA. PA is calculated based on the MET values, weekly frequency, and duration of each PA using the formula: PA (MET-min/week) = MET × weekly frequency × duration17. A PA value of 0 indicates participants who are not engaged in any PA. Finally, the participants were classified into different groups based on their PA level. The groups were defined as NPA: Non-physical activity (0 MET-min/week), Low PA(LPA) (Q1), Medium PA (MPA) (Q2) and High PA (HPA) (Q3)18.

Assessment of covariates

The following information was gathered using a household structured questionnaire: age (20–44, 45–64, ≥ 65 years), gender (male and female), race (Mexican American, Other Hispanic, Non-Hispanic white, Non-Hispanic black and Other race), level of education (less than high school, high school, more than high school), Examination data included BMI (kg/m2). Laboratory data covered blood urea nitrogen (BUN, mmol/L)19, total calcium (Ca, mmol/L)20, phosphorus (P, mmol/L)21, triglycerides (TG, mmol/L)22, uric acid (UA, µmol/L)23, and total cholesterol (TC, mmol/L)24. As a final point, questionnaire data included information on smoking behavior (Yes/No)25, alcohol consumption (Yes/No)26, hypertension (Yes/No)27, and diabetes (Yes/No)28.

Statistical analysis

We used weighted samples and considered stratification and clustering in the design to generate nationally representative estimates for the US population. The NHANES cycles and weights were constructed according to the guidelines for continuous NHANES analysis29. The participants’ baseline characteristics were classified based on their physical activity level and OA status. A descriptive statistical analysis was performed on basic demographic information (age, gender, race, level of education), physical activity types (vigorous work activity, vigorous recreational activities, moderate work activity, moderate recreational activities, walk or bicycle, physical activity levels), body measurements (BMI), and diseases (hypertension, diabetes).Laboratory data included blood urea nitrogen, total calcium, phosphorus, triglycerides, uric acid, and total cholesterol. The categorical variables were reported in frequency (%), with the chi-square test for unordered categorical variables and the chi-square trend test for ordered multi-categorical variables.

STATA version 17.0 was used for statistical analysis (Stata Corp LP, College Station, TX, USA). R version 4.2.3 was used to create the forest graphs. A p-value under 0.05 was considered statistically significant.

Result

Among the 22,901 participants (41.5% male, 58.4% female), the overall prevalence of OA was 11.8%. Table 1 shows the baseline characteristics of the participants stratified by OA status. In general, OA was more common in middle-aged (45–64 years) and elderly (≥ 65 years) individuals, females and non-Hispanic whites. Higher education levels of above high school, BMI ≥ 25 kg/m2, heavy smoking, hypertension diagnosis and diabetes diagnosis were related to higher risk of OA. Compared to vigorous work activity and moderate work activity, OA prevalence was lower in those engaging in vigorous recreational activity and moderate recreational activity (Table 1). Females and participants with higher BMI were more likely to perform non-physical activity, and had higher risks of OA, hypertension or diabetes. Participants who received education higher than high school tended to have low physical activity level. Those with lower BMI and non-whites were more likely to have medium physical activity level. Younger individuals (20–44 years) and males were more inclined to have high physical activity level, and experienced lower risks of OA, hypertension, and diabetes (Table 2).

Table 1 Baseline characteristics of the research population with and without OA.
Table 2 Baseline Characteristics of Participants by Physical Activity Levels Status: NHANES,2007–2020.

Figure 2 shows the relationship between different types of physical activity and risk of OA. The three models were defined as follows: Model 1 had no covariates adjusted; Model 2 adjusted for gender and age; Model 3 further adjusted for a comprehensive set of variables including age, gender, race, level of education, BMI, smoking status, alcohol consumption, hypertension, diabetes, blood urea nitrogen, total calcium, phosphorus, triglycerides, uric acid, and total cholesterol. We found that all types of physical activity in this study were associated with OA. In Model 1, moderate work activity was a risk factor (OR: 1.20, 95% CI: 1.20–1.40), while vigorous work activity, vigorous recreational activity and walking or bicycling were protective factors. Additionally, in particular, there was the strongest statistically significant inverse relationship between vigorous recreational activity and OA risk (OR: 0.31, 95% CI: 0.21–0.48). In Model 3 after controlling for age, gender, race, education level, BMI, smoking status, alcohol consumption, hypertension, diabetes, as well as laboratory data (blood urea nitrogen, total calcium, phosphorus, triglycerides, uric acid, total cholesterol) (in the fully adjusted model) high-intensity recreational activity (OR: 0.61, 95% CI: 0.41–0.92) and moderate-intensity recreational activity (OR: 0.79, 95% CI: 0.68–0.92) were both protective factors against OA, reducing the risk by 39% and 21%, respectively. High-intensity recreational activity showed the strongest protective effect against OA.

Fig. 2
figure 2

The association between physical activity types and OA. Model 1: no covariates were adjusted; Model 2 adjusted for gender, age; Model 3 adjusted for age, gender, race, level of education, BMI, smoking status, alcohol consumption, hypertension, diabetes, blood urea nitrogen, total calcium, phosphorus, triglycerides, uric acid, total cholesterol.

Figure 3 illustrates that in Model 1, as the physical activity level increases, the risk of OA gradually decreases. Compared to the population who does not conduct physical activity, the risk of OA decreased by 25% (OR: 0.75, 95% CI: 0.65–0.87), 36% (OR: 0.64 95% CI: 0.54–0.76) 和32% (OR: 0.68 95% CI: 0.58–0.81) for low, medium, and high physical activity level, respectively. In Model 3, among all levels of physical activity, individuals with medium physical activity level showed the strongest protective factor against OA (OR: 0.72, 95% CI: 0.60–0.87). The dose–response relationship between total physical activity level and the risk of OA varied with age and gender. Regardless of the physical activity level, no significant correlation was found with the risk of OA in age groups of ≥ 65 years (Additional file S1C). Medium physical activity, however, demonstrated preventive benefits against OA in the age groups of 20–44 years and 45–64 years and (OR:0.46, 95% CI:0.26–0.81, OR:0.73, 95% CI:0.54–0.99), reducing the risk of OA (Additional files S1A and B). In gender-stratified analyses, medium physical activity level was the strongest protective factor against OA in males (OR:0.70,95% CI:0.52–0.95), while no significant correlation was found with the risk of OA in females (Additional files S2A and B).

Fig. 3
figure 3

The association between physical activity levels and OA. Model 1: no covariates were adjusted; Model 2 adjusted for age, gender; Model 3 adjusted for age, gender, race, level of education, BMI, smoking status, alcohol consumption, hypertension, diabetes, blood urea nitrogen, total calcium, phosphorus, triglycerides, uric acid, total cholesterol.

As displayed in Figs. 4A–C, across five different types of physical activity, namely vigorous work activity, vigorous recreational activity, moderate work activity, moderate recreational activity, and walking or bicycling, the impact of three physical activity levels on the risk of OA varies compared to non-physical activity. Individuals engaged in vigorous work activity with different levels of physical activity shows no significant relationship to OA risk (Fig. 4A). In addition, for vigorous recreational activity type, participating in vigorous recreational activities at a medium physical activity level reduces the incidence of OA by 59% (OR: 0.41, 95% CI: 0.17 -0.95) (Fig. 4A). In the moderate work activity type, low physical activity level is positively related to OA, raising the risk by 29%, respectively (OR:1.29, 95% CI: 1.01 -1.66). In moderate recreational activity type, a high physical activity level lowers the risk of OA by 29% (OR: 0.71, 95% CI: 0.55–0.91) (Fig. 4B). Regarding walking or bicycling related to transportation, there is no significant association between different levels of physical activity and the risk of OA (Fig. 4C).

Fig. 4
figure 4

The association between physical activity levels and OA. Model 1: no covariates were adjusted; Model 2 adjusted for gender, age; Model 3 adjusted for age, gender, race, level of education, BMI, smoking status, alcohol consumption, hypertension, diabetes, blood urea nitrogen, total calcium, phosphorus, triglycerides, uric acid, total cholesterol.

Discussion

In this nationally representative study, vigorous recreational activities and moderate recreational activities serve as protective factors against OA, with vigorous recreational activities showing the greatest protective effect. Subdividing weekly metabolic equivalents revealed that engaging in weekly exercise at a medium physical activity level reduces the risk of OA, although this does not apply to all kinds of physical activities, including work-related activities and vigorous recreational activity. After categorizing different physical activity types, it was found that a medium physical activity level of vigorous recreational activity can reduce the risk of OA, while a high physical activity level of moderate recreational activity can lower the risk of OA. However, a low physical activity level of moderate work activity increases the risk of OA. Furthermore, the correlation between physical activity levels and the risk of OA varies with gender and age.

The relation between physical activity and OA has long been a focus of discussion. Physical activity, defined as movement of the body and limbs causing energy expenditure through muscle contraction, is a primary therapeutic approach against the health impacts of a sedentary lifestyle. It is also a frontline treatment for OA, with a benefit/risk ratio far exceeding that of drug therapy, supported by high levels of evidence and expert consensus, especially for lower limb OA, as confirmed by numerous studies and meta-analyses30. However, in a systematic review in 2016, 18 studies, mainly small to moderate-sized quasi-experimental (non-randomized) clinical ones were identified, providing conflicting evidence on the preventive role of physical activity in knee OA following anterior cruciate ligament injury31. Another study demonstrated that moderate physical activity may increase the sGAG content in cartilage and potentially prevent OA. However, when sGAG content in the cartilage is depleted, exercise may induce OA32. Therefore, the type, frequency, duration, and intensity of physical activity can have different impacts on OA.

Regarding the impact of work-related activities on OA, this study found that a low physical activity level of moderate work activity increases the risk of OA, mainly due to prolonged standing, sitting, and repetitive movements in occupational activities. From a biomechanical perspective, standing places continuous compressive loads on knee cartilage, which induces wear-promoting lubrication modes33. Long-term engagement in occupations that involve squatting, kneeling, or climbing stairs can lead to chronic wear of joint cartilage, eventually causing OA34,35. A comparative study found that the duration of squatting is positively associated with tibiofemoral OA36. For individuals who sit for extended periods, prolonged sitting can lead to increased pressure within the knee joint, muscle atrophy, and joint stiffness37. As muscle strength declines, joint stability is reduced, accelerating cartilage degeneration38,39. On the other hand, when joints are deprived of regular loading, cartilage degradation may also occur.

In recreational activities, high physical activity levels of vigorous recreational activity and medium physical activity levels of moderate recreational activity both reduce the risk of OA. Various activities, including exercise and fitness have been validated for their protective effects against OA40,41,42. Some activities, including Tai Chi43, walking44, water sports45, and cycling42, stand out in both the prevention and treatment of OA. However, according to surveys, certain traumatic physical activities, including American football, soccer, weightlifting, and ice hockey, may increase the occurrence of OA46,47,48. Nevertheless, it is the rising rate of knee joint injuries in these groups rather than participation in physical activity itself that elevates the risk of OA.

Among different physical activity levels, a medium physical activity level has a protective effect against OA. In a study, it was found that individuals with the highest level of physical activity were characterized by an elevated risk of knee OA outcomes49. Yang et al.50demonstrated through histological and protein analysis that moderate-intensity treadmill exercise can improve OA by downregulating TRAIL. They also found that adaptive intensity exercise can lower the sensitivity of cartilage cells to inflammation, while excessive exercise can lead to progressive damage, hinder matrix synthesis, and stimulate the production of matrix-degrading enzymes, thereby increasing the risk of OA51,52. Research has revealed that normal joint cartilage and chondrocytes can benefit from mechanical stimuli, but when the mechanical stimuli are excessive or prolonged and cannot be tolerated, it can cause damage42. Regular participation in physical activity has been demonstrated to be dramatically conducive to the treatment of knee OA. In addition, if patients cannot maintain activity and become sedentary, it may accelerate joint mechanical damage, potentially resulting in further deterioration of cartilage and aggravating OA53. Therefore, maintaining a medium physical activity level and achieving a balance in physical activity is recommended. After conducting subgroup analyses of medium physical activity levels, this study found that the protective effect of medium physical activity levels diminished in individuals aged 65 and above. This could be due, on one hand, to the decline in joint cartilage repair capacity associated with aging54. On the other hand, the weakening of muscles that accompanies aging makes joints more prone to damage, meaning that medium physical activity levels might no longer offer effective protection55. Additionally, the increased likelihood of comorbidities with advancing age may further impact the effectiveness of physical activity56.

The subgroup analysis in this study suggests that low and medium physical activity levels are protective against osteoarthritis risk in men, while no significant association was observed in women. This may be related to sex differences in joint cartilage57. Jones et al.58 found that men tend to have more knee cartilage compared to women. Faber et al. 59noted that, compared to men, women have 19.9% less patellar cartilage and 46.6% less medial tibial cartilage. Additionally, Pachowsky et al.60 discovered that women have poorer cartilage quality and lower healing capacity compared to men. Moreover, these differences could be related to metabolic variations influenced by hormone levels61.

Our study collected data from a large-scale, nationally representative sample of U.S. adults with long-term longitudinal follow-up and adjusted for most potential confounding factors. We also used appropriate NHANES sample weights to analyze the data, enhancing the reliability and generalizability of our results. The results indicated that vigorous recreational activity and moderate recreational activity are protective factors against OA. Specifically, at a medium physical activity level, vigorous recreational activity can reduce the risk of OA, while at a high physical activity level, moderate recreational activity can also reduce the risk. However, moderate work activity at a low physical activity level is associated with an increased risk of OA. Based on these findings, we can develop _targeted exercise recommendations for individuals at higher risk of OA, guide individual physical activity interventions, and improve public health strategies to reduce the incidence of OA. Finally, multiple subgroup analyses were performed to confirm the robustness of the results.

Despite these strengths, our study has certain limitations. Firstly, it is a cross-sectional study, and NHANES data are typically cross-sectional, providing health and nutrition information at a specific point in time without revealing trends or causal relationships. Secondly, This study relied on self-reported data to determine the presence of OA, which may introduce reporting bias and result in some individuals with actual OA being incorrectly classified as non-OA. This could lead to an underestimation of the true prevalence of OA and potentially impact the statistical analysis results. However, despite the possibility of misclassification of OA cases, the large sample size used in this study enhances the generalizability of the findings. Additionally, the use of multiple models and subgroup analyses to adjust for confounding factors increases the reliability of the results. Therefore, the findings of this study still maintain a certain degree of validity and robustness. Thirdly, despite adjusting for various covariates, some potential confounding factors might still have been overlooked.

Conclusion

To conclude, this study demonstrates the cross-sectional impact of physical activity on OA. Among all types of physical activity, vigorous recreational activity and moderate recreational activity are found to have a preventive effect on OA. In terms of physical activity levels, low physical activity levels of moderate work activity are associated with an increased risk of OA, while moderate physical activity levels are confirmed to have a protective effect against OA in the age groups of 20–44 and 45–64. However, gender-stratified analyses reveal that both low and moderate physical activity levels provide protection against OA in males, with moderate physical activity levels showing a more significant protective effect.