Identification of candidate biomarkers and therapeutic agents for heart failure by bioinformatics analysis

Heart failure (HF) is a heterogeneous clinical syndrome and affects millions of people all over the world. HF occurs when the cardiac overload and injury, which is a worldwide complaint. The aim of this study was to screen and verify hub genes involved in developmental HF as well as to explore active drug molecules. The expression profiling by high throughput sequencing of GSE141910 dataset was downloaded from the Gene Expression Omnibus (GEO) database, which contained 366 samples, including 200 heart failure samples and 166 non heart failure samples. The raw data was integrated to find differentially expressed genes (DEGs) and were further analyzed with bioinformatics analysis. Gene ontology (GO) and REACTOME enrichment analyses were performed via ToppGene; protein-protein interaction (PPI) networks of the DEGs was constructed based on data from the HiPPIE interactome database; modules analysis was performed; target gene - miRNA regulatory network and target gene - TF regulatory network were constructed and analyzed; hub genes were validated; molecular docking studies was performed. A total of 881 DEGs, including 442 up regulated genes and 439 down regulated genes were observed. Most of the DEGs were significantly enriched in biological adhesion, extracellular matrix, signaling receptor binding, secretion, intrinsic component of plasma membrane, signaling receptor activity, extracellular matrix organization and neutrophil degranulation. The top hub genes ESR1, PYHIN1, PPP2R2B, LCK, TP63, PCLAF, CFTR, TK1, ECT2 and FKBP5 were identified from the PPI network. Module analysis revealed that HF was associated with adaptive immune system and neutrophil degranulation. The target genes, miRNAs and TFs were identified from the target gene - miRNA regulatory network and target gene - TF regulatory network. Furthermore, receiver operating characteristic (ROC) curve analysis and RT-PCR analysis revealed that ESR1, PYHIN1, PPP2R2B, LCK, TP63, PCLAF, CFTR, TK1, ECT2 and FKBP5 might serve as prognostic, diagnostic biomarkers and therapeutic target for HF. The predicted targets of these active molecules were then confirmed. The current investigation identified a series of key genes and pathways that might be involved in the progression of HF, providing a new understanding of the underlying molecular mechanisms of HF.