Abstract
BACKGROUND: Cardiovascular complications are a major cause of death and disability in patients with diabetes mellitus, but how such complications arise is unclear. METHODS: Weighted gene correlation network analysis (WGCNA) was performed on gene expression profiles from healthy controls, individuals with diabetes mellitus, and individuals with diabetes mellitus-associated coronary artery disease (DMCAD). Phenotypically related module genes were analyzed for enrichment in Gene Ontology (GO) terms and Kyoto Gene and Genome Encyclopedia (KEGG) pathways. Predicted biological functions were validated using gene set enrichment analysis (GSEA) and ClueGo analysis. Based on the TRRUST v2 database and hypergeometric tests, a global network was built to identify transcription factors (TFs) and downstream target genes potentially involved in DMCAD. RESULTS: WGCNA identified three modules associated with progression from diabetes mellitus to DMCAD. The module genes were significantly involved in biological processes related to interferon and viral infection, while GSEA of DMCAD samples suggested involvement in viral myocarditis, chemokine signaling and phagosomes. RUNX1 was identified as a potential TF regulating these module genes. Analysis of the global regulatory network of TFs and their targets suggested that CCL3 may be a key regulator in DMCAD. CONCLUSION: We found bioinformatic evidence that CCL3 may be a key regulator and RUNX1 a key TF in DMCAD.