The mRNA Expression Levels of TIMP2 in Human Tumors
We first compared mRNA expression level of TIMP2 in different cancers and normal tissues in the TCGA dataset using TIMER database. The result revealed that it expressed at relatively higher levels in Cholangio carcinoma, Head and Neck squamous cell carcinoma, Esophageal carcinoma, Kidney renal papillary cell carcinoma, and low expression in Bladder Urothelial Carcinoma, Breast invasive carcinoma, Colon adenocarcinoma, Kidney Chromophobe, Lung cancer, Prostate adenocarcinoma and Rectum adenocarcinoma (Figure 1A). Similar results were found in the TNMPlot and GEPIA database (Figure 1B and C).
The Correlation Between TIMP2 Expression and Cancer Patient’s Prognosis
We used the PrognoScan database to analyze whether TIMP2 expression was associated with prognosis in cancer patients. TIMP2 expression significantly affects prognosis in 5 cancer types, including colorectal, breast, head and neck, bladder, and ovarian cancers (Figure 2A-J). Two cohorts (GSE17536 and GSE17537) (26) including 177 and 55 colorectal cancer samples with different stages were analyzed, and the results showed that TIMP2 expression was correlated with prognosis (Figure 2E-H). To further investigate the prognostic potential of TIMP2 in different cancers, we evaluated the prognostic value of TIMP2 using the Kaplan-Meier plotter database (Figure 2K-T). The results indicated that high expression of TIMP2 is associated with poor prognosis in gastric cancer (overall survival [OS] hazard ratio [HR] = 1.38, 95% confidence interval [CI]: [1.16–1.63]; P=0.0002) and ovarian cancers (OS HR=1.23, 95% CI [1.08-1.4]; P=0.0013) (Figure 2M and S), and low expression of TIMP2 is association with poor prognosis in lung cancer (OS HR=0.83, 95% CI [0.73-0.94]; P=0.0034) (Figure 2Q). However, TIMP2 expression has less influence on breast cancer and liver cancer (Figures 2K and O). These results suggest the TIMP2 expression has an impact on the prognosis of ovarian, lung and gastric cancer.
The Expression of TIMP2 in Gastric and Colon Cancer Patients
To determine TIMP2 expression in gastric and colon cancer, we analyzed its mRNA expression using the TNMplot database. The results indicated that TIMP2 expression was decreased in colon cancer (Figures 3 A-C) and increased in gastric cancer (Figures 3 D-F). In order to verify the protein level of TIMP2, western blotting was used to detect the protein level in gastric cancer tissues and adjacent tissues of 5 patients with gastric cancer, and the results showed that the expression of TIMP2 in cancer tissues was significantly higher than that in adjacent tissues (Figures 4 A). The expression of TIMP2 protein was next analyzed by immunohistochemical staining in 50 paired gastric cancer tissues and corresponding normal tissues. Similar results were obtained from the assessment of TIMP2 protein expression (Figures 4 B). These results indicate that TIMP2 expression is decreased in colon cancer tissues and its mRNA and protein levels are increased in gastric cancer tissues.
TIMP2 Expression Correlated with Clinicopathological Parameters for Colon and Gastric Cancer Patients
In order to explore the relationship between TIMP2 expression and the clinicopathology of colon and gastric cancer, we first analyzed the correlation between TIMP2 expression and clinicopathological features of colon cancer in the TCGA database. As shown in Table 1, the low expression of TIMP2 is associated with the depth of tumor invasion, lymphatic metastasis, and tumor stage (P=0.037, P=0.001 and P=0.049, respectively). Furthermore, we investigated the correlation between TIMP2 expression and clinicopathological factors of the 50 paired gastric cancer and adjacent normal samples.
The results showed that TIMP2 expression was significantly correlated with Histologic grade in gastric cancer patients (P=0.007, Supplementary Table 1). In order to better understand the correlation of TIMP2 expression in cancer and the potential mechanism, we used Kaplan-Meier plotter database to investigated the relationship between TIMP2 expression and clinical characteristics in patients with gastric cancer. In stage 3 and 4 gastric cancer patients, high TIMP2 expression was associated with poorer OS and PFS, but not with OS and PFS in stage 1 (OS HR=1.96, P= 0.29; PFS HR=0.43, P= 0.19) and stage 2 (OS HR=1.59, P= 0.12; PFS HR=1.47, P= 0.2) and stage N0 patients (OS HR = 1.6, P= 0.35; PFS HR= 1.56, P= 0.38) of gastric cancer patients (Table 2). These results suggest that the expression level of TIMP2 may affect the prognosis of gastric cancer patients with lymph node metastasis.
The Expression of TIMP2 is Associated with Immune Infiltration in Gastric and Colon Cancer
Tumor infiltrate lymphocytes as an independent predictor of sentinel lymph node status and survival for cancer(27-29). Furthermore, we analyzed the relationship between TIMP2 expression and tumor-infiltrating immune cells in gastric cancer and colon cancer in the TIMER database. There was a significant positive correlation between TIMP2 expression and different type of immune cells, including CD4+ T cells (P = 3.62e-30, cor = 0.527), CD8+ T cells (P = 1.00e-09, cor = 0.297), B cells (P = 7.75e-03, cor =0.132), macrophages (P = 1.09e-52, cor = 0.664), Neutrophils (P = 1.18e-36, cor = 0.575), and dendritic cells (P = 1.83e-44, cor = 0.622) in COAD (Figure 5 A). Similarly, there were positive correlations with infiltrating levels of CD4+ T cells (P = 7.50e-13, cor = 0.363), CD8+ T cells (P = 5.45e-04, cor = 0.179), macrophages (P = 6.62e-53, cor = 0.687), Neutrophils (P = 1.99e-05, cor = 0.288), and dendritic cells (P = 1.67e-18, cor = 0.434) in STAD (Figure 5 A). In addition, TIMP2 expression has negative correlations with tumor purity in COAD (P = 4.60e-20, cor = -0.433) and STAD (P = 2.85e-03, cor = -0.153). Next, we analyzed the effect of TIMP2 expression and tumor-infiltrating immune cells on patient survival in COAD and STAD. As shown in Figure 5 B, the immune infiltration of tumor cells had no significant effect on the survival rate of COAD patients. However, high expression of TIMP2 in CD4+ T cells, macrophages and Neutrophils lead to a poor prognosis in STAD patients (P=0.035, P=0.002 and P=0.046, respectively). In addition, we found that several immune cell infiltration levels appeared to be associated with TIMP2 copy numbers changes, including B cells, CD8+T cells, neutrophils, and dendritic cells, especially in STAD (Figure 5 C). These results indicate that TIMP2 expression plays an important role in tumor immune cell infiltration.
Correlation Between TIMP2 Expression and Immune Markers
To further investigate the potential relationship between TIMP2 expression and immune infiltrating cells in STAD and COAD, we analyzed the correlations between TIMP2 and immune markers of various immune cells in the TIMER and GEPIA databases. These immune cell markers mainly include CD8+ T cells, T cells (general), B cells, monocytes, TAMs, M1 and M2 macrophages, neutrophils, NK cells and Dendritic cell. In addition, we also analyzed T cells with different functions, such as Th1 cells, Th2 cells, Tfh cells, Th17 cells, Treg cells and exhausted T cells. After adjusting for correlation by purity, the results showed that the expression level of TIMP2 was significantly correlated with the majority of immune marker sets in various immune cells and different T cells in STAD and COAD (Table 3).
According to the results in Table 3, the expression levels of most monocyte markers, TAMs and M2 macrophages and Dendritic cell are strongly correlated with the expression of TIMP2 in STAD and COAD. Specifically, we showed CD86, CSF1R of Monocyte, CCL2, CD68, IL10 of TAMs, CD163, VSIG4 and MS4A4A of M2 phenotype, HLA-DPB1, HLA-DQB1, HLA-DRA, HLA-DPA1, CD1C, NRP1 and ITGAX of Dendritic cell are significantly correlate with TIMP2 expression in STAD and COAD (P <0.001, Figure 6 A and B). Then, we further analyzed the correlation between TIMP2 expression and the above mononuclear markers and TAMs in STAD and COAD by GEPIA database. In COAD, the correlation between TIMP2 and mononuclear markers and TAM is similar to that of TIMER. In STAD, the correlation between TIMP2 and Dendritic cell are slightly different from those of TIMER (Supplementary Table 2). These results suggest that TIMP2 may regulate macrophage polarization, but there is no close relationship between TIMP2 and Dendritic cell infiltration in STAD and COAD.