A MODEL TO PREDICT CARDIOVASCULAR EVENTS IN PATIENTS WITH NEWLY DIAGNOSED WEGENER’S GRANULOMATOSIS AND MICROSCOPIC POLYANGIITIS (original) (raw)

. Author manuscript; available in PMC: 2015 Jul 22.

Published in final edited form as: Arthritis Care Res (Hoboken). 2011 Apr;63(4):588–596. doi: 10.1002/acr.20433

Abstract

Objectives

To create a prognostic tool to quantify the 5 year cardiovascular (CV) risk in patients with newly diagnosed Wegener’s granulomatosis (WG) and microscopic polyangiitis (MPA) without pre-morbid CV disease.

Methods

We reviewed CV outcomes during the long term follow up of patients in the first 4 European Vasculitis Study Group (EUVAS) trials of WG and MPA. CV events were defined as: CV-death, stroke, myocardial infarction, coronary artery bypass graft, or percutaneous coronary intervention. Logistic regression was performed to create a model to predict the absolute risk of a CV event. The model was tested using the Wegener’s Granulomatosis Etanercept Trial (WGET) cohort.

Results

74/535 (13.8%) of the patients with 5 years of follow up from the EUVAS trials had at least one CV event: 33/281 (11.7%) WG vs. 41/254 (16%) MPA. The independent determinants of CV outcomes were; older age [OR 1.45 (95%CI 1.11 – 1.90)]; diastolic hypertension [OR 1.97 (95%CI 0.98 – 3.95)], and positive PR3 ANCA status [OR 0.39 (95%CI 0.20 – 0.74)]. The model was validated using the WGET cohort (Area under ROC curve = 0.80).

Conclusion

Within 5 years of diagnosis of WG or MPA, 14% of patients will have a cardiovascular event. We have constructed and validated a tool to quantify the risk of a cardiovascular event based on age, diastolic hypertension and PR3 ANCA status in patients without prior CV disease. In patients with vasculitis, PR3 ANCA is associated with reduced cardiovascular risk compared to MPO ANCA or negative ANCA status.

Keywords: Vasculitis, Cardiovascular, Wegener’s granulomatosis, microscopic polyangiitis, ANCA, myocardial infarction, stroke

Introduction

Wegener’s granulomatosis (WG) and microscopic polyangiitis (MPA) are the two most common types of small and medium vessel vasculitis with an estimated combined prevalence of 49–254 per million in Europe and North America.[1] Untreated these diseases are fatal, but modern therapy has dramatically improved survival.[28] However, patients continue to experience long-term morbidity and mortality from persistent low grade activity and permanent damage caused by the acute phase of vasculitis or its treatment.[9] An important component of this increased long-term morbidity and mortality is from cardiovascular (CV) disease which was highlighted in a retrospective review of the Danish National Hospital Register, where patients with WG showed an increased rate of a myocardial infarction within the 1st 5 years after a diagnosis of WG compared to the general population (hazard ratio 3.6).[10] In addition, a retrospective study showed that patients with WG and MPA when matched for renal function and other traditional risk factors, had double the rate of CV events.[11] Theories for the increased CV events in vasculitis include systemic inflammation and endothelial dysfunction,[12, 13] factors associated with increased cardiovascular risk in other inflammatory diseases.[1416]

In 1995 the European Vasculitis Study Group (EUVAS) launched randomised controlled trials for the treatment of WG and MPA with unified trial protocols and data collection procedures; 4 of these trials have been published.[4, 68, 17] Patient enrolment in the 4 trials was stratified by the degree of renal involvement. Subsequent to the completion of these trials, EUVAS has performed a long-term follow-up study on patients enrolled in these original trials. Together, this cohort represents the largest ever prospectively studied group of patients with WG and MPA.

The aims of this study were to 1) review the cardiovascular events in the first 5 years from the long-term follow-up of the 4 EUVAS trials; 2) use this data to create a prognostic tool aimed at predicting the 5 year risk of cardiovascular events in patients with newly diagnosed WG and MPA who have no cardiovascular disease at diagnosis; 3) validate the tool in a second separate cohort of vasculitis patients.

Patients and Methods

The patients and methods of the first 4 EUVAS randomised controlled trials included in the long-term follow-up (LTFU) study have previously been described.[4, 68, 17] In summary, all patients had a new diagnosis of WG, MPA or renal limited vasculitis (RLV). Patients were stratified into the 4 studies by the severity of renal involvement, listed here in order: 1) MEPEX: Comparison of Plasma Exchange to Methyl prednisone in patients with severe renal disease (Cr >500mmol/L or requiring dialysis).[7] 2) CYCAZAREM: Maintenance therapy with azathioprine vs. cyclophosphamide in patients with a renal manifestation (Cr<500mmol/L) or with generalised life threatening disease.[4] 3) CYCLOPS: Induction therapy comparing daily oral cyclophosphamide with pulsed intravenous cyclophosphamide in patients with generalised disease (Creatinine <500mmol/L).[8] 4) NORAM: A comparison of methotrexate to cyclophosphamide for induction treatment in patients with a Cr<150mmol/L and without critical organ threatening disease.[6] The study methods, data collection procedures and disease scoring were consistent between the 4 trials.

For the validation set we used anonymised data from patients enrolled in the Wegener’s granulomatosis etanercept trial (WGET): 180 patients with active Wegener’s granulomatosis who were randomized to receive etanercept or placebo in addition to standard maintenance therapy.[18] Long-term follow-up of this cohort is now available including cardiovascular events. Data from 136 patients without pre-existing CV disease was available for testing our regression model. 16/136 patients in the WGET cohort had a CV event: 14 with a non fatal cardiac event (myocardial infarction, coronary artery bypass graft, or percutaneous coronary intervention), 1 non fatal stroke and 1 CV death. Baseline demographics for the EUVAS and WGET cohorts are shown in table 1.

Table 1.

Baseline demographics in the EUVAS long-term follow-up and WGET cohorts

EUVAS LTFU (n=535) WGET (n=180) P*
Age in years, mean (SD) 58 (14) 47 (16) <0.001
Male gender, n (%) 288 (54) 108 (60) 0.17
Wegener’s granulomatosis, n (%) 281 (53) 180 (100) <0.001
Microscopic polyangiitis, n (%) 254 (47) 0 (0) <0.001
New diagnosis of vasculitis, n (%) 535 (100) 80 (44) <0.001
Existing diagnoses of vasculitis, n (%) 0 (0) 100 (56) <0.001
Serum creatinine μmol/L, mean(SD) 341 (321) 153 (177) <0.001
BVAS v2 (new/worse) score, mean (SD) 16.9 (9.2) Not available -
BVAS/WG score, mean (SD) Not available 6.9 (3.4) -
MPO ANCA positive, n(%) 190 (35.5) 21 (12) <0.001
PR3 ANCA positive, n (%) 302 (56.5) 131 (73) <0.001
ESR, mean (SD) 76 (36) Not available -
CRP, mean (SD) 88 (139) Not available -

Ethics approval was obtained for each participating site as per local requirements. A questionnaire on vital status and CV events was completed on all eligible patients by clinical investigators at each participating site. Information was collected at 5 years from entry into the original study, and also at the last available follow-up. The databases containing baseline information for the 4 trials were merged with the long-term follow-up results. All patients with an entry diagnosis of RLV were considered to have MPA for this analysis.

The primary outcome measure of interest in this sub-analysis was a CV event within 5 years after enrolment into the original EUVAS trial. We defined a CV event as death from any CV cause, non-fatal stroke, non-fatal myocardial infarction (MI), and coronary artery bypasses graft (CABG) or percutaneous coronary intervention (PCI). The identification of cardiovascular death was based on the local investigator reporting death and cause of death. The cause of death was adjudicated by an independent panel. The observed CV death rate in the EUVAS trials was compared to the predicted CV death rate (adjusted for age and country of origin). The relevant European and Mexican CV death rates for the year 2002 based on the World Health Organisation Statistical Information System (WHOSIS) database (http://www.who.int/whosis/en/) are provided in Appendix A.

Determination of non fatal MI required evidence of ECG changes or cardiac enzyme elevation. Non fatal stroke was defined as focal neurological deficit present for at least 3 months (i.e. transient ischemic events were not included). All items recorded at baseline from entry into each trial such as patient demographics, clinical features, blood results, disease activity scores (Birmingham Vasculitis Activity Score (BVAS 1&2))[19] and disease damage scores (Vasculitis Damage Index (VDI))[20] including individual components of these instruments were available as predictor variables. Hypertension in the development cohort (EUVAS trials) was defined as a diastolic blood pressure > 95 mmHg at the time of entry into the trial which was attributable to active vasculitis (irrespective of whether the patient had a previous history of hypertension or were on antihypertensive treatment). This definition of hypertension was chosen because of how the data was captured in the original trials (based on a BVAS item). For the validation sample (WGET cohort), hypertension was defined as a diastolic blood pressure >95mmHg at trial entry irrespective of previous history of hypertension or antihypertensive therapy. Attribution of the diastolic hypertension to vasculitis in the WGET cohort was not included in the definition because it was not captured in the original trial.

A points based Cox Framingham cardiovascular risk prediction tool,[21] which uses body mass index instead of total and HDL cholesterol was used to test how well a general population risk model performed in patients with vasculitis. We simplified Systolic Blood Pressure to < 140 mmHg (0 points on the scoring system) versus ≥140 mmHg (3 points women, 2 points men) because of the way in which hypertension was recorded in the dataset.

We used two methods to validate our new regression model: 1) Bootstrapping - an internal validation technique used to obtain a bias-corrected estimate in the development sample.[22]; and 2) Testing the model in the validation cohort using receiver operator characteristics.

Data analysis

Data was analysed using the statistical software package Stata, release 10.1 (College Station, Texas, USA). Univariate logistic regression modelling was used to examine the association between baseline variables and cardiovascular events. Fractional polynomial regression modelling was used to model non-linear relationships for continuous variables. A multivariate logistic regression model was then fitted including all predictor variables, regardless whether or not they were statistically associated with a CV event in univariate analysis. Having fitted the full multivariate model, a backwards selection process was used to exclude variables that did not improve model fit. Likelihood ratio tests were used to compare model fit. The results of complete case analyses can be biased because the cumulative effect of missing data in several variables often leads to exclusion of a substantial proportion of the original sample, causing a loss of precision and power. [23] This bias can be overcome by using multiple imputation methods which we did using the ICE procedure in Stata.[2326] For logistic regression, 10 times as many outcome events as predictor variables are required to avoid model over fit. Using clinical judgment, we initially created a model using traditional risk factors (age, gender, body mass index, smoking status, hypertension and diabetes), and then a model with only disease specific items (Birmingham Vasculitis Activity Score (BVAS), proteinase 3 (PR3) and myeloperoxidase (MPO) anti neutrophil cytoplasm antibody (ANCA) status, baseline estimated glomerular filtration rate (eGFR), age and gender). Regression diagnostics were performed to identify outlying data that may overly influence the model. Cholesterol measurements were not recorded at baseline therefore could not be included in the traditional risk factor model.

Performance of the predictive model was assessed in terms of calibration and discrimination. Calibration measures how closely predicted risk agrees with observed risk. This was assessed for each tenth of predicted risk ensuring 10 equally sized groups, and a Hosmer-Lemeshow goodness of fit test performed.[27] Discrimination is the ability of the model to differentiate between patients who experienced a CV event during the 5 years of follow-up in this study and those who did not. We tested discrimination by calculating the area under the receiver operating characteristic (AUROC) curve in a) the original dataset, b) by bootstrapping and c) testing in the WGET cohort.[18] To assess the goodness of fit of a logistic model we used McKelvey-Zavoina pseudo R2.[28, 29]

In addition, CART (Classification and Regression Tree) analysis, a binary recursive partitioning method was performed in the R statistical software package (R Foundation, Vienna, Austria), using the ‘party’ package. In each node of the tree a significance test was made between any of the covariates and the response, and a split established when the p-value was < 0.05.

Results

A total of 535 patients (281 WG, 254 MPA) including patients who died, had 5 years of follow-up; 74/535 (13.8%) of the patients had at least one CV event; 33/281 (11.7%) WG vs. 41/254 (16.1%) MPA. There were 32 (6.0%) CV deaths, 25 (4.7%) non fatal strokes and 42 (7.9%) had a non fatal myocardial infarction, coronary artery bypass graft or a percutaneous coronary intervention. The observed age standardised CV death rate for the EUVAS cohort was 699 compared to the predicted rate of 190 per 100,000 population per annum.

When developing the predictive model, 31 patients were excluded due to missing baseline values; 32 due to missing outcome records, and 45 due to CV disease prior to trial entry. The remaining 427 patients, 50 with a CV event (19/237 WG vs. 31/190 MPA, P<0.01) were used to develop the model. A summary of the baseline variables analysed for the 427 patients is provided in table 2. Smoking status was available in 292/427 patients in the study: 7.3% were current smokers, 23.7% ex-smokers and 37.5% had never smoked. For analysis we combined the current and ex-smokers as an ever smoked group. Comparing those that had never smoked to ever smoked gave an OR 0.61 (95%CI 0.32 – 1.18), P = 0.14, for any CV event within 5 years. Re-analysing the data comparing never smoked to current smokers gave OR 0.90 (95%CI 0.24 – 3.36), P = 0.87.

Table 2.

Baseline variables for patients included in the development of the cardiovascular risk prediction model. (N=427)

Total number of patients Number of patients with at least one cardiovascular event (%) Multiple Imputation [n=427]Crude OR (95%CI) P-value
Trial name:
CYCAZAREM 131 14 (10.7%) 1.00
CYCLOPS 107 14 (13.1%) 1.26 (0.57 – 2.77) 0.57
MEPEX 99 18 (18.2%) 1.86 (0.87 – 3.95) 0.11
NORAM 90 4 (4.4%) 0.39 (0.12 – 1.22) 0.11
Diagnosis:
MPA 190 31 (16.3%) 1.00
WG 237 19 (8.0%) 0.45 (0.24 – 0.82) 0.009
ANCA PR3:
No 168 33 (19.6%) 1.00
Yes 257 17 (6.6%) 0.29 (0.16 – 0.54) <0.001
ANCA MPO:
No 279 22 (7.9%) 1.00
Yes 146 28 (19.2%) 2.78 (1.52 – 5.06) 0.001
Age*: 1.58 (1.22 – 2.06) 0.001
Sex:
Female 207 23 (11.1%) 1.00
Male 220 27 (12.3%) 1.12 (0.62 – 2.02) 0.71
BMI: 1.03 (0.93 – 1.14) 0.57
Ever smoked:
No 160 15 (9.4%) 0.61 (0.32 – 1.18) 0.14
Yes 132 20 (15.2%) 1.00
Diabetes:
No 339 40 (11.8%) 1.00
Yes 14 3 (21.4%) 2.06 (0.57 – 7.40) 0.27
New hypertension:
No 350 33 (9.4%) 1.00
Yes 72 16 (22.2%) 2.74 (1.41 – 5.30) 0.003
Haemoglobin: 426 0.86 (0.75 – 1.00) 0.051
log WBC*: 427 1.20 (0.56 – 2.59) 0.63
Platelets**: 427 0.97 (0.95 – 0.99) 0.011
eGFR* (MDRD): 426 0.83 (0.74 – 0.93) 0.001
BVAS score: 422 1.02 (0.99 – 1.05) 0.23
VDI 6-month score: 238 1.38 (1.09 – 1.74) 0.007

All of the traditional and disease-specific risk factors were considered for inclusion in the final model. To avoid over-fitting we considered the effect of each predictor independently, and assessed the effect of confounding for that predictor with each risk factor. Variables that were significantly associated with cardiovascular events in univariate analysis, or were important confounders were included in the full multivariable regression model. Older age was associated with an increased risk of a CV event, however the effect was slightly attenuated due to confounding by eGFR. Crude analyses suggested that increasing eGFR was associated with reduced CV risk, but this was attenuated after adjustment for age and PR3 ANCA. As MPO ANCA and PR3 ANCA were strongly inversely associated with each other, only PR3 ANCA was considered for inclusion in the full model because it was a statistically stronger predictor. PR3 ANCA was associated with reduced CV risk and MPO with an increased CV risk, in both cases this was attenuated slightly after adjusting for age, eGFR and diastolic hypertension. Patients with diastolic hypertension were at increased risk of CV events. This effect was attenuated after adjusting for age, eGFR and PR3 ANCA, but strengthened if adjusting for gender and smoking status. Gender, BMI, smoking status, baseline BVAS score and diabetes were not associated with cardiovascular events.

Age, gender, eGFR, PR3 ANCA status, smoking status, and diastolic hypertension at baseline, were included in creating the full multivariable regression model. After backward selection: age, PR3 ANCA status and diastolic hypertension remained in the final model. The final model is shown in Figure 1. A cut off of <10% as low risk, 10–20% as moderate risk and >20% as high risk were chosen to provide memorable cut-offs. In the EUVAS cohort: 245 (57.4%) of patients fell into the low risk, 106 (24.8%) into the moderate risk, and 76 (17.8%) into the high risk group using our new 5 year risk prediction model. Observed CV events were: 15/245 (6.1%) for the low risk group, 14/106 (13.2%) for the medium risk group and 21/76 (27.6%) for the high risk group.

Figure 1.

Figure 1

Final model to predict the risk of a cardiovascular event in the first 5 years from diagnosis with microscopic polyangiitis or Wegener’s granulomatosis

Five year CV risk is shown as: green = low risk (<10%), orange = moderate risk (10–20%), Red = high risk (>20%).

The Hosmer-Lemeshow goodness of fit test suggested that the model was well calibrated (P = 0.55). The model demonstrated good discrimination with an AUROC of 0.73 for the original dataset; bootstrapped analysis gave a bias-corrected AUROC of 0.72; and validation with the WGET cohort demonstrated very good discrimination with an area under the ROC of 0.80; Figure 2. The pseudo R2 (McKelvey-Zavonia) for the final model was 18.5%. The model was better at predicting coronary events than strokes; AUROC 0.85 (95%CI 0.77 – 0.92) with a pseudo R2 of 49.7% vs. AUROC 0.70 (95%CI 0.60 – 0.80) and a pseudo R2 of 19.6% respectively for the EUVAS cohort. In comparison, the AUROC for any cardiac event using the points based Cox Framingham model was 0.65 (95%CI 0.57 – 0.72), Figure 2.

Figure 2.

Figure 2

Receiver operating characteristic curves comparing our new disease specific model with the traditional points based Cox-Framingham model

Figure 3 shows the conditional inference trees we created using CART analysis to predict CV events. In this tree, PR3 ANCA status was the most discriminative starting point. Thereafter, age, eGFR, and BVAS items were important predictors of the risk of a CV event. Using this inference tree, a patient that was PR3 ANCA positive, had eGFR > 32.3, and no items present on the Cardiovascular section of the BVAS, had the lowest risk (1.4%) of a CV event in the following 5 years. In the EUVAS cohort, 146 individuals fitted these criteria and 2 of them (1.4%) had a CV event. Conversely, 12 of 47 (44.4%) of patients that were PR3 ANCA negative and older than 72.4 years had a CV event.

Figure 3.

Figure 3

Conditional inference tree showing approximate risk of a cardiovascular event within 5 years based on baseline clinical features

% shown is the risk of a cardiovascular event within 5 years of diagnosis of Wegener’s granulomatosis or microscopic polyangiitis.

N=Cardiovascular events/number of individuals in Group

Discussion

Within 5 years of diagnosis of WG or MPA, approximately 14% of patients will have a major cardiovascular event. The age standardised annual CV mortality rate of 699 per 100,000 in this cohort is 3.7 (95% CI 3.2–4.3) times higher than we would expect in the background population. The result of our study adds further weight to the growing body of evidence that individuals with a diagnosis of WG or MPA are at significantly increased risk of cardiovascular morbidity and mortality.[10, 11]

There are many potential reasons for the increased CV risk. Endothelial dysfunction which is a recognised risk factor for CV disease[16, 3032] has been shown to be present in ANCA associated vasculitis and is independent of disease activity or renal involvement.[33] Renal dysfunction, which occurs frequently in patients with MPA and WG, is an established contributor to cardiovascular disease by affecting metabolic, inflammatory and hemodynamic pathways.[34] In addition, vasculitis represents a chronic inflammatory state, and other inflammatory diseases, such as rheumatoid arthritis and systemic lupus erythematosis have strong associations with CV disease.[35, 36] Animal models suggest that vessels following arteritis are more prone to atherosclerotic change.[37] Further, corticosteroids, which are a routine part of treatment for vasculitis, present an interesting conflict: they increase cardiovascular risk by accelerating the development of diabetes,[38] dyslipidemia[39] and hypertension,[40] but may also have a protective role in vasculitis by reducing systemic inflammation and improving endothelial dysfunction.[13, 15, 41]

We have developed 2 complementary models (logistic regression and CART) to predict the 5-year cardiovascular risk at the time of first presentation with WG and MPA. The regression model was validated in the WGET cohort and performed better than existing generic risk tools available for the general population. Both models indentified a positive PR3 ANCA status as an important determinant of lower risk among those with vasculitis. However, almost all patients in the EUVAS trials were either MPO or PR3 ANCA positive implying that those with a positive MPO ANCA are associated with an increased risk of CV disease. Interestingly, neither model identified gender as a predictor of risk, suggesting that vasculitis may remove the CV risk benefit usually observed in women. A potential explanation for this is that a large proportion of the women enrolled in the EUVAS trials were in the peri- or post-menopausal age group (mean age was 58, SD 14, i.e. 84% were over the age of 46 at baseline) and others may have been pushed into premature menopause by the use of cyclophosphamide. In the regression model, a positive PR3 ANCA status decreased the risk of a CV event, older age and the presence of diastolic hypertension at time of enrolment into the study increased the risk. The discrimination of this model compares very favourably with the QRISK score (area under ROC curve of 0.76 and 0.79 for males and females respectively) and the 1991 Framingham equation (area under ROC curve 0.74 for male and 0.76 for females).[4244] Our second model used CART analysis to create an inference tree to determine the cardiovascular risk which provides insight into possible interactions between variables. For example, the inference tree shows that patients who are PR3 ANCA positive and have good renal function are generally at low risk for cardiovascular events, however, if these patients have any cardiovascular system involvement at baseline then their risk of a major cardiovascular event is 9-fold higher.

The main reason for developing risk algorithms is not just prognostication, but prognostication that allows clinicians to develop and test preventative strategies. For vasculitis we do not yet know whether treating traditional risk factors such as smoking and dyslipidemia would change the cardiovascular outcome or whether it is more important to treat the inflammation and renal disease. Our model could be used to do the power calculation to determine how many patients would be needed in an interventional study to detect a 20% reduction in CV events. For example, you would need X high risk patients vs. Y low risk patients or Z all patients.

There are limitations to our study. Lipid results, glycosylated haemoglobin levels and specific cardiovascular medications associated with prevention of CV disease (e.g. aspirin, statins, ACE inhibitors) were absent from the baseline dataset which may have omitted important predictor variables. The original EUVAS trials were not designed to evaluate cardiovascular risk and therefore a study designed to look specifically at this issue may have resulted in an even stronger model. This is something to aspire to for future studies. However, ours is the only disease specific model available for vasculitis currently. When interpreting our statistical model it is important to remember that we are trying to predict a specific outcome and the key objective is to predict the risk accurately. There are potentially a large number of variables that have overlapping contributions to the cardiovascular risk but if a few variables are able to predict the risk with the same accuracy, then whether or not we include the other variables in the equation should not matter. For example ‘hypertension’ in our model would very likely have an overlapping contribution to cardiovascular risk with renal impairment. The patients included in our development dataset comprised a larger proportion of patients with mild WG (NORAM trial), and more patients with severe MPA (MEPEX trial), therefore patient selection bias may have unduly influenced the CV outcomes despite our best efforts to correct for renal function and disease severity in our statistical modelling. In addition, there are some major differences between the development and validation groups; all patients had newly diagnosed WG or MPA in the EUVAS trials whereas the WGET cohort only included WG but comprised patients with new and existing disease. We acknowledge that a large mixed cohort with both WG and MPA patients would have been the best cohort for testing and validation. However, there are currently no other suitable cohorts available (i.e. with sufficient number of patients, long enough duration of follow-up and with accurate recording of cardiovascular events). Therefore we have used the best alternative possible. There was only 1 stroke in the WGET cohort so validation of this outcome needs to be interpreted with caution. Despite these potential biases and cohort differences our model still works very well in both groups, confirming the strength of our model and reinforcing that age, diastolic hypertension and PR3 ANCA status are strong predictors of cardiovascular outcomes in patients with ANCA associated vasculitis. A further consideration is that renal function my change dramatically from baseline to later in the disease, therefore its use as a predictive variable may depend on when it is measured. The effect of disease flares, cumulative dose of steroids and changes in renal function will need to be taken into account when evaluating cardiovascular risk at different time points in the disease course.

In conclusion, we have shown that the risk of a cardiovascular event in the first 5 years after the diagnosis of WG or MPA is raised. To quantify this risk for an individual patient we have created and validated a statistical model using baseline clinical features. Identifying those at highest risk may help target those who require closer monitoring and further intervention.

Table 3.

Predictor variables included in developing the final regression model.

Complete case Multiple Imputation [n=427] Multiple Imputation [n=427] Multiple Imputation [n=427]
UnivariableOR (95%CI) UnivariableOR (95%CI) Pseudo R2 MultivariableOR (95%CI)R2 = 22.1% Reduced modelOR (95%CI)R2 = 18.5%
Age [n=427]* 1.58 (1.22, 2.06) 1.58 (1.22, 2.06) 11.7% 1.35 (1.02, 1.79) 1.45 (1.11, 1.90)
Gender [n=427]:
Female 1.00 1.00 1.00 -
Male 1.12 (0.62, 2.02) 1.12 (0.62, 2.02) 0.1% 1.10 (0.56, 2.15) -
BMI [n=196] 1.03 (0.93, 1.14) 1.03 (0.93, 1.14) 1.3% - -
BVAS1 [n=422] 1.02 (0.99, 1.05) 1.02 (0.99, 1.05) 1.0% - -
eGFR [n=426]* 0.83 (0.74, 0.93) 0.83 (0.74, 0.93) 11.3% 0.93 (0.82, 1.06) -
ANCA MPO [n=425]
No 1.00 1.00 - -
Yes 2.77 (1.52, 5.05) 2.78 (1.52, 5.06) 6.7% - -
ANCA PR3 [n=425]
No 1.00 1.00 1.00 1.00
Yes 0.29 (0.16, 0.54) 0.29 (0.16, 0.54) 10.0% 0.41 (0.21, 0.80) 0.39 (0.20, 0.74)
Ever smoked [n=427]:
Yes 1.00 1.00 1.00 -
No 0.58 (0.28, 1.18) 0.61 (0.32, 1.18) 2.3% 0.62 (0.29, 1.31) -
New Hypertension [n=422]:
No 1.00 1.00 1.00 1.00
Yes 2.74 (1.42, 5.32) 2.74 (1.41, 5.30) 4.0% 1.92 (0.93, 3.95) 1.97 (0.98, 3.95)
Previous diabetes [n=353]:
No 1.00 1.00 - -
Yes 2.04 (0.55, 7.62) 2.06 (0.57, 7.40) 0.9% - -

Acknowledgments

Funding

The long term follow up study of the EUVAS trials was funded by a project grant from the European League Against Rheumatism (EULAR). RS was supported by the Rose Hellaby Medical Scholarship, New Zealand. AJ, RB, MKJ and RL were supported by Oxford NIHR BRU Musculoskeletal Research Group, University of Oxford. DJ was supported by the Cambridge Biomedical Research Centre. The WGET was supported by grants from the National Institute of Arthritis and Musculoskeletal and Skin Diseases, the National Center for Research Resources, and the Office of Rare Diseases Research 1 U54-RR01949, U54-AR057319, and U01 AR1874, NIH N01-AR92240 and the Office of Orphan Products, FDA (grant FD-R-001652), General Clinical Research Center Grants M01-RRO-00533 (Boston University), M01-RRO-0042 (The University of Michigan), MO1-RR-30 (Duke University), and M01-RRO-2719 (Johns Hopkins University School of Medicine), from the National Center for Research Resources/NIH. Drs. Merkel, Stone, and St. Clair were supported by NIAMS grants K24 AR049185-01, K24 AR2224-01A1, and K24 AR02126-04.

APPENDIX A. A comparison of this EUVAS cohort’s age standardised cardiovascular mortality rate with selected European countries and Mexico

Country Age-standardized mortality rate for cardiovascular diseases (per 100 000 population)* Country of origin of patients in EUVAS cohort as percentage of whole group (n=535)
Belgium 162 3.9
Czech Republic 315 7.5
Denmark 182 1.7
Finland 201 2.6
France 118 4.7
Germany 211 10.7
Ireland 214 3.4
Italy 174 7.7
Lithuania 391 1.1
Mexico 163 0.9
Netherlands 171 8.2
Spain 137 6.2
Sweden 176 9.7
United Kingdom 182 31.6
Predicted CV mortality rate for EUVAS cohort adjusted for country of origin 190
Observed CV mortality rate for EUVAS cohort 699

APPENDIX B. The physicians involved in the EUVAS trials

Abramowicz D, Abuzakouk M, Adu D, Andrassy K, Bacon, P, Ballarin J, Bataille P, Blockmans D, Bruchfeld A, Burns A, Calero F, Chabova V, Chizzolini C, Cohen Tervaert JW, Confalonieri R, Dadoniene J, de Groot K, de Lind van Wijngaarden R, de Souza R, Dhaene M, Ekstrand A, Eriksson P, Esnault VL, Feehally J, Feighery C, Ferrario F, Flores Suarez L, Garibotto G, Gaskin G, Germanis G, Gibson J, Gregorini G, Grönhagen C, Gross WL, Guillevin L, Hagen EC, Harper L, Haubitz M, Heigl Z, Heimburger M, Hergesell O, Jayne DR, Kallenberg CGM, Kieley P, Leblau J, Lesavre P, Lundberg I, Luqmani R, Madhoun P, Mason P, Mathieson P, Mirapeix E, Natusch A, Nowack R, Olveira D, Oudejans I, Pettersson E, Poisetti G, Poveda R, Puechal X, Pusey C, Rasmussen N, Reinhold-Keller E, Riska, Rupprecht H, Rychlik I, Santostefano M, Savage C, Schmitt W, Schneider M, Segelmark M, Selga D, Sennesael J, Serra A, Siegert C, Sinico A, Specker C, Stahl-Hallengren C, Stegeman CM, Sterner G, Stevens J, Stolear J, Svenungsson E, Tesar V, Theander E, Tidman M, Tomson C, Turner N, Valles M, Valles M, van der Woude F, van Gurp E, Vanhille P, Verburgh C, Vischedyck M, Watts R, Weber P, Weidner S, Westman K, Wiik A, Williams A, Wissing M.

Appendix C. The WGET research group

WGET Chairman: Stone J (The Johns Hopkins Vasculitis Center), WGET Co-Chairman: Hoffman G (The Cleveland Clinic Foundation Center for Vasculitis Research and Care). The Johns Hopkins University Center for Clinical Trials (co-ordinating center): Holbrook JT, Meiner CL, Doge J, Donithan J, Min N, Murrow L, Smith J, Tibbs A, Van Natta M. The Beth Israel Medical Center, New York: Spiera R, Berman R, Enuha S. Boston University: Merkel PA, Gelbard R, Nuite M, Schiller A. The Cleveland Clinic Foundation: Hoffman G, Blumenthal D, Bork D, Clark T, Crook S, Calabrese L, Farkas S, Sridharan S, Strom K, Wilke. Duke University: St. Clair EW, Allen N, Rodin K, Scarlett E. Johns Hopkins University: Stone J, Hellman D, Moore A, Pinachos L, Regan M, Uhlfelder M. The Mayo Clinic: Specks U, Bradt K, Carlson K, Fisher S, Hammel B, Mieras K, Ytterberg S. University of California, San Francisco: Davis JC, Fitzpatrick M, Fye K, Lund S. University of Michigan: McCune J, Coomer B, Gilson B, Haftel H, Morrel-Samuels A, Neckel S. The Johns Hopkins University Immune Diseases Laboratory: Rose N, Burek CL, Barin J, Talor M. Data and Safety Monitoring Board: Canner P (Maryland Medical Research Institute), Conn DL (Emory University), Klippel JH (Arthritis Foundation), Landis JR (University of Pennsylvania).

Footnotes

Competing interests: none

References