TIGIT was significantly up-regulated in advanced human gliomas.
When analyzing the expression of TIGIT depending on clinical and histological parameters available in TCGA database (Table 1), the results of non-parametric t-test showed that TIGIT was linked to sex, age, high grade and IDH wild type status (p = 0.0145; p = 0.0006; p < 0,0001; p = 0,0041 respectively) but not with the overall survival status (p = 0.1896).
According to the TCGA database, and depending on the 2016’ WHO classification of histological status (Fig. 1.A), TIGIT was highly expressed in GBM compared to astrocytoma (GII & GIII) (p < 0.0001), oligoastrocytoma (GII & GIII) (p < 0.0001) and oligodendroglioma (GII & GIII) (p = 0.0016). In contrast, no significant differences were detected when comparing the other subtypes (astrocytoma vs oligoastrocytoma; p = 0.6269), (astrocytoma vs oligodendroglioma; p = 0.1618) and (oligoastrocytoma vs oligodendroglioma p = 0.5209). These data indicated that high expression of TIGIT was positively associated with the most malignant subtypes of human gliomas.
ANOVA test revealed a significant difference between the means of the three available grades of Astrocytoma in TCGA database (GII, GIII and GIV/GBM) (p < 0.0001) (Table 1). Also, t-test was performed to evaluate the expression of TIGIT depending on the grades from TCGA database. Results showed higher expression in glioma Grade IV (GBM) compared to grade II (Astrocytoma, Oligoastrocytoma and Oligodendroglioma) (p < 0.0001) and grade III (Astrocytoma, Oligoastrocytoma and Oligodendroglioma) (p < 0.0001). In contrast, no statistical difference was observed between grade II and grade III (p = 0.1639) (Fig. 1B). These results corroborated those reported in (Fig. 1A) and showed that elevated levels of TIGIT were associated to the most advanced grade of human gliomas.
Table 1
Expression of TIGIT depending on distinct patients’ pathophysiological parameters in the TCGA cohort.
Clinical parameters
|
TCGA Cohorts (n)
|
Median of TIGIT mRNA
|
p value
|
Sex
|
Women (n = 285)
|
2.284
|
0.0145
|
Man (n = 382)
|
2.557
|
Age
|
< 40 (n = 262)
|
2.179
|
0.0006
|
> 40 (n = 393)
|
2.585
|
IDH mutation
|
Wild Type (n = 169)
|
2.972
|
0.0041
|
Mutant (n = 99)
|
2.362
|
Overall Survival Status
|
Living (n = 441)
|
2.337
|
0.1896
|
Deceased (n = 225)
|
2.585
|
Grade
|
Low grade (GII-GII) (n = 514)
|
2.221
|
< 0.0001
|
High grade (GIV) (n = 152)
|
3.000
|
Expression Subtype (GBM)
|
Mesenchymal (n = 49)
|
3.700
|
0,0101
|
Classical (n = 39)
|
2.585
|
Neural (n = 26)
|
2.989
|
Proneural (n = 29)
|
2.322
|
G-CIMP (n = 8)
|
2.953
|
Histological Diagnosis
|
Astrocytoma GII (n = 63)
|
1.958
|
< 0,0001
|
Astrocytoma GIII (n = 131)
|
2.312
|
Astrocytoma GIV (GBM) (n = 152)
|
3.170
|
Oligoastrocytoma GII (n = 74)
|
2.350
|
0.1223
|
Oligoastrocytoma GIII (n = 55)
|
2.219
|
Oligodendroglioma GII (n = 112)
|
2.274
|
0.9708
|
Oligodendroglioma GIII (n = 79)
|
2.323
|
*GBM: Glioblastoma multiform
*p value < 0.05 indicates a significant statistical test
The second cohort used in this study was related to the Moroccan patients, and which was composed of 53 human glioma patients (Table 2). In this cohort, TIGIT showed similar trends as those observed in the TCGA cohort. However, TIGIT expression did not reach significant difference when comparing groups of patients depending on age and sex (p = 0.8734; p = 0,1268 respectively).
Table 2
Expression of TIGIT depending on distinct patients’ pathophysiological parameters in the Moroccan cohort.
Clinical data
|
Cohorts (n)
|
Median of TIGIT mRNA
|
p value
|
Sex
|
Female (n = 23)
|
1.480
|
0.1268
|
Male (n = 30)
|
1.955
|
Age
|
< 40 (n = 28)
|
1.705
|
0,8734
|
> 40 (n = 24)
|
1.820
|
Grade
|
Low grade* (n = 24)
|
1.310
|
0.0423
|
High grade* (n = 29)
|
2.060
|
*Low grade (Grade I-Grade II)
*High grade (Grade III-Grade IV)
*p value < 0.05 indicates a significant statistical test
Interestingly, when exploring the 53 Moroccan glioma biopsies (Fig. 1C), TIGIT expression was significantly higher in advanced grades (GIII & GIV, n = 29) compared to lower grades of gliomas (GI & GII, n = 24) (p = 0.0423), confirming data from the TCGA cohort.
TIGIT expression was higher in IDH-wildtype glioma and was associated to mesenchymal-molecular subtype.
TIGIT expression was assessed according to the glioma molecular status. The expression of TIGIT was first evaluated depending on IDH gene status, IDH wildtype versus mutant in patients from TCGA database (Table 1). The results showed that IDH-wildtype patients displayed higher expressions of TIGIT compared to patients with IDH-mutant (p = 0.0041). This analysis could not be achieved in the second cohort (the Moroccan cohort) because of the small sample size (Wildtype, n = 3; Mutant, n = 4).
Subsequently, Kruskal-Wallis test was used to assess the difference in the expression of GBM molecular subtypes (Mesenchymal, Classical, Neural, Proneural and G-CIMP) and reported a significant difference between these groups (p = 0,0101) (Supplementary Fig. 1). Furthermore, when performing several t-tests between the five subtypes above, we found that TIGIT was significantly higher in the mesenchymal molecular subtype compared to Proneural and Classical (p = 0,0073; p = 0,0043 respectively) (Supplementary table 1).
Higher TIGIT expression was linked to higher expression of immune T cell markers and was positively correlated with FoxP3.
As previously mentioned in the introduction, several studies confirmed a higher expression of TIGIT on T-cells in advanced grades of different tumor sites. While exploring TCGA database, we wanted to assess whether a high expression of TIGIT could be associated with a higher infiltration of CD4 + T-cells, CD8 + T-cells and regulatory T-cells (via FoxP3) in glioma (Fig. 2). Our results showed that indeed, markers of CD4, CD8 and regulatory T-cells were upregulated in patients with high expression of TIGIT versus cases with lower expression of TIGIT in human glioma (p < 0.0001). Furthermore, we evaluated the correlation of TIGIT with the expression levels of the three cell markers (Table 3), and found that high expression of TIGIT was positively correlated with FoxP3, CD8 and CD4 expression (r = 0.4269, p < 0.0001; r = 0.2565, p < 0.0001; r = 0.1926, p = 0.0004). This suggested that high expression of TIGIT, whose expression was found to be associated with advanced subtypes of glioma from previous results, also was positively correlated with high infiltration of T cells.
Table 3: Correlation between the expression of T-cells markers with high TIGIT expression profile depending on TCGA database.
|
High TIGIT
|
r
|
p
|
CD4 T-cells
|
0.1926
|
0.0004
|
CD8 T-cells
|
0.2565
|
< 0,0001
|
FoxP3 Tregs
|
0.4269
|
< 0,0001
|
*p value < 0.05 indicates a significant statistical test
Subsequently, we wanted to check if the data from the Moroccan cohort were similar to those observed with the TCGA database. 18 patients having available mRNA expression data for both TIGIT and Foxp3 were analyzed to evaluate their correlation in the context of human glioma (Fig. 3). The result joined those of Table 3 and corroborated that TIGIT was positively correlated to FoxP3 (r = 0.534; p = 0.0225).
TIGIT was associated with an immunosuppressive microenvironment.
Since the regulatory T cell marker (FoxP3) exhibited the highest rate of correlation with TIGIT, we checked whether higher TIGIT expression would correlate with an immunosuppressive microenvironment in human glioma, globally. We therefore analyzed the level of expression of other known and potent regulatory T cell secreting immunosuppressive cytokines (Fig. 4), and found that the expression of both IL-10 and TGF-beta was upregulated in the case of high TIGIT compared to low TIGIT conditions (p < 0.0001 and p < 0.0001, respectively). This suggested that TIGIT highly contributes to the setting up of a potent immunosuppressive microenvironment in human glioma.
TIGIT was significantly upregulated along with other immune checkpoints in human glioma.
To assess the immunosuppressive status of human gliomas according to TIGIT expression level, three immune checkpoints (PD-1, VISTA, Tim-3) were assessed. Interestingly, and according to TCGA database (Fig. 5), we found that in the context of low TIGIT, the immune checkpoints PD-1, VISTA and Tim-3 exhibited lower expression compared to the context of high TIGIT expression where all the three immune checkpoints were upregulated (p < 0.0001). This showed that in human glioma’s microenvironment, the expression level of the three inhibitory molecules PD-1, VISTA and Tim-3 goes along with the level of expression of TIGIT, supporting the important role of TIGIT as a potential potent immunosuppressive immune checkpoint in human glioma.