Skip to main content

Comprehensive evaluation of the influence of sex differences on composite disease activity indices for rheumatoid arthritis: results from a nationwide observational cohort study

Abstract

Background

The effects and their magnitudes of sex on disease activity indices for rheumatoid arthritis are not clear. We aimed to comprehensively evaluate the influence of sex on disease activity indices in the real-world setting using a large observational database.

Methods

We analyzed 14,958 patients registered in the National Database of Rheumatic Diseases in Japan (NinJa) in 2017. We evaluated the sex differences in the 28-joint disease activity score (DAS28) using erythrocyte sedimentation rate (ESR), DAS28 using C-reactive protein (DAS28-CRP), Simplified Disease Activity Index (SDAI), and Clinical Disease Activity Index by disease activity category using Cliff’s delta and regression analysis. Differences in the share of components of indices were evaluated using permutational multivariate analysis of variance. Correction equations were constructed to estimate the number of misclassification in male patients who achieve DAS28-ESR remission.

Results

DAS28-ESR showed higher values in female patients than male patients in remission despite no obvious difference in other indices or disease activity categories. Among the components of DAS28-ESR, only ESR was higher in female patients than male patients in remission. In DAS28-CRP and SDAI, 28-tender joint count was higher and CRP was lower in female patients than male patients. In addition, the profiles in the components were different between female and male patients, especially among those with high disease activity. Using correction equations, almost 12% of male patients with DAS28-ESR remission were estimated to be misclassified, mainly due to differences in ESR.

Conclusion

Among the disease activity indices, significant sex difference was observed only in DAS28-ESR remission. The degree of misclassification in DAS28-ESR remission would be unignorable.

Peer Review reports

Introduction

The concept of “Treat to Target” (T2T), which involves setting a goal and determining an appropriate treatment, has improved outcomes in patients with rheumatoid arthritis (RA) [1,2,3]. Assessing the disease activity is an important part of T2T. Disease activity is usually evaluated using scoring systems, such as 28-joint disease activity score (DAS28) using erythrocyte sedimentation rate (ESR), DAS28 using C-reactive protein (DAS28-CRP), Clinical Disease Activity Index (CDAI), and Simplified Disease Activity Index (SDAI) [4,5,6,7]. Although these indices are widely used and recommended to assess RA disease activity [8], they are influenced by sex, age, body mass index, and other factors [9,10,11,12,13,14,15,16,17,18,19,20,21]. These effects and their magnitude on disease activity indices are not clear; therefore, disease activity should be carefully interpreted based on the properties of each index and individual patient factors.

Sex differences in RA have been studied from multifaceted aspects, for example, incidence, phenotype, comorbidities, treatment response, and prognosis [22]. In addition, sex differences in disease activity indices have been studied, mainly in DAS28-ESR. Many studies have reported that DAS28-ESR is lower in male patients compared to female patients [10,11,12,13,14,15,16], and ESR is considered to contribute to the sex differences [10, 12, 13]. The association of sex difference with CRP level and discrepancy between DAS28-ESR and DAS28-CRP have been studied [16, 23,24,25,26,27], but few studies have evaluated the effect of sex differences on DAS28-CRP [16]. Although the effects of sex differences on CDAI and SDAI have not been thoroughly evaluated, these indices are influenced by pain perception and sex [9]. Therefore, the effect of sex differences on these indices cannot be ignored. Nevertheless, sex differences are not taken into account while assessing disease activity and it is not clear how sex differences influence the composite measure indices, which may lead to biased interpretation of these indices. The current problem is a lack of large-scale systematic evaluation of the effects of sex differences on disease activity indices in the era of biologics and molecular targeting therapy. Most previous studies that evaluated the effects of sex differences on disease activity indices were conducted in the 2000s and involved less than 1000 patients. Since then, the profile of drugs used in RA patients has changed. However, the effects of sex differences on disease activity indices are not well clear in patients treated with biologics or molecular targeting therapy. Furthermore, most previous studies only performed a simple comparison of the disease activity indices without stratification, which may introduce bias in the results due to the combined analysis of patients with varying disease activity. Thus, it is unclear whether the previous studies accurately evaluated the effects of sex differences on disease activity indices or the results were obtained due to differences in disease activity between the groups. For the aforementioned reasons, the effects of sex differences on disease activity indices, while taking into account the different drugs used, require a comprehensive evaluation. Therefore, we aimed to evaluate the influence of sex differences on composite measure indices and their clinical impact by analyzing the impact on each disease activity category using a large nationwide observational database. The results would enhance our understanding and allow more appropriate use of composite measure indices by taking into account the sex differences.

Patients and methods

Study population

We collected data from a nationwide observational cohort database of RA in Japan (National Database of Rheumatic Diseases in Japan; NinJa) [26, 28] in 2017. Forty-nine hospitals and institutions from all over Japan participated in the NinJa project in 2017. NinJa included RA patients diagnosed according to the standard diagnostic criteria for RA [29,30,31,32], regardless of disease duration, onset age, and treatment. Once a year, NinJa collects information about important events (e.g., hospitalization, surgical operation, malignancy, tuberculosis, herpes zoster, or childbirth) and data are arbitrarily collected at one point in the year for each patient, including 28-tender joint count (TJC28), 28-swollen joint count (SJC28), disease activity indices, Health Assessment Questionnaire-Disability Index (HAQ-DI), drug use, and joint destruction. The collected data are curated in National Hospital Organization Sagamihara National Hospital and verified in case of doubt about accuracy. Of a total of 15,185 patients registered in NinJa in 2017, 15,056 had onset of RA at age > 16 years. We analyzed 14,958 out of 15,056 patients, thereby excluding 98 patients in whom the drug used was unknown (94 patients) or not approved for RA in Japan (4 patients).

Measures and disease activity categories

We evaluated the influence of sex differences on DAS28-ESR, DAS28-CRP, CDAI, and SDAI. Patients were classified into remission, low disease activity (LDA), moderate disease activity (MDA), and high disease activity (HDA) based on DAS28-ESR, DAS28-CRP, CDAI, and SDAI, in accordance with the updated American College of Rheumatology recommendations [8].

Statistical analysis

Statistical analyses were performed using R version 4.0.3 software (R Foundation for Statistical Computing). The URLs and/or references of R packages and function are listed in the Additional file 1. The values and graphs are the results of available-case analysis (pairwise deletion), unless stated otherwise. Figures were generated using ggplot2 or car package or geom_flat_violin function. To compare continuous variables, differences in 25% trimmed mean between female and male patients were calculated and its 95% confidence intervals (CIs) were obtained by the percentile bootstrap method (5000 iterations) using simpleboot and boot packages. Cliff’s delta, non-parametric effect size, and its 95% CIs were calculated using effsize package for comparing disease activity indices and their components between female and male patients. The magnitude of effect size based on Cliff’s delta for disease activity indices are not established, thus thresholds of magnitude for the absolute values were assessed that less than 0.147, 0.147 or more and less than 0.330, 0.330 or more and less than 0.474, and 0.474 or more corresponded to negligible, small, medium, and large according to the threshold values proposed by Romano et al. [33].

Permutational multivariate analysis of variance (PERMANOVA) was performed using vegan package for examining the difference in the share of components to disease activity indices between female and male patients in each disease activity category (number of permutations: 1000). The share of components is considered to be compositional data containing essential zeros, which cannot calculate Aitchison’s distance. Thus, Bray–Curtis dissimilarity was applied for the analysis of PERMANOVA. The threshold of p-value was not defined because this study used observational database without prespecified analysis plan and sample size design to control for type I error.

The results of regression analysis were confirmed by generalized linear model (GLM) with gamma distribution using identity link and quantile (median) regression (QR) model. Although the canonical link function of GLM with gamma distribution provides a reciprocal link, we used identity link because coefficients can be interpreted in terms of the effects of independent variables on disease activity indices at the original scale. Dependent variable in GLM with gamma distribution should not contain zero or negative values. In case that the remission dataset to be analyzed in GLM contained patients whose disease activity indices were zero, a small value (1.0 × 10–15) was added to the value for all patients in the dataset. QR was performed using quantreg package.

In regression models adjusting for patient-related factors, analysis using stacked dataset imputed by chained equations were also performed in addition to the available-case analysis to confirm the robustness of the results. Multivariate imputation by chained equations were conducted using mice package. Imputation method and models are described in the Supplementary Methods (see Additional file 1). As conventional procedure of multiple imputation, the coefficients of regression models were pooled by Rubin’s rule. However, imbalance of cases between imputed datasets occurs by stratifying disease activity categories due to imputation of the disease activity indices; therefore, the results cannot be combined by Rubin’s rule. Alternatively, we analyzed the stacked imputed dataset, which was deemed as the complete data. The point estimates calculated by the stacked method are unbiased, but the confidence intervals are invalid [34]. Thus, we present the coefficients and their confidence intervals in available-case analysis as the main results.

Results

Patient characteristics

Among a total of 14,958 patients analyzed in this study, 11,916 were female patients (79.7%) and 3042 were male patients (20.3%). Table 1 presents the patient characteristics. Male patients were older and had a higher age at onset, whereas female patients had longer disease duration, higher HAQ-DI score, and higher values of disease activity indices. Male patients had higher remission rate than female patients in all disease activity indices, especially DAS28-ESR. Of the 14,958 patients, 1043 were not treated with disease modifying anti-rheumatic drugs (DMARDs), 9593 were treated with conventional synthetic DMARDs (csDMARDs) alone, 2100 were treated with tumor necrotizing factor inhibitors (TNFi), 1144 were treated with interleukin 6 receptor inhibitors (IL-6i), 772 were treated with cytotoxic T lymphocyte-associated antigen 4 immunoglobulin (CTLA-4-Ig), and 306 were treated with Janus kinase inhibitors (JAKi). Additional file 1: Tables S1–S6 present the patient characteristics by sex for each treatment type (see Additional file 1).

Table 1 Patient characteristics

Sex differences in DAS28-ESR, DAS28-CRP, CDAI, and SDAI by disease activity category

A comparison of the distributions of DAS28-ESR, DAS28-CRP, CDAI, and SDAI values between female and male patients by disease activity category showed no difference, except for remission in DAS28-ESR (Fig. 1A). Cliff’s delta also showed sex difference only in DAS28-ESR remission, with higher values in female patients than male patients (Fig. 1B). Furthermore, we evaluated the sex difference in DAS28-ESR components and explored whether the components of these indices that showed no sex differences were less sensitive to sex difference. Cliff’s delta showed that ESR was higher in female patients than male patients in remission, while other components of DAS28-ESR and ESR in other disease activity categories showed no obvious sex differences (Cliff’s delta for each component of DAS28-ESR are shown in Fig. 1C and that of other indices are shown in Additional file 1: Fig. S1; the values of Cliff’s delta are shown in Additional file 1: Table S7). Therefore, we concluded that sex difference in DAS28-ESR remission was mainly caused by sex difference in ESR. Cliff’s delta was also calculated for each treatment group (Additional file 1: Table S7), which showed higher DAS28-ESR values and ESR for female patients in the DMARDs free, csDMARDs, TNFi, and IL-6i groups compared to male patients in remission; however, the results for CTLA-4-Ig and JAKi were equivocal.

Fig. 1
figure 1

Comparison of DAS28-ESR, DAS28-CRP, CDAI, and SDAI between female and male patients. A Distributions of disease activity indices by sex. Jitter, box, and violin plots are depicted in each disease activity category by sex (red, female; blue, male). Dashed lines are drawn on the cutoff values for the disease activity categories (DAS28-ESR and DAS28-CRP: 2.6, 3.2, and 5.1. CDAI: 2.8, 10, and 22. SDAI: 3.3, 11, and 26). B, C Cliff’s delta for sex difference in disease activity indices. Positive values of Cliff’s delta indicate that the values of the indices and their components were higher in female patients compared to male patients, whereas negative values indicate the opposite. Black line is drawn at the value of 0.000. Blue, green, and red dashed lines are drawn on the values of 0.147, 0.330, and 0.474, respectively. Points and bars indicate the estimates and 95% CIs of Cliff’s delta, respectively. Rem, remission; LDA, low disease activity; MDA, moderate disease activity; HDA, high disease activity

Sex differences were also observed in TJC28 and CRP in DAS28-CRP and SDAI (Additional file 1: Fig. S1, Table S7). In DAS28-CRP and SDAI, TJC28 tended to be higher in female patients compared to male patients in LDA to MDA and showed obvious difference in HDA. In contrast, CRP tended to be higher in male patients compared to female patients in remission to MDA and showed obvious difference in HDA in DAS28-CRP and SDAI. No obvious differences were observed in the other components of DAS28-CRP and SDAI. In each treatment group, similar results were observed in DMARDs free and csDMARDs groups (Additional file 1: Table S7). However, in the other treatment groups, the sex differences in TJC28 and CRP were equivocal or could not be evaluated due to small sample size.

Profile in the share of components of disease activity indices

Considering the contrasting dynamics of TJC28 and CRP as well as no obvious differences in the other components of DAS28-CRP and SDAI, DAS28-CRP and SDAI had no sex differences because TJC28 and CRP cancel each other. Therefore, the share of each component of the composite measures would be different between female and male patients, even if the values of the composite measures are the same. We compared the share of components between female and male patients. The share of each component to the total value of the composite measures was calculated using the procedure described by Radovits et al. [10], after excluding patients whose disease activity index value was zero. The share of components of each composite measure by disease activity category is presented in Fig. 2 and was similar between female and male patients, except for HDA in DAS28-CRP and SDAI. In HDA, the share of TJC28 in DAS28-CRP and SDAI was slightly lower in male patients compared to female patients. Distributions of the share supported that the share of TJC28 tended to be higher in female patients compared to male patients and that of CRP was the opposite in the two indices (Additional file 1: Fig. S2). Furthermore, we performed PERMANOVA to quantify the impact of sex on the share (Additional file 1: Table S8). In HDA, the R2 values in DAS28-CRP and SDAI were relatively higher than others, while sex only partly explained the variance (DAS28-CRP, R2 = 0.06577; SDAI, R2 = 0.02520). Thus, we concluded that the composition was similar between female and male patients, but the component profile in HDA slightly differed in DAS28-CRP and SDAI.

Fig. 2
figure 2

Share of the components of disease activity indices by sex. The stacked bar plot shows the share in each patient. The horizontal axis indicates the number of patients and the vertical axis indicates percentage of share of each component. Rem, remission; LDA, low disease activity; MDA, moderate disease activity; HDA, high disease activity

Adjustment of patient-related factors by regression models

We performed multivariable regression analysis by disease activity category to adjust patient-related factors using GLM and QR model. In both GLM and QR model, sex differences were only seen in remission in DAS28-ESR after adjustment of patient-related factors (Table 2 shows only the results of “Male” variable; the results of all variables are shown in Additional file 1: Table S9). The regression analysis using stacked dataset imputed by chained equations showed similar results (Additional file 1: Table S10).

Table 2 Partial regression coefficients of sex in GLM and QR models for adjustment of patient-related factors

Estimation of the impact of sex difference on remission rate

To quantify the clinical impact of sex difference in DAS28-ESR remission, we constructed correction equations using GLM and QR model. Previous studies suggest that ESR levels change in a sex specific age-dependent manner [10, 35]. Assuming that tender and swollen joint count is surrogate index which accurately reflects disease activity under the condition ESR and DAS28-ESR may not accurately reflect disease activity, we used regression model with DAS28-ESR as a dependent variable and age, 0.56 × √(TJC28) + 0.28 × √(SJC28), sex, interaction term of sex and age, and interaction term of sex and 0.56 × √(TJC28) + 0.28 × √(SJC28) as independent variables for patients with DAS28-ESR remission. We calculated the coefficients using both GLM and QR model (Table 3). Using the results of regression analysis, we constructed the following correction equations (Eqs. 1 and 2).

$$0.825 - 0.008 \times {\text{age}} - 0.177 \times \left( {0.56 \times \surd \left( {{\text{TJC}}28} \right) + 0.28 \times \surd \left( {{\text{SJC}}28} \right)} \right)$$
(1)
$$1.021 - 0.011 \times {\text{age}} - 0.202 \times \left( {0.56 \times \surd \left( {{\text{TJC}}28} \right) + 0.28 \times \surd \left( {{\text{SJC}}28} \right)} \right)$$
(2)
Table 3 Regression analysis for construction of correction equations

If the values of equations are above zero, they were added to DAS28-ESR. If the values are zero or below, the original DAS28-ESR value was used. These equations were only applied to male patients. After correction using the above equations, sex difference was deemed to be negligible (Fig. 3A and Additional file 1: Fig. S3). Cliff’s delta indicated negligible sex difference after correction [Eq. (1), − 0.0530 (95% CI − 0.0887 to − 0.0172); Eq. (2), − 0.0446 (95% CI − 0.0800 to − 0.0090)]. Of the total 1183 male patients in DAS28-ESR remission (available-case analysis), 143 (almost 12.1%) and 141 (almost 11.9%) male patients in DAS28-ESR remission were classified as LDA by correction using Eqs. 1 and 2, respectively (Fig. 3B). Therefore, almost 12% of male patients who achieve DAS28-ESR remission criteria are estimated to be misclassified.

Fig. 3
figure 3

The effect of correction by equations on DAS28-ESR. A Relationship between DAS28-ESR, age, and joint findings before and after corrections. Three-dimensional scatter plot with fitted surface are depicted by sex (red, female; blue, male). “Joint” represents 0.56 × √(TJC28) + 0.28 × √(SJC28). B Distributions of DAS28-ESR before and after corrections. Jitter, box, and violin plots are depicted for each disease activity category by sex (red, female; blue, male). Dashed lines are drawn at the value of 2.6 in DAS28-ESR

Discussion

This study confirmed that among DAS28-ESR, DAS28-CRP, CDAI, and SDAI, sex difference in the composite measures is observed only in remission based on DAS28-ESR. In addition, sex difference in DAS28-ESR remission is mainly due to difference in ESR, which is consistent with previous studies [12, 13, 36]. However, our results showed no meaningful sex differences in other components of DAS28-ESR, which is not consistent with previous studies [12,13,14, 36]. In previous studies, tender joint count and patient global assessment values, as well as ESR, were consistently higher in female patients compared to male patients. The reason for this discrepancy is uncertain, but may result from differences in stratification and study population (e.g., remission rate or pain intensity). The previous studies compared female and male patients without stratification or with stratification by swollen joint count, while we stratified patients by disease activity based on each composite measure. Therefore, the populations analyzed in the previous studies may have included patients with various disease activities than our analysis. Both tender joint count and patient global assessment are affected by pain [13, 37]. Although the mechanisms underlying sex difference in pain perception are not fully clear, female sex is more sensitive to pain than male sex [38, 39]. Thus, pain may contribute to sex differences in composite measures. However, adjusting for individual differences in pain intensity is challenging because the most frequently used measures, Visual Analog Scale and Numeric Rating Scale, are subjective and dependent on individual pain tolerance and other factors [40].

We found that the values of composite measures do not differ between female and male patients, but profile in the share of components differ between the sexes in composite measures including CRP. Several studies reported that sex difference is not observed in CRP level in RA patients [13, 36], and it is widely believed that CRP is less sensitive to sex difference than ESR. However, some studies reported that the CRP level is higher in male patients compared to female patients, similar to the results of our study [14, 16, 23]. Although, this discrepancy needs to be evaluated in further studies, our findings suggest that results of clinical trials for patients whose disease activity is assessed as HDA by DAS28-CRP and SDAI may be carefully interpreted because of the possibility that study population is not homogeneous between female and male patients. Compared to DAS28-ESR, sex differences in DAS28-CRP, CDAI, and SDAI have not been well studied. This study revealed behaviors of DAS28-CRP, CDAI, and SDAI in sex difference using stratification by disease activity category and evaluation by effect size at the level of components. The importance of our results is that they clarify whether or not disease activity indices are affected by sex differences in each disease activity category. Evidence for robustness to sex differences is necessary for the appropriate selection of indices to be used in clinical trials and routine clinical practice. Based on the findings of our study, CDAI is the most robust to sex difference among the four indices. Our study provides a cautionary implication that composite measure indices should be carefully interpreted in light of confounders, such as sex difference, before making clinical decisions based on the values of composite measure indices.

In this study, we also found that almost 12% of male patients may be misclassified as being in remission based on DAS28-ESR, suggesting that the criterion overestimates the remission rate in male patients. Generally, remission rate is lower in female patients compared to male patients, and male sex is a predictor of remission [41, 42]. In our study, sex difference had a greater effect on remission rate defined using the DAS28-ESR criteria compared to other criteria; however, to a lesser extent, the remission rates were higher in all disease activity indices in male patients compared to female patients. Previous studies showed that the remission rate was lower in female patients compared to male patients using DAS28-ESR, but not when other criteria were used, thereby suggesting that the sex difference in remission rate based on DAS28-ESR is due to sex difference in ESR [36, 43]. Although further study is needed to determine whether sex difference in remission rate may be explained by factors other than bias caused by composite measures themselves, e.g., pathophysiological or psychosocial differences, most part of the sex difference in remission rate based on DAS28-ESR would be explained by ESR [12, 36]. Overestimation of the remission rate due to sex difference in ESR is important in the context of determining treatment outcomes. In early RA, achieving remission confer favorable radiographic, quality of life and functional outcome compared to LDA [44]. Thus, misclassification of LDA into remission would have great impact on clinical decision making in T2T era, while further studies are needed to investigate the effects of misclassification of male patients with LDA into remission due to ESR on the radiographic appearance, quality of life, and functional outcome in both early and established RA using longitudinal data.

The present study had several limitations. First, we did not analyze all the factors that affect ESR and other components of the indices, such as hematocrit level, alcohol consumption, race, ethnicity, and history of fibromyalgia, hypergammaglobulinemia, Sjögren's syndrome, or other comorbidities. Second, the correction equations were not validated in an independent dataset. Therefore, correction using these equations is not yet suitable for clinical use, and estimated misclassification rates may fluctuate especially in applying to non-Japanese populations. Moreover, the equations and estimations are based on the assumption that tender and swollen joint count accurately reflects disease activity because no unbiased gold standard measure of disease activity has been established. The definitions of disease activity and remission are complex and controversial with respect to whether they should be defined solely by the state of inflammation (e.g., joint findings, ESR, or CRP) or by the combination of patients’ overall health status (e.g., patient global assessment) and state of inflammation [45, 46]. The dual target strategy, which separately manages inflammation measures as disease activity and patient-reported outcomes as disease impact, may avoid overtreatment and improve the quality of life of patients [46, 47]. Patient-reported outcomes are not only proposed as the pillar of the dual target strategy, but are also associated with functional outcomes and sustained remission [48,49,50]. Thus, it is unclear whether the assumption that joint findings are a proxy for disease activity is optimal. However, to the best of our knowledge, this is the first report to quantitatively determine the impact of bias due to sex difference on remission using correction equations, and these are candidates to correct sex difference in remission based on DAS28-ESR.

Conclusion

Large-scale and detailed analysis, including the drug type used, of the disease activity indices showed significant sex differences only in DAS28-ESR remission, mainly due to differences in ESR. In DAS28-CRP and SDAI, the values of the composite measures showed no significant sex differences, but TJC28 was higher and CRP was lower in female patients than male patients. This sex difference in the components indicated that the profiles of male and female patients were different, especially in those with high disease activity. Furthermore, almost 12% of male patients with DAS28-ESR remission were considered to be equivalent to LDA using equations to correct the effects of sex and age differences on ESR. Our results will help to understand the properties of composite measures of disease activity and allow the appropriate selection of indices based on the sex differences.

Availability of data and materials

The datasets generated and/or analysed during the current study are not publicly available because they contain patient information and unrestricted access to them is not permitted by the ethics committee. But they are available from the corresponding author on reasonable request subject to ethical approval and with permission of all investigators participating this study.

Abbreviations

CDAI:

Clinical Disease Activity Index

CIs:

Confidence intervals

CRP:

C-reactive protein

csDMARDs:

Conventional synthetic disease modifying anti-rheumatic drugs

CTLA-4-Ig:

Cytotoxic T lymphocyte-associated antigen 4 immunoglobulin

DAS28:

28-joint disease activity score

DMARDs:

Disease modifying anti-rheumatic drugs

ESR:

Erythrocyte sedimentation rate

GLM:

Generalized linear model

HAQ-DI:

Health Assessment Questionnaire-Disability Index

HDA:

High disease activity

IL-6i:

Interleukin 6 receptor inhibitors

JAKi:

Janus kinase inhibitors

LDA:

Low disease activity

MDA:

Moderate disease activity

NinJa:

National Database of Rheumatic Diseases in Japan

PERMANOVA:

Permutational multivariate analysis of variance

QR:

Quantile regression

RA:

Rheumatoid arthritis

SDAI:

Simplified Disease Activity Index

SJC28:

28-swollen joint count

T2T:

Treat to Target

TJC28:

28-tender joint count

TNFi:

Tumor necrotizing factor inhibitors

References

  1. Stoffer MA, Schoels MM, Smolen JS, et al. Evidence for treating rheumatoid arthritis to target: results of a systematic literature search update. Ann Rheum Dis. 2016;75:16–22.

    Article  Google Scholar 

  2. Smolen JS, Aletaha D, Bijlsma JWJ, et al. Treating rheumatoid arthritis to target: recommendations of an international task force. Ann Rheum Dis. 2010;69:631–7.

    Article  PubMed  Google Scholar 

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

    Article  PubMed  Google Scholar 

  4. Prevoo MLL, van'T Hof MA, Kuper HH, et al. Modified disease activity scores that include twenty-eight-joint counts development and validation in a prospective longitudinal study of patients with rheumatoid arthritis. Arthritis Rheum. 1995;38:44–8.

    Article  CAS  PubMed  Google Scholar 

  5. Fransen J, Welsing PMJ, de Keijzer RMH, et al. Disease activity scores using C-reactive protein: CRP may replace ESR in the assessment of RA disease activity. Ann Rheum Dis. 2004;62(Suppl 1):151.

    Google Scholar 

  6. Smolen JS, Breedveld FC, Schiff MH, et al. A simplified disease activity index for rheumatoid arthritis for use in clinical practice. Rheumatology (Oxford). 2003;42:244–57.

    Article  CAS  PubMed  Google Scholar 

  7. Aletaha D, Nell VPK, Stamm T, et al. Acute phase reactants add little to composite disease activity indices for rheumatoid arthritis: validation of a clinical activity score. Arthritis Res Ther. 2005;7:796–806.

    Article  Google Scholar 

  8. England BR, Tiong BK, Bergman MJ, et al. 2019 update of the American College of Rheumatology recommended rheumatoid arthritis disease activity measures. Arthritis Care Res (Hoboken). 2019;71:1540–55.

    Article  PubMed  Google Scholar 

  9. Rintelen B, Haindl PM, Maktari A, et al. SDAI/CDAI levels in rheumatoid arthritis patients are highly dependent on patient’s pain perception and gender. Scand J Rheumatol. 2008;37:410–3.

    Article  CAS  PubMed  Google Scholar 

  10. Radovits BJ, Fransen J, van Riel PLCM, et al. Influence of age and gender on the 28-joint Disease Activity Score (DAS28) in rheumatoid arthritis. Ann Rheum Dis. 2008;67:1127–31.

    Article  CAS  PubMed  Google Scholar 

  11. Tengstrand B, Ahlmén M, Hafström I. The influence of sex on rheumatoid arthritis: a prospective study of onset and outcome after 2 years. J Rheumatol. 2004;31:214–22.

    PubMed  Google Scholar 

  12. Sokka T, Toloza S, Cutolo M, et al. Women, men, and rheumatoid arthritis: analyses of disease activity, disease characteristics, and treatments in the QUEST-RA study. Arthritis Res Ther. 2009;11(1):R7.

    PubMed  PubMed Central  Google Scholar 

  13. Leeb BF, Haindl PM, Maktari A, et al. Disease activity score-28 values differ considerably depending on patient’s pain perception and sex. J Rheumatol. 2007;34:2382–7.

    CAS  PubMed  Google Scholar 

  14. Ahlmén M, Svensson B, Albertsson K, et al. Influence of gender on assessments of disease activity and function in early rheumatoid arthritis in relation to radiographic joint damage. Ann Rheum Dis. 2010;69:230–3.

    Article  PubMed  Google Scholar 

  15. Shin S, Park EH, Kang EH, et al. Sex differences in clinical characteristics and their influence on clinical outcomes in an observational cohort of patients with rheumatoid arthritis. Joint Bone Spine. 2021;88:105124.

    Article  PubMed  Google Scholar 

  16. Iikuni N, Sato E, Hoshi M, et al. The influence of sex on patients with rheumatoid arthritis in a large observational cohort. J Rheumatol. 2009;36:508–11.

    Article  PubMed  Google Scholar 

  17. Matsui T, Kuga Y, Nishino J, et al. Comparison of composite disease activity indices for rheumatoid arthritis. Mod Rheumatol. 2011;21:134–43.

    Article  PubMed  Google Scholar 

  18. Jawaheer D, Olsen J, Lahiff M, et al. Gender, body mass index and rheumatoid arthritis disease activity: results from the QUEST-RA study. Clin Exp Rheumatol. 2010;28:454–61.

    CAS  PubMed  PubMed Central  Google Scholar 

  19. George MD, Giles JT, Katz PP, et al. Impact of obesity and adiposity on inflammatory markers in patients with rheumatoid arthritis. Arthritis Care Res (Hoboken). 2017;69:1789–98.

    Article  PubMed  Google Scholar 

  20. Ranzolin A, Brenol JCT, Bredemeier M, et al. Association of concomitant fibromyalgia with worse disease activity score in 28 joints, health assessment questionnaire, and short form 36 scores in patients with rheumatoid arthritis. Arthritis Care Res (Hoboken). 2009;61:794–800.

    Article  Google Scholar 

  21. Ruiz-Medrano E, Espinosa-Ortega HF, Arce-Salinas CA. The effect of concomitant hand osteoarthritis on pain and disease activity in patients with rheumatoid arthritis. Clin Rheumatol. 2019;38:2709–16.

    Article  PubMed  Google Scholar 

  22. Favalli EG, Biggioggero M, Crotti C, et al. Sex and management of rheumatoid arthritis. Clin Rev Allergy Immunol. 2019;56:333–45.

    Article  PubMed  Google Scholar 

  23. Siemons L, Ten Klooster PM, Vonkeman HE, et al. How age and sex affect the erythrocyte sedimentation rate and C-reactive protein in early rheumatoid arthritis. BMC Musculoskelet Disord. 2014;15:1–7.

    Article  Google Scholar 

  24. Hamann PDH, Shaddick G, Hyrich K, et al. Gender stratified adjustment of the DAS28-CRP improves inter-score agreement with the DAS28-ESR in rheumatoid arthritis. Rheumatology (Oxford). 2019;58:831–5.

    Article  Google Scholar 

  25. Castrejón I, Ortiz AM, García-Vicuña R, et al. Are the C-reactive protein values and erythrocyte sedimentation rate equivalent when estimating the 28-joint disease activity score in rheumatoid arthritis? Clin Exp Rheumatol. 2008;26:769–75.

    PubMed  Google Scholar 

  26. Matsui T, Kuga Y, Kaneko A, et al. Disease Activity Score 28 (DAS28) using C-reactive protein underestimates disease activity and overestimates EULAR response criteria compared with DAS28 using erythrocyte sedimentation rate in a large observational cohort of rheumatoid arthritis patients. Ann Rheum Dis. 2007;66:1221–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Siemons L, Vonkeman HE, Ten Klooster PM, et al. Interchangeability of 28-joint disease activity scores using the erythrocyte sedimentation rate or the C-reactive protein as inflammatory marker. Clin Rheumatol. 2014;33:783–9.

    Article  PubMed  Google Scholar 

  28. Saeki Y, Matsui T, Saisho K, et al. Current treatments of rheumatoid arthritis: from the “NinJa” registry. Expert Rev Clin Immunol. 2012;8(5):455–65.

    Article  CAS  PubMed  Google Scholar 

  29. Arnett FC, Edworthy SM, Bloch DA, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis. Arthritis Rheum. 1988;31(3):315–24.

    Article  CAS  PubMed  Google Scholar 

  30. Aletaha D, Neogi T, Silman AJ, et al. 2010 rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative. Arthritis Rheum. 2010;62(9):2569–81.

    Article  PubMed  Google Scholar 

  31. Funovits J, Aletaha D, Bykerk V, et al. The 2010 American College of Rheumatology/European league against rheumatism classification criteria for rheumatoid arthritis: methodological report phase I. Ann Rheum Dis. 2010;69(9):1589–95.

    Article  PubMed  Google Scholar 

  32. Neogi T, Aletaha D, Silman AJ, et al. The 2010 American College of Rheumatology/European League Against Rheumatism classification criteria for rheumatoid arthritis: phase 2 methodological report. Arthritis Rheum. 2010;62(9):2582–91.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Romano J, Kromrey JD, Coraggio J, et al.: Exploring methods for evaluating group differences on the NSSE and other surveys: are the t-test and Cohen’s d indices the most appropriate choices? Annu Meet South Assoc Institutional Res 2006:1–51.

  34. van Buuren S. Flexible imputation of missing data. 2nd ed. Boca Raton, FL: Chapman and Hall/CRC Press; 2018.

    Book  Google Scholar 

  35. Miller A, Green M. Simple rule for calculating normal erythrocyte sedimentation rate. Br Med J (Clin Res Ed). 1983;286(6361):266.

    Article  Google Scholar 

  36. Mäkinen H, Hannonen P, Sokka T. Sex: A major predictor of remission as measured by 28-joint Disease Activity Score (DAS28) in early rheumatoid arthritis? Ann Rheum Dis. 2008;67:1052–3.

    Article  PubMed  Google Scholar 

  37. Ward MM, Guthrie LC, Dasgupta A. Direct and indirect determinants of the patient global assessment in rheumatoid arthritis: differences by level of disease activity. Arthritis Care Res (Hoboken). 2017;69:323–9.

    Article  PubMed  Google Scholar 

  38. Fillingim RB, King CD, Ribeiro-Dasilva MC, et al. Sex, gender, and pain: a review of recent clinical and experimental findings. J Pain. 2009;10:447–85.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Nasser SA, Afify EA. Sex differences in pain and opioid mediated antinociception: modulatory role of gonadal hormones. Life Sci. 2019;237: 116926.

    Article  CAS  PubMed  Google Scholar 

  40. Wagemakers SH, van der Velden JM, Gerlich AS, et al. A systematic review of devices and techniques that objectively measure patients’ pain. Pain Physician. 2019;22:1–13.

    Google Scholar 

  41. Forslind K, Hafström I, Ahlmén M, et al. Sex: A major predictor of remission in early rheumatoid arthritis? Ann Rheum Dis. 2007;66:46–52.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  42. Yu C, Jin S, Wang Y, et al. Remission rate and predictors of remission in patients with rheumatoid arthritis under treat-to-target strategy in real-world studies: a systematic review and meta-analysis. Clin Rheumatol. 2019;38:727–38.

    Article  PubMed  Google Scholar 

  43. Maynard C, Mikuls TR, Cannon GW, et al. Sex differences in the achievement of remission and low disease activity in rheumatoid arthritis. Arthritis Care Res (Hoboken). 2020;72:326–33.

    Article  PubMed  Google Scholar 

  44. Nikiphorou E, Norton SJ, Carpenter L, et al. Remission vs low disease activity: function, quality of life and structural outcomes in the Early Rheumatoid Arthritis Study and Network. Rheumatology (Oxford). 2020;59:1272–80.

    Article  PubMed  Google Scholar 

  45. Nikiphorou E, Radner H, Chatzidionysiou K, et al. Patient global assessment in measuring disease activity in rheumatoid arthritis: a review of the literature. Arthritis Res Ther. 2016;18:251.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Ferreira RJO, Duarte C, Ndosi M, et al. Suppressing inflammation in rheumatoid arthritis: Does patient global assessment blur the target? A practice-based call for a paradigm change. Arthritis Care Res (Hoboken). 2018;70:369–78.

    Article  PubMed  Google Scholar 

  47. Ferreira RJO, Ndosi M, de Wit M, et al. Dual target strategy: a proposal to mitigate the risk of overtreatment and enhance patient satisfaction in rheumatoid arthritis. Ann Rheum Dis. 2019;78: e109.

    Article  PubMed  Google Scholar 

  48. Studenic P, Felson D, de Wit M, et al. Testing different thresholds for patient global assessment in defining remission for rheumatoid arthritis: Are the current ACR/EULAR Boolean criteria optimal? Ann Rheum Dis. 2020;79:445–52.

    Article  PubMed  Google Scholar 

  49. Doumen M, de Cock D, Pazmino S, et al. Psychosocial burden predicts sustained remission in early rheumatoid arthritis: unraveling the complex interplay of wellbeing and disease activity. Arthritis Care Res (Hoboken). 2021. https://doi.org/10.1002/acr.24847.

    Article  PubMed  Google Scholar 

  50. Pazmino S, Lovik A, Boonen A, et al. New indicator for discordance between patient-reported and traditional disease activity outcomes in patients with early rheumatoid arthritis. Rheumatology (Oxford). 2022. https://doi.org/10.1093/rheumatology/keac213.

    Article  Google Scholar 

Download references

Acknowledgements

The authors would like to acknowledge all investigators in the National Database of Rheumatic Diseases in Japan (NinJa) in 2017. We also acknowledge Akiko Komiya, who curated the NinJa database, and Satomi Hanawa, who assisted administrative work. Investigators in NinJa 2017: Kenji Ichikawa (National Hospital Organization (NHO) Hokkaido Medical Center), Yukitomo Urata (Tsugaru General Hospital), Fuminori Hirano (NHO Asahikawa Medical Center), Noriyuki Chiba (NHO Morioka Medical Center), Shintaro Yano (Maebashi Hirosegawa Clinic), Takayuki Sumida (Faculty of Medicine, University of Tsukuba), Kensuke Koyama (University of Yamanashi Hospital), Toshihide Mimura (Saitama Medical University Hospital), Yuji Akiyama (Japanese Red Cross Ogawa Hospital), Takao Sugiyama (NHO Shimoshizu National Hospital), Noboru Hagino (Teikyo University Chiba Medical Center), Toshikazu Kano (Kohnodai Hospital, National Center for Global Health and Medicine), Sakae Tanaka (The University of Tokyo Hospital), Keishi Fujio (The University of Tokyo Hospital), Hisaji Ohshima (NHO Tokyo Medical Center), Naoto Yokogawa (Tokyo Metropolitan Tama Medical Center), Tetsuji Sawada (Tokyo Medical University Hospital), Fumio Hirano (Tokyo Medical And Dental University Hospital), Akiko Aoki (Tokyo Medical University Hachioji Medical Center), Takuji Nishikawa (Tokyo Metropolitan Bokutoh Hospital), Hisanori Nakayama (Soshigaya Okura Clinic), Miho Ohshima (Tama-Hokubu Medical Center), Yasuhiko Kita (Yokohama Rosai Hospital), Atsushi Ihata (NHO Yokohama Medical Center), Satoshi Ito (Niigata Rheumatic Center), Hideshi Yamazaki (Marunouchi Hospital), Hiroshi Tsutani (NHO Awara Hospital), Isao Matsushita (University of Toyama, Faculty of Medicine), Atsushi Kaneko (NHO Nagoya Medical Center), Masao Katayama (NHO Nagoya Medical Center), Yukihiko Saeki (NHO Osaka Minami Medical Center), Koichiro Takahi (NHO Osaka Toneyama Medical Center), Kiyoshi Matsui (The Hospital of Hyogo College of Medicine), Satoshi Kashiwagi (Amagasaki Health Co-operative Hospital), Akira Okamoto (NHO Himeji Medical Center), Yasuhika Yoshinaga (Kurashiki Medical Center), Tetsu Oyama (Oyama Clinic), Eiji Sugiyama (Hiroshima University Hospital), Akinori Matsumori (NHO Kochi National Hospital), Eiichi Suematu (NHO Kyushu Medical Center), Shigeru Yoshizawa (NHO Hospital), Yojiro Kawabe (NHO Ureshino Medical Center), Kenji Kumagai (NHO Nagasaki Medical Center), Yasuo Suenaga (NHO Beppu Medical Center), Koichiro Saisho (NHO Miyakonojo Medical Center), Hirokazu Takaoka (Kumamoto Shinto General Hospital), Kenichiro Tokunaga (Japanese Red Cross Kumamoto Hospital), Hideo Otsubo (Japanese Red Cross Kagoshima Hospital), Issaku Toyohara (Northern Okinawa Medical Center)

Funding

This work was supported by the Health and Labor Sciences Research Grant from the Ministry of Health, Labor, and Welfare of Japan [19FE1003]. The funding body played no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

Authors and Affiliations

Authors

Contributions

TN and TM designed the study, interpreted data and wrote the manuscript. TN conceived the study and performed the data analysis. TM contributed to project administration and obtained funding. AH designed the study and contributed to project administration. ST contributed to the conception of the study, acquisition of data, and project administration. All authors approved the final version of the manuscript.

Corresponding author

Correspondence to Toshihiro Matsui.

Ethics declarations

Ethics approval and consent to participate

This study was conducted in accordance with the Ethical Guidelines for Medical and Health Research Involving Human Subjects. The ethics committee of National Hospital Organization Sagamihara National Hospital approved that this study is appropriate for the opt-out method according to ethical guidelines and allowed to waive the need of written informed consent. The NinJa project was reviewed and approved by the ethics committee of National Hospital Organization Sagamihara National Hospital (2014031816) and the ethics committees of participating institutions. The list of participating institutions: National Hospital Organization Tokyo National Hospital, National Hospital Organization Sagamihara National Hospital, National Hospital Organization Hokkaido Medical Center, Tsugaru General Hospital, National Hospital Organization Asahikawa Medical Center, National Hospital Organization Morioka Medical Center, Maebashi Hirosegawa Clinic, Faculty of Medicine, University of Tsukuba, University of Yamanashi Hospital, Saitama Medical University Hospital, Japanese Red Cross Ogawa Hospital, National Hospital Organization Shimoshizu National Hospital, Teikyo University Chiba Medical Center, Kohnodai Hospital, National Center for Global Health and Medicine, The University of Tokyo Hospital, National Hospital Organization Tokyo Medical Center, Tokyo Metropolitan Tama Medical Center, Tokyo Medical University Hospital, Tokyo Medical And Dental University Hospital, Tokyo Medical University Hachioji Medical Center, Tokyo Metropolitan Bokutoh Hospital, Soshigaya Okura Clinic, Tama-Hokubu Medical Center, Yokohama Rosai Hospital, National Hospital Organization Yokohama Medical Center, Niigata Rheumatic Center, Marunouchi Hospital, National Hospital Organization Awara Hospital, University of Toyama, Faculty of Medicine, National Hospital Organization Nagoya Medical Center, National Hospital Organization Osaka Minami Medical Center, National Hospital Organization Osaka Toneyama Medical Center, The Hospital of Hyogo College of Medicine, Amagasaki Health Co-operative Hospital, National Hospital Organization Himeji Medical Center, Kurashiki Medical Center, Oyama Clinic, Hiroshima University Hospital, National Hospital Organization Kochi National Hospital, National Hospital Organization Kyushu Medical Center, National Hospital Organization Fukuoka Hospital, National Hospital Organization Ureshino Medical Center, National Hospital Organization Nagasaki Medical Center, National Hospital Organization Beppu Medical Center, National Hospital Organization Miyakonojo Medical Center, Kumamoto Shinto General Hospital, Japanese Red Cross Kumamoto Hospital, Japanese Red Cross Kagoshima Hospital, Northern Okinawa Medical Center.

Consent for publication

Not applicable.

Competing interests

TN and AH have no competing interests to disclose. ST has received grants from Pfizer Japan Inc., Mitsubishi-Tanabe Pharm Co., Daiichi Sankyo Co., Ltd., Chugai Pharmaceutical Co., Ltd., and AbbVie Japan Co., Ltd. and honorariums for speakers from Pfizer Japan Inc., Ono Pharmaceutical Co., Ltd., Chugai Pharmaceutical Co., Ltd., Takeda Pharmaceutical Co., Ltd., Astellas Pharm Inc., Ayumi Pharmaceutical, QIAGEN K. K., Mitsubishi-Tanabe Pharm Co., and AsahiKASEI Co., Ltd. TM has received grants from AbbVie Japan Co., Ltd., Asahi Kasei Pharma, Astellas Pharma Inc., Chugai Pharmaceutical Co., Ltd., Eisai Co., Ltd., Pfizer Inc., Mitsubishi-Tanabe Pharma Co., and Taisho Pharmaceutical Co., Ltd. and honorarium for lectures from AbbVie Japan Co., Ltd., Asahi Kasei Pharma, Astellas Pharma Inc., Chugai Pharmaceutical Co., Ltd., Daiichi Sankyo Co., Ltd., Eisai Co., Ltd., Eli Lilly Japan K.K., Mitsubishi-Tanabe Pharma Co., Ono Pharmaceutical, Pfizer Inc., Takeda Pharmaceutical Co., Ltd., and UCB Japan Co., Ltd.

Additional information

Publisher's Note

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

Supplementary Information

Additional file 1: Methods S1.

Imputation method and models for multiple imputation by chained equations. Methods S2. Predictor matrix. Fig. S1. Cliff’s delta for sex difference in each component in available-case analysis. Fig. S2. Distributions of the share of each component. Fig. S3. Relationships between DAS28-ESR and age or joint findings before and after correction. Table S1. Patient characteristics (DMARDs free). Table S2. Patient characteristics (csDMARDs). Table S3. Patient characteristics (TNFi). Table S4. Patient characteristics (IL-6i). Table S5. Patient characteristics (CTLA-4-Ig). Table S6. Patient characteristics (JAKi). Table S7. Cliff’s delta for sex difference in disease activity indices and their components by treatment type. Table S8. PERMANOVA for sex difference in the share of components. Table S9. Adjustment for patient-related factors with GLM and QR in available-case analysis. Table S10. Adjustment for patient-related factors with GLM and QR using stacked dataset imputed by chained equations.

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

Nishino, T., Hashimoto, A., Tohma, S. et al. Comprehensive evaluation of the influence of sex differences on composite disease activity indices for rheumatoid arthritis: results from a nationwide observational cohort study. BMC Rheumatol 7, 4 (2023). https://doi.org/10.1186/s41927-023-00328-9

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s41927-023-00328-9

Keywords