Comprehensive protein tyrosine phosphatase mRNA profiling identifies new regulators in the progression of glioma
© The Author(s). 2016
Received: 13 July 2016
Accepted: 19 August 2016
Published: 1 September 2016
The infiltrative behavior of diffuse gliomas severely reduces therapeutic potential of surgical resection and radiotherapy, and urges for the identification of new drug-targets affecting glioma growth and migration. To address the potential role of protein tyrosine phosphatases (PTPs), we performed mRNA expression profiling for 91 of the 109 known human PTP genes on a series of clinical diffuse glioma samples of different grades and compared our findings with in silico knowledge from REMBRANDT and TCGA databases. Overall PTP family expression levels appeared independent of characteristic genetic aberrations associated with lower grade or high grade gliomas. Notably, seven PTP genes (DUSP26, MTMR4, PTEN, PTPRM, PTPRN2, PTPRT and PTPRZ1) were differentially expressed between grade II-III gliomas and (grade IV) glioblastomas. For DUSP26, PTEN, PTPRM and PTPRT, lower expression levels correlated with poor prognosis, and overexpression of DUSP26 or PTPRT in E98 glioblastoma cells reduced tumorigenicity. Our study represents the first in-depth analysis of PTP family expression in diffuse glioma subtypes and warrants further investigations into PTP-dependent signaling events as new entry points for improved therapy.
KeywordsGlioblastoma Astrocytoma EGFR Oligodendroglioma IDH1 DUSP26 MTMR4 PTEN PTP PTPRM PTPRN2 PTPRT PTPRZ1 Malignancy
Gliomas arise from glial (precursor) cells and represent the most frequent type of primary brain tumor. The vast majority is classified as diffuse gliomas, reflecting their infiltrative growth into the brain parenchyma along myelinated axon bundles and blood vessels [1, 2]. Diffuse gliomas are histopathologically graded into WHO grades II-IV, with grade IV astrocytoma (glioblastoma) being most malignant . Their diffuse infiltrative character severely complicates glioma treatment. Whereas surgery and radiotherapy mostly address the central part of the tumor, cells that have migrated into the surrounding brain tissue remain relatively unaffected and give rise to recurrences. Although glioblastoma patient survival time has slightly improved over the past decades, the prospect with current treatment is only a median 15 months following diagnosis , urging for the development of novel drug modalities to combat these tumors.
Diffuse gliomas can also be subdivided into several classes based on molecular criteria. Lower grade (grade II-III) gliomas and secondary grade IV glioblastomas (arising from lower grade gliomas) display mutations in isocitrate dehydrogenase genes (IDH1/2) in 80 % of the cases . Primary glioblastomas (originating de novo) are wild type for IDH1 and IDH2 . Presence of these mutations is now often used as an independent predictor of prolonged survival [7, 8]. Glioblastomas have been further subdivided on the basis of transcriptional profiles as cancers of the pro-neural, neural, mesenchymal and classical subtype [9, 10]. Pro-neural and neural glioblastomas display IDH1/2 mutations, whereas classical and mesenchymal gliomas often show loss of chromosome 10 (containing the tumor suppressor gene PTEN at 10q23), and/or amplification of chromosome 7 (containing the proto-oncogene EGFR at 7p12) . Constitutively active mutant versions of receptor tyrosine kinases (RTKs), e.g. EGFR (EGFRvIII), PDGFRA or MET (METΔ7–8 ), are commonly encountered in primary glioblastomas . As a result, in about 90 % of primary glioblastomas the RTK/RAS/PI3K signaling pathway is affected by mutations , underscoring deregulated phosphotyrosine-based signaling as a major driver in glioma etiology.
Protein tyrosine phosphorylation is controlled via balanced actions of protein tyrosine kinases and protein tyrosine phosphatases (PTPs) . Tumor cell biological research has mostly focused on the druggable kinases, and many tyrosine kinase genes are listed as proto-oncogenes (e.g. EGFR, MET, ERBB2, and PDGFRA) . Consequently, investigations on the counteracting, hence potentially tumor suppressive, PTPs [13, 14] lag behind. A well-known exception is PTEN [15, 16], the PTP family member that is mutated or deleted in 41 % of primary glioblastomas  and whose tumor suppressive action is linked to its phospholipid phosphatase activity . Although seemingly counter-intuitive, it has become clear that some PTPs (e.g. SHP2, DUSP6 and PTPRZ1) behave as oncogenes [18, 19]. Thus, the relevance of specific PTPs in certain cancer types is emerging [19, 20] but PTPome studies in diffuse glioma specimens are lacking so far.
We performed expression profiling of the PTP gene family in diffuse glioma samples of different grades and genetic backgrounds, and found that overall PTP expression levels decrease with increasing glioma grade and appear independent of characteristic genetic aberrations. Furthermore, for two of the PTPs whose expression levels differed consistently and significantly between lower grade (II-III) and high grade (IV, glioblastoma) gliomas and correlated with patient survival, we could demonstrate tumor suppressive activity. These findings highlight PTP impact on glioma tumorigenicity.
Materials and methods
Tumor samples and cell lines
Analyses of patient material were performed in accordance with the guidelines of the local medical ethical committee of the Radboud University Medical Center, Nijmegen, The Netherlands. Patient glioma samples (n = 83), comprising 19 grade II, 15 grade III and 49 grade IV tumors, were obtained from the archives of the Radboud University Medical Center [see Additional file 1]. Histologically normal brain control tissue (temporal neocortex and white matter) was obtained from surgery on epilepsy patients at the VU University Medical Center, Amsterdam, The Netherlands.
HEK293FT cells were purchased from Invitrogen, LN-229 glioblastoma cells were from the American Type Culture Collection, and U-251 MG cells were kindly provided by Joost Schalkwijk (Radboudumc). Derivation and use of the xenograft-derived glioblastoma cell model E98 has been described elsewhere . U-251 MG, LN-229, HEK293FT and E98 cells were cultured in high glucose Dulbecco’s Modified Eagle Medium (Life Technologies; cat.no.11960–044) supplemented with 10 % Fetal Bovine Serum (PAA laboratories; Cat.no. A15–101), 4 mM glutamine and 1 mM pyruvate (Life Technologies), at 37 °C in a humidified incubator under 7.5 % CO2. For some experiments 1 μM Gefitinib (Selleck Chemicals), 4 mM 2-hydroxyglutarate (Sigma H8378), or DMSO (solvent control) was added 48 h prior to isolation.
Total RNA was isolated from frozen tumor material or from cell lines. Ten 4 μm cryosections were homogenized in 1.3 ml RNA-Bee (Tel-Test Inc., cs104B) through mechanical dissociation. Next, 260 μl chloroform was added and samples were vortexed and incubated on ice for 15 min before being centrifuged at maximum speed for 15 min at 4 °C. The top clear supernatant was transferred into a new tube and RNA was precipitated with 1 volume of isopropanol. The RNA pellet was dissolved in 400 μl NSE (50 mM NaAc, pH 5.2; 0.2 % SDS; 2 mM EDTA in RNase-free water) and re-precipitated with 1 ml ethanol and collected by centrifugation (30 min at 4 °C, 14,000 g). The RNA pellet was briefly air-dried and dissolved in 30 μl of RNase-free water. Concentrations were determined on a NanoVue plus spectrophotometer (GE Healthcare).
Reverse transcriptase reactions were performed using the iScripttm cDNA synthesis kit (Bio-Rad). RNA concentrations that were used ranged from 20 ng/μl to 2 μg/μl, and applied RNA and RT-mix volumes were according to the manufacturers’ specifications. Synthesis of cDNA was performed by incubating for 5 min at 25 °C, followed by 30 min at 42 °C and 5 min at 85 °C. Resulting cDNA samples were diluted in distilled water (1:3 – 1:5) and stored at −80 °C until further use.
Validated qPCR primer sets for mRNA transcripts from 91 PTPs and housekeeping genes ACTB, CANX and PSMB (described in ) were purchased from Qiagen and SABiosciences. Transcript abundance was determined real-time in 10 μl reactions using SYBR GREEN (Bio-Rad) and 3 μl of the diluted cDNA samples on a CFX96™ system using the C1000™ Thermal Cycler (Bio-Rad). PCR reactions were initialized at 95 °C for 15 min and followed by 40 cycles of 15 s at 95 °C and 40 s at 60 °C. To monitor primer specificity, at the end of the last cycle a melting curve with increments of 0.5 °C was recorded between 60 °C and 95 °C. PTP transcript amounts were normalized to housekeeping gene expression levels (ΔCt). Normalization to any of three housekeeping genes gave comparable results, and values normalized towards ACTB are used here. The ΔCt values are presented relative to the transcript levels in control brain tissue, according to the ΔΔCt method of Livak and Schmittgen . The average expression level for a given PTP was considered to be meaningfully different between tumor grades when p < 0.05 and │ΔΔCt│ > 1.5.
Tissue micro-array staining
Immunohistochemical stainings were performed on 4-μm sections of tissue micro-arrays (TMAs) containing formalin-fixed paraffin-embedded (FFPE) glioma (lower grade, n = 16; high grade, n = 47). Slides were deparaffinized and rehydrated according to standard methods. Heat-induced epitope retrieval (HIER) was performed in TRIS/EDTA buffer, pH 9 (Klinipath), for 10 min in a microwave oven at 180 W. After cooling down, endogenous peroxidase activity was blocked using 3 % H2O2 in PBS for 10 min at room temperature. After washing with PBS, primary antibodies against PTPRZ1 (BD Transduction Laboratories; #610179), DUSP26 (Atlas antibodies; HPA018221), MTMR4 (Pierce; #Pa5–13711), or PTPRT (Ab frontier; #LFMA0345 clone T20–3C7) were diluted in Normal Antibody Diluent (ImmunoLogic) and incubated for 1 h at room temperature. Slides were then washed with PBS and incubated with Poly-HRP-GAM/R/R IgG (ImmunoLogic) as secondary antibody for 30 min at room temperature prior to detection with 3,3-diaminobenzidine staining (Power-DAB, ImmunoLogic). All sections were counter-stained with haematoxylin and mounted in QuickD Mounting Medium (Klinipath). Sections were scored by two observers (AB and KV), after instructions and partly supervised by a neuropathologist (PW), for intensity and percentage of cells stained. Samples were binned into three classes: high staining (high intensity staining, in the majority of tumor cells), medium staining (lower intensity staining, in a considerable amount of tumor cells) and low or no staining (low intensity or no staining, in a minor part of the tumor cells). For DUSP26 we defined a sample as low staining if there were positive cells in the sample and as no staining if not a single stained tumor cell was present.
Generation of cells with aberrant EGFR expression
Third generation lentiviral constructs that drive expression of wildtype (EGFR WT) or variant III (EGFRvIII) epidermal growth factor receptor were generated as follows. A NotI linker (self-annealed 5′- CTA GTC CGG GCG GCC GCC CGG A–3′) was ligated into the cleaved XbaI site in PRK5-EGFR-WT (a mammalian expression vector harboring wild-type human EGFR cDNA, kindly provided by Reiner Lammers, Tübingen, Germany). The resulting plasmid was digested using NotI and HindIII and the EGFR WT-encoding fragment was ligated into NotI/HindIII-digested pENTR/NotI-XhoI . EGFRvIII cDNA was obtained using primers fw: 5′- GAT ATC ACC ATG CGA CCC TCC GGG–3′ and rev: 5′- CAG CGG CCG CTC ATG CTC CAA TAA ATT CAC TG–3′ in a PCR reaction with an EGFRvIII-positive glioblastoma cDNA pool as template. The amplicon was digested with NotI and resulting NotI/blunt-end EGFRvIII cDNA was ligated into HincII/NotI-digested vector pENTR-1A (Invitrogen). Resulting pENTR plasmids were subjected to an LR reaction (Gateway® LR clonasetm II Enzyme mix; Invitrogen # 11791–020) with pLenti6.2/V5-DEST as destination vector, and pENTR/NotI-XhoI served to generate the empty vector control.
Lentiviruses were generated in HEK293FT cells as described earlier . Small samples of virus-containing medium were diluted 1:2 with methanol to monitor viral titers by dot blot on nitrocellulose membranes (0.45 μM Schleicher&Schuell #401196). After blocking for 30 min with 3 % non-fat dry milk in TBST (10 mM Tris–HCl, pH 8.0; 150 mM NaCl; 0.1 % Tween-20 (Sigma-Aldrich)) membranes were incubated with anti-VSV monoclonal antibody P5D4  in blocking buffer for 1 h. Membranes were then washed three times for 10 min in TBST and subsequently incubated with IRDye® 800CW-conjugated goat-anti-mouse IgG (Li-COR Biosciences, #926–32350) for 1 h at room temperature in the dark. Blots were analyzed using the Odyssey imaging system (LI-COR Biosciences).
U-251 MG glioblastoma cells were transduced with lentiviral expression constructs for EGFR WT, EGFRvIII or empty vector control, and cultured in blasticidin-containing medium for 2 weeks. Resulting cell pools were cultured for 48 h in presence or absence of Gefitinib before being harvested for RNA and protein isolation. RNA was used for qRT-PCR as described above, to quantify transcript levels of three genes; PTPRT, DUSP16 and PTPRG. Proteins from cell lysates were size-separated on 8 % SDS-PAGE gels and immunoblotted as described previously  except that for detection of phosphorylation-dependent epitopes 5 % Bovine Serum Albumin in TBST was used as blocking and incubation buffer. Primary antibodies were directed against EGFR (rabbit 1:2000 SC03; #1005), pTyr (mouse 1:1000 pY20; sc508), and GAPDH (rabbit 1:5000; Cell Signaling Technology #2118). Secondary antibodies were IRDye-680- or −800-conjugated goat-anti-mouse and goat-anti-rabbit antibodies. Detection was performed on the Odyssey imaging system.
Generation of glioblastoma cells with PTEN deletion
PTEN-deficient LN-229 cells were created using CRISPR/Cas9 technology and a previously described PTEN targeting sequence  aimed at the start of the coding sequence in exon 1 of PTEN. Oligonucleotides that contained the target sequence and overlapped with sequences in the receiving guide RNA expression plasmid (PTEN-CRISPR-1-Fw: 5′- TTT CTT GGC TTT ATA TAT CTT GTG GAA AGG ACG AAA CAC CGG ATC GTT AGC AGA AAC AAA -3′; PTEN-CRISPR-1-Rev: 5′- GAC TAG CCT TAT TTT AAC TTG CTA TTT CTA GCT CTA AAA CTT TGT TTC TGC TAA CGA TCC -3′) were annealed and extended using Klenow DNA Polymerase fragment. The PTEN CRISPR-1 heteroduplex was introduced into the AflII-linearized gRNA cloning vector (a gift from George Church; Addgene plasmid # 41824 ) using Gibson assembly (New England Biolabs) according to manufacturer’s instructions. PTEN wildtype LN-229 glioma cells were transfected with a 1:1:0.1 mixture of the PTEN-targeting gRNA plasmid, plasmid hCas (a gift from George Church; Addgene plasmid # 41815 ), and a GFP expression plasmid (pLenti6/shSCR-GFP ), respectively, using jetPRIME (Polyplus; Cat.no 114–07). Individual clones were derived via limiting dilution and tested for successful CRISPR/Cas9-mediated PTEN gene editing using a T7 endonuclease I assay. Briefly, cells were lysed (100 mM Tris–HCl, pH 7.5; 5 mM EDTA; 0.2 % SDS; 200 mM NaCl; 100 μg/mL Proteinase K) overnight at 50 °C. Subsequently, samples were cleared by centrifugation (5 min, 8,000 rpm) and DNA in the supernatant was precipitated using an equal volume of isopropanol, washed once with 70 % ethanol, air-dried and finally dissolved in T10E0.1 (10 mM Tris–HCl, pH 8.0; 0.1 mM EDTA). Genomic DNA was then used in a PCR reaction with Q5 polymerase (NEB) containing Q5 GC enhancer according to manufacturer’s instructions. Primers used were: PTEN-T7-fw: 5′-GGC TGG GAA CGC CGG AGA GT-3′ and PTEN-T7-rev: 5′-CTA TCC TAC ATC ACC CCA TAG GG-3′. Amplicons were purified using the NucleoSpin® Gel and PCR Clean-up kit (Macherey-Nagel) and subjected to a T7 endonuclease I digestion assay (NEB, M0302S) for 15 min according to manufacturer’s instructions. Clones displaying edited PTEN alleles were further tested for PTEN status through qRT-PCR and immunoblot analyses as described above. Primary antibodies used were: rabbit anti-PTEN (1:2,000, CST#9559), mouse anti-tubulin monoclonal E7 (1:5,000; DSHB, University of Iowa), and rabbit anti-pAKT-Ser473 (1:2,000; CST #4058).
PTPRT and DUSP26 expression constructs
PTPRT cDNA encompassing the full-length coding sequence was isolated from plasmid pGW-PTPRT (a kind gift from Jae-Ran Lee, Daejeon, Korea) using BamHI and NotI, and ligated into the corresponding sites in vector pENTR/NotI-XhoI . Similarly, plasmid pEF-HA-DUSP26 (a kind gift from Rafael Pulido, Barakaldo, Spain) was digested using EcoRI and XbaI and the DUSP26 open reading frame was ligated into pENTR/NotI-XhoI. An empty vector control (EV) was constructed as described above. An EGFP vector control was generated by PCR (Fw: 5′- TAT AGG ATC CAT GGT GAG CAA GGG CGA GG -3′; Rev: 5′- TAT ACT CGA GTT ACT TGT ACA GCT CGT CC -3′) using pEGFP as template, and the obtained amplicon was digested with BamHI and XhoI and ligated into pENTR/NotI-XhoI. All resulting plasmids were used in LR Gateway® reactions (Invitrogen) using pLenti6/PGK-DEST-TagRFP  as destination vector. Obtained expression plasmids were used in HEK293FT cells to produce lentiviruses, as described above. E98 human glioblastoma cells  were lentivirally transduced twice, and 3 days later further processed for glioma cell behavior assays.
Proliferation and migration assays
E98 cells were grown on collagen I-coated (10 μg/cm2; Invitrogen) coverslips to 60–80 % confluency over 48 h, and incubated for 1 h with culture medium containing 10 μM EdU (5-ethynyl-2′-deoxyuridine). EdU incorporation was visualized using the click-iT® EdU Imaging kit (Thermo Fisher Scientific, #C10086) via the manufacturers’ instruction. Coverslips were mounted on microscope slides in DAPI-containing Mowiol (Sigma-Aldrich) and images were collected on a Leica DMRA Fluorescence microscope, equipped with a DFC340 FX CCD camera, using 40x and 63x objectives. DAPI- and EdU-positive nuclei were counted automatically using FIJI software .
Migration of E98 cells was assessed in spheroid outgrowth assays as follows. E98 spheroids were generated in hanging drops using methylcellulose (12 mg/mL; Sigma, M6385) in DMEM supplemented with 10 % FCS (2500 cells per spheroid). The next day, individual spheroids were seeded in a 96-well imaging culture dish (BD Falcon, #353219) on top of a confluent mouse astrocyte layer in Matrigel-coated (30 μg/mL PBS; BD Biosciences, #356237) culture wells. 24 h later, cells were fixed and fluorescent (tagRFP) images were collected. Average migration distance of cells from spheroids (n > 37), calculated as change in radius of the spheroid over 24 h, were analyzed semi-automatically using FIJI software.
In silico analyses
REpository for Molecular BRAin Neoplasia DaTa (REMBRANDT) was used to determine Kaplan-Meier curves based on PTP transcript levels (as determined via micro-array analyses) using the online Project Betastasis representation tool (http://betastasis.com/glioma/rembrandt/). Threshold was set at median expression for the respective PTP transcript. RNA-seq data from two different TCGA datasets (Brain Lower Grade Glioma  and Glioblastoma Multiforme; Provisional) were downloaded using cBioPortal (http://www.cbioportal.org/index.do) [30, 31]. These represented RNA-seq data of 166 IDH1 WT glioblastoma samples, 9 IDH1 mutant glioblastoma samples, 221 IDH1 mutant lower grade glioma samples and 309 IDH1 WT lower grade glioma samples.
Statistical significance was tested using non-paired two-tailed Student’s t-test or non-parametric Kruskal-Wallis one-way analysis of variance with post-hoc Dunn’s Multiple Comparison Test. Patient survival data was automatically calculated by Project Betastasis using a log-rank test. Heat map was generated in R (version 3.0.1).
Differential expression of PTP transcripts between lower grade and high grade glioma samples
PTP expression patterns correlate with disease outcome
Characteristic genetic aberrations in gliomas do not drive PTP expression
PTEN-inactivating deletions or mutations represent another common aberration in primary glioblastoma samples and could potentially explain the pattern observed on the PTP transcripts. To genocopy PTEN loss, the PTEN WT glioblastoma cell line LN-229 was subjected to CRISPR/Cas9-mediated genome editing and resulting clones were inspected for PTEN levels and activity, respectively, as witnessed by PTEN and phospho-AKT levels under low-serum conditions (Fig. 5b). Although we successfully generated PTEN-deficient LN-229 derivates and appropriate controls, no significant changes in expression levels were observed for the three representative reporter PTPs that were tested (Fig. 5d).
We also monitored effects of the oncometabolite 2-hydroxyglutarate (2-HG) that is produced in lower grade gliomas and secondary glioblastomas as a result of the characteristic IDH1R132H mutation and inhibits DNA and histone demethylation. Tumor-relevant levels of 2-HG were added for 48 h to U-251 MG (IDH1 wild type) glioblastoma cells and potential changes in DUSP16, PTPRG and PTPRT expression levels were monitored. Addition of 2-HG significantly reduced the mRNA levels for two out of the three PTPs measured (Fig. 5e). These results are in line with epigenetic regulation of DUSP16 and PTPRG but do not explain the generally higher PTP mRNA levels in IDH-mutant low-grade gliomas. RNAseq data from WT or IDH1 R132H-containing glioma xenografts E434 and E478 [27, 33] also do not point to an IDH-mutation associated difference in expression for DUSP16, PTPRG and PTPRT (WPJL, unpublished data). Collectively, these data make it rather unlikely that genetic alterations characterizing lower grade glioma and glioblastoma subgroups are major determinants of the general PTP expression pattern observed (Fig. 1a,b).
Overexpression of DUSP26 or PTPRT in E98 glioblastoma cells results in reduced tumorigenicity
To also monitor potential DUSP26 or PTPRT effects on glioblastoma cell migration, we turned to a previously established spheroid outgrowth assay [4, 28]. Homogeneously-sized spheroids of lentivirally transduced E98 cells, generated using the hanging-drop method, were placed on a layer of mouse astrocytes and the migratory performance of outgrowing fluorescent E98 cells was calculated as the change in radius of the spheroid over 24 h. As compared to GFP-expressing control cells, the expression of either PTPRT or DUSP26 significantly reduced migration of E98 cells from the spheroid edge (Fig. 6e-f). DUSP26 and PTPRT thus impinge on both growth and motility of glioblastoma cells. Collectively, these data underscore their involvement as tumor suppressors in glioma, in line with the mRNA profiling data (Figs. 1, 3, 4).
Diffuse gliomas are invariably fatal tumors, with most glioblastoma patients surviving only up to 15 months after diagnosis , stressing the need for more effective treatment modalities. In this study we focused on the protein tyrosine phosphatase gene family as alternative entry points for treatment strategies. Seven out of ninety-one investigated PTPs displayed significantly different mRNA expression levels when comparing lower and high grade gliomas. Furthermore, for four out of these seven, expression levels correlated with patient survival probability, underscoring their role in tumor behavior. Importantly, overexpression of the two PTPs that showed the largest expression difference in our qPCR and in silico analyses, DUSP26 and PTPRT, resulted in reduced glioblastoma cell proliferation and migration, supportive of a tumor suppressive role.
In general, PTP expression levels in lower grade gliomas resembled those in normal brain tissue whereas levels in glioblastomas were generally lower. This overall reduction in PTP expression with increasing grade may relate to differences in the cell type of origin. Alternatively, it may be imposed by characteristic genetic alterations acquired during glioma development. We found that experimental introduction of EGFR hyperactivity, or of PTEN inactivation, in a glioblastoma cell line had no effect on overall PTP expression patterns. Furthermore, the transcriptional impact of 2-HG on PTP expression in IDH1 wildtype glioblastoma cells was at best opposite to the expected, and not able to tone down PTP mRNA levels to those in lower grade gliomas or normal brain tissue. These findings suggest that other differences in tumor pathogenesis between lower and high grade tumors are the major determinants for the PTP expression pattern. The overall reduced PTP levels in glioblastomas as compared to histologically normal brain and lower grade tumor levels, points to PTP agonists as a potential option for glioblastoma treatment.
In our screening rationale, stringent criteria led to a selection of seven PTPs that differed significantly in their expression between glioblastomas and lower grade gliomas. Our data revealing that PTEN expression was low in high grade gliomas concords with literature data on PTEN inactivation or absence in over one-third of high grade gliomas, resulting in increased proliferation and survival through activation of the PI3K/AKT pathway [10, 35]. PTPRM has also been proposed as a glial tumor suppressor (e.g. [19, 32]) and is consistently down-regulated in high grade samples in our screen. Together with its subfamily members PTPRU and PTPRK, they are amongst the most frequently mutated PTPs in a number of cancer types [16, 19, 20, 32, 36]. Some PTPs that had previously been implicated in primary glioblastoma etiology  were not selected by our screen, although some only just failed to reach our criteria (e.g. PTPN11, encoding SHP2). This may reflect that such PTPs could govern processes equally important to lower grade and high grade tumors.
We found that PTPRZ1 expression is higher in lower grade tumors than in high grade tumors, but in online RNA-seq data the inverse was observed. This apparent discrepancy may result from histopathological differences in the tumor collections. In the RNA-seq dataset, a third of the lower grade samples represent astrocytomas while in our sample cohort lower grade tumors consisted mostly of oligodendroglial neoplasms with 1p/19q co-deletion. PTPRZ1 is highly expressed in oligodendrocyte precursor cells and is involved in normal oligodendrocyte development [37, 38], thus providing rationale for the high PTPRZ1 expression in our sample cohort. Irrespective, PTPRZ1 expression is consistently up-regulated in all glioma specimens. We and others have shown that the encoded transmembrane PTP acts as an oncogene in gliomas [21, 39, 40] and, interestingly, two HLA-presented peptides derived from PTPRZ1’s extracellular part are included in an experimental glioma vaccine that is currently explored in phase II clinical trials .
We additionally identified PTPRT, DUSP26, PTPRN2 and MTMR4 in our screen for glioma-relevant PTPs. PTPRT belongs to the same cell adhesion molecule-like PTP subfamily as PTPRM, PTPRU and PTPRK. Although compelling evidence has put forward PTPRT as a colon cancer susceptibility gene [20, 42–44], thus far no reports have linked PTPRT to gliomas. PTPRT mRNA levels are quite low in normal brain tissue and lower grade gliomas but are still on average 40-fold higher than those in high grade glioblastomas. In line with a tumor suppressive role, we found that PTPRT overexpression significantly suppressed glioblastoma cell proliferation and migration. The finding that the closely related PTPRT and PTPRM show distinct expression patterns in lower versus high grade gliomas (PTPRU just failed to reach the criteria in the first cohort) and correlate with patient survival, warrants further studies on their role in glioma-associated signaling pathways.
DUSP26 represents a dual-specificity phosphatase, i.e. it is able to dephosphorylate not only phosphotyrosines but also phospho-serine and -threonine residues in substrate proteins. In neuroblastoma cells DUSP26 dephosphorylates and inhibits the tumor suppressor p53 . Furthermore, DUSP26 amplification has been found in thyroid cancer, and DUSP26 knockdown impeded growth of anaplastic thyroid cancer cells . Collectively, this indicates that DUSP26 acts as an oncogene, and contrasts with our observation of low DUSP26 levels in glioblastomas in comparison with lower grade glioma samples and normal brain tissue. Importantly, we found that DUSP26 overexpression in E98 glioblastoma cells resulted in decreased growth and motility. Other investigators also reported on down-regulated DUSP26 levels in glioblastoma  and Tanuma and colleagues additionally provided evidence that DUSP26 facilitates catenin/cadherin delivery to cell-cell junction sites, hence cell-cell adhesiveness . It is becoming increasingly clear that diffuse glioma cells form functional networks via cell-cell contacts  and in such a context DUSP26 protein may indeed have tumor suppressive activity. Together with PTPRM and PTPRT, DUSP26 therefore holds promise for therapeutic intervention.
Intriguingly, PTPRN2 and MTMR4 both encode PTPs that, like PTEN, have phosphoinositides as substrates. PTPRN2 preferentially dephosphorylates phosphatidylinositol 4,5-bisphosphate (PI(4,5)P)  and MTMR4 prefers phosphatidylinositol 3-phosphate (PI(3)P) . Whereas the tumor suppressive impact of PTEN’s phospholipid phosphatase activity is beyond doubt , only limited and contradicting data exist regarding PTPRN2 involvement in cancer. The hypermethylation of the PTPRN2 promoter region in glioblastomas and lung adenocarcinomas [53, 54] suggests tumor suppressor roles. In contrast, high PTPRN2 expression has been reported to correlate with poor clinical outcomes in breast cancer cases . Furthermore, experimental down-regulation of PTPRN2 in metastatic breast cancer cells inhibited migratory potential and yielded smaller tumors in vivo , rather pointing to oncogenic potential. The tumor biological impact of PTPRN2 and MTMR4 may thus be context-dependent, and how these PTPs feed into glioma biology needs further investigation.
Taken together, our study – to our knowledge the first comprehensive PTP mRNA profiling in diffuse gliomas – highlights seven PTPs (DUSP26, MTMR4, PTEN, PTPRM, PTPRN2, PTPRT and PTPRZ1) of which expression decreases with increasing malignancy. Largest expression differences between lower grade and high grade diffuse gliomas were found for DUSP26 and PTPRT. Their expression correlated with patient survival probability, and overexpression inhibited glioblastoma cell growth and motility. Together with PTPs that impinge upon cellular contacts (DUSP26, PTPRZ1) and phospholipid signaling (MTMR4, PTEN and PTPRN2) they provide novel cues to explore and design glioma treatment options.
We thank Sandra Boots-Sprenger, Judith Jeuken and Sarmistha Guha for assistance during the starting phase of the study, Marieke Willemse, Huib Croes and Jack Fransen for Microscopical Imaging Centre core support, and Bé Wieringa for useful discussions and critical reading of the manuscript. Rafael Pulido, Reiner Lammers and Jai-Ran Lee are gratefully acknowledged for providing cDNA constructs, Joost Schalkwijk for supply of U-251 MG cells, and Pavlo Grytsenko for providing the astrocyte feeder layer protocol. We also thank Mark ter Laan and Radboudumc coworkers for supportive tumor specimens biobanking activities. This work was supported in part by a Marie Curie Research Training Network grant (PTPNET/MRTN-CT-2006–035830) and a UMC St Radboud research grant (R0000761).
Conception and design: PW, WPJL, WJAJH. Development of methodology: AMB, JTGS, KV, ME. Acquisition of data: AMB, ACN, JTGS, KV, JAFP, CBCP. Analysis and interpretation of data: AMB, WPJL, WJAJH. Writing, review and/or revision of the manuscript: AMB, RHH, PW, WPJL, WJAJH. Study supervision: WPJL, WJAJH. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Claes A, Idema AJ, Wesseling P. Diffuse glioma growth: a guerilla war. Acta Neuropathol. 2007;114:443–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Kohler BA, Ward E, McCarthy BJ, Schymura MJ, Ries LAG, Eheman C, Jemal A, Anderson RN, Ajani UA, Edwards BK. Annual report to the nation on the status of cancer, 1975–2007, featuring tumors of the brain and other nervous system. J Natl Cancer Inst. 2011;103:714–36.View ArticlePubMedPubMed CentralGoogle Scholar
- Louis DN, Ohgaki H, Wiestler OD, Cavenee WK, Burger PC, Jouvet A, Scheithauer BW, Kleihues P. The 2007 WHO classification of tumours of the central nervous system. Acta Neuropathol. 2007;114:97–109.View ArticlePubMedPubMed CentralGoogle Scholar
- Stupp R, Mason WP, van den Bent MJ, Weller M, Fisher B, Taphoorn MJ, Belanger K, Brandes AA, Marosi C, Bogdahn U, Curschmann J, Janzer RC, Ludwin SK, Gorlia T, Allgeier A, Lacombe D, Cairncross JG, Eisenhauer E, Mirimanoff RO; European Organisation for Research and Treatment of Cancer Brain Tumor and Radiotherapy Groups; National Cancer Institute of Canada Clinical Trials Group. Radiotherapy plus concomitant and adjuvant temozolomide for glioblastoma. N Engl J Med. 2005;352:987–96.Google Scholar
- Balss J, Meyer J, Mueller W, Korshunov A, Hartmann C, von Deimling A. Analysis of the IDH1 codon 132 mutation in brain tumors. Acta Neuropathol. 2008;116:597–602.View ArticlePubMedGoogle Scholar
- Ohgaki H, Kleihues P. The definition of primary and secondary glioblastoma. Clin Cancer Res. 2013;19:764–72.View ArticlePubMedGoogle Scholar
- Boots-Sprenger SHE, Sijben A, Rijntjes J, Tops BBJ, Idema AJ, Rivera AL, Bleeker FE, Gijtenbeek AM, Diefes K, Heathcock L, Aldape KD, Jeuken JWM, Wesseling P. Significance of complete 1p/19q co-deletion, IDH1 mutation and MGMT promoter methylation in gliomas: use with caution. Mod Pathol. 2013;26:922–9.View ArticlePubMedGoogle Scholar
- Van Den Bent MJ, Dubbink HJ, Marie Y, Brandes AA, Taphoorn MJB, Wesseling P, Frenay M, Tijssen CC, Lacombe D, Idbaih A, Van Marion R, Kros JM, Dinjens WNM, Gorlia T, Sanson M. IDH1 and IDH2 mutations are prognostic but not predictive for outcome in anaplastic oligodendroglial tumors: a report of the european organization for research and treatment of cancer brain tumor group. Clin Cancer Res. 2010;16:1597–604.View ArticlePubMedGoogle Scholar
- Verhaak RGW, Hoadley KA, Purdom E, Wang V, Wilkerson MD, Miller CR, Ding L, Golub T, Jill P, Alexe G, Lawrence M, Kelly MO, Tamayo P, Weir BA, Gabrie S, Winckler W, Gupta S, Feiler HS, Hodgson JG, James CD, Sarkaria JN, Brennan C, Kahn A, Spellman PT, Wilson RK, Terence P, Gray JW, Meyerson M, Getz G, Perou CM, et al. An integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR and NF1. Cancer Cell. 2011;17:1–25.Google Scholar
- Brennan C, Verhaak R, McKenna A, Campos B, Noushmehr H, Salama S, Zheng S, Chakravarty D, Sanborn J, Berman S, Beroukhim R, Bernard B, Chin L. The somatic genomic landscape of glioblastoma. Cell. 2013;155:462–77.View ArticlePubMedPubMed CentralGoogle Scholar
- Navis AC, van Lith SAM, van Duijnhoven SMJ, de Pooter M, Yetkin-Arik B, Wesseling P, Hendriks WJAJ, Venselaar H, Timmer M, van Cleef P, van Bergen En Henegouwen P, Best MG, Wurdinger TD, Tops BBJ, Leenders WPJ. Identification of a novel MET mutation in high-grade glioma resulting in an auto-active intracellular protein. Acta Neuropathol. 2015;130:131–44.View ArticlePubMedPubMed CentralGoogle Scholar
- Sugawa N, Ekstrand AJ, Jamest CD, Collins VP. Identical splicing of aberrant epidermal growth factor receptor transcripts from amplified rearranged genes in human glioblastomas. Proc Natl Acad Sci. 1990;87:8602–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Tonks NK. Protein tyrosine phosphatases: from genes, to function, to disease. Nat Rev Mol Cell Biol. 2006;7:833–46.View ArticlePubMedGoogle Scholar
- Hendriks WJAJ, Elson A, Harroch S, Stoker AW. Protein tyrosine phosphatases: functional inferences from mouse models and human diseases. FEBS J. 2008;275:816–30.View ArticlePubMedGoogle Scholar
- Verhoeff JJC, van Tellingen O, Claes A, Stalpers LJA, van Linde ME, Richel DJ, Leenders WPJ, van Furth WR. Concerns about anti-angiogenic treatment in patients with glioblastoma multiforme. BMC Cancer. 2009;9:444.View ArticlePubMedPubMed CentralGoogle Scholar
- Kaur H, Burden-Gulley SM, Phillips-Mason PJ, Basilion JP, Sloan AE, Brady-Kalnay SM. Protein tyrosine phosphatase mu regulates glioblastoma cell growth and survival in vivo. Neuro Oncol. 2012;14:561–73.View ArticlePubMedPubMed CentralGoogle Scholar
- Pulido R, Stoker AW, Hendriks WJAJ. PTPs emerge as PIPs: protein tyrosine phosphatases with lipid-phosphatase activities in human disease. Hum Mol Genet. 2013;22:R66–76.View ArticlePubMedGoogle Scholar
- Hendriks WJAJ, Pulido R. Protein tyrosine phosphatase variants in human hereditary disorders and disease susceptibilities. Biochim Biophys Acta. 1832;2013:1673–96.Google Scholar
- Navis AC, van den Eijnden M, Schepens JTG, Hooft Van Huijsduijnen R, Wesseling P, Hendriks WJAJ. Protein tyrosine phosphatases in glioma biology. Acta Neuropathol. 2010;119:157–75.View ArticlePubMedGoogle Scholar
- Julien SG, Dubé N, Hardy S, Tremblay ML. Inside the human cancer tyrosine phosphatome. Nat Rev Cancer. 2011;11:35–49.View ArticlePubMedGoogle Scholar
- Bourgonje AM, Navis AC, Schepens JTG, Verrijp K, Hovenstad L, Hilhorst R, Harroch S, Wesseling P, Leenders WPJ, Hendriks WJAJ. Intracellular and extracellular domains of protein tyrosine phosphatase PTPRZ-B differentially regulate glioma cell growth and motility. Oncotarget. 2014;5:8690–702.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmidt F, van den Eijnden M, Pescini Gobert R, Saborio GP, Carboni S, Alliod C, Pouly S, Staugaitis SM, Dutta R, Trapp B, Hooft Van Huijsduijnen R. Identification of VHY/Dusp15 as a regulator of oligodendrocyte differentiation through a systematic genomics approach. PLoS One. 2012;7, e40457.View ArticlePubMedPubMed CentralGoogle Scholar
- Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(−Delta Delta C(T)) Method. Methods. 2001;25:402–8.View ArticlePubMedGoogle Scholar
- Kreis TE. Microinjected antibodies against the cytoplasmic domain of vesicular stomatitis virus glycoprotein block its transport to the cell surface. EMBO J. 1986;5:931–41.PubMedPubMed CentralGoogle Scholar
- Xue W, Chen S, Yin H, Tammela T, Papagiannakopoulos T, Joshi NS, Cai W, Yang G, Bronson R, Crowley DG, Zhang F, Anderson DG, Sharp PA, Jacks T. CRISPR-mediated direct mutation of cancer genes in the mouse liver. Nature. 2014;514:3–7.View ArticleGoogle Scholar
- Mali P, Yang L, Esvelt K, Aach J, Guell M, DiCarlo J, Norville J, Church G. RNA-guided human genome engineering via Cas9. Science. 2013;339(80-):823–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Claes A, Schuuring J, Boots-Sprenger S, Hendriks-Cornelissen S, Dekkers M, van der Kogel AJ, Leenders WP, Wesseling P, Jeuken JW. Phenotypic and genotypic characterization of orthotopic human glioma models and its relevance for the study of anti-glioma therapy. Brain Pathol. 2008;18:423–33.View ArticlePubMedGoogle Scholar
- Schindelin J, Arganda-Carreras I, Frise E, Kaynig V, Longair M, Pietzsch T, Preibisch S, Rueden C, Saalfeld S, Schmid B, Tinevez J-Y, White DJ, Hartenstein V, Liceiri K, Tomancak P, Cardona A. Fiji: an open source platform for biological image analysis. Nat Methods. 2012;9:676–82.View ArticlePubMedGoogle Scholar
- The Cancer Genome Atlas Research Network. Comprehensive, integrative genomic analysis of diffuse lower- grade gliomas. N Engl J Med. 2015;372:2481–98.View ArticlePubMed CentralGoogle Scholar
- Cerami E, Gao J, Dogrusoz U, Gross BE, Sumer SO, Aksoy BA, Jacobsen A, Byrne CJ, Heuer ML, Larsson E, Antipin Y, Reva B, Goldberg AP, Sander C, Schultz N. The cBio cancer genomics portal: an open platform for exploring multidimensional cancer genomics data. Cancer Discov. 2012;2:401–4.View ArticlePubMedGoogle Scholar
- Gao J, Aksoy BA, Dogrusoz U, Dresdner G, Gross B, Sumer SO, Sun Y, Jacobsen A, Sinha R, Larsson E, Cerami E, Sander C, Schultz N. Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci Signal. 2013;6:l1.View ArticleGoogle Scholar
- Craig SEL, Brady-Kalnay SM. Regulation of development and cancer by the R2B subfamily of RPTPs and the implications of proteolysis. Semin Cell Dev Biol. 2015;37C:108–18.View ArticleGoogle Scholar
- Navis AC, Niclou SP, Fack F, Stieber D, van Lith S, Verrijp K, Wright A, Stauber J, Tops B, Otte-Holler I, Wevers RA, van Rooij A, Pusch S, von Deimling A, Tigchelaar W, van Noorden CJF, Wesseling P, Leenders WPJ. Increased mitochondrial activity in a novel IDH1-R132H mutant human oligodendroglioma xenograft model: in situ detection of 2-HG and α-KG. Acta Neuropathol Commun. 2013;1:18.View ArticlePubMedPubMed CentralGoogle Scholar
- Navis AC, Bourgonje A, Wesseling P, Wright A, Hendriks W, Verrijp K, van der Laak JA WM, Heerschap A, Leenders WPJ. Effects of dual targeting of tumor cells and stroma in human glioblastoma xenografts with a tyrosine kinase inhibitor against c-MET and VEGFR2. PLoS One. 2013;8, e58262.View ArticlePubMedPubMed CentralGoogle Scholar
- Pulido R, Hooft van Huijsduijnen R. Protein tyrosine phosphatases: dual-specificity phosphatases in health and disease. FEBS J. 2008;275:848–66.View ArticlePubMedGoogle Scholar
- Burgoyne AM, Palomo JM, Phillips-Mason PJ, Burden-Gulley SM, Major DL, Zaremba A, Robinson S, Sloan AE, Vogelbaum MA, Miller RH, Brady-Kalnay SM. PTPmu suppresses glioma cell migration and dispersal. Neuro Oncol. 2009;11:767–78.View ArticlePubMedPubMed CentralGoogle Scholar
- Kuboyama K, Fujikawa A, Masumura M, Suzuki R, Matsumoto M, Noda M. Protein tyrosine phosphatase receptor type z negatively regulates oligodendrocyte differentiation and myelination. PLoS One. 2012;7:e48797.View ArticlePubMedPubMed CentralGoogle Scholar
- Lamprianou S, Chatzopoulou E, Thomas J-L, Bouyain S, Harroch S. A complex between contactin-1 and the protein tyrosine phosphatase PTPRZ controls the development of oligodendrocyte precursor cells. Proc Natl Acad Sci U S A. 2011;108:17498–503.View ArticlePubMedPubMed CentralGoogle Scholar
- Müller S, Kunkel P, Lamszus K, Ulbricht U, Lorente GA, Nelson AM, von Schack D, Chin DJ, Lohr SC, Westphal M, Melcher T. A role for receptor tyrosine phosphatase zeta in glioma cell migration. Oncogene. 2003;22:6661–8.View ArticlePubMedGoogle Scholar
- Ulbricht U, Eckerich C, Fillbrandt R, Westphal M, Lamszus K. RNA interference targeting protein tyrosine phosphatase zeta/receptor-type protein tyrosine phosphatase beta suppresses glioblastoma growth in vitro and in vivo. J Neurochem. 2006;98:1497–506.View ArticlePubMedGoogle Scholar
- Dutoit V, Herold-Mende C, Hilf N, Schoor O, Beckhove P, Bucher J, Dorsch K, Flohr S, Fritsche J, Lewandrowski P, Lohr J, Rammensee HG, Stevanovic S, Trautwein C, Vass V, Walter S, Walker PR, Weinschenk T, Singh-Jasuja H, Dietrich PY. Exploiting the glioblastoma peptidome to discover novel tumour-associated antigens for immunotherapy. Brain. 2012;135:1042–54.View ArticlePubMedGoogle Scholar
- Lee JW, Jeong EG, Lee SH, Nam SW, Kim SH, Lee JY, Yoo NJ, Lee SH. Mutational analysis of PTPRT phosphatase domains in common human cancers. APMIS. 2007;115:47–51.View ArticlePubMedGoogle Scholar
- Scott A, Wang Z. Tumour suppressor function of protein tyrosine phosphatase receptor-T. Biosci Rep. 2011;31:303–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Wu C-W, Kao H-L, Li AF-Y, Chi C-W, Lin W-C. Protein tyrosine-phosphatase expression profiling in gastric cancer tissues. Cancer Lett. 2006;242:95–103.View ArticlePubMedGoogle Scholar
- Shang X, Vasudevan SA, Yu Y, Ge N, Ludwig AD, Wesson CL, Wang K, Burlingame SM, Zhao Y-J, Rao PH, Lu X, Russell HV, Okcu MF, Hicks MJ, Shohet JM, Donehower LA, Nuchtern JG, Yang J. Dual-specificity phosphatase 26 is a novel p53 phosphatase and inhibits p53 tumor suppressor functions in human neuroblastoma. Oncogene. 2010;29:4938–46.View ArticlePubMedGoogle Scholar
- Yu W, Imoto I, Inoue J, Onda M, Emi M, Inazawa J. A novel amplification target, DUSP26, promotes anaplastic thyroid cancer cell growth by inhibiting p38 MAPK activity. Oncogene. 2007;26:1178–87.View ArticlePubMedGoogle Scholar
- Patterson KI, Brummer T, Daly RJ, O’Brien PM. DUSP26 negatively affects the proliferation of epithelial cells, an effect not mediated by dephosphorylation of MAPKs. Biochim Biophys Acta. 1803;2010:1003–12.Google Scholar
- Tanuma N, Nomura M, Ikeda M, Kasugai I, Tsubaki Y, Takagaki K, Kawamura T, Yamashita Y, Sato I, Sato M, Katakura R, Kikuchi K, Shima H. Protein phosphatase Dusp26 associates with KIF3 motor and promotes N-cadherin-mediated cell-cell adhesion. Oncogene. 2009;28:752–61.View ArticlePubMedGoogle Scholar
- Osswald M, Jung E, Sahm F, Solecki G, Venkataramani V, Horstmann H, Liao Y, Syed M, Blaes J, Berghoff A, Ratliff M, Lemke D, Gömmel M, Pauli M, Häring P, Herl V, Steinhäuser C, Miletic H, Preusser M, Weiss S, Liu H, Huber PE, Kuner T, von Deimling A, Wick W, Winkler F. Brain tumor cells interconnect to a functional and resistant network. Nature. 2015;528:93–8.PubMedGoogle Scholar
- Caromile LA, Oganesian A, Coats SA, et al. The neurosecretory vesicle protein phogrin functions as a phosphatidylinositol phosphatase to regulate insulin secretion. J Biol Chem. 2010;285:10487–96.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhao R, Qi Y, Chen J, Zhao ZJ. FYVE-DSP2, a FYVE domain-containing dual specificity protein phosphatase that dephosphorylates phosphotidylinositol 3-phosphate. Exp Cell Res. 2001;265:329–38.View ArticlePubMedGoogle Scholar
- Hendriks WJAJ, Elson A, Harroch S, Pulido R, Stoker A, den Hertog J. Protein tyrosine phosphatases in health and disease. FEBS J. 2013;280:708–30.View ArticlePubMedGoogle Scholar
- Lee E-J, Rath P, Liu J, Ryu D, Pei L, Noonepalle SK, Shull AY, Feng Q, Litofsky NS, Miller DC, Anthony DC, Kirk MD, Laterra J, Deng L, Xin H-B, Wang X, Choi J-H, Shi H. Identification of global DNA methylation signatures in glioblastoma-derived cancer stem cells. J Genet genomics. 2015;42:355–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Selamat SA, Galler JS, Joshi AD, Fyfe MN, Campan M, Siegmund KD, Kerr KM, Laird-Offringa IA. DNA methylation changes in atypical adenomatous hyperplasia, adenocarcinoma in situ, and lung adenocarcinoma. PLoS One. 2011;6:e21443.View ArticlePubMedPubMed CentralGoogle Scholar
- Sorokin A, Nair B, Wei Y, Aziz K, Evdokimova V, Hung M, Chen J. Aberrant expression of proPTPRN2 in cancer cells confers resistance to apoptosis. Cancer Res. 2015;75:1846–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Sengelaub CA, Navrazhina K, Ross JB, Halberg N, Tavazoie SF. PTPRN2 and PLCb1 promote metastatic breast cancer cell migration through PI(4,5)P2 -dependent actin remodeling. EMBO J. 2016;35:62–76.View ArticlePubMedGoogle Scholar