Skip to main content
Download PDF

Correlate tumor mutation burden with immune signatures in human cancers

BMC Immunology volume 20, Article number: 4 (2019) Cite this article

Abstract

Background

Tumor mutation burden (TMB) has been associated with cancer immunotherapeutic response and cancer prognosis. Although many explorations have revealed that high TMB may yield many neoantigens to incite antitumor immune response, a systematic exploration of the correlation between TMB and immune signatures in different cancer types is lacking.

Results

We classified cancer into the lower-TMB subtype and the higher-TMB subtype for each of 32 cancer types based on their somatic mutation data from the Cancer Genome Atlas (TCGA), and compared the expression levels of immune-related genes and gene-sets between both subtypes of cancers in each cancer type. In some cancer types most of the immune signatures analyzed were upregulated in the lower-TMB subtype, while in some other cancer types the immune signatures were prone to be upregulated in the higher-TMB subtype. However, the regulatory T cells, immune cell infiltrate, tumor-infiltrating lymphocytes, and cytokine signatures tended to be upregulated in the lower-TMB subtype, and the cancer-testis antigen (CTA) and pro-inflammatory signatures were inclined to be upregulated in the higher-TMB subtype. Importantly, high TMB was associated with elevated expression of PD-L1 in diverse prevailing cancers. Furthermore, we found that higher TMB was associated with better survival prognosis in numerous cancer types while was associated with worse prognosis in a few cancer types.

Conclusions

High TMB may inhibit immune cell infiltrations while promote CTAs expression and inflammatory response in cancer. In many common cancer types, higher TMB may respond favorably to anti-PD-1/PD-L1 immunotherapy. Our data implicate that higher-TMB patients could gain a more favorable prognosis in diverse cancer types if treated with immunotherapy, otherwise would have a poorer prognosis compared to lower-TMB patients.

Background

Cancer immunotherapy is becoming increasingly noteworthy for its effectiveness in treating advanced and refractory cancers [1]. Particularly, the immune checkpoint blockade is being clinically used for treating diverse malignancies, such as melanoma [2, 3] and lung cancer [4]; the chimeric antigen receptor T cell therapy has been successfully utilized to treat refractory leukemia and lymphoma [5]. Nevertheless, these immunotherapies are beneficial to only 20% of cancer patients [6]. Thus, many efforts have been devoted to discovering the molecular determinants of immunotherapeutic responsiveness [7]. Some well-recognized molecular determinants include PD-L1 expression on tumor [8], DNA mismatch-repair deficiency [9], neoantigen load [10], and tumor-infiltrating lymphocytes [11]. Besides, to improve the efficacy of cancer immunotherapy, the combination of different immunotherapeutic methods [12, 13], or the combination of immunotherapy with other therapeutic approaches [14,15,16,17] have been explored.

A number of studies have explored the association between tumor mutation burden (TMB) and immunotherapy response [2,3,4, 18, 19]. These studies demonstrated that higher nonsynonymous mutation burden in tumors is inclined to form more neoantigens that make tumors to have higher immunogenicity, and thus result to improved clinical response to immunotherapy [4]. Nevertheless, a systematic exploration of the correlation between TMB and tumor immune activities in different cancer types remains lacking. To explore the association of TMB with tumor immunity in different cancer types, we compared the expression levels of immune-related genes and gene-sets between the lower-TMB subtype and the higher-TMB subtype of 32 cancer types based on the Cancer Genome Atlas (TCGA) data (https://portal.gdc.cancer.gov/). We tried to address several questions, including: Is the immune activity of the higher-TMB subtype different from that of the lower-TMB subtype of cancers? Are there any immune-related genes or gene-sets which are differentially expressed between the lower-TMB subtype and the higher-TMB subtype of cancers and whose expression is associated with clinical outcomes in cancer? Is the TMB itself associated with clinical outcomes in cancer?

Results

Association of TMB with regulatory T cell marker genes expression in human cancers

Regulatory T (Treg) cells play an important role in the maintenance of tumor immunosuppression [20]. We compared expression levels of 70 tumor-infiltrating Treg gene signatures [21] between the lower-TMB and the higher-TMB cancers in each of the 32 cancer types. We found that seven genes (ADPRH, IL1R1, KSR1, SOCS2, JAK1, NFAT5, and SSH1) were more highly expressed in the lower-TMB subtype than in the higher-TMB subtype of more than 10 cancer types (Additional file 1: Table S1). Of note, ADPRH had significantly higher expression levels in the lower-TMB subtype of 14 cancer types. The expression levels of the Treg gene-set were significantly higher in the lower-TMB subtype of 12 cancer types (HNSC, STAD, CHOL, UVM, PRAD, ACC, THCA, LUSC, ESCA, DLBC, KIRP, and LIHC) while were significantly higher in the higher-TMB subtype of 1 cancer type (THYM) (Wilcoxon rank-sum test, P < 0.05) (Fig. 1a). Interestingly, 27 Treg genes were more highly expressed in lower-TMB LIHC versus 1 more highly expressed in higher-TMB LIHC (Fisher’s exact test, P = 1.2*10− 8, OR = 42.42). In contrast, 14 Treg genes were more highly expressed in lower-TMB THYM versus 28 more highly expressed in higher-TMB THYM (Fisher’s exact test, P = 0.016, OR = 0.38). These results suggest that the relatedness between TMB and Treg cells infiltration degree depends on cancer types, whereas the lower-TMB subtype is likely to have stronger Treg cells infiltration than the higher-TMB subtype in diverse cancers.

Fig. 1
figure 1

Comparison of the expression levels of Treg, immune checkpoint, and immune cell infiltrate gene signatures between the lower-TMB and the higher-TMB subtypes of cancers. a The cancer types in which the Treg gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P < 0.05). b The cancer types in which the immune checkpoint gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P < 0.05). c Heat-map for the expression levels of immune cell infiltrate genes in lower-TMB and higher-TMB LIHC. d The cancer types in which the immune cell infiltrate (ICI) gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P-value< 0.05)

Full size image

Association of TMB with immune checkpoint genes expression in human cancers

Immune checkpoint molecules are important for tumor immune evasion [22]. We compared expression levels of 47 immune checkpoint genes [21] between the lower-TMB subtype and the higher-TMB subtype of cancers. We found that 12 genes had significantly higher expression levels in the higher-TMB subtype than in the lower-TMB subtype of at least 6 cancer types (Additional file 1: Table S2). The 12 genes included LAG3, CD80, TNFSF9, IDO1, CD70, KIR3DL1, CTLA4, PD-1, PD-L1, PD-L2, TIGIT, and TNFRSF9. Notably, LAG3 had higher expression levels in the higher-TMB subtype of 10 cancer types versus 2 cancer types of which LAG3 showed higher expression levels in the lower-TMB subtype. Interestingly, many immune checkpoint genes, which are established or promising targets for immune checkpoint blockade therapy, had significantly higher expression levels in the higher-TMB subtype of various cancers, such as CTLA4, PD-1, PD-L1, PD-L2, LAG3, IDO1 and TIGIT. In contrast, 16 immune checkpoint genes (C10orf54, CD200, CD40LG, ADORA2A, TNFSF14, BTLA, CD160, CD44, CD48, CD28, VTCN1, CD200R1, NRP1, TMIGD2, ICOS, and TNFSF15) had significantly higher expression levels in the lower-TMB subtype than in the higher-TMB subtype of at least 6 cancer types.

The expression levels of the immune checkpoint gene-set were significantly higher in the lower-TMB subtype than in the higher-TMB subtype of TGCT, KIRC, HNSC, ACC, THCA, LIHC, and THYM, while were significantly higher in the higher-TMB subtype of CESC, COAD, UCEC, and BLCA (Wilcoxon rank-sum test, P < 0.05) (Fig. 1b). Of the 47 immune checkpoint genes, 22 were more highly expressed in higher-TMB CESC versus 3 more highly expressed in lower-TMB CESC (Fisher’s exact test, P = 1.3*10− 5, OR = 12.55). In addition, 20 immune checkpoint genes were more highly expressed in higher-TMB COAD versus 1 more highly expressed in lower-TMB COAD (Fisher’s exact test, P = 1.98*10− 6, OR = 33). In contrast, no any immune checkpoint gene was more highly expressed in higher-TMB LIHC versus 22 more highly expressed in lower-TMB LIHC (Fisher’s exact test, P = 1.88*10− 8), and the similar result was observed in THCA. These data indicate that the association between TMB and the immune checkpoint activity is cancer type dependent, with in some cancers TMB being positively correlated with the immune checkpoint activity while in some other cancers they showing an inverse correlation.

Association of TMB with immune cell infiltration in human cancers

We compared the infiltration densities of 16 different immune cell subpopulations [23] between the lower-TMB subtype and the higher-TMB subtype of cancers. We found that 11 immune cell subpopulation marker genes had significantly higher expression levels in the lower-TMB subtype than in the higher-TMB subtype of at least 6 cancer types (Additional file 1: Table S3). The 11 genes included ENG (blood vessels), CD45RO (memory T cells), CD20 (B cells), CD1A (immature Dendritic Cells (iDCs)), TPSB2 (mast cells), PDPN (lymph vessels), CXCR5 (follicular helper T (Tfh cells)), IL3RA (plasmacytoid Dendritic Cells (pDCs)), CSF3 (neutrophils), CD3 (T cells), and CD57 (activated T or NK cells). In contrast, only CD68 (macrophages) had significantly higher expression levels in the higher-TMB subtype than in the lower-TMB subtype of at least 6 cancer types. It suggests that most of these immune cells have stronger infiltration in the lower-TMB subtype than in the higher-TMB subtype of cancers. Typically, 10 of the 16 immune cell subpopulation marker genes were more highly expressed in lower-TMB LIHC, but no one was more highly expressed in higher-TMB LIHC (Fig. 1c), indicating that heavy mutation load may inhibit immune cell infiltration in LIHC. Furthermore, we found 12 cancer types in which the expression levels of the immune cell subpopulation gene-set being significantly higher in the lower-TMB subtype, versus 1 cancer type in which the expression levels of this gene-set being significantly higher in the higher-TMB subtype (Fig. 1d). Again, this suggests that high TMB tends to inhibit immune cell infiltration in cancer.

Association of TMB with tumor-infiltrating lymphocytes (TILs) infiltration in human cancers

We compared expression levels of 120 TILs gene signatures [24] between the lower-TMB subtype and the higher-TMB subtype of cancers. We found that 90 genes had significantly higher expression levels in the lower-TMB subtype of at least 6 cancer types versus 5 having significantly higher expression levels in the higher-TMB subtype of at least 6 cancer types (Fisher’s exact test, P < 2.2*10− 16, OR = 67.11) (Additional file 1: Table S4). Notably, GIMAP6 had significantly higher expression levels in the lower-TMB subtype than in the higher-TMB subtype of 15 cancer types. GIMAP6 encodes the immunity-associated nucleotide 6 protein, a member of the GTPases of immunity-associated proteins family. This gene has been shown to be downregulated in several cancer types such as NSCLC [25] and LIHC [26]. Our results showed that this gene was more lowly expressed in higher-TMB LIHC than in lower-TMB LIHC. The expression levels of the TILs gene-set were significantly higher in the lower-TMB subtype of 9 cancer types (TGCT, KIRC, DLBC, HNSC, STAD, ACC, THCA, LIHC and THYM), while were significantly higher in the higher-TMB subtype of 2 cancer types (CESC and UCEC) (Fig. 2a). These results indicated that although the association between TMB and TILs infiltration was cancer type dependent, high TMB tended to inhibit TILs infiltration in various cancer types. Strikingly, 99 of the 120 TILs genes were more highly expressed in lower-TMB THYM compared to zero showing higher expression levels in higher-TMB THYM (Fisher’s exact test, P < 2.2*10− 16) (Fig. 2b). In addition, 88 TILs genes were more highly expressed in lower-TMB LIHC versus 3 more highly expressed in higher-TMB LIHC (Fisher’s exact test, P < 2.2*10− 16, OR = 104.46) (Fig. 2c).

Fig. 2
figure 2

Comparison of the expression levels of tumor-infiltrating lymphocytes (TILs) genes between the lower-TMB and the higher-TMB subtypes of cancers. a The cancer types in which the TILs gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P-value< 0.05). b Heat-map for the expression levels of TILs genes in lower-TMB and higher-TMB THYM. c Heat-map for the expression levels of TILs genes in lower-TMB and higher-TMB LIHC

Full size image

Association of TMB with cancer-testis antigen genes expression in human cancers

Cancer-testis antigens (CTAs) are a group of immunogenic proteins that are aberrantly activated in a variety of cancer types, and thus are important targets for developing cancer immunotherapy [27]. We compared expression levels of 276 CTA genes [28] between the lower-TMB subtype and the higher-TMB subtype of cancers. We found that 28 CTA genes had significantly higher expression levels in the lower-TMB subtype of at least 6 cancer types, while 63 had significantly higher expression levels in the higher-TMB subtype of at least 6 cancer types (Fisher’s exact test, P = 8.23*10− 5, OR = 0.38) (Additional file 1: Table S5). Of note, 6 CTA genes CEP55, KIF2C, TTK, OIP5, CASC5, and NUF2 had higher expression levels in the higher-TMB subtype of at least 15 cancer types, while had higher expression levels in the lower-TMB subtype of at most 3 cancer types. Interestingly, a number of genes encoding CTAs that are potentially useful for developing cancer vaccines were in the list of 63 CTA genes with higher expression levels in the higher-TMB subtype, such as MAGEA (MAGEA-1, 2, 3, 4, 6, 8, 9B, 10, 11, 12), NY-ESO-1, and PRAME.

The expression levels of the CTA gene-set were significantly higher in the higher-TMB subtype of 13 cancer types, including BRCA, LUAD, LIHC, SKCM, CESC, BLCA, THYM, LAML, LGG, HNSC, LUSC, ACC, and SARC, and were significantly higher in the lower-TMB subtype of 4 cancer types (COAD, UVM, THCA, and UCEC) (Fig. 3a). In many cancer types, the number of CTA genes with higher expression levels in the higher-TMB subtype far exceeded that of CTA genes with higher expression levels in the lower-TMB subtype. For example, 111 CTA genes were more highly expressed in higher-TMB LUAD versus 20 more highly expressed in lower-TMB LUAD (Fisher’s exact test, P < 2.2*10− 16, OR = 8.58); 102 CTA genes were more highly expressed in higher-TMB BRCA versus 29 more highly expressed in lower-TMB BRCA (Fisher’s exact test, P = 1.83*10− 13, OR = 4.98). These results suggest that high TMB is associated with elevated expression of many CTAs in cancer.

Fig. 3
figure 3

Comparison of the expression levels of cancer-testis antigen (CTA) and HLA genes between the lower-TMB and the higher-TMB subtypes of cancers. a The cancer types in which the CTA gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P-value< 0.05). b The cancer types in which the HLA gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P-value< 0.05). c Heat-maps for the expression levels of HLA genes in lower-TMB and higher-TMB ACC and LUAD. d Heat-maps for the expression levels of HLA genes in lower-TMB and higher-TMB CESC and COAD

Full size image

Association of TMB with HLA genes expression in human cancers

HLA (human leukocyte antigen) plays an important role in regulating the immune system in humans [29]. We compared expression levels of 24 HLA genes (with expression values available in the TCGA RNA-Seq data) between the lower-TMB subtype and the higher-TMB subtype of cancers. We found that 6 genes (HLA-J, DOA, DOB, DPB1, DQA1, and DQB2) had significantly higher expression levels in the lower-TMB subtype of at least 5 cancer types, while no any HLA gene showed significantly higher expression levels in the higher-TMB subtype of at least 5 cancer types (Additional file 1: Table S6). The expression levels of the HLA gene-set were significantly higher in the lower-TMB subtype of LIHC, HNSC, ACC, PRAD, and LUAD, and were significantly higher in the higher-TMB subtype of CESC, COAD, LGG, and SKCM (Fig. 3b). In LIHC, HNSC, ACC, PRAD, and LUAD, there were much more HLA genes showing higher expression levels in the lower-TMB subtype than those showing higher expression levels in the higher-TMB subtype (10 versus 1, 5 versus 0, 14 versus 0, 10 versus 0, and 16 versus 0 for LIHC, HNSC, ACC, PRAD, and LUAD, respectively) (Fig. 3c). In contrast, in CESC, COAD, and SKCM, there were much more HLA genes with higher expression levels in the higher-TMB subtype than those with higher expression levels in the lower-TMB subtype (12 versus 0, 12 versus 0, and 6 versus 0 for CESC, COAD, and SKCM, respectively) (Fig. 3d). These results suggest that the association between TMB and HLA expression is cancer type dependent.

Association of TMB with cytokine-related genes expression in human cancers

Cytokines are important components of the tumor immune microenvironment (TIM) [30]. Of 261 cytokine and cytokine receptor (CCR) genes [31], 93 showed significantly higher expression levels in the lower-TMB subtype versus 31 showing significantly higher expression levels in the higher-TMB subtype of at least 6 cancer types (Fisher’s exact test, P = 1.77*10− 10, OR = 4.10) (Additional file 1: Table S7). Notably, TNFAIP8L3, CCL14, CX3CR1, CCL21, IL1R1, and IL33 had higher expression levels in the lower-TMB subtype of at least 13 cancer types, while had higher expression levels in the higher-TMB subtype of at most 2 cancer types. In contrast, ILF2 showed higher expression levels in the higher-TMB subtype of 13 cancer types, while showed higher expression levels in the lower-TMB subtype of 1 cancer type.

Interestingly, the expression levels of the CCR gene-set were significantly higher in the lower-TMB subtype of 12 cancer types (LUSC, DLBC, UVM, TGCT, PRAD, LUAD, KIRC, ACC, HNSC, THCA, STAD, and LIHC), while were significantly higher in the higher-TMB subtype of THYM solely (Fig. 4a). It suggests that high TMB may lead to depressed cytokine activity in diverse cancers. Indeed, 98 CCR genes were more highly expressed in lower-TMB LIHC versus 7 more highly expressed in higher-TMB LIHC (Fisher’s exact test, P < 2.2*10− 16, OR = 21.71). In addition, 86 CCR genes were more highly expressed in lower-TMB THCA versus 2 more highly expressed in higher-TMB THCA (Fisher’s exact test, P < 2.2*10− 16, OR = 63.35).

Fig. 4
figure 4

Comparison of the expression levels of cytokine-related and pro-inflammatory genes between the lower-TMB and the higher-TMB subtypes of cancers. a The cancer types in which the cytokine and cytokine receptor (CCR) gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P-value< 0.05). b The cancer types in which the pro-inflammatory gene-set is differentially expressed between the lower-TMB and the higher-TMB subtypes (Wilcoxon rank-sum test, P-value< 0.05). c Heat-map for the expression levels of pro-inflammatory genes in lower-TMB and higher-TMB THCA. d Heat-map for the expression levels of pro-inflammatory genes in lower-TMB and higher-TMB COAD

Full size image

Association of TMB with pro-inflammatory genes expression in human cancers

Inflammatory responses play important roles in regulating the TIM [32]. We compared expression levels of 15 pro-inflammatory genes [33] between the lower-TMB subtype and the higher-TMB subtype of cancers. We found that 10 pro-inflammatory genes were more highly expressed in the lower-TMB subtype versus 4 more highly expressed in the higher-TMB subtype of at least 5 cancer types (Additional file 1: Table S8). Notably, CXCL9, CXCL10, and IFNG had significantly higher expression levels in the higher-TMB subtype of at least 10 cancer types, while had significantly higher expression levels in the lower-TMB subtype of at most 4 cancer types. The expression levels of the pro-inflammatory gene-set were significantly higher in the lower-TMB subtype of 6 cancer types (KIRC, TGCT, ACC, HNSC, KIRP, and THCA), and were significantly higher in the higher-TMB subtype of 10 cancer types (CESC, COAD, UCEC, BLCA, LUAD, BRCA, STAD, LGG, SKCM, and SARC) (Fig. 4b). Interestingly, 12 of the 15 pro-inflammatory genes were more highly expressed in lower-TMB THCA versus zero more highly expressed in higher-TMB THCA (Fig. 4c). In contrast, 11 pro-inflammatory genes were more highly expressed in higher-TMB COAD versus zero more highly expressed in lower-TMB COAD (Fig. 4d). These results indicate that the association between TMB and pro-inflammatory activity is cancer type dependent, with in some cancer types high TMB enhancing pro-inflammatory activity while in some other cancer types high TMB inhibiting pro-inflammatory activity in cancer.

Associations among TMB, immune signatures, and cancer prognosis

To explore the association among TMB, immune signatures, and cancer prognosis, we compared survival prognosis (overall survival (OS) and disease free survival (DFS)) between immune gene-set higher-expression-level and lower-expression-level lower-TMB cancers and higher-TMB cancers, respectively. We found that some cancer types showed a significant correlation between immune gene-set expression and survival prognosis in the lower-TMB subtype but not in the higher-TMB subtype. For example, higher expression levels of the Treg, immune checkpoint, immune cell infiltrate, TILs, and CCR gene-sets were consistently associated with worse DFS in lower-TMB GBM (log-rank test, P < 0.05), but there was no any immune gene-set whose expression was associated with survival prognosis in higher-TMB GBM (Fig. 5a). In contrast, some cancer types showed a significant correlation between immune gene-set expression and survival prognosis in the higher-TMB subtype but not in the lower-TMB subtype. For example, higher expression levels of the immune checkpoint, immune cell infiltrate, TILs, CCR, HLA, and pro-inflammatory gene-sets were consistently associated with better OS and/or DFS in higher-TMB SARC, but there was no any immune gene-set whose expression was associated with survival prognosis in lower-TMB SARC (Fig. 5b). These results indicate that there exist significant associations among TMB, immune signatures, and survival prognosis in some cancer types.

Fig. 5
figure 5

Association among TMB, immune signatures, and cancer prognosis. a Kaplan-Meier survival curves show that the elevated expression of Treg, immune checkpoint, immune cell infiltrate, TILs, and CCR gene-sets are consistently associated with worse disease free survival (DFS) prognosis in lower-TMB GBM (log-rank test, P < 0.05). b The elevated expression of immune checkpoint, TILs, HLA, CCR, and pro-inflammatory gene-sets are consistently associated with better overall survival (OS) and/or DFS prognosis in higher-TMB SARC (log-rank test, P < 0.05)

Full size image

Discussion

We analyzed the associations between TMB and diverse immune signatures in 32 human cancer types. We found that these associations were generally cancer type dependent. For example, most of the immune signatures were upregulated in the lower-TMB subtype of HNSC, ACC, THCA, and LIHC, while were downregulated in the lower-TMB subtype of CESC relative to their higher-TMB subtypes. However, the Treg cells, immune cell infiltrate, TILs, and CCR signatures were inclined to be upregulated in the lower-TMB subtype of various cancer types, suggesting that high TMB may inhibit immune cell infiltration in the TIM. In contrast, the CTA and pro-inflammatory signatures tended to be upregulated in the higher-TMB subtype of various cancer types, suggesting that high TMB may promote CTA expression and tumor inflammatory response. Interestingly, HNSC, ACC, THCA, and LIHC were the cancer types in which TMB and the immune gene-set expression alteration had a significant correlation for almost all the immune signatures analyzed (Table 1). It implies that TMB could have a significant impact on the TIM in these cancer types. In fact, when we compared survival prognosis between the lower-TMB subtype and the higher-TMB subtype of cancers, we found that the lower-TMB subtype had better OS and/or DFS prognosis than the higher-TMB subtype in three of the four cancer types including HNSC, ACC, and LIHC (Fig. 6a). These data suggest that TMB is associated with survival prognosis in some cancer types, and that the mechanism underlying this association could lie in the marked differences in immune cell infiltration densities and immune activities between the lower-TMB subtype and the higher-TMB subtype of these cancers.

Table 1 Comparison of immune activities between the lower-TMB and the higher-TMB subtypes of cancers
Full size table
Table 2 32 TCGA cancer types used in this study
Full size table
Fig. 6
figure 6

Association of TMB with survival prognosis and PD-L1 expression in cancers. a Kaplan-Meier survival curves show that TMB is associated with survival prognosis in diverse cancer types (log-rank test, P < 0.05). b PD-L1 is differentially expressed between the lower-TMB and the higher-TMB subtypes of diverse cancer types (Student’s t test, P < 0.05, fold change > 1.5)

Full size image

Interestingly, we found that high TMB was associated with elevated pro-inflammatory immune activity while depressed immune cell infiltration in diverse cancers. These findings appear to be contradictory and disagree with the established notion that high TMB may yield numerous neoantigens that incite anti-tumor immune response [4]. The possible explanations are that high TMB is often associated with genome instability that may inhibit anti-tumor immune response [34], and that the increased pro-inflammatory immune activity could be attributed to the higher percent of tumor necrosis component elicited by gene mutations in the higher-TMB cancer.

A recent study explored the landscape of TMB across 100,000 cancer cases in over 100 tumor types, and revealed that there were a substantial number of high-TMB cases across nearly every cancer type [35]. It justifies that the comparisons of high-TMB with low-TMB cases within each cancer type in the current study are sensible, although the methods in measuring the high-TMB and low-TMB are distinct between both studies. In another recent study [19], Goodman et al. revealed that higher TMB was associated with better clinical outcome in diverse cancers treated with various immunotherapies. However, our analyses showed that higher TMB was associated with better survival (OS and/or DFS) prognosis in SKCM, OV, and BLCA, while was associated with worse prognosis in ACC, ESCA, THYM, LIHC, LGG, HNSC, KICH, and READ based on the TCGA data (Fig. 6a). Seemingly, the revealed associations between TMB and cancer prognosis exhibit discrepancy between both studies. The main reason behind the discrepancy could be that most of the TCGA patients were not treated with immunotherapy. Indeed, for the TCGA cases likely with immunotherapy such as melanoma, higher TMB was associated with better prognosis, consistent with the conclusion drawn in [19]. Our data, together with the data from [19], implicate that higher-TMB patients could gain a more favorable prognosis compared to lower-TMB patients in diverse cancers if treated with immunotherapy, otherwise higher-TMB patients would have an unfavorable prognosis compared to lower-TMB patients.

Previous studies have shown that PD-L1 expression is positively associated with response to anti-PD-L1 immunotherapy [36]. Our data showed that PD-L1 had significantly higher expression levels in the higher-TMB subtype than in the lower-TMB subtype of 7 cancer types (STAD, COAD, BLCA, SKCM, PAAD, LUAD, and CESC), compared to 3 cancer types (KIRP, GBM, and UVM) of which PD-L1 was more lowly expressed in the higher-TMB subtype (Student’s t test, P < 0.05, fold change > 1.5) (Fig. 6b). It indicates that in the diverse prevailing malignancies, higher TMB could respond favorably to anti-PD-1/PD-L1 immunotherapy, bolstering the conclusion drawn in [19].

In this study, we defined the higher-TMB and the lower-TMB subtypes of cancers on the basis of the TMB scores in the individual cancer types instead of across all cancer types that was a different approach from prior studies [4, 19, 35]. We used this approach mainly considering the intertumor heterogeneity. Indeed, there were marked gaps in TMB among different cancer types, with some cancers having high TMB such as SKCM, LUSC, LUAD, PAAD, ESCA, and BLCA while some other cancers having low TMB such as LAML, PCPG, TGCT, THCA, and UVM in general (Fig. 7). If we defined the higher-TMB and the lower-TMB subtypes of cancers based on the TMB scores across all cancer types, the sample size of the lower-TMB subtype in the highly-mutated cancer types, and the sample size of the higher-TMB subtype in the lowly-mutated cancer types would be too small to perform effective comparisons in these cancer types. As a result, it would be difficult to look into the correlations between TMB and immune activities in these cancer types. Interestingly, a recent pan-cancer analysis showed that the anti-tumor immune response could be more effective when the immune system robustly responded against a few antigens other than diversely responded against numerous different antigens [37]. It implicates that even in some cancer types with low TMB, anti-tumor immune response may take effect if the immune system heavily targets a few proper neoantigens.

Fig. 7
figure 7

Distribution of TMB scores across 32 cancer types

Full size image

There are several limitations in the present study. First, we have used a relatively loose significance level (P < 0.05) to identify differentially expressed genes (or gene-sets) when tens or hundreds of genes were tested. We did not use the method of adjusting for multiple tests such as the false discovery rate (FDR) [38] to define a more stringent significance threshold considering that the numbers of multiple tests were less than 100 in most of the cases, and were only eight in the comparisons of gene-set expression. Second, the cutoff for defining higher-TMB samples (the samples with TMB scores of upper quartile) and lower-TMB samples (the samples with TMB scores of lower quartile) was somewhat arbitrary. Nevertheless, when we reanalyzed the data using two different cutoffs to define higher-TMB samples versus lower-TMB samples (upper third vs. lower third, and upper half vs. lower half, respectively), we obtained the consistent results (Additional file 2: Table S9).

Conclusions

High TMB may inhibit immune cell infiltrations while promote CTAs expression and inflammatory response in cancer. High TMB was associated with elevated expression of PD-L1 in diverse prevailing cancers. Higher TMB was associated with better clinical outcomes in SKCM, OV, and BLCA, while was associated with worse prognosis in ACC, ESCA, THYM, LIHC, LGG, HNSC, KICH, and READ. Our data provide insights into the correlation between TMB and tumor immune response in different types of cancers, and have potential clinical implications for cancer immunotherapy.

Methods

Comparisons of expression levels of genes and gene-sets between two classes of samples

The TCGA RNA-Seq gene expression data (Level 3) were normalized by base-2 log transformation. Student’s t test was used to compare expression levels of a single gene between two groups of samples, and Wilcoxon rank-sum test was used to compare expression levels of a gene-set between two groups of samples.

Definition of the higher-TMB and the lower-TMB subtypes of cancers

Consistent with a method we proposed previously [39], the TMB score of a tumor sample was calculated as follows:

total number of truncating mutations * 2.0 + total number of non-truncating mutations * 1.0

Nonsense, frame-shift deletion or insertion, and splice-site mutations were included in the truncating mutation category, and missense, in-frame deletion or insertion, and nonstop mutations were included in the non-truncating mutation category. The higher-TMB subtype (samples with TMB scores higher than the third quartile value) and the lower-TMB subtype (samples with TMB scores lower than the first quartile value) were defined in each individual cancer type based on the TMB scores of its tumor samples.

Survival analyses

Besides the higher-TMB and the lower-TMB subtypes, we also defined two classes of cancer patients based on gene-set expression levels that indicate the activity of immune signatures. The expression level of a gene-set was evaluated as the average expression value of all the genes in the gene-set. The Kaplan-Meier method was used to compare survival (OS and DFS) between two classes of cancer patients (higher-TMB versus lower-TMB, and gene-set higher-expression-level (upper half) versus gene-set lower-expression-level (bottom half)). The log-rank test was used to evaluate the significance of survival-time differences between two classes of cancer patients.

Statistical and computational analyses

For all statistical tests, a two-tailed P < 0.05 was considered statistically significant. All the statistical and computational analyses were performed using R programming (version 3.2.3, https://www.r-project.org/).

Abbreviations

ACC:

Adrenocortical carcinoma

BLCA:

Bladder urothelial carcinoma

BRCA:

Breast invasive carcinoma

CCR:

Cytokine and cytokine receptor

CESC:

Cervical squamous-cell carcinoma and endocervicaladeno-carcinoma

CHOL:

Cholangiocarcinoma

COAD:

Colon adenocarcinoma

CTA:

Cancer-testis antigen

CTLA4:

Cytotoxic T-lymphocyte-associated protein 4

DFS:

Disease free survival

DLBC:

Lymphoid neoplasm diffuse large B-cell lymphoma

ESCA:

Esophageal carcinoma

FDR:

False discovery rate

GBM:

Glioblastoma multiforme

HLA:

Human leukocyte antigen

HNSC:

Head and neck squamous cell carcinoma

iDCs:

Immature Dendritic Cells

KICH:

Kidney chromophobe

KIRC:

Kidney renal clear cell carcinoma

KIRP:

Kidney renal papillary cell carcinoma

LAML:

Acute myeloid leukemia

LGG:

Brain lower-grade glioma

LIHC:

Liver hepatocellular carcinoma

LUAD:

Lung adenocarcinoma

LUSC:

Lung squamous cell carcinoma

NSCLC:

Non-small cell lung cancer

OR:

Odds ratio

OS:

Overall survival

OV:

Ovarian serous cystadenocarcinoma

PAAD:

Pancreatic adeno-carcinoma

PD-1:

Programmed cell death protein 1

pDCs:

Plasmacytoid Dendritic Cells

PD-L1:

Programmed cell death 1 ligand

PRAD:

Prostate adenocarcinoma

READ:

Rectum adenocarcinoma

SARC:

Sarcoma

SKCM:

Skincutaneous melanoma

STAD:

Stomach adenocarcinoma

TCGA:

The Cancer Genome Atlas

Tfh:

Follicular helper T

TGCT:

Testicular germ-cell tumors

THCA:

Thyroid carcinoma

THYM:

Thymoma

TILs:

Tumor-infiltrating lymphocytes

TIM:

Tumor immune microenvironment

TMB:

Tumor mutation burden

Treg:

Regulatory T

UCEC:

Uterine corpus endometrial carcinoma

UCS:

Uterine carcino-sarcoma

UVM:

Uveal melanoma

References

  1. Hoos A. Development of immuno-oncology drugs - from CTLA4 to PD1 to the next generations. Nat Rev Drug Discov. 2016;15(4):235–47.

    CAS  Article  Google Scholar 

  2. Van Allen EM, Miao D, Schilling B, Shukla SA, Blank C, Zimmer L, Sucker A, Hillen U, Geukes Foppen MH, Goldinger SM, et al. Genomic correlates of response to CTLA-4 blockade in metastatic melanoma. Science. 2015;350(6257):207–11.

    Article  Google Scholar 

  3. Snyder A, Makarov V, Merghoub T, Yuan J, Zaretsky JM, Desrichard A, Walsh LA, Postow MA, Wong P, Ho TS, et al. Genetic basis for clinical response to CTLA-4 blockade in melanoma. N Engl J Med. 2014;371(23):2189–99.

    Article  Google Scholar 

  4. Rizvi NA, Hellmann MD, Snyder A, Kvistborg P, Makarov V, Havel JJ, Lee W, Yuan J, Wong P, Ho TS, et al. Cancer immunology. Mutational landscape determines sensitivity to PD-1 blockade in non-small cell lung cancer. Science. 2015;348(6230):124–8.

    CAS  Article  Google Scholar 

  5. June CH, O'Connor RS, Kawalekar OU, Ghassemi S, Milone MC. CAR T cell immunotherapy for human cancer. Science. 2018;359(6382):1361–5.

    CAS  Article  Google Scholar 

  6. Braun DA, Burke KP, Van Allen EM. Genomic approaches to understanding response and resistance to immunotherapy. Clin Cancer Res : an official journal of the American Association for Cancer Research. 2016;22(23):5642–50.

    CAS  Article  Google Scholar 

  7. Gibney GT, Weiner LM, Atkins MB. Predictive biomarkers for checkpoint inhibitor-based immunotherapy. The Lancet Oncology. 2016;17(12):e542–51.

    CAS  Article  Google Scholar 

  8. Patel SP, Kurzrock R. PD-L1 expression as a predictive biomarker in Cancer immunotherapy. Mol Cancer Ther. 2015;14(4):847–56.

    CAS  Article  Google Scholar 

  9. Le DT, Uram JN, Wang H, Bartlett BR, Kemberling H, Eyring AD, Skora AD, Luber BS, Azad NS, Laheru D, et al. PD-1 blockade in tumors with mismatch-repair deficiency. N Engl J Med. 2015;372(26):2509–20.

    CAS  Article  Google Scholar 

  10. Efremova M, Finotello F, Rieder D, Trajanoski Z. Neoantigens generated by individual mutations and their role in Cancer immunity and immunotherapy. Front Immunol. 2017;8:1679.

    Article  Google Scholar 

  11. Zito Marino F, Ascierto PA, Rossi G, Staibano S, Montella M, Russo D, Alfano R, Morabito A, Botti G, Franco R. Are tumor-infiltrating lymphocytes protagonists or background actors in patient selection for cancer immunotherapy? Expert Opin Biol Ther. 2017;17(6):735–46.

    Article  Google Scholar 

  12. Postow MA, Chesney J, Pavlick AC, Robert C, Grossmann K, McDermott D, Linette GP, Meyer N, Giguere JK, Agarwala SS, et al. Nivolumab and ipilimumab versus ipilimumab in untreated melanoma. N Engl J Med. 2015;372(21):2006–17.

    Article  Google Scholar 

  13. Wang B, Zhang W, Jankovic V, Golubov J, Poon P, Oswald EM, Gurer C, Wei J, Ramos I, Wu Q, et al. Combination cancer immunotherapy targeting PD-1 and GITR can rescue CD8(+) T cell dysfunction and maintain memory phenotype. Science Immunol. 2018;3(29).

  14. Song W, Shen L, Wang Y, Liu Q, Goodwin TJ, Li J, Dorosheva O, Liu T, Liu R, Huang L. Synergistic and low adverse effect cancer immunotherapy by immunogenic chemotherapy and locally expressed PD-L1 trap. Nat Commun. 2018;9(1):2237.

    Article  Google Scholar 

  15. Bommareddy PK, Aspromonte S, Zloza A, Rabkin SD, Kaufman HL. MEK inhibition enhances oncolytic virus immunotherapy through increased tumor cell killing and T cell activation. Sci Transl Med. 2018;10(471).

  16. Goel S, DeCristo MJ, Watt AC, BrinJones H, Sceneay J, Li BB, Khan N, Ubellacker JM, Xie S, Metzger-Filho O, et al. CDK4/6 inhibition triggers anti-tumour immunity. Nature. 2017;548(7668):471–5.

    CAS  Article  Google Scholar 

  17. Li M, Liu Z, Wang X. Exploration of the combination of PLK1 inhibition with immunotherapy in cancer treatment. Journal of Oncology. 2018;2018(Article ID 3979527):13.

    Google Scholar 

  18. Hugo W, Zaretsky JM, Sun L, Song C, Moreno BH, Hu-Lieskovan S, Berent-Maoz B, Pang J, Chmielowski B, Cherry G, et al. Genomic and transcriptomic features of response to anti-PD-1 therapy in metastatic melanoma. Cell. 2016;165(1):35–44.

    CAS  Article  Google Scholar 

  19. Goodman AM, Kato S, Bazhenova L, Patel SP, Frampton GM, Miller V, Stephens PJ, Daniels GA, Kurzrock R. Tumor mutational burden as an independent predictor of response to immunotherapy in diverse cancers. Mol Cancer Ther. 2017;16(11):2598–608.

    CAS  Article  Google Scholar 

  20. Tanaka A, Sakaguchi S. Regulatory T cells in cancer immunotherapy. Cell Res. 2017;27(1):109–18.

    CAS  Article  Google Scholar 

  21. De Simone M, Arrigoni A, Rossetti G, Gruarin P, Ranzani V, Politano C, Bonnal RJP, Provasi E, Sarnicola ML, Panzeri I, et al. Transcriptional landscape of human tissue lymphocytes unveils uniqueness of tumor-infiltrating T regulatory cells. Immunity. 2016;45(5):1135–47.

    Article  Google Scholar 

  22. Topalian SL, Drake CG, Pardoll DM. Immune checkpoint blockade: a common denominator approach to cancer therapy. Cancer Cell. 2015;27(4):450–61.

    CAS  Article  Google Scholar 

  23. Bindea G, Mlecnik B, Tosolini M, Kirilovsky A, Waldner M, Obenauf AC, Angell H, Fredriksen T, Lafontaine L, Berger A, et al. Spatiotemporal dynamics of intratumoral immune cells reveal the immune landscape in human cancer. Immunity. 2013;39(4):782–95.

    CAS  Article  Google Scholar 

  24. Massink MPG, Kooi IE, Martens JWM, Waisfisz Q, Meijers-Heijboer H. Genomic profiling of CHEK2*1100delC-mutated breast carcinomas. BMC Cancer. 2015;15:877.

    Article  Google Scholar 

  25. Shiao YM, Chang YH, Liu YM, Li JC, Su JS, Liu KJ, Liu YF, Lin MW, Tsai SF. Dysregulation of GIMAP genes in non-small cell lung cancer. Lung Cancer. 2008;62(3):287–94.

    Article  Google Scholar 

  26. Huang Z, Zhang W, Gao C, Ji B, Chi X, Zheng W, Wang HL. Dysregulation of GTPase IMAP family members in hepatocellular cancer. Mol Med Rep. 2016;14(5):4119–23.

    CAS  Article  Google Scholar 

  27. Caballero OL, Chen Y-T. Cancer/testis (CT) antigens: potential targets for immunotherapy. Cancer Sci. 2009;100(11):2014–21.

    CAS  Article  Google Scholar 

  28. Almeida LG, Sakabe NJ, de Oliveira AR, MCC S, Mundstein AS, Cohen T, Chen Y-T, Chua R, Gurung S, Gnjatic S, et al. CTdatabase: a knowledge-base of high-throughput and curated data on cancer-testis antigens. Nucleic Acids Res. 2009;37(Database issue):D816–9.

    CAS  Article  Google Scholar 

  29. Munro A, Bright S. Products of the major histocompatibility complex and their relationship to the immune response. Nature. 1976;264(5582):145–52.

    CAS  Article  Google Scholar 

  30. Lippitz BE. Cytokine patterns in patients with cancer: a systematic review. Lancet Oncol. 2013;14(6):e218–28.

    CAS  Article  Google Scholar 

  31. Wong HS-C, Chang C-M, Liu X, Huang W-C, Chang W-C. Characterization of cytokinome landscape for clinical responses in human cancers. Oncoimmunology. 2016;5(11):e1214789.

    Article  Google Scholar 

  32. Danaher P, Warren S, Lu R, Samayoa J, Sullivan A, Pekker I, Wallden B, Marincola FM, Cesano A. Pan-cancer adaptive immune resistance as defined by the tumor inflammation signature (TIS): results from the Cancer genome atlas (TCGA). J Immunother Cancer. 2018;6(1):63.

    Article  Google Scholar 

  33. Bedognetti D, Hendrickx W, Marincola FM, Miller LD. Prognostic and predictive immune gene signatures in breast cancer. Curr Opin Oncol. 2015;27(6):433–44.

    CAS  Article  Google Scholar 

  34. Jiang Z, Liu Z, Li M, Chen C, Wang X. Immunogenomics analysis reveals that TP53 mutations inhibit tumor immunity in gastric Cancer. Transl Oncol. 2018;11(5):1171–87.

    Article  Google Scholar 

  35. Chalmers ZR, Connelly CF, Fabrizio D, Gay L, Ali SM, Ennis R, Schrock A, Campbell B, Shlien A, Chmielecki J, et al. Analysis of 100,000 human cancer genomes reveals the landscape of tumor mutational burden. Genome medicine. 2017;9(1):34.

    Article  Google Scholar 

  36. Herbst RS, Soria JC, Kowanetz M, Fine GD, Hamid O, Gordon MS, Sosman JA, McDermott DF, Powderly JD, Gettinger SN, et al. Predictive correlates of response to the anti-PD-L1 antibody MPDL3280A in cancer patients. Nature. 2014;515(7528):563–7.

    CAS  Article  Google Scholar 

  37. Thorsson V, Gibbs DL, Brown SD, Wolf D, Bortone DS, Ou Yang TH, Porta-Pardo E, Gao GF, Plaisier CL, Eddy JA, et al. The immune landscape of Cancer. Immunity. 2018;48(4):812–30 e814.

    CAS  Article  Google Scholar 

  38. Benjami Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Statist Soc B. 1995;57:289–300.

    Google Scholar 

  39. Liu Z, Li M, Jiang Z, Wang X. A comprehensive immunologic portrait of triple-negative breast Cancer. Transl Oncol. 2018;11(2):311–29.

    Article  Google Scholar 

Download references

Acknowledgements

We thank the anonymous reviewers for their valuable comments which greatly helped us to improve the manuscript. We thank Dr. Sen Zhou from China Pharmaceutical University for providing statistical help in this work.

Funding

This work was supported by the China Pharmaceutical University [Grant number 3150120001].

Availability of data and materials

We downloaded RNA-Seq gene expression profiles (Level 3), gene somatic mutations (Level 3), and clinical data for 32 cancer types from the TCGA data portal (https://portal.gdc.cancer.gov/). For survival analyses, we used clinical data from FireBrowse (http://gdac.broadinstitute.org/) for the TCGA data. The numbers of cancer samples and normal samples for each cancer type are listed in Table 2.

Author information

Authors and Affiliations

  1. Biomedical Informatics Research Lab, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 211198, China

    Xiaosheng Wang & Mengyuan Li

  2. Cancer Genomics Research Center, School of Basic Medicine and Clinical Pharmacy, China Pharmaceutical University, Nanjing, 211198, China

    Xiaosheng Wang & Mengyuan Li

  3. Big Data Research Institute, China Pharmaceutical University, Nanjing, 211198, China

    Xiaosheng Wang & Mengyuan Li

Authors
  1. Xiaosheng Wang
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Mengyuan Li
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

XW conceived of the research, performed the research, and wrote the manuscript. ML performed the research, and helped prepare for the manuscript. Both authors read and approved the final manuscript.

Corresponding author

Correspondence to Xiaosheng Wang.

Ethics declarations

Ethics approval and consent to participate

Ethical approval was waived since we used only publicly available data and materials in this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Comparison of regulatory T cell marker genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S2. Comparison of immune checkpoint genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S3. Comparison of immune cell infiltrate genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S4. Comparison of tumor-infiltrating lymphocytes genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S5. Comparison of cancer-testis antigen genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S6. Comparison of HLA genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S7. Comparison of cytokine and cytokine receptor genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. Table S8. Comparison of pro-inflammation genes expression levels between the lower-TMB and the higher-TMB subtypes of cancers. (XLSX 130 kb)

Additional file 2:

Table S9. Comparison of the results obtained by using different cutoffs to define the higher-TMB and the lower-TMB tumor samples. (XLSX 21 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Wang, X., Li, M. Correlate tumor mutation burden with immune signatures in human cancers. BMC Immunol 20, 4 (2019). https://doi.org/10.1186/s12865-018-0285-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12865-018-0285-5

Keywords

  • Tumor mutation burden
  • Immune signatures
  • Tumor immune microenvironment
  • Cancer immunotherapy
  • Cancer prognosis

BMC Immunology

ISSN: 1471-2172

Contact us