Gastric cancer heterogeneity represents a barrier to disease management. We generated a comprehensive single-cell atlas of gastric cancer (>200,000 cells) comprising 48 samples from 31 patients across clinical stages and histologic subtypes. We identified 34 distinct cell-lineage states including novel rare cell populations. Many lineage states exhibited distinct cancer-associated expression profiles, individually contributing to a combined tumor-wide molecular collage. We observed increased plasma cell proportions in diffuse-type tumors associated with epithelial-resident KLF2 and stage-wise accrual of cancer-associated fibroblast subpopulations marked by high INHBA and FAP coexpression. Single-cell comparisons between patient-derived organoids (PDO) and primary tumors highlighted inter- and intralineage similarities and differences, demarcating molecular boundaries of PDOs as experimental models. We complemented these findings by spatial transcriptomics, orthogonal validation in independent bulk RNA-sequencing cohorts, and functional demonstration using in vitro and in vivo models. Our results provide a high-resolution molecular resource of intra- and interpatient lineage states across distinct gastric cancer subtypes. Significance: We profiled gastric malignancies at single-cell resolution and identified increased plasma cell proportions as a novel feature of diffuse-type tumors. We also uncovered distinct cancer-associated fibroblast subtypes with INHBA–FAP-high cell populations as predictors of poor clinical prognosis. Our findings highlight potential origins of deregulated cell states in the gastric tumor ecosystem. This article is highlighted in the In This Issue feature, p. 587

Almost all gastric cancers are adenocarcinomas, which have considerable heterogeneity among patients. We sought to identify subtypes of gastric adenocarcinomas with particular biological properties and responses to chemotherapy and targeted agents. We compared gene expression patterns among 248 gastric tumors; using a robust method of unsupervised clustering, consensus hierarchical clustering with iterative feature selection, we identified 3 major subtypes. We developed a classifier for these subtypes and validated it in 70 tumors from a different population. We identified distinct genomic and epigenomic properties of the subtypes. We determined drug sensitivities of the subtypes in primary tumors using clinical survival data, and in cell lines through high-throughput drug screening. We identified 3 subtypes of gastric adenocarcinoma: proliferative, metabolic, and mesenchymal. Tumors of the proliferative subtype had high levels of genomic instability, TP53 mutations, and DNA hypomethylation. Cancer cells of the metabolic subtype were more sensitive to 5-fluorouracil than the other subtypes. Furthermore, in 2 independent groups of patients, those with tumors of the metabolic subtype appeared to have greater benefits with 5-fluorouracil treatment. Tumors of the mesenchymal subtype contain cells with features of cancer stem cells, and cell lines of this subtype are particularly sensitive to phosphatidylinositol 3-kinase−AKT−mTOR inhibitors in vitro. Based on gene expression patterns, we classified gastric cancers into 3 subtypes, and validated these in an independent set of tumors. The subgroups have differences in molecular and genetic features and response to therapy; this information might be used to select specific treatment approaches for patients with gastric cancer.

Many solid cancers are known to exhibit a high degree of heterogeneity in their deregulation of different oncogenic pathways. We sought to identify major oncogenic pathways in gastric cancer (GC) with significant relationships to patient survival. Using gene expression signatures, we devised an in silico strategy to map patterns of oncogenic pathway activation in 301 primary gastric cancers, the second highest cause of global cancer mortality. We identified three oncogenic pathways (proliferation/stem cell, NF-κB, and Wnt/β-catenin) deregulated in the majority (>70%) of gastric cancers. We functionally validated these pathway predictions in a panel of gastric cancer cell lines. Patient stratification by oncogenic pathway combinations showed reproducible and significant survival differences in multiple cohorts, suggesting that pathway interactions may play an important role in influencing disease behavior. Individual GCs can be successfully taxonomized by oncogenic pathway activity into biologically and clinically relevant subgroups. Predicting pathway activity by expression signatures thus permits the study of multiple cancer-related pathways interacting simultaneously in primary cancers, at a scale not currently achievable by other platforms. Gastric cancer is the second leading cause of global cancer mortality. With current treatments, less than a quarter of patients survive longer than five years after surgery. Individual gastric cancers are highly disparate in their cellular characteristics and responses to standard chemotherapeutic drugs, making gastric cancer a complex disease. Pathway based approaches, rather than single gene studies, may help to unravel this complexity. Here, we make use of a computational approach to identify connections between molecular pathways and cancer profiles. In a large scale study of more than 300 patients, we identified subgroups of gastric cancers distinguishable by their patterns of driving molecular pathways. We show that these identified subgroups are clinically relevant in predicting survival duration and may prove useful in guiding the choice of targeted therapies designed to interfere with these molecular pathways. We also identified specific gastric cancer cell lines mirroring these pathway subgroups, which should facilitate the pre-clinical assessment of responses to targeted therapies in each subgroup.

Despite multidisciplinary treatment for patients with advanced gastric cancer, their prognosis remains poor. Therefore, the development of novel therapeutic strategies is urgently needed, and immunotherapy utilizing anti‐programmed death 1/‐programmed death ligand‐1 mAb is an attractive approach. However, as there is limited information on how programmed death ligand‐1 is upregulated on tumor cells within the tumor microenvironment, we examined the mechanism of programmed death ligand‐1 regulation with a particular focus on interferon gamma in an in vitro setting and in clinical samples. Our in vitro findings showed that interferon gamma upregulated programmed death ligand‐1 expression on solid tumor cells through the JAK‐signal transducer and activator of transcription pathway, and impaired the cytotoxicity of tumor antigen‐specific CTL against tumor cells. Following treatment of cells with anti‐programmed death ligand‐1 mAb after interferon gamma‐pre‐treatment, the reduced anti‐tumor CTL activity by interferon gamma reached a higher level than the non‐treatment control targets. In contrast, programmed death ligand‐1 expression on tumor cells also significantly correlated with epithelial‐mesenchymal transition phenotype in a panel of solid tumor cells. In clinical gastric cancer samples, tumor membrane programmed death ligand‐1 expression significantly positively correlated with the presence of CD8‐positive T cells in the stroma and interferon gamma expression in the tumor. The results suggest that gastric cancer patients with high CD8‐positive T‐cell infiltration may be more responsive to anti‐programmed death 1/‐programmed death ligand‐1 mAb therapy.

Background & Aims Gastric cancer (GC) is a heterogeneous disease comprising multiple subtypes that have distinct biological properties and effects in patients. We sought to identify new, intrinsic subtypes of GC by gene expression analysis of a large panel of GC cell lines. We tested if these subtypes might be associated with differences in patient survival times and responses to various standard-of-care cytotoxic drugs. Methods We analyzed gene expression profiles for 37 GC cell lines to identify intrinsic GC subtypes. These subtypes were validated in primary tumors from 521 patients in 4 independent cohorts, where the subtypes were determined by either expression profiling or subtype-specific immunohistochemical markers (LGALS4, CDH17). In vitro sensitivity to 3 chemotherapy drugs (5-fluorouracil, cisplatin, oxaliplatin) was also assessed. Results Unsupervised cell line analysis identified 2 major intrinsic genomic subtypes (G-INT and G-DIF) that had distinct patterns of gene expression. The intrinsic subtypes, but not subtypes based on Lauren's histopathologic classification, were prognostic of survival, based on univariate and multivariate analysis in multiple patient cohorts. The G-INT cell lines were significantly more sensitive to 5-fluorouracil and oxaliplatin, but more resistant to cisplatin, than the G-DIF cell lines. In patients, intrinsic subtypes were associated with survival time following adjuvant, 5-fluorouracil–based therapy. Conclusions Intrinsic subtypes of GC, based on distinct patterns of expression, are associated with patient survival and response to chemotherapy. Classification of GC based on intrinsic subtypes might be used to determine prognosis and customize therapy.