Application of Weighted Gene Coexpression Network Analysis to Identify Key Modules and Hub Genes in Systemic Juvenile Idiopathic Arthritis

Abstract

Systemic juvenile idiopathic arthritis (sJIA) is a severe autoinflammatory disorder with a still not clearly defined molecular mechanism. To better understand the disease, we used scattered datasets from public domains and performed a weighted gene coexpression network analysis (WGCNA) to identify key modules and hub genes underlying sJIA pathogenesis. Two gene expression datasets, GSE7753 and GSE13501, were used to construct the WGCNA. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were applied to the genes and hub genes in the sJIA modules. Cytoscape was used to screen and visualize the hub genes. We further compared the hub genes with the genome-wide association study (GWAS) genes and used a consensus WGCNA to verify that our conclusions were conservative and reproducible across multiple independent datasets. A total of 5,414 genes were obtained for WGCNA, from which highly correlated genes were divided into 17 modules. The red module demonstrated the highest correlation with the sJIA module ( $r=0.8$ , $p=3e^{-29}$ ), whereas the green-yellow module was found to be closely related to the non-sJIA module ( $r=0.62$ , $p=1e^{-14}$ ). Functional enrichment analysis demonstrated that the red module was mostly enriched in the activation of immune responses, infection, nucleosomes, and erythrocytes, and the green-yellow module was mostly enriched in immune responses and inflammation. Additionally, the hub genes in the red module were highly enriched in erythrocyte differentiation, including ALAS2, AHSP, TRIM10, TRIM58, and KLF1. The hub genes from the green-yellow module were mainly associated with immune responses, as exemplified by the genes KLRB1, KLRF1, CD160, and KIRs. We identified sJIA-related modules and several hub genes that might be associated with the development of sJIA. Particularly, the modules may help understand the mechanisms of sJIA, and the hub genes may become biomarkers and therapeutic targets of sJIA in the future.