Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
GASTRIC CANCER TUMOR MICROENVIRONMENTS
Document Type and Number:
WIPO Patent Application WO/2022/192399
Kind Code:
A1
Abstract:
Techniques for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk of having gastric cancer, whether the subject is likely to respond to an immunotherapy. The techniques include: obtaining RNA expression data for the subject; generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in a plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression data; identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy.

Inventors:
KUDRYASHOVA OLGA (RU)
SHAH MANISH (US)
KOTLOV NIKITA (RU)
MELIKHOVA DARIA (RU)
GUSAKOVA MARIIA (RU)
SAMARINA NAIRA (RU)
PODSVIROVA SVETLANA (RU)
TYCHININ DMITRY (RU)
Application Number:
PCT/US2022/019546
Publication Date:
September 15, 2022
Filing Date:
March 09, 2022
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
BOSTONGENE CORP (US)
BOSTONGENE LLC (RU)
WEILL CORNELL MEDICAL COLLEGE (US)
UNIV CORNELL (US)
KUDRYASHOVA OLGA (RU)
SHAH MANISH (US)
KOTLOV NIKITA (RU)
MELIKHOVA DARIA (RU)
GUSAKOVA MARIIA (RU)
SAMARINA NAIRA (RU)
PODSVIROVA SVETLANA (RU)
TYCHININ DMITRY (RU)
International Classes:
C12Q1/6886
Foreign References:
CN112133365A2020-12-25
US20200273543A12020-08-27
Other References:
ZENG DONGQIANG ET AL: "Tumor Microenvironment Characterization in Gastric Cancer Identifies Prognostic and Immunotherapeutically Relevant Gene Signatures", CANCER IMMUNOLOGY RESEARCH, vol. 7, no. 5, 1 May 2019 (2019-05-01), US, pages 737 - 750, XP055930069, ISSN: 2326-6066, Retrieved from the Internet DOI: 10.1158/2326-6066.CIR-18-0436
CHO JUNHUN ET AL: "Four distinct immune microenvironment subtypes in gastric adenocarcinoma with special reference to microsatellite instability", ESMO OPEN : CANCER HORIZONS, vol. 3, no. 3, 1 January 2018 (2018-01-01), London, pages e000326, XP055931360, ISSN: 2059-7029, DOI: 10.1136/esmoopen-2018-000326
O'DONNELL JAKE S ET AL: "Cancer immunoediting and resistance to T cell-based immunotherapy", NATURE REVIEWS CLINICAL ONCOLOGY, NATURE, NY, US, vol. 16, no. 3, 6 December 2018 (2018-12-06), pages 151 - 167, XP036706114, ISSN: 1759-4774, [retrieved on 20181206], DOI: 10.1038/S41571-018-0142-8
MA ET AL., ONCOL LETT, vol. 11, no. 5, May 2016 (2016-05-01), pages 2959 - 2964
HENNEQUIN ET AL., ONCOLMMUNOLOGY, vol. 5, 2016, pages 2
VAUGHT ET AL., CANCER EPIDEMIOL BIOMARKERS PREV, vol. 21, no. 2, February 2012 (2012-02-01), pages 253 - 5
VAUGHTHENDERSON, IARC SCI PUBL, no. 163, 2011, pages 23 - 42
NICOLAS L BRAYHAROLD PIMENTELPALL MELSTEDLIOR PACHTER: "Near-optimal probabilistic RNA-seq quantification", NATURE BIOTECHNOLOGY, vol. 34, 2016, pages 525 - 527
BIOINFORMATICS, vol. 20, no. 3, 12 February 2004 (2004-02-12), pages 307 - 15
RITCHIE MEPHIPSON BWU DHU YLAW CWSHI WSMYTH GK: "limma powers differential expression analyses for RNA-sequencing and microarray studies", NUCLEIC ACIDS RES., vol. 43, no. 7, 20 April 2015 (2015-04-20), pages e47
WAGNER ET AL., THEORY BIOSCI, vol. 131, 2012, pages 281 - 285
BARBIE ET AL., NATURE, vol. 462, no. 7269, 5 November 2009 (2009-11-05), pages 108 - 112
SAUTES-FRIDMAN ET AL., NAT REV CANCER, vol. 19, 2019, pages 307 - 325
Attorney, Agent or Firm:
RUDOY, Daniel, G. et al. (US)
Download PDF:
Claims:
CLAIMS What is claimed is: 1. A method for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk of having gastric cancer, whether the subject is likely to respond to an immunotherapy, the method comprising: using at least one computer hardware processor to perform: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; (c) identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and (d) identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy. 2. The method of claim 1, wherein the RNA expression data further indicates second RNA expression levels for at least some genes in each group of at least some of a second plurality of gene groups listed in Table 2, wherein the GC TME signature further comprises second gene group scores for respective gene groups of the at least some of the second plurality of gene groups, and wherein the generating further comprises determining the second gene group scores using the second RNA expression levels. 3. The method of claim 1, wherein the RNA expression data further indicates second RNA expression levels for at least some genes in each group of at least some of a second plurality of gene groups listed in Table 2, and the method further comprises: (e) generating a second GC TME signature for the subject using the RNA expression data, the second GC TME signature comprising second gene group scores for respective gene groups in the at least some of the second plurality of gene groups, the generating comprising: determining the second gene group scores using the second RNA expression levels. 4. The method of any one of claims 1 to 3, wherein obtaining the RNA expression data for the subject comprises obtaining sequencing data previously obtained by sequencing a biological sample obtained from the subject. 5. The method of claim 4, wherein the sequencing data comprises at least 1 million reads, at least 5 million reads, at least 10 million reads, at least 20 million reads, at least 50 million reads, or at least 100 million reads. 6. The method of claim 4 or 5, wherein the sequencing data comprises whole exome sequencing (WES) data, bulk RNA sequencing (RNA-seq) data, single cell RNA sequencing (scRNA-seq) data, or next generation sequencing (NGS) data. 7. The method of claim 4 or 5, wherein the sequencing data comprises microarray data. 8. The method of any one of claims 1 to 7, further comprising: normalizing the RNA expression data to transcripts per million (TPM) units prior to generating the GC TME signature. 9. The method of any one of claims 1 to 8, wherein obtaining the RNA expression data for the subject comprises sequencing a biological sample obtained from the subject. 10. The method of claim 9, wherein the biological sample comprises gastrointestinal tissue of the subject, optionally wherein the biological sample comprises tumor tissue of the subject. 11. The method of any one of claims 1 to 10, wherein the RNA expression levels comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 12. The method of any one of claims 1 to 11, wherein the RNA expression levels comprise RNA expression levels for each of the genes from each of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 13. The method of any one of claims 1 to 12, wherein the RNA expression levels comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA- DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93.

14. The method of any one of claims 1 to 13, wherein determining the gene group scores comprises: determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 15. The method of any one of claims 1 to 14, wherein determining the gene group scores comprises: determining a respective gene group score for each of the following gene groups, using, for each gene group, RNA expression levels for each of the genes in each gene group to determine the gene group score for each particular group, the gene groups including: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 16. The method of claim 15, wherein determining the gene group scores comprises determining a first score of a first gene group using a single-sample GSEA (ssGSEA) technique from RNA expression levels for at least some of the genes in one of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 17. The method of any one of claims 1 to 16, wherein determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 18. The method of any one of claims 1 to 17, wherein determining the gene group scores comprises: determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA- DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 19. The method of claim 18, wherein determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for at least some of the genes in each one of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA- DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. 20. The method of any one of claims 1 to 19, wherein determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA- DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93.

21. The method of any one of claims 1 to 20, wherein generating the GC TME signature further comprises normalizing the gene group scores, wherein the normalizing comprises applying rank estimation and/or median scaling to the gene group scores. 22. The method of any one of claims 1 to 21, wherein the plurality of GC TME types is associated with a respective plurality of GC TME signature clusters, wherein identifying, using the GC TME signature, and from among a plurality of GC TME types, the GC TME type for the subject comprises: associating the GC TME signature of the subject with a particular one of the plurality of GC TME signature clusters; and identifying the GC TME type for the subject as the GC TME type corresponding to the particular one of the plurality of GC TME signature clusters to which the GC TME signature of the subject is associated. 23. The method of claim 22, further comprising generating the plurality of GC TME signature clusters, the generating comprising: obtaining multiple sets of RNA expression data by sequencing biological samples from multiple respective subjects, each of the multiple sets of RNA expression data indicating RNA expression levels for at least some genes in each of the at least some of the plurality of gene groups listed in Table 1; generating multiple GC TME signatures from the multiple sets of RNA expression data, each of the multiple GC TME signatures comprising gene group scores for respective gene groups in the plurality of gene groups, the generating comprising, for each particular one of the multiple GC TME signatures, determining the GC TME signature by determining the gene group scores using the RNA expression levels in the particular set of RNA expression data for which the particular one GC TME signature is being generated; and clustering the multiple GC signatures to obtain the plurality of GC TME signature clusters. 24. The method of claim 23, wherein the clustering comprises dense clustering, spectral clustering, k-means clustering, hierarchical clustering, and/or an agglomerative clustering.

25. The method of claim 24, wherein the hierarchical clustering is performed using a Louvain community detection algorithm. 26. The method of any one of claims 23 to 25, further comprising: updating the plurality of GC TME signature clusters using the GC TME signature of the subject, wherein the GC TME signature of the subject is one of a threshold number GC TME signatures for a threshold number of subjects, wherein when the threshold number of GC TME signatures is generated the GC TME signature clusters are updated. 27. The method of claim 26, wherein the threshold number of GC TME signatures is at least 50, at least 75, at least 100, at least 200, at least 500, at least 1000, or at least 5000 GC TME signatures. 28. The method of claim 27, wherein the updating comprises applying dense clustering, spectral clustering, k-means clustering, hierarchical clustering, and/or agglomerative clustering. 29. The method of claim 28, wherein the hierarchical clustering is performed using a Louvain community detection algorithm. 30. The method of any one of claims 23 to 29, further comprising: determining an GC TME type of a second subject, wherein the GC TME type of the second subject is identified using the updated GC TME signature clusters, wherein the identifying comprises: determining an GC TME signature of the second subject from RNA expression data obtained by sequencing a biological sample obtained from the second subject; associating the GC TME signature of the second subject with a particular one of the plurality of the updated GC TME signature clusters; and identifying the GC TME type for the second subject as the GC TME type corresponding to the particular one of the plurality of updated GC TME signature clusters to which the GC TME signature of the second subject is associated.

31. The method of any one of claims 1 to 30, wherein the plurality of a plurality of GC TME types comprises: GC TME type A, GC TME type B, GC TME type C, GC TME type D, and GC TME type E. 32. The method of any one of claims 1 to 31, further comprising: identifying the subject as having tertiary lymphoid structures (TLS) when the subject is identified as having GC TME type E. 33. The method of any one of claims 1 to 32, further comprising: identifying the subject as having an increased likelihood of having a good prognosis, optionally, as measured by overall survival (OS) or progression-free survival (PFS), when the subject is identified as having GC TME type E. 34. The method of any one of claims 3 to 33, wherein the RNA expression levels for genes in the second plurality of gene groups comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2.

35. The method of claim 34, wherein the RNA expression levels for genes in the second plurality of gene groups comprise RNA expression levels for each of the genes from each of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. 36. The method of claim 34 or 35, wherein determining the gene group scores comprises determining a first score of a first gene group using a single-sample GSEA (ssGSEA) technique from second RNA expression levels for at least some of the genes in one of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. 37. The method of any one of claims 34 to 36, wherein determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each one of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. 38. The method of any one of claims 34 to 37, wherein the identifying in (c) is performed using the second GC TME signature. 39. The method of any one of claims 34 to 38, wherein the identifying in (d) is performed using the second GC TME signature. 40. The method of any one of claims 1 to 39, further comprising: identifying the subject as having an increased likelihood of responding to an immunotherapy when the subject is assigned GC TME type B.

41. The method of any one of claims 1 to 39, further comprising: identifying the subject as having an increased likelihood of responding to an immunotherapy when the subject is assigned GC TME type E. 42. The method of claim 40 or 41, wherein the subject is Microsatellite Stable (MSS). 43. The method of any one of claims 40 to 42, wherein the subject is Epstein Barr virus (EBV) low or EBV negative. 44. The method of any one of claims 1 to 43 further comprising: administering an immunotherapy to the subject. 45. The method of any one of claims 1 to 44, wherein the immunotherapy comprises a PD1 inhibitor, optionally pembrolizumab. 46. A method for predicting disease control rate (DCR) to an immunotherapy of a subject, the method comprising: using at least one computer hardware processor to perform: (i) identifying the GC TME of the subject using the method of any one of the preceding claims; and (ii) identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D. 47. The method of claim 46, wherein the immunotherapy comprises a PD1 inhibitor, optionally pembrolizumab. 48. The method of claim 46 or 47, wherein the DCR of a subject having GC TME type B is greater than 50%.

49. The method of claim 46 or 47, wherein the DCR of a subject having GC TME type E is greater than 50%. 50. The method of any one of claims 46 to 49, wherein the subject is Microsatellite Stable (MSS). 51. The method of any one of claims 46 to 50, wherein the subject is Epstein Barr virus (EBV) low or EBV negative. 52. The method of any one of claims 46 to 51, further comprising administering an immunotherapy to the subject. 53. A system, comprising: at least one computer hardware processor; and at least one computer-readable storage medium storing processor-executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk, of having a gastric cancer, whether the subject is likely to respond to an immunotherapy, the method comprising: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; (c) identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and (d) identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy.

54. At least one computer-readable storage medium storing processor-executable instructions that, when executed by at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk, of having a gastric cancer, whether the subject is likely to respond to an immunotherapy, the method comprising: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; (c) identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and (d) identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy. 55. A system, comprising: at least one computer hardware processor; and at least one non-transitory computer-readable storage medium, storing processor executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for predicting disease control rate (DCR) to an immunotherapy of a subject, the method comprising: (i) identifying the GC TME of the subject using the method of any one of the preceding claims; and (ii) identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D.

56. At least one non-transitory computer-readable storage medium, storing processor executable instructions that, when executed by at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for predicting disease control rate (DCR) to an immunotherapy of a subject, the method comprising: (i) identifying the GC TME of the subject using the method of any one of the preceding claims; and (ii) identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D.

Description:
B CELL-ENRICHED TUMOR MICROENVIRONMENTS CROSS REFERENCE TO RELATED APPLICATIONS This application claims the benefit under 35 U.S.C. § 119(e) of the filing date of U.S. Provisional Application No.63/158,818, filed March 9, 2021, entitled “B CELL-ENRICHED TUMOR MICROENVIRONMENTS,” the entire contents of which are incorporated by reference herein. BACKGROUND Correctly characterizing the type or types of cancer a patient or subject has and, potentially, selecting one or more effective therapies for the patient can be crucial for the survival and overall wellbeing of that patient. Advances in characterizing cancers, predicting prognoses, identifying effective therapies, and otherwise aiding in personalized care of patients with cancer are needed. SUMMARY Aspects of the disclosure relate to techniques for characterizing subjects having gastric cancers. The disclosure is based, in part, on methods for identifying the tumor microenvironment (TME) of a subject having gastric cancer by using RNA expression data obtained from the subject to produce a gastric cancer (GC) tumor microenvironment (TME) signature that allows for assignment of a GC TME type to the subject. In some embodiments, the GC TME type of a subject is indicative of one or more characteristics of the subject (or the subject’s cancer), for example the likelihood a subject will have a good prognosis and/or that the subject is likely to respond to treatment with certain therapeutic agents, such as immunotherapeutic agents (e.g., immune checkpoint inhibitors). Accordingly, in some aspects, the disclosure provides a method for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk of having gastric cancer, whether the subject is likely to respond to an immunotherapy, comprising: using at least one computer hardware processor to perform: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; (c) identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and (d) identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy. Aspects of the disclosure include a method for predicting disease control rate (DCR) to an immunotherapy of a subject, comprising: using at least one computer hardware processor to perform: (i) identifying the GC TME of the subject using the method of any one of the preceding claims; and (ii) identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D. Aspects of the disclosure include a system, comprising at least one computer hardware processor; and at least one computer-readable storage medium storing processor-executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk, of having a gastric cancer, whether the subject is likely to respond to an immunotherapy, the method comprising: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; (c) identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and (d) identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy. Aspects of the disclosure include at least one computer-readable storage medium storing processor-executable instructions that, when executed by at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for identifying, based at least in part on a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk, of having a gastric cancer, whether the subject is likely to respond to an immunotherapy, comprising: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; (c)identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject; and (d) identifying, using the GC TME type of the subject, whether or not the subject is likely to respond to the immunotherapy. Aspects of the disclosure include a system, comprising at least one computer hardware processor; and at least one non-transitory computer-readable storage medium, storing processor executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for predicting disease control rate (DCR) to an immunotherapy of a subject, the method comprising: (i) identifying the GC TME of the subject using the method of any one of the preceding claims; and (ii) identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D. Aspects of the disclosure include at least one non-transitory computer-readable storage medium, storing processor executable instructions that, when executed by at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for predicting disease control rate (DCR) to an immunotherapy of a subject, comprising: (i) identifying the GC TME of the subject using the method of any one of the preceding claims; and (ii) identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D. In some embodiments, the RNA expression data further indicates second RNA expression levels for at least some genes in each group of at least some of a second plurality of gene groups listed in Table 2, wherein the GC TME signature further comprises second gene group scores for respective gene groups of the at least some of the second plurality of gene groups, and wherein the generating further comprises determining the second gene group scores using the second RNA expression levels. In some embodiments, the RNA expression data further indicates second RNA expression levels for at least some genes in each group of at least some of a second plurality of gene groups listed in Table 2, and the method further comprises: (e) generating a second GC TME signature for the subject using the RNA expression data, the second GC TME signature comprising second gene group scores for respective gene groups in the at least some of the second plurality of gene groups, the generating comprising: determining the second gene group scores using the second RNA expression levels. In some embodiments, obtaining the RNA expression data for the subject comprises obtaining sequencing RNA data previously obtained by sequencing a biological sample obtained from the subject. In some embodiments, the sequencing data comprises at least 1 million reads, at least 5 million reads, at least 10 million reads, at least 20 million reads, at least 50 million reads, or at least 100 million reads. In some embodiments, the sequencing data comprises whole exome sequencing (WES) data, bulk RNA sequencing (RNA-seq) data, single cell RNA sequencing (scRNA-seq) data, or next generation sequencing (NGS) data. In some embodiments, the sequencing data comprises microarray data. In some embodiments, the method described herein further comprises normalizing the RNA expression data to transcripts per million (TPM) units prior to generating the GC TME signature. In some embodiments, obtaining the RNA expression data for the subject comprises sequencing a biological sample obtained from the subject. In some embodiments, the biological sample comprises gastrointestinal tissue of the subject. In some embodiments, the biological sample comprises tumor tissue of the subject. In some embodiments, the RNA expression levels comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, the RNA expression levels comprise RNA expression levels for each of the genes from each of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, the RNA expression levels comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (a) MHC I group: HLA- A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA- DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises: determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises: determining a respective gene group score for each of the following gene groups, using, for each gene group, RNA expression levels for each of the genes in each gene group to determine the gene group score for each particular group, the gene groups including: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining a first score of a first gene group using a single-sample GSEA (ssGSEA) technique from RNA expression levels for at least some of the genes in one of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises: determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA- DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for at least some of the genes in each one of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA- DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, generating the GC TME signature further comprises normalizing the gene group scores, wherein the normalizing comprises applying rank estimation and/or median scaling to the gene group scores. In some embodiments, the plurality of GC TME types is associated with a respective plurality of GC TME signature clusters, wherein identifying, using the GC TME signature, and from among a plurality of GC TME types, the GC TME type for the subject comprises: associating the GC TME signature of the subject with a particular one of the plurality of GC TME signature clusters; and identifying the GC TME type for the subject as the GC TME type corresponding to the particular one of the plurality of GC TME signature clusters to which the GC TME signature of the subject is associated. In some embodiments, the method described herein further comprises generating the plurality of GC TME signature clusters, comprising: obtaining multiple sets of RNA expression data by sequencing biological samples from multiple respective subjects, each of the multiple sets of RNA expression data indicating RNA expression levels for at least some genes in each of the at least some of the plurality of gene groups listed in Table 1; generating multiple GC TME signatures from the multiple sets of RNA expression data, each of the multiple GC TME signatures comprising gene group scores for respective gene groups in the plurality of gene groups, comprising: for each particular one of the multiple GC TME signatures, determining the GC TME signature by determining the gene group scores using the RNA expression levels in the particular set of RNA expression data for which the particular one GC TME signature is being generated; and clustering the multiple GC signatures to obtain the plurality of GC TME signature clusters. In some embodiments, the clustering comprises dense clustering, spectral clustering, k- means clustering, hierarchical clustering, and/or an agglomerative clustering. In some embodiments, the hierarchical clustering is performed using a Louvain community detection algorithm. In some embodiments, the method described herein further comprises updating the plurality of GC TME signature clusters using the GC TME signature of the subject, wherein the GC TME signature of the subject is one of a threshold number GC TME signatures for a threshold number of subjects, wherein when the threshold number of GC TME signatures is generated the GC TME signature clusters are updated. In some embodiments, the threshold number of GC TME signatures is at least 50, at least 75, at least 100, at least 200, at least 500, at least 1000, or at least 5000 GC TME signatures. In some embodiments, the updating comprises applying dense clustering, spectral clustering, k-means clustering, hierarchical clustering, and/or agglomerative clustering. In some embodiments, the hierarchical clustering is performed using a Louvain community detection algorithm. In some embodiments, the method described herein further comprises determining an GC TME type of a second subject, wherein the GC TME type of the second subject is identified using the updated GC TME signature clusters, wherein the identifying comprises: determining an GC TME signature of the second subject from RNA expression data obtained by sequencing a biological sample obtained from the second subject; associating the GC TME signature of the second subject with a particular one of the plurality of the updated GC TME signature clusters; and identifying the GC TME type for the second subject as the GC TME type corresponding to the particular one of the plurality of updated GC TME signature clusters to which the GC TME signature of the second subject is associated. In some embodiments, the plurality of a plurality of GC TME types comprises: GC TME type A, GC TME type B, GC TME type C, GC TME type D, and GC TME type E. In some embodiments, the method described herein further comprises identifying the subject as having tertiary lymphoid structures (TLS) when the subject is identified as having GC TME type E. In some embodiments, the method described herein further comprises identifying the subject as having an increased likelihood of having a good prognosis. In some embodiments, the increased likelihood of having a good prognosis is as measured by overall survival (OS) or progression-free survival (PFS), when the subject is identified as having GC TME type E. In some embodiments, the RNA expression levels for genes in the second plurality of gene groups comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. In some embodiments, the RNA expression levels for genes in the second plurality of gene groups comprise RNA expression levels for each of the genes from each of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. In some embodiments, the determining the gene group scores comprises determining a first score of a first gene group using a single-sample GSEA (ssGSEA) technique from second RNA expression levels for at least some of the genes in one of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. In some embodiments, the determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each one of the following gene groups: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressor group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) Myc targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2. In some embodiments, the identifying in (c) of the method described herein is performed using the second GC TME signature. In some embodiments, the identifying in (d) of the method described herein is performed using the second GC TME signature. In some embodiments, the method described herein further comprises identifying the subject as having an increased likelihood of responding to an immunotherapy when the subject is assigned GC TME type B. In some embodiments, the method described herein further comprises identifying the subject as having an increased likelihood of responding to an immunotherapy when the subject is assigned GC TME type E. In some embodiments, the subject is Microsatellite Stable (MSS). In some embodiments, the subject is Epstein Barr virus (EBV) low or EBV negative. In some embodiments, the method described herein further comprises administering an immunotherapy to the subject. In some embodiments, the immunotherapy comprises a PD1 inhibitor. In some embodiments, the PD1 inhibitor comprises pembrolizumab. In some embodiments, the DCR of a subject having GC TME type B is greater than 50%. In some embodiments, the DCR of a subject having GC TME type E is greater than 50%. Aspects of the disclosure relate to methods, systems, and computer-readable storage medium that can be used for identifying a gastric cancer (GC) tumor microenvironment (TME) type for a subject. In some aspects, the disclosure provides a method for identifying a gastric cancer (GC) tumor microenvironment (TME) type for a subject, comprising: using at least one computer hardware processor to perform: (a) obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; (b) generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; and (c) identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject. Aspects of the disclosure include a system, comprising: at least one computer hardware processor; and at least one computer-readable storage medium storing processor-executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for identifying a gastric cancer (GC) tumor microenvironment (TME) type for a subject, the method comprising: obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; and identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject. Aspects of the disclosure include at least one computer-readable storage medium storing processor-executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform a method for identifying a gastric cancer (GC) tumor microenvironment (TME) type for a subject, the method comprising: obtaining RNA expression data for the subject, the RNA expression data indicating RNA expression levels for at least some genes in each group of at least some of a plurality of gene groups listed in Table 1; generating a GC TME signature for the subject using the RNA expression data, the GC TME signature comprising gene group scores for respective gene groups in the at least some of the plurality of gene groups, the generating comprising: determining the gene group scores using the RNA expression levels; and identifying, using the GC TME signature and from among a plurality of GC TME types, a GC TME type for the subject. In some embodiments, the subject has, is suspected of having, or is at risk of having gastric cancer. In some embodiments, the method described herein further comprises identifying the subject as having GC TME type E; and when the subject is identified as having the GC TME type E, administering an immunotherapy to the subject. In some embodiments, obtaining the RNA expression data for the subject comprises obtaining sequencing RNA data previously obtained by sequencing a biological sample obtained from the subject. In some embodiments, the sequencing data comprises at least 1 million reads, at least 5 million reads, at least 10 million reads, at least 20 million reads, at least 50 million reads, or at least 100 million reads. In some embodiments, the sequencing data comprises whole exome sequencing (WES) data, bulk RNA sequencing (RNA-seq) data, single cell RNA sequencing (scRNA-seq) data, or next generation sequencing (NGS) data. In some embodiments, the sequencing data comprises microarray data. In some embodiments, the method described herein further comprises normalizing the RNA expression data to transcripts per million (TPM) units prior to generating the GC TME signature. In some embodiments, obtaining the RNA expression data for the subject comprises sequencing a biological sample obtained from the subject. In some embodiments, the biological sample comprises gastrointestinal tissue of the subject. In some embodiments, the biological sample comprises tumor tissue of the subject. In some embodiments, the RNA expression levels comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, the RNA expression levels comprise RNA expression levels for each of the genes from each of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, the RNA expression levels comprise RNA expression levels for at least three genes from each of at least two of the following gene groups: (a) MHC I group: HLA- A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA- DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises: determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises: determining a respective gene group score for each of the following gene groups, using, for each gene group, RNA expression levels for each of the genes in each gene group to determine the gene group score for each particular group, the gene groups including: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining a first score of a first gene group using a single-sample GSEA (ssGSEA) technique from RNA expression levels for at least some of the genes in one of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each of the following gene groups: (i) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (ii) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (iii) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (iv) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (v) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (vi) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (vii) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (viii) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises: determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA- DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for at least some of the genes in each one of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the gene group scores comprises determining the gene group scores, using a single-sample GSEA (ssGSEA) technique, from RNA expression levels for each of the genes in each of the following gene groups: (a) MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA- DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, generating the GC TME signature further comprises normalizing the gene group scores. In some embodiments, the normalizing comprises applying rank estimation and/or median scaling to the gene group scores. In some embodiments, the plurality of GC TME types is associated with a respective plurality of GC TME signature clusters, wherein identifying, using the GC TME signature and from among a plurality of GC TME types, the GC TME type for the subject comprises: associating the GC TME signature of the subject with a particular one of the plurality of GC TME signature clusters; and identifying the GC TME type for the subject as the GC TME type corresponding to the particular one of the plurality of GC TME signature clusters to which the GC TME signature of the subject is associated. In some embodiments, the method described herein further comprises generating the plurality of GC TME signature clusters, comprising: obtaining multiple sets of RNA expression data by sequencing biological samples from multiple respective subjects, each of the multiple sets of RNA expression data indicating RNA expression levels for at least some genes in each of the at least some of the plurality of gene groups listed in Table 1; generating multiple GC TME signatures from the multiple sets of RNA expression data, each of the multiple GC TME signatures comprising gene group scores for respective gene groups in the plurality of gene groups, the generating comprising, for each particular one of the multiple GC TME signatures, determining the GC TME signature by determining the gene group scores using the RNA expression levels in the particular set of RNA expression data for which the particular one GC TME signature is being generated; and clustering the multiple GC signatures to obtain the plurality of GC TME signature clusters. In some embodiments, the clustering comprises dense clustering, spectral clustering, k- means clustering, hierarchical clustering, and/or an agglomerative clustering. In some embodiments, the hierarchical clustering is performed using a Louvain community detection algorithm. In some embodiments, the method described herein further comprises updating the plurality of GC TME signature clusters using the GC TME signature of the subject, wherein the GC TME signature of the subject is one of a threshold number GC TME signatures for a threshold number of subjects, wherein when the threshold number of GC TME signatures is generated the GC TME signature clusters are updated. In some embodiments, the threshold number of GC TME signatures is at least 50, at least 75, at least 100, at least 200, at least 500, at least 1000, or at least 5000 GC TME signatures. In some embodiments, the updating comprises applying a dense clustering, spectral clustering, k-means clustering, hierarchical clustering, and/or agglomerative clustering. In some embodiments, the hierarchical clustering is performed using a Louvain community detection algorithm. In some embodiments, the method described herein further comprises determining an GC TME type of a second subject, wherein the GC TME type of the second subject is identified using the updated GC TME signature clusters, wherein the identifying comprises: determining an GC TME signature of the second subject from RNA expression data obtained by sequencing a biological sample obtained from the second subject; associating the GC TME signature of the second subject with a particular one of the plurality of the updated GC TME signature clusters; and identifying the GC TME type for the second subject as the GC TME type corresponding to the particular one of the plurality of updated GC TME signature clusters to which the GC TME signature of the second subject is associated. In some embodiments, the plurality of a plurality of GC TME types comprises: GC TME type A, GC TME type B, GC TME type C, GC TME type D, and GC TME type E. In some embodiments, the method described herein further comprises identifying the subject as having tertiary lymphoid structures (TLS) when the subject is identified as having GC TME type E. In some embodiments, the method described herein further comprises identifying the subject as having an increased likelihood of having a good prognosis. In some embodiments, increased likelihood of having a good prognosis is as measured by overall survival (OS) or progression-free survival (PFS) when the subject is identified as having GC TME type E. In some embodiments, the method described herein further comprises administering an immunotherapy to the subject. The immunotherapy may be administered when the subject is identified as having GC TME type E. The immunotherapy may be administered when the subject is identified as having TLS. In some embodiments, the immunotherapy comprises a PD1 inhibitor. In some embodiments, the PD1 inhibitor comprises pembrolizumab. BRIEF DESCRIPTION OF DRAWINGS FIG.1 is a diagram depicting a flowchart of an illustrative process for determining a gastric cancer (GC) tumor microenvironment (TME) type for a subject having, suspected of having, or at risk of having a GC, according to some embodiments of the technology as described herein. FIG.2 is a diagram depicting a flowchart of an illustrative process for processing sequencing data to obtain RNA expression data, according to some embodiments of the technology as described herein. FIG.3 is a diagram depicting an illustrative technique for determining gene group scores, according to some embodiments of the technology as described herein. FIG.4 is a diagram depicting an illustrative technique for identifying a GC TME type for a subject using a GC TME signature, according to some embodiments of the technology as described herein. FIG.5A provides an exemplary heatmap of GC samples classified into five distinct GC TME type (A, B, C, D, E) using example GC TME signatures comprising 20 gene group scores, according to some aspects of the technology described herein. Gene groups are listed vertically on the left side of the heatmap. Each column represents one sample. Panel on the top corresponds to the sample annotation: MFP – Gastric cancer type; Cohort - ID of dataset; Lauren - histological subtype of stomach cancer; Stage - tumor stage. FIG.5B provides a heatmap of pairwise gene group score correlation for the gene groups listed in Table 1. FIG.6 provides an example of gene group score comparison across GC TME types, according to some embodiments of the technology described herein. FIG.7 provides an exemplary heatmap of GC TME type gene group scores and cell deconvolution across the samples, according to some embodiments of the technology described herein. The columns represent single samples and are ordered by GC TME types as follows: B, E, D, A, and C. Panel on the top corresponds to the sample annotation: MFP – GC TME type; Lauren - histological subtype of stomach cancer; Cohort - ID of the dataset. The bottom bar plots represent the cell composition fraction. Each color corresponds to one cell type (according to the legend on the right). Boxplots on the right represent comparison of cell composition across GC TME types. FIGs.8A-8F provides representative data indicating The Cancer Genome Atlas (TCGA) histological data supports GC TME types identified by methods descried herein, according to some embodiments of the technology described herein. FIG.8A shows GC TME type A (Mesenchymal, EMT) and GC TME type C (Fibrotic) histology slides indicate a high level of fibrous tissue and fibroblasts, with lower purity (e.g., malignant cell content) (left); GC TME type B (Immune Enriched, Non-Fibrotic) histology slides indicate lymphoid infiltration and an inflamed immunophenotype (second from left); GC TME type D (Depleted) histology slides indicate high purity, a desert of immune cells, and low fibroblast numbers (second from right); GC TME type E (B cell Enriched) histology slides indicate lymphoid infiltration (inflamed immunophenotype) and the presence of germinative centers (e.g., high level of B cells, indicated by arrow) (right). FIG.8B shows a representation of the relative cell type content of GC TME type A. FIG.8C shows a representation of the relative cell type content of GC TME type B. FIG.8D shows a representation of the relative cell type content of GC TME type C. FIG.8E shows a representation of the relative cell type content of GC TME type D. FIG.8F shows a representation of the relative cell type content of GC TME type E. FIG.9 provides representative data indicating survival and progression analysis across different GC TME types, according to some embodiments of the technology described herein. FIG.10 provides representative data indicating that most microsatellite- unstable/instability (MSI) samples classified as having a GC TME type B, while most microsatellite stable (MSS) samples were classified as having GC TME type A, in accordance with some embodiments of the technology described herein. FIGs.11A and 11B provide representative data of mutation burden across GC TME types. FIG.11A shows that GC TME type B (i.e., immune enriched type) has the highest mutation load (ML). FIG.11B shows that the highest purity (e.g., cellularity or amount of malignant cells) (left) and the highest ploidy(right) were in the GC TME type D based on the ACGR cohort data. FIG.12 provides representative data indicating strong relationship between EBV positive status and GC TME B and E types (immune enriched and B cell enriched).1.0 = EBV positive; 00.0 = EBV negative. FIG.13 depicts an illustrative implementation of a computer system that may be used in connection with some embodiments of the technology described herein. FIG.14 provides representative data showing a comparison of Th17 T cell and Gamma- delta γδ T cell signatures across gastric cancer TME types, according to some embodiments of the technology described herein. FIG.15 provides a representative heatmap of gastric cancer samples classified into five distinct GC TME types (B, E, D, A, C) based on unsupervised dense clustering of GC TME signatures comprising 27 gene group scores, according to some embodiments of the technology described herein. Gene groups are listed vertically on the left side of the heatmap. Each column represents one sample. Panel on the top corresponds to the sample annotation: MFP – GCT; Cohort – ID of dataset; Lauren – histological subtype of gastric cancer; Stage – tumor stage. FIG.16 provides representative data for landscape of TCGA cohort and percentage of EBV positive/MSI samples across GC TME types, according to some embodiments of the technology described herein. Each column represents one sample. Top panels correspond to GC TME types; presence of EBV - black EBV positive and white EBV negative; MSI or MSS status; number of somatic mutations; GC TME type processes. On the right - the percentage of samples in each of the GC TME types with their characteristics of EBV and MSI/MSS. FIG.17 provides representative data for the number of responders classified as having each GC TME type, according to some embodiments of the technology described herein. PD = progressive disease, SD = stable disease, PR = partial response, CR = complete response. FIG.18 provides a representative heatmap of GC TME types, process signatures, and cell deconvolution across the samples in the Cabazi cohort, according to some embodiments of the technology described herein. The columns represent single samples and are ordered by the GC TME types as follows: B, E, D, A, and C. Panel on the top corresponds to the sample annotation: MFP - GC TME types; Recist (Response Evaluation Criteria in Solid Tumors) to cabazitaxel treatment; Lauren - histological subtype of gastric cancer; Biopsy location - location of the biopsy; Specific site - additional biopsy location for samples. The bottom bar plots represent the cell composition fraction. Each color corresponds to one cell type (according to the legend on the right). FIG.19 provides representative data for a comparison of GC TME types in pre- and post- treatment samples. On the left panel, similar to FIG.17, pre- and post- treated samples are presented. On the right panel - dynamic change in GC TME type in response to immunotherapeutic treatment was identified using principal component analysis. DETAILED DESCRIPTION Aspects of the disclosure relate to methods for characterizing subjects having certain cancers, for example gastric cancers. The disclosure is based, in part, on methods for determining the tumor microenvironment (TME) type of a subject’s gastric cancer. In some embodiments the methods comprise identifying a subject as having a particular gastric cancer (GC) TME type based upon a GC TME signature computed for the subject from their RNA expression data. The GC TME signature may comprise gene group scores for gene groups that are associated with malignant (e.g., cancer cells, tumor cells, etc.) and certain immune cells (e.g., T cells, B cells, Treg cells, etc.). The GC TME type identified for the subject may have various prognostic, diagnostic, and/or therapeutic applications. For example, in some embodiments, methods developed by the inventors and described herein are useful for identifying a subject’s prognosis based upon the GC TME type identified for the subject. In some embodiments, identification of a subject’s GC TME type is useful for determining whether or not a subject is likely to respond to treatment with certain therapeutic agents (e.g., an immunotherapy). Gastric cancer is the fifth most common cancer in the world. The 5-year overall survival rate for patients having advanced gastric cancer is approximately ~20%. Gastric cancer tumors are highly heterogeneous, and therefore present significant therapeutic challenges. Methods of classifying gastric cancer generally include histological analysis and/or molecular analysis. Histological analysis comprises obtaining a gastric tumor biopsy sample and analyzing the tissue to ascertain the presence of certain cell morphologies, and characterizing the tumor based upon the analysis. One example of histological analysis of gastric cancers is Lauren classification, for example as described by Ma et al. Oncol Lett.2016 May;11(5):2959-2964. Histological analysis is generally viewed as expensive, laborious, and time-consuming. Additionally, histological classification is typically reliant on a subjective determination made by laboratory personnel, and therefore the possibility of inaccuracy or human-error exists. Efforts have also been made to classify gastric cancers using molecular techniques in order to address the limitations of histological analysis. However, many studies performing molecular analysis on gastric cancer tissue are limited to gene expression studies of the tumor microenvironment, in particular tumor infiltrating immune lymphocytes (TIILs), for example as described by Hennequin et al.2016 OncoImmunology, 5:2, DOI: 10.1080/2162402X.2015.1054598. Accordingly, the inventors have recognized that there is a need to develop methods for molecular characterization of GC types specifically based upon the underlying biology of both the tumor microenvironment and malignant cells, rather than more broadly defined cancer biomarkers. Aspects of the disclosure relate to statistical techniques for analyzing expression data (e.g., RNA expression data), which was obtained from a biological sample obtained from a subject that has gastric cancer (GC), is suspected of having GC, or is at risk of developing GC, in order to generate a GC tumor microenvironment (TME) signature for the subject (termed an “GC TME signature” herein) and use this GC TME signature to identify a particular GC TME type that the subject may have. The inventors have recognized that a combination of certain gene group scores determined using RNA expression data of a subject (e.g., gene group scores for at least some of the gene groups listed in Table 1 and/or Table 2) may be combined to form a GC TME signature that characterizes patients having GC more accurately than previously developed methods. A GC TME signature comprising a combination of gene group scores from gene groups associated with the tumor microenvironment, and gene group scores from gene groups associated with malignant cells, in turn, may be used to identify the subject as having a particular gastric cancer (GC) tumor microenvironment (TME) type. In some embodiments, such GC TME types are useful for identifying the prognosis and/or likelihood that a subject will respond to particular therapeutic interventions (e.g., immunotherapy agents). The use of GC TME signatures comprising the combinations of gene group scores described by the disclosure represents an improvement over previously described GC molecular biomarkers or tumor microenvironment analyses because the specific groups of genes used to produce the GC TME signatures described herein better reflect the molecular tumor microenvironments of GC because these gene groups are associated with 1) gastric cancer tumor cells, and 2) gastric cancer tumor microenvironment. These focused combinations of gene groups (e.g., gene groups consisting of some or all of the gene group genes listed in Table 1 and/or Table 2) are unconventional, and differ from previously described molecular signatures, which attempt to incorporate expression data from either very large numbers of genes, or only account for TIILs. The GC TME typing methods described herein have several utilities. For example, identifying a subject’s GC TME type using methods described herein may allow for the subject to be diagnosed as having (or being at a high risk of developing) an aggressive form of GC at a timepoint that is not possible with previously described GC characterization methods. Earlier detection of aggressive GC types, enabled by the GC TME signatures described herein, improve the patient diagnostic technology by enabling earlier chemotherapeutic intervention for patients than currently possible for patients tested for GC using other methods (e.g., histological analysis). As described herein, the inventors have also determined that subjects identified by methods described herein as having GC TME type E are characterized as being more likely to have tertiary lymphoid structures (TLS), and thus having a good prognosis and/or an increased likelihood of responding to therapeutic treatments (e.g., an immunotherapy). As described further in the Examples, subjects identified as having GC TME type B were also observed to respond to immunotherapy treatment, regardless of MSS/MSI or EBV status. The recognition of GC TME type E and GC TME type B patients having a good prognosis or being likely to respond to immunotherapy independent of MSS/MSI or EBV status is surprising because it allows selection of certain patients (e.g., patients that are MSS or EBV-low) who have an increased likelihood of responding to the immunotherapy but typically are characterized as potential non-responders when other GC characterization methods are used (e.g., MSS/MSI or EBV status). Conversely, the inventors have determined that identifying a subject as having GC TME type A or GC TME type C using methods described herein, are less likely to have TLS and/or are more likely to have an aggressive or treatment-resistant form of GC. Thus, the techniques developed by the inventors and described herein improve patient treatment and associated outcomes by increasing patient comfort, and avoiding toxic side effects of chemotherapy that is not expected to be effective for the subject. In some aspects, the disclosure relates to producing one or more additional GC TME signatures (e.g., in addition to a GC TME signature comprising gene group scores from some or all of the gene groups listed in Table 1) of a gastric cancer subject. For example, the disclosure provides methods of determining a second GC TME signature comprising one or more gene group scores for some or all of the gene groups listed in Table 2. In some embodiments, a second GC signature comprises a Th17 signature gene group score and/or a Gamma-delta (γΔ) T cell gene group score. In some embodiments, a Th17 signature gene group score is produced using RNA expression levels for one or more (e.g., 1, 2, 3, 4, 5, 6, 7, or 8) genes in a TH17 gene group comprising the following genes: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, and CCR6. In some embodiments, a γΔ T cell gene group score is produced using RNA expression levels for one or more (e.g., 1, 2, 3, 4, or 5) genes in a γΔ T cell gene group comprising the following genes: TRGC1, TRDC, TRDV2, TRGC2, and TRDV1. Without wishing to be bound by any particular theory, calculation of a Th17 signature gene group score and/or a γΔ T cell gene group score may assist, either alone or in combination with a GC TME type in identifying dynamic changes of the GC TME types correlated with response to certain therapeutic agents. Gastric Cancers Aspects of the disclosure relate to methods of determining the gastric cancer (GC) TME type of a subject having, suspected of having, or at risk of having GC. As used herein, a subject may be a mammal, for example a human, non-human primate, rodent (e.g., rat, mouse, guinea pig, etc.), dog, cat, horse etc. In some embodiments, the subject is a human. The terms “individual” or “subject” may be used interchangeably with “patient.” As used herein, “gastric cancer” or “GC” or “stomach cancer” refers to any gastric or gastrointestinal cancer, for example, gastric or gastrointestinal adenocarcinoma, or any other type of malignancy caused by one or more various genetic mutations in the body that affect cells (originally present in or metastasized to) the stomach and/or intestine of a subject. As used herein, “cancer” refers to any malignant and/or invasive growth or tumor caused by abnormal cell growth in a subject, including solid tumors, blood cancer, bone marrow or lymphoid cancer, etc. Examples of gastric cancers include but are not limited to esophageal cancer, stomach cancer, liver cancer, pancreatic cancer, colorectal cancer, anal cancer, adenocarcinoma of the stomach, fungating (polypoid) stomach cancer, ulcerating stomach cancer, superficial spreading stomach cancer, diffusely spreading stomach cancer, malignant lymphoma of the stomach, liposarcoma, fibrosarcoma, carcinosarcoma, and gastrointestinal stromal tumor (GST). A subject having GC may exhibit one or more signs or symptoms of GC, for example the presence of cancerous cells (e.g., tumor cells), fever, swelling, bleeding, nausea and vomiting, heartburn, and weight loss. In some embodiments, a subject having GC does not exhibit one or more signs or symptoms of GC. In some embodiments, a subject having GC has been diagnosed by a medical professional (e.g., a licensed physician) as having GC based upon one or more assays (e.g., clinical assays, molecular diagnostics, etc.) that indicate that the subject has GC, even in the absence of one or more signs or symptoms. A subject suspected of having GC typically exhibits one or more signs or symptoms of GC. In some embodiments, a subject suspected of having GC exhibits one or more signs or symptoms of GC but has not been diagnosed by a medical professional (e.g., a licensed physician) and/or has not received a test result (e.g., a clinical assay, molecular diagnostic, etc.) indicating that the subject has GC. A subject at risk of having GC may or may not exhibit one or more signs or symptoms of GC. In some embodiments, a subject at risk of having GC comprises one or more risk factors that increase the likelihood that the subject will develop GC. Examples of risk factors include the presence of pre-cancerous cells in a clinical sample, having one or more genetic mutations that predispose the subject to developing cancer (e.g., GC), taking one or more medications that increase the likelihood that the subject will develop cancer (e.g., GC), family history of GC, and the like. FIG.1 is a flowchart of an illustrative process 100 for determining a GC TME signature for a subject, using the determined GC TME signature to identify the GC TME type for the subject, and using the GC TME type of the subject to identify whether or not the subject is likely to respond to immunotherapy. Various (e.g., some or all) acts of process 100 may be implemented using any suitable computing device(s). For example, in some embodiments, one or more acts of the illustrative process 100 may be implemented in a clinical or laboratory setting. For example, one or more acts of the process 100 may be implemented on a computing device that is located within the clinical or laboratory setting. In some embodiments, the computing device may directly obtain RNA expression data from a sequencing apparatus located within the clinical or laboratory setting. For example, a computing device included in the sequencing apparatus may directly obtain the RNA expression data from the sequencing apparatus. In some embodiments, the computing device may indirectly obtain RNA expression data from a sequencing apparatus that is located within or external to the clinical or laboratory setting. For example, a computing device that is located within the clinical or laboratory setting may obtain expression data via a communication network, such as Internet or any other suitable network, as aspects of the technology described herein are not limited to any particular communication network. Additionally or alternatively, one or more acts of the illustrative process 100 may be implemented in a setting that is remote from a clinical or laboratory setting. For example, the one or more acts of process 100 may be implemented on a computing device that is located externally from a clinical or laboratory setting. In this case, the computing device may indirectly obtain RNA expression data that is generated using a sequencing apparatus located within or external to a clinical or laboratory setting. For example, the expression data may be provided to computing device via a communication network, such as Internet or any other suitable network. It should be appreciated that, in some embodiments, not all acts of process 100, as illustrated in FIG.1, may be implemented using one or more computing devices. For example, the act 120 of administering one or more immunotherapies to the subject may be implemented manually (e.g., by a clinician) Process 100 begins at act 102 where sequencing data for a subject is obtained. In some embodiments, the sequencing data may be obtained by sequencing a biological sample (e.g., stomach biopsy and/or tumor tissue) obtained from the subject using any suitable sequencing technique. The sequencing data may include sequencing data of any suitable type, from any suitable source, and be in any suitable format. Examples of sequencing data, sources of sequencing data, and formats of sequencing data are described herein including in the section called “Obtaining RNA Expression Data”. As one illustrative example, in some embodiments, the sequencing data may comprise bulk sequencing data. The bulk sequencing data may comprise at least 1 million reads, at least 5 million reads, at least 10 million reads, at least 20 million reads, at least 50 million reads, or at least 100 million reads. In some embodiments, the sequencing data comprises bulk RNA sequencing (RNA-seq) data, single cell RNA sequencing (scRNA-seq) data, or next generation sequencing (NGS) data. In some embodiments, the sequencing data comprises microarray data. Next, process 100 proceeds to act 104, where the sequencing data obtained at act 102 is processed to obtain RNA expression data. This may be done in any suitable way and may involve normalizing bulk sequencing data to transcripts-per-million (TPM) units (or other units) and/or log transforming the RNA expression levels in TPM units. Converting the data to TPM units and normalization are described herein including with reference to FIG.2. Next, process 100 proceeds to act 106, where a gastric cancer (GC) tumor microenvironment (TME) signature is generated for the subject using the RNA expression data generated at act 104 (e.g., from bulk-sequencing data, converted to TPM units and subsequently log-normalized, as described herein including with reference to FIG.2). As described herein, in some embodiments, a GC TME signature comprises two or more (e.g., 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, etc.) gene group scores. In some embodiments, the two or more gene group scores comprise gene group scores (which may also be referred to as gene group enrichment scores or gene group expression scores) for some or all of the gene groups shown in Table 1. Accordingly, act 106 comprises: act 108 where the gene group scores are determined, act 110 where the GC TME signature is determined, and act 112 where the GC TME type is determined by using GC TME signature. In some embodiments, determining the gene group scores comprises determining, for each of multiple (e.g., some or all of the) gene groups listed in Table 1, a respective gene group score. In some embodiments, determining the gene group scores comprises determining respective gene group scores for 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20 gene groups (e.g., gene groups listed in Table 1). In some embodiments, determining the gene group scores further comprises determining respective gene group scores for 1, 2, 3, 4, 5, 6, or 7 gene groups listed in Table 2. The gene group score for a particular gene group may be determined using RNA expression levels for at least some of the genes in the gene group (e.g., the RNA expression levels obtained at act 104). The RNA expression levels may be processed using a gene set enrichment analysis (GSEA) technique to determine the score for the particular gene group. For example, in some embodiments, determining the GC TME signature comprises: determining gene group scores using the RNA expression levels for at least three genes from each of at least two of the gene groups, the gene groups including: NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining the GC TME gene signature comprises: determining gene group scores using the RNA expression levels for at least three genes from each of at least two of the gene groups, the two gene groups including: (a) MHC I group: HLA-A, HLA-B, HLA- C, B2M, TAP1, TAP2, NLRC5, TAPBP; (b) MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; and (c) Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; (d) Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; (e) T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; (f) NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; (g) T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; (h) B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; (i) M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; (j) Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; (k) Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; (l) Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; (m) Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; (n) MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; (o) M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; (p) Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; (q) Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; (r) Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and (s) Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. Aspects of determining the gene group scores are described herein, including with reference to FIG.3 and in the Section titled “Gene Expression Signatures”. As described above, at act 110, the GC TME signature is generated. In some embodiments, the GC TME signature consists of only gene group scores for one or more (e.g., all) of the gene groups listed in Table 1. In some embodiments, the GC TME signature comprises gene group scores for at least 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20 gene groups listed in Table 1. In some embodiments, each gene group score for a particular gene group is determined using RNA expression levels of some or all (e.g., at least 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, etc.) of the genes of each gene group listed in Table 1. In other embodiments, the GC TME signature includes one or more other gene group scores in addition to the gene group scores listed in Table 1. In some embodiments, the GC TME signature further comprises gene group scores from at least 2, 3, 4, 5, 6, or 7, gene groups listed in Table 2. In some embodiments each gene group score is determined using RNA expression levels of some or all (e.g., at least 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, etc.) of the genes of each gene group listed in Table 2. In some embodiments, act 110 comprises generating a GC TME signature using only gene group scores from some or all of the gene groups listed in Table 1. In some embodiments, act 110 further comprises a step of generating a second GC TME signature comprising gene group scores for some or all of the gene groups listed in Table 2. In some embodiments, the second GC TME signature comprises a Th17 signature gene group score (e.g., a gene group score generated using RNA expression levels of some or all of the genes of the Th17 signature gene group listed in Table 2) and/or a γΔ T cell gene group score (e.g., a gene group score generated using RNA expression levels of some or all of the genes of the γΔ T cell gene group listed in Table 2). In some embodiments, a second GC TME signature is not used at act 112 to identify a GC TME type of a subject. In some embodiments, a second GC TME signature is indicative of changes in a subject’s GC TME type after treatment with an immunotherapy. Next, process 100 proceeds to act 112, where a GC TME type is identified for the subject using the GC TME signature generated at act 110. This may be done in any suitable way. For example, in some embodiments, the each of the possible GC TME types is associated with a respective plurality of GC TME signature clusters. In such embodiments, a GC TME type for the subject may be identified by associating the GC TME signature of the subject with a particular one of the plurality of GC TME signature clusters; and identifying the GC TME type for the subject as the GC TME type corresponding to the particular one of the plurality of GC TME signature clusters to which the GC TME signature of the subject is associated. Examples of GC TME types are described herein. Aspects of identifying a GC TME type for a subject are described herein including in the section below titled “Generating GC TME Signature and Identifying TME Type”. As described above, a subject’s GC TME type is identified at act 112. In some embodiments, the GC TME type of a subject is identified to be one of the following GC TME types: GC TME type A, GC TME type B, GC TME type C, GC TME type D, or GC TME type E. Next, process 100 proceeds to act 118, where the subject’s likelihood of responding to an immunotherapy is identified using the GC TME type identified at act 112. In some embodiments, when a subject is identified as having a GC TME type E or GC TME type B at act 112, the subject is identified as having an increased likelihood of responding to an immunotherapy (e.g., a PD1 antibody, such as pembrolizumab) relative to a subject having other GC TME types, at act 118. Aspects of identifying whether or not a subject is likely to respond to an immunotherapy are described herein including in the section below titled “Therapeutic Indications”. In some embodiments, process 100 completes after act 118 completes. In some such embodiments, the determined GC TME signature and/or identified GC TME type, and/or the identified likelihood the subject will respond to an immunotherapy may be stored for subsequent use, provided to one or more recipients (e.g., a clinician, a researcher, etc.), and/or used to update the GC TME signature clusters (as described hereinbelow). However, in some embodiments, one or more other acts are performed after act 112.. For example, in the illustrated embodiment of FIG.1, process 100 may include one or more of optional acts 114, 116, and 120 shown using dashed lines in FIG.1. For example, at act 114, a subject may be identified as having tertiary lymphoid structures (TLS) based on the GC TME type determined for the subject. For example, a subject may be identified (at act 114) as having a TLS when the subject is identified as having GC TME type E. In another example, at act 116, a prognosis may be identified for the subject. In some embodiments, when a subject is identified as having GC TME type E at act 112 and/or identified as having TLS at act 114, the subject is then identified as having an increased likelihood of having a good prognosis (e.g., as measured by overall survival (OS) or progression-free survival (PFS) at act 116. In another example, when a subject is identified as having GC TME type B or GC TME type E at act 112, and/or identified as having an increased likelihood of responding to immunotherapy at act 118, the subject is administered one or more immunotherapies at act 120. Examples of immunotherapies are provided herein. It should be appreciated that although acts 114, 116, and 120 are indicated as optional in the example of FIG.1, in other embodiments, one or more other acts may be optional (in addition to or instead of acts 114, 116, and 120). For example, in some embodiments, acts 102 and 104 may be optional (e.g., when the sequencing data is obtained and processed to obtain RNA expression data previously, process 100 may begin at act 106 by accessing the previously obtained RNA expression data). In some embodiments, the process 100 may comprise acts 102, 104, 106, 118 and 120, without acts 114 and 116. In some embodiments, the process 100 may comprise acts 102, 104, 106, 116, 114, 118, and 120 without act 116. Biological Samples Aspects of the disclosure relate to methods for determining a GC TME type of a subject by obtaining sequencing data from a biological sample that has been obtained from the subject. The biological sample may be from any source in the subject’s body including, but not limited to, any fluid such as blood (e.g., whole blood, blood serum, or blood plasma), lymph node, stomach, small intestine. Other source in the subject’s body may be from saliva, tears, synovial fluid, cerebrospinal fluid, pleural fluid, pericardial fluid, ascitic fluid, and/or urine, hair, skin (including portions of the epidermis, dermis, and/or hypodermis), oropharynx, laryngopharynx, esophagus, bronchus, salivary gland, tongue, oral cavity, nasal cavity, vaginal cavity, anal cavity, bone, bone marrow, brain, thymus, spleen, appendix, colon, rectum, anus, liver, biliary tract, pancreas, kidney, ureter, bladder, urethra, uterus, vagina, vulva, ovary, cervix, scrotum, penis, prostate, testicle, seminal vesicles, and/or any type of tissue (e.g., muscle tissue, epithelial tissue, connective tissue, or nervous tissue). The biological sample may be any type of sample including, for example, a sample of a bodily fluid, one or more cells, one or more pieces of tissue(s) or organ(s). In some embodiments, the biological sample comprises gastrointestinal tissue sample of the subject. Examples of gastrointestinal tissue samples include but are not limited to mucosal tissue, submucosal tissue, muscular layer tissue, and serous layer tissue (also referred to as serosa tissue). In some embodiments, a gastrointestinal tissue sample comprises one or more cell types derived from a stomach (e.g., mucous cells, parietal cells, chief cells, endocrine cells, etc.). In some embodiments, a gastrointestinal tissue sample comprises one or more cell types derived from gastrointestinal tissue, for example enterocytes, Paneth cells, goblet cells, neuroendocrine cells, etc. In some embodiments, a gastrointestinal tissue sample may be obtained from a subject using a surgical procedure (e.g., laparoscopic surgery, microscopically controlled surgery, or endoscopy), bone marrow biopsy, punch biopsy, endoscopic biopsy, or needle biopsy (e.g., a fine-needle aspiration, core needle biopsy, vacuum-assisted biopsy, or image-guided biopsy). A sample of lymph node or blood, in some embodiments, refers to a sample comprising cells, e.g., cells from a blood sample or lymph node sample. In some embodiments, the sample comprises non-cancerous cells. In some embodiments, the sample comprises pre-cancerous cells. In some embodiments, the sample comprises cancerous cells. In some embodiments, the sample comprises blood cells. In some embodiments, the sample comprises lymph node cells. In some embodiments, the sample comprises lymph node cells and blood cells. A sample of blood may be a sample of whole blood or a sample of fractionated blood. In some embodiments, the sample of blood comprises whole blood. In some embodiments, the sample of blood comprises fractionated blood. In some embodiments, the sample of blood comprises buffy coat. In some embodiments, the sample of blood comprises serum. In some embodiments, the sample of blood comprises plasma. In some embodiments, the sample of blood comprises a blood clot. In some embodiments, a sample of blood is collected to obtain the cell-free nucleic acid (e.g., cell-free DNA) in the blood. In some embodiments, the sample may be from a cancerous tissue or an organ or a tissue or organ suspected of having one or more cancerous cells. In some embodiments, the sample may be from a healthy (e.g., non-cancerous) tissue or organ. In some embodiments, a sample from a subject (e.g., a biopsy from a subject) may include both healthy and cancerous cells and/or tissue. In certain embodiments, one sample will be taken from a subject for analysis. In some embodiments, more than one (e.g., 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, or more) samples may be taken from a subject for analysis. In some embodiments, one sample from a subject will be analyzed. In certain embodiments, more than one (e.g., 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, or more) samples may be analyzed. If more than one sample from a subject is analyzed, the samples may be procured at the same time (e.g., more than one sample may be taken in the same procedure), or the samples may be taken at different times (e.g., during a different procedure including a procedure 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 days; 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 weeks; 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 months, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 years, or 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 decades after a first procedure). A second or subsequent sample may be taken or obtained from the same region (e.g., from the same tumor or area of tissue) or a different region (including, e.g., a different tumor). A second or subsequent sample may be taken or obtained from the subject after one or more treatments, and may be taken from the same region or a different region. As a non-limiting example, the second or subsequent sample may be useful in determining whether the cancer in each sample has different characteristics (e.g., in the case of samples taken from two physically separate tumors in a patient) or whether the cancer has responded to one or more treatments (e.g., in the case of two or more samples from the same tumor prior to and subsequent to a treatment). Any of the biological samples described herein may be obtained from the subject using any known technique. See, for example, the following publications on collecting, processing, and storing biological samples, each of which is incorporated by reference herein in its entirety: Biospecimens and biorepositories: from afterthought to science by Vaught et al. (Cancer Epidemiol Biomarkers Prev.2012 Feb;21(2):253-5), and Biological sample collection, processing, storage and information management by Vaught and Henderson (IARC Sci Publ. 2011;(163):23-42). Any of the biological samples from a subject described herein may be stored using any method that preserves stability of the biological sample. In some embodiments, preserving the stability of the biological sample means inhibiting components (e.g., DNA, RNA, protein, or tissue structure or morphology) of the biological sample from degrading until they are measured so that when measured, the measurements represent the state of the sample at the time of obtaining it from the subject. In some embodiments, a biological sample is stored in a composition that is able to penetrate the same and protect components (e.g., DNA, RNA, protein, or tissue structure or morphology) of the biological sample from degrading. As used herein, degradation is the transformation of a component from one form to another form such that the first form is no longer detected at the same level as before degradation. In some embodiments, the biological sample is stored using cryopreservation. Non- limiting examples of cryopreservation include, but are not limited to, step-down freezing, blast freezing, direct plunge freezing, snap freezing, slow freezing using a programmable freezer, and vitrification. In some embodiments, the biological sample is stored using lyophilisation. In some embodiments, a biological sample is placed into a container that already contains a preservant (e.g., RNALater to preserve RNA) and then frozen (e.g., by snap-freezing), after the collection of the biological sample from the subject. In some embodiments, such storage in frozen state is done immediately after collection of the biological sample. In some embodiments, a biological sample may be kept at either room temperature or 4 o C for some time (e.g., up to an hour, up to 8 h, or up to 1 day, or a few days) in a preservant or in a buffer without a preservant, before being frozen. Non-limiting examples of preservants include formalin solutions, formaldehyde solutions, RNALater or other equivalent solutions, TriZol or other equivalent solutions, DNA/RNA Shield or equivalent solutions, EDTA (e.g., Buffer AE (10 mM Tris·Cl; 0.5 mM EDTA, pH 9.0)) and other coagulants, and Acids Citrate Dextronse (e.g., for blood specimens). In some embodiments, special containers may be used for collecting and/or storing a biological sample. For example, a vacutainer may be used to store blood. In some embodiments, a vacutainer may comprise a preservant (e.g., a coagulant, or an anticoagulant). In some embodiments, a container in which a biological sample is preserved may be contained in a secondary container, for the purpose of better preservation, or for the purpose of avoid contamination. Any of the biological samples from a subject described herein may be stored under any condition that preserves stability of the biological sample. In some embodiments, the biological sample is stored at a temperature that preserves stability of the biological sample. In some embodiments, the sample is stored at room temperature (e.g., 25 °C). In some embodiments, the sample is stored under refrigeration (e.g., 4 °C). In some embodiments, the sample is stored under freezing conditions (e.g., -20 °C). In some embodiments, the sample is stored under ultralow temperature conditions (e.g., -50 °C to -800 °C). In some embodiments, the sample is stored under liquid nitrogen (e.g., -1700 °C). In some embodiments, a biological sample is stored at -60 o C to -8- o C (e.g., -70 o C) for up to 5 years (e.g., up to 1 month, up to 2 months, up to 3 months, up to 4 months, up to 5 months, up to 6 months, up to 7 months, up to 8 months, up to 9 months, up to 10 months, up to 11 months, up to 1 year, up to 2 years, up to 3 years, up to 4 years, or up to 5 years). In some embodiments, a biological sample is stored as described by any of the methods described herein for up to 20 years (e.g., up to 5 years, up to 10 years, up to 15 years, or up to 20 years). Obtaining RNA Expression Data Aspects of the disclosure relate to methods of determining a GC TME type of a subject using sequencing data or RNA expression data obtained from a biological sample from the subject. The RNA expression data used in methods described herein typically is derived from sequencing data obtained from the biological sample. The sequencing data may be obtained from the biological sample using any suitable sequencing technique and/or apparatus. In some embodiments, the sequencing apparatus used to sequence the biological sample may be selected from any suitable sequencing apparatus known in the art including, but not limited to, Illumina TM , SOLid TM , Ion Torrent TM , PacBio TM , a nanopore-based sequencing apparatus, a Sanger sequencing apparatus, or a 454 TM sequencing apparatus. In some embodiments, sequencing apparatus used to sequence the biological sample is an Illumina sequencing (e.g., NovaSeq TM , NextSeq TM , HiSeq TM , MiSeq TM , or MiniSeq TM ) apparatus. After the sequencing data is obtained, it is processed in order to obtain the RNA expression data. RNA expression data may be acquired using any method known in the art including, but not limited to whole transcriptome sequencing, whole exome sequencing, total RNA sequencing, mRNA sequencing, targeted RNA sequencing, RNA exome capture sequencing, next generation sequencing, and/or deep RNA sequencing. In some embodiments, RNA expression data may be obtained using a microarray assay. In some embodiments, the sequencing data is processed to produce RNA expression data. In some embodiments, RNA sequence data is processed by one or more bioinformatics methods or software tools, for example RNA sequence quantification tools (e.g., Kallisto) and genome annotation tools (e.g., Gencode v23), in order to produce expression data. The Kallisto software is described in Nicolas L Bray, Harold Pimentel, Páll Melsted and Lior Pachter, Near- optimal probabilistic RNA-seq quantification, Nature Biotechnology 34, 525–527 (2016), doi:10.1038/nbt.3519, which is incorporated by reference in its entirety herein. In some embodiments, microarray expression data is processed using a bioinformatics R package, such as “affy” or “limma”, in order to produce expression data. The “affy” software is described in Bioinformatics.2004 Feb 12;20(3):307-15. doi: 10.1093/bioinformatics/btg405. “affy--analysis of Affymetrix GeneChip data at the probe level” by Laurent Gautier 1, Leslie Cope, Benjamin M Bolstad, Rafael A Irizarry PMID: 14960456 DOI: 10.1093/bioinformatics/btg405, which is incorporated by reference herein in its entirety. The “limma” software is described in Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK "limma powers differential expression analyses for RNA-sequencing and microarray studies." Nucleic Acids Res.2015 Apr 20;43(7):e47.20. https://doi.org/10.1093/nar/gkv007 PMID: 25605792, PMCID: PMC4402510, which is incorporated by reference herein its entirety. In some embodiments, sequencing data and/or expression data comprises more than 5 kilobases (kb). In some embodiments, the size of the obtained RNA data is at least 10 kb. In some embodiments, the size of the obtained RNA sequencing data is at least 100 kb. In some embodiments, the size of the obtained RNA sequencing data is at least 500 kb. In some embodiments, the size of the obtained RNA sequencing data is at least 1 megabase (Mb). In some embodiments, the size of the obtained RNA sequencing data is at least 10 Mb. In some embodiments, the size of the obtained RNA sequencing data is at least 100 Mb. In some embodiments, the size of the obtained RNA sequencing data is at least 500 Mb. In some embodiments, the size of the obtained RNA sequencing data is at least 1 gigabase (Gb). In some embodiments, the size of the obtained RNA sequencing data is at least 10 Gb. In some embodiments, the size of the obtained RNA sequencing data is at least 100 Gb. In some embodiments, the size of the obtained RNA sequencing data is at least 500 Gb. In some embodiments, the expression data is acquired through bulk RNA sequencing. Bulk RNA sequencing may include obtaining expression levels for each gene across RNA extracted from a large population of input cells (e.g., a mixture of different cell types.) In some embodiments, the expression data is acquired through single cell sequencing (e.g., scRNA-seq). Single cell sequencing may include sequencing individual cells. In some embodiments, bulk sequencing data comprises at least 1 million reads, at least 5 million reads, at least 10 million reads, at least 20 million reads, at least 50 million reads, or at least 100 million reads. In some embodiments, bulk sequencing data comprises between 1 million reads and 5 million reads, 3 million reads and 10 million reads, 5 million reads and 20 million reads, 10 million reads and 50 million reads, 30 million reads and 100 million reads, or 1 million reads and 100 million reads (or any number of reads including, and between). In some embodiments, the expression data comprises next-generation sequencing (NGS) data. In some embodiments, the expression data comprises microarray data. Expression data (e.g., indicating expression levels) for a plurality of genes may be used for any of the methods or compositions described herein. The number of genes which may be examined may be up to and inclusive of all the genes of the subject. In some embodiments, expression levels may be determined for all of the genes of a subject. As a non-limiting example, four or more, five or more, six or more, seven or more, eight or more, nine or more, ten or more, eleven or more, twelve or more, 13 or more, 14 or more, 15 or more, 16 or more, 17 or more, 18 or more, 19 or more, 20 or more, 21 or more, 22 or more, 23 or more, 24 or more, 25 or more, 26 or more, 27 or more, 28 or more, 29 or more, 30 or more, 35 or more, 40 or more, 50 or more, 60 or more, 70 or more, 80 or more, 90 or more, 100 or more, 125 or more, 150 or more, 175 or more, 200 or more, 225 or more, 250 or more, 275 or more, or 300 or more genes may be used for any evaluation described herein. As another set of non-limiting examples, the expression data may include, for each gene group listed in Table 1, expression data for at least 5, at least 10, at least 15, at least 20, at least 25, at least 35, at least 50, at least 75, at least 100 genes selected from each gene group. In some embodiments, RNA expression data is obtained by accessing the RNA expression data from at least one computer storage medium on which the RNA expression data is stored. Additionally or alternatively, in some embodiments, RNA expression data may be received from one or more sources via a communication network of any suitable type. For example, in some embodiment, the RNA expression data may be received from a server (e.g., a SFTP server, or Illumina BaseSpace). The RNA expression data obtained may be in any suitable format, as aspects of the technology described herein are not limited in this respect. For example, in some embodiments, the RNA expression data may be obtained in a text-based file (e.g., in a FASTQ, FASTA, BAM, or SAM format). In some embodiments, a file in which sequencing data is stored may contains quality scores of the sequencing data. In some embodiments, a file in which sequencing data is stored may contain sequence identifier information. Expression data, in some embodiments, includes gene expression levels. Gene expression levels may be detected by detecting a product of gene expression such as mRNA and/or protein. In some embodiments, gene expression levels are determined by detecting a level of a mRNA in a sample. As used herein, the terms “determining” or “detecting” may include assessing the presence, absence, quantity and/or amount (which can be an effective amount) of a substance within a sample, including the derivation of qualitative or quantitative concentration levels of such substances, or otherwise evaluating the values and/or categorization of such substances in a sample from a subject. FIG. 2 shows an exemplary process 104 for processing sequencing data to obtain RNA expression data from sequencing data. Process 104 may be performed by any suitable computing device or devices, as aspects of the technology described herein are not limited in this respect. For example, process 104 may be performed by a computing device part of a sequencing apparatus. In other embodiments, process 104 may be performed by one or more computing devices external to the sequencing apparatus. Process 104 begins at act 200, where sequencing data is obtained from a biological sample obtained from a subject. The sequencing data is obtained by any suitable method, for example, using any of the methods described herein including in the Section titled “Biological Samples”. In some embodiments, the sequencing data obtained at act 104 comprises RNA-seq data. In some embodiments, the biological sample comprises blood or tissue. In some embodiments, the biological sample comprises one or more tumor cells, for example, one or more GC tumor cells. Next, process 104 proceeds to act 202 where the sequencing data obtained at act 200 is normalized to transcripts per kilobase million (TPM) units. The normalization may be performed using any suitable software and in any suitable way. For example, in some embodiments, TPM normalization may be performed according to the techniques described in Wagner et al. (Theory Biosci. (2012) 131:281–285), which is incorporated by reference herein in its entirety. In some embodiments, the TPM normalization may be performed using a software package, such as, for example, the gcrma package. Aspects of the gcrma package are described in Wu J, Gentry RIwcfJMJ (2021). “gcrma: Background Adjustment Using Sequence Information. R package version 2.66.0.”, which is incorporated by reference in its entirety herein. In some embodiments, RNA expression level in TPM units for a particular gene may be calculated according to the following formula: Next, process 104 proceeds to act 204, where the RNA expression levels in TPM units (as determined at act 202) may be log transformed. Process 104 is illustrative and there are variations. For example, in some embodiments, one or both of acts 202 and 204 may be omitted. Thus, in some embodiments, the RNA expression levels may not be normalized to transcripts per million units and may, instead, be converted to another type of unit (e.g., reads per kilobase million (RPKM) or fragments per kilobase million (FPKM) or any other suitable unit). Additionally or alternatively, in some embodiments, the log transformation may be omitted. Instead, no transformation may be applied in some embodiments, or one or more other transformations may be applied in lieu of the log transformation. RNA expression data obtained by process 104 can include the sequence data generated by a sequencing protocol (e.g., the series of nucleotides in a nucleic acid molecule identified by next-generation sequencing, sanger sequencing, etc.) as well as information contained therein (e.g., information indicative of source, tissue type, etc.) which may also be considered information that can be inferred or determined from the sequence data. In some embodiments, expression data obtained by process 104 can include information included in a FASTA file, a description and/or quality scores included in a FASTQ file, an aligned position included in a BAM file, and/or any other suitable information obtained from any suitable file. Gene Expression Signatures Aspects of the disclosure relate to processing of expression data to determine one or more gene expression signatures (e.g., a GC TME signature). In some embodiments, expression data (e.g., RNA expression data) is processed using a computing device to determine the one or more gene expression signatures. In some embodiments, the computing device may be operated by a user such as a doctor, clinician, researcher, patient, or other individual. For example, the user may provide the expression data as input to the computing device (e.g., by uploading a file), and/or may provide user input specifying processing or other methods to be performed using the expression data. In some embodiments, expression data may be processed by one or more software programs running on computing device. In some embodiments, methods described herein comprise an act of determining a GC TME signature comprising gene group scores for respective gene groups in a plurality of gene groups. In some embodiments, a GC TME signature comprises gene group scores for at least one (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20) of the gene groups listed in Table 1. In some embodiments, a GC TME signature comprises gene group scores for at least one (e.g., 1, 2, 3, 4, 5, 6, or 7,) of the gene groups listed in Table 2. In some embodiments, a GC TME signature comprises gene group scores for at least one (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20) of the gene groups listed in Table 1 and gene group scores for at least one (e.g., 1, 2, 3, 4, 5, 6, or 7,) of the gene groups listed in Table 2. The number of genes in a gene group used to determine a gene group score may vary. In some embodiments, all RNA expression levels for all genes in a particular gene group may be used to determine a gene group score for the particular gene group. In other embodiments, RNA expression data for fewer than all genes may be used (e.g., RNA expression levels for at least two genes, at least three genes, at least five genes, between 2 and 10 genes, between 5 and 15 genes, between 3 and 30 genes, or any other suitable range within these ranges). In some embodiments, a GC TME signature comprises a gene group score for the MHC I gene group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, or at least seven genes) in the MHC I gene group, which is defined by its constituent genes: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP. In some embodiments, a GC TME signature comprises a gene group score for the MHC II gene group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes) in the MHC II gene group, which is defined by its constituent gene: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA- DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA. In some embodiments, a GC TME signature comprises a gene group score for the Coactivation molecules group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, or at least ten genes) in the Coactivation molecules group, which is defined by its constituent genes: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70. In some embodiments, a GC TME signature comprises a gene group score for the Effector cells group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than ten genes) in the Effector cells group, which is defined by its constituent genes: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B. In some embodiments, a GC TME signature comprises a gene group score for the T cell traffic group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, or at least eight genes) in the T cell traffic group, which is defined by its constituent genes: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6. In some embodiments, a GC TME signature comprises a gene group score for the NK cells group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than ten genes) in the NK cells group, which is defined by its constituent genes: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3. In some embodiments, a GC TME signature comprises a gene group score for the T cells group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, or at least ten genes) in the T cells group, which is defined by its constituent genes: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, PLK1, CCNB1, MCM2, and MCM6. In some embodiments, a GC TME signature comprises a gene group score for the B cells group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than ten genes) in the B cells group, which is defined by its constituent genes: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1. In some embodiments, a GC TME signature comprises a gene group score for the M1 signature group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, or at least eight genes) in the M1 signature group, which is defined by its constituent genes: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A. In some embodiments, a GC TME signature comprises a gene group score for the Antitumor cytokines group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, or at least five genes) in the Antitumor cytokines group, which is defined by its constituent genes: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21. In some embodiments, a GC TME signature comprises a gene group score for the Checkpoint inhibition group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than ten genes) in the Checkpoint inhibition group, which is defined by its constituent genes: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR. In some embodiments, a GC TME signature comprises a gene group score for the Treg group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, or at least six genes) in the Treg group, which is defined by its constituent genes: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2. In some embodiments, a GC TME signature comprises a gene group score for the Neutrophil signature group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, or at least nine genes) in the Neutrophil group, which is defined by its constituent genes: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1. In some embodiments, a GC TME signature comprises a gene group score for the MDSC group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, or at least six genes) in the MDSC group, which is defined by its constituent genes: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6. In some embodiments, a GC TME signature comprises a gene group score for the M2 signature group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, or at least seven genes) in the M2 group, which is defined by its constituent genes: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68. In some embodiments, a GC TME signature comprises a gene group score for the Cancer associated fibroblast (CAF) group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than 10 genes) in the Cancer associated fibroblast (CAF) group, which is defined by its constituent genes: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX. In some embodiments, a GC TME signature comprises a gene group score for the angiogenesis group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than ten genes) in the angiogenesis group, which is defined by its constituent genes: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5. In some embodiments, a GC TME signature comprises a gene group score for the Proliferation rate group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, or at least six genes) in the Proliferation rate group, which is defined by its constituent genes: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6. In some embodiments, a GC TME signature comprises a gene group score for the Lgr5 ISC group. In some embodiments, this gene group score may be calculated using RNA expression levels of at least three genes (e.g., at least three genes, at least four genes, at least five genes, at least six genes, at least seven genes, at least eight genes, at least nine genes, at least ten genes, or more than 10 genes) in the Lgr5 ISC group, which is defined by its constituent genes: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. In some embodiments, determining a GC TME signature comprises determining a respective gene group score for each of at least two of the following gene groups, using, for a particular gene group, RNA expression levels for at least three genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: MHC I group: HLA-A, HLA-B, HLA-C, B2M, TAP1, TAP2, NLRC5, TAPBP; MHC II group: HLA-DRA, HLA-DRB1, HLA-DMA, HLA-DPA1, HLA-DPB1, HLA-DMB, HLA-DQB1, HLA-DQA1, CIITA; Coactivation molecules group: CD28, CD40, TNFRSF4, ICOS, TNFRSF9, CD27, CD80, CD86, CD40LG, CD83, TNFSF4, ICOSLG, TNFSF9, CD70; Effector cells group: IFNG, GZMA, GZMB, PRF1, GZMK, ZAP70, GNLY, FASLG, TBX21, EOMES, CD8A, CD8B; T cell traffic group: CXCL9, CXCL10, CXCL11, CX3CL1, CCL3, CCL4, CX3CR1, CXCL16, CXCR6; NK cells group: NKG7, CD160, CD244, NCR1, KLRC2, KLRK1, CD226, GZMH, GNLY, IFNG, KIR2DL4, EOMES, GZMB, FGFBP2, KLRF1, SH2D1B, NCR3; T cells group: TBX21, ITK, CD3D, CD3E, CD3G, TRAC, TRBC1, TRBC2, CD28, CD5, TRAT1; B cells group: CD19, MS4A1, TNFRSF13C, CR2, TNFRSF17, TNFRSF13B, CD22, CD79A, CD79B, BLK, FCRL5, PAX5, STAP1; M1 signature group: NOS2, TNF, IL1B, SOCS3, CMKLR1, IRF5, IL12A, IL12B, IL23A; Antitumor cytokines group: TNF, IFNB1, IFNA2, CCL3, TNFSF10, IL21; Checkpoint inhibition group: PDCD1, CD274, CTLA4, LAG3, PDCD1LG2, BTLA, HAVCR2, TIGIT, VSIR; Treg group: FOXP3, CTLA4, IL10, TNFRSF18, CCR8, IKZF4, IKZF2; Neutrophil signature group: MPO, ELANE, PRTN3, CTSG, CXCR1, CXCR2, FCGR3B, CD177, FFAR2, PGLYRP1; MDSC group: IDO1, ARG1, IL10, CYBB, PTGS2, IL4I1, IL6; M2 signature group: IL10, MRC1, MSR1, CD163, CSF1R, IL4I1, SIGLEC1, CD68; Cancer associated fibroblast (CAF) group: LGALS1, COL1A1, COL1A2, COL5A1, ACTA2, FAP, LRP1, CD248, COL6A1, COL6A2, COL6A3, COL11A1, CXCL12, FBLN1, LUM, MFAP5, MMP3, MMP2, PDGFRB, PDGFRA, FN1, COL1A1, COL1A2, COL4A1, COL3A1, VTN, LGALS7, LGALS9, LAMA3, LAMB3, LAMC2, TNC, COL5A1, COL11A1, LGALS3, CA9, MMP9, MMP2, MMP1, MMP3, MMP12, MMP7, MMP11, PLOD2, ADAMTS4, ADAMTS5, LOX; Angiogenesis group: VEGFA, VEGFB, VEGFC, PDGFC, CXCL8, CXCR2, FLT1, PGF, KDR, ANGPT1, ANGPT2, TEK, VWF, CDH5; Proliferation rate group: MKI67, ESCO2, CETN3, CDK2, CCND1, CCNE1, AURKA, AURKB, E2F1, MYBL2, BUB1, CCNB1, MCM2, MCM6; and Lgr5 ISC group: ABTB2, AFAP1L1, APCDD1, ARHGEF4, ARNT2, AXIN2, BCL2, BEX1, BEX2, CAP2, CCDC46, CYP2E1, DGKG, DLGAP1, DTL, DYNC2H1, EPHA4, FAM64A, FGFR4, FMNL2, FSTL1, GRAMD1A, GRK4, IGF1R, IGFBP4, IL17RD, KIF12, KIF26B, KLHL13, LDHB, LGR5, LIFR, LOC285141, MDFIC, MPP3, NPNT, PITPNC1, PLP1, RASSF4, RNF157, SCN2B, SEPT6, SERTAD4, SLC1A2, SLC38A4, SLCO3A1, SLIT2, SOAT1, SORBS2, SOX4, TACC1, TMEM182, TNFRSF19, UTRN, ZNF141, ZNF273, ZNF493, ZNF626, ZNF678, ZNF680, ZNF714, ZNF85, ZNF92, ZNF93. A list of gene groups is provided in Table 1 below: Table 1: List of Gene Groups, the left column providing the name of the Gene Group and the right column providing examples of genes in the Gene Group. In some embodiments, a GC TME signature may include scores for two or more of the gene groups in this table. In some embodiments, a GC TME signature may include scores for the gene groups in this table that are denoted by bold text. It should be noted that the names of the gene groups shown in Table 1 appear in some of the figures (e.g., Figs.5A, 5B, 7, 8A, 15, 16, 18, and 19) with “_” instead of spaces due to how the graphics were generated. For example, “T cells” appears as “T_cells” in those figures. Also “M1 signature” is labeled “M1_signatures” in the figures.

In some embodiments, expression data from a subject is further analyzed after a GC TME type has been assigned. The disclosure is based, in part, on secondary GC TME signatures, which can used for identifying a shift in a subject’s TME type and/or subject’s shift in response to an immunotherapy. In some embodiments, a second GC TME signature is generated using RNA expression levels for genes in the plurality of gene groups listed in Table 2. In some embodiments, the second GC TME signature comprises RNA expression levels for at least three (e.g., at least 2, 3, 4, 5, between 2 and 10 and any integer therebetween, or more than 10) genes from each of at least two of the following gene groups: (i) Th17 signature group; (ii) Gamma-delta (γΔ) T cells group; (iii) Tumor suppressor group; (iv) Oncogenes group (v) Myc targets group; (vi) Th2 signature group; and (vii) Th1 signature group. In some embodiments, the second GC TME signature comprises RNA expression levels for each of the genes from each of the following gene groups: (i) Th17 signature group; (ii) Gamma-delta (γΔ) T cells group; (iii) Tumor suppressor group; (iv) Oncogenes group (v) Myc targets group; (vi) Th2 signature group; and (vii) Th1 signature group. The gene group scores of the second GC TME signature may be determined using gene set enrichment (GSEA) (e.g., using a single sample GSEA technique). For example, in some embodiments, ssGSEA is performed on RNA expression data using one or more (e.g., 1, 2, 3, 4, 5, 6, or 7) gene groups listed in Table 2. In some embodiments, each gene group of Table 2 comprises one or more (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, or more) genes listed in Table 2. In some embodiments, the additional analysis comprises performing ssGSEA on the Th17 gene group and gamma-delta T cells (γΔ T cells) gene group to produce a Th17 signature and a Gamma-delta (γΔ) T cell signature. Aspects of GSEA and ssGSEA are described herein. In some embodiments, the additional analysis comprises performing cell deconvolution analysis and/or histological analysis (e.g., by assessing Lauren subtypes) of the expression data. In some embodiments, determining one or more additional signatures (e.g., a second GC TME signature) comprises determining a respective gene group score for each of at least one of the following gene groups, using, for a particular gene group, RNA expression levels for at least two genes in the particular gene group to determine the gene group score for the particular group, the gene groups including: (i) Th17 signature group: IL17F, IL17A, IL26, IL22, RORC, BATF, IL23R, CCR6; (ii) Gamma-delta (γΔ) T cells group: TRGC1, TRDC, TRDV2, TRGC2, TRDV1; (iii) Tumor suppressors group: LATS1, LATS2, FBXW7, BAP1, KLF4, SASH1, CASP3, MLKL, SOCS1; (iv) Oncogenes group: KRAS, AKT1, EGFR, PIK3CA, PIK3CB, MYC, VAV3, MET, CCND1, INPPL1, CHD1L, RACGAP1; (v) MYC targets group: AIMP2, BYSL, CBX3, CDK4, DCTPP1, DDX18, DUSP2, EXOSC5, FARSA, GNL3, GRWD1, HK2, HSPD1, HSPE1, IMP4, IPO4, LAS1L, MAP3K6, MCM4, MCM5, MPHOSPH10, MRTO4, MYBBP1A, MYC, NDUFAF4, NIP7, NOC4L, NOLC1, NOP16, NOP2, NOP56, NPM1, PA2G4, PES1, PHB, PLK1, PLK4, PPAN, PPRC1, PRMT3, PUS1, RABEPK, RCL1, RRP12, RRP9, SLC19A1, SLC29A2, SORD, SRM, SUPV3L1, TBRG4, TCOF1, TFB2M, TMEM97, UNG, UTP20, WDR43, WDR74; (vi) Th2 signature group: IL4, IL5, IL13, IL10, GATA3, CCR4; and (vii) Th1 signature group: IFNG, IL2, CD40LG, IL21, TBX21, STAT4, IL12RB2 Table 2: List of Gene Groups, the left column providing the name of the Gene Group and the right column providing examples of genes in the Gene Group. It should be noted that the names of the gene groups shown in Table 2 appear in some of the figures (e.g., Figs.15, 16, 18, and 19) with “_” instead of spaces due to how the graphics were generated. For example, “Th17 signature” appears as “Th17_signature” in those figures. Also “MYC targets” is labeled “MYC_targets_2” in the figures. As described above, aspects of the disclosure relate to determining an GC TME signature for a subject. That signature may include gene group scores (e.g., gene group scores generated using RNA expression data for gene groups listed in Table 1 and/or Table 2). Aspects of determining of GC TME signatures is described next with reference to FIG.3. In some embodiments, a GC TME signature comprises gene group scores generated using a gene set enrichment analysis (GSEA) technique to determine a gene group score for one or more (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20) gene groups listed in Table 1. In some embodiments, a GC TME signature comprises gene group scores generated using a gene set enrichment analysis (GSEA) technique to determine a gene group score for eight or more (e.g., 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20) gene groups listed in Table 1. In some embodiments, each gene group score is generated using a gene set enrichment analysis (GSEA) technique, using RNA expression levels of at least some genes in the gene group. In some embodiments, using a GSEA technique comprises using single-sample GSEA. Aspects of single sample GSEA (ssGSEA) are described in Barbie et al. Nature.2009 Nov 5; 462(7269): 108–112, the entire contents of which are incorporated by reference herein. In some embodiments, ssGSEA is performed according to the following formula: where r i represents the rank of the ith gene in expression matrix, where N represents the number of genes in the gene set (e.g., the number of genes in the first gene group when ssGSEA is being used to determine a gene group score for the first gene group using expression levels of the genes in the first gene group), and where M represents total number of genes in expression matrix. Additional, suitable techniques of performing GSEA are known in the art and are contemplated for use in the methods described herein without limitation. In some embodiments, a GC TME signature is calculated by performing ssGSEA on expression data from a plurality of subjects, for example expression data from one or more cohorts of subjects, such as GSE15459, GSE34942, GSE26253, GSE62254, GSE13861, GSE26901, GSE29272, GSE84437, GSE26899, GSE28541, GSE113255, PRJEB25780, SRP219269, and TCGA – stomach samples (e.g., STAD TCGA project), in order to produce a plurality of enrichment scores. FIG.3 depicts an illustrative example of how gene group scores may be determined as part of act 108 of process 100. As shown in the example of FIG.3, a “GC TME signature” comprises multiple gene group scores 320 determined for respective multiple gene groups. Each gene group score, for a particular gene group, is computed by performing GSEA 310 (e.g., using ssGSEA) on RNA expression data for one or more (e.g., at least two, at least three, at least four, at least five, at least six, etc., or all) genes in the particular gene group 300. For example, as shown in FIG.3, a gene group score (labelled “Gene Group Score 1”) for gene group 1 (e.g., the Treg group) is computed from RNA expression data for one or more genes in gene group 1. As another example, a gene group score (labelled “Gene Group Score 2”) for gene group 2 (e.g., the T cells group) is computed from RNA expression data for one or more genes in gene group 2. As another example, a gene group score (labelled “Gene Group Score 3”) for gene group 3 (e.g., the NK cells group) is computed from RNA expression data for one or more genes in gene group 3. As another example, a gene group score (labelled “Gene Group Score 4”) for gene group 4 (e.g., the B cells group) is computed from RNA expression data for one or more genes in gene group 4. As another example, a gene group score (labelled “Gene Group Score 5”) for gene group 5 (e.g., the MDSC group) is computed from RNA expression data for one or more genes in gene group 5. As another example, a gene group score (labelled “Gene Group Score 6”) for gene group 6 (e.g., the CAF group) is computed from RNA expression data for one or more genes in gene group 6. As another example, a gene group score (labelled “Gene Group Score 7”) for gene group 7 (e.g., the Proliferation rate group) is computed from RNA expression data for one or more genes in gene group 7. As another example, a gene group score (labelled “Gene Group Score 8”) for gene group 8 (e.g., the Lgr5 ISC group) is computed from RNA expression data for one or more genes in gene group 8. Although the example of FIG.3 shows that the gene expression group score includes eight gene group scores for a respective set of eight gene groups, it should be appreciated that in other embodiments, the first gene expression signature may include scores for any suitable number of groups (e.g., not just 8; the number of groups could be fewer or greater than 8). As indicated by the vertical ellipsis in FIG.3, determining gene group scores of a GC TME signature may comprise determining gene group scores for 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, or more gene groups using RNA expression data from one or more respective genes in each respective gene group, as aspects of the technology described herein are not limited in this respect. In another example, a GC TME signature may include scores for only a subset of the gene groups listed in Table 1 above. As another example, the gene expression group score may include one or more scores for one or more gene groups other than those gene groups listed in Table 1 (either in addition to the score(s) for the groups in Table 1 or instead of one or more of the scores for the groups in Table 1, for example the gene groups listed in Table 2). In some embodiments, RNA expression levels for a particular gene group may be embodied in at least one data structure having fields storing the expression levels. The data structure or data structures may be provided as input to software comprising code that implements a GSEA technique (e.g., the ssGSEA technique) and processes the expression levels in the at least one data structure to compute a score for the particular gene group. The number of genes in a gene group used to determine a gene group score may vary. In some embodiments, all RNA expression levels for all genes in a particular gene group may be used to determine a gene group score for the particular gene group. In other embodiments, RNA expression data for fewer than all genes may be used (e.g., RNA expression levels for at least two genes, at least three genes, at least five genes, between 2 and 10 genes, between 5 and 15 genes, or any other suitable range within these ranges). In some embodiments, RNA expression levels for a particular gene group may be embodied in at least one data structure having fields storing the expression levels. The data structure or data structures may be provided as input to software comprising code that is configured to perform suitable scaling (e.g., median scaling) to produce a score for the particular gene group. In some embodiments, ssGSEA is performed on expression data comprising three or more (e.g., 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20) gene groups set forth in Table 1. In some embodiments, each of the gene groups separately comprises one or more (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72, 73, 74, 75, 76, 77, 78, 79, 80, or more) genes listed in Table 1. In some embodiments, a GC TME signature is produced by performing ssGSEA on all 20 of the gene groups in Table 1, each gene group including all listed genes in Table 1. In some embodiments, one or more (e.g., a plurality) of enrichment scores are normalized in order to produce a GC TME signature for the expression data (e.g., expression data of the subject or of a cohort of subjects). In some embodiments, the enrichment scores are normalized by median scaling. In some embodiments, the enrichment scores are normalized by rank estimation and median scaling. In some embodiments, median scaling comprises clipping the range of enrichment scores, for example clipping to about -1.0 to about +1.0, -2.0 to about +3.0, -3.0 to about +3.0, -4.0 to +4.0, -5.0 to about +5.0. In some embodiments, median scaling produces a GC TME signature of the subject. In some embodiments, a GC TME signature of a subject processed using a clustering algorithm to identify a GC tumor microenvironment type (e.g. a GC TME type). In some embodiments, the clustering comprises unsupervised clustering. In some embodiments, the unsupervised clustering comprises a dense clustering approach. In some embodiments, the unsupervised clustering comprises a hierarchical clustering approach. In some embodiments, clustering comprises calculating intersample similarity (e.g., using a Pearson correlation coefficient that, for example, may take on values in the range of [-1,1]), converting the distance matrix into a graph where each sample forms a node and two nodes form an edge with a weight equal to their Pearson correlation coefficient, removing edges with weight lower than a specified threshold, and applying a Louvain community detection algorithm to calculate graph partitioning into clusters. In some embodiments, the optimum weight threshold for observed clusters was calculated by employing minimum DaviesBouldin, maximum Calinski-Harabasz, and Silhouette techniques. In some embodiments, separations with low-populated clusters (< 5% of samples) are excluded. In some embodiments, a GC TME signature of a subject is compared to pre-existing clusters of GC TME types and assigned a GC TME type based on that comparison. Some aspects of determining gene group scores for gene groups are also described in U.S. Patent Publication No.2020-0273543, entitled “SYSTEMS AND METHODS FOR GENERATING, VISUALIZING AND CLASSIFYING MOLECULAR FUNCTIONAL PROFILES”, the entire contents of which are incorporated by reference herein. Generating GC TME Signature and Identifying TME Type As described herein, FIG.1 illustrates the determination of a subject’s GC TME signature, identification of the subject’s GC TME type using the GC TME signature, and identification of whether the subject is likely to respond to an immunotherapy based on the identified GC TME type. As described herein, in some embodiments, one of a plurality of different GC TME types may be identified for the subject using the GC TME signature determined for the subject using the techniques described herein. In some embodiments, the GC TME type comprises GC TME type A, GC TME type B, GC TME type C, GC TME type D, and GC TME type E, as described herein and further below. In some embodiments, each of the plurality of GC TME types is associated with a respective GC TME signature cluster in a plurality of GC TME signature clusters. The GC TME type for a subject may be determined by: (1) associating the GC TME signature of the subject with a particular one of the plurality of GC TME signature clusters; and (2) identifying the GC TME type for the subject as the GC TME type corresponding to the particular one of the plurality of GC TME signature clusters to which the GC TME signature of the subject is associated. FIG.4 shows an illustrative GC TME signature 400. In some embodiments, the GC TME signature comprises at least eight gene group scores for the following gene groups: NK cell group, T cell group, B cell group, Treg cells group, MDSC group, CAF group, Proliferation rate group, and the Lgr5 ISC group. However, it should be appreciated, that a GC TME signature may include fewer scores than the number of scores shown in FIG.4 (e.g., by omitting scores for one or more of the gene groups listed in Table 1) or more scores than the number of scores shown in FIG.4 (e.g., by including scores for one or more other gene groups in addition to or instead of the gene groups listed in Table 1). In some embodiments, a GC TME signature may be embodied in at least one data structure comprising fields storing the gene group scores part of the GC TME signature. In some embodiments, the GC TME signature clusters may be generated by: (1) obtaining GC TME signatures (using the techniques described herein) for a plurality of subjects; and (2) clustering the GC TME signatures so obtained into the plurality of clusters. Any suitable clustering technique may be used for this purpose including, but not limited to, a dense clustering algorithm, spectral clustering algorithm, k-means clustering algorithm, hierarchical clustering algorithm, and/or an agglomerative clustering algorithm. For example, intersample similarity may be calculated using a Pearson correlation. A distance matrix may be converted into a graph where each sample forms a node and two nodes form an edge with a weight equal to their Pearson correlation coefficient. Edges with weight lower than a specified threshold may be removed. A Louvain community detection algorithm may be applied to calculate graph partitioning into clusters. To mathematically determine the optimum weight threshold for observed clusters minimum DaviesBouldin, maximum Calinski- Harabasz, and Silhouette techniques may be employed. Separations with low-populated clusters (< 5% of samples) may be excluded. Accordingly, in some embodiments, generating the GC TME signature clusters involves: (A) obtaining multiple sets of RNA expression data obtained by sequencing biological samples from multiple respective subjects, each of the multiple sets of RNA expression data indicating RNA expression levels for genes in a first plurality of gene groups (e.g., one or more of the gene groups in Table 1); (B) generating multiple GC TME signatures from the multiple sets of RNA expression data, each of the multiple GC TME signatures comprising gene group scores for respective gene groups, the generating comprising, for each particular one of the multiple TME signatures: (i) determining the GC TME signature by determining the gene group scores using the RNA expression levels in the particular set of RNA expression data for which the particular one GC TME signature is being generated, and (ii) clustering the multiple GC signatures to obtain the plurality of GC TME signature clusters. The resulting GC TME signature clusters may each contain any suitable number of GC TME signatures (e.g., at least 10, at least 100, at least 500, at least 500, at least 1000, at least 5000, between 100 and 10,000, between 500 and 20,000, or any other suitable range within these ranges), as aspects of the technology described herein are not limited in this respect. The number of GC TME signature clusters in this example is five. And although, in some embodiments, it may be possible that the number of clusters is different, it should be appreciated that an important aspect of the present disclosure is the inventors’ discovery that GC may be characterized into five types based upon the generation of GC TME signatures using methods described herein. For example, as shown in FIG.4, a subject’s GC TME signature 400 may be associated with one of five GC TME clusters: 402, 404, 406, 408, and 410. Each of the clusters 402, 404, 406, 408and 410 may be associated with respective GC TME type. In this example, the GC TME signature 400 is compared to each cluster (e.g., using a distance-based comparison or any other suitable metric) and, based on the result of the comparison, the GC TME signature 400 is associated with the closest GC signature cluster (when a distance-based comparison is performed, or the “closest” in the sense of whatever metric or measure of distance is used). In this example, GC TME signature 400 is associated with GC TME Type Cluster 5410 (as shown by the consistent shading) because the measure of distance D5 between the GC TME signature 400 and (e.g., a centroid or other point representative of) cluster 410 is smaller than the measures of the distance D1, D2, D3, and D4 between the GC TME signature 400 and (e.g., a centroid or other point(s) representative of) clusters 402, 404, 406, and 408, respectively. In some embodiments, a subject’s GC TME signature may be associated with one of five GC TME signature clusters by using a machine learning technique (e.g., such as k-nearest neighbors (KNN) or any other suitable classifier) to assign the GC TME signature to one of the four GC TME signature clusters. The machine learning technique may be trained to assign GC TME signatures on the meta-cohorts represented by the signatures in the clusters. In some embodiments, GC TME types include GC TME type A, GC TME type B, GC TME type C, GC TME type D, and GC TME type E. The GC TME types described herein may be described by qualitative characteristics, for example high signals for certain gene expression signatures or scores or low signals for certain other gene expression signatures or scores. In some embodiments, a “high” signal refers to a gene expression signal or score (e.g., an enrichment score) that is at least 1-fold, 2-fold, 3-fold, 4-fold, 5-fold, 6-fold, 7-fold, 8-fold, 9- fold, 10-fold, 20-fold, 50-fold, 100-fold, 1000-fold, or more increased relative to the score of the same gene or gene group in a subject having a different type of GC. In some embodiments, a “low” signal refers to a gene expression signal or score (e.g., an enrichment score,) that is at least 1-fold, 2-fold, 3-fold, 4-fold, 5-fold, 6-fold, 7-fold, 8-fold, 9-fold, 10-fold, 20-fold, 50-fold, 100-fold, 1000-fold, or more decreased relative to the score of the same gene or gene group in a subject having a different type of GC TME. Without wishing to be bound by any theory, the tumor microenvironment of GC may contain variable numbers of immune cells, stromal cells, blood vessels and extracellular matrix. In some embodiments, “GC TME type A” is characterized by “Mesenchymal, EMT” type. In some embodiments, GC TME type A (mesenchymal) type is characterized by having the highest LGR5+ stem cell signature relative to any other GC TME types. In some embodiments, GC TME type A GC subjects also have high stromal compartment gene group signatures, for example cancer-associated fibroblasts (CAF) and angiogenesis, relative to other GC TME types. In some embodiments, GC TME type A subject are characterized by high levels of pro- tumor cytokines and medium immune compartment signatures. In some embodiments, GC TME type A subjects have low tumor proliferation rates (relative to other GC TME types). In some embodiments, type A GC tumors are diffuse (by e.g. as assessed by Lauren classification). GC TME type A has been observed, in some embodiments, to be the most aggressive GC TME type, and subjects having GC TME type A have been observed in some embodiments to have worse prognosis (e.g., by overall survival (OS)) relative to other GC TME types. In some embodiments, “GC TME type B” is also referred to as “Immune-enriched, non- fibrotic (IE)” type. In some embodiments, GC TME type B (IE) is characterized by high levels of immune-active infiltrate with a significant number of effector cells and NK cells relative to other GC TME types. In some embodiments, GC TME type B subjects have the highest tumor mutational burden (TMB). In some embodiments, subjects having GC TME type B have been observed to have better prognosis (e.g., by overall survival (OS)) relative to other GC TME types. In some embodiments, “GC TME type C” is also referred to as “Fibrotic” or “immune non-inflamed” type. In some embodiments, GC TME type C is characterized by a highly fibrotic tumor microenvironment with dense collagen formation. In some embodiments, GC TME type C is characterized by minimal lymphocyte infiltration with elevated angiogenesis, relative to other GC TME types. In some embodiments, Cancer-associated fibroblasts (CAFs) signatures are abundant in GC TME type C. In some embodiments, GC TME type C is characterized by a high level of both protumor and antitumor cytokines. In some embodiments, subjects having GC TME type C have a poorer prognosis (e.g., as measured by OS) than subjects having other GC TME types. In some embodiments, GC TME type D is also referred to as “Immune desert” type. In some embodiments, GC TME type D is characterized by the highest malignant cell percentage relative to other GC TME types. In some embodiments, GC TME type D is characterized by minimal or complete absence of leukocyte/lymphocyte infiltration (e.g., relative to other GC TME types). In some embodiments, GC TME type D is characterized by a high level of tumor proliferation rate relative to other GC TME types. In some embodiments, GC TME type E is also referred to as “B-cell enriched” type. In some embodiments, GC TME type E is characterized by high levels of immune infiltrate with a significant number of B-cells, relative to other GC TME types. In some embodiments, GC TME type D is characterized by having a low proliferation rate relative to other GC TME types. In some embodiments, subjects having GC TME type E have been observed to have better prognosis (e.g., by overall survival (OS)) relative to other GC TME types. Tables 3, 4, 5, 6, and 7 below describe examples of GC TME signatures and gene group scores produced by ssGSEA analysis and normalization (e.g., median scaling) of expression data from one or more GC subjects. Table 3: Statistics of GC TME signatures of samples having GC TME type A. The statistics in the table show, for each gene group score in the GC TME signature, the mean, minimum, maximum, 25th percentile, 50th percentile, and 75% percentile values.

Table 4: Statistics of GC TME signatures of samples having GC TME type B. The statistics in the table show, for each gene group score in the GC TME signature, the mean, minimum, maximum, 25th percentile, 50th percentile, and 75% percentile values.

Table 5: Statistics of GC TME signatures of samples having GC TME type C. The statistics in the table show, for each gene group score in the GC TME signature, the mean, minimum, maximum, 25th percentile, 50th percentile, and 75% percentile values. Table 6: Statistics of GC TME signatures of samples having GC TME type D. The statistics in the table show, for each gene group score in the GC TME signature, the mean, minimum, maximum, 25th percentile, 50th percentile, and 75% percentile values. Table 7: Statistics of GC TME signatures of samples having GC TME type E. The statistics in the table show, for each gene group score in the GC TME signature, the mean, minimum, maximum, 25th percentile, 50th percentile, and 75% percentile values.

In some embodiments, the present disclosure provides methods for identifying a subject having, suspected of having, or at risk of having GC as having an increased likelihood of having a good prognosis (e.g., as measured by overall survival (OS) or progression-free survival (PFS). In some embodiments, the method comprises determining a GC TME type of the subject as described herein. In some embodiments, the methods comprise identifying the subject as having a decreased risk of GC progression relative to other GC TME types when the subject is assigned GC TME type E. In some embodiments, “decreased risk of GC progression” may indicate better prognosis of GC or decreased likelihood of having advanced disease in a subject. In some embodiments, “decreased risk of GC progression” may indicate that the subject who has GC is expected to be more responsive to certain treatments. For instance, “decreased risk of GC progression” indicates that a subject is at least 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or 100% likely to experience a progression-free survival event (e.g., relapse, retreatment, or death) than another GC patient or population of GC patients (e.g., patients having GC, but not the same GC TME type as the subject). In some embodiments, the methods further comprise identifying the subject as having an increased risk of GC progression relative to other GC TME types when the subject is assigned a GC TME type other than GC TME type E. In some embodiments, “increased risk of GC progression” may indicate less positive prognosis of GC or increased likelihood of having advanced disease in a subject. In some embodiments, “increased risk of GC progression” may indicate that the subject who has GC is expected to be less responsive or unresponsive to certain treatments and show less or no improvements of disease symptoms. For instance, “increased risk of GC progression” indicates that a subject is at least 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, or 100% more likely to experience a progression-free survival event (e.g., relapse, retreatment, or death) than another GC patient or population of GC patients (e.g., patients having GC, but not the same GC TME type as the subject). The disclosure is based, in part, on the recognition that subjects having certain GC TME types (e.g., GC TME type E) are characterized as having (or being more likely to have) tertiary lymphoid structures (TLS), compared to subjects having other types of GC TME. Tertiary lymphoid structures are ectopic lymphoid organs that develop in non-lymphoid tissues at sites of chronic inflammation including tumors, for example as described by Sautès-Fridman et al. Nat Rev Cancer 19, 307–325 (2019). In the context of certain cancers (e.g., gastric cancers) the presence of TLS is associated with an improved subject prognosis. Accordingly, in some embodiments, subjects identified as having GC TME type E have an increased likelihood of having PFS or increased overall survival (OS) relative to subjects having other GC TME types. In some embodiments, the methods described herein comprise the use of at least one computer hardware processor to perform the determination. In some embodiments, the present disclosure provides a method for providing a prognosis, predicting survival, or stratifying patient risk of a subject suspected of having, or at risk of having GC. In some embodiments, the method comprises determining a GC TME type of the subject as described herein. Updating GC TME Clusters Based on New Data Techniques for generating GC TME clusters are described herein. It should be appreciated that the GC TME clusters may be updated as additional GC TME signatures are computed for patients. In some embodiments, the GC TME signature of the subject is one of a threshold number GC TME signatures for a threshold number of subjects. In some embodiments, when the threshold number of GC TME signatures is generated the GC TME signature clusters are updated. For example, once a threshold number of new GC TME signatures are obtained (e.g., 1 new signature, 10 new signatures, 100 new signatures, 500 new signatures, any suitable threshold number of signatures in the range of 10-1,000 signatures), the new signatures may be combined with the GC TME signatures previously used to generate the GC TME clusters and the combined set of old and new GC TME signatures may be clustered again (e.g., using any of the clustering algorithms described herein or any other suitable clustering algorithm) to obtain an updated set of GC TME signature clusters. In this way, data obtained from a future patient may be analyzed in a way that takes advantage of information learned from patients whose GC TME signature was computed prior to that of the future patient. In this sense, the machine learning techniques described herein (e.g., the unsupervised clustering machine learning techniques) are adaptive and learn with the accumulation of new patient data. This facilitates improved characterization of the GC TME type that future patients may have and may improve the selection of treatment for those patients. Therapeutic Indications Aspects of the disclosure relate to methods of identifying or selecting a therapeutic agent for a subject based upon determination of the subject’s GC TME type. The disclosure is based, in part, on the recognition that subjects having GC TME type B or GC TME type E have an increased likelihood of responding to certain therapies (e.g., an immunotherapy) relative to subjects having other GC TME types. The inventors have surprisingly discovered that determination of GC TME type E or GC TME type B patients are likely to respond to immunotherapy independent of their MSS/MSI or EBV status. In some aspects, the disclosure provides a method for predicting disease control rate (DCR) to an immunotherapy of a subject, the method comprising determining the GC TME type of the subject, and identifying the subject as likely to have an increased DCR when the subject is assigned GC TME type B or GC TME type E relative to subjects having GC TME type A, GC TME type C, or GC TME type D. As used herein, “disease control rate” or “DCR” refers to the percentage of patients with advanced or metastatic cancer who have achieved complete response, partial response, and stable disease to a therapeutic intervention. In some embodiments, a subject determined to have GC TME type B or GC TME type E has a DCR greater than 50% (e.g., at least 50%, 60%, 70%, 80%, 90%, 95%, or 99%) for an immunotherapy. In some embodiments, the immunotherapy comprises a PD1 inhibitor. Examples of PD1 inhibitors include but are not limited to pembrolizumab, nivolumab, cemilimab, JTX-4014, camrelizumab, sintilimab, tisleizumab, toripalimab, dostarlimab, AMP-224, and AMP-514. In some embodiments, the immunotherapy comprises pembrolizumab. The disclosure is based, in part, on the recognition that certain subjects having GC TME type B or GC TME type E may have a higher likelihood of responding to an immunotherapy (e.g., pembrolizumab) even if that subject has previously been characterized as Epstein-Barr virus (EBV)-low and/or microsatellite stable (MSS). This is surprising because such patients have previously been observed to be low responders to immunotherapies. Aspects of the disclosure relate to methods of treating a subject having (or suspected or at risk of having) GC based upon a determination of the GC TME type of the subject. In some embodiments, the methods comprise administering one or more (e.g., 1, 2, 3, 4, 5, or more) therapeutic agents to the subject. In some embodiments, the therapeutic agent (or agents) administered to the subject are selected from small molecules, peptides, nucleic acids, radioisotopes, cells (e.g., CAR T-cells, etc.), and combinations thereof. Examples of therapeutic agents include chemotherapies (e.g., cytotoxic agents, etc.), immunotherapies (e.g., immune checkpoint inhibitors, such as PD-1 inhibitors, PD-L1 inhibitors, etc.), antibodies (e.g., anti- HER2 antibodies), cellular therapies (e.g. CAR T-cell therapies), gene silencing therapies (e.g., interfering RNAs, CRISPR, etc.), antibody-drug conjugates (ADCs), and combinations thereof. In some embodiments, a subject is administered an effective amount of a therapeutic agent. “An effective amount” as used herein refers to the amount of each active agent required to confer therapeutic effect on the subject, either alone or in combination with one or more other active agents. Effective amounts vary, as recognized by those skilled in the art, depending on the particular condition being treated, the severity of the condition, the individual patient parameters including age, physical condition, size, gender and weight, the duration of the treatment, the nature of concurrent therapy (if any), the specific route of administration and like factors within the knowledge and expertise of the health practitioner. These factors are well known to those of ordinary skill in the art and can be addressed with no more than routine experimentation. It is generally preferred that a maximum dose of the individual components or combinations thereof be used, that is, the highest safe dose according to sound medical judgment. It will be understood by those of ordinary skill in the art, however, that a patient may insist upon a lower dose or tolerable dose for medical reasons, psychological reasons, or for virtually any other reasons. Empirical considerations, such as the half-life of a therapeutic compound, generally contribute to the determination of the dosage. For example, antibodies that are compatible with the human immune system, such as humanized antibodies or fully human antibodies, may be used to prolong half-life of the antibody and to prevent the antibody being attacked by the host's immune system. Frequency of administration may be determined and adjusted over the course of therapy, and is generally (but not necessarily) based on treatment, and/or suppression, and/or amelioration, and/or delay of a cancer. Alternatively, sustained continuous release formulations of an anti-cancer therapeutic agent may be appropriate. Various formulations and devices for achieving sustained release are known in the art. In some embodiments, dosages for an anti-cancer therapeutic agent as described herein may be determined empirically in individuals who have been administered one or more doses of the anti-cancer therapeutic agent. Individuals may be administered incremental dosages of the anti-cancer therapeutic agent. To assess efficacy of an administered anti-cancer therapeutic agent, one or more aspects of a cancer (e.g., tumor microenvironment, tumor formation, tumor growth, or GC TME types, etc.) may be analyzed. Generally, for administration of any of the anti-cancer antibodies described herein, an initial candidate dosage may be about 2 mg/kg. For the purpose of the present disclosure, a typical daily dosage might range from about any of 0.1 µg/kg to 3 µg /kg to 30 µg /kg to 300 µg /kg to 3 mg/kg, to 30 mg/kg to 100 mg/kg or more, depending on the factors mentioned above. For repeated administrations over several days or longer, depending on the condition, the treatment is sustained until a desired suppression or amelioration of symptoms occurs or until sufficient therapeutic levels are achieved to alleviate a cancer, or one or more symptoms thereof. An exemplary dosing regimen comprises administering an initial dose of about 2 mg/kg, followed by a weekly maintenance dose of about 1 mg/kg of the antibody, or followed by a maintenance dose of about 1 mg/kg every other week. However, other dosage regimens may be useful, depending on the pattern of pharmacokinetic decay that the practitioner (e.g., a medical doctor) wishes to achieve. For example, dosing from one-four times a week is contemplated. In some embodiments, dosing ranging from about 3 µg /mg to about 2 mg/kg (such as about 3 µg /mg, about 10 µg /mg, about 30 µg /mg, about 100 µg /mg, about 300 µg /mg, about 1 mg/kg, and about 2 mg/kg) may be used. In some embodiments, dosing frequency is once every week, every 2 weeks, every 4 weeks, every 5 weeks, every 6 weeks, every 7 weeks, every 8 weeks, every 9 weeks, or every 10 weeks; or once every month, every 2 months, or every 3 months, or longer. The progress of this therapy may be monitored by conventional techniques and assays and/or by monitoring GC TME types as described herein. The dosing regimen (including the therapeutic used) may vary over time. When the anti-cancer therapeutic agent is not an antibody, it may be administered at the rate of about 0.1 to 300 mg/kg of the weight of the patient divided into one to three doses, or as disclosed herein. In some embodiments, for an adult patient of normal weight, doses ranging from about 0.3 to 5.00 mg/kg may be administered. The particular dosage regimen, e.g., dose, timing, and/or repetition, will depend on the particular subject and that individual's medical history, as well as the properties of the individual agents (such as the half-life of the agent, and other considerations well known in the art). For the purpose of the present disclosure, the appropriate dosage of an anti-cancer therapeutic agent will depend on the specific anti-cancer therapeutic agent(s) (or compositions thereof) employed, the type and severity of cancer, whether the anti-cancer therapeutic agent is administered for preventive or therapeutic purposes, previous therapy, the patient's clinical history and response to the anti-cancer therapeutic agent, and the discretion of the attending physician. Typically, the clinician will administer an anti-cancer therapeutic agent, such as an antibody, until a dosage is reached that achieves the desired result. Administration of an anti-cancer therapeutic agent can be continuous or intermittent, depending, for example, upon the recipient's physiological condition, whether the purpose of the administration is therapeutic or prophylactic, and other factors known to skilled practitioners. The administration of an anti-cancer therapeutic agent (e.g., an anti-cancer antibody) may be essentially continuous over a preselected period of time or may be in a series of spaced dose, e.g., either before, during, or after developing cancer. As used herein, the term “treating” refers to the application or administration of a composition including one or more active agents to a subject, who has a cancer, a symptom of a cancer, or a predisposition toward a cancer, with the purpose to cure, heal, alleviate, relieve, alter, remedy, ameliorate, improve, or affect the cancer or one or more symptoms of GC, or the predisposition toward GC. Alleviating GC includes delaying the development or progression of the disease, or reducing disease severity. Alleviating the disease does not necessarily require curative results. As used therein, “delaying” the development of a disease (e.g., a cancer) means to defer, hinder, slow, retard, stabilize, and/or postpone progression of the disease. This delay can be of varying lengths of time, depending on the history of the disease and/or individuals being treated. A method that “delays” or alleviates the development of a disease, or delays the onset of the disease, is a method that reduces probability of developing one or more symptoms of the disease in a given time frame and/or reduces extent of the symptoms in a given time frame, when compared to not using the method. Such comparisons are typically based on clinical studies, using a number of subjects sufficient to give a statistically significant result. “Development” or “progression” of a disease means initial manifestations and/or ensuing progression of the disease. Development of the disease can be detected and assessed using clinical techniques known in the art. Alternatively, or in addition to the clinical techniques known in the art, development of the disease may be detectable and assessed based on other criteria. However, development also refers to progression that may be undetectable. For purpose of this disclosure, development or progression refers to the biological course of the symptoms. “Development” includes occurrence, recurrence, and onset. As used herein “onset” or “occurrence” of a cancer includes initial onset and/or recurrence. Examples of the antibody anti-cancer agents include, but are not limited to, alemtuzumab (Campath), trastuzumab (Herceptin), Ibritumomab tiuxetan (Zevalin), Brentuximab vedotin (Adcetris), Ado-trastuzumab emtansine (Kadcyla), blinatumomab (Blincyto), Bevacizumab (Avastin), Cetuximab (Erbitux), ipilimumab (Yervoy), nivolumab (Opdivo), pembrolizumab (Keytruda), atezolizumab (Tecentriq), avelumab (Bavencio), durvalumab (Imfinzi), and panitumumab (Vectibix). Examples of an immunotherapy include, but are not limited to, a PD-1 inhibitor or a PD- L1 inhibitor, a CTLA-4 inhibitor, adoptive cell transfer, therapeutic cancer vaccines, oncolytic virus therapy, T-cell therapy, and immune checkpoint inhibitors. Examples of radiation therapy include, but are not limited to, ionizing radiation, gamma- radiation, neutron beam radiotherapy, electron beam radiotherapy, proton therapy, brachytherapy, systemic radioactive isotopes, and radiosensitizers. Examples of a surgical therapy include, but are not limited to, a curative surgery (e.g., tumor removal surgery), a preventive surgery, a laparoscopic surgery, and a laser surgery. Examples of the chemotherapeutic agents include, but are not limited to, R-CHOP, Carboplatin or Cisplatin, Docetaxel, Gemcitabine, Nab-Paclitaxel, Paclitaxel, Pemetrexed, and Vinorelbine. Additional examples of chemotherapy include, but are not limited to, Platinating agents, such as Carboplatin, Oxaliplatin, Cisplatin, Nedaplatin, Satraplatin, Lobaplatin, Triplatin, Tetranitrate, Picoplatin, Prolindac, Aroplatin and other derivatives; Topoisomerase I inhibitors, such as Camptothecin, Topotecan, irinotecan/SN38, rubitecan, Belotecan, and other derivatives; Topoisomerase II inhibitors, such as Etoposide (VP-16), Daunorubicin, a doxorubicin agent (e.g., doxorubicin, doxorubicin hydrochloride, doxorubicin analogs, or doxorubicin and salts or analogs thereof in liposomes), Mitoxantrone, Aclarubicin, Epirubicin, Idarubicin, Amrubicin, Amsacrine, Pirarubicin, Valrubicin, Zorubicin, Teniposide and other derivatives; Antimetabolites, such as Folic family (Methotrexate, Pemetrexed, Raltitrexed, Aminopterin, and relatives or derivatives thereof); Purine antagonists (Thioguanine, Fludarabine, Cladribine, 6-Mercaptopurine, Pentostatin, clofarabine, and relatives or derivatives thereof) and Pyrimidine antagonists (Cytarabine, Floxuridine, Azacitidine, Tegafur, Carmofur, Capacitabine, Gemcitabine, hydroxyurea, 5-Fluorouracil (5FU), and relatives or derivatives thereof); Alkylating agents, such as Nitrogen mustards (e.g., Cyclophosphamide, Melphalan, Chlorambucil, mechlorethamine, Ifosfamide, mechlorethamine, Trofosfamide, Prednimustine, Bendamustine, Uramustine, Estramustine, and relatives or derivatives thereof); nitrosoureas (e.g., Carmustine, Lomustine, Semustine, Fotemustine, Nimustine, Ranimustine, Streptozocin, and relatives or derivatives thereof); Triazenes (e.g., Dacarbazine, Altretamine, Temozolomide, and relatives or derivatives thereof); Alkyl sulphonates (e.g., Busulfan, Mannosulfan, Treosulfan, and relatives or derivatives thereof); Procarbazine; Mitobronitol, and Aziridines (e.g., Carboquone, Triaziquone, ThioTEPA, triethylenemalamine, and relatives or derivatives thereof); Antibiotics, such as Hydroxyurea, Anthracyclines (e.g., doxorubicin agent, daunorubicin, epirubicin and relatives or derivatives thereof); Anthracenediones (e.g., Mitoxantrone and relatives or derivatives thereof); Streptomyces family antibiotics (e.g., Bleomycin, Mitomycin C, Actinomycin, and Plicamycin); and ultraviolet light. In some aspects, the disclosure provides a method for treating gastric cancer, the method comprising administering one or more therapeutic agents (e.g., one or more anti-cancer agents, such as one or more chemotherapeutic agents) to a subject identified as having a particular GC TME type, wherein the GC TME type of the subject has been identified by method as described by the disclosure. In some embodiments, a subject has been identified as having GC TME type A, GC TME type B, GC TME type C, GC TME type D, or GC TME type E. In some embodiments, a subject has been identified as having GC TME type B or GC TME type E. The disclosure is based, in part, on the inventors’ recognition that subjects having certain GC TME types are likely to respond well to certain immunotherapies (e.g., immune checkpoint inhibitors, such as pembrolizumab). In some embodiments, the therapeutic agent comprises pembrolizumab when the subject has been identified as having GC TME type E or GC TME type B. Dosages of pembrolizumab are well known, for example 200 mg every 3 weeks or 400 mg every 6 weeks, by infusion over 30 minutes. Aspects of the disclosure are based on the inventors’ recognition that subjects having certain GC TME types are unlikely to respond well to certain GC therapies, such as immunotherapeutic agents. Thus, in some embodiments, the therapeutic agent comprises a therapeutic agent other than an immunotherapy when the subject has been identified as having a GC TME type C or GC TME type A. Examples of other gastric cancer therapies include but are not limited to Cyramza (Ramucirumab), Docetaxel, Doxorubicin Hydrochloride, Enhertu (Fam- Trastuzumab Deruxtecan-nxki), 5-FU (Fluorouracil Injection), Fam-Trastuzumab Deruxtecan- nxki, Fluorouracil Injection, Herceptin (Trastuzumab), Keytruda (Pembrolizumab), Lonsurf (Trifluridine and Tipiracil Hydrochloride), Mitomycin, Ramucirumab, Taxotere (Docetaxel), Trastuzumab, and Trifluridine and Tipiracil Hydrochloride. Reports In some aspects, methods disclosed herein comprise generating a report for assisting with the preparation of recommendation for prognosis and/or treatment. The generated report can provide summary of information, so that the clinician can identify the GC TME type or suitable therapy. The report as described herein may be a paper report, an electronic record, or a report in any format that is deemed suitable in the art. The report may be shown and/or stored on a computing device known in the art (e.g., handheld device, desktop computer, smart device, website, etc.). The report may be shown and/or stored on any device that is suitable as understood by a skilled person in the art. In some embodiments, methods disclosed herein can be used for commercial diagnostic purposes. For example, the generated report may include, but is limited to, information concerning expression levels of one or more genes from any of the gene groups described herein, clinical, and pathologic factors, patient’s prognostic analysis, predicted response to the treatment, classification of the GC TME environment (e.g., as belonging to one of the types described herein), the alternative treatment recommendation, and/or other information. In some embodiments, the methods and reports may include database management for the keeping of the generated reports. For instance, the methods as disclosed herein can create a record in a database for the subject (e.g., subject 1, subject 2, etc.) and populate the specific record with data for the subject. In some embodiments, the generated report can be provided to the subject and/or to the clinicians. In some embodiments, a network connection can be established to a server computer that includes the data and report for receiving or outputting. In some embodiments, the receiving and outputting of the date or report can be requested from the server computer. Computer Implementation An illustrative implementation of a computer system 1300 that may be used in connection with any of the embodiments of the technology described herein (e.g., such as the method of FIG.1) is shown in FIG.13. The computer system 1300 includes one or more processors 1310 and one or more articles of manufacture that comprise non-transitory computer- readable storage media (e.g., memory 1320 and one or more non-volatile storage media 1330). The processor 1310 may control writing data to and reading data from the memory 1320 and the non-volatile storage device 1330 in any suitable manner, as the aspects of the technology described herein are not limited to any particular techniques for writing or reading data. To perform any of the functionality described herein, the processor 1310 may execute one or more processor-executable instructions stored in one or more non-transitory computer-readable storage media (e.g., the memory 1320), which may serve as non-transitory computer-readable storage media storing processor-executable instructions for execution by the processor 1310. Computing device 1300 may also include a network input/output (I/O) interface 1340 via which the computing device may communicate with other computing devices (e.g., over a network), and may also include one or more user I/O interfaces 1350, via which the computing device may provide output to and receive input from a user. The user I/O interfaces may include devices such as a keyboard, a mouse, a microphone, a display device (e.g., a monitor or touch screen), speakers, a camera, and/or various other types of I/O devices. The above-described embodiments can be implemented in any of numerous ways. For example, the embodiments may be implemented using hardware, software, or a combination thereof. When implemented in software, the software code can be executed on any suitable processor (e.g., a microprocessor) or collection of processors, whether provided in a single computing device or distributed among multiple computing devices. It should be appreciated that any component or collection of components that perform the functions described above can be generically considered as one or more controllers that control the above-discussed functions. The one or more controllers can be implemented in numerous ways, such as with dedicated hardware, or with general purpose hardware (e.g., one or more processors) that is programmed using microcode or software to perform the functions recited above. In this respect, it should be appreciated that one implementation of the embodiments described herein comprises at least one computer-readable storage medium (e.g., RAM, ROM, EEPROM, flash memory or other memory technology, CD-ROM, digital versatile disks (DVD) or other optical disk storage, magnetic cassettes, magnetic tape, magnetic disk storage or other magnetic storage devices, or other tangible, non-transitory computer-readable storage medium) encoded with a computer program (i.e., a plurality of executable instructions) that, when executed on one or more processors, performs the above-discussed functions of one or more embodiments. The computer-readable medium may be transportable such that the program stored thereon can be loaded onto any computing device to implement aspects of the techniques discussed herein. In addition, it should be appreciated that the reference to a computer program which, when executed, performs any of the above-discussed functions, is not limited to an application program running on a host computer. Rather, the terms computer program and software are used herein in a generic sense to reference any type of computer code (e.g., application software, firmware, microcode, or any other form of computer instruction) that can be employed to program one or more processors to implement aspects of the techniques discussed herein. The foregoing description of implementations provides illustration and description but is not intended to be exhaustive or to limit the implementations to the precise form disclosed. Modifications and variations are possible in light of the above teachings or may be acquired from practice of the implementations. In other implementations the methods depicted in these figures may include fewer operations, different operations, differently ordered operations, and/or additional operations. Further, non-dependent blocks may be performed in parallel. It will be apparent that example aspects, as described above, may be implemented in many different forms of software, firmware, and hardware in the implementations illustrated in the figures. Further, certain portions of the implementations may be implemented as a “module” that performs one or more functions. This module may include hardware, such as a processor, an application-specific integrated circuit (ASIC), or a field-programmable gate array (FPGA), or a combination of hardware and software. EXAMPLES Example 1: Identification of Gastric Cancer Tumor Microenvironment (TME) This example describes an illustrative technique for generating an GC TME signature for a subject from RNA expression data for the subject, according to some embodiments of the technology described herein. The produced GC TME signature reflects and/or indicates the abundance of both the malignant and microenvironment (TME) cell subpopulations and the activity of tumor-promoting and tumor-suppressive processes occurring within a tumor, and constitutes a personalized tumor map. The generated GC TME signature for the subject is used to identify an GC TME type for the subject from among five GC TME types: GC TME type A, GC TME type B, GC TME type C, GC TME type D, and GC TME type E. Aspects of some of the steps of the process described in this example are described in further detail herein including with reference to FIGs.1-4 above. RNA expression data (including both RNA-seq and microarray expression data) were obtained from multiple public databases. Data were subjected to basic quality control (QC) measures. For example, outlier samples and samples with signs of RNA degradation were excluded. Preprocessing of expression data included normalization and log-transformation. For microarrays normalization is performed automatically using gcrma package. RNA-seq data was subsequently normalized to TPM (transcript per million) units. TPM normalization techniques are described in Wagner et al. (Theory Biosci. (2012) 131:281–285), which is incorporated by reference herein in its entirety. TPM normalization may be performed using a software package, such as, for example, the gcrma package. Aspects of the gcrma package are described in Wu J, Gentry RIwcfJMJ (2021). “gcrma: Background Adjustment Using Sequence Information. R package version 2.66.0.”, which is incorporated by reference in its entirety herein. In some embodiments, RNA expression level in TPM units for a particular gene may be calculated according to: The GC TME type determined for a subject is determined using a GC TME signature. The GC TME signature includes gene group scores (e.g., one or more of the gene groups described in Table 1) obtained using ssGSEA. The gene group scores in the GC TME signatures were calculated from log-transformed RNA expression valued. The ssGSEA was performed according to the following formula: where ri represents the rank of the ith gene in expression matrix, where N represents the number of genes in the gene set (e.g., the number of genes in the first gene group when ssGSEA is being used to determine a score for the first gene group using expression levels of the genes in the first gene group), and where M represents total number of genes in expression matrix. Additional, suitable techniques of performing GSEA are known in the art and are contemplated for use in the methods described herein without limitation. In some embodiments, a GC signature is calculated by performing ssGSEA on expression data from a plurality of subjects, for example expression data from one or more cohorts of subjects. After calculation, the scores were scaled using rank estimation and median-scaling, which was important for removing undesirable batch effects and to enable all the datasets to be combined together. Median scaling consisted of estimating median and MAD (median absolute deviation) for each signature within each dataset, and applying the formula xi- median(x)/MAD(x). After GC TME signatures were calculated according to the above for multiple patients, unsupervised clustering was performed to generate GC TME clusters. To classify a new sample, it is grouped together with the dataset used to get the GC TME type. Scores are calculated for the sample and scaled together with the selected cohort. After that the sample type, can be predicted by applying a machine learning model (e.g., K-nearest neighbor, “knn”) trained on the scaled metacohort. In this example, inter-sample similarity was calculated using Pearson correlation. The distance matrix was converted into a graph where each sample formed a node and two nodes formed an edge with a weight equal to their Pearson correlation coefficient. Edges with weight lower than a specified threshold were removed. The Louvain community detection algorithm was applied to calculate graph partitioning into clusters. To mathematically determine the optimum weight threshold for observed clusters, minimum DaviesBouldin, maximum Calinski- Harabasz, and Silhouette techniques were employed. Separations with low-populated clusters (< 5% of samples) were excluded. This analysis resulted in the generation of five (5) GC TME signature clusters, corresponding to five (5) GC TME types. Example of Characterizing GC TME types Using the aforementioned approach on several publicly available cancer data sets, five distinct types of GC were observed (FIG.5A). To identify GC TME types, a meta-cohort of gene expression data was collected from public datasets. GC TME signatures comprising twenty gene group scores were used to estimate different biological processes in each of the samples. Genes included in each gene group are described in Table 1. To overcome batch effect from different datasets, rank estimation and median scaling transformation were used prior to clustering, which led to identification of the five GC TME types. Examples of processes for identifying GC TME types are described in FIG. 1. Using unsupervised clustering, five stable GC TME types were identified (FIG.5A): Immune Enriched/Non fibrotic (type B), Depleted (type D), Mesenchymal/EMT (type A), Fibrotic (type C), and “B cell-enriched” (type E). FIG.5B provides a heatmap of pairwise signature correlation for the gene groups listed in Table 1. Data indicate the gene group scores for the NK cells group, T cells group, B cells group, Treg group, MDSC group, CAF group, Proliferation rate group, and Lgr5 ISC group may facilitate typing the GC TME of a subject. Accordingly, in some embodiments, a GC TME signature may include gene group scores for these gene groups. Though, in other embodiments, one or more gene scores for one or more other gene groups may be used in addition to or instead of these gene group scores. GC TME type E is characterized by the following features: low or medium stromal/vascularized component (e.g., lower than GC TME types A and C), medium or high immune component (e.g., similar to GC TME type B), highest B cell signal (e.g., of all GC TME types), and low or medium neoplasm properties (e.g., low expression level of Proliferation rate and Oncogenes signatures, higher level of Tumor suppressors signature, relative to other GC TME types). In some embodiments, GC TME type is prognostic for patient outcome. An analysis between GC TME types was performed and compared to TCGA histological data. FIG.6 provides selected process signature comparisons across GC TME types. Data indicates that while GC TME types A and C have the highest Angiogenesis and stromal (CAF = Cancer- associated fibroblasts) signals, GC TME type E represents the highest B cells signal and similar T cell signals. FIG.7 provides an exemplary heatmap of gastric cancer TME type gene group signatures and cell deconvolution across the samples. For RNA-seq samples, high B cell content was also supported by a cell deconvolution algorithm-based analysis. This algorithm allows reconstructing cell composition from bulk RNA-seq data and estimating the percentage of different cell types (fibroblasts, B cells, T cells, macrophages, etc.). GC TME type E samples proved the highest B cells percentage (FIG.5A). Medians and cumulative distribution function (CDF) for each signature in each GC TME type are provided below in Tables 8 and 9. Table 8: Median signature values by TME types in Gastric Cancer

Table 9: Mean cumulative distribution function (CDF) of signatures values by TME types in Gastric Cancer

Using TCGA histology data, immune infiltration and stromal compartments were compared in samples belonging to different GC TME types (FIG.8A). Samples of GC TME types A and C had high stromal content. GC TME type D samples had desert immune, low fibroblast composition, and high tumor cellularity. GC TME type B and E samples, on the contrary, showed high immune infiltration. Remarkably, GC TME type E samples also showed the existence of germinative centers/tertiary lymphoid structures (TLSs) - ectopic lymphoid organs that develop in non-lymphoid tissues (shown by the arrow in FIG.8A). FIGs.8B-8F show graphic representations of relative cell type content of different GC TME types. FIG.8B shows a representation of the relative cell type content of GC TME type A, which is characterized by WNT-activation, and the prevalence of LGR5+ stem cells.; type A was also associated with a low tumor proliferation rate. Cancer-associated fibroblasts are often observed. Signs of epithelial-mesenchymal transition (EMT) were present. Subjects having GC TME type A have been observed to have a poor prognosis. FIG.8C shows a representation of the relative cell type content of GC TME type B, which is characterized by high levels of tumor-infiltrating immune cells, including cytotoxic effector cells. This type has the highest tumor mutational burden (TMB). High PD-L1 expression is commonly observed. This GC TME type is often responsive to immunotherapy and associated with good prognosis. FIG.8D shows a representation of the relative cell type content of GC TME type C, which is associated with high fibrosis and dense collagen formation. This type is also characterized by minimal lymphocyte infiltration (non-inflamed) with intense angiogenesis. Cancer-associated fibroblasts (CAF) are abundant. High levels of both protumor and antitumor cytokines have been observed. This GC TME type is associated with poor prognosis. FIG.8E shows a representation of the relative cell type content of GC TME type D, which is characterized by the highest percentage of malignant cells, while leukocyte/lymphocyte infiltration is only minimal or completely absent. A high tumor proliferation rate has been observed. This GC TME type is associated with good prognosis. FIG.8F shows a representation of the relative cell type content of GC TME type E, which is characterized by the presence of tertiary lymphoid structures (TLS), and high levels of immune infiltrate with a significant number of B cells. A low tumor proliferation rate has been observed. This GC TME type is associated with good prognosis. Certain GC TME types were found to be associated with malignant properties: Mesenchymal A type comprised the majority of the studied metastatic samples and was characterized by the highest signal of LGR5+ stem cell signature. TLSs have been observed in various cancers: colorectal cancer, lung cancer, clear-cell renal cell carcinoma, sarcoma, urothelial bladder cancer, and others. The existence of TLS may provide additional opportunities in treatment decision making. For example, it has been observed that TLS-positive samples in primary gastrointestinal stromal tumors demonstrated better postoperative outcomes. Analysis of overall survival (OS) and progression-free survival (PFS) revealed better prognosis for GC TME type E compared with stromal enriched and WNT activated GC TME types A and C (FIG.9). Samples of GC TME types B, E and D had the best overall survival rates. Example 2: Additional embodiments Association of STAD TME types with EBV, MSI, TMB Status GC is histologically classified into intestinal, diffuse and the mixed types, and into four molecular types based on genetic profiling (i.e., microsatellite instable (MSI), EBV positive, chromosomal instable, and genomically stable). Here, samples with various statuses of EBV, MSI, and tumor mutational burden (TMB) across the GC TME types were examined. FIG.10 shows comparison of MSS and MSI samples in different GC TME types. As shown in FIG.10, most MSI samples were located in GC TME type B (i.e., immune enriched), indicating that subjects having GC TME type B may have a higher likelihood of responding to therapies such as immunotherapies. FIGs.11A-11B show mutation burden across the GC TME types in the TCGA cohort. As shown in FIG.11A, immune enriched type has the highest mutation load (ML) in comparison to all other GC TME types. However, in the ACGR cohort, the highest malignant cell content (i.e., purity or cellularity) and ploidy were observed in the immune depleted type, GC TME type D (FIG.11B, left). Additional studies were conducted to examine the EBV status of the TCGA cohorts across all GC TME types. FIG.12 shows strong relationships between EBV-positive status (i.e., “1.0") and GC TME type B and E (immune enriched and B cell enriched, respectively). The overall results confirmed that subjects who have GC TME types B and E would be determined to have increased likelihood of having a good prognosis and/or progression-free survival (PFS). Example 3: Th17 and γδ T cells signatures for GC typing After identification of the five GC TME types, additional gene group scores were evaluated, including Th17 (T helper 17) and γδ T cells gene group scores; and gene group scores for Tumor suppressors, Oncogenes, MYC targets, Th2, and Th1. The additional gene groups used to produce the additional GC TME signatures are described in Table 2. Comparative analysis revealed higher γδ T cells signal in immune enriched GC TME types (B and E) and the lowest signal in immune desert subtype. Th17 was higher in all these three types and has the lowest signal in Mesenchymal, WNT-activated subtype A (FIG.14). The combination of the novel Th17 and γδ T cells signatures with the GC TME types of the TCGA cohort is shown in FIG.15. Dynamic changes of the GC TME types may correlate with response to chemotherapy. Example 4: Immunotherapy Immunotherapy for gastric cancer has shown promising results, indicating better response in Epstein-Barr virus (EBV)-positive samples or in microsatellite-unstable/instability (MSI) samples (or samples with high mutation burden) than other cancer types. The GC typing (particularly for types E and B) may identify responders to immunotherapy among a large population of patients having conventional biomarkers associated with an adverse response to immunotherapy (e.g., microsatellite stable (MSS) and EBV-negative). Analysis of the TCGA cohort revealed that the majority of EBV positive samples belonged to E and B types; MSI samples were also mostly located in type B (FIG.16). Indeed, 54 samples from a total 69 B-type samples in the TCGA cohort were EBV positive or MSI with high mutation burden. Fifteen samples are EBV negative and microsatellite stable (MSS) but still these samples may be better responders for immunotherapy. In type E more than a half of samples are EBV negative and MSS (FIG.16). Using a public dataset, response to immunotherapy was investigated in metastatic gastric cancer (2nd or 3rd line of pembrolizumab) across different GC TME types, considering their EBV and MSS status. It was observed that patients from types B and E had a higher Disease Control Rate (DCR): 9/9 for type B and 3/4 for type E, while other subtypes had lower: type A - 3 /6, type C - 4 /11 and type D – 1/4. Next, all patients with well-known markers for immunotherapy treatment were excluded: MSI samples and EBV positive samples. Interestingly, B and E types still included higher percentages of responders (FIG.17). All metrics were calculated for different groups and are presented in Table 10. Table 10

According to the data in FIGs.16 and 17, and Table 10, GC TME typing serves as a predictive parameter independent from MSI/MSS and EBV status. Although MSI and EBV-high patients are consistent responders to immunotherapy (Table 10), a significant percentage of MSS and EBV-low patients could still show better response rates if they belong to E- and B- subtypes (FIG.17). Since more than half of E-type patients from TCGA cohort were MSS and EBV- negative (FIG.16) GC typing may allow covering this large group of patients who previously escaped from the responder category. E- and B-type metastatic gastric cancer patients have a higher disease control rate (DCR) in response to immunotherapy (2nd or 3rd line of pembrolizumab), regardless of MSI and EBV status. All the five described subtypes were found in Cabazi-cohort (FIG.18). All subtypes have similar distribution of samples from gastric, GEJ and esophagus. The cohort comprises samples of different biopsy sites, including metastatic samples. Interestingly, the most metastatic samples (5/9) were identified as Mesenchymal, WNT-activated type (type A). In addition, these five type A metastatic samples have the highest signal of LGR5+ stem cell signature, hence the highest tumor stemness potential. Metastatic samples were also observed in types E and D. For nine samples from the Cabazi-cohort both pre- and post-treatment biopsies were sequenced in order to identify their GC TME type and analyze the impact of treatment on type- switching. Two of these samples retained their initial GC TME type after treatment (both samples had Immune desert D type). The other seven samples showed various GC TME type switching. Interestingly, the only one partial responder to cabazitaxel in this group changed fibrotic type C to Immune enriched type B (FIG.19). Table 11: Exemplary NCBI Accession Numbers for genes listed in Table 1.

Table 12: Exemplary NCBI Accession Numbers for genes listed in Table 2.

EQUIVALENTS Having thus described several aspects and embodiments of the technology set forth in the disclosure, it is to be appreciated that various alterations, modifications, and improvements will readily occur to those skilled in the art. Such alterations, modifications, and improvements are intended to be within the spirit and scope of the technology described herein. For example, those of ordinary skill in the art will readily envision a variety of other means and/or structures for performing the function and/or obtaining the results and/or one or more of the advantages described herein, and each of such variations and/or modifications is deemed to be within the scope of the embodiments described herein. Those skilled in the art will recognize or be able to ascertain using no more than routine experimentation many equivalents to the specific embodiments described herein. It is, therefore, to be understood that the foregoing embodiments are presented by way of example only and that, within the scope of the appended claims and equivalents thereto, inventive embodiments may be practiced otherwise than as specifically described. In addition, any combination of two or more features, systems, articles, materials, kits, and/or methods described herein, if such features, systems, articles, materials, kits, and/or methods are not mutually inconsistent, is included within the scope of the present disclosure. The above-described embodiments can be implemented in any of numerous ways. One or more aspects and embodiments of the present disclosure involving the performance of processes or methods may utilize program instructions executable by a device (e.g., a computer, a processor, or other device) to perform, or control performance of, the processes or methods. In this respect, various inventive concepts may be embodied as a computer readable storage medium (or multiple computer readable storage media) (e.g., a computer memory, one or more floppy discs, compact discs, optical discs, magnetic tapes, flash memories, circuit configurations in Field Programmable Gate Arrays or other semiconductor devices, or other tangible computer storage medium) encoded with one or more programs that, when executed on one or more computers or other processors, perform methods that implement one or more of the various embodiments described above. The computer readable medium or media can be transportable, such that the program or programs stored thereon can be loaded onto one or more different computers or other processors to implement various ones of the aspects described above. In some embodiments, computer readable media may be non-transitory media. The terms “program” or “software” are used herein in a generic sense to refer to any type of computer code or set of computer-executable instructions that can be employed to program a computer or other processor to implement various aspects as described above. Additionally, it should be appreciated that according to one aspect, one or more computer programs that when executed perform methods of the present disclosure need not reside on a single computer or processor, but may be distributed in a modular fashion among a number of different computers or processors to implement various aspects of the present disclosure. Computer-executable instructions may be in many forms, such as program modules, executed by one or more computers or other devices. Generally, program modules include routines, programs, objects, components, data structures, etc. that perform particular tasks or implement particular abstract data types. Typically the functionality of the program modules may be combined or distributed as desired in various embodiments. Also, data structures may be stored in computer-readable media in any suitable form. For simplicity of illustration, data structures may be shown to have fields that are related through location in the data structure. Such relationships may likewise be achieved by assigning storage for the fields with locations in a computer-readable medium that convey relationship between the fields. However, any suitable mechanism may be used to establish a relationship between information in fields of a data structure, including through the use of pointers, tags or other mechanisms that establish relationship between data elements. When implemented in software, the software code can be executed on any suitable processor or collection of processors, whether provided in a single computer or distributed among multiple computers. Further, it should be appreciated that a computer may be embodied in any of a number of forms, such as a rack-mounted computer, a desktop computer, a laptop computer, or a tablet computer, as non-limiting examples. Additionally, a computer may be embedded in a device not generally regarded as a computer but with suitable processing capabilities, including a Personal Digital Assistant (PDA), a smartphone, a tablet, or any other suitable portable or fixed electronic device. Also, a computer may have one or more input and output devices. These devices can be used, among other things, to present a user interface. Examples of output devices that can be used to provide a user interface include printers or display screens for visual presentation of output and speakers or other sound generating devices for audible presentation of output. Examples of input devices that can be used for a user interface include keyboards, and pointing devices, such as mice, touch pads, and digitizing tablets. As another example, a computer may receive input information through speech recognition or in other audible formats. Such computers may be interconnected by one or more networks in any suitable form, including a local area network or a wide area network, such as an enterprise network, and intelligent network (IN) or the Internet. Such networks may be based on any suitable technology and may operate according to any suitable protocol and may include wireless networks, wired networks or fiber optic networks. Also, as described, some aspects may be embodied as one or more methods. The acts performed as part of the method may be ordered in any suitable way. Accordingly, embodiments may be constructed in which acts are performed in an order different than illustrated, which may include performing some acts simultaneously, even though shown as sequential acts in illustrative embodiments. All definitions, as defined and used herein, should be understood to control over dictionary definitions, definitions in documents incorporated by reference, and/or ordinary meanings of the defined terms. The indefinite articles “a” and “an,” as used herein in the specification and in the claims, unless clearly indicated to the contrary, should be understood to mean “at least one.” The phrase “and/or,” as used herein in the specification and in the claims, should be understood to mean “either or both” of the elements so conjoined, i.e., elements that are conjunctively present in some cases and disjunctively present in other cases. Multiple elements listed with “and/or” should be construed in the same fashion, i.e., “one or more” of the elements so conjoined. Other elements may optionally be present other than the elements specifically identified by the “and/or” clause, whether related or unrelated to those elements specifically identified. Thus, as a non-limiting example, a reference to “A and/or B”, when used in conjunction with open-ended language such as “comprising” can refer, in one embodiment, to A only (optionally including elements other than B); in another embodiment, to B only (optionally including elements other than A); in yet another embodiment, to both A and B (optionally including other elements); etc. As used herein in the specification and in the claims, the phrase “at least one,” in reference to a list of one or more elements, should be understood to mean at least one element selected from any one or more of the elements in the list of elements, but not necessarily including at least one of each and every element specifically listed within the list of elements and not excluding any combinations of elements in the list of elements. This definition also allows that elements may optionally be present other than the elements specifically identified within the list of elements to which the phrase “at least one” refers, whether related or unrelated to those elements specifically identified. Thus, as a non-limiting example, “at least one of A and B” (or, equivalently, “at least one of A or B,” or, equivalently “at least one of A and/or B”) can refer, in one embodiment, to at least one, optionally including more than one, A, with no B present (and optionally including elements other than B); in another embodiment, to at least one, optionally including more than one, B, with no A present (and optionally including elements other than A); in yet another embodiment, to at least one, optionally including more than one, A, and at least one, optionally including more than one, B (and optionally including other elements); etc. In the claims, as well as in the specification above, all transitional phrases such as “comprising,” “including,” “carrying,” “having,” “containing,” “involving,” “holding,” “composed of,” and the like are to be understood to be open-ended, i.e., to mean including but not limited to. Only the transitional phrases “consisting of” and “consisting essentially of” shall be closed or semi-closed transitional phrases, respectively. The terms “approximately,” “substantially,” and “about” may be used to mean within ±20% of a target value in some embodiments, within ±10% of a target value in some embodiments, within ±5% of a target value in some embodiments, within ±2% of a target value in some embodiments. The terms “approximately,” “substantially,” and “about” may include the target value.