Abstract
Objective. Systemic autoimmune rheumatic diseases (SARD) are associated with an increased risk of premature cardiovascular disease (CVD) and all-cause mortality. We examined the potential survival benefit of statin use among patients with SARD in a general population setting.
Methods. We conducted an incident user cohort study using a UK general population database. Our population included patients with a SARD as determined by Read code diagnoses of systemic lupus erythematosus, systemic sclerosis, Sjögren syndrome, dermatomyositis, polymyositis, mixed connective tissue disease, Behçet disease, or antineutrophil cytoplasmic antibodies-associated vasculitis between January 1, 2000, and December 31, 2014. We compared propensity score–matched cohorts of statin initiators and noninitiators within 1-year cohort accrual blocks to account for potential confounders, including disease duration, body mass index, lifestyle factors, comorbidities, and medication use.
Results. Of 2305 statin initiators, 298 died during the followup period (mean 5.1 yrs), whereas among 2305 propensity score–matched noninitiators, 338 died during the followup period (mean 4.8 yrs). This corresponded to mortality rates of 25.4/1000 and 30.3/1000 person-years, respectively. Statin initiation was associated with reduced all-cause mortality (HR 0.84, 95% CI 0.72–0.98). When we compared the unmatched cohorts, the statin initiators (n = 2863) showed increased mortality (HR 1.85, 95% CI 1.58–2.16) compared with noninitiators (n = 2863 randomly selected within 1-year cohort accrual blocks) because of confounding by indication.
Conclusion. In this general population–based study, statin initiation was shown to reduce overall mortality in patients with SARD after adjusting for relevant determinates of CVD risk.
Systemic autoimmune rheumatic diseases (SARD) are associated with premature mortality, largely through an increased risk of premature cardiovascular disease (CVD)1,2,3,4,5,6. These conditions share the commonalities of chronic systemic inflammation as well as treatment strategies often involving glucocorticoids and other chronic immunosuppressive therapies. Moreover, they are typically managed longitudinally by rheumatologists who often take an active role in the management of comorbidities. In systemic lupus erythematosus (SLE), many studies have shown an increased risk of CVD, including a doubled risk of stroke and myocardial infarction (MI) in the Nurse’s Health Study and other large studies7,8,9,10,11,12, as well as doubled risk of overall mortality13. A large UK population-based study in systemic sclerosis (SSc) found 2.6 times the risk of stroke and 1.8 times the risk of MI relative to peers14. Similarly, in granulomatosis with polyangiitis, a large population-based study in British Columbia, Canada, found a 1.5-times increased risk of stroke and a 1.8-times increased risk of MI15, and a recent metaanalysis of observational studies found 2.7-fold increased mortality associated with antineutrophil cytoplasmic antibody-associated vasculitis (AAV)16. In myositis, overall mortality is increased, and mounting evidence suggests premature CVD risk, as a previous large population-based study found 1.7-times increased risk of stroke and 3.9-times increased risk of MI in dermatomyositis (DM) and polymyositis (PM)17,18,19. With substantial evidence for a collectively higher risk of CVD and mortality for patients with SARD, investigations into strategies to mitigate this risk are urgently needed.
In the general population, a major strategy for reducing the risk of CVD and mortality is treatment with statins. A wide body of literature supports the efficacy of statins for primary and secondary prevention of CVD and associated mortality in the general population. Briefly, a Cochrane review of statin use for CVD prevention in the general population found a 25% reduction in CVD events, and a previous metaanalysis of statin use in the general population found a 9–14% reduction in overall mortality20,21. However, given the rarity of these diseases, evidence regarding the potential mortality benefit of statins among patients with individual SARD is scarce. To help address this gap in knowledge, we conducted a proof-of-concept analysis of the potential survival benefit of statins collectively among individuals with SARD from a large general population database.
MATERIALS AND METHODS
Data source
The setting for our study was The Health Improvement Network (THIN), an electronic medical record database that represents 6.2% of the UK population, including over 11 million patients. The THIN database is representative of the general UK population in demographics and the prevalence of common medical diagnoses22. Health information included in this database includes demographics, lifestyle factors, information from general practitioner (GP) visits, diagnoses from hospital admissions and specialists, medications, and laboratory results. The specific diagnoses are recorded by the Read code classification system, which is the standard nomenclature of clinical terms used by the National Health Service in describing clinical diagnoses23. Medication prescriptions are recorded by the Multifunctional Standardized Lexicon for European Community Language (Multilex) classification system24.
Study design and cohort definition
We conducted a population-based incident user cohort study. The classification of SARD has been previously described2,25. We identified subjects with SARD, defined by having at least 1 Read code for one of the included conditions: SLE, SSc, Sjögren syndrome (SS), DM, PM, mixed connective tissue disease, AAV, or Behçet disease2,25. We did not include individuals with diagnoses of giant cell arteritis or polymyalgia rheumatica to improve comparability within the SARD cohort, because those conditions exclusively affect an older age group26. Read code diagnoses have been previously validated in a similar UK database, and diagnoses of connective tissue diseases were found to be 80% accurate27.
Within this defined population with SARD, we identified subjects who initiated statin medications between January 1, 2000, and December 31, 2014. To account for potential secular trends in statin prescribing in SARD, subjects were divided into 1-year accrual blocks during this time frame. Within each accrual block, statin initiators were defined as patients with new statin prescriptions during that 1-year period. To be considered a statin initiator, a subject also had to be enrolled in the THIN database for at least 1 year prior to the first recorded statin prescription. The comparators were subjects with SARD matched 1:1 by the 1-year accrual blocks of entry into the cohort and who did not initiate statin medications during the study period using a 5-to-1 digit “greedy matching” algorithm28,29. The index date was the date of statin initiation for statin users and a random date within the accrual year for noninitiators30. Individuals were excluded if they were current or prior statin users (i.e., not incident users) or if they had incomplete records of covariates.
Because we expected there would be systematic differences in the baseline characteristics of subjects receiving statin therapy compared to those not prescribed statins, we performed propensity score matching to adjust for these potential imbalances. We calculated propensity scores (the predicted probability of statin initiation) within each 1-year accrual block by logistic regression. The variables included in the propensity score estimation comprised demographics (age, sex), body mass index (BMI), alcohol and tobacco use, healthcare use (as measured by the number of GP visits, specialist referrals, and hospitalizations), specific SARD diagnosis, SARD duration prior to the index date, medication use, comorbid conditions (determined by Read code diagnoses), and cholesterol levels assessed within 1 year prior to the index date of entry into the cohort (Table 1).
Baseline characteristics of statin initiators and noninitiators in the propensity score–matched and unmatched cohorts. Values are % unless otherwise specified.
Assessment of outcome
The primary outcome of interest was all-cause mortality. This was assessed by the death date automatically recorded in the THIN database when death is registered in the Personal Demographics Service database, which contains demographic data for all patients registered with the National Health Service in the UK. This automatic update has been shown to be an accurate reflection of national death rates in the UK22.
Statistical analysis
Descriptive statistics were calculated for baseline characteristics of the statin initiators and noninitiator comparators in both the propensity score–matched and unmatched cohorts. We calculated person-years (PY) of followup for each subject from the index date until either death, end of the study period, or disenrollment from the THIN database. Subjects remained grouped as “statin initiators” regardless of continuation of the statin prescription. This is analogous to an intention-to-treat analysis used in clinical trials and provides conservative estimates for the target effect. Nevertheless, to examine the potential effect of discontinuation of statins among initiators over time, we performed analyses with the followup time truncated at 1, 2, 3, and 4 years for all subjects. We generated survival plots with estimates of cumulative mortality over time. We then used Cox proportional hazard models to estimate the effect of statin initiation on mortality, stratified by 1-year cohort accrual blocks. To examine the potential effect of individual SARD on our effect estimation, we repeated our analysis by excluding 1 subset of all SARD at a time.
This study was approved by the Partners Human Research Committee, approval number 2017P000399. Informed consent was waived because all data were anonymous.
RESULTS
Unmatched analysis
We identified 2863 statin initiators and 2863 noninitiators in the unmatched cohort with complete covariate information (Table 1). An additional 475 statin initiators were excluded because of missing covariates; their baseline covariates did not differ from the included patients (Supplementary Table 1, available with the online version of this article). Statin initiators had higher baseline cholesterol levels and were older, more often male, and had greater comorbidities. There were 232 deaths among statin initiators over a mean followup time of 3.1 years. The overall mortality rate was 26.0 deaths per 1000 PY. In the comparator group, there were 123 deaths over a mean of 3.1 years, with an overall mortality rate of 13.8 deaths per 1000 PY. In this unmatched analysis without controlling for confounding by indication, statin initiation was associated with an overall HR for mortality of 1.85 (95% CI 1.58–2.16) and greater cumulative mortality as depicted in Figure 1.
Cumulative mortality over time in statin initiators and noninitiators.
Propensity score–matched analysis
In the propensity score–matched cohorts (2305 statin initiators and 2305 noninitiators), subjects were well-balanced by age, sex, BMI, SARD diagnosis, disease duration, alcohol and tobacco use, comorbidities, medication use, and cholesterol levels at baseline (Table 1). SLE represented the largest subset of all SARD, followed by SS and SSc. There were 298 deaths among statin initiators over a mean followup time of 5.1 years. The overall mortality rate was 25.4 deaths per 1000 PY. In the comparator group, there were 338 deaths over a mean of 4.8 years, with an overall mortality rate of 30.3 deaths per 1000 PY. Figure 1A depicts the cumulative overall mortality for the 2 exposure groups. The HR for overall mortality associated with statin initiation was 0.84 (95% CI 0.72–0.98) over the complete followup period. With the truncated followup of 1, 2, 3, and 4 years, the overall mortality HR remained < 1 for statin initiation, trending toward significance for 1 and 2 years of followup, and with HR 0.78 (95% CI 0.62–0.98) for 3 years of followup and 0.77 (95% CI 0.63–0.95) for 4 years of followup (Table 2).
Association between statin initiation and all-cause mortality in propensity score–matched cohorts.
Subgroup and sensitivity analysis
In our subgroup analyses, mortality HR associated with statin initiation was 0.81 (95% CI 0.69–0.96) among patients ≥ 55 years and 0.89 (95% CI 0.45–1.66) among patients < 55 years. Females, comprising the majority (76%), had overall mortality HR of 0.77 (95% CI 0.64–0.92) among statin initiators, whereas among males it was 1.05 (95% CI 0.77–1.43). The corresponding HR were 0.82 (95% CI 0.66–1.01) among patients without a prior history of CVD (e.g., coronary artery disease, stroke, or peripheral vascular disease) at baseline and 0.81 (95% CI 0.65–1.02) among patients with such conditions. For users of nonsteroidal antiinflammatory drugs and nonusers, the corresponding HR were 0.90 (95% CI 0.67–1.21) and 0.80 (95% CI 0.67–0.97), respectively (data not shown).
In our analyses per specific SARD subgroups, the mortality HR were 0.83 (95% CI 0.66–1.06) for patients with SLE, 0.87 (95% CI 0.66–1.16) for those with SS, and 0.63 (95% CI 0.42–0.94) for those with SSc. The numbers of subjects with DM, PM, AAV, and Behçet disease were small (Table 1), limiting meaningful individual subgroup analyses. With the exclusion of 1 subset of all SARD at a time from the propensity score–matched cohorts, the effect estimates of the all-cause mortality reduction associated with statin initiation were similar (HR 0.82, 95% CI 0.70–0.96 excluding vasculitides; HR 0.82, 95% CI 0.69–0.96 excluding myositis; HR 0.83, 95% CI 0.68–1.02 excluding SLE; HR 0.82, 95% CI 0.68–0.99 excluding SS, and HR 0.88, 95% CI 0.75–1.05 excluding SSc; Supplementary Table 2, available with the online version of this article).
DISCUSSION
In this large-scale UK general population–based cohort study of subjects with SARD, statin initiation was associated with a 16% reduction of all-cause mortality, compared with propensity score–matched patients with SARD who were not treated with statins. This reduction in overall mortality became significant at 3 years of followup and was sustained during the rest of the study followup. The magnitude of the protective association between statin use and mortality risk reduction was similar to that seen in the general population and appeared to be similar across the included individual SARD conditions as well as for primary and secondary CVD prevention indications, although subgroup analysis was limited by sample size20,21.
The survival benefit associated with statin initiation in patients with SARD may be related to the pleotropic antiinflammatory effects of statins, in addition to their cholesterol-lowering and antiatherosclerotic effects31. Statins have immunomodulatory properties such as altering the function of antigen-presenting cells and T cells, which may provide a direct benefit in these autoimmune diseases32. They have been shown to reduce C-reactive protein levels in patients with SLE33 and may play a role in reducing proteinuria34. In SSc, statins have been implicated in reducing endothelial scarring and reducing the development of digital ulcers35. One study found simvastatin to reduce neutrophil degranulation in patients with ANCA-associated vasculitis36, suggesting a potential role in reducing vascular damage in this disease as well.
Additionally, as newer treatments have improved our ability to control disease activity for patients with SLE, AAV, and other SARD, mortality because of the underlying diseases has declined over time16,37. However, these patients continue to experience multiple risk factors for premature atherosclerosis, including side effects from chronic glucocorticoid exposure38,39 and an increase in traditional CVD risk factors such as smoking and a sedentary lifestyle40, in addition to longstanding chronic inflammation because of the SARD. Therefore, this population at heightened risk of developing CVD may benefit from the effects of statins on reducing atherosclerosis as well as broader antiinflammatory effects41,42.
Previous randomized trials of patients with SLE (N ranging from 60 to 221) have investigated the potential effect of statins on surrogate imaging endpoints of atherosclerosis [i.e., carotid plaque, carotid intima media thickness, and coronary artery calcium (CAC)] with conflicting results41. While 3 trials did not demonstrate an effect of statin therapy on the subclinical atherosclerosis endpoints43,44,45, 1 study showed a reduced progression of CAC in the atorvastatin-treated group46. These studies had limited sample sizes for clinical endpoints such as mortality. A previous Taiwanese study has assessed mortality risk among SLE patients with hyperlipidemia, but the study did not use an incident user design (as was done in ours)47, which is critically important to ensure validity in pharmacoepidemiology48. Further, that study’s calculation of person-time in relation to SLE onset and duration is highly suggestive of immortal time bias49, which would lead to a very protective effect as reported (i.e., mortality HR of high-dose statin use compared with nonusers = 0.45, 95% CI 0.25–0.84)47. We are not aware of any other studies investigating the mortality effect of statins in patients with SLE or other SARD. Thus our study provides data to suggest the efficacy of statins in reducing mortality for patients with SARD.
Our study has several strengths and limitations that warrant recognition. We used a large general population database, which made our study generalizable and allowed us to identify large numbers of individuals with these rare conditions who have sufficient followup to investigate the outcome of mortality. In considering the individual SARD diagnoses together in 1 analysis, we improved our power to detect overall differences in outcomes. However, our study was not powered to specify the relative effect of statins on mortality in each specific SARD condition. To that end, confirming our findings in individual SARD in a much larger dataset or pooled analyses with multiple datasets would be valuable. In particular, further data on the safety of statins among patients with DM and PM would be useful. Until then, the potential concern of concomitant statin myopathy may outweigh the benefit other than in high cardiovascular risk–dyslipidemic patients with these conditions. Nonetheless, our findings demonstrating the benefit of statins in reducing mortality among a collective SARD cohort provides important preliminary evidence. We did not demonstrate the same benefit among males. While this may imply a subgroup effect by sex, the proportion of males was small, limiting meaningful analysis. Further studies would be needed in a larger cohort of males with SARD. Our source population derives from a GP-driven medical record database, and it does not include accurate information regarding disease activity. Therefore, we could not determine the duration of disease flares or details of specific disease manifestations such as Raynaud phenomenon or digital ulcers among patients with SSc. However, it does include numbers of GP visits, hospital admissions, and specialist referrals, which we used to measure healthcare use as a surrogate of illness severity. For ascertainment of the exposure, we relied upon medication prescribing data and could not assess compliance with the statin medications. However, this would tend to provide a conservative estimate of the effect of statin treatment. We considered different potency statins and dosing regimens together in the analysis because of sample size constraints. Future studies could assess the effect of these differences on mortality in patients with SARD. We focused our analysis on all-cause mortality because this outcome has been shown to be highly reliable in our source database, and we had incomplete data on cause-specific mortality22. Regardless, the demonstration of reduced all-cause mortality is itself critically important, as this endpoint reflects the net health outcome of the overall benefits and risks associated with statin use.
We have shown that statin initiation was associated with a reduction in overall mortality for patients with SARD, with a magnitude of effect similar to that previously seen in the general population. The proper use of statins may be beneficial in reducing premature mortality for patients with SARD.
ONLINE SUPPLEMENT
Supplementary material accompanies the online version of this article.
Acknowledgment
The work was presented at the 2016 American College of Rheumatology Annual Meeting.
Footnotes
This project was supported in part by the Ruth L. Kirschstein Institutional National Research Service Award T32-AR-007258 and National Institutes of Health Grant P60-AR-047785.
- Accepted for publication May 11, 2018.