Protein tyrosine phosphatase receptor delta acts as a neuroblastoma tumor suppressor by destabilizing the aurora kinase a oncogene

Background Protein tyrosine phosphatase receptor delta (PTPRD) is a member of a large family of protein tyrosine phosphatases which negatively regulate tyrosine phosphorylation. Neuroblastoma is a major childhood cancer arising from precursor cells of the sympathetic nervous system which is known to acquire deletions and alterations in the expression patterns of PTPRD, indicating a potential tumor suppressor function for this gene. The molecular mechanism, however, by which PTPRD renders a tumor suppressor effect in neuroblastoma is unknown. Results As a molecular mechanism, we demonstrate that PTPRD interacts with aurora kinase A (AURKA), an oncogenic protein that is over-expressed in multiple forms of cancer, including neuroblastoma. Ectopic up-regulation of PTPRD in neuroblastoma dephosphorylates tyrosine residues in AURKA resulting in a destabilization of this protein culminating in interfering with one of AURKA's primary functions in neuroblastoma, the stabilization of MYCN protein, the gene of which is amplified in approximately 25 to 30% of high risk neuroblastoma. Conclusions PTPRD has a tumor suppressor function in neuroblastoma through AURKA dephosphorylation and destabilization and a downstream destabilization of MYCN protein, representing a novel mechanism for the function of PTPRD in neuroblastoma.

Neuroblastoma is derived from primitive cells of the sympathetic nervous system, and is the most common extracranial solid tumor in children accounting for 15% of all childhood cancer deaths [13]. These tumors are particularly noted for extensive heterogeneity in clinical behaviour, ranging from spontaneous regression to aggressive clinical course and death from disease. Notably, amplification of the MYCN transcription factor is one of the most powerful adverse prognostic factors in neuroblastoma [14] and we have previously demonstrated that PTPRD is expressed at significantly lower levels in MYCN amplified neuroblastoma relative to non-MYCN amplified tumors [10]. In addition, PTPRD mRNA expression is higher in normal adrenal fetal neuroblasts, the cell of origin of neuroblastoma, relative to unfavourable neuroblastoma tumors, indicating that PTPRD down-regulation might be an important step in the development of these tumors [7,10]. Multiple mechanisms appear to exist for the down-regulation of PTPRD in neuroblastoma, including intragenic microdeletions which can include coding sequence, or in some instances be restricted to non-coding exons of an extended 5' UTR [5]. Aberrant splicing of the 5' UTR also has been noted in neuroblastoma cell lines and primary tumors, which could potentially cause destabilization of the mRNA sequence [10].
In this report, we demonstrate for the first time that experimental up-regulation of PTPRD in neuroblastoma cell lines significantly decreases cell growth and increases apoptosis. Moreover, we identify aurora kinase A, a serine/threonine kinase oncogene that is up-regulated in many forms of cancer, including high risk neuroblastoma [15], as an interaction partner of PTPRD. We further demonstrate that PTPRD has a tumor suppressor function in neuroblastoma through dephosphorylating and destabilizing AURKA, leading to a downstream decrease of MYCN protein. Our findings represent a novel mechanism of action for the function of PTPRD in neuroblastoma.

PTPRD functions as a tumor suppressor in neuroblastoma
In order to further examine the possibility that PTPRD acts as a tumor suppressor gene in neuroblastoma, we initially analyzed the levels of PTPRD mRNA transcripts in a set of 88 neuroblastoma tumors using the R2: microarray analysis and visualization platform (http://r2. amc.nl) (University of Amsterdam). Lower than median PTPRD mRNA levels were significantly associated with both poor relapse free and overall patient survival, consistent with PTPRD having a tumor suppressor function ( Figure 1A and 1B).
We then sought to determine the effects of PTPRD over-expression in neuroblastoma cell lines that have only minimally detectable PTPRD mRNA transcripts. Kelly cells (MYCN amplified) have a homozygous deletion in the PTPRD 5' UTR [5], while SHSY-5Y (MYCN non-amplified) has a 5' UTR that was generated by exon skipping [10]. Although there is no evidence for deletion or aberrant splicing of PTPRD in CHP212 (MYCN amplified), endogenous mRNA transcripts are nevertheless at very low levels. These three cell lines were transfected with a wild type PTPRD V5-tagged cDNA clone [6] to determine the impact of over-expression on cell growth. Ectopic expression of PTPRD at mRNA and protein level was validated by TaqMan qPCR and by Western blot analysis (for a V5 epitope tag) (Additional File 1). Ectopic up-regulation of PTPRD resulted in a significant decline in cell growth in all three cell lines (Figure 2A-C). To determine if the phosphatase activity of PTPRD is required to have anti-proliferative effects on cell lines, a mutant form of PTPRD containing a cancer specific mutation, Q1481X, that results in a truncated protein product lacking a functional C-terminal phosphatase domain, was transfected into Kelly cells. Over-expression of this PTPRD mutant did not lead to a significant decrease in cell growth compared to empty vector transfected cells ( Figure 2D), indicating that a functional phosphatase domain is required. FACS analysis of Kelly cells transfected with either wild type PTPRD or empty vector, followed by staining with propidium iodide and staining for Annexin V, revealed that PTPRD up-regulation resulted in a median 1.7 fold increase in apoptosis relative to empty vector control ( Figure 2E, F).

Identification of proteins that interact with PTPRD
Taken together, all available evidence indicates that PTPRD is a neuroblastoma tumor suppressor, but the molecular mechanism(s) is unknown, as only a small number of PTPRD interacting partners have been identified. PTPRD is known to interact with the metastasis suppressor protein 1 (MTSS1), liprin-alpha 1 (LAR), and signal transducer and activator of transcription 3 (STAT3) [6,16,17]. In order to discover new potential substrates for PTPRD in the context of neuroblastoma, we used a commercially available protein-protein interaction array containing 9,400 proteins. A biotinylated PTPRD protein was used to probe the protein microarray to detect potential interaction partners. This experiment was conducted in duplicate (biological replicates as two different biotinylated PTPRD baits were used). The consensus interacting proteins with above-background signal are listed in Table 1.
One of the top scoring proteins on the PTPRD interaction arrays was aurora kinase A (AURKA), a serine/ threonine kinase that is involved in microtubule formation and stabilization at the spindle pole during chromosome segregation [18]. In neuroblastoma, high AURKA expression is associated with poor patient survival [15], so that the interaction of a tumor suppressor with AURKA was of immediate interest. In order to confirm an in vivo interaction between PTPRD and AURKA in neuroblastoma, protein extract isolated from Kelly cells transfected with either the V5-tagged PTPRD cDNA or empty vector was immunoprecipitated with a V5 epitope antibody. As illustrated in Figure 3A, western blotting demonstrated that V5-tagged PTPRD was immunoprecipitated from cells, with a clear band of approximately 85 kDa being detected compared with empty vector control. The PTPRD pro-protein (220 kDa) is cleaved intracellularly into two subunits of 150 kDa (E-subunit) and 85 kDa (P-subunit) [19]. This process is called "ectodomain shedding" and affects a number of transmembrane proteins. The 150-kDa E-subunit is predicted to contain the three Ig-like and eight FN-III-like domains, whereas the 85-kDa P-subunit contains a short ectodomain segment, the transmembrane peptide, and the intracellular PTPase domains and is expressed in the membrane of the cell. The 150 kDa band is not visualized because it has lost the V5 epitope tag. Western blot probed with AURKA antibody indicates that AURKA co-immunoprecipitates with the V5 epitope pulled down with the V5 antibody, confirming a PTPRD-AURKA interaction in vivo ( Figure 3B).

PTPRD destabilizes AURKA by dephosphorylating tyrosine sites
Having confirmed a PTPRD-AURKA interaction, we next investigated whether AURKA is a potential substrate for PTPRD. Within the catalytic domain of AURKA, there are a number of conserved tyrosine phosphorylation sites (Y148, Y199, Y197, Y212) [20]). Immunoprecipitation of tyrosine phosphorylated proteins with a pan-phosphotyrosine agarose conjugate and subsequent western blotting with an AURKA antibody revealed an abolishment of AURKA tyrosine phosphorylation in PTPRD cDNA transfected cells compared to phosphatase dead mutant or empty vector control transfectants, confirming the mechanism of action of PTPRD ( Figure 3C). In addition to deactivating AURKA, ectopic over-expression of PTPRD also leads to a decrease in protein stability of AURKA. After 4 hours of cycloheximide treatment, there is a significant difference in the decay of AURKA in the PTPRD transfected cells compared to empty vector (EV) or phosphatase dead mutant ( Figure 3D, E). qPCR analysis of AURKA mRNA at 48 hours post PTPRD transfection indicated no significant change relative to negative control, confirming that the AURKA decrease is post-transcriptional (Additional File 2A). A critical function of AURKA in neuroblastoma is to regulate MYCN protein stability by sequestering MYCN from ubiquitin-mediated proteolytic degradation mediated by the FBXW7 ubiquitin ligase [21], leading us to hypothesize that a down-stream effect of ectopic over-expression PTPRD in neuroblastoma cells would be the destabilization of MYCN. As illustrated in Figure  3F, MYCN protein levels are also decreasing as a consequence of PTPRD ectopic over-expression. The decrease in MYCN protein is due to a post-transcriptional mechanism, as there is no change in MYCN mRNA (Additional File 2B).

Discussion
In this report we demonstrate a novel in vivo interaction between the tumor suppressor PTPRD and the oncogenic AURKA protein in neuroblastoma. Furthermore, to our best knowledge, this is the first report to demonstrate dephosphorylation of tyrosine residues in AURKA leading to its subsequent destabilization. The identification of AURKA as a novel target of PTPRD mediated dephosphorylation adds to a growing list of substrates for this tumor suppressor. Activity of AURKA is also regulated by phosphorylation of Thr288 in a cell cycle dependent manner [22]. Constitutive phosphorylation on Ser51 of AURKA has also been shown to be associated with its over-expression and protein stabilization, preventing protein degradation mediated by the ubiquitin-proteasome system [23] and it is likely that constitutive phosphorylation of tyrosine residues in AURKA also prevent its degradation. AURKA destabilization also has the down-stream effect of destabilizing MYCN, consistent with prior findings by Otto et al. [21]. Thus, PTPRD is likely to have a tumor suppressor function in any type of cancer with MYCN amplification, including other pediatric cancers such as medulloblastoma and rhabdomyosarcoma, or in adult tumors that are dependent on high AURKA levels.
MYCN directly up-regulates AURKA by binding to its promoter region, forming an auto regulatory feedback loop [21]. As previously mentioned, PTPRD expression is significantly lower in MYCN amplified tumors relative to non-MYCN amplified tumors [10], indicating that MYCN might be directly repressing PTPRD in order to maintain AURKA levels. However, analysis of a previously published data set on MYCN binding sites in neuroblastoma cell lines indicated that MYCN does not bind to the PTPRD promoter region [24,25], revealing that down-regulation of PTPRD in MYCN amplified tumors is an indirect effect.
Along with AURKA, a number of other proteins scored highly in our protein-protein interaction arrays, including SLAIN2, and two serine/threonine kinases, NEK2 and PLK1. SLAIN2 has recently been shown to play an important role in the microtubule complex by controlling microtubule dynamics and organization [26]. NEK2 is an integral component of the mitotic spindle-assembly checkpoint which is necessary for proper chromosome segregation during metaphase-anaphase transition [27], while PLK1 is an early trigger for G2/M transition [28]. Intriguingly, one of the primary functions of AURKA is the activation of PLK1 by direct phosphorylation of Thr210 [29], and high expression of PLK1 is also significantly associated with high-risk neuroblastoma and unfavourable patient outcome [30]. Inhibition of PLK1 protein has an anti-proliferative effect on neuroblastoma cell proliferation [30]. Although the interaction of PTPRD with these proteins requires validation, it seems likely that over-expression of PTPRD should have a pleiotrophic tumor suppressor effect by associating with, and dephosphorylating a number of interacting partners. The recent development of AURKA inhibitors, along with impressive results in pre-clinical models of pediatric cancer, indicate that AURKA could be a therapeutic target of great value in the treatment of neuroblastoma [31,32]. PTPRD, as a natural antagonist of AURKA, might be of therapeutic value in cancers where its inactivation is epigentically reversible, such as in instances of promoter region hypermethylation.

Conclusions
Over-expression of PTPRD in neuroblastoma cell lines results in decreased cell viability through the activation of apoptosis. We further demonstrate a novel PTPRD-AURKA protein interaction, and that PTPRD overexpression results in the dephosphorylation and the destabilization of AURKA, along with a down-stream destabilization of MYCN protein. These findings, taken together, indicate a novel mechanism of action for PTPRD as a neuroblastoma tumor suppressor through the inhibition of these important oncogene products.

Methods
Cell Culture, In Vitro Growth and Apoptosis Assays SHSY-5Y and CHP-212 neuroblastoma cell lines were obtained from the American Type Culture Collection, while Kelly was obtained from European Collection of Animal Cell Cultures. All cell lines were validated for the presence of previously published genomic imbalances using aCGH. Culture media was supplemented with 10% FBS.
Cells were transfected with 2 μg PTPRD cDNA, empty vector (pcDNA3.1-V5) or phosphatase dead mutant PTPRD (Q14811) using lipofectamine 2000 reagent according to the manufacturer's instructions (Invitrogen, Carlsbad, CA) and were analysed by a Cyquant assay (Invitrogen) at 24, 48 and 72 h. PTPRD induced apoptosis was demonstrated by FACS analysis of annexin-V staining. Kelly cells were transfected with empty vector or PTPRD cDNA. Cells were harvested 72 h post transfection and apoptosis was quantified using the FITC Annexin-V Apoptosis Detection Kit I (BD Pharmingen, San Diego, CA, USA). Experiments were performed in triplicate to quantify apoptosis by phosphatidylserine (PS) externalization.

Immunoprecipitation experiments
Kelly cells (80-90% confluence) were transfected with 2 μg PTPRD cDNA or empty vector (pcDNA3.1-V5) in 100 mm dishes. After 24 h transfection, cells were incubated in a RIPA buffer (Sigma-Aldrich, Arklow) containing complete protease inhibitor mixture. The lysates were cleared by centrifugation for 20 min at 14,000 × g and then subjected to immunoprecipitation with either anti-V5 agarose affinity beads (Sigma-Aldrich) for co-IP experiments or 4 G10 ® Platinum, anti-phosphotyrosine agarose conjugate (Millipore) for pan-tyrosine IP experiments, and then incubated at 4 C for 2 h. The beads were washed five times with RIPA buffer, incubated in laemmli sample buffer and heated to 95 C for five mins. Resultant samples were separated by SDS-PAGE and were analysed by western blotting.

Western blotting
Cells were washed with PBS and lysed at 4 C for 20 min using RIPA buffer containg phosphatase and protease inhibitors. Lysates were centrifuged for 20 min at 14,000 × g and resultant protein extracts were separated by SDS-PAGE and blotted onto PVDF membranes. Antibodies used were as follows: anti-V5 (Invitrogen), anti-AURKA (Sigma), anti-MYCN (Santa Cruz) and antialpha tubulin (Abcam). Secondary antibodies were horseradish peroxidise-conjugated rabbit anti-mouse (Abcam) or goat anti-rabbit (Cell Signaling Technology). Reactive proteins were visualized by chemiluminescence with ECL (Thermo).

PTPRD Interaction arrays
Purified PTPRD protein (99% pure, mass spectrum verified) (Clone ID 8613c5BCt10p1, protein sequence 1201-1495) was biotinylated and quantitated in duplicate using the Invitrogen FluoReporter ® Biotin Quantitation Assay Kit. The array control protein (biotinylated calmodulin kinase) supplied as part of the Invitrogen protein interaction kit was used to probe the Protoarray ® control protein microarray to validate assay conditions. Biotinylated PTPRD baits (two biological replicates) were then used to probe the Invitrogen ProtoArray ® Protein Microarray Kit for biotinylated proteins Version 4 in order to detect potential interaction partners. The results were analyzed using the Protoarray prospector ® package. The consensus hits (hits that are above the back ground signal in both the experiments) along with their Z-scores are listed in Table 1. A Z-Score indicates how far and in what direction the value of an individual data point in a population falls from the mean in units of standard deviations. ProtoArray ® Prospector collects all the signals from the proteins on a microarray, calculates the mean value and standard deviation, and then calculates a Z-Score corresponding to the signal from each feature.

Cycloheximide half-life experiments
Kelly cells (80-90% confluence) were transfected with 2 μg PTPRD cDNA, empty vector (pcDNA3.1-V5) or PTPRD mutant (Q1481X) in 100 mm dishes. After 72 h transfection, cells were incubated in a RIPA buffer (Sigma-Aldrich, Arklow) containing complete protease inhibitor mixture. Cells were treated with 1 μl of a 20 mg/ml stock solution of cycloheximide (Sigma-Aldrich) 5 min prior to starting the indicated time course, and cells were collected at the indicated points. Quantification of western blots was performed by densitometry using GelEval (FrogDance Software). Normalization for loading differences was achieved by dividing the densitometry values for individual bands by the densitometry values for alpha-tubulin in the same lane.

Additional material
Additional File 1: Expression of PTPRD (A) mRNA and (B) protein following transfection of Kelly cells with increasing concentrations of PTPRD cDNA. mRNA was extracted at 24 hours and qPCR was performed. Lysates were harvested at 48 h and subjected to SDS PAGE and western blot analysis with a monoclonal antibody to the V5 epitope tag or alpha tubulin. All experiments were performed in triplicate.