Objective To examine the risk of serious infection conveyed by tumour necrosis factor α (TNFα) inhibitors in the treatment of rheumatoid arthritis (RA).
Methods Data from patients with RA enrolled in the German biologics register RABBIT were used for analysis. Baseline patient characteristics, time-varying risk factors (treatment changes, functional capacity) and selection processes caused by dropout, death or switching to non-anti-TNF treatment were taken into account to estimate the adjusted incidence rate ratios (IRRadj) of serious infection during treatment with TNF inhibitors compared with non-biological disease-modifying antirheumatic drug treatment.
Results Data were available on 5044 patients, in whom 392 serious infections occurred. The crude rates of serious infections in patients treated with TNF inhibitors declined over the first 3 years of observation (from 4.8 to 2.2/100 patient-years). This decline was driven by (1) treatment termination or loss to follow-up in patients at increased risk and (2) a risk reduction through decreasing glucocorticoid doses and improvement in function. Adjusted for selection processes and time-varying risk factors, the following parameters assessed at baseline (age, chronic diseases) or at follow-up prior to the infection were significantly associated with an increased risk: age >60 years, chronic lung or renal disease, low functional capacity, history of serious infections, treatment with glucocorticoids (7.5–14 mg/day, IRRadj 2.1 (95% CI 1.4 to 3.2); ≥15 mg/day, IRRadj 4.7 (95% CI 2.4 to 9.4)) and treatment with TNFα inhibitors (IRRadj 1.8 (95% CI 1.2 to 2.7)).
Conclusion Reasons for the decline in infection rates observed at the group level were identified. The results enable expected infection rates to be calculated in individual patients based on their risk profiles.
This paper is freely available online under the BMJ Journals unlocked scheme, see http://ard.bmj.com/info/unlocked.dtl
Statistics from Altmetric.com
A decade ago tumour necrosis factor (TNF) inhibitors were approved for the treatment of rheumatoid arthritis (RA). Uncertainties regarding the specific risks of these new agents led to intensified efforts to investigate their safety. Meta-analyses of randomised controlled trials (RCTs) were undertaken, and various data sources including claims databases were used to assess the risk of serious adverse events (eg, serious infections) possibly associated with these agents.1,–,11 One major innovation was the establishment of biologics registers to evaluate the safety and effectiveness of these drugs under the conditions of daily rheumatological care.12 Upon analysis of these data, it became increasingly clear that the question of how a drug affects the risk of a specific adverse event is far more complex than evaluating the efficacy of the drug. Some of the methodological difficulties we faced while analysing real-world data will be described in this paper using the example of examining how anti-TNF therapy affects the risk of serious infection.
An increased rate of serious infections with anti-TNF therapy compared with conventional disease-modifying antirheumatic drug (DMARD) therapy was found in a meta-analysis of RCTs,1 in an analysis of claims data8 and in observational data.3 6 9 Other findings were in contrast to these results2 10 11 or reported a decline in the infection risk over time in patients treated with TNF inhibitors,5 7 8 13 and therefore raised the question whether or not the risk is increased only during the first months of treatment.5 7 13
The first aim of our study therefore was to determine whether or not there is a methodological or clinical explanation for this decline in risk and, if so, what the relative contributions might be.
Our second question was how the risk reduction seen in the cohorts is reflected at the level of the individual patient. Patient demographics, clinical features and follow-up information such as treatment response and patient use of additional medications were studied to calculate the expected incidence rates of serious infections for defined subgroups of patients. The aim of this approach was to enable the treating physician to assess the magnitude of infection risk that he or she imposes on a patient when making specific treatment decisions.
The analysis was based on patients with RA enrolled in the German biologics register RABBIT, an ongoing prospective cohort study, at the start of treatment with a biological agent or a conventional DMARD, between 1 May 2001 and 31 December 2006. Patients were followed up independent of any change in their treatment regimes. Observation time following a start of treatment with a non-anti-TNF biological agent was excluded. Patients treated with anakinra at baseline (n=89) and those for whom only baseline data were available (n=141) were excluded.
At baseline and at predefined points of follow-up (at 3, 6 months and thereafter every 6 months), rheumatologists assessed the clinical status of the patient including the components of the disease activity score based on 28-joint counts (DAS28),14 reported treatment details and adverse events. Patients assessed, among other items, their functional capacity as a percentage of full function by means of the Hannover Functional Status Questionnaire (Funktionsfragebogen Hannover (FFbH)).15 16 Reasons for dropout from the study and causes of death were ascertained by contacting health authorities. Follow-up data (including serious infections) reported prior to 1 November 2009 were included. An infection was attributed to anti-TNF treatment when the treatment was either ongoing or terminated ≤3 months prior to infection. A patient who switched from anti-TNF treatment to treatment with non-biological DMARDs contributed to anti-TNF exposure time until 3 months after switching and to DMARD exposure time thereafter. A DMARD-treated patient who switched to anti-TNF treatment also contributed time to both exposure groups (for further details see Strangfeld et al17 or Listing et al3).
Three different statistical models developed in a complementary fashion were used. In model A we adjusted the data only for differences between the anti-TNF and DMARD groups at the start of treatment (confounding by indication). A propensity score, as described earlier,18 was used for adjustment. Further adjustment was made for risk factors for infection (table 3). In model B we also considered changes in glucocorticoid dosage and changes in functional status during follow-up as time-varying covariables. The differences between models A and B in the estimates of the trend parameter were used to calculate the relative contribution of clinical improvement to the decreased risk of infection. In the fully adjusted model C, generalised estimation equations (GEEs) were applied to a weighted sample of patients. In addition to model B, this model takes into account confounding by anti-TNF treatment termination, the impact of dropout and previous serious infections.19 We defined 'dropout' as either loss to follow-up, death or switching to non-anti-TNF therapy.
For each patient and each 6-month follow-up period we calculated the probability of being in a certain status. Using logistic regression we estimated the probabilities of being a dropout with and without accounting for gender and patient characteristics at follow-up (occurrence of serious infections, DAS28, functional status). We calculated weights as ratios of both probabilities. The weights for different time periods were combined as described in detail by Molenberghs et al.20 Again, for each 6-month period the likelihood of being treated with anti-TNF agents was calculated by taking and not taking into account covariables (age, number of DMARD failures and the patient characteristics listed above). This resulted in weights which account for treatment decisions.19 20 Both types of weights (for dropout and treatment groups) were uncorrelated (|r|<0.1) and were then combined by multiplication to a final weight. By using this inverse probability weighting we obtained balanced samples of patients treated with or without TNFα inhibitors and of patients who continued or discontinued follow-up.
The application of multivariate GEE models requires a sufficient number of observed serious infections. We aimed to detect a ≥1.5-fold increase in infection risk caused by TNF inhibitors or a trend in risk of similar size (≤0.67=1/1.5) that would be comparable to other observations.2 The number of serious infections observed in the first and second years would therefore ensure sufficient power (80%) for this analysis, whereas the total number of infections observed in the third year (n=83) was too low to achieve sufficient power. Because multiple infections in individual patients are not independent of each other, the GEE models require an estimate of intrapatient correlations. The number of these correlations increased considerably when the data were analysed for all 3 years. Therefore, a robust estimate of the 3-year infection risk was not possible. For both reasons (the small number of events and the high number of intracorrelations), we restricted the multivariate analysis to the first 2 years.
Owing to strict monitoring, information about DMARD or anti-TNF treatment exposure was complete in >99% of all patients. For all follow-up time points within the first 2 years (in dropouts before their last visit), 5.6% of DAS28 scores, 3.7% of FFbH scores and 3.6% of glucocorticoid dosages were missing. These values were replaced using statistical imputation methods based on the expectation maximation algorithm.20 SAS software version 9.2 (PROC GLIMMIX, PROC MI) was used for the computations. p Values <0.05 were considered to be statistically significant.
Patient characteristics and treatment status at study entry
Between 1 May 2001 and 31 December 2006, 5274 patients were enrolled in RABBIT of whom 5044 met the inclusion criteria for this study. The average follow-up time was 2.6 years. Upon enrolment, patients treated with TNFα inhibitors differed significantly from those receiving conventional DMARDs with regard to age, disease duration, disease activity (DAS28), functional capacity (FFbH) or number of failed previous DMARD treatments (table 1, columns 2–5). The proportion of patients receiving higher doses of glucocorticoids was also significantly different. Inverse probability weighting was used for the final analysis and allowed for balanced anti-TNF and DMARD samples (table 1, columns 6–8).
Trends in crude infection rates
The crude rates of serious infections decreased significantly from the first to the second and third years in patients treated with TNF inhibitors, whereas these rates remained stable in the cohort of DMARD-treated patients (table 2). The decrease in the absolute rates of infection in patients treated with TNF inhibitors led to a decrease in the relative rates or incidence rate ratios (IRRs).
Changes in patient risk profiles over time
Overall, 1893 patients dropped out because of loss to follow-up, death or changing to non-anti-TNF biological agents within 3 years (figure 1). Dropouts experienced a significantly greater number of serious infections during the 6 months prior to dropping out than the patients who remained in the study. This association remained significant after adjustment for other risk factors of dropping out (male gender, elevated DAS28 scores, poor functional capacity) and repeated significance testing. The adjusted OR for dropping out after a serious infection calculated for each 6-month period ranged from 2.7 (95% CI 1.2 to 6.2) to 4.7 (95% CI 2.5 to 8.8).
Furthermore, patients who developed serious infections were more likely to switch from anti-TNF treatment to DMARD treatment. Since patients with prior serious infection were at increased risk of developing a further serious infection (see below), this switching led to changes in the risk profiles of the treatment groups (data not shown).
Another factor leading to a change in patient risk profile was the decreasing number of patients who received higher doses of glucocorticoids. The decline in the percentage of patients who needed this additional treatment was far more pronounced in the anti-TNF group (figure 2).
To estimate the relative contributions of changes in the clinical status of the patients and of selection processes to the time-varying risk of infection, we applied three different statistical models (see Methods section). Based on model A which adjusted only for the baseline characteristics of the patients, we observed a statistically significant 31% decrease in the IRRs from year 1 to year 2 in the anti-TNF treatment group compared with the DMARD group (table 3). This decrease was comparable to the decrease in crude IRRs (table 2). However, model A disregards the considerable changes in glucocorticoid dose during the follow-up period (figure 2) and therefore misclassifies exposure at follow-up. Taking into account the time-varying change of glucocorticoid dose and the improvement in functional capacity (model B), approximately one-third (1 − (1 − 0.79)/(1 − 0.69)) of the decrease observed in model A can be attributed to the efficacy of TNF inhibitors. Dropout and treatment changes (either to non-anti-TNF biological agents or from anti-TNF to conventional DMARD treatment) in patients at risk of infection are responsible for the remaining two-thirds of the decline.
The additional adjustment for dropout processes and treatment adaptations in model C led to changes in the IRRs compared with model B which were observed in the expected direction (table 3). Overall, these changes resulted in an IRR for trend of 1.0, suggesting that there is no time-dependent decline in the risk of infection with TNF inhibition if all risk factors are controlled for (table 3).
We found a significantly increased risk of serious infection in patients who had already developed a serious infection earlier in the observation period. Furthermore, older age, chronic lung disease, chronic renal disease and a high number of previous treatment failures increased the risk, whereas a good functional status reduced the risk. The risk was significantly increased in patients treated with glucocorticoids in a dose of ≥15 mg/day (IRR=4.7) or 7.5–14 mg/day (IRR=2.1) and in patients treated with TNF inhibitors (IRR=1.8). No significant increase in risk was observed for treatment with lower doses of glucocorticoids (IRR=1.1 (0.8; 1.7)) or in patients with other co-morbid conditions (eg, diabetes, IRR=1.3 (0.8; 2.0)). Higher disease activity was not directly associated with an increased infection risk but indirectly via the use of glucocorticoids and decline in function.
Absolute risks in different groups of patients
The fully adjusted model C allows estimation of relative risks (IRRs) and also the calculation of expected incidence rates of serious infection for individual patients depending on their risk profiles. As the IRRs of model C are not time-dependent, these incidence rates are also not time-dependent. They only change in a patient when the risk factors she or he is exposed to change.
The rates shown in figure 3 (and in file 1 in the online supplement) reflect how the infection risk in an individual patient is influenced by treatment with different doses of glucocorticoids or by treatment with anti-TNF agents compared with DMARDs. The rates further describe how the risk increases if a patient has one, two or three risk factors for serious infection (greater age, comorbid conditions such as chronic lung disease or chronic renal disease, or history of serious infection). The risk of infection increases steadily with the number of risk factors and with the dose of glucocorticoids, both in the anti-TNF and DMARD treatment groups.
The number of patients with two or more risk factors and treatment with glucocorticoids in a dose of ≥15 mg/day was rather low (70 exposed to anti-TNF, 39 treated with DMARDs). Nevertheless, the estimated incidence rates from model C were in accordance with the observed incidence rates in these patients. Taking into account the duration of exposure to high-dose glucocorticoids, the specific risk factors and the functional capacity of these patients, we would have expected 12.8 (95% CI 6.0 to 27.6) serious infections in patients exposed to anti-TNF agents. This modelling aligned well with the 13 serious infections that were actually observed in these 70 patients. A similar result was found for the 39 DMARD-treated patients in whom three serious infections were observed and 3.7 (95% CI 1.6 to 8.5) were predicted by the model.
Long-term safety data from observational studies reflect time-dependent changes in the risk profiles of individual patients as well as changes in the composition of the cohorts. Our study shows how variable the results from these studies can be, depending on the statistical model applied to adjust for differences between groups.
Using our statistical model A, we saw a similar decline in the absolute and relative risk of infection in the anti-TNF group in the second year of treatment, as has been reported from other observational data.5 7 13 Taking into account changes in clinical status (model B) and the effect of treatment changes and dropouts (model C), we were able to explain this decrease and to estimate the relative contributions of both processes. We found that approximately one-third of the decrease in risk was caused by improvement in the clinical status and reduction of concomitant glucocorticoid therapy (mirroring a reduction of risk in individual patients). The remaining two-thirds of the decrease in risk could be explained by selective switching of patients who were at increased risk of infection. Both processes (clinical improvement and depletion of susceptible patients) led to a 'healthy drug survivor effect' in the anti-TNF cohort—that is, those patients who did well and responded to treatment remained under therapy. Of note, such a 'healthy drug survivor effect' is also present in extension studies of RCTs and might explain the contrasting results of two meta-analyses.1 2 With regard to the question as to what are the risk factors of serious infection, we found that age, functional status and comorbidity contributed significantly to the overall infection risk, in addition to treatment with TNF inhibitors or glucocorticoids.
The adjusted relative risk of treatment with TNF inhibitors compared with that of conventional DMARDs was 1.8. Furthermore, the adjusted relative risk of glucocorticoids was 2.1 for a dose of ≥7.5 mg/day and 4.7 for a dose of ≥15 mg/day. This result underlines the high risk conveyed by glucocorticoids that has also been addressed by other researchers.8,–,11 21 22
The risk imposed on a patient by one treatment must be balanced against the risk conferred by an alternative treatment. For example, consider a patient with highly active RA, aged 65, treated with methotrexate and 10 mg/day glucocorticoids. This patient has an 'average' risk of developing 0.032 serious infections per year. This risk increases to 0.058 if anti-TNF treatment is started, but it may decrease to 0.027 if glucocorticoid treatment can be reduced below 7.5 mg/day due to the response to anti-TNF treatment. The number needed to harm in this example decreases from 1 out of 31 with DMARD treatment plus 10 mg glucocorticoids to 1 out of 17 after start of anti-TNF and increases again to 1 out of 37 after stopping glucocorticoid treatment and maintaining TNF inhibition.
According to our data, improvement in function as a result of effective treatment can reduce the risk of serious infection significantly and to a greater extent than improvement in disease activity measured by the DAS28. This result is clinically plausible when one considers that immobility is a strong risk factor for developing comorbidities such as pneumonia or urogenital infections. Cardiovascular risk factors, as discussed prepviously,7 10 are linked to age and were therefore not identified in our data.
Our analysis has some limitations. We investigated only the overall infection risk and were therefore not able to consider site-specific risk factors. There was also limited statistical power to analyse the data from the third year. We did not investigate the three TNF inhibitors separately. Although different effects on the reactivation of tuberculosis23 and viral infections17 have been reported for adalimumab, infliximab and etanercept, the overall infection risk appears to be similar.3 5 13
The important message from this study to the practising physician is that the time-dependent decline in risk seen in cohorts of patients treated with anti-TNF agents is partially caused by changes in the case mix of these cohorts and cannot directly be transferred to individual patients. Independent of treatment duration, TNF inhibition imposes an increased risk of serious infection on a patient which must be balanced against other risks resulting from co-medication, in particular higher doses of glucocorticoids and active disease.
Our data provide, for the first time, detailed information about the expected absolute risk of serious infection in subgroups of patients. This work therefore helps rheumatologists to identify patients at increased risk and to avoid combinations of treatments that confer a very high risk in predisposed patients.
The authors acknowledge the invaluable contributions of all participating consultant rheumatologists. In particular, they would like to thank those who enrolled 25 patients or more: Ulrich von Hinüber, Hildesheim; Andreas Krause, Immanuel Hospital Berlin; Karin Babinsky, Halle; Arnold Bussmann, Geilenkirchen; Hans Peter Tony, Medizinische Poliklinik der Universität Würzburg; Katja Richter, Universitätsklinikum Carl Gustav Carus, Dresden; Anett Grässler, Pirna; Elke Wilden, Köln; Michael Hammer, St Josef-Stift Sendenhorst; Edmund Edelmann, Bad Aibling; Christina Eisterhues, Braunschweig; Wolfgang Ochs, Bayreuth; Thomas Karger, Eduardus-Krankenhaus Köln-Deutz; Michael Bäuerle, Universität Erlangen, Erlangen; Herbert Kellner, München; Silke Zinke, Berlin; Angela Gause, Elmshorn; Lothar Meier, Hofheim; Karl Alliger, Zwiesel; Martin Bohl-Bühler, Brandenburg; Carsten Stille, Hannover; Susanna Späthling-Mestekemper and Thomas Dexel, München; Harald Tremel, Hamburg; Stefan Schewe, Medizinische Poliklinik der Ludwig-Maximilians-Universität München; Helmut Sörensen, Krankenhaus Waldfriede Berlin; Florian Schuch, Erlangen; Klaus Krüger, München; Andreas Teipel, Leverkusen; Kirsten Karberg, Berlin; Gisela Maerker-Alzer and Dorothea Pick, Holzweiler; Volker Petersen, Hamburg; Kerstin Weiss, Lichtenstein; Werner Liman, Ev. Krankenhaus Hagen-Haspe; Kurt Gräfenstein, Johanniter-Krankenhaus im Fläming, Treuenbrietzen; Jochen Walter, Rendsburg; Werner A Biewer, Saarbrücken; Roland Haux, Berlin; Wolfgang Gross, Lübeck; Michael Zänker, Evangelisches Freikirchliches Krankenhaus Eberswalde; Gerhard Fliedner, Osnabrück; Thomas Grebe, Ev. Krankenhaus Kredenbach; Karin Leumann, Riesa; Jörg-Andres Rump, Freiburg; Joachim Gutfleisch, Biberbach; Michael Schwarz-Eywill, Evangelisches Krankenhaus Oldenburg; Kathrin Fischer, Greifswald; Monika Antons, Köln. The authors also acknowledge the significant contributions of Bernhard Manger, University of Erlangen, Jörn Kekow, University of Magdeburg and Peter Herzer, Munich for their service as members of the advisory board. Their work on the RABBIT advisory board is honorary and without financial compensation. The authors gratefully recognise the substantial contributions of Christina Bungartz, Ulrike Kamenz, Susanna Zernicke and all employees of the German Rheumatism Research Centre, Berlin in the study monitoring and support of the data analyses. The authors would like to thank Geert Molenberghs (University of Hasselt, Belgium) for sharing his experiences in using weighted GEE models and James Galloway (University of Manchester, UK) for his helpful comments related to this manuscript.
Funding RABBIT has been supported by an unconditional joint grant from Essex Pharma/MSD (since 2001), Wyeth Pharma/Pfizer (since 2001), Amgen (2003–8), Swedish Orphan Biovitrum (since 2009) and Abbott (since 9/2003). Since 2007 it has also been supported by Bristol-Myers Squibb and Roche. Since 2009 it has received support from UCB. The principal investigators and their team have full academic freedom in study design and conduct, data analysis and publication of results. These stipulations are delineated in their contract with the sponsors. However, for the purposes of information, all seven funding companies receive every manuscript 30 days prior to submission. The data interpretation, drafting, critical revision and approval of the final manuscript were performed solely by the authors.
Competing interests None.
Patient consent Obtained.
Ethics approval This study was approved in 2001 by the ethics committee of the Charité University Medicine, Berlin.
Provenance and peer review Not commissioned; externally peer reviewed.
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.