Research Paper Volume 16, Issue 13 pp 10860—10867

Identification of differentially expressed genes to predict the risk of heart failure in older patients with hypertrophic cardiomyopathy

Hao Dong1, *, , Chufan Yin2, *, , Dongping Xiao3, , Yong Tang1, ,

  • 1 Department of Cardiology, The Second Hospital of Nanjing, Nanjing University of Chinese Medicine, Nanjing 210000, China
  • 2 Education Section, The Second Hospital of Nanjing, Nanjing University of Chinese Medicine, Nanjing 210000, China
  • 3 Department of Cardiology, The First Hospital of Nanchang, Nanchang 330000, China
* Equal contribution

Received: January 31, 2024       Accepted: May 29, 2024       Published: June 20, 2024      

https://doi.org/10.18632/aging.205956
How to Cite

Copyright: © 2024 Dong et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Aim: Hypertrophic cardiomyopathy (HCM) is a common heart disease. Old people with HCM are at high risk of heart failure (HF). This study aimed to identify differentially expressed genes (DEGs) to evaluate the risk of HF in older patients with HCM.

Methods: GSE89714 and GSE116250 were downloaded from Gene Expression Omnibus (GEO) database, and DEGs were identified by using limma R package with P < 0.05 and logFC> 1 as cut off. Protein-protein interaction (PPI) network, Genome Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed for the identified DEGs. NetworkAnalyst online tool was applied for Gene Set Enrichment Analysis (GSEA) analysis.

Results: We identified 124 overlap DEGs from the 2 datasets. PPI network showed that COL1A1, COL3A1, COL1A2, BGN, COL5A1, LUM, TGFB2, FMOD, ASPN, and COL14A1 were the top ten genes related to HCM and HF compared with control. Functional and pathway analyses showed that the overlap genes were mainly related to ECM-receptor interaction, ECM organization, Focal adhesion, PI3K-Akt signaling, TGF-beta signaling, and Platelet activation signaling and aggregation. Among the overlap genes, COL5A1 and LUM were significantly upregulated, while TGFB2, FMOD, ASPN, and COL14A1 were significantly downregulated in HF dataset compared with HCM dataset.

Conclusions: Bioinformatics-based analysis revealed potential genes associated with HCM and HF, which could be utilized to evaluate the risk of HF in older patients with HCM.

Introduction

Hypertrophic cardiomyopathy (HCM) is a common heart disease characterized by the hypertrophy of left ventricle (LV) [1]. Older patients with HCM have a high risk of heart failure (HF) and sudden cardiac death (SCD) [2].

The development of HF in HCM patients remains unclear [3]. In most patients with HCM, HF shows the phenotype of HF with preserved ejection fraction (HFpEF), but few patients with HCM develop HF with reduced ejection fraction (HFrEF) at a late stage [4]. Several studies showed that the prevalence of HF in patients with HCM could be from 50% to 67% [5, 6]. The high morbidity and mortality of older HCM patients indicate that it is important to investigate the mechanism by which HCM patients develop HF.

The development of bioinformatics tools to explore gene expression data provides the opportunity to understand the molecular mechanism of a variety of diseases [7, 8]. This study aimed to identify differentially expressed genes (DEGs) to evaluate the risk of HF in older patients with HCM based on bioinformatics analysis of gene expression data in patients with HCM and HF.

Materials and Methods

High-throughput sequencing data and differentially expressed genes analysis

The high-throughput sequencing datasets GSE89714 and GSE116250 were downloaded from Gene Expression Omnibus (GEO) database. GSE89714 dataset formed the study of “Differential gene expressions in the heart of hypertrophic cardiomyopathy patients”, and GSE116250 formed the study of “RNA-seq of heart failure in human left ventricles”. The DEGs were analyzed using linear models for microarray analysis R program with P < 0.05 and logFC > 1 as cut off. The heatmap and the hub genes were created using the NetworkAnalyst program (https://www.networkanalyst.ca/) [9].

Network establishment

The overlap of HCM and HF datasets’ DEGs were outlined with Venn diagram (https://bioinfogp.cnb.csic.es/tools/venny/index.html). HCM-overlap genes-HF network was established using Cytoscape 3.8.0 software, and protein-protein interactions (PPIs) of the DEGs were analyzed using STRING database (https://string-db.org/).

Gene functional and pathway enrichment analysis

Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) database (https://david.ncifcrf.gov/). NetworkAnalyst online tool was applied for Gene Set Enrichment Analysis (GSEA) analysis. Briefly, GO items for the overlap genes were collected and imported into GSEA software to identify significantly enriched GO items.

Statistical analysis

Data were analyzed by using GraphPad Prism 8.0 software. Two-tailed Student’s t-test was used to compare the data in two groups. P < 0.05 was deemed significant.

Availability of data and materials

All data are available from the corresponding author on request.

Results

Identification of DEGs in HCM and HF patients

LV tissues high-throughput sequencing data datasets GSE89714 and GSE116250 were downloaded from GEO database. The volcano plot and heatmap of DEGs showed a total of 655 DEGs with 437 upregulated and 218 downregulated in HCM patients (Figure 1), and a total of 859 DEGs with 490 upregulated and 369 downregulated in HF patients (Figure 2).

Bioinformatic analysis of DEGs in LV tissue of normal control (NC) and hypertrophic cardiomyopathy (HCM) patients. (A) The volcano plot of DEGs in LV tissue between NC group and HCM group. (B) Heatmaps of DEGs in LV tissue of NC group and HCM group. Red color indicated high expression while blue color indicated low expression.

Figure 1. Bioinformatic analysis of DEGs in LV tissue of normal control (NC) and hypertrophic cardiomyopathy (HCM) patients. (A) The volcano plot of DEGs in LV tissue between NC group and HCM group. (B) Heatmaps of DEGs in LV tissue of NC group and HCM group. Red color indicated high expression while blue color indicated low expression.

Bioinformatic analysis of DEGs in LV tissue of normal control (NC) and heart failure (HF) patients. (A) The volcano plot of DEGs in LV tissue between NC group and HF group. (B) Heatmaps of DEGs in LV tissue of NC group and HF group. Red color indicated high expression while blue color indicated low expression.

Figure 2. Bioinformatic analysis of DEGs in LV tissue of normal control (NC) and heart failure (HF) patients. (A) The volcano plot of DEGs in LV tissue between NC group and HF group. (B) Heatmaps of DEGs in LV tissue of NC group and HF group. Red color indicated high expression while blue color indicated low expression.

The protein-protein interaction analysis of the HCM and HF overlap genes

To further analyze the relationship between HCM and HF, the interaction of HCM and HF DEGs was shown in a Venn diagram. We found 124 genes closely related to HF progression (Figure 3A). In addition, we established an HCM-overlap gene-HF network (Figure 3B), highlighting that HCM may deteriorate to HF through these genes.

HCM and HF datasets overlap DEGs. (A) Venn diagram of HCM and HF DEGs. (B) The network of HCM, HF, and all the overlap DEGs. Yellow nodes represented the overlap DEGs.

Figure 3. HCM and HF datasets overlap DEGs. (A) Venn diagram of HCM and HF DEGs. (B) The network of HCM, HF, and all the overlap DEGs. Yellow nodes represented the overlap DEGs.

PPI analysis of the 124 overlap genes established a network: the top 10 hub genes were highlighted in red, and other genes connected with hub genes were shown as blue nodes (Figure 4A). Hub genes such as COL1A1, COL3A1, COL1A2, BGN, and COL5A1 are involved in the biological process of HF (Figure 4B).

The interaction network of HCM and HF overlap genes. (A) PPI network of the overlap genes. The hub genes were represented as red and yellow nodes. A deeper red color indicated more connections. (B) Bar plot of the number of hub gene links.

Figure 4. The interaction network of HCM and HF overlap genes. (A) PPI network of the overlap genes. The hub genes were represented as red and yellow nodes. A deeper red color indicated more connections. (B) Bar plot of the number of hub gene links.

Functional enrichment analysis of overlap genes

GO enrichment analysis showed that extracellular matrix (ECM) structural constituent and collagen binding (Figure 5A), ECM and extracellular region (Figure 5B), and collagen fibril organization and ECM organization (Figure 5C) were enriched and could play an important role in the biological process of HF.

GO enrichment analysis of HCM and HF overlap genes. (A) The plot of enriched molecular functions. (B) The plot of enriched cellular components. (C) The plot of enriched biological processes. The number of genes enriched in each GO term was shown as the circle size, and the p-value was shown as different colors.

Figure 5. GO enrichment analysis of HCM and HF overlap genes. (A) The plot of enriched molecular functions. (B) The plot of enriched cellular components. (C) The plot of enriched biological processes. The number of genes enriched in each GO term was shown as the circle size, and the p-value was shown as different colors.

In addition, eight KEGG pathways were enriched (Figure 6A) with genes mainly distributed in ECM-receptor interaction pathway (Figure 6B). Furthermore, GSEA analysis of the overlap genes showed that the DEGs were significantly correlated with ECM organization, ECM receptor interaction, Focal adhesion, TGF-beta signaling, PI3K-Akt signaling, and Platelet activation signaling and aggregation (Figure 7).

KEGG enrichment analysis of the HCM and HF overlap genes. (A) KEGG annotation of overlap genes. The number of genes enriched in each KEGG term was shown as the circle size, and the p-value was shown as different colors. (B) The overlap genes were mainly distributed in the ECM-receptor interaction pathway. Arrows represented activation effect, T-arrows represented inhibition effect and segments showed activation effect or inhibition effect. The red nodes were the intersection genes.

Figure 6. KEGG enrichment analysis of the HCM and HF overlap genes. (A) KEGG annotation of overlap genes. The number of genes enriched in each KEGG term was shown as the circle size, and the p-value was shown as different colors. (B) The overlap genes were mainly distributed in the ECM-receptor interaction pathway. Arrows represented activation effect, T-arrows represented inhibition effect and segments showed activation effect or inhibition effect. The red nodes were the intersection genes.

GSEA of the HCM and HF overlap genes. The enrichment plot of the ECM receptor interaction (A), Extracellular matrix organization (B), Focal adhesion (C), PI3K-Akt signaling pathway (D), TGF-beta signaling pathway (E), and Platelet activation signaling and aggregation (F). The green curve represented the enrichment profile, and the black vertical line represented the gene hits.

Figure 7. GSEA of the HCM and HF overlap genes. The enrichment plot of the ECM receptor interaction (A), Extracellular matrix organization (B), Focal adhesion (C), PI3K-Akt signaling pathway (D), TGF-beta signaling pathway (E), and Platelet activation signaling and aggregation (F). The green curve represented the enrichment profile, and the black vertical line represented the gene hits.

The changes of overlap hub gene expression in HCM and HF datasets

Finally, we selected ten hub genes COL1A1, COL3A1, COL1A2, BGN, COL5A1, LUM, TGFB2, FMOD, ASPN, and COL14A1 to confirm their expression patterns in HCM and HF datasets (Figure 8). Among them, COL5A1 and LUM were significantly increased in the HF dataset (P < 0.05, Figure 8E, 8F), but TGFB2, FMOD, ASPN, and COL14A1 were significantly decreased in the HF dataset (P < 0.05, Figure 8G8J).

The overlap genes expression pattern in HCM and HF datasets. The expression pattern of COL1A1 (A), COL3A1 (B), COL1A2 (C), BGN (D), COL5A1 (E), and LUM (F), TGFB2 (G), FMOD (H), ASPN (I), COL14A1 (J) in different datasets. Values were normalized to NC group and represented as Mean ± SD (n = 3 in HCM dataset, n = 13 in HF dataset). *P **P

Figure 8. The overlap genes expression pattern in HCM and HF datasets. The expression pattern of COL1A1 (A), COL3A1 (B), COL1A2 (C), BGN (D), COL5A1 (E), and LUM (F), TGFB2 (G), FMOD (H), ASPN (I), COL14A1 (J) in different datasets. Values were normalized to NC group and represented as Mean ± SD (n = 3 in HCM dataset, n = 13 in HF dataset). *P < 0.05, **P< 0.01 vs. HCM disease group.

Discussion

The rate of HF complications in HCM patients ranges from 2.4% to 20%, and increases significantly in older patients [10, 11]. However, the mechanistic link between HCM and HF remains unclear. Therefore, in this study we utilized gene expression datasets of HCM and HF to identify 124 overlap genes. Among them, top ten hub genes included COL1A1, COL3A1, COL1A2, BGN, COL5A1, LUM, TGFB2, FMOD, ASPN, and COL14A1. These results revealed that DEGs were mainly involved in myocardial fibrosis, and HCM may deteriorate to HF through these genes.

Myocardial fibrosis in pathological cardiac remodeling contributes to HCM and HF [12]. COL5A1 is an alpha chain for one of the low abundance fibrillar collagens, and regulates the function of connective tissues [13]. As a member of the leucine-rich proteoglycan family, LUM regulates the assembly of collagen fibers, and both COL1A1 and LUM are involved in heart fibrosis [14].

To gain insight into the biological function of the overlap genes, we performed GO enrichment analysis. The overlap genes were enriched in the following functional categories including, ECM structural constituent and Collagen binding, cellular components as ECM and Extracellular region, and molecular function as Collagen fibril organization and ECM organization. KEGG pathway analysis showed that the DEGs were enriched in ECM organization, ECM receptor interaction, Focal adhesion, TGF-beta signaling, PI3K-Akt signaling, and Platelet activation signaling and aggregation. ECM-receptor interaction and Focal adhesion have been shown to be important for HCM and HF progression [15]. A previous study showed that platelet aggregation was positively correlated with LV hypertrophy [16]. TGFβ signaling is known to play a role in hypertrophy [17]. Notably, a recent study identified six TGFβ related genes involved in both cardiac hypertrophy and HF based on single-cell RNA sequencing [18].

The limitations of this study should be pointed out. This study is in silico and our results should be confirmed by further in vivo studies. Moreover, we only screened the DEGs and pathways involved in HCM and HF. Further dissection of these DEGs and pathways will help elucidate the mechanism underlying the progression of HCM and HF.

In conclusion, we performed bioinformatics analysis to identify potential targets involved in HCM and HF, which could be utilized to evaluate the risk of HF in older patients with HCM.

Author Contributions

YT designed the study. HD, CY, and DX performed the study and analyzed the data. All authors wrote and approved the manuscript.

Conflicts of Interest

The authors have no conflicts of interest to declare.

Funding

This study received no funding.

References

  • 1. Maron BJ, Maron MS. Hypertrophic cardiomyopathy. Lancet. 2013; 381:242–55. https://doi.org/10.1016/S0140-6736(12)60397-3 [PubMed]
  • 2. Maron BJ, Rowin EJ, Udelson JE, Maron MS. Clinical Spectrum and Management of Heart Failure in Hypertrophic Cardiomyopathy. JACC Heart Fail. 2018; 6:353–63. https://doi.org/10.1016/j.jchf.2017.09.011 [PubMed]
  • 3. Olivotto I, Cecchi F, Gistri R, Lorenzoni R, Chiriatti G, Girolami F, Torricelli F, Camici PG. Relevance of coronary microvascular flow impairment to long-term remodeling and systolic dysfunction in hypertrophic cardiomyopathy. J Am Coll Cardiol. 2006; 47:1043–8. https://doi.org/10.1016/j.jacc.2005.10.050 [PubMed]
  • 4. Seferović PM, Polovina M, Bauersachs J, Arad M, Ben Gal T, Lund LH, Felix SB, Arbustini E, Caforio AL, Farmakis D, Filippatos GS, Gialafos E, Kanjuh V, et al. Heart failure in cardiomyopathies: a position paper from the Heart Failure Association of the European Society of Cardiology. Eur J Heart Fail. 2019; 21:553–76. https://doi.org/10.1002/ejhf.1461 [PubMed]
  • 5. Charron P, Elliott PM, Gimeno JR, Caforio AL, Kaski JP, Tavazzi L, Tendera M, Maupain C, Laroche C, Rubis P, Jurcut R, Calò L, Heliö TM, et al, and EORP Cardiomyopathy Registry Investigators. The Cardiomyopathy Registry of the EURObservational Research Programme of the European Society of Cardiology: baseline data and contemporary management of adult patients with cardiomyopathies. Eur Heart J. 2018; 39:1784–93. https://doi.org/10.1093/eurheartj/ehx819 [PubMed]
  • 6. Maron BJ, Rowin EJ, Casey SA, Link MS, Lesser JR, Chan RH, Garberich RF, Udelson JE, Maron MS. Hypertrophic Cardiomyopathy in Adulthood Associated With Low Cardiovascular Mortality With Contemporary Management Strategies. J Am Coll Cardiol. 2015; 65:1915–28. https://doi.org/10.1016/j.jacc.2015.02.061 [PubMed]
  • 7. Shen H, Ren S, Wang W, Zhang C, Hao H, Shen Q, Duan Y, Wang Z, Ge W. Profiles of immune status and related pathways in sepsis: evidence based on GEO and bioinformatics. Biocell 2020; 44:583–9. https://doi.org/10.32604/biocell.2020.011345
  • 8. Gao W, Ma Q, Tang C, Zhan Y, Duan Y, Ni H, Xu Y. Microenvironment and related genes predict outcomes of patients with cervical cancer: evidence from TCGA and bioinformatic analysis. Biocell 2020; 44:597–605. https://doi.org/10.32604/biocell.2020.011328
  • 9. Zhou G, Soufan O, Ewald J, Hancock RE, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. 2019; 47:W234–41. https://doi.org/10.1093/nar/gkz240 [PubMed]
  • 10. Biagini E, Coccolo F, Ferlito M, Perugini E, Rocchi G, Bacchi-Reggiani L, Lofiego C, Boriani G, Prandstraller D, Picchio FM, Branzi A, Rapezzi C. Dilated-hypokinetic evolution of hypertrophic cardiomyopathy: prevalence, incidence, risk factors, and prognostic implications in pediatric and adult patients. J Am Coll Cardiol. 2005; 46:1543–50. https://doi.org/10.1016/j.jacc.2005.04.062 [PubMed]
  • 11. Melacini P, Basso C, Angelini A, Calore C, Bobbo F, Tokajuk B, Bellini N, Smaniotto G, Zucchetto M, Iliceto S, Thiene G, Maron BJ. Clinicopathological profiles of progressive heart failure in hypertrophic cardiomyopathy. Eur Heart J. 2010; 31:2111–23. https://doi.org/10.1093/eurheartj/ehq136 [PubMed]
  • 12. Raman B, Ariga R, Spartera M, Sivalokanathan S, Chan K, Dass S, Petersen SE, Daniels MJ, Francis J, Smillie R, Lewandowski AJ, Ohuma EO, Rodgers C, et al. Progression of myocardial fibrosis in hypertrophic cardiomyopathy: mechanisms and clinical implications. Eur Heart J Cardiovasc Imaging. 2019; 20:157–67. https://doi.org/10.1093/ehjci/jey135 [PubMed]
  • 13. Roulet M, Ruggiero F, Karsenty G, LeGuellec D. A comprehensive study of the spatial and temporal expression of the col5a1 gene in mouse embryos: a clue for understanding collagen V function in developing connective tissues. Cell Tissue Res. 2007; 327:323–32. https://doi.org/10.1007/s00441-006-0294-1 [PubMed]
  • 14. Weber KT. Cardioreparation in hypertensive heart disease. Hypertension. 2001; 38:588–91. https://doi.org/10.1161/01.hyp.38.3.588 [PubMed]
  • 15. Perestrelo AR, Silva AC, Oliver-De La Cruz J, Martino F, Horváth V, Caluori G, Polanský O, Vinarský V, Azzato G, de Marco G, Žampachová V, Skládal P, Pagliari S, et al. Multiscale Analysis of Extracellular Matrix Remodeling in the Failing Heart. Circ Res. 2021; 128:24–38. https://doi.org/10.1161/CIRCRESAHA.120.317685 [PubMed]
  • 16. Cambronero F, Marín F, Roldán V, Hernández-Romero D, Valdés M, Lip GY. Biomarkers of pathophysiology in hypertrophic cardiomyopathy: implications for clinical management and prognosis. Eur Heart J. 2009; 30:139–51. https://doi.org/10.1093/eurheartj/ehn538 [PubMed]
  • 17. Bagnall RD, Tsoutsman T, Shephard RE, Ritchie W, Semsarian C. Global microRNA profiling of the mouse ventricles during development of severe hypertrophic cardiomyopathy and heart failure. PLoS One. 2012; 7:e44744. https://doi.org/10.1371/journal.pone.0044744 [PubMed]
  • 18. Huang K, Wu H, Xu X, Wu L, Li Q, Han L. Identification of TGF-β-related genes in cardiac hypertrophy and heart failure based on single cell RNA sequencing. Aging (Albany NY). 2023; 15:7187–218. https://doi.org/10.18632/aging.204901 [PubMed]