ORIGINAL RESEARCH article

Front. Public Health, 28 April 2025

Sec. Aging and Public Health

Volume 13 - 2025 | https://doi.org/10.3389/fpubh.2025.1518572

Global shifts in osteoarthritis subtype trends among older adults due to elevated BMI: an age-period-cohort analysis based on the global burden of disease database

Zhiqiang Li
Zhiqiang Li*Yaxiang ChenYaxiang ChenZhao Shen
Zhao Shen*
  • Department of Orthopedics, Cangzhou People’s Hospital, Cangzhou, Hebei, China

Purpose: It remains unclear whether elevated Body Mass Index(BMI)has a similar impact on the disease burden of osteoarthritis subtypes in older adults. This study aims to compare the long-term trends of osteoarthritis subtypes caused by high BMI across different gender groups globally from 1990 to 2021.

Methods: We obtained cross-sectional data from the Global Burden of Disease, Injuries, and Risk Factors Study (GBD) 2021 (https://vizhub.healthdata.org/gbd-results/). The disease burden of osteoarthritis subtypes in older adults attributable to high BMI was quantified using Years Lived with Disability (YLDs). Linear regression and the Age-Period-Cohort (APC) method were employed to calculate the trends in Age-standardized Years lived with disability rate (ASYR), adjusting for age, period, and cohort effects.

Results: The ASYR of osteoarthritis attributable to high BMI in older adults globally has shown a continuous upward trend over the past 32 years, with an Estimated Average Percentage Change (EAPC) of 0.96 (95% CI: 0.94 to 0.99). Specifically, the EAPC for hip osteoarthritis was 0.73 (95% CI: 0.70 to 0.76), while for knee osteoarthritis, it was 0.99 (95% CI: 0.96 to 1.02). China recorded the highest number of osteoarthritis YLDs globally, reaching 0.59 million (95% UI: −0.05 to 1.71). The United States had one of the highest ASYR rates for osteoarthritis at 410.85 per 100,000 (95% UI: −44.47 to 1,083.52), while India exhibited the highest EAPC for osteoarthritis worldwide at 2.74 (95% CI: 2.70 to 2.79), with hip osteoarthritis at 3.36 (95% CI: 3.25 to 3.48) and knee osteoarthritis at 2.70 (95% CI: 2.65 to 2.75). The local drift curves indicated a slow upward trend in the annual percentage change of YLDs for both hip and knee osteoarthritis attributable to high BMI across all age groups. In terms of gender, males exhibited a higher rate and risk of YLDs associated with high BMI.

Conclusion: Our findings provide strong evidence that the ASYR associated with high BMI globally have continuously increased over the past 32 years, with consistent patterns of change observed across different osteoarthritis subtypes. This highlights the critical role of BMI control in effectively alleviating the burden of osteoarthritis in older adults.

Introduction

Osteoarthritis is one of the most prevalent chronic musculoskeletal diseases among older adults globally, leading to long-term pain, functional impairment, and a decline in quality of life (13). With the global rise in ageing and obesity rates, the burden of osteoarthritis is expected to continue to increase. High BMI is widely recognised as a major risk factor for osteoarthritis, particularly in weight-bearing joints such as the knee and hip (4). According to findings from the GBD 2021 study, obesity is one of the top five modifiable health risk factors worldwide, and the contribution of high BMI to the osteoarthritis burden has become increasingly evident over the past few decades (5).

Although the association between high BMI and osteoarthritis is well-established, previous studies, including those from the GBD project, have yet to provide a comprehensive estimate of the long-term impact of high BMI on osteoarthritis and its subtypes (6, 7). Most existing research has focused solely on the overall burden of osteoarthritis, without thoroughly exploring the differences between osteoarthritis subtypes in different anatomical sites, such as knee osteoarthritis and hip osteoarthritis, nor considering the impact of high BMI across different age groups (810). Therefore, assessing the burden trends of various osteoarthritis subtypes caused by high BMI is crucial for developing targeted prevention strategies.

In this study, we utilised data obtained from the GBD 2021 database to describe the YLDs due to osteoarthritis and its subtypes attributable to high BMI. We calculated the EAPC to reveal the global trends in the burden of osteoarthritis subtypes in older adults caused by high BMI from 1990 to 2021. Additionally, we employed an APC model to analyse the effects of age, period, and cohort on the prevalence of osteoarthritis subtypes.

Methods

Study population and data collection

The data for this study were sourced from the GBD 2021 dataset, which provides comprehensive information on the global and regional burden of 371 diseases and injuries, along with 88 risk factors, across 204 countries and territories from 1990 to 2021. The data was accessed through the Global Health Data Exchange (GHDx) platform,1 where detailed information on the disease burden indicators, including YLDs for osteoarthritis in older adults aged 60 and above, was downloaded. The dataset also provides data on various demographic groups, osteoarthritis subtypes (hip and knee osteoarthritis), and BMI-related factors. According to the World Health Organization (WHO) standards, high BMI is defined as a BMI of 25 kg/m2 or higher, which includes overweight (BMI 25–29.9 kg/m2) and obesity (BMI ≥ 30 kg/m2) (11). In contrast to the WHO definition, which applies to all age groups, the GBD group uses different definitions for children and adults based on specific age cutoffs. In this study, we applied the GBD-specific definition of high BMI, as outlined in the GBD methods (source: https://www.healthdata.org/sites/default/files/methods_appendices/2021/Cogen_metab_bmi_writeup_gbd2020_updated01-31%E2%80%932024.pdf). Data on socio-demographic index (SDI) were also collected to assess the impact of socioeconomic factors, which could provide a context for understanding the influence of regional disparities on the disease burden. As this study utilised a publicly available database containing de-identified data, ethical approval was not required. This study adhered to the Guidelines for Accurate and Transparent Health Estimates Reporting (GATHER) (12).

Statistical analysis

To examine the global distribution, regional disparities, and population-level status of the osteoarthritis burden and its subtypes among older adults, we quantified the burden attributable to high BMI using the ASYR obtained from the GBD 2021 dataset. Global mapping and regional comparative analyses were conducted to visualise the distribution of osteoarthritis burden. Data were aggregated according to the geographical regions defined by the GBD study, with maps created using R (version 4.3.2) and the ‘ggplot2’ and ‘sf’ packages. Linear regression was chosen as a statistical method due to its effectiveness in modeling continuous data trends over time, which allows for the evaluation of changes in osteoarthritis burden in relation to BMI over the 32-year period. The APC model was specifically employed because it accounts for the impact of age, period, and cohort effects independently, making it suitable for analyzing long-term trends influenced by various temporal factors, such as population aging, medical advancements, and societal changes [38913425]. The EAPC was employed to assess the time trends of ASYR from 1990 to 2021, with EAPC values calculated through linear regression models and outputs processed using the ‘broom’ package. This method is widely used for studying disease burden trends and has been validated in previous studies for its capacity to quantify annual percentage changes in disease rates [39,435,408, 38,745,964].

To eliminate the impact of differences in population age structure across regions or time periods on disease burden, this study uses the direct standardization method to calculate the Age-Standardized Rate (ASR). The specific age group rate data is selected, and the age categories from the GBD database are matched with the global standard population data to obtain the standard population weights for each age group. The weighted average rate is then calculated using the following formula: (Val is the raw rate value for the specific age group, Population_weight is the standard population weight for that age group).

A S R = val × Population _ weight Population _ weight

The APC model was employed to assess the impact of age, period, and cohort on YLD rates for osteoarthritis subtypes attributable to high BMI. APC analysis minimises the interaction effects among these three factors, providing more precise estimates for each. The age effect reflects the influence of population ageing, while the period effect captures changes in YLD rates over time that affect all age groups, often due to advancements in disease screening, medical technologies, or reclassification of diseases. The cohort effect highlights long-term trends in disease incidence and mortality, shaped by cohort-specific lifestyles, environmental changes, and exposure to risk factors (13, 14).

In this study, within the APC framework, we evaluated the following parameters: net drift, which represents the overall annual percentage change in the attributable YLD rate; local drift, indicating the annual percentage change in YLD rates for each age group; the longitudinal age curve, which reflects the age effect by showing age-specific rates after adjusting for drift; period rate ratios (PRR), which represent the period effect and the relative risk of YLDs across different time periods; and cohort rate ratios (CRR), which reflect the cohort effect and the relative risk of YLDs across different birth cohorts. The APC analysis considered age groups starting from 60 years, segmented into 5-year intervals (e.g., 60–64 years, 65–69 years, etc.). Trends were assessed for statistical significance using 95% confidence intervals (CI), and descriptive statistics for all key variables were presented as mean values with 95% uncertainty intervals (UI). All analyses were conducted using R version 4.3.2.

Results

Global and SDI regional trends in YLDs

In 1990, the YLDs attributable to high BMI for osteoarthritis in older adults globally were 84.02 ten thousands (95% UI: −7.23 to 240.44), corresponding to an ASYR of 175.06 per 100,000 (95% UI: −14.94 to 500.89). By 2021, the global YLDs had significantly increased to 2.55 million (95% UI: −0.24 to 7.18), with the ASYR rising to 234.51 per 100,000 (95% UI: −21.52 to 659.64). Although the 95% UI for YLDs and ASYR include negative values, it is important to note that the presence of these negative intervals reflects the inherent uncertainty in the data and does not invalidate the overall trends. These intervals are a statistical artifact of the model’s estimation process, particularly when estimating changes over long periods across diverse global regions. The EAPC was 0.96 (95% CI: 0.94 to 0.99). While the 95% CI for the EAPC is narrow, indicating statistical significance, it is essential to acknowledge that the UI ranges for YLDs and ASYR suggest that the estimates are subject to uncertainty, which should be considered when interpreting the results. The YLDs for hip osteoarthritis increased from 9.04 ten thousands (95% UI: −0.81 to 24.50) in 1990 to 25.82 (95% UI: −2.44 to 68.07) in 2021. The global ASYR for hip osteoarthritis rose from 19.37 per 100,000 (95% UI: −1.71 to 52.62) in 1990 to 23.93 per 100,000 (95% UI: −2.25 to 63.13) in 2021, with an EAPC of 0.73 (95% CI: 0.70 to 0.76). The YLDs for knee osteoarthritis increased significantly from 0.75 million (95% UI: −0.64 to 2.15) in 1990 to 2.30 million (95% UI: −0.21 to 6.46) in 2021. In 1990, the global ASYR for knee osteoarthritis was 155.69 per 100,000 (95% UI: −13.22 to 446.14), rising to 210.57 per 100,000 (95% UI: −19.24 to 593.06) by 2021. The EAPC for ASYR of knee osteoarthritis was 0.99 (95% CI: 0.96 to 1.02), showing a growth trend closely aligned with overall osteoarthritis, slightly higher than the growth rate for hip osteoarthritis (Figures 1A-C, 2A, Table 1).

Figure 1
www.frontiersin.org

Figure 1. The EAPC of the ASYR for osteoarthritis (A), hip osteoarthritis (B), and knee osteoarthritis (C) among individuals aged 60 years and above, stratified by GBD regions and SDI quintiles, from 1990 to 2021. EAPC, Estimated Average Percentage Changes; GBD, Global Burden of Disease; SDI, Socio-demographic Index; ASYR, Age-standardized Years lived with disability rate.

Figure 2
www.frontiersin.org

Figure 2. The trends in osteoarthritis, hip osteoarthritis, and knee osteoarthritis among individuals aged 60 years and above from 1990 to 2021 (A), as well as the trends for osteoarthritis (B), hip osteoarthritis (C), and knee osteoarthritis (D) classified by the SDI quintiles. ASYR, Age-standardized Years lived with disability rate; YLDs, Years Lived with Disability; SDI, Socio-demographic Index.

Table 1
www.frontiersin.org

Table 1. The global YLDs and ASYR, as well as those of the SDI regions and 21 specific regions in 1990 and 2021, along with the EAPC of ASYR from 1990 to 2021.

Among the five distinct SDI regions, older adults in high-SDI areas consistently exhibited the highest ASYR for osteoarthritis and its two subtypes. However, the EAPC values in high-SDI regions were markedly lower than those in other SDI regions for both overall osteoarthritis and its subtypes. Over the past 32 years, the ASYR in all five SDI regions has demonstrated a consistently increasing trend (Figures 1AC, 2B–D, Table 1).

Regional trends in YLDs

In 2021, among the 21 regions, High-income North America recorded the highest ASYR of osteoarthritis, at 390.27 per 100,000 (95% UI: −41.82 to 1031.22). while Southeast Asia had the lowest ASYR, at 111.66 per 100,000 (95% UI: −8.43 to 330.80). Over the past 30 years, South Asia had the highest EAPC at 2.64 (95% CI: 2.59 to 2.69), ranking first among the 21 regions, while High-income North America had the lowest EAPC at 0.38 (95% CI: 0.21 to 0.54). Regarding osteoarthritis subtypes, in 2021, High-income North America had the highest ASYR for hip osteoarthritis, reaching 64.46 per 100,000 (95% CI: −6.93 to 166.74), while South Asia had the lowest ASYR for hip osteoarthritis at 8.15 per 100,000 (95% CI: −0.57 to 22.94). South Asia also had the highest EAPC for hip osteoarthritis at 2.64 (95% CI: 2.59 to 2.69), whereas Oceania had the lowest EAPC at 0.71 (95% CI: 0.63 to 0.78). For knee osteoarthritis, Australia recorded the highest ASYR in 2021, at 329.85 per 100,000 (95% UI: −32.74 to 890.47), while Southeast Asia had the lowest ASYR for knee osteoarthritis, at 102.63 per 100,000 (95% UI: −7.83 to 302.24). Once again, South Asia showed the highest EAPC for knee osteoarthritis at 2.61 (95% CI: 2.56 to 2.66), while High-income North America had the lowest EAPC at 0.25 (95% CI: 0.05 to 0.44) (Figures 1AC, 3A–C, Table 1).

Figure 3
www.frontiersin.org

Figure 3. The ASYR of osteoarthritis (A), hip osteoarthritis (B), and knee osteoarthritis (C) among individuals aged 60 years and above across 21 regions in 1990 and 2021. ASYR, Age-standardized Years lived with disability rate.

National trends in YLDs

Among 204 countries, China had the highest number of YLDs for osteoarthritis in older adults in 2021, with a total of 0.59 million (95% UI: −0.05 to 1.71). The United States recorded the highest YLDs for hip osteoarthritis at 5.32 ten thousand (95% UI: −0.58 to 13.71). For knee osteoarthritis, China led globally with 0.56 million YLDs (95% UI: −0.05 to 1.63). In terms of ASYR for overall osteoarthritis and its two subtypes, the United States had the highest rates globally, with an ASYR of 410.85 per 100,000 (95% UI: −44.47 to 1,083.52) for osteoarthritis, 67.23 per 100,000 (95% UI: −7.31 to 173.39) for hip osteoarthritis, and 343.62 per 100,000 (95% UI: −36.94 to 917.26) for knee osteoarthritis. On the other hand, the Democratic Republic of Timor-Leste had the lowest ASYR, with ASYRs for osteoarthritis, hip osteoarthritis, and knee osteoarthritis being 13.66 (−0.90 to 42.16), 0.90 (−0.05 to 2.58), and 12.76 (−0.85 to 39.37), respectively. Regarding EAPC, India showed the highest global rates for osteoarthritis and its two subtypes, with an EAPC of 2.74 (95% CI: 2.70 to 2.79) for osteoarthritis, 3.36 (95% CI: 3.25 to 3.48) for hip osteoarthritis, and 2.70 (95% CI: 2.65 to 2.75) for knee osteoarthritis. On the other hand, Georgia had the lowest EAPC for osteoarthritis, at 0.38 (95% CI: 0.38 to 0.39), Kingdom of Denmark had the lowest EAPC for hip osteoarthritis, at 0.21 (95% CI: −0.00 to 0.42), and the People’s Republic of Bangladesh had the lowest EAPC for knee osteoarthritis, at 3.07 (95% CI: 2.96 to 3.19). (Figures 4AC, Supplementary Figures 1A–C, 2A–C, Supplementary Table 1).

Figure 4
www.frontiersin.org

Figure 4. The EAPC of the ASYR for osteoarthritis (A), hip osteoarthritis (B), and knee osteoarthritis (C) among individuals aged 60 years and above across 204 countries and territories from 1990 to 2021. ASYR, Age-standardized Years lived with disability rate; EAPC, Estimated Average Percentage Changes.

Net drift and local drift across different age groups

The net drift for both males and females was positive, indicating an upward trend in the global ASYR for osteoarthritis among individuals aged 60 and above from 1990 to 2021. This suggests a progressively increasing burden of osteoarthritis in the older adult population worldwide. The local drift curves for both genders also showed an upward trend, signifying that the proportion of the osteoarthritis burden increased with age (Figure 5A).

Figure 5
www.frontiersin.org

Figure 5. The local drift of YLD rates for osteoarthritis (A), hip osteoarthritis (B), and knee osteoarthritis (C) across eight age groups, stratified by sex, from 1992 to 2021. The dots and shaded areas denote the local drift (ie, annual percentage change of age-specific incidence, % per year) and their corresponding 95% CIs. YLD, Year Lived with Disability; CIs, Confidence Intervals.

From 1990 to 2021, the net drift for hip osteoarthritis in both males and females remained above zero, indicating a continuous rise in ASYR globally. Notably, the net drift value was higher in males, reflecting a faster rate of increase in the incidence of hip osteoarthritis among men. The local drift curve for males exhibited a U-shape, with significant increases in the 60–69 age group and after 85 years, suggesting an elevated burden of hip osteoarthritis in these age brackets. In contrast, the local drift for females showed a declining curve, indicating a slower progression of hip osteoarthritis and a reduced burden among women aged 60 and above (Figure 5B).

The global net drift for knee osteoarthritis was positive, meaning the ASYR for knee osteoarthritis in individuals aged 60 and above continued to rise. Similar to hip osteoarthritis, the rate of increase was faster in males. The local drift curves for both males and females demonstrated that the burden of knee osteoarthritis experienced a brief decline between the ages of 60 and 84, followed by stabilisation (Figure 5C).

Impact of APC on osteoarthritis and its subtypes in older adults

In terms of the CRR for osteoarthritis and its two subtypes, the CRR values for both males and females were close to 1, indicating that there was no significant difference in the incidence of osteoarthritis and its subtypes across different generations from 1895 to 1959. However, the global overall CRR, as well as that for males and females, showed an upward trend, reflecting an increasing burden of osteoarthritis in older adults, with men consistently having a higher risk than women (Figures 6AC).

Figure 6
www.frontiersin.org

Figure 6. The APC model demonstrates the impact of age, period, and birth cohort on the YLD rates of osteoarthritis, hip arthritis, and knee arthritis among the older adult. The cohort effect is evidenced by the cohort-specific relative risk, with each cohort being compared to the reference cohort (1925–1929), while adjusting for the influence of age and period (A–C). The PRR of each period group is compared with the reference period (2002–2006), with adjustments made for age and cohort effects (D–F). The age effect is described through the longitudinal rate of a specific age, adjusted for variations across different birth cohorts and accounting for period-specific biases (G–I). The dots and shaded areas denote the incidence rates or rate ratios and their corresponding 95% CIs. YLD, Year Lived with Disability; CIs, Confidence Intervals; PRR, Period Rate Ratios.

The PRR demonstrated a continuous increase in the relative risk of osteoarthritis and its two subtypes for the global population, males, and females since 1990, indicating a rising disease burden (Figures 6DF).

The YLD rates for hip osteoarthritis in older adults showed a consistent upward trend, with a rapid increase in recent years. For osteoarthritis and knee osteoarthritis, the YLD rates for both genders and overall population exhibited an increasing trend in the 60–80 age group. In osteoarthritis, the YLD rates for females and the overall population remained stable after the age of 80, while males showed a gradual increase. For knee osteoarthritis, the YLD rates for females and the overall population showed a slight decline after the age of 80, in contrast to the previous trend, while males experienced a more stable and moderate increase (Figures 6GI).

Discussion

Osteoarthritis poses a significant global public health challenge, presenting particular risks to high-risk populations, including older adults, who may experience considerable joint pain, functional limitations, and a decline in quality of life. Osteoarthritis is a chronic degenerative disease that often leads to severe joint damage, affecting patients’ daily activities and increasing the risk of disability and healthcare burden among the older adult population. With the ageing population, the incidence of osteoarthritis continues to rise, particularly among individuals aged 60 and over, resulting in substantial economic and social burdens (7, 10, 15, 16). This vulnerability may also lead to a decline in self-care abilities, an increased risk of falls, and further health deterioration, making it a critical issue in global public health. Therefore, a comprehensive understanding of the epidemiological trends of osteoarthritis and its associated risk factors in older adults is essential for developing effective intervention strategies and achieving the World Health Organization’s 2030 Sustainable Development Goals. Accurate assessment of osteoarthritis trends within this population is crucial for disease evaluation and health target achievement. This study outlines the trends of osteoarthritis and its two subtypes using YLDs and ASYR, and explores their temporal dynamics using the APC model.

Overall trends

Our findings indicate a clear upward trend in ASYR, period, and cohort for the two subtypes of osteoarthritis in older adults attributable to high BMI globally, from 1990 to 2021, with the greatest increase observed in knee osteoarthritis. Notably, the ASYR for hip osteoarthritis in older adults showed a sharp increase with age, whereas the ASYR for knee osteoarthritis exhibited a significant rise between the ages of 60 and 80, with no substantial increase observed in the population over 80 years. Specifically, the YLDs for knee osteoarthritis experienced the largest growth, increasing from 749,802.64 in 1990 to 2,295,374 in 2021, demonstrating that the knee joint is more susceptible to the negative effects of high BMI. As the largest weight-bearing joint in the body, the knee endures significant pressure, and high BMI increases the joint’s burden, leading to cartilage wear and a rising incidence of osteoarthritis. In contrast, while hip osteoarthritis also exhibited an upward trend, its increase was slightly lower than that of knee osteoarthritis, possibly due to the anatomical structure and functional characteristics of the hip joint. Previous guidelines and studies have shown that reducing BMI can effectively prevent the onset and progression of both knee and hip osteoarthritis (11, 17, 18). Research by Li JS et al. demonstrated that weight loss has been proven to significantly improve knee pain and functional impairment (19). A large longitudinal study by Gabby B. Joseph et al. also suggested that weight loss may prevent both radiographic and symptomatic knee osteoarthritis (20). Since 1980, the global average BMI in adults has continued to rise, and the problems of overweight and obesity have become increasingly severe (21).

Differences across SDI regions

Among the five different SDI regions, older adults in high-SDI areas consistently had the highest ASYR for osteoarthritis and its two subtypes. This may be attributed to a higher degree of population ageing, better healthcare leading to higher diagnostic rates, and lifestyle factors in high-SDI regions. However, the EAPC in high-SDI regions was significantly lower than in other SDI regions, suggesting some success in controlling the rising burden of osteoarthritis. Possible reasons for this include more effective public health interventions, the promotion of healthier lifestyles, advances in early diagnosis, and improved treatment techniques. In contrast, the burden of osteoarthritis has been growing faster in low- and middle-SDI regions, with India having the highest global EAPC for osteoarthritis and its two subtypes. This rapid increase could be related to lifestyle changes brought on by fast economic development (e.g., the westernisation of diets and reduced physical activity), limited healthcare resources, and the accelerated pace of ageing in these countries. These findings indicate the need for enhanced osteoarthritis prevention and management in low-SDI countries, as well as the development of targeted public health strategies to address the rapidly growing disease burden.

Regions-level trend analysis

In 2021, notable regional disparities in the burden of osteoarthritis were identified, with variations closely linked to socio-demographic factors, particularly the SDI and obesity trends. High-income regions such as North America and Western Europe, with higher SDI and rising obesity rates, exhibited a higher osteoarthritis burden, as reflected in the highest ASYR of osteoarthritis in High-income North America, reaching 390.27 per 100,000 (95%UI: −41.82 to 1031.22). These regions have experienced significant increases in obesity rates over the past three decades, which are strongly correlated with a rise in osteoarthritis cases, particularly in weight-bearing joints such as the knees and hips. In contrast, regions with lower SDI, such as South Asia and parts of Sub-Saharan Africa, have shown slower increases in obesity prevalence. Despite these regions having lower obesity rates, South Asia has demonstrated the highest EAPC of osteoarthritis at 2.64 (95%CI: 2.59 to 2.69) over the past 30 years, signaling a rapid increase in osteoarthritis cases. This is likely due to a combination of factors such as rapid urbanization, changes in lifestyle, and an aging population. Although obesity prevalence is lower in South Asia compared to high-income regions, other factors such as a growing urban population and improved healthcare access may contribute to the increase in osteoarthritis burden. The positive correlation between BMI and osteoarthritis remains evident, and as South Asia urbanizes and adopts more Westernized lifestyles, obesity rates are expected to rise, leading to further increases in osteoarthritis prevalence. The differences between regions highlight the complex interplay between socio-economic development, urbanization, and obesity trends, suggesting that a multifactorial approach is required to understand the burden of osteoarthritis. Regions with high SDI are experiencing a significant burden due to obesity and lifestyle factors, while regions with lower SDI, though starting with lower obesity rates, may see increases in osteoarthritis cases due to urbanization and demographic shifts.

Country-level trend analysis

At the country level, China and India represent contrasting examples of osteoarthritis burden and trends, influenced by their SDI, obesity trends, and demographic factors. China ranks among the highest globally in terms of YLDs for osteoarthritis, reaching 0.59 million (95% UI: −0.05 to 1.71) in 2021, with knee osteoarthritis accounting for a substantial portion of this burden. China’s rapid urbanization and aging population, combined with a steady increase in obesity rates, are contributing to the rising prevalence of osteoarthritis. Although China’s obesity rates are not as high as those in high-income countries, they have been rising steadily over the past three decades, leading to an increase in osteoarthritis cases, particularly in older populations. This trend underscores the positive correlation between BMI and osteoarthritis, as higher obesity rates lead to more weight-bearing joint stress, contributing to higher osteoarthritis prevalence. The United States, with its high SDI and a high proportion of individuals with elevated BMI, also exhibits one of the highest osteoarthritis burdens globally, especially in the hip and knee osteoarthritis subtypes. The availability of extensive healthcare resources and diagnostic capabilities in the United States has led to higher diagnostic rates, further amplifying the observed osteoarthritis burden. The rising rates of obesity in the US have directly contributed to the increasing prevalence of osteoarthritis, reflecting the strong association between BMI and the disease. In contrast, India has the fastest-growing burden of osteoarthritis, as indicated by the highest EAPC at 2.64 (95%CI: 2.59 to 2.69). While India has a lower SDI and obesity rates compared to high-income countries, rapid urbanization, changing diets, and an aging population are driving the increase in osteoarthritis cases. The interaction between urbanization and rising obesity rates in India suggests that, as the population adopts more sedentary behaviors and less healthy diets, the prevalence of osteoarthritis is likely to continue to rise. Addressing obesity through public health interventions that promote active lifestyles and healthy diets will be crucial in managing the osteoarthritis burden in India. These country-level analyses highlight the importance of considering regional and national socio-economic factors in understanding the burden of osteoarthritis. Countries with higher SDI and obesity rates face an increasing burden due to the direct link between BMI and osteoarthritis, while countries like India, which are undergoing rapid socio-economic transitions, may experience a growing osteoarthritis burden as obesity rates rise alongside urbanization and demographic shifts.

Differences in gender and age groups

This study found that the net drift for both males and females was positive, indicating a continuous increase in osteoarthritis among individuals aged 60 and above, with a faster rate of burden increase in males. The incidence of hip osteoarthritis in males showed particularly notable growth, especially in the 60–69 and over 85 age groups. This may be related to gender differences in physical activity patterns, workload, and physiological factors in these age brackets. For knee osteoarthritis, the YLD rates increased significantly between the ages of 60 and 80 but tended to stabilise or slightly decline after 80 years. This could be associated with changes in survival rates, health management in osteoarthritis patients, and the influence of other comorbidities in the older adult. Additionally, the burden of hip osteoarthritis in females aged 60 and above was significantly lower than in males, possibly due to differences in bone density, hormone levels, and higher health awareness and more proactive health behaviours among women.

Cohort and period effects

Cohort and period effect analyses showed that the CRR for the global population, as well as for males and females, exhibited an upward trend, indicating an increasing incidence of osteoarthritis across different generations. This trend was particularly evident in populations with high BMI, likely due to the rising obesity rates, reduced physical activity, and unhealthy dietary patterns associated with modern lifestyles. Additionally, the continuous increase in PRR suggests that the relative risk of osteoarthritis and its subtypes has been growing globally, further emphasising its significance as a global public health issue.

Limitations

Despite providing a comprehensive analysis of the global and regional trends in osteoarthritis burden, this study has several limitations. First, the reliance on data from the GBD study means that data quality and comparability may vary due to differences in national reporting systems. Data from some low-income countries may be incomplete, potentially introducing bias into the results. Second, the focus on high BMI as the sole risk factor limits the scope of the analysis, as other important risk factors, such as genetics, occupational exposure, and lifestyle, were not considered. This may underestimate the overall burden of osteoarthritis. Third, while YLDs reflect the impact of osteoarthritis on quality of life, they do not fully capture the broader socio-economic consequences of the disease, such as healthcare costs and productivity loss. Finally, the complexity of cohort and period effects poses challenges for interpretation. Although the APC model can reveal temporal dynamics, its results are influenced by multiple factors, introducing uncertainty and complexity to the analysis.

Conclusion

Overall, this study reveals a significant increase in the burden of osteoarthritis among older adults attributable to high BMI over the past 32 years, both globally and across different SDI regions, particularly in countries such as China and India. Although high-SDI regions have a higher baseline burden, their slower growth rate indicates the effectiveness of existing management and prevention measures. The findings highlight the importance of controlling obesity, strengthening osteoarthritis prevention and management, and provide valuable insights for global public health policy development. Further research is needed to explore the multi-dimensional impact of osteoarthritis and optimise public health interventions to reduce the burden of this chronic disease among the global older adult population.

Data availability statement

Publicly available datasets were analyzed in this study. This data can be found here: data used for the analyses are publicly available from the Institute of Health Metrics and Evaluation (https://vizhub.healthdata.org/gbd-results/).

Author contributions

ZL: Supervision, Visualization, Writing – review & editing. YC: Data curation, Formal analysis, Resources, Software, Writing – original draft, Writing – review & editing. ZS: Data curation, Formal Analysis, Resources, Writing – original draft, Writing – review & editing.

Funding

The author(s) declare that no financial support was received for the research and/or publication of this article.

Acknowledgments

As an unofficial GBD collaborator, we would like to thank the Institute for Health Metrics and Evaluation and the countless individuals who have contributed to GBD Study 2021 in various capacities. This analysis fully complies with the conditions of use of the data in the GBD databases.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Gen AI was used in the creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpubh.2025.1518572/full#supplementary-material

SUPPLEMENTARY FIGURE 1 | In 2021, the ASYR of osteoarthritis (A), hip osteoarthritis (B), and knee osteoarthritis (C) among individuals aged 60 years and above across 204 countries and territories. ASYR, Age-standardized Years lived with disability rate.

SUPPLEMENTARY FIGURE 2 | In 2021, the YLDs of osteoarthritis (A), hip osteoarthritis (B), and knee osteoarthritis (C) among individuals aged 60 years and above across 204 countries and territories. YLDs, Years Lived with Disability

Abbreviations

APC, Age-Period-Cohort; ASR, Age-standardized Rate; ASYR, Age-standardized Years lived with disability rate; BMI, Body Mass index; CRR, Cohort rate ratios; CI, Confidence interval; EAPC, Estimated average percentage change; GBD, Global burden of disease; GHDx, Global health data exchange; PRR, Period rate ratios; SDI, Socio-demographic Index; UI, Uncertainty intervals; YLDs, Years lived with disability.

Footnotes

References

1. Zheng, L, Lu, Z, Xu, G, Niu, X, and Zhao, J. Dual-targeted disease-modifying therapies for osteoarthritis. Lancet. (2024) 403:2591. doi: 10.1016/S0140-6736(24)00475-6

PubMed Abstract | Crossref Full Text | Google Scholar

2. Jiang, T, Weng, Q, Liu, K, He, H, Zhang, Y, Zhang, W, et al. Bidirectional association identified between synovitis and knee and hand osteoarthritis: a general population-based study. Lancet Reg Health West Pac. (2024) 50:101169. doi: 10.1016/j.lanwpc.2024.101169

PubMed Abstract | Crossref Full Text | Google Scholar

3. Katz, JN, Arant, KR, and Loeser, RF. Diagnosis and treatment of hip and knee osteoarthritis: a review. JAMA. (2021) 325:568–78. doi: 10.1001/jama.2020.22171

PubMed Abstract | Crossref Full Text | Google Scholar

4. Dogruoz, F, Yapar, A, Buyukarslan, V, Egerci, OF, Etli, I, and Kose, O. Circumferential patellar denervation does not reduce anterior knee pain in total knee arthroplasty without patellar resurfacing; a prospective comparison. J Orthop Surg Res. (2024) 19:653. doi: 10.1186/s13018-024-05161-5

PubMed Abstract | Crossref Full Text | Google Scholar

5. GBD 2019 Risk Factors Collaborators. Global burden of 87 risk factors in 204 countries and territories, 1990-2019: a systematic analysis for the global burden of disease study 2019. Lancet. (2020) 396:1223–49. doi: 10.1016/S0140-6736(20)30752-2

PubMed Abstract | Crossref Full Text | Google Scholar

6. Ho, J, Mak, C, Sharma, V, To, K, and Khan, W. Mendelian randomization studies of lifestyle-related risk factors for osteoarthritis: a PRISMA review and Meta-analysis. Int J Mol Sci. (2022) 23:1906. doi: 10.3390/ijms231911906

PubMed Abstract | Crossref Full Text | Google Scholar

7. Wang, Y, Wang, J, Chen, S, Li, B, Lu, X, and Li, J. Different changing patterns for stroke subtype mortality attributable to high sodium intake in China during 1990 to 2019. Stroke. (2023) 54:1078–87. doi: 10.1161/STROKEAHA.122.040848

PubMed Abstract | Crossref Full Text | Google Scholar

8. GBD 2021 Osteoarthritis Collaborators. Global, regional, and national burden of osteoarthritis, 1990-2020 and projections to 2050: a systematic analysis for the global burden of disease study 2021. Lancet Rheumatol. (2023) 5:e508–22. doi: 10.1016/S2665-9913(23)00163-7

PubMed Abstract | Crossref Full Text | Google Scholar

9. Cao, F, Xu, Z, Li, XX, Fu, ZY, Han, RY, Zhang, JL, et al. Trends and cross-country inequalities in the global burden of osteoarthritis, 1990-2019: a population-based study. Ageing Res Rev. (2024) 99:102382. doi: 10.1016/j.arr.2024.102382

PubMed Abstract | Crossref Full Text | Google Scholar

10. GBD 2019 Australia Collaborators. Pre-COVID life expectancy, mortality, and burden of diseases for adults 70 years and older in Australia: a systematic analysis for the global burden of disease 2019 study. Lancet Reg Health West Pac. (2024) 47:101092. doi: 10.1016/j.lanwpc.2024.101092

PubMed Abstract | Crossref Full Text | Google Scholar

11. Moseng, T, Vliet Vlieland, TPM, Battista, S, Beckwée, D, Boyadzhieva, V, Conaghan, PG, et al. EULAR recommendations for the non-pharmacological core management of hip and knee osteoarthritis: 2023 update. Ann Rheum Dis. (2024) 83:730–40. doi: 10.1136/ard-2023-225041

PubMed Abstract | Crossref Full Text | Google Scholar

12. GBD 2019 MSK in Adolescents Collaborators. Global pattern, trend, and cross-country inequality of early musculoskeletal disorders from 1990 to 2019, with projection from 2020 to 2050. Med. (2024) 5:943–962.e6. doi: 10.1016/j.medj.2024.04.009

PubMed Abstract | Crossref Full Text | Google Scholar

13. Stevens, GA, Alkema, L, Black, RE, Boerma, JT, Collins, GS, et al. Guidelines for accurate and transparent health estimates reporting: the GATHER statement. Lancet. (2016) 388:e19–23. doi: 10.1016/S0140-6736(16)30388-9

PubMed Abstract | Crossref Full Text | Google Scholar

14. Gao, Y, Liu, K, and Fang, S. Changing patterns of stroke and subtypes attributable to high systolic blood pressure in China from 1990 to 2019. Stroke. (2024) 55:59–68. doi: 10.1161/STROKEAHA.123.044505

PubMed Abstract | Crossref Full Text | Google Scholar

15. Araujo-Castillo, RV, Culquichicon, C, and Solis Condor, R. Burden of disease due to hip, knee, and unspecified osteoarthritis in the Peruvian social health insurance system (EsSalud), 2016. F1000Res. (2020) 9:238. doi: 10.12688/f1000research.22767.2

PubMed Abstract | Crossref Full Text | Google Scholar

16. Allen, KD, Thoma, LM, and Golightly, YM. Epidemiology of osteoarthritis. Osteoarthr Cartil. (2022) 30:184–95. doi: 10.1016/j.joca.2021.04.020

PubMed Abstract | Crossref Full Text | Google Scholar

17. Kolasinski, SL, Neogi, T, Hochberg, MC, Oatis, C, Guyatt, G, Block, J, et al. 2019 American College of Rheumatology/Arthritis Foundation guideline for the Management of Osteoarthritis of the hand, hip, and knee. Arthritis Care Res. (2020) 72:149–62. doi: 10.1002/acr.24131

PubMed Abstract | Crossref Full Text | Google Scholar

18. Jevsevar, DS, Brown, GA, Jones, DL, Matzkin, EG, Manner, PA, Mooar, P, et al. The American Academy of Orthopaedic Surgeons evidence-based guideline on: treatment of osteoarthritis of the knee, 2nd edition. J Bone Joint Surg Am. (2013) 95:1885–6. doi: 10.2106/00004623-201310160-00010

PubMed Abstract | Crossref Full Text | Google Scholar

19. Li, J, Tsai, TY, Clancy, MM, Li, G, Lewis, CL, and Felson, DT. Weight loss changed gait kinematics in individuals with obesity and knee pain. Gait Posture. (2019) 68:461–5. doi: 10.1016/j.gaitpost.2018.12.031

PubMed Abstract | Crossref Full Text | Google Scholar

20. Joseph, GB, McCulloch, C, Nevitt, MC, Lynch, J, Lane, NE, and Link, TM. Effects of weight change on knee and hip radiographic measurements and pain over four years: data from the osteoarthritis initiative. Arthritis Care Res. (2023) 75:860–8. doi: 10.1002/acr.24875

PubMed Abstract | Crossref Full Text | Google Scholar

21. Dai, H, Alsalhe, TA, Chalghaf, N, Riccò, M, Bragazzi, NL, and Wu, J. The global burden of disease attributable to high body mass index in 195 countries and territories, 1990-2017: an analysis of the global burden of disease study. PLoS Med. (2020) 17:e1003198. doi: 10.1371/journal.pmed.1003198

PubMed Abstract | Crossref Full Text | Google Scholar

Keywords: osteoarthritis, disease burden, BMI, older adults, Sociodemographic status, GBD 2021

Citation: Li Z, Chen Y and Shen Z (2025) Global shifts in osteoarthritis subtype trends among older adults due to elevated BMI: an age-period-cohort analysis based on the global burden of disease database. Front. Public Health. 13:1518572. doi: 10.3389/fpubh.2025.1518572

Received: 28 October 2024; Accepted: 31 March 2025;
Published: 28 April 2025.

Edited by:

Assunta Pozzuoli, University of Padua, Italy

Reviewed by:

Lichun Qiao, Xi’an Jiaotong University, China
Chen Chao, Albert Einstein College of Medicine, United States
Meaad Elbashir, Jazan University, Saudi Arabia
Keshu Cai, The First Affiliated Hospital of Nanjing Medical University, China

Copyright © 2025 Li, Chen and Shen. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Zhao Shen, MTQ1MDA4MThAcXEuY29t; Zhiqiang Li, ZG9jdG9ybGkxMDAzQDE2My5jb20=

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.