Integrated analysis of gene expression changes associated with coronary artery disease

This study investigated the pathways and genes involved in coronary artery disease (CAD) and the associated mechanisms.

Two array data sets of GSE19339 and GSE56885 were downloaded. The limma package was used to analyze the differentially expressed genes (DEGs) in normal and CAD specimens. Examination of DEGs through Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment and Gene Ontology annotation was achieved by Database for Annotation, Visualization and Integrated Discovery (DAVID). The Cytoscape software facilitated the establishment of the protein-protein interaction (PPI) network and Molecular Complex Detection (MCODE) was performed for the significant modules.

We identified 413 DEGs (291 up-regulated and 122 down-regulated). Approximately 256 biological processes, only 1 cellular component, and 21 molecular functions were identified by GO analysis and 10 pathways were enriched by KEGG. Moreover, 264 protein pairs and 64 nodes were visualized by the PPI network. After the MCODE analysis, the top 4 high degree genes, including interleukin 1 beta (IL1B, degree = 29), intercellular adhesion molecule 1 (ICAM1, degree = 25), Jun proto-oncogene (JUN, degree = 23) and C-C motif chemokine ligand 2 (CCL2, degree = 20) had been identified to validate in RT-PCR and Cox proportional hazards regression between CAD and normals.

The relative expression of IL1B, ICAM1 and CCL2 was higher in CAD than in normal controls (P < 0.05-0.001), but only IL1B and CCL2 genes were confirmed after testing the gene expression in blood and/or analyzing in Cox proportional hazards regression (P < 0.05-0.001), and the proper mechanism may involve in the AGE-RAGE signaling pathway, fluid shear stress, the tumor necrosis factor (TNF) and cytokine-cytokine receptor interaction.