Abstract
Background:
Lung adenocarcinoma (LUAD) is the most common diagnostic histologic subtype of non-small cell lung cancer, but the role of receptor-type tyrosine-protein phosphatase-like N (PTPRN) in LUAD has not been studied.
Methods:
We conducted a bioinformatic analysis to identify the expression of PTPRN on LUAD data from the Cancer Genome Atlas (TCGA) and the relationship between PTPRN and overall survival of LUAD patients. The effects of PTPRN on the migration ability of LUAD cells and the underlying mechanisms were investigated by in vitro and in vivo assays (i.e., wound healing assay, transwell assay, western blotting, xenograft model, and immunohistochemistry). Gene-set enrichment analysis and computational resource were used to analyze the correlation between PTPRN and different tumor-infiltrating immune cells (TIICs). Lactate dehydrogenase assay and Enzyme-linked immunosorbent assay were conducted to examine natural killer (NK) cell cytotoxicity.
Results:
In our study, we found that PTPRN was up-regulated in LUAD and related to metastasis of LUAD patients. Besides, PTPRN was correlated with poor prognosis in the TCGA-LUAD dataset. PTPRN overexpression promoted LUAD cell migration and the expression of EMT markers by influencing MEK/ERK and PI3K/AKT signaling. Moreover, PTPRN expression was significantly associated with TIICs, especially NK cells. A549 and H1299 cells overexpressed PTPRN inhibited NK cell cytotoxicity.
Conclusion:
Taken together, these findings demonstrated that PTPRN might be a potential and novel therapeutic target modulating antitumor immune response in treatment of LUAD.
Introduction
Lung cancer, the main cause of cancer incidence and mortality, accounted for nearly a fifth of cancer deaths in 2018 (). Besides, lung cancer is predicted to cause around a quarter of all cancer deaths in America in 2020 (). More than 80% of lung cancer cases are non-small cell lung cancer (NSCLC; ). Lung adenocarcinoma (LUAD), the most common histologic type of NSCLC, is responsible for almost 40% of deaths related to lung cancer (; ). Although there are many remarkable improvements in the treatment of LUAD, the overall survival (OS) of LUAD remains unsatisfactory. The 5-year survival rate of LUAD patients is no more than 20% and has not increased in the past few years (; ). Many patients die of its metastatic relapse after surgery (). Therefore, it is vital to discover new prognostic biomarkers and therapeutic strategies for patients who succumb to metastatic characteristics of LUAD.
Carcinogenesis is proved to be a long-term and sophisticated process that metastasis is a cardinal feature of OS (). Epithelial–mesenchymal transition (EMT) plays a key role in embryonic development, fibrosis, and cancer metastasis (). EMT is regarded as the initial event related to metastasis of cancers. However, the mechanism of EMT is not fully elucidated. Thus, focusing on regulated molecules related to EMT might contribute to identify potential biomarkers of prognosis in cancers.
Receptor-type tyrosine-protein phosphatase-like N (PTPRN), also called islet cell autoantigen 512 (ICA512/IA-2), is a tyrosine phosphatase-like intrinsic membrane protein which participates in the biological process of insulin secretory granules in pancreatic islet β-cells (). It was reported to have a high expression in neuroendocrine cells (; ) and is related to diabetes (). Several studies reported that PTPRN plays a potential role in solid tumors (; ; ; Xu et al., 2018; Zhangyuan et al., 2018; Yin et al., 2019). PTPRN was overexpressed in a cluster of glioblastoma patients and indicated poor survival (). Furthermore, Wen et al. and Xu et al. (2018) showed that PTPRN was a biomarker for predicting the prognosis of glioblastoma and could estimate the OS of glioblastoma patients (Yin et al., 2019). PTPRN was also an independent prognostic factor of hepatocellular carcinoma (Zhangyuan et al., 2018). Abnormal DNA hypermethylation of PTPRN was related to OS of patients with ovarian cancer (). PTPRN mRNA and protein expression levels were increased in breast cancer cells under hypoxia (). Although PTPRN could be a new marker for detecting neuroendocrine differentiation in lung cancer, the potential relationship between PTPRN and on EMT is still unclear ().
The tumor microenvironment (TME) promotes immune escape which is regarded as one of cancer hallmarks (). The relationship between immune cells in TME and cancer cells is continuous and participates from tumorigenesis initiation to tumor metastasis (). A study indicated that EMT was closely related to TME (). NK cells are recruited to TME (). NK cells play an important role in early tumors and metastasis elimination (). The relationship between PTPRN, TME, and NK cells in LUAD needs to be further elucidated.
In the present study, we demonstrated that the expression of PTPRN was highly up-regulated in LUAD and associated with metastasis and poor OS of LUAD patients. PTPRN overexpression promoted LUAD cell metastasis capacity by influencing the MEK/ERK and PI3K/AKT pathways. Our study provided functional links between PTPRN and LUAD metastasis. Furthermore, our study reported that PTPRN influenced TME compositions and regulated NK cell regulation.
Materials and Methods
Publicly Available Data Analysis
Gene expression profiles and clinical data of LUAD patients were downloaded from the Cancer Genome Atlas (TCGA) database1. The samples of LUAD were divided into two groups according to the cutoff of the median PTPRN expression. The differences in N-cadherin, Slug, Snail, Twist, and Vimentin between the PTPRN low- and high-expression groups were analyzed. The Kaplan–Meier Plotter was used to analyze the correlation between PTPRN expression and survival time of lung cancer patients2.
The ESTIMATE algorithm by feat of R-package “estimate” was used to calculate immune scores that infer the proportion of stromal and immune cells of tumor samples (Yoshihara et al., 2013). It also exhibited three kinds of scores called ImmuneScore, StromalScore, and ESTIMATEScore, which were related to immune, stromal cells and the sum of them in TME.
Evaluation of Tumor-Infiltrating Immune Cells
Computational resource (CIBERSORTx)3, an analytical tool based on gene expression profiles, predicts the ratios of member cell types in mixed cells (). Thus, we can use this process to estimate tumor-infiltrating immune cells (TIICs). Immune response of 22 TIICs was measured to assess the association with PTPRN expression in LUAD. We established a dataset of gene expression by using standard annotation files and a defaulted signature matrix at 1000 permutations. CIBERSORTx predicted a P-value for deconvolution through Monte Carlo sampling to establish the confidence level in the outcomes. In order to analyze the impact of PTPRN in TME of LUAD, 534 tumor samples were utilized in the TCGA-LUAD dataset. We set P-value < 0.05 as the criterion to define the lymphocyte types affected by PTPRN.
Furthermore, the “correlation” module of GEPIA4 was used to determine the association between PTPRN expression and possible gene markers of TIICs. The gene markers included markers of CD8+ T cells, T cells (general), B cells, monocytes, tumor-associated macrophage (TAM), T-helper 1 (Th1) cells, T-helper 2 (Th2) cells, follicular helper T (Tfh) cells, T-helper 17 (Th17) cells, regulatory T cells (Treg), exhausted T cells, neutrophils, NK cells, and Mast cells. The correlation heat map was conducted to detect the relationship between every two different immune cells in 22 types of immune cells in TME.
Gene-Set Enrichment Analysis
The TCGA-LUAD dataset with the gene sets (hallmark and C5 gene set v7.1 collections) was analyzed by Gene-set enrichment analysis (GSEA) in order to obtain biological processes enriched by PTPRN. The samples were divided into two groups: the top 50% of them were defined as the high PTPRN expression group, and the bottom 50% were defined as the low PTPRN expression group.
NK Cell Isolation
Primary NK cells were purified from fresh peripheral blood samples from healthy volunteers. Peripheral blood mononuclear cells were isolated by RosetteSepTM Human NK Cell Enrichment Cocktail (STEMCELL Technologies).
Flow Cytometry Analysis
We stained NK cells with 5 μL/test FITC-CD3 (Biolegend) and 5 μL/test PE-CD56 (Biolegend) at 4°C for 30 min. NK cells were washed with PBS and suspended in 400 μL of PBS after staining. Flow cytometry analysis was performed on a BD FACSVia Flow Cytometer (BD Biosciences, United States).
Cell Culture and Transfection
The H1299 and A549 lung cancer cell lines were obtained from the Institute of Biochemistry and Cell Biology at the Chinese Academy of Sciences (Shanghai, China). H1299 was maintained in the RPMI-1640 medium (HyClone) with 10% fetal bovine serum (FBS, PAA, Germany), 100 mg/ml penicillin (Invitrogen), and 100 U/mL streptomycin (Invitrogen). A549 was maintained in F-12K medium (Gibco) with 10% fetal bovine serum (FBS, PAA, Germany), 100 mg/ml penicillin (Invitrogen), and 100 U/mL streptomycin (Invitrogen). Purified NK cells were cultured in NK cell culture medium (Cell Science and Technology Institute Inc.). NK cells and H1299 and A549 cell lines were cultured at 37°C in a 5% CO2 and 95% air incubator. PTPRN-pcDNA3.1 and pcDNA3.1 vectors were purchased from Sangon Co., Ltd. (Shanghai, China). The empty vector was obtained from Sangon. The cells were transfected with plasmids using Lipofectamine 3000 (Invitrogen) according to the manufacturer’s protocol. To obtain LLC cells which overexpressed PTPRN stably, LLC cells were electroporated with the BIO-RAD GenePulser Xcell system (BIO-RAD, United States). The condition of different transfection parameters (pulse settings) was tested in order to optimize the parameters (). The electroporation pulse voltage and width were adjusted. We found that the optimal condition of voltage pulse was 248 V and the optimal condition of pulse width was 30 ms.
Wound-Healing Assay
For detecting cell migration, wound assay was conducted; A549 and H1299 cells with transection of each group were seeded into six-well plates. When cells reached around 90% confluency, a pipette tip was used vertically to make scratches on the cell layer. Cells were washed with PBS to clear scratched cell debris and cultured with a serum-free medium in the incubator. Cells were observed and photographed using a microscope (Nikon, Japan) after cultured for 0 h and 24 h.
Transwell Migration Assay
Cell migration was assessed in a 24-well transwell chamber (3422, Corning, United States) according to the manufacturer’s protocol. H1299 and A549 cells were seeded at a density of 4 × 104 cells/well in the upper chambers with a serum-free RPMI-1640 medium. The lower chambers contained the RPMI-1640 medium with 20% FBS. After 24 h, cells which migrated out of the inserts were immobilized by methanol and stained with 0.1% crystal violet (Solarbio, China). Cells were observed using the microscope and photographed (magnification, ×100). Cells were counted by five fields per well at random.
Western Blotting
Cells were lysed by radioimmunoprecipitation assay lysate buffer (Beyotime, Jiangsu, China) on ice for 30 min. The supernatants of cells were collected. Bicinchoninic acid protein quantitation (Takara, Japan) was used to determine protein concentrations. Proteins from cell lysate were transferred to PVDF membranes (Millipore, Bedford, MA, United States) after SDS-PAGE electrophoresis. The primary and secondary antibodies are shown in Table 1.
TABLE 1
| Antibody | Company/provider | Dilution ratio |
| Anti-human PTPRN | Proteintech | 1:1000 |
| Anti-human Vimentin | Cell signaling technology | 1:1000 |
| Anti-human E-cadherin | Cell signaling technology | 1:1000 |
| Anti-human N-cadherin | Cell signaling technology | 1:1000 |
| Anti-human β-actin | BOSTER | 1:1000 |
| Anti-human ERK1/2 | Cell signaling technology | 1:1000 |
| Anti-human p-ERK1/2 | Cell signaling technology | 1:1000 |
| Anti-human AKT | Proteintech | 1:1000 |
| Anti-human p-ERK1/2 | Proteintech | 1:1000 |
| Anti-human Cyclin D1 | Cell signaling technology | 1:1000 |
| Goat anti-mouse IgG | Absin bioscience Inc | 1:10000 |
| Goat anti-rabbit IgG | Absin bioscience Inc | 1:10000 |
| Anti-human PTPRN (for IHC) | Proteintech | 1:100 |
| Anti-human Vimentin (for IHC) | Cell signaling technology | 1:100 |
| Anti-mouse KLRA1 (for IHC) | Bioss | 1:100 |
Antibodies used for WB and IHC.
IHC and Quantification Analysis
Mouse tumor sections and LUAD tissue samples were dewaxed by xylene and rehydrated by graded alcohol series. We used citrate buffer (pH = 6.0) to retrieve the antigen under high pressure. We used Ultra-sensitive S-P kit (Maixin-Bio, China) to block the endogenous peroxidase activity of sections. The sections were incubated with primary antibodies (shown in Table 1) at 4°C overnight. Then, sections were incubated with a secondary antibody. DAB reagent was used to visualize the sections.
Every slide was examined five times, and cells were observed at 400 × magnification. The expression of PTPRN was estimated by staining intensity and percentage of positive cells. Staining intensity was defined as 0 (negative), 1 (weak), 2 (moderate), and 3 (strong). The percentage of positive cells was defined as 0 (<5%), 1 (6–25%), 2 (26–50%), 3 (51–75%), and 4 (>76%). The final staining score of Immunohistochemistry (IHC) was calculated by multiplying the staining intensity levels with the positive percentage staining scores (Zhao et al., 2018).
Lactate Dehydrogenase Cytotoxicity Assay
The cytotoxicity of NK cells was conducted by lactate dehydrogenase (LDH) cytotoxicity assay (Beyotime, China), according to the manufacturer’s instructions. Briefly, we collected cell supernatants which cocultured NK cells and A549 and H1299 cells for 4 h. Then, we performed LDH assay. The formula of cytotoxicity was calculated: Cytotoxicity (%) = [(ODExperimental value–ODSpontaneous value)/(ODMaximum value–ODSpontaneous value)] × 100 ().
Enzyme-Linked Immunosorbent Assay
IFN-γ and TNF-α in cell supernatants were detected using human IFN-γ precoated Enzyme-linked immunosorbent assay (ELISA) kit and human TNF-α precoated ELISA kit (DAKEWE, China) according to the manufacturer’s instructions. A549 and H1299 cells were seeded in 96-well plates at a density of 1.0 × 104 cells per well. NK cells were added into the 96-well plates at a density of 1.0 × 104 cells per well. After 4 h of incubation, cell supernatants were collected to perform ELISA assay. Absorbance was read at Infinite 200 PRO microplate reader (TECAN, Switzerland).
Xenograft Model
To study the influence of PTPRN in LUAD, LLC cells (1 × 106) were suspended in 100 μL of PBS and injected into mammary fat pads of 4- to 5-week-old BALB/c(nu/nu) mice (Hua Fukang Biological Technologies Inc, Beijing, China). Mice were randomized into the following four groups (n = 6 per group): NC-cDNA, PTPRN-cDNA, NC-cDNA with TGF-β1, and PTPRN-cDNA with TGF-β1. TGF-β1 (Novoprotein; 4 μg/kg bodyweight) was injected intraperitoneally (i.p.) every 5th day after cell inoculation as previously described to establish TGF-β-induced cancer cell invasion (). The tumor diameter and weights of mice were measured every other day. Tumors were removed from xenograft models and weighed until tumors appeared to collapse. Tumor volume (mm3) was measured using a digital caliper and calculated according to the following equation: (width)2 × (length/2). All mice were bred at pathogen-free conditions in the Animal Centre of China Medical University. All animal studies were performed according to the National Institute of Health Guide for the Care and Use of Laboratory Animals.
Statistical Analysis
Quantitative data were expressed as the means ± SD of at least three independent experiments. GraphPad Prism 7.0 (GraphPad Software, La Jolla, CA, United States) was used to evaluate all experimental values. Student’s independent t test was used to perform statistical analysis between two experimental groups, while analysis of variance (ANOVA) was used to perform analyses among three experimental groups. Pearson’s correlation coefficient was used to analyze the correlation between PTPRN and ImmuneScore, StromalScore, and ESTIMATEScore. P value < 0.05 was considered statistically significant in all cases.
Results
Detection of PTPRN Protein Expression Is Analyzed by IHC
Sectioning experiments of 40 LUAD clinical samples have been done to elucidate the relationship between PTPRN and OS. Representative results for PTPRN protein in LUAD tissue are shown in Figure 1A. According to the ROC curve (AUC = 0.727, P = 0.0019), IHC scores of PTPRN were classified as low expression group (≤2) or high expression group (>2; data not shown). For survival analysis, the LUAD patients were divided into two groups with low or high PTPRN expression by IHC scores. The Kaplan–Meier survival curve demonstrated that LUAD with PTPRN high expression had an unfavorable OS (P = 0.0257; Figure 1B).
FIGURE 1
PTPRN Influences LUAD Cell Migration in vitro
To identify the effect of PTPRN on migration of LUAD cells, PTPRN was overexpressed in A549 and H1299 with low expression of PTPRN. The overexpression efficiency of PTPRN was evaluated by western blot. After transfection with PTPRN, PTPRN expression was upregulated in comparison with the negative control group (Figure 2A). The migration ability of the control group and PTPRN-overexpressing cells was examined. The wound-healing assay significantly revealed that the cell-covered area was increased in the PTPRN-overexpressed group compared with the negative control group in Figure 2B. Furthermore, we found that the number of migrated cells was increased in the PTPRN-overexpressed group compared with the control group in transwell migration assay (Figure 2C).
FIGURE 2
PTPRN Is Related to Metastasis of LUAD Patients
To elucidate whether PTPRN affected metastasis of lung cancer, we evaluated the expression level of PTPRN in TCGA datasets. The TCGA RNA Seq data demonstrated that PTPRN was significantly overexpressed in lung cancer compared with noncancerous tissue samples (Figure 3A). Besides, PTPRN was highly expressed in NSCLC patients with distant metastasis and lymph-node metastasis compared with non-metastatic NSCLC patients (Figure 3B). PTPRN expression in LUAD was distinctively higher than that of normal tissues (Figure 3C). PTPRN was also highly expressed in metastatic LUAD patients compared with non-metastatic LUAD patients (Figure 3D). Meanwhile, PTPRN was also significantly expressed in LUSC compared with normal tissues (Figure 3E). There were no differences in PTPRN expression between metastatic LUSC patients and non-metastatic patients (Figure 3F). Then, we performed bioinformatic analysis of data mining with the Kaplan–Meier plotter. We found that PTPRN high expression was relevant to a shorter OS in NSCLC patients and LUAD patients (Figures 3G,H). However, OS in LUSC patients showed no correlation between high PTPRN expression and low PTPRN expression (Figure 3I). These results demonstrated that PTPRN was related to metastasis of LUAD and affected OS of LUAD patients, not in LUSC patients.
FIGURE 3
High PTPRN Expression Was an Independent Prognostic Factor for Poor OS in LUAD Patients
In order to elucidate whether PTPRN was involved in LUAD development, we analyzed the relationship between the expression of PTPRN and clinical pathology factors. Univariate analysis (Table 2) revealed that the advanced pathologic stage (III–IV), tumor stage (T3–T4), regional lymph node (N1–N3), metastasis (M1), and high PTPRN expression were relevant to remarkably shorter OS in LUAD patients (Table 2). Then, the multivariate analysis verified that the high expression of PTPRN was an independent prognostic factor for poor OS (Table 2).
TABLE 2
| Univariate analysis | Multivariate analysis | |||||
| HR | p value | 95% CI | HR | p value | 95% CI | |
| Age | ||||||
| >66 versus ≤66 | 1.218 | 0.268 | 0.859–1.726 | |||
| Gender | ||||||
| Male versus female | 1.076 | 0.679 | 0.760–1.523 | |||
| Smoking history | ||||||
| 2/3/4/5 versus 1 | 1.075 | 0.772 | 0.659–1.753 | |||
| Pathologic stage | ||||||
| III–IV versus I–II | 2.503 | <0.001 | 1.741–3.599 | 1.009 | 0.979 | 0.514–1.980 |
| Tumor stage | ||||||
| T3–T4 versus T1–T2 | 2.586 | <0.001 | 1.682–3.977 | 2.194 | 0.002 | 1.344–3.583 |
| Regional lymph node | ||||||
| N1–N3 versus N0 | 2.438 | <0.001 | 1.719–3.457 | 2.104 | 0.015 | 1.158–3.823 |
| Metastasis | ||||||
| M1 versus M0 | 1.878 | 0.032 | 1.056–3.341 | 1.502 | 0.188 | 0.820–2.750 |
| PTPRN | ||||||
| High versus Low | 1.710 | 0.003 | 1.200–2.436 | 1.470 | 0.036 | 1.026–2.107 |
Univariate and multivariate analyses of OS in patients with LUAD in TCGA.
PTPRN Is Associated With EMT of LUAD Through MEK/ERK and PI3K/AKT Signaling Pathways
To illustrate how PTPRN influenced the metastatic phenotype of LUAD, LUAD samples of TCGA were divided into PTPRN high and low expression groups by the median of PTPRN expression. GSEA was conducted to choose significant enriched pathways according to the normalized enrichment score (NES). Then, we found that PTPRN was related to EMT (Figure 4A and Table 3). EMT is regarded as a vital step in metastasis progression and promotes cancer cell migration and invasion. Then, we analyzed the relationship between PTPRN and EMT markers. The expressions of EMT markers including N-cadherin, Slug, Snail, Twist, and Vimentin were significantly higher in PTPRN-high patients than PTPRN-low patients (Figures 4B–F).
FIGURE 4
TABLE 3
| Description | Gene markers | LUAD | |||
| Tumor | Normal | ||||
| R | P | R | P | ||
| CD8+ T cell | CD8A | -0.018 | 0.7 | -0.27 | 0.04 |
| CD8B | 0.0036 | 0.94 | -0.22 | 0.098 | |
| T cell (general) | CD3D | -0.055 | 0.22 | -0.12 | 0.35 |
| CD3E | -0.091 | 0.046 | -0.12 | 0.38 | |
| CD2 | -0.011 | 0.013 | 0.018 | 0.89 | |
| B cell | CD19 | 0.052 | 0.26 | 0.053 | 0.69 |
| CD27 | 0.027 | 0.56 | -0.15 | 0.26 | |
| Monocyte | CD86 | 0.085 | 0.061 | -0.19 | 0.16 |
| CD115 | 0.062 | 0.18 | 0.051 | 0.70 | |
| TAM | CCL2 | 0.15 | 0.00076 | -0.065 | 0.63 |
| CD68 | 0.16 | 0.00046 | -0.19 | 0.16 | |
| IL-10 | 0.079 | 0.085 | -0.11 | 0.39 | |
| Th1 | IFN-γ (IFNG) | 0.059 | 0.19 | 0.17 | 0.20 |
| STAT1 | 0.18 | 6.3e-05 | 0.17 | 0.21 | |
| T-bet (TBX21) | -0.058 | 0.21 | 0.14 | 0.29 | |
| TNF-α (TNF) | 0.022 | 0.62 | 0.15 | 0.25 | |
| Th2 | STAT6 | -0.24 | 6e-08 | 0.062 | 0.64 |
| GATA3 | 0.0014 | 0.97 | -0.01 | 0.94 | |
| IL13 | -0.0085 | 0.85 | -0.038 | 0.77 | |
| STAT5A | -0.051 | 0.27 | 0.054 | 0.69 | |
| Tfh | BCL6 | -0.057 | 0.21 | 0.044 | 0.74 |
| Th17 | STAT3 | 0.0018 | 0.97 | 0.078 | 0.56 |
| IL17A | -0.028 | 0.54 | 0.014 | 0.92 | |
| Treg | CCR8 | 0.015 | 0.75 | 0.16 | 0.22 |
| STAT5B | 0.013 | 0.78 | 0.14 | 0.27 | |
| FOXP3 | 0.043 | 0.35 | 0.15 | 0.27 | |
| T cell exhaustion | LAG3 | 0.11 | 0.014 | 0.19 | 0.16 |
| CTLA4 | 0.021 | 0.64 | 0.044 | 0.74 | |
| PD-1 (PDCD1) | 0.078 | 0.085 | -0.099 | 0.46 | |
| TIM-3 | 0.059 | 0.19 | -0.21 | 0.10 | |
| Neutrophils | CD11b | 0.027 | 0.55 | -0.15 | 0.26 |
| CCR7 | -0.13 | 0.003 | 0.013 | 0.92 | |
| Natural killer cell | KIR2DL1 | 0.088 | 0.052 | 0.037 | 0.78 |
| KIR2DL3 | 0.1 | 0.022 | 0.039 | 0.77 | |
| KIR2DL4 | 0.24 | 7e-08 | -0.11 | 0.42 | |
| KIR3DL1 | 0.05 | 0.27 | -0.0078 | 0.95 | |
| KIR3DL2 | 0.077 | 0.09 | 0.0047 | 0.97 | |
| KIR3DL3 | 0.13 | 0.0036 | -0.084 | 0.53 | |
| CD56 | 0.14 | 0.0026 | 0.15 | 0.24 | |
| Mast cell | TPSB2 | -0.094 | 0.039 | -0.05 | 0.7 |
| TPSAB1 | -0.08 | 0.079 | 0.033 | 0.8 | |
| CPA3 | -0.092 | 0.044 | -0.037 | 0.78 | |
| MS4A2 | -0.15 | 0.0082 | -0.0059 | 0.96 | |
| HDC | -0.099 | 0.029 | -0.088 | 0.51 | |
A correlation analysis using the “Correlation” module of GEPIA.
Tumor, correlation analysis in LUAD tumor tissue of TCGA. Normal, correlation analysis in LUAD normal tissue of TCGA.
Furthermore, we also found that PTPRN was related to the PI3K/AKT/mTOR signaling pathway by GSEA (Figure 4G). Western blot assay showed that Vimentin and N-cadherin expression levels were increased in PTPRN-overexpressing A549 and H1299 cells (Figure 4H). However, E-cadherin expression level was decreased in PTPRN-overexpressing A549 and H1299 cells (Figure 4H). The protein expression of the phosphorylated ERK1/2, AKT, and cell cycle-related protein Cyclin D1 was increased by overexpressing PTPRN in A549 and H1299 cells (Figure 4I). However, PTPRN barely affected the total protein expression of ERK1/2 and AKT (Figure 4I). These results above added increasing evidence that PTPRN was related to EMT characteristics and regulated MEK/ERK and PI3K/AKT signaling pathways.
PTPRN Is Related to Immune Infiltrates
Gene-set enrichment analysis also found that high PTPRN expression was enriched in inflammatory response, TNFα signaling via NF-κB and TGFβ signaling pathways (Figures 5A–C and Table 3), indicating that PTPRN might be associated with immune response in LUAD to promote EMT. Furthermore, we demonstrated several immune-related pathways by GSEA (Figures 5D–F and Table 4), including negative regulation of immune response, regulation of immune effector process, and regulation production involved in immune response. These results showed that PTPRN was involved in immune response. Higher scores calculated in ImmuneScore, StromalScore, and ESTIMATEScore represented a higher proportion of immune and stromal components in TME. The previous study showed that these increased scores of LUAD inferred favorable survival of LUAD patients (). PTPRN was negatively correlated with ImmuneScore (P = 0.0108, Figure 5G). However, PTPRN expression had no correlation with StromalScore and ESTIMATEScore (P = 0.4799 and P = 0.0690, respectively, Figures 5H,I). These results reminded us that PTPRN might influence immune cell infiltration.
FIGURE 5
TABLE 4
| MsigDB collection | Gene set name | NES | NOM p-val | FDR q-val |
| h.all.v7.1.symbols.gmt | HALLMARK_EPITHELIAL_MESENCHYMAL_TRANSITION | −2.30 | 0.000 | 0.000 |
| HALLMARK_PI3K_AKT_mTOR | −1.89 | 0.000 | 0.000 | |
| HALLMARK_TNFA_SIGNALING_VIA_NFKB | −1.93 | 0.004 | 0.004 | |
| HALLMARK_INFLAMMATORY_RESPONSE | −1.89 | 0.002 | 0.004 | |
| HALLMARK_TGF_BETA_SIGNALING | −1.63 | 0.031 | 0.038 | |
| c5.all.v7.1.symbols.gmt | GO_NEGATIVE_REGULATION_OF_IMMUNE_RESPONSE | −1.93 | 0.000 | 0.000 |
| GO_REGULATION_OF_IMMUNE_EFFECTOR_PROCESS | −2.08 | 0.000 | 0.000 | |
| GO_REGULATION_OF_CYTOKINE_PRODUCTION_INVOLVED_IN_IMMUNE_RESPONSE | −2.02 | 0.000 | 0.000 |
Gene set enriched in LUAD samples with high-expression PTPRN.
NES: normalized enrichment score; NOM: nominal; FDR: false discovery rate.
The Expression of PTPRN Is Correlated With Proportion of Immune Infiltrating Cells
Studies have shown that tumor-infiltrating lymphocytes are an independent factor predicting survival in cancer patients (; ). Different cancers had various immune infiltrations, which affects cancer cell proliferation and metastasis. Therefore, we further investigated the correlation between PTPRN expression and immune infiltrating cells in the microenvironment of LUAD patients. The data of 534 LUAD tumor tissue samples were downloaded from TCGA and divided into high and low expression by the median expression of PTPRN. We used an established CIBERSORTx to explore the gene expression profiles of LUAD samples to infer the different fractions of 22 types of immune cells between high expression group and low expression group. Figure 6 shows the results of 22 immune cell subtypes. As shown in Figure 6A, immune cells such as naive B cells, T cells CD8, resting memory T cells CD4, activated memory T cells CD4, resting NK cells, monocytes, M0 macrophages, resting mast cells, activated mast cells, and neutrophils were affected by PTPRN expression. The fractions of activated memory CD4 T cells (P = 0.002), resting NK cells (P = 0.024), macrophages M0 (P < 0.001), activated mast cells (P = 0.004), and neutrophils (P < 0.001) were up-regulated in the high-PTPRN expression group. Meanwhile, the proportions of naive B cells (P = 0.004), CD8 T cells (P = 0.05), resting memory T cells (P < 0.001), monocytes (P = 0.001), and resting mast cells (P < 0.001) were down-regulated in high-expression PTPRN. Then, relationships between 22 immune cell subsets were analyzed. The correlation heat map (Figure 6B) reflected that there were weak and moderate correlations between the proportions of different TIICs subpopulations.
FIGURE 6
The correlation between PTPRN expression and different immune cell type markers was further investigated in LUAD through the “Correlation” module of GEPIA. We used Spearman correlation to evaluate the correlation coefficient. Our study reminded us that gene markers affected by PTPRN expression included CD3E and CD2 of T cell (general); CCL2 and CD68 of TAM; STAT1 of Th1; STAT6 of Th2; LAG3 of T cell exhaustion; CCR7 of neutrophils; KIR2DL3, KIR2DL4, KIR3DL3, and CD56 of NK cell; and TPSB2, CPA3, MS4A2, and HDC of mast cell (Table 4). Results between PTPRN expression and markers of T cells, macrophages, neutrophils, mast cells, and NK cells were similar to the results of CIBERSORTx. According to the results of CIBERSORTx and GEPIA, these findings indicated that PTPRN might play a key role in regulating T cells, macrophage, neutrophils, mast cells, and NK cells of the tumor-infiltrating immune environment. According to GEPIA analysis, PTPRN was associated with NK cell inhibitory markers, indicating that PTPRN might influence NK cell activation.
PTPRN Is Related to Inhibitory Markers of NK Cells and Inhibits NK Cell Cytotoxicity
To identify how PTPRN influences NK cells, we analyzed the relationship between PTPRN and NK cell inhibitory receptors in the TCGA-LUAD dataset. We found that NK cell inhibitory markers such as LIRB1 and KLRC1 were expressed higher in the PTPRN high-expression group compared with the PTPRN low-expression group (Figures 7A,B). Peripheral blood NK cells were isolated from healthy donors. In peripheral blood, human NK cells accounted for approximately 10–15% of cells. NK cells were identified according to the CD3–CD56+ profile. NK cells were determined by flow cytometry on Day 14. We obtained 90.97 ± 2.55% CD3–CD56+ NK cells after purification of NK cells (Figure 7C). NK cell cytotoxic activity was measured by LDH assay. As shown in Figure 7D, stimulated NK cells decreased cytotoxicity when A549 and H1299 overexpressed PTPRN. Furthermore, we found that IFN-γ and TNF-α levels in cell culture supernatant were reduced when A549 and H1299 were overexpressed PTPRN (Figures 7E,F). These results suggested that PTPRN overexpression in LUAD cells might inhibit NK cell cytotoxicity.
FIGURE 7
PTPRN Influences EMT and NK Cell Inhibitory Marker in vivo
To identify the effect of PTPRN on EMT and NK cell inhibition, LLC cells transfected with NC-cDNA or PTPRN-cDNA were subcutaneously injected into nude mice. Then, TGF-β1 was injected i.p. post cell inoculation. A study showed that TGF-β1 is one of the strongest inducers of EMT (). Tumor size was significantly increased in the PTPRN overexpression group compared with negative control. Besides, TGF-β1 had the same effect on increasing tumor size (Figures 8A,B). Consistent with in vitro experiments, PTPRN promoted the expression of vimentin and KLRA1 (NK cells inhibitory marker of mice) detected by IHC (). TGF-β1 had the same effect on increasing vimentin expression compared with PTPRN (Figure 8C). Furthermore, TGF-β1 and PTPRN had a synergistic effect on the increase in vimentin expression (Figure 8C). These results indicated that PTPRN promoted EMT and NK cell inactivation in vivo.
FIGURE 8
Discussion
Protein phosphatase-like N was a membrane protein associated with diabetes. Mziaut and colleagues found that the depletion of ICA512 reduced insulin release stimulated by glucose through downregulating F-actin modifier villin in pancreatic islet β-cells (). Besides, PTPRN was reported to be overexpressed in several cancers such as glioblastoma (Xu et al., 2018; Yin et al., 2019) and hepatocellular carcinoma (Zhangyuan et al., 2018). In our study, data mining by the TCGA dataset was performed to demonstrate that PTPRN is up-regulated in LUAD tissues and high PTPRN expression is related to poor prognosis in LUAD patients. Surprisingly, we also found that high PTPRN expression is associated with metastasis of LUAD patients. The molecular mechanism of PTPRN promoting LUAD metastasis is still unclear.
Cancer metastasis is a complex process where cancer cells migrate from the primary site to another tissue for colonization (; Zada et al., 2019), causing most cancer-related deaths (). Metastasis is divided into four steps: (1) detachment from the primary tumor; (2) blood migration; (3) stromal tissue invasions; and (4) proliferation in distant location (). EMT with characteristics of loss of epithelial properties and increase of mesenchymal properties is crucial to initiation of metastasis (; Zada et al., 2019).
Furthermore, GSEA found that PTPRN was related to EMT. The TCGA-LUAD dataset showed that PTPRN is associated with EMT markers such as N-cadherin, Slug. Then, we transfected PTPRN into A549 and H1299 cells. The cells with PTPRN high expression expressed more EMT markers such as N-cadherin at mRNA and protein levels. Wound healing assay and transwell cell assay were conducted, suggesting that PTPRN promotes LUAD cell invasion and migration by regulating the EMT pathway. Then, our results indicated that the expressions of Vimentin, N-cadherin, p-ERK, and p-AKT were up-regulated in the PTPRN overexpression group compared with the negative control group. Some studies have shown that EMT was related to the MEK/ERK and PI3K/AKT signaling pathways (; ; ).
The immune cells in TME are vital to the development of tumors and could be an indicator of prognostic and diagnostic assessment of cancer patients (). Then, we discussed whether PTPRN influences immune infiltration in TME. GSEA found that PTPRN influenced several immune-related pathways, suggesting that PTPRN might participate in EMT and regulated immune infiltration. The results of CIBERSORTx and GEPIA indicated that the expression of PTPRN is correlated with type markers of different immune cells, especially NK cells.
The antitumor function of innate immunity has gained great attention. NK cells, a subset of lymphoid cells triggering innate immunity, exert natural cytotoxicity against cancer cells and inhibit metastasis to different tissues (). The main reason for cancer-associated deaths is metastasis, and immune escape leads to invasion-metastasis cascade (). NK cells which target tumor cells rapidly are suppressed by molecules derived from tumors, stromal cells educated by tumors, and cancer cells, leading to initiation and multiple metastatic processes of tumors (). In our study, we found that the PTPRN high expression group expressed more NK cell inhibitory markers such as LIRB1 and KLRC1 in the TCGA-LUAD dataset. LDH assay indicated that A549 and H1299 cells with a PTPRN high expression level inhibit NK cell activation.
There are several strengths in our study. First, we found the impact of PTPRN in lung cancer, especially in LUAD. Second, we further identified the role of PTPRN in promoting migration and invasion of LUAD cells by regulating MEK/ERK and PI3K/AKT signaling pathways. Moreover, PTPRN was associated with immune cell infiltration in TME and inhibited NK cell activation. These results indicated that PTPRN could be a prognostic marker in LUAD patients and predict metastatic status of LUAD patients. PTPRN served as a potential antitumor immunotherapy strategy for LUAD patients. However, there are several limitations in our study. First, the way that PTPRN affected MEK/ERK and PI3K/AKT signaling pathways needed further study. Besides, biological mechanisms on how PTPRN influenced NK cell activation needed to be elucidated.
In conclusion, our study found that PTPRN can be a prognostic marker of LUAD patients and promotes the migration of LUAD cells in vitro and in vivo. Besides, PTPRN is correlated with immune infiltration of LUAD and inhibited NK cells activation. Our study identified the role of PTPRN in LUAD. Furthermore, PTPRN may play a key role in TME of LUAD by regulating immune cells infiltration especially NK cells, indicating that PTPRN as a potential and novel therapeutic target modulates antitumor immune response in treatment of LUAD.
Conclusion
In conclusion, the study demonstrated that PTPRN was up-regulated in LUAD and promoted metastasis of LUAD patients. Furthermore, high PTPRN expression indicated poor OS in LUAD patients. Importantly, PTPRN promoted cell migration by regulating MEK/ERK and PI3K/AKT signaling pathways and influenced TIICs, especially inhibiting NK cell activation. The study reminds us that PTPRN might be a prognostic marker of LUAD patients and therapeutic target regulating antitumor immune response in LUAD treatment.
Statements
Data availability statement
The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author/s.
Ethics statement
The studies involving human participants were reviewed and approved by the study was supported by Ethics Committee of China Medical University and they claimed that they approved this study. All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and national research committee and with the Helsinki declaration and its later amendments or comparable ethical standards. The plans for the use of the data were approved. The informed consent (written) is on file at China Medical University. The patients/participants provided their written informed consent to participate in this study. The animal study was reviewed and approved by all animal experiments were approved by the Institutional Animal Care and Use Committee of China Medical University.
Author contributions
XS, LZ, and MW conceived and designed the project. XS, XJ, MZ, and LJ designed and supervised experiments conducted in the laboratories. HY, LY, LC, and MJ performed experiments and/or data analyses. QW, LZ, and MW contributed reagents/analytic tools and/or grant support. XS and LW wrote the manuscript. XJ revised the manuscript. All authors discussed the results and commented on the manuscript.
Funding
This work was supported by grants from the National Natural Science Foundation of China (Nos. 82073281, 81573462, 81903658, 81703560, 81703756, and 82073884), Liaoning Revitalization Talents Program (No. XLYC1807201), Major Special S&T Projects in Liaoning Province (2019JH1/10300005), Shenyang S&T Projects (Nos. 19-109-4-09 and 20-204-4-22), Program for Shenyang High Level Talent Innovation and Entrepreneurship Team (2019-SYRCCY-B-01), and NSFC joint fund for regional innovation and development (No. U20A20413).
Conflict of interest
MW was employed by Shenyang Kangwei Medical Laboratory Analysis Co. LTD. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Abbreviations
- LUAD
Lung adenocarcinoma
- PTPRN
Receptor-type tyrosine-protein phosphatase-like N
- TCGA
The Cancer Genome Atlas
- EMT
Epithelial–mesenchymal transition
- NK cells
Natural killer cells
- NSCLC
Non-small cell lung cancer
- OS
Overall survival
- TME
Tumor microenvironment
- TIICs
Tumor-infiltrating immune cells
- TAM
Tumor-associated macrophage
- Th1 cells
T-helper 1 cells
- Th2 cells
T-helper 2 cells
- Tfh cells
Follicular helper T cells
- Th17 cells
T-helper 17 cells
- Treg cells
Regulatory T cells
- GSEA
Gene-set enrichment analysis
- PBMC
Peripheral blood mononuclear cells
- LDH
Lactate dehydrogenase
- ELISA
Enzyme-linked immunosorbent assay
- LLC
Lewis lung carcinoma
- IHC
Immunohistochemistry.
Footnotes
1.^https://portal.gdc.cancer.gov/
2.^http://kmplot.com/analysis/
References
1
AzimiF.ScolyerR. A.RumchevaP.MoncrieffM.MuraliR.McCarthyS. W.et al (2012). Tumor-infiltrating lymphocyte grade is an independent predictor of sentinel lymph node status and survival in patients with cutaneous melanoma.J. Clin. Oncol.302678–2683. 10.1200/JCO.2011.37.8539
2
BaldT.KrummelM. F.SmythM. J.BarryK. C. (2020). The NK cell-cancer cycle: advances and new challenges in NK cell-based immunotherapies.Nat. Immunol.21835–847. 10.1038/s41590-020-0728-z
3
BauerschlagD. O.AmmerpohlO.BrautigamK.SchemC.LinQ.WeigelM. T.et al (2011). Progression-free survival in ovarian cancer is reflected in epigenetic DNA methylation profiles.Oncology8012–20. 10.1159/000327746
4
BiK. W.WeiX. G.QinX. X.LiB. (2020). BTK has potential to be a prognostic factor for lung adenocarcinoma and an indicator for tumor microenvironment remodeling: a study based on TCGA data mining.Front. Oncol.10:424. 10.3389/fonc.2020.00424
5
BianJ.HeL.WuY.LiuW.MaH.SunM.et al (2020). Anterior gradient 2-derived peptide upregulates major histocompatibility complex class I-related chains A/B in hepatocellular carcinoma cells.Life Sci.246:117396. 10.1016/j.lfs.2020.117396
6
BrayF.FerlayJ.SoerjomataramI.SiegelR. L.TorreL. A.JemalA. (2018). Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries.CA Cancer J. Clin.68394–424. 10.3322/caac.21492
7
CaoZ.LiaoQ.SuM.HuangK.JinJ.CaoD. (2019). AKT and ERK dual inhibitors: the way forward?Cancer Lett.45930–40. 10.1016/j.canlet.2019.05.02
8
ChenY.WangD. D.WuY. P.SuD.ZhouT. Y.GaiR. H.et al (2017). MDM2 promotes epithelial-mesenchymal transition and metastasis of ovarian cancer SKOV3 cells.Br. J. Cancer1171192–1201. 10.1038/bjc.2017.265
9
ChenZ.FillmoreC. M.HammermanP. S.KimC. F.WongK. K. (2014). Non-small-cell lung cancers: a heterogeneous set of diseases.Nat. Rev. Cancer14535–546. 10.1038/nrc3775
10
EcclesS. A.WelchD. R. (2007). Metastasis: recent discoveries and novel treatment strategies.Lancet3691742–1757. 10.1016/S0140-6736(07)60781-8
11
FennellD. A.SummersY.CadranelJ.BenepalT.ChristophD. C.LalR.et al (2016). Cisplatin in the modern era: the backbone of first-line chemotherapy for non-small cell lung cancer.Cancer Treat. Rev.4442–50. 10.1016/j.ctrv.2016.01.003
12
Glatzel-PlucinskaN.PiotrowskaA.DziegielP.Podhorska-OkolowM. (2019). The role of SATB1 in tumour progression and metastasis.Int. J. Mol. Sci.20:4156. 10.3390/ijms20174156
13
GomiH.Kubota-MurataC.YasuiT.TsukiseA.ToriiS. (2013). Immunohistochemical analysis of IA-2 family of protein tyrosine phosphatases in rat gastrointestinal endocrine cells.J. Histochem. Cytochem.61156–168. 10.1369/0022155412466872
14
HanahanD.WeinbergR. A. (2011). Hallmarks of cancer: the next generation.Cell144646–674. 10.1016/j.cell.2011.02.013
15
IngegnereT.MariottiF. R.PelosiA.QuintarelliC.De AngelisB.TuminoN.et al (2019). Human CAR NK cells: a new non-viral method allowing high efficient transfection and strong tumor cell killing.Front. Immunol.10:957. 10.3389/fimmu.2019.00957
16
KasinskiA. L.SlackF. J. (2012). miRNA-34 prevents cancer initiation and progression in a therapeutically resistant K-ras and p53-induced mouse model of lung adenocarcinoma.Cancer Res.725576–5587. 10.1158/0008-5472.CAN-12-2001
17
LambertA. W.PattabiramanD. R.WeinbergR. A. (2017). Emerging biological principles of metastasis.Cell168670–691. 10.1016/j.cell.2016.11.037
18
LiC. M.ChenG.DaytonT. L.Kim-KiselakC.HoerschS.WhittakerC. A.et al (2013). Differential Tks5 isoform expression contributes to metastatic invasion of lung adenocarcinoma.Genes Dev.271557–1567. 10.1101/gad.222745.113
19
LiJ.LiX.ZhangC.ZhangC.WangH. (2020). A signature of tumor immune microenvironment genes associated with the prognosis of nonsmall cell lung cancer.Oncol. Rep.43795–806. 10.3892/or.2020.7464
20
LiangH.XiaoJ.ZhouZ.WuJ.GeF.LiZ.et al (2018). Hypoxia induces miR-153 through the IRE1alpha-XBP1 pathway to fine tune the HIF1alpha/VEGFA axis in breast cancer angiogenesis.Oncogene371961–1975. 10.1038/s41388-017-0089-8
21
LiuX.SunL.ZhangS.ZhangS.LiW. (2020). GINS2 facilitates epithelial-to-mesenchymal transition in non-small-cell lung cancer through modulating PI3K/Akt and MEK/ERK signaling.J. Cell. Physiol.2357747–7756. 10.1002/jcp.29381
22
MehlenP.PuisieuxA. (2006). Metastasis: a question of life or death.Nat. Rev. Cancer6449–458. 10.1038/nrc1886
23
MziautH.MulliganB.HobothP.OttoO.IvanovaA.HerbigM.et al (2016). The F-actin modifier villin regulates insulin granule dynamics and exocytosis downstream of islet cell autoantigen 512.Mol. Metab.5656–668. 10.1016/j.molmet.2016.05.015
24
NewmanA. M.SteenC. B.LiuC. L.GentlesA. J.ChaudhuriA. A.SchererF.et al (2019). Determining cell type abundance and expression from bulk tissues with digital cytometry.Nat. Biotechnol.37773–782. 10.1038/s41587-019-0114-2
25
NietoM. A.HuangR. Y.JacksonR. A.ThieryJ. P. (2016). EMT: 2016.Cell16621–45. 10.1016/j.cell.2016.06.028
26
OrrM. T.LanierL. L. (2011). Inhibitory Ly49 receptors on mouse natural killer cells.Curr. Top. Microbiol. Immunol.35067–87. 10.1007/82_2010_85
27
PastushenkoI.BlanpainC. (2019). EMT transition states during tumor progression and metastasis.Trends Cell Biol.29212–226. 10.1016/j.tcb.2018.12.001
28
Perez-RuizE.MeleroI.KopeckaJ.Sarmento-RibeiroA. B.Garcia-ArandaM.de Las RivasJ. (2020). Cancer immunotherapy resistance based on immune checkpoints inhibitors: targets, biomarkers, and remedies.Drug Resist. Updates53:100718. 10.1016/j.drup.2020.100718
29
SchaferA.HaenigB.ErupathilJ.StricknerP.SabatoD.WelfordR. W. D.et al (2021). Inhibition of endothelin-B receptor signaling synergizes with MAPK pathway inhibitors in BRAF mutated melanoma.Oncogene401659–1673. 10.1038/s41388-020-01628-x
30
ShergalisA.BankheadA.IIILuesakulU.MuangsinN.NeamatiN. (2018). Current challenges and opportunities in treating glioblastoma.Pharmacol. Rev.70412–445. 10.1124/pr.117.014944
31
SiegelR. L.MillerK. D.JemalA. (2020). Cancer statistics, 2020.CA Cancer J. Clin.707–30. 10.3322/caac.21590
32
SongY. H.ZhangC. Q.ChenF. F.LinX. Y. (2018). Upregulation of neural precursor cell expressed developmentally downregulated 4-1 is associated with poor prognosis and chemoresistance in lung adenocarcinoma.Chin. Med. J.13116–24. 10.4103/0366-6999.221262
33
SpiraA.EttingerD. S. (2004). Multidisciplinary management of lung cancer.N. Engl. J. Med.350379–392. 10.1056/NEJMra035536
34
StankovicB.BjorhovdeH. A. K.SkarshaugR.AamodtH.FrafjordA.MullerE.et al (2018). Immune cell composition in human non-small cell lung cancer.Front. Immunol.9:3101. 10.3389/fimmu.2018.03101
35
SteinbrennerH.NguyenT. B.WohlrabU.ScherbaumW. A.SeisslerJ. (2002). Effect of proinflammatory cytokines on gene expression of the diabetes-associated autoantigen IA-2 in INS-1 cells.Endocrinology1433839–3845. 10.1210/en.2002-220583
36
TakeyamaN.AnoY.WuG.KubotaN.SaekiK.SakudoA.et al (2009). Localization of insulinoma associated protein 2, IA-2 in mouse neuroendocrine tissues using two novel monoclonal antibodies.Life Sci.84678–687. 10.1016/j.lfs.2009.02.012
37
ToledoP. L.TorkkoJ. M.MullerA.WegbrodC.SonmezA.SolimenaM.et al (2019). ICA512 RESP18 homology domain is a protein-condensing factor and insulin fibrillation inhibitor.J. Biol. Chem.2948564–8576. 10.1074/jbc.RA119.007607
38
WangL.TongX.ZhouZ.WangS.LeiZ.ZhangT.et al (2018). Circular RNA hsa_circ_0008305 (circPTK2) inhibits TGF-beta-induced epithelial-mesenchymal transition and metastasis by controlling TIF1gamma in non-small cell lung cancer.Mol. Cancer17:140. 10.1186/s12943-018-0889-7
39
WangW.RenS.WangZ.ZhangC.HuangJ. (2020). Increased expression of TTC21A in lung adenocarcinoma infers favorable prognosis and high immune infiltrating level.Int. Immunopharmacol.78:106077. 10.1016/j.intimp.2019.106077
40
WuS. Y.FuT.JiangY. Z.ShaoZ. M. (2020). Natural killer cells in cancer biology and therapy.Mol. Cancer19:120. 10.1186/s12943-020-01238-x
41
XieH.NotkinsA. L.LanM. S. (1996). IA-2, a transmembrane protein tyrosine phosphatase, is expressed in human lung cancer cell lines with neuroendocrine phenotype.Cancer Res.562742–2744.
42
XuP.YangJ.LiuJ.YangX.LiaoJ.YuanF.et al (2018). Identification of glioblastoma gene prognosis modules based on weighted gene co-expression network analysis.BMC Med. Genomics11:96. 10.1186/s12920-018-0407-1
43
YinW.TangG.ZhouQ.CaoY.LiH.FuX.et al (2019). Expression profile analysis identifies a novel five-gene signature to improve prognosis prediction of glioblastoma.Front. Genet.10:419. 10.3389/fgene.2019.00419
44
YoshiharaK.ShahmoradgoliM.MartinezE.VegesnaR.KimH.Torres-GarciaW.et al (2013). Inferring tumour purity and stromal and immune cell admixture from expression data.Nat. Commun.4:2612. 10.1038/ncomms3612
45
ZadaS.HwangJ. S.AhmedM.LaiT. H.PhamT. M.KimD. R. (2019). Control of the epithelial-to-mesenchymal transition and cancer metastasis by autophagy-dependent SNAI1 degradation.Cells8:129. 10.3390/cells8020129
46
ZhangyuanG.YinY.ZhangW.YuW.JinK.WangF.et al (2018). Prognostic value of phosphotyrosine phosphatases in hepatocellular carcinoma.Cell. Physiol. Biochem.462335–2346. 10.1159/000489625
47
ZhaoX.ZhouL. L.LiX.NiJ.ChenP.MaR.et al (2018). Overexpression of KIF20A confers malignant phenotype of lung adenocarcinoma by promoting cell proliferation and inhibiting apoptosis.Cancer Med.74678–4689. 10.1002/cam4.1710
Summary
Keywords
PTPRN, LUAD, EMT, tumor-infiltrating immune cells, NK cells
Citation
Song X, Jiao X, Yan H, Yu L, Jiang L, Zhang M, Chen L, Ju M, Wang L, Wei Q, Zhao L and Wei M (2021) Overexpression of PTPRN Promotes Metastasis of Lung Adenocarcinoma and Suppresses NK Cell Cytotoxicity. Front. Cell Dev. Biol. 9:622018. doi: 10.3389/fcell.2021.622018
Received
27 October 2020
Accepted
06 April 2021
Published
02 June 2021
Volume
9 - 2021
Edited by
Zhengfei Zhu, Fudan University, China
Reviewed by
Avisek Majumder, University of California, San Francisco, United States; Pier Giorgio Petronini, University of Parma, Italy
Updates
Copyright
© 2021 Song, Jiao, Yan, Yu, Jiang, Zhang, Chen, Ju, Wang, Wei, Zhao and Wei.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Minjie Wei, weiminjiecmu@163.comLin Zhao, zl_cmu@163.com; lzhao@cmu.edu.cn
†These authors have contributed equally to this work
This article was submitted to Molecular and Cellular Oncology, a section of the journal Frontiers in Cell and Developmental Biology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.