Rong Yu

Background and Aims: This study aimed to develop a prognostic model based on DNA methylation-driven genes for patients with early-stage gastric cancer and to examine immune infiltration and function across varying risk levels. Methods: We analyzed data from stage I/II gastric cancer patients in The Cancer Genome Atlas which included clinical details, mRNA expression profiles, and level 3 DNA methylation array data. Using the empirical Bayes method of the limma package, we identified differentially expressed genes (DEGs), and the MethylMix package facilitated the identification of DNA methylation-driven genes (DMGs). Univariate Cox regression and LASSO (least absolute shrinkage and selector operation) analyses were utilized to pinpoint critical genes. A risk score prediction model was formulated using two genes that demonstrated the most significant hazard ratios (HRs). Model performance was evaluated within the initial cohort and verified in the GSE84437 cohort; a nomogram was also constructed based on these genes. We further examined 50 methylation sites associated with three CpG islands in C1orf35 and 14 methylation sites linked to one CpG island in FAAH. The CIBERSORT package was employed to identify immune cell clusters in the prediction model. Results: A total of 176 DNA methylation-driven genes were refined down to a four-gene signature (ZC3H12A was hypermethylated; GATA3, C1orf35, and FAAH were hypomethylated), which exhibited a significant correlation with overall survival (OS), as evidenced by p -values below 0.05 following univariate Cox regression and LASSO analysis. Specifically, for the risk score prediction model, C1orf35, which had the highest hazard ratio (HR = 2.035, p = 0.028), and FAAH, with the lowest hazard ratio (HR = 0.656, p = 0.012), were selected. The Kaplan–Meier analysis demonstrated distinct survival outcomes between the high-risk and low-risk score groups. The model’s predictive accuracy was confirmed with an area under the curve (AUC) of 0.611 for 3-year survival and 0.564 for 5-year survival. Notably, the hypomethylation of the three CpG islands in C1orf35 and the single CpG island in FAAH was significantly different in stage I/II gastric cancer patients compared to normal tissues. Additionally, the high-risk score group showed a notable association with resting CD4 memory T cells. Conclusion: Promoter hypomethylation of C1orf35 and FAAH in early-stage gastric cancer underscores their potential as biomarkers for accurate diagnosis and treatment. The developed predictive model employing genes affected by DNA methylation serves as a crucial independent prognostic factor in early-stage gastric cancer.

According to the World Health Organization (WHO), gastric cancer (GC) is the fourth leading cause of tumor-related mortality globally and one of the most prevalent malignant tumors. To better understand the role of tumor-infiltrating B cells (TIBs) in GC, this work used single-cell RNA sequencing (scRNA-Seq) and bulk RNA sequencing (bulk RNA-Seq) data to identify candidate hub genes. Both scRNA-Seq and bulk RNA-Seq data for stomach adenocarcinoma (STAD) were obtained from the GEO and TCGA databases, respectively. Using scRNA-seq data, the FindNeighbors and FindClusters tools were used to group the cells into distinct groups. Immune cell clusters were sought in the massive RNA-seq expression matrix using the single-sample gene set enrichment analysis (ssGSEA). The expression profiles were used in Weighted Gene Coexpression Network Analysis (WGCNA) to build TCGA’s gene coexpression networks. Next, univariate Cox regression, LASSO regression, and Kaplan–Meier analyses were used to identify hub genes in scRNA-seq data from sequential B-cell analyses. Finally, we examined the correlation between the hub genes and TIBs utilizing the TISIDB database. We confirmed the immune-related markers in clinical validation samples using reverse transcriptase polymerase chain reaction (RT-PCR) and immunohistochemistry (IHC). 15 cell clusters were classified in the scRNA-seq database. According to the WGCNA findings, the green module is most associated with cancer and B cells. The intersection of 12 genes in two separate datasets (scRNA and bulk) was attained for further analysis. However, survival studies revealed that increased C-X-C motif chemokine receptor 4 (CXCR4) expression was linked to worse overall survival. CXCR4 expression is correlated with active, immature, and memory B cells in STAD were identified. Finally, RT-PCR and IHC assays verified that in GC, CXCR4 is overexpressed, and its expression level correlates with TIBs. We used scRNA-Seq and bulk RNA-Seq to study STAD’s cellular composition. We found that CXCR4 is highly expressed by TIBs in GC, suggesting that it may serve as a hub gene for these cells and a starting point for future research into the molecular mechanisms by which these immune cells gain access to tumors and potentially identify therapeutic targets.