Article Text

Download PDFPDF

Original research
Impact of systemic lupus erythematosus disease activity, hydroxychloroquine and NSAID on the risk of subsequent organ system damage and death: analysis in a single US medical centre
  1. Deanna D Hill1,
  2. Amanda M Eudy2,
  3. Peter J Egger3,
  4. Qinggong Fu1 and
  5. Michelle A Petri4
  1. 1Real World Evidence, Epidemiology, GlaxoSmithKline, Collegeville, Pennsylvania, USA
  2. 2Duke University School of Medicine, Durham, North Carolina, USA
  3. 3Real World Evidence, Epidemiology, GlaxoSmithKline, Uxbridge, UK
  4. 4Department of Medicine, Johns Hopkins University School of Medicine, Baltimore, Maryland, USA
  1. Correspondence to Dr Michelle A Petri; mpetri{at}jhmi.edu

Abstract

Objective To assess the impact of mild-moderate systemic lupus erythematosus (SLE) disease activity during a 12-month period on the risk of death or subsequent organ system damage.

Methods 1168 patients with ≥24 months of follow-up from the Hopkins Lupus Cohort were included. Disease activity in a 12-month observation period was calculated using adjusted mean Safety of Estrogens in Lupus Erythematosus National Assessment (SELENA) version of the SLE Disease Activity Index (SLEDAI), defined as the area under the curve divided by the time interval. Damage accrual in the follow-up period was defined as change in Systemic Lupus International Collaborating Clinics/American College of Rheumatology Damage Index (SDI) score ≥1 among patients without prior damage. Patients visited the clinic quarterly and had SELENA-SLEDAI and SDI assessed at every visit.

Results During follow-up (median 7 years), 39% of patients accrued new damage in any organ system (7% cardiovascular and 3% renal) and 8% died. In adjusted models, an increased SELENA-SLEDAI score increased the risk of death (HR=1.22, 95% CI 1.13 to 1.32, p<0.001), renal damage (HR=1.24, 95% CI 1.08 to 1.42, p=0.003) and cardiovascular damage (HR=1.17, 95% CI 1.07 to 1.29, p<0.001). Hydroxychloroquine use reduced the risk of death (HR=0.46, 95% CI 0.29 to 0.72, p<0.05) and renal damage (HR=0.30, 95% CI 0.13 to 0.68, p<0.05). Non-steroidal anti-inflammatory drug use increased the risk of cardiovascular damage (HR=1.66, 95% CI 1.04 to 2.63, p<0.05). Without prior damage, an increased adjusted mean SELENA-SLEDAI score increased the risk of overall damage accrual (HR=1.09, 95% CI 1.04 to 1.15, p<0.001).

Conclusions Each one-unit increase in adjusted mean SELENA-SLEDAI during a 12-month observation period was associated with an increased risk of death and developing cardiovascular and renal damage.

  • lupus erythematosus
  • systemic
  • anti-inflammatory agents
  • non-steroidal
  • antirheumatic agents

Data availability statement

Data are available on reasonable request.

http://creativecommons.org/licenses/by-nc/4.0/

This is an open access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited, appropriate credit is given, any changes made indicated, and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/.

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.

Key messages

What is already known about this subject?

  • Previous studies have demonstrated that systemic lupus erythematosus (SLE) disease activity and damage accrual are strong predictors of SLE prognosis and survival.

What does this study add?

  • This study involved patients with SLE from the large and racially diverse Hopkins Lupus Cohort, composed of 55% white patients and 39% patients of Black African Ancestry. This prospective analysis demonstrated that increases in adjusted mean Safety of Estrogens in Lupus Erythematosus National Assessment-SLE Disease Activity Index (SELENA-SLEDAI) score during a 12-month observational period were associated with an increased risk of death or developing renal or cardiovascular damage.

  • Chronic non-steroidal anti-inflammatory drug (NSAID) use resulted in significantly increased risk of developing cardiovascular damage.

  • SLE may progress to irreversible damage in selected organ systems among patients with stable, mild-to-moderate SLE disease activity over 7 years.

How might this impact on clinical practice or future developments?

  • This analysis demonstrates that, even in patients with mild-to-moderate disease activity, an increase in adjusted mean SELENA-SLEDAI score and NSAID use during a 12-month period increases the risk of death and organ specific damage, highlighting the need for more active measures to manage SLE disease activity over time and to limit NSAID use.

Introduction

Systemic lupus erythematosus (SLE) is an autoimmune disease characterised by a complex and multifactorial aetiology. It affects multiple organ systems, which contributes to the complicated patient burden of SLE. Some patients with SLE experience a waxing and waning pattern of disease activity during the course of disease, whereas others have continuous activity.1 Cumulative SLE disease activity2 and SLE therapies3 4 impact the risk of developing organ system damage and survival.5–7 The 5-year survival among individuals with SLE has increased dramatically in the last four decades to 95%,7 with similar trends observed for 10-year survival.7 8 The impact of irreversible organ system damage in the prognosis of SLE remains a major concern because patients who develop damage are more likely to accrue additional damage and die.9 Published research from prospective cohorts has demonstrated that SLE disease activity is a strong predictor of SLE prognosis and survival.10 11

The objective of this analysis was to investigate whether the pattern of SLE disease activity (as measured by the Safety of Estrogens in Lupus Erythematosus National Assessment (SELENA) version of the SLE Disease Activity Index (SLEDAI)) during a 12-month period after enrolment in a racially diverse cohort of patients prospectively followed quarterly by protocol impacted the risk of developing damage (overall or by specific organ system) and death over time.

Methods

Patient population

Demographic, clinical and laboratory data were collected prospectively in enrolled patients in the Hopkins Lupus Cohort from 1987 to 2010.12–15 Cohort patients with a minimum of 24-month follow-up were included in this analysis. Patients were seen by a single rheumatologist in the Hopkins Lupus Center (MAP) and had a clinical diagnosis of SLE according to the revised American College of Rheumatology (ACR) or Systemic Lupus International Collaborating Clinics (SLICC) criteria.16–18 The study was approved by the Johns Hopkins University School of Medicine Institutional Review Board. All patients provided informed written consent. Patients with a SLICC/ACR Damage Index (SDI) score ≥3, end-stage kidney disease, history of major organ transplant or malignant neoplasm at cohort entry were excluded from this analysis.

Patient follow-up

Follow-up time in the cohort was then divided into three discrete periods: (a) background (background period: first 12 months after cohort entry), (b) observation (observation period: second 12 months after cohort entry) and (c) follow-up (follow-up period: remainder of time under observation until damage occurred, death or end of available data). SLE disease duration was defined as the time interval from SLE clinical diagnosis (revised ACR or SLICC criteria) until the end of the follow-up period or death.

Study variables

Key variables included in the observation period were demographic characteristics, SLE clinical characteristics, laboratory parameters, medication use and comorbid medical conditions (hypertension, diabetes mellitus, obesity, etc.). Hypertension was defined as systolic blood pressure greater than 140 mm Hg or diastolic blood pressure greater than 90 mm Hg or requiring antihypertensive medication. Diabetes mellitus was defined as requiring an oral hypoglycaemic agent or insulin therapy. Immunosuppressive therapies excluded steroids and included the following: leflunomide, mycophenolate, cyclophosphamide, tacrolimus, azathioprine, methotrexate or rituximab (biologic). Antimalarial therapies were almost exclusively comprised of hydroxychloroquine and non-steroidal anti-inflammatory drug (NSAID) therapy included several drugs (prescription and over-the-counter) in this class, but almost exclusively naproxen. For the aforementioned SLE therapies, use was defined as ever prescribed during the observation period (yes/no). Oral steroid (prednisone) use in the observational period was stratified into the following groups: ever prescribed any dose (yes) or never prescribed (no). Among those prescribed oral steroids, daily dose during the observation period was stratified into two categories: (a) prescribed >7.5 mg/day at least once during the observation period or (b) only prescribed ≤7.5 mg/day during the observation period. Antinuclear antibody positivity (ANA+) was defined as a titre ≥1:80 and antidouble stranded DNA positivity (anti-dsDNA+) was defined as a titre ≥1:10 on the Crithidia luciliae indirect immunofluorescence test.

Disease activity was measured at every clinic visit (quarterly by protocol or more often if clinically warranted) using the SELENA-SLEDAI.19 To describe disease activity over time in this analysis, adjusted mean SELENA-SLEDAI20 21 was calculated from the cumulative area under the curve for each period divided by the time interval under evaluation and reported in months. At least one clinic visit per period was required to estimate adjusted mean SELENA-SLEDAI for each period. Damage accrual was assessed at each visit by the SDI.22 23

Data analysis

All statistical analyses were completed in SAS V.9.13 (SAS Institute, Cary, North Carolina, USA).

Continuous variables were reported as median values unless otherwise noted. Categorical variables were described as frequency counts. SLE outcomes of interest included (a) development of overall organ system damage, (b) development of specific organ system damage or (c) death during the follow-up period. The outcome of death was defined as all-cause mortality that may or may not have been clinically attributable to SLE. Damage accrual analyses were restricted to patients without prior damage in the organ system of interest through the end of the observation period. Accrued damage was defined as at least a one (≥1) unit increase in SDI score between the observation period and last recorded visit in the follow-up period.

In this analysis, organ specific system damage was restricted to renal, cardiovascular, peripheral vascular, neuropsychiatric, pulmonary, musculoskeletal, seizure or stroke. Damage due to malignancy was not explored as it was an exclusion criterion. Damage accrual in organ systems not commonly associated with death (ocular, gastrointestinal, skin, or gonadal failure) was not evaluated in this analysis.

Adjusted Mean SELENA-SLEDAI was our key explanatory variable. Laboratory parameters were not considered in models evaluating the impact of adjusted mean SELENA-SLEDAI on the outcomes of interest since these variables are components of the SELENA-SLEDAI score. Time-varying variables that were measured during the observation period and SLE clinical history variables were individually evaluated in univariate, unadjusted models for the outcomes of death or damage accrual (overall or specific organ systems). When more than one variable was available to evaluate comorbidity status, the variable that was most relevant to the time period of interest was used. If variables were comparable, then both were retained in the univariate analyses and the results from the univariate analysis determined which variable was considered in the multivariable model (when applicable). Variables with a p<0.10 in the univariate analysis were considered as potential covariates in the multivariable models for the outcomes of death or damage accrual. Cox proportional hazard models were used to estimate the impact of adjusted mean SELENA-SLEDAI as a time-dependent variable on the risk of death or developing any new organ damage over time. Kaplan–Meier survival curves were also generated to explore the impact of adjusted mean SELENA-SLEDAI in the observation period on the risk of death over time during the follow-up period.

For each outcome, a multivariable model was created to adjust for well-established potential confounding factors that included age, gender, race, SLE duration, overall SDI score at the start of the follow-up period and ever prescribed oral prednisone >7.5 mg/day in the observation period. With these factors forced in the model, other covariates identified in the univariate analyses were then evaluated for inclusion in the model using backward elimination methods with the criteria of p<0.10 for entry into the model and p<0.05 for retention. To evaluate racial and gender differences in the association of adjusted mean SELENA-SLEDAI and organ system damage, models were stratified by race and gender (data not shown).

Patient and public involvement

We did not involve patients and/or the public in this work.

Results

Characteristics

Overall, 1168 adult patients with SLE were included in this study. Patients were predominantly women (93%), 55% were white and 39% were of Black African Ancestry (table 1). During the background period, the median (range) adjusted mean SELENA-SLEDAI was 3 (0–16) and 55% of patients had mild-to-moderate disease activity, defined by an adjusted mean SELENA-SLEDAI <3 (table 1). Median (range) adjusted mean SELENA-SLEDAI in the observation period was 2 (0–13) and remained the same during the follow-up period. The distribution of adjusted mean SELENA-SLEDAI, stratified by period, is shown in figure 1.

Figure 1

Distribution of adjusted mean SELENA-SLEDAI stratified by period (N=1168). SELENA-SLEDAI, Safety of Estrogens in Lupus Erythematosus National Assessment-Systemic Lupus Erythematosus Disease Activity Index.

Table 1

Demographic and clinical characteristics of the Hopkins Lupus Cohort analytical cohort (N=1168)

Death

Ninety-two of 1168 patients (8%) died during the follow-up period (table 2). When exploring the association between adjusted mean SELENA-SLEDAI in the observation period and risk of death in the follow-up period, an adjusted mean SELENA-SLEDAI ≥5 significantly increased the risk of death after 7 years of follow-up (Log rank p<0.001; figure 2). In the multivariable model, each one-unit adjusted mean SELENA-SLEDAI increase in the observation period was significantly associated with a 22% increased risk of death in the follow-up period (HR=1.22, 95% CI 1.13 to 1.32, p<0.001; table 2). Use of hydroxychloroquine in the observation period was significantly associated with a 54% lower risk of death in the follow-up period (HR=0.46, 95% CI 0.29 to 0.72, p<0.05; table 2). The following variables were considered but not retained in the model: have a history of smoking at cohort entry, hypertension and diabetes.

Figure 2

Kaplan–Meier curve: time to death by adjusted mean SELENA-SLEDAI score (N=1168). AMS, adjusted mean SELENA-SLEDAI; SELENA-SLEDAI, Safety of Estrogens in Lupus Erythematosus National Assessment-Systemic Lupus Erythematosus Disease Activity Index.

Table 2

Final models for adjusted mean SELENA-SLEDAI in the observation period as predictors of death/organ system damage in the follow-up period

Overall damage accrual

After restricting the analysis (univariate and multivariable) to patients without a prior history of overall damage (SDI=0; n=888) at the start of the follow-up period, more than one-third (39%) of patients developed damage (SDI≥1) in any organ system by the end of the follow-up period (table 2). In multivariable models restricted to patients without prior damage in the organ system at the start of the follow-up period, adjusted for the effects of age, race, gender, SLE duration and ever prescribed >7.5 mg/day oral prednisone during the observation period, a one-unit increase in adjusted mean SELENA-SLEDAI during the observation period was significantly associated with an increased risk of accruing subsequent overall damage (HR=1.09, 95% CI 1.04 to 1.15, p<0.001; table 2). Among women without prior damage at the start of the follow-up period, adjusted mean SELENA-SLEDAI during the observation period was associated with an increased risk of accruing subsequent overall damage (HR=1.08, 95% CI 1.03 to 1.13, p=0.002; data not shown). No significant differences were observed for the outcomes of interest when stratified by race (HR=1.07, 95% CI 0.99 to 1.16, p=0.089 in white patients and HR=1.04, 95% CI 0.96 to 1.12, p=0.343 in patients of Black African Ancestry; data not shown).

Renal damage accrual

Approximately 3% of patients without renal damage in the observational period developed renal damage (SDI≥1; n=1147) in the follow-up period (table 2). Each one-unit increase in adjusted mean SELENA-SLEDAI during the observation period significantly increased the risk of renal damage by 24% (HR=1.24, 95% CI 1.08 to 1.42, p=0.003; table 2) in the follow-up period. Ever use of hydroxychloroquine in the observation period was significantly associated with lowering the risk of renal damage in the follow-up period by 70% (HR=0.30, 95% CI 0.13 to 0.68, p<0.05; table 2). When we restricted the analysis to women, the results did not change (HR=1.23, 95% CI 1.06 to 1.42, p=0.005). It was not possible to evaluate renal damage among men due to only one man having developed renal damage.

Cardiovascular damage accrual

Subsequent cardiovascular damage (SDI≥1; n=1135) was observed in approximately 7% of patients during the follow-up period (table 2). A one-unit increase in adjusted mean SELENA-SLEDAI during the observation period was associated with a 17% (HR=1.17, 95% CI 1.07 to 1.29, p<0.001; table 2) significant increased risk of cardiovascular damage accrual in the adjusted model. Patients who had been prescribed NSAID therapies in the observation period had a 66% (HR=1.66, 95% CI 1.04 to 2.63; table 2) increased risk of cardiovascular damage accrual. On further exploration, patients with any antihypertensive use in the observation period had an 81% (HR=1.81, 95% CI 1.09 to 3.02, p<0.05) significant increased risk of cardiovascular damage accrual.

Other organ system or subsystem damage accrual

Adjusted mean SELENA-SLEDAI during the observation period did not impact the risk of developing damage in the following organ systems or subsystems in our analysis: peripheral vascular (HR=1.12, 95% CI 0.96 to 1.32, p=0.154), pulmonary (HR=0.96, 95% CI 0.87 to 1.07, p=0.484), neuropsychiatric (HR=1.01, 95% CI 0.92 to 1.10, p=0.880), musculoskeletal (HR=1.05, 95% CI 0.98 to 1.12, p=0.194), stroke or seizure (data not shown).

Discussion

This analysis described the impact of disease activity during a 12-month period on the subsequent risk of death or organ system damage in a large, prospective cohort of racially diverse patients with SLE. This study included a large proportion of patients of Black African Ancestry (39.4%), which contrasts with previous Spanish studies (eg, RELESSER cohort)24–26 where the majority of patients were white (>90%), and Latin American studies (eg, GLADEL),27–29 which had lower percentages of African-Latin American patients (<13%). In this cohort of patients with SLE and mild-to-moderate disease activity at cohort entry, adjusted mean SELENA-SLEDAI measured during a prior 12-month period (corresponding to the length of a typical Phase III trial) significantly impacted the risk of renal and cardiovascular damage accrual (SDI≥1) and risk of death in the follow-up period. In previous studies, patients of Black African Ancestry were reported to experience a higher prevalence of SLE, greater disease severity, higher risk of cardiovascular events, greater organ damage and higher mortality rates, compared with white patients.30–37 In the current study, although non-white (Black African Ancestry and other race) patients had higher disease activity in the follow-up period compared with white patients, no significant racial differences were observed for the outcomes of interest. Therefore, results from this cohort provide an important contribution in characterising non-white patient populations. In adjusted Cox regression models, a one-unit increase in adjusted mean SELENA-SLEDAI during the 12-month observation period was associated with a 22% increased risk of death and a 17% and 24% increased risk of subsequent cardiovascular and renal damage accrual, respectively. Adjusted mean SELENA-SLEDAI in the observation period was associated with overall damage accrual after adjustment for age, gender, race, SLE duration, ever prescribed oral prednisone therapy >7.5 mg/day in observation period and SDI at the start of the follow-up period. Hydroxychloroquine use in the observation period was associated with a 54% and 70% decreased risk of death and renal damage, respectively, and NSAID use in the observation period increased the risk for cardiovascular damage accrual by approximately 70%.

Adjusted mean SELENA-SLEDAI in a 12-month prior period increased the risk of developing new renal and cardiovascular damage accrual (SDI≥1) and risk of death, which corroborates findings from another published Hopkins Lupus Cohort analysis with different analytical methods.38 Other studies have demonstrated that high disease activity increased the risk of poor SLE outcomes including accrual of overall organ damage20 39–44 and risk of death.7 40 42 45 The findings in this analysis corroborate the influence of disease activity for renal and cardiovascular damage accrual and death and also extend the findings to patients with SLE and mild-to-moderate disease activity.

Although 61.5% of the cohort remained free of any organ system damage through the end of the follow-up period, a reasonable proportion of patients with SLE with mild-to-moderate disease activity at cohort entry and through the start of the follow-up period accrued damage (renal and cardiovascular) (9.2%) in a relatively short time (median 7 years). Similarly, 8% of patients died during follow-up despite overall mild-to-moderate disease activity in the 24 months after cohort entry.

Detailed methods from real-world, large, prospective SLE cohorts have been described elsewhere10–13 46 47 and have contributed to SLE disease understanding over the past four decades. It is well recognised that damage accrual, due to active inflammation, comorbidities, previous SLE disease activity and/or exposure to SLE therapies (particularly corticosteroids), is an important prognostic factor for death and has a clear impact on burden of disease for patients,15 20 40 48 49 but the time for interval for damage to manifest has varied somewhat between different studies. Rahman et al found that patients with SLE exhibiting damage within the first year of admission to the clinic had a higher mortality rate after 10 years as compared with patients with no early damage,50 whereas results from Gladman et al suggested that damage accrued gradually over the 15 years of follow-up.2 In a study by Becker-Merok and Nossent, damage accrual occurred in 54% of patients with SLE in a linear fashion over the first 10 years of the disease.42 Similarly, a Swedish study also observed that 54% of patients with SLE had damage accrual in the first 5 years after diagnosis.51 Our findings demonstrate that SLE may progress to accrual of irreversible damage in selected organ systems among patients with stable, mild-to-moderate disease activity over 7 years. Our study confirms the results from a Norwegian SLE cohort that showed disease activity at baseline was a predictor of accrued organ damage that occurred over 2 years of follow-up in a small cohort of patients with relatively stable disease activity52 but extends them to specific organ damage subtypes. Unlike the findings from the SLICC inception cohort, we did not observe a pattern of decreasing disease activity over the first 5 years in conjunction with organ damage accrual.10 This may be in part due to the combination of inception and prevalent patients with SLE at the start of the observation period in this Hopkins Lupus Cohort analysis.

Subsequent accrual of organ system damage attributed to active SLE disease has been observed in several SLE prospective cohorts of patients with SLE receiving care from rheumatology specialty centres.5 10 11 20 48 53–55 One strength of this analysis is that Hopkins Lupus Cohort patients had clinic visits on average every 3 months per registry protocol. Ibañez et al recently demonstrated that adjusted mean SELENA-SLEDAI derived from quarterly clinic visits was more reliable than adjusted mean SELENA-SLEDAI based on more infrequent clinic visits.56

We recognise the limitations of this analysis and similar evaluations in prospective SLE cohorts. One limitation is that all patients in this analysis received care at a single tertiary medical centre, under the care of a single provider, and their SLE clinical characteristics, and hence our findings, may not extend to all patients with SLE. Furthermore, although we attempted to adjust for the effects of known risk factors and potential confounders, unmeasured or residual confounding factors may have influenced our findings.

We observed that adjusted mean SELENA-SLEDAI during a prior 12-month period impacted the risk of death and developing damage in a racially diverse cohort of patients with SLE, who had on average mild-to-moderate disease activity during a median follow-up of 7 years, after adjusting for potential confounders. Exposure to hydroxychloroquine during a prior 12-month period decreased the risk of death and developing renal damage. The London University College Hospital Lupus Cohort reported that more than half of patients with SLE were prescribed hydroxychloroquine during a 12-month observation period and this exposure was also associated with a decreased risk of death (50% reduction) and renal damage accrual (47% reduction) in adjusted models.11 A new finding in our analysis was the 66% increased risk of cardiovascular organ damage accrual associated with NSAID use in the previous 12 months, after controlling for the effects of other covariates in the model. It has been reported that NSAIDs may negatively affect the cardiovascular system, yet they also decrease inflammation, which is an independent risk factor in cardiovascular pathology.57 In our analysis, we observed that chronic NSAID use, commonly taken by patients with SLE to alleviate musculoskeletal pain, resulted in a significantly increased risk of cardiovascular damage accrual. NSAID use is also correlated with an increase in blood pressure, which may have an effect on cardiovascular damage.58 Indeed, in our study, we observed a significant increase in cardiovascular damage accrual with antihypertensive use. This may suggest that the known cardiovascular risk of NSAIDs in the general population is also applicable to patients with SLE59–62 and highlights the importance of assessing cardiovascular risk in this patient population.

In summary, our findings corroborate other published data that demonstrated that cumulative SLE disease activity over time impacted risk of developing organ damage (ie, SLE prognosis) and was associated with an increased risk of death. Furthermore, the findings suggest routine clinical care and SLE disease management, characterised by minimal fluctuations in disease activity, may not prevent development of organ damage or risk of premature death. An increase in adjusted mean SELENA-SLEDAI score during a 12-month period increased the risk of death and developing renal and cardiovascular organ system damage, even in patients with mild-to-moderate disease severity, which underscores the need for active measures to manage SLE disease activity over time.

Data availability statement

Data are available on reasonable request.

Ethics statements

Ethics approval

The study was approved by the Johns Hopkins University School of Medicine Institutional Review Board (Study number NA_00039294).

Acknowledgments

Medical writing support was provided by Helen Taylor, of Fishawack Indicia Ltd, UK, part of Fishawack Health, and was funded by GSK.

References

Footnotes

  • Contributors Conception or design: DDH, PJE, QF, MAP and AME. Acquisition of data: MAP. Data analysis or interpretation: DDH, PJE, QF, MAP and AME.

  • Funding The Hopkins Lupus Cohort is supported by grants R01AR043727 and R01AR069572 from the National Institutes of Health (NIH). This work was funded by GSK Study Number: WEUKBRE4566. Johns Hopkins University School of Medicine received a research grant to support this study. AME is currently supported by NIH NCATS Award Number 1KL2TR002554.

  • Competing interests DDH, PJE and QF were paid employees of GlaxoSmithKline (GSK) with stock options and conducted the study as part of their employment using GSK resources. AME was a paid contractor for GSK. MAP was a paid consultant to GSK.

  • Patient and public involvement Patients and/or the public were not involved in the design, or conduct, or reporting, or dissemination plans of this research.

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

  • Author note DDH current affiliation is Center for Observational and Real-World Evidence, Pharmacoepidemiology, Merck & Co, Inc, Kenilworth, New Jersey, USA.