- Primary research
- Open Access
The expression patterns and the diagnostic/prognostic roles of PTPN family members in digestive tract cancers
Cancer Cell International volume 20, Article number: 238 (2020)
Non-receptor protein tyrosine phosphatases (PTPNs) are a set of enzymes involved in the tyrosyl phosphorylation. The present study intended to clarify the associations between the expression patterns of PTPN family members, and diagnosis as well as the prognosis of digestive tract cancers.
Oncomine and Ualcan were used to analyze PTPN expressions. Data from The Cancer Genome Atlas (TCGA) were downloaded through UCSC Xena for validation and to explore the relationship of the PTPN expression with diagnosis, clinicopathological parameters and survival of digestive tract cancers. Gene ontology enrichment analysis was conducted using the DAVID database. The gene–gene interaction network was performed by GeneMANIA and the protein–protein interaction (PPI) network was built using STRING portal coupled with Cytoscape. The expression of differentially expressed PTPNs in cancer cell lines were explored using CCLE. Moreover, by histological verification, the expression of four PTPNs in digestive tract cancers were further analyzed.
Most PTPN family members were associated with digestive tract cancers according to Oncomine, Ualcan and TCGA data. Several PTPN members were differentially expressed in digestive tract cancers. For esophageal carcinoma (ESCA), PTPN1 and PTPN12 levels were correlated with incidence; PTPN20 was associated with poor prognosis. For stomach adenocarcinoma (STAD), PTPN2 and PTPN12 levels were correlated with incidence; PTPN3, PTPN5, PTPN7, PTPN11, PTPN13, PTPN14, PTPN18 and PTPN23 were correlated with pathological grade; PTPN20 expression was related with both TNM stage and N stage; PTPN22 was associated with T stage and pathological grade; decreased expression of PTPN5 and PTPN13 implied worse overall survival of STAD, while elevated PTPN6 expression indicated better prognosis. For colorectal cancer (CRC), PTPN2, PTPN21 and PTPN22 levels were correlated with incidence; expression of PTPN5, PTPN12, and PTPN14 was correlated with TNM stage and N stage; high PTPN5 or PTPN7 expression was associated with increased hazards of death. CCLE analyses showed that in esophagus cancer cell lines, PTPN1, PTPN4 and PTPN12 were highly expressed; in gastric cancer cell lines, PTPN2 and PTPN12 were highly expressed; in colorectal cancer cell lines, PTPN12 was highly expressed while PTPN22 was downregulated. Results of histological verification experiment showed differential expressions of PTPN22 in CRC, and PTPN12 in GC and CRC.
Members of PTPN family were differentially expressed in digestive tract cancers. Correlations were found between PTPN genes and clinicopathological parameters of patients. Expression of PTPN12 was upregulated in both STAD and CRC, and thus could be used as a diagnostic biomarker. Differential expression of PTPN12 in GC and CRC, and PTPN22 in CRC were presented in our histological verification experiment.
Tyrosyl phosphorylation is a dynamic and reversible process which plays a pivotal part in many cellular signaling pathways . The dephosphorylation of tyrosine residues are catalyzed by a series of enzymes named protein tyrosine phosphatases (PTPs) . Encoded by 103 genes, PTPs are sorted into four major superfamily classes , and every single PTP was denominated an official gene name by The Human Genome Organization in Nomenclature Committee. According to the latter system, 17 non-receptor PTPs which belong to the biggest family class I, are designated PTPN, followed by a number . There is mounting evidence suggesting that the cross-talk of the PTPN gene family members is involved in extensive physiological processes, such as cell proliferation, survival, immune response, migration, and metabolism [3,4,5,6]. Previously published study stated that PTPN family members play an essential part in numerous diseases. For example, the expression of PTPN6 with the loss of pSTAT3 expression could be chosen as a biomarker for the prognosis of peripheral-T cell lymphoma . One research published recently noted that a deficiency of PTPN2 could enhance anti-tumor immunity and the therapeutic efficacy of CAR T cells to solid cancers . PTPN22 plays an important part in regulating autophagy and NLR family pyrin domain containing 3 inflammasome activation . It is hopeful that PTPN genes have potential to be served as prognostic and diagnostic indicators , and even therapeutic targets.
As main malignancies of gastrointestinal tract, esophagus cancer (EC), gastric cancer (GC) and colorectal cancer (CRC) are responsible for a large portion of cancer-related deaths worldwide, of which the data are 5.3%, 8.2% and 9.0% respectively . All of the three digestive tract cancers are ranked in the top 10 for incidence rates of tumors and have poor prognosis . Although a couple of diagnostic biomarkers have been observed, robust biomarkers to predict clinical outcomes are still in urgent need . Former investigations reported that the expression patterns of individual PTPNs and its correlations with patients in various digestive tract neoplasms, but research to date has not yet observed the whole picture of the entire PTPN family, from the aspects of diagnosis, prognosis and expression characteristics. Personalized therapies based upon the genetics of individual cancer will be prior to other treatments in the near future. This research therefore is intended to illuminate the clinical value of different PTPN genes to support potential biomarkers and new individualized targets for patients. In the current study, we analyzed the expression status of different PTPN members and their diagnostic and prognostic values to comprehensively evaluate the role of PTPNs in various digestive tract cancers.
Oncomine database analysis for the expression patterns of PTPN family in digestive tract cancers
Oncomine database is a cancer microarray database (https://www.oncomine.org/) which shows the expression information of genes in cancer and normal samples . Oncomine provides both microarray information from 715 datasets and a set of online data-mining functions. The expression levels of individual PTPN family members in different types of cancer were obtained from the Oncomine database. Student’s t test was applied to calculate the P value for expression differences of PTPN family genes between normal controls and cancer samples. The threshold parameters of P value and fold change were demarcated as 0.05 and 2 respectively.
Ualcan database analysis for the validation
Ualcan is a publicly available web-portal (http://ualcan.path.uab.edu) that offers online analysis of data from The Cancer Genome Atlas (TCGA) . In this study, we used it to analyze the relative expression of PTPN family genes in esophageal carcinoma (ESCA), stomach adenocarcinoma (STAD), colon adenocarcinoma (COAD), rectum adenocarcinoma (READ) and normal samples. The expression level of PTPN members was normalized as transcript per million reads, and a P value of no more than 0.01 conducted through Student’s t test was considered to be significant.
TCGA data obtainment through UCSC Xena
The Cancer Genome Atlas (TCGA) is a cancer genomics database (https://cancergenome.nih.gov/) which contains genomic information of over 2000 primary neoplasms and matched normal samples [15, 16]. In our study, case information concerning the mRNA expression profiles and clinical features were obtained from TCGA database, which was downloaded through the University of California Santa Cruz Xena (UCSC Xena; https://xena.ucsc.edu/) platform . Data of 184 ESCA and 11 matched normal samples, 415 STAD and 35 matched normal samples, and 380 CRC and 51 matched normal samples were extracted for further analysis.
Analysis of PTPN members in diagnosis and prognosis of digestive tract cancers
Receiver operating characteristic (ROC) models were established to evaluate the diagnostic values of the differentially expressed PTPNs according to Oncomine and Ualcan analysis. SPSS software version 23.0 (IBM, SPSS, and Chicago, IL, USA) was applied to plot ROC curves. P < 0.01 was considered significant. The Pearson X2 test was performed to evaluate the relationship between PTPN expression status and clinicopathological parameters including TNM stage and grade. All of the patients were differentiated into low expression and high expression groups using the third quartile of mRNA expressions as the cut-off value. The effects of different PTPN expressions on overall survival (OS) were estimated through univariate and multivariate Cox proportional hazards models with or without adjustment for confounding factors. Variables including gender, TNM stage, and age were further adjusted during the evaluation. The medical data obtained was managed by R language (Version 3.6.1). P value < 0.05 was considered statistically significant.
GO and PPI analysis for function and interaction of PTPN family
Enrichment analysis of Gene Ontology (GO) of PTPN genes was explored using the Database for Annotation, Visualization and Integrated Discovery (DAVID; v.6.8; https://david.ncifcrf.gov/home.jsp; accessed on November 20, 2019) . The gene–gene interaction network was structured using the Gene Multiple Association Network Integration Algorithm (GeneMANIA; https://www.genemania.org/; accessed on November 21, 2019)  and the Search Tool for the Retrieval of Interacting Genes Database (STRING v.10.0; https://string-db.org/; accessed on November 23, 2019) was used to create a protein–protein interaction (PPI) network . The Cytoscape software was applied to visualize network diagrams for PPI analysis .
Evaluation of PTPN expression in cancer cell lines by CCLE analysis
The Broad Institute Cancer Cell Line Encyclopedia (CCLE; http://www.broadinstitute.org/ccle) is a database, which aims to provide available access to detailed genomic information, computational analyses and visualization for over 1100 cell lines representing 37 cancer types . The expression levels of PTPN1, PTPN2, PTPN4, PTPN12, PTPN21 and PTPN22 in different types of cancer cell lines were investigated through CCLE database.
Collection of specimens and histological verification experiment
All clinical specimens used to measure mRNA levels of PTPN1, PTPN2, PTPN12 and PTPN22 were collected from the First Hospital of China Medical University, Shenyang, China. The study was approved by the ethics committee of the First Hospital of China Medical University. The written informed consent has been obtained from each participant before specimen collection. Specimens were collected in accordance with the Declaration of Helsinki and legal regulations.
Total RNA was extracted from 30 pairs of patient tissues (10 EC, 10 GC and 10 CRC) and adjacent non-tumor tissues. Relative mRNA expression levels were detected by quantitative real-time PCR (qRT-PCR).The qRT-PCR experiment was performed with a real-time PCR 480 system and SYBR-green PCR Master Mix. All of the qRT-PCR curves were with single peak. P < 0.05 was considered to imply significant results. Primer sequences for PTPN members and beta-actin were tabulated in Additional file 1: Table S1.
PTPN genes expression patterns in digestive tract cancers
Before expression profiling analysis, we refined the chromosome location of all PTPN members through published literature review. The detailed information was summarized in Table 1. All of the PTPN family members have been analyzed in our study. According to the analysis results of Oncomine, the expression of PTPN genes were different in all types of cancer and its matched normal tissues (Fig. 1). For EC samples, PTPN1, PTPN4, PTPN12, PTPN18 and PTPN21 were over-expressed, while PTPN3, PTPN11, PTPN13 and PTPN21 were downregulated. In GC tissues, the expression of PTPN5 and PTPN13 was decreased while at the same time PTPN2, PTPN12 and PTPN22 were highly expressed in patients. As for CRC patients, the expression of PTPN1, PTPN3, PTPN7, PTPN11, PTPN12, PTPN13 and PTPN14 was upregulated, while PTPN2, PTPN18, PTPN20, PTPN21 and PTPN22 were expressed in a lower level in tissues. Comprehensively, the expression of PTPN12 was higher in EC, GC and CRC samples.
Then we used Ualcan database to further verify the expression of PTPNs among ESCA, STAD, COAD, READ and their matched normal tissues in the TCGA datasets. As shown in the results of Ualcan (Figs. 2, 3, 4, 5), there was a clear trend that the expression levels of PTPN1, PTPN2, PTPN4, PTPN6, PTPN7, PTPN9, PTPN12, PTPN22 and PTPN23 were statistically higher in samples of ESCA than in normal tissues; in STAD, PTPN1, PTPN2, PTPN3, PTPN4, PTPN6, PTPN7, PTPN9, PTPN11, PTPN12, PTPN18, PTPN22, and PTPN23 were significantly increased compared to normal samples; in COAD, PTPN1, PTPN2, PTPN4, PTPN11, PTPN12, PTPN13, and PTPN23 were overexpressed compared with normal samples, while PTPN9, PTPN18, PTPN21 and PTPN22 were lower than normal tissues; in READ, the expression levels of PTPN4, PTPN12, and PTPN23 were higher than normal samples while PTPN21 and PTPN22 were decreased in patients compared to normal groups. Taken together, the results above indicated that expression status of PTPN4, PTPN12 and PTPN23 was upregulated in ESCA, STAD, COAD and READ. PTPN21 and PTPN22 were upregulated in ESCA and STAD but downregulated in COAD and READ.
The diagnostic efficiency of PTPN genes in digestive tract cancers
According to the expression analysis of Oncomine and Ualcan above, some PTPN genes were over-expressed in digestive tract cancers, while others were downregulated. Therefore, we computed ROC curves to further analyze the diagnostic efficiency of these differentially expressed PTPNs. The results revealed that upregulation of PTPN1 (Area = 0.784, 95% CI 0.675–0.893, and P < 0.002) and PTPN12 (Area = 0.860, 95% CI 0.768–0.951, and P < 0.001) was correlated with ESCA incidence; upregulation of PTPN2 (Area = 0.839, 95% CI 0.771–0.906, and P < 0.001) and PTPN12 (Area = 0.775, 95% CI 0.680–0.871, and P < 0.001) was associated with STAD incidence; and upregulation of PTPN12 (Area = 0.946, 95% CI 0.921–0.970, and P < 0.001), and downregulation of PTPN21 (Area = 0.986, 95% CI 0.977–0.995, and P < 0.001) and PTPN22 (Area = 0.764, 95% CI 0.704–0.825, and P < 0.001) were associated with CRC incidence (Table 2). Detailed results of the ROC curves are shown in Additional file 2. Our results above suggest that levels of these genes could be exploited as effective biomarkers for the diagnosis of specific digestive tract cancers. Especially, PTPN12 might serve as a useful diagnostic biomarker for both ESCA, STAD and CRC.
The correlations between PTPN genes and clinicopathological parameters
Association of PTPN genes expression status with different clinicopathological features of ESCA, STAD and CRC were analyzed in the study. For ESCA, no statistically association was found between PTPN members and clinicopathological parameters (Additional file 1: Table S2). For STAD, grade was significantly correlated with PTPN3, PTPN11, PTPN13, PTPN14, and PTPN18 (P = 0.015, 0.021, 0.005, 0.002 and 0.034, resp.; Additional file 1: Table S3). In addition, PTPN5, PTPN7 and PTPN23 showed statistically association with pathological grade of STAD with a p value of no more than 0.001. Besides, PTPN20 expression was related with TNM stage (P = 0.002) and N stage (P = 0.038); PTPN22 was significantly related with grade (P < 0.001) and T stage (P = 0.039) (Additional file 1: Table S3). For CRC, TNM stage was correlated with expression of PTPN5, PTPN12, and PTPN14 (P = 0.049, 0.027, 0.002, resp.; Additional file 1: Table S4); N stage was associated with expression of PTPN5, PTPN12, and PTPN14 (P = 0.034, 0.018 and 0.004, resp.; Additional file 1: Table S4). According to the results above, PTPN5 and PTPN14 were associated with the clinicopathological parameters of both STAD and CRC.
The prognostic roles of PTPN genes in digestive tract cancers
The analysis results of the relationship between PTPN family genes and the prognosis of digestive tract cancers showed that for ESCA, elevated expression of PTPN20 was associated with a worse overall survival (OS) of ESCA in multivariate model (P = 0.038, adjusted HR = 1.982, 95% CI 1.039–3.780; Table 3). For STAD, decreased PTPN5 expression was associated with worse OS according to univariate analysis (P = 0.006, HR = 1.680, 95% CI 1.162–2.429) and multivariate analysis (P = 0.004, adjusted HR = 1.753, 95% CI 1.192–2.577), while elevated expression of PTPN6 indicated longer OS based on the results of univariate analysis (P = 0.028, HR = 0.633, 95% CI 0.421–0.951) and multivariate analysis (P = 0.017, adjusted HR = 0.606, 95% CI 0.402–0.915) (Table 4). Besides, univariate analysis exhibited that there was a significant association between decreased PTPN13 expression and decreased OS (P = 0.036, HR = 1.471, 95% CI 1.026–2.109). CRC patients with high PTPN5 expression showed increased hazards of death in univariate model (P = 0.021, HR = 1.947, 95% CI 1.105–3.431; Table 5). Elevated expression of PTPN7 was significantly associated with the unfavorable OS for every CRC patients according to multivariate analysis (P = 0.013, adjusted HR = 2.043, 95% CI 1.164–3.584; Table 5). From the results above we can see that the expression of PTPN5 was associated with the prognosis of both STAD and CRC.
Function and interaction of PTPN family genes
GO analysis was basically grouped into three terms including molecular function groups, cellular component groups and biological process groups. The top 5 enriched categories obtained from the analysis results of each group were showed in Fig. 6a. GO analysis revealed that PTPN proteins were mainly related to cytoplasmic side of plasma membrane. PTPN genes exert their functions primarily on peptidyl-tyrosine dephosphorylation and protein tyrosine phosphatase activity. Further, the interaction analysis of PTPN genes at the gene level was performed by GeneMANIA to clarify the correlations among colocalization, shared protein domains, co-expression, prediction and pathways (Fig. 6b). As the protein–protein interaction network of STRING analysis revealed, interrelationships among PTPN gene family members were intricate (Fig. 6c).
Differential expression of PTPNs in cancer cell lines
We used CCLE to confirm the expression of PTPN1, PTPN2, PTPN4, PTPN12, PTPN21 and PTPN22 in digestive tract cancer cell lines. All of these six PTPNs were differentially expressed in cancer and control groups according to both Oncomine and Ualcan database analysis results. As the Additional file 3 showed, levels of PTPN1, PTPN4 and PTPN12 were higher in esophageal cancer cell lines; PTPN2 and PTPN12 expressions were higher in stomach cancer cell lines; higher PTPN12 and decreased PTPN22 were identified in colorectal cancer cell lines. These results were consistent with our analyses results of Oncomine and Ualcan database. Different from these results, PTPN22 was downregulated in stomach cancer cell lines and PTPN21 were increased in colorectal cancer cell lines, which were inconsistent with the conclusions of Oncomine and Ualcan database. The expression profiles of PTPN21 in colorectal cancer and PTPN22 in stomach cancer warrant further examination.
Verification of four PTPN expressions in digestive tract cancer tissues
To confirm our analysis results, we then did tissue verification with qRT-PCR experiment. The results revealed that in GC, PTPN12 expression were higher than normal tissues (P = 0.037); in CRC, higher PTPN12 expression and lower PTPN22 expression were observed (P = 0.007 and 0.009, resp.); and marginal difference was found with the expression of PTPN1 in EC (P = 0.074) (Table 6).
Being significant members of protein tyrosine phosphatases, PTPN regulates tyrosine phosphorylation and dephosphorylation in cellular signal transduction together with protein-tyrosine kinases . There are abundant investigations focusing on the links between individual PTPN family members and diverse neoplasms. However, no report has offered an overview of how PTPN family genes tie up with diverse gastrointestinal tumors up to now. This research for the first time investigated the expression patterns of PTPN genes, and illustrated the association between PTPN family genes and the diagnosis as well as prognosis of digestive tract cancers, which provides deeper insights into the clinical values of all the PTPN genes in various digestive tract cancers.
We used online database Oncomine and Ualcan to explore the PTPN genes expression patterns in digestive tract cancers, ROC models were further established to investigate their diagnostic values. CCLE analysis and human tissue experiment by qRT-PCR were conducted to verify our results. As revealed by our study, expression levels of PTPN5 and PTPN13 were decreased, while PTPN2, PTPN12, and PTPN22 were increased in human STAD. PTPN2 and PTPN12 were further exploited as STAD diagnostic biomarkers. According to previous studies, PTPN2 antagonizes cytokine signaling needed for T cell differentiation, homeostasis, and function by regulating the dephosphorylation and inactivation of Janus-activated kinase (JAK)-1 and JAK-3, and their target substrates signal transducer and activator of transcription (STAT)-1, STAT-3 and STAT-5 in a cell context-dependent manner [23,24,25,26]. Recent researches demonstrated that PTPN2 deletion can enhance the effectiveness of anti-tumor immunity [8, 27]. Thus, PTPN2 expression is likely to be exploited as an attractive immunotherapy target in the treatment of cancer, such as STAD. For CRC patients, PTPN3, PTPN7, PTPN11, PTPN12, PTPN13 and PTPN14 were upregulated while PTPN2, PTPN18, PTPN21, and PTPN22 expressions were downregulated according to Oncomine. PTPN12, PTPN21, and PTPN22 were further verified in Ualcan, and proved to be correlated with CRC incidence. PTPN3 was demonstrated to be served as a tumor suppressor gene in CRC . While another study illustrated that elevated PTPN3 expression promotes tumor recurrence and is detrimental to the prognosis of intrahepatic cholangiocarcinoma patients . These findings indicate that regulation mechanisms of PTPN3 in various cancers are complicated. Consistent with our results, Slattery et al.  reported that PTPN11 expression is upregulated in CRC. PTPN11 was suggested to impact the tumorigenesis and metastasis of CRC  and termed as a potential prognostic marker . The data in one report showed that in patients with gastric B cell non-Hodgkin’s lymphoma, PTPN21 is over expressed . In ESCA, the expression levels of PTPN1, PTPN4, and PTPN12 were increased. Differential expression of PTPN12 in GC and CRC was observed in our histological experiment. ROC results suggested that PTPN12 could serve as a diagnostic biomarker gastrointestinal cancers. In several researches, PTPN12 is characterized as a tumor suppressor which antagonizes EGFR/HER2 signaling [34, 35], which is contrary to our findings. In hepatocellular carcinoma cells, PTPN12 regulates epithelial–mesenchymal transition which contributes to chemoresistance and metastasis . All of these indicate that functions of PTPN12 are different in various cancer types, and comprehensive studies are required to clarify its specific mechanisms.
There may also exist associations between differentially expressed PTPNs and prognosis of digestive tract cancers. To further figure out the associations, we then analyzed the interrelationships between PTPN expressions and the clinical outcomes of digestive tract cancers. As revealed by the current study, PTPN5, PTPN13 and PTPN22 were also associated with clinicopathological parameters of STAD. And decreased expression of PTPN5 and PTPN13 indicated worse OS of STAD patients, while high PTPN6 expression was associated with a favorable STAD OS. It is worth noting that SHP-1 protein encoded by PTPN6 mediates the tumor-suppressive function of TMEFF2 in STAD , which indicates expression of PTPN6 might influence the carcinogenesis of STAD patients. PTPN13 has been reported to regulate the resistance of human lung fibroblasts to Fas-induced apoptosis in previous study . Mutated PTPN13 was suggested to be a tumor suppressor gene in colorectal cancer . PTPN22 gene encodes an enzyme called lymphoid-specific tyrosine phosphatase, which functions as a master regulator in the biological process of relevant immune responses . Upregulation of PTPN22 could result in impairment of regulatory T-cell differentiation in patients with non-ST-segment elevation acute coronary syndromes . As for CRC, expression levels of PTPN5, PTPN12, and PTPN14 were correlated with clinicopathological parameters. Besides, elevation of PTPN5 and PTPN7 was significantly associated with increased death hazards of CRC patients respectively. One study presents PTPN12 as a novel candidate that contributes to the heterogeneous susceptibility to colorectal cancer . PTPN5 could regulate tyrosine dephosphorylation needed for the activation of BAK, a noteworthy cell-death mediator in apoptosis . It was once reported that PTPN14 regulated phosphorylation of p130Cas Y128 plays a crucial role in colorectal carcinogenesis . And mutated PTPN14 is suggested to be a tumor suppressor gene for colorectal cancer, regulating cellular pathways that are appropriate for therapeutic intervention . A recent study proved that PTPN14 suffices to inhibit migration and invasion of metastatic cancer cells . In a xenograft breast cancer model, PTPN14 acts as a suppressor of metastasis of triple-negative breast cancer cells . In our analysis, PTPN5 and PTPN14 expression was associated with the clinicopathological parameters of both STAD and CRC. However, differential expressions of PTPN5 in GC and controls, and PTPN14 in CRC and controls were observed only in Oncomine database, but were not further verified in Ualcan database. Taken together, we hypothesized that PTPN5 and PTPN14 mainly function in the process of tumor progression, and its correlation with the risk of cancer requires further evaluation and experiments. As for ESCA, our analysis indicated that only PTPN20 expression was related with a worse OS of patients. However, former research indicated that PTPN12 may serve as a potential prognostic indicator for esophagus cancer patients . In our study, some PTPNs were associated with overall survival of cancer patients, but not associated with clinicopathological parameters. As we all know, variables including patient-level, tumor-level and environment-level factors all exert great influence on survival outcomes of cancer patients, which have been illustrated by multiple laboratory evidence and clinical experience. Thus it can be ratiocinated that intrinsic interactions may exist between PTPNs and other clinical factors and consequently affect the prognosis of cancer. In order to understand the functions of PTPNs in gastrointestinal cancers better, investigations concerning more clinicopathological parameters are requisite in the future.
In this study, GO analysis revealed that PTPN proteins were mainly related to cytoplasmic side of plasma membrane. PTPN genes exert their functions primarily on peptidyl-tyrosine dephosphorylation and protein tyrosine phosphatase activity, which was illustrated by dozens of publications. Substrate-trapping and biochemical analyses demonstrated that PTPN22 mainly dephosphorylates the tyrosine residues of SFKs and SFK substrates including the E3ubiquitin ligase c-Cbl57-59 and the T cell signaling protein Zap70 [47,48,49]. PTPN5 was reported to mediate internalization and phosphorylation of AMPA receptors after metabotropic glutamate receptor stimulation . In addition, the results of interaction network analysis at gene and protein levels further indicated that PTPN members and other genes comprehensively interact with each other. One experiment suggested that both PTPN1 and PTPN2 target protein mediator of IRF3 activation for dephosphorylation at Y245 . PTPN2 and PTPN22 were demonstrated negatively regulating T cell receptor signaling by dephosphorylating lymphocyte-specific protein tyrosine kinase [47, 52]. Multiple efforts demonstrated that through binding and inactivating the mitogen activated protein kinase Erk2 and p38, PTPN5 and PTPN7 could negatively regulate cell proliferation as well as differentiation [53,54,55,56]. Promoter hypermethylation of PTPN6 and PTPN13 was reported to inhibit the progression of diffuse large B cell lymphomas . In our study, PTPN5 and PTPN7 were found to be correlated with prognosis of colon adenocarcinoma patients, while PTPN6 and PTPN13 were statistically associated with the prognosis of STAD. All of these implied that PTPN members could function through alliance mechanisms in many diseases, including gastrointestinal cancers.
In summary, findings of our study illustrated the expression status as well as diagnostic and prognostic values of PTPN members in digestive tract cancers. The results indicated that several PTPN members were differentially expressed and related to clinical outcomes of patients with digestive tract cancers. Especially, upregulation of PTPN12 was correlated with the incidence of ESCA, STAD and CRC. Differential expression of PTPN12 in GC and CRC, and PTPN22 in CRC were presented in our histological verification experiment. Future well-designed investigations are required to elucidate the significance of our findings and thus develop the clinical utility of PTPNs.
Availability of data and materials
The datasets analyzed during the current study are available in the TCGA repository (https://cancergenome.nih.gov/). The authors declare that the data supporting the findings of this study are available within the article.
Protein tyrosine phosphatases
Non-receptor protein tyrosine phosphatase
The Cancer Genome Atlas
Receiver operating characteristics
Database for Annotation, Visualization and Integrated Discovery
Gene Multiple Association Network Integration Algorithm
Search Tool for the Retrieval of Interacting Genes Database
Broad Institute Cancer Cell Line Encyclopedia
Signal transducer and activator of transcription
Li X, Wilmanns M, Thornton J, Kohn M. Elucidating human phosphatase-substrate networks. Science signaling. 2013;6(275):10.
Hunter T. Tyrosine phosphorylation: thirty years and counting. Curr Opin Cell Biol. 2009;21(2):140–6.
Alonso A, Sasin J, Bottini N, Friedberg I, Friedberg I, Osterman A, et al. Protein tyrosine phosphatases in the human genome. Cell. 2004;117(6):699–711.
Tonks NK. Protein tyrosine phosphatases: from genes, to function, to disease. Nat Rev Mol Cell Biol. 2006;7(11):833–46.
Julien SG, Dube N, Hardy S, Tremblay ML. Inside the human cancer tyrosine phosphatome. Nat Rev Cancer. 2011;11(1):35–49.
Yu ZH, Zhang ZY. Regulatory mechanisms and novel therapeutic targeting strategies for protein tyrosine phosphatases. Chem Rev. 2018;118(3):1069–91.
Han JJ, O’Byrne M, Stenson MJ, Maurer MJ, Wellik LE, Feldman AL, et al. Prognostic and therapeutic significance of phosphorylated STAT3 and protein tyrosine phosphatase-6 in peripheral-T cell lymphoma. Blood Cancer J. 2018;8(11):110.
Wiede F, Lu KH, Du X, Liang S, Hochheiser K, Dodd GT, et al. PTPN2 phosphatase deletion in T cells promotes anti-tumour immunity and CAR T-cell efficacy in solid tumours. EMBO J. 2019;39:e103637.
Spalinger MR, Lang S, Gottier C, Dai X, Rawlings DJ, Chan AC, et al. PTPN22 regulates NLRP3-mediated IL1B secretion in an autophagy-dependent manner. Autophagy. 2017;13(9):1590–601.
Zhang ZY. Drugging the undruggable: therapeutic potential of targeting protein tyrosine phosphatases. Acc Chem Res. 2017;50(1):122–9.
Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J Clin. 2018;68(6):394–424.
Taniguchi H, Moriya C, Igarashi H, Saitoh A, Yamamoto H, Adachi Y, et al. Cancer stem cells in human gastrointestinal cancer. Cancer Sci. 2016;107(11):1556–62.
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 (New York, NY). 2004;6(1):1–6.
Chandrashekar DS, Bashel B, Balasubramanya SAH, Creighton CJ, Ponce-Rodriguez I, Chakravarthi B, et al. UALCAN: a portal for facilitating tumor subgroup gene expression and survival analyses. Neoplasia (New York, NY). 2017;19(8):649–58.
Weinstein JN, Collisson EA, Mills GB, Shaw KR, Ozenberger BA, Ellrott K, et al. The Cancer Genome Atlas Pan-Cancer analysis project. Nat Genet. 2013;45(10):1113–20.
Akbani R, Ng PK, Werner HM, Shahmoradgoli M, Zhang F, Ju Z, et al. A pan-cancer proteomic perspective on The Cancer Genome Atlas. Nat Commun. 2014;5:3887.
Goldman M, Craft B, Swatloski T, Cline M, Morozova O, Diekhans M, et al. The UCSC Cancer Genomics Browser: update 2015. Nucleic Acids Res. 2015;43(Database issue):D812–7.
da Huang W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
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.
Szklarczyk D, Morris JH, Cook H, Kuhn M, Wyder S, Simonovic M, et al. The STRING database in 2017: quality-controlled protein-protein association networks, made broadly accessible. Nucleic Acids Res. 2017;45(D1):D362–8.
Su G, Morris JH, Demchak B, Bader GD. Biological network exploration with Cytoscape 3. Curr Protoc Bioinform. 2014;47:8.13.1-24.
Barretina J, Caponigro G, Stransky N, Venkatesan K, Margolin AA, Kim S, et al. The cancer cell line encyclopedia enables predictive modelling of anticancer drug sensitivity. Nature. 2012;483(7391):603–7.
ten Hoeve J, de Jesus Ibarra-Sanchez M, Fu Y, Zhu W, Tremblay M, David M, et al. Identification of a nuclear Stat1 protein tyrosine phosphatase. Mol Cell Biol. 2002;22(16):5662–8.
Simoncic PD, Lee-Loy A, Barber DL, Tremblay ML, McGlade CJ. The T cell protein tyrosine phosphatase is a negative regulator of janus family kinases 1 and 3. Curr Biol CB. 2002;12(6):446–53.
Wiede F, Dudakov JA, Lu KH, Dodd GT, Butt T, Godfrey DI, et al. PTPN2 regulates T cell lineage commitment and alphabeta versus gammadelta specification. J Exp Med. 2017;214(9):2733–58.
Wiede F, Sacirbegovic F, Leong YA, Yu D, Tiganis T. PTPN2-deficiency exacerbates T follicular helper cell and B cell responses and promotes the development of autoimmunity. J Autoimmun. 2017;76:85–100.
LaFleur MW, Nguyen TH, Coxe MA, Miller BC, Yates KB, Gillis JE, et al. PTPN2 regulates the generation of exhausted CD8(+) T cell subpopulations and restrains tumor immunity. Nat Immunol. 2019;20(10):1335–47.
Wang Z, Shen D, Parsons DW, Bardelli A, Sager J, Szabo S, et al. Mutational analysis of the tyrosine phosphatome in colorectal cancers. Science (New York, NY). 2004;304(5674):1164–6.
Gao Q, Zhao YJ, Wang XY, Guo WJ, Gao S, Wei L, et al. Activating mutations in PTPN3 promote cholangiocarcinoma cell proliferation and migration and are associated with tumor recurrence in patients. Gastroenterology. 2014;146(5):1397–407.
Slattery ML, Herrick JS, Mullany LE, Samowitz WS, Sevens JR, Sakoda L, et al. The co-regulatory networks of tumor suppressor genes, oncogenes, and miRNAs in colorectal cancer. Genes Chromosom Cancer. 2017;56(11):769–87.
Zhang Q, Li Y, Zhao R, Wang X, Fan C, Xu Y, et al. The gain-of-function mutation E76K in SHP2 promotes CAC tumorigenesis and induces EMT via the Wnt/β-catenin signaling pathway. Mol Carcinog. 2018;57(5):619–28.
Yu SJ, Yu JK, Ge WT, Hu HG, Yuan Y, Zheng S. SPARCL1, Shp2, MSH2, E-cadherin, p53, ADCY-2 and MAPK are prognosis-related in colorectal cancer. World J Gastroenterol. 2011;17(15):2028–36.
Plani-Lam JH, Chow TC, Fan YH, Garcia-Bloj B, Cheng L, Jin DY, et al. High expression of PTPN21 in B-cell non-Hodgkin’s gastric lymphoma, a positive mediator of STAT5 activity. Blood Cancer J. 2016;6:e388.
Sun T, Aceto N, Meerbrey KL, Kessler JD, Zhou C, Migliaccio I, et al. Activation of multiple proto-oncogenic tyrosine kinases in breast cancer via loss of the PTPN12 phosphatase. Cell. 2011;144(5):703–18.
Zheng Y, Zhang C, Croucher DR, Soliman MA, St-Denis N, Pasculescu A, et al. Temporal regulation of EGF signalling networks by the scaffold protein Shc1. Nature. 2013;499(7457):166–71.
Kodama T, Newberg JY, Kodama M, Rangel R, Yoshihara K, Tien JC, et al. Transposon mutagenesis identifies genes and cellular processes driving epithelial-mesenchymal transition in hepatocellular carcinoma. Proc Natl Acad Sci USA. 2016;113(24):E3384–93.
Sun T, Du W, Xiong H, Yu Y, Weng I, Ren L. Correction: TMEFF2 deregulation contributes to gastric carcinogenesis and indicates poor survival outcome. Clin Cancer Res. 2015;21(15):3572.
Bamberg A, Redente EF, Groshong SD, Tuder RM, Cool CD, Keith RC, et al. Protein tyrosine phosphatase-N13 promotes myofibroblast resistance to apoptosis in idiopathic pulmonary fibrosis. Am J Respir Crit Care Med. 2018;198(7):914–27.
Vang T, Liu WH, Delacroix L, Wu S, Vasile S, Dahl R, et al. LYP inhibits T-cell activation when dissociated from CSK. Nat Chem Biol. 2012;8(5):437–46.
Flego D, Severino A, Trotta F, Previtero M, Ucci S, Zara C, et al. Increased PTPN22 expression and defective CREB activation impair regulatory T-cell differentiation in non-ST-segment elevation acute coronary syndromes. J Am Coll Cardiol. 2015;65(12):1175–86.
de Voer RM, Hahn MM, Weren RD, Mensenkamp AR, Gilissen C, van Zelst-Stams WA, et al. Identification of novel candidate genes for early-onset colorectal cancer susceptibility. PLoS Genet. 2016;12(2):e1005880.
Fox JL, Ismail F, Azad A, Ternette N, Leverrier S, Edelmann MJ, et al. Tyrosine dephosphorylation is required for Bak activation in apoptosis. EMBO J. 2010;29(22):3853–68.
Zhang P, Guo A, Possemato A, Wang C, Beard L, Carlin C, et al. Identification and functional characterization of p130Cas as a substrate of protein tyrosine phosphatase nonreceptor 14. Oncogene. 2013;32(16):2087–95.
Diaz-Valdivia NI, Diaz J, Contreras P, Campos A, Rojas-Celis V, Burgos-Ravanal RA, et al. The non-receptor tyrosine phosphatase type 14 blocks caveolin-1-enhanced cancer cell metastasis. Oncogene. 2020;39(18):3693–709.
Belle L, Ali N, Lonic A, Li X, Paltridge JL, Roslan S, et al. The tyrosine phosphatase PTPN14 (Pez) inhibits metastasis by altering protein trafficking. Science signaling. 2015;8(364):ra18.
Cao X, Li Y, Luo RZ, He LR, Yang J, Zeng MS, et al. Tyrosine-protein phosphatase nonreceptor type 12 is a novel prognostic biomarker for esophageal squamous cell carcinoma. Ann Thorac Surg. 2012;93(5):1674–80.
Cloutier JF, Veillette A. Cooperative inhibition of T-cell antigen receptor signaling by a complex between a kinase and a phosphatase. J Exp Med. 1999;189(1):111–21.
Wu S, Bottini M, Rickert RC, Mustelin T, Tautz L. In silico screening for PTPN22 inhibitors: active hits from an inactive phosphatase conformation. ChemMedChem. 2009;4(3):440–4.
Gjörloff-Wingren A, Saxena M, Williams S, Hammi D, Mustelin T. Characterization of TCR-induced receptor-proximal signaling events negatively regulated by the protein tyrosine phosphatase PEP. Eur J Immunol. 1999;29(12):3845–54.
Zhang Y, Venkitaramani DV, Gladding CM, Zhang Y, Kurup P, Molnar E, et al. The tyrosine phosphatase STEP mediates AMPA receptor endocytosis after metabotropic glutamate receptor stimulation. J Neurosci. 2008;28(42):10561–6.
Xia T, Yi XM, Wu X, Shang J, Shu HB. PTPN1/2-mediated dephosphorylation of MITA/STING promotes its 20S proteasomal degradation and attenuates innate antiviral response. Proc Natl Acad Sci USA. 2019;116(40):20063–9.
Wiede F, Shields BJ, Chew SH, Kyparissoudis K, van Vliet C, Galic S, et al. T cell protein tyrosine phosphatase attenuates T cell signaling to maintain tolerance in mice. J Clin Investig. 2011;121(12):4758–74.
Francis DM, Różycki B, Koveal D, Hummer G, Page R, Peti W. Structural basis of p38α regulation by hematopoietic tyrosine phosphatase. Nat Chem Biol. 2011;7(12):916–24.
Muñoz JJ, Tárrega C, Blanco-Aparicio C, Pulido R. Differential interaction of the tyrosine phosphatases PTP-SL, STEP and HePTP with the mitogen-activated protein kinases ERK1/2 and p38alpha is determined by a kinase specificity sequence and influenced by reducing agents. Biochem J. 2003;372(Pt 1):193–201.
Francis DM, Kumar GS, Koveal D, Tortajada A, Page R, Peti W. The differential regulation of p38α by the neuronal kinase interaction motif protein tyrosine phosphatases, a detailed molecular study. Structure. 2013;21(9):1612–23.
Francis DM, Koveal D, Tortajada A, Page R, Peti W. Interaction of kinase-interaction-motif protein tyrosine phosphatases with the mitogen-activated protein kinase ERK2. PLoS ONE. 2014;9(3):e91934.
Wang W, Wang J, Li Z, Zhu M, Zhang Z, Wang Y, et al. Promoter hypermethylation of PTPL1, PTPN6, DAPK, p16 and 5-azacitidine inhibits growth in DLBCL. Oncol Rep. 2016;35(1):139–46.
The authors gratefully acknowledge financial support from the National Natural Science Foundation of China (Ref No. 31800124).
This work was supported by grants from the National Natural Science Foundation of China (Ref No. 31800124).
Ethics approval and consent to participate
All participants involved in this study have provided written informed consent and the study protocols were permitted by the Ethics Committee of the First Hospital of China Medical University and performed in accordance with the Declaration of Helsinki.
Consent for publication
The authors declare that they have no competing interests in this work.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
The list of primers; Table S2. Association of PTPN family genes expression with clinicopathological parameters of ESCA; Table S3. Association of PTPN family genes expression with clinicopathological parameters of STAD; Table S4. Association of PTPN family genes expression with clinicopathological parameters of CRC.
Diagnostic value of differentially expressed PTPNs in digestive tract cancers. a ROCs of PTPN1/4/12 levels in ESCA, showing that elevated PTPN1/12 levels were correlated with ESCA incidence. b ROCs of PTPN2/12/22 levels in STAD, showing that elevated PTPN2/12 levels were correlated with STAD incidence. c ROCs of PTPN12/21/22 levels in CRC, showing that elevated PTPN12 levels and decreased PTPN21/22 levels were correlated with CRC incidence.
Expression profiles of certain PTPNs in cancer cell lines. (a) PTPN1 (b) PTPN2 (c) PTPN4 (d) PTPN12 (e) PTPN21 (f) PTPN22.
About this article
Cite this article
Chen, J., Zhao, X., Yuan, Y. et al. The expression patterns and the diagnostic/prognostic roles of PTPN family members in digestive tract cancers. Cancer Cell Int 20, 238 (2020). https://doi.org/10.1186/s12935-020-01315-7
- PTPN family members
- Digestive tract cancers
- Clinical features