ABSTRACT
Background Evidence for the association between glucocorticoid dose and cardiovascular risk is weak for moderate and low doses. To quantify glucocorticoid dose-dependent cardiovascular risk in people with six immune-mediated inflammatory diseases.
Methods and Findings Population-based cohort analysis of medical records from 389 primary care practices contributing data to the UK Clinical Practice Research Datalink, linked to hospital admissions and deaths in 1998-2017. There were 87,794 patients with giant cell arteritis and/or polymyalgia rheumatica (n=25,581), inflammatory bowel disease (n=27,739), rheumatoid arthritis (n=25,324), systemic lupus erythematosus (n=3951), and/or vasculitis (n=5199); and no prior cardiovascular disease (CVD). Mean age was 56 years and 34.1% were men. Median follow-up time was 5.0 years. Time-variant daily and cumulative glucocorticoid prednisolone-equivalent dose-related risks and hazard ratios of first all-cause and type-specific CVD.
We found 13,426 (15.3%) people with incident CVD, including 6,013 atrial fibrillation, 7,727 heart failure and 2,809 acute myocardial infarction events. At 1 and 5 years, the cumulative risks of all-cause CVD increased from 1.5% in periods of non-use to 9.1% for a daily prednisolone-equivalent dose of ≥25.0mg, and from 7.6% to 29.9%, respectively. We found strong dose-dependent estimates for all immune-mediated diseases (hazard ratio [HR] for <5.0mg daily dose vs. non-use=1.74, 95%CI 1.64-1.84; range 1.52 for polymyalgia rheumatica and/or giant cell arteritis to 2.82 for systemic lupus erythematosus), all cardiovascular outcomes, regardless of disease activity level. The highest estimates were for heart failure and acute myocardial infarction.
Conclusions We estimated glucocorticoid dose-dependent cardiovascular risk in six immune-mediated diseases. Results highlight the importance of prompt and regular monitoring of cardiovascular risk and use of primary prevention treatment at all glucocorticoid doses.
INTRODUCTION
Patients with immune-mediated inflammatory diseases often receive long-term courses of oral glucocorticoids to reduce disease activity and inflammation during the initial episode and subsequent episodic flares. Prolonged glucocorticoid treatment often causes adverse events, including cardiovascular diseases (CVDs)[1-3]. Glucocorticoids can increase cardiovascular risk through direct and indirect metabolic syndrome enhancement[4-6] and mineralocorticoid effects, including cellular membrane electrolyte-mediated efflux[7-9]. However, the anti-inflammatory and immune-suppressive effects of glucocorticoids could also lower or neutralise the atherosclerotic and vascular injury effects of chronic inflammatory diseases[10]. Demonstration of cost-effectiveness of newly licenced glucocorticoid-sparing drugs, such as biologics, is critical to guide their introduction for the treatment of immune-mediated inflammatory diseases in routine healthcare. These studies require the accurate estimation of glucocorticoid dose-response relationships to quantify cost-savings associated with the toxicity profile of new drugs.
Evidence of the relationship between glucocorticoids and CVDs comes primarily from studies of associations with current baseline medication use or dose[2, 3, 7, 11-16], ignoring the doses previously administered and their changes over time, as well as the concomitant use of other common medications that can affect the risk of CVDs (e.g. non-steroidal anti-inflammatory drugs). Many also have failed to adjust for important cardiovascular risk factors, such as smoking[1, 2, 7, 13, 16]. These studies have reported a dose-dependent risk of CVD with weaker associations for daily prednisolone-equivalent doses lower than 5-10mg[2, 3, 13, 17].
Our study aimed to estimate daily and cumulative dose-dependent oral glucocorticoid cardiovascular disease risk accurately in people diagnosed with six common immune-mediated inflammatory diseases in England.
METHODS
Setting and data sources
We analysed linked electronic health records from people registered at family practices in the Clinical Practice Research Datalink (CPRD) between 1st Jan 1998 and 15th March 2017. CPRD contains demographic and lifestyle data, diagnoses (e.g. stroke), prescribed medication and results of laboratory tests and clinical examinations, prospectively recorded during primary care contacts[18].
Previous validation studies have provided evidence of the accuracy of diagnostic and prescribing data[18]. Patients are broadly representative of the UK population with regard to age, sex and ethnicity[18]. CPRD data were linked to hospital records and the mortality registry (Text 1 in S1 File). Hospital records from the Hospital Episode Statistics (HES) (www.hscic.gov.uk/hes) contain diagnoses recorded during elective and emergency hospital admission across all National Health Service hospitals in England. Mortality data from the Office of National Statistics (ONS) (https://www.ons.gov.uk/atoz?query=mortality&size=10) were used to identify dates and causes of death.
Ethical considerations
The study was approved by the Independent Scientific Advisory Committee for Medicines and Healthcare products Regulatory Agency database research (ISAC), reference 16_146.
Study design and follow-up
This was a cohort study including all patients continuously registered in a CPRD practice for 1 year or more, aged ≥18 years and free of CVD, who had been diagnosed with at least one of 6 immune-mediated inflammatory diseases commonly treated with oral glucocorticoids at or before the start of follow-up. These were polymyalgia rheumatica, giant cell arteritis (GCA), systemic lupus erythematosus, rheumatoid arthritis, vasculitis, and inflammatory bowel disease (Figure 1 in S1 File). Diagnostic codes used to identify patients with each immune-mediated inflammatory disease are shown in Table 1 in S1 File. For each patient, the follow-up started when they first became eligible. It ended on the earliest of the following dates: occurrence of the outcome analysed (e.g. stroke), leaving the family practice, death, or last data collection date.
Note: AF, atrial fibrillation; CI, confidence interval; HF, heart failure; HR, Hazard ratios from Cox proportional imputed models adjusted for baseline age, sex, index of multiple deprivation, smoking status, ethnicity, body mass index, type of immune-mediated inflammatory disease, comorbidities (diabetes, diagnosed hypertension, cancer, asthma, chronic obstructive pulmonary disease, and renal disease), biomarkers (total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, c-reactive protein, creatinine), number of hospital admissions in last year, and prescribed non-oral glucocorticoids; and time-variant use of disease-modifying anti-rheumatic drugs and non-steroidal anti-inflammatory drugs; the practice identifier was included as a random intercept to account for clustering effect.
Oral glucocorticoid exposure
For each prescription of oral glucocorticoids issued to the patients between one year before the start and the end of the follow-up dates, we derived the daily dose from the recorded product name, which included information on product strength (e.g. 2mg), directions given (e.g. one tablet once a day) and quantity prescribed (e.g. 28 tablets). We then estimated the duration of each oral glucocorticoid prescription dividing the quantity of tablets prescribed by the daily dose. Given the variation in relative anti-inflammatory effects of different types of glucocorticoids, for each prescription we finally converted the daily dosage into milligrams of prednisolone-equivalent dose (Table 2 in S1 File).
Associations between time variant oral glucocorticoid prednisolone-equivalent dose and incident all-cause cardiovascular diseases by immune-mediated inflammatory disease
We defined several time-variant glucocorticoid variables to quantify current and cumulative drug exposure: i) ever use from one year prior to follow-up start (binary variable); ii) current daily use (i.e. whether or not the patient was prescribed glucocorticoids at a given time point [binary variable]); iii) current daily dose per 5 mg/day with zero value when medication was not prescribed (continuous and categorical variables: non-use, 1-4.9 mg, 5.0-14.9 mg, 15.0-24.9 mg, ≥ 25.0 mg/day); iv) cumulative dose since one year prior to follow-up start per 1,000 mg (i.e. sum of the total dose prescribed up to that point divided by 1,000; considered as continuous and categorical variables: non-use, 1-959 mg, 960-3,054 mg, 3,055-7,299 mg, and ≥7,300 mg; as defined previously[19, 20]).
Outcome measures
The primary outcome was the first occurrence of a composite of fatal and non-fatal cardiovascular diseases (all-cause CVD). Secondary outcomes were the first occurrence of the following common types of CVDs: atrial fibrillation, heart failure, myocardial infarction, cerebrovascular disease, peripheral arterial disease and abdominal aortic aneurysm. Diagnostic codes[21] used to define the outcomes are listed in Table 3 in S1 File and have been validated and used in multiple previous studies[22-26].
Confounding variables
We considered the following variables as a priory confounders: baseline age, sex, ethnicity, socioeconomic status (index of multiple deprivation[24, 27], area-based indicator linked through the patient’s home postcode), smoking status, body mass index (BMI), biomarkers (total, high and low density lipoprotein-cholesterol, systolic blood pressure, c-protein reactive protein, creatinine), underlying disease (e.g. rheumatoid arthritis), comorbidities recorded in primary or hospital care (diabetes, diagnosed hypertension, cancer, asthma, chronic obstructive pulmonary disease [COPD], renal disease); prescribed non-oral glucocorticoid medication (inhaled, nasal, parenteral/intra-articular, topical and rectal), and the number of hospital visits 1 year before baseline. We also considered time-variant prescribed medication (disease-modifying antirheumatic drugs and non-steroidal anti-inflammatory drugs) during follow-up. Detailed definition of covariates are shown in Text 1 in S1 File
Statistical analysis
We replaced missing daily dose of oral glucocorticoids (i.e. during tapering periods) and confounders through multiple imputation with chained equations with generation of 25 datasets (Text 1 in S1 File). Models for dose imputation included patient demographics (i.e. age, sex, and index of multiple deprivation[28]), underlying immune-mediated inflammatory disease, time between follow-up start and prescription, type of oral glucocorticoid (e.g. prednisone) and diagnosed comorbidities (e.g. diabetes).
We used standard descriptive statistics to describe baseline patient characteristics. We estimated cumulative probabilities of CVD outcomes using Kaplan-Meier methods. We calculated incident rates with 95% confidence intervals (95% CI) dividing the number of patients with incident CVD by the total number of person-years of follow-up.
We assessed the association between the outcomes and each of the oral glucocorticoid exposures using Cox proportional hazards models adjusted for the a priori confounders, with the practice identifier included as a random intercept to account for clustering effect. No interaction terms were included. The proportional hazards assumption was assessed using Schoenfeld residuals tests. The primary analysis was based on covariate imputed data. We generated models for each of the 25 imputed datasets and pooled estimates and accompanying 95% CIs following Rubin’s rules. We used 2-sided tests and considered significant at p<0.05. We performed the data management in Stata (StataCorp LP, College Station, USA; version 15) and analyses in R (http://cran.r-project.org/; version 3.3.1).
In secondary analyses, we modelled cardiovascular risk separately for men and women, for each of the 6 immune-mediated inflammatory diseases studied, and according to duration of these diseases at the start of follow-up (newly diagnosed/incident, within 2 years and over 2 years since diagnosis).
In sensitivity analyses, we obtained estimates from complete case models (i.e. restricted to patients with complete covariate data), from models including covariates with a separate category for missing data and from models unadjusted for biomarker data with a level of missingness >60%. We also additionally adjusted the models for the level of disease activity. We defined periods of active disease based on c-reactive protein and erythrocyte sedimentation rate levels (≥10 mg/mL and ≥30 mm/h, respectively) and the glucocorticoid daily dose (increase in prednisolone-equivalent dose by >5 or 10 mg that was sustained for over 3 weeks) (Text 1 in S1 File).
RESULTS
Patient characteristics
The study included 87,794 adults from 389 general practices with at least one immune-mediated inflammatory disease diagnosed; 25,581 with polymyalgia and/or giant cell arteritis, 27,739 (31.6%) with inflammatory bowel disease, 25,324 (28.8%) with rheumatoid arthritis, 5,199 (5.9%) with vasculitis, and 3,951 (4.5%) with systemic lupus erythematosus (Table 1). The overall mean age was 56 years (SD 18.3), 29,935 (34.1%) were men and 21,264 (24.2%) were current smokers. At baseline, the mean duration since immune-mediated disease diagnosis was 9.6 years (SD=8.7; range from 6.8 years for polymyalgia and/or giant cell arteritis and 11.7 years for inflammatory bowel disease). The most common patient comorbidities were hypertension (25.1%), asthma (14.6%) and diabetes (6.4%).
In the year prior to follow-up start, 14,356 (16.4%) patients were prescribed oral glucocorticoids, 13,331 (15.2%) inhaled or nasal glucocorticoids and 39,690 (45.2%) non-steroidal anti-inflammatory drugs. During follow-up 18,877 (21.5%; range 4.2% for polymyalgia and/or giant cell arteritis to 48.0% for rheumatoid arthritis) received disease-modifying antirheumatic drugs.
Incidence of fatal and non-fatal cardiovascular diseases
The median time of follow-up per patient was 5.0 (IQR 2.0-6.2) years. A total of 13,426 incident cardiovascular events occurred (15.3% of patients) over 541,655 person-years of follow-up (Table 4 in S1 File), including 6,013 episodes of atrial fibrillation, 4,727 of heart failure and 2,809 of acute myocardial infarction. The incidence of all-cause CVD was 24.8 per 1,000 person-years (95%CI 24.4-25.2). It increased from 18.5 (95%CI 18.1-18.9) per 1,000 person-years for periods of non-glucocorticoid use to 45.6 (95%CI 42.1-49.2) for periods of ≥25mg daily dose; and from 19.9 (95%CI 19.3-20.5) for unexposed periods to 26.4 (95%CI 25.5-27.2) per 1,000 person-years for ≥7,300 mg cumulative those). A total of 7,940 cardiovascular events happened during periods of non-exposure.
Cumulative probabilities of cardiovascular diseases
The Kaplan-Meier estimates of all-cause CVD at 1 year increased from 1.5% (95%CI 1.4%-1.6%) for periods of non-use, through 3.8% (95%CI 3.4%-4.3%) for <5 mg, to 9.1% (95%CI 7.6%-10.7%) for ≥25.0 mg daily dose; and from 7.6% (95%CI 7.4%-7.8%) for unexposed periods to 10.3% (95%CI 9.7%-10.9%) for ≥7,300 mg cumulative dose (Table 5 in S1 File). We found higher dose-response estimates in men than in women, in women with systemic lupus erythematosus, men with rheumatoid arthritis, and for atrial fibrillation and heart failure (Tables 6&7 in S1 File).
Relationship between glucocorticoid dose and cardiovascular diseases
The increase in the hazard of all-cause CVD per 5 mg increase in daily dose was 1.08 (95%CI 1.07-1.10 per 5 mg/day), ranging from 1.07 (95%CI 1.06-1.09) for inflammatory bowel disease to 1.30 (95%CI 1.22-1.38) for systemic lupus erythematosus (Table 2). We found strong dose-response estimates for current daily doses of <5.0mg for all immune-mediated diseases (HR=1.74, 95%CI 1.64-1.84; range 1.52 for polymyalgia and/or giant cell arteritis to 2.82 for systemic lupus erythematosus), for all cardiovascular outcomes and for daily and cumulative dose (Figures 1-3; Figures 2-13 in S1 File). The highest glucocorticoid dose-response estimates were for heart failure and for acute myocardial infarction. We found similar patterns in sensitivity analyses, including restriction to patients with complete covariate data (Tables 8&10 in S1 File). Daily and cumulative dose-response estimates were generally higher amongst patients with longer underlying inflammatory disease duration and in those newly diagnosed (Tables 11&14 in S1 File). Further, adjustment for the level of disease activity generally decreased the dose-response estimates, but associations remained statistically significant (Tables 15&16 in S1 File).
Note: AMI, acute myocardial infarction; CI, confidence interval; HR, Hazard ratios from Cox proportional imputed models adjusted for baseline age, sex, index of multiple deprivation, smoking status, ethnicity, body mass index, type of immune-mediated inflammatory disease, comorbidities (diabetes, diagnosed hypertension, cancer, asthma, chronic obstructive pulmonary disease, and renal disease), biomarkers (total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, c-reactive protein, creatinine), number of hospital admissions in last year, and prescribed non-oral glucocorticoids; and time-variant use of disease-modifying anti-rheumatic drugs and non-steroidal anti-inflammatory drugs; the practice identifier was included as a random intercept to account for clustering effect. PAD, peripheral arterial disease.
Note: AAA, abdominal aortic aneurysm; CI, confidence interval; CVA, cerebrovascular disease; HR, Hazard ratios from Cox proportional imputed models adjusted for baseline age, sex, index of multiple deprivation, smoking status, ethnicity, body mass index, type of immune-mediated inflammatory disease, comorbidities (diabetes, diagnosed hypertension, cancer, asthma, chronic obstructive pulmonary disease, and renal disease), biomarkers (total cholesterol, high-density lipoprotein cholesterol, low-density lipoprotein cholesterol, c-reactive protein, creatinine), number of hospital admissions in last year, and prescribed non-oral glucocorticoids; and time-variant use of disease-modifying anti-rheumatic drugs and non-steroidal anti-inflammatory drugs; the practice identifier was included as a random intercept to account for clustering effect.
DISCUSSION
In this longitudinal study of 87,794 adults diagnosed with at least one of six common immune-mediated inflammatory diseases, we quantified oral glucocorticoid dose-dependent risks of all-cause and type-specific CVDs taking into account changes in prescribed medication over time. At 1 year, the cumulative risk of all-cause CVD increased from 1.5% during periods without medication, through 3.8% for a daily prednisolone-equivalent dose <5mg, to 9.1% for periods with a daily dose of ≥25.0mg. We found strong dose-dependent increases in hazards of all-cause CVD, atherosclerotic diseases, heart failure, atrial fibrillation and abdominal aortic aneurysm, regardless of the underlying immune-mediated disease, its activity and duration. The cardiovascular risk profile of the study patients showed high prevalence of modifiable risk factors, including current smoking (24.2% of patients), BMI ≥30 kg/m2 (24.5%) and hypertension (25.1%).
Previous studies have reported increased risk of composites of cardiovascular[2, 3, 11, 12] or coronary heart disease[1], myocardial infarction[2, 3, 11, 14, 15], heart failure[2, 3, 11], stroke[2, 3, 16] and atrial fibrillation[7, 13] in current glucocorticoid users. Some found an increased risk of CVD only for daily doses of 5-10mg or higher[13, 16, 17]. Estimates from previous studies were based on current, baseline medication use[2, 3, 7, 11-16] or dose, or average glucocorticoid dose in the last 6-12 months[1, 3], without consideration of previously administered doses and changes in dose or medication use over time. Consistent with our findings, other studies assessing the relationship between glucocorticoid use and the risk of different types of CVDs reported stronger associations for heart failure than for other cardiovascular outcomes[2, 3].
The high absolute risk of CVD in patients receiving high doses of glucocorticoids, of similar magnitude to that of patients with diabetes or established CVD, warrants the need to implement and evaluate intensive lifestyle modification interventions to this high-risk group. The dose-dependent increased risk of CVDs, including atherosclerotic diseases, heart failure and atrial fibrillation, observed in our study, supports the need for close monitoring of cardiovascular risk in patients diagnosed with immune-mediated inflammatory diseases during glucocorticoid treatment and in the period after therapy discontinuation. Of cardiovascular risk scores currently used to guide decision on when to start primary cardiovascular prophylaxis, only the QRISK3[29] considers whether the patient is currently taking glucocorticoids (as a binary ‘yes/no’ predictor that ignores dose and recent exposure) and whether he/she is diagnosed with rheumatoid arthritis or systemic lupus erythematosus. Further refinements of this risk prediction tool, taking into account cumulative and/or current dose, might therefore improve its performance to identify patients in need for primary cardiovascular prevention. Our findings also emphasise the importance of rapid glucocorticoid dose tapering and discontinuation as soon as disease control is achieved, as well as the importance of evaluating the safety profile of alternative therapeutic options for patients with autoimmune-mediated inflammatory diseases.
The estimation of drug dose-response risks in this population-based cohort of all people with immune-mediated inflammatory diseases with different levels of activity and duration minimised the introduction of selection bias and increased the generalisability of the results. The use of linked health data from primary care and hospital facilities and the mortality registry, and diagnostic codes extensively used and validated for cardiovascular research[23-26, 30, 31], increased ascertainment of all the outcomes assessed. Information on prescribed medication is prospectively collected and includes all prescriptions issued in primary care, where patients with immune-mediated diseases are primarily treated. We derived the dose of oral glucocorticoids and the duration of prescribed medication from the directions given to patients on how to take their treatment. During periods of dose tapering, when these directions were unspecific (e.g. written ‘as directed’), we used the longitudinal doses prescribed to the patients to impute the dose taken. The lack of data on hospital prescribed medication and on drug adherence is likely to have resulted in underestimation of the dose taken when specialists treated the patients and might have overestimated the dose taken in periods of low disease activity for some patients. The resulting misclassification is likely to have reduced the size of dose-response estimates. We minimised time-related bias through use of time-variant medication variables (both exposure and confounders) and a start of patient follow-up that was unrelated to the start or use of glucocorticoid therapy. Although the main purpose of the study was to provide estimates of oral glucocorticoid dose-response for patients with the inflammatory diseases studied without making aetiological inferences, we examined the effect of confounding by indication and ascertainment bias through adjustment by periods of disease activity and disease-modifying antirheumatic drugs use, and performing analyses according to duration of the underlying disease. Dose-response associations remained strong and statistically significant. We adjusted estimates of risk for established cardiovascular risk factors (e.g. smoking, hypertension and diabetes), concomitant use of medications (e.g. time-variant disease-modifying antirheumatic drugs, non-steroidal anti-inflammatory drugs, non-oral glucocorticoid use). In the primary analysis, we used multiple imputation to handle missing baseline biomarkers and smoking data. We found similar patterns of dose-response in sensitivity, including analyses in which we used a separate category for missing covariate data, in those restricted to individuals with complete covariate data and in those unadjusted for baseline biomarker information with high level of missingness (>60% of patients).
In conclusion, we reported improved estimates of dose-dependent risks of CVDs. Our findings highlight the importance of implementing and evaluating targeted intensive cardiovascular risk factor modification interventions; promptly and regularly monitor patient cardiovascular risk, beyond diagnosis of inflammatory arthropathies and systemic lupus erythematosus, even when prescribing low prednisolone-equivalent doses; and the need for refining existing risk prediction tools for primary prevention of CVDs. Furthermore, the estimates of risk can be used to conduct cost-effectiveness and benefit-harm evaluations that guide the introduction of newly licensed glucocorticoid-sparing drugs for the treatment of immune-mediated inflammatory diseases. These estimates are to be complemented by future work on the estimation of risk of cardiovascular events beyond the first occurrence of CVD that are considered in calculations of glucocorticoid-associated costs.
Data Availability
Access to raw data can be requested from the CPRD (https://cprd.com)
SUPPORTING INFORMATION
S1 File. Supplemental material
ACKNOWLEDGEMENTS
AWM is supported by the Medical Research Council TARGET Partnership Grant (Treatment According to Response in Giant cEll arTeritis) (MR/N011775/1), the National Institute for Health Research (NIHR) infrastructure @ Leeds, including the NIHR Biomedical Research Centre (IS-BRC-1215-20015), and NIHR Medtech and In vitro Diagnostics Co-operatives (MIC-2016-015). JW is supported by the NIHR infrastructure @ Leeds. The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. The study funders had no role in the study design, data collection, analysis or interpretation, in the writing of the paper or in the decision to submit the paper for publication.
REFERENCES
