- Open Access
Comprehensive analysis of expression, prognosis and immune infiltration for TIMPs in glioblastoma
BMC Neurology volume 21, Article number: 447 (2021)
Tissue inhibitors of metalloproteinase (TIMP) family proteins are peptidases involved in extracellular matrix (ECM) degradation. Various diseases are related to TIMPs, and the primary reason is that TIMPs can indirectly regulate remodelling of the ECM and cell signalling by regulating matrix metalloproteinase (MMP) activity. However, the link between TIMPs and glioblastoma (GBM) is unclear.
This study aimed to explore the role of TIMP expression and immune infiltration in GBM.
Oncomine, GEPIA, OSgbm, LinkedOmics, STRING, GeneMANIA, Enrichr, and TIMER were used to conduct differential expression, prognosis, and immune infiltration analyses of TIMPs in GBM.
All members of the TIMP family had significantly higher expression levels in GBM. High TIMP3 expression correlated with better overall survival (OS) and disease-specific survival (DSS) in GBM patients. TIMP4 was associated with a long OS in GBM patients. We found a positive relationship between TIMP3 and TIMP4, identifying gene sets with similar or opposite expression directions to those in GBM patients. TIMPs and associated genes are mainly associated with extracellular matrix organization and involve proteoglycan pathways in cancer. The expression levels of TIMPs in GBM correlate with the infiltration of various immune cells, including CD4+ T cells, macrophages, neutrophils, B cells, CD8+ T cells, and dendritic cells.
Our study inspires new ideas for the role of TIMPs in GBM and provides new directions for multiple treatment modalities, including immunotherapy, in GBM.
Glioblastoma (GBM) is the most aggressive of the four grades of gliomas classified by WHO  and accounts for 17% of all primary brain tumours . When radiotherapy was combined with temozolomide, the survival time of GBM patients reached 14.6 months, and the 2-year OS and 5-year OS rates were 27.2 and 9.8%, respectively . Matrix metalloproteinases (MMPs) are a family of zinc-dependent endopeptidases implicated in degrading components of the extracellular matrix (ECM) , tumour cell invasion and angiogenesis, and suppression of antitumour immune surveillance . Previous studies have found that MMPs are associated with different tumour types, such as papillary thyroid carcinomas  and lung adenocarcinoma . Several studies have indicated that MMPs are related to increased invasiveness  and a poor prognosis  in GBM. Tissue inhibitors of metalloproteinases (TIMPs) comprise four natural endogenous secreted proteins that inhibit the activity of MMPs . Therefore, TIMPs regulate cellular processes such as migration, proliferation, and survival by controlling ECM degradation via interaction with MMPs . In addition to this originally described feature , over the past decade, many additional relatively independent biological activities have been identified, including cell growth , apoptosis , differentiation , angiogenesis , and proliferation .
The TIMP gene family has been studied in various tumours, including GBM, but no systematic analysis has investigated the relationship between their expression and prognosis in glioblastoma. With the development of microarray technology, immune infiltration research, and the progress of single-cell technology, using bioinformatics analysis is vital to examine the expression pattern and prognostic value of the TIMP gene family in glioblastoma and to use cutting-edge technology and related research data for further analysis. This process will also lead to an update of our understanding of this ancient family.
Materials and methods
Oncomine (https://www.oncomine.org/resource/login.html) is a platform with powerful analysis functions that compute gene expression signatures, clusters, and gene-set modules from the data . We analysed the mRNA expression of TIMPs and compared the difference between the normal and tumour groups by Student’s t-test in various tumours on this platform. A p value < 0.05 and fold change ≥2 for the p value indicated that the difference between the normal and tumour groups was statistically significant.
GEPIA2 (http://gepia.cancer-pku.cn/index.html) is an updated version of Gene Expression Profiling Interactive Analysis (GEPIA) that processes mRNA data from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) projects . The website has various customized functions, such as similar genetic testing, correlation analysis, and dimensionality reduction analysis. In our study, we performed differential expression analysis of family genes as a validation set.
OSgbm (http://bioinfo.henu.edu.cn/GBM/GBMList.jsp.) is a web server that performs online survival analysis of glioblastoma . This website comprises 684 samples with transcriptome profiles and clinical information from TCGA, Gene Expression Omnibus (GEO), and Chinese Glioma Genome Atlas (CGGA). We used it to explore the impact of TIMPs on the prognosis of GBM patients. The survival analysis results are presented as Kaplan–Meier (KM) plots with hazard ratios (HRs) and log-rank plots with p values.
LinkedOmics (http://www.linkedomics.orglogin.php) is a public portal for multiomics data that mainly includes TCGA cancer types . The mutual combination of different modules allows the analysis of multiomics data from various cancers. Here, we explored the expression relationships between TIMPs and different genes in GBM.
Gliovis (http://gliovis.bioinfo.cnio.es/) is a user-friendly website tool that includes brain tumour research results containing more than 6,500 tumour samples, approximately 50 expression datasets, and mainly gliomas . We explored GBM data from the TCGA database at this server to investigate the correlation between the expression of TIMP family gene members and GBM subgroups according to different classifications.
STRING (https://string-db.org/) is a website to perform protein interactions. Its goal is to build a comprehensive, objective global network that can perform direct (physical) and indirect (functional) interactions of proteins. We used it to generate protein–protein interaction (PPI) networks of TIMPs.
GeneMANIA (https://genemania.org/) is a server that performs gene association and similar gene exploration . It can extend the submitted gene list by analysing a large amount of association data in the background and the gene association function in the list. We found the associated genes of TIMPs here and explored their different functions with GeneMANIA.
Enrichr (https://maayanlab.cloud/Enrichr/) is a search engine that performs multiple functional enrichment analyses of input genes by querying many collections of annotated genes [24,25,26]. We performed functional analysis here for the set of genes generated in GeneMANIA.
cBioPortal (http://cbioportal.org) is an open portal to explore and analyse multidimensional data from different cancers and combine molecular data and clinical data with visualization. We used this database to explore the mutation of TIMP family genes in GBM and link them to prognosis.
TIMER (http://timer.cistrome.org/) is a user-friendly portal that evaluates the molecular characterization of tumour interactions with immunity . It calculates the level of immunity in various infiltrating tumours in advance and provides six main analysis modules so that the relationship between tumours and immunity can be interactively explored. We conducted a correlative analysis of TIMPs at this site for measures of immune infiltration, including tumour purity and six immune cells in GBM.
Transcriptional levels of TIMPs in patients with GBM
Four TIMP factors have been identified in humans. We compared the transcriptional levels of TIMPs in GBM with those in normal samples using the filter of differential analysis in Glioblastoma vs. Normal and Brain Glioblastoma vs. Normal in ONCOMINE databases (Fig. 1). The mRNA expression levels of TIMP1 were significantly upregulated in patients with GBM in five datasets (Table 1). In detail, the expression level of TIMP1 in the tumour group was dramatically augmented with a fold change of 18.922 in Bredel Brain 2’s dataset . Additionally, the lowest fold change of 2.419 was found in the differential comparison between GBM and normal samples of neural stem cells in Lee Brain’s dataset . Regarding TIMP2, increased expression occurred in both TCGA Brain and Liang Brain  datasets (Table 1). The samples of these two datasets were brain vs. brain GBM, brain and cerebellum vs. GBM. The fold change of TIMP2 in TCGA Brain was 2.218, and the latter was 2.304. In Bredel Brain 2’s  and Lee Brain’s datasets , TIMP3 was overexpressed with a fold change≥2 (Table 1). TIMP4 was overexpressed in two types of samples in the TCGA Brain’s dataset: brain GBM vs. brain and GBM vs. brain. Their fold changes, 8.998 and 8.887, respectively, were both higher (Table 1). The lowest fold change of 2.774 appeared in Shai Brain’s dataset , and the sample of this dataset was GBM vs. White Matter.
Differential expression of TIMPs in GBM
In the GEPIA (Gene Expression Profiling Interactive Analysis) dataset (http://gepia.cancer-pku.cn/), we compared the mRNA expression level of TIMPs between GBM and normal samples. The gene expression levels of TIMP1, TIMP2, TIMP3, and TIMP4 were all higher in GBM than in normal samples (Fig. 2A–E).
The prognostic values of TIMPs in glioblastoma and co-expression genes of significant gene
We investigated the prognostic value of TIMP1, TIMP2, TIMP3 and TIMP4 in GBM data from TCGA using OSgbm . The TIMP3 overexpression group was associated with better overall survival (OS) and disease-specific survival (DSS) (Fig. 3A). TIMP4 overexpression was associated with better OS but without significance in DSS (Fig. 3A). Neither TIMP1 nor TIMP2 showed differences in OS and DSS (Fig. 3A). Next, we found a positive correlation between TIMP3 and TIMP4 expression in the LinkedOmics dataset (R = 0.2168, p <0.01) (Fig. 3B) .
Given the prognostic value of TIMP3 and TIMP4, using the functional module of LinkedOmics, we explored their associated genes in the TCGA_GBM cohort. The genes positively (dark red dots) and negatively (dark green dots) correlated with TIMP3 and TIMP4 are shown in Fig. 3C. The top 50 significant genes positively and negatively correlated with TIMP3 and TIMP4 are shown in the heat map (Fig. 3D). TIMP3 expression showed a positive association with the expression of TOB2 (positive rank #1; r = 0.513; p = 1.275e-11), MTMR3 (r = 0.511; p = 1.448e-11), and EIF4ENIF1 (r = 0.510; p = 1.832e-11). The top three genes with a positive correlation with TIMP4 were SCRG1 (r = 0.544; p = 3.778e-13), SHISA4 (r = 0.540; p = 5.826e-13) and COMMD2 (r = 0.526; p = 2.890e-12). These genes are more likely to be protective genes similar to TIMP3 and TIMP4.
Validation of the expression of TIMPs in different subgroups of GBM
We divided the GBM data from TCGA into different subgroups according to different criteria in Gliovis. Next, we investigated whether TIMPs were differentially expressed between different subgroups (Fig. 4). The four TIMP family members were not differentially expressed between sexes.
The expression levels of both TIMP1 and TIMP2 were significantly decreased in patients with IDH1 mutations. However, no association was found between TIMP3 and TIMP4 expression and IDH1 mutation. Among the TIMP family, only TIMP1 expression correlated with methylation of the MGMT promoter. We focused on evaluating the differential expression of TIMPs between different subtypes of GBM. The four genes showed significant expression differences among classical, mesenchymal and proneural GBM. TIMP1 and TIMP2 expression was highest in the mesenchymal region, but TIMP3 and TIMP4 expression lower in this region. Specific data can be viewed in the Supplementary file.
Protein-protein interaction analyses and functional enrichment analysis of TIMPs
We performed protein–protein interaction PPI network analysis of TIMPs using STRING to explore the potential interactions among them. Eventually, we obtained a protein interaction network comprising 9 nodes and 29 edges (Fig. 5A). We also used the GeneMANIA database to identify the genes associated with TIMPs, revealing that the functions of TIMPs with 20 related molecules (such as RECK, MMP1, MMP14, MMP3, MMP2, AGTR2, and PCSK5) were primarily associated with extracellular matrix organization, extracellular structure organization, negative regulation of peptidase activity and endopeptidase activity (Fig. 5B). To further analyze the functions of TIMPs and related genes as well as to make a comparative control with the functions presented on GeneMANIA, we performed functional enrichment analysis of these genes at the Enrichr website. The results showed that the top 10 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways for TIMPs and their 20 correlated genes were mainly associated with proteoglycans in cancer, the IL-17 signalling pathway, the TNF signalling pathway, the GnRH signalling pathway, and others (Fig. 5C) . These pathways may involve mechanisms of progression in GBM. Next, using the Enrichr tool, GO analysis was performed using TIMPs and their correlated genes to analyse functions in biological processes, molecular functions, and cellular components. TIMPs and correlated genes were mainly associated with extracellular matrix disassembly in biological processes (Fig. 5D), tertiary granule lumen in cellular components (Fig. 5E), and metalloendopeptidase activity in molecular functions (Fig. 5F).
Genetic alteration of TIMPs in patients with GBM
We analysed the genetic changes of TIMPs in the database. The TIMP gene was genetically altered in 23% of cases. The alteration frequency was 14% for TIMP1, 1.5% for TIMP2, 8% for TIMP3 and 5% for TIMP4 (Fig. 6A). However, these data primarily represent expression changes. We explored the mutation status of TIMPs in the four databases Glioblastoma (TCGA, Nature 2008), Glioblastoma Multiforme (TCGA, Firehose Legacy), Glioblastoma Multiforme (TCGA, PanCancer Atlas) and Glioblastoma (TCGA, Cell 2013) and found that the highest mutation frequency did not exceed 1.4% (Fig. 6B). We also performed correlation analysis between patients with genetic alterations in TIMPs and prognostic conditions in the database and found that the patients were significantly associated with better OS and DFS (Fig. 6C, D).
Relationship between TIMPs expression and immune infiltration in GBM
Immune infiltration is a component of the tumour microenvironment and is highly correlated with the diagnosis, progression, and prognosis of tumours. We used the TIMER database in this study to reveal the correlation between the gene expression levels of TIMP family members and immune cell infiltration in patients with GBM (Fig. 7A). TIMP1 expression was negatively associated with tumour purity but positively correlated with dendritic cells. Additionally, this correlation was also observed for TIMP2. TIMP3 expression was associated with B cells, CD4+ T cells, macrophages, and neutrophils. However, no negative results are shown for TIMP3. Finally, TIMP4 expression was positively associated with the infiltration of B cells, CD8+ T cells, and macrophages in GBM patients. We also explored the relationship between the copy number variation (CNV) of TIMPs and immune infiltration level (Fig. 7B). The association of CNV with immune infiltration varied among the four genes. The levels of CNS in TIMP1 and TIMP2 were associated with more immune infiltration indicators, and , their arm-level deletion and arm-level gain were closely related to neutrophils. Regarding the remaining two genes, TIMP3 and TIMP4, that affected the prognosis, the CNV level of TIMP3 was the least significantly associated with immune infiltration indicators but was the only member of the family whose CNV affected the immune infiltration level of B cells. The CNV of TIMP4 was significantly correlated with both CD4+ T cells and dendritic cells, a finding similar to that of TIMP2. However, the CNV showed no correlation with neutrophils.
Previous studies have found that TIMP gene family members are differentially expressed in multiple cancers [35,36,37]. Because MMPs are highly correlated with the prognosis of multiple cancers, TIMPs, at least theoretically, may affect prognosis in some cancer cases as specific inhibitors. In brain tumours, including glioma, many reports have investigated the differential expression of TIMPs [38,39,40]. However, previous reports on TIMPs in GBM are largely fragmented, and systematic analysis of the TIMP family in GBM is rare. In the present study, we performed systemic analysis concerning TIMPs for the prognosis of patients with GBM.
As shown in the results (Fig. 1), TIMPs were differentially expressed in brain and CNS cancer, colorectal cancer, and breast cancer, a finding that is also consistent with previous reports [35, 41]. Interestingly, the expression levels of these genes were almost all upregulated in brain and CNS cancer, but TIMP expression was downregulated in normal individuals (Fig. 2). TIMP-1 demonstrated both antitumour and antimetastatic functions in multiple genetic models [42, 43]. However, in our study, TIMP1 expression was not associated with the prognosis of GBM (Fig. 3A). Most of the studies on TIMP1 gene expression levels showed a negative relationship with prognosis [44, 45]. High TIMP-1 expression is strongly associated with a poor prognosis in almost all known cancer types . Many explanations exist for this contradictory phenomenon, such as the interaction between TIMP1 and CD63 to activate MAPK signalling and participation in the HGF/c-Met pathway to promote cancer progression. Additionally, MMPs have demonstrated potential cancer-inhibiting effects [47, 48]. Our results were not significant, a finding similar to the findings of Paulina Vaitkiene . However, the high expression group of TIMP1 showed a poor prognosis, similar to the findings of Charlotte Aaberg-Jessen . This phenomenon may be due to the insufficient number of patients and requires further work.
TIMP-2 is a unique member of the TIMP family because it directly inhibits endothelial cell proliferation instead of metalloproteinase activity . TIMP2 has recently made adequate progress among biomarkers of renal injury [51, 52]. Regarding cancer, the current findings show that the role of TIMP2 in prognosis depends on the tissue of origin of the tumor . However, epigenetic silencing of TIMP2 can promote invasion of ovarian cancer (OC) . Concerning renal cell carcinoma (RCC), TIMP2 is inversely correlated with markers of tumour progression . However, TIMP2 is a poor prognostic marker in gastric cancer , in contrast to its role in colorectal cancer . TIMP2 is involved in glioma inhibition by different drugs as an intermediate link in the upregulated expression in different studies [58, 59]. However, regarding GBM, laboratory evidence has shown that TIMP2 promotes MMP2 activation and GBM cell invasion . This finding indicates that the effect of TIMP2 on GBM is more complex than previously thought in the real environment. We showed that TIMP2 expression was positive, although not significantly, associated with the prognosis of GBM patients (Fig. 3A). This result may also be a manifestation of its complex role.
TIMP3 is the only TIMP family gene with an affinity for proteoglycans in the ECM, and it has the broadest range of substrates . TIMP3 frequently leads to a poor prognosis because of silencing in tumours . In addition to being a marker for predicting cancer progression, TIMP3 can also be used as a therapeutic target for cancer and a marker of tumour sensitivity to drugs [62,63,64]. In conclusion, TIMP3 has been shown to play a role in inhibiting cancer progression in numerous studies. A similar role has been reported for TIMP3 in glioma. MicroRNA 21 (miR-21) enhances the malignant degree of glioma by inhibiting TIMP3 expression . In another study, upregulated expression of TIMP3 emerged in euxanthone-inhibited GBM cell lines . TIMP3 can suppress angiogenesis by binding to vascular endothelial growth factor receptor 2 (VEGFR2) and disrupting the interaction of VEGF with VEGFR2 and downstream signalling . Additionally, TIMP3 induces apoptosis and inhibits cell proliferation . These functions all underlie TIMP3 as an inhibitor in various cancers, including GBM. Our findings also reflect these characteristics of TIMP3. TIMP3 overexpression was associated with better prognosis in GBM patients for both OS and DSS (Fig. 3A). This finding provides evidence for the positive impact of TIMP3 on the prognosis of GBM patients and validates the potential of TIMP3 as an intermediate link in drug therapy for various malignancies, including GBM, for further development. We also list genes associated with differential expression of TIMP3 in GBM that may develop similar application values to TIMP3 because of its similar role in GBM (Fig. 3D). Existing studies also demonstrated similar potential for these genes. For example, TOP2 is a member of the TOB/BTG protein family and inhibits cell proliferation and stimulates cell differentiation . The dysregulation of TOP2 is associated with the development of hepatocellular carcinoma (HCC) . MTMR3 belongs to the myotubularin-related protein family and is a target gene of various miRNAs that regulate tumour development [71, 72]. MTMR3 modulates the proliferative phenotype and apoptosis of glioma cells by mediating the Wnt/β-catenin pathway and correlates with glioma grade . EIF4ENIF1 encodes a nucleocytoplasmic shuttle protein for the translation initiation factor eIF4E that constitutes a repressor complex involved in neurogenesis and confers resistance to DNA damaging agents through interaction with eIF4 [74,75,76].
TIMP4 is expressed only in specific tissues of the human body, and this restricted expression implies its functional specificity compared with other family members . Thus, relatively few studies have investigated the role of TIMP4 in tumours, making its role in some tumours unclear in published studies. In studies on renal cell carcinomas (RCCs), TIMP4 showed opposite expression patterns in clear cell and papillary cell carcinomas . regarding breast cancer, in contrast to Wang’s study in which TIMP4 played a role in stimulating tumorigenesis , Jiang’s study showed that TIMP4 inhibited tumour growth and invasion . Recent studies have revealed that TIMP4 is involved in the inhibitory effect of various substances on cancer, showing a negative regulatory role in the progression of various tumours. Low TIMP4 expression induces resistance to chemotherapeutic agents in non-small cell lung cancer (NSCLC) . Stable TIMP-4 knockdown promotes the proliferation and migration of well-differentiated liposarcoma (WDLS) . In bladder and prostate tumours, TIMP4 is downregulated by chemokine (C-X-C motif) ligand 1 (CXCL1) as an intermediate link to maintain tumour growth . Regarding the relationship between TIMP4 and glioma progression, the directivity of the available findings is inconsistent. TIMP4 binds CRN2 (an actin filament binding protein) and promotes perivascular invasion of GBM cells . High TIMP4/CD63 coexpression contributes to the poor prognosis of GBM patients . TIMP4 has been shown in Pullen’s study to have an inverse expression relationship with MMP1, which performs enhanced tumorigenic function in GBM . TIMP4 is also downregulated during human herpesvirus 6 (HHV-6) promotion of glioma development and progression . In our study, high TIMP4 expression was associated with better OS in GBM patients. Although the relationship between the expression level of TIMP4 and DSS was not significant, the effect of TIMP4 expression on DSS was positive (Fig. 3A). This result may be explained by the insufficient number of enrolled patients. Notably, a significantly positive relationship was found between TIMP4 and TIMP3 expression (Fig. 3B), warranting further investigation and providing evidence for the effect of TIMP4 on the prognosis of GBM patients. Finally, a list of genes highly correlated with TIMP4 expression is also provided for exploration (Fig. 3D). Recent studies provide a basis for the further exploration of these genes. For example, SCRG1 (Scrapie Responsive Gene 1) is involved in neurodegeneration and autophagy and is associated with transmissible spongiform encephalopathy. It can also maintain the stemness of stem cells through the SCRG1/BST1 axis [88,89,90]. SHISA4 belongs to the SHISA family and mediates synaptic transmission in the central nervous system by constituting an auxiliary subunit of the AMPA receptor . COMMD2 is the second member of the copper metabolism MURR1 domain protein family, in which COMMD1 inhibits GBM cell proliferation and promotes lung cancer cell apoptosis [92, 93].
Based on the clinical and molecular pathological features of GBM, multiple classification methods exist. Our study demonstrated that the expression of TIMP family genes did not show a correlation with sex. The gene encoding isocitrate dehydrogenase (IDH) distinguishes GBM patients from IDH-mutant (IDH-mt) and wild-type (IDH-wt) GBM patients . Existing studies have shown that the IDH mutant body type has a better prognosis and therapeutic effect than the IDH wild-type, accounting for approximately 10% of GBM [95, 96]. TIMP1 expression was markedly elevated in wild-type GBM, corresponding to patients with high TIMP1 expression showing a worse prognosis, although not significantly. TIMP2 expression did not differ greatly between the two groups, corresponding to the close results of the two prognostic subgroups. The MGMT gene is associated with DNA repair, and its methylation produces gene silencing, allowing glioblastoma patients receiving chemotherapy to achieve a better prognosis . We observed even lower TIMP1 expression in MGMT-methylated patients, a finding that may also be responsible for the better prognosis of patients with downregulated TIMP1 expression. GBM was originally divided into four subtypes according to gene expression profiles—classical, mesenchymal, proneural, and neural . However, some studies have found that the neural subtype is produced because normal samples are contaminated . Mutations of EGFR often occur in the classical subtype, and it and the mesenchymal subtype are more likely to be affected by aggressive treatment . GBM patients with proneural subtypes frequently have rare genetic mutations, particularly IDH mutations with better prognosis . Our study lays the foundation to evaluate the mechanism and how differential expression of TIMPs affects the survival of GBM patients with central subtypes.
The tumour microenvironment (TME) has gained increasing attention because it regulates tumour progression and significantly affects the treatment response [100, 101]. The infiltration of immune cells in the TME promotes or antagonizes tumours . We analysed the effect of TIMP family members on immune cell infiltration in GBM with the TIMER database and found that the expression and CNV of TIMPs are closely related to the infiltration of six immune cells. These four members can be identified in GBM because of differences in multiple microenvironment indicators, including tumour purity. Our results help to further reveal the molecular mechanism of GBM progression and are expected to use this gene family as a link in immunotherapy to specifically affect the GBM tumour microenvironment and improve the treatment of GBM.
Our study has some limitations, primarily because our data were all obtained from public databases, making some heterogeneity in the source and use of the data inevitable. Our results mainly provide a direction and basis for studies targeting TIMPs in GBM. Further experimental and clinical validation for translational purposes is warranted.
We conducted a systematic study on the expression, prognosis, and immune infiltration of TIMPs in GBM. Our results confirm that the expression of all four members of the TIMP family is upregulated in GBM. High expression of both TIMP3 and TIMP4 can be used as a predictor of better survival in GBM patients and a positive relationship was found between them. These findings led to the identification of related genes for further analysis. Finally, the expression of this gene family was associated with the infiltration of immune cells, including CD4+ T cells, macrophages, neutrophils, B cells, CD8+ T cells, and dendritic cells, to varying degrees. This effect provides new exploration directions for research on the TIMP family and GBM.
Availability of data and materials
Publicly available datasets were analyzed in this study. Institutional review board approval was not demanded in our study for all database is publicly available; This data can be found in following website. Oncomine: http://www.oncomine.org; GEPIA2: http://gepia2.cancer-pku.cn/; OSgbm: http://bioinfo.henu.edu.cn/GBM/GBMList.jsp; LinkedOmics: http://www.linkedomics.orglogin.php; STRING: https://string-db.org; GeneMANIA: http://www.genemania.org; Enrichr: https://maayanlab.cloud/Enrichr/; KEGG: https://www.kegg.jp/kegg/kegg1.html; cBioPortal: http://cbioportal.org; TIMER: https://cistrome.org/.
The tissue inhibitors of metalloproteinases
Disease special survival
Neural stem cells
Basic fibroblast growth factor
Gene Expression Profiling Interactive Analysis
The Cancer Genome Atlas
Gene Expression Omnibus
Chinese Glioma Genome Atlas
Kyoto Encyclopedia of Genes and Genomes
Copy number variation
Renal cell carcinoma
Vascular endothelial growth factor receptor 2
Renal cell carcinomas
Non-small cell lung cancer
Human herpesvirus 6
Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, Jouvet A, et al. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol. 2007;114(2):97–109.
Ostrom QT, Bauchet L, Davis FG, Deltour I, Fisher JL, Langer CE, et al. The epidemiology of glioma in adults: a “state of the science” review. Neuro Oncol. 2014;16(7):896–913.
Stupp R, Hegi ME, Mason WP, van den Bent MJ, Taphoorn MJ, Janzer RC, et al. Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet Oncol. 2009;10(5):459–66.
Lu P, Takai K, Weaver VM, Werb Z. Extracellular matrix degradation and remodeling in development and disease. Cold Spring Harb Perspect Biol. 2011;3(12).
Ahir BK, Engelhard HH, Lakka SS. Tumor development and angiogenesis in adult brain tumor: glioblastoma. Mol Neurobiol. 2020;57(5):2461–78.
Wang JR, Li XH, Gao XJ, An SC, Liu H, Liang J, et al. Expression of MMP-13 is associated with invasion and metastasis of papillary thyroid carcinoma. Eur Rev Med Pharmacol Sci. 2013;17(4):427–35.
Cai X, Zhu H, Li Y. PKCζ, MMP-2 and MMP-9 expression in lung adenocarcinoma and association with a metastatic phenotype. Mol Med Rep. 2017;16(6):8301–6.
Kenig S, Alonso MB, Mueller MM, Lah TT. Glioblastoma and endothelial cells cross-talk, mediated by SDF-1, enhances tumour invasion and endothelial proliferation by increasing expression of cathepsins B, S, and MMP-9. Cancer Lett. 2010;289(1):53–61.
Turpeenniemi-Hujanen T. Gelatinases (MMP-2 and -9) and their natural inhibitors as prognostic indicators in solid cancers. Biochimie. 2005;87(3-4):287–97.
Arpino V, Brock M, Gill SE. The role of TIMPs in regulation of extracellular matrix proteolysis. Matrix Biol. 2015;44-46:247–54.
Ries C. Cytokine functions of TIMP-1. Cell Mol Life Sci. 2014;71(4):659–72.
Lambert E, Dassé E, Haye B, Petitfrère E. TIMPs as multifacial proteins. Crit Rev Oncol Hematol. 2004;49(3):187–98.
Yamashita K, Suzuki M, Iwata H, Koike T, Hamaguchi M, Shinagawa A, et al. Tyrosine phosphorylation is crucial for growth signaling by tissue inhibitors of metalloproteinases (TIMP-1 and TIMP-2). FEBS Lett. 1996;396(1):103–7.
Singla DK, McDonald DE. Factors released from embryonic stem cells inhibit apoptosis of H9c2 cells. Am J Physiol Heart Circ Physiol. 2007;293(3):H1590–5.
Smith MF, Ricke WA, Bakke LJ, Dow MP, Smith GW. Ovarian tissue remodeling: role of matrix metalloproteinases and their inhibitors. Mol Cell Endocrinol. 2002;191(1):45–56.
Cruz-Munoz W, Khokha R. The role of tissue inhibitors of metalloproteinases in tumorigenesis and metastasis. Crit Rev Clin Lab Sci. 2008;45(3):291–338.
Fassina G, Ferrari N, Brigati C, Benelli R, Santi L, Noonan DM, et al. Tissue inhibitors of metalloproteases: regulation and biological activities. Clin Exp Metastasis. 2000;18(2):111–20.
Rhodes DR, Yu J, Shanker K, Deshpande N, Varambally R, Ghosh D, et al. ONCOMINE: a cancer microarray database and integrated data-mining platform. Neoplasia. 2004;6(1):1–6.
Tang Z, Kang B, Li C, Chen T, Zhang Z. GEPIA2: an enhanced web server for large-scale expression profiling and interactive analysis. Nucleic Acids Res. 2019;47(W1):W556–w60.
Dong H, Wang Q, Li N, Lv J, Ge L, Yang M, et al. OSgbm: an online consensus survival analysis web server for glioblastoma. Front Genet. 2019;10:1378.
Vasaikar SV, Straub P, Wang J, Zhang B. LinkedOmics: analyzing multi-omics data within and across 32 cancer types. Nucleic Acids Res. 2018;46(D1):D956–d63.
Bowman RL, Wang Q, Carro A, Verhaak RG, Squatrito M. GlioVis data portal for visualization and analysis of brain tumor expression datasets. Neuro Oncol. 2017;19(1):139–41.
Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, et al. The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010;38(Web Server issue):W214–20.
Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14:128.
Kuleshov MV, Jones MR, Rouillard AD, Fernandez NF, Duan Q, Wang Z, et al. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016;44(W1):W90–7.
Xie Z, Bailey A, Kuleshov MV, Clarke DJB, Evangelista JE, Jenkins SL, et al. Gene set knowledge discovery with enrichr. Curr Protoc. 2021;1(3):e90.
Li T, Fan J, Wang B, Traugh N, Chen Q, Liu JS, et al. TIMER: a web server for comprehensive analysis of tumor-infiltrating immune cells. Cancer Res. 2017;77(21):e108–e10.
Bredel M, Bredel C, Juric D, Harsh GR, Vogel H, Recht LD, et al. Functional network analysis reveals extended gliomagenesis pathway maps and three novel MYC-interacting genes in human gliomas. Cancer Res. 2005;65(19):8679–89.
Lee J, Kotliarova S, Kotliarov Y, Li A, Su Q, Donin NM, et al. Tumor stem cells derived from glioblastomas cultured in bFGF and EGF more closely mirror the phenotype and genotype of primary tumors than do serum-cultured cell lines. Cancer Cell. 2006;9(5):391–403.
Liang Y, Diehn M, Watson N, Bollen AW, Aldape KD, Nicholas MK, et al. Gene expression profiling reveals molecularly and clinically distinct subtypes of glioblastoma multiforme. Proc Natl Acad Sci U S A. 2005;102(16):5814–9.
Shai R, Shi T, Kremen TJ, Horvath S, Liau LM, Cloughesy TF, et al. Gene expression profiling identifies molecular subtypes of gliomas. Oncogene. 2003;22(31):4918–23.
Murat A, Migliavacca E, Gorlia T, Lambiv WL, Shay T, Hamou MF, et al. Stem cell-related “self-renewal” signature and high epidermal growth factor receptor expression associated with resistance to concomitant chemoradiotherapy in glioblastoma. J Clin Oncol. 2008;26(18):3015–24.
Sun L, Hui AM, Su Q, Vortmeyer A, Kotliarov Y, Pastorino S, et al. Neuronal and glioma-derived stem cell factor induces angiogenesis within the brain. Cancer Cell. 2006;9(4):287–300.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Böckelman C, Beilmann-Lehtonen I, Kaprio T, Koskensalo S, Tervahartiala T, Mustonen H, et al. Serum MMP-8 and TIMP-1 predict prognosis in colorectal cancer. BMC Cancer. 2018;18(1):679.
Ozden F, Saygin C, Uzunaslan D, Onal B, Durak H, Aki H. Expression of MMP-1, MMP-9 and TIMP-2 in prostate carcinoma and their influence on prognosis and survival. J Cancer Res Clin Oncol. 2013;139(8):1373–82.
Azevedo Martins JM, Rabelo-Santos SH, do Amaral Westin MC, Zeferino LC. Tumoral and stromal expression of MMP-2, MMP-9, MMP-14, TIMP-1, TIMP-2, and VEGF-A in cervical cancer patient survival: a competing risk analysis. BMC Cancer. 2020;20(1):660.
Guo H, Sun Z, Wei J, Xiang Y, Qiu L, Guo L, et al. Expressions of matrix metalloproteinases-9 and tissue inhibitor of metalloproteinase-1 in pituitary adenomas and their relationships with prognosis. Cancer Biother Radiopharm. 2019;34(1):1–6.
Mashayekhi F, Saberi A, Mashayekhi S. Serum TIMP1 and TIMP2 concentration in patients with different grades of meningioma. Clin Neurol Neurosurg. 2018;170:84–7.
Aaberg-Jessen C, Christensen K, Offenberg H, Bartels A, Dreehsen T, Hansen S, et al. Low expression of tissue inhibitor of metalloproteinases-1 (TIMP-1) in glioblastoma predicts longer patient survival. J Neurooncol. 2009;95(1):117–28.
Wu ZS, Wu Q, Yang JH, Wang HQ, Ding XD, Yang F, et al. Prognostic significance of MMP-9 and TIMP-1 serum and tissue expression in breast cancer. Int J Cancer. 2008;122(9):2050–6.
Khokha R. Suppression of the tumorigenic and metastatic abilities of murine B16-F10 melanoma cells in vivo by the overexpression of the tissue inhibitor of the metalloproteinases-1. J Natl Cancer Inst. 1994;86(4):299–304.
Khokha R, Waterhouse P, Yagel S, Lala PK, Overall CM, Norton G, et al. Antisense RNA-induced reduction in murine TIMP levels confers oncogenicity on Swiss 3T3 cells. Science. 1989;243(4893):947–50.
Würtz S, Schrohl AS, Sørensen NM, Lademann U, Christensen IJ, Mouridsen H, et al. Tissue inhibitor of metalloproteinases-1 in breast cancer. Endocr Relat Cancer. 2005;12(2):215–27.
Grunnet M, Mau-Sørensen M, Brünner N. Tissue inhibitor of metalloproteinase 1 (TIMP-1) as a biomarker in gastric cancer: a review. Scand J Gastroenterol. 2013;48(8):899–905.
Jackson HW, Defamie V, Waterhouse P, Khokha R. TIMPs: versatile extracellular regulators in cancer. Nat Rev Cancer. 2017;17(1):38–53.
Grünwald B, Vandooren J, Gerg M, Ahomaa K, Hunger A, Berchtold S, et al. Systemic ablation of MMP-9 triggers invasive growth and metastasis of pancreatic cancer via deregulation of IL6 expression in the bone marrow. Mol Cancer Res. 2016;14(11):1147–58.
Jobin PG, Butler GS, Overall CM. New intracellular activities of matrix metalloproteinases shine in the moonlight. Biochim Biophys Acta Mol Cell Res. 2017;1864(11 Pt A):2043–55.
Vaitkiene P, Urbanaviciute R, Grigas P, Steponaitis G, Tamasauskas A, Skiriutė D. Identification of astrocytoma blood serum protein profile. Cells. 2019;9(1).
Seo DW, Li H, Guedez L, Wingfield PT, Diaz T, Salloum R, et al. TIMP-2 mediated inhibition of angiogenesis: an MMP-independent mechanism. Cell. 2003;114(2):171–80.
Andreucci M, Faga T, Pisani A, Perticone M, Michael A. The ischemic/nephrotoxic acute kidney injury and the use of renal biomarkers in clinical practice. Eur J Intern Med. 2017;39:1–8.
Waskowski J, Pfortmueller CA, Schenk N, Buehlmann R, Schmidli J, Erdoes G, et al. (TIMP2) x (IGFBP7) as early renal biomarker for the prediction of acute kidney injury in aortic surgery (TIGER). A single center observational study. PLoS One. 2021;16(1):e0244658.
Eckfeld C, Häußler D, Schoeps B, Hermann CD, Krüger A. Functional disparities within the TIMP family in cancer: hints from molecular divergence. Cancer Metastasis Rev. 2019;38(3):469–81.
Yi X, Guo J, Guo J, Sun S, Yang P, Wang J, et al. EZH2-mediated epigenetic silencing of TIMP2 promotes ovarian cancer migration and invasion. Sci Rep. 2017;7(1):3568.
Lyu J, Zhu Y, Zhang Q. An increased level of MiR-222-3p is associated with TMP2 suppression, ERK activation and is associated with metastasis and a poor prognosis in renal clear cell carcinoma. Cancer Biomark. 2020;28(2):141–9.
Wang W, Zhang Y, Liu M, Wang Y, Yang T, Li D, et al. TIMP2 is a poor prognostic factor and predicts metastatic biological behavior in gastric cancer. Sci Rep. 2018;8(1):9629.
Jin Z, Zhang G, Liu Y, He Y, Yang C, Du Y, et al. The suppressive role of HYAL1 and HYAL2 in the metastasis of colorectal cancer. J Gastroenterol Hepatol. 2019;34(10):1766–76.
Wang Z, Wang B, Shi Y, Xu C, Xiao HL, Ma LN, et al. Oncogenic miR-20a and miR-106a enhance the invasiveness of human glioma stem cells by directly targeting TIMP-2. Oncogene. 2015;34(11):1407–19.
Mei PJ, Chen YS, Du Y, Bai J, Zheng JN. PinX1 inhibits cell proliferation, migration and invasion in glioma cells. Med Oncol. 2015;32(3):73.
Lu KV, Jong KA, Rajasekaran AK, Cloughesy TF, Mischel PS. Upregulation of tissue inhibitor of metalloproteinases (TIMP)-2 promotes matrix metalloproteinase (MMP)-2 activation and cell invasion in a human glioblastoma cell line. Lab Invest. 2004;84(1):8–20.
Fan D, Kassiri Z. Biology of tissue inhibitor of metalloproteinase 3 (TIMP3), and its therapeutic implications in cardiovascular pathology. Front Physiol. 2020;11:661.
Span PN, Lindberg RL, Manders P, Tjan-Heijnen VC, Heuvel JJ, Beex LV, et al. Tissue inhibitors of metalloproteinase expression in human breast cancer: TIMP-3 is associated with adjuvant endocrine therapy success. J Pathol. 2004;202(4):395–402.
Deng X, Bhagat S, Dong Z, Mullins C, Chinni SR, Cher M. Tissue inhibitor of metalloproteinase-3 induces apoptosis in prostate cancer cells and confers increased sensitivity to paclitaxel. Eur J Cancer. 2006;42(18):3267–73.
Yu JL, Lv P, Han J, Zhu X, Hong LL, Zhu WY, et al. Methylated TIMP-3 DNA in body fluids is an independent prognostic factor for gastric cancer. Arch Pathol Lab Med. 2014;138(11):1466–73.
Gabriely G, Wurdinger T, Kesari S, Esau CC, Burchard J, Linsley PS, et al. MicroRNA 21 promotes glioma invasion by targeting matrix metalloproteinase regulators. Mol Cell Biol. 2008;28(17):5369–80.
Li W, Du H, Zhou G, Song D. Euxanthone represses the proliferation, migration, and invasion of glioblastoma cells by modulating STAT3/SHP-1 signaling. Anat Rec. 2020.
Qi JH, Ebrahem Q, Moore N, Murphy G, Claesson-Welsh L, Bond M, et al. A novel function for tissue inhibitor of metalloproteinases-3 (TIMP3): inhibition of angiogenesis by blockage of VEGF binding to VEGF receptor-2. Nat Med. 2003;9(4):407–15.
Baker AH, George SJ, Zaltsman AB, Murphy G, Newby AC. Inhibition of invasion and induction of apoptotic cell death of cancer cell lines by overexpression of TIMP-3. Br J Cancer. 1999;79(9-10):1347–55.
Winkler GS. The mammalian anti-proliferative BTG/Tob protein family. J Cell Physiol. 2010;222(1):66–72.
Shen H, Li W, Tian Y, Xu P, Wang H, Zhang J, et al. Upregulation of miR-362-3p modulates proliferation and anchorage-independent growth by directly targeting Tob2 in hepatocellular carcinoma. J Cell Biochem. 2015;116(8):1563–73.
Lin Y, Zhao J, Wang H, Cao J, Nie Y. miR-181a modulates proliferation, migration and autophagy in AGS gastric cancer cells and downregulates MTMR3. Mol Med Rep. 2017;15(5):2451–6.
Gong Y, He T, Yang L, Yang G, Chen Y, Zhang X. The role of miR-100 in regulating apoptosis of breast cancer cells. Sci Rep. 2015;5:11650.
Yan Y, Yan H, Wang Q, Zhang L, Liu Y, Yu H. MicroRNA 10a induces glioma tumorigenesis by targeting myotubularin-related protein 3 and regulating the Wnt/β-catenin signaling pathway. FEBS J. 2019;286(13):2577–92.
Martínez A, Sesé M, Losa JH, Robichaud N, Sonenberg N, Aasen T, et al. Phosphorylation of eIF4E confers resistance to cellular stress and DNA-damaging agents through an interaction with 4E-T: a rationale for novel therapeutic approaches. PLoS One. 2015;10(4):e0123352.
Dostie J, Ferraiuolo M, Pause A, Adam SA, Sonenberg N. A novel shuttling protein, 4E-T, mediates the nuclear import of the mRNA 5’ cap-binding protein, eIF4E. EMBO J. 2000;19(12):3142–56.
Räsch F, Weber R, Izaurralde E, Igreja C. 4E-T-bound mRNAs are stored in a silenced and deadenylated form. Genes Dev. 2020;34(11-12):847–60.
Melendez-Zajgla J, Del Pozo L, Ceballos G, Maldonado V. Tissue inhibitor of metalloproteinases-4. The road less traveled. Mol Cancer. 2008;7:85.
Hagemann T, Gunawan B, Schulz M, Füzesi L, Binder C. mRNA expression of matrix metalloproteases and their inhibitors differs in subtypes of renal cell carcinomas. Eur J Cancer. 2001;37(15):1839–46.
Wang M, Liu YE, Greene J, Sheng S, Fuchs A, Rosen EM, et al. Inhibition of tumor growth and metastasis of human breast cancer cells transfected with tissue inhibitor of metalloproteinase 4. Oncogene. 1997;14(23):2767–74.
Jiang Y, Wang M, Celiker MY, Liu YE, Sang QX, Goldberg ID, et al. Stimulation of mammary tumorigenesis by systemic tissue inhibitor of matrix metalloproteinase 4 gene delivery. Cancer Res. 2001;61(6):2365–70.
Zhang YW, Zheng Y, Wang JZ, Lu XX, Wang Z, Chen LB, et al. Integrated analysis of DNA methylation and mRNA expression profiling reveals candidate genes associated with cisplatin resistance in non-small cell lung cancer. Epigenetics. 2014;9(6):896–909.
Shrestha M, Ando T, Chea C, Sakamoto S, Nishisaka T, Ogawa I, et al. The transition of tissue inhibitor of metalloproteinases from -4 to -1 induces aggressive behavior and poor patient survival in dedifferentiated liposarcoma via YAP/TAZ activation. Carcinogenesis. 2019;40(10):1288–97.
Miyake M, Furuya H, Onishi S, Hokutan K, Anai S, Chan O, et al. Monoclonal antibody against CXCL1 (HL2401) as a novel agent in suppressing IL6 expression and tumoral growth. Theranostics. 2019;9(3):853–67.
Solga R, Behrens J, Ziemann A, Riou A, Berwanger C, Becker L, et al. CRN2 binds to TIMP4 and MMP14 and promotes perivascular invasion of glioblastoma cells. Eur J Cell Biol. 2019;98(5-8):151046.
Rorive S, Lopez XM, Maris C, Trepant AL, Sauvage S, Sadeghi N, et al. TIMP-4 and CD63: new prognostic biomarkers in human astrocytomas. Mod Pathol. 2010;23(10):1418–28.
Pullen NA, Anand M, Cooper PS, Fillmore HL. Matrix metalloproteinase-1 expression enhances tumorigenicity as well as tumor-related angiogenesis and is inversely associated with TIMP-4 expression in a model of glioblastoma. J Neurooncol. 2012;106(3):461–71.
Gu B, Li M, Zhang Y, Li L, Yao K, Wang S. DR7 encoded by human herpesvirus 6 promotes glioma development and progression. Cancer Manag Res. 2019;11:2109–18.
Chosa N, Ishisaki A. Two novel mechanisms for maintenance of stemness in mesenchymal stem cells: SCRG1/BST1 axis and cell-cell adhesion through N-cadherin. Jpn Dent Sci Rev. 2018;54(1):37–44.
Dron M, Bailly Y, Beringue V, Haeberlé AM, Griffond B, Risold PY, et al. Scrg1 is induced in TSE and brain injuries, and associated with autophagy. Eur J Neurosci. 2005;22(1):133–46.
Dron M, Bailly Y, Beringue V, Haeberlé AM, Griffond B, Risold PY, et al. SCRG1, a potential marker of autophagy in transmissible spongiform encephalopathies. Autophagy. 2006;2(1):58–60.
Abdollahi Nejat M, Klaassen RV, Spijker S, Smit AB. Auxiliary subunits of the AMPA receptor: The Shisa family of proteins. Curr Opin Pharmacol. 2021;58:52–61.
Mu P, Akashi T, Lu F, Kishida S, Kadomatsu K. A novel nuclear complex of DRR1, F-actin and COMMD1 involved in NF-κB degradation and cell growth suppression in neuroblastoma. Oncogene. 2017;36(41):5745–56.
Fernández Massó JR, Oliva Argüelles B, Tejeda Y, Astrada S, Garay H, Reyes O, et al. The antitumor peptide CIGB-552 increases COMMD1 and inhibits growth of human lung cancer cells. J Amino Acids. 2013;2013:251398.
Tan AC, Ashley DM, López GY, Malinzak M, Friedman HS, Khasraw M. Management of glioblastoma: state of the art and future directions. CA Cancer J Clin. 2020;70(4):299–312.
Ceccarelli M, Barthel FP, Malta TM, Sabedot TS, Salama SR, Murray BA, et al. Molecular profiling reveals biologically discrete subsets and pathways of progression in diffuse glioma. Cell. 2016;164(3):550–63.
Brennan CW, Verhaak RG, McKenna A, Campos B, Noushmehr H, Salama SR, et al. The somatic genomic landscape of glioblastoma. Cell. 2013;155(2):462–77.
Hegi ME, Diserens AC, Gorlia T, Hamou MF, de Tribolet N, Weller M, et al. MGMT gene silencing and benefit from temozolomide in glioblastoma. N Engl J Med. 2005;352(10):997–1003.
Kristensen BW, Priesterbach-Ackley LP, Petersen JK, Wesseling P. Molecular pathology of tumors of the central nervous system. Ann Oncol. 2019;30(8):1265–78.
Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, et al. Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010;17(1):98–110.
Hinshaw DC, Shevde LA. The tumor microenvironment innately modulates cancer progression. Cancer Res. 2019;79(18):4557–66.
Wu T, Dai Y. Tumor microenvironment and therapeutic response. Cancer Lett. 2017;387:61–8.
Lei X, Lei Y, Li JK, Du WX, Li RG, Yang J, et al. Immune cells within the tumor microenvironment: biological functions and roles in cancer immunotherapy. Cancer Lett. 2020;470:126–33.
The present study was supported by the National Natural Science Foundation of China (Grant Number 8167050561). The institutions that provided funding did not participate in the data analysis, article writing, and submission decision of this paper.
Ethics approval and consent to participate
All methods were performed in accordance with the relevant guidelines and regulations.
Consent for publication
The authors declare that there are no competing interests associated with the study.
About this article
Cite this article
Han, J., Jing, Y., Han, F. et al. Comprehensive analysis of expression, prognosis and immune infiltration for TIMPs in glioblastoma. BMC Neurol 21, 447 (2021). https://doi.org/10.1186/s12883-021-02477-1
- Bioinformatics analysis
- Tissue inhibitors of metalloproteinases family
- Immune infiltration