Predictors for major cardiovascular outcomes in stable ischaemic heart disease (PREMAC): statistical analysis plan for data originating from the CLARICOR (clarithromycin for patients with stable coronary heart disease) trial
Diagnostic and Prognostic Research volume 1, Article number: 10 (2017)
The purpose of the predictors for major cardiovascular outcomes in stable ischaemic heart disease (PREMAC) study is exploratory and hypothesis generating. We want to identify biochemical quantities which—conditionally on the values of available standard demographic, anamnestic, and biochemical data—may improve the prediction of cardiovascular outcomes and/or death in patients suffering from stable ischaemic heart disease. The candidate biochemical quantities include N-terminal pro-B-type natriuretic peptide, YKL-40, osteoprotegerin, high-sensitive assay cardiac troponin T (hs-cTnT), pregnancy-associated plasma protein-A (PAPP-A), cathepsin B, cathepsin S, soluble TNF receptor 1 and 2, neutrophil gelatinase-associated lipocalin, endostatin, and calprotectin. As an extra objective, we also want to assess if skewness in these predictors may explain why the clarithromycin for patients with stable coronary heart disease (CLARICOR) trial found increased all-cause and cardiovascular (CV) mortality on a brief clarithromycin regimen compared with placebo.
Baseline data were obtained from the hospital files at five cardiology clinics covering the Copenhagen area. The CLARICOR trial included data from 4372 stable coronary artery disease patients recruited among such patients alive and diagnosed with acute myocardial infarction or unstable angina pectoris during 1993 to 1999 in Copenhagen and randomised during October 1999 to April 2000 to the CLARICOR trial of 14 days clarithromycin versus placebo.
Initial follow-up lasted for 2.6 years, during which outcomes were collected through hospital and death registries and assessed by an adjudication committee. Corresponding register data later showed to produce similar results. The adjudicated outcomes were therefore replaced and augmented by register data on outcomes to cover 10 years of follow-up. Biochemical marker data were obtained from analysis of serum from the CLARICOR bio-bank collected at randomisation and stored at −80° C.
Using Cox proportional hazard method, we will identify among the candidate biochemical quantities those which are significant predictors when used alone and in combination with the standard predictors as defined in the present study.
Patients who became stable during the period 1993 to 1999 and died before October 1999 are missing. The data from the placebo patients are nevertheless useful to identify new prognostic biomarkers in patients with stable coronary artery disease, and data from both trial groups are useful to assess important potential skewness between randomised groups. However, due to the potential selection bias, we do not feel that it is advisable to try to rank identified biochemical predictors relative to each other nor to use the results for predictive purposes.
Date of registration 13 July 2005
Date of enrolment of first participant 12 October 1999
Cardiovascular diseases, and ischaemic heart disease in particular, affect large fractions of the elderly population and constitute one of the two dominant causes of death [1,2,3,4,5]. Identifying high-risk patients would allow one to assess more aggressive measures for treatment of cardiovascular disease . The Prognosis Research Strategy (PROGRESS) series introduced a framework of four themes [7,8,9], including identification of specific factors (such as biomarkers and treatment modalities) that are associated with prognosis (prognostic factor research). This theme is the focus of the present “predictors for major cardiovascular outcomes in stable ischaemic heart disease” (PREMAC) study.
Objectives of the PREMAC study
Our primary objective is to identify biochemical predictors of a combined outcome including acute myocardial infarction (AMI), unstable angina pectoris (UAP), cerebro-vascular disease (CeVD), cardiovascular mortality (CV death), and all-cause mortality in non-hospitalised patients with stable coronary artery disease.
The patient group and data material
The patient group comprises the patients who participated in the clarithromycin for patients with stable coronary heart disease (CLARICOR) trial [10, 11], and the data material includes the baseline data collected at randomisation augmented by biomarker data obtained by analysing the bio-bank material collected at baseline and outcome data prospectively obtained from public registers .
The patient group is defined by the inclusion and the exclusion criteria of the patients who, during the winter 1999–2000, entered the CLARICOR trial.
The CLARICOR trial was an investigator-initiated, randomised, placebo-controlled, multicentre superiority trial including outpatients with stable coronary artery disease (CAD), using central 1:1 randomisation and blinding of all parties. All patients discharged from wards or outpatient clinics in the Copenhagen area were available in an existing database. We invited all 13,702 patients who were alive and aged 18–85 years in 1999 and identified with a diagnosis of myocardial infarction or unstable angina pectoris during the years 1993–1999 to visit one of five cardiology centres in the Copenhagen area. Six thousand one hundred sixteen (44.6%) patients accepted the invitation, and of these, 4372 (71.5%) were randomised, while 1567 (25.6%) were excluded, and 177 (2.9%) refused to participate. Exclusion criteria included AMI or UAP within the previous 3 months, percutaneous transluminal coronary angioplasty and coronary bypass surgery within the previous 6 months, impaired renal or hepatic function, congestive heart failure (New York Heart Association (NYHA) IV classification of heart failure), active malignancy, incapacity to manage own affairs, breast feeding, and possible pregnancy. Between October 1999 and April 2000, the 4372 patients were randomised to receive oral clarithromycin 500 mg once daily for 2 weeks versus matching placebo to assess the effects on the risk of major cardiovascular outcomes and death.
Data material collected from the 4372 patients include (1) demographic and anamnestic (hospital) data gathered prior to the randomisation, (2) values of biochemical quantities measured in plasma specimens obtained from the patients at randomisation, and (3) vital data and diagnostic information on first occurrence of cardiovascular outcomes covering the period from the start of the trial until December 31, 2009, and obtained from public registers.
Demographic and anamnestic data
Clinical data were obtained during enrolment interviews (smoking status, current medication, and known hypertension or diabetes), while information concerning sex, age, and history of myocardial infarction or unstable angina pectoris were obtained from the local hospital files.
Biochemical measurements on plasma collected at entry
Blood samples on trial participants were obtained just before randomisation and stored at −80° C. The following quantities were later measured: lipoproteins (total cholesterol, HDL cholesterol, LDL cholesterol, triglycerides, apoprotein A1, apoprotein B), high-sensitivity c-reactive protein (hs-CRP) , YKL-40 , high-sensitive assay cardiac troponin T (hs-cTnT) , pregnancy associated plasma protein-A (PAPP-A) , N-terminal pro-B-type natriuretic peptide , cathepsin B , endostatin , cathepsin S , soluble TNF receptor 1 and 2 (sTNFR1 and sTNFR2) , neutrophil gelatinase-associated lipocalin (NGAL) , calprotectin , osteoprotegerin , and glomerular filtration rate using creatinine (GFR) .
Vital data and diagnostic information on cardiovascular and other vascular outcomes
Vital status was monitored via the Danish Central Civil Register. Information about the underlying cause of death was obtained from the National Register of Causes of Death (RCD) . The Danish National Patient Register (NPR), covering all somatic hospital admissions, provided hospitalisation data . A blinded adjudication committee, consisting of three cardiologists, working in randomised rotation, assessed each hospital admission or death, using standard diagnostic criteria during the first 2.6 years of follow-up [10, 29]. These results were later compared to the results obtained when only registry data were used and similar conclusions were reached [29, 30]. The follow-up period was therefore extended to 10 years, and the analysis of the full 10 years period was solely based on registry data , as explained in the following.
The Danish 10-digit central person registration (CPR) number is used at all contacts with the health care system. Somatic hospital contact cannot be completed without a diagnosis based upon the International Statistical Classification of Diseases, 10th revision (ICD-10) and subsequent notification of the NPR. Each department must issue at least an action diagnosis (A diagnosis), describing the main reason for the admission. Other important diagnoses may be recorded as B diagnoses. All registers have coverage close to 100%. Based on this material, we used the following algorithm to transform the registry information into CV outcomes: for each A code in discharge notifications and, in case of death, each ‘underlying cause of death’ code (in the official terminology of the Registry of Causes of Death), we classified the outcome according to the ICD-10 coding system into a list of disjoint and exhaustive categories: acute myocardial infarction (AMI) (I21.0–23.9), unstable angina pectoris (UAP) (I20.0, I24.8–24.9), cerebro-vascular disease (CeVD) (I60.0–64.9 and G45.0–46.8), peripheral vascular disease (PVD) (I70.2–70.9), and non-cardiovascular disease (A00.0–T98.3 except the codes already covered).
For each patient and each of the outcomes AMI, UAP, CeVD, CV death, and all-cause mortality, we scanned the discharge codes in chronological order. We began with the codes of the first admission following randomisation and recorded the period from the date of randomisation until the date of the first occurrence of the outcome or (in case of a non-fatal outcome) until the date of death or until the date of censoring (December 31, 2009) whichever came first.
Classification of the predictors
We have chosen collectively to classify as ‘standard predictors’ the below mentioned groups of predictors, i.e. (1) clinical predictors, (2) current medical treatment, and (3) standard biochemical predictors. The term ‘standard predictors’ is only a collective term used by us in this particular study to refer to those baseline quantities that were available to us during the CLARICOR trial and which are either established prognostic predictors or proxies of such predictors not available to us [31, 32]. The biochemical predictors we have measured in addition for the present PREMAC study are labelled ‘advanced biochemical predictors’. The two groups of predictors, i.e. ‘standard predictors’ and ‘advanced biochemical predictors’ are defined below.
Sex, age, smoking history, history of myocardial infarction compared to angina only, hypertension, and diabetes.
Current medical treatment
Aspirin (Yes/No), beta-blocker (Yes/No), calcium-antagonist (Yes/No), ACE-inhibitor (Yes/No), long lasting nitrate (Yes/No), diuretic (Yes/No), digoxin (Yes/No), statin (Yes/No), and anti-arrhythmic drugs (Yes/No). The current medical treatment was included as proxy predictors because information about post-infarction heart failure and post-infarction angina pectoris are not available to us.
Standard biochemical predictors
High-sensitivity-reactive protein, glomerular filtration rate (GFR) estimated by creatinine, and lipoproteins (total cholesterol, HDL cholesterol, LDL cholesterol, triglycerides, apoprotein A1, and apoprotein B).
Advanced biochemical predictors
N-terminal pro-B-type natriuretic peptide, YKL-40, osteoprotegerin, hs-cTnT, PAPP-A, cathepsin B, cathepsin S, sTNFR1, sTNFR2, NGAL, endostatin, and calprotectin. Tests for all these quantities are commercially available.
Background information on the advanced biochemical predictors
Patients with coronary artery disease (CAD) are usually only evaluated using clinical variables. But some CAD patients have a high incidence of CV insults which are difficult to predict . Measurements of biomarkers may potentially help identifying CAD patients at high risk of such CV insults.
Previously, we found that increased serum N-terminal pro-B-type natriuretic peptide (NT-pro-BNP), a marker of left ventricular dysfunction, and heart failure was a stronger predictor of myocardial infarction (MI), CV death, and non-cardiovascular death than high-sensitive assay C-reactive protein (hs-CRP) in patients with CAD during a 2.6 year follow up . Similar results were found for YKL-40 in the same CAD patients . In this context, it is of pathogenetic interest that YKL-40 is expressed by arteriosclerotic plaque macrophages, particularly macrophages which have infiltrated into the lesion . The highest expression of YKL-40 is found in macrophages in the early atherosclerotic lesion. Another promising marker is high-sensitive assay cardiac troponin T (hs-cTnT) indicating myocardial ischaemia which when combined with the NT-pro-BNP results was found to be significantly associated with all-cause mortality, CV death, and MI after adjustment for traditional risk factors and NT-pro-BNP . Pregnancy-associated plasma protein-A (PAPP-A), a marker of vulnerable plaques in coronary arteries, has also been found to be predictive of CV insults and death in CAD patients . The glycoprotein osteoprotegerin (OPG), which is positively related to coronary calcification, vascular stiffness, and the presence of unstable atherosclerotic plaques , is included among the candidate predictors because we found it to be an independent predictor of mortality in CAD patients . The cathepsins are a group of proteinases that have been suggested to be causally involved in the different stages of the atherosclerotic process, from the early stages such as foam cell formation  to the later stages, such as destabilisation of the fibrous cap . Endostatin is an endogenous angiogenesis inhibitor where circulating levels have been suggested to mirror an increased neovascularisation induced by vascular or myocardial ischaemia , but endostatin has also been suggested to be a marker for an increased extracellular matrix remodelling . Inflammation is a key underlying factor in the atherosclerotic process  and tumour necrosis factor receptor alpha TNF-α, and its soluble receptors sTNFR1 and sTNFR2 are inflammatory markers that have been suggested to portray information about a systemic inflammatory state that is independent of other more established inflammatory markers such as CRP or IL-6 . Previous studies report higher levels of cathepsins, endostatin, and TNF-receptors in patients with atherosclerosis and/or CAD [41,42,43], but whether the proteins are relevant risk markers in these patients remains to be established. Calprotectin, neutrophil gelatinase-associated lipocalin (NGAL), and myeloperoxidase are all released from neutrophils when the cells are activated. Circulating levels of neutrophils and their activation products have been shown to be markers for plaque instability in both primary and secondary prevention of cardiovascular diseases [44, 45]. We have also previously shown that U-NGAL is associated with cardiovascular mortality . We now aim to explore if calprotectin is a risk marker in these patients.
Each of the abovementioned advanced biochemical predictors will be studied individually on the below mentioned outcomes provided that the advanced quantity concerned has shown promise in the primary (‘gatekeeper’) analysis. This primary ‘gatekeeper’ analysis will examine a composite outcome, defined to be present if at least one of the events AMI, UAP, CeVD, and all-cause death has occurred in a patient. The analyses will be conducted using the Cox proportional hazard model supplemented by Breslow estimation of the baseline hazard. We will use SAS 9.4, and the analysis will be based on the data from the placebo group because it was previously shown that clarithromycin had a significant and harmful effect on the prognosis of the trial patients [10,13,, 12–14].
The proportional hazards assumption
The joint proportional hazard property, covering all covariates included in a Cox analysis and the chosen functional forms of quantitative covariates, will be tested using cumulative sums of martingale-based residuals over follow-up time and/or covariate values . The test statistic is a Kolmogorow-type supremum test, using for each outcome a P threshold = 0.05 and Bonferroni adjusted P values.
Flow of analyses
For each advanced biochemical quantity, we will examine if the hazard ratio of the biochemical quantity is different from 1.00 with a P value <0.01 when used alone and when used in combination with the standard predictors (forced to be included in the analyses whether they are significant or not). Using centre ID as a class variable, we will stratify on centre ID implying that the observations are conditionally independent within centres and that the coefficients of the covariates are the same across centres.
For each outcome, we will test for interaction between sex and each clinical predictor, each standard biochemical predictor, and each advanced biochemical predictor using Bonferroni adjusted P values and 0.05 as a threshold. The interactions corresponding to adjusted P values below 0.05 will be included in the model. Explanation will be sought for any marked discrepancies between the two settings (candidate used alone and candidate used in combination with standard predictors); but we will let us be guided by the results obtained when the latter setting is used.
Those advanced biochemical quantities which have a hazard ratio different from 1.00 and a P value <0.01 when used in combination with the standard predictors in the initial ‘gatekeeper’ analysis, will be included in subsequent and similar analyses of each individual outcome of the composite outcome.
We will use the P values of the analyses as a data-reducing device where a threshold of P < 0.01 is used at each decision node on our way to identify a candidate biochemical predictor among the advanced biochemical predictors. We regard the set of selected advanced biochemical predictors as the primary result of this essentially hypothesis-generating effort. However, clinically implausible results will be discussed and commented on.
Illustration of predictive impact
Nine, six and three-year estimated incident risks will be calculated for each patient using the Cox-Breslow procedure. We will illustrate discrimination by the reclassification indices, reported separately for those with and those without the event in question as recommended by Kerr et al. . For example the risk estimates could be categorised as less than 25% (‘negative’) or 25% or more (‘positives’). With only two categories, this amounts to reporting the changes in the true-positive and false-positive rates.
For the mortality models, we present measures of risk differentiation, obtained by imagining preventive measures administered under budget constraints. The C statistic (Harrell’s concordance statistic)  represents the proportion of all patient pairs where the predicted survival is better for the patient who survived longer. The fixed budget constraint counts ‘lives saveable [by a hypothetical intervention in high-risk patients]’ .
There are virtually no missing data from entry information in the CLARICOR trial [10, 12]. We expect missing laboratory data to be MCAR due to the fact that we know that some serum samples are missing either due to missing blood sample in the first place or due to a vial being damaged (see discussion). We will use Little’s test to decide whether a multiple imputation or a complete case analysis should be conducted . Thus, each analysis of an advanced predictor will be initiated by Little’s test which will include all covariates with missing values (the standard biochemical predictors plus the advanced biochemical predictor to be assessed) plus all other variables in the model.
Additional supplementary analyses include (1) an analysis of the impact of alternative GFR definitions and (2) an assessment of whether chance imbalance of new biochemical predictors may dramatically change the original findings of the CLARICOR trial.
The creatinine-based GFR will be used as our quantitative measure of kidney function in this study. Within the framework of the model that we will develop, we will in a separate report compare the predictive power of different formulas to estimate GFR (creatinine-based GFR, cystatin C-based GFR, or a combined creatinine/cystatin GFR formula). We will divide participants in GFR categories (GFR >60, 30–60, and <30 mL/min/1.73 m2) according to each specific GFR equation. We will then study the impact of reclassification across the GFR categories by the different GFR equations, by assessing the risk of adverse outcomes in those who were reclassified to a higher GFR category or to a lower GFR category, compared to those who were not reclassified.
Finally, we will assess if skewness (chance imbalance) in any predictors may explain the increased cardiovascular and all-cause mortality in participants who received clarithromycin as compared with placebo participants, i.e. we will repeat the analysis of the effect of clarithromycin but this time include the biochemical prognostic markers identified in the above described study.
The present PREMAC study is important mainly for two reasons. First, new prognostic indicators are needed to improve preventive interventions in patients with ischaemic heart disease. Second, if known or new prognostic indicators show skewness between the two intervention groups at entry into the CLARICOR trial, this may help us better understand the harmful effects of clarithromycin. However, in previous publications we found no skewness of hs-CRP , N-terminal-pro-B-type natriuretic peptide , PAPP-A , YKL-40 , osteoprotegerin , and hs-cTnT .
The strengths of the CLARICOR trial are the considerable size of the patient population, long duration of follow-up, very few losses to follow-up (0.5%), the ethnic homogeneity of the patient population, and rarity of missing values. However, the values of 20 biochemical quantities have been added to the CLARICOR data and these data have missing values. The set of placebo patients with one or more missing biochemical quantity values includes 210 patients out of 2199 placebo treated patients giving 210/2199 = 9.5% patients having one or more biochemical quantity values missing. The percentage missing for the individual biochemical quantities ranges between 1.6 and 5.7%.
Potential weaknesses of the present cohort include the lack of information about left ventricle function, body mass index, blood pressure, and changes in medications during the follow-up period, and the fact that the prognosis of this type of patients today may have changed somewhat from the year 2000. Information about post-infarction heart failure and post-infarction angina pectoris was not available to us. Therefore, we added information about the medication at entry into CLARICOR as proxy information instead. The lack of information about left ventricular ejection fraction may be partially or completely compensated as Solomon et al.  found that age, sex, hypertension, prior AMI, creatinine, diuretics, digoxin, and mineralocorticoid receptor antagonist were related to left ventricular ejection fraction, all quantities that we have included within the group referred to as ‘standard predictors’.
Selection bias is also a potential weakness of the study. Only 2.6% of the patients selected for the CLARICOR trial refused to participate. So far that problem seems negligible. However, among those 7586 patients who declined our invitation to visit a cardiology centre, many must have been eligible for the CLARICOR trial, and we do not know how they looked and fared. So selection bias is a possibility.
The design of the CLARICOR trial may also contribute to selection bias. In the CLARICOR trial, the rate of cardiovascular events and deaths was studied in patients who had contracted an AMI or UAP during the period January 1993 to August 1999, were alive during the randomisation period (October 1999 to April 2000) and fulfilled the inclusion criteria and none of the exclusion criteria as defined in the CLARICOR trial .
Normally, to assess the rates in stable CAD patients, one would have started the study in January 1993 by including consecutively each patient who contracted an AMI or UAP, monitor the patient at regular intervals, say monthly from then on, note when the patient fulfilled the stability criteria (i.e. entered the stable state) and when the patient left the stable state again, and noting the characteristics of the new state (e.g. death). The patients included in such a design would differ from those of the CLARICOR trial in that they would include patients who entered the stable state and then died before October 1999. These patients are missing in the present PREMAC study. Thus, using the data to predict the prognosis of patients with stable CAD (as defined by the CLARICOR group), the results may very well be biased. This is the reason why we emphasise that the PREMAC study is an exploratory and hypothesis generating study.
We confined the net reclassification indices to the increase in the true-positive rate and the decrease in false-positives by including only two risk groups because when there are more than two risk categories, the indices do not adequately account for clinically important differences in shifts among risk categories . We did not use category-free net classification indices because they suffer from many of the same problems as similar measures such as, e.g. the area under the receiver operating characteristic curve does. Besides, it can mislead investigators by overstating the incremental value of a biomarker, even in independent data .
A single number summary of the prediction increment would be the improvement in net benefit. This would require us to know the ratio between the cost of falsely classifying an event as a non-event and the cost of classifying a non-event as an event. (If this ratio was 3:1, it would dictate a predicted risk of 25% as a decision threshold; the cut at 25% mentioned above was foreseen to be an illustrative cut and was not suggested to be a clinical decision threshold.) We do not know which value of this ratio will be appropriate; anyhow, we may decide to adapt the Vickers-Elkin net benefit diagram  to the PREMAC context.
Clearly, the prognosis for the occurrence of all-cause mortality or CV mortality as predicted by the baseline values may very well change once it becomes known that one of the non-fatal vascular outcomes has occurred. Thus, this valuable outcome information should be studied not only in its role as being an outcome as described above but also in its role as being a predictor of all-cause mortality and CV mortality. We plan to analyse these time-series aspects in a subsequent study using intervening MI and UAP as predictive events.
Acute myocardial infarction
Clarithromycin for patients with stable coronary heart disease
Central personal registration number
Glomerular filtration rate
High-sensitive assay cardiac troponin T
International Statistical Classification of Diseases
Neutrophil gelatinase-associated lipocalin
National Patient Register
New York Heart Association
Pregnancy-associated plasma protein-A
Predictors for major cardiovascular outcomes in stable ischaemic heart disease
Prognosis Research Strategy
Register of Causes of Death
Soluble TNF receptor 1
Soluble TNF receptor 2
Unstable angina pectoris
A 40 kDa heparin- and chitin-binding glycoprotein also known as human cartilage glycoprotein 39
Nichols M, Townsend N, Scarborough P, Rayner M. Cardiovascular disease in Europe 2014: epidemiological update. Eur Heart J. 2014;35:2929.
Cooper R, Cutler J, Desvigne-Nickens P, Fortmann SP, Friedman L, Havlik R, et al. Trends and disparities in coronary heart disease, stroke, and other cardiovascular diseases in the United States: findings of the national conference on cardiovascular disease prevention. Circulation. 2000;102:3137–47.
Schmidt M, Jacobsen JB, Lash TL, Botker HE, Sorensen HT. 25 year trends in first time hospitalisation for acute myocardial infarction, subsequent short and long term mortality, and the prognostic impact of sex and comorbidity: a Danish nationwide cohort study. BMJ. 2012;344:e356.
Wilson PWF, Cannon CP, Downey BC. Overview of the risk equivalents and established risk factors for cardiovascular disease. UpToDate. Post TW, Editors, UpToDate Waltham; 2015. http://www.uptodate.com/contents/Overview-of-the-risk-equivalents-and-established-risk-factors-for-cardiovascular-disease. Accessed 2 Mar 2017.
WHO. 2015. http://www.who.int/cardiovascular_diseases/en/. Accessed 2 Mar 2017.
Vilarino-Rico J, Pita-Fernandez S, Segura-Iglesias RJ. Clinical predictors of major adverse cardiovascular events during long-term follow-up after carotid endarterectomy. Ann Vasc Surg. 2015;29:419–25.
Collins GS, Reitsma JB, Altman DG, Moons KG. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ. 2015;350:g7594. doi:10.1136/bmj.g7594. Review.PubMedPMID: 25569120.
Hemingway H, Croft P, Perel P, Hayden JA, Abrams K, Timmis A, Briggs A, Udumyan R, Moons KG, Steyerberg EW, et al. Prognosis research strategy (PROGRESS) 1: a framework for researching clinical outcomes. BMJ. 2013;346:e5595.
Wynants L, Collins GS, Van Calster B. Key steps and common pittfalls in developing and validating risk models. BJOG. 2016. doi:10.1111/1471-0528.14170.
Jespersen CM, Als-Nielsen B, Damgaard M, Hansen JF, Hansen S, Helo OH, et al. Randomised placebo controlled multicentre trial to assess short term clarithromycin for patients with stable coronary heart disease: CLARICOR trial. BMJ. 2006;332:22–7.
Hansen S, Als-Nielsen B, Damgaard M, Helø OH, Petersen L, Jespersen CM, et al. Intervention with clarithromycin in patients with stable coronary heart disease: the CLARICOR trial design. Heart Drug. 2001;1:14–9.
Winkel P, Hilden J, Hansen JF, Kastrup J, Kolmos HJ, Kjoller E, et al. Clarithromycin for stable coronary heart disease increases all-cause and cardiovascular mortality and cerebrovascular morbidity over 10 years in the CLARICOR randomised, blinded clinical trial. Int J Cardiol. 2015;182:459–65.
Winkel P, Hilden J, Fischer Hansen J, Hildebrandt P, Kastrup J, Kolmos HJ, et al. Excess sudden cardiac deaths after short-term clarithromycin administration in the CLARICOR trial: why is this so, and why are statins protective? Cardiology. 2011;118:63–7.
Gluud C, Als-Nielsen B, Damgaard M, Fischer Hansen J, Hansen S, Helo OH, et al. Clarithromycin for 2 weeks for stable coronary heart disease: 6-year follow-up of the CLARICOR randomized trial and updated meta-analysis of antibiotics for coronary heart disease. Cardiology. 2008;111:280–7.
Harutyunyan MJ, Mathiasen AB, Winkel P, Gotze JP, Hansen JF, Hildebrandt P, et al. High-sensitivity C-reactive protein and N-terminal pro-B-type natriuretic peptide in patients with stable coronary artery disease: a prognostic study within the CLARICOR trial. Scand J Clin Lab Invest. 2011;71:52–62.
Kastrup J, Johansen JS, Winkel P, Hansen JF, Hildebrandt P, Jensen GB, et al. High serum YKL40 concentration is associated with cardiovascular and all-cause mortality in patients with stable coronary artery disease. Eur Heart J. 2009;30:1066–72.
Lyngbaek S, Winkel P, Gotze JP, Kastrup J, Gluud C, Kolmos HJ, et al. Risk stratification in stable coronary artery disease is possible at cardiac troponin levels below conventional detection and is improved by use of N-terminal pro-B-type natriuretic peptide. Eur J Prev Cardiol. 2014;21:1275–84.
Iversen KK, Teisner B, Winkel P, Gluud C, Kjoller E, Kolmos HJ, et al. Pregnancy associated plasma protein-A as a marker for myocardial infarction and death in patients with stable coronary artery disease: a prognostic study within the CLARICOR Trial. Atherosclerosis. 2011;214:203–8.
Wang N, Bai X, Jin B, Han W, Sun X, Chen X. The association of serum cathepsin-B concentration with age-related cardiovascular-renal subclinical state in a healthy Chinese population. Arch Gerontol Geriatr. 2016;65:146–55.
Ärnlöv J, Ruge T, Ingelsson E, Sundstrom J, Lind L. Serum endostatin and risk of mortality in the elderly: findings from 2 community-based cohorts. Arterioscler Thromb Vasc Biol. 2013;33:2689–95.
Jobs E, Ingelsson E, Risérus U, Nerpin E, Jobs M, Sundström J, et al. Association between serum cathepsin S and mortality in older adults. JAMA. 2011;306:1113–21.
Carlsson AC, Larsson TE, Helmersson-Karlqvist J, Larsson A, Lind L, Ärnlöv J. Soluble TNF receptors and kidney dysfunction in the elderly. J Am Soc Nephrol. 2014;25:1313–20.
Helmersson-Karlqvist J, Larsson A, Carlsson AC, Venge P, Sundstrom J, Ingelsson E, et al. Urinary neutrophil gelatinase-associated lipocalin (NGAL) is associated with mortality in a community-based cohort of older Swedish men. Atherosclerosis. 2013;227:408–13.
Larsen SB, Grove EL, Pareek M, Kristensen SD, Hvas AM. Calprotectin and platelet aggregation in patients with stable coronary artery disease. PLoS ONE. 2015;10(5):e0125992.
Bjerre M, Hilden J, Kastrup J, Skoog M, Hansen JF, Kolmos HJ, et al. Osteoprotegerin independently predicts mortality in patients with stable coronary artery disease: the CLARICOR trial. Scand J Clin Lab Invest. 2014;74:657–64.
Inker LA, Schmid CH, Tighiouart H, Eckfeldt JH, Feldman HI, Greene T, et al. Estimating glomerular filtration rate from serum creatinine and cystatin C. N Engl J Med. 2012;367:20–9.
Helweg-Larsen K. The Danish register of causes of death. Scand J Public Health. 2011;39(7 30.Suppl):26–9.
Lynge E, Sandegaard JL, Rebolj M. The Danish National Patient Register. Scand J Public Health. 2011;39(7 suppl):30–3.
Kjoller E, Hilden J, Winkel P, Frandsen NJ, Galatius S, Jensen G, et al. Good interobserver agreement was attainable on outcome adjudication in patients with stable coronary heart disease. J Clin Epidemiol. 2012;65:444–53.
Kjoller E, Hilden J, Winkel P, Galatius S, Frandsen NJ, Jensen GB, et al. Agreement between public register and adjudication committee outcome in a cardiovascular randomized clinical trial. Am Heart J. 2014;168:197–204.
Wong ND. Epidemiological studies of CHD and the evolution of preventive cardiology. Nat Rev Cardiol. 2014;11:276–89.
Solomon SD, Claggett B, Desai AS, Packer M, Zile M, Swedberg K, et al. Influence of ejection fraction on outcomes and efficacy of sacubitril/valsartan (LCZ696) in heart failure with reduced ejection fraction. Circ Heart Fail. 2016;9:e00274.
Boot RG, van Achterberg TAE, van Aken BE, Renkema GH, Jacobs MHM, Aerts JMFG, de Vries CJM. Strong induction of members of the chitinase family of proteins in atherosclerosis. Chitotriosidase and human cartilage gp-39 expressed in lesion macrophages. Arterioscler Thromb Vasc Biol. 1999;19:687–94.
Vik A, Mathiesen EB, Noto AT, Sveinbjornsson B, Brox J, Hansen JB. Serum osteoprotegerin is inversely assiociated to carotid plaque echogenecity in humans. Atherosclerosis. 2007;191:128–34.
Lindstedt L, Lee M, Oörni K, Brömme D, Kovanen PT. Cathepsins F and S block HDL3-induced cholesterol efflux from macrophage foam cells. Biochem Biophys Res Commun. 2003;312(4):1019–24.
Rodgers KJ, Watkins DJ, Miller AL, Chan PY, Karanam S, Brissette WH, Long CJ, Jackson CL. Destabilizing role of cathepsin S in murine atherosclerotic plaques. Arterioscler Thromb Vasc Biol. 2006;26(4):851–6.
Sluimer JC, Gasc JM, van Wanroij JL, Kisters N, Groeneweg M, Sollewijn Gelpke MD, Cleutjens JP, van den Akker LH, Corvol P, Wouters BG, Daemen MJ, Bijnens AP. Hypoxia, hypoxia-inducible transcription factor, and macrophages in human atherosclerotic plaques are correlated with intraplaque angiogenesis. J Am Coll Cardiol. 2008;51(13):1258–65.
Heljasvaara R, Nyberg P, Luostarinen J, Parikka M, Heikkilä P, Rehn M, Sorsa T, Salo T, Pihlajaniemi T. Generation of biologically active endostatin fragments from human collagen XVIII by distinct matrix metalloproteases. Exp Cell Res. 2005;307(2):292–304.
Libby P. Inflammation in atherosclerosis. Arterioscler Thromb Vasc Biol. 2012;32:2045-51.
Carlsson AC, Juhlin CC, Larsson TE, Larsson A, Ingelsson E, Sundström J, Lind L, Ärnlöv J. Soluble tumor necrosis factor receptor 1 (sTNFR1) is associated with increased total mortality due to cancer and cardiovascular causes—findings from two community based cohorts of elderly. Atherosclerosis. 2014;237(1):236–42.
Liu J, Ma L, Yang J, Ren A, Sun Z, Yan G, Sun J, Fu H, Xu W, Hu C, Shi GP. Increased serum cathepsin S in patients with atherosclerosis and diabetes. Atherosclerosis. 2006;186(2):411–9.
Seko Y, Fukuda S, Nagai R. Serum levels of endostatin, vascular endothelial growth factor (VEGF) and hepatocyte growth factor (HGF) in patients with acute myocardial infarction undergoing early reperfusion therapy. Clin Sci (Lond). 2004;106(5):439–42.
Blankenberg S, McQueen MJ, Smieja M, Pogue J, Balion C, Lonn E, Rupprecht HJ, Bickel C, Tiret L, Cambien F, Gerstein H, Münzel T, Yusuf S, HOPE Study Investigators. Comparative impact of multiple biomarkers and N-Terminal pro-brain natriuretic peptide in the context of conventional risk factors for the prediction of recurrent cardiovascular events in the Heart Outcomes Prevention Evaluation (HOPE) Study. Circulation. 2006;114:201–8.
Carbone F, Mach F, Montecucco F. Update on the role of neutrophils in atherosclerotic plaque vulnerability. Curr Drug Targets. 2015;16:321–33.
Elneihoum AM, Falke P, Hedblad B, Lindgärde F, Ohlsson K. Leukocyte activation in atherosclerosis: correlation with risk factors. Atherosclerosis. 1997;131:79–84.
Lin DY, Wei LJ, Ying Z. Checking the Cox model with cumulative sums of martingale-based residuals. Biometrika. 1993;80:557–72.
Kerr KF, Wang Z, Janes H, McClelland RL, Psaty BM, Pepe MS. Net reclassification indices for evaluating risk-prediction instruments: a critical review. Epidemiology. 2014;25:114–21.
Newson RB. Comparing the predictive powers of survival models using Harrell’s C or Somers’ D. Stata J. 2010;10:339–58.
Little R. A test of missing completely at random for multivariate data with missing values. J Am Stat Assoc. 1988;83:1198–1202.
Pepe MS, Fan J, Feng Z, Gerds T, Hilden J. The Net Reclassification Index (NRI): a misleading measure of prediction improvement even with independent test data sets. Stat Biosci. 2015;7:282–95. Epub 2014 Aug 23.
Vickers AJ, Elkin EB. Decision curve analysis: a novel method for evaluating prediction models. Med Decis Making. 2006;26:565–74. doi:10.1177/0272989X06295361.
We thank the CLARICOR trial participants. We thank the investigators and other staff involved in the first phases of the CLARICOR trial (for full list of names please see references 10 and 13). We thank our funders of the CLARICOR trial. The Copenhagen Trial Unit, Centre for Clinical Intervention Research, is thanked for providing monetary support for part of biochemical analyses for the PREMAC study as well as wages for Per Winkel, Janus C. Jacobsen, and Christian Gluud.
This study was funded by the Copenhagen Trial Unit, Centre for Clinical Intervention Research; original funders of the CLARICOR trial and The Swedish Research Council, Swedish Heart-Lung foundation; Thuréus foundation; Marianne and Marcus Wallenberg Foundation, Dalarna University; and Uppsala University.
Availability of data and materials
All pertinent anonymised data will be uploaded at ZENODO (http://zenodo.org/) when the individual manuscripts have been published.
JCJ, CG, and PW contributed substantially to the concept and design and drafted the manuscript. GBJ, EK, JK, HJK, AS, PW, AL, and JÄ reviewed and summarized prognostic and patho-physiological literature relevant for the candidate prognostic markers. JH and TL critically revised the statistical analysis section. All authors revised the manuscript critically for important intellectual content, gave final approval of version to be published, and agreed to be accountable for all aspects of the work in assuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.
The authors declare that they have no competing interests.
Ethics approval and consent to participate was given by
VEKKF01-076/99; Danish Medicines Agency 2612-975; Danish Data Protection Agency 1999-1200-174; VEK H-B-2009-015.
About this article
Cite this article
Winkel, P., Jakobsen, J.C., Hilden, J. et al. Predictors for major cardiovascular outcomes in stable ischaemic heart disease (PREMAC): statistical analysis plan for data originating from the CLARICOR (clarithromycin for patients with stable coronary heart disease) trial. Diagn Progn Res 1, 10 (2017). https://doi.org/10.1186/s41512-017-0009-y