Identification and interaction analysis of molecular markers in myocardial infarction by bioinformatics and next-generation sequencing data analysis

Abstract

Background Cardiovascular diseases are prevalent worldwide with any age, and it is characterized by sudden blockage of blood flow to heart and permanent damage to the heart muscle, whose cause and underlying molecular mechanisms are not fully understood. This investigation aimed to explore and identify essential genes and signaling pathways that contribute to the progression of MI. Methods The aim of this investigation was to use bioinformatics and next-generation sequencing (NGS) data analysis to identify differentially expressed genes (DEGs) with diagnostic and therapeutic potential in MI. NGS dataset (GSE132143) was downloaded from the Gene Expression Omnibus (GEO) database. DEGs between MI and normal control samples were identified using the DESeq2 R bioconductor tool. The gene ontology (GO) and REACTOME pathway enrichment analyses of the DEGs were performed using g:Profiler. Next, four kinds of algorithms in the protein–protein interaction (PPI) were performed to identify potential novel biomarkers. Next, miRNA-hub gene regulatory network analysis and TF-hub gene regulatory network were constructed by miRNet and NetworkAnalyst database, and Cytoscape software. Finally, the diagnostic effectiveness of hub genes was predicted by receiver operator characteristic curve (ROC) analysis and AUC more than 0.800 was considered as having the capability to diagnose MI with excellent specificity and sensitivity. Results A total of 958 DEGs were identified, consisting of 480 up-regulated genes and 478 down-regulated genes. The enriched GO terms and pathways of the DEGs include immune system, neuronal system, response to stimulus and multicellular organismal process. Ten hub genes (namely cftr, cdk1, rps13, rps15a, rps27, notch1, mrpl12, nos2, ccdc85b and atn1) were obtained via protein–protein interaction analysis results. MiRNA-hub gene regulatory network and TF-hub gene regulatory network showed that hsa-mir-409-3p, hsa-mir-3200-3p, creb1 and tp63 might play an important role in the MI. Conclusions Analysis of next-generation sequencing dataset combined with global network information and validation presents a successful approach to uncover the risk hub genes and prognostic markers of MI. Our investigation identified four risk- and prognostic-related gene signatures, including cftr, cdk1, rps13, rps15a, rps27, notch1, mrpl12, nos2, ccdc85b and atn1. This gene sets contribute a new perspective to improve the diagnostic, prognostic, and therapeutic outcomes of MI.