Identification of the susceptible genes and mechanism underlying the comorbid presence of coronary artery disease and rheumatoid arthritis: a network modularization analysis

Abstract

Abstract Objective The comorbidities of coronary artery disease (CAD) and rheumatoid arthritis (RA) are mutual risk factors, which lead to higher mortality, but the biological mechanisms connecting the two remain unclear. Here, we aimed to identify the risk genes for the comorbid presence of these two complex diseases using a network modularization approach, to offer insights into clinical therapy and drug development for these diseases. Method The expression profile data of patients CAD with and without RA were obtained from the GEO database (GSE110008). Based on the differentially expressed genes (DEGs), weighted gene co-expression network analysis (WGCNA) was used to construct a gene network, detect co-expression modules, and explore their relation to clinical traits. The Zsummary index, gene significance (GS), and module membership (MM) were utilized to screen the important differentiated modules and hub genes. The GO and KEGG pathway enrichment analysis were applied to analyze potential mechanisms. Result Based on the 278 DEGs obtained, 41 modules were identified, of which 17 and 24 modules were positively and negatively correlated with the comorbid occurrence of CAD and RA (CAD&RA), respectively. Thirteen modules with Zsummary < 2 were found to be the underlying modules, which may be related to CAD&RA. With GS ≥ 0.5 and MM ≥ 0.8, 49 hub genes were identified, such as ADO, ABCA11P, POT1, ZNF141, GPATCH8, ATF6 and MIA3, etc. The area under the curve values of the representative seven hub genes under the three models (LR, KNN, SVM) were greater than 0.88. Enrichment analysis revealed that the biological functions of the targeted modules were mainly involved in cAMP-dependent protein kinase activity, demethylase activity, regulation of calcium ion import, positive regulation of tyrosine, phosphorylation of STAT protein, and tissue migration, etc. Conclusion Thirteen characteristic modules and 49 susceptibility hub genes were identified, and their corresponding molecular functions may reflect the underlying mechanism of CAD&RA, hence providing insights into the development of clinical therapies against these diseases.