Background: Ovarian metastasis (OM) poses a major threat to the outcome of gastric cancer (GC) patients. Recently, immunotherapy emerged as a novel promising therapeutic strategy to treat late-stage GC, whereas its efficacy is influenced by tumor immune microenvironment (TIME). M2 macrophage, a key subset within TIME, plays dual immunosuppressive and pro-tumorigenic roles in cancer progression and is recognized as a potential therapeutic target. However, molecular mechanisms underlying OM remain elusive and the TIME-related prognostic and immunotherapeutic index for these patients is yet to establish.
Methods: Differential expressed genes (DEGs) between paired normal mucosa, primary GC and OM of patients from Fudan University Shanghai Cancer Center (FUSCC) cohort (n = 6) were identified by transcriptome sequencing, followed by the functional annotation of enriched hallmark pathways of DEGs between them. CIBERSORT was used to profile the relative expression level of 22 immune cell subsets in normal tissues, primary and metastatic tumors, followed by weighted gene coexpression network analysis (WGCNA) uncovering immune cell-correlated gene sets. The intersected genes between DEGs and M2 macrophage-related genes were processed by least absolute shrinkage and selection operator (LASSO) regression analysis to construct a predictive signature, M2GO, which was further validated by training set and test set of The Cancer Genome Atlas-Stomach Adenocarcinoma (TCGA-STAD), GSE62254 and GSE84437 cohorts. GC patients were divided into M2GO-high and -low subgroup according to the optimal cutoff value of the M2GO score. Furthermore, the clinical, molecular and immune features between M2GO-high and -low subgroups were analyzed. Clinical cohorts of immunotherapy were used to validate the predictive value of M2GO in regard to immunotherapy effectiveness.
Results: Transcriptomic sequencing and follow-up analyses of triple-matched normal tissues, primary and ovarian metastatic tumors identified distinctive sets of DEGs and enriched immune-, cancer- and metastasis-related pathways between them. Of note, M2 macrophage, a major immunosuppressive and pro-tumorigenic component within TIME, was significantly up-regulated in OMs. WGCNA and LASSO regression were applied to establish a novel OM- and M2 macrophage-related predictive signature, M2GO, based on M2 macrophage-related prognostic genes including GJA1, MAGED1 and SERPINE1. M2GO served as an independent prognostic factor of GC patients. Comprehensive molecular and immune characterization of M2GO-based subgroups uncovered their distinctive features in terms of enriched functional pathways, tumor mutation burden, key immune checkpoints, major regulators of natural immune cGAS-STING pathway, infiltrated subsets of immune cells and tumor immune exclusion/dysfunction (TIDE) score. Notably, the M2GO score was significantly lower in responsive group than non-responsive group (P < 0.05) in clinical cohort of metastatic GC patients undergoing immunotherapy.
Conclusion: Transcriptomic characterization of paired normal mucosae, primary and ovarian metastatic tumors revealed their unique molecular and immune features. Follow-up analyses established a novel OM- and M2 macrophage-related signature, M2GO, which served as a promising prognostic and immunotherapeutic biomarker to distinguish the clinical outcome, molecular and immune features of GC patients and predict their differential responses to immunotherapy.
Keywords: Gastric cancer; Immunotherapy; M2 macrophage; Ovarian metastasis; Prognostic signature; Transcriptome sequencing.
© 2022. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.