Eht 439
Eht 439
Received 26 July 2013; revised 17 September 2013; accepted 26 September 2013; online publish-ahead-of-print 14 October 2013
Aims Hypertrophic cardiomyopathy (HCM) is a leading cause of sudden cardiac death (SCD) in young adults. Current risk algo-
rithms provide only a crude estimate of risk and fail to account for the different effect size of individual risk factors. The aim
of this study was to develop and validate a new SCD risk prediction model that provides individualized risk estimates.
.....................................................................................................................................................................................
Methods The prognostic model was derived from a retrospective, multi-centre longitudinal cohort study. The model was devel-
and results oped from the entire data set using the Cox proportional hazards model and internally validated using bootstrapping. The
cohort consisted of 3675 consecutive patients from six centres. During a follow-up period of 24 313 patient-years
(median 5.7 years), 198 patients (5%) died suddenly or had an appropriate implantable cardioverter defibrillator
(ICD) shock. Of eight pre-specified predictors, age, maximal left ventricular wall thickness, left atrial diameter, left ven-
tricular outflow tract gradient, family history of SCD, non-sustained ventricular tachycardia, and unexplained syncope
were associated with SCD/appropriate ICD shock at the 15% significance level. These predictors were included in the
final model to estimate individual probabilities of SCD at 5 years. The calibration slope was 0.91 (95% CI: 0.74, 1.08),
C-index was 0.70 (95% CI: 0.68, 0.72), and D-statistic was 1.07 (95% CI: 0.81, 1.32). For every 16 ICDs implanted in
patients with ≥4% 5-year SCD risk, potentially 1 patient will be saved from SCD at 5 years. A second model with the
data set split into independent development and validation cohorts had very similar estimates of coefficients and perform-
ance when externally validated.
.....................................................................................................................................................................................
Conclusion This is the first validated SCD risk prediction model for patients with HCM and provides accurate individualized estimates
for the probability of SCD using readily collected clinical parameters.
-----------------------------------------------------------------------------------------------------------------------------------------------------------
Keywords Hypertrophic cardiomyopathy † Sudden cardiac death † Implantable cardioverter defibrillator † Risk prediction
model
* Corresponding author. Tel: +44 203 456 4801, Fax: +44 207 456 4901, Email: [email protected]
Published on behalf of the European Society of Cardiology. All rights reserved. & The Author 2013. For permissions please email: [email protected]
Risk prediction model for SCD in HCM 2011
parameters that reflect the severity of the underlying myocardial physicians). Deaths were assessed without knowledge of the candidate
disease. The presence or absence of these risk factors is then used predictor variables. The primary endpoint was SCD or an equivalent
to guide clinical decision-making with respect to prophylactic ICD event. Sudden cardiac death was defined as witnessed sudden death
implantation.1,2 Although observational cohort studies show that with or without documented ventricular fibrillation or death within 1 h
of new symptoms or nocturnal deaths with no antecedent history of wor-
this approach identifies patients with the greatest risk of SCD, valid-
sening symptoms.15 Aborted SCD during follow-up and appropriate ICD
ation of current algorithms suggests that they overestimate risk,
shock therapy were considered equivalent to SCD.16 – 21 Implantable car-
resulting in inappropriate prophylactic ICD implantation in a substan-
dioverter defibrillator shocks were considered appropriate if the treated
tial number of patients.3 The aim of this study was to derive and val- tachyarrhythmia was ventricular in origin in an identical manner to previ-
idate a new sudden death risk model that can be used to generate ous studies.16 – 21 All ICDs had the capacity to store intracardiac electro-
individualized risk estimates for SCD and improve the targeting of grams.
ICD therapy in patients with HCM.
Selection of predictor variables and coding
techniques based on chained equations.36 All predictors of missingness SCD.43,44 A value of 0.5 for C-index indicates no discrimination and a
were included in the multiple imputation model, together with the value equal to 1 indicates perfect discrimination. The D-statistic quanti-
outcome, all pre-specified predictors of the risk model, and the estimate fies the observed separation between subjects with low and high pre-
of the cumulative hazard function.37 A total of 25 imputed data sets were dicted risks as predicted by the model and can be interpreted as the
generated and the estimates were combined using Rubin’s rules.38 log hazard ratio for having SCD between the low and high risk groups
Patients with .50% missing predictors were excluded from model devel- of patients. A model with no discriminatory ability will produce a value
opment. of 0 for D-statistic, with increasing values indicating greater separation.
Sensitivity analyses seven of the eight predictors: maximal wall thickness in 11 (0.3%)
A model to adjust for centre effect was developed and validated using patients, fractional shortening in 286 (7.8%) patients, LA diameter
bootstrapping. An additional sensitivity analysis was carried out by devel- in 92 (2.5%) patients, left ventricular outflow gradient in 94 (2.5%)
oping a model using only patients with complete data, with and without patients, NSVT in 453 (12.3%) patients, family history of SCD in 7
centre. The same modelling strategy and validation measures as the ones (0.2%) patients, and unexplained syncope in 65 (1.8%) patients. Miss-
employed to develop and validate the SCD risk prediction model on the ingness was associated with New York Heart Association functional
entire data set were used. class III/IV, amiodarone treatment at baseline, ICD implantation
during follow-up, cardiovascular deaths not secondary to SCD/
Comparison with conventional risk factors heart failure/stroke, non-cardiovascular deaths, and the year of exit
Clinical practice guidelines consider severe hypertrophy (maximal wall from the study.
thickness ≥30 mm), non-sustained ventricular tachycardia (NSVT),
family history of SCD, and unexplained syncope as conventional risk
factors of SCD.1,2 To compare the risk model developed in this article Model development
All The Heart A Coruña 1St Department of Institute of University Hospital Monaldi
Hospital, UK University Cardiology, University Cardiology, Virgen de la Hospital,
Hospital, Spain of Athens, Greece Bologna, Italy Arrixaca, Spain Italy
.............................................................................................................................................................................................................................................
Variables are expressed as mean + standard deviation (SD), median and interquartile range (IQR) or counts and percentages as appropriate.
NYHA, New York Heart Association; ICD, implantable cardioverter defibrillator; AF, atrial fibrillation; NSVT, non-sustained ventricular tachycardia; LA, left atrium; LVOTG, left ventricular outflow tract gradient at rest or Valsalva; LVedd, left
ventricular end diastolic dimension; MWT, maximal wall thickness; FS, fractional shortening; FHSCD, family history of sudden cardiac death; SCD, sudden cardiac death.
a
Three patients were excluded from model development because .50% of predictors were missing. One of the three excluded patients had suffered SCD/appropriate ICD shock.
b
Range: 1%; 99% centiles: age: 16.9; 81.4 years, LA diameter: 28; 67 mm, LVOTGmax: 2; 154 mmHg, MWT: 10; 36 mm, FS: 15; 64%.
c
207 patients (5%) with FS , 27%, 489 patients (13%) with FS.50%.
C. O’Mahony et al.
Risk prediction model for SCD in HCM 2015
Table 4 Summary of the characteristics of patients with sudden cardiac death endpoints and univariable Cox regression
models
Predictor variable SCD group characteristics Hazard ratio 95% confidence P-value
(n 5 198)a,b interval
...............................................................................................................................................................................
Age (years) 42.5 + 15 0.988 0.979, 0.997 0.007
Maximal wall thickness (mm) 21.5 + 6 1.048 1.025, 1.071 ,0.001
Fractional shortening (%) 41.0 + 10 0.992 0.977, 1.008 0.344
Left atrial diameter (mm) 46.2 + 9 1.035 1.018, 1.052 ,0.001
Left ventricular outflow gradient (mmHg) 18 (6 –58) 1.005 1.001, 1.008 0.005
Family history of sudden cardiac death 73 (37%) 1.760 1.318, 2.350 ,0.001
Non-sustained ventricular tachycardia 62 (31%) 2.533 1.849, 3.469 ,0.001
Unexplained syncope 52 (26%) 2.326 1.693, 3.195 ,0.001
Variables are expressed as mean + standard deviation (SD), median and interquartile range (IQR) or counts and percentages as appropriate.
a
Range of values (minimum; maximum) in SCD group: age: 16.3; 77.4 years, maximal wall thickness: 9; 37 mm, fractional shortening: 15; 62%, left atrial diameter: 28; 70 mm, maximal left
ventricular outflow tract gradient: 2;190 mmHg.
b
Missing data in SCD group: maximal wall thickness: 3%, fractional shortening: 11%, left atrial diameter: 6%, left ventricular outflow tract gradient: 3%, non-sustained ventricular
tachycardia: 19%, unexplained syncope: 2%.
2016 C. O’Mahony et al.
Table 5 Sudden cardiac death risk prediction model and sensitivity analysis for centre effect
Predictor variable SCD risk prediction model Sensitivity analysis: model with centre
................................................................. ..................................................................
Hazard ratio (95% confidence interval) P-value Hazard ratio (95% confidence interval) P-value
...............................................................................................................................................................................
Age (years) 0.98 (0.97, 0.99) 0.001 0.98 (0.97, 0.99) ,0.001
Maximal wall thickness (mm) 1.17 (1.01, 1.37) 0.042 1.15 (0.99, 1.35) 0.069
Maximal wall thickness2 (mm2) 0.997 (0.99, 1.0003) 0.078 0.997 (0.99, 1.001) 0.116
Left atrial diameter (mm) 1.03 (1.01, 1.05) 0.006 1.02 (1.01, 1.04) 0.008
LV outflow gradient (mmHg) 1.004 (1.001, 1.01) 0.021 1.004 (1.0004, 1.01) 0.031
Family history SCD 1.58 (1.18, 2.13) 0.002 1.62 (1.20, 2.19) 0.002
NSVT 2.29 (1.64, 3.18) ,0.001 2.14 (1.53, 2.99) ,0.001
A 5-year SCD risk of ≥4% identified 60 (71%) of 84 SCD end- The usefulness of this model lies in providing accurate prognostic
points with 896 (30%) ICD implants in patients without SCD at 5 information that aids clinical decision making. The model achieved
years (Figure 2D). For every 16 ICD implantations in patients with this by showing good agreement between the predicted and
≥4% 5-year SCD risk, 1 patient can potentially be saved from SCD observed hazards of SCD and by demonstrating the ability to separ-
at 5 years. Patients not reaching the SCD endpoint at 5 years ate patients with regard to their 5-year risk of SCD.42 The C-indices
(n ¼ 2982) had a mean predicted 5-year SCD risk of 3.7% (IQR: indicated that the proposed risk prediction model has superior dis-
1.8 –4.5%), while the corresponding figures for those reaching the crimination compared with the model of conventional risk factors
SCD endpoint (n ¼ 84) were 7.3% (IQR: 3.4 –10%). used in contemporary clinical practice. The sensitivity analysis
demonstrated that the relationship between the predictors and
SCD remains unchanged with the inclusion of centre in the model
and the risk model without centre is proposed for general clinical use.
Discussion The risk prediction model has the potential to improve the man-
This is the first validated risk prediction model for SCD in patients agement of patients with a solitary and multiple risk factors by simul-
with HCM. The model was derived from a large, diverse, and well- taneously reducing unnecessary and potentially harmful ICD
characterized population of patients followed at six different Euro- implants in patients who do not suffer SCD and correctly identifying
pean centres and provides accurate, individualized estimates for the majority of those who suffer SCD and are most likely to benefit
the probability of SCD using readily collected clinical parameters. from an ICD. Currently, patients without conventional risk factors
The broad patient inclusion criteria of the study mean that the are reassured and reassessed and are not routinely offered ICD
model can be used in the majority of adult patients with HCM, includ- therapy.1,2 However, approximately one-third of all SCD come
ing those with mild disease identified during family screening. from this subgroup of patients, and contemporary management strat-
Current clinical guidelines for HCM in the USA and Europe recom- egies fail to address this problem. The risk prediction model may help
mend SCD risk algorithms based on a simple summation of a limited identify a small proportion of SCD in this group which represents an
number of binary clinical parameters (NSVT, severe hypertrophy, un- improvement when compared with current clinical practice, but the
explained syncope, family history of SCD, and abnormal BPRE).1,2 performance of the model in this patient subgroup is not optimal.
Even though this approach has been used in clinical practice for The probability of SCD at 5 years is derived from a range of readily
more than a decade, it provides only a very crude estimate of relative available clinical parameters, each with a unique contribution to risk.
risk of SCD and fails to account for the different effect size of individ- For example, consider the management of two patients with NSVT,
ual risk factors.3 Moreover, some risk factors such as hypertrophy are maximal wall thickness of 23 mm, and LA diameter of 44 mm. One is
considered as binary variables when in fact they are associated with a aged 24 years with a resting LVOTGmax of 64 mmHg and the other is
continuous increase in SCD risk.25 As a result, existing algorithms 64 years old with an LVOTGmax of 36 mmHg. Current guidelines treat
have a low positive predictive accuracy for SCD that results in the un- these two patients identically as they each have a single risk factor
necessary treatment of patients who are at intrinsically low risk.3 (NSVT). By applying the clinical risk prediction model in this clinical
Risk prediction model for SCD in HCM 2017
Table 6 The simulated effect of using different thresholds of 5-year sudden cardiac death risk to implant an implantable
cardioverter defibrillator in patient groups defined using conventional risk factors as recommended in current treatment
guidelines
Patient groups defined using four Estimate for 5-year SCD risk Patients reaching SCD Patients without SCD
conventional risk factors used to implant ICD (%) endpointa endpointb
................................. ...................................
ICD No ICD No ICD ICD
implant Implant implant implant
...............................................................................................................................................................................
Two or more conventional risk factors: 412 patients ≥3 32 (100%) 0 9 (2%) 371 (98%)
with 32 SCD endpoints at 5 years ≥4 32 (100%) 0 32 (8%) 348 (92%)
≥5 32 (100%) 0 69 (18%) 311 (82%)
≥6 31 (97%) 1 (3%) 131 (34%) 249 (66%)
a
Denominator for %¼ number of SCD endpoints.
b
Denominator for %¼ number of patients without SCD endpoints.
c
During the study period of 1972 –2011, of the 3066 patients with the necessary data required to calculate the 5-year SCD risk, 49% with ≥2 risk factors, 19% with a single risk factor,
and 5% with no conventional risk factors had an ICD.
rate (Supplementary material online, Table S2) which is reflected in Hypertrophic Cardiomyopathy
the nonlinear relation between maximal wall thickness and 5-year
SCD risk. This subgroup was also small in size, and pending further Outcomes Investigators (www.
analysis, the model should also be used cautiously in patients with HCMRisk.org)
maximal wall thickness≥ 35 mm.
This work is not an epidemiological study exploring new associations The following were also investigators in the study: Martin Ortiz-
with SCD. Instead, previously described associations were used to Genga MD (A Coruña University Hospital, Spain), Xusto Fernandez
develop a prognostic model. Predictors were only included if shown MD (A Coruña University Hospital, Spain), Vassiliki Vlagouli MD
to be independently associated with SCD in published multivariable (University of Athens, Greece), Chris Stefanadis MD (University of
survival analyses. This strategy achieved a high event-to-predictor Athens, Greece), Fabio Coccolo MD (University of Bologna, Italy),
ratio which improves the generalizability and accuracy of predictions. Maria Jose Oliva Sandoval MD (Murcia, Spain), Giuseppe Pacileo
Although prospective cohort studies are desirable,46 the low rate MD (Monaldi Hospital, Italy), Daniele Masarone MD (Monaldi
of SCD in HCM made a prospective design impractical. All participat- Hospital, Italy), Antonios Pantazis MD (The Heart Hospital, UK),
ing centres assess HCM patients in an almost identical manner, which Maite Tome-Esteban MD (The Heart Hospital, UK), Shaughan Dickie
allowed the use of retrospectively collected data from all sites. The (The Heart Hospital, UK), Pier D Lambiase PhD (The Heart Hospital,
prediction model will be further improved through the incorporation UK), and Shafiq Rahman PhD (University College London, UK).
of additional cohorts, thus increasing the number of events that will
reduce model optimism.42 Authors’ contributions
C.O’M designed the study, collected and interpreted the data, carried
out the descriptive statistical analysis, and wrote the article. P.M.E.
Conclusions designed the study, interpreted the data, and wrote the article.
The risk prediction model proposed in this study provides accurate R.Z.O. was involved in study design, led the statistical aspects of
prognostic information regarding SCD in HCM. The model provides the risk modelling and wrote the article. F.J. carried out the statistical
the basis for a rational and informed approach to treatment that analysis. Dr Rahman evaluated some of the measures used for model
empowers patients by helping them to understand the relative validation as part of his PhD thesis with Prof Omar. M.P. carried out
merits of prophylactic therapies including implantable cardioverter the statistical analysis in relation to the external validation of the
defibrillators. model. LM., A.A., E.B., J.R.G., G.L., and C.R. collected and interpreted
Risk prediction model for SCD in HCM
Downloaded from https://academic.oup.com/eurheartj/article/35/30/2010/467191 by guest on 01 October 2024
Figure 2 The impact of the risk prediction model on clinical decision making. Conventional risk factors (wall thickness ≥30 mm, NSVT, family history of SCD and unexplained syncope) were used to
classify patients in those with multiple (A), single (B), no risk factors (C). Each figure shows the implications of implanting an ICD in all (bottom bar) or none of the patients (top bar), while the rest of the bars
show the impact of using a threshold of ≥3, ≥4, ≥5, ≥6 5-year SCD risk to implant an ICD. The impact on the whole cohort is shown in (D). The raw data are shown in Table 6. Only patients with all data
necessary for the 5-year SCD risk calculation were included (n ¼ 3066).
2019
2020 C. O’Mahony et al.
the data and critically reviewed the manuscript. Ortiz-Genga, Fernan- 4. Barriales-Villa R, Centurion-Inda R, Fernandez-Fernandez X, Ortiz MF,
Perez-Alvarez L, Rodriguez GI, Hermida-Prieto M, Monserrat L. Severe cardiac con-
dez, Vlagouli, Sandoval, Pacileo, Pantazis, Tome-Esteban, Dickie,
duction disturbances and pacemaker implantation in patients with hypertrophic car-
Lambiase, Stefanadis, Coccolo, and Masarone were involved in data diomyopathy. Rev Esp Cardiol 2010;63:985 – 988.
collection and interpretation. W.J.M. was involved in the drafting of 5. Di Salvo G, Pacileo G, Limongelli G, Baldini L, Rea A, Verrengia M, D’Andrea A,
the article and revising it critically for important intellectual content. Russo MG, Calabro R. Non sustained ventricular tachycardia in hypertrophic cardio-
myopathy and new ultrasonic derived parameters. J Am Soc Echocardiogr 2010;23:
581 –590.
6. Oliva-Sandoval MJ, Ruiz-Espejo F, Monserrat L, Hermida-Prieto M, Sabater M,
Ethics approval Garcia-Molina E, Ortiz M, Rodriguez-Garcia MI, Nunez L, Gimeno JR,
Castro-Beiras A, Valdes M. Insights into genotype-phenotype correlation in hyper-
Patients at A Coruña University Hospital (Spain), Unit of Inherited trophic cardiomyopathy. findings from 18 Spanish families with a single mutation in
Cardiovascular diseases at the 1st Department of Cardiology MYBPC3. Heart 2010;96:1980 –1984.
(Athens), University Hospital Virgen de la Arrixaca (Spain), and 7. Monserrat L, Lopez B, Gonzalez A, Hermida M, Fernandez X, Ortiz M,
Barriales-Villa R, Castro-Beiras A, Diez J. Cardiotrophin-1 plasma levels are asso-
Monaldi Hospital (Italy) provided written informed consent. The ciated with the severity of hypertrophy in hypertrophic cardiomyopathy. Eur
24. Spirito P, Bellone P, Harris KM, Bernabo P, Bruzzi P, Maron BJ. Magnitude of left ven- 34. Olivotto I, Cecchi F, Casey SA, Dolara A, Traverse JH, Maron BJ. Impact of atrial fib-
tricular hypertrophy and risk of sudden death in hypertrophic cardiomyopathy. N rillation on the clinical course of hypertrophic cardiomyopathy. Circulation 2001;104:
Engl J Med 2000;342:1778 –1785. 2517 –2524.
25. Elliott PM, Gimeno B Jr, Mahon NG, Poloniecki JD, McKenna WJ. Relation between 35. Peduzzi P, Concato J, Feinstein AR, Holford TR. Importance of events per independ-
severity of left-ventricular hypertrophy and prognosis in patients with hypertrophic ent variable in proportional hazards regression analysis. II. Accuracy and precision of
cardiomyopathy. Lancet 2001;357:420–424. regression estimates. J Clin Epidemiol 1995;48:1503 –1510.
26. Kofflard MJM, Ten Cate FJ, van der Lee C, van Domburg RT. Hypertrophic cardio- 36. van Buuren S, Boshuizen HC, Knook DL. Multiple imputation of missing blood pres-
myopathy in a large community-based population: clinical outcome and identifica- sure covariates in survival analysis. Stat Med 1999;18:681 – 694.
tion of risk factors for sudden cardiac death and clinical deterioration. J Am Coll 37. White IR, Royston P, Wood AM. Multiple imputation using chained equations: issues
Cardiol 2003;41:987 –993. and guidance for practice. Stat Med 2011;30:377 –399.
27. Autore C, Bernabo P, Barilla CS, Bruzzi P, Spirito P. The prognostic importance of left 38. Rubin D. Multiple Imputation for Nonresponse in Surveys. New York: John Wiley and
ventricular outflow obstruction in hypertrophic cardiomyopathy varies in relation to Sons; 1987.
the severity of symptoms. J Am Coll Cardiol 2005;45:1076 –1080. 39. Ambler G, Seaman S, Omar RZ. An evaluation of penalised survival methods for
developing prognostic models with rare events. Stat Med 2012;31:1150 –1161.
28. Olivotto I, Maron MS, Adabag AS, Casey SA, Vargiu D, Link MS, Udelson JE, Cecchi F,
40. Schoenfeld DAVI. Partial residuals for the proportional hazards regression model.
Maron BJ. Gender-related differences in the clinical presentation and outcome of
Biometrika 1982;69:239 – 241.
hypertrophic cardiomyopathy. J Am Coll Cardiol 2005;46:480 –487.