Development and Validation of a Robust Immune-Related Prognostic Signature for Gastric Cancer

Development and verification of a manganese metabolism- and immune-related genes signature for prediction of prognosis and immune landscape in gastric cancer

Background

Gastric cancer (GC) poses a global health challenge due to its widespread prevalence and unfavorable prognosis. Although immunotherapy has shown promise in clinical settings, its efficacy remains limited to a minority of GC patients. Manganese, recognized for its role in the body's anti-tumor immune response, has the potential to enhance the effectiveness of tumor treatment when combined with immune checkpoint inhibitors.

Methods

Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases was utilized to obtain transcriptome information and clinical data for GC. Unsupervised clustering was employed to stratify samples into distinct subtypes. Manganese metabolism- and immune-related genes (MIRGs) were identified in GC by univariate Cox regression and least absolute shrinkage and selection operator (LASSO) regression analysis. We conducted gene set variation analysis, and assessed the immune landscape, drug sensitivity, immunotherapy efficacy, and somatic mutations. The underlying role of NPR3 in GC was further analyzed in the single-cell RNA sequencing data and cellular experiments.

Results

GC patients were classified into four subtypes characterized by significantly different prognoses and tumor microenvironments. Thirteen genes were identified and established as MIRGs, demonstrating exceptional predictive effectiveness in GC patients. Distinct enrichment patterns of molecular functions and pathways were observed among various risk subgroups. Immune infiltration analysis revealed a significantly greater abundance of macrophages and monocytes in the high-risk group. Drug sensitivity analysis identified effective drugs for patients, while patients in the low-risk group could potentially benefit from immunotherapy. NPR3 expression was significantly downregulated in GC tissues. Single-cell RNA sequencing analysis indicated that the expression of NPR3 was distributed in endothelial cells. Cellular experiments demonstrated that NPR3 facilitated the proliferation of GC cells.

Conclusion

This is the first study to utilize manganese metabolism- and immune-related genes to identify the prognostic MIRGs for GC. The MIRGs not only reliably predicted the clinical outcome of GC patients but also hold the potential to guide future immunotherapy interventions for these patients.

DNA damage repair molecular subtype derived immune signature applicable for the prognosis and immunotherapy response prediction in colon cancer

Background

The DNA damage repair (DDR) pathway is one of the pathways of tumor pathogenesis, but its relationship with the immunophenotype has not been clarified in colon cancer (CC).

Methods

We identified the differentially expressed immune-related genes (DEIRGs) between two DDR molecular subtypes, namely, C1 and C2, and used univariate Cox analysis and least absolute shrinkage and selection operator (LASSO) penalized Cox regression analysis to construct the risk score in the training cohort [n=1,009, a combination of The Cancer Genome Atlas (TCGA) and GSE39582]. Regarding the median risk score as the unified cutoff to classify the patients into high- and low-risk groups. Two independent cohorts (GSE17538, n=232; GSE38832, n=122) were used for external validation of the prognostic value of the risk score. The IMvigor210 cohort (n=348) was used to test the predictive value of the risk score for immunotherapy response. Gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were performed to discover the underlying mechanism. Immune cell infiltration was quantified by the single sample gene set enrichment analysis (ssGSEA) algorithm.

Results

The high-risk group showed significantly reduced overall survival (OS), disease-specific survival (DSS), disease-free survival (DFS), progression-free survival (PFS), and relapse-free survival (RFS) compared to the low-risk group, and the two groups differed significantly in lymphatic invasion, American Joint Committee on Cancer (AJCC) TNM stage, preoperative carcinoembryonic antigen (CEA) level, etc. The enrichment levels of pathways related to colorectal cancer, epithelial-mesenchymal transition (EMT), angiogenesis, hypoxia, P53, TGF-β, KRAS signaling, etc., were upregulated in the high-risk group, but DDR-related pathways were defective in the high-risk group. The immunophenotypes of the high-risk group tended to be desert and excluded, and the risk score of patients who responded to immunotherapy was significantly lower than that of patients who did not respond to immunotherapy. The higher the infiltration levels of gamma delta T cells (γδ T cells), immature dendritic cells, and T follicular helper (Tfh) cells, the more significant adverse impact on the prognosis of CC patients was exhibited and an obviously positive correlation with the risk score was showed.

Conclusions

An immune gene risk score associated with the DDR molecular subtype was built and verified herein; that is applicable to the prognosis and immunotherapy response prediction in CC.

Integrated bioinformatics analysis for conducting a prognostic model and identifying immunotherapeutic targets in gastric cancer

BACKGROUND

Gastric cancer is the third leading cause of death from cancer worldwide and has a poor prognosis. Practical risk scores and prognostic models for gastric cancer are lacking. While immunotherapy has succeeded in some cancers, few gastric cancer patients benefit from immunotherapy. Immune genes and the tumor microenvironment (TME) are essential for cancer progression and immunotherapy response. However, the roles of immune genes and the tumor microenvironment in immunotherapy remain unclear. The study aimed to construct a prognostic prediction model and identify immunotherapeutic targets for gastric cancer (GC) patients by exploring immune genes and the tumor microenvironment.

RESULTS

An immune-related risk score (IRRS) model, including APOH, RNASE2, F2R, DEFB126, CXCL6, and CXCL3 genes, was constructed for risk stratification. Patients in the low-risk group, which was characterized by elevated tumor mutation burden (TMB) have higher survival rate. The risk level was remarkably correlated with tumor-infiltrating immune cells (TIICs), the immune checkpoint molecule expression, and immunophenoscore (IPS). CXCL3 and CXCL6 were significantly upregulated in gastric cancer tissues compared with normal tissues using the UALCAN database and RT-qPCR. The nomogram showed good calibration and moderate discrimination in predicting overall survival (OS) at 1-, 3-, and 5- year for gastric cancer patients using risk-level and clinical characteristics.

CONCLUSION

Our findings provided a risk stratification and prognosis prediction tool for gastric cancer patients and further the research into immunotherapy in gastric cancer.

HIF1A overexpression predicts the high lymph node metastasis risk and indicates a poor prognosis in papillary thyroid cancer

Objective

To investigate the value of Hypoxia-inducible factor 1 A (HIF1A) in predicting lymph node metastasis (LNM) stage and clinical outcomes of papillary thyroid cancer (PTC) patients.

Materials and methods

The HIF1A gene expression analysis in PTC was performed by bioinformatics approaches followed by evaluating its protein level using immunohistochemistry analysis. The role of HIF1A in predicting the LNM stage was evaluated by logistic regression analysis, nomogram construction, and receiver operating characteristic (ROC) analysis. We performed survival analyses to determine its prognostic value. Enrichment analysis was conducted, and immune cell infiltration and stromal content were evaluated to examine the underlying mechanism of HIF1A in PTC.

Results

HIF1A transcription and protein levels were significantly high in PTC tissue (P < 0.05). Its overexpression predicted high LNM risk and unfavorable prognosis for PTC patients (P < 0.05). Cox regression analysis revealed HIF1A as an independent prognostic biomarker for the disease-free interval (DFI) (P < 0.01). In addition, HIF1A was positively related to tumor-suppressive immunity but was negatively correlated with anti-tumor immunity. HIF1A upregulation was also associated with increased stromal content.

Conclusions

HIF1A overexpression is an independent predictor for worse DFI in PTC. The HIF1A expression may affect the prognosis of PTC patients through immune- and stroma-related pathways. Our study provides new insight into the role of HIF1A in PTC biology and clinical management.

An immune-related prognostic signature associated with immune landscape and therapeutic responses in gastric cancer

Immune-related genes (IRGs) have attracted attention in recent years as therapeutic targets in various tumors. However, the role of IRGs in gastric cancer (GC) has not been clearly elucidated. This study presents a comprehensive analysis exploring the clinical, molecular, immune, and drug response features characterizing the IRGs in GC. Data were acquired from the TCGA and GEO databases. The Cox regression analyses were performed to develop a prognostic risk signature. The genetic variants, immune infiltration, and drug responses associated with the risk signature were explored using bioinformatics methods. Lastly, the expression of the IRS was verified by qRT-PCR in cell lines. In this manner, an immune-related signature (IRS) was established based on 8 IRGs. According to the IRS, patients were divided into the low-risk group (LRG) and high-risk group (HRG). Compared with the HRG, the LRG was characterized by a better prognosis, high genomic instability, more CD8+ T cell infiltration, greater sensitivity to chemotherapeutic drugs, and greater likelihood of benefiting from the immunotherapy. Moreover, the expression result showed good consistency between the qRT-PCR and TCGA cohort. Our findings provide insights into the specific clinical and immune features underlying the IRS, which may be important for patient treatment.

A nine–consensus–prognostic –gene–based prognostic signature, recognizing the dichotomized subgroups of gastric cancer patients with different clinical outcomes and therapeutic strategies

Background: The increasing prevalence and mortality of gastric cancer (GC) has promoted the urgent need for prognostic signatures to predict the long-term risk and search for therapeutic biomarkers.

Methods and materials: A total of 921 GC patients from three GEO cohorts were enrolled in the current study. The GSE15459 and GSE62254 cohorts were used to select the top prognostic gene via the evaluation of the area under the receiver operating characteristic (ROC) curve (AUC) values. The GSE84437 cohort was used as the external validation cohort. Least absolute shrinkage and selector operation (LASSO) regression analysis was applied to reduce the feature dimension and construct the prognostic signature. Furthermore, a nomogram was constructed by integrating the independent prognostic analysis and validated by calibration plot, decision curve analysis and clinical impact curve. The molecular features and response to chemo-/immunotherapy among risk subgroups were evaluated by the “MOVICS” and “ESTAMATE” R packages and the SubMap algorithm. Lauren classification and ACRG molecular subtype were obtained to compare with the risk model.

Results: Forty-four prognosis-associated genes were identified with a preset cutoff AUC value of 0.65 in both the GSE62254 and GSE15459 cohorts. With the 10-fold cross validation analysis of LASSO, nine genes were selected to construct the nine-consensus-prognostic-gene signature. The signature showed good prognostic value in the GSE62254 (p < 0.001, HR: 3.81, 95% CI: 2.44–5.956) and GSE15459 (p < 0.001, HR: 2.65, 95% CI: 1.892–3.709) cohorts and the external validation GSE84437 cohort (p < 0.001, HR: 2.06, 95% CI: 1.554–2.735). The nomogram constructed based on two independent predictive factors, tumor stage and the signature, predicted events tightly consistent with the actual (Hosmer–Lemeshow p value: 1-year, 0.624; 3-years, 0.795; 5-years, 0.824). For the molecular features, we observed the activation of apical junction, epithelial mesenchymal transition, and immune pathways in the high-risk group, while in the low-risk group, cell cycle associated G2M, E2F and MYC target pathways were activated. Based on the results we obtained, we indicated that gastric patients in the low-risk group are more suitable for 5-fluorouracil therapy, while high-risk group patients are more suitable for anti-CTLA4 immunotherapy, these results need more support in the further studies. After compare with proposed molecular subtypes, we realized that the nine-consensus prognostic gene signature is a powerful addition to identify the gastric patients with poor prognosis.

Conclusion: In summary, we constructed a robust nine-consensus-prognostic-gene signature for the prediction of GC prognosis, which can also predict the personalized treatment of GC patients.

Comprehensive analysis of metabolic pathway activity subtypes derived prognostic signature in hepatocellular carcinoma

OBJECTIVE

Metabolic reprogramming is one of the hallmarks of cancer, but metabolic pathway activity-related subtypes of hepatocellular carcinoma (HCC) have not been identified.

METHODS

Based on the quantification results of 41 metabolic pathway activities by gene set variation analysis, the training cohort (n = 609, merged by TCGA and GSE14520) was clustered into three subtypes (C1, C2, and C3) with the nonnegative matrix factorization method. Totally 1371 differentially expressed genes among C1, C2, and C3 were identified, and an 8-gene risk score was established by univariable Cox regression analysis, least absolute shrinkage and selection operator method, and multivariable Cox regression analysis.

RESULTS

C1 had the strongest metabolic activity, good prognosis, the highest CTNNB1 mutation rate, with massive infiltration of eosinophils and natural killer cells. C2 had the weakest metabolic activity, poor prognosis, was younger, was inclined to vascular invasion and advanced stage, had the highest TP53 mutation rate, exhibited a higher expression level of immune checkpoints, accompanied by massive infiltration of regulatory T cells. C3 had moderate metabolic activity and prognosis, the highest LRP1B mutation rate, and a higher infiltration level of neutrophils and macrophages. Internal cohorts (TCGA, n = 370; GSE14520, n = 239), external cohorts (ICGC, n = 231; GSE116174, n = 64), and clinical subgroup validation showed that the risk score was applicable for patients with diverse clinical features and was effective in predicting the prognosis and malignant progression of patients with HCC. Compared with the low-risk group, the high-risk group had a poor prognosis, enhanced cancer stem cell characteristics, activated DNA damage repair, weakened metabolic activity, cytolytic activity, and interferon response.

CONCLUSION

We identified HCC subtypes from the perspective of metabolism-related pathway activity and proposed a robust prognostic signature for HCC.

A Five-Gene Signature Associated With DNA Damage Repair Molecular Subtype Predict Overall Survival for Hepatocellular Carcinoma

Background: DNA damage repair (DDR) is an important mechanism for the occurrence and development of hepatocellular carcinoma (HCC), but its impact on prognosis has not been fully understood.

Materials and methods: A total of 904 HCC patients were included in our study, TCGA (n = 370) and GSE14520 (n = 239) were merged into a large-sample training cohort (n = 609). The training cohort was clustered into C1 and C2 based on prognostic DDR-related genes, the differentially expressed genes (DEGs) between C1 and C2 were identified by the Wilcoxon signed-rank test referred to criteria (|log2FC|≥1 and FDR&lt; 0.05). The univariate Cox analysis was used to screen the prognostic-related DEGs, and Lasso penalized Cox regression analysis was used to construct the risk score. The patients were clarified into high- and low-risk groups based on the median risk score. ICGC (n = 231) and GSE116174 (n = 64) cohorts were used for external validation of the risk score’s prognostic value.

Results: The Kaplan–Meier survival analysis showed that the high-risk group had a significantly reduced overall survival (OS) compared to the low-risk group in the three independent cohorts, and the time-dependent ROC curve showed that the five-gene (STMN1, PON1, PLOD2, MARCKSL1, and SPP1) risk score with a high accuracy in predicting OS. The patients with AFP &gt;300 ng/ml, tumor poor differentiation (grade 3–4), micro and macro vascular tumor invasion, advanced stage (AJCC III-IV, BCLC stage B-C, and CLIP score &gt;2) exhibited a higher risk score. Subgroup survival analysis found that the risk score was applicable to patients with different clinical characteristics. GO and KEGG functional enrichment analysis revealed that cell cycle, p53 signaling, TNF signaling-related pathways were upregulated in the high-risk group. The higher infiltration level of activated CD4 T cell, CD56 bright natural killer cell, plasmacytoid dendritic cell, and type 2 T helper cells were found to lead an unfavorable impact on the OS of HCC patients, and these four kinds of immune cells exhibited a higher infiltration level in the high-risk group.

Conclusion: The five-gene risk score proposed in the research may provide new insights into the individualized evaluation of HCC prognosis.

A Ferroptosis and Pyroptosis Molecular Subtype-Related Signature Applicable for Prognosis and Immune Microenvironment Estimation in Hepatocellular Carcinoma

Background: Due to the heterogeneity of tumors and the complexity of the immune microenvironment, the specific role of ferroptosis and pyroptosis in hepatocellular carcinoma (HCC) is not fully understood, especially its impact on prognosis.

Methods: The training set (n = 609, merged by TCGA and GSE14520) was clustered into three subtypes (C1, C2, and C3) based on the prognosis-related genes associated with ferroptosis and pyroptosis. The intersecting differentially expressed genes (DEGs) among C1, C2, and C3 were used in univariate Cox and LASSO penalized Cox regression analysis for the construction of the risk score. The median risk score served as the unified cutoff to divide patients into high- and low-risk groups.

Results: Internal (TCGA, n = 370; GSE14520, n = 239) and external validation (ICGC, n = 231) suggested that the 12-gene risk score had high accuracy in predicting the OS, DSS, DFS, PFS, and RFS of HCC. As an independent prognostic indicator, the risk score could be applicable for patients with different clinical features tested by subgroup (n = 26) survival analysis. In the high-risk patients with a lower infiltration abundance of activated B cells, activated CD8 T cells, eosinophils, and type I T helper cells and a higher infiltration abundance of immature dendritic cells, the cytolytic activity, HLA, inflammation promotion, and type I IFN response in the high-risk group were weaker. The TP53 mutation rate, TMB, and CSC characteristics in the high-risk group were significantly higher than those in the low-risk group. Low-risk patients have active metabolic activity and a more robust immune response. The high- and low-risk groups differed significantly in histology grade, vascular tumor cell type, AFP, new tumor event after initial treatment, main tumor size, cirrhosis, TNM stage, BCLC stage, and CLIP score.

Conclusion: The ferroptosis and pyroptosis molecular subtype-related signature identified and validated in this work is applicable for prognosis prediction, immune microenvironment estimation, stem cell characteristics, and clinical feature assessment in HCC.

Identification and validation of a novel immune-related signature associated with macrophages and CD8 T cell infiltration predicting overall survival for hepatocellular carcinoma

BACKGROUND

Although the effects of macrophages and CD8 T cell infiltration on clinical outcome have been widely reported, the association between immunity-associated gene with them for hepatocellular carcinoma (HCC) remains unclear.

MATERIALS AND METHODS

The ssGSEA served for quantifying the macrophages as well as CD8 T cell infiltration in the HCC samples obtained from TCGA database. Kaplan-Meier (KM) survival assay was used to determine the associations between macrophages and CD8 T cell infiltration with OS. LASSO Cox regressive method assisted in developing an immune gene signature as well as building a risk score. The performance was evaluated by the time-dependent ROC together with the KM survival analysis. The ICGC database were adopted for external verification. CIBERSORT was applied to the correlation analysis on the immune-related signature and the immunocyte infiltration. GSEA were employed exploring the underlying molecular mechanisms.

RESULTS

Increased CD8+ T cell infiltration was associated with longer OS, whereas a greater infiltration of macrophages was related to shorter OS. There were 398 differential expression genes (DEGs) between the high- and low infiltration groups with the "edgeR" package. An prognostic signature consisted of 10 immune genes was built in TCGA and examined in ICGC. The uniform cutoff (0.927) was adopted for separating sufferers into the high-risk (HR) and low-risk (LR) groups. The ROC curves revealed that the AUC data for this signature predicting 1, 2, 3, 4 and 5 year were all above 0.7 in both TCGA and ICGC cohort and patients in the HR group exhibited an evidently weaker prognostic results compared with the LR group. The HR group presented evidently greater Tregs and Macrophage M0 relative to the LR group, whereas the LR group saw the enrichment of CD8 T cells.

CONCLUSION

The immune signature associated with macrophages as well as CD8 T cell infiltration has reliable prognostic and predictive value for HCC patients.