Comparative risk of infections between JAK inhibitors versus TNF inhibitors among patients with rheumatoid arthritis: a cohort study

Abstract

Background

To compare infectious risk between JAK inhibitors (JAKis) versus TNF inhibitors (TNFis) among rheumatoid arthritis (RA) patients in Korea.

Methods

Using 2009–2019 Korea National Health Insurance Service database, we conducted a cohort study on RA patients initiating a JAKi or TNFi. The primary outcomes were herpes zoster (HZ), serious bacterial (SBI), and opportunistic infections (OI). Propensity-score fine-stratification (PSS) and weighting were applied to adjust for > 70 baseline covariates. Hazard ratios (HRs) and 95% confidence intervals (CIs) were estimated using Cox proportional hazard models comparing JAKi versus TNFi users.

Results

We included 2963 JAKi initiators PSS-weighted on 5169 TNFi initiators. During a follow-up of 1.16 years, the most frequent type of infections was HZ with incidence rate (IR) per 100 person-years of 11.54 and 4.88 in JAKi and TNFi users, respectively. The IR of SBI was 1.39 and 1.32, respectively. The OI was rare with a majority being tuberculosis and showed an IR of 0.11 and 0.49 in JAKi and TNFi users, respectively. The PSS-weighted HR (95% CI) for individual types of infections was 2.37 (2.00–2.80) for HZ, 1.04 (0.71–1.52) for SBI, and 0.25 (0.09–0.73) for OI.

Conclusions

This population-based cohort study on RA patients treated with JAKi or TNFi in Korea showed an exceptionally high IR of HZ in both treatment groups compared to that from Western countries, with an approximately doubled risk associated with JAKi versus TNFi use. The risk of SBI was comparable, but the risk of OI, particularly tuberculosis, was less among JAKi than TNFi initiators.

Introduction

Rheumatoid arthritis (RA) is a systemic inflammatory disease in which synovial joints are the primary target of autoimmunity [1]. However, chronic inflammation causes not only the joint failure but also a wide spectrum of comorbidities [2]. Therefore, international guidelines endorse the treat-to-target strategy to achieve remission or low disease activity [3, 4]. With the introduction of biologic and targeted synthetic disease-modifying anti-rheumatic drugs (bDMARDs and tsDMARDs, respectively), clinical outcome of RA refractory to conventional DMARDs (cDMARDs) has dramatically improved [5]. In particular, Janus kinase inhibitors (JAKis) of the tsDMARD class have shown impressive efficacy against RA [6]. According to the international guidelines, the b/tsDMARDs are used as monotherapy or in combination with DMARDs including methotrexate (MTX) to treat moderate-to-severe RA [3, 4]

Infection is one of the most common treatment-emergent adverse events in RA patients due to disease-associated immune alteration and/or treatment-related immune suppression. Overall, there is a twofold risk of serious infections among RA patients compared to the non-RA population [7]. Also, there has been a particular concern for infection among RA patients treated with high efficacy DMARDs including bDMARDs or JAKis [8,9,10]. The incidence rate (IR) of serious infections in randomized controlled trials (RCTs) investigating tofacitinib, a JAKi, was similar to that in RCTs evaluating bDMARDs including TNF inhibitors (TNFis) in patients with RA [10]. However, unlike bacterial infections, the risk of herpes zoster (HZ) with tofacitinib was significantly higher than that with TNFis [11, 12]. In particular, higher susceptibility to developing HZ in Asian patients has been suggested [12].

Despite such backgrounds, population-based studies have been few in the real-world setting that directly compared the risk of infections of JAKis versus bDMARDs users among RA patients of Asian ancestry [13]. To meet this end, we compared the risk of HZ, serious bacterial infections (SBI), and opportunistic infections (OI) among RA patients treated with JAKis versus TNFis using the nationally representative Korea National Health Insurance Service (KNHIS) database.

Methods

Data source

We used the 2009–2019 KNHIS database. The KNHIS database contains longitudinal patient data including demographics, International Classification of Diseases Tenth Revision (ICD10) diagnosis codes, procedures, prescription records (drug names, prescription and dispensing dates, days’ supply, dose, and route of administration), and type of medical utilization (outpatient, inpatient, or emergency department) of all Korean citizens [14]. The Institutional Review Board of the Seoul National University Bundang Hospital approved the study protocol (X-2207–770-901) and waived the need for written patient consent based on de-identified database. This study was conducted in accordance with the principles of the Declaration of Helsinki and Good Clinical Practice Guidelines.

Study population

RA patients with at least two ICD10 codes of RA and any DMARD (MTX, leflunomide, hydroxychloroquine, sulfasalazine, tacrolimus, cyclosporine, mizoribine, azathioprine, infliximab, adalimumab, golimumab, etanercept, abatacept, tocilizumab, rituximab, tofacitinib, or baricitinib; certolizumab not available in Korea) were eligible [15]. We also applied a V-code (V223) that indicates a copayment beneficiary (90% reduction of the cost) associated with seropositive RA in Korea. Thus, all included study participants were eventually seropositive for rheumatoid factor and/or anti-citrullinated protein antibody.

Among the above patients, we first selected those who initiated a JAKi (tofacitinib or baricitinib) or TNFi (infliximab, adalimumab, golimumab, or etanercept) with the first dispensing date of the corresponding index drug being defined as the index date. In Korea, bDMARD was first introduced in the market in 2004 to treat RA refractory to cDMARDs, and JAKi was first approved in 2015 as a 3rd line treatment after bDMARD failure and then approved as a 2nd line treatment after cDMARD failure from 2017. The approved dose of tofacitinib and baricitinib in Korea to treat RA does not exceed twice daily dose of 5 mg and once daily dose of 4 mg, respectively. We only included new users of individual study drugs by excluding patients who had prior use of the given study drugs during the 365-day pre-index period (= baseline period). A significant proportion of JAKi initiators used bDMARDs during the baseline period. Therefore, we allowed bDMARD use in both treatment groups as long as the bDMARDs used at baseline were not the index drug in order to utilize as many as JAKi initiators without losing comparability between the two groups (Fig. 1). But we required the TNFi group be free of JAKi during the baseline period. Other exclusion criteria included those with dialysis services or human immunodeficiency virus infection at baseline, and those hospitalized due to serious infection within 30 days prior to the index date to avoid re-hospitalization from the same episode of infection or residual effect of previous serious infection on patients’ general health.

Fig. 1

Study cohort selection process. HIV, human immunodeficiency virus; ICD, International Classification of Diseases Tenth Revision; JAKi, JAK inhibitor; KNHIS, Korea National Health Insurance Service; PSS, propensity score fine stratification; RA, rheumatoid arthritis

Outcomes

The primary outcomes were HZ, SBI, and OI. HZ cases were identified with either inpatient primary diagnosis codes or outpatient diagnosis codes plus use of anti-viral medications (i.e., acyclovir, valacyclovir, or famciclovir) within ± 7 days of the diagnosis for HZ [16]. SBI included meningitis, encephalitis, cellulitis, endocarditis/myocarditis, pneumonia, pyelonephritis, septic arthritis, osteomyelitis, and septicemia/bacteremia, and OI included tuberculosis, non-tuberculous mycobacterial infection, and systemic fungal infections (cryptococcosis, or aspergillosis). Such algorithm for SBI and OI has shown a positive predictive value of 80% using their inpatient diagnosis codes in the primary position [17].

Secondary outcomes included (1) serious HZ infection defined as the hospitalized cases where inpatient diagnosis was HZ in the primary position and (2) tuberculosis among our OI cases.

Covariates

During the 365-day pre-index baseline period, we measured > 70 variables including patients’ demographics, index calendar years, RA medications including the number and class of bDMARD used, non-RA medications, potential risk factors for infection, use of anti-microbials, and markers of health care utilizations (listed in Table 1). We also estimated a Charlson-Deyo score for multi-morbidities [18].

Table 1 Baseline characteristics of study participants

Statistical analysis

In our primary as-treated analysis, patients were followed from the day after the index date to the first occurrence of the following events: outcome occurrence, disenrollment, death, discontinuation of the index treatment, or adding any other DMARDs over the index treatment. Switching between different TNFis or between different JAKis was not a censoring event. Drug discontinuation was defined as no dispensing within 90 days from the expected refill date. Patients who discontinued the study medication were followed up until the last expected refill date plus 30-day grace period. The expected refill date was calculated by adding days’ supply to the last dispensing date of the study medication. The days’ supply of individual TNFi was 56 days for infliximab, 14 days for adalimumab, 28 days for golimumab, and 7 days for etanercept. For the secondary intention-to-treat (ITT) analysis, we followed patients up to 365 days after the index date without censoring on drug switching, adding, or discontinuation.

For confounding adjustment, we used propensity score (PS) fine stratification and weighting to account for > 70 baseline covariates listed in Table 1 [19]. A multivariable logistic regression model incorporated all of these covariates including the index year to estimate a PS, which was defined as the predicted probability of a patient initiating a JAKi versus TNFi given aforementioned baseline covariates. After trimming patients in the non-overlapping areas of PS, we created 50 strata based on the distribution of PS among the exposed (i.e., JAKi treatment); patients in the comparator therapy group were weighted proportionally to the distribution of the JAKi group within each of the 50 PS strata. The covariate balance between the two groups among the PS stratification (PSS) weighted study cohort was evaluated by standardized mean differences: a balanced covariate distribution was considered achieved with a standardized mean difference of < 0.1 [20]. PSS-weighted IRs of primary and secondary outcomes were calculated per 100-person-years. We used a Cox proportional hazard model estimating the hazard ratio (HR) and 95% confidence intervals (CIs). Proportional hazard assumptions were not violated in any of the models when tested using the interaction term between the exposure and follow-up time [21]. All analyses were performed using the SAS 9.4 (SAS Institute) software.

Subgroup analysis

To investigate the risk factors associated with infections of interest, we performed analyses on subgroups stratified by age (≥ and < 60 years), concurrent use of MTX, and concurrent use of steroids. The PSS-weighting was separately done for each subgroup analysis. The interaction between treatment and individual stratifying factors was also tested using the Cox model.

Results

Baseline patient characteristics

We identified 2967 JAKi initiators and 13,152 TNFi initiators from the database and generated 2963 JAKi initiators (78.8% tofacitinib, 21.2% baricitinib) PSS-weighted on 5169 TNFi initiators (Fig. 1). The baseline patient characteristics before and after PSS-weighting were summarized in Table 1.

Before PSS-weighting, bDMARD use (40.2 vs 17.3%) prior to the index date was more common among JAKi than TNFi initiators, and the average number of bDMARDs used (0.5 vs 0.2) was also higher in JAKi initiators. This is expected since JAKi had been initially approved as a 3rd-line treatment after bDMARD failure in Korea and later approved as a 2nd-line treatment. TNFi than JAKi initiators showed a higher mean cumulative steroid dose (1641 vs 1437 mg of prednisone equivalent dose) and more frequently used analgesics (both non-steroidal anti-inflammatory drugs and opioids) at baseline. The comorbidity profile was comparable in general between the two groups. In particular, 7–8% of patients already experienced anti-HZ treatment before the index date. Healthcare service use was more common among TNFi initiators than JAKi: 35.8 vs. 30.9% of hospitalization, 18.3 vs. 15.1% of emergency department visits, 20.2 vs.12.0% of HbA1c test ordered, 92.0 vs 53.8% of serum creatinine test ordered, and 86.9 vs 50.6% of lipid profile test ordered.

After PSS-weighting, all of the baseline covariates were well-balanced according to a standardized difference of < 0.1. MTX use at baseline was observed in 87.1% of both treatment groups. Although MTX dose was not included in the PS-estimating logistic model, we observed a well-balanced distribution between the two groups (standardized differences of < 0.1) of mean index and maximal doses of MTX (standard deviation) during follow-up: the index dose of 9.5 (6.9) mg and maximal dose of 10.2 (6.6) mg among JAKi initiators and 9.6 (7.7) mg and 10.7 (7.2) mg among TNFi initiators.

Comparative risk of infections between JAK inhibitor and TNF inhibitor users

The most frequent type of infections was HZ among others (Table 2). During a mean follow-up of 1.16 years, 582 cases of HZ occurred with the IR of HZ per 100 person-years of 11.54 and 4.88 in JAKi and TNFi users, respectively. The PSS-weighted HR (95% CI) of HZ was 2.37 (2.00–2.80) comparing JAKi and TNFi users. Among all HZ cases in JAKi users, 17.2% were serious, requiring hospitalizations. The risk of such serious HZ infection was even higher among JAKi users than TNFi with the PSS-weighted HR (95% CI) of 7.43 (3.91–14.11).

Table 2 Infectious risk comparing JAK inhibitor and TNF inhibitor users

The IR of SBI was similar between the two groups: 1.39 and 1.32 per 100 person-years in JAKi and TNFi users, respectively, with the PSS-weighted HR (95% CI) of 1.04 (0.71–1.52). Upper urinary tract infection (47.7% of all SBI) was most common, followed by septic arthritis (22.0%), cellulitis (21.1%), and pneumonia (9.2%) (Supplemental Table 1).

The OI was rare with a majority (22/27, 81.5%) being tuberculosis. Four OI cases among JAKi users were tuberculosis (n = 2) and systemic fungal infections (n = 1 for systemic candidiasis, n = 1 for aspergillosis) while 20 out of 23 OI cases among TNFi users were tuberculosis with the rest being non-tuberculous mycobacterial infection. The PSS-weighted HR (95% CI) for OI was 0.25 (0.09–0.73) comparing JAKi versus TNFi users. The ITT results were similar.

Subgroup analyses

We performed subgroups analyses stratified by age (≥ and < 60 years), concurrent use of MTX, and concurrent use of steroids (Supplemental Table 2). The IR of HZ, SBI, and OI was substantially higher among those aged ≥ than < 60 years. However, the combined use of MTX or steroid did not alter the IR of these infections. Although IRs of individual types of infections within each subgroup were affected by stratifying factors, the PSS-weighted HR between JAKi versus TNFi users was similar to that of the main PSS-weighted cohort (Table 3). There was no significance when interaction between treatment and stratifying factors was tested.

Table 3 Comparative infection risk in subgroups

Discussion

In this population-based cohort study, we confirmed that the IR of HZ is exceptionally high among patients with RA in Korea compared to patients from Western countries [12, 22, 23]. The risk of HZ was 2.37-fold higher in JAKi users than TNFi in our study but the IR of HZ was substantially high even in our TNFi users. The risk of SBI was comparable between the two treatments. The risk of opportunistic infection, with the majority of which cases being tuberculosis, was lower among JAKi users than TNFi.

An approximately doubled risk of HZ among JAKi users than TNFi has been shown in the ORAL Surveillance study or US claims/registry studies [22,23,24,25]. Despite the similar comparative risk of HZ between the two treatments, the IR of HZ in our study was exceptionally high among both JAKi and TNFi users, compared to that from the Western countries [26, 27]: in a study by Winthrop et al. [27] on 6192 RA patients treated with tofacitinib for a median of 3.4 years from the global clinical trials, the IR of HZ was 3.9 (95% CI 3.6 to 4.2). Yet, the IRs of HZ varied across regions, from 2.4 (95% CI 2.0–2.9) in Eastern Europe to 8.0 (95% CI 6.6–9.6) in Japan and 8.4 (95% CI 6.4–10.9) in Korea. In our real-world setting, the IR of HZ exerted by JAKi was far higher (11.54, 95% CI 10.47–12.72). Moreover, the IR of TNFi users in our study was as high as 4.88 (95% CI 4.29–5.55), also higher than the IR of patients from the Western countries treated with a JAKi [26, 27]. Surprisingly, the risk of HZ was not elevated in the Chinese or Taiwanese compared to that from a pooled analysis on the global trials [28, 29]. This finding suggests that heterogeneity exists even within Asian ethnicities. Of note, the IR (95% CI) of serious or hospitalized HZ was 1.81 (1.42–2.31) among JAKi users, which is 7.43 times than among TNFi users (IR 0.24, 95% CI 0.13–0.42).

According to the ORAL Surveillance trial [22] and the US-based registry study [25], there was no difference in the risk of serious infections comparing JAKi versus TNFi or bDMARD users. Similarly, we did not find any difference in the risk of SBI between JAKi users and TNFi. Also, pneumonia, cellulitis, and urinary tract infections were most common types of SBI as in other studies [22, 25].

The IR of OI was low as in previous RCTs. The pooled analysis of the global clinical trials and long-term extension studies showed that the IR per 100 person-years was 0.3 for OI other than tuberculosis and 0.2 for tuberculosis [26]. In our study, the IR of OI including tuberculosis was 0.11–0.49 per 100 person-years, with the majority of OI being tuberculosis. Because the case definition of OI in our study required hospitalization, our findings suggest that TNFi might confer a higher risk of severe tuberculosis than JAKi. However, this interpretation needs cautions due to the rarity of cases and further research is needed with a larger cohort of patients.

Old age, Asian ethnicity, and/or steroid use were found to be baseline risk factors for HZ in previous studies [13, 24, 27]. We found in our subgroup analyses that old age increased the risk of HZ not only among JAKi but also in TNFi users and that it was a risk factor for non-HZ infections as well. However, we did not find any significant increase in the IR of infections associated with concomitant use of steroid unlike studies on the Western populations [24, 27]. Such lack of association is similar to the Japanese multicenter study [13]. Concomitant use of MTX was not associated with a higher risk of infection, as in previous studies [13, 24, 27]. Overall, the comparative risk of infections associated with JAKi versus TNFi use was consistent across in all subgroups.

The strengths of this study are as follows: first, we used rigorous pharmacoepidemiologic methods to reduce confounding between comparator groups. The new user design with active comparator is a powerful tool to cope with both measured and unmeasured confounding [30]. In addition to that, we used PSS-based weighting to further adjust for > 70 covariates at baseline. Second, this is one of the few studies that provide a head-to-head comparison on infections between JAKis and TNFis in the real-world setting, and to the best of our knowledge, the first population-based study on Asians. Third, we used a nationally representative data to ensure high generalizability. Fourth, we used validated algorithms to define outcomes and investigated different types of infections ranging from viral, bacterial to opportunistic [16, 17]. Fourth, we performed relevant subgroup analyses to identify high risk subsets.

There are also limitations. First, inherent to any observational studies, there is a concern for residual or unmeasured confounding particularly due to lack of direct information on RA duration or activity at baseline. Nonetheless, to minimize such limitation, we used the active comparator design and further accounted for many proxies for RA activity such as the number of DMARD used, individual DMARDs used, and steroids use and their cumulative dose. Second, the number of patients receiving JAKi was small even in this nation-wide database, leading to limited power for rare outcomes. The incidence of OI was small, and more research is needed to study the comparative risk of OI to compare JAKi and TNFi. Also, the small number of outcomes in certain subgroup analyses resulted in lack of precision and wide 95% CIs. Third, the KNHIS data do not provide information on vaccination status for herpes zoster since vaccines are dispensed at individuals’ own expense.

Conclusion

Taken together, despite the similar comparative risk of HZ among JAKi users versus TNFi to that in other ethnic groups, the IR of HZ was far higher among RA patients in Korea. We were able to confirm that the risk of SBI was very similar between the two treatments. We also cautiously suggest that TNFi might confer a higher risk for tuberculosis that requires hospitalizations. When we performed a subgroup analysis to identify risk factors associated with infections in JAKi and TNFi users, we found that the elderly was more susceptible for all types of infections but concomitant use of either steroid or MTX did not increase the risk of infections.