Skip to main content

Secular trend of mortality and incidence of rheumatoid arthritis in global ,1990–2019: an age period cohort analysis and joinpoint analysis

Abstract

Background

Rheumatoid arthritis (RA) is a major public health problem. Unfortunately, there is a scarcity of comprehensive and up-to-date information regarding the burden of RA and its dynamic trends in subsequent years. To examine the changing trends in the global burden of RA and forecast for 2044, which will facilitate the development of strategies tailored to RA burden and provide reference for the development of effective treatment guidelines.

Methods

Following the general analytical strategy used the Global Burden of Disease Study (GBD) 2019, which included 204 countries, the age-standardized incidence rate (ASIR), age-standardized mortality rate (ASMR) and age-standardized disability adjusted of life year (DALY) rate for RA were analyzed.

Results

The ASIR, ASMR and age-standardized DALY rate for RA in 2019 were 13.001/100,000 (95% UI, 11.833 ~ 14.274), 0.574/100,000 (95% UI, 0.356 ~ 0.793) and 39.565/100,000 (95% UI, 49.529 ~ 30.508), respectively. America had the highest ASIR [18.578(95% UI, 17.147 ~ 20.148)] and age-standardized DALY rate [53.676(95% UI, 40.106 ~ 67.968)] in 2019. Asia had the highest ASMR [0.681(95% UI, 0.802 ~ 0.480)] in 2019. From 1990 to 2019, a significant average annual percentage change (AAPC) in the ASIR was observed in both males [0.237% (95% CI, 0.216 ~ 0.259%)] and females [0.197% (95% CI, 0.141 ~ 0.254%)], AAPC in the ASMR was observed in both males [-0.398% (95% CI, -0.605~-0.191%)] and females [-0.295% (95% CI, -0.424~-0.65%)]. Age effects indicated that the relative risk (RR) of RA-associated incidence and mortality rates increased with age among males and females. The RR of RA increased over time and started to gradually increase from 1990. Cohort effects showed decreases in incidence, mortality and DALY rates in successive birth cohorts. The global incidence of RA would continue to increase in the future, while mortality would continue to decrease.

Conclusion

The increased risk of RA is dominantly influenced by age effects and period effects and the ethnic area. The results suggest that early identification and treatment of RA is important for reducing the ongoing burden with age, and targeted health education and specific intervention programs should be promoted to control middle-elderly population.

Peer Review reports

Background

Rheumatoid arthritis (RA) is a systemic autoimmune disease with symmetrical inflammatory polyarthritis as the main clinical manifestation [1]. It usually begins in the small joints of the hands and feet, then spreads to the larger joints, eventually leading to disabling joint destruction and premature mortality [2, 3]. Despite medicine advances, the disease burden of RA remains high. RA also causes many health care costs, as it has been an incurable disease until now. The world’s population is aging, with the World Health Organization (WHO) estimating that the number of individuals aged 65 years or older will increase from 524 million in 2010 to 1.5 billion in 2050 [4]. As a result of this increased longevity and a rising frequency of chronic disease risk factors, chronic disease prevalence is projected to rise steadily. The consequences are more severe and include death and disability, reduced quality of life, and increased health care utilization and costs [5]. This is particularly obvious in RA, a disease that may accelerate aging, portend poor long-term outcomes, and have significant financial implications [6,7,8]. It follows that RA is a major public health problem, associated with a substantial burden of functional disability.

Based on the analysis of previous disease burden data, forecasting the future incidence rate can better guide the prevention of RA and promote the rational treatment. To this end, this study described and analyzed global RA incidence, mortality and disability-adjusted of life years (DALYs) from 1990 to 2019. Meanwhile, an age-period-cohort model with an intrinsic estimator (IE) is established based on these data, and the incidence and mortality rates of RA in the next 30 years are predicted, which can provide a reference for the adjustment of global intervention and treatment policies for RA.

Materials and methods

Data sources

Contributed by the Institute for Health Metrics and Evaluation (IHME: http://www.healthdata.org/gbd/2019), the Global Burden of Disease 2019 (GBD 2019) aims to quantify the comparative magnitude of health loss due to diseases, injuries, and risk factors by age, sex, and geographies for specific points in a series of times. Estimates such as incidence, mortality and DALYs were conducted among groups stratified by age, sex and year on the GBD [9, 10]. The study provides an improved standardized methodology and comprehensive assessment of incidence, mortality and DALYs for 369 diseases and injuries and 87 risk factors in 204 countries and territories.

RA events were diagnosed and classified based on the International Statistical Classification of Diseases and Related Health Problems, Tenth Revision (ICD-10). For RA, ICD-10 codes are M05.3, M05.8, M05.9, M06.0, M06.8, M06.9 and M08.0. The modelling strategy for estimating RA incidence, mortality and DALYs was generalized by grouping countries based on quality and types of data available [11]. Our data on RA were obtained from the global health data exchange (GHDx) section result tools (https://vizhub.healthdata.org/gbd-results/). Ethical approval was not needed for this study because there was no direct involvement of human subjects.

Data analysis

Join point regression analysis

The annual percentage change (APC) and the average annual percentage change (AAPC) for each segment were estimated by join point regression, which focused on estimating the temporal trends in the age-standardized incidence rate (ASIR), age-standardized mortality rate (ASMR) and age-standardized DALY rate of RA. We applied the join point regression program version 4.8.0.1 from the Statistical Research and Applications Branch of the Surveillance Research Program of the U.S. National Cancer Institute.

APC analysis

This study used the APC framework to assess the incidence, mortality and DALY rates of RA globally and to assess the potential impact of age, period, and cohort effects on these trends. The APC model was developed based on the Poisson distribution, and can be generally expressed as follows [12].

Y = log(M) =µ + α agej + β periodj + γ cohortj + εi, where M is the incidence, mortality, or DALY rate of the corresponding age group; α, β, and γ denote the coefficients of age, period, and cohort of the model, respectively; µ denotes the intercept of the model; and εi denotes the residual of the model.

The time series from 1990 to 2019 was chopped into six consecutive intervals. Successive 5-year age groups from 0 ~ 14 to 89 ~ 94 years were divided. Individuals younger than 9 and older than 95 were excluded. The risk ratio (RR) was the exponential value of the coefficient, which denoted the risk in a particular age, period, or birth cohort relative to the reference groups. The specific value of the middle group was selected to represent the whole group. All analyses were conducted by Stata 15 software (Stata Corp, College Station, TX, USA). The Wald chi-square test was utilized to calculate the significance of the estimable parameters and functions. All statistical tests were 2-sided, and the significance level considered was 5%. Model fitting was evaluated by deviance, SE coefficients, Akaike’s information criterion (AIC), and the Bayesian information criterion (BIC).

Incidence and mortality rate predicted based on bayesian age-period cohorts

This study predicted the rate of incidence and mortality due to RA from 2019 to 2044 by running an APC analysis by sex using the Nordpred package in R, taking into account the changing rates and changing population structure, which has been fully demonstrated and recognized in previous studies.

Results

Global level

The incidence rate

The ASIR of RA showed a relatively stable trend from 12.208/100,000 (95% UI, 11.126 ~ 13.376) in 1990 to 13.001/100,000 (95% UI, 11.833 ~ 14.274) in 2019. Both sexes followed a similar pattern, but the rates were higher in females than in males across all groups. The ASIR increased from 7.590/100,000 (95% UI, 6.904 ~ 8.336) in 1990 to 8.133/100,000 (95% UI, 7.380 ~ 8.944) in 2019 for males, and from 16.79% (95% UI, 15.355 ~ 18.343%) in 1990 to 17.80/100,000 (95% UI, 16.218 ~ 19.523) in 2019 for females.

Before 2015, the ASIR of RA showed a slighter upward trend globally, and an opposite trend was observed after 2015. The ASIR decreased from 8.223/100,000 (95% UI, 7.488 ~ 9.014) in 2015 to 8.133/100,000 (95% UI, 7.380 ~ 8.944) in 2019 in males, and decreased from 18.272/100,000 (95% UI, 16.672 ~ 19.934) in 2015 to 17.715/100,000 (95% UI, 16.133 ~ 19.433) in 2018 in females (Fig. 1A).

Fig. 1
figure 1

The trend of the RA in global between 1990 and 2019. (A): ASIR; (B): ASMR; (C): Age-standardized DALY rate

The mortality rate

The ASMR of RA showed a fluctuating downward trend, from 0.627/100,000 (95% UI, 0.491 ~ 0.775) in 1990 to 0.566/100,000 (95% UI, 0.431 ~ 0.660) in 2019. The ASMR decreased from 0.452/100,000 (95% UI, 0.354 ~ 0.499) in 1990 to 0.404/100,000 (95% UI, 0.350 ~ 0.457) in 2019 in males, and from 0.759/100,000 (95% UI, 1.029 ~ 0.553) in 1990 to 0.697/100,000 (95% UI, 0.465 ~ 0.859) in 2019 in females (Fig. 1B).

The DALY rate

From 1990 to 2019, the age-standardized DALY rate showed a slightly fluctuating upward trend, from 39.116/100,000 (95% UI, 30.132 ~ 48.565) in 1990 to 39.565/100,000 (95% UI, 49.529 ~ 30.508) in 2019, which increased from 1.103% for males and 1.699% for females (Fig. 1C).

Regional level

At the regional-level, the ASIR of RA was found to be highest in America, with 18.578 (95% UI, 17.147 ~ 20.148) in 2019. Whereas, Africa showed the lowest ASIR, with 9.351 (95% UI, 8.374 ~ 10.393) in 2019 (Table 1). The ASMR was found to be highest in Asia, with 0.681 (95% UI, 0.802 ~ 0.480) in 2019. However, Africa showed the lowest rate, with 0.348 (95% UI, 0.428 ~ 0.212) in 2019 (Table 1). The age-standardized DALY rate was found to be highest in America, with 53.676 (95% UI, 40.106 ~ 67.968). Africa showed the lowest age-standardized DALY rate, with 27.024 (95% UI, 20.471 ~ 34.619) (Table 1).

Table 1 Rate of RA in 1990 and 2019 and age-standardized rates by region (95% uncertainty interval)

The ASIR showed an increasing trend (Fig. 2A), while the ASMR showed a decreasing trend (Fig. 2B). Except for Europe, the age-standardized DALY rate showed an upward trend (Fig. 2C). The percentage change in regional-level ASIR, ASMR and age-standardized DALY rate estimates for all GBD 2019 regions are presented between 1990 and 2019 in online supplemental Figure S1.

Fig. 2
figure 2

The trend of the RA in regional between 1990 and 2019. (A): ASIR; (B): ASMR; (C): Age-standardized DALY rate

Join point regression analysis

Generally, the ASIR shown an overall fluctuating upward trend, but decreased for both males and females from 2014 to 2017. The AAPC in the whole period was 0.237% (95% CI, 0.216 ~ 0.259%) and 0.197% (95% CI, 0.141 ~ 0.254%) for males and females, respectively. The ASMR decreased from 2002 to 2019, except for a flat trend among males, during the periods of 2002 to 2008 and 2012 to 2019. While in females, the ASMR substantially decreased from 2002 to 2010 and then maintained a flat trend from 2010 to 2014. Throughout the whole period, the AAPC was − 0.398% (95% CI, -0.605~-0.191%) and 0.295% (95% CI, -0.424~-0.165%) for males and females, respectively.

The age-standardized DALY rate substantially increased from 1996 to 2004 and then maintained a flat trend during the period from 2004 to 2014 in males. However, in females, the age-standardized DALY rate substantially increased from 1997 to 2003 and then maintained a flat trend during the period from 2003 to 2011. Throughout this period, the AAPC was 0.034% (95% CI, -0.072~-0.039%) and 0.058% (95% CI, 0.011 ~ 0.105%) for males and females, respectively. Join point regression analysis of the sex-specific, ASIR, ASMR and age-standardized DALY rate for RA in global from 1990 to 2019 is shown in Table 2; Fig. 3.

Table 2 AAPC and APC of Incidence, Mortality and DALY rate of RA among males and females
Fig. 3
figure 3

Join point regression analysis for RA between 1990 and 2019. (A): ASIR for male; (B): ASIR for female; (C): ASMR for male; (D): ASMR for female; (E): Age-standardized DALY rate for male; (F): Age-standardized DALY rate for female. *: Indicates that the APC is significantly different from zero at the alpha = 0.05 level

APC analysis

Age effect

The age effect on the ASIR of RA is shown in Table 3. The minimum age of the male RR was 0.137 (95% CI, 0.054 ~ 0.345) in the 4 ~ 9-year-old group, and the maximum age of the male RR was 2.630 (95% CI, 2.164 ~ 3.197) in the 65 ~ 69- year-old group. The highest RR was 2.799 (95% CI, 2.418 ~ 3.239) in the 60 ~ 64 age group, which was 24.129-fold higher than that in the 4 ~ 9 age group (Fig. 4A).

Table 3 Specific relative risk of RA incidence, mortality, and DALY rate due to age, period, and cohort effects
Fig. 4
figure 4

APC Analysis on RA for both sexes: Age effect on RA for both sexes. (A): ASIR; (B): ASMR; (C): Age-standardized DALY rate;Period effect on RA for both sexes. (D): ASIR; (E): ASMR; (F): Age-standardized DALY rate; Birth cohort effect on RA for both sexes. (G): ASIR; (H): ASMR; (I): Age-standardized DALY rate

The effect of age on RA mortality showed no significant difference in ASMR before 64 years of age. The age effect exhibited an increasing trend from 5 ~ 89 age in both sexes and reached the highest in the 85 ~ 89 age group of males, in which the RR was 16.086 (95% CI, 1.716 ~ 150.776), and reached the highest in the 90 ~ 94 age group of females, in which the RR was 12.807 (95% CI, 1.397 ~ 117.42) (Fig. 4B).

The trend on age-standardized DALY rate of the age effect was similar in both sexes. From the 9 ~ 14 age group to the 64 ~ 69 age group, which showed an upward trend, while from the 74 ~ 79 to 89 ~ 94 age groups, it dropped slightly (Figure 4C).

Period effect

The period effect on incidence, mortality and DALYs is shown in Fig. (4D-4F). RR for period effects was estimated separately after controlling for the other two effects. Both sexes had similar patterns of period effects in incidence, mortality and DALYs. Compared with 1994, the incidence rate in 2019 increased 1.361-fold for males and 1.259-fold for females. The effect of period on RA mortality rate showed no significant difference on ASMR. The age-standardized DALY rate of RR in 2019 was 1.277 (95% CI, 1.174 ~ 1.390) for males and 1.275 (95% CI, 1.200 ~ 1.354) for females.

Birth cohort effect

After controlling for the other two effects, the birth cohort RR revealed a slowly downward trend from the 1990 ~ 1995 cohort to the 2010 ~ 2014 cohort (Fig. 4G and I). In total, the incidence rate of the RA showed a slow downward trend, both in males and females, but with a lower level in males than in females after birth cohort 1945. However, the effect of cohort on the incidence rate of RA has not been significantly different since 1930. The cohort effect on mortality rate in RA is not significantly different across the entire interval. The influence of male and female cohorts on DALY rate of RA was similar in the whole interval.

Morbidity and mortality projections

We predicted that the future incidence rate of RA would continue to increase while the mortality rate would continue to decrease globally, especially among males. Although the trends of incidence and mortality rates in females were somewhat similar to those in males, the incidence and mortality rates were substantially higher than those in males. Furthermore, the incidence rate was predicted to relatively stabilize after 2030, and the mortality rate was predicted to remain at a relatively low level, largely due to progress in a series of therapeutic progress and population increased awareness of the disease (Fig. 5).

Fig. 5
figure 5

Forecasting global trends in RA. (A): ASIR; (B): ASMR

Discussion

The study explores the long-term trends in global RA ASIR, ASMR and age-standardized DALY rates of RA from 1990 to 2019 using the age-period-cohort framework based on data from the GBD 2019 study. The ASIR of RA showed a relatively stable trend from 12.208/100,000 (95% UI, 11.126 ~ 13.376) in 1990 to 13.001/100,000 (95% UI, 11.833 ~ 14.274) in 2019. The ASMR of RA showed a fluctuating downward trend. The age-standardized DALY rate showed a slight fluctuating upward trend. These data were consistent with the findings reported by W. Yang et al. [13]. Furthermore, we employed a join point regression analysis and developed an APC model to further investigate the incidence and DALY rate of RA. Both sexes followed a similar pattern, but the rates were higher in females than in males across all groups by join point regression analysis. APC with IE was established to analyze the influence of age, period and cohort effect on RA. The highest RR on the ASIR for males was 2.630 (95% CI, 2.164 ~ 3.197) in the 65 ~ 69 year-old-group. The future incidence of RA will increase, while the mortality rate will continue to decrease globally according to the model of prediction. These findings are of great help to epidemiologists and clinical experts in controlling the incidence trend of RA.

The minor decrease in the mortality is similar to that reported in a previous study in Global that investigated the period from 1990 to 2010 [14]. This may be the result of clear improvement in treatment strategies and more profound knowledge of the disease course. However, a previous global study examined trends in RA incidence and DALYs from 1990 to 2017 with inconsistent results [15]. The ASIR of RA increased by 8.2% between 1990 and 2017. Additionally, the age-standardized DALY rate of RA decreased by 3.6% between 1990 and 2017. This may be because we used the 2019 GBD database and they used the 2017 database. The updated data of two years will affect the previous algorithm.

In this study, the ASIR was higher among females than among males in 2019. This was consistent with previous studies [16]. The ASIR in females was 2.22 times the ASIR in males in 2015, which indicates that more females are affected by RA than males worldwide. This is consistent with previous studies [17]. The ASMR in females was also higher than that in males, which is in line with a higher incidence in females than in males, indicating that the female predominance is not simply due to the longer lifespan of females.

At the regional level, differences in incidence, mortality and DALY rates in regions should be noted. America has the highest incidence and DALY rates of RA, yet Asia has the highest mortality rate, and Africa has the lowest incidence, mortality and DALY rates. This might be due to regional, ethnic and economic differences. Therefore, more individualized guidelines for diagnosis and management may be needed.

According to the join point regression analysis, an increase in the age-standardized rate of incidence and DALYs related to RA represents an ongoing challenge. This suggests the importance of increased investment in the prevention and treatment of rheumatoid arthritis. However, there was an obvious downward trend from 2014 to 2017. It may be related to major events that had an impact on population changes during this period, such as the Ukrainian crisis in 2014, the European refugee crisis in 2015, and the Rohingya crisis in 2017, which resulted in many people being displaced and population migrations. The Ebola outbreak in 2014 also led to thousands of infections and deaths. Additionally, the approval of new therapies can contribute to changes in the ASIR and age-standardized DALY rates. For example, the 2012 American College of Rheumatology (ACR) guidelines [18] added three biologic agents, certolizumab, tocilizumab, and golimumab, and emphasized the importance of early disease control for preventing irreversible joint damage and maintaining long-term function. The 2015 ACR guidelines [19] for the first time included tofacitinib, a new small molecule drug for the first time. Meanwhile, the decreased in age-standardised rate of mortality related to RA suggests effective control. This may be related to effective treatment and early diagnosis [20]. Moreover, we also found that males had a higher incidence of AAPCs than females in the current study. This may reveal a potential increasing burden of RA in males. While males had lower AAPCs in the ASMR due to RA than females. This may reveal two points. First, the disease may progress more rapidly and be more aggressive in females; Second, male patients may respond better to treatment.

The age effect analysis showed that RA incidence consistently increased with advancing age from 9 to 64 years in females and from 9 to 69 years in males.The peak age at RA onset in both sexes was 49 ~ 69 years age group and 54 ~ 74 years age group, respectively. However, a previous study showed that worldwide mean age of RA onset was 44 ± 14 years (95% CI, 44 ~ 45) [21], which is slightly lower than our age. This finding indicates that RA incidence significantly increased with advancing age, mainly in older people, which might be attributable to the aging population globally. The age effect analysis showed that ASMR of RA consistently increased with advancing age from 9 to 94 years in both sexes. As we all known, the most common causes of mortality among RA patients are cardiovascular diseases, respiratory diseases, and infections [22]. The risk of cardiovascular disease and infection in RA patients increases greatly with age. This may be one reason for the increase in mortality with age.

The period effect analysis showed that ASIR, ASMR and age-standardized DALY rates of RA consistently increased over time. It is usually due to a series of complex historical events and environmental factors. This may be due to environmental changes and a more rapid pace of life. Lifestyle factors have changed substantially globally, including the adoption of unhealthy eating habits [23]. Previous studies have evaluated for several factors associated with the risk of RA in prior studies, including obesity [24], smoking, periodontitis [25], vitamin D deficiency [26] breast feeding and oral contraceptive use [27, 28].

The cohort effect on the ASIR, ASMR and age-standardized DALY rates due to RA continuously decreased in later birth cohorts in both males and females. The possible reason was that the later birth cohorts received better education and had a greater awareness of health and disease prevention than earlier birth cohorts [29, 30], as well as, the continuous update and improvement of treatment strategies.

The estimated increases in the incidence rate of RA and the numbers of new cases and deaths are expected to continue to increase globally during the next 25 years due to population growth and aging. Although the incidence and mortality trends were somewhat similar between both sexes, they were much higher in females than in males. Mortality rate is expected to remain at a relatively low level, mainly due to advances in treatment and improved the understanding of the disease among people.

However, there are some limitations that should be noted. First, although the GBD 2019 adapted numerous adjustments and corrections to the source, collation, and evaluation of the RA incidence, mortality and DALY rates to enhance data accuracy and comparability, it was undeniable that certain deviations in the completeness and accuracy of the GBD data were inevitable. Second, our study lacked an analysis of the comparison of RA incidence, mortality and DALY rates in the different national and different disease subtypes. Third, similar to other studies based on a population level, the APC-IE framework concluded the results at the population level. Ecological fallacy might occur because the study might not focus on the individual level. Therefore, more relevant studies focusing on these limitations should be carried out in the future to reveal more accurate analysis results.

Conclusions

In summary, the ASIR of RA increased in globally for both males and females during the 1990 ~ 2019 period, with the peak age in the group aged 49 ~ 74 years. Due to the increase in high-risk behaviors and demographic changes in the coming decades, the incidence and mortality of RA worldwide will continue to increase. According to the regional level trend, the age-standardized incidence and DALY rates are increasing, especially in America. Therefore, an enhanced understanding of risk profiles and morbidity patterns associated with RA could facilitate the early identification of individuals, thereby reducing the future burden of this disease and better monitoring of disease burden and health outcomes. Meanwhile, more individualized guidelines for diagnosis and management may be needed.

Data Availability

The datasets analyzed during the current study are available in the (Global Health Data Exchange) query tool produced by the IHME repository, [http://ghdx.healthdata.org/gbd-results-tool].

Abbreviations

AAPC:

Average Annual Percent change

APC:

Annual Percent change

ASIR:

Age Standardized Incidence Rate

ASMR:

Age Standardized Mortality Rate

DALY:

Disability-adjusted life year

GBD:

Global Burden of Disease (study)

GHDx:

Global Health Data exchange

ICD:

International Classification of Disease

RR:

Relative Risk

RA:

Rheumatoid arthritis

WHO:

World Health Organization

AIC:

Akaike’s information criterion

BIC:

Bayesian information criterion

IE:

Intrinsic estimator

UI:

Uncertainty interval

References

  1. Sherrer YS, Bloch DA, Mitchell DM, Young DY, Fries JF. The development of disability in rheumatoid arthritis. Arthritis Rheum. 1986;29:494–500.

    Article  CAS  PubMed  Google Scholar 

  2. England BR, Sayles H, Michaud K, et al. Cause-specific mortality in male US veterans with rheumatoid arthritis. Arthritis Care Res. 2016;68:36–45.

    Article  Google Scholar 

  3. Sparks JA, Chang SC, Liao KP, et al. Rheumatoid arthritis and mortality among women during 36 years of prospective follow-up: results from the nurses’ health study. Arthritis Care Res. 2016;68:753–62.

    Article  Google Scholar 

  4. Fitzmaurice C, Abate D, Abbasi N. Et al;for the global burden of Disease Cancer Collaboration.Global, Regional, and National Cancer incidence, mortality, years of Life Lost, Years lived with disability, and disability-adjusted life-years for 29 Cancer groups, 1990 to 2017: a systematic analysis for the global burden of Disease Study. JAMA Oncol. 2019;5(12):1749–68.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Salive ME. Multimorbidity in older adults. Epidemiol Rev. 2013;35:75–83.

    Article  PubMed  Google Scholar 

  6. Crowson CS, Liang KP, Therneau TM, Kremers HM, Gabriel SE. Could accelerated aging explain the excess mortality in patients with seropositive rheumatoid arthritis? Arthritis Rheum. 2010;62:NA–82.

    Article  Google Scholar 

  7. Goronzy JJ, Shao L, Weyand CM. Immune aging and rheumatoid arthritis. Rheum Dis Clin North Am. 2010;36:297–310.

    Article  PubMed  PubMed Central  Google Scholar 

  8. Birnbaum H, Pike C, Kaufman R, Marynchenko M, Kidolezi Y, Cifaldi M. Societal cost of rheumatoid arthritis patients in the US. Curr Med Res Opin. 2010;26:77–90.

    Article  PubMed  Google Scholar 

  9. Gregory Roth LB et al. Marczak,Kalkidan Hassen Abate, Global, regional, and national age-sex-specific mortality for 282 causes of death in 195 countries and territories, 1980–2017: a systematic analysis for the Global Burden of Disease Study 2017. Lancet. (2018) 392:1736–88.

  10. Wu D, Wong P, Guo C, Tam LS, Gu J. Pattern and trend of five major musculoskeletal disorders in China from 1990 to 2017: findings from the global burden of Disease Study 2017. BMC Med. 2021;19(1):34.

    Article  PubMed  PubMed Central  Google Scholar 

  11. Jin Z, Wang D, Zhang H, et al. Incidence trend of five common musculoskeletal disorders from 1990 to 2017 at the global, regional and national level: results from the global burden of disease study 2017. Ann Rheum Dis. 2020;79(8):1014–22.

    Article  PubMed  Google Scholar 

  12. Dhamnetiya D, Patel P, Jha RP, Shri N, Singh M, Bhattacharyya K. Trends in incidence and mortality of tuberculosis in India over past three decades: a joinpoint and age-period-cohort analysis. BMC Pulm Med. 2021;21:375.

    Article  PubMed  PubMed Central  Google Scholar 

  13. Yang W, Xu Y, Lin J et al. Global, regional and national burden of rheumatoid arthritis, and attributable risk factors from 1990 to 2019: update from the global burden of Disease 2019 study. Clin Exp Rheumatol. 2022 Dec 5.

  14. Cross M, Smith E, Hoy D, et al. The global burden of rheumatoid arthritis: estimates from the global burden of disease 2010 study. Ann Rheum Dis. 2014;73(7):1316–22.

    Article  PubMed  Google Scholar 

  15. Safiri S, Kolahi AA, Hoy D, et al. Global, regional and national burden of rheumatoid arthritis 1990–2017: a systematic analysis of the global burden of disease study 2017. Ann Rheum Dis. 2019;78:1463–71.

    Article  PubMed  Google Scholar 

  16. Finckh A, Gilbert B, Hodkinson B, et al. Global epidemiology of rheumatoid arthritis Nat Rev Rheumatol. 2022;18(10):591–602.

    PubMed  Google Scholar 

  17. Megha Arora E, Gakidou Y, He, et al. Global, regional, and national age-sex-specific mortality and life expectancy, 1950–2017: a systematic analysis for the global burden of Disease Study 2017. Lancet. 2018;392:1684–735.

    Article  Google Scholar 

  18. Singh JA, Furst DE, Bharat A, et al. 2012 update of the 2008 American College of Rheumatology recommendations for the use of disease-modifying antirheumatic drugs and biologic agents in the treatment of rheumatoid arthritis. Arthritis Care Res (Hoboken). 2012;64(5):625–39.

    Article  CAS  PubMed  Google Scholar 

  19. Singh JA, Saag KG, Bridges SL Jr, et al. 2015 American College of Rheumatology Guideline for the treatment of rheumatoid arthritis. Arthritis Rheumatol. 2016;68(1):1–26. https://doi.org/10.1002/art.39480. Epub 2015 Nov 6.

    Article  PubMed  Google Scholar 

  20. Aletaha D, Smolen JS. Diagnosis and management of rheumatoid arthritis: a review. JAMA. 2018;320(13):1360–72.

    Article  PubMed  Google Scholar 

  21. GEO-RA Group. Latitude gradient influences the age of onset of rheumatoid arthritis: a worldwide survey. Clin Rheumatol. 2017;36(3):485–97.

    Article  Google Scholar 

  22. Otón T, Carmona L. The epidemiology of established rheumatoid arthritis. Best Pract Res Clin Rheumatol. 2019;33(5):101477.

    Article  PubMed  Google Scholar 

  23. Myasoedova E, Davis J, Matteson EL, et al. Is the epidemiology of rheumatoid arthritis changing? Results from a population-based incidence study, 1985–2014. Ann Rheum Dis. 2020;79(4):440–4.

    Article  PubMed  Google Scholar 

  24. Myasoedova E, Davis J, Matteson EL, Crowson CS. Contribution of obesity to the rise in incidence of rheumatoid arthritis. Arthritis Care Res. 2013;65:71–7.

    Article  Google Scholar 

  25. Chou YY, Lai KL, Chen DY, Lin CH, Chen HH. Rheumatoid arthritis risk associated with periodontitis exposure: a nationwide, population-based cohort study. PLoS ONE. 2015;10:e0139693.

    Article  PubMed  PubMed Central  Google Scholar 

  26. Merlino LA, Curtis J, Mikuls TR, Cerhan JR, Criswell LA, Saag KG. Vitamin D intake is inversely associated with rheumatoid arthritis: results from the Iowa women’s health study. Arthritis Rheum. 2004;50:72–7.

    Article  CAS  PubMed  Google Scholar 

  27. Karlson EW, Mandl LA, Hankinson SE, Grodstein F. Do breast-feeding and other reproductive factors influence future risk of rheumatoid arthritis? Results from the nurses’ health study. Arthritis Rheum. 2004;50:3458–67.

    Article  PubMed  Google Scholar 

  28. Doran MF, Crowson CS, O’Fallon WM, Gabriel SE. The effect of oral contraceptives and estrogen replacement therapy on the risk of rheumatoid arthritis: a population based study. J Rheumatol. 2004;31:207–13.

    CAS  PubMed  Google Scholar 

  29. Jack Fox SL, James KJ, Krohn K, Lau B, Miltz, Christopher JL, Murray. Global, regional, and national incidence, prevalence, and years lived with disability for 354 diseases and injuries for 195 countries and territories, 1990–2017: a systematic analysis for the global burden of Disease Study 2017. Lancet. 2018;392:1789–858.

    Article  Google Scholar 

  30. Megha Arora M, Biehl S, Deiparine C, Ikeda GR, Kemp. Hmwe Kyu. Global, regional, and national disability-adjusted life-years (DALYs) for 359 diseases and injuries and healthy life expectancy (HALE) for 195 countries and territories, 1990–2017: a systematic analysis for the global burden of Disease Study 2017. Lancet. 2018;392:1859–922.

    Article  Google Scholar 

Download references

Acknowledgements

We would like to thank the Institute for Health Metrics and Evaluation staff and its collaborators who prepared these publicly available data.

Funding

Bill and Melinda Gates Foundation funded the GBD study, but did not fund this study. This study was supported by National Key R&D Program of China (2020YFA0710800), the Major International (Regional) Joint Research Project (No.81720108020), National Natural Science Foundation of China (No. 81871283), and Nanjing Medical Science and technique Development Foundation (JQX20004).

Author information

Authors and Affiliations

Authors

Contributions

FY, MT and LG designed the study. MT analysed the data and performed the statistical analyses. FY and QL drafted the initial manuscript. ZJ, HC, LS and LG helped in revising the article. All authors contributed to the interpretation of the data, and critically revised all versions of the manuscript and approved the final version.

Corresponding authors

Correspondence to Linyu Geng or Lingyun Sun.

Ethics declarations

Ethics approval and consent to participate

This is a database study that reuses data from GBD study to solve a research problem. It is observational and does not require ethical approval and consent.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Electronic supplementary material

Below is the link to the electronic supplementary material.

Supplementary Material 1

Rights and permissions

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

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Yu, F., Chen, H., Li, Q. et al. Secular trend of mortality and incidence of rheumatoid arthritis in global ,1990–2019: an age period cohort analysis and joinpoint analysis. BMC Pulm Med 23, 356 (2023). https://doi.org/10.1186/s12890-023-02594-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12890-023-02594-2

Keywords