Analysis and identification of ferroptosis-related genes in ulcerative colitis

Previous studies have shown that ferroptosis is associated with the pathogenesis of ulcerative colitis (UC). Therefore, this study aimed to identify key ferroptosis-related genes (FRGs) associated with the diagnosis of UC.

UC-related expression datasets were downloaded from the Gene Expression Omnibus (GEO) database. First, Weighted Gene Co-expression Network Analysis (WGCNA) was used to identify UC-related genes (UCRGs). Differentially expressed genes (DEGs) between normal and UC groups were screened in GSE87466, and DEGs were subjected to an intersection analysis with FRGs and UCRGs to obtain ferroptosis-related DEGs (FR DEGs). Then a protein-protein interaction (PPI) network was constructed for FR DEGs. The hub genes were extracted based on the degree, Maximum Neighborhood Component (MNC), closeness, and Maximal Clique Centrality (MCC). Biomarkers with diagnostic values were screened by support vector machine (SVM) and the least absolute shrinkage and selection operator (LASSO) algorithms. Next, the infiltration of immune cells was compared between UC and normal groups, and the correlation between different immune cells and diagnostic genes was analyzed. The biological functions, classical pathways, and intermolecular interaction networks of diagnostic genes were characterized utilizing ingenuity pathway analysis (IPA). Finally, a TF-mRNA network was constructed and potential small-molecule compounds were screened.

Thirty-six FR DEGs were obtained, and these were enriched in biological processes such as positive regulation of cytokine production, cytokine-mediated signalling pathway, long-chain fatty acid-CoA ligase activity, etc. Among 18 hub genes, five genes (ALOX5, TIMP1, TNFAIP3, SOCS1, DUOX2) were captured with diagnostic values for UC, and they displayed significant differences between UC and normal groups. Sixteen immune cell infiltrates were significantly different between UC and normal groups, such as activated dendritic cells and resting dendritic cells. TNFAIP3 and ALOX5 were positively correlated with neutrophils, and TIMP1, SOCS1, ALOX5, and DUOX2 were negatively correlated with M2 macrophages. IPA showed that diagnostic genes were related to 43 function modules and activated 17 pathways. The constructed TF-mRNA regulatory network comprised three diagnostic genes and 17 differentially expressed TFs. Potential small-molecule compounds including helveticoside and cymarin were identified.

Our findings yielded several promising FRGs for UC, providing a scientific reference for further studies on the pathogenesis of UC.