Prognostic Model and Immune Infiltration of Ferroptosis Subcluster-Related Modular Genes in Gastric Cancer

Gastric cancer (GC) is one of the gastrointestinal tumors with the highest mortality rate. The number of GC patients is still high. As a way of iron-dependent programmed cell death, ferroptosis activates lipid peroxidation and accumulates large reactive oxygen species. The role of ferroptosis in GC prognosis was underrepresented. The objective was to investigate the role of ferroptosis-related genes (FRGs) in the prognosis and development of GC.

Datasets of GC patients were obtained from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) database that include clinical information and RNA seq data. Through nonnegative matrix factorization (NMF) clustering, we identified and unsupervised cluster analysis of the expression matrix of FRGs. And we constructed the co-expression network between genes and clinical characteristics by consensus weighted gene co-expression network analysis (WGCNA). The prognostic model was constructed by univariate and multivariate regression analysis. The potential mechanisms of development and prognosis in GC were explored by Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, gene ontology (GO), tumor immune microenvironment (TIME), and tumor mutation burden (TMB).

Two molecular subclusters with different expression patterns of FRGs were identified, which have significantly different survival states. Ferroptosis subcluster-related modular genes were identified by WGCNA. Based on 8 ferroptosis subcluster-related modular genes (collagen triple helix repeat containing 1 (CTHRC1), podoplanin (PDPN), procollagen-lysine,2-oxoglutarate 5-dioxygenase 2 (PLOD2), glutamine-fructose-6-phosphate transaminase 2 (GFPT2), ATP-binding cassette subfamily A member 1 (ABCA1), G protein-coupled receptor 176 (GPR176), serpin family E member 1 (SERPINE1), dual specificity phosphatase 1 (DUSP1)) and clinicopathological features, a nomogram was constructed and validated for their predictive efficiency on GC prognosis. Through receiver operating characteristic (ROC) analysis, the results showed that the area under the curve (AUC) of 1-, 3-, and 5-year survival were 0.721, 0.747, and 0.803, respectively, indicating that the risk-scoring model we constructed had good prognosis efficacy in GC. The degree of immune infiltration in high-risk group was largely higher than low-risk group. It indicated that the immune cells have a good response in high-risk group of GC. The TMB of high-risk group was higher, which could generate more mutations and was more conducive to the body's resistance to the development of cancer.

The risk-scoring model based on 8 ferroptosis subcluster-related modular genes has shown outstanding advantages in predicting patient prognosis. The interaction of ferroptosis in GC development may provide new insights into exploring molecular mechanisms and targeted therapies for GC patients.