Article Text

Extended report
The comparative risk of serious infections among rheumatoid arthritis patients starting or switching biological agents
  1. Jeffrey R Curtis1,2,
  2. Fenglong Xie1,
  3. Lang Chen1,
  4. John W Baddley1,
  5. Timothy Beukelman1,
  6. Kenneth G Saag1,2,
  7. Claire Spettell3,
  8. Raechele M McMahan4,
  9. Joaquim Fernandes3,
  10. Kevin Winthrop5,
  11. Elizabeth Delzell2
  1. 1Division of Clinical Immunology and Rheumatology, University of Alabama at Birmingham, Birmingham, Alabama, USA
  2. 2Department of Epidemiology, University of Alabama at Birmingham, Birmingham, Alabama, USA
  3. 3Aetna Informatics, Aetna, Blue Bell, Pennsylvania, USA
  4. 4Aetna Specialty Pharmacy, Orlando, Florida, USA
  5. 5Oregon Health and Science University, Portland, Oregon, USA
  1. Correspondence to Jeffrey R Curtis, 510 20th Street South, FOT 805D, Birmingham, AL 35294, USA; jcurtis{at}uab.edu

Abstract

Background It is unclear whether anti-tumour necrosis factor alpha and biological agents with different mechanisms of action have similar safety. This study evaluated the incidence of hospitalised infections among rheumatoid arthritis (RA) patients starting or switching various biological agents.

Methods Using a database from a large US healthcare organisation from January 2005 to August 2009, the authors identified enrollees with RA and their treatment episodes entailing the new use of a biological agent, stratified by no biological use in the previous year (‘biological-free’) or switching from a different biological agent (‘switchers’). Outcomes were hospitalised infections identified using previously validated algorithms. Proportional hazards models estimated the hazard ratio of hospitalised infections, comparing each biological agent with infliximab.

Results Among 7847 biological treatment episodes, 63% were for biological-free patients and 37% for switchers. There were 364 hospitalised infections. Rates of hospitalised infection among biological-free patients and switchers were 4.6 and 7.0 per 100 person-years, respectively (p<0.0001). After multivariable adjustment controlling for biological-free/switcher status and other infection-related factors and compared with infliximab, users of abatacept (HR 0.68, 95% CI 0.48 to 0.96), adalimumab (HR 0.52, 0.39 to 0.71), etanercept (HR 0.64, 0.49 to 0.84) and rituximab (HR 0.81, 0.55 to 1.20) had lower rates of hospitalised infection. Patient risk factors contributed more to the risk of infection than did the risk associated with specific biological therapies.

Conclusion The rate of hospitalised infections among RA patients was highest for infliximab. Most of the variability in patients' risk of infection was driven by factors other than biological agent exposure.

Statistics from Altmetric.com

Request Permissions

If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.

Numerous factors contribute to the decision to choose biological agents for the treatment of rheumatoid arthritis (RA). Appropriately powered randomised studies directly comparing biological agents are lacking, but for the average patient, the tumour necrosis factor alpha antagonist (anti-TNF) agents are thought to have efficacy that is relatively comparable to one another and perhaps to other biological therapies with different mechanisms of action.1,,4 Given efficacy that many consider relatively comparable, other considerations such as safety may be important in choosing among the biological agents for an individual patient.

Some previous reports have suggested that anti-TNF therapies may differ in their respective safety profiles.5,,8 A relatively small trial comparing abatacept and infliximab suggested safety-related differences between these two agents.9 Few other data exist with which to compare the anti-TNF biological agents with newer agents including abatacept and rituximab. One potentially confounding factor is that most patients receive these biological agents with non-tumour necrosis factor alpha mechanisms of action only after failing one or more anti-TNF therapies, consistent with some national guidelines.10 Generally, these more refractory patients are thought to be at higher risk of serious infections due to their age, comorbidities or concomitant therapies such as glucocorticoid use. However, current knowledge is inadequate to determine whether patients who switch biological therapies might be at higher risk of infections compared with those who are biological agent naive.

In light of these evidence gaps, we evaluated the drug-specific absolute and relative incidence of hospitalised bacterial infections among RA patients. We compared the individual biological agents with one another, and also compared the patients starting these agents who were free of recent biological exposure with those who were switching from one biological agent to another. Finally, because of uncertainty regarding physicians' comfort with restarting the same or a different biological agent after a hospitalised infection, we descriptively characterised whether patients were started on the same or a different biological agent in the 6 months after a hospitalised infection.

Methods

Eligible patient population and observation period

We used the administrative medical and pharmacy databases of Aetna, one of the largest health insurers in the USA that provides medical and/or pharmacy benefits to more than 18 million individuals through numerous employer sponsored and individual plans. We identified RA patients on the basis of International Classification of Disease, 9th Edition (ICD-9) codes. Patients were eligible to be under observation after they had at least two ICD-9 codes for RA (ICD-9 714.X) from a physician office visit separated by over 7 days and occurring within 1 year, or with one or more physician RA diagnosis code followed by a prescription for methotrexate or a biological agent. Patients with malignancy (excluding non-melanoma skin cancer) were excluded.

After meeting eligibility criteria, patients must also have had at least one prescription or infusion for an anti-TNF or non-anti-TNF biological agent that they had not received in the preceding 12 months, consistent with a new user design.11 The date of the biological prescription/infusion defined the ‘index date’ and began a new treatment episode. Patients must have had medical and pharmacy benefits in the 12 months before the index date and throughout follow-up. Individuals aged 65 years and older were included only if they were enrolled in the health plan through Medicare Advantage (with a pharmacy benefit). Any other type of health plan coverage (eg, commercial) for these older individuals is often a secondary, concomitant payor to governmental insurance (ie, Medicare), and enrollees' claims profile in the health plan data may be incomplete. The 12-month ‘baseline’ period was used for covariate assessment for all non-time-varying covariates.

Biological exposures

At the start of each new treatment episode, patients were assigned into one of several biological exposure groups, which included infliximab, etanercept, adalimumab, abatacept and rituximab. The use of golimumab or certolisumab was insufficient to study these agents. We classified patients as ‘biological-free’, if they had no exposure to any biological medication in the preceding 1 year baseline period, or as ‘switchers’ if they had received a different biological agent in the baseline period. Because of left censoring, some biological-free patients may have received a biological agent even before the baseline period and thus might not be truly ‘biological-naive’.

Exposure was defined as the current use of the biological agent. We specified days of exposure for injectable biological agents as the days supply plus a 90-day ‘extension’ period. The extension was added because patients who become ill often stop medications, and using an extension may capture exposure relevant to outcomes.12 For infused drugs, in which days supply is inapplicable, we assigned days of exposure as 30 days to abatacept, 180 days to rituximab and 56 days to infliximab; the same 90-day extension was added to each. As part of a sensitivity analysis, we shortened this extension to 30 days.

Outcome of interest

The outcome of interest was the first hospitalisation with at least one infection. We identified infections using a comprehensive set of ICD-9 codes developed as part of a systematic literature review and two validation studies in RA patients. The validation studies compared cases identified using these codes with cases confirmed using hospital medical records abstracted across the USA13 and with cases confirmed by abstracting records from a university hospital.14

Controlling for confounders through an infection risk score

To control for confounding, we derived and validated an infection-specific disease risk score. A disease risk score15 allows for multivariable reduction into a single composite measure, and is similar to a propensity score16 17 except that it models the risk for the outcome rather than the risk for exposure. The WHO FRAX calculator18 and several well-known cardiovascular risk calculators (eg, Framingham)19 are examples of disease risk scores. Disease risk scores achieve similar control for confounding compared with traditional multivariable adjustment, but are more efficient when outcome data are sparse. The infection risk score was used to control for infection-related confounding for all factors unrelated to biological therapy and biological-free/switcher status by computing a predicted infection score for each subject (see supplementary appendix).

Statistical analysis and covariate selection

For each treatment episode, patients were censored at the time of their first hospitalised infection, after they stopped taking the biological medication, lost medical or pharmacy benefits in the health plan (which would impede ascertaining their exposures and outcomes), started a new biological treatment episode or reached the end of the study. The study period was 1 January 2005 to 31 August 2009. The characteristics of the patients contributing biological-free and switch treatment episodes were compared. We used Cox proportional hazards models to estimate the adjusted hazard risk (HR) of infection. We used infliximab as the reference exposure group because it had sufficient data to serve as the comparator group for both the biological-free and switcher subcohorts. A sensitivity analysis used etanercept as the referent. The Cox models were stratified by biological-free versus switcher status, allowing the baseline infection hazard for these two groups to differ while still providing a combined HR. The proportional hazards assumption was verified using the SAS statement ASSESS to inspect the martingale residuals visually over the follow-up time. Control for confounding was achieved through adjustment for the infection risk score, modelled in quintiles. In a sensitivity analysis, we controlled for the infection risk score in deciles. Additional sensitivity analyses censored all observation time at 6 months and at 12 months, in order to account for the possibility that the risk of infection varied between agents over time.

In order to describe the discrimination of the statistical models in relation to the contribution of patients' infection risk score and their biological use, we calculated concordance probabilities from the Cox models.20 Concordance probabilities are similar to an area under the receiver-operator curve and allow for estimation of the incremental discrimination of adding variables to a statistical model, similar to a change in a c statistic.

Finally, we identified all biological-treated patients who developed a hospitalised infection and had at least 6 months of follow-up in the health plan after hospital discharge. Over that 6-month period, we described the cumulative proportion of people who restarted the same biological agent and those who changed to a different biological agent. Analyses were performed using SAS version 9.2. The study was approved by the university institutional review board.

Results

A total of 4916 treatment episodes occurred among biological-free patients, and an additional 2931 treatment episodes occurred in biological switchers. The mean age of the biological-free and switcher groups was 49 years, and three-quarters were women (table 1). A total of 823 patients was represented in both the biological-free cohort (17.0% of cohort) and the switcher cohort (38.9% of cohort). Fewer than 1% of patients over the age of 65 years were censored due to not being fully observable in the health plan. Compared with biological-free patients, switchers were more likely to have comorbidities such as chronic obstructive pulmonary disease and diabetes, to use narcotics and to use prednisone and at higher doses. Most of the abatacept and rituximab treatment episodes were assigned to the switcher group.

Table 1

Characteristics* of biological treatment episodes for biological-free and switch patients

Using the infection risk score, we estimated for each biological-treated patient contributing to the analysis the 1-year predicted risk of a hospitalised infection associated with components included in the infection risk score (which did not include biological agent use). As shown in figure 1, the 1-year predicted risk for hospitalisation with infection due to patient factors (and independent of biological agent use) varied across a wide range. In addition, the switcher group had a somewhat higher predicted risk of infection than the biological-free group (mean prediction infection risk score 4.2% vs 5.2%, p<0.0001).

Figure 1

Distribution of predicted infection risk among biological-free versus switchers. The x-axis is the 1-year predicted risk of hospitalised infection, expressed as a percentage ranging from 0 to 100. The y-axis is the cumulative distribution of infection risk in the biological-free and switcher groups; the area under each of the two curves sums to 100%, representing all patients in that group. The x-axis was right-truncated at 15% to facilitate visual inspection; 97% of the biological-free patients had a predicted 1-year infection risk of 15% or less; 94% of the switcher patients had a predicted 1-year infection risk of 15% or less.

After applying censoring criteria, including censoring patients once they stopped taking the biological agent, the median (IQR) follow-up was 7.7 months (4.0–14.1) after the index date; the median (IQR) follow-up was slightly greater among the biological-free patients (8.4 months (4.3–15.4), 4611 person-years of exposure) than for the switcher group (6.6 months (3.4–12.1), 2205 person-years of exposure). A total of 364 unique hospitalisations with at least one infection occurred. The site and type of infections are shown in table 2. Pneumonia, skin and soft tissue and septicemia were the most common types of infections.

Table 2

Type and number of hospitalised infections experienced by biological agent users (n=364 unique hospitalisations)*

Overall, the mean rate of hospitalised infection among biological-free patients was 4.6 per 100 patient-years and was lower than the mean rate of infection for switchers (7.0 per 100 patient-years, p<0.0001 vs biological-free patients). The rate of infection for each specific biological agent was consistently one to three per 100 person-years higher among the switcher versus the biological-free patients except for rituximab. Among rituximab users, the rate of hospitalised infection was 10.4 (6.3, 17.2) among biological-free patients and 7.1 (4.6, 10.9) among switchers.

The main results of the study are described in figure 2 and table 3. Compared with infliximab and after adjusting for the infection risk score, the rate of infection was significantly lower for most other drugs (figure 2) in both the biological-free and switcher subcohorts. The combined HR shown in figure 2 that pooled results for the biological-free and switcher cohorts are shown numerically in table 3. The differences in the HR yielded much smaller differences in the absolute rate of infection for low-risk patients compared with high-risk patients. For example, among low-risk biological-free patients and comparing infliximab with adalimumab users (combined HR 0.52), there was a −2.1 per 100 person-years lower rate of infections among adalimumab users for low-risk patients and a −8.2 per 100 patient-years difference for high-risk patients. The quintile of the predicted infection risk score had a strong dose-dependent relationship with the rate of infection (bottom half of table 3). Indeed, the HR associated with the infection risk score were much greater in magnitude than the differences between the biological agents.

Figure 2

Adjusted* HR of hospitalised infection for biological-free, switcher and combined cohorts, compared with infliximab. *Adjusted for quintile of the infection risk score, described in the supplementary appendix (available online only).

Table 3

Multivariable adjusted* rates of hospitalised infection and associated rate differences for biological-free low and high-risk patient subgroups

The results of the sensitivity analyses for table 3 that defined current use with a 30-day extension period and of the analysis that controlled for the infection risk score in deciles were similar to the primary results (data not shown). Additional sensitivity analyses that censored patients at 6 months and at 12 months yielded similar magnitude HR and statistical significance to the primary results. In the sensitivity analysis that used etanercept rather than infliximab as the referent group in the adjusted Cox models, there was no significant difference in the rate of infection between etanercept and any of the other biological agents except for infliximab.

The concordance probability for infection including only the predicted infection risk score and biological-free/switcher status in the model was 0.67; after adding the specific biological drugs to the model, it increased minimally to 0.68. In contrast, the concordance probability with only the specific biological drugs and biological-free/switcher status in the model was 0.59.

Among the 240 patients with a hospitalised infection who had at least 6 months of follow-up, a majority of patients (69.2%) were re-started on the same biological agent within 6 months and did not switch biological agents. An additional 6.7% of patients restarted the same biological agent but switched within 6 months to a different agent. Ten per cent of patients switched to a different biological agent within 6 months after hospital discharge, and the remaining 14.2% of patients did not use any biological agent. Data were insufficient to characterise the risk for subsequent hospitalised infection after re-starting or changing biological therapy.

Discussion

In this analysis of RA patients enrolled in a large US health plan, the rate of hospitalised infections for specific biological agents varied across a relatively small range. Patients switching biological agents had an increased rate of infection in the range of one to three per 100 patient-years compared with those not receiving biologic agents in the previous 12 months. The rate of hospitalised infections for patients treated with infliximab was higher than for other agents. However, to put these findings into context, patient characteristics that were part of the infection risk score contributed to much wider variability in the rate of serious infections than the risk associated with specific biological drugs.

A substantial body of literature supports a modest increase in the risk of serious infection for patients treated with anti-TNF agents compared with non-biological disease-modifying antirheumatic drugs (DMARD).21,,23 The rate of infections for biological-treated patients compared with those receiving non-biological DMARD reported in past studies is typically a difference of one to three per 100 person-years. We found rate differences similar to this magnitude for patients switching to a new biological agent. A higher infection risk among infliximab users compared with some other anti-TNF therapies has previously been observed(5), although may be most different for opportunistic infections (6–8). Consistent with earlier findings,24 pneumonia and soft tissue infections were the most common infections observed in our cohort.

The clinical importance of RR in the range of approximately 0.5–0.8, as we found in table 3, is a matter of judgement. For patients at low risk of infection, this translates into a relatively small rate difference of approximately two per 100 patient-years. For patients at higher baseline infection risk, the magnitude of the rate difference is greater, in the range of five to eight per 100 patient-years. An earlier survey of the membership of the American College of Rheumatology suggested that rate differences less than six to eight per 100 person-years may not be very important to clinical decision-making.25 The magnitude of the incremental risk related to biological treatments compared with patients' baseline infection risk may impact both patients and physicians' perceptions of the risk of starting biological agents and should be examined prospectively in relation to treatment preferences for RA medications.

Our results provide useful information to put differences in the safety profiles of biological agents into context. The contribution of demographics, comorbidities and glucocorticoid use on infection risk (summarised as a single infection risk score) contributed substantially greater variability to infection risk compared with smaller differences in infection risk between individual biological agents. Similarly, the improvement in the discrimination of the infection risk model was substantial after adding the infection risk score, and was minimal after adding the biological agents. While some of these infection risk factors are not modifiable (eg, age), others might be to a variable extent (eg, glucocorticoid use, some comorbidities), and may offset any increased infection risk associated with use of specific biologics.

The strengths of our study include a large population of patients receiving care in a real-world setting. Unlike clinical trials and even some registries, patients were included in these safety-focused analyses irrespective of age or concomitant comorbidities. Despite these strengths, the results of our study must be interpreted in light of its design. Administrative data provide limited information on RA disease severity and comorbidity. This yields potential residual confounding. We also lacked hospital medical records to confirm the infections, although our claims-based algorithms have been validated in two large studies.13 14 We observed a modest number of outcome events for some exposures, especially for rituximab, which led to CI that were wider than for other agents. We also observed higher rates of infection for biological-free rituximab users, which probably reflects channelling of higher risk patients to this agent. Although we had a somewhat limited ability to examine risk in specific time windows after biological agent initiation,5 12 most of our follow-up time occurred within the first year after initiation, and sensitivity analyses that censored observation time at 6 or 12 months supported our primary findings. Finally, enrollees in the health plan were relatively young and healthy, and these results may not be generalisable to older patients or those with a higher burden of comorbidities.

In conclusion, we found wide variability in patients' risk of serious infections related to demographics, comorbidities and other patient-specific factors, and much less variability related to the individual biological agent or switching patterns. A majority of patients were continued on the same biological therapy despite having been hospitalised with an infection, despite little information on the safety of this practice. More attention needs to be focused on patient-specific, potentially modifiable risk factors for infection such as glucocorticoid use and immunisations. The ultimate importance of perceived safety risks associated with biological agents in relation to patients' other infection risk factors, and the impact of these perceptions on treatment selection, may be useful to examine in future studies.

Appendix

Derivation and performance of infection risk score in original and bootstrap datasets

To derive and validate the infection risk score, we identified rheumatoid arthritis (RA) patients who were new users of methotrexate (MTX), leflunomide (LEF), sulfasalazine (SSZ) or hydroxychloroquine (HCQ) that were not exposed to any biologic in a subsequent 1-year period. The SSZ and HCQ patients must also have been on MTX in the preceding year so as to identify patients with greater RA severity. The outcome was hospitalised infection in the first year of follow-up (n=202 hospitalised infections used to develop the infection risk score). Logistic regression was used to derive the infection risk score equation by estimating the coefficients for each of the factors included in the score. Factors in the model were selected based upon clinical interest and review of the literature.13 26 These factors included age and sex; extra-articular RA (Felty's syndrome or other systemic involvement); comorbidity including prior outpatient infection, prior hospitalised infection, COPD, dementia, peptic ulcer disease, diabetes, obesity, alcoholism, liver disease, ischemic cardiovascular disease, CIRAS,27 the Charlson Comorbidity Index; health services utilisation including number of outpatient physician visits, number of hospitalisations for conditions other than infections, number of liver function lab tests, any rheumatoid factor lab test, mammography; and medications including steroid use/dose (none, 0 to <5 mg/day, 5 to <10 mg/day, ≥10 mg/day) and narcotics.

Because any prediction model is subject to overfitting, we used bootstrapping to assess performance of the infection risk score, randomly selecting 200 bootstrap samples of equal size to the original dataset. Discrimination and calibration were averaged across the bootstrap samples (figure).

The c statistic (area under the receiver operator curve) was 0.74 for the original data; the median c statistic was 0.74 in the bootstrap data. Error bars indicate the SD of the infection risk score. Note: All factors were measured in the 1-year baseline period before starting MTX/LEF/SSZ/HCQ.

References

Supplementary materials

Footnotes

  • Funding This work was supported by the Agency for Healthcare Research and Quality (R01HS018517) and the Doris Duke Charitable Foundation. Some of the investigators receive support from the National Institutes of Health (AR053351).

  • Competing interests JRC: (research grants and consulting) UCB, Roche/Genentech, CORRONA, Centocor, Amgen, BMS, Abbott. ED: (research grants) Amgen. Remainder of co-authors: none.

  • Ethics approval This study was conducted with the approval of the University of Alabama at Birmingham institutional review board.

  • Provenance and peer review Not commissioned; externally peer reviewed.