- Research article
- Open Access
- Open Peer Review
Tartrate-resistant acid phosphatase (TRAP/ACP5) promotes metastasis-related properties via TGFβ2/TβR and CD44 in MDA-MB-231 breast cancer cells
BMC Cancervolume 17, Article number: 650 (2017)
Tartrate-resistant acid phosphatase (TRAP/ACP5), a metalloenzyme that is characteristic for its expression in activated osteoclasts and in macrophages, has recently gained considerable focus as a driver of metastasis and was associated with clinically relevant parameters of cancer progression and cancer aggressiveness.
MDA-MB-231 breast cancer cells with different TRAP expression levels (overexpression and knockdown) were generated and characterized for protein expression and activity levels. Functional cell experiments, such as proliferation, migration and invasion assays were performed as well as global phosphoproteomic and proteomic analysis was conducted to connect molecular perturbations to the phenotypic changes.
We identified an association between metastasis-related properties of TRAP-overexpressing MDA-MB-231 breast cancer cells and a TRAP-dependent regulation of Transforming growth factor (TGFβ) pathway proteins and Cluster of differentiation 44 (CD44). Overexpression of TRAP increased anchorage-independent and anchorage-dependent cell growth and proliferation, induced a more elongated cellular morphology and promoted cell migration and invasion. Migration was increased in the presence of the extracellular matrix (ECM) proteins osteopontin and fibronectin and the basement membrane proteins collagen IV and laminin I. TRAP-induced properties were reverted upon shRNA-mediated knockdown of TRAP or treatment with the small molecule TRAP inhibitor 5-PNA. Global phosphoproteomics and proteomics analyses identified possible substrates of TRAP phosphatase activity or signaling intermediates and outlined a TRAP-dependent regulation of proteins involved in cell adhesion and ECM organization. Upregulation of TGFβ isoform 2 (TGFβ2), TGFβ receptor type 1 (TβR1) and Mothers against decapentaplegic homolog 2 (SMAD2), as well as increased intracellular phosphorylation of CD44 were identified upon TRAP perturbation. Functional antibody-mediated blocking and chemical inhibition demonstrated that TRAP-dependent migration and proliferation is regulated via TGFβ2/TβR, whereas proliferation beyond basal levels is regulated through CD44.
Altogether, TRAP promotes metastasis-related cell properties in MDA-MB-231 breast cancer cells via TGFβ2/TβR and CD44, thereby identifying a potential signaling mechanism associated to TRAP action in breast cancer cells.
Tartrate-resistant acid phosphatase (TRAP/ACP5) is a metalloenzyme of the category of acid phosphatases  that is synthesized as a monomeric proenzyme (TRAP 5a, 35 kDa) . A disulfide linked heterodimer (TRAP 5b) with an N-terminal fragment of 20–23 kDa joined to the 16 to 17 kDa C-terminal part originates from post-translational cleavage of the monomeric form, which significantly increases phosphatase activity . TRAP derived from different mammalian sources reveals almost identical homology at the amino acid sequence level and identical biochemical properties [4,5,6,7]. During bone resorption, TRAP is secreted into the resorption lacuna of active osteoclasts, where it dephosphorylates the bone matrix protein osteopontin (OPN), thereby promoting osteoclast detachment and migration . Additionally, TRAP has been suggested to regulate OPN bioactivity in autoimmune conditions [9,10,11]. The isoform TRAP 5b was proposed as a serum marker for bone metastases in various types of primary cancers [12,13,14,15,16]. Interestingly, TRAP has also been detected in several cancer cells and tissues (breast, ovarian, cervical cancer and malignant melanoma) and its expression level correlates with the severity of the tumor [17,18,19]. Moreover, high TRAP expression correlates with reduced tumor- and metastasis-free survival in malignant melanoma , and with decreased overall survival and increased incidence of metastasis in hepatocellular cancer . In gastric cancer, elevated TRAP expression is an independent risk factor for peritoneal dissemination and is associated with shorter patient survival . In lung cancer, patients with high TRAP expression had a significantly lower overall survival than the patients with low TRAP expression .
Altogether, previous studies underscore the potential clinical relevance of TRAP to monitor cancer development and progression; nevertheless, underlying cellular and molecular processes remain unclear.
TRAP was shown to interact intracellularly with the Transforming growth factor β (TGFβ) receptor interacting protein-1 (TRIP-1), thereby activating TGFβ receptor type II (TβR2) and osteoblast differentiation through the Mothers against decapentaplegic homolog 2/3 (SMAD2/3) pathway at sites of prior bone resorption . Furthermore, TRIP-1 knock-down abrogates osteoblast differentiation and proliferation . TRAP 5a interaction with TRIP-1 has also been demonstrated in mouse pre-adipocytes .
TGFβ ligands exist in three highly homologous isoforms, TGFβ1, TGFβ2 and TGFβ3 and are part of a large family of structurally related secreted cytokines [27,28,29]. Upon ligand binding to the constitutively active serine/threonine kinase TβR2 , the latter forms a hetero-oligomeric complex with the type I receptor (TβR1). TβR2 then trans-phosphorylates TβR1, ultimately leading to the transcription of various target genes via both SMAD- and non-SMAD mediated pathways (reviewed in [29, 31,32,33,34]). TGFβ acts as tumor repressor in early stages of tumorigenesis and as an oncogene in late stages . For instance, in breast cancer patients expression of TGFβ was increased in tumor tissue [35, 36] and was associated with disease progression [36, 37]. TGFβ2 has been proposed as a predictive marker for breast cancer, as high levels of TGFβ2 correlate with advanced tumor stage and shortened survival . Additionally, TGFβ2 was reported as a catalyzer of TGFβ signaling through an autocrine loop . Finally, TβR1 contain a single cytoplasmic binding site for Cluster of differentiation 44 (CD44) , suggesting a potential interaction between the TGFβ pathway and CD44, a cancer-associated glycoprotein previously reported as an OPN receptor . Activated CD44 stimulates the serine/threonine kinase activity of TβR1, which in turns increases SMAD2/3 phosphorylation .
Aim of this study was to delineate by proof-of-concept, how TRAP promotes cellular properties related to metastasis in breast cancer cells at advanced state. As there is only limited knowledge about the molecular perturbations and possible substrates of TRAP, global phosphoproteomic and proteomic analysis was applied to connect possible signaling mechanisms to the TRAP-dependent phenotypic changes.
Matrix proteins: bovine milk OPN, previously purified by our group [42, 43], murine Engelbrecht-holm-swarm Laminin-1 (Lam I, Sigma-Aldrich, #L2020), human plasma Fibronectin (FN, Life Technologies, #PHE0023), human recombinant natural Vitronectin (VN, Life Technologies, #PHE0011), Cultrex® Rat Collagen I, (Col I, Trevigen, Cat# 3440–100-01), human placenta Collagen IV (Col IV, Merck Millipore, Cat# CC076) (10 μg/ml).
Antibodies: rabbit antibody serum against total TRAP (raised by immunization of New Zealand rabbits , 1:1000, TRAP 5a: 37 kDa, TRAP 5b: 16 and 25 kDa), mouse anti-β-Actin (1:1000, 42 kDa; Cat# 8224, Abcam), rabbit anti-TGFβ2 (Western blotting 1:1000, 48 kDa; ICC 1:200; Cat# 113670, Abcam), secondary donkey anti-mouse (Licor IRDye® 800CW, Cat# 925–32,212; 1:15,000); donkey anti-rabbit (Licor IRDye® 680RD, Cat# 926–68,073; 1:15,000); goat anti-rabbit (Licor IRDye® 680RD Cat# 926–68,071; 1:15,000); secondary goat anti-rabbit Alexa 488 (1:100, Cat#A11008, Life Technologies).
Recombinant protein, blocking antibodies and inhibitory compounds: rat monoclonal IgG2A Isotype Control (respective similar concentration, Novus Biologicals, MAB006, 54,447), rat monoclonal anti-CD44 (10 μg/mL, Novus Biologicals, Hermes-1, NBP2–22530), rabbit polyclonal anti-TGFβ2 (0,25 μg/mL, R&D systems, AB-12-NA); TGFβ receptor type I/type II kinases inhibitor LY2109761 (2 μM, Santa Cruz; CAS 700874–71-1, PubChem CID 11655119), Human recombinant TGFβ1 (10 ng/mL, R&D systems, 240-B), TRAP inhibitor 5-PNA (200 μM, 5-phenylnicotinic acid, Maybridge code CC24201, Sigma-Aldrich code CDS013984; PubChem CID 346160, previously characterized by our group ).
Cell transfection and culture
MDA-MB-231 breast cancer cell line was obtained from American Type Culture Collection (Manassas, U.S., ATCC® Number: HTB-26™). MDA-MB-231 were previously stably transfected with the eukaryotic expression vector pcI-neo containing the full size rat TRAP  and different subpopulations maintained in complete medium (RPMI 1640, 10% fetal bovine serum, 0.1 mg/mL Gentamicin) (Life technologies, Carlsbad, CA, U.S.) at 37 °C in a 5% CO2 humidified atmosphere. Cells were continuously tested for contamination with the MycoAlert™ mycoplasma detection kit (Lonza, Cat# LT07).
Knockdown of rat TRAP was achieved by the use of different custom cloned MISSION shRNA constructs within the lentivirus plasmid vector pLKO.1-puro containing ampicillin and puromycin antibiotic resistance genes (Sigma Aldrich, St. Louis, MO, U.S.). Transfections were done with Escort II transfection reagent according to the manufacturer’s instructions. 1 μg of purified DNA was complexed with 5 μl transfection reagent and applied to the cells for 24 h. Cells transfected with plasmids encoding antibiotic resistance were selected by culture of complete medium supplemented with 1 μg/mL puromycin (Sigma Aldrich, St. Louis, MO, U.S.). TRAP gene expression in the cells bearing a knockdown was quantified according to a detailed description in the Additional file 1: Material and Methods.
Tartrate-resistant acid phosphatase activity assay
Cell lysates and conditioned media were prepared according to a detailed description in the Additional file 1: Material and Methods. TRAP activity in cell lysate and medium was measured under optimal enzyme conditions using 10 mM p-nitrophenyl-phosphate as a synthetic substrate (Sigma Aldrich) as previously described . Enzymatic activity was calculated in enzyme Units according to Lambert-Beer and normalized to total protein concentration for lysate or per 106 cells per 24 h for medium.
Cell lysates and conditioned media were prepared according to a detailed description in the Additional file 1: Material and Methods. 25–50 μg of total protein and normalized volumes of corresponding medium were subjected to SDS-PAGE (Mini-PROTEAN®TGX™ precast gel, Biorad) and transferred to PVDF membrane (Trans-Blot turbo mini PVDF packs, Biorad) according to the manufacturer’s instructions. Unspecific binding was blocked by 3% bovine serum albumin (BSA) in PBS for 1 h at room temperature (R.T.) and protein bands detected by subsequent incubation with respective primary antibodies in each 3% BSA overnight (o.n.) at 4 °C and fluorescently labelled secondary antibodies for 1 h at R.T.. The membranes were washed with TBST (20 mM Tris- HCl pH 7.5, 500 mM NaCl, 0.05% Tween-20) after antibody incubations and visualized in the Licor Odyssey Fc Imager and quantified by densitometry with the Licor Image Studio software 3.1.4 (Licor Biosciences, Lincoln, NE, U.S.) upon normalization to β-Actin expression for lysates.
6-well plates were precoated with a bottom layer of 0.8% low melting Agarose solution (Sigma Aldrich, A9414) prepared in complete medium and gel formation allowed for at least 1 h at R.T.. After gelling, 10,000 cells were incorporated in a 0.35% Agarose solution and overlaid on the previous high Agarose-containing lower layer. During regular feeding, cells were allowed to form colonies for 3 weeks and agarose gels fixed in formaldehyde (Solvecco). Colonies were stained in 0.01% Crystal violet solution in formaldehyde o.n. and 1X images taken. Images were changed to binaries and threshold adjusted on all images using Image J 1.48 software. Colony number and sizes were quantified by the “count particles” function with a threshold of 50 pixels as minimum detection limit and inclusion of colonies with circularity from 0 to 1.
Cell growth assay
Fifty thousand cells were seeded in 48-well plates in complete medium and let adhere and grow for 24 h or 48 h. The cells were washed with PBS, fixed in formaldehyde and stained for 5 min in 0.1% Toluidine blue. Excess color was washed away with PBS and color dissolved for 5 min in 50% EtOH/50 mM HCl. Absorbance was measured at 630 nm in a PowerWave HT Microplate Spectrophotometer (Biotek, Winooski, VT, U.S.).
Fifty thousand cells were seeded into glass 8-well chambers (Labtek II, 154534) and let adhere o.n.. After a PBS wash, medium was changed, respectively, conditionally containing blocking antibodies or inhibiting compounds or deficient in serum and maintained for another 24 h. Following, the cells were pulsated with 10 μM EdU for 1 h and subsequently fixed in 4% formaldehyde for 10 min. Fluorescent staining was performed according to the manufacturer’s protocol with the Click-iT® Plus EdU Alexa Fluor 488 Imaging Kit (Cat# C10637, Invitrogen, Life Technologies Europe BV, Stockholm, Sweden). Confocal images were acquired in the Nikon A1+ confocal laser microscope system and image batch analysis performed in the Nis Elements Advance research imaging software 4.1.0 (Nikon). For automated cell counting, thresholding was adapted to experimental controls and nuclei with circularity greater than 0.1 and size (diameter) greater than 5 μm taken into account.
Bright field phase contrast images of cells cultured in complete medium were taken at equal confluence and at different time points in a Nikon Eclipse TE300 Inverted microscope equipped with a DS-Fi1 digital microscope camera and a DS-U2 camera control unit (Nikon, Tokyo, Japan). Ratios of cell length to cell width were acquired for all cells using the Straight line tool in Image J 1.48 program. A minimum of 200 cells was randomly measured per time point and ratio medians calculated.
Live cell wound migration
Wound migration experiments were performed and analyzed as previously described . Medium was respectively enriched with blocking antibodies, chemical inhibitors or compounds.
Transwell migration assay
Transwell assays membranes (Corning Incorporated Costar-Transwell CLS3422-48EA) were precoated on the lower side of the insert with 10 μg/mL of matrix proteins in PBS under humidified atmosphere at 4 °C for 24 h. Following, 200,000 cells were seeded into the upper chamber of the Transwell and both wells filled with serum-free medium. Migration was allowed for two different time periods during linear increase (OPN: 8 h, 22 h; FN: 2 h, 4 h; VN: 6 h, 8 h; Col I: 2 h, 4 h; Col IV: 4 h, 8 h; Lam I: 2 h, 4 h). Cells that had migrated on the lower membrane surface were fixed in 4% formaldehyde and stained in 0.1% Crystal violet solution (200 mM Borat, pH 9) for 10 min. Cells left on the upper side of the Transwell membrane were removed and excess color washed away with water. Images were taken at 20 X magnification. Color was dissolved in 10% acetic acid for 5 min and read at 600 nm.
Transwell invasion assay
CytoSelect™ 96-Well Cell Invasion Assays (Cellbiolabs, Cat#CBA-112, San Diego, CA, U.S.;) were performed according to the manufacturer’s protocol as previously described .
Twenty thousand cells were allowed to grow for 48 h in complete medium in an 8-well chambered slide (Labtek II, 154534). Cells were washed in serum-free medium, fixed with 4% formaldehyde (Solvecco) and washed in TBST (25 mM Tris pH 7.4, 150 mM NaCl, 0.1% Tween 20) followed by permeabilization in 0.1% Triton X-100 (Sigma-Aldrich) for each 10 mins at R.T.. Unspecific binding was blocked in 1% BSA/TBST (Sigma-Aldrich) at R.T. for 1 h followed by incubation with primary and secondary antibodies for each 1 h at R.T.. Cells were after antibody incubations washed twice with TBST and finally stained with Hoechst (1:7500 dilutions, Life technologies) for 3 mins. The cells were washed, mounted in fluorescent mounting media (DAKO) and imaged using a Nikon A1+ confocal laser microscope system equipped with Nis Elements Advance research imaging software 4.3.0 (Nikon, Sweden, Stockholm) at 60 X. For quantification fluorescent signal intensity in 100–400 cells per experiment was measured with respect to a threshold set on experimental controls.
Peptide level high-resolution isoelectric focusing (HiRIEF)
SILAC labeled peptide samples were prepared for quantitative phosphoproteomics and proteomics analysis of control cells and TRAP3high cells, while TMT labeled peptide samples were prepared for scrambled and knockdown cells (sh2 and sh3 + 4) according to a detailed description in the Additional file 1: Material and Methods.
HiRIEF was performed as described previously for analysis of both SILAC and TMT labeled samples . Immobilized pH gradient (IPG) gel strips (GE Healthcare Bio-Sciences AB, Uppsala, Sweden) with linear pH ranges of 2.5–3.7 (“ultra-acidic” range) and of 3–10 (“wide range”) were employed for standard phosphoproteomics and for proteomics analyses, respectively. Strips were divided into 72 fractions (fraction numbering proceeds from the acidic end towards the basic end of the strips), and extracted to V-bottom 96-well plates with a liquid handling robot (GE Healthcare prototype modified from Gilson liquid handler 215). Plates were lyophilized in a Speedvac prior to liquid chromatography-mass spectrometry (LC-MS) analysis.
LC-MS data statistical analyses
The LC-MS proteomics data have been deposited to the ProteomeXchange Consortium via the PRIDE partner repository with the dataset identifier PXD006430. (https://www.ebi.ac.uk/pride/archive/login). Details regarding the LC-MS analysis and the calculations of protein and phospho-site ratios are described in Additional file 1: Material and Methods.
To define significantly regulated events in both SILAC-based or TMT-based quantitative analysis, a median absolute deviation (MAD) value for all the ratios in each experimental condition was calculated as the median absolute value of the differences between each individual log2(ratio) and the median log2(ratio). All log2 transformed ratios were median absolute deviation (MAD) scaled by calculating robust z-scores (MADs away from the median), as described before [48, 49], to correct for their skewed distributions. Significantly regulated events were those with robust z-scores of at least −/+ 2.5 and −/+ 3 in both replicates for SILAC-based phosphoproteomics and standard proteomics analysis respectively. For TMT-based standard proteomics analysis, significantly regulated events were those with average robust z-scores of at least −/+ 3 and t-test p-value <0.01 for both sh2 and sh3 + 4 samples.
Plots were generated using RStudio. Gene assignment to different classes was based on information retrieved from publicly available databases for protein kinases , protein phosphatases , transcription factors [52,53,54] and enzymes of the ubiquitin and ubiquitin-like (UBL) conjugation systems . The list of human phosphorylation sites with a previously reported function was obtained from the PhosphositePlus database, released 2017–02-16 . Gene ontology (GO) enrichment analysis was performed with the web service GOrilla by selecting the “Two unranked lists of genes” option and setting a p-value threshold of 10−3. The target set included all the genes that were significantly regulated in either of the proteomics and phosphoproteomics analyses (842 genes), while the background set was the list of all the identified genes (9570 genes).
Network analysis of the phosphorylated proteins was performed using the Cytoscape software platform. Significantly regulated genes from the GO analysis belonging to biological adhesion or ECM organization processes were employed to generate the network. Direct interactions between those genes were extracted from the STRING database (version 10.0, 2016–04-16)  and visualized with Cytoscape.
Results visualized in columnar graphs were expressed as mean values ± standard deviation. Boxplots display the 25th and 75th percentile with medians as vertical line. Whiskers represent the minimum and maximum. Analysis and statistical comparison of the mean of at least three biological replicates was compared with GraphPad Prism 6 Software and each experiment based on one to five technical replicates. The respective amount of biological experiments was denoted more specifically in each figure legend. If not denoted specifically, TRAP-overexpressing cells were compared to control cells and TRAP knockdown cells compared to scrambled cells, respectively. In case of normal distribution, statistics was done by parametric two sample t-test (2 groups) or ANOVA test (> 2 groups). Otherwise, non-parametric Mann- Whitney test (2 groups) or Kruskal- Wallis one- way ANOVA on ranks (> 2 groups) was applied. A p-value <0.05 was considered significant (*); values p < 0.01 (**), values p < 0.001 (***), values p < 0.0001 (****) were marked in the graphs, respectively.
See Additional file 1: Material and Methods for description of gene expression analysis, cell lysis, SILAC labeling, protein extraction for mass spectrometric analyses, protein digestion, Tandem Mass Tag labeling, liquid chromatography tandem mass spectrometric analyses and proteomics database search, protein and phospho-peptide ratios calculation.
Generation and characterization of TRAP-overexpressing and TRAP knockdown MDA-MB-231 cells
MDA-MB-231 cells, an invasive breast cancer cell line that has previously not been tested with regard to the role of TRAP was employed to test whether TRAP overexpression could further enhance invasive capabilities of these cells. As this cell line expresses low levels of TRAP, processing within the different cellular compartments, secretion, as well as signaling of TRAP in the overexpression system should be similar to the parental one. Four clonal populations of cells stably transfected with a full-length rat TRAP and a control cell population transfected with a mock insert (control) , were characterized by comparing TRAP protein expression level and enzymatic activity in cell lysates and media (Fig. 1a-c). As previously reported, in TRAP-overexpressing cells  both isoforms TRAP 5a and TRAP 5b are present in cell lysates, whereas only the monomeric isoform 5a is present in the media (Fig. 1a). TRAP-overexpressing cell populations possess intermediate (TRAP1low and TRAP2low) and high (TRAP3high and TRAP4high) TRAP protein levels and enzymatic activity in lysate and medium as compared to control cells (Fig. 1b, c).
Additionally, TRAP3high cells were transfected with two different shRNA sequences targeting rat TRAP (sh2 and sh3 + 4) to generate TRAP knockdown cells. TRAP3high cells transfected with a scrambled shRNA sequence (scr) were used as control (Fig. 1d-g). TRAP mRNA (Fig. 1d) and protein expression (Fig. 1e, f) were decreased by around 80% in TRAP knockdown cell lysates and media. TRAP enzymatic activity was decreased in both scrambled cells and TRAP knockdown cells compared to TRAP3high cells, however significantly lower in TRAP knockdown cells compared to scrambled cells (Fig. 1g).
Overexpression of TRAP increases cell growth and proliferation
Anchorage-independent growth of TRAP3high cells was assayed by their ability to form colonies in soft agar (Fig. 2a, one representative well per condition). Numbers of colonies were increased by 58% in TRAP3high cells compared to control cells (Fig. 2b) and significantly larger than the ones formed by control cells (Fig. 2c). Anchorage-dependent growth of TRAP3high and TRAP4high cells compared to control cells were increased by 44% and 50%, respectively, after 48 h of culture in serum-supplemented medium (Fig. 2d). Proliferation, determined by the incorporation of the thymidine analogue EdU during S-phase, was increased from 31% in the control cells to 37%–43% in the TRAP-overexpressing cells (Fig. 2e). Accordingly, TRAP knockdown cells (sh3 + 4) displayed significantly lower proliferation than scrambled cells (Fig. 2f). Additionally, inhibition of TRAP by the small molecule inhibitor 5-PNA normalized proliferation of the TRAP3high cells down to the level of the control cells, without affecting proliferation of the control cells (Fig. 2g). Interestingly, TRAP3high cells also displayed an increased capacity to proliferate in serum-free medium after 48 h (Fig. 2h), where around 9.5% of the TRAP3high cells remained cycling compared to 7% of the control cells.
TRAP enhances an elongated morphology and metastasis-related hallmarks
TRAP affects cellular morphology as cells with increasing TRAP protein levels showed increasingly elongated phenotypes. Cellular length-to-width ratios were increased in the TRAP-overexpressing (TRAP1low, TRAP3high, TRAP4high) cells compared to control cells. Additionally, TRAP3high cells were more elongated than the TRAP1low cells (Fig. 3a). This morphological phenotype was reverted upon TRAP downregulation, as shown by decreased length-to-width ratios in sh2 and sh3 + 4 cells compared to scrambled cells (Fig. 3b). Furthermore, transwell migration was higher in TRAP3high cells than in control cells when the wells were coated with fibronectin (1.4-fold), collagen IV (1.6-fold) and laminin I (1.4-fold) (Fig. 3c). Interestingly, a very prominent increase in transwell migration was observed in the presence of the phosphorylated extracellular matrix (ECM) protein osteopontin (OPN) (2.6-fold), a well-known substrate for TRAP . No difference in transwell migration was detected in the presence of vitronectin or collagen I. Furthermore, wound migration was assessed by live cell imaging over 30 h. TRAPhigh cells but not TRAPlow cells displayed significantly increased migration velocity compared to the control cells; additionally migration velocity was significantly increased in TRAP3high cells compared to TRAP1low cells, indicating a level-dependent regulation of migration by TRAP (Fig. 3d). Accordingly, migration velocity was decreased in sh2 and sh3 + 4 cells compared to scrambled cells (Fig. 3e). Finally, TRAP also promoted transwell invasion through a basement membrane layer, as indicated by a 2.3–3.3-fold increase in invasion in TRAP-overexpressing cells as compared to control cells (Fig. 3f). Importantly, treatment with the TRAP inhibitor 5-PNA was recently shown to revert TRAP-dependent promotion of cell wound migration and transwell invasion .
Proteomics and phosphoproteomics profiling of TRAP perturbed cells
In order to investigate the molecular changes underlying the phenotype of TRAP-overexpressing cells, a large-scale quantitative proteomics and phosphoproteomics analyses using high-resolution isoelectric focusing (HiRIEF) fractionation  coupled to LC-MS was performed. Analyses were conducted in biological duplicates and relative quantifications are expressed as ratios of TRAP3high (Heavy SILAC labeled) relative to control (Light SILAC labeled) samples.
Since TRAP possesses phosphatase activity, phosphoproteomics analysis upon TRAP overexpression was employed to identify putative TRAP targets. For phosphoproteomics analysis, HiRIEF fractionation using “ultra-acidic” pH range (2.5–3.7) IPG strips was employed to enrich for phosphorylated peptides, as the addition of phosphate groups decreases the peptide pI. This approach led to the identification of 6266 and 5771 unique peptides in replicate 1 and replicate 2, respectively, of which 3325 and 3892 were phosphorylated (Table 1). Across both replicates, 3290 unique phosphorylations sites corresponding to 1059 genes were identified (Fig. 4a).
Notably, these results demonstrate that HiRIEF fractionation alone can be used to perform phosphoproteomics profiling with moderate analytical depth, independently than other enrichment methods, as previously reported . This approach is particularly suitable to identify multiply phosphorylated peptides, as about 50% of the identified phospho-peptides carry two or more phosphorylations. Multiply phosphorylated peptides are enriched in the first 30 fractions (acidic end) of the IPG 2.5–3.7 strip, while singly phosphorylated peptides are identified mostly in the more basic strip fractions (Fig. 4b; Additional file 2: Figure S1).
Phospho-site ratios (normalized to the total protein levels) and protein ratios have Pearson coefficients of correlation between replicates of 0.60 and 0.84 respectively (Additional file 3: Figure S2). One hundred ninety-three and 353 proteins and 116 and 58 phosphorylation sites were significantly down- and upregulated in both TRAP3high replicates (Fig. 4c, d; Additional file 4: Table S1). Among the top upregulated phosphorylation sites, S179, S180 and S183 of CD44 displayed a much higher magnitude of regulation than any other sites, with an increase of more than 11-fold in TRAP3high cells compared to control cells. These three sites located in the intracellular portion of CD44 do not have any previously reported function. The 116 phosphorylation sites downregulated in response to TRAP overexpression represent putative targets of TRAP phosphatase activity. Among the significantly downregulated phosphorylation sites, 11 sites were previously annotated to be functional, including S183 of protein phosphatase 1G (PPM1G); S131 and S137 of lysine-specific histone demethylase 1A (KDM1A); S957 and S966 of structural maintenance of chromosomes protein 1A (SMC1A); S27 of isoform 2 of X-ray repair cross-complementing protein 6 (XRCC6); and S453 and T455 of isoform 3 of Mediator of DNA damage checkpoint protein 1 (MDC1). All these sites are involved in promotion of DNA damage response. Finally, the list of genes corresponding to the significantly regulated phosphorylation sites is enriched in GO Molecular Function terms “cell adhesion molecule binding” and “cadherin binding”, as well as Cellular Component terms related to cellular junctions (Additional file 5: Figure S3).
Additionally quantitative proteomics analysis by tandem mass tags (TMT) of scrambled and TRAP knockdown (sh2 and sh3 + 4) cells was performed in three biological replicates each (Fig. 5a), identifying and quantifying 9848 proteins corresponding to 9189 genes. Quantifications for each TMT channel are expressed as ratios relative to the average of the three scrambled cell samples. Hierarchical clustering based on Euclidian distance and complete linkage separates samples based on experimental condition, illustrating good reproducibility of replicate measurements (Fig. 5b). 174 and 144 proteins were significantly down- and upregulated respectively, upon TRAP knockdown in both sh2 and sh3 + 4 samples (Additional file 6: Figure S4; marked with red row-side colors in Fig. 5b). TRAP protein expression level was significantly reduced by at least 1.5-fold in sh2 and sh3 + 4 cells.
Proteomics and phosphoproteomics analyses of TRAP-perturbed cells reveal regulation of cell adhesion and extracellular matrix organization network
Enrichment of gene ontology (GO) processes was evaluated comparing the list of all genes regulated upon TRAP perturbation in either of the phosphoproteomics or proteomics analyses (Additional file 7: Table S2) with the list of all identified genes. Significantly regulated genes included 21 protein kinases, 9 protein phosphatases, 25 transcription factors and 24 ubiquitin and ubiquitin-like (UBL)-conjugating system enzymes (Additional file 7: Table S2).
Biological adhesion and ECM organization processes showed a high degree of enrichment based on significance and fold enrichment (Fig. 6a). This result is in agreement with our previous observations of increased migration and invasion upon TRAP overexpression, as those processes are tightly interconnected. Additionally, we observed enrichment in processes related to mitochondrial translational termination; examination of the genes included in these processes shows that several mitochondrial ribosomal proteins were upregulated upon TRAP knockdown. As mitochondrial ribosomes translate exclusively mitochondrial encoded mRNA (13 genes, mainly electron transport chain proteins), this result might indicate increased mitochondrial protein synthesis in the mitochondria upon TRAP knockdown, possibly suggesting that TRAP affects cellular metabolism.
To examine the regulation of the genes specifically involved in modulating biological adhesion and ECM organization upon TRAP perturbation, we generated a protein-protein interaction network of TRAP regulated proteins classified into those GO terms, using the STRING database  (Fig. 6b). Proteins in the network involved with cell adhesion include protein kinases ephrin type-A receptor 2 (EPHA2), peripheral plasma membrane protein (CASK) and ROCK1 as well as receptor-type tyrosine-protein phosphatases mu, F and U (PTPRM, PTPRF and PTPRU). Interestingly, several proteins known to be involved in the migration and invasion processes were included in the network, such as matrix metalloproteases (MMPs), collagens (core components of the ECM), integrins, CD44, TGFβ2 and TβR1.
Based on previous publications reporting a connection of TRAP with TGFβ signaling through TRIP-1 and the significant upregulation of CD44 phosphorylation sites and the TGFβ pathway-associated proteins TGFβ2, TβR1 and SMAD2 measured in the SILAC proteomics analysis (Table 2), we hypothesized that an activation of TGFβ signaling or regulation of CD44 upon TRAP overexpression is responsible for the observed metastasis-related properties.
TRAP effects on proliferation and migration are mediated via TGFβ pathway-associated proteins
To consolidate our hypothesis, protein expression of the most upregulated TGFβ pathway protein, TGFβ2, was quantified in TRAP3high cells and compared to control cells (Fig. 7). Concordantly with the global analysis, TGFβ2 protein was upregulated in TRAP3high cells as shown by both immunocytochemistry (Fig. 7a, b) and Western blotting (Fig. 7c, d). Additionally, TGFβ2 expression was reduced to control cell level upon inhibition with the TRAP inhibitor 5-PNA (Fig. 7c, d).
To examine the functional impact of TGFβ2/TβR in TRAP3high cells, we assayed cell proliferation and migration upon treatment with blocking antibodies against TGFβ2 (Fig. 8a-c) or with a small molecule inhibiting TβR1/2 (LY2109761; Fig. 8d-f). Blocking of either TGFβ2 or TβR1/2 kinase activity reduced the proliferation of TRAP3high cells to the level of control cells (Fig. 8a, d). Furthermore, interfering with the TGFβ pathway-associated proteins decreased migration (Fig. 8b, c, e, f) in TRAP3high cells to the levels detected in the control cells. Migration curves over 30 h are shown in Fig. 8c and f. Control cell proliferation or migration were not affected when interfering with the TGFβ pathway associated proteins by using blocking antibodies against TGFβ2 or the small molecule inhibiting TβR1/2 activity (Fig. 8a-f). Treatment with the closely related ligand TGFβ1 did not affect proliferation and migration of control and TRAP3high cells (Additional file 8: Figure S5).
As discovered by phosphoproteomics analysis, phosphorylation of S179, S180 and S183 of CD44 was increased more than 11-fold in TRAP3high compared to control cells, representing the top upregulated phospho-events in TRAP3high cells. These three phosphorylation sites are located in the intracellular domain of CD44 and do not have any previously reported function. No antibodies targeting those phosphorylation sites are currently available, thus we probed the role of CD44 in determining the TRAP-dependent phenotype by inhibiting CD44 using a blocking antibody (Fig. 8g-i). Upon anti-CD44 treatment, the proliferation rate of TRAP3high cells was reduced to a level lower than that of control cells (Fig. 8g), while the TRAP-dependent increase in migration was not affected (Fig. 8h, i).
TRAP has been shown to be associated with tumor progression in several types of cancer and suggested to be clinically relevant as a marker for peritoneal dissemination in gastric cancer [20,21,22]. Recently, TRAP was identified as a pro-invasion oncogene and a prognostic marker in melanoma . Moreover, TRAP’s functional role in invasion, cell motility and metastasis was suggested to be mediated through phosphorylation of focal adhesion complexes .
In this study, we investigated the effects of TRAP expression on cell properties related to the development of metastasis and on the cellular signaling of the MDA-MB-231 breast cancer cell line at the system level. The small molecule 5-phenylnicotinic acid (5-PNA), recently identified as a specific inhibitor of TRAP activity and TRAP-dependent migration and invasion, was employed as a tool to investigate the molecular events mediated by TRAP [45, 59]. We demonstrate that TRAP regulates metastasis-related features such as anchorage-independent and –dependent cell growth, proliferation, migration and invasion. Global proteomics and phosphoproteomic analyses showed that regulated events upon TRAP perturbation are mainly proteins and phosphorylation sites involved with cellular adhesion and extracellular matrix (ECM) organization. Based on these analyses and on the literature, we hypothesized that TGFβ pathway-associated proteins and three previously unreported intracellular phosphorylation sites of CD44 mediate the observed TRAP-dependent cellular phenotypic properties.
TRAP has been proposed as a differentiation and growth factor for cells of hematopoietic origin . Effects of TRAP expression on cell transformation and tumor progression have been clinically validated in melanoma, as well as proven in vitro by an anchorage-independent growth assay . TRAP overexpression correlates with increased tumor size and poor differentiation in hepatocellular cancer . We consolidated and expanded these findings by using the breast cancer cell line MDA-MB-231, showing that TRAP overexpression increases cell growth as well as colony formation and cell proliferation above basal levels. Moreover, a higher number of TRAP-overexpressing cells compared to control cells were actively proliferating after 48 h serum starvation, indicating a lower requirement for exogenous growth stimulation.
Morphological changes such as cell rounding or cell spreading upon perturbation of TRAP expression have been reported in melanoma cells . Modulation of TRAP expression was shown to impact the migration and invasion of melanoma and hepatocellular carcinoma cells both in vitro and in vivo, when either non-invasive cancer cells, expressing low amounts of TRAP or metastatic high-TRAP expressing cells were subjected to upregulation or knockdown, respectively [20, 21]. Moreover, also in non-malignant epithelial cells TRAP expression was linked to a regulation of cell migration . This study demonstrates that TRAP overexpression enhances the elongated phenotype, migration and invasion capabilities of invasive breast cancer cells. Importantly, the elongated morphology and migration were regulated by TRAP in a dose-dependent manner. The presence of ECM proteins and basement membrane proteins Collagen IV and Laminin I increased transwell migration of TRAP3high cells as compared to control cells, underscoring the role of TRAP during the invasive process. Transwell migration was particularly increased in the presence of osteopontin (OPN), a highly phosphorylated ECM protein previously suggested to be a physiological substrate for TRAP , and involved in the progression of TRAP-related pathologies such as the immuno-osseous disorder Spondyloenchondrodysplasia [9, 62]. OPN has been reported as a ligand to the CD44 receptor  and was shown to increase osteoclast migration , which is blunted upon antibody-mediated blocking of CD44 . Inhibition of TRAP by the small molecule inhibitor 5-PNA was previously reported to decrease TRAP3high cells migration and invasion ; here we showed that also proliferation of TRAP3high cells is reduced to basal levels upon treatment with 5-PNA, altogether providing evidence that the above mentioned phenotypes of TRAP-overexpressing MDA-MB-231 cells are attributable to the overexpression of TRAP.
In parallel, global proteomics analysis of TRAP3high cells revealed regulation of various proteins belonging to the GO terms “biological adhesion” and “ECM organization”. Coherently, an increase in migration and invasion on various ECM and basement membrane proteins was observed in the TRAP3high cells. Enrichment in closely related GO terms, such as “cell adhesion molecule binding” and “cell junction”, was noted when analyzing phosphosites regulated in TRAP-overexpressing cells compared to control cells, further substantiating the involvement of TRAP in these functions. The list of 119 phosphorylation sites downregulated upon TRAP overexpression represent an inventory of putative targets of TRAP phosphatase activity or possible signaling intermediates; among those, eight sites with known regulatory function are involved in DNA damage response, another hallmark of cancer.
Most importantly, we identified a regulation of the TGFβ pathway-associated proteins TGFβ2, TβR1 and SMAD2, as well as a highly significant upregulation of previously unreported phosphorylation sites of CD44 upon TRAP perturbation in the MDA-MB-231 breast cancer cell line. Quantification of expression levels by several methodological approaches confirmed the upregulation of the ligand TGFβ2, which could be reverted by treatment with the TRAP inhibitor 5-PNA. Functional blocking of TGFβ2 or inhibition of TβR1/2 kinase activity restrained the increase in migration and proliferation promoted by TRAP. Antibody-mediated inhibition of CD44 reduced proliferation beyond the basal level of control cells.
Several reports support the concept that TGFβ promotes or restrains cell proliferation depending on the context . TGFβ inhibited the proliferation of most epithelial cells and its growth inhibitory effect could be partially reverted by treatment with a specific inhibitor [64,65,66]. Poorly differentiated prostate cancer cells were resistant to TGFβ growth inhibitory effect in vivo . Upon functional blocking of TGFβ2 or inhibition of TβR1/2 kinase activity, we observed decreased proliferation of TRAP-overexpressing cells but not of control MDA-MB-231 cells. Such lack of response in the control cells suggests that the malignant MDA-MB-231 cells are resistant to the growth-inhibitory effect of TGFβ, and that TRAP-dependent increase in proliferation of TRAP-overexpressing cells is TGFβ2-mediated. Furthermore, the invasive capacity of malignant breast cancer cells is enhanced by TGFβ1  and inactivation of TGFβ signaling inhibited invasiveness in vitro and in vivo for colon carcinoma cell lines . Additionally, here we detected a TRAP-dependent promotion of migration mediated by TGFβ2.
We could exclude a possible contribution of TGFβ1 because treatment with this TGFβ ligand modulated neither migration nor proliferation, as reported previously, despite constitutive expression of its receptors . Thereby we substantiated that the TGFβ2 isoform is crucial for the TRAP-mediated effects. In support to this notion, TGFβ2 was previously attributed a dominant role in predicting the outcome of breast cancers  and aberrant expression of the TGFβ2 isoform exclusively was induced through an autocrine loop in glioma .
We were also interested in the regulation of CD44, as it has previously been reported to be connected to TGFβ pathway [70,71,72] and to be phosphorylated by TRβ1 . Our analysis identified three phosphorylation sites of CD44 as the top upregulated phosphorylation events in TRAP3high cells. Additionally, CD44 is a receptor for OPN , which is the protein that prominently increased transwell migration in TRAP-overexpressing cells in this study. Functional blocking of CD44 resulted in a decrease of proliferation of TRAP3high cells beyond the basal level of control cells, but had no effect on cell migration. This suggests that cell migration is regulated by TRAP via TGFβ independently of CD44, but that the basal and TRAP-dependent proliferative activity of the TRAP3high cells is regulated through both TGFβ signaling and CD44.
As this study is performed in vitro and uses a metastatic cell line, similar experiments in other cancer cell lines at different stages of metastatic progression would be crucial to fully delineate the effects of TRAP, possibly revealing a role in promoting metastasis of non-invasive cell lines as well [20, 21]. Additionally, knockdown of TRAP in invasive cancer cells expressing high TRAP levels might be employed to test whether their invasive phenotype is TRAP-dependent. Limitations include also a full dissection of the respective TGFβ pathway and a reconfirmation, to allow for a generalization to other TRAP expressing cell lines.
This is the first study investigating the effects of TRAP perturbation in cancer cells on a global scale and identifying possible substrates and signaling intermediates. Herein, we show that TRAP promotes metastasis-related cellular properties, such as cancer cell proliferation, migration and invasion beyond basal levels in malignant breast cancer cells and regulates cell adhesion and extracellular matrix (ECM) organization. TRAP-dependent migration and proliferation can be neutralized upon inhibition of TGFβ2/TRβ and CD44 or inhibition of TRAP by the small molecule 5-PNA. Altogether, this data provide the basis for further studies investigating TRAP signaling and a novel possible targeting strategy for the treatment of tumors with high TRAP expression.
Cluster of differentiation 44
Transforming growth factor β
Tartrate-resistant acid phosphatase
Transforming growth factor receptor β
Oddie GW, Schenk G, Angel NZ, Walsh N, Guddat LW, de Jersey J, et al. Structure, function, and regulation of tartrate-resistant acid phosphatase. Bone [Internet]. 2000;27:575–84. Available from: http://www.ncbi.nlm.nih.gov/pubmed/11062342
Ljusberg J, Ek-Rylander B, Andersson G. Tartrate-resistant purple acid phosphatase is synthesized as a latent proenzyme and activated by cysteine proteinases. Biochem J [Internet]. 1999;343 Pt 1:63–9. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1220524&tool=pmcentrez&rendertype=abstract
Uppenberg J, Lindqvist F, Svensson C, Ek-Rylander B, Andersson G. Crystal structure of a mammalian purple acid phosphatase. J Mol Biol. 1999;290:201–11.
Ek-Rylander B, Bill P, Norgård M, Nilsson S, Andersson G. Cloning, sequence, and developmental expression of a type 5, tartrate-resistant, acid phosphatase of rat bone. J Biol Chem [Internet]. 1991;266:24684–9. Available from: http://www.ncbi.nlm.nih.gov/pubmed/1722212
Lord DK, Cross NC, Bevilacqua MA, Rider SH, Gorman PA, Groves AV, et al. Type 5 acid phosphatase. Sequence, expression and chromosomal localization of a differentiation-associated protein of the human macrophage. Eur J Biochem [Internet]. 1990; [cited 2015 Oct 18];189:287–93. Available from: http://www.ncbi.nlm.nih.gov/pubmed/2338077
Hayman AR, Dryden AJ, Chambers TJ, Warburton MJ. Tartrate-resistant acid phosphatase from human osteoclastomas is translated as a single polypeptide. Biochem. J. [Internet]. 1991;277(Pt 3):631–4. Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1151288&tool=pmcentrez&rendertype=abstract
Ling P, Roberts RM. Uteroferrin and intracellular tartrate-resistant acid phosphatases are the products of the same gene. J. Biol. Chem. [Internet]. 1993; [cited 2015 Oct 18];268:6896–902. Available from: http://www.ncbi.nlm.nih.gov/pubmed/8463220
Ek-Rylander B, Andersson G. Osteoclast migration on phosphorylated osteopontin is regulated by endogenous tartrate-resistant acid phosphatase. Exp Cell Res. 2010;316:443–51. Elsevier B.V.
An J, Briggs TA, Dumax-Vorzet A, Alarcón-Riquelme ME, Belot A, Beresford M, et al. Tartrate-Resistant Acid Phosphatase Deficiency in the Predisposition to Systemic Lupus Erythematosus. Arthritis Rheumatol. 2016; [cited 2016 Jul 22]; (Hoboken, N.J.) [Internet]. Available from: http://www.ncbi.nlm.nih.gov/pubmed/27390188
Andersson G, Ek-Rylander B, Hollberg K, Ljusberg-SjöLander J, Lång P, Norgård M, et al. TRACP as an Osteopontin Phosphatase. J Bone Miner Res [Internet]. 2003;18:1912–5. [cited 2017 Apr 4]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/14584906
Lausch E, Janecke A, Bros M, Trojandt S, Alanay Y, De Laet C, et al. Genetic deficiency of tartrate-resistant acid phosphatase associated with skeletal dysplasia, cerebral calcifications and autoimmunity. Nat Genet [Internet]. 2011;43:132–7. [cited 2011 Sep 15]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/21217752
Capeller B, Caffier H, Sütterlin MW, Dietl J. Evaluation of tartrate-resistant acid phosphatase (TRAP) 5b as serum marker of bone metastases in human breast cancer. Anticancer Res [Internet]. 2003;23:1011–5. [cited 2015 Oct 15]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/12820340
Chao TY, Wu YY, Janckila AJ. Tartrate-resistant acid phosphatase isoform 5b (TRACP 5b) as a serum maker for cancer with bone metastasis. Clin Chim Acta. 2010;411:1553–64. Elsevier B.V.
Lyubimova NV, Pashkov MV, Tyulyandin SA, Gol’dberg VE, Kushlinskii NE. Tartrate-resistant acid phosphatase as a marker of bone metastases in patients with breast cancer and prostate cancer. Bull Exp Biol Med. 2004;138:77–9.
Tauchert S, di Liberto A, Cordes T, Thill M, Salehin D, Friedrich M. Tartrate-resistant acid phosphatase (TRAP) as a serum marker for bone resorption in breast cancer patients with bone metastases. Clin Exp Obs Gynecol [Internet]. 2009;36:219–25. [cited 2011 Nov 28]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/20101852
Yao NS, Wu YY, Janckila AJ, Ku CH, Hsieh AT, Ho CL, et al. Serum tartrate-resistant acid phosphatase 5b (TRACP5b) activity as a biomarker for bone metastasis in non-small cell lung cancer patients. Clin Chim Acta. 2011;412:181–5. Elsevier B.V.
Adams LM, Warburton MJ, Hayman AR. Human breast cancer cell lines and tissues express tartrate-resistant acid phosphatase (TRAP). Cell Biol Int. 2007;31:191–5.
Honig A, Rieger L, Kapp M, Krockenberger M, Eck M, Dietl J, et al. Increased tartrate-resistant acid phosphatase (TRAP) expression in malignant breast, ovarian and melanoma tissue: an investigational study. BMC Cancer [Internet]. 2006;6:199. [cited 2015 Oct 20]; Available from: http://www.biomedcentral.com/1471-2407/6/199
Zenger S, He W, Ek-Rylander B, Vassiliou D, Wedin R, Bauer H, et al. Differential expression of tartrate-resistant acid phosphatase isoforms 5a and 5b by tumor and stromal cells in human metastatic bone disease. Clin Exp Metastasis [Internet]. 2011;28:65–73. [cited 2011 Sep 8]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/20967488
Scott KL, Nogueira C, Heffernan TP, van Doorn R, Dhakal S, Hanna JA, et al. Proinvasion Metastasis Drivers in Early-Stage Melanoma Are Oncogenes. Cancer Cell [Internet]. 2011;20:92–103. [cited 2015 Oct 22]; Available from: http://www.sciencedirect.com/science/article/pii/S1535610811001954
Xia L, Huang W, Tian D, Chen Z, Zhang L, Li Y, et al. ACP5, a direct transcriptional target of FoxM1, promotes tumor metastasis and indicates poor prognosis in hepatocellular carcinoma. Oncogene [Internet]. 2014;33:1395–406. Available from: http://www.ncbi.nlm.nih.gov/pubmed/23604121
Kawamura M, Tanaka K, Toiyama Y, Okugawa Y, Okigami M, Yasuda H, et al. Clinical significance of tartrate-resistant acid phosphatase type-5 expression in human gastric cancer. Anticancer Res. [Internet]. 2014;34:3425–9. Available from: http://www.ncbi.nlm.nih.gov/pubmed/24982350
Gao Y-L, Liu M-R, Yang S-X, Dong Y-J, Tan X-F. Prognostic significance of ACP5 expression in patients with lung adenocarcinoma. Clin Respir J [Internet]. 2017 [cited 2017 Apr 18]; Available from: http://doi.wiley.com/10.1111/crj.12637
Sheu T-J, Schwarz EM, Martinez DA, O’Keefe RJ, Rosier RN, Zuscik MJ, et al. A phage display technique identifies a novel regulator of cell differentiation. J. Biol. Chem. [Internet]. 2003;278:438–43. American Society for Biochemistry and Molecular Biology; [cited 2015 Oct 23]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/12403789
Metz-Estrella D, Jonason JH, Sheu TJ, Mroczek-Johnston RM, Puzas JE. TRIP-1: A regulator of osteoblast function. J Bone Miner Res [Internet]. 2012;27:1576–84. [cited 2017 Jan 24];Available from: http://www.ncbi.nlm.nih.gov/pubmed/22460930
Patlaka C, Mai HA, Lång P, Andersson G. The growth factor-like adipokine tartrate-resistant acid phosphatase 5a interacts with the rod G3 domain of adipocyte-produced nidogen-2. Biochem Biophys Res Commun. 2014;454:446–52.
Akhurst RJ, Derynck R. TGF-beta signaling in cancer-a double-edged sword. Trends Cell Biol [Internet]. 2001;11:S44–51. [cited 2016 Jun 8]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/11684442
Derynck R, Akhurst RJ, Balmain A. TGF-beta signaling in tumor suppression and cancer progression. Nat. Genet. [Internet]. 2001;29:117–29. Nature Publishing Group; [cited 2016 Jun 8]; Available from: http://www.nature.com/doifinder/10.1038/ng1001-117
de Caestecker MP. Role of Transforming Growth Factor-beta Signaling in Cancer. J Natl Cancer Inst [Internet]. 2000;92:1388–402. [cited 2016 Jun 8]; Available from: http://jnci.oxfordjournals.org/content/92/17/1388.full
Wrana JL, Attisano L, Wieser R, Ventura F, Massagué J. Mechanism of activation of the TGF-beta receptor. [Internet] Nature. 1994:341–7. [cited 2016 Jun 8] Available from: http://www.nature.com/doifinder/10.1038/370341a0
Huang F, Chen Y-G. Regulation of TGF-β receptor activity. Cell Biosci [Internet]. 2012;2:9. [cited 2016 Jun 8]; Available from: http://www.cellandbioscience.com/content/2/1/9
Massagué J. TGFβ Signal Transduction. Annu Rev Biochem [Internet]. 1998;67:753–91. [cited 2016 Jun 8]; Available from: http://www.annualreviews.org/doi/abs/10.1146/annurev.biochem.67.1.753
Morikawa M, Derynck R, Miyazono K. TGF-β and the TGF-β Family: Context-Dependent Roles in Cell and Tissue Physiology. Cold Spring Harb Perspect Biol [Internet]. 2016;8:a021873. [cited 2016 May 17]; Available from: http://cshperspectives.cshlp.org/content/8/5/a021873.long
Wrighton KH, Lin X, Feng X-H. Phospho-control of TGF-beta superfamily signaling. Cell Res [Internet]. 2009;19:8–20. [cited 2016 Jun 8]; Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=2929013&tool=pmcentrez&rendertype=abstract
MacCallum J, Bartlett JMS, Thompson AM, Keen JC, Dixon JM, Miller WR. Expression of transforming growth factor beta mRNA isoforms in human breast cancer. Br J Cancer [Internet]. 1994;69:1006–9. [cited 2016 Jun 10]; Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1969445/?tool=pmcentrez&report=abstract
Walker RA, Dearing SJ. Transforming growth factor beta 1 in ductal carcinoma in situ and invasive carcinomas of the breast. Eur J Cancer [Internet]. 1992;28:641–4. [cited 2016 Jun 10]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/1317202
Gorsch SM, Memoli VA, Stukel TA, Gold LI, Arrick BA. Immunohistochemical staining for transforming growth factor beta 1 associates with disease progression in human breast cancer. Cancer Res [Internet]. 1992;52:6949–52. [cited 2016 Jun 10];Available from: http://www.ncbi.nlm.nih.gov/pubmed/1458485
Dave H, Trivedi S, Shah M, Shukla S. Transforming growth factor beta 2: a predictive marker for breast cancer. Indian J Exp Biol [Internet]. 2011;49:879–87. [cited 2016 Nov 25]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/22126020
Rodón L, Gonzàlez-Juncà A, Del Mar IM, Sala-Hojman A, Martínez-Sáez E, Seoane J. Active CREB1 promotes a malignant TGFb2 autocrine loop in glioblastoma. Cancer Discov [Internet]. 2014;4:1230–41. [cited 2016 Nov 4];Available from: http://cancerdiscovery.aacrjournals.org/cgi/doi/10.1158/2159-8290.CD-14-0275
Bourguignon LYW, Singleton PA, Zhu H, Zhou B. Hyaluronan promotes signaling interaction between CD44 and the transforming growth factor beta receptor I in metastatic breast tumor cells. J. Biol. Chem. [Internet]. 2002;277:39703–12. [cited 2017 Mar 2]; Available from: http://www.jbc.org/cgi/content/short/277/42/39703
Weber GF, Ashkar S, Glimcher MJ, Cantor H. Receptor-Ligand Interaction Between CD44 and Osteopontin (Eta-1). Science (80-. ). [Internet]. 1996;271:509–12. Available from: http://science.sciencemag.org/content/271/5248/509
Ek-Rylander B, Flores M, Wendel M, Heinegard D, Andersson G. Dephosphorylation of osteopontin and bone sialoprotein by osteoclastic tartrate-resistant acid phosphatase. Modulation of osteoclast adhesion in vitro. J. Biol. Chem. [Internet]. 1994;269:14853–6. Available from: http://www.ncbi.nlm.nih.gov/pubmed/8195113
Franzén A, Heinegård D. Isolation and characterization of two sialoproteins present only in bone calcified matrix. Biochem. J. [Internet]. 1985;232:715–24. [cited 2015 Oct 15]; Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1152943&tool=pmcentrez&rendertype=abstract
Ek-Rylander B, Barkhem T, Ljusberg J, Ohman L, Andersson KK, Andersson G. Comparative studies of rat recombinant purple acid phosphatase and bone tartrate-resistant acid phosphatase. Biochem. J. [Internet]. 1997;321(2):305–11. Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC1218069/?tool=pmcentrez
Krumpel M, Reithmeier A, Senge T, Baeumler TA, Frank M, Nyholm PG, et al. The small chemical enzyme inhibitor 5-phenylnicotinic acid/CD13 inhibits cell migration and invasion of tartrate-resistant acid phosphatase/ACP5-overexpressing MDA-MB-231 breast cancer cells. Exp Cell Res [Internet]. 2015;339:154–62. [cited 2015 Oct 13]; Available from: http://www.sciencedirect.com/science/article/pii/S0014482715301014
Zenger S, Ek-Rylander B, Andersson G. Biogenesis of tartrate-resistant acid phosphatase isoforms 5a and 5b in stably transfected MDA-MB-231 breast cancer epithelial cells. Biochim Biophys Acta - Mol Cell Res. 2010;1803:598–607. Elsevier B.V.;
Branca RMM, Orre LM, Johansson HJ, Granholm V, Huss M, Pérez-Bercoff Å, et al. HiRIEF LC-MS enables deep proteome coverage and unbiased proteogenomics. Nat Methods [Internet]. 2014;11:59–62. Available from: http://www.nature.com/doifinder/10.1038/nmeth.2732
Muller JW. Possible advantages of a robust evaluation of comparisons. J Res Natl Inst Stand Technol [Internet]. 2000;105:551. [cited 2016 Dec 12]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/27551622
Hsu PP, Kang SA, Rameseder J, Zhang Y, Ottina KA, Lim D, et al. The mTOR-Regulated Phosphoproteome Reveals a Mechanism of mTORC1-Mediated Inhibition of Growth Factor Signaling. Science (80-. ). [Internet]. 2011;332:1317–22. [cited 2016 Dec 12]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/21659604
Manning G, Whyte DB, Martinez R, Hunter T, Sudarsanam S. The Protein Kinase Complement of the Human Genome. Science (80-. ). [Internet]. 2002;298:1912–34. Available from: http://www.sciencemag.org/content/298/5600/1912.abstract
Duan G, Li X, Köhn M. The human DEPhOsphorylation database DEPOD: A 2015 update. Nucleic Acids Res. 2015;43:D531–5.
Barrera LA, Vedenko A, Kurland JV, Rogers JM, Gisselbrecht SS, Rossin EJ, et al. Survey of variation in human transcription factors reveals prevalent DNA binding changes. Science (80-. ). [Internet]. 2016;351:1450–4. [cited 2017 Apr 4]; Available from: http://science.sciencemag.org/content/351/6280/1450
Chen H-F, Wang J-K. The databases of transcription factors. Yi Chuan. 2010;32:1009–17.
Vaquerizas JM, Kummerfeld SK, Teichmann SA, Luscombe NM. A census of human transcription factors: function, expression and evolution. Nat Rev Genet. 2009;10:252–63.
Gao T, Liu Z, Wang Y, Cheng H, Yang Q, Guo A, et al. UUCD: A family-based database of ubiquitin and ubiquitin-like conjugation. Nucleic Acids Res. [Internet]. 2013;41. Available from: http://uucd.biocuckoo.org/
Hornbeck PV, Zhang B, Murray B, Kornhauser JM, Latham V, Skrzypek E. PhosphoSitePlus, 2014: Mutations, PTMs and recalibrations. Nucleic Acids Res. 2015;43:D512–20.
Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res [Internet]. 2015;43:D447–52. [cited 2016 Aug 3]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/25352553
Panizza E, Branca RMM, Oliviusson P, Orre LM, Lehtiö J. Isoelectric point-based fractionation by HiRIEF coupled to LC-MS allows for in-depth quantitative analysis of the phosphoproteome. Sci Rep. 2017;7:4513.
Feder D, Hussein WM, Clayton DJ, Kan MW, Schenk G, McGeary RP, et al. Identification of Purple Acid Phosphatase Inhibitors by Fragment-Based Screening: Promising New Leads for Osteoporosis Therapeutics. Chem Biol Drug Des. 2012;80:665–74.
Bazer FW, Worthington-White D, Fliss MF, Gross S. Uteroferrin: a progesterone-induced hematopoietic growth factor of uterine origin. Exp Hematol [Internet]. 1991;19:910–5. [cited 2015 Oct 15]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/1893968
Simpson KJ, Selfors LM, Bui J, Reynolds A, Leake D, Khvorova A, et al. Identification of genes that regulate epithelial cell migration using an siRNA screening approach. Nat Cell Biol [Internet]. 2008;10:1027–38. [cited 2011 Jul 5]; Available from: http://www.ncbi.nlm.nih.gov/entrez/query.fcgi?cmd=Retrieve&db=PubMed&dopt=Citation&list_uids=19160483
Briggs TA, Rice GI, Daly S, Urquhart J, Gornall H, Bader-Meunier B, et al. Tartrate-resistant acid phosphatase deficiency causes a bone dysplasia with autoimmunity and a type I interferon expression signature. Nat Genet [Internet]. 2011;43:127–31. [cited 2017 Sept 4]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/21217755
Chellaiah M, Hruska K. The integrin alpha(v)beta(3) and CD44 regulate the actions of osteopontin on osteoclast motility. Calcif Tissue Int [Internet]. 2003;72:197–205. [cited 2011 Oct 4]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/12469249
Böttinger EP, Jakubczak JL, Roberts ISD, Mumy M, Hemmati P, Bagnall K, et al. Expression of a dominant-negative mutant TGF-β type II receptor in transgenic mice reveals essential roles for TGF-β in regulation of growth and differentiation in the exocrine pancreas. EMBO J [Internet]. 1997;16:2621–33. EMBO Press; [cited 2016 Jun 9];Available from: http://emboj.embopress.org/cgi/doi/10.1093/emboj/16.10.2621
Halder SK, Beauchamp RD, Datta PK. A specific inhibitor of TGF-beta receptor kinase, SB-431542, as a potent antitumor agent for human cancers. Neoplasia [Internet]. 2005;7:509–21. [cited 2016 Mar 8]; Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=1501161&tool=pmcentrez&rendertype=abstract
Liu X, Sun Y, Constantinescu SN, Karam E, Weinberg RA, Lodish HF. Transforming growth factor beta-induced phosphorylation of Smad3 is required for growth inhibition and transcriptional induction in epithelial cells. Proc Natl Acad Sci U S A [Internet]. 1997;94:10669–74. [cited 2016 Jun 8]; Available from: https://www.ncbi.nlm.nih.gov/pmc/articles/PMC23442/
Morton DM, Barrack ER. Modulation of transforming growth factor beta 1 effects on prostate cancer cell proliferation by growth factors and extracellular matrix. Cancer Res [Internet]. 1995;55:2596–602. [cited 2016 Jun 8]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/7780974
Farina AR, Coppa A, Tiberio A, Tacconelli A, Turco A, Colletta G, et al. Transforming growth factor-beta1 enhances the invasiveness of human MDA- MB-231 breast cancer cells by up-regulating urokinase activity. Int J Cancer. Wiley Subscription Services, Inc., A Wiley Company; 1998;75:721–30.
Oft M, Heider K-H, Beug H. TGFβ signaling is necessary for carcinoma cell invasiveness and metastasis. Curr Biol [Internet]. 1998;8:1243–52. [cited 2016 Jun 8]; Available from: http://www.cell.com/article/S0960982207005337/fulltext
Mima K, Okabe H, Ishimoto T, Hayashi H, Nakagawa S, Kuroki H, et al. CD44s regulates the TGFbeta-mediated mesenchymal phenotype and is associated with poor prognosis in patients with hepatocellular carcinoma. Cancer Res. [Internet]. 2012;72:3414–23. American Association for Cancer Research; [cited 2016 Jul 25]; Available from: http://www.ncbi.nlm.nih.gov/pubmed/22552294
Ouhtit A, Madani S, Gupta I, Shanmuganathan S, Abdraboh ME, Al-Riyami H, et al. TGFbeta2: A novel target of CD44-promoted breast cancer invasion. J Cancer [Internet]. 2013;4:566–72. [cited 2016 Apr 7]; Available from: http://www.pubmedcentral.nih.gov/articlerender.fcgi?artid=3753531&tool=pmcentrez&rendertype=abstract
Zhang Y, Tripathi V, Sixt K, Heller M. TGF-beta Regulates Alternative Splicing of CD44 by Inducing Smad3 Binding to CD44 pre-mRNA. FASEB J [Internet]. 2015;29:562.15. [cited 2016 Apr 7]; Available from: http://www.fasebj.org/cgi/content/long/29/1_Supplement/562.15
The authors would like to thank Shuvojit Moulik, Melanie Weiss and Joseph Azem for their assistance in laboratory preparations.
This work was funded by the Swedish Research Council and supported by grants supplied from the Swedish Cancer Society, the Cancer Society Stockholm and the Swedish Foundation for Strategic Research. Funding bodies did not have any influence in the design of the study and collection, analysis and interpretation of data or in writing the manuscript.
Availability of data and materials
The datasets generated and/or analyzed during the current study are available in the ProteomeXchange/ PRIDE database: “Proteomic and phosphoproteomics profiling of breast cancer cells with altered expression level of the protein phosphatase TRAP”, dataset identifier PXD006430 (https://www.ebi.ac.uk/pride/archive/login). Other raw data generated and/or analyzed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Consent for publication
The authors declare that they do not have any competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional explanations on Material and Methods including descriptions of gene expression analysis, cell lysis, SILAC labeling, protein extraction for mass spectrometric analyses, protein digestion, Tandem Mass Tag labeling, liquid chromatography tandem mass spectrometric analyses and proteomics database search, protein and phospho-peptide ratios calculation. (DOCX 32 kb)
Distribution of unique peptides across fractions generated by HiRIEF separation on IPG 2.5–3.7 strips. Number of unique peptides and phospho-peptides identified across fractions in each biological replicate (A) and (B), displayed by number of phosphorylations. Fraction numbering proceeds from the acidic end towards the basic end of the strips. Multiply phosphorylated peptides are identified prevalently in the first 30 fractions in the IPG 2.5–3.7 strip. (PDF 148 kb)
Quantitative reproducibility of biological replicates employed for phosphoproteomics and proteomics analysis of TRAP3high (Heavy SILAC labeled) and control (Light SILAC labeled) MDA-MB-231 cells. Correlation of log2 transformed ratio (H/L) values for replicate pairs for phosphoproteomics analysis (A) and standard proteomics analysis (B); Pearson correlation coefficient is displayed. (PDF 2872 kb)
Phosphorylation sites significantly regulated in TRAP3high cells compared to control (ctrl) cells. Unique phosphorylation sites are defined by a sequence window of 15 amino acids centered at the phosphorylated residue. The column “Protein class” describes to which of the examined protein class a protein belongs to. Examined proteins classes include protein kinases, protein phosphatases, ubiquitin and UBL system enzymes and transcription factors. The column “Functional site” describes whether the phospho-site has a known regulatory function based on the information available on the database PhosphoSitePlus. (XLSX 48 kb)
Volcano plots representing the selection of proteins significantly regulated upon TRAP knockdown in TRAP3high MDA-MB-231 cells. The average log2 transformed ratio of the three replicates is plotted for each experimental condition. Significantly regulated events are defined by log2 transformed ratios at least −/+3 MAD away from the median and p-values lower than 0.01. TRAP sh2 samples versus scrambled (A). TRAP sh3 + 4 samples vs scrambled (B). (PDF 5171 kb)
GO enrichment analyses of genes corresponding to phosphorylation sites significantly regulated upon TRAP overexpression in MDA-MB-231 cells. GO terms enrichment was evaluated for the set of genes corresponding to phospho-sites significantly regulated in TRAP3high cells (99 genes) compared to the set of genes identified across all proteomics analyses (9570 genes). Enriched GO terms are displayed by significance of the enrichment (q-value, multiple testing corrected) and fold enrichment. GO Molecular Function enrichment analysis (A). GO Cell Component enrichment analysis (B). (PDF 86 kb)
Proteins significantly regulated upon TRAP perturbation in any of the phosphoproteomics or proteomics analyses (TRAP3high vs ctrl, phospho and standard proteomics analyses; TRAP sh2 or sh3 + 4 vs scrambled shRNA, standard proteomics analysis). The column “Enriched GO term” specifies whether the protein belong to any of the GO terms that were found to be enriched in the list of regulated proteins. The column “Protein class” describes to which of the examined protein class a protein belongs to. Examined proteins classes include protein kinases, protein phosphatases, ubiquitin and UBL system enzymes and transcription factors. (XLSX 79 kb)
Functional analysis after treatment with TGFβ1 (10 μg/mL, A, B, n = 5) and respective controls. Live cell imaging was performed over 30 h in serum-free medium containing blocking antibody or small molecule inhibitor and wound migration velocity (A) and migration curves are compared. Cell proliferation was assessed after 24 h treatment in complete medium (B). Statistical comparison was performed on biological replicates by ANOVA test. Groups are generally compared to their respective controls (ctrl, untreated) or indicated with brackets and significance is annotated with an asterisk (*). ns = non-significant. “n=” indicates the number of replicates. (PDF 952 kb)