The expression patterns and the diagnostic/prognostic roles of PTPN family members in digestive tract cancers

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.


Background
Tyrosyl phosphorylation is a dynamic and reversible process which plays a pivotal part in many cellular signaling pathways [1]. The dephosphorylation of tyrosine residues are catalyzed by a series of enzymes named protein tyrosine phosphatases (PTPs) [2]. Encoded by 103 genes, PTPs are sorted into four major superfamily classes [3], 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 [3]. 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 [7]. 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 [8]. PTPN22 plays an important part in regulating autophagy and NLR family pyrin domain containing 3 inflammasome activation [9]. It is hopeful that PTPN genes have potential to be served as prognostic and diagnostic indicators [10], 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 [11]. All of the three digestive tract cancers are ranked in the top 10 for incidence rates of tumors and have poor prognosis [11]. Although a couple of diagnostic biomarkers have been observed, robust biomarkers to predict clinical outcomes are still in urgent need [12]. 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.oncom ine.org/) which shows the expression information of genes in cancer and normal samples [13]. 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://ualca n.path.uab.edu) that offers online analysis of data from The Cancer Genome Atlas (TCGA) [14]. 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 ://cance rgeno me.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 [17]. Data of 184 ESCA and 11 matched normal samples, 415 STAD and 35 matched normal 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.
Keywords: PTPN family members, Digestive tract cancers, Expression, Diagnosis, Prognosis, Clinical features 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 X 2 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 .ncifc rf.gov/home.jsp; accessed on November 20, 2019) [18]. The gene-gene interaction network was structured using the Gene Multiple Association Network Integration Algorithm (GeneMANIA; https :// www.genem ania.org/; accessed on November 21, 2019) [19] and the Search Tool for the Retrieval of Interacting Genes Database (STRING v.10.0; https ://strin g-db. org/; accessed on November 23, 2019) was used to create a protein-protein interaction (PPI) network [20]. The Cytoscape software was applied to visualize network diagrams for PPI analysis [21].

Evaluation of PTPN expression in cancer cell lines by CCLE analysis
The Broad Institute Cancer Cell Line Encyclopedia (CCLE; http://www.broad insti tute.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 [22]. 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 nontumor 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  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).

Discussion
Being significant members of protein tyrosine phosphatases, PTPN regulates tyrosine phosphorylation and dephosphorylation in cellular signal transduction together with protein-tyrosine kinases [4]. 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  [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 [28]. While another study illustrated that elevated PTPN3 expression promotes tumor recurrence and is detrimental to the prognosis of intrahepatic cholangiocarcinoma patients [29]. These findings indicate that regulation mechanisms of PTPN3 in various cancers are complicated. Consistent with our results, Slattery et al. [30] reported that PTPN11 expression is upregulated in CRC. PTPN11 was suggested to impact the tumorigenesis and metastasis of CRC [31] and termed as a potential prognostic marker [32]. The data in one report showed that in patients with gastric B cell non-Hodgkin's lymphoma, PTPN21 is over expressed [33]. 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 [36]. 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 [37], 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 [38]. Mutated PTPN13 was suggested to be a tumor suppressor gene in colorectal cancer [28]. PTPN22 gene encodes an enzyme called lymphoid-specific tyrosine phosphatase, which functions as a master regulator in the biological process of relevant immune responses [39]. Upregulation of PTPN22 could result in impairment of regulatory T-cell differentiation in patients with non-ST-segment elevation acute coronary syndromes [40]. 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 [41]. PTPN5 could regulate tyrosine dephosphorylation needed for the activation of BAK, a noteworthy cell-death mediator in apoptosis [42]. It was once reported that PTPN14 regulated phosphorylation of p130Cas Y128 plays a crucial role in colorectal carcinogenesis [43]. And mutated PTPN14 is suggested to be a tumor suppressor gene for colorectal cancer, regulating cellular pathways that are appropriate for therapeutic intervention [28]. A recent study proved that PTPN14 suffices to inhibit migration and invasion of metastatic cancer cells [44]. In a xenograft breast cancer model, PTPN14 acts as a suppressor of metastasis of triple-negative breast cancer cells [45]. 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 [46]. 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 patientlevel, 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 [50]. 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 [51]. PTPN2 and PTPN22 were demonstrated negatively regulating T cell receptor signaling by dephosphorylating lymphocytespecific 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 [57]. 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.

Conclusions
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 welldesigned investigations are required to elucidate the significance of our findings and thus develop the clinical utility of PTPNs.
Additional file 1: Table S1. 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.