Disease Markers
Volume 2020 (2020), Article ID 8847635, 7 pages
https://doi.org/10.1155/2020/8847635

Research Article

A Novel Gene Signature to Predict Survival Time and Incident Ventricular Arrhythmias in Patients with Dilated Cardiomyopathy

Chenliang Ge1,2 and Yan He¹

¹The First Affiliated Hospital of Guangxi Medical University, Nanning 530021, China
²The First Affiliated Hospital of University of South China, Hengyang 421001 ., China

Correspondence should be addressed to Yan He; hyxjwxy@126.com

Received 28 June 2020; Accepted 1 September 2020; Published 15 September 2020

Academic Editor: Ibadete Bytyçi

Copyright © 2020 Chenliang Ge and Yan He. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The mortality in nonischaemic dilated cardiomyopathy (NIDCM) patients is still at a high level; sudden death in NIDCM can be caused by ventricular tachycardia. It is necessary to explore the pathogenesis of ventricular arrhythmias (VA) in NIDCM. Differentially expressed genes (DEGs) were identified by comparing the gene expression of NIDCM patients with or without VA in the gene expression profile of GSE135055. A total of 228 DEGs were obtained, and 3 genes were screened out to be significantly related to the survival time of NIDCM patients. We established a prediction model on two-gene (TOMM22, PPP2R5A) signature for the survival time of NIDCM patients. The area under the curve (AUC) was 0.75 calculated by the ROC curve analysis. These risk genes are probably new targets for exploring the pathogenesis of NIDCM with VA; the prediction model for survival time and incident ventricular arrhythmias is useful in clinical decision making for individual treatment.

1. Introduction

Nonischaemic dilated cardiomyopathy (NIDCM) is one of the most common inherited cardiomyopathy and is considered to be one of the main causes of heart failure and sudden cardiac death. Heart transplantation is usually needed [1, 2]. One-third of patients with NIDCM may represent arrhythmogenic phenotypes and have an increased risk of arrhythmias during follow-up; ventricular arrhythmia (VA) is a major cause of clinical deterioration and demise in patients with NIDCM [3–6]. In order to prevent these situations, it is recommended that prophylactic intervention with implantable cardioverter-defibrillators (ICD) in patients with heart failure and left ventricular ejection in current international guidelines [7, 8]. Although NIDCM is characterized by genetic and clinical heterogeneity, few studies have explored the pathogenesis of VA in patients with NIDCM. In order to reveal the inherent molecular mechanism, NIDCM patients were divided into sinus rhythm (SR) group and VA group and differentially expressed genes (DEGs) were identified via bioinformatic methods.

2. Materials and Methods

2.1. Data Source

We download the microarray expression dataset (GSE135055) from the Gene Expression Omnibus (GEO) database [9]. The multilevel transcriptional data of GSE135055 were generated from the heart tissues of 21 heart failure (HF) patients and 9 healthy donors. Among these HF patients, there are 18 NIDCM patients, during which 6 patients suffered VA, including ventricular premature beats, ventricular tachycardia, and ventricular fibrillation [10]. The characteristics of these patients are listed in Table 1 (supplemental file Table S1). We divided the NIDCM patients into two groups, with 12 patients in the sinus rhythm (SR) group and 6 patients in the VA group. We used the “limma” package [11] to quantile the datasets and screened out DEGs between the SR group and VA group with 5 and were considered statistically significant. A volcano map and heat map were used to show the DEGs.

Table 1: Baseline parameters in NIDCM patients with and without ventricular arrhythmia.

Clinical features
SR (12)
VA (6)
value

Age (y)
30.33 (16.87)
39.33 (16.17)
0.296

Male sex ( (%))
7 (58.3)
3 (50.0)
0.563

LA (mm)
41.08 (13.53)
35.33 (17.78)
0.454

LVEDD (mm)
68.00 (11.72)
58.33 (28.74)
0.318

LVEF (%)
27.00 (11.50)
24.00 (13.74)
0.631

AF (%)
2 (16.7)
2 (33.3)
0.407

NYHA functional class

        II
2 (16.7)
0 (0.0)
0.806

        III
4 (33.3)
2 (33.3)

        IV
6 (50.0)
4 (66.7)

Diabetes mellitus (%)
1 (8.3)
1 (16.7)
0.569

Smoking history (%)
2 (16.7)
1 (16.7)
0.755

Drug therapy

        Amiodarone (%)
1 (8.3)
2 (33.3)
0.245

        β-Blockers (%)
8 (66.7)
4 (66.7)
0.691

        Digoxin (%)
10 (83.3)
4 (66.7)
0.407

        ACEI/ARB (%)
10 (83.3)
4 (66.7)
0.407

        CCB (%)
2 (16.7)
1 (16.7)
0.755

        Diuretic (%)
11 (91.7)
5 (83.3)
0.569

Values are , (%); LA: left atrium; LVEDD: left ventricular diastolic diameter; LVEF: left ventricular ejection fraction; AF: Atrial fibrillation; NYHA: New York Heart Association; ACEI: angiotensin-converting enzyme inhibitors; ARB: angiotensin receptor blocker; CCB: calcium channel blocker.

2.2. Construction of Survival Time Prediction Model

17 NIDCM patients (6 patients suffered VA included) with complete survival time (time period of each patient from symptoms to heart transplantation) information were selected for the study. By using univariate Cox regression analysis, we screened out DEG significantly related to survival time () [12]. Then, a gene signature and a prediction model for survival time were constructed via multivariate cox analysis, survival analysis, and random survival forest algorithm.

2.3. Evaluation of Survival Time Prediction Model

Based on the regression coefficient and the expression value of each selected gene obtained by the multivariate Cox regression model, the risk score of each patient was calculated, then we separated 17 patients into high-risk and low-risk groups using the median risk score as the cutoff. A high-risk score indicates poor survival time for the NIDCM patients. The accuracy of the prediction model was evaluated by time-dependent ROC analysis.

2.4. Enrichment Analysis by Metascape

The pathway enrichment of DEGs was analyzed using Metascape. With the species limited to “Homo sapiens,” functional enrichment analysis was performed based on Gene Ontology, Kyoto Encyclopedia of Genes and Genomes, Reactive pathways, and Canonical Pathways [13]. All genes in the genome have been used as the enrichment background. Terms with a value < 0.01, a minimum count of 3, and an enrichment are collected and grouped into clusters based on their membership similarities.

3. Results

3.1. Identification of Differentially Expressed Genes (DEGs)

A total of 228 DEGs were obtained via the microarray data analysis by Limma, including 86 upregulated and 142 downregulated DEG (supplemental file Table S2). A DEG expression volcano map and heat map were shown in Figure 1.

(a)

(b)

Figure 1: DEGs between SR and VA. (a) Volcano map of DEGs between SR and VA. Red represents upregulated differential genes, black represents no significant difference genes, and green represents downregulated differential genes. (b) Heat map of all DEGs between SR and VA. Each column represents a tissue sample, and each row represents a DEG. The gradual color change from green to red indicates the changing process from downregulation to upregulation.

3.2. Identification of Survival-Related DEGs and Construction of Survival Time Prediction Model

3 DEGs were screened out to be significantly related to the survival time of NIDCM patients via univariate Cox regression analysis including TOMM22, TMSB4XP8, and PPP2R5A () (Table 2). The 3 genes were used to establish a model for predicting survival time of NIDCM patients. Multivariate Cox proportional hazard regression analysis was performed to confirm the optimal model and a forest map was shown in Figure 2. PPP2R5A and TOMM22 were identified as risk genes in the survival time model. We confirmed PPP2R5A and TOMM22 as high-risk genes in the survival time model which are negatively correlated with patient prognosis (Table 3).

Table 2: DEGs significantly related to survival time of NIDCM patients by univariate Cox regression analysis.

Gene symbol
HR
HR.95L
HR.95H
value

TOMM22
1.42163
1.063455
1.900439
0.017531

TMSB4XP8
1.119692
1.019037
1.230289
0.018656

PPP2R5A
1.231818
1.009141
1.503629
0.040422

Figure 2: Forrest plot of the multivariate Cox regression analysis in NIDCM.

Table 3: Characteristics of risk DEGs in the prognostic risk models.

Gene symbol
Coef
HR
HR.95L
HR.95H
value

PPP2R5A
0.198883
1.220039
0.999961
1.488553
0.050045

TOMM22
0.350995
1.42048
1.053479
1.915333
0.021358

3.3. Testing the Survival Time Prediction Model

The prognostic risk scores for the 17 NIDCM patients were calculated using the gene expression values and the regression coefficients of the risk genes according to the prediction model: . Patients were subdivided into high and low-risk groups for survival time based on the median risk scores. Survival analysis revealed that high-risk scores were significantly related to poor survival time (symptoms to heart transplantation). The 2-year and 5-year survival rates for the high-risk patients were 37.5% and 25%, whereas the 8-year and 11-year survival rates for the low-risk patients were 55.6% and 44.4%. We then measured the predictive performance of the prognostic risk models using the time-dependent receiver operating characteristic (ROC) curves. The area under the curve (AUC) was 0.75, which indicated superior predictive accuracy in NIDCM patients for survival time (Figure 3).

(a)

(b)

(c)

Figure 3: Validation of the survival time prognostic risk models in NIDCM patients. (a) Survival analysis of symptoms to heart transplantation in the high-risk (red line) and low-risk (blue line) NIDCM patients. (b) Time-dependent ROC curve analyses show AUC values for time from symptoms to heart transplantation in NIDCM patients. (c) Expression levels of risk genes in the high-risk and low-risk patients.

3.4. Function Enrichment Analysis by Metascape

The Metascape function enrichment results for DEGs including 204 Gene Ontology terms, 53 Reactome pathways, and 6 Kyoto Encyclopedia of Genes and Genomes pathways. The risk gene PPP2R5A is mainly involved in Influenza Infection, protein localization to membrane, PID IGF1 PATHWAY, CTLA4 inhibitory signaling, Platelet activation, signaling and aggregation, and Hemostasis. TOMM22 is mainly enriched in Influenza Infection, establishment of protein localization to organelle, establishment of protein localization to membrane, protein targeting, protein localization to membrane, mitochondrion organization, autophagy, process utilizing autophagic mechanism, macroautophagy, regulation of mRNA metabolic process, and mitochondrial membrane organization (Figure 4).

(a)

(b)

Figure 4: Enrichment analysis by Metascape. (a) Enriched terms of PPP2R5A. (b) Enriched terms of TOMM22.

4. Discussion

Due to high incidence, poor therapeutic effect, and prognosis, NIDCM is a major health problem that threatens human health, it is necessary to explore its pathogenesis and establish an accurate prediction prognosis model [14–16]. NIDCM is characterized as a genetically determined disease; the genetic basis of NIDCM highlights the importance of screening biomarkers for diagnosis and prognosis [17–19]. Although the survival time of patients with NIDCM is quite different [20, 21], there are few methods to predict the survival time after being diagnosed as NIDCM. We first established a novel two-gene signature to predict survival in patients with NIDCM. To our knowledge, the two-gene signature related prognostic model for NIDCM has not been reported previously. The risk score was based on mRNA expression but not somatic mutations or methylation status of only two prognostic genes, which could be more routine and cost-effective in practice.

VA is common in patients with NIDCM. Ventricular fibrillation and cardiac arrest are important reasons for death in patients with NIDCM [22, 23], while there have been few studies focusing solely on VA in NIDCM. Although sudden cardiac death rates have decreased, implantation of ICD for primary prevention in NIDCM patients does not provide overall survival benefits [24]. Establishing a prediction model for incident VA in patients with NIDCM helps make personalized treatment to improve the survival, which can be used as a reference index of whether patients need to be implanted with ICD or not. In this study, we established a novel two-gene signature (including TOMM22 and PPP2R5A) to predict the occurrence of VA in patients. Consistent with our research, some studies support PPP2R5A as a novel target for the treatment of arrhythmia. PPP2R5A, protein phosphatase 2 regulatory subunit B’alpha, encodes an alpha isoform of the regulatory subunit B56 subfamily. Downregulated B56α myocytes are insensitive to isoproterenol-induced induction of arrhythmogenic Na⁺ channel late component. Voltage-gated Na⁺ channel 1.5 is critical for normal cardiac excitability, PP2A-B56α complex proved to interact with the primary cardiac voltage-gated Na⁺ channel 1.5 [25]. TOMM22, Translocase of outer mitochondrial membrane 22, is responsible for the recognition and translocation of synthesized mitochondrial precursor proteins; phosphorylation of TOMM22 is a critical switch for mitophagy [26, 27]. Combined with the functional enrichment analysis results above, TOMM22 may regulate mitochondrial autophagy in NIDCM, which needs to be confirmed by further research.

Tachycardia-induced cardiomyopathy (TIC) is characterized by diverse tachyarrhythmias, including supraventricular arrhythmias (such as atrial fibrillation) and ventricular arrhythmias [28, 29]. PPP2R5A and TOMM22 show significantly different expression levels between the VA group and SR group in NIDCM patients, which suggests they may be related to the pathogenesis of TIC. More study needs to be done with PPP2R5A and TOMM22 to determine its exact role in the pathogenesis of TIC.

However, there are some limitations in our research. We studied a small number of patients, although the ROC value (0.75) indicated superior predictive accuracy of the prediction model for survival time in NIDCM patients. In future research, we will carry out a large sample, randomized, and follow-up studies. Besides, the time period of each patient from symptoms to heart transplantation still can represent the prognosis of patients with NIDCM to some extent, overall survival data should be more accurate, and we will carry out patients’ long-term follow-up to collect more detail clinical information.

5. Conclusions

Our study identified 3 new genes that were significantly related to the survival time of NIDCM and established a novel two-gene signature to predict survival time and incident VA of NIDCM. Although the sample size of our study is relatively small, these risk genes are probably new targets for exploring the pathogenesis of NIDCM with VA, and the prediction model is useful in clinical decision making for individual treatment.

Abbreviations

NIDCM: | Nonischaemic dilated cardiomyopathy
VA: | Ventricular arrhythmias
GEO: | Gene Expression Omnibus
DEGs: | Differentially expressed genes
SR: | Sinus rhythm
ROCs: | Receiver operating characteristics
AUC: | Area under the curve
HF: | Heart failure.

Data Availability

The dataset with mRNA expression profiling was taken from GEO: GSE135055. The associated corresponding clinical information of patients in this study was downloaded from BMC Medicine (doi:10.1186/s12916-019-1469-4).

Ethical Approval

The databases are publicly available and open access, and the present study followed the data access policy and publishing guidelines of these databases. Therefore, no local ethics committee is required to approve this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

Yan He did the conception and design; Yan He did the administrative support; Chenliang Ge did the provision of study materials or patients; Chenliang Ge did the collection and assembly of data; All authors did the data analysis and interpretation; All authors did the manuscript writing. All authors did the final approval of the manuscript.

Acknowledgments

This research was supported by the National Natural Science Foundation of China (No. 81760061) and the Open Project of Guangxi Key Laboratory of Precision Medicine in Cardio-cerebrovascular Diseases Control and Prevention (GXXNXG202005).

Supplementary Materials

Supplementary 1. Supplementary Table 1: summary of the basic information of DCM patients included in the study.

Supplementary 2. Supplementary Table 2: differentially expressed genes.

References

Y. M. Pinto, P. M. Elliott, E. Arbustini et al., “Proposal for a revised definition of dilated cardiomyopathy, hypokinetic non-dilated cardiomyopathy, and its implications for clinical practice: a position statement of the ESC working group on myocardial and pericardial diseases,” European Heart Journal, vol. 37, no. 23, pp. 1850–1858, 2016.

A. G. Japp, A. Gulati, S. A. Cook, M. R. Cowie, and S. K. Prasad, “The diagnosis and evaluation of dilated cardiomyopathy,” Journal of the American College of Cardiology, vol. 67, no. 25, pp. 2996–3010, 2016.

A. Spezzacatene, G. Sinagra, M. Merlo et al., “Arrhythmogenic phenotype in dilated cardiomyopathy: natural history and predictors of life-threatening arrhythmias,” Journal of the American Heart Association, vol. 4, article e002149, 2015.

F. Sedaghat-Hamedani, J. Haas, F. Zhu et al., “Clinical genetics and outcome of left ventricular non-compaction cardiomyopathy,” European Heart Journal, vol. 38, pp. 3449–3460, 2019.

M. Akhtar and P. M. Elliott, “Risk stratification for sudden cardiac death in non-Ischaemic dilated cardiomyopathy,” Current Cardiology Reports, vol. 21, no. 12, p. 155, 2019.

M. Ebert, S. Richter, B. Dinov, K. Zeppenfeld, and G. Hindricks, “Evaluation and management of ventricular tachycardia in patients with dilated cardiomyopathy,” Heart Rhythm, vol. 16, no. 4, pp. 624–631, 2019.

C. W. Yancy, B. Jessup, M. Fau-Bozkurt et al., “2013 ACCF/AHA guideline for the management of heart failure: executive Summary,” Circulation, vol. 128, no. 16, pp. 1810–1852, 2013.

C. P. S. Fau-Blomström-Lundqvist and C. Blomström-Lundqvist, “European Society of Cardiology Guidelines for the management of patients with ventricular arrhythmias and the prevention of sudden cardiac death summarized by co-chairs,” European Heart Journal, vol. 36, no. 2015, pp. 2757–2759, 2015.

T. Barrett, P. W. S. Fau-Ledoux, C. L. P. Fau-Evangelista et al., “NCBI GEO: archive for functional genomics data sets--update,” Nucleic Acids Research, vol. 41, pp. D991–D995, 2013.

X. Hua, Y. Y. Wang, P. Jia et al., “Multi-level transcriptome sequencing identifies COL1A1 as a candidate marker in human heart failure progression,” BMC Medicine, vol. 18, no. 1, p. 2, 2020.

M. E. Ritchie, B. Phipson, D. Wu et al., “Limma powers differential expression analyses for RNA-sequencing and microarray studies,” Nucleic Acids Research, vol. 43, article e47, 2015.

B. Wan, B. Liu, G. Yu, Y. Huang, and C. Lv, “Differentially expressed autophagy-related genes are potential prognostic and diagnostic biomarkers in clear-cell renal cell carcinoma,” Aging (Albany NY), vol. 11, no. 20, pp. 9025–9042, 2019.

Y. Zhou, B. Zhou, L. Pache et al., “Metascape provides a biologist-oriented resource for the analysis of systems-level datasets,” Nature Communications, vol. 10, no. 1, p. 1523, 2019.

M. Merlo, B. P. A. Fau-Pinamonti, D. P. B. Fau-Stolfo et al., “Long-term prognostic impact of therapeutic strategies in patients with idiopathic dilated cardiomyopathy: changing mortality over the last 30 years,” European Journal of Heart Failure, vol. 16, no. 3, pp. 317–324, 2014.

M. Merlo, A. Cannatà, C. Pio Loco et al., “Contemporary survival trends and aetiological characterization in non-ischaemic dilated cardiomyopathy,” European Journal of Heart Failure, 2020.

A. A.-O. Cannatà, G. De Angelis, A. Boscutti et al., “Arrhythmic risk stratification in non-ischaemic dilated cardiomyopathy beyond ejection fraction,” Heart Failure Clinics, vol. 106, pp. 656–664, 2020.

B. P. Halliday, J. G. F. Cleland, J. J. Goldberger, and S. K. Prasad, “Personalizing risk stratification for sudden death in dilated cardiomyopathy: the past, present, and future,” Circulation, vol. 136, no. 2, pp. 215–231, 2017.

P. Charron, E. A. M. Fau-Arbustini, C. A. E. Fau-Basso et al., “Genetic counselling and testing in cardiomyopathies: a position statement of the European Society of Cardiology Working Group on Myocardial and Pericardial Diseases,” European Heart Journal, vol. 31, pp. 2715–2726., 2010.

R. E. Hershberger, M. M. Givertz, C. Y. Ho et al., “Genetic evaluation of cardiomyopathy-A Heart Failure Society of America practice guideline,” Journal of Cardiac Failure, vol. 24, no. 5, pp. 281–302, 2018.

D. A.-O. Reichart, C. A.-O. Magnussen, T. Zeller, and S. Blankenberg, “Dilated cardiomyopathy: from epidemiologic to genetic phenotypes: a translational review of current literature,” Journal of Internal Medicine, vol. 286, no. 4, pp. 362–372, 2019.

R. G. Weintraub, C. Semsarian, and P. Macdonald, “Dilated cardiomyopathy,” The Lancet, vol. 390, no. 10092, pp. 400–414, 2017.

S. R. Piers, K. Everaerts, R. J. van der Geest et al., “Myocardial scar predicts monomorphic ventricular tachycardia but not polymorphic ventricular tachycardia or ventricular fibrillation in nonischemic dilated cardiomyopathy,” Heart Rhythm, vol. 12, no. 10, pp. 2106–2114, 2015.

D. Stolfo, N. Ceschia, M. Zecchin et al., “Arrhythmic risk stratification in patients with idiopathic dilated cardiomyopathy,” The American Journal of Cardiology, vol. 121, no. 12, pp. 1601–1609, 2018.

L. Køber, J. J. Thune, J. C. Nielsen et al., “Defibrillator implantation in patients with nonischemic systolic heart failure,” The New England Journal of Medicine, vol. 375, no. 13, pp. 1221–1230, 2016.

M. El Refaey, H. Musa, N. P. Murphy et al., “Protein phosphatase 2A regulates cardiac Na(+) channels,” Circulation Research, vol. 124, no. 5, pp. 737–746, 2019.

B. Kravic, A. B. Harbauer, V. Romanello et al., “In mammalian skeletal muscle, phosphorylation of TOMM22 by protein kinase CSNK2/CK2 controls mitophagy,” Autophagy, vol. 14, no. 2, pp. 311–335, 2018.

G. Bertolin, M. F.-M. R. Fau-Jacoupy, S. J. M. Fau-Traver et al., “The TOMM machinery is a molecular switch in PINK1 and PARK2/PARKIN-dependent mitochondrial clearance,” Autophagy, vol. 9, pp. 1801–1817, 2013.

A. Khasnis, G. J. K. Fau-Abela, S. A. G. Fau-Veerareddy, V. V. S. Fau-Reddy, R. R. V. Fau-Thakur, and R. Thakur, “Tachycardia-induced cardiomyopathy: a review of literature,” Pacing and Clinical Electrophysiology, vol. 28, no. 7, pp. 710–721, 2005.

D. A.-O. Kim, S. A.-O. Kim, and K. A.-O. Ryu, “Tachycardia induced cardiomyopathy,” Vnitr̆ní Lékar̆ství, vol. 61, pp. 56–59, 2015.