Integrated bioinformatics analysis identifies microRNA-376a-3p as a new microRNA biomarker in patient with coronary artery disease

Introduction: Coronary artery disease (CAD) is a major global health problem with high incidence and mortality. Despite many advances in treatment, the prognosis of patients with CAD still remains poor. Therefore, this study aimed to identify potential biomarkers and targets associated with the progression of CAD.

Methods: Two gene expression profile datasets (GSE20681 and GSE12288), and two microRNA (miRNA) expression profile datasets (GSE59421 and GSE105449) were downloaded from the Gene Expression Omnibus (GEO) database; R language was used to screen out the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs), respectively. In addition, five online bioinformatics tools (miRWalk et al.) were used to identify the target genes of DEMs, and miRNA-gene network was constructed using Cytoscape software. Moreover, CCK-8, flow cytometry assays were used to detect the cell proliferation and apoptosis in human umbilical vein endothelial cells (HUVECs). Meanwhile, the dual luciferase reporter system assay was used to explore the interaction of miR-376a-3p and NRIP1 in HUVECs.

Results: In the present study, 150 common DEGs and 5 common DEMs were screened using a Venn diagram in R language. First, a total of 6812 target genes were identified from the overlapping DEMs. Second, 26 overlapping dysregulated genes from 150 overlapping DEGs and 6812 miRNA target genes were identified. Meanwhile, 43 miRNA-gene regulatory pairs were obtained between the 5 common DEMs and 26 dysregulated genes. Downregulation of miR-376a-3p significantly inhibited the proliferation of HUVECs via inducing apoptosis. Moreover, overexpression of miR-376a-3p markedly inhibited the growth of HUVECs via downregulating NRIP1.

Conclusion: In this study, miR-376a-3p-NRIP1 pair might involve in the progression of CAD, implying that miR-376a-3p may be a therapeutic target for the treatment of CAD.