Abstract
Background and Aims
Pathogenic desmoplakin (DSP) gene variants are associated with the development of a distinct form of arrhythmogenic cardiomyopathy known as DSP cardiomyopathy. Patients harbouring these variants are at high risk for sustained ventricular arrhythmia (VA), but existing tools for individualized arrhythmic risk assessment have proven unreliable in this population.
Methods
Patients from the multi-national DSP-ERADOS (Desmoplakin SPecific Effort for a RAre Disease Outcome Study) Network patient registry who had pathogenic or likely pathogenic DSP variants and no sustained VA prior to enrolment were followed longitudinally for the development of first sustained VA event. Clinically guided, step-wise Cox regression analysis was used to develop a novel clinical tool predicting the development of incident VA. Model performance was assessed by c-statistic in both the model development cohort (n = 385) and in an external validation cohort (n = 86).
Results
In total, 471 DSP patients [mean age 37.8 years, 65.6% women, 38.6% probands, 26% with left ventricular ejection fraction (LVEF) < 50%] were followed for a median of 4.0 (interquartile range: 1.6–7.3) years; 71 experienced first sustained VA events {2.6% [95% confidence interval (CI): 2.0, 3.5] events/year}. Within the development cohort, five readily available clinical parameters were identified as independent predictors of VA and included in a novel DSP risk score: female sex [hazard ratio (HR) 1.9 (95% CI: 1.1–3.4)], history of non-sustained ventricular tachycardia [HR 1.7 (95% CI: 1.1–2.8)], natural logarithm of 24-h premature ventricular contraction burden [HR 1.3 (95% CI: 1.1–1.4)], LVEF < 50% [HR 1.5 (95% CI: .95–2.5)], and presence of moderate to severe right ventricular systolic dysfunction [HR 6.0 (95% CI: 2.9–12.5)]. The model demonstrated good risk discrimination within both the development [c-statistic .782 (95% CI: .77–.80)] and external validation [c-statistic .791 (95% CI: .75–.83)] cohorts. The negative predictive value for DSP patients in the external validation cohort deemed to be at low risk for VA (<5% at 5 years; n = 26) was 100%.
Conclusions
The DSP risk score is a novel model that leverages readily available clinical parameters to provide individualized VA risk assessment for DSP patients. This tool may help guide decision-making for primary prevention implantable cardioverter-defibrillator placement in this high-risk population and supports a gene-first risk stratification approach.
Keywords: DSP, VA, DSP cardiomyopathy, ICD, Sudden cardiac death, Risk prediction, ACM, ARVC
Structured Graphical Abstract
See the editorial comment for this article ‘Sudden cardiac death risk stratification in heritable cardiomyopathies: phenotype-based to genotype-based risk scores’, by R. Tadros, https://doi.org/10.1093/eurheartj/ehae381.
Introduction
Desmoplakin (DSP) is a cardiac desmosomal protein that plays a critical role in myocardial force transmission,1 and pathogenic variants in its encoding gene (DSP) are an important cause of arrhythmogenic cardiomyopathies.2,3 While traditionally associated with the development of arrhythmogenic right ventricular (RV) cardiomyopathy (ARVC), DSP-related disease is characterized by several unique features and has become increasingly recognized as a distinct clinical entity.4–6 Desmoplakin variant carriers demonstrate high rates of left ventricular (LV) pathology including both LV systolic dysfunction and LV myocardial fibrosis, suffer from episodes of inflammatory myocardial injury, and are at particularly increased risk for life-threatening ventricular arrhythmias (VA).4–7 Implantable cardioverter-defibrillators (ICDs) are thus an essential tool for preventing DSP-related sudden cardiac death (SCD).8 However, while current guidelines agree upon the benefits of offering an ICD for secondary prevention of sustained VA,8–10 the indications for prophylactic ICD use in patients with P/LP DSP variants are less clear.
The 2019 ARVC risk calculator is a validated, patient-specific risk prediction tool designed to assist clinicians in decision-making regarding the use of primary prevention ICDs in ARVC.11 However, while this tool reliably predicts sustained VA events in patients with classical ARVC,12–14 it has demonstrated poor performance in both patients with left-sided disease15,16 and those with P/LP DSP variants.17,18 The inability to provide dependable risk assessment therefore represents a major deficiency in the clinical care of these patients. While genotype-tailored strategies have demonstrated benefit in other forms of arrhythmic cardiomyopathy,19–21 the rarity of patients with DSP variants described in the current literature has thus far limited the development of DSP-specific risk models.
With this goal in mind, we identified DSP patients without sustained VA prior to enrolment from within the DSP-ERADOS (Desmoplakin SPecific Effort for a RAre Disease Outcome Study) cohort of patients with P/LP DSP variants. Using readily available clinical variables, we developed and validated the DSP risk score, a clinical tool providing personalized VA risk assessment for DSP patients being considered for primary prevention ICD implantation.
Methods
This study conformed to the Declaration of Helsinki and was approved by local ethics and/or institutional review boards. To maintain patient confidentiality, data and study materials have not been made publicly available, but a limited data set may be made available upon reasonable request.
Study population
This was a retrospective analysis of patients from the DSP registry maintained by the DSP-ERADOS Network group.5 This international collaboration includes 26 academic institutions across nine countries, with each institution functioning as an independent, prospective, observational patient registry. Individuals seen at one of these institutions before 1 March 2023 were entered in the registry if they harboured a P/LP variant in the DSP gene in accordance with the American College of Medical Genetics and Genomics (ACMGG) criteria.22 The study population thus includes both proband patients (defined as the first affected individual in a family seeking medical attention) and relatives in whom the P/LP DSP variant was identified for diagnostic purposes or cascade screening. Registry patients were included in the current study if they (i) underwent ambulatory cardiac monitoring at the time of their initial evaluation, (ii) underwent either echocardiographic or cardiac magnetic resonance (CMR) imaging at the time of their initial evaluation, and (iii) had clinical follow-up available for longitudinal outcome ascertainment. Patients were excluded from this study if they had a history of sustained VA prior to enrolment or if sustained VA was their presenting symptom. Given the gene-first nature of this cohort, we refer to participants of this study as DSP patients throughout the manuscript.
The overall cohort was divided into two groups along institutional lines for model development and external model validation. Briefly, contributing centres were selected as part of the external model validation cohort such that this cohort (i) comprised ∼20% of the overall patient cohort and (ii) reflected the diverse geographic locations of the DSP-ERADOS Network. Specific institutions comprising each of the two cohorts are shown in Figure 1.
Data collection
Data were collected independently at each institution according to a set of standardized definitions (see Supplementary data online, Table S1). Enrolment (time = 0) of DSP patients was defined as the first clinical evaluation of the patient at which initial work-up including echocardiography or CMR imaging and ambulatory cardiac monitor was either available or obtained within the next 90 days. Available demographics, patient medical history, genetic test results, and baseline cardiac instrumental exams, including 12-lead electrocardiogram (ECG), echocardiography, CMR, and ambulatory cardiac monitor, were retrieved for each patient, if available. All DSP genetic variants initially considered P/LP underwent centralized expert review by specialists in cardiac genetics (B.M. and C.A.J.), in accordance with the ACMGG guidelines and a previously published arrhythmogenic cardiomyopathy variant protocol.22,23 Moderate to severe RV systolic dysfunction was defined as the presence of RV ejection fraction (RVEF) < 40% on CMR or RV fractional area change < 33% on echocardiography (if CMR was not available); further RV functional categorizations are presented in Supplementary data online, Methods S1. While protocolized quantitation of late gadolinium enhancement (LGE) burden on CMR was not available in this multinational cohort, high-risk LGE pattern was defined as previously published24 according to the distribution of LGE involvement: (i) sub-epicardial, (ii) transmural, or (iii) involving both interventricular septum and LV free wall.
Outcomes
The primary outcome of this study was the development of a first sustained VA event, defined as a composite of sustained ventricular tachycardia (VT) (lasting ≥30 s at ≥100 b.p.m. or with haemodynamic compromise requiring cardioversion), ventricular fibrillation/flutter, sudden cardiac arrest (SCD), or appropriate ICD interventions including both ICD shock or anti-tachycardia pacing for sustained VA. We additionally defined life-threatening VA (LTVA) events as a composite of SCD, aborted SCD, ventricular fibrillation, or fast sustained VT (≥250 b.p.m.). Secondary outcomes used for competing risk analysis included heart transplantation and overall mortality. Outcomes were adjudicated locally at each centre via review of ECG, ICD interrogations, and patient medical records.
Statistics
Analyses were performed in Python version 3.9.13 using the open-source Pandas and Lifelines statistical code libraries. Baseline characteristics were presented as mean ± standard deviation for continuous variables and as proportions for categorical variables; continuous variables that did not follow normal distributions were presented as median with interquartile range (IQR). Unpaired Student’s t-test, proportional Z-test, and Mann–Whitney U test were used to test differences between patients who did and did not develop sustained VA for normally distributed continuous, categorical, and non-normally distributed continuous variables, respectively.
Missing data were considered to be missing at random and imputed using multiple imputation with chained equations; a more detailed description of imputation and missingness is presented in Supplementary data online, Methods - Data Inputation .
Model development
The DSP risk score was developed in accordance with the Transparent Reporting of a multivariable Prediction model for individual prognosis or Diagnosis (TRIPOD) statement.25 The association between potential arrhythmic risk factors in patients with DSP variants and the primary outcome of first sustained VA was assessed using Cox regression, and both univariable and multivariable hazard ratios (HRs) are reported. Proportional hazard assumptions were verified using scaled Schoenfeld residual testing. Conversion of variables from continuous to dichotomous values was decided based upon graphical assessment of the hazard for sustained VA over different clinically relevant ranges of those variables. Candidate variables known to have strong associations with VA based on previously published work4,5,7,8,11,12 were considered iteratively for addition to the prediction model (a complete list is included in Supplementary data online, Methods S2). The continued predictive benefit of each added variable was assessed using log-likelihood ratio testing, and model development was finalized once the test P-value exceeded .1 or the addition of five predictive variables (allowing for ∼10 sustained VA events per variable), whichever came first. Given the high missingness in CMR-derived data, non-CMR variables were preferentially selected for inclusion. While this study was not adequately powered to develop and externally validate an independent model for the prediction of LTVA events, we assessed the proportion of predicted arrhythmic risk that was comprised of LTVA events using Cox regression over the combined DSP cohort with the previously developed DSP risk score prognostic index as the sole predictor of LTVA.
Model performance was assessed using concordance-based c-statistics in both the development and external validation cohorts. Here, 95% confidence intervals (CIs) were assessed using five-fold cross-validation. Optimism correction was not used. Model calibration was assessed graphically as a comparison between predicted and observed rates of sustained VA within the overall cohort. We also assessed the clinical utility of the model for prospective risk stratification and ICD decision-making by dividing patients into the following risk groups: (i) low-risk, defined as predicted 5-year VA risk < 5%, (ii) intermediate-risk, defined as predicted 5-year VA risk between 5% and 20%, and (iii) high-risk, defined as predicted 5-year VA risk > 20%. The cumulative survival free from sustained VA was estimated using the Kaplan–Meier method. Event rates are reported as averages over the 5-year period following the initial patient evaluation. Follow-up duration was calculated from the date of initial evaluation to the date of the study outcome or censoring, which was defined as death from any other cause, heart transplantation, or the most recent follow-up visit at which the outcomes could be ascertained. Log-rank testing was used to assess differences in survival curves between subgroups. Competing risk analysis was also performed using an Aalen–Johansen risk estimator and incident VA, heart transplantation, and mortality as competing risks. We calculated sensitivity, specificity, positive (PPV), and negative (NPV) predictive values for sustained VA events at 5 years at these clinically relevant risk-group cut-offs based on previously published methods for estimating these metrics in survival data.26
Results
Patient characteristics and outcomes
The overall cohort comprised 471 DSP patients followed for a median of 4.0 (IQR: 1.6–7.3) years. The average age of these patients was 37.8 ± 17.2 years, 309 (65.6%) were women, and 182 (38.6%) were probands. There were 211 different P/LP DSP variants (see Supplementary data online, Table S2); variants were truncating in 408 patients (86.6%) and missense in 58 patients (12.3%). Echocardiography was available for 424 (90.0%) patients, and CMR evaluation was performed in 387 (82.2%) patients. Baseline patient characteristics stratified by the presence or absence of VA outcome are presented in Table 1. Baseline characteristics for patients included in the model development vs. external model validation cohorts, and by contributing centres of the DSP-ERADOS Network, are presented in Supplementary data online, Tables S3 and Table S4, respectively.
Table 1.
Variable | All patients (N = 471) |
VA free survival (n = 400) |
Incident VA (n = 71) |
P-value |
---|---|---|---|---|
Age (years) | 37.8 (±17.2) | 37.9 (±17.4) | 37.3 (±15.8) | .78 |
Female Sex | 309 (65.6%) | 254 (63.5%) | 55 (77.5%) | .02 |
Non-White ethnicity | 31 (6.6%) | 25 (6.2%) | 6 (8.5%) | .49 |
History of cardiac syncope | 46 (9.8%) | 36 (9.0%) | 10 (14.1%) | .18 |
Symptomatic dyspnoea | 98 (20.8%) | 78 (19.5%) | 20 (28.2%) | .10 |
Proband status | 182 (38.6%) | 140 (35.0%) | 42 (59.2%) | .0001 |
Missense variant | 58 (12.3%) | 49 (12.2%) | 9 (12.7%) | .92 |
Truncating variant | 408 (86.6%) | 346 (86.5%) | 62 (87.3%) | .85 |
Non-sense | 194 (41.2%) | 162 (40.5%) | 32 (45.1%) | .47 |
Frameshift | 183 (38.9%) | 158 (39.5%) | 25 (35.2%) | .49 |
Splice | 31 (6.6%) | 26 (6.5%) | 5 (7.0%) | .87 |
Deletion/duplication | 5 (1.1%) | 5 (1.2%) | 0 (.0%) | .34 |
Number of ECG leads with T-wave inversion | 1 [0–3] | 1 [0–2] | 3 [0–4] | <.0001 |
NSVT on ambulatory monitor | 81 (17.2%) | 53 (13.2%) | 28 (39.4%) | <.0001 |
24-h PVC count | 500 [12–2346] | 300 [5–1699] | 2000 [798–4519] | <.0001 |
PVC count/24 h > 500 | 235 (49.9%) | 174 (43.5%) | 61 (85.9%) | <.0001 |
LVEF | 53.8 (±11.9) | 55.0 (±11.1) | 47.2 (±13.7) | <.0001 |
LVEF < 50% | 121 (25.9%) | 83 (21.0%) | 38 (53.5%) | <.0001 |
Moderate to severe RV dysfunction | 12 (2.6%) | 3 (.8%) | 9 (12.7%) | <.0001 |
CMR evaluation | 387 (82.2%) | 331 (82.8%) | 56 (78.9%) | .43 |
Any LGE | 245 (63.3%) | 200 (60.4%) | 45 (80.4%) | .038 |
High-risk LGE pattern | 185 (47.8%) | 145 (43.8%) | 40 (71.4%) | .001 |
Right ventricular LGE | 32 (8.3%) | 23 (6.9%) | 9 (16.1%) | .03 |
Variables are expressed as frequency (%), mean ± standard deviation, or median [IQR], as appropriate.
CMR, cardiac magnetic resonance; ECG, electrocardiogram; LGE, late gadolinium enhancement; LVEF, left ventricular ejection fraction; NSVT, non-sustained ventricular tachycardia; PVC, premature ventricular contraction.
Sustained VA was experienced by 71 (15.1%) patients at a median of 3.9 (IQR: 1.6–6.7) years following enrolment. Initial sustained VA events included 10 aborted SCD events (14.1%), 13 ICD interventions for fast VA (VF or VT ≥ 250 b.p.m.) (18.3%), 28 ICD interventions for slower VT (<250 b.p.m.) (39.4%), and 20 spontaneous sustained VT events (<250 b.p.m. or self-terminating after 30 s) (28.2%). Life-threatening VA occurred in 26 DSP patients, including three patients in which a slower sustained VA event occurred prior to the LTVA (see Supplementary data online, Figure S1). Five (1.1%) patients died during follow-up. Fourteen (3%) patients underwent cardiac transplantation. Within the development and external model validation cohorts, 53 (13.8%) and 18 (20.9%) patients experienced sustained VA, respectively. Overall, 2.6% (95% CI: 2.0–3.5) of patients experienced a sustained VA event per year (Figure 2). Compared with patients who did not experience sustained VA events during follow-up, patients who did were more likely to be female (77.5% vs. 63.5%, P = .02), were more likely to be probands (59.2% vs. 38.6%, P = .0001), had greater numbers of T-wave inversions (TWIs) in inferior and anterior ECG leads (3 (IQR: 0–4) vs. 1 (IQR: 0–2), P < .0001], had greater burdens of both non-sustained ventricular tachycardia (NSVT) (13.2% vs. 39.4%, P < .0001) and premature ventricular contractions (PVCs) [2000 (IQR: 800–4500) vs. 300 (IQR: 5–1700), P < .0001], and had both lower LV ejection fraction (LVEF) (47 ± 14% vs. 55 ± 11%, P < .0001) and higher rates of moderate to severe RV dysfunction (12.7% vs. .8%, P < .0001). At the end of follow-up, ICDs were implanted in 184 (39.1%) patients, including 57 (80.3%) patients who did and 127 (31.8%) patients who did not experience a first sustained VA event.
To inform model development and confirm the relevance of the most important risk factors, relationships between clinically relevant cut-offs of continuous variables and HRs for sustained VA events were assessed graphically (Figure 3). Hazard for sustained VA events increased linearly with both the logarithm of 24-h PVC count and the number of inferior and anterior ECG leads with TWI. Hazard for sustained VA events increased sharply below an LVEF of 50% but remained relatively constant even for more severely reduced LVEF. Both moderate and severe RV dysfunction demonstrated increased hazards for sustained VA. Based on these results, we performed univariable and multivariable Cox regression between potential arrhythmic risk factors and incident VA (Table 2). Here, female sex, history of NSVT, presence of moderate or severe RV dysfunction, and the log of the burden of PVCs were strongly associated with incident VA on both univariable and multivariable Cox regression. Both LVEF < 50% and the number of inferior and anterior ECG leads with TWI were significantly associated with incident VA on univariable analysis. However, while there was a trend towards an association between LVEF < 50% and incident VA on multivariable analysis (P = .07), there was no association between number of inferior and anterior ECG leads with TWI and incident VA after adjustment for other risk factors. A competing risk assessment was performed as a sensitivity analysis using heart transplant and mortality as competing risks but did not impact results.
Table 2.
Variable | Univariable CPH (95% CI) | P-value | Multivariable CPH (95% CI) | P-value |
---|---|---|---|---|
Female sex | 1.905 (1.091–3.324) |
.023 | 1.914 (1.080–3.394) |
.03 |
NSVT on ambulatory monitor | 2.801 (1.738–4.514) |
<.0001 | 1.798 (1.073–3.011) |
.03 |
LVEF < 50% | 2.786 (1.742–4.456) |
<.0001 | 1.591 (.958–2.643) |
.07 |
Moderate or severe RV dysfunction | 6.503 (3.199–13.216) |
<.0001 | 6.543 (2.866–14.939) |
<.0001 |
ln(24-h PVC count) | 1.363 (1.211–1.534) |
<.0001 | 1.266 (1.117–1.436) |
<.0001 |
Number of ECG leads with T-wave inversion | 1.150 (1.055–1.254) |
.002 | .978 (.875–1.092) |
.69 |
Here, only those variables ultimately selected for the DSP risk score were included in multivariable modelling. Presence of LGE on CMR was excluded due to collinearity with LV systolic dysfunction (LVEF < 50%) and higher missingness in CMR data. Hazards are shown along with 95% CIs.
CI, confidence interval; CPH, Cox proportional hazard; ECG, electrocardiogram; LVEF, left ventricular ejection fraction; NSVT, non-sustained ventricular tachycardia; RV, right ventricular.
The desmoplakin risk score
Based upon a step-wise addition of variables to the model, the finalized DSP risk score (available at www.DSP-risk.com) included five clinical and imaging-based risk factors: female sex, history of NSVT, LVEF < 50%, moderate or severe RV dysfunction (RVEF < 40% on CMR or RV fractional area change < 33% on echocardiography), and the natural logarithm of the total 24-h PVC count. Analysis of other potential arrhythmic risk predictors is presented in Supplementary data online, Methods S2 and for LVEF as a continuous risk predictor in Supplementary data online, Methods S3. Here, the risk of sustained VA events at 5 years from the time of initial evaluation is calculated according to Equation (1).
(1) |
Here, PI represents the prognostic index, which is calculated according to Equation (2). Here, dichotomous predictors are defined numerically as 1 (if the predictor is present) or 0 (if the predictor is absent).
(2) |
Other risk factors considered for model addition included the number of TWI on inferior and anterior ECG leads, patient age, and the presence of high-risk LGE pattern on CMR, but these did not meet our pre-specified statistical criteria (P < .1 during log-likelihood ratio testing and >10 events per predictive variable).
Sub-analysis of the DSP risk score examining the proportion of arrhythmic risk comprised of LTVA yielded Equation (3).
(3) |
The DSP risk score demonstrated good discrimination, with c-statistics of .782 (95% CI: .765–.799) during model development and .791 (95% CI: .751–.830) during external validation. Prediction of LTVA over the combined DSP cohort using the DSP risk score prognostic index also had strong discrimination, with c-statistic of .752 (95% CI: .735–.768). Desmoplakin risk score calibration in the development and external validation cohorts (see Supplementary data online, Figure S2) and after stratification of proband/family member status (see Supplementary data online, Figure S3) was good. Model performance was good for both male [n = 162; c-statistic .780 (95% CI: .753–.807)] and female patients [n = 309; c-statistic .765 (95% CI: .760–.770)]; differences in patient characteristics by sex are presented in Supplementary data online, Table S5.
The DSP risk score was able to effectively stratify patient risk (Structured Graphical Abstract). Observed rates of sustained VA events were .4%/year (95% CI: .1–1.6) in those predicted to be at low-risk (n = 186, 39.5%), 2.8%/year (95% CI: 1.9–4.1) among those predicted to be intermediate-risk (n = 217, 46.1%), and 6.6%/year (95% CI: 4.5–9.4) among those predicted to be at high risk (n = 68, 14.4%). Log-rank testing showed significant differences between all predicted risk groups (P < .0001). Within the external testing cohort (n = 87), we assessed the clinical utility of using these risk groups to determine a patient’s hypothetical appropriateness for ICD implantation (Table 3). No DSP patients prospectively determined to be low risk (n = 26, 30.2%) had sustained VA events at 5 years, resulting in an NPV of 100%. Specific characteristics of patients in each of the three risk groups are presented in Supplementary data online, Table S6.
Table 3.
5% predicted 5-year predicted risk for VA | 20% predicted 5-year predicted risk for VA | |
---|---|---|
Number below risk threshold | 26 (30.2%) | 72 (83.7%) |
Sensitivity (95% CI) |
92% (77–100) |
43% (15–71) |
Specificity (95% CI) |
37% (24–50) |
93% (86–100) |
PPV (95% CI) |
24% (11–38) |
49% (18–80) |
NPV (95% CI) |
100% (100–100) |
90% (83–98) |
Positive likelihood ratio |
1.47 | 5.93 |
Negative likelihood ratio |
.21 | .61 |
CI, confidence interval; NPV, negative predictive value; PPV, positive predictive value; VA, ventricular arrhythmia.
We performed a sub-analysis examining the role of high-risk LGE pattern for the prediction of sustained VA events in the subset of patients with available CMR data (n = 387). High-risk LGE pattern was correlated with other risk predictors included in the DSP risk score (see Supplementary data online, Figure S4). While high-risk LGE pattern was strongly associated with sustained VA events on univariable Cox regression [HR 2.9 (95% CI: 1.6–5.5), P = .0009) within the development cohort (n = 321), this association was no longer significant after controlling for other risk factors [HR 1.54 (95% CI: .82–2.92), P = .18]; addition of high-risk LGE to the model could not be justified based on log-likelihood ratio testing (P = .86). In the low-risk group (5-year risk for VA < 5%; n = 149), only 1 of 51 patients (2%) with high-risk LGE pattern experienced a sustained VA event. This patient also had PVC/24 h > 1000 and their VA event occurred 6.1 years following their initial evaluation. Given the signal for additive prognostic information, we performed an additional Bayesian analysis of high-risk LGE pattern in the low- and intermediate-risk groups (5-year VA risk < 20%; n = 335) (see Supplementary data online, Methods S4 and Table S7); positive and negative likelihood ratios were 1.7 and .5, respectively, for patients with or without high-risk LGE.
Discussion
In this study, we present the first clinical tool for assessing individualized sustained VA risk in DSP patients without prior VA events and validate its performance in a distinct external cohort. This tool, which we term the DSP risk score, was able to reliably distinguish between those patients with incident-sustained VA during 5-year follow-up and those who survived free of sustained VA based upon five readily available clinical parameters: female sex, history of NSVT, LVEF < 50%, moderate or severe RV dysfunction, and the natural logarithm of the total 24-h PVC count. For those DSP patients deemed to have low baseline risk (<5% risk for sustained VA at 5 years), the DSP risk score demonstrated excellent NPV.
Compared with existing risk assessment tools, the DSP risk score demonstrated marked improvement in risk stratification for DSP patients.5,17 We believe that this model has the potential to be used as part of everyday shared decision-making discussions of the use of primary prevention ICDs in this important, clinically distinct form of arrhythmogenic cardiomyopathy. To facilitate its use, we have made our model available online at www.DSP-risk.com. The DSP risk score represents an important step forward in arrhythmic risk assessment that is tailored not just to a patient’s clinical phenotype but also to his or her underlying genetic substrate.
Desmoplakin-specific risk factors for sustained ventricular arrhythmia
Our findings reinforce the unique nature of DSP-specific VA risk that has been observed by other groups. As in other arrhythmogenic cardiomyopathies, electrophysiologic derangement, including both increased PVC burden and the presence of NSVT on ambulatory monitoring, appears to be important predictors of subsequent VA events. The magnitude of increased VA risk conferred by these findings appears to be very similar across DSP-related disease, ARVC,11 and phospholamban (PLN)-related disease.21 While RV dysfunction is another known risk factor for sustained VA events in ARVC, the magnitude of risk conferred by its presence was even more pronounced in DSP (HR DSP: 6.0; HR ARVC: 1.8).27 Importantly, this risk factor was much less common in DSP patients (moderate to severe RV dysfunction present in 2.6% of patients). Moderate to severe RV dysfunction is thus an infrequent manifestation of DSP,5 and, when paired with the fact that it was identified only in the presence of multiple other VA risk factors, this finding suggests that RV dysfunction may be a marker of more advanced disease [analogous to LV involvement in plakophilin (PKP2)-mediated ARVC]. Unlike for patients with other cardiomyopathies in whom LVEF < 30%–35% usually implies increased VA risk (and, as a result, is often used as a cut-off for ICD implantation decisions), we found that DSP patients had elevated VA risk at even low levels of LV systolic dysfunction. Here, VA risk increased below an LVEF threshold of 50% and remained stably elevated even for more severely reduced LVEFs. As suggested by prior studies demonstrating the predilection of DSP variants to cause LV fibrosis and subsequent systolic dysfunction,4–7 the importance of LV morphologic and functional evaluation is another distinguishing feature of VA risk assessment in these patients. The ARVC risk calculator, which was derived in a predominantly PKP2 variant or gene-elusive patient cohort, does not incorporate any direct assessments of LV involvement (either the presence of LGE or LVEF). This has proven to be a limitation for risk stratification of patients with less common ARVC genotypes in which LV involvement is frequent (e.g. DSP), and multiple external validation studies have demonstrated unreliable performance of the ARVC risk calculator for patients with LV involvement.15,18 Finally, in stark contrast to other cardiomyopathies,28,29 we identified female sex, rather than male sex, as a strong independent risk factor for VA. The underlying mechanism of this association is yet to be fully decoded, but differences in sex hormone levels and history of pregnancy are suspected to play roles.30
Of note, while DSP patients experiencing incident-sustained VA events were more likely to be probands than those who did not, defining proband status for an individual patient can be difficult from a practical standpoint. Family history may be unknown or ambiguous, particularly given the possibility of SCD as the first manifestation of DSP-related disease. In the era of increasingly widespread, community-based genetic testing,31 there is also likely to be a shift towards initial identification of phenotype-negative patients or patients with mild phenotype within carrier families. Proband status was thus purposefully excluded as a predictor during the development of the DSP RISK SCORE. Importantly, the DSP risk score demonstrated well-calibrated risk prediction in both proband and family member subgroups (see Supplementary data online, Figure S3), suggesting that much of the predictive value of proband status is accounted for by the other included risk predictors. Likewise, we did not include LGE as a quantitative predictor due to the practical difficulties in achieving robust adjudication of this variable across sites. Since focal sub-epicardial LGE occurs early in the disease process with subsequent progression to circumferential LGE,4 prioritization of LGE as a continuous variable in future studies may help to further refine risk assessment in these patients.
Value of genotype-tailored arrhythmic risk assessment
While patients with P/LP DSP variants often fulfil the criteria for the diagnosis of ARVC (due to the presence of a P/LP DSP variant fulfilling a major criterion), the poor performance of the ARVC risk calculator within this population highlights the need for genotype-tailored strategies for risk assessment.17,18 These strategies have demonstrated improved risk prediction in other genetic cardiomyopathies,32 including those related to filamin C (FLNC),19 lamin A (LMNA),20 and PLN.21 Development of the DSP risk score is thus part of a larger cultural shift away from basing arrhythmic risk stratification on groups of patients with loosely similar clinical phenotype and towards a ‘gene-first’ approach.33 However, the way in which this new paradigm impacts the use of ICDs in patients with arrhythmogenic cardiomyopathies remains to be seen. Unlike in hypertrophic cardiomyopathy where specific risk cut-offs derived from analogous risk models have been recommended for deciding ICD eligibility,10,34 ICD management in ARVC has historically leveraged risk scores as a tool for guiding shared decision-making conversations,11,12 but regional differences exist.35 Recent guidelines from the European Society of Cardiology (ESC) on VA and the prevention of SCD10 highlight a limited set of high-risk genes in dilated cardiomyopathy (including FLNC, LMNA, and PLN) for which mutations plus LVEF < 50% and an additional risk factor should trigger consideration of ICD implantation. While DSP was not recognized as a high-risk gene in these guidelines, updated ESC guidelines for the management of cardiomyopathies do recognize DSP cardiomyopathy as a high-risk subgroup of both dilated cardiomyopathy and the newly defined phenotype of non-dilated LV cardiomyopathy and suggest consideration of primary prevention ICD in these patients.36 Our findings likewise demonstrate that DSP patients, particularly those with manifestations of disease such as left- or right-sided systolic dysfunction, high PVC burden, or history of NSVT, are at similarly elevated risk for VA compared with these other high-risk genotypes.
The DSP risk score likely represents only the first step in individualized risk assessment for DSP patients. As for ARVC patients, many important modifiers to VA risk will require further study. Athletic activity is strongly associated with worsening risk for VA in ARVC,37,38 but to what extent this remains true in DSP patients remains unclear. Incident arrhythmic risk in ARVC has also been shown to decrease over time,27 but whether this holds true in patients with DSP-related disease requires further study with longitudinal risk modelling. In particular, episodes of myocardial injury are observed in around 9% of DSP patients over the course of their disease, and these episodes are followed by a dramatic increase in subsequent arrhythmic risk.4,5,36,39 The extent to which these episodes parallel the ‘hot phases’ of increased ectopy seen in ARVC,40 reflect rapid periods of otherwise normal disease progression, or represent some other distinct mechanism leading to increased arrhythmic burden is not clear at this time. Finally, specialized imaging such as CMR for assessment of LGE and invasive procedures such as programmed ventricular stimulation for assessment of arrhythmic inducibility41 and ablation42 have demonstrated improved risk stratification in intermediate-risk ARVC patients for whom the decision to implant ICD is ambiguous. While we identified a signal suggesting the potential for high-risk LGE pattern to further stratify VA risk predictions provided by the DSP risk score, due to relatively high missingness in CMR (∼20%), this study was not powered to fully elucidate the role of LGE assessment in DSP patients and further study is required. Although the incidence of sustained VA in low-risk patients with high-risk LGE pattern on CMR was low in the present cohort (only 1 of 51 patients experienced VA), VA was still possible in this important subset of patients. Nuanced clinical judgement that considers high-risk LGE pattern and other potential risk markers not adequately captured by the DSP risk score (e.g. syncope deemed to be very likely arrhythmic) thus remains an essential part of shared decision-making around ICD implantation even in low-risk patients.
Clinical implications
Desmoplakin patients have a higher risk for arrhythmia even relative to patients with other forms of arrhythmogenic cardiomyopathy,15,17,18 but not all DSP patients will go on to have VA events. Distinguishing low-risk patients unlikely to benefit from ICD placement from high-risk patients for whom ICD placement is critical for reducing SCD risk is therefore of paramount importance. The DSP risk score offers prognostic information to help guide decision-making around the placement of primary prevention ICDs. This type of decision is highly personal, however, and patients may have different tolerance for risk based upon their age, sex, lifestyle, or cultural background. For this reason, the DSP risk score assesses VA risk as a continuous variable rather than simply grouping patients into low-, intermediate-, or high-risk categories. The DSP risk score thus seeks to provide clinicians with the data necessary to have more nuanced, data-driven, shared decision-making conversations with their patients.
While there is no concrete risk threshold below which DSP patients should not receive an ICD, our model demonstrated excellent NPV in those patients for whom 5-year VA risk was estimated to be below 5%. In fact, no low-risk patients within the external validation cohort (n = 26) went on to have incident VA events in the 5 years following their initial evaluation. Thus, patients and clinicians may feel reassured that a decision to withhold or defer ICD placement based upon 5-year DSP risk score prediction below 5% is reasonable.
Limitations
Our study population was drawn from the DSP-ERADOS Network, an international collaboration of academic centres from North America, Europe, and Australia. The method for identifying patients with P/LP DSP variants was defined locally by each contributing centre, and there may be between-centre differences in the clinical indications of enrolled patients. These patients were predominantly White, and our results should consequently be extrapolated with caution to patients from other racial or ethnic backgrounds. Our ascertainment from primarily tertiary care settings may also have created a referral bias that could lead to overestimation of VA risk in a community-derived population. This was a retrospective cohort study with variable patient follow-up time. Some patients may have been lost to follow-up prior to an initial VA event. While this study represents the largest cohort of primary prevention DSP patients published to date, statistical power to analyse potential arrhythmic risk predictors was limited to approximately five variables (based on our requirement for ∼10 VA events per added predictor); additional arrhythmic risk predictors may exist and warrant further study. We did not identify a significant association between sustained VA and syncopal events; this could reflect the challenge of retrospective adjudication of the aetiology of syncopal events (e.g. vagal/orthostatic vs. arrhythmic/unexplained). Additionally, family history of SCD was not systematically collected and therefore was not tested as a potential risk factor for sustained VA.
As in similar studies, we used a surrogate composite endpoint that included appropriate ICD therapy to infer risk of SCD. While most clinicians agree that ICD-treated VA represents a severe event, ICD therapies are an imperfect substitute for SCD.43 Further study examining the utilization of primary prevention ICDs in DSP patients and their ultimate impact on adverse clinical arrhythmic outcomes is warranted. Finally, the DSP risk score provides VA risk estimation from the single time of a patient’s initial clinical evaluation. Arrhythmic risk in ARVC is known to change over time,27 and thus DSP patients should be re-evaluated periodically. Longitudinal studies are needed to identify the frequency with which clinical re-evaluation of DSP patients should be performed.
Conclusion
In this study, we present the DSP risk score, a novel model that leverages readily available clinical parameters to generate individualized assessments of VA risk in patients with P/LP DSP variants. The DSP risk score demonstrated good ability to predict incident-sustained VA by 5 years during both model development and subsequent external validation. It has the potential to guide the implantation of primary prevention ICDs in this rare but important form of arrhythmogenic heart disease. This clinical tool thus represents an important step forward in arrhythmic risk assessment that is tailored not just to a patient’s clinical phenotype but also to his or her underlying genetic.
Supplementary Material
Acknowledgements
We are grateful to the DSP patients and families who have made this work possible. For the purpose of open access, the authors have applied a Creative Commons Attribution (CC BY) license to any Author Accepted Manuscript version arising.
Contributor Information
Richard T Carrick, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Alessio Gasperetti, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA; Department of Genetics, University Medical Center Utrecht, University of Utrecht, Utrecht, The Netherlands; Department of Cardiology, University Medical Center Utrecht, Utrecht University, Utrecht, The Netherlands.
Alexandros Protonotarios, Inherited Cardiovascular Diseases Unit, St Bartholomew's Hospital, UCL Institute of Cardiovascular Science, London, UK.
Brittney Murray, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Mikael Laredo, Institut de Cardiologie, Sorbonne Université, AP-HP, IHU-ICAN, Groupe Hospitalier Pitié-Salpêtrière, Institut de Cardiologie, Paris, France.
Iris van der Schaaf, Department of Cardiology, University Medical Center Utrecht, Utrecht University, Utrecht, The Netherlands.
Dennis Dooijes, Department of Cardiology, University Medical Center Utrecht, Utrecht University, Utrecht, The Netherlands.
Petros Syrris, Inherited Cardiovascular Diseases Unit, St Bartholomew's Hospital, UCL Institute of Cardiovascular Science, London, UK.
Douglas Cannie, Inherited Cardiovascular Diseases Unit, St Bartholomew's Hospital, UCL Institute of Cardiovascular Science, London, UK.
Crystal Tichnell, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Nisha A Gilotra, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Chiara Cappelletto, Division of Cardiology, Cardiothoracovascular Department, Azienda Sanitaria Universitaria Giuliano Isontina, University of Trieste, Trieste, Italy.
Kristen Medo, Department of Medicine, Division of Cardiology, University of Colorado Cardiovascular Institute, University of Colorado Anschutz Medical Campus, Aurora, CO, USA.
Ardan M Saguner, Department of Cardiology, Arrhythmia Unit, University Heart Center, University Hospital Zurich, Zurich, Switzerland.
Firat Duru, Department of Cardiology, Arrhythmia Unit, University Heart Center, University Hospital Zurich, Zurich, Switzerland.
Robyn J Hylind, Center for Cardiovascular Genetics, Boston Children’s Hospital, Boston, MA, USA.
Dominic J Abrams, Center for Cardiovascular Genetics, Boston Children’s Hospital, Boston, MA, USA.
Neal K Lakdawala, Center for Advanced Heart Disease, Brigham and Women's Hospital Cardiovascular Medicine, Boston, MA, USA.
Julia Cadrin-Tourigny, Cardiovascular Genetics Center, Montreal Heart Institute, Université de Montréal, Montréal, QC, Canada.
Mattia Targetti, Department of Experimental and Clinical Medicine, University of Florence, Meyer Children Hospital, Careggi University Hospital, Florence, Italy.
Iacopo Olivotto, Department of Experimental and Clinical Medicine, University of Florence, Meyer Children Hospital, Careggi University Hospital, Florence, Italy.
Maddalena Graziosi, Cardiology Unit, Cardiac Thoracic and Vascular Department, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy.
Moniek Cox, Department of Cardiology, University Medical Centre Groningen, Groningen, The Netherlands.
Elena Biagini, Cardiology Unit, Cardiac Thoracic and Vascular Department, IRCCS Azienda Ospedaliero-Universitaria di Bologna, Bologna, Italy.
Philippe Charron, Institut de Cardiologie, Sorbonne Université, AP-HP, IHU-ICAN, Groupe Hospitalier Pitié-Salpêtrière, Institut de Cardiologie, Paris, France.
Paolo Compagnucci, Cardiology and Arrhythmology Clinic, University Hospital ‘Ospedali Riuniti’, Ancona, Italy.
Michela Casella, Cardiology and Arrhythmology Clinic, University Hospital ‘Ospedali Riuniti’, Ancona, Italy; Department of Clinical, Special and Dental Sciences, Marche Polytechnic University, Ancona, Italy.
Giulio Conte, Department of Cardiology, Cardiocentro Ticino Institute, Ente Ospedaliero Cantonale, Lugano, Switzerland.
Claudio Tondo, Department of Clinical Electrophysiology and Cardiac Pacing, Centro Cardiologico Monzino, Cen, IRCCS, University of Milan, Milan, Italy; Department of Biochemical, Surgical and Dentist Sciences, University of Milan, Milan, Italy.
Momina Yazdani, Department of Cardiology, National Heart and Lung Institute and and MRC London Institute of Medical Sciences, London, United Kingdom; Royal Brompton & Harefield Hospitals, Guy’s and St. Thomas’ NHS Foundation Trust, London, UK.
James S Ware, Department of Cardiology, National Heart and Lung Institute and and MRC London Institute of Medical Sciences, London, United Kingdom; Royal Brompton & Harefield Hospitals, Guy’s and St. Thomas’ NHS Foundation Trust, London, UK.
Sanjay K Prasad, Department of Cardiology, National Heart and Lung Institute and and MRC London Institute of Medical Sciences, London, United Kingdom; Royal Brompton & Harefield Hospitals, Guy’s and St. Thomas’ NHS Foundation Trust, London, UK.
Leonardo Calò, Department of Cardiology, Policlinico Casilino, Rome, Italy.
Eric D Smith, Department of Internal Medicine, Division of Cardiovascular Medicine, University of Michigan, Ann Arbor, USA.
Adam S Helms, Department of Internal Medicine, Division of Cardiovascular Medicine, University of Michigan, Ann Arbor, USA.
Sophie Hespe, Centre for Population Genomics, Garvan Institute of Medical Research, UNSW Sydney, Sydney, Australia.
Jodie Ingles, Centre for Population Genomics, Garvan Institute of Medical Research, UNSW Sydney, Sydney, Australia.
Harikrishna Tandri, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Flavie Ader, UF de Cardiogénétique et Myogénétique Moléculaire et Cellulaire, APHP Sorbonne Université, DMU BioGem, 75013 Paris, France; Université Paris Cité, UFR de Pharmacie, UP Biochimie, 75006 Paris, France.
Giovanni Peretto, Department of Cardiac Electrophysiology and Arrhythmology, IRCCS San Raffaele Hospital, Milan, Italy.
Stacey Peters, Department of Cardiology, Royal Melbourne Hospital, Melbourne, Victoria 3050, Australia.
Ari Horton, Department of Cardiology, Royal Melbourne Hospital, Melbourne, Victoria 3050, Australia.
Jessica Yao, Department of Cardiology, Royal Melbourne Hospital, Melbourne, Victoria 3050, Australia.
Eric Schulze-Bahr, Department of Cardiovascular Medicine, Institute for Genetics of Heart Diseases, University Hospital Münster, Münster, Germany.
Sven Dittman, Department of Cardiovascular Medicine, Institute for Genetics of Heart Diseases, University Hospital Münster, Münster, Germany.
Eric D Carruth, Department of Translational Data Science and Informatics, Geisinger, Danville, PA, USA.
Katelyn Young, Department of Translational Data Science and Informatics, Geisinger, Danville, PA, USA.
Maria Qureshi, Department of Translational Data Science and Informatics, Geisinger, Danville, PA, USA.
Chris Haggerty, Department of Translational Data Science and Informatics, Geisinger, Danville, PA, USA; The Heart Institute, Geisinger, Danville, PA, USA.
Victoria N Parikh, Stanford Center for Inherited Cardiovascular Disease, Division of Cardiovascular Medicine, Department of Medicine, Stanford University School of Medicine, Stanford, CA, USA.
Matthew Taylor, Department of Medicine, Division of Cardiology, University of Colorado Cardiovascular Institute, University of Colorado Anschutz Medical Campus, Aurora, CO, USA.
Luisa Mestroni, Department of Medicine, Division of Cardiology, University of Colorado Cardiovascular Institute, University of Colorado Anschutz Medical Campus, Aurora, CO, USA.
Arthur Wilde, Department of Cardiology, Amsterdam UMC Location University of Amsterdam, Meibergdreef 9, Amsterdam, The Netherlands; Amsterdam Cardiovascular Sciences, Heart Failure and Arrhythmias, Amsterdam UMC Location University of Amsterdam, Amsterdam, The Netherlands.
Gianfranco Sinagra, Division of Cardiology, Cardiothoracovascular Department, Azienda Sanitaria Universitaria Giuliano Isontina, University of Trieste, Trieste, Italy.
Marco Merlo, Division of Cardiology, Cardiothoracovascular Department, Azienda Sanitaria Universitaria Giuliano Isontina, University of Trieste, Trieste, Italy.
Estelle Gandjbakhch, Institut de Cardiologie, Sorbonne Université, AP-HP, IHU-ICAN, Groupe Hospitalier Pitié-Salpêtrière, Institut de Cardiologie, Paris, France.
J Peter van Tintelen, Department of Genetics, University Medical Center Utrecht, University of Utrecht, Utrecht, The Netherlands; Netherlands Heart Institute, Utrecht, The Netherlands.
Anneline S J M te Riele, Department of Cardiology, University Medical Center Utrecht, Utrecht University, Utrecht, The Netherlands; Netherlands Heart Institute, Utrecht, The Netherlands.
Perry Elliott, Inherited Cardiovascular Diseases Unit, St Bartholomew's Hospital, UCL Institute of Cardiovascular Science, London, UK.
Hugh Calkins, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Katherine C Wu, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Cynthia A James, Division of Cardiology, Department of Medicine, Johns Hopkins University, 601 North Caroline St., Baltimore, MD 21287, USA.
Supplementary data
Supplementary data are available at European Heart Journal online.
Declarations
Disclosure of Interest
A.M.S. received educational grants through his institution from Abbott, Bayer Healthcare, Biosense Webster, Biotronik, Boston Scientific, BMS/Pfizer, and Medtronic and speaker/advisory board/consulting fees from Bayer Healthcare, Biotronik, Daiichi Sankyo, Medtronic, Novartis, Pfizer, and StrideBio. H.C. is a consultant for Medtronic Inc., Biosense Webster, Pfizer, StrideBio, and Abbott, Tenaya. H.C. receives research support from Boston Scientific. C.T. and C.A.J. receive salary support from this grant. C.A.J. receives research support from Lexeo, Inc. C.T. and A.G. receive salary support from this grant. C.A.J. is a consultant for Pfizer, Lexeo, and StrideBio. J.C.T. is consultant for Tenaya, Lexeo, Bayer. and BMS/Pfizer. N.A.G. is a consultant for Kiniksa Pharmaceuticals. L.M. is a consultant for Tenaya and received research funding from Pfizer, Greenstone Bio, and Tenaya. A.S.J.M.t.R. and J.P.v.T. are consultants for StrideBio and Tenaya. G.S. is a consultant for Novartis, Impulse Dynamics, Novo Nordisk, and Biotronik and a speaker and honoraria for Novartis, Bayer, AstraZeneca, Boston Scientific, Vifor Pharma, Menarini, and Akcea Therapeutics. M.M. is a consultant for Novartis, AstraZeneca, Novo Nordisk, Pfizer, and Vitor Pharma. C.C. is a consultant for Novo Nordisk and AstraZeneca. J.S.W. has received research support from Bristol Myers Squibb and advisory board/consultancy fees from Bristol Myers Squibb, Foresite Labs, Pfizer, and Health Lumen. Authors for whom no specific disclosures are listed have nothing to disclose.
Data Availability
To maintain patient confidentiality, data and study materials have not been made publicly available, but a limited data set may be made available upon reasonable request.
Funding
R.T.C. is funded by the National Institutes of Health (T32HL007227 and L30HL165535) and is a recipient of the Semyon and Janna Friedman Fellowship award. The Zurich ARVC Program is supported by the Georg und Bertha Schwyzer-Winiker Foundation, Baugarten Foundation, USZ Foundation (Dr Wild Grant), Swiss Heart Foundation grant no. FF17019 and FF21073 to A.M.S., Swiss National Science Foundation grant No. 160327 to F.D. J.C.-T. receives salary support as a FRQS research scholar. L.M. and M.R.T. are founded by NIH/NHLBI (R01HL164634 and R01HL147064). A.S.J.M.t.R. is funded by ZonMW (Off Road grant 2021). This work was supported by the Netherlands Cardiovascular Research Initiative with the support of the Dutch Heart Foundation (PREDICT2 2018-30). J.S.W., M.Y., and S.K.P. are supported by the Medical Research Council (UK), British Heart Foundation (RE/18/4/34215 and FS/CRTF/23/24448), the NIHR Imperial College Biomedical Research Centre. the NIHR Royal Brompton Biomedical Research Centre, the Sir Jules Thorn Charitable Trust (21JTA), and Alexander Jansons Myocarditis UK, Rosetrees Trust. The views expressed in this work are those of the authors and not necessarily those of the funders.
Ethical Approval
This study conformed to the Declaration of Helsinki and was approved by local ethics and/or institutional review boards.
Pre-registered Clinical Trial Number
None supplied.
References
- 1. Smith EA, Fuchs E. Defining the interactions between intermediate filaments and desmosomes. J Cell Biol 1998;141:1229–41. 10.1083/jcb.141.5.1229 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 2. Rampazzo A, Nava A, Malacrida S, Beffagna G, Bauce B, Rossi V, et al. . Mutation in human desmoplakin domain binding to plakoglobin causes a dominant form of arrhythmogenic right ventricular cardiomyopathy. Am J Hum Genet 2002;71:1200–6. 10.1086/344208 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 3. Marcus FI, McKenna WJ, Sherrill D, Basso C, Bauce B, Bluemke DA, et al. . Diagnosis of arrhythmogenic right ventricular cardiomyopathy/dysplasia: proposed modification of the task force criteria. Eur Heart J 2010;31:806–14. 10.1093/eurheartj/ehq025 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 4. Smith ED, Lakdawala NK, Papoutsidakis N, Aubert G, Mazzanti A, McCanta AC, et al. . Desmoplakin cardiomyopathy, a fibrotic and inflammatory form of cardiomyopathy distinct from typical dilated or arrhythmogenic right ventricular cardiomyopathy. Circulation 2020;141:1872–84. 10.1161/CIRCULATIONAHA.119.044934 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 5. Gasperetti A, Carrick RT, Protonotarios A, Murray B, Laredo M, van der Schaaf I, et al. . Clinical features and outcomes of 815 patients harboring desmoplakin pathogenic variants: gene-specific evaluation of a distinct clinical entity. Heart Rhythm 2024;21:S20. 10.1016/j.hrthm.2024.03.279 [DOI] [Google Scholar]
- 6. Di Lorenzo F, Marchionni E, Ferradini V, Latini A, Pezzoli L, Martino A, et al. . DSP-related cardiomyopathy as a distinct clinical entity? Emerging evidence from an Italian cohort. Int J Mol Sci 2023;24:2490. 10.3390/ijms24032490 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 7. Wang W, Murray B, Tichnell C, Gilotra NA, Zimmerman SL, Gasperetti A, et al. . Clinical characteristics and risk stratification of desmoplakin cardiomyopathy. Europace 2022;24:268–77. 10.1093/europace/euab183 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 8. Towbin JA, McKenna WJ, Abrams DJ, Ackerman MJ, Calkins H, Darrieux FCC, et al. . 2019 HRS expert consensus statement on evaluation, risk stratification, and management of arrhythmogenic cardiomyopathy. Heart Rhythm 2019;16:e301–72. 10.1016/j.hrthm.2019.05.007 [DOI] [PubMed] [Google Scholar]
- 9. Al-Khatib SM, Stevenson WG, Ackerman MJ, Bryant WJ, Callans DJ, Curtis AB, et al. . 2017 AHA/ACC/HRS guideline for management of patients with ventricular arrhythmias and the prevention of sudden cardiac death: a report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines and the Heart Rhythm Society. J Am Coll Cardiol 2018;72:e91–220. 10.1016/j.hrthm.2017.10.036 [DOI] [PubMed] [Google Scholar]
- 10. Zeppenfeld K, Tfelt-Hansen J, de Riva M, Winkel BG, Behr ER, Blom NA, et al. . 2022 ESC guidelines for the management of patients with ventricular arrhythmias and the prevention of sudden cardiac death. Eur Heart J 2022;43:3997–4126. 10.1093/eurheartj/ehac262 [DOI] [PubMed] [Google Scholar]
- 11. Cadrin-Tourigny J, Bosman LP, Nozza A, Wang W, Tadros R, Bhonsale A, et al. . A new prediction model for ventricular arrhythmias in arrhythmogenic right ventricular cardiomyopathy. Eur Heart J 2019;40:1850–8. 10.1093/eurheartj/ehz103 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 12. Jordà P, Bosman LP, Gasperetti A, Mazzanti A, Gourraud JB, Davies B, et al. . Arrhythmic risk prediction in arrhythmogenic right ventricular cardiomyopathy: external validation of the arrhythmogenic right ventricular cardiomyopathy risk calculator. Eur Heart J 2022;43:3041–52. 10.1093/eurheartj/ehac289 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 13. Zhang N, Wang C, Gasperetti A, Song Y, Niu H, Gu M, et al. . Validation of an arrhythmogenic right ventricular cardiomyopathy risk-prediction model in a Chinese cohort. J Clin Med 2022;11:1973. 10.3390/jcm11071973 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 14. Baudinaud P, Laredo M, Badenco N, Rouanet S, Waintraub X, Duthoit G, et al. . External validation of a risk prediction model for ventricular arrhythmias in arrhythmogenic right ventricular cardiomyopathy. Can J Cardiol 2021;37:1263–6. 10.1016/j.cjca.2021.02.018 [DOI] [PubMed] [Google Scholar]
- 15. Aquaro GD, De Luca A, Cappelletto C, Raimondi F, Bianco F, Botto N, et al. . Prognostic value of magnetic resonance phenotype in patients with arrhythmogenic right ventricular cardiomyopathy. J Am Coll Cardiol 2020;75:2753–65. 10.1016/j.jacc.2020.04.023 [DOI] [PubMed] [Google Scholar]
- 16. Casella M, Gasperetti A, Gaetano F, Busana M, Sommariva E, Catto V, et al. . Long-term follow-up analysis of a highly characterized arrhythmogenic cardiomyopathy cohort with classical and non-classical phenotypes-a real-world assessment of a novel prediction model: does the subtype really matter. Europace 2020;22:797–805. 10.1093/europace/euz352 [DOI] [PubMed] [Google Scholar]
- 17. Protonotarios A, Bariani R, Cappelletto C, Pavlou M, García-García A, Cipriani A, et al. . Importance of genotype for risk stratification in arrhythmogenic right ventricular cardiomyopathy using the 2019 ARVC risk calculator. Eur Heart J 2022;43:3053–67. 10.1093/eurheartj/ehac235 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 18. Gasperetti A, Carrick R, Protonotarios A, Laredo M, van der Schaaf I, Syrris P, et al. . Long-term arrhythmic follow-up and risk stratification of patients with desmoplakin-associated arrhythmogenic right ventricular cardiomyopathy. JACC Adv 2024;3:100832. 10.1016/j.jacadv.2024.100832 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 19. Gigli M, Stolfo D, Graw SL, Merlo M, Gregorio C, Nee Chen S, et al. . Phenotypic expression, natural history, and risk stratification of cardiomyopathy caused by filamin C truncating variants. Circulation 2021;144:1600–11. 10.1161/CIRCULATIONAHA.121.053521 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 20. Nishiuchi S, Makiyama T, Aiba T, Nakajima K, Hirose S, Kohjitani H, et al. . Gene-based risk stratification for cardiac disorders in LMNA mutation carriers. Circ Cardiovasc Genet 2017;10:e001603. 10.1161/CIRCGENETICS.116.001603 [DOI] [PubMed] [Google Scholar]
- 21. Verstraelen TE, van Lint FHM, Bosman LP, de Brouwer R, Proost VM, Abeln BGS, et al. . Prediction of ventricular arrhythmia in phospholamban p.Arg14del mutation carriers-reaching the frontiers of individual risk prediction. Eur Heart J 2021;42:2842–50. 10.1093/eurheartj/ehab294 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 22. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. . Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. Genet Med 2015;17:405–24. 10.1038/gim.2015.30 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 23. van Lint FHM, Murray B, Tichnell C, Zwart R, Amat N, Lekanne Deprez RH, et al. . Arrhythmogenic right ventricular cardiomyopathy-associated desmosomal variants are rarely de novo. Circ Genom Precis Med 2019;12:e002467. 10.1161/CIRCGEN.119.002467 [DOI] [PubMed] [Google Scholar]
- 24. Di Marco A, Brown PF, Bradley J, Nucifora G, Claver E, de Frutos F, et al. . Improved risk stratification for ventricular arrhythmias and sudden death in patients with nonischemic dilated cardiomyopathy. J Am Coll Cardiol 2021;77:2890–905. 10.1016/j.jacc.2021.04.030 [DOI] [PubMed] [Google Scholar]
- 25. Collins GS, Reitsma JB, Altman DG, Moons KGM. Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement. BMJ 2015;350:g7594. 10.1136/bmj.g7594 [DOI] [PubMed] [Google Scholar]
- 26. Heagerty PJ, Lumley T, Pepe MS. Time-dependent ROC curves for censored survival data and a diagnostic marker. Biometrics 2000;56:337–44. 10.1111/j.0006-341X.2000.00337.x [DOI] [PubMed] [Google Scholar]
- 27. Carrick RT, Te Riele ASJM, Gasperetti A, Bosman L, Muller SA, Pendleton C, et al. . Longitudinal prediction of ventricular arrhythmic risk in patients with arrhythmogenic right ventricular cardiomyopathy. Circ Arrhythm Electrophysiol 2022;15:e011207. 10.1161/CIRCEP.122.011207 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 28. Bosman LP, Sammani A, James CA, Cadrin-Tourigny J, Calkins H, van Tintelen JP, et al. . Predicting arrhythmic risk in arrhythmogenic right ventricular cardiomyopathy: a systematic review and meta-analysis. Heart Rhythm 2018;15:1097–107. 10.1016/j.hrthm.2018.01.031 [DOI] [PubMed] [Google Scholar]
- 29. Sammani A, Kayvanpour E, Bosman LP, Sedaghat-Hamedani F, Proctor T, Gi WT, et al. . Predicting sustained ventricular arrhythmias in dilated cardiomyopathy: a meta-analysis and systematic review. ESC Heart Fail 2020;7:1430–41. 10.1002/ehf2.12689 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 30. Akdis D, Saguner AM, Shah K, Wei C, Medeiros-Domingo A, von Eckardstein A, et al. . Sex hormones affect outcome in arrhythmogenic right ventricular cardiomyopathy/dysplasia: from a stem cell derived cardiomyocyte-based model to clinical biomarkers of disease outcome. Eur Heart J 2017;38:1498–508. 10.1093/eurheartj/ehx011 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 31. Carruth ED, Beer D, Alsaid A, Schwartz MLB, McMinn M, Kelly MA, et al. . Clinical findings and diagnostic yield of arrhythmogenic cardiomyopathy through genomic screening of pathogenic or likely pathogenic desmosome gene variants. Circ Genom Precis Med 2021;14:e003302. 10.1161/CIRCGEN.120.003302 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 32. Paldino A, Dal Ferro M, Stolfo D, Gandin I, Medo K, Graw S, et al. . Prognostic prediction of genotype vs phenotype in genetic cardiomyopathies. J Am Coll Cardiol 2022;80:1981–94. 10.1016/j.jacc.2022.08.804 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 33. James CA, Gasperetti A. Genotype-based risk stratification can outperform phenotype-based practice for inherited cardiomyopathies: not all paths are equal. J Am Coll Cardiol 2022;80:1995–7. 10.1016/j.jacc.2022.09.018 [DOI] [PubMed] [Google Scholar]
- 34. Elliott PM, Anastasakis A, Borger MA, Borggrefe M, Cecchi F, Charron P, et al. . 2014 ESC guidelines on diagnosis and management of hypertrophic cardiomyopathy: the task force for the diagnosis and management of hypertrophic cardiomyopathy of the European Society of Cardiology (ESC). Eur Heart J 2014;35:2733–79. 10.1093/eurheartj/ehu284 [DOI] [PubMed] [Google Scholar]
- 35. Carrick RT, De Marco C, Gasperetti A, Bosman LP, Gourraud JB, Trancuccio A, et al. . Implantable cardioverter defibrillator use in arrhythmogenic right ventricular cardiomyopathy in North America and Europe. Eur Heart J 2024;45:538–48. 10.1093/eurheartj/ehad799 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 36. Arbelo E, Protonotarios A, Gimeno JR, Arbustini E, Barriales-Villa R, Basso C, et al. . 2023 ESC guidelines for the management of cardiomyopathies. Eur Heart J 2023;44:3503–626. 10.1093/eurheartj/ehad194 [DOI] [PubMed] [Google Scholar]
- 37. Gasperetti A, Dello Russo A, Busana M, Dessanai M, Pizzamiglio F, Saguner AM, et al. . Novel risk calculator performance in athletes with arrhythmogenic right ventricular cardiomyopathy. Heart Rhythm 2020;17:1251–9. 10.1016/j.hrthm.2020.03.007 [DOI] [PubMed] [Google Scholar]
- 38. James CA, Bhonsale A, Tichnell C, Murray B, Russell SD, Tandri H, et al. . Exercise increases age-related penetrance and arrhythmic risk in arrhythmogenic right ventricular dysplasia/cardiomyopathy-associated desmosomal mutation carriers. J Am Coll Cardiol 2013;62:1290–7. 10.1016/j.jacc.2013.06.033 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 39. Poller W, Haas J, Klingel K, Kühnisch J, Gast M, Kaya Z, et al. . Familial recurrent myocarditis triggered by exercise in patients with a truncating variant of the desmoplakin gene. J Am Heart Assoc 2020;9:e015289. 10.1161/JAHA.119.015289 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 40. Gasperetti A, Cappelletto C, Carrick R, Targetti M, Tichnell C, Martino A, et al. . Association of premature ventricular contraction burden on serial Holter monitoring with arrhythmic risk in patients with arrhythmogenic right ventricular cardiomyopathy. JAMA Cardiol 2022;7:378–85. 10.1001/jamacardio.2021.6016 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 41. Gasperetti A, Carrick RT, Costa S, Compagnucci P, Bosman LP, Chivulescu M, et al. . Programmed ventricular stimulation as an additional primary prevention risk stratification tool in arrhythmogenic right ventricular cardiomyopathy: a multinational study. Circulation 2022;146:1434–43. 10.1161/CIRCULATIONAHA.122.060866 [DOI] [PMC free article] [PubMed] [Google Scholar]
- 42. Gasperetti A, Peretto G, Muller SA, Hasegawa K, Compagnucci P, Casella M, et al. . Catheter ablation for ventricular tachycardia in patients with desmoplakin cardiomyopathy. JACC Clin Electrophysiol 2024;10:487–98. 10.1016/j.jacep.2023.11.017 [DOI] [PubMed] [Google Scholar]
- 43. Ellenbogen KA, Levine JH, Berger RD, Daubert JP, Winters SL, Greenstein E, et al. . Are implantable cardioverter defibrillator shocks a surrogate for sudden cardiac death in patients with nonischemic cardiomyopathy? Circulation 2006;113:776–82. 10.1161/CIRCULATIONAHA.105.561571 [DOI] [PubMed] [Google Scholar]
Associated Data
This section collects any data citations, data availability statements, or supplementary materials included in this article.
Supplementary Materials
Data Availability Statement
To maintain patient confidentiality, data and study materials have not been made publicly available, but a limited data set may be made available upon reasonable request.